振型分解反应谱法.docx

1、附录一振型分解反应谱法振型分解反应谱法作为弹性多自由体系的主要分析方法，很有必要对振型分解反应谱法有充分的了解。 本文仅作为大家参考之用，如有理解上的错误或者不当，敬请谅解。1 、单自由度体系在地震作用下的运动如图（ 1）所示，根据达朗贝尔原理有：f c f I f s 0（1）也即： mucu kumug（2）方程两边同时除以m ，可化为：图（ 1）u 2 u2 uug（ 3）式中，2k / m ，令c，为体系阻尼比。2m2 、多自由度体系在地震作用下的运动类似于单自由度体系分析过程，体系运动方程为： m u c u k u mug（4）无阻尼体系自由振动时，ug0， c0 ，上式即为：m

2、uk u 0（5）根据方程解的特征，设其解的形式为： u sin(t)（ 6）代入（ 5）式有：( k2 m)sin(t) 0（ 7）由于 sin(t) 0则 ( k2 m) 0（ 8）另外， 0 ，故特征方程为：k2 0（ 9）m由（ 9）式可以求出2 ，进而可以求得各阶振型对应的圆频率2i ，再代入（ 8）式可求对应于各个2的特征向量i ，即为振型。i振型：多自由度体系自由振动时，各质点在任意时刻位移比值是一定的，不随时间变化，即体系自由振动过程中形状保持不变。振型是结构形状保持不变的振动形式，唯一的。 N 个自由度的体系具有N 个振型。则结构的变形总可以表示成这N 个振型的线性组合：Nu

3、qi ii 1其中 qi 称为正则坐标。3、振型的正交性由于 k 2 0m则 k r r2 m r 0（ 12）式两边同时左乘 n T ， (nr ) ，得到：n T kr r2 n T mr 同理， r T kn n2 r T mn，该式两边同时转置一次，得到：n T kr n2 n T mr （ 13），（14）两式左右对应相减，得到：（ r2n2 )n T mr 0(rn)22因为 rn所以 n T m r 0(rn)同理亦有 n T k r 0(rn)即所说的振型关于质量和刚度矩阵满足正交性质。对于阻尼：根据瑞雷阻尼的基本假定，若用矩阵形式表达，即：camb k由于该式是线性表达式，根

4、据前面推导的振型正交性质，可以得出：n T c r 0(r n)但要注意的是， 体系振型关于质量和刚度矩阵满足正交性质是无条件的，振型的形状是（10）（ 11）（12）（13）（ 14）（15）（ 16）（17）（ 18）（19）而振型关于阻尼矩阵满足正交性质却是有一定条件的，阻尼不满足正交的情况下就不能在理论上严格的对结构进行振型分解来求解。3.1振型正交性质的物理意义振型关于质量矩阵的正交性：第 n 阶振型的惯性力在经历第r 阶振型时所做的功为0；振型关于刚度矩阵的正交性：与第 n 阶振型位移有关的弹性力在经历第r 阶振型时所做的功为0；总体来说就是各个振型按照自己的规律振动，而相互之间没

5、有干扰，从而为把方程求解分解成按各个振型分别求解提供了可能。4、具有经典阻尼的多自由度体系在P(t)sp(t) 激励下的反应这种荷载形式是一个固定的分布向量乘以常数， 具有振型类似的特征， 但是实际结构是很少会受到这种形式的荷载作用， 但它对于后面要进行的地震作用下的响应分析是有重要作用的。4.1运动方程求解（后面u 与等参数都是表示向量）如图（ 2）所示，得到体系的平衡方程： mu cuk usp(t)（ 20）由公式（ 10），把 ui(t) iqi (t ) 代人，则：NNN m(i qi )c(i qi ) k (i qi ) s p(t )（ 21）i1i1i 1方程前乘 nT，根据

6、振型的正交性质，上式可化为：M n qn Cn qnK n qnnT sp(t)（22）其中 M nnT m n称为广义质量CnnTcn称为广义阻尼K nnTkn称为广义刚度上式两边同时除以M n ，由2K n，并另 nC n得到：n2 n M nM nqn2nn qnn2 qnnT s p(t)M n令nnTs，称为振型参与系数。M n则（ 23）式化为：qn2 nn qnn2qnn p(t )n 是一个与振型正则化方式(即的取值 )及荷载分布向量有关的一个量， 它反映了外部激励对振型的影响程度。比如 s具有 m i的形式，根据正交性n 0(i n)或1(i n)，说明这种荷载只会引起第 i

