Διαγράμματα ελέγχου FEMA440 + ATC40

In [12]:
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sns
sns.set_style("whitegrid")

from streng.codes.usa.fema440.raw.ch3 import csm as fema440_csm
from streng.codes.usa.atc40.raw.ch8 import csm as atc40_csm

Τeq vs μ (FEMA440 Figure 3-3)

In [2]:
μs = np.linspace(0,10, 101)
plt.plot(μs, [fema440_csm.Teq(1.0, x, 0.00) for x in μs], label = 'α = 0.00')
plt.plot(μs, [fema440_csm.Teq(1.0, x, 0.03) for x in μs], label = 'α = 0.03')
plt.plot(μs, [fema440_csm.Teq(1.0, x, 0.05) for x in μs], label = 'α = 0.05')
plt.ylabel('$T_{eq}/T_0$')
plt.xlabel('$μ$') #'μθpl'
plt.title('$Τ_{eq}$ vs $μ$')
plt.legend()
plt.axis([0.,10.0,0.,4.])
plt.show()
../../_images/jupyters_tools_capacity_spectrum_method_charts_3_0.png

κ vs β0 (ATC40 Figure 8-15)

In [3]:
β0s = np.linspace(0, 0.6, 101)
plt.plot(β0s, [atc40_csm.κ(x, 'A') for x in β0s], label = 'Type A')
plt.plot(β0s, [atc40_csm.κ(x, 'B') for x in β0s], label = 'Type B')
plt.plot(β0s, [atc40_csm.κ(x, 'C') for x in β0s], label = 'Type C')
plt.ylabel('$κ$')
plt.xlabel('$β_0$') #'μθpl'
plt.title('$κ$ vs $β_0$')
plt.legend()
plt.axis([0.,0.6,0.,1.1])
plt.show()
../../_images/jupyters_tools_capacity_spectrum_method_charts_5_0.png

κ vs μ (FEMA440 Figure 3-4)

In [4]:
μs = np.linspace(0.,10., 101)
β0s = [fema440_csm.β0(μ, 0.0) for μ in μs]
plt.plot(μs, [atc40_csm.κ(x, 'A') for x in β0s], label = 'Type A')
plt.plot(μs, [atc40_csm.κ(x, 'B') for x in β0s], label = 'Type B')
plt.plot(μs, [atc40_csm.κ(x, 'C') for x in β0s], label = 'Type C')
plt.ylabel('$κ$')
plt.xlabel('$μ$') #'μθpl'
plt.title('$κ$ vs $μ$')
plt.legend()
plt.axis([0.,10,0.,1.2])
plt.show()
../../_images/jupyters_tools_capacity_spectrum_method_charts_7_0.png

β0 vs μ

In [5]:
μs = np.linspace(0.,10, 101)
β0s = [fema440_csm.β0(μ, 0.0) for μ in μs]

plt.plot(μs, [fema440_csm.β0(μ, 0.0) for μ in μs], label = 'α = 0.00')
plt.plot(μs, [fema440_csm.β0(μ, 0.3) for μ in μs], label = 'α = 0.03')
plt.plot(μs, [fema440_csm.β0(μ, 0.5) for μ in μs], label = 'α = 0.05')

plt.ylabel('$β_{0}$')
plt.xlabel('$μ$') #'μθpl'
plt.title('$β_{0}$ vs $μ$')
plt.legend()
plt.axis([0.,10.0,0.,1.])
plt.show()
../../_images/jupyters_tools_capacity_spectrum_method_charts_9_0.png

βeff vs β0 (ATC40 Figure 8-16)

In [6]:
β0s = np.linspace(0, 0.6, 101)
plt.plot(β0s, [atc40_csm.βeff(0.05, x, 'A') for x in β0s], label = 'Type A')
plt.plot(β0s, [atc40_csm.βeff(0.05, x, 'B') for x in β0s], label = 'Type B')
plt.plot(β0s, [atc40_csm.βeff(0.05, x, 'C') for x in β0s], label = 'Type C')
plt.ylabel('$β_{eff}$')
plt.xlabel('$β_0$') #'μθpl'
plt.title('$β_{eff}$ vs $β_0$')
plt.legend()
plt.axis([0.,0.6,0.,0.5])
plt.show()
../../_images/jupyters_tools_capacity_spectrum_method_charts_11_0.png

βeff vs μ (FEMA440 Figure 3-5)

In [7]:
μs = np.linspace(0.,10, 101)
β0s = [fema440_csm.β0(μ, 0.0) for μ in μs]

