pygmm.tavakoli_pezeshk_2005.TavakoliPezeshk05¶

class pygmm.tavakoli_pezeshk_2005.TavakoliPezeshk05(scenario)[source]¶

Tavakoli and Pezeshk (2005, [Tavakoli and Pezeshk, 2005]) model.

Developed for the Eastern North America with a reference velocity of 2880 m/s.

Parameters:: scenario (pygmm.model.Scenario) – earthquake scenario

NAME = 'Tavakoli and Pezeshk (2005)'¶: Long name of the model

ABBREV = 'TP05'¶: Short name of the model

V_REF = 2880.0¶

COEFF = rec.array([(0. , 1.14 , 0.623, -0.0483, -1.81, -0.652, 0.446, -2.93e-05, -0.00405, 0.00946 , 1.41, -0.961, 0.000432, 0.000133, 1.21, -0.111 , 0.409), (0.05, 1.82 , 0.533, -0.0475, -1.63, -0.567, 0.454, 7.77e-03, -0.00491, -0.00314 , 0.98, -0.939, 0.000512, 0.00093 , 1.22, -0.108 , 0.441), (0.08, 0.683, 0.743, -0.0293, -1.71, -0.756, 0.46 , -9.68e-04, -0.00494, -0.0055 , 1.13, -0.916, 0.000482, 0.000733, 1.22, -0.108 , 0.449), (0.1 , 0.869, 0.607, -0.0474, -1.52, -0.704, 0.449, -6.19e-03, -0.0047 , -0.00424 , 1.04, -0.913, 0.000411, 0.000358, 1.23, -0.108 , 0.456), (0.15, 2.38 , 0.501, -0.0642, -1.73, -0.976, 0.414, 6.60e-03, -0.0048 , 0.00393 , 1.51, -0.865, 0.000364, 0.000684, 1.24, -0.108 , 0.464), (0.2 , -0.548, 0.857, -0.0262, -1.68, -0.861, 0.433, 2.79e-03, -0.00365, -0.00202 , 1.64, -0.925, 0.000161, 0.000643, 1.24, -0.108 , 0.469), (0.3 , -0.513, 0.667, -0.0443, -1.42, -0.47 , 0.468, 1.08e-02, -0.00541, 0.00644 , 1.52, -0.915, 0.000432, 0.000287, 1.26, -0.109 , 0.479), (0.5 , 0.24 , 0.611, -0.0789, -1.55, -0.844, 0.414, 7.89e-03, -0.00365, -0.000265, 1.59, -0.859, 0.000277, 0.000146, 1.28, -0.173 , 0.505), (0.75, -0.679, 0.666, -0.083 , -1.48, -0.734, 0.435, 9.53e-03, -0.00337, -0.00119 , 1.55, -0.784, 0.000245, 0.000547, 1.28, -0.105 , 0.522), (1. , -1.55 , 0.764, -0.0859, -1.49, -0.941, 0.424, -5.84e-03, -0.00209, 0.0033 , 1.52, -0.757, 0.000117, 0.000759, 1.28, -0.103 , 0.537), (1.5 , -2.3 , 0.794, -0.0884, -1.45, -0.886, 0.412, 8.30e-03, -0.00327, 0.00251 , 1.71, -0.769, 0.000233, 0.000166, 1.27, -0.0999, 0.551), (2. , -2.7 , 0.805, -0.0929, -1.44, -0.923, 0.408, 2.06e-02, -0.00214, 0.0023 , 1.43, -0.755, 0.000214, 0.000391, 1.26, -0.0978, 0.562), (3. , -2.42 , 0.801, -0.108 , -1.65, -0.898, 0.437, 1.67e-02, -0.00203, 0.00358 , 1.93, -0.818, 0.000116, 0.000398, 1.26, -0.0952, 0.573), (4. , -3.69 , 0.817, -0.118 , -1.46, -0.845, 0.425, 1.13e-02, -0.00172, -0.00334 , 1.69, -0.737, 0.00011 , 0.000359, 1.25, -0.0926, 0.589)], dtype=[('period', '<f8'), ('c_1', '<f8'), ('c_2', '<f8'), ('c_3', '<f8'), ('c_4', '<f8'), ('c_5', '<f8'), ('c_6', '<f8'), ('c_7', '<f8'), ('c_8', '<f8'), ('c_9', '<f8'), ('c_10', '<f8'), ('c_11', '<f8'), ('c_12', '<f8'), ('c_13', '<f8'), ('c_14', '<f8'), ('c_15', '<f8'), ('c_16', '<f8')])¶

