有限元编程算例.docxVIP

有限元编程算例（Fortran） 本程序通过 Fortran 语言编写，程序在 In tel Parallel Studio XE 2013 with VS2013中成功运行，程序为《计算力学》（龙述尧等编）一书中的源程序，仅作 研究学习使用，省去了敲写的麻烦。 3.7.4算例 例工夕设谦聲祇受沟布戴荷、如图3-3叙叮所示「假進￡=」，泊松比世二①17?不计 容重+厚度才=1 为平面应力问題，因对称取半边结构计真?結枸支承，跟兀划分，节点 绸号如图3,16（b）JPfziip试矗出及y = S m截啲的竖向位移图応=总mflg面的靳 应力分布图， ⑴ （山 图玉開燮均布栽荷的简支漂梁 源程序为： !Page149 C0MM0N/X1/NJ,NE,NZ,NDD,NPJ,IND,NJ2,E0,UN,GAMA,TE,AE COMMON/X2/JM(100,3),NZC(50),CJZ(100,2),PJ(100,2),B(3,6),D(3,3),S(3,6),TKZ(200,20),EKE(6,6),P(200) 0PEN(5,FILE=DATAIN) !OPEN(6,FILE=DATAOUT,STATUS=NEW) CALL DATA IF(IND.EQ.0)GOTO 10 EO=EO/(1.0-UN*UN) UN=UN/(1.0-UN) 10 CALL TOTSTI CALL LOAD CALL SUPPOR CALL SOLVEQ CALL STRESS PAUSE !STOP END SUBROUTINE DATA COMMON/X1/NJ,NE,NZ,NDD,NPJ,IND,NJ2,EO,UN,GAMA,TE,AE COMMON/X2/JM(100,3),NZC(50),CJZ(100,2),PJ(100,2),B(3,6),D(3,3),S(3,6),TKZ(200,20),EKE(6,6),P(200) READ(5,*)NJ,NE,NZ,NDD,NPJ,IND NJ2=NJ*2 NPJ1=NPJ+1 READ(5,*)EO,UN,GAMA,TE READ(5,*)((JM(I,J),J=1,3),I=1,NE) READ(5,*)((CJZ(I,J),J = 1,2),I=1,NJ) !Page150 READ(5,*)(NZC(I),I=1,NZ) READ(5,*)((PJ(I,J),J=1,2),I=1,NPJ1) WRITE(6,10)(I,(CJZ(I,J),J=1,2),I=1,NJ) 10 FORMA T(4X,2HNO,6X,1HX,6X,1HY/(I6,2X,F7.2,F7.2)) RETURN END SUBROUTINE ELEST(MEO,IASK) COMMON/X1/NJ,NE,NZ,NDD,NPJ,IND,NJ2,EO,UN,GAMA,TE,AE COMMON/X2/JM(100,3),NZC(50),CJZ(100,2),PJ(100,2),B(3,6),D(3,3),S(3,6),TKZ(200,20),EKE(6,6),P(200) IE=JM(MEO,1) JE=JM(MEO,2) ME=JM(MEO,3) CM=CJZ(JE,1)-CJZ(IE,1) BM=CJZ(IE,2)-CJZ(JE,2) CJ=CJZ(IE,1)-CJZ(ME,1) BJ=CJZ(ME,2)-CJZ(IE,2) AE=(BJ*CM-BM*CJ)/2.0 IF(IASK.LE.1) GOTO 50 DO 10 I=1,3 DO 10 J=1,6 B(I,J)=0.0 10 CONTINUE B(1,1)=-BJ-BM B(1,3)=BJ B(1,5)=BM B(2,2)=-CJ-CM B(2,4)=CJ B(2,6)=CM B(3,1)=B(2,2) B(3,2)=B(1,1) B(3,3)=B(2,4) B(3,4)=B(1,3) B(3,5)=B(2,6) !Page151 B(3,6)=B(1,5) DO 20 I=1,3 DO 20 J=1,6 B(I,J)=B(I,J)/(2.0*AE) 20 CONTINUE D(1,1)=EO/(1.0-UN*UN) D(1,2)=EO*UN/(1.0-UN*UN) D(2,1)=D(1,2) D(2,2)=D(1,1) D(1,3)=0.0 D(2,3)=0.0 D(3,1)=0.0 D(3,2)=0.0 D(3,3)=EO/(2.0*(1.0+UN)) DO 30 I=1,3 DO 30 J=1,6 S(I,J)=0.0 DO 30 K=1,3 S(I,J)=S(I,J)+D(I,K)*B(K,J) 30 CONTINUE IF(IASK.LE.2) G