pygmm.campbell_2003.Campbell2003¶

class pygmm.campbell_2003.Campbell2003(scenario)[source]¶

Campbell (2003, [Campbell, 2003]) model.

This model was developed for the Eastern US.

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

NAME = 'Campbell (2003)'¶: Long name of the model

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

V_REF = 2800.0¶

COEFF = rec.array([(0.01 , 0.0305, 0.633, -0.0427, -1.591, -0.00428, 4.83e-04, 0.683, 0.416, 1.14 , -0.873, 1.03 , -0.086 , 0.414), (0.02 , 1.3535, 0.63 , -0.0404, -1.787, -0.00388, 4.97e-04, 1.02 , 0.363, 0.851, -0.715, 1.03 , -0.086 , 0.414), (0.03 , 1.186 , 0.622, -0.0362, -1.691, -0.00367, 5.01e-04, 0.922, 0.376, 0.759, -0.922, 1.03 , -0.086 , 0.414), (0.05 , 0.3736, 0.616, -0.0353, -1.469, -0.00378, 5.00e-04, 0.63 , 0.423, 0.771, -1.239, 1.042, -0.0838, 0.443), (0.075, -0.0395, 0.615, -0.0353, -1.383, -0.00421, 4.86e-04, 0.491, 0.463, 0.955, -1.349, 1.052, -0.0838, 0.453), (0.1 , -0.1475, 0.613, -0.0353, -1.369, -0.00454, 4.60e-04, 0.484, 0.467, 1.096, -1.284, 1.059, -0.0838, 0.46 ), (0.15 , -0.1901, 0.616, -0.0478, -1.368, -0.00473, 3.93e-04, 0.461, 0.478, 1.239, -1.079, 1.068, -0.0838, 0.469), (0.2 , -0.4328, 0.617, -0.0586, -1.32 , -0.0046 , 3.37e-04, 0.399, 0.493, 1.25 , -0.928, 1.077, -0.0838, 0.478), (0.3 , -0.6906, 0.609, -0.0786, -1.28 , -0.00414, 2.63e-04, 0.349, 0.502, 1.241, -0.753, 1.081, -0.0838, 0.482), (0.5 , -0.5907, 0.534, -0.1379, -1.216, -0.00341, 1.94e-04, 0.318, 0.503, 1.166, -0.606, 1.098, -0.0824, 0.508), (0.75 , -0.5429, 0.48 , -0.1806, -1.184, -0.00288, 1.60e-04, 0.304, 0.504, 1.11 , -0.526, 1.105, -0.0806, 0.528), (1. , -0.6104, 0.451, -0.209 , -1.158, -0.00255, 1.41e-04, 0.299, 0.503, 1.067, -0.482, 1.11 , -0.0793, 0.543), (1.5 , -0.9666, 0.441, -0.2405, -1.135, -0.00213, 1.19e-04, 0.304, 0.5 , 1.029, -0.438, 1.099, -0.0771, 0.547), (2. , -1.4306, 0.459, -0.2552, -1.124, -0.00187, 1.03e-04, 0.31 , 0.499, 1.015, -0.417, 1.093, -0.0758, 0.551), (3. , -2.2331, 0.492, -0.2646, -1.121, -0.00154, 8.40e-05, 0.31 , 0.499, 1.014, -0.393, 1.09 , -0.0737, 0.562), (4. , -2.7975, 0.507, -0.2738, -1.119, -0.00135, 7.40e-05, 0.294, 0.506, 1.018, -0.386, 1.092, -0.0722, 0.575)], 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')])¶

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

INDICES_PSA = array([ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15])¶: 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_PGA = None¶: Index of the peak ground acceleration

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).