平面桁架杆单元fortran程序.docVIP

PROGRAM BAR3D DIMENSION MEA(100,4),JZ(50,2),AAI(10),AEU(10,2),R(2,6),RT(6,2) DIMENSION AKEP(2,2),AKE(6,6),P(300),X(100),Y(100),Z(100) DIMENSION LD(300),IS(6),A(50000),UPE(2),PPE(2),UE(6),RA(6,2) OPEN(1,FILE=FNAME1_t.txt) OPEN(2,FILE=FNAME2_t.txt) READ(1,*) NN,NE,NM,NA,NC write(*,*) NN,NE,NM,NA,NC NF=3 ND=2 NFD=NF*ND N=NN*NF DO 50 I=1,NN 50 READ(1,*)K,X(I),Y(I),Z(I),(P(NF*(I-1)+J),J=1,NF) READ(1,*)((JZ(I,J),J=1,2),I=1,NC) DO 100 I=1,NE 100 READ(1,*) IE,(MEA(I,J),J=1,4) READ(1,*) (AAI(I),I=1,NA) READ(1,*) ((AEU(I,J),J=1,2),I=1,NM) CLOSE(1) WRITE(2,460) NN,NE,NM,NA,NC WRITE(2,465) (I,X(I),Y(I),Z(I),P(3*I-2),P(3*I-1),P(3*I),I=1,NN) WRITE(2,470) ((JZ(I,J),J=1,2),I=1,NC) write(2,475) (I,(Mea(i,j),j=1,4),i=1,Ne) WRITE(2,480) (I,AAI(I),I=1,NA) WRITE(2,485) (I,(AEU(I,J),J=1,2),I=1,NM) 400 FORMAT(5I5) 405 FORMAT(I5,6F10.0) 410 FORMAT(2I5) 415 FORMAT(5I5) 420 FORMAT(F10.0) 425 FORMAT(2F10.0) 460 FORMAT(/5X,THE INPUT OF NN,NE,NM,NA,NC//(5X,5I5)) 465 FORMAT(/5X,THE INPUT OF X,Y,Z,P//(5X,I5,6F10.2)) 470 FORMAT(/5X,THE INPUT OF JZ//(5X,I5,5X,I5)) 475 FORMAT(/5X,THE INPUT OF MEA//(5X,I5,5X,4I5)) 480 FORMAT(/5X,THE INPUT OF AAI//(5X,I5,5X,F20.6)) 485 FORMAT(/5X,THE INPUT OF AEU//(5X,I5,5X,2E15.6)) CALL FLD(NN,NE,MEA,NF,ND,N,NT,LD) DO 500 I=1,NT 500 A(I)=0.0 DO 600 IE=1,NE CALL KEP(NN,NA,IE,MEA,AEU,AAI,X,Y,Z,AKEP) ! CALL CR(NN,IE,MEA,X,Y,Z,R) CALL TRAN (2,6,R,RT) CALL DOT(6,2,2,RT,AKEP,RA) CALL DOT(6,2,6,RA,R,AKE) CALL FIS(IE,MEA,NF,ND,NFD,IS) DO 560 I=1,NFD DO 560 J=1,NFD IF (IS(I)-IS(J)) 560,520,520 520 NI=IS(I) IJ= LD(NI)-NI+IS(J) A(ij)=A(ij)+AKE(I,J) 560 CONTINUE 600 CONTINUE CALL FCC(NC,N,NT,NF,JZ,LD,A) CALL DECOM(N,NT,A,P,LD) WRITE(2,850) (I,P(3*I-2),P(3*I-1),P(3*I),I=1,NN) WRITE(2,860) DO 800 IE=1,NE CALL KEP(NN,NA,IE,MEA,AEU,AAI,X,Y,Z,AKEP) CALL CR(NN,IE,MEA,X,Y,Z,R) CALL FIS(IE,MEA,NF,ND,NFD,IS) DO 750 I=1,NFD NI=IS(I) UE(I)=P(NI) 750 CONT