底部剪力法和振型分解法比较分析.docx

1、底部剪力法和振型分解法比较分析 SHANGHAI UNIVERSITY土木工程数值方法及应用课程论文COURSE THESIS题 目:底部剪力法和振型分解法比较分析学 院 土木工程系专 业 土木工程学 号 *学生姓名 曹宝雅指导教师 朱杰江 教授日 期 2011/03/18上海大学20102011学年冬季学期研究生课程论文课程名称： 土木工程数值方法及应用 课程编号： 181101903 论文题目: 底部剪力法和振型分解法比较分析 研究生姓名: 曹宝雅 学 号: * 论文评语:成 绩: 任课教师: 朱杰江 教授 评阅日期: 底部剪力法和振型分解法比较分析曹宝雅（上海大学 土木工程系）摘要：为验

2、证抗震设计规范对底部剪力法和振型分析法的适用条件和范围的规定，本文用Fortran编写Jocabi法程序，对三层均匀、三层非均匀、四层均匀、四层非均匀、八层均匀和十一层均匀结构进行了比较分析，证实了规范的正确性，并给出了新的规律。关键词：抗震；底部剪力法；振型分解法；Fortran；Jocabi法Comparative analysis between bottom shear method and modal analysis methodCao Baoya(Department of Civil Engineering, Shanghai Univercity)Abstract: In o

3、rder to vertify the rule of application conditions and range about bottom shear method and modal analysis method in earthquake resistant design code, this article has compared and analysed three-storey uniform, three-storey nonuniform, four-storey uniform, four-storey nonuniform, eight-storey unifor

4、m and eleven storey uniform structure by Jocabi method through Fortran program. The result has confirmed correctness of earthquake resistant design code. And it also give some new regularity. Keywords: seismic resistance; bottom shear method; model analysis method; Fortran; Jocabi method 1 三种抗震分析方法的

5、比较11.1底部剪力法适用范围：对于重量和刚度沿高度分布比较均匀、高度不超过40m，并以剪切变形为主（房屋高宽比小于4时）的结构，振动时具有以下特点；(1)位移反应以基本振型为主；(2)基本振型接近直线。 基本原理：在振型分解反应谱法的基础上，针对某些建筑物的特定条件做进一步简化，而得到的一种近似计算水平地震作用的方法：将多自由度体系简化成单自由度体系，计算出结构总的地震作用（即结构底部剪力），再将其按倒三角形原则分配到各个楼层，计算结构内力。1.2振型分解反应谱法适用范围：除上述底部剪力法外的建筑结构。基本原理：利用振型分解法的概念，把多自由度体系分解成若干个单自由度体系振动的组合，并利用单

6、自由度体系的反应谱理论计算各个振型振动的地震作用，最后将各个振型计算出的地震效应按一定的规则组合起来，求出总的地震响应。1.3时程分析法适用范围：抗震规范规定，重要的工程结构，例如：大跨桥梁，特别不规则建筑、甲类建筑，高度超出规定范围的高层建筑应采用时程分析法进行补充计算。基本原理：时程分析法是对结构物的运动微分方程直接进行逐步积分求解的一种动力分析方法。由时程分析可得到各质点随时间变化的位移、速度和加速度动力反应，并进而可计算出构件内力的时程变化关系。2实例分析某8层钢筋混凝土框架（如图1所示），集中于楼盖和屋盖处的重力荷载代表值为G1=G2=G12=1200kN，柱的截面尺寸为800mmx

7、800mm，采用C40的混凝土，E=32.5kN/mm2，梁的刚度EI=，建筑场地为类，抗震设防烈度7度，设计地震分组为第二组，设计基本地震加速度为0.1g，结构阻尼比=0.05。图1 某8层钢筋混凝土框架计算简图2.1振型分解反应谱法求结构的各周期和振型2.1.1建立振动微分方程（1）柔度法：A= ：特征值 ：特征向量N层均匀结构柔度矩阵通式：（2）刚度法： B= :特征值 ：特征向量N层结构刚度矩阵通式：Kn= 若k1=k2=k3=kn2.1.2经典Jacobi法计算周期和振型2只能用于实对称矩阵计算（1）经典Jacobi法计算基本步骤:选定Ak-1非主对角线元素的最大值确定计算Ak计算X

