资源描述：
《有限单元法程序报告.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、目录程序一：平面刚架静力分析程序（PF.FOR）1程序二：平面三节点有限元程序17程序三：四节点矩形薄板单元程序24程序一：平面刚架静力分析程序（PF.FOR）已知各杆截面均为矩形，柱截面宽0.4m,高0.5m,梁截面宽0.4m,高0.4m，各杆E=3.65×104MPa。图2节点、单元编号如下图3，1.2.3…..为节点号，①②③……为单元号：图3总共有13个节点，13个单元。计算源程序如下：!PF.FOR(Aprogramforanalysisofplaneframe)!Version6.32004!Mainprogramreadsthecontrolda
2、ta&organizesthewhole!calculationbycallingsubroutines.DIMENSIONW(80000)CHARACTERIDFN*20,TITLE(5)*72READ(*,'(A12)')IDFNOPEN(3,FILE=IDFN,STATUS='OLD')READ(3,'(A72)')(TITLE(M),M=1,5)READ(3,*)E,NM,NJ,NS,NLCL1=1L2=L1+NML3=L2+NML4=L3+NML11=L4+NML12=L11+NJL21=L12+NJL22=L21+NSL31=L22+NSL41=L
3、31+6*NMCALLIOMJS(TITLE,E,NM,NJ,NS,NLC,W(L1),W(L2),W(L3),&W(L4),W(L11),W(L12),W(L21),W(L22))CALLLCVCT(NM,W(L1),W(L2),W(L31),NJ,N)CALLLCDIA(NM,N,W(L31),W(L41),W(L41),W(L41),MAXBDW,NA)L51=L41+NL52=L51+36L53=L52+NA*2L54=L53L61=L54+N*2NW=L61+6*NM-1WRITE(*,1)NA,NW1FORMAT(/40X,'(NA=',I6,')
4、',/40X,'(NW=',I6,')')CALLFORMA(E,NM,NJ,N,NA,W(L1),W(L2),W(L3),W(L4),&W(L11),W(L12),W(L31),W(L51),W(L41),W(L52))CALLAS(NS,N,NA,W(L21),W(L41),W(L52))CALLLDLT(N,NA,W(L41),W(L52),W(L53))DO100LC=1,NLCREAD(3,*)NLJL62=L61+NLJL63=L62+NLJL64=L63+NLJL71=L61L81=L71+6*NMCALLB0(LC,N,NLJ,W(L54))I
5、F(NLJ.NE.0)CALLIOLJB(N,NLJ,W(L61),W(L62),&W(L63),W(L64),W(L54))READ(3,*)NLML82=L81+NLML83=L82+NLML84=L83+NLMCALLF0(NLM,NM,W(L71))IF(NLM.NE.0)CALLIOLMFB(NM,NJ,N,NLM,W(L81),&W(L82),W(L83),W(L84),W(L1),W(L2),W(L11),&W(L12),W(L31),W(L71),W(L54))CALLBS(NS,N,W(L21),W(L22),W(L54))CALLSLVEQ
6、(N,NA,MAXBDW,W(L41),W(L52),W(L54))CALLOJD(NJ,N,W(L54))CALLCOTF(E,NM,NJ,N,W(L1),W(L2),W(L3),W(L4),&W(L11),W(L12),W(L31),W(L54),W(L71))NW=L84+NLM-1WRITE(*,1)NA,NW100CONTINUEWRITE(*,'(/)')ENDSUBROUTINEIOMJS(TITLE,E,NM,NJ,NS,NLC,IST,IEN,&AR,RI,X,Y,IS,VS)!Readdataofmembers,joints,support
7、s&printthemDIMENSIONIST(NM),IEN(NM),AR(NM),RI(NM),&X(NJ),Y(NJ),IS(NS),VS(NS)CHARACTERTITLE(5)*72WRITE(*,'(/)')WRITE(*,1)(TITLE(M),M=1,5)1FORMAT(1X,A72)WRITE(*,2)E,NM,NJ,NS,NLC2FORMAT(/13X,'TheInputData'&//10X,'TheGeneralInformation'&//6X,'E',9X,'NM',5X,'NJ',5X,'NS',5X,'NLC'&/1X,1PE1
8、0.3,4I7)READ(3,*)(I
