3.5. User-Defined Average Spectral Acceleration

imcalculator.get_saavg_user_defined(periods_list)[source]

Computes the geometric mean of spectral accelerations for a user-defined list of periods.

Parameters:: periods_list (list or numpy.ndarray) – List of user-defined periods (s) at which spectral accelerations are computed.
Returns:: sa_avg – Geometric mean of spectral accelerations (g) over the user-defined periods.
Return type:: float

Example

import numpy as np
from openquake.vmtk.imcalculator import imcalculator

acc = np.loadtxt("openquake/vmtk/tests/test_data/acceleration.txt")
im = imcalculator(acc, dt=0.005)

periods = np.linspace(0.1, 1.0, 10)
avg_sa = im.get_saavg_user_defined(periods=periods)
print(f"User-defined AvgSa = {avg_sa:.4f} g")