振动的测量傅里叶变换duhamel积分反应谱

1、8振动的测量8.1 前言有的时候，一些微小的、不显著的振动，会与结构，或者结构的某一部分产生共振，从而将振动放大。共振也会发生在人的身上， 人体的自振频率大概为 7.5Hz,因此次声（20Hz） 会对人体造成伤害。所以说，对于结构来说，利用合适的装置或者设计来减小这样的共振是非常有必要的。那么，想要研究如何减小共振，我们首先要知道将要发生的振动的参数。想要知道这些参数，我们就需要一些仪器来测量，这些仪器就是我们这章要了解的。首先来看一下一些概念。 在结构工程中常常进行运动量 （位移、速度或加速度）的测量， 例如地震动时程的测量；振动台试验中结构模型的动力反应的测量； 脉动作用下结构物的振 动的

2、测量；大桥、超高层结构风振的测量等。用于测量振动量的仪器（拾振仪）主要有三种：加速度位移计：测量加速度的时程（强震仪）。位移计：测量位移时程（地震仪） 。速度计：测量速度。8.2 理论8.2.1 运动方程的建立D' Alembert理：在质点系的运动的任意瞬间，如果除了实际作用于每一质点的主动力和约束反力外，再加上假想的惯性力，则在该瞬间质点系将处于假想的平衡状态，称之为动力平衡状态。记？? ?希别为质点？?所受的主动力、惯性力和约束反力，则D' Alembert原理可表示为通常主动力？泡括外荷载、阻尼力和弹性恢复力。上图质量块m所受的主动力为?(?= ?(?) c?- ?惯性

3、力为?= -?由于该体系是约束反力不做功的理想约束体系，故列运动方程时仅考虑运动方向上的受力，此时的约束反力是没有的。将上面两式代入D' Alembert原理表达式，有?+ c?+ ?= ?(?)当然，建立运动方程的方法有多种，除了上面介绍的D' Alembert原理之外，还有虚位移原理、Hamilton原理和Lagrange方程，这四种方法对建立运动方程是完全等同的，可以 推得完全相同的运动方程。8.2.2 Fourier变化法(频域分析法)最简单的测量仪器模型是一单自由度弹簧-质点-阻尼体系，被封闭在一个刚性盒子里面,如图所示单自由度体系运动方程为：?+ c?+ ?= -?

4、(? (1)其中：c= 2?/？加则(1)式可以写为：?+ 2?r?+ ?2?= -?(? (2)使用傅里叶变换法(之后补上介绍)，正变换，把问题从时间域(自变量为t)转变到频域(自变量为？，可得：-?2?+ ?2?(? + ?2? = ? ? (3)?苗？?7?) = ?)z二=?才？2? ? (4)?(1 -法)+ ?2旗?|乃| 二 2?1？田、2（1 - 帝）172_? 2+ (2? (5)_ 一 . 一 ？ ？.卜面给出了 |赤？|与？?关于频率比的图像: （ , ?-0.00 '£+ 0.25, 0.70 ，-1.001_ 0.00 1 .,0 0.50 -11*

5、-彳如。.了0.,Lno-?（?内复频反应函数，也叫传递函数。相角？2的含义，在动力荷载作用下，有阻尼体系的动力反应（位移、速度、加速度）一定要滞后动力荷载一段时间，即存在反应滞后现 象。相角？校际是反映结构体系位移相对于动力荷载的反应滞后时间，从下图可以发现，频 率比越大，即外荷载作用的越快，动力反应的滞后时间越长。补充：傅立叶变换Fourier变换的定义为oo?1? = / ?-正变换-81?=?(? ?-?-逆变换2? -oo其中，？称为位移？?刑Fourier谱。根据Fourier变换的性质，速度和加速度的Fourier变换为oo/ 7?1? ?(?) -8? ? -? 2 ?(?)-

6、ooOO8.3振动测量仪器8.3.1 加速度计(强震仪)加速度计测量的是加速度，在基底加速度作用下仪器质点的运动方程为式(1):?+ c?+ ?= -?(?)包含在一定的频段分布的一系列简谐分量,实际要测量的加速度运动时程是任意变化的， 我们可以先分析仪器对一个简谐分量的测量。仪器基底加速度时程:?)?= ? ? ? ? ? ?为简谐运动，其中？?？为地面运动加速度的振幅。带入上式，可得仪器质点的相对位移？?为?籥？1?=k一 sin (? ?/2 2°V1 - (?) + 2?,?2?=-9? (? ?)? (8)下面仅讨论？?勺振幅？?=9?？?)? (9)作为一种测量仪器，其测

7、量的对象的单位可以不同，但一定要成比例。由上式可知，仪器记录值 ？与被测量的地面加速度值之间的关系变化由动力放大系数?决定，而动力放大系数 ？?是由？?勺函数，观察不同阻尼比时动力放大系数曲线图(下图)J;Q_9上J.654 $21 J Qo_0o.o.o.oro.o.o.4:第味KKS留可以发现，当？?= 0.7时，在频段0W?W0.5范围内，？?=1为常数。< ；0.60 f = 0.65 <=g7Q < -0,75即在以上的频段范围内，仪器反应的振幅？，或者说是对于不同的简谐运动，仪器的记录与仪器要测量的加速的振幅成线性关系。这时可以用仪器反应？来度量要测量的加速的。为

8、了保证加速的不同简谐分量的频率？?TB满足0W?W0.5,可以采用提高？?痢方法来实现，因为？?= v?/?>, 一般情况下，质量？环变，可以采用提高加速度计中的弹簧刚度？的方法来实现提高？?甥目的。因此，加速度计或强震仪中的弹簧刚度比较大，仪器是比较刚性的。8.3.2 位移计(地震仪)位移计是用来测量仪器基底的位移量，设任意一个简谐位移为?= ?0? ? ? ? ? ? (10)其中？粉为地面运动位移振幅，将？?？= -? 2?0?沌?嬷动方程式（1）,解得仪器的相对位移反应为：?= ?夕?0? ?= I / f ： < ： W ' 1? 2(亚?0? ? ? (11)相

9、对位移反应的振幅为:?=（1）2?3北?0? ? ? ? ? （12）?与加速度计的原理相同，我们希望在所测量的振动频率范围内，对于不同的频率分量均 有?？和？?江间的比例系数（？?笳2?接近常量。下图给出了（？?？为2?抽频率（？?？）变化的曲线。由图可以看出，当阻尼比 ？?= 0.5时，在频率范围？?？?> 1时，（？）2?獗近于常数,因此可以用位移计来测量频率范围在?> ?>撷段范围内的位移量。为了保证被测位移的频 率满足?> ?r?可通过降低仪器自振频率 ？?甥方法来实现。实际中采取降低弹簧刚度 ？或者 增大质量？硒方法来实现。因此，位移计一般都是比较柔的。第十

10、章反应谱10.1 前言地震动引起地面的运动，并通过地面的运动，是结构也产生振动。因此，在地震中，结构上所受的荷载是由于其支座的运动而产生的。地面的运动有三个平动分量和三个转动分量，但是由于测量水平的限制，转动分量很难测得，而平动的分量可以由加速度计测得。相对于平动分量来说，转动分量很小，因此在对结构进行抗震分析的时候，转动分量忽略不计。对一个结构来说，在弹性范围内，它的响应是由于地面的一个平动分量产生的，对于一个结构体系来说，就是这些分量的和。由于结构的自振频率是未知的，在设计的时候需要多次迭代才能求出来。所以，结构设计者就需要反应谱的帮助了。10.2 傅立叶谱将振动的信号（或任意变化的函数）

11、分解为简谐振动（三角函数）的过程称为傅里叶分 解。得到振幅和相位随频率变化的关系称为傅里叶谱，包括振幅谱和相位谱，统称为频谱，完成分解的运算称为傅立叶变化。傅里叶谱全面描述了地震动过程的频谱特征，包括了各频率分量的相位及幅值信息。因此，从两个傅里叶谱可以反推出地震动的时程，而功率谱和反应谱则不行。傅立叶谱是复数，由实部和虚部组成，他的模称为幅值谱，幅角为相位谱。oo? = / 7V?)?7? (13)-8假设地面运动的加速度在 ？?(咤，则上式可以变为?/? = / ? ?)? ? ? ? ? (14)00则强震下傅立叶振幅谱和相位谱可以用下式定义?2?2|?| = v / ?+ / ? ?

12、? ? ? (15)00?.匚？?？?？?？d, Y /, , V - - 72? = -?-0?j? (16)，？?？?？?？力 Y /, 110.3反应谱单自由度体系在给定的地震作用下某个最大反应与体系自振周期的关系曲线称为该反 应的地震反应谱。10.3.1 Duhamel 积分(1)单位脉冲反应函数单位脉冲：作用时间很短，冲量等于1的荷载，实际上就是数学中的特殊函数一？函数。？函数的定义为8 , ? ? ? = 0,其他？ (17)oo? ? ? 10? (18)在?= ?时刻一个单位脉冲？?= ?乍用在单自由度体系上，使结构的质点获得一个单位冲量，在脉冲结束后，质点获得一个初速度，即?

13、+?+?(19)? ? = f ? r ? i ? Y r*/ /, , o/ ?当？?0时，1?=刃由于脉冲作用时间很短，当 ？A0时，由单位脉冲引起质点的位移为零，即?= 0求体系在单位脉冲下的反应，即是求解单位脉冲作用后的自由振动问题。将初始条件1?=不和？?= 0代入单自由度体系自由振动(低阻尼体系)一般解式?= ?,?z?s(0)?+>( ?0) + ?；梆?(0)? ? ? ? (20)所以有阻尼体系单位脉冲反应函数为1»?(? ?= ?= ?-?(?-?(? ?) ?俗? ? ? ? (21)?(2)对任意荷载的反应我们可以把荷载分解为一系列脉冲，获得每一个脉冲下

14、结构的反应，最后叠加得到结构的总反应。如果已经将作用于结构体系的外荷载？?离散成一系列脉冲，首先计算其中任一脉冲?动力反应。则该脉冲下结构的反应为? (22)?= ? ?,?> ?在任意时间？结构的反应，就是在？之前所有脉冲作用下的反应之和为? (23)?= r ? r ?(? ? y %/ %/ a z f将(21)式代入上式，可以求解阻尼体系动力反应的Duhamel积分公式1?= ?夕?9?-? ?p? ? ? ? ? ? (24)?L?C? 0其中？?= ?/1- ?为阻尼体系的自振频率。?球冲引您的反应时刻脓同引起的反应总反应10.3.2反应谱法前面介绍了 一些结构动力反应分析的

15、方法，可以对结构地震反应问题展开分析计算。当地震动较小时，结构处于线弹性范围，可以采用时域的 Duhamel积分法，或频域的Fourier变换方法获得地震下结构的反应，并根据得到的结构最大变形最大内力进行抗震设 计。当地震动较强时，结构反应可能进入塑性，需要用到时域逐步积分法进行弹塑性反应 分析。我们仅讨论结构线弹性地震反应问题，采用Duhamel积分法介绍地震反应谱。地震作用的特点是地震动过程非常复杂，随时间不规则、快速变化。设地震加速度时 程为？没?，其特点为：第一阶段，振幅快速增长；第二阶段，相对稳定；第三阶段，震荡 衰减。地震作用下结构的运动方程为?+ c?+ ?= -?(?)地震等效

16、荷载为? = -?)应用Duhamel积分，结构地震的位移反应为-1?Ci”?= 一 / 2然取可?2 f ? ?)? ? ? ? ? ? (25)? 0'?'?观察上式可以发现，对于给定的地震动？?？结构的地震反应仅与结构的阻尼比？和自 振频率？3有关。换句话说，对于大小尺寸不同的结构，当结构阻尼比和自振频率相同时, 对同一个地震的反应完全相同。当阻尼比较小时，？'=?)?,则结构地震反应式可以简化为1?= XTT / 2找邛?1?"?-%?1? ?)? ? ? ? ? ? (26) ? 0则最大位移反应为：.一1? . . X 一一?= I? ?= ? |

17、 / ?"" ? ?)? ? (27), ? 0?质点相对于地面的速度为：?= = r ?2?夕?？?-?2(? ? ?QQQQ ?' ?,?0? ? ”? ? (28)? ? ? (29)+ - / ?,?(?-?(? ?)? ? ? 0则质点相对于地面的最大速度反应为:?= I?= I / ? " " " " " " .?0? 7? . a 一+ ?夕?3d？ ？?？?(？2 ?!)? ?马？ 0在实际工程中，我们比较看重的是结构的绝对加速度?+ ?*?)从单指点体系的在地震作用下的运动方程可以知道:?

18、+ ?超?？= -2?？- ?2? 一？ , 一。 一+ - / ?宠？?？?>4?(? ?)? ? (30) ? 0将上面求得的？?？、？?)?弋入上式:?+ ?>?= 2?f ?夕?>4 ?-?彳？ ？?! ?，?02?， ?- / ?夕?;(?-?(?2 ?p?0?2 ?+ ?彳？?夕?4?-?2 ?)? ? ? ? (31) ?. ,? 0由于低阻尼体系，？?"?方 则上式可以简化为?+ ?(?= ?/ ?>，(?-?(?- ?)? ? ? (32)0则点相对于地面的最大加速度反应为：皿?= |?)?+ ?"?|?= ?| /?/?即?-?；

19、(? ?)? ? (33)0?ELC波以上，我们求得了 ？洒？私?他们具有如下关系?= ?3北?二当阻尼比给定时，结构对任一地震的最大相对位移反应、最大相对速度和最大绝对加速度反应仅由？?）决定。工程中，一般习惯采用结构的自振周期？?= 2?/?弋替圆频率，因而工程谱中使用的反应谱一般以自振周期为自变量。EL CENTRO 1940(N5七=% 02C三0【。N秒）A 2 =326gil，msEkv<iir<，1940 (N Si 地震 心液料加速度反应谱绝对加忌度反应诣加通速“算法这囤尼汨阻出国二首相同有国后时拟和遨度段必牌小于胸加速度及皮猎，在怖昵时这 神差京旅不E月显，高阳同

20、时比兢阻界T而种反南晋的M津店都隔朋 尼增大而涉h且反曲a的波动雌、蝇剧解，隋延尼比的用 大而暹设.速度-位将景大时，速度为零，此固在覆值上二者相等言的自据周期里番t m增大.报贡田道 不同阻尼上匕之间的差舁 逐渐减少，到任周搦睡 在重合.随着T的增大，绝对加速度反应唱不回 阻尼此间的差异逐渐遍小.在Tf 左 右於近1零而到T?4S的长周期朦谶 差异又开始增大最后趋于平铝h并且高 阻尼的痂:岩花低俎后凶反应沿之上.1曲的作周期琳，西柏辰府酒苗Ml谋席都盼明刚汴而弑小绝对加速度谱相对速度反谱口 (原来)5 Tf秒）位移反应谱10.3.3影响地震反应谱的因素影响地震反应谱的因素有两点，一是体系阻尼

21、比，二是地震动。(1)体系阻尼比(2)地震动不同的地震动将会有不同的地震反应谱，大家都知道地震动特性的三要素： 振幅、频谱、 持时。振幅越大，地震反应谱值越大， 但是地震动振幅仅对地震反应谱值有影响，且呈线性比 例关系。地面运动的各种频率与加速度幅值的对应关系场地越软，地震震中距越大，地震反应谱峰值对应的周期也越长。所以说，地震动频谱对地震反应谱的 形状有影响。持时对地震反应谱影响不大。10.3.4地震反应谱的特点(1)阻尼比对反应谱影响很大，不仅能降低结构反应的幅值，而且可以削平不少峰 值，使反应谱曲线变得平缓。(2)对于加速度反应谱，当结构的周期小于某个值时，幅值随周期急剧增大，大于某 个

22、值时，快速下降，当?> 3?时，加速度反应谱值下降缓慢。(3)对于速度反应谱，当结构周期小于某个值时幅值随周期增大，随后趋于常数。(4)对于位移反应谱，幅值随周期增大。(5) 土质条件对反应谱形状有很大的影响，土质越松软，地震反应谱峰值对应的周期 也越长。小结：结构的最大地震反应，对于高频结构来说，主要取决于地面运动的加速度；对 于中频结构来说，主要取决于地面运动的最大速度；对于低频结构来说，只要取决于地面 运动的最大位移。10.3.5设计反应谱地震反应谱直接用于结构的抗震设计有一定的困难，所以需要有专门的可供结构抗震设计用的反应谱，我们称之为设计反应谱。?= ?於？ = ?|?盟? ?

23、纵？?|?猾？?？?(? ? ? (33)(1)地震系数？?= 1?可以将地震动振幅对地震反应谱的影响分离出来。基本烈度67n地震系数k(LO5OJO (045 )()2()(OJO )烈度每增加一度，地震系数大致增加一倍。(2)体系最大加速度的放大系数?=?(?I?初?= 2?1?I / ?蜜?2?(?-?)?2?(? ? ? ? ? (33) |?,?F? 0,?那么，对于给定的加速度记录 ？?？和结构白阻尼比?那么可以用上式计算出对应不 同的结构自振周期？的动力系数？?直。?= |?+ ?(?|?= ?| /?-?(? ?)? ? ? (34)0?对比两个公式，可以看出，地面最大加速度？

24、?寸于给定的地震是一个常数，所以 ？?-?的曲线形式与拟加速度反应谱的曲线的形状是完全一致的，只是纵坐标的数值不同。？?-?岫线的纵坐标为石口，而拟加速度反应谱的纵坐标是？私L"? ?回顾一下，地震反应谱是现阶段计算地震作用的基础，通过反应谱把随时程变化的地震作用转化为等效侧向力。对于单自由度体系，把惯性力看作反映地震对结构体系影响的等效力，用它对结构进行抗震验算。叫?？ ？?(？?= | ?(?= ?= ?一?'?一=?=? ? ? ? ? (35)?1?彳”?= ?水平地震影响系数。上面讲过，？? ?的曲线形式与拟加速度反应谱的曲线的形状是完全一致的，只是纵坐标的数值不同，由于对于给定的地震(或者设防烈度)，地震系数？为常数，所以呢，？? ?的曲线形式与拟加速度反应谱的曲线的形状是完全一致的，其纵坐标为?I?但是，由于地震的随机性，同一地点、相同地震烈度的前后两次地震记录的