有限元程序求解弹性力学平面问题.doc

计算力学 课程设计报告 有限元程序求解弹性力学平面问题 专 业： 班 级： 姓 名： 学 号： 指导教师： 有限元程序求解弹性力学平面问题 设计目的: 学习有限元程序求解弹性力学平面问题的方法； 学习有限元程序编写技巧； 加深对有限元方法的理解； 锻炼处理复杂弹性力学问题的能力。 题一： 例3.9 设深梁承受均布荷载，如下图（a）所示。假定E=1，泊松比，不计容重，厚度t=1m，为平面应力问题。因对称去半边结构进行计算，结构支承、单元划分、节点编号如图（b）所示。试画出y=0及y=6m截面的竖向位移图，x=3m截面的应力分布图。 1、有限元Fortran源程序如下： 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) 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 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) PRINT*,INPUT:NJ,NE,NZ,NDD,NPJ,IND READ(5,*)NJ,NE,NZ,NDD,NPJ,IND NJ2=NJ*2 NPJ1=NPJ+1 PRINT*,INPUT:EO,UN,GAMA,TE READ(5,*)EO,UN,GAMA,TE PRINT*,INPUT:JM READ(5,*)((JM(I,J),J=1,3),I=1,NE) PRINT*,INPUT:CJZ READ(5,*)((CJZ(I,J),J=1,2),I=1,NJ) PRINT*,INPUT:NZC READ(5,*)(NZC(I),I=1,NZ) PRINT*,INPUT:PJ READ(5,*)((PJ(I,J),J=1,2),I=1,NPJ1) OPEN(100,FILE=1.TXT) WRITE(100,10)(I,(CJZ(I,J),J=1,2),I=1,NJ) 10 FORMAT(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