单质点地震作用计算的计算方法

1、单质点地震作用计算的计算方法所谓单质点弹性体质，是指可以将结构参与振动的全部质量集中于一点，用无重量的弹性直 杆支承于地面上的结构.例如水塔、单层房屋等建筑物，由于它们的质量大部分集中于结构 的顶部，所以通常将这些结构简化成单质点体系.目前，计算弹性体系的反应时，一般假定 地基不产生转动，而把地基的运动分解为一个竖向和两个水平向的分量，然后分别计算这些 运动分量对结构的影响.主要内容：1单自由度弹性体系地震反应分析，主要是运动方程解的一般形式及水 平地震作用的基本公式及计算方法。2.计算水平地震作用关键在于求出地震系数k和动力系数卩。一、地震概述 地震是一种地质现象，就是人们常说的地动，它主要

2、是由于地球的内力作用而产 生的一种地壳振动现象。据统计，地球上每年约有15 万次以上或大或小的地震。人们 能感觉到的地震平均每年达三千次，具有很大破坏性的达100 次。每次中等程度的地震 就会造成重大损失和人员伤亡，研究地震的危害和抗震的方法极有必要，目前已经研究 到了多质点体系地震作用和整体结构的地震作用，但这些研究都离不开单质点地震作用 的计算，我们组准备理论研究并在现有的计算基础上做一点拓展。地震危害直接2005年2月15日新疆乌什发生6.2级地震，经济损失达15757.43万元，主要是土木结 构的房屋破坏严重。近期，云南普洱发生严重的地震，震中位于人口稠密的县城，造成严重 的财产损失和

3、人员伤亡。目前，因灾受伤群众为3 0 0余人，其中3人死亡。全县各乡（镇） 房屋受损严重，土木结构房屋墙体倒塌较多，砖混结构房屋普遍出现墙体开裂，承重柱移位。 作为将来的结构工程师，抗震是我们拦路虎，必须加以重视，那我们先从基础理论着手。三、单质点弹性体系的地震反应目前，我国和其他许多国家的抗震设计规范都采用反应谱理论来确定地震作用。这种计算理论是根据地震时地面运动的实测纪录，通过计算分析所绘制的加速度（在计算中通常采 用加速度相对值）反应谱曲线为依据的。所谓加速度反应谱曲线，就是单质点弹性体系在一 定地震作用下，最大反应加速度与体系自振周期的函数曲线。如果已知体系的自振周期，那 么利用加速度

4、反应谱曲线或相应公式就可以很方便地确定体系的反应加速度，进而求出地震 作用。应用反应谱理论不仅可以解决单质点体系的地震反应计算问题，而且，在一定假设条件 下，通过振型组合的方法还可以计算多质点体系的地震反应。1运动方程的建立 为了研究单质点弹性体系的地震反应，我们首先建立体系在地震作用下的运动方程。图 2-1表示单质点弹性体系的计算简图。HL图2-1单质点弹性体系计算简图由结构动力学方法可得到单质点弹性体系运动方程: mx(t)+cx(t)+kx(t) = mx (t)(2-3)g其中x (t)表示地面水平位移，是时间t的函数，它的变化规律可自地震时地面运动实 测记录求彳得x (t)表示质点对

5、于地面的相对弹性位移或相对位移反应，它也是时间t的函 数，是待求的未知量。若将式(2-3)与动力学中单质点弹性体系在动荷载F(t)作用下的运动方程 mx(t)+cx(t)+kx(t)=F(t)(2-4)进行比较，不难发现两个运动方程基本相同，其区别仅在于式(2-3)等号右边为地震时地面 运动加速度与质量的乘积；而式(2-4)等号右边为作用在质点上的动荷载。由此可见，地面 运动对质点的影响相当于在质点上加一个动荷载，其值等于加兀(t)，指向与地面运动加速 度方向相反。因此，计算结构的地震反应时，必须知道地面运动加速度Xg (t)的变化规律， g而X (t)可由地震时地面加速度记录得到。而X (t

6、)可由地震时地面加速度记录得到。为了使方程进一步简化，设k2 =m(2-5)ccz2、km 2 m(2-6)将上式代入式(2-3)，经简化后得：x(t) +2Z x(t) + 2 x(t)= x (t)g(2-7)式(2-7)就是所要建立的单质点弹性体系在地震作用下的运动微分方程。2运动方程的解答式(2-7)是一个二阶常系数线性非齐次微分方程，它的解包含两个部分：一个是对应于为求方程(2-7)的全部解答，先讨论齐次方程x(t) +2Z x(t) + 2 x(t) = 0(2-8)的通解。由微分方程理论可知，其通解为：x(t) = e z A cos 眉 + B sin 眉)(2-9)式中心=叭