8、k2.1.3 Fortran程序计算3（1）程序内容: DIMENSION A(8,8),V(8,8),P(8) DOUBLE PRECISION A,V,P N=8DATA P/8*1/ DO 10 I=1,N-1 A(I,I)=P(I)+P(I+1)4 A(I,I+1)=-P(I+1) A(I+1,I)=A(I,I+1) A(N,N)=P(N)10 CONTINUE EPS=0.0001 CALL CJCBI(A,N,EPS,V,L) IF (L.NE.0) THEN WRITE(*,20) (A(I,I),I=1,N) END IF20 FORMAT(1X,8D13.4) WRITE(*

9、,*) WRITE(*,20)(V(I,J),J=1,N),I=1,N) END SUBROUTINE CJCBI(A,N,EPS,V,L) DIMENSION A(N,N),V(N,N) DOUBLE PRECISION A,V,FM,CN,SN,OMEGA,X,Y INTEGER P,Q L=1 DO 20 I=1,N V(I,I)=1.0 DO 10 J=1,N IF (I.NE.J) V(I,J)=0.010 CONTINUE20 CONTINUE25 FM=0.0 DO 30 I=2,N DO 30 J=1,I-1 IF (ABS(A(I,J).GT.FM) THEN FM=ABS(

10、A(I,J) P=I Q=J END IF30 CONTINUE IF (FM.LT.EPS) THEN L=1 RETURN END IF IF (L.GT.1000) THEN L=0 RETURN END IF L=L+1 X=-A(P,Q) Y=(A(Q,Q)-A(P,P)/2.0 OMEGA=X/SQRT(X*X+Y*Y) IF (Y.LT.0.0) OMEGA=-OMEGA SN=1.0+SQRT(1.0-OMEGA*OMEGA) SN=OMEGA/SQRT(2.0*SN) CN=SQRT(1.0-SN*SN) FM=A(P,P) A(P,P)=FM*CN*CN+A(Q,Q)*SN

11、*SN+A(P,Q)*OMEGA A(Q,Q)=FM*SN*SN+A(Q,Q)*CN*CN-A(P,Q)*OMEGA A(P,Q)=0.0 A(Q,P)=0.0 DO 60 J=1,N IF (J.NE.P).AND.(J.NE.Q) THEN FM=A(P,J) A(P,J)=FM*CN+A(Q,J)*SN A(Q,J)=-FM*SN+A(Q,J)*CN END IF60 CONTINUEDO 70 I=1,N IF (I.NE.P).AND.(I.NE.Q) THEN FM=A(I,P) A(I,P)=FM*CN+A(I,Q)*SN A(I,Q)=-FM*SN+A(I,Q)*CN END

12、IF70 CONTINUE DO 80 I=1,N FM=V(I,P) V(I,P)=FM*CN+V(I,Q)*SN V(I,Q)=-FM*SN+V(I,Q)*CN80 CONTINUE GOTO 25 END（2）对8层钢筋混凝土框架结构程序的运行结果如下图所示：图2 该8层结构的特征值、特征向量程序计算结果2.1.4振型组合取前三阶振型进行计算T1= 0.80783 T2= 0.27234 T3= 0.16719表1 前三阶振型的特征向量X10.089120.17520.25540.32680.38710.43420.46650.483X20.25540.43420.4830.38710.

13、1752-0.08912-0.3268-0.4666X30.38710.46650.1752-0.2554-0.483-0.32690.089120.4342结构j振型i质点的水平地震作用标准值，按下式计算4式中Fjij振型i质点的水平地震作用标准值；j相应于j振型自振周期的地震影响系数；Xjij振型i质点的水平相对位移；rjj振型的参与系数。 水平地震作用效应 图3 地震影响系数表2 水平地震作用效应组合层数F1iF2iF3iV1iV2iV3iV/2M111.8964 20.8920 18.0992 349.3551 69.7506 22.7720 178.4886 446.2214 223

14、.3844 35.5320 21.8116 337.4587 48.8586 4.6728 170.5046 426.2616 334.0889 39.5255 8.1963 314.0743 13.3265 -17.1388 157.4119 393.5297 443.6188 31.6777 -11.9368 279.9854 -26.1989 -25.3351 141.1737 352.9343 551.6672 14.3372 -22.5831 236.3665 -57.8766 -13.3984 121.8589 304.6472 657.9538 -7.2954 -15.2845

15、184.6993 -72.2138 9.1847 99.2636 248.1589 762.2783 -26.7431 4.1678 126.7455 -64.9183 24.4692 72.2454 180.6134 864.4672 -38.1752 20.3014 64.4672 -38.1752 20.3014 38.8121 97.0303 2.2底部剪力法（1）基本原理4:表3 顶点附加地震系数顶点附加地震作用系数Tg(s)T11.4TgT11.4Tg0.350.08T1+0.070.00.550.08T1-0.02（2）该8层钢筋混凝土框架结构计算结果：表4 该8层结构底部剪力法

16、计算结果FiViVi/2M5.0225206.6112103.3056258.264010.0450201.5887100.7944251.985915.0675191.543795.7719239.429720.0900176.476388.2381220.595325.1124156.386378.1932195.482930.1349131.273965.6369164.092335.1574101.138950.5695126.423740.179965.981532.990882.47692.3该8层结构底部剪力法和振型分解法比较表5 两种方法比较八层均匀（ki=k,i=1,8；m）

17、层号振型分解反应谱法（KN）底部剪力法（KN）相对误差(%)1178.49 173.39 -2.86 2170.50 168.93 -0.92 3157.41 160.02 1.65 4141.17 146.65 3.88 5121.86 128.82 5.71 699.26 106.53 7.32 772.25 79.79 10.45 838.81 48.59 25.20 3其它比较对结构层数发生变化，结构均匀性发生变化的情况下进行了类似的计算比较，如表6-表9所示：表6 三层均匀、非均匀三层均匀(k1=k2=k3=k；m)三层非均匀（k1=k2=k3=k；m）层号振型分解反应谱法（KN）底

18、部剪力法（KN）相对误差(%)振型分解反应谱法（KN）底部剪力法（KN）相对误差（%）1132.08 122.40 -7.33 136.48 122.40 -10.32 2105.79 102.75 -2.87 103.95 102.70 -1.20 360.17 63.45 5.45 60.15 63.31 5.24 表7 四层均匀、非均匀四层均匀(k1=k2=k3=k4=k；m)四层非均匀（k1=k2=k3=k4=k；m）层号振型分解反应谱法（KN）底部剪力法（KN）相对误差(%)振型分解反应谱法（KN）底部剪力法（KN）相对误差(%)1161.82 153.16 -5.35 186.25

19、 163.20 -12.37 2141.66 138.52 -2.21 153.19 147.50 -3.71 3106.31 109.25 2.77 106.84 116.11 8.68 458.44 65.34 11.81 55.56 69.03 24.24 表8 八层均匀八层均匀（ki=k,i=1,8；m）层号振型分解反应谱法（KN）底部剪力法（KN）相对误差(%)1178.49 173.39 -2.86 2170.50 168.93 -0.92 3157.41 160.02 1.65 4141.17 146.65 3.88 5121.86 128.82 5.71 699.26 106.

20、53 7.32 772.25 79.79 10.45 838.81 48.59 25.20 表9 十一层均匀十一层均匀(ki=k,k-1,12；m)层号振型分解反应谱法（KN）底部剪力法（KN）相对误差(%)1187.71 181.74 -3.18 2182.02 179.26 -1.52 3172.53 174.29 1.02 4161.40 166.83 3.36 5149.54 156.89 4.92 6136.69 144.46 5.69 7122.47 129.55 5.78 8106.41 112.15 5.39 987.22 92.26 5.78 1063.10 69.89 10

21、.76 1133.42 45.04 34.76 4数据分析与结论4.1数据分析（1）结构越不均匀，两种计算方法误差越大，如表6、表7所示。（2）高度越高，误差越大，如图4所示。（3）在三层（15m）以下，振型分解法算出的内力较大；三层以上，基底剪力法算出的内力较大；三层时，误差最小。如图4所示。图4 高度不同的均匀结构计算结果比较4.2 结论从以上分析结果可以看出，对于高度小于40m（本例中八层以下）的均匀结构，底部剪力法和振型分解法计算结果具有较好的吻合性；而对于高度高于40m或非均匀结构，两种算法结果差别较大。因此，规范对高度小于40m的均匀规则结构采用底部剪力法，而对于一般的高层或非均匀

22、结构采用振型分解法的规定是经济合理的。此外，我们还可以看到，对均匀结构，在15m（三层）的高度时，两种方法误差最小；在小于15m时，振型分解法算出的内力较大；而在大于15m时，基底剪力法算出的内力较大。问题：本文采用的是经典Jocabi法进行的数值分析，经典Jocabi法只用用于对称矩阵的计算，对于质量变化的非均匀结构，刚度法公式中的将是非对称的，这就要采用其它的数值分析方法来计算，如QR算法等。参 考 文 献1金春福.浅谈建筑结构抗震设计方法J.中国科技信息,2006,(1).2同济大学计算数学教研室.现代数值数学和计算M.上海：同济大学出版社，2004.219-234.3程誉，王惇，翁卫兵.Fortran语言程序设计M.北京：中国电力出版社，2002.4GB 50011-2010 建筑抗震设计规范S.