地震动反映谱的数值计算（含fortran程序）.docVIP

PAGE PAGE 1 地震动反映谱的数值算 冯宏光 （云南大学资源环境与地球科学学院，云南 昆明 650091） 摘要：本文从单质点系入手，通过解方程而得到直接积分法中的线性加速度法的递推公式，通过递推公式以及线性加速度发的主要思想，利用fortran编写程序得到地震动加速度反应谱的计算程序进而计算出地震动加速度度反应谱。 关键词：单质点系、线性加速度法、地震动反应谱 Compute the Date of Earthquake Response Spectrum Feng hongguang （Yunnan University, School of Resources Environment and Earth Scitnces,Kunming Yunnan China 650091） Abstract: Thesis is base on the single quality point department, and square distance in solution got the formula of linear acceleration method, and make use of the main thought of linear acceleration method, using the procedure of fortran got the earthquake response spectrum. Key words：single quality point department、linear acceleration method、earthquake response spectrum 序言 地震是一项难以预报的自然灾害，随着人类文明的发展地震所带来的危害被更多人重视，其一次性所带来的灾害往往是难以用简单的经济数据来描述的，而地震中危害最大的就是房屋的倒塌所带来的损失。进入二十一世纪以来地震安全性评估就提出了小震不坏，中震可修，大震不倒的宗旨。因此伴随着计算机技术的急速发展利用计算机编程计算地震动反应谱成为了普遍方法。 通过研究地震动反应谱我们可以的道地震动与地面建筑响应的关系，通过了解这种关系就可以指导我们修建地面建筑，因此反应谱的研究无疑是为我们地面建筑的修建指明了一条光明大道。 因此我们可以看到研究地震动反应谱的必要性。 第一章：单质点系的震动 我们知道地震波是弹性波，弹性波最简单的模型便是单质点系的振动，同时反应谱本身是以单质点系振动理论为基础建立的，单质点系的振动对于研究反应谱有很重要的意义。 我们采用Kelvin Voit模型来研究单质点振动，此模型为一弹簧与一个活塞并联。我们设质点的质量为m，弹簧的弹性系数为k，阻尼器的阻尼系数为c，并且设该模型只可以在一个方向上是自由的。 如上图所示x+y为质点的绝对位移，x为质点的相对位移。那么根据牛顿第三定律我们可以得到： （1.1） 1.1考虑了模型中同时存在阻力与粘滞力。当仅仅是考虑简单的一维无阻尼振动时我们可以得到：，我们可以解得： （1.2） 根据初始条件： 将A、B带入1.2式我们可以得到： （1.3） 我们便可以得到一维无阻尼振动的特解 当我们引入阻尼后我们便可以得到： （1.4） 通过解1.4式我们可以解得： （1.5） 1.5式是对于有阻尼的情况的通解，令；：我们便可以得到： （1.6） 引入初始条件： 我们可以得到： （1.7） 第二章：地震动反应 第一章的结论只考虑了单质点震动的问题，但是在实际情况中以单摆式地震仪为例，第一章的结果并没有考虑仪器的响应问题，也就是说仪器对于单质点产生的地震动会做何反应，大而言之我们可以认为地面建筑对于地震动的反应是如何的。 在解决了单质点有阻尼的振动后我们便要考虑到波动响应的问题。对第一章中的（1.1）便是考虑了场地反应的问题。已知：，其中是由地震动引起的惯性力。该式的物理意义可以理解为，假定地面在空间中固定不动，而惯性力代替了地面的运动。惯性力与系统的弹性恢复力、阻尼力构成平衡力系，从而把运动状态变成了静止状态。 在第一章中已知， 那么1.1式可以成为： （2.1） 惯性力是作用于质点上的随时间而变化的力，因此可令： 于是（1.1）式为：