《有限单元法基本原理和数值方法》一书的源程序讲述.docVIP

《有限单元法基本原理和数值方法》一书的源程序 -------------------------------------------------------------------------------- !********************************************************************* ! * ! PL ---- PROGRAM OF PLANE PROBLEM 96.1 * ! * !********************************************************************* C C-------- 输入数据顺序-------- C 1.NG 1整型 C NG 结构结点总数 C NG=0 则停止运行 C 2.NE,MC,NX,NB,ND,EO,VO,T 5整型3实型 C NE 结构单元总数 C MC 计算控制类型参数 C MC=0 平面应力 C =1 平面应变 C NX 作用载荷组数 C NB 给定位移个数 C ND 结构刚度矩阵的半带宽 C EO 弹性模量 C VO 泊松比 C T 单元(结构)的厚度 C 3.NWA NWE,NWK,NWP,NWD 5整型 C 输出控制参数 C =1 输出 C =0 不输出 C NWA 单元参数的输出控制参数 C NWE 单元刚度矩阵的输出控制参数 C NWK 结构刚度矩阵的输出控制参数 C NWP 载荷向量的输出控制参数 C NWD 结点位移的输出控制参数 C 4.IJM(3,NE) 单元结点编码数组 3×NE整型 C IJM(1,I); IJM(2,I); IJM(3,I) C 第I个三角形单元的结点编号,按结点编号顺序填写 C 5.XY(2,NG) 结构结点坐标数组 2×NG实型 C 6.MB(2,N,ZB(N 2×NB整型,NB实型 C MB(1,I)---第I个给定位移所在的结点号 C MB(2,I)=1--给定X方向位移 C =0--给定Y方向位移 C ZB(I)----给定位移值(以坐标正向为正) C 7.NF,NP 2整型 C NF-----作用在结点上的集中载荷(坐标方向)的个数 C NP-----作用均布侧压的单元边数 C 若 NF0 则填写 C 8.MF(2,NF),ZF(NF) 2×NF整型,NF实型 C MF(1,I)---第I个集中载荷所在的结点号 C MF(2,I)=1--给定X方向集中力 C =0--给定Y方向集中力 C ZF(I)-----作用的集中力值 C 若 NP0 则填写 C 9.MP(2,NP),ZP (NP) 2×NP整型,NP实型 C MP(1,I)----第I个载荷作用边的起始结点号 C MP(2,I)----第I个载荷作用边的起始结点号 C ZP(I)------第I个均布载荷值 C C 若 NX1 重复 7.-9. (NX-1) 次 C C 最后 NG=0 表示数据结束 C !-------- 输出数据顺序-------- C 1.IJM(3,NE) 单元结点编码数组 3×NE整型 C IJM(1,I); IJM(2,I); IJM(3,I) C 第I个三角形单元的结点编号,按结点编号顺序填写 C 2.XY(2,NG) 结构结点坐标数组 2×NG实型 C C 若NWA=1,则输出 C 3.I,B(7) 单元参数 NE行,1×NE整型,7×NE实型 C 每行结构为:NE=+单元号+Bi+Bj+Bm+Ci+Cj+Cm+A C C 若NWE=1,则输出 C 4.IO,EK(6×6)单元刚度阵 NE行,1×NE整型,6×6×NE实型 C 每行结构为:NE=+单元号+EK(单元刚度阵) C C 若NWK=1,则输出 C 5.SK(NT,ND) 结构刚度矩阵 NT×ND=2NG×ND实型 C C 若NWD=1,则输出 C 6.I,B 结点位移数据 NG行,1×NG整型,2×NG实型 C 每行结构为:单元号+U+V C C 7.S1,S2,S3,X1,X2,CTA C 单元应力数据 6×NE实型 C 分别代表σx,σy,τxy,σ1,σ2和主应力方向 C C 若 NX1 重复 6.-7. (NX-1) 次 C !--------可调数组分配-------- C C 实型数组 C(100000) 整型数组 IA(100000) C C(