7、 个振型的反应，而不会引起其他振型的反应。令： D nqn ，则（ 24）式可变为：图（ 2）（23）（24）图（ 3）nD n 2 n n D nn2 D n p(t)（25）该方程是如图（3）的一个标准的单自由度体系的运动微分方程，求解微分方程得到Dn (t ) 。4.2求激励下的效应对应于此刻第n 阶振型的位移 un (t ) nqn (t )nnD n (t) ，则如图（4）所示，对应于此时变形的弹性恢复力为：fn (t ) kun (t) k nn D n (t )n2 mnnD n (t )n m n n2D n (t )（26）N令 snn m n ，称之为振型贡献，且ssn ，

8、则（ 26)式可写为：n1f n (t)snn2 D n (t )（ 27）sn 是一个 n 阶列向量， 是个常量， 与结构特性及荷载作用方式都有关，而与振型的正则化方式无关。那么f n (t) 为一个静力部分与动力部分的组合，作用在自由度上的f n (t)由 sn 与n2D n (t ) 的乘积决定，则在f n (t) 作用下的效应也可以表示为由sn 作用下的效应经动力放大后的总效应，而这个弹性恢复力是用来计算结构各种反应效应的直接量。图（ 4）如图（ 5）所示，st 为在s作用下对应于结构需要求的作用效应值，可以是弯矩、剪nn力，位移等，则经动力放大后得到该效应的动力反应时程，可表示为：n

9、 (t )st2（ 28）nn D n (t )那么总的效应可以表示为：N(t )n (t)（ 29）n 1总静力效应：Nstst（ 30）nn1st图（ 5）令：n（ 31）nstn 称为振型贡献系数，则第n 阶振型的效应可表示为n (t )stn n2Dn (t ) 。4.2.1 对n 性质的讨论 n 是个无量纲的物理量，对于不同的反应量取值是不同；Ts n 与正则化方式无关， snn m n ，n，上下有两个n ， n 的影响nnT mn被约去；Nn1 ；i 1 n 有正有负，所以只计算前几个振型并不一定都是小于总效应的，也可能偏大的估计了荷载的作用效应。4.2.2 振型贡献系数与振型参

10、与系数的区别n 是一个与振型正则化方式无关而与结构特性和激励分布特性有关的量，是衡量各振型对反应量最直观的参数， 虽然它不包括动力部分， 不能精确反映振型对反应的贡献， 但已经能够表现各振型反应的相对大小。4.2.3 位移也是f n 作用下的一种效应st212112D n (t )n n D n (t )（ 32）un (t ) unn D n (t ) k sn n D n (t )2 mn m n nn这与之前的 un (t )n qn (t ) nn D n (t ) 的表达式是相同的，说明效应的表达式具有广泛的适用性。5、具有经典阻尼的多自由度体系在地震激励下的反应如图（ 6）所示，以

11、地面为参考系，则各质点受到的惯性力为 m ug ，它的分布与质量有关， 代表地面发生单位位移时各质点发生该方向上的位移，比如转动以及竖向自由度上的位移为 0，对于一般受水平方向作用的多自由度体系，常见的简化后的结构计算模型， 为1。5.1求解平衡方程及参数定义如图（ 6）所示，以地面为参考系，各力分量为：f I M 1u g ， f ccu ， f s ku由牛顿第二定律，可得：f Ifc fs mu ，（ 33）图（ 6）相对于地面， u 为相对于地面的加速度，则上式也可记为： mu cu kum1 ug（34）以地面为参考系，右边是体系受到的地震激励，即该参考系下的f I 。我们可以看到地

12、震激励虽然是一个复杂的随机过程，体现在u g ，但它的分布函数是 m1 ，是一个只与结构特性有关的量。我们联想到之前介绍的n ， sn 等都与激励分布函数有关，那么这部分的影响也一起整合到结构的特性里去了。类似于前面的分析过程：M n qnCn qnK n qnnT m1ug（35）两边同时除以 M n ，得到：qn 2 nn qnn2qnnT m1ug(36)M n令 nnT m 1，这个值与结构特性以及振型正则化方式有关，但却与外部激励无关。M n具有性质 :Nnn1（37）n 1则（ 36）式可记为：qn2 n n qnn2qnn ug（38）与前面相同，令D nqnn上式进一步可以化为

13、：Dn2D2 Dugn n nn n这里可以通过各种数值积分的方法求解该运动的微分方程，可以得到弹性力：fn (t) ku(t )n2 mu(t )n2 m n qn (t)n2 m n n D n (t )n mnn2 Dn (t )令 snn mnsn 具有的性质：Nm 1snn1则上式可以写为：fn (t ) snn2 D n (t)snAn (t)式中 An (t ) 称为伪加速度反应。nst 为 sn 作为静力作用下的结构反应。n (t )stAn (t )n（39）D n (t) 。（40）（41）（42）（43）（44）stst振型贡献系数：n (t )nn（45）Nststn

