大学课件_建筑结构抗震设计2.4 振型分解反应谱法.pptx

1、第二章结构抗震基本知识2.1 计算原则2.2 地震作用2.3 设计反应谱2.4 振型分解反应谱法2.5 底部剪力法2.6 结构竖向地震作用2.7 结构抗震验算2.4 振型分解反应谱法大多数结构物都应简化为多质点体系分析。而振型分解反应谱法是弹性体系地震反应的基本方法。其基本概念：(1)假定建筑结构是线弹性多自由度体系；(2)利用振型分解，变为求解n个独立的等效单自由度弹性体系的最大地震反应，从而求得每一振型的作用效应；(3)按SRSS或CQC法则进行作用效应组合。振型分解法只需考虑前几阶振型，减小计算量。一、不考虑扭转影响时结构的地震作用和作用效应对大多数质量和刚度分布比较均匀和对转影响，可在

2、建筑物的两个主轴方向分别考虑水平地震作用进行验算。称的结构不需要考虑水平地震作用下的扭()()()()()()1.多自由度弹性体系的运动方程()图2-12多自由度弹性体系位移1 2 n M diag m m m=L,在n质点即n个自由度的弹性体系：M为质量矩阵，一般采用集中质量阵形式：K刚度矩阵，nn阶对角矩阵，如果只考虑层间剪切变形的层间剪切结构K为三对角矩阵&Mx(t)+Cx(t)+Kx(t)=MIxg(t)（2-24）（2-25）C阻尼矩阵C=M+K瑞利阻尼形式（2-26）定；I单位列向量。其中 、由1、2、1 、2 确 2多自由度弹性体系的自由振动将式（2-24）略去阻尼项和右端项，振

3、动方程：Mx(t)+Kx(t)=0&（2-27）设（2-27）式的解为：x(t)=Xsin(t+)&初相角。（2-28）x(t)=2Xsin(t+)=2x(t)（2-29）其中：(K M)X=0将式(2-28)、(2-29)代入式(2-27)，得：2(2-30)为了体系振动，必须是非零解，则：XK2M=0(2-31)该方程的n个根2 21、2 Ln即是体系的n个自振频率，一般有：2则n个自振周期：T1 T2 L Tj L Tn将所求的 j 依次代回（2-30），可得到与之相对应的Xj，即为振型。=0K22 m2 X 2mK11 1 12KK22 m2K21可解出 一个两自由度体系：体系的自由振

