EAK2000 Φάσματα

[1]:
from streng.codes.greek.eak2000.raw.ch2.seismic_action import spectra
import numpy as np
import matplotlib
import matplotlib.pyplot as plt
[2]:
ground_type = 'Β'
[3]:
T1 = spectra.T1(ground_type)
T2 = spectra.T2(ground_type)
print(f'T1={T1}')
print(f'T2={T2}')
T1=0.15
T2=0.6
[4]:
x = np.linspace(1e-10, 2, 201)
y = spectra.Φd(T=x,
               α=0.16,
               γI=1.0,
               T1=T1,
               T2=T2,
               q=3.5,
               η=1.0,
               θ=1.0,
               β0=2.5)

plt.plot(x, y, lw=2, label='EAK spectrum')
plt.axis([0, 2, 0, 1.2 * y.max()])
plt.ylabel('$S_{e}$ (g)')
plt.xlabel('$T$ (sec)') #'μθpl'
plt.title('EEAK spectrum')
plt.legend()
plt.show()
../../../_images/jupyters_codes_greek_eak_spectrum_4_0.png
[5]:
print(f'{x[3]:.2f}')
0.03