文档仅供参考
L=I-K+1
C=AK(K,L)/AK(K,1)
IW=N_BAND-L+1
DO 10 J=1,IW M=J+I-K
10 AK(I,J)=AK(I,J)-C*AK(K,M)
20 RESULT_N(I)=RESULT_N(I)-C*RESULT_N(K)
RESULT_N(N_DOF)=RESULT_N(N_DOF)/AK(N_DOF,1)
DO 40 I1=1,N_DOF-1 I=N_DOF-I1
IF(N_BAND.GT.N_DOF-I-1)JQ=N_DOF-I+1 IF(N_BAND.LE.N_DOF-I-1)JQ=N_BAND DO 30 J=2,JQ K=J+I-1
30 RESULT_N(I)=RESULT_N(I)-AK(I,J)*RESULT_N(K) 40 RESULT_N(I)=RESULT_N(I)/AK(I,1) WRITE(8,50) 50
FORMAT(/12X,'* * * * * RESULTS BY FEM2D * * * * *',//8X,
OF
NODE--'//5X,'NODE
NO',8X,'X-
&'--DISPLACEMENT DISP',8X,'Y-DISP')
DO 60 I=1,N_NODE
60 WRITE(8,70) I,RESULT_N(2*I-1),RESULT_N(2*I)
21
2020年4月19日
文档仅供参考
70 FORMAT(8X,I5,7X,2E15.6) RETURN END
c
c calculate the stress components of element and node
SUBROUTINE CAL_STS(N_ELE,N_NODE,N_DOF,PE,PR,IJK_ELE,X,Y,RESULT_N, &STS_ELE,STS_ND)
DIMENSION
IJK_ELE(500,3),X(N_NODE),Y(N_NODE),DD(3,3),BB(3,6), &SS(3,6),RESULT_N(N_DOF),DISP_E(6) DIMENSION STS_ELE(500,3),STS_ND(500,3)
WRITE(8,10)
10 FORMAT(//8X,'--STRESSES OF ELEMENT--') CALL CAL_DD(PE,PR,DD) DO 50 IE=1,N_ELE
CALL CAL_BB(IE,N_NODE,N_ELE,IJK_ELE,X,Y,AE,BB) DO 20 I=1,3 DO 20 J=1,6 SS(I,J)=0.0
DO 20 K=1,3
20 SS(I,J)=SS(I,J)+DD(I,K)*BB(K,J)
2020年4月19日
22
文档仅供参考
DO 30 I=1,3
DO 30 J=1,2 IH=2*(I-1)+J
IW=2*(IJK_ELE(IE,I)-1)+J
30 DISP_E(IH)=RESULT_N(IW) STX=0
STY=0 TXY=0 DO 40 J=1,6
STX=STX+SS(1,J)*DISP_E(J) STY=STY+SS(2,J)*DISP_E(J) 40 TXY=TXY+SS(3,J)*DISP_E(J) STS_ELE(IE,1)=STX STS_ELE(IE,2)=STY
STS_ELE(IE,3)=TXY
50 WRITE(8,60)IE,STX,STY,TXY
60 FORMAT(1X,'ELEMENT NO.=',I5/18X,'STX=',E12.6,5X,'STY=', &E12.6,2X,'TXY=',E12.6)
c the following part is to calculate stress components of node
WRITE(8,55)
55 FORMAT(//8X,'--STRESSES OF NODE--') DO 90 I=1,N_NODE
23
2020年4月19日
文档仅供参考
A=0. B=0. C=0. II=0
DO 70 K=1,N_ELE DO 70 J=1,3
IF(IJK_ELE(K,J).EQ.I) THEN
II=II+1
A=A+STS_ELE(K,1) B=B+STS_ELE(K,2) C=C+STS_ELE(K,3)
END IF
70 CONTINUE STS_ND(I,1)=A/II STS_ND(I,2)=B/II STS_ND(I,3)=C/II
WRITE(8,75)I,STS_ND(I,1),STS_ND(I,2),STS_ND(I,3)
75 FORMAT(1X,'NODE NO.=',I5/18X,'STX=',E12.6,5X,'STY=', &E12.6,2X,'TXY=',E12.6) 90 CONTINUE RETURN
END
2020年4月19日
24
