FEAP目录FEAP研制者：RobertL.Taylor教授美国加州大学伯克利分校土木与环境工程系研制目的：教学与科研配套教材：O.C.Zienkiewica,R.L.Taylor,TheFiniteElementMethod,McGraw-Hill,London,4theditionFEAP界面FEAP求解过程建立有限元模型单元类型节点载荷及边界条件材料模型其它求解求解算法FEAP-运行FEAP程序将提示用户键入输入数据文件名（包含扩展名）和输出数据文件名（存放计算结果）。FEAP–输入数据文件各行最大长度为255个字符每行最多可包括16个数据项，每个数据项最长为15个字符（14个字符，1个分隔符）数据分隔符：空格或逗号自由格式常数3.56,-12.37,1.34e+5,-4.36d-05参数（参数名最多为两个字符，且以字符开头，如a,aa,a2,er等）pi=3.1415926,E=200.0e9,mu=0.3FEAP–输入数据文件表达式（由常数、参数、函数和运算符组成）FEAP–输入数据文件函数abs,exp,int,log,sqrt,sin,cos,tan,atan,asin,acos,sind,cosd,tand,atand,asind,acosd,cosh,sinh,tanh以d结尾的三角函数和反三角函数的自变量的单位为度，其它的为弧度。函数的自变量可以是参数，但不能含有括号。如：p=4.*atan(1),q=1./(3.+a),q=tan(q)建立有限元网格FEAP**King-posttrussanalysis！总信息331222！（节点总数，单元总数，材料组数，维数，节点最大自由度，单元最大节点数）MATErial1！定义材料TRUSSELASticisotropic2.0e11CROSssection21.98e-4COORdinates！定义节点100.00.0！节点号，增量，x，y，z202.00.0301.01.0ELEMents！定义节点10112！单元号，增量号，材料号，节点1，节点22012330131BOUNdaryrestraints！定义边界1011！节点，增量，自由度1，自由度22001FORCe！定义荷载300.-10e3！节点，增量，x向力，y向力ENDINTESTOPFEAP–网格控制信息FEAPtitleNUMNP,NUMEL,NUMMAT,NDM,NDF,NENNUMNP–节点总数（可以0代替，程序自动计算）NUMEL–单元总数（可以0代替，程序自动计算）NUMMAT–材料性质组数NDM–空间维数(1、2或3)NDF–每个节点的最大自由度数目NEN–每个单元的最大节点数目END(网格数据输入结束)PRINt/NOPRint(输出/不输出网格数据)FEAP–网格COORdinates命令（输入节点坐标）COORnode1,ngen1,(x(i),i=1,ndm)node2,ngen2,(x(i),i=1,ndm)…以空行结束node1–节点编号ngen1–节点自动生成增量。生成的节点序列为node1,node1+ngen1,node1+2*ngen1,…,node2ngen1等于0时不自动生成节点。x(i)–节点坐标FEAP–网格例子–圆盘FEAP–网格FEAP**Example1.CircularDisk:Basicinputs19121224FEAP–网格ELEMent命令（输入单元数据）ELEMnelm1,ngen1,matl1,(ix(i),i=1,nen)nelm2,ngen2,matl2,(ix(i),i=1,nen)…<以空行结束>nelm1–单元号ngen1–单元自动生成节点号增量。matl1–材料组号ix(i)–单元节点号单元数据必须按单元编号次序排列。FEAP–网格FEAP–边界条件BOUNdary命令（指定边界约束条件）BOUNnode1,ngen1,(id(i),i=1,ndf)node2,ngen2,(id(i),i=1,ndf)…<以空行结束>id(i)=0：在该自由度上位移不约束(缺省)id(i)≠0：在该自由度上位移被约