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

SHANGHAI UNIVERSITY土木工程数值方法及应用课程论文COURSE THESIS题 目:底部剪力法和振型分解法比较分析学 院 土木工程系专 业 土木工程学 号 10721988学生姓名 曹宝雅指导教师 朱杰江 教授日 期 2011/03/18上海大学20102011学年冬季学期研究生课程论文课程名称： 土木工程数值方法及应用 课程编号： 181101903 论文题目: 底部剪力法和振型分解法比较分析 研究生姓名: 曹宝雅 学 号: 10721988 论文评语:成 绩: 任课教师: 朱杰江 教授 评阅日期: 底部剪力法和振型分解法比较分析曹宝雅（上海大学 土木工程系）摘要：为验证抗震设计规范对底部剪力法和振型分析法的适用条件和范围的规定，本文用Fortran编写Jocabi法程序，对三层均匀、三层非均匀、四层均匀、四层非均匀、八层均匀和十一层均匀结构进行了比较分析，证实了规范的正确性，并给出了新的规律。关键词：抗震；底部剪力法；振型分解法；Fortran；Jocabi法Comparative analysis between bottom shear method and modal analysis methodCao Baoya(Department of Civil Engineering, Shanghai Univercity)Abstract: In order 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 uniform 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 三种抗震分析方法的比较11.1底部剪力法适用范围：对于重量和刚度沿高度分布比较均匀、高度不超过40m，并以剪切变形为主（房屋高宽比小于4时）的结构，振动时具有以下特点；(1)位移反应以基本振型为主；(2)基本振型接近直线。 基本原理：在振型分解反应谱法的基础上，针对某些建筑物的特定条件做进一步简化，而得到的一种近似计算水平地震作用的方法：将多自由度体系简化成单自由度体系，计算出结构总的地震作用（即结构底部剪力），再将其按倒三角形原则分配到各个楼层，计算结构内力。1.2振型分解反应谱法适用范围：除上述底部剪力法外的建筑结构。基本原理：利用振型分解法的概念，把多自由度体系分解成若干个单自由度体系振动的组合，并利用单自由度体系的反应谱理论计算各个振型振动的地震作用，最后将各个振型计算出的地震效应按一定的规则组合起来，求出总的地震响应。1.3时程分析法适用范围：抗震规范规定，重要的工程结构，例如：大跨桥梁，特别不规则建筑、甲类建筑，高度超出规定范围的高层建筑应采用时程分析法进行补充计算。基本原理：时程分析法是对结构物的运动微分方程直接进行逐步积分求解的一种动力分析方法。由时程分析可得到各质点随时间变化的位移、速度和加速度动力反应，并进而可计算出构件内力的时程变化关系。2实例分析某8层钢筋混凝土框架（如图1所示），集中于楼盖和屋盖处的重力荷载代表值为G1=G2=G12=1200kN，柱的截面尺寸为800mmx800mm，采用C40的混凝土，E=32.5kN/mm2，梁的刚度EI=，建筑场地为类，抗震设防烈度7度，设计地震分组为第二组，设计基本地震加速度为0.1g，结构阻尼比=0.05。图1 某8层钢筋混凝土框架计算简图2.1振型分解反应谱法求结构的各周期和振型2.1.1建立振动微分方程（1）柔度法：M-12IX=0A= M =12：特征值 X：特征向量N层均匀结构柔度矩阵通式：A=l324EI11111 12222 12333 123n-1n-1 123n-1n （2）刚度法： K-2MX=0 M-1K-2X=0 B= M-1K =2:特征值 X：特征向量N层结构刚度矩阵通式：Kn= 若k1=k2=k3=knKn=24EIl32-10 -12-1 -12 -1 2-1 0-112.1.2经典Jacobi法计算周期和振型2只能用于实对称矩阵计算（1）经典Jacobi法计算基本步骤:选定Ak-1非主对角线元素的最大值apq(k-1)确定sin和cosX=-apq(k-1)Y=12aqq(k-1)-app(k-1)=sign(Y)XX2+Y2sin2=sin=2+1+1-2cos=1-sin2计算Akapp(k)=app(k-1)cos2+aqq(k-1)sin2+apq(k-1)sin2aqq(k)=app(k-1)sin2+aqq(k-1)cos2-apq(k-1)sin2apq(k)=aqp(k)=aqq(k-1)-app(k-1)sincos+apq(k-1)cos2api(k)=aip(k)=aip(k-1)cos+aiqk-1sin(ip,q)aqi(k)=aiq(k)=-aip(k-1)sin+aiqk-1cos(ip,q)计算XkXip(k)=Xipk-1cos+Xiqk-1sin(i=1,2,n)Xiq(k)=-Xipk-1sin+Xiqk-1cos(i=1,2,n)Xijk=Xijk-1i=1,2,n;jp,q2.1.3 Fortran程序计算3（1）程序内容: DIMENSION A(8,8),V(8,8),P(8)DOUBLE PRECISION A,V,PN=8DATA P/8*1/DO 10 I=1,N-1A(I,I)=P(I)+P(I+1)4A(I,I+1)=-P(I+1)A(I+1,I)=A(I,I+1)A(N,N)=P(N)10 CONTINUEEPS=0.0001CALL CJCBI(A,N,EPS,V,L)IF (L.NE.0) THEN WRITE(*,20) (A(I,I),I=1,N)END IF20FORMAT(1X,8D13.4) WRITE(*,*)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,YINTEGER P,QL=1DO 20 I=1,N V(I,I)=1.0 DO 10 J=1,N IF (I.NE.J) V(I,J)=0.010 CONTINUE20CONTINUE25FM=0.0DO 30 I=2,NDO 30 J=1,I-1 IF (ABS(A(I,J).GT.FM) THEN FM=ABS(A(I,J) P=I Q=J END IF30CONTINUEIF (FM.LT.EPS) THEN L=1 RETURNEND IFIF (L.GT.1000) THEN L=0 RETURNEND IFL=L+1X=-A(P,Q)Y=(A(Q,Q)-A(P,P)/2.0OMEGA=X/SQRT(X*X+Y*Y)IF (Y.LT.0.0) OMEGA=-OMEGASN=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*SN+A(P,Q)*OMEGAA(Q,Q)=FM*SN*SN+A(Q,Q)*CN*CN-A(P,Q)*OMEGAA(P,Q)=0.0A(Q,P)=0.0DO 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 IF60CONTINUEDO 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 IF70CONTINUEDO 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)*CN80CONTINUEGOTO 25END（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.1752-0.08912-0.3268-0.4666X30.38710.46650.1752-0.2554-0.483-0.32690.089120.4342结构j振型i质点的水平地震作用标准值，按下式计算4Fji=jjXjiGi（i=1,2,3,n；j=1,2,3，n）j=i=1nXjiGi/i=1nXji2Gi式中Fjij振型i质点的水平地震作用标准值；j相应于j振型自振周期的地震影响系数；Xjij振型i质点的水平相对位移；rjj振型的参与系数。 水平地震作用效应SEk=Sj2 图3 地震影响系数表2 水平地震作用效应组合层数F1iF2iF3iV1iV2iV3iV/2M111.8964 20.8920 18.0992 349.3551 69.7506 22.7720 178.4886 446.2214 223.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 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:FEk=lGeqFi=GiHii=1nGjHjFEk1-n(i=1,2,n)Fn=nFEk表3 顶点附加地震系数n顶点附加地震作用系数Tg(s)T11.4TgT11.4Tg0.350.08T1+0.070.00.550.08T1-0.02（2）该8层钢筋混凝土框架结构计算结果：表4 该8层结构底部剪力法计算结果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）层号振型分解反应谱法（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）底部剪力法（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 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.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