用伪代码方式写有限元求解步骤.docVIP

作业：一 1 用伪代码方式写出有限元求解步骤 解：有限元分析流程图： Ⅰ.平面问题程序使用说明 输入数据格式依次如下： 1.NG NG:结构的结点总数 NE,MC,NX,NB,ND,EO,VO,T NE:结构单元总数 MC:计算控制类型参数=0平面应力 =1平面应变 NX:作用荷载组数 NB：给定位移的个数 NB:给定位移的个数 ND：结构刚度阵的半带宽 EO:弹性模量 VO：泊松比 T：单元的厚度 NWA,NWE,NWK,NWP,NWD NWA：单元参数的输出控制参数 NWE:单元刚度矩阵的输出控制参数 NWK：结构刚度矩阵的输出控制参数 NWP：荷载向量的输出控制参数 NWD：结点位移的输出控制参数 输出控制参数=1输出 =0不输出 IJM(3，NE):单元结点编码数组 IJM(1，I),IJM(2，I),IJM(3，I)，第I个三角形单元的节点号，按单元编号顺序填写。 XY(2，NG)：结点坐标数组 XY(1，I):第I个结点的坐标，XY(2，I):第I个结点的Y坐标 按结点编号顺序填写。 MB(2，NB),ZB(NB):给定位移约束的信息数组与值数组 MB(1，I):第I个给定位移所在的结点号 NB(2，I)=1：给定X方向位移 =2：给定Y方向位移 ZB(NB):给定位移值（以坐标正向为正) 7. NF,NP NF:作用于结点上的集中荷载的个数 NP:作用于均布侧压的单元边数 若NF0，填8 MF(2，NF),ZF(NF):作用于结点上集中荷载的信息组与值数组 MF(1，I):第I个集中荷载作用的结点号 MF(2，I)=1:作用于x方向的集中力 =2:作用y方向的集中力 ZF(NF):作用的集中力值 若NP0，填9 MP(2，NP),ZP(NP):作用于单元边上均布荷载的信息数组与值数组 MF(1,I):第I个均布荷载作用边的起始结点号 MF(2,I):第I个均布荷载作用边的终止结点号，逆时针排列 ZP(NP）：第I个均布荷载值 输入数据格式，建立数据文件，文件名小于12个字符： NG NE,MC,NX,NB,ND,E,P,T NWA,NWE,NWK,NWP,NWD IJM(3,NE) XY(2,NG) MB(2,NB),ZB(NB) NF,NP 当NF0,则填8 MF(2,NF),ZF(NF) 当NP0,则填9 MP(2,,NP），ZP(NP) 若NX1,即多组荷载情况，重复7～９,若计算多个结构则重复１～９。 结束 平面问题主程序框图 输入子程序INPUT伪代码 DIMENSIONG IJM(3,NE),XY(2,NG),MB(2,NB),ZB(NB) READ(5,*)((IJM(I,L),I=1,3),L=1,NE) READ(5,*)((XY(I,J),I=1,2),J=1,NG) READ(5,*)(MB(I,L),I=1,2),L=1,NB),(ZB(L),L=1,NB) WRITE(6，10) 10 FORMAT(/10X,‘ELEMENT CODE BLOCKS IJM(2，NE)’,/) WRITE(6，20)((IJM(M,I)M=1，3),I=1，NE) 20 FORMAT(1X,4(3I4，3X),3I4) WRITE(6，30) 30 FORMATE(/,10X,‘COORDINATES OF NODES XY(2,NG)’,/) WRITE(6,40)((XY(M,I),M=1,2),I=1,NG) 40 FORMATE(1X,6E12.5) WRITE(6,50) 50 FORMATE(/10X,INRMATION AND VALUES OF GIVEN DISPLACEMENTS,/,20X,‘MB’,15X,‘ZB’,/) DO 70 I=1,NB WRITE(6，60)(MB(K,I),K=1,2),ZB(I) 60 FORMATE(1OX,2I5，5X,E14.6) 70 CONGTINUE RETURN END 2,计算单元参数子程序伪代码 X(2，5)当前计算单元的结点坐标 B(7) 当前计算单元的单元参数 BCA(7，NE) 将B（7）按单元存储于该数组中 DIMENSION IJM(3，NE),XY(2,NG),BCA(7,NE),X(2,5),B(7) IF(NEWA.EQ.1) WRITE(6，5) 5 FORMAT(/10X,’PARAMETERS OF ELEMENTS BCA(7，NE)’/) DO 80 I=1,NE DO 10 K=1,3