振型分解反应谱法子程序.doc

SUBROUTINE DZL(M,LD,T,L2,IP,LE,TU,JD,CT,PW,JW,DHY) INTEGER TU,T,JW(50) REAL S(50,50),A(50,50),PW(50) REAL*8 W(450) COMMON N/C3/IJ(300,2)/C6/P(450)/C7 /W/SA/S,A/C17/JDW(150,3) DO 10 I=1,JD ‘JD是振动位移个数 10 JW(I)=JDW(JW(I),1) ‘JW是振动位移对应的节点号，计算完后是振动位移对应的自由度编号 NN=MAX0(JDW(N,1),JDW(N,2),JDW(N,3))‘N是总可动节点数，JWD是节点定位向量 C 形 成 柔 度 矩 阵 DO 50 II=1,JD DO 35 I1=1,NN 35 P(I1)=0 P(JW(II))=1‘在第JW(II)个自由度上施加单位荷载 CALL JFC(LD,2) DO 40 L=1,JD 40 A(L,II)=W(JW(L))‘W（）是位移，取出运动自由度上的位移，存入柔度系数矩阵A（） 50 CONTINUE C 雅 可 比 法 求 解 CALL JACOBI(JD,PW) DO 500 II=1,IP‘对振型循环 C 求 各 振 型 的 地 震 荷 载 A(II,2)=A(II,2)*CT‘CT：周期折减系数 IF(LE.EQ.7) THEN(根据设防烈度LE，求水平地震影响系数最大值AX) IF DHY=1 THEN AX=0.08 多遇地震 IF DHY=2 THEN AX=0.5 罕遇地震 ENDIF IF(LE.EQ.8) THEN IF DHY=1 THEN AX=0.16 IF DHY=2 THEN AX=0.9 ENDIF IF(LE.EQ.9) THEN IF DHY=1 THEN AX=0.32 IF DHY=2 THEN AX=1.4 ENDIF IF(A(II,2).GE.6) THEN(取出自振周期，规范只能最大考虑6秒，大于6秒按6秒计) PP=6 ELSE PP=A(II,2) ENDIF IF(PP.LE.0.1)THEN’计算地震影响系数α AF=(5.5*PP+.45)*AX GOTO 310 ENDIF IF(TU.EQ.1) THEN ‘TU是场地类型 IF(IC.EQ.1)AF=(.2/PP)**.9*AX 近震（老抗震规范，分近震和远震） IF(IC.EQ.2)AF=(.3/PP)**.9*AX 远震 ENDIF IF(TU.EQ.2) THEN IF(IC.EQ.1)AF=(.3/PP)**.9*AX IF(IC.EQ.2)AF=(.4/PP)**.9*AX ENDIF IF(TU.EQ.3) THEN IF(IC.EQ.1)AF=(.4/PP)**.9*AX IF(IC.EQ.2)AF=(.55/PP)**.9*AX ENDIF IF(TU.EQ.4) THEN IF(IC.EQ.1)AF=(.65/PP)**.9*AX IF(IC.EQ.2)AF=(.8/PP)**.9*AX ENDIF 300 IF (AF.GT.AX) AF=AX 310 SD=0 XD=0 DO 320 I=1,JD SD=SD+S(I,II)*PW(I) ‘S：振型向量，PW：重力代表值 320 XD=XD+S(I,II)*S(I,II)*PW(I) XD=SD/XD（振型参与系数） DO 330 I=1,NN 330 P(I)=0 DO 340 I=1,JD 340 P(JW(I))=AF*XD*S(I,II)*PW(I)（） CALL JFC(LD,2)解方程 CALL DYL(M,0,T,L2,LE)求单元力 CALL