PERIODS = array([0. , 0.05, 0.08, 0.1 , 0.15, 0.2 , 0.3 , 0.5 , 0.75, 1. , 1.5 , 2. , 3. , 4. ])¶: Indices of the periods

INDEX_PGA = 0¶: Index of the peak ground acceleration

INDICES_PSA = array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13])¶: Indices for the spectral accelerations

PARAMS = [<pygmm.model.NumericParameter object>, <pygmm.model.NumericParameter object>]¶: Model parameters

__init__(scenario)[source]¶

Initialize the model.

Parameters:: scenario (Scenario)

INDEX_PGD = None¶: Index of the peak ground displacement

INDEX_PGV = None¶: Index of the peak ground velocity

LIMITS = {}¶: Limits of model applicability

PGD_SCALE = 1.0¶: Scale factor to apply to get PGD in cm

PGV_SCALE = 1.0¶: Scale factor to apply to get PGV in cm/sec

interp_ln_spec_accels(periods, kind='linear')¶

Interpolate the spectral acceleration.

Interpolation of the spectral acceleration is done in natural log space.

Parameters:

periods (array_like) – spectral periods to interpolate the response.
kind (str, optional) – see scipy.interpolate.interp1d() for description of kind. Options include: ‘linear’ (default), ‘nearest’, ‘zero’, ‘slinear’, ‘quadratic’, and ‘cubic’

Returns:

ln_spec_accels – interpolated spectral accelerations

Return type:

np.ndarray

interp_ln_stds(periods, kind='linear')¶

Interpolate the logarithmic standard deviation.

Interpolate the logarithmic standard deviation (\(\sigma_{\ln}\)) of spectral acceleration at the provided damping at specified periods.

Parameters:

periods (array_like) – spectral periods to interpolate the response.
kind (str, optional) – see scipy.interpolate.interp1d() for description of kind. Options include: ‘linear’ (default), ‘nearest’, ‘zero’, ‘slinear’, ‘quadratic’, and ‘cubic’

Returns:

ln_stds – interpolated logarithmic standard deviations

Return type:

np.ndarray

interp_spec_accels(periods, kind='linear')¶

Interpolate the spectral acceleration.

Interpolation of the spectral acceleration is done in natural log space.

Parameters:

periods (array_like) – spectral periods to interpolate the response.
kind (str, optional) – see scipy.interpolate.interp1d() for description of kind. Options include: ‘linear’ (default), ‘nearest’, ‘zero’, ‘slinear’, ‘quadratic’, and ‘cubic’

Returns:

spec_accels – interpolated spectral accelerations

Return type:

np.ndarray

property ln_std_pga: float¶: Peak ground accelaration log-standard deviation.

property ln_std_pgd: float¶: Peak ground displacement log-standard deviation.

property ln_std_pgv: float¶: Peak ground velocity log-standard deviation.

property ln_stds: ndarray¶: Pseudo-spectral accelerations log-standard deviation.

property periods: ndarray¶: Periods specified by the model.

property pga: float¶: Peak ground acceleration (PGA) computed by the model (g).

property pgd: float¶: Peak ground displacement (PGD) computed by the model (cm).

property pgv: float¶: Peak ground velocity (PGV) computed by the model (cm/sec).

property scenario¶

property spec_accels: ndarray¶: Pseudo-spectral accelerations computed by the model (g).