-
Notifications
You must be signed in to change notification settings - Fork 0
/
Stress_Strain_2.py
46 lines (34 loc) · 1.12 KB
/
Stress_Strain_2.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
import numpy as np
import matplotlib.pylab as plt
from matplotlib.backends.backend_pdf import PdfPages
import subprocess
output_filename = 'Stress_Strain_2.pdf'
OpenSees_files = [
'Output/story1_SS_EC_Comp2.txt',
'Output/story1_SS_CC_Comp2.txt',
'Output/story1_SS_ES_Comp2.txt',
'Output/story1_SS_EC_Ten2.txt',
'Output/story1_SS_CC_Ten2.txt',
'Output/story1_SS_ES_Ten2.txt',
]
pp = PdfPages(output_filename)
fig = plt.figure()
count = 1
for file_name in OpenSees_files:
figure = fig.add_subplot(3,2,count)
count = count + 1
print file_name
OpenSees = np.loadtxt(file_name)
plt.minorticks_on()
plt.grid(b=True, which='major', alpha=1)
plt.tight_layout()
plt.plot(OpenSees[:, 2], OpenSees[:, 1], lw=1, label='Fiber 1')
# plt.xlim(-5, 5)
# plt.ylim(-800, 800)
# plt.legend(loc=0, shadow=True, numpoints=1)
# plt.xlabel("Normalized axial displacement ($\Delta/h_w$), %")
# plt.ylabel(r"Normalized Axial Core Stress, $P_c / (A_\mathrm{core} \, f'_c)$")
# plt.grid(b=True, which='minor', alpha=0.2)
pp.savefig(fig)
pp.close()
subprocess.Popen([output_filename], shell=True)