14、1nn (t)stn An (t)（46）st 对各阶振型计算都相同，而n 反应了各阶振型对静力部分贡献的相对大小，只与结构特性有关，只有An (t ) 是不确定的。5.2一般结构振型贡献系数的特征基本振型，即第一振型的振型贡献系数最大，结构越高，自由度越多，高阶振型的影响越大；高阶振型对力的影响要大于位移， 力在位移的基础上有2的系数， 而高阶振型的2nn很大， 加大了对力的影响， 更细致的关系是顶层剪力的影响要高于底部剪力的影响，高于顶点位移以及倾覆力矩的影响；框架梁端对柱的约束越大， 高阶振型影响越小， 即高阶振型对弯曲型结构的影响越大，这可能是底部剪力法考虑一个振型为什么只针对剪切型结

15、构的原因。6、振型分解反应谱法6.1伪加速度反应由（ 46）式看出结构的各阶反应时程与伪加速度反应的时程相关，如果我们求各阶振型的最大反应， 那么我们就可以根据伪加速度反应谱来求， 伪加速度反应谱提供了各个周期点的伪加速度的最大反应 An ，如图（ 7）所示。图（ 7）这样得到各阶振型的最大反应：stAn（ 47）nn6.2振型组合通过反应谱得到各阶振型的最大反应，但并不知道各阶振型最大反应的具体时刻，而它们的最大反应时刻几乎没有可能会相同，所以直接相加这些最大反应必然会过大的估计结构的最大反应。 所以我们需要合适的最大反应组合方式来合理的估计结构的最大反应，基于地震激励是个平稳的随机过程的基

