用常应变三角形单元解弹性力学平面问题的程序

用常应变三角形单元解弹性力学平面问题

的程序

******************************************************************* * ANALYSIS PROGTAM OF FINITE ELEMENT METHOD * * FOR PLANE STRESS/STRAIN OF TRIANGULAR ELEMENT * * ----- FEMT3.FOR ----- * *------------------------------------------------------------- * * Subroutines: 1-SDATA, 2-STE, 3-ATE, 4-DTE, 5-BTE, 6-STIFF * * 7-EQUPE, 8-INSCD, 9-BGSMT, 10-SIGME * ******************************************************************* DIMENSION LND(50,3),X(100),Y(100),JR(20,3),PJ(20,3),P(200) REAL KS(200,100)

OPEN(5,FILE='FEMT3.DAT')

OPEN(6,FILE='FEMT3.OUT',STATUS='NEW')

READ(5,*) NJ,NE,NS,NPJ,IPS(结点、单元、支承、荷载、类型) WRITE(6,*)' FINITE ELEMENT ANALYSIS IN PLANE PROBLEM' WRITE(6,*)' SOURCE DATA OUTPUT'

WRITE(6,20) NJ,NE,NS,NPJ,IPS

20 FORMAT(4X,'NJ',3X,'NE',3X,'NS',3X,'NPJ',2X,'IPS'/1X,5I5) IF(IPS.EQ.0) WRITE(6,*)' PLANE STRESS PROBLEM' IF(IPS.EQ.1) WRITE(6,*)' PLANE STRAIN PROBLEM'

CALL SDATA(NJ,NE,NS,NW,NPJ,IPS,E,PR,T,V,LND,X,Y,JR,PJ) NJ2=2*NJ

WRITE(6,50) NJ2

50 FORMAT(/1X,'DEGREES OF FREEDOM=',I5)

WRITE(6,60) NW

60 FORMAT(1X,'BAND WIDTH=',I5)

CALL STIFF(NJ,NE,NJ2,NW,LND,X,Y,E,PR,T,KS)(总刚6) CALL EQUPE(NJ,NE,NPJ,NJ2,T,V,LND,X,Y,PJ,P)({P}7) CALL INSCD(NS,NW,NJ2,JR,KS,P)(引入支承条件8) CALL BGSMT(NJ,NJ2,NW,KS,P)（解方程9）

CALL SIGME(NE,NJ,NJ2,E,PR,LND,X,Y,P)（求应力10） CLOSE(5)

CLOSE(6)

END

*--------------------------------------------------------

C SUBPROGRAM-1

C INPUT STRUCTURAL DATA

SUBROUTINE SDATA(NJ,NE,NS,NW,NPJ,IPS,E,PR,

* T,V,LND,X,Y,JR,PJ)

DIMENSION LND(NE,3),X(NJ),Y(NJ),JR(NS,3),PJ(NPJ,3) READ(5,*) E,PR,T,V(弹性模量、泊松比、厚度、容重)