数值分析底部剪力法与振型分解反应谱法对比分析

底部剪力法与振型分解反应谱法对比分析摘要：建筑结构抗震设计是建筑结构设计中必不可少，也是非常重要的一部分。

结构抗震在建筑结构的总成本中占有相当大的比例。

建筑抗震设计规范中有关于结构抗震计算的方法以及适用范围，水平地震力的计算方法主要是底部剪力法和振型分解反应谱法，底部剪力法适用于质量和刚度沿高度分布比较均匀的结构，而振型分解反应谱法能反应结构的真实情况，对一般结构都适用。

本文通过对五层、八层、十层，质量和刚度分布均匀和不均匀框架结构的各层剪力计算，来比较两种方法的计算结果，验证底部剪力法的适用范围以及有效性。

本文对结构特征周期的计算是用广义Jacobi方法，通过Fortran语言编程实现的。

关键词：底部剪力法；振型分解反应谱法；Jacobi方法；Fortran语言Comparative Analysis between Equivalent Base Shear Method and ModalAnalysis MethodAbstract: Seismic design plays an essential and important part in the structure design. It also makes up a significant proportion of the total cost. About the horizontal seismic force, the code has detailed specification of the calculation principle and applicable scope. The calculation method for horizontal seismic force mainly has the equivalent base shear method and modal analysis method. The equivalent base shear method is suitable for mass and stiffness along the height of structure with uniform distribution, and the modal analysis method reflects the true action of the structure and has a wide usage. By calculating the shear of five-story, eight-story and ten-story framework with mass uniform or non-uniform distribution, this paper verified the scope and the effectiveness of the equivalent base shear method. The eigenperiod of the structure is calculated by generalized Jacobi method though the Fortran language programming.Key words: Equivalent Base Shear Method; Modal Analysis Method; Jacobi Method; Fortran Language引言实际的建筑结构其质量一般均是连续分布的，因此，严格的说，其动力自由度均是无限的。

但采用无限自由度模型一方面计算过于复杂；另一方面也没有这种必要，因为选用有限自由度模型的计算结果已能充分满足一般情况下工程设计的精度要求[1]。

因此，对于多高层房屋，一般均采用多自由度模型。

对于建筑结构而言，一般每层楼面及屋面可作为一个质点，而楼面与楼面（屋面）之间墙、柱的质量则分别向上、向下集结到楼面及屋面质点处。

这种多自由度模型即为工程上的层间模型。

图1所示为这种层间模型的计算简图，这种简图即所谓的“具有n个质量的悬臂柱”，本文所讨论的即“悬臂柱具有n个质量，比较基底剪力法和振型分解反应谱法的结果”。

关于建筑结构水平地震力的弹性计算方法，建筑抗震设计规范中[2]的5.1.2条规定，各类建筑结构的抗震计算，应采用的方法中：（1）高度不超过40m、以剪切变形为主且质量和刚度沿高度分布比较均匀的结构，以及近似于单质点体系的结构，可采用底部剪力法等简化方法。

（2）除（1）款外的建筑结构，宜采用振型分解反应谱法。

那么，规范为什么要如此规定呢，具体有哪些方面的考虑。

本文将通过对框架结构的两种方法计算的结果进行对比分析，说明规范规定的原因及考虑，这里只考虑规则结构，也就是不计结构扭转耦联。

1 底部剪力法对于一般的建筑结构，应采用振型分解反应谱法计算其地震作用效应，但当房屋结构满足下述条件时，可采用更为简便的底部剪力法计算其地震作用效应[1]：（1）结构的质量和刚度沿高度分布比较均匀；（2）房屋的总高度不超过40m；（3）建筑结构在地震作用下的变形以剪切应变为主；（4）建筑结构在地震作用时的扭转效应可忽略不计。

满足上述条件的结构在地震作用下其反应通常以第一振型为主，且近似为直线。

图2 结构水平地震作用计算简图建筑抗震设计规范文献5.2.1条规定：采用底部剪力法时，各楼层可仅取一个自由度，结构的水平地震作用标准值，应取下列公式确定（图2）：1EK eqF Gα=（1）1(1)i ii EK nni ijG HF FG Hδ==-∑（i=1,2,3…n）（2）n n EK F F δ∆= （3）其中，EK F 是结构总水平地震作用标准值；1α是相应于结构基本自振周期的水平地震影响系数，根据规范按图3进行计算；eq G 是结构等效重力荷载，单质点取总重力荷载代表值，多质点可取总重力荷载代表值的85%；i F 是质点i 的水平地震作用标准值；i G 是集中于i 质点的重力荷载代表值；iH 是i 质点的计算高度；n δ是顶部附加地震作用系数按表1；n F ∆是顶部附加水平地震作用。