7、:口 ； A和B为常数，其值可由问题的初始条件确定。当阻尼力为0时，式 (2-9)变为：(2-10)x(t) = A cos t + B sin t(2-10)式(2-io)为无阻尼单质点体系自由振动的通解，表示质点做简谐振动，这里=ikrm为无 阻尼自振频率。对比式(2-9)和式(2-10)可知，有阻尼单质点体系的自由振动为按指数函数 衰减的简谐振动，其振动频率为心=J-Z2，称为有阻尼的自振频率。根据初始条件t=0可以确定常数A和B，将t=0和x(t) = x(0)代入式(2-9)得：A = x (0)为确定常数B，对时间t求一阶导数，并将t=0，X(t)= X(0)代入，得：x (0)

8、+Zx (0)B =将 A、B 值代入式(2-9)得：x(t)=ex (0) cos x(t)=ex (0) cos rot +x() +3班0儿说 rot(2-11)上式就是式(2-8)在给定的初始条件时的解答。由ro=和Z = c /2mro可以看出，有阻尼自振频率ro随阻尼系数c增大而减小， 即阻尼愈大，自振频率愈慢。当阻尼系数达到某一数值c时，即rc = c = 2mro = 2km(2-12)r时，则ro = 0，表示结构不再产生振动。这时的阻尼系数c称为临界阻尼系数。它是由结构 r的质量m和刚度k决定的，不同的结构有不同的阻尼系数。而ccZ =(2-13)2mro cr上式表示结构

9、的阻尼系数c与临界阻尼系数c的比值，所以Z称为临界阻尼比,简称阻尼比。r在建筑抗震设计中，常采用阻尼比Z表示结构的阻尼参数。由于阻尼比Z的值很小，它 的变化范围在0.010.1之间，因此，有阻尼自振频率ro = rol-z2和无阻尼自振频率ro很 接近，因此计算体系的自振频率时，通常可不考虑阻尼的影响。(2) 地震作用下运动方程的特解 进一步考察运动方程(2-7) x(t) +2Zro x(t) +ro2x(t) = x (t)g可以看到，方程与单位质量的弹性体系在单位质量扰力作用下的运动方程基本相同，区别仅 在于方程等号右端为地震地面加速度 xg 在于方程等号右端为地震地面加速度 xg(t)

10、， g是随时间而变化的单位质量的“扰力” 为了便于求方程(2-7)的特解，我们将 “扰力”x g (t)看作是无穷多个连续作用的微分脉冲，如图2-2所示。现在讨论任一 微分脉冲的作用。设它在t = T - dT开始作 用，作用时间为dT，此时微分脉冲的大小 为xg (T)dT。显然，体系在微分脉冲作用所以，在求方程的解答时，可将xg (t)看作g后仅产生自由振动。这时，体系的位移可按式(2-3)确定。但式中的x(0)和X(0)应为微分 脉冲作用后瞬时的位移和速度值。根据动量定理：根据动量定理：x(0) = x (T)dT(2-14)g将x(0) =0和X(0)的值代入式(2-3)，即可求得时间

11、T作用的微分脉冲所产生的位移反应X (T)dx = e zm-t) tsin (t - i)di(2-15)将所有组成扰力的微分脉冲作用效果叠加，就可得到全部加载过程所引起的总反应。因此，将式(2-15)积分，可得时间为t的位移1 x(t ) = 一 f X (T)e z(t-t) sin ro(t - i)di(2-16)O o g上式就是非齐次线性微分方程(2-7)的特解，通称杜哈梅(Duhamel)积分。它与齐次微分方 程(2-8)的通解之和就是微分方程(2-7)的全解。但是，由于结构阻尼的作用，自由振动很快 就会衰减，公式(2-9)的影响通常可以忽略不计。分析运动方程及其解答可以看到：

