FEAPpv Linux编译和运行方法.pdf

FEAPpv Linux 编译及运行方法 整理日期：2016-03-06 问题反馈：师访，QQ ：1549221758，pomato157300@126.com 1. FEAPpv 简介 FEAP 是由加州大学伯克利分校 （University of California, Berkeley ）开发的通用有限元程序， FEAPpv 是其对应的开放 源代码的个人版本。程序最初是由 Robert L. Taylor 教授 （/~rlt/ ）为满足教学和科研目的而编写的。与 FEAPpv 对应的教材为 The Finite Element Method, 6th edition。 FEAPpv 的功能包括： （1）有限的单元库 （与付费版本相比）和用户自定义单元。 ①结构分析： 弹簧-质点单元（pointnd.f ）、一维-二维-三维桁架单元（trussnd.f ）、二维梁单元 （frame2d.f）、三维梁单元（frame3d.f）、二维壳单元（shell2d.f）、三维壳单元（shell3d.f）、二维（平 面应力-平面应变-轴对称）结构单元（solid2d.f ）、三维实体单元（solid3d.f ）、三维四边形膜单元 （membr3d ）。 ②热传导分析：二维热传导表面单元（convec2d.f）、三维热传导表面单元 （convec3d.f）、二维热 传导单元（therm2d.f ）、三维热传导单元（therm3d.f ）。 ③用户自定义单元接口：可通过ver31\user 文件夹下的elmt01.f 至elmt05.f 实现。 （2 ）少量的材料本构（与付费版本相比）和用户自定本构。 ①弹性、横观各向同性、正交各向异性、超弹性（Neo-Hookean ）、粘弹性、塑性 （包括可用于金 属材料模拟的Mises 和岩土材料模拟的Drucker-Prager 屈服准则和相关流动法则）、和热传导本构。这 些材料相关代码都在material.f 文件中。 ②用户自定义材料本构关系接口：可通过ver31\user 文件夹下的umati1.f 至umati5.f 和umatl1.f 至umatl5.f 实现。 （3）云图输出功能。 相关程序在ver31\plot 文件夹下。 FEAPpv 的用法：用户编写输入文件，在终端中通过feappv 命令激活程序，并读入输入文件，程 序运行完毕后保存输出文件、并根据输入文件中的plot 语句绘制相关图形。以下内容为官网上的输入 1 / 5 文件示例 （Iblock ，该示例计算一个二维平面应变块体受拉变形问题）： feap * * test on solution !项目名称定义 0 0 0 2 2 4 !问题描述，6 个数字的含义分别是：节点数目、单元数目、材料种类、维度、每个节点的 !未知量数目、每个单元的最大节点数目。由于程序可自动识别，所以1-3 个数字可设为0 block !建模及单元划分 （规则网格可直接用block 命令形成） cart 4 4 !笛卡尔坐标系下4x4=16 单元模型 1 0 0 !模型几何坐标定义 2 5 0 3 5 5 4 0 5 ebou !边界条件定义 1 0 1 0 !给x=0 的节点施加x 方向的位移约束 2 0 0 1 !给y=0 的节点施加y 方向的位移约束 forc !节点载荷(施加到24,25 号节点上) 24 0 0 10 25 0 0 10 mate 1 !材料参数定义 solid