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

一、得到地震波數(shù)據(jù)1二、地震波分析22.1時程曲線繪制22.2傅里葉譜的繪制32.3反應(yīng)譜結(jié)果分析8參考文獻17附件說明17 一、得到地震波數(shù)據(jù)訪問/ 網(wǎng)站，下載相應(yīng)的附錄說明文件，在這些文件中可以找到相關(guān)地震的信息，便于有效的規(guī)劃檢索詞，可以看到題目需要的地震波RSN 編號為6,；圖1.1查詢相應(yīng)的編號根據(jù)所查到的RSN編號直接檢索，如圖1.2所示輸入RSN號即可檢索圖1.2檢索條件得到檢索結(jié)果下載即可得到所需的地震波，如圖1.3所示圖1.3檢索結(jié)果具體地震波數(shù)據(jù)見附件1.二、地震波分析2.1時程曲線繪制圖2.1.1豎向地震加速度時程曲線根據(jù)得到的地震波數(shù)據(jù)，進行MATLAB程序編制，繪制的到相應(yīng)的豎向地震、180度地震程曲線，270度地震的加速度時程曲線，分別為圖2.1.1，2.1.2,2.1.3. 圖2.1.2 180度地震加速度時程曲線圖2.1.3 270度地震加速度時程曲線2.2傅里葉譜的繪制 根據(jù)離散傅里葉的變換準則可以得到Ak（k=0,1，,N/2）、Bk(k=1,.N/2-1)，進而計算得到相應(yīng)的參數(shù)T2*Ak2+Bk2； 式2.1=acrtan(-BkAk); 式2.2f=kT 式2.3Sk=TCk2 式2.4式中N=1024，T=10.24s;在具體計算時采用了MATLAB中的FFT函數(shù)，并且對參考文獻中1中的例題波進行了試算得到了理想的結(jié)果；故可證明程序中所使用的算法是沒有問題的；得到的結(jié)果見下。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反應(yīng)譜結(jié)果分析 根據(jù)地震波是可以求得相應(yīng)的不同周期的單自由度體系的各種反應(yīng)譜的，這包括：位移反應(yīng)譜，速度反應(yīng)譜，準速度反應(yīng)譜，加速度反應(yīng)譜，絕對加速度反應(yīng)譜，具體的原理是根據(jù)杜哈梅積分對單自由度進行積分獲得，相應(yīng)的公式如下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計算程序的相應(yīng)算法進行說明，在實際計算反應(yīng)譜時使用的是復(fù)合的辛普森積分公式3qf=h6fa+4i=1n-1fxi+12+2i=2n-1fxi+f(b) 式2.8因此在計算中是根據(jù)t=0.02進行數(shù)值積分運算的，這樣選擇的理由有二，1.辛普森公式具有高階的精度，2由于時間間隔較小是可以捕捉到最大值的，不至于引起太大誤差又可以減小計算量，大大較少程序運行時間。同時在計算中結(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 豎向地震波準速度反應(yīng)譜圖2.3.4 豎向地震波加速度反應(yīng)譜圖2.3.5 豎向地震波絕對加速度反應(yīng)譜圖2.36 180度震波位移反應(yīng)譜圖2.3.8 180地震波速度反應(yīng)譜圖2.3.9 180地震波準速度反應(yīng)譜圖2.3.9 180地震波加速度反應(yīng)譜 圖2.3.10 180地震波絕對加速度反應(yīng)譜圖2.3.11 270地震波位移反應(yīng)譜圖2.3.12 270地震波速度反應(yīng)譜圖2.3.13 270地震波準速度反應(yīng)譜圖2.3.14 270地震波加速度反應(yīng)譜圖2.3.15 270地震波絕對反應(yīng)譜參考文獻1大崎順彥，地震動的譜分析入門，地震出版社，1980；附件說明（1） 附件一為題