有限单元法程序报告.doc

目录 程序一：平面刚架静力分析程序（PF.FOR） 1 程序二：平面三节点有限元程序 17 程序三：四节点矩形薄板单元程序 24  程序一：平面刚架静力分析程序（PF.FOR） 已知各杆截面均为矩形，柱截面宽0.4m,高0.5m, 梁截面宽0.4m,高0.4m，各杆E=3.65×104 MPa。 图2 节点、单元编号如下图3，1.2.3…..为节点号，①②③……为单元号： 图3 总共有13个节点，13个单元。 计算源程序如下： ! PF.FOR (A program for analysis of plane frame) ! Version 6.3 2004 ! Main program reads the control data organizes the whole ! calculation by calling subroutines. DIMENSION W(80000) CHARACTER IDFN*20,TITLE(5)*72 READ(*,(A12))IDFN OPEN(3,FILE=IDFN,STATUS=OLD) READ(3,(A72))(TITLE(M),M=1,5) READ(3,*)E,NM,NJ,NS,NLC L1=1 L2=L1+NM L3=L2+NM L4=L3+NM L11=L4+NM L12=L11+NJ L21=L12+NJ L22=L21+NS L31=L22+NS L41=L31+6*NM CALL IOMJS(TITLE,E,NM,NJ,NS,NLC,W(L1),W(L2),W(L3), W(L4),W(L11),W(L12),W(L21),W(L22)) CALL LCVCT (NM,W(L1),W(L2),W(L31),NJ,N) CALL LCDIA (NM,N,W(L31),W(L41),W(L41),W(L41),MAXBDW,NA) L51=L41+N L52=L51+36 L53=L52+NA*2 L54=L53 L61=L54+N*2 NW=L61+6*NM-1 WRITE (*,1)NA,NW 1 FORMAT(/40X,( NA=,I6, ),/40X,( NW=,I6, )) CALL FORMA (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)) CALL AS(NS,N,NA,W(L21),W(L41),W(L52)) CALL LDLT(N,NA,W(L41),W(L52),W(L53)) DO 100 LC=1,NLC READ (3,*)NLJ L62=L61+NLJ L63=L62+NLJ L64=L63+NLJ L71=L61 L81=L71+6*NM CALL B0(LC,N,NLJ,W(L54)) IF(NLJ.NE.0) CALL IOLJB(N,NLJ,W(L61),W(L62), W(L63),W(L64),W(L54)) READ(3,*)NLM L82=L81+NLM L83=L82+NLM L84=L83+NLM CALL F0(NLM,NM,W(L71)) IF(NLM.NE.0) CALL IOLMFB(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)) CALL BS(NS,N,W(L21),W(L22),W(L54)) CALL SLVEQ(N,NA,MAXBDW,W(L41),W(L52),W(L54)) CALL OJD(NJ,N,W(L54)) CALL COTF(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-1 WRITE(*,1)NA,NW 100 CONTINUE WRITE(*,(/)) END SUBROUTINE IOMJS(TITL