12、地面运动加速度Xg (t)直接影响体系地震反应的大g 小；而不同频率(或周期)的单自由度体系，在相同的地面运动下会有不同的地震反应；阻 尼比Z对体系的地震反应有直接的影响，阻尼比愈大则弹性反应愈小。四、单质点弹性体系水平地震作用1水平地震作用基本公式由结构力学可知，作用在质点上的惯性力等于质量m乘以它的绝对加速度，方向与加 速度的方向相反，即 F(t)= mXg(t)+X(t)(2-17)g式中F(t)为作用在质点上的惯性力。其余符号意义同前。 如果将式(2-3) 代入式(2-17) ，并考 虑到cX(t)远小于kx (t)而略去不计，则得：F (t) = kX(t) = m2X(t)(2-1

13、8)由上式可以看到，相对位移X(t)与惯性力F(t)成正比，因此，可以认为在某瞬时地震作用 使结构产生相对位移是该瞬时的惯性力引起的。也就是为什么可以将惯性力理解为一种能反 应地震影响的等效载荷的原因。将式(2-16)代入式(2-18)，并注意到和的微小差别，令=，则得：F(t) = mro ftX (T)e z(t-t) sin ro(t - T)dT(219)0g由上式可见，水平地震作用是时间t的函数，它的大小和方向随时间t而变化。在结构 抗震设计中，并不需要求出每一时刻的地震作用数值，而只需求出水平作用的最大绝对值 设F表示水平地震作用的最大绝对值，由式(2-19)得：F = mro f

14、x (T)e zw-t)sin(t - F = mro fx (T)e zw-t)sin(t - T)dT0g(2-20)max这里F = mSa(2-21)S = ro fx (T)e zro (t-t)sin ro(t - T)dT0g(2-22)maxgmax代入式(2-21)并以FEk代替F ,Xgmax则得：FEk=mkpg = kpG(2-23)地震动峰值加速度; max标准值)。式中 F 水平地震作用标准值； S 地震动峰值加速度; max标准值)。k 地震系数；卩一动力系数；G 建筑的重力荷载代表值式(2-23)就是计算水平地震作用的基本公式。由此可见求作用在质点上的水平地震作

15、 用F ，关键在于求出地震系数k和动力系数卩。Ek2地震系数 k地震系数k是地震动峰值加速度与重力加速度之比，即(2-24)k =-(2-24)g也就是以重力加速度为单位的地震动峰值加速度。显然，地面加速度愈大，地震的影响 就愈强烈，即地震烈度愈大。所以，地震系数与地震烈度有关，都是地震强烈程度的参 数。3动力系数卩动力系数卩是单质点弹性体系在地震作用下反应加速度与地面最大加速度之比，即S卩=a(2-25)|Xg|max 也就是质点最大反应加速度对地面最大加速度放大的倍数。4.地震影响系数为了简化计算，将上述地震系数k和动力系数卩的乘积用a来表示，并称为地震影响系 数。a = kp(2-26)

16、这样，式(2-23)可以写成F =aG(2-27)Ek图2-3图2-3地震影响系数曲线因为a = k Ba = k B=XgSSaaxgmax(2-28)所以，地震影响系数a就是单质点弹性体系在地震时最大反应加速度（以重力加速度g为单 位）。另一方面，若将式（2-27）写成a = F /G，则可以看出，地震影响系数乃是作用在质 Ek 点上的地震作用与结构重力荷载代表值之比。抗震规范就是以地震影响系数a作为抗震设计依据的，其数值应根据烈度、场地类 别、设计地震分组以及结构自振周期和阻尼比确定。这时水平地震影响系数曲线按图 2-3 确定，形状参数和阻尼调整系数应按教材规定调 整。五，相关抗震的措施抗震就是达到“小震没事，中震可修，大震不倒”。从上面的分析过程可以看出，首先 就是从消耗地震带来的巨大能量入手，那么要增加体系的固有周期和它的柔度，现在通俗讲 的就是柔度方案。关键是用哪一种材料，然后再从结构本身考虑。（一）：目前隔震支座广泛 应用于建筑，桥梁领域。它的原理是在建筑层间设置隔震装置，大大降低隔震层的水平刚度， 延长隔震体系的自震周期，减小结构的地震反应。常用的叠层橡胶支座具有较好的减震效果， 但同时也会带来橡胶的老化等问题。（二）；现在又研究一种盆式橡胶支座，一方面它的钢盆 具有较大的竖向承载能力，另一