有限元程序.docVIP

学院：工程学院 姓名：XX 学号：000000 算例一 ：平面刚桁架静力分析程序（PF.FOR） 如图1所示刚架材料为钢筋混凝土，杆件截面为矩形，柱截面宽0.4m，高0.4m；梁截面宽0.4m，高0.7m。 图1 解：建立坐标系，确定节点，划分单元。单元和节点编号如上图所示。 计算源程序如下： ! 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(TITLE,E,NM,NJ,NS,NLC,IST,IEN, AR,RI,X,Y,IS,VS) ! Read data of members, joints, supports print them DIMENSION IST(NM),I