PFEM程序输入说明.doc

PFEM,EXE平面问题有限元程序的数据文件输入格式：基本输入数据文件名：*.DAT计算结果输出数据文件名：*.OUT中间参数输出数据文件名：*.OU1下面按在文件中的前后顺序进行说明：1 基本控制参数信息：NG,NE,MC,NX,NB,EO,VO,DENSITY ,T(共计5个整形数，4个实型数) NG：结构的结点总数； NE：结构的单元总数； MC：平面问题的类型，MC=0，为平面应力，MC=1，为平面应变； NX：荷载工况数； NB：支承位移数； EO：材料弹性模量(Pa)； VO：材料泊松比； DENSITY ：容重(N/m3) T ：材料厚度(m)；2 打印输出控制参数：NWA,NEW,NWK,NWP(4个整形数) 等于1时，输出，否则不输出。3 单元结点信息：(K,(IJK(I,K),I=1,3),K=1,NE) (每行4个整形数，共计NE行) K：单元号； IJK(1,K)：K单元I结点编号； IJK(2,K)：K单元J结点编号； IJK(3,K)：K单元K结点编号；4 结点坐标信息：(K,XY(1,K),XY(2,K),K=1,NG)(每行3个整形数，共计NG行) K：结点号 XY(1,K)：K结点X坐标； XY(2,K)：K结点Y坐标；5 支承信息：(K,MB(1,K),MB(2,K),ZB(K),K=1,NB)(每行3个整形数，1个实型数，共计NB行) K：支承位移序号； MB(1,K)：第K个支承位移所在的结点号； MB(2,K)：第K个支承位移的坐标方向； ZB(K)： 第K个支承位移的数值；6 按NX荷载工况数输入荷载信息：每一荷载工况如下 (1) NF,NP,NPX,NPY,NM(3个整型数) NF：集中荷载个数； NP：均布侧压分布荷载个数； NPC：坐标轴方向侧压分布荷载个数； NM：计自重单元数； (2) 若NF0，则输入下面数据 K,MF(1,K),MF(2,K),ZF(K)(每行3个整形数，1个实型数，共计NF行) K：集中荷载序号； MF(1,K)：第K个集中荷载作用的结点号； MF(2,K)：第K个集中荷载的坐标方向； ZF(K)： 第K个集中荷载的数值； (3) 若NP0，则输入下面数据 K,MP(1,K),MP(2,K),ZP(K)(每行3个整形数，1个实型数，共计NP行) K：分布荷载序号； MP(1,K)：第K个均布侧压分布荷载作用单元边的首结点号； MP(2,K)：第K个均布侧压分布荷载作用单元边的末结点号； ZP(K)： 第K个均布侧压分布荷载的数值；注意：从首结点到末结点必须按照单元结点编号的逆时针方向。 (4) 若NPC0，则输入下面数据 K,MPC(1,K),MPC(2,K) ,MPC(3,K,ZPC1(K), ZPC2(K) (每行4个整形数，2个实型数，共计NPC行) K：分布荷载序号； MPC(1,K)：第K个坐标轴方向侧压分布荷载作用单元边的首结点号； MPC(2,K)：第K个坐标轴方向侧压分布荷载作用单元边的末结点号； MPC(3,K)：第K个坐标轴方向侧压分布荷载作用的坐标轴方向，X方向其值为1，Y方向其值为0； ZPC1(K)： 第K个坐标轴方向侧压分布荷载首结点处的数值； ZPC2(K)： 第K个坐标轴方向侧压分布荷载末结点处的数值；注意：从首结点到末结点必须按照单元结点编号的逆时针方向。 (5) 若NM0，则输入下面数据 若NMNE，则表示计所有单元的自重，不需输入计自重的单元号； 若NMNE，则需要输入计自重的单元号；例题验证1 (选自徐芝伦编弹性力学简明教程(第二版)设有对角受压的正方形薄板，荷载沿厚度均匀分布，为2N/m，由于xz和yz面均为该薄板的对称面，所以只需取1/4部分作为计算对象，图1(a,b)，将该对象划分为4