振型分解反应谱法子程序.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
显示全部