def test_get_response_spectrum_2(self):
"""
"""
rss = read_rs(os.path.join(DATA_DIR, '15279_V.txt'))
tgt_02_rs, tgt_05_rs, _, _, _, tgt_30_rs = rss
periods = tgt_02_rs.periods
acc = read_acc(os.path.join(DATA_DIR, '15279_V.cor.acc'))
cal_02_rs = acc.get_response_spectrum(
periods, damping=0.02, gmt='acc', pgm_frequency=150)
cal_05_rs = acc.get_response_spectrum(
periods, damping=0.05, gmt='acc', pgm_frequency=150)
cal_30_rs = acc.get_response_spectrum(
periods, damping=0.30, gmt='acc', pgm_frequency=150)
labels, colors, widths = ['tgt', 'cal'], ['r', 'b'], [3, 1]
fs1 = os.path.join(OUTPUT_DIR,
'recordings.accelerogram.get_response_spectrum.2.02.png')
fs2 = os.path.join(OUTPUT_DIR,
'recordings.accelerogram.get_response_spectrum.2.05.png')
fs3 = os.path.join(OUTPUT_DIR,
'recordings.accelerogram.get_response_spectrum.2.30.png')
plot_response_spectra([tgt_02_rs, cal_02_rs], labels=labels,
colors=colors, widths=widths, unit='m/s2', filespec=fs1)
plot_response_spectra([tgt_05_rs, cal_05_rs], labels=labels,
colors=colors, widths=widths, unit='m/s2', filespec=fs2)
plot_response_spectra([tgt_30_rs, cal_30_rs], labels=labels,
colors=colors, widths=widths, unit='m/s2', filespec=fs3)
def test_get_response_spectrum_2(self):
"""
"""
tgt_rss = read_rs(os.path.join(DATA_DIR, '15279_V.txt'))
tgt_02_rs, tgt_05_rs, tgt_07_rs, tgt_10_rs, tgt_20_rs, tgt_30_rs = tgt_rss
periods = tgt_02_rs.periods
acc = read_acc(os.path.join(DATA_DIR, '15279_V.cor.acc'))
cal_02_rs = acc.get_response_spectrum(periods, damping=0.02, gmt='acc')
cal_05_rs = acc.get_response_spectrum(periods, damping=0.05, gmt='acc')
cal_30_rs = acc.get_response_spectrum(periods, damping=0.30, gmt='acc')
labels = ['tgt', 'cal']
colors = ['r', 'b']
widths = [2, 0.5]
fs1 = os.path.join(FIG_DIR, 'recordings.accelerogram.get_response_spectrum.2.02.png')
fs2 = os.path.join(FIG_DIR, 'recordings.accelerogram.get_response_spectrum.2.05.png')
fs3 = os.path.join(FIG_DIR, 'recordings.accelerogram.get_response_spectrum.2.30.png')
plot.response_spectra([tgt_02_rs, cal_02_rs], labels=labels, colors=colors, widths=widths, unit='m/s2', png_filespec=fs1)
plot.response_spectra([tgt_05_rs, cal_05_rs], labels=labels, colors=colors, widths=widths, unit='m/s2', png_filespec=fs2)
plot.response_spectra([tgt_30_rs, cal_30_rs], labels=labels, colors=colors, widths=widths, unit='m/s2', png_filespec=fs3)
def test_read_rs(self):
"""
"""
rs = read_rs(os.path.join(DATA_DIR, '15279_V.txt'))[1]
periods = """
0.000 0.010 0.020 0.030 0.040 0.050 0.075 0.100 0.110 0.120 0.130
0.140 0.150 0.160 0.170 0.180 0.190 0.200 0.220 0.240 0.260 0.280
0.300 0.320 0.340 0.360 0.380 0.400 0.420 0.440 0.460 0.480 0.500
0.550 0.600 0.650 0.700 0.750 0.800 0.850 0.900 0.950 1.000 1.100
1.200 1.300 1.400 1.500 1.600 1.700 1.800 1.900 2.000 2.200 2.400
2.600 2.800 3.000 3.200 3.400 3.600 3.800 4.000 4.200 4.400 4.600
4.800 5.000 5.500 6.000 6.500 7.000 7.500 8.000 8.500 9.000 9.500
10.00
"""
periods = list(np.array(periods.split(), dtype=float))
f_responses = [
3.9679000E-001, 4.0477622E-001, 4.0532070E-001, 5.5488199E-001,
7.5385058E-001, 7.2152662E-001, 7.3033845E-001, 8.3194721E-001,
7.8486401E-001, 9.5955211E-001, 1.0797322E+000, 1.2527586E+000,
1.4042230E+000, 1.2237327E+000, 1.2319688E+000, 1.0519416E+000,
]
l_responses = [
4.2007200E-003, 3.7803600E-003, 3.5943200E-003, 3.4261000E-003,
3.2630000E-003, 3.3139600E-003, 3.0978400E-003, 2.7881100E-003,
2.1490700E-003, 1.9131700E-003, 1.7104300E-003, 1.5357800E-003,
1.3851500E-003, 1.2545400E-003, 1.1407000E-003, 1.0410500E-003,
]
self.assertListEqual(list(rs.periods), periods)
self.assertListEqual(list(rs.responses[:+16]), f_responses)
self.assertListEqual(list(rs.responses[-16:]), l_responses)
self.assertEqual(rs.unit, 'm/s2')
self.assertEqual(rs.damping, 0.05)
def test_read_rs(self):
"""
"""
rs = read_rs(os.path.join(DATA_DIR, '15279_V.txt'))[1]
periods = """
0.000 0.010 0.020 0.030 0.040 0.050 0.075 0.100 0.110 0.120 0.130
0.140 0.150 0.160 0.170 0.180 0.190 0.200 0.220 0.240 0.260 0.280
0.300 0.320 0.340 0.360 0.380 0.400 0.420 0.440 0.460 0.480 0.500
0.550 0.600 0.650 0.700 0.750 0.800 0.850 0.900 0.950 1.000 1.100
1.200 1.300 1.400 1.500 1.600 1.700 1.800 1.900 2.000 2.200 2.400
2.600 2.800 3.000 3.200 3.400 3.600 3.800 4.000 4.200 4.400 4.600
4.800 5.000 5.500 6.000 6.500 7.000 7.500 8.000 8.500 9.000 9.500
10.00
"""
periods = list(np.array(periods.split(), dtype=float))
f_responses = [
3.9679000E-001, 4.0477622E-001, 4.0532070E-001, 5.5488199E-001,
7.5385058E-001, 7.2152662E-001, 7.3033845E-001, 8.3194721E-001,
7.8486401E-001, 9.5955211E-001, 1.0797322E+000, 1.2527586E+000,
1.4042230E+000, 1.2237327E+000, 1.2319688E+000, 1.0519416E+000,
]
l_responses = [
4.2007200E-003, 3.7803600E-003, 3.5943200E-003, 3.4261000E-003,
3.2630000E-003, 3.3139600E-003, 3.0978400E-003, 2.7881100E-003,
2.1490700E-003, 1.9131700E-003, 1.7104300E-003, 1.5357800E-003,
1.3851500E-003, 1.2545400E-003, 1.1407000E-003, 1.0410500E-003,
]
self.assertListEqual(list(rs.periods), periods)
self.assertListEqual(list(rs.responses[:+16]), f_responses)
self.assertListEqual(list(rs.responses[-16:]), l_responses)
self.assertEqual(rs.unit, 'm/s2')
self.assertEqual(rs.damping, 0.05)