plt.plot(μs, [atc40_csm.βeff(0.05, x, 'A') for x in β0s], label = 'Type A')
plt.plot(μs, [atc40_csm.βeff(0.05, x, 'B') for x in β0s], label = 'Type B')
plt.plot(μs, [atc40_csm.βeff(0.05, x, 'C') for x in β0s], label = 'Type C')

plt.ylabel('$β_{eff}$')
plt.xlabel('$μ$') #'μθpl'
plt.title('$β_{eff}$ vs $μ$')
plt.legend()
plt.axis([0.,10.0,0.,0.5])
plt.show()
../../_images/jupyters_tools_capacity_spectrum_method_charts_13_0.png

SRA vs β0 (ATC40 Figure 8-17)

In [8]:
β0s = np.linspace(0, 0.6, 101)
plt.plot(β0s, [atc40_csm.SRA(atc40_csm.βeff(0.05, x, 'A'),'A')  for x in β0s], label = 'Type A')
plt.plot(β0s, [atc40_csm.SRA(atc40_csm.βeff(0.05, x, 'B'),'B')  for x in β0s], label = 'Type B')
plt.plot(β0s, [atc40_csm.SRA(atc40_csm.βeff(0.05, x, 'C'),'C')  for x in β0s], label = 'Type C')
plt.ylabel('$SR_{A}$')
plt.xlabel('$β_0$') #'μθpl'
plt.title('$SR_{A}$ vs $β_0$')
plt.legend()
plt.axis([0.,0.6,0.,1.0])
plt.show()
../../_images/jupyters_tools_capacity_spectrum_method_charts_15_0.png

SRV vs β0 (ATC40 Figure 8-18)

In [9]:
β0s = np.linspace(0, 0.6, 101)
plt.plot(β0s, [atc40_csm.SRV(atc40_csm.βeff(0.05, x, 'A'),'A')  for x in β0s], label = 'Type A')
plt.plot(β0s, [atc40_csm.SRV(atc40_csm.βeff(0.05, x, 'B'),'B')  for x in β0s], label = 'Type B')
plt.plot(β0s, [atc40_csm.SRV(atc40_csm.βeff(0.05, x, 'C'),'C')  for x in β0s], label = 'Type C')
plt.ylabel('$SR_{V}$')
plt.xlabel('$β_0$') #'μθpl'
plt.title('$SR_{V}$ vs $β_0$')
plt.legend()
plt.axis([0.,0.6,0.,1.0])
plt.show()
../../_images/jupyters_tools_capacity_spectrum_method_charts_17_0.png

SRA vs μ (FEMA440 Figure 3-6)

In [10]:
μs = np.linspace(0.,10, 101)
β0s = [fema440_csm.β0(μ, 0.0) for μ in μs]

plt.plot(μs, [atc40_csm.SRA(atc40_csm.βeff(0.05, x, 'A'),'A')  for x in β0s], label = 'Type A')
plt.plot(μs, [atc40_csm.SRA(atc40_csm.βeff(0.05, x, 'B'),'B')  for x in β0s], label = 'Type B')
plt.plot(μs, [atc40_csm.SRA(atc40_csm.βeff(0.05, x, 'C'),'C')  for x in β0s], label = 'Type C')

plt.ylabel('$SR_{A}$')
plt.xlabel('$μ$') #'μθpl'
plt.title('$SR_{A}$ vs $μ$')
plt.legend()
plt.axis([0.,10.0,0.,1.2])
plt.show()
../../_images/jupyters_tools_capacity_spectrum_method_charts_19_0.png

SRV vs μ (FEMA440 Figure 3-7)

In [11]:
μs = np.linspace(0.,10, 101)
β0s = [fema440_csm.β0(μ, 0.0) for μ in μs]

plt.plot(μs, [atc40_csm.SRV(atc40_csm.βeff(0.05, x, 'A'),'A')  for x in β0s], label = 'Type A')
plt.plot(μs, [atc40_csm.SRV(atc40_csm.βeff(0.05, x, 'B'),'B')  for x in β0s], label = 'Type B')
plt.plot(μs, [atc40_csm.SRV(atc40_csm.βeff(0.05, x, 'C'),'C')  for x in β0s], label = 'Type C')

plt.ylabel('$SR_{V}$')
plt.xlabel('$μ$') #'μθpl'
plt.title('$SR_{V}$ vs $μ$')
plt.legend()
plt.axis([0.,10.0,0.,1.2])
plt.show()
../../_images/jupyters_tools_capacity_spectrum_method_charts_21_0.png