资源描述：
《振型分解反应谱法子程序》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、SUBROUTINEDZL(M,LD,T,L2,IP,LE,TU,JD,CT,PW,JW,DHY)INTEGERTU,T,JW(50)REALS(50,50),A(50,50),PW(50)REAL*8W(450)COMMONN/C3/IJ(300,2)/C6/P(450)/C7/W/SA/S,A/C17/JDW(150,3)DO10I=1,JD‘JD是振动位移个数10JW(I)=JDW(JW(I),1)‘JW是振动位移对应的节点号，计算完后是振动位移对应的自由度编号NN=MAX0(JDW(N,1),JDW(N,2)
2、,JDW(N,3))‘N是总可动节点数，JWD是节点定位向量C形成柔度矩阵DO50II=1,JDDO35I1=1,NN35P(I1)=0P(JW(II))=1‘在第JW(II)个自由度上施加单位荷载CALLJFC(LD,2)DO40L=1,JD40A(L,II)=W(JW(L))‘W（）是位移，取出运动自由度上的位移，存入柔度系数矩阵A（）50CONTINUEC雅可比法求解CALLJACOBI(JD,PW)DO500II=1,IP‘对振型循环C求各振型的地震荷载A(II,2)=A(II,2)*CT‘CT：周期折减系数
3、IF(LE.EQ.7)THEN(根据设防烈度LE，求水平地震影响系数最大值AX)IFDHY=1THENAX=0.08多遇地震IFDHY=2THENAX=0.5罕遇地震ENDIFIF(LE.EQ.8)THENIFDHY=1THENAX=0.16IFDHY=2THENAX=0.9ENDIFIF(LE.EQ.9)THENIFDHY=1THENAX=0.32IFDHY=2THENAX=1.4ENDIFIF(A(II,2).GE.6)THEN(取出自振周期，规范只能最大考虑6秒，大于6秒按6秒计)PP=6ELSEPP=A(II
4、,2)ENDIFIF(PP.LE.0.1)THEN’计算地震影响系数αAF=(5.5*PP+.45)*AXGOTO310ENDIFIF(TU.EQ.1)THEN‘TU是场地类型IF(IC.EQ.1)AF=(.2/PP)**.9*AX近震（老抗震规范，分近震和远震）IF(IC.EQ.2)AF=(.3/PP)**.9*AX远震ENDIFIF(TU.EQ.2)THENIF(IC.EQ.1)AF=(.3/PP)**.9*AXIF(IC.EQ.2)AF=(.4/PP)**.9*AXENDIFIF(TU.EQ.3)THENIF(
5、IC.EQ.1)AF=(.4/PP)**.9*AXIF(IC.EQ.2)AF=(.55/PP)**.9*AXENDIFIF(TU.EQ.4)THENIF(IC.EQ.1)AF=(.65/PP)**.9*AXIF(IC.EQ.2)AF=(.8/PP)**.9*AXENDIF300IF(AF.GT.AX)AF=AX310SD=0XD=0DO320I=1,JDSD=SD+S(I,II)*PW(I)‘S：振型向量，PW：重力代表值320XD=XD+S(I,II)*S(I,II)*PW(I)XD=SD/XD（振型参与系数）DO
6、330I=1,NN330P(I)=0DO340I=1,JD340P(JW(I))=AF*XD*S(I,II)*PW(I)（）CALLJFC(LD,2)解方程CALLDYL(M,0,T,L2,LE)求单元力CALLWRITEDYL(M,T,L2,IP,II)输出单元力500CONTINUECALLETQ(M,T,L2,IP)END
