资源描述：
《有限元程序求解弹性力学平面问题》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、计算力学课程设计计算力学课程设计报告有限元程序求解弹性力学平面问题专业：班级：姓名：学号：指导教师：计算力学课程设计有限元程序求解弹性力学平面问题设计目的:1、学习有限元程序求解弹性力学平面问题的方法；2、学习有限元程序编写技巧；3、加深对有限元方法的理解；4、锻炼处理复杂弹性力学问题的能力。题一：例3.9设深梁承受均布荷载，如下图（a）所示。假定E=1，泊松比，不计容重，厚度t=1m，为平面应力问题。因对称去半边结构进行计算，结构支承、单元划分、节点编号如图（b）所示。试画出y=0及y=6m截面的竖向位移图，x=3m截面的应
2、力分布图。(b)(a)1、有限元Fortran源程序如下：COMMON/X1/NJ,NE,NZ,NDD,NPJ,IND,NJ2,EO,UN,GAMA,TE,AECOMMON/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)CALLDATAIF(IND.EQ.0)GOTO10EO=EO/(1.0-UN*UN)UN=UN/(1.0-UN)10CALLTOTSTICALLLOADCAL
3、LSUPPORCALLSOLVEQCALLSTRESSSTOPENDSUBROUTINEDATACOMMON/X1/NJ,NE,NZ,NDD,NPJ,IND,NJ2,EO,UN,GAMA,TE,AECOMMON/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,INDNJ2=N
4、J*2NPJ1=NPJ+1PRINT*,"INPUT:EO,UN,GAMA,TE"READ(5,*)EO,UN,GAMA,TEPRINT*,"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(10
5、0,FILE='1.TXT')WRITE(100,10)(I,(CJZ(I,J),J=1,2),I=1,NJ)10FORMAT(4X,2HNO,6X,1HX,6X,1HY/(I6,2X,F7.2,F7.2))RETURNENDSUBROUTINEELEST(MEO,IASK)COMMON/X1/NJ,NE,NZ,NDD,NPJ,IND,NJ2,EO,UN,GAMA,TE,AECOMMON/X2/JM(100,3),NZC(50),CJZ(100,2),PJ(100,2),B(3,6),D(3,3),*S(3,6),TKZ(20
6、0,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.0IF(IASK.LE.1)GOTO50DO10I=1,3DO10J=1,6B(I,J)=0.010CONTINUEB(1,1)=-BJ-BMB(1,3)=BJB(1,5)=BMB(2,2
7、)=-CJ-CMB(2,4)=CJB(2,6)=CMB(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)B(3,6)=B(1,5)DO20I=1,3DO20J=1,6B(I,J)=B(I,J)/(2.0*AE)20CONTINUED(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.0D(2,3)=0.0D(3,1)=0.0D(3,2)=0.0D
8、(3,3)=EO/(2.0*(1.0+UN))DO30I=1,3DO30J=1,6S(I,J)=0.0DO30K=1,3S(I,J)=S(I,J)+D(I,K)*B(K,J)30CONTINUEIF(IASK.LE.2)GOTO50DO40I=1,6DO40J=1,6计算力学
