桥梁结构.doc

174 页 DIMENSION IA(NMTOT),IC(NMTOT) C一一 SET STORAGE TO ZERO DO 100 I=1, NMTOT !给整型矩阵冲零 IA(I)=0 100 IC(I)=0. 0 RETURN END C-----------------------RFILE SUBROUTINE RFILE(LM,ASA,RF,AMSA,ND) IMPLICIT INTEGER*4 (I-N) IMPLICIT REAL*8 (A-H,O-Z) DIMENSION ASA(ND,ND),LM(ND),RF(ND),AMSA(ND,ND) COMMON/IOLIST/ NTM, NTR, NIN, NOT, NSP, NFL, NT7, NT8, NT9 C------------------------- !从临时文件 K3 中读入数据 READ(NT7)N,ND, (LM(I), I=1,ND),((ASA(I,J),J=1, ND), I=1, ND). 1 (RF(I) , I=1,ND).((AMSA(I,J),J=1, ND) , I=1, ND) RETURN END C一一一—一 RFILEI SUBROUTINE RFILEI ( LM, GMSA, ND ) IMPLICIT INTEGER*4 (I-N) IMPLICIT REAL*8 (A-H,O-Z) DIMENSION GMSA(ND,ND),LM(ND) COMMON/IOLIST/ NTM, NTR, NIN, NOT, NSP, NFL, NT7, NT8, NT9 C一一一一一一一一一一—!从临时文件 N3 中读入数据 READ (NT9)(LM(I),I=1, ND),(( GMSA(I,J)J=1, ND) , I=1,ND) RETURN END C------------------------------------------FORMID SUBROUTINE FORMID(ID,NUMNP,NRJ,NEQ) IMPLICIT INTEGER*4 (I-N) IMPLICIT REAL*8 (A-H, O-Z ) DIMENSION ID(NUMNP , 3), II(3),NM(3) COMMON/IOLIST/ NTM,NTR,NIN,NOT,NSP,NFL,NT7,NT8,NT9 C-------------SET ID ARRAY------------------------- DO 100 N =1,NUMNP DO 100 I=1,3 175 页 ID(N,I)=0 100 CONTINUE NM(1)=0 C一一 READ RESTRAINT INFORMATION 一一!读入约束信息 CALL FIND( ′ REST ′ ,KL)! 寻找 SYST 分隔符, KL =0表示找到该分隔,否则暂停 IF(KL. EQ. 0) GOTO 200! WRITE(NOT,1999) WRITE(NTM,1999) STOP 200 WRITE(NOT,1000) WRITE(NTM,1000) DO 105 I=1, 3 II( I)=0 NM(I)=0 105 CONTINUE 600 CALL FREEa CALL FREEPT CALL FREEI( ′′ ,NM(1),3)! 将三个整型数读入 NM(1),NM(2),NM(3) IF( NM (1) . EQ. 0) GOTO201 !当 NM(1)=0 时表示数据段已读完 NL=NM(1) !NL 为具有相同约束信息的一串结点号中最小结点号码 NH=NM(2) !为这一串结点号中最大结点号码 INC=NM(3) !读入约束信息,为这一串结点号的增量 IF(INC. LE. 0)INC=1 IF( NH. LE. 0) NH = NL CALL FREEI( ′ R′ , II,3)!读入引导字 R后的三个整型数 DO 160 L=NL, NH , INC DO 160 I=1,3 IF(II(I).) ID(L,I)=1 !形成初始 ID表 160 CONTINUE GO TO 600 C一一 EVALUATION OF EQUATION NUMBERS —!计算有效自由度总数 NEQ 201 NEQ=0 DO 400 N=1,NUMMP DO 350 I=1, 3 IF(ID(N, I). NE. 0) GOTO 355 !形成最终 ID表 NEO=NEQ 十 1!计算 NEQ 总数 ID( N, I)= NEQ GO TO 350 176 页 355 ID(N, I)=0 350 CONTINUE 400 CONTINUE