基于matlab语言的地震波动力特性分析方法.pdfVIP

№. 1 　　　　　　　　　　　　　　　 　　　　陕西科技大学学报 　　　　　　　　　　　　　　　　　　Feb . 2009 　　　　　　J OU RN AL O F SHAAN XI UN IV ER SIT Y O F SCIEN CE T EC HN OL O GY 　　　　　Vol . 27 ·130 · ( ) 　文章编号 2009 0 10 13005 基于 MATL AB 语言的地震波动力特性分析方法 1 1 1 2 马乐为 , 钟骁珺 , 谢异同 , 张同亿 ( 1 西安建筑科技大学土木工程学院 , 陕西 西安 　7 10055 ;2 中元国际工程设计研究院 , 北京 　100089) 摘 　要 :地震波频谱组成的研究主要是利用 Fourier 分析或小波理论 ,而地震波对建筑物动力 反应的研究则是通过拟加速度反应谱来进行. 现行分析方法的数值化过程必须通过计算机编 程来实现 ,而采用常规的编程语言实现起来往往较为复杂. 基于此 ,作者提出了基于 MA T L AB 平台的分析方法 ,可以大大提高地震波动力特性分析的效率 ,便于分析结果的图形化输 出 ,为工程结构抗震计算提供有益的参考. 关键词 :MA TL AB 语言 ; 地震波频谱组成 ; 拟加速度反应谱 ; Fourier 分析 ; 小波理论 中图分类号 :P3 15 . 3 + 1 ; TP3 12 　　　　　　文献标识码 :A 0 　引言 地震波的动力特性是指其频谱组成以及对不同建筑物的动力反应特点 ,地震波的这些特性对建筑结 构抗震计算是非常重要的. 一般而言 ,在地震波三要素振幅、持时和频率中以第三者最为复杂. 目前 ,有关 地震波频谱特性分析或动力反应研究的理论主要包括时域或频域内的分析方法 , 即时程分析法 、Fourier ( ) 变换法 、拉普拉斯变换法 L ap lace t ran sfer 以及小波分析法等 ;而地震波对建筑物动力反应的研究则往 往是通过拟加速度反应谱来进行的. 结构地震作用计算中的反应谱理论是各国规范普遍采用的计算原则 ,其中以地震波的加速度反应谱 最为有效. 绘制加速度反应谱曲线的常规方法是时域内的 Duhamel 积分 ,其数值积分可利用线性加速度 θ 法或 Will sion 法等. 上述分析方法的数值化过程必须通过计算机编程来实现 ,而采用常规的编程语言如 Fort ran 或 C 语言实现起来往往较为复杂. 其实 ,动力方程的求解也可以通过拉普拉斯变换来实现 ,而且 利用 MA TL AB 提供的函数可使求解过程大大简化. 作为 目前科技人员普遍使用的计算工具 ,MA TL AB 最大的优点在于其极为丰富的函数库以及几乎 ( ) [ 1 ] 能够满足各学科要求的工具箱 Toolbo x ,人们可以用最少的语句完成最复杂的计算 . 本文上述的 Fou rier 变换 、拉普拉斯变换以及小波分析等方法在 MA TL AB 中均可以快速实现 ,从而大大提高地震波动力 特性的分析效率 ,便于分析结果的图形化输出 ,为工程结构抗震计算提供有益的参考. 1 　基于 Matlab 的地震拟加速度反应谱 单自由度线性系统的运动微分方程可表示为 : ( ) ( ) mx¨+ cx+ kx = f t 1 ( )