地震工程学傅里叶变换反应谱计算.docx
- 文档编号:30399482
- 上传时间:2023-08-14
- 格式:DOCX
- 页数:18
- 大小:922.79KB
地震工程学傅里叶变换反应谱计算.docx
《地震工程学傅里叶变换反应谱计算.docx》由会员分享,可在线阅读,更多相关《地震工程学傅里叶变换反应谱计算.docx(18页珍藏版)》请在冰豆网上搜索。
地震工程学傅里叶变换反应谱计算
一、得到地震波数据1
二、地震波分析2
2.1时程曲线绘制2
2.2傅里叶谱的绘制3
2.3反应谱结果分析8
[参考文献]17
附件说明17
一、得到地震波数据
访问http:
//ngawest2.berkeley.edu/网站,下载相应的附录说明文件,在这些文件中可以找到相关地震的信息,便于有效的规划检索词,可以看到题目需要的地震波RSN编号为6,;
根据所查到的RSN编号直接检索,如图1.2所示输入RSN号即可检索
得到检索结果下载即可得到所需的地震波,如图1.3所示
具体地震波数据见附件1.
二、地震波分析
2.1时程曲线绘制
根据得到的地震波数据,进行MATLAB程序编制,绘制的到相应的竖向地震、180度地震程曲线,270度地震的加速度时程曲线,分别为图2.1.1,2.1.2,2.1.3.
2.2傅里叶谱的绘制
根据离散傅里叶的变换准则可以得到Ak(k=0,1,…,N/2)、Bk(k=1,….N/2-1),进而计算得到相应的参数
;式2.1
;式2.2
式2.3
式2.4
式中N=1024,T=10.24s;
在具体计算时采用了MATLAB中的FFT函数,并且对参考文献中[1]中的例题波进行了试算得到了理想的结果;故可证明程序中所使用的算法是没有问题的;得到的结果见下。
10
2.3反应谱结果分析
根据地震波是可以求得相应的不同周期的单自由度体系的各种反应谱的,这包括:
位移反应谱,速度反应谱,准速度反应谱,加速度反应谱,绝对加速度反应谱,具体的原理是根据杜哈梅积分对单自由度进行积分获得,相应的公式如下[2]
式2.5
式2.6
其中,
式中ε为阻尼比,在本题中分别取为0和0.05。
此处对MATLAB计算程序的相应算法进行说明,在实际计算反应谱时使用的是复合的辛普森积分公式[3]
式2.8
因此在计算中是根据△t=0.02进行数值积分运算的,这样选择的理由有二,1.辛普森公式具有高阶的精度,2由于时间间隔较小是可以捕捉到最大值的,不至于引起太大误差又可以减小计算量,大大较少程序运行时间。
同时在计算中结构的角速度的取值参照了振型分解的反应谱法中的取值取了相应的值为
具体源代码见附件
下面即为相应的反应谱结果:
[参考文献]
[1]大崎顺彦,《地震动的谱分析入门》,地震出版社,1980;
附件说明
(1)附件一为题目中所用的地震波数据;
(2)附件二为相应的MATLAB程序,其中tfp.m为相应的时程分析,傅里叶幅值、相位、功率谱分析绘图程序;resp.m为反应谱分析程序.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 地震 工程学 傅里叶变换 反应 计算