4、动方程为：(即式(2-30)K12 X12K11 m1 2 K21(2-32)(2-33)频率方程(即式(2-31)2=022 21，2，将之带回(2-32)式K11 1 m2式（2-32）是齐次方程组，两个方程线性相关2 2第一振型：第二振型：X11X12X 21X 22值：K122 K12K11 2m1=每一振型的幅值之比都是常数，不随时间而变。(2-34a)(2-34b)X MXk3振型的正交性(1)振型关于质量矩阵的正交性：其矩阵表达式为：=0Tjj k(2-35)式(2-35)是根据功的互等定理推导而来，式中Xj ，Xk 分别为体系第 j、k 振型的振幅向量。物理意义 ：某一振型在振

5、动过程中所引起的惯性力不在其它振型上作功，说明某一个振型的动能不会转移到其它振型上去，也就是体系按某一振型作自由振动时不会激起该体系其它振型的振动。X KXk =0(2)振型关于刚度矩阵的正交性其矩阵表达式为(由功的互等定理而来)：Tj(j k)(2-36)2物理意义：该体系按k振型振动引起得弹性恢复力在j振型位移所作的功之和等于零，也即体系按某一振型振动时，它得位移不会转移到其它振型上去。X CXkC =X CXj(3)振型关于阻尼矩阵的正交性令：C=M+K（2-26）则有：=0Tj(j k)当j=k时，j振型的广义阻尼为：Tj j（2-37）（2-38）()()()X1 1 =0.618X

6、 1 2 =1.0X2 2 =1.0X2 1 =-1.6 18图2-13两质点弹性体系的振型（a）两质点弹性体系；(b)第一振型；(c)第二振型m K例题2-3已知某两个质点的弹性体系（如图），其结构参数为：1=m2=m，1=K2=K 。验算质量矩阵和刚度矩阵的正交性。M=K=K X=质量矩阵和刚度矩阵分别为：m00m K2K K0.618 1 1.618 1 X2第一振型和第二振型分别为：X MX2=0.618 1.000解：（1）验算振型关于质量矩阵的正交性由式（2-61）得：01.618m 1.000 m0T11.618m 1.000m X KX2=0.618 1.000（2）验算振型关

7、于刚度矩阵的正交性由式（2-68）得：K1.618K 1.000 2K KT1 4.236K 2.618K =x x xX1 11 12 13=x x xX2 21 22 234.振型分解n个自由度的弹性体系具有n个独立振型。例如：3个自由度体系。TTTXji:j振型i自由度的相对位移。（2-39）x=21x23由(2-39)形成矩阵：A=X1 X2 X3x11 x12 x13x22x31 x32 x3333（2-40）xi(t)=xji j(t)按照振型叠加原理，弹性结构体系中每一个自由度在振动过程中的位移 可以表示为：qnj=1q式中：j(t)为j振型的广义坐标,即以振型作为坐标系的位移值

8、。(2-41)xi(t)t(=1L n =整个结构体系的位移列向量、速度列向量和加速度列向量可分别表示为:X2 x 1(t)x2(t)X X M xn(t)&x(t)=Aq(t)x(t)=Aq(t)q 1q2X A qM qn(2-42a)(2-42b)(2-42c)对方程两端左乘 得：A5.计算水平地震作用的振型分解反应谱法地震作用下的运动微分方程为：&Mx(t)+Cx(t)+Kx(t)=MIxg(t)(2-24)将(2-42a)、(2-42b)、(2-42c)式代入上式并T ()()A M A q t A C A q t+&()g A K A A M I x t+=&()()j j j j

9、 j j M t C X q t+&XX X ()()j g j j j K q t M t+=&XX X IT TT Tq(t)利用正交性，对式(2-43)化简，展开后可得n个独立的二阶微分方程，对于第j振型有：T TT T(2-43)(2-44)M =XjM j*X其中：广义质量：TjTjTj则有：Tj j j j(2-45)和 的关系：KC =2 j jM K j =jM j 、*j仿照单自由度体系，M*jC*j(2-46)j j 2 其中，j 表示第j振型的阻尼的圆频率。将上式代入式(2-45)，并用j振型的广义质量除等式的两端，得：T2j TX(j =1,2,L n)(2-47)j

10、j M要)，表示第i振型对结构总响应贡献的大小。其中：T 表示j振型参与系数(很重X X分析式(2-45)：相当于单自由度弹性体系的运动方程，区别在于：自由度弹性体系的运动方程进行解耦。此广义坐标 qj(t)作为未知量而不是x(t)；方程右端项多1个 j，经过坐标变换将n0 ()x e&()t t j ()sin ()()j g j j j q t t d t =参照单自由度体系，Duhamel积分，方程(2-49)式的解为：j j j j(t)为 j 、j 单自由度弹性体系的位移。(j=1,2,L,n)（2-48）xi(t)=j j(t)X ji&x&i(t)=j&j(t)X jij=1多自

11、由度弹性体系i质点相对于基础的位移和加速度为：nnj=1jnj=1X ji =1q qX ji 为振型幅值，j(t)、&(t)分别见式（2-41）和式（2-42）。由结构力学得：（2-49a）（2-49b）（2-50）x第i自由度t时刻的水平地震作用Fi(t)=作用在i质点的惯性力：(t)=m (t)+&nj=1则j振型i自由度得水平作用Fi(t)为：&x Fji(t)=mi j j(t)X ji+jX ji&g(t)（2-52）为单自由度体系的其中：是j振型自振周期的maxFji =Fji(t)&=mi jX jixg(t)+j(t)maxmaxx&g(t)+j(t)Sa(j,j)；j =S

12、a(j,j)/g影响系数。=j jX jiGi (i=1,2,L,n;j=1,2,L,n)(2-53)V ji =Fjk6.地震作用效应有了Fji后 计算j振型水平地震作用产生的作用效应。对于层间剪切结构，j振型地震作用下各楼层水平地震层间剪力：nk=i(i=1,2,Ln)(2-54)图2-15 j振型水平地震作用式(2-53)、(2-54)是整个地震作用过程中各振型作用效应的最大值。而最大值并不出现在同一时刻，不能简单叠加 存在一个各振型作用效应组合问题。规范规定了两种组合方法：平方和开方SRSS；完整的二次项组合法(CQC法)。SRSS法适用范围：平稳随机过程；各振型之间独立。jS 2S

13、=(2-55)例2-4.试用振型分解反应谱法计算人如图所示的三层框架，在多遇地震时的层间地震剪力。已知抗震设防烈度为8度，设计地震分组为第二组，类场地，阻尼比取0.05。321图2-16 三层框架示意第一振型 000.1 667.0 334.0 1 =X第二振型 000.1 666.0 667.0 2 =X 000.1 035.3 019.4 =X第三振型解：（1）求解结构体系的周期和振型由矩阵迭代法（或雅可比法）可计算出结构体系的三个自振周期和振型分别为：TT1=0.467sTT2=0.208sT3T3=0.134s1=0.40 2max =（2）计算各振型的地震影响系数查得多遇地震时设防烈

14、度为8度的 max =0.16查得类场地、设计地震分组为第二组Tg =0.40s当阻尼比 =0.05 时，2=1.0 =0.9第一振型，因Tg T1 5Tg ，所以：0.91.00.16=0.139 0.467Tg T第二振型2=max =0.16 第三振型 3=max =0.16X1764 000.1 2646 667.0 2646 334.0 +1 =1764 000.1 2646 667.0 2646 334.0 +1=iX 2 0.428 =第二振型：X i G2=1.363（3）计算各振型的参与系数由式（2-84）计算各振型的振型参与系数第一振型：31iGii=13 2 2 2X1i

15、Gi32iGii=1322ii=1 X G(4)计算各振型各楼层的水平地震作用各振型各楼层的水平地震作用由式(2-96)计算：第三振型：23iGiX3i i3i=13i=1 3=0.063第一振型：F1i =1 1X1iGiF11=0.1391.3630.3342646=167.4kNF12=0.1391.3630.6672646=334.4kNF13=0.1391.3631.0001764=334.2kN第二振型：F2i =2 2X 2iGiF21=120.9kN第三振型：F31=107.2kNF22=120.7kNF3i =3 3X 3iGiF32=80.9kNF23=120.8kNF33

16、=17.8kNF1k(5)计算各振型的层间剪力nk=i第一振型：V1i =V11=167.4+334.4+334.2=836.0kNV12=334.4+334.2=668.6kNV13=334.2kNV2i =F2k第三振型：V3i =F3k第二振型：V21=120.8kNnk=iV22=0.1kNV23=120.8kNnk=iV31=44.1kN V32=63.1kNV33=17.8kN（6）计算水平地震作用效应各层层间剪力由式（2-97）计算各层层间剪力2 2V1=V11+V21+V31=845.8kNV2=671.6kNV3=355.8kN图2-17 各振型的地震剪力(kN)（a）第一振

17、型地震剪力；（b）第二振型地震剪力；（c）第三振型地震剪力和组合后各层地震剪力(kN)(短边)1.15(长边)1.05扭转K小时乘数不小于1.3二 、估计水平地震作用扭转影响的结构地震作用和作用效应规则结构:由于施工、使用等原因所产生的偶然偏心有扭转现象；规范规定：规则结构不进行扭转耦联计算时，平行于地震作用方向的两个边榀，其地震作用效应增大系数：平面复杂、不规则、质量刚度明显不均匀、不对称的多高层建筑大量出现震作用下的扭转影响；过去采用“偏心矩法”，人为主观因素较多，物理概念不明确。规范：采用振型分解反应谱法计算平移-扭转耦联的多高层建筑的水平地震作用。规范规定：对这类建筑应考虑水平地要解决

18、以下三个问题：（1）求解平移扭转耦联体系的自由振动；（2）计算各振型水平地震作用标准值的表达式；（3）各振型地震作用效应的组合方法。1 平移扭转耦联体系的自由振动 楼板自身平面内绝对刚性，平面外刚度忽略不计；各榀抗侧力结构(框架或剪力墙)自身平面内刚度很大，平面外忽略；所有构件不考虑自身的抗扭作用；振动计算时，将质量(包括柱、墙的质量)集中到各层楼板处。基本假定：集中在每一楼层的质量有三个自由度，两个正交水平移动和一个转角，坐标原点一般选在各楼层的质心处，如图2-19所示，坐标轴为一折线形轴，运动方程：MD(t)+CD(t)+KD(t)&=MDg(t)&（2-56）图2-19计算简图Ji=mi

19、(a +b )/12其中：M 为广义质量阵（3n3n）M=diagm m Jm=diagm1 m2Lmn其中mi是第 i楼层的质量；J=diagJ1 J2LJn2 2上式 Ji表示第i层质量对本楼层质心的转动惯量。a、b分别是本楼层的长短边。（2-57）（2-58）（2-59）K=0KYYKYKxx=KxsK yy=KrrnK zpKK为广义侧移刚度阵：K zz 0K zpT KY（2-60）上式表明两个平动之间无耦联；nys=1其中ny 为平行于轴框架的榀数Kzs为平行于轴第s榀框架的刚度矩阵xr=1（2-61）（2-62）y x xs s r ynyK =Ks=1Ys=diagy1s y2

20、sLyns其中 y1s为第i层第s榀方向框架的y向坐标（2-63）（2-64）x图2-20第i层平面图K x=rrrnxK yr=1Xr=diagx1r x2rLxnr（2-65）（2-66）其中 xir为第i层第r榀y方向框架的x向坐标ny nxT TK Ks=1 r=1（2-67）D =u1其中D为广义位移，且：Tu2 L un v 1 v2 L vn 1 2 L n&Dg(t)地面运动水平加速度时间历程函数：T T T Tg g D Dd&g(t)为地面运动加速度的时间历程；D 为地面运动方向与x轴的夹角。其中：先求解平动扭转耦联体系的自由振动，其方程为：MD+KD=0&结构体系的自由度

21、为3n（n为结构的层数），常用雅可比法求解。（2-69）2 .结构体系考虑扭转影响的水平地震作用仍按振型分解法：&D(t)=Aq(t)D(t)=Aq(t)D(t)=Aq(t)A为振型矩阵。（2-70）&j&将式(2-70)代入运动方程（2-56），并利用振型正交性原理，可将方程（2-56）分解成为3n个相互独立的二阶微分方程，其通式为：q j+2 j jq j+2q j =jd g(t)(j=1,2,Ln,L3n)（2-71）Fyji =j tjYjiGi 经过与单向平移振动时相类似的推导，可以得到考虑扭转地震效应时第j振型第i层的水平地震作用标准值计算公式：Fxji =j tjX jiGi

22、Ftji =j tjri2 jiGi(i=1,2,Ln;j=1,2,Lm)（2-72）tj =X jiGi当仅考虑x方向地震时：ni=12 2 2ni=12 2 2(X ji+Yji +jiri2)Gi（2-73）当仅考虑y方向地震时：n n(X ji+Yji +jiri2)Gi （2-74）i=1 i=1当仅考虑与x方向斜交的地震时：（2-75）其中 xj、yj 分别为由式（273）和式（274）求得的参与系数3 考虑扭转作用时的振型组合问题采用用振型分解反应谱法，用式（272）一定的规则组合，以获得总的地震作用效应。计算各振型的水平地震作用 再计算每一振型水平地震作用产生的作用效应 各效应

23、按考虑扭转影响时体系振动的特点：扭转分量的影响不一定随高振型而衰减 增加振型的数可取915个。联严重 考虑振型之间的相关性；自由度增加 频率分布密集，振型耦1)+4 j k(1+T)T规范规定按下式组合（单向水平地震作用的扭转效应）：jkS jSkSEk =m mj=1 k=11.8 j k(1+T)T52 2 2T jk =（2-76）（2-77）SEk =S +(0.85Sx)关于 jk：当j=k时，T 1，jk 1，完全自相关。T增加，减小，说明低振型 与高振型之间相关性很小。双向水平地震作用时扭转效应，取maxxSEk =S 2+(0.85S y)222ySx或、S y 分别为x、y单向水平地震作用按式(276)计算的扭转效应。（2-78）（2-79）