振型分解反应谱法子程序.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,NN35 P(I1)=0 P(JW(II)=1在第JW(II)个自由度上施加单位荷载 CALL JFC(LD,2) DO 40 L=1,JD40 A(L,II)=W(JW(L)W（）是位移，取出运动自由度上的位移，存入柔度系数矩阵A（）50 CONTINUEC 雅 可 比 法 求 解 CALL JACOBI(JD,PW) DO 500 II=1,IP对振型循环C 求 各 振 型 的 地 震 荷 载 A(II,2)=A(II,2)*CTCT：周期折减系数 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 ENDIF300 IF (AF.GT.AX) AF=AX310 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,NN330 P(I)=0 DO 340 I=1,JD340 P(JW(I)=AF*XD*S(I,II)*PW(I)（） CALL JFC(LD,2)解方程 CALL DYL(M,0,T