ruau[k] = ruau_rec
Ruau_rec[i] = ruau
# Original Acceleration file
try:
earthquake = cH.calleqdata(test_id[test_name.index(i)][0])
acceleration = earthquake['Data']*earthquake['Scale_Factor']*Magnitude
DT_a = earthquake['DT']
except:
acceleration = np.array([0, 0, 0, 0])
DT_a = 0.05
eq_time_a = np.arange(offset_eq, len(acceleration)*DT_a+offset_eq, DT_a)
orig_earthquake[i] = [acceleration, eq_time_a, DT_a]
# design spectra
[T, dSA, dSD, eta]=design_spectra()
#####################################################################################
# Figure stuff
fontsize = 8
rcParams['axes.labelsize'] = fontsize
rcParams['xtick.labelsize'] = fontsize
rcParams['ytick.labelsize'] = fontsize
rcParams['legend.fontsize'] = fontsize
rcParams['axes.titlesize'] = fontsize
colors = cm.rainbow(np.linspace(0.05, 0.95, len(test_name)))
colors2 = cm.rainbow(np.linspace(0.05, 0.95, 2))
colors3 = cm.rainbow(np.linspace(0.05, 0.95, len(rec)+1))
colors4 = cm.rainbow(np.linspace(0.05, 0.95, 4))
midiHi = mi*di*Hi
deff = sum(midi2)/sum(midi)
meff = sum(midi)/deff
Heff = sum(midiHi)/sum(midi)
sig_yield = 0.0004*bay_length/0.25*100 #0.0004*eps_y*L_b/h_b *100
diy = Hi*sig_yield/100
midiy = mi*diy
midiy2 = mi*diy**2
dyeff = sum(midiy2)/sum(midiy)
mu = deff/dyeff
# design spectra
[T, dSA, dSD, eta] = design_spectra(mu=mu, Z=0.4, Lambda=Lambda, Sp=1, alpha=0.5)
if PT == 1:
eta = 1
# design values
Teff = np.interp(deff, dSD*eta, T)
Teff5 = np.interp(deff, dSD, T)
Keff = 4*np.pi**2*meff/Teff**2
Vb = Keff*deff
FACTOR = 0.45 # LATEX text width
RATIO = (np.sqrt(5) - 1.0) / 2.0 # Height : width
fig_dims = [FACTOR*6.2677, RATIO*FACTOR*6.2677] # fig dims as a list
fontsize = 8
rcParams['axes.labelsize'] = fontsize
plt.title('ADRS')
plt.plot(avged, avge, 'r', lw=2)
plt.plot(target_disp, target, 'k', lw=2)
plt.plot()
plt.show()
# shear force - displacement plot
if plotter == 7:
for record in Anadata.listNodes(HDFgroupname, classname='Table'):
if record._v_name == motion:
shear_force = -record[30]['data']
base_acceleration = record[7]['data']
drift = (record[2]['data'])/3.2
[T, dSA, dSD, RSA, RSDU] = design_spectra(mu=3.7)
design_force = RSA*5.99*9.81 # 5.99 is effective building mass
scaled_dSD = RSDU*(scale)/(0.834*3.2) # 0.834*H is effective building height
plt.figure(plotter)
plt.plot(drift, shear_force, lw=2)
plt.plot(scaled_dSD, design_force, lw=2)
plt.title('ADRS plot of model building')
plt.xlabel('Building drift')
plt.ylabel('Base shear [kN]')
plt.ylim(0, 70)
plt.xlim(0,0.08)
plt.show()
# peak floor-beam displacement
if plotter == 8:
