fortran語(yǔ)言編寫(xiě)及有限元結(jié)構(gòu)程序-資料下載頁(yè)
2025-08-21 01:52本頁(yè)面
【正文】 DIMENSION JE(2,NE),X(NJ),Y(NJ) REAL*8 BL,SI,CO,DX,DY J1=JE(1,M) J2=JE(2,M) DX=X(J2)X(J1) DY=Y(J2)Y(J1) BL=DSQRT(DX*DX+DY*DY) SI=DY/BL CO=DX/BL RETURN ENDC (八)形成單元?jiǎng)偠染仃? SUBROUTINE ESM(M,NE,BL,EA,EI,KD) DIMENSION EA(NE),EI(NE) REAL*8 KD(6,6),BL,S,G,G1,G2,G3 G=EA(M)/BL G1=*EI(M)/BL G2=*G1/BL G3=*G2/BL DO 10 I=1,6 DO 10 J=1,610 KD(I,J)= KD(1,1)=G KD(1,4)=G KD(4,4)=G KD(2,2)=G3 KD(5,5)=G3 KD(2,5)=G3 KD(2,3)=G2 KD(2,6)=G2 KD(3,5)=G2 KD(5,6)=G2 KD(3,3)=*G1 KD(6,6)=*G1 KD(3,6)=G1 DO 20 I=1,5 I1=I+1 DO 30 J=I1,630 KD(J,I)=KD(I,J)20 CONTINUE RETURN ENDC (九)形成單元坐標(biāo)轉(zhuǎn)換矩陣 SUBROUTINE CTM(SI,CO,T) REAL*8 T(6,6),SI,CO DO 10 I=1,6 DO 10 J=1,610 T(I,J)= T(1,1)=CO T(1,2)=SI T(2,1)=SI T(2,2)=CO T(3,3)= DO 20 I=1,3 DO 20 J=1,320 T(I+3,J+3)=T(I,J) RETURN ENDC (十)形成單元固端力 SUBROUTINE EFX(I,NPF,BL,PF,FO) DIMENSION PF(4,NPF) REAL*8 FO(6),A,B,C,G,Q,S,BL IND=PF(2,I) A=PF(3,I) Q=PF(4,I) C=A/BL G=C*C B=BLA DO 5 J=1,65 FO(J)= GO TO (10,20,30,40,50,60,70),IND10 S=Q*A* FO(2)=S*(*G+C*G) FO(5)=S*G*() S=S*A/ FO(3)=S*(*C+*G) FO(6)=S*C*(*C) GO TO 10020 S=B/BL FO(2)=Q*S*S*(+*C) FO(5)=Q*G*(+*S) FO(3)=Q*S*S*A FO(6)=Q*B*G GO TO 10030 S=B/BL FO(2)=*Q*C*S/BL FO(5)=FO(2) FO(3)=Q*S*(*S) FO(6)=Q*C*(*C) GO TO 10040 S=Q*A* FO(2)=S*(*G+*G*C) FO(5)=S*G*(*C) S=S*A FO(3)=S*(*C+*G)/ FO(6)=S*C*(*C) GO TO 10050 FO(1)=Q*A*(*C) FO(4)=*Q*C*A GO TO 10060 FO(1)=Q*B/BL FO(4)=Q*C GO TO 10070 S=B/BL FO(2)=Q*G*(*S+C) FO(5)=FO(2) S=S*B/BL FO(3)=Q*S*A FO(6)=Q*G*B100 RETURN END
規(guī)章制度相關(guān)推薦
鄂ICP備17016276號(hào)-1