PFEM程序输入说明.docVIP

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) 若NF≠0，则输入下面数据 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) 若NP≠0，则输入下面数据 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) 若NPC≠0，则输入下面数据 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) 若NM≠0，则输入下面数据 若NM≥NE，则表示计所有单元的自重，不需输入计自重的单元号； 若NMNE，则需要输入计自重的单元号； 例题验证 1 (选自徐芝伦编《弹性力学简明教程》(第二版))设有对角受压的正方形薄板，荷载沿厚度均匀分布，为2N/m，由于xz和yz面均为该薄板的对称面，所以只需取1/4部分作为计算对象，图1(a,b)，将该对象划分为4个单元，共有6个结点。对称面上的结点没有垂直于对称面的位移分量，因此，在1、2、4三个结点上设置了水平链杆支座，在4、5、6三个结点设置了铅直链杆支座，如图1c所示。