振型分解反应谱法子程序.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)
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/S A/S, A/C17/JDW(150,3)
DO 10 1=1,JD 4JD是振动位移个数
10 JW⑴=JDW(JW(1),1) 4JW 是振动位移对应的节
点号，计算完后是振动位移对应的自由度编号
NN=MAX0(JDW(N,l),JDW(N,2),JDW(N,3)) 'N
是总可动节点数，JWD是节点定位向量
C 形成柔度矩阵
DO 50 II=1,JD
DO 35 11=1,NN
35 P(Il>0
P(JW(II))=1 '在第JW(II)个自由度上施加单位 荷载
CALL JFC(LD,2)
DO 40 L=1,JD
40 A(L,II)=W(JW(L)) 'WO 是位移，取出运动自
由度上的位移，存入柔度系数矩阵A ()
50 CONTINUE
雅可比法求解
CALL JACOBI(JD,PW)
DO 500 11=1,IP，对振型循环
C 求各振型的地震荷载
A(II,2)=A(II,2)*CT 'CT：周期折减系数
IF() THEN(根据设防烈度LE,求水平 地震影响系数最大值AX)
IFDHY=1 THEN AX= 多遇地震
IF DHY=2 THEN AX= 罕遇地震
ENDIF
IF() THEN
AX=
IF DHY=1 THEN
IF DHY=2 THEN AX=
ENDIF
IF() THEN
IF DHY=1 THEN AX=
IF DHY=2 THEN AX=
ENDIF
IF(A(II,2).) THEN(取出自振周期，规范只
能最大考虑6秒，大于6秒按6秒计)
PP=6
ELSE
PP=A(II,2)
ENDIF
IF()THEN,计算地震影响系数a
AF=(*PP+.45)*AX
GOTO 310
ENDIF
IF() THEN 'TU 是场地类型
IF()AF=(.2/PP)**.9*AX 近震(老抗震
规范，分近震和远震)
IF()AF=(.3/PP)**.9*AX 远震
ENDIF
IF() THEN
IF()AF=(.3/PP)**.9*AX
IF()AF=(.4