16、本假定，得出了CQC 和 SRSS 的基本组合原则。6.2.1完全二次组合（ CQC 法）计算表达式如下：NNn )1/2(ini（48）i1n 11 ，当 i n 时，所有可能的组in 称为相关系数，反应了两个振型的相关程度，in合都要计算到。Der kiureghian 提出的相关系数表达式：8in (inin ) in3/2（49）in(1n2 )24 i n in (1in2) 4( i 2n2 ) in2其中inin如果在 in时：则in82 (1in )in3/2（50）(1in2 )242in (1in )2图（ 8）表示不同的下的inin关系：图（ 8）如图（ 8）所示，可以看

17、出的特征是随着阻尼比的减小，相关系数随周期比的增大退化的很快，在常见的结构阻尼比0.05 的情况下，周期相差两倍，相关系数已经接近于0，但如果结构的周期相当密集，则各振型的相关项不能忽略。一般情况下， 包括各种大跨结构，质量、刚度以及形状不规则的结构要考虑平扭耦合的影响可能会使得结构周期比较密集。6.2.2 平方开方法（SRSS 法）在非高阻尼化的周期较分散的情况下，则式（48）可写为：NNN(2in i n )1/2n（51）n 1i 1n 1i n则后面相关系数项近似为0，则（ 51）式可化为：Nn2 )1/2(（52）n 1本式就是平方开方法（SRSS），其实是一种退化以后的CQC 法。

18、6.2.3 振型组合的一些注意事项 SRSS 以及 CQC 法是在随机振动理论下得到的组合关系，是具有一定概率意义的对结构最大反应的包络， 从统计观点来看， 能够反应出结构反应量的大小水平。而对于某一次的地震激励， 并不一定实际的结果会与组合的结果相当，有可能偏大或偏小于组合结果，但一般情况下能够满足工程上的精度要求。求某一个反应量不能通过间接地来求，比如不能通过一种内力的峰值去求另外一种内力的峰值， 不能通过各层的位移峰值来求层间位移的峰值，而应该针对每一个反应量分别求出它各阶振型的峰值，然后组合得到该反应量的总反应的峰值。另外就是完全组合的时候，是全部的组合， 例如（ 48）式中的 i1,

19、n2 ，i2, n1都应该分别组合一次。6.3结合中国规范的振型分解反应谱法在实际的结构设计时， 我们不可能使用单条地震波的反应谱来计算， 而是使用综合了场地条件、地震环境等因素的有一定概率保证的设计反应谱。中国的振型分解反应谱法没有提出振型贡献系数的概念， 振型贡献系数是对反应量相对贡献大小最直观的一种表达， 是考虑多少个振型最有效地衡量标准， 中国规范没有给出定量的衡量所需振型个数的标准。6.3.1 公式推导为与中国规范表述一致，振型参与系数n 用i 表示，意义相同。同（ 39）式的动力微分方程：Di2 i i Dii2 Diug（53）可计算出：Di (t )ui (t )i qi i

20、i D i (t )（54）ni 与正则化方式有关，但根据（37）式具有i i1 的性质，证明见文献5 。那么第 i 个质点的位移，表达为各振型第 ii1质点的位移之和。nui (t )j D j (t ) ji（55）j 1这里的 ui (t ) 与（ 54）式中的 ui(t) 不相同，一个是表示第i 个质点，一个是表示第 i个振型。（ 55）式中的 ji 表示的是第 j 个振型的第 i 个质点的值。则相对加速度：n&&&（56）j D j (t ) jiui (t)j 1由（ 37）式的关系：&&&&n&&(j(t )

21、（57）ug (t) 1 ug (t)ji )ugj1则任一时刻的惯性力为：&&&&fimi (ui (t ) ug (t)n&&n&&mi (j D j (t)ji(t)j ji ugj 1j 1n&&&&nmij ji ( D j(t)(t )f jiugj 1j1则质点 i 在第 j个振型的惯性力为：&&&f jimijji(D j (t)ug (t)则质点 i 在第 j个振型的最大惯性力为：&&F jif jimaxmijji D j (t)ug (t

22、 ) maxmijjiSa(Tj)mi g jjijGijjijj 为地震影响系数，祥见文献5 。6.3.2 中国建筑抗震设计规范 2010版规定在不考虑扭转耦联的情况下：按式Fjijj X ji Gi其中nnjX ji Gi /X 2ji Gii1i 1（ 58）（ 59）（ 60）(61)把式（ 60）计算的地震作用加到结构上求作用效应（弯矩、剪力、轴向力和变形），当相邻振型的周期比小于0.85 时，按SRSS 法组合。当不满足时，规范5.2.2 条的条文说明中建议使用CQC法组合。当要考虑扭转耦联时，请参照建筑抗震设计规范5.2.3 条的规定。抗震规范在考虑需要计算多少个振型虽然在5.2

23、.2 条的条文说明里说一般取振型参与质量达到总质量的90%（有关振型参与质量将在后面介绍），但没有给出具体的计算式，只是在规范中说使用SRSS 法时，一般只取前23 个振型； 当基本自振周期大于1.5s 或房屋高宽比大于5 时，振型个数适当增加，而在抗震规范5.2.3 条中使用CQC组合方法时建议考虑915 个振型，但都没有给出具体的量化标准。6.3.3 相关问题讨论中国规范是按照加速度来推导惯性力，然后加到结构上来求效应，而之前是基于位移状态来推导内力。两则的出发点有所不同，结构的弹性恢复力并不与结构的惯性力相等，地震激励下弹性恢复力应该等于结构的惯性力和阻尼力之和。我们的理解是， 每阶振型

24、的最大反应时刻是位移的峰值点，也就是说位移的导数是0，即该振型结构的反应速度为0，根据瑞雷阻尼的假定，阻尼力也为0，从这个方面来说二者似乎是等效的，这方面的准确解释有待商榷。Sa(Tj ) 是加速度谱值，从（59）式的推导来看是真谱的概念，但之前的推导是用的伪谱，两者的出发点不同， 如果如果我们承认上述两种方法都是适用的， 那么也从另外一个侧面说明了真谱与伪谱的差别不大的事实，特别是加速度谱的短周期段。对于规范给出当自振周期大于1.5 或高宽比大于5 时，振型个数适当的增加的规定，虽然随着结构周期的加长，高阶振型的参与系数也加大，但从公式（47）可以看出结构真正的贡献大小不仅仅与静力部分的贡献

25、系数有关，还与An 有关，下面就讨论一下动力部分的影响。从图（ 9）与图（ 10）中可以看出随着结构周期的增大，高阶振型的反应量的动力部分相应的增大，关键点在图中标示出来，只要高阶振型的周期值跨过了该点进入平台段以后，则以后周期更短的 An 就不会再大于首次跨过该点的那个振型，也就是说我们考察动力部分影响的重点在没有跨过该特征点的振型上面。在我们看出， 随着周期的加大，结构高阶振型的动力部分会放大该振型的贡献，另外从图（10），图（ 11）看出结构的周期越密集，则高阶振型的动力部分影响也变显著，这也从另外一个侧面解释了周期加长以及周期密集应该多考虑一些振型的原因，虽然不一定是主要的因素。6.4

26、底部剪力法底部剪力法是仅考虑结构的基本振型的一种方法，不存在振型组合的问题。基于基本振型按着倒三角形的分布基本假定，推导了基本计算公式，详见抗震规范5.2.1 条。但之前介绍过随着结构周期的加长，高阶振型的影响变得更显著，特别是对于顶层的剪力，所以采用在顶部附加一个随周期加长而增大的作用力，来考虑高阶振型的影响。7、振型参与的相关参数商业软件一般在满足工程精度要求的基础上尽可能的减少所需的计算时间， 也就是减少要计算的振型数， 据了解 Etabs是根据结构的有效参与质量的大小来判断需要多少个振型的。图（ 12）如图（ 12）的三个自由度的体系为例。对于基底剪力，以第一阶振型为例：sts11s12 s13 M 1(62)1从（