地震工程学-傅里叶变换-反应谱计算教案分析.docx

 TOC \o 1-3 \h \z \u  HYPERLINK \l _Toc420097238 一、得到地震波数据  PAGEREF _Toc420097238 \h 1  HYPERLINK \l _Toc420097239 二、地震波分析  PAGEREF _Toc420097239 \h 2  HYPERLINK \l _Toc420097240 2.1时程曲线绘制  PAGEREF _Toc420097240 \h 2  HYPERLINK \l _Toc420097241 2.2傅里叶谱的绘制  PAGEREF _Toc420097241 \h 3  HYPERLINK \l _Toc420097242 2.3反应谱结果分析  PAGEREF _Toc420097242 \h 8  HYPERLINK \l _Toc420097243 [参考文献]  PAGEREF _Toc420097243 \h 17  HYPERLINK \l _Toc420097244 附件说明  PAGEREF _Toc420097244 \h 17    PAGE \* Arabic \* MERGEFORMAT17 / 17 一、得到地震波数据 访问 HYPERLINK / / 网站，下载相应的附录说明文件，在这些文件中可以找到相关地震的信息，便于有效的规划检索词，可以看到题目需要的地震波RSN 编号为6,； 图1.1查询相应的编号 根据所查到的RSN编号直接检索，如图1.2所示输入RSN号即可检索 图1.2检索条件 得到检索结果下载即可得到所需的地震波，如图1.3所示 图1.3检索结果 具体地震波数据见附件1. 二、地震波分析 2.1时程曲线绘制 图2.1.1竖向地震加速度时程曲线 根据得到的地震波数据，进行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.2 f=kT 式2.3 Sk=T×Ck2 式2.4 式中N=1024，T=10.24s; 在具体计算时采用了MATLAB中的FFT函数，并且对参考文献中[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反应谱结果分析 根据地震波是可以求得相应的不同周期的单自由度体系的各种反应谱的，这包括：位移反应谱，速度反应谱，准速度反应谱，加速度反应谱，绝对加速度反应谱，具体的原理是根据杜哈梅积分对单自由度进行积分获得，相应的公式如下[2] xt=-1ωd0taτe-εω(t-τ)sin?(ωd(t-τ))dτmax 式2.5 xt=-ωωd0taτe-εω(t-τ)cos?(ωdt-τ+α)dτ max 式2.6 xt=-ω2ωd0taτe-εωt-τsinωdt-τ+2αdτ max 式2.7 其中，ωd=ω(1-ε2),式中ε为阻尼比，在本题中分别取为0和0.05。此处对MATLAB计算程序的相应算法进行说明，在实际计算反应谱时使用的是复合的辛普森积分公式[3] qf=h6fa+4i=1n-1fxi+12+2i=2n-1fxi+f(b)