fortran語言編寫及有限元結構程序(編輯修改稿)
2025-09-17 01:52 本頁面
【文章內(nèi)容簡介】 ,3),J=1,NJ) READ (5,*) ((JE(I,J),I=1,2),EA(J),EI(J),J=1,NE) IF () READ (5,*) ((PJ(I,J),I=1,2),J=1,NPJ) IF () READ (5,*) ((PF(I,J),I=1,4),J=1,NPF) WRITE (6,10) NE,NJ,N,NW,NPJ,NPF WRITE (6,20) (J,X(J),Y(J),(JN(I,J),I=1,3),J=1,NJ) WRITE (6,30) (J,(JE(I,J),I=1,2),EA(J),EI(J),J=1,NE) IF () WRITE (6,40) ((PJ(I,J),I=1,2),J=1,NPJ) IF () WRITE (6,50) ((PF(I,J),I=1,4),J=1,NPF)10 FORMAT(/6X,39。NE=39。,I5,2X,39。NJ=39。,I5,2X,39。N=39。,I5,2X,39。NW=39。,I5,2X, $39。NPJ=39。,I5,2X,39。NPF=39。I5)20 FORMAT(/7X,39。NODE39。,7X,39。X39。,11X,39。Y39。,12X,39。XX39。,8X,39。YY39。,8X,39。ZZ39。/ $(1X,I10,3I10))30 FORMAT(/4X,39。ELEMENT39。,4X,39。NODEI39。,4X,39。NODEJ39。,11X,39。EA39。,13X,39。EI39。/ $(1X,3I10,))40 FORMAT(/7X,39。CODE39。,7X,39。PXPYPM39。/(1X,))50 FORMAT(/4X,39。ELEMENT39。,7X,39。IND39。,10X,39。A39。,14X,39。Q39。,/ $(1X,))C (二)形成總結點荷載向量 DO 55 I=1,N55 P(I)= IF() GO TO 65 DO 60 I=1,NPJ L=PJ(1,I)60 P(L)=PJ(2,I)65 IF() GO TO 90 DO 70 I=1,NPF M=PF(1,I) CALL SCL(M,NE,NJ,BL,SI,CO,JE,X,Y) CALL EFX(I,NPF,BL,PF,FO) CALL CTM(SI,CO,T) CALL EJC(M,NE,NJ,JE,JN,JC) DO 75 L=1,6 S= DO 80 K=1,680 S=ST(K,L)*FO(K) F(L)=S75 CONTINUE DO 85 J=1,6 L=JC(J) IF() GO TO 85 P(L)=P(L)+F(J)85 CONTINUE70 CONTINUEC (三)形成整體剛度矩陣90 DO 95 I=1,N DO 100 J=1,NW100 KB(I,J)=95 CONTINUE DO 105 M=1,NE CALL SCL(M,NE,NJ,BL,SI,CO,JE,X,Y) CALL CTM(SI,CO,T) CALL ESM(M,NE,BL,EA,EI,KD) CALL EJC(M,NE,NJ,JE,JN,JC) DO 110 I=1,6 DO 115 J=1,6 S= DO 120 L=1,6 DO 125 K=1,6125 S=S+T(L,I)*KD(L,K)*T(K,J)120 CONTINUE KE(I,J)=S115 CONTINUE110 CONTINUE DO 130 L=1,6 I=JC(L) IF() GO TO 130 DO 135 K=1,6 J=JC(K) IF()GO TO 135 JJ=
規(guī)章制度相關推薦
鄂ICP備17016276號-1