图3 地震影响系数曲线[2]2 振型分解反应谱法适用范围：除上述可以采用底部剪力法以外的建筑结构，宜采用振型分解反应谱法。

建筑抗震设计规范5.2.2条规定：采用振型分解反应谱法时，不进行扭转耦联计算的结构，应按以下规定计算其地震作用和作用效应。

表1 顶部附加作用影响系数[2]注：1T 为结构基本自振周期。

（1）结构第j 振型i 质点的水平地震作用标准值，应按下列公式计算：ji j j ji i F X G αγ= （i=1,2,…n,j=1,2,…m ）（4）211=/n nj ji i ji i i i X G X G γ==∑∑ （5）其中，ji F 是j 振型i 质点的水平地震作用标准值；j α是相应于第j 自振周期的地震影响系数；ji X 是j 振型i 质点的水平相对位移；j γ是j 振型的参与系数。

（2）水平地震作用效应（弯矩、剪力、轴向力和变形），当相邻振型的周期比小于0.85时，可按下式确定：EK S =其中，EK S 是水平地震作用标准值的效应；j S 是j 振型水平地震作用标准值的效应，可只取前2~3个振型，当基本自振周期大于1.5s 或是房屋高宽比大于5时，振型个数应当适当增加。

3 用Jacobi 方法求解基本振型3.1 经典Jacobi 方法的计算基本步骤经典Jacobi 方法就是用一系列平面转换矩阵，从实对称矩阵A 出发，逐次作相似变换而使之对角化，则此对角阵的对角元素就是矩阵A 的全部特征值，而且所作正交变换的正交矩阵乘积的各个列向量即为矩阵A 的相应的特征向量[3]。

经典Jacobi 方法的基本步骤如下： （1）选定A k -1的非对角元素的绝对值最大的元素-1k pq a （）；（2）计算sinθ和cosθ的值-1k pq x a=-（）（6）-1(1)1()2k k qqpp y a a -=-（） （7）(sign y ω= （8）sin 2θω= （9）sin θ=（10）于是平面旋转变换矩阵为：（3）计算变换之后的A k由于对A 进行平面旋转变换[]Q 时，只是改变A 的p 行p 列和q 行q 列，所以只需计算A k 的p 行p 列和q 行q 列元素即可。

()(1)2(1)2(1)cos sin sin 2k k k k ppppqqpqaaaaθθθ---=++ （11）()(1)2(1)2(1)sin cos sin 2k k k k qq pp qq pq a a a a θθθ---=+- （12）()()(1)(1)(1)()sin cos cos2k k k k k pq qp qq pp pq a a a a a θθθ---==-+ （13）()()(1)(1)cos sin (,)k k k k pi ip ip iq a a a a i p q θθ--==+≠ （14）()()(1)(1)sin cos (,)k k k k qi iq ip iq a a a a i p q θθ--==-+≠ （15）（4）计算特征向量X k()(1)(1)cos sin (1,2...)k k k ip ip iq x x x i n θθ--=+=（16）()(1)(1)sin cos (1,2...)k k k iq ip iq x x x i n θθ--=-+= （17）()(1)(,1,2...;,,)k k ij ij x x i j n i j p q -==≠（18）3.2 广义Jacobi 方法的计算基本步骤经典Jacobi 方法只适用于求解A 为实对称矩阵的标准特征值问题Ax x λ=的特征值和特征向量。

而广义Jacobi 方法可以用来求解广义特征值问题Ax Bx λ=的全部特征值和特征向量。

在广义特征值问题中，A 为实对称矩阵，B 为正定近似对角矩阵[4]。

本例中，需要求解质量分布不均匀悬臂柱的特征值和特征向量，所以必须采用广义Jacobi 方法来做。

广义Jacobi 方法的基本原理和经典Jacobi 方法类似：（1）按次序（1，2），（1，3），…，（1，n ），（2，3），…，（n -1，n ）的顺序选取（i ，j ）的值。

用广义Jacobi 矩阵代替经典Jacobi 方法中的平面转换矩阵作合同的矩阵序列。

[]1cos sin sin cos q 1p Q θθθθ⎡⎤⎢⎥⎢⎥⎢⎥←⎢⎥=⎢⎥⎢⎥-←⎢⎥⎢⎥⎢⎥⎣⎦行行1111k i P jαβ⎡⎤⎢⎥⎢⎥⎢⎥←⎢⎥=⎢⎥⎢⎥←⎢⎥⎢⎥⎢⎥⎣⎦1T k kk k A P A P +=，1Tk kk k B P B P +=，0,1,2,k=（2）令()()()()k k k k ii ij ii ij a a b b a =-，()()()()k k k k ii jj ii jj b a b b a =-，()()()()k k k k jj ij jj ij c a b b a =-，(2b d sign b =+。