1、 一、得到地震波數據1二、地震波分析22.1時程曲線繪制22.2傅里葉譜的繪制32.3反應譜結果分析8參考文獻17附件說明17 一、得到地震波數據訪問/ 網站，下載相應的附錄說明文件，在這些文件中可以找到相關地震的信息，便于有效的規劃檢索詞，可以看到題目需要的地震波RSN 編號為6,；圖1.1查詢相應的編號根據所查到的RSN編號直接檢索，如圖1.2所示輸入RSN號即可檢索圖1.2檢索條件得到檢索結果下載即可得到所需的地震波，如圖1.3所示圖1.3檢索結果具體地震波數據見附件1.二、地震波分析2.1時程曲線繪制圖2.1.1豎向地震加速度時程曲

2、線根據得到的地震波數據，進行MATLAB程序編制，繪制的到相應的豎向地震、180度地震程曲線，270度地震的加速度時程曲線，分別為圖2.1.1，2.1.2,2.1.3. 圖2.1.2 180度地震加速度時程曲線圖2.1.3 270度地震加速度時程曲線2.2傅里葉譜的繪制 根據離散傅里葉的變換準則可以得到Ak（k=0,1，,N/2）、Bk(k=1,.N/2-1)，進而計算得到相應的參數T2*Ak2+Bk2； 式2.1=acrtan(-BkAk); 式2.2f=kT 式2.3Sk=T×Ck2 式2.4式中N=1024，T=10.24s;在具體計算時采用了MATLAB中的FFT函數，并且對

3、參考文獻中1中的例題波進行了試算得到了理想的結果；故可證明程序中所使用的算法是沒有問題的；得到的結果見下。10圖2.2.1 豎向地震波的傅里葉幅值曲線圖2.2.2 180度地震波的傅里葉幅值曲線圖2.2.3 270度地震波的傅里葉幅值曲線圖2.2.4 豎向地震波的傅里葉相位曲線圖2.2.5 180度地震波的傅里葉相位曲線圖2.2.6 270度地震波的傅里葉相位曲線圖2.2.7豎向地震波的功率譜曲線圖2.2.8 180度地震波的功率譜曲線圖2.2.9 270度地震波的功率譜曲線2.3反應譜結果分析 根據地震波是可以求得相應的不同周期的單自由度體系的各種反應譜的，這包括：位移反應譜，速度反應譜，準

4、速度反應譜，加速度反應譜，絕對加速度反應譜，具體的原理是根據杜哈梅積分對單自由度進行積分獲得，相應的公式如下2xt=-1d0tae-(t-)sin(d(t-)dmax 式2.5xt=-d0tae-(t-)cos(dt-+)d max 式2.6xt=-2d0tae-t-sindt-+2d max 式2.7其中，d=(1-2),式中為阻尼比，在本題中分別取為0和0.05。此處對MATLAB計算程序的相應算法進行說明，在實際計算反應譜時使用的是復合的辛普森積分公式3qf=h6fa+4i=1n-1fxi+12+2i=2n-1fxi+f(b) 式2.8因此在計算中是根據t=0.02進行數值積分運算的，這

5、樣選擇的理由有二，1.辛普森公式具有高階的精度，2由于時間間隔較小是可以捕捉到最大值的，不至于引起太大誤差又可以減小計算量，大大較少程序運行時間。同時在計算中結構的角速度的取值參照了振型分解的反應譜法中的取值取了相應的值為=2*pi/linspace(0.2,659);,具體源代碼見附件下面即為相應的反應譜結果：圖2.3.1 豎向地震波位移反應譜圖2.3.2 豎向地震波速度反應譜圖2.3.3 豎向地震波準速度反應譜圖2.3.4 豎向地震波加速度反應譜圖2.3.5 豎向地震波絕對加速度反應譜圖2.36 180度震波位移反應譜圖2.3.8 180地震波速度反應譜圖2.3.9 180地震波準速度反應譜圖2.3.9 180地震波加速度反應譜 圖2.3.10 180地震波絕對加速度反應譜圖2.3.11 270地震波位移反應譜圖2.3.12 270地震波速度反應譜圖2.3.13 270地震波準速度反應譜圖2.3.14 270地震波加速度反應譜圖2.3.15 270地震波絕對反應譜參考文獻1大崎順彥，地震動的譜分析入門，地震出版社，1980；附件說明（1） 附件一為題