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

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 AV/SA/S ,A/C 17/JDW( 150,3) DO 10 1=1,JD 4JD是振动位移个数 10 JW(I)=JDW(JW(I),1) 4JW是振动位移对应的节 点号，计算完后是振动位移对应的自由度编号 NN=MAX0(JDW(N, 1),JDW(N,2),JDW(N,3)) N 是总可动节点数，JWD是节点定位向量 C 形成柔度矩阵 DO 50 II=1,JD DO 35 11=1,NN 35 P(I1)=O 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 AX=0.08多遇地震AX=0.5 AX=0.08多遇地震 AX=0.5罕遇地震 AX=0.16 AX=0.9 AX=0.32 AX=1.4 雅可比法求解 CALL JACOBI(JD,PW) DO 500 11=1,IP 对振型循环 C 求各振型的地震荷载 A(II,2)=A(II,2)*CT CT：周期折减系数