地震工程學(xué)-傅里葉變換-反應(yīng)譜計(jì)算

一、得到地震波數(shù)據(jù) /17一、得到地震波數(shù)據(jù)訪問/網(wǎng)站，下載相應(yīng)的附錄說明文件，在這些文件中可以找到相關(guān)地震的信息，便于有效的規(guī)劃檢索詞，可以看到題目需要的地震波RSN編號(hào)為6,；圖圖1.1查詢相應(yīng)的編號(hào)根據(jù)所查到的RSN編號(hào)直接檢索，如圖1.2所示輸入RSN號(hào)即可檢索圖圖1.2檢索條件得到檢索結(jié)果下載即可得到所需的地震波，如圖1.3所示圖圖2.2.3270度地震波的傅里葉幅值曲線圖圖2.2.4豎向地震波的傅里葉相位曲線圖圖2.2.5180度地震波的傅里葉相位曲線圖圖2.2.6270度地震波的傅里葉相位曲線圖圖2.2.7豎向地震波的功率譜曲線圖圖2.2.8180度地震波的功率譜曲線圖圖2.2.9270度地震波的功率譜曲線2.3反應(yīng)譜結(jié)果分析根據(jù)地震波是可以求得相應(yīng)的不同周期的單自由度體系的各種反應(yīng)譜的，這包括：位移反應(yīng)譜，速度反應(yīng)譜，準(zhǔn)速度反應(yīng)譜，加速度反應(yīng)譜，絕對加速度反應(yīng)譜，具體的原理是根據(jù)杜哈梅積分對單自由度進(jìn)行積分獲得，相應(yīng)的公式如下[2]xt=-xt=-x其中，ωd=ω(1-ε2),式中ε為阻尼比，在本題中分別取為0qf=h6f因此在計(jì)算中是根據(jù)△t=0.02進(jìn)行數(shù)值積分運(yùn)算的，這樣選擇的理由有二，1.辛普森公式具有高階的精度，2由于時(shí)間間隔較小是可以捕捉到最大值的，不至于引起太大誤差又可以減小計(jì)算量，大大較少程序運(yùn)行時(shí)間。同時(shí)在計(jì)算中結(jié)構(gòu)的角速度的取值參照了振型分解的反應(yīng)譜法中的取值取了相應(yīng)的值為ω=2*pi/linspace(0.2,659);,具體源代碼見附件下面即為相應(yīng)的反應(yīng)譜結(jié)果： 圖圖2.3.1豎向地震波位移反應(yīng)譜圖圖2.3.2豎向地震波速度反應(yīng)譜圖圖2.3.3豎向地震波準(zhǔn)速度反應(yīng)譜圖圖2.3.4豎向地震波加速度反應(yīng)譜圖圖2.3.5豎向地震波絕對加速度反應(yīng)譜圖圖2.36180度震波位移反應(yīng)譜圖圖2.3.8180地震波速度反應(yīng)譜圖圖2.3.9180地震波準(zhǔn)速度反應(yīng)譜 圖圖2.3.9180地震波加速度反應(yīng)譜圖圖2.3.10180地震波絕對加速度反應(yīng)譜 圖圖2.3.11270地震波位移反應(yīng)譜圖圖2.3.12270地震波速度反應(yīng)譜圖圖2.3.13270地震波準(zhǔn)速度反應(yīng)譜 圖圖2.3.14270地震波加速度反應(yīng)譜 圖圖2.3.15270地震波絕對反應(yīng)譜 [參考文獻(xiàn)][1]大崎順彥，《地震動(dòng)的譜分析入門》，地震出版社，1980；附件說明附件一為題目中所用的地震波數(shù)據(jù)；附件二為相應(yīng)的