ANSYS建模apdl命令流实例应用

结合自身经验，谈ANSYS中的APDL命令(二)发表时间：2009-5-10 作者: 倪欣来源: e-works关键字: ANSYS APDL 命令流在ANSYS中，命令流是由一条条ANSYS的命令组成的一个命令组合，这些命令按照一定顺序排布，能够完成一定的ANSYS功能，本文是作者结合自身经验所总结的一些命令。

1.1 /prep7(进入前处理)定义几何图形：关键点、线、面、体（1）.csys,kcnkcn , 0 迪卡尔坐标系1 柱坐标2 球4 工作平面5 柱坐标系（以Y轴为轴心）n 已定义的局部坐标系（2）.numstr, label, value 设置以下项目编号的开始nodeelemkplineareavolu注意：vclear, aclear, lclear, kclear 将自动设置节点、单元开始号为最高号，这时如需要自定义起始号，重发numstr（3）.K, npt, x,y,z, 定义关键点Npt：关键点号，如果赋0，则分配给最小号（4）.Kgen,itime,Np1,Np2,Ninc,Dx,Dy,Dz,kinc,noelem,imoveItime：拷贝份数Np1,Np2,Ninc：所选关键点Dx,Dy,Dz：偏移坐标Kinc：每份之间节点号增量noelem: “0” 如果附有节点及单元，则一起拷贝。

“1”不拷贝节点和单元imove：“0” 生成拷贝“1”移动原关键点至新位置，并保持号码，此时（itime,kinc,noelem）被忽略注意：MAT,REAL,TYPE 将一起拷贝，不是当前的MAT,REAL,TYPE（5）.A, P1, P2, ……… P18 由关键点生成面（6）.AL, L1,L2, ……,L10 由线生成面面的法向由L1按右手法则决定，如果L1为负号，则反向。

（线需在某一平面内坐标值固定的面内）（7）.vsba, nv, na, sep0,keep1,keep2 用面分体（8）.vdele, nv1, nv2, ninc, kswp 删除体kswp: 0 只删除体1 删除体及面、关键点（非公用）（9）.vgen, itime, nv1, nv2, ninc, dx, dy, dz, kinc, noelem, imove 移动或拷贝体itime: 份数nv1, nv2, ninc：拷贝对象编号dx, dy, dz ：位移增量kinc: 对应关键点号增量noelem,：0：同时拷贝节点及单元1：不拷贝节点及单元imove：0：拷贝体1：移动体（10）.cm, cname, entity 定义组元，将几何元素分组形成组元cname: 由字母数字组成的组元名entity: 组元的类型（volu, area, line, kp, elem, node）（11）.cmgrp, aname, cname1, ……,cname8 将组元分组形成组元集合aname: 组元集名称cname1……cname8: 已定义的组元或组元集名称1.2 定义几个所关心的节点，以备后处理时调用节点号。

ANSYS计算结果无难事，APDL经典命令让你的模型舞起来1、让你的ANSYS模型'舞'起来ANSYS计算结果的动画可采用ANTIME、ANMODE、ANCNTR、ANHARM等自动生成动画(AVI格式)，使结果展示更加生动直观，相信使用ANSYS的都会制作。

然而，几何模型或有限元模型则无动画显示功能，有时为展示模型本身，会从多个角度截取图片。

那么，模型能否也可制作动画呢？答案是肯定的。

利用ANSYS的图形存储命令/SEG可以实现此功能，让你的模型动起来。

具体过程详见命令流中及其注释，动画上传总是失败，自己生成不要观看吧。

Finish$/clear$/prep7!简单的创建几何模型以减少篇幅blc4,0,0,4,2,5cyl4,2,4,1,,2,,4!关闭图例信息/plopts,info,off!以下开始制作模型动画!删除当前储存的图形/seg,dele/seg,multi,jhdh,1 !独立存储且不覆盖，文件名为jhdh/auto,1 !自动计算与图形区合适显示方式!正视/view,1,0,0,1$vplot!侧视/view,1,1$vplot!俯视/view,1,,1$vplot!D视图/view,1,1,1,1$vplot!循环36次，每次改变10度视角*do,i,1,36$/ang,1,10,ys,1$/replot$*enddo!关闭图形存储操作，保存为jhdh.avi文件/seg,off$/anfile,save,jhdh,avi其实比较简单，一旦进入模型动画制作过程，所有的xPLOT(x=KLAVNE)绘制的图形都将进入动画序列，按显示过程形成一部连续的动画。

2、用一个命令解决ANSYS数据列表分页早年初学ANSYS时，经常用到xLIST(如NLIST、ELIST、KLIST、LLIST、ALIST、VLIST等命令)和PRxSOL（如PRNSOL、PRESOL、PRRSOL、PRETAB、PRPATH）等列表命令，并希望将这些内容保存到TXT文件中，然后再导入EXCEL中处理。

在ANSYS中，命令流是由一条条ANSYS的命令组成的一个命令组合，这些命令按照一定顺序排布，能够完成一定的ANSYS功能，这些功能一般来说通过菜单操作也能够实现（而那些命令流能够实现，菜单操作实现不了的单个命令比较少见）。

以下命令是结合我自身经验，和前辈们的一些经验而总结出来的，希望对大家有帮助。

（1）．Lsel, type, item, comp, vmin, vmax, vinc, kswp 选择线type: s 从全部线中选一组线r 从当前选中线中选一组线a 再选一部线附加给当前选中组aunoneu(unselect)inve: 反向选择item: line 线号loc 坐标length 线长comp: x,y,zkswp: 0 只选线1 选择线及相关关键点、节点和单元（2）．Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点type: S: 选择一组新节点（缺省）R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0”使用正负号“1”仅用绝对值（3）．Esel, type, item, comp, vmin, vmax, vinc, kabs 选择一组单元type: S: 选择一组单元（缺省）R: 在当前组中再选一部分作为一组A: 为当前组附加单元U: 在当前组中不选一部分单元All: 选所有单元None: 全不选Inve: 反向选择当前组Stat: 显示当前选择状态Item：Elem: 单元号Type: 单元类型号Mat: 材料号Real: 实常数号Esys: 单元坐标系号（4）. mp, lab, mat, co, c1,…….c4 定义材料号及特性lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens）ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号（缺省为当前材料号）c : 材料特性值，或材料之特性，温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数（5）. 定义DP材料：首先要定义EX和泊松比：MP，EX，MAT，……MP，NUXY，MAT，……定义DP材料单元表（这里不考虑温度）：TB，DP，MAT进入单元表并编辑添加单元表：TBDATA，1，CTBDATA，2，ψTBDATA，3，……如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：MP，EX，1，1E8MP，NUXY，1，0.3TB，DP，1TBDATA，1，27TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,deg（6）. 根据需要耦合某些节点自由度cp, nset, lab,,node1,node2,……node17nset: 耦合组编号lab: ux,uy,uz,rotx,roty,rotz ,allnode1-node17: 待耦合的节点号。

ANSYS 数据文件读写的APDL命令详解及实例作者：huright一 FORTRAN数据格式I格式(又叫整数格式)一般形式：Iw 或：Iw.m其中：w 一个数据占的位数宽度（又称“字段宽度”），m 需要输出的最少数字位数。

例1：（1）数字在指定的区域内向右端靠齐，如果数字位数比指定的字段宽度w小，则左边补以空格。

负数的符号也包含在字段宽度内。

（2）如果数字的位数超过了规定的字段宽度w，则不输出有效数据，而在该字段宽度范围内充满“*”符号。

（3）如果数字的位数超过了m,则按实际应输出的位数输出（但条件是不能超过w)。

m 不包括负号所占的一列。

F格式(又叫小数型格式)一般形式：Fw.dw 各数值占的总位数 d 输出数据的小数位数（小数点后的位数）。

例1：（1）数字在指定的区域内向右端靠齐，如果数字位数（含小数点和符号位）比指定的字段宽度w小，则左边补以空格；如果数字的位数超过了规定的字段宽度w，则不输出有效数据，而在该字段宽度范围内充满“*”符号。

（2）如果数据的小数位数比指定的小数位数d小，则在小数右边补0以凑足d位；如果小数位数大于d位，则输出时多于的小数位数按“四舍五入”规则舍去。

（3）假设b为数据整数部分的位数，则应使w≥b+d+1(小数点占一列），如果输出负数，则应保证w≥b+d+2(小数点和负号各占一列）。

（4）用F格式输出时应注意，由于难以事先确切估计出数据的大小，输出大的数时容易产生“宽度不够”的错误（由于w不够大），输出小的数时会出现丢掉有用数字的情况（由于d不够大而将后面的数字截去），这就是“大数印错，小数印丢”。

E格式(又叫指数型格式)一般形式：w.dw 各数值占的总位数，d 输出数据的小数位数（小数点后的位数）。

例1：（1）采取标准化的指数形式输出一个实数，d为以指数形式出现的数据的数字部分的小数位数。

（2）指数部分一般占4列，其中字母“E”和指数的符号各占一列，指数2列。

ANSYS 数据文件读写的APDL命令详解及实例作者：huright一 FORTRAN数据格式I格式(又叫整数格式)一般形式：Iw 或：Iw.m其中：w 一个数据占的位数宽度（又称“字段宽度”），m 需要输出的最少数字位数。

例1：（1）数字在指定的区域内向右端靠齐，如果数字位数比指定的字段宽度w小，则左边补以空格。

负数的符号也包含在字段宽度内。

（2）如果数字的位数超过了规定的字段宽度w，则不输出有效数据，而在该字段宽度范围内充满“*”符号。

（3）如果数字的位数超过了m,则按实际应输出的位数输出（但条件是不能超过w)。

m 不包括负号所占的一列。

F格式(又叫小数型格式)一般形式：Fw.dw 各数值占的总位数 d 输出数据的小数位数（小数点后的位数）。

例1：（1）数字在指定的区域内向右端靠齐，如果数字位数（含小数点和符号位）比指定的字段宽度w小，则左边补以空格；如果数字的位数超过了规定的字段宽度w，则不输出有效数据，而在该字段宽度范围内充满“*”符号。

（2）如果数据的小数位数比指定的小数位数d小，则在小数右边补0以凑足d位；如果小数位数大于d位，则输出时多于的小数位数按“四舍五入”规则舍去。

（3）假设b为数据整数部分的位数，则应使w≥b+d+1(小数点占一列），如果输出负数，则应保证w≥b+d+2(小数点和负号各占一列）。

（4）用F格式输出时应注意，由于难以事先确切估计出数据的大小，输出大的数时容易产生“宽度不够”的错误（由于w不够大），输出小的数时会出现丢掉有用数字的情况（由于d不够大而将后面的数字截去），这就是“大数印错，小数印丢”。

E格式(又叫指数型格式)一般形式：w.dw 各数值占的总位数，d 输出数据的小数位数（小数点后的位数）。

例1：（1）采取标准化的指数形式输出一个实数，d为以指数形式出现的数据的数字部分的小数位数。

（2）指数部分一般占4列，其中字母“E”和指数的符号各占一列，指数2列。

ANSYS 分析的命令流（APDL语言）前面分析的过程都是基于ANSYS用户图形界面(GUI命令流)形式，GUI 形式非常的直观明了，分析过程中通常会采用此种的分析方式。

但是这种分析方式会有一个缺点，就是当我们操作失误时候，没有后退的功能，因此我们不得不重新操作，这样就会给分析带来很多的麻烦。

但如果以命令流的形式进行计算的话就能很容易的减少这些麻烦。

下面就是本次点机前处理和分析计算过程的命令流/UNITS,SI !采用国际单位制/TITLE,2D DIANJI Static Analysis !定义分析名称KEYW,MAGNOD,1*AFUN,DEG !指定角度单位为度/PREP7 !进入前处理器ET,1,PLANE53 !设定单元类型为plane53EMUNIT,MKS !电磁单位LOCAL,12,0,0,0,0,90 !定义局部坐标系12--17LOCAL,13,0,0,0,0,30LOCAL,14,0,0,0,0,330LOCAL,15,0,0,0,0,270LOCAL,16,0,0,0,0,210LOCAL,17,0,0,0,0,150MP,MURX,1,1.0 ! 定义第一种材料的相对磁导率MP,MURX,2,1.0 ! 定义第二种材料的相对磁导率TB,BH,3TBPT,,35.03,0.1 ! 第三种材料的B-H磁化特性TBPT,,46.97,0.2TBPT,,57.32,0.3TBPT,,66.08,0.4TBPT,,74.04,0.5TBPT,,82.01,0.6TBPT,,91.56,0.7TBPT,,103.5,0.8TBPT,,117.83,0.9TBPT,,136.15,1TBPT,,160.83,1.1TBPT,,203.03,1.2TBPT,,286.62,1.3TBPT,,461.78,1.4TBPT,,955.41,1.5TBPT,,2547.8,1.6!TBPLOT,BH,3 !绘制B-H曲线TB,BH,4TBPT,,130,0.1 ! 定义第四种材料的B-H磁化特性TBPT,,170,0.2TBPT,,197,0.3TBPT,,218,0.4TBPT,,250,0.5TBPT,,290,0.6TBPT,,338,0.7TBPT,,400,0.8TBPT,,472,0.9TBPT,,570,1TBPT,,682,1.1TBPT,,810,1.2TBPT,,970,1.3TBPT,,1600,1.4TBPT,,2520,1.5TBPT,,3520,1.6!TBPLOT,BH,4 !绘制B-H曲线MP,MURX,5,1.17 ! 定义第五种材料的相对磁导率MP,MGXX,5,4.4E+005 ! 定义矫顽力矢量为X方向分量MP,MURX,6,1.17 ! 定义第六种材料的相对磁导率MP,MGXX,6,4.4E+005 ! 定义矫顽力矢量为X方向分量MP,MURX,7,1.17 ! 定义第七种材料的相对磁导率MP,MGXX,7,4.4E+005 ! 定义矫顽力矢量为X方向分量MP,MURX,8,1.17 ! 定义第八种材料的相对磁导率MP,MGXX,8,4.4E+005 ! 定义矫顽力矢量为X方向分量MP,MURX,9,1.17 ! 定义第九种材料的相对磁导率MP,MGXX,9,4.4E+005 ! 定义矫顽力矢量为X方向分量MP,MURX,10,1.17 ! 定义第十种材料的相对磁导率MP,MGXX,10,4.4E+005 ! 定义矫顽力矢量为X方向分量RECTNG,-0.00275,0.00275,0.0485,0.062 ! 创建一个矩形RECTNG,-0.00025,0.00025,0.048,0.0485 ! 再创建一个矩形AADD,ALL ! 将两个矩形融合在一起，形成槽结构CSYS,1 !转动工作平面到柱坐标下AGEN,36,3,38,1,,10,0 ! 复制槽/PNUM,AREA,ALL ! 给所有槽编号PCIRC,0.048,0.074 ! 创建定子圆环AOVLAP,ALL !对定子圆环和槽进行交迭布尔操作使得在连接线上共节点PCIRC,0.037,0.047,78,102 ! 创建永磁体极CSYS,1 !转到柱坐标下AGEN,6,3,8,1,,60,0 ! 复制6个磁极PCIRC,0.023,0.047 ! 创建转子圆环AOVLAP,3,10,45,51,57,63,69 ! 对转子圆环与磁极进行交迭操作PCIRC,0.047,0.0475 !创建两层气隙PCIRC,0.0475,0.048AOVLAP,ALLAGLUE,ALL ! 将所有交接面黏在一起NUMCMP,AREA !压缩面编号ASEL,S,AREA,,1,35 !选中槽ASEL,A,AREA,,42AATT,1 !材料定义为1号材料ASEL,S,AREA,,45ASEL,A,AREA,,43 !选中气隙AATT,2 !材料定义为2号材料ASEL,S,AREA,,46 !选中定子铁芯AATT,3 !材料定义为3号材料ASEL,S,AREA,,44 !选中转子磁轭AATT,4 !材料定义为4号材料ASEL,S,AREA,,40 !选中第1-6个永磁体磁极AATT,5,,,12 !材料定义为5-10号材料ASEL,S,AREA,,39AATT,6,,,13ASEL,S,AREA,,38AATT,7,,,14ASEL,S,AREA,,41AATT,8,,,15ASEL,S,AREA,,37AATT,9,,,16ASEL,S,AREA,,36AATT,10,,,17/PNUM,MAT,1 !打开材料编号APLOT !重新显示ALLSEL,ALL !选择所有实体SMRTSIZE,1 !指定智能网格划分等级L AMESH,ALL !智能网格划分!ESEL,S,MAT,,1!CM,ARM,ELEM!FMAGBC,'ARM'!SA VE,DIANJI_2D_MESH.db !保存!SA VEFINISH !完成网格划分/SOLU !进入求解器ANTYPE,STATIC !选择静态磁场分析NROPT,AUTO !采用牛顿---拉夫森迭代方法!ASEL,S,AREA,,1,35 !选中槽!ASEL,A,AREA,,42!BFE,ALL,JS,1,,,15000000 !施加电流密度LOCAL,11,1 ! 定义柱坐标系FLST,2,8,4,ORDE,4FITEM,2,289FITEM,2,-292FITEM,2,561FITEM,2,-564DL,P51X, ,AZ,0, !加边界条件ALLSEL,ALLMAGSOLV!SA VE,EMAGE_2D_MESH.db !保存!SA VEFINISH/POST1PLF2D,27,0,10,1 !显示磁力线图FMAGSUM,'ARM' ！对电磁力求和plvect,B,,,,VECT,ELEM,ON,0 !显示磁感应强度矢量plvect,H,,,,VECT,ELEM,ON,0 !显示磁场强度矢量/GRAPHICS,POWERA VRES,2PLNSOL,B,SUM !显示磁通密度等值云图flishESEL,S,MAT,,2CM,ARM,ELEM!生成一个组件FMAGSUM,'ARM'！以下命令流是保存图片时使用，使其背景为白，结合PlotCtrlsHard ﹥Copy命令可设置格式/GRAPHICS,POWER/RGB,INDEX,100,100,100, 0/RGB,INDEX, 80, 80, 80,13/RGB,INDEX, 60, 60, 60,14/RGB,INDEX, 0, 0, 0,15/REPLOT !背景反白显示/SHOW,JPEG,,0JPEG,QUAL,75,JPEG , ORIENT, HORIZJPEG , Color, 2JPEG , TMOD, 1/GFILE,600,EPLOT/SHOW,CLOSE/DEVICE,VECTOR,0。

APDL命令流：将ansys分析结果输出为tecplot格式/post1file,Jobname,rst !指明从哪⼀个结果⽂件中读取数据*get,NodeNum,NODE,0,COUNT !得到模型的所有节点数*get,nd,NODE,0,NUM,MIN !得到模型的最⼩节点编号*dim,nodes,array,NodeNum !定义⼀个存储节点的数组*dim,xyz,array,NodeNum,3 !定义⼀个存储节点坐标的数组*dim,NodeTemp,array,NodeNum !定义⼀个存储节点的数组set,4,1 !读⼊第四载荷步、第⼀个⼦步的结果*do,i,1,NodeNum,1nodes(i)=nd !将节点编号存储在nodes数组中xyz(i,1)=NX(nd) !在xyz(i,1)中存储节点的X坐标xyz(i,2)=NY(nd) !在xyz(i,1)中存储节点的Y坐标xyz(i,3)=NZ(nd) !在xyz(i,1)中存储节点的Z坐标*Get,NodeTemp(i),NODE,nodes(i),TEMPnd=NDNEXT(nd) !取得下⼀个节点编号*enddo*get,ElementNum,ELEM,0,COUNT !得到模型的所有单元数*get,el,ELEM,0,NUM,MIN !得到模型的最⼩单元编号*dim,elems,array,ElementNum !定义⼀个存储单元的数组*dim,ndlst,array,ElementNum,8 !定义⼀个存储节点坐标的数组*do,i,1,ElementNum,1elems(i)=el !将单元编号存储在elems数组中*do,j,1,8,1ndlst(i,j)=NELEM(el,j) !将单元对应的节点编号存储在ndlst数组中*enddoel=ELNEXT(el) !取得下⼀个单元编号*enddo*cfopen,mesh,dat !输出的⽂件名为：mesh.dat!写TecPlot⽂件的⽂件头*vwrite(’TITLE=Ansys Temperature Analysis’)!写变量名*vwrite(’VARIABLES=”X”,”Y”,”Z”,”Temp”‘)!写节点数和单元数*vwrite,NodeNum,ElementNum(’ZONE n=’f6.1,’ e=’f6.1,’ f=fepoint’,' et=brick’)!写节点的坐标和温度值*vwrite,xyz(1,1),xyz(1,2),xyz(1,3),NodeTemp(1)(4f12.6)!写⼀个空⾏*vwrite(’ ‘)!写各单元对应的节点号*vwrite,ndlst(1,1),ndlst(1,2),ndlst(1,3),ndlst(1,4),ndlst(1,5),ndlst(1,6),ndlst(1,7),ndlst(1,8)(8f6.0)!关闭⽂件*CFCLOSEfinish。

ANSYS 数据文件读写的APDL命令详解及实例源自机知网作者：huright一 FORTRAN数据格式I格式(又叫整数格式)一般形式：Iw 或：Iw.m其中：w 一个数据占的位数宽度（又称“字段宽度”），m 需要输出的最少数字位数。

例1：（1）数字在指定的区域内向右端靠齐，如果数字位数比指定的字段宽度w小，则左边补以空格。

负数的符号也包含在字段宽度内。

（2）如果数字的位数超过了规定的字段宽度w，则不输出有效数据，而在该字段宽度范围内充满“*”符号。

（3）如果数字的位数超过了m,则按实际应输出的位数输出（但条件是不能超过w)。

m 不包括负号所占的一列。

F格式(又叫小数型格式)一般形式：Fw.dw 各数值占的总位数 d 输出数据的小数位数（小数点后的位数）。

例1：（1）数字在指定的区域内向右端靠齐，如果数字位数（含小数点和符号位）比指定的字段宽度w小，则左边补以空格；如果数字的位数超过了规定的字段宽度w，则不输出有效数据，而在该字段宽度范围内充满“*”符号。

（2）如果数据的小数位数比指定的小数位数d小，则在小数右边补0以凑足d位；如果小数位数大于d位，则输出时多于的小数位数按“四舍五入”规则舍去。

（3）假设b为数据整数部分的位数，则应使w≥b+d+1(小数点占一列），如果输出负数，则应保证w≥b+d+2(小数点和负号各占一列）。

（4）用F格式输出时应注意，由于难以事先确切估计出数据的大小，输出大的数时容易产生“宽度不够”的错误（由于w不够大），输出小的数时会出现丢掉有用数字的情况（由于d不够大而将后面的数字截去），这就是“大数印错，小数印丢”。

E格式(又叫指数型格式)一般形式：w.dw 各数值占的总位数，d 输出数据的小数位数（小数点后的位数）。

例1：（1）采取标准化的指数形式输出一个实数，d为以指数形式出现的数据的数字部分的小数位数。

（2）指数部分一般占4列，其中字母“E”和指数的符号各占一列，指数2列。

ANSYS房屋框架结构模型命令流APDL!********************提示********************************************!第一步：建立框架结构模型!第二步：框架结构施加重力荷载!第三步：框架结构施加活荷载!第四步：框架结构施加风荷载!第五步：框架结构荷载组合运算!第六步：框架结构模态分析!第七步：框架结构屈曲特征值分析!第八步：框架结构地震时程弹性分析!第九步：结束!采用地震波时间文件为:TIME.TXT!采用地震波文件分别为:AC_X.TXT/AC_Y.TXTFINISH/CLEAR/FILENAME,FRAME_SHEAR_WALL/TITLE, FRAME_SHEAR_WALL!************************************************************** ************ !************************第一步******************************************** !********************框架结构建模型***************************************/PREP7!采用单位为N/m/S 国际单位制!定义单元类型ET,1,BEAM4ET,2,SHELL63!定义实常数R,1,0.25,0.0052,0.0052,0.5,0.5R,2,0.2025,0.0034,0.0034,0.45,0.45 R,3,0.06,0.0002,0.00045,0.3,0.2 R,4,0.25,0.25,0.25,0.25R,5,0.1,0.1,0.1,0.1!定义材料参数MP,EX,1,3.0E10MP,NUXY,1,0.2MP,DENS,1,2500!关键点*DO,II,1,1K, 2+(II-1)*20, 6, 0, (II-1)*6K, 3+(II-1)*20, 12, 0, (II-1)*6K, 4+(II-1)*20, 17, 0, (II-1)*6K, 5+(II-1)*20, 0, 5, (II-1)*6K, 6+(II-1)*20, 6, 5, (II-1)*6K, 7+(II-1)*20, 12, 5, (II-1)*6K, 8+(II-1)*20, 17, 5, (II-1)*6K, 9+(II-1)*20, 0, 10, (II-1)*6K, 10+(II-1)*20, 6, 10, (II-1)*6K, 11+(II-1)*20, 12, 10, (II-1)*6 K, 12+(II-1)*20, 17, 10, (II-1)*6 K, 14+(II-1)*20, 6, 1.875, (II-1)*6 K, 16+(II-1)*20, 6, 3.125, (II-1)*6 K, 18+(II-1)*20, 6, 1.875, 3.9K, 20+(II-1)*20, 6, 3.125, 3.9*ENDDO*DO,II,1,10K, 1+II*20, 0, 0, (II-1)*3+6K, 2+II*20, 6, 0, (II-1)*3+6K, 3+II*20, 12, 0, (II-1)*3+6K, 4+II*20, 17, 0, (II-1)*3+6K, 5+II*20, 0, 5, (II-1)*3+6K, 6+II*20, 6, 5, (II-1)*3+6K, 7+II*20, 12, 5, (II-1)*3+6K, 8+II*20, 17, 5, (II-1)*3+6K, 9+II*20, 0, 10, (II-1)*3+6K, 10+II*20, 6, 10, (II-1)*3+6K, 11+II*20, 12, 10, (II-1)*3+6K, 12+II*20, 17, 10, (II-1)*3+6K, 13+II*20, 0, 1.875, (II-1)*3+6K, 14+II*20, 6, 1.875, (II-1)*3+6K, 15+II*20, 0, 3.125, (II-1)*3+6K, 16+II*20, 6, 3.125, (II-1)*3+6K, 17+II*20, 0, 1.875, (II-1)*3+6+1.5 K, 18+II*20, 6, 1.875, (II-1)*3+6+2 K, 19+II*20, 0, 3.125, (II-1)*3+6+1.5 K, 20+II*20, 6, 3.125, (II-1)*3+6+2 *ENDDO*DO,II,11,11K, 1+II*20, 0, 0, (II-1)*3+6K, 2+II*20, 6, 0, (II-1)*3+6K, 3+II*20, 12, 0, (II-1)*3+6K, 4+II*20, 17, 0, (II-1)*3+6K, 6+II*20, 6, 5, (II-1)*3+6K, 7+II*20, 12, 5, (II-1)*3+6K, 8+II*20, 17, 5, (II-1)*3+6K, 9+II*20, 0, 10, (II-1)*3+6K, 10+II*20, 6, 10, (II-1)*3+6K, 11+II*20, 12, 10, (II-1)*3+6K, 12+II*20, 17, 10, (II-1)*3+6K, 13+II*20, 0, 1.875, (II-1)*3+6K, 14+II*20, 6, 1.875, (II-1)*3+6K, 15+II*20, 0, 3.125, (II-1)*3+6K, 16+II*20, 6, 3.125, (II-1)*3+6K, 17+II*20, 0, 1.875, (II-1)*3+6+1.8 K, 18+II*20, 6, 1.875, (II-1)*3+6+2.5 K, 19+II*20, 0, 3.125, (II-1)*3+6+1.8 K, 20+II*20, 6, 3.125, (II-1)*3+6+2.5 *ENDDO*DO,II,12,12K, 1+II*20, 0, 0, 39.5K, 2+II*20, 6, 0, 39.5K, 3+II*20, 12, 0, 39.5K, 4+II*20, 17, 0, 39.5K, 5+II*20, 0, 5, 39.5K, 6+II*20, 6, 5, 39.5K, 7+II*20, 12, 5, 39.5K, 8+II*20, 17, 5, 39.5K, 9+II*20, 0, 10, 39.5K, 10+II*20, 6, 10, 39.5K, 11+II*20, 12, 10, 39.5K, 12+II*20, 17, 10, 39.5*ENDDO*DO,II,2,13L, (II-1)*20+1, (II-1)*20+2L, (II-1)*20+2, (II-1)*20+3L, (II-1)*20+3, (II-1)*20+4L, (II-1)*20+5, (II-1)*20+6L, (II-1)*20+6, (II-1)*20+7L, (II-1)*20+7, (II-1)*20+8L, (II-1)*20+9, (II-1)*20+10L, (II-1)*20+10, (II-1)*20+11L, (II-1)*20+11, (II-1)*20+12L, (II-1)*20+5, (II-1)*20+9L, (II-1)*20+6, (II-1)*20+10 L, (II-1)*20+3, (II-1)*20+7 L, (II-1)*20+4, (II-1)*20+8 L, (II-1)*20+8, (II-1)*20+12 *ENDDO*DO,II,1,12L, (II-1)*20+1, II*20+1L, (II-1)*20+2, II*20+2L, (II-1)*20+3, II*20+3L, (II-1)*20+4, II*20+4L, (II-1)*20+5, II*20+5L, (II-1)*20+6, II*20+6L, (II-1)*20+7, II*20+7L, (II-1)*20+8, II*20+8L, (II-1)*20+9, II*20+9L, (II-1)*20+10, II*20+10 L, (II-1)*20+11, II*20+11 L, (II-1)*20+12, II*20+12*ENDDO*DO,II,1,1L, (II-1)*20+14, (II-1)*20+18 L, (II-1)*20+16, (II-1)*20+20 *ENDDO*DO,II,2,12L, (II-1)*20+13, (II-1)*20+17 L, (II-1)*20+15, (II-1)*20+19 L, (II-1)*20+14, (II-1)*20+18 L, (II-1)*20+16, (II-1)*20+20 *ENDDOA, 1, 21, 25, 5*DO,II,2,12A, (II-1)*20+1, II*20+1, II*20+5, (II-1)*20+5, (II-1)*20+15, (II-1)*20+19, (II-1)*20+17, (II-1)*20+13 *ENDDO*DO,II,1,12A, (II-1)*20+2, II*20+2, II*20+6, (II-1)*20+6, (II-1)*20+16, (II-1)*20+20, (II-1)*20+18, (II-1)*20+14 *ENDDO*DO,II,1,12A, II*20+1, II*20+2, II*20+6, II*20+5A, II*20+2, II*20+3, II*20+7, II*20+6A, II*20+5, II*20+6, II*20+10,II*20+9A, II*20+6, II*20+7, II*20+11,II*20+10A, II*20+7, II*20+8, II*20+12,II*20+11*ENDDO!立柱网格划分LSEL,S,,,181,192 !底层立柱LATT,1,1,1LESIZE,ALL,,,10LMESH,ALLLSEL,S,,,192,324 !2～12层立柱LATT,1,2,1LESIZE,ALL,,,5LMESH,ALL!梁网格划分LSEL,S,,,1,180LATT,1,3,1,LESIZE,ALL,,,8LMESH,ALL!剪力墙网格划分LSEL,S,,,372 !外剪力墙两端、内剪力墙顶端LSEL,A,,,413LSEL,A,,,461LESIZE,ALL,,,8LSEL,s,,,418 !内剪力墙底端LSEL,A,,,420LESIZE,ALL,,,3LSEL,s,,,325 !内剪力墙底端内侧LSEL,A,,,326LESIZE,ALL,,,6LSEL,S,,,375 !剪力墙空洞顶部LSEL,A,,,379LSEL,A,,,383LSEL,A,,,387LSEL,A,,,391LSEL,A,,,395LSEL,A,,,399LSEL,A,,,403LSEL,A,,,407LSEL,A,,,415LSEL,A,,,419LSEL,A,,,423LSEL,A,,,427LSEL,A,,,431LSEL,A,,,435LSEL,A,,,439LSEL,A,,,443LSEL,A,,,447LSEL,A,,,451LSEL,A,,,455LSEL,A,,,459LSEL,A,,,463LESIZE,ALL,,,2!剪力墙网格划分ASEL,S,,,1,24AATT,1,4,2AMESH,ALL!楼板网格划分ASEL,S,,,25,96AATT,1,5,2AMESH,ALLNSEL,S,LOC,Z,0 ！选取模型底端节点D,ALL,ALL ！施加位移约束ALLSEL,ALL ！重新选取所有节点/eshape,1.0/VIEW,1,1,1,1/ANG,1,270,XM,0/REPlotFINISH!************************************************************** ************ !************************第二步******************************************** !********************框架结构施加重力荷载*********************************ANTYPE,STATICNSEL,S,LOC,Z,0 ！选取模型底端节点D,ALL,ALL ！施加位移约束ALLSEL,ALL ！重新选取所有节点ACEL,0,0,9.8SOLVEFINISH/POST1SET,FIRSTPLNSOL,U,Z,0,1FINISH!************************************************************** ************ !************************第三步******************************************** !*************框架结构施加楼面活荷载D=3KN/M^2***************************** !框架结构施加楼面活荷载D=3KN/M^2FINISH/SOLUANTYPE,STATICNSEL,S,LOC,Z,0 ！选取模型底端节点D,ALL,ALL ！施加位移约束ALLSEL,ALL ！重新选取所有节点ACEL,0,0,9.8SOLVEFINISH/POST1SET,FIRSTPLNSOL,U,Z,0,1FINISH!************************************************************** ************ !************************第四步******************************************** !***********框架结构施加风荷载(基本风压=0.25KN/M^2,体形系数为1.0********** ANTYPE,STATIC*DIM,LOAD_1,ARRAY,12LOAD_1(1)=3.78,2.16,2.39,2.57,2.72,2.84,2.95,3.17,3.20,3.29, 3.39,3.51*DIM,LOAD_2,ARRAY,12LOAD_2(1)=7.56,4.32,4.78,5.14,5.44,5.68,5.90,6.34,6.40,6.58, 6.76,7.02*DIM,LOAD_3,ARRAY,12LOAD_3(1)=6.93,3.96,4.38,4.71,4.99,5.21,5.41,5.81,5.87,6.03, 6.20,6.44*DIM,LOAD_4,ARRAY,12LOAD_4(1)=3.15,1.80,1.99,2.14,2.27,2.37,2.46,2.64,2.67,2.74, 2.82,2.93*DIM,LOAD_A,ARRAY,12LOAD_A(1)=3.15,1.80,1.99,2.14,2.27,2.37,2.46,2.64,2.67,2.74, 2.82,2.93*DIM,LOAD_B,ARRAY,12LOAD_B(1)=6.30,3.60,3.98,4.28,4.54,4.74,4.92,5.28,5.34,5.48, 5.64,5.86*DIM,LOAD_C,ARRAY,12LOAD_C(1)=3.15,1.80,1.99,2.14,2.27,2.37,2.46,2.64,2.67,2.74, 2.82,2.93*DO,II,1,12FK,1+20*II,FY,LOAD_1(II)*ENDDO*DO,II,1,12FK,2+20*II,FY,LOAD_2(II)*ENDDO*DO,II,1,12FK,3+20*II,FY,LOAD_3(II)*ENDDO*DO,II,1,12FK,4+20*II,FY,LOAD_4(II)*ENDDOSOLVEFKDELE,ALL,ALL !第二次施加风荷载,删除第一次的。

ANSYS中的APDL命令ANSYS中的APDL命令(一)（1）．Lsel, type, item, comp, vmin, vmax, vinc, kswp 选择线 type: s 从全部线中选一组线r 从当前选中线中选一组线a 再选一部线附加给当前选中组aunoneu(unselect)inve: 反向选择item: line 线号loc 坐标length 线长comp: x,y,zkswp: 0 只选线1 选择线及相关关键点、节点和单元（2）．Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点type: S: 选择一组新节点（缺省）R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0”使用正负号“1”仅用绝对值（3）．Esel, type, item, comp, vmin, vmax, vinc, kabs 选择一组单元type: S: 选择一组单元（缺省）R: 在当前组中再选一部分作为一组A: 为当前组附加单元U: 在当前组中不选一部分单元All: 选所有单元None: 全不选Inve: 反向选择当前组Stat: 显示当前选择状态Item： Elem: 单元号Type: 单元类型号Mat: 材料号Real: 实常数号Esys: 单元坐标系号（4）. mp, lab, mat, co, c1,…….c4定义材料号及特性lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens）ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号（缺省为当前材料号）c : 材料特性值，或材料之特性，温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数（5）. 定义DP材料：首先要定义EX和泊松比：MP，EX，MAT，……MP，NUXY，MAT，……定义DP材料单元表（这里不考虑温度）：TB，DP，MAT进入单元表并编辑添加单元表：TBDATA，1，CTBDATA，2，ψTBDATA，3，…… 如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：MP，EX，1，1E8MP，NUXY，1，0.3TB，DP，1TBDATA，1，27TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,deg（6）. 根据需要耦合某些节点自由度cp, nset, lab,,node1,node2,……node17nset: 耦合组编号lab: ux,uy,uz,rotx,roty,rotz ,allnode1-node17: 待耦合的节点号。

ANSYS-APDL的实例学习10.1有限元模型的建立/TITLE，3D analysis on shield tunnel in Metro ! 确定分析标题/NOPR !菜单过滤设置/PMETH，OFF，0KEYW，PR_SET，1KEYW，PR_STRUC，1 !保留结构分析部分菜单/COM，/COM，Preferences for GUI filtering have been set to display: /COM，Structural1．材料、实常数和单元类型定义/clear !更新数据库/prep7 !进入前处理器et,1,solid45 !设置单元类型et,2,mesh200,6save !保持数据（2）定义模型中的材料参数。

!土体材料参数mp,ex,1,3.94e6 !地表层土弹性模量mp,prxy,1,0.35 !地表层土泊松比mp,dens,1,1828 !地表层土密度mp,ex,2,20.6e6 !盾构隧道所在地层参数mp,prxy,2,0.30mp,dens,2,2160mp,ex,3,500e6 !基岩地层参数mp,prxy,3,0.33mp,dens,3,2160!管片材料参数，管片衬砌按各向同性计算mp,ex,4,27.6e9 !管片衬砌弹性模量mp,prxy,4,0.2 !管片衬砌泊松比mp,dens,4,2500 !管片衬砌密度!注浆层，参数按水泥土取值mp,ex,5,1e9 !注浆层弹性模量mp,prxy,5,0.2 !注浆层泊松比mp,dens,5,2100 !注浆层密度save !保持数据第1章大型有限元软件ANSYS简介2 2．建立平面内模型并划分单元（1）在隧道中心线定义局部坐标，便于后来的实体选取。

local,11,0,0,0,0 !局部笛卡儿坐标local,12,1,0,0,0 !局部极坐标csys,11 !将当前坐标转换为局部坐标wpcsys,-1 !同时将工作平面转换到局坐标cyl4,,,,,2.7,90 !画部分圆半径为2.7cyl4,0,0,2.7,0,3,90 !画管片层部分圆cyl4,0,0,3,0,3.2,90 !画注浆层部分圆rectng,0,4.5,0,4.5 !画外边界矩形aovlap,all !做面递加nummrg,all !合并所有元素numcmp,all !压缩所有元素编号rectng,4.5,31.5,0,4.5 !画矩形面nummrg,all !合并所有元素numcmp,all !压缩所有元素编号save !保持数据（2）划分单元，如图10-1所示。

在ANSYS中，命令流是由一条条ANSYS的命令组成的一个命令组合，这些命令按照一定顺序排布，能够完成一定的ANSYS功能，这些功能一般来说通过菜单操作也能够实现（而那些命令流能够实现，菜单操作实现不了的单个命令比较少见）。

以下命令是结合我自身经验，和前辈们的一些经验而总结出来的，希望对大家有帮助。

（1）．Lsel, type, item, comp, vmin, vmax, vinc, kswp 选择线type: s 从全部线中选一组线r 从当前选中线中选一组线a 再选一部线附加给当前选中组aunoneu(unselect)inve: 反向选择item: line 线号loc 坐标length 线长comp: x,y,zkswp: 0 只选线1 选择线及相关关键点、节点和单元（2）．Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点type: S: 选择一组新节点（缺省）R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0”使用正负号“1”仅用绝对值（3）．Esel, type, item, comp, vmin, vmax, vinc, kabs 选择一组单元type: S: 选择一组单元（缺省）R: 在当前组中再选一部分作为一组A: 为当前组附加单元U: 在当前组中不选一部分单元All: 选所有单元None: 全不选Inve: 反向选择当前组Stat: 显示当前选择状态Item：Elem: 单元号Type: 单元类型号Mat: 材料号Real: 实常数号Esys: 单元坐标系号（4）. mp, lab, mat, co, c1,…….c4 定义材料号及特性lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens）ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号（缺省为当前材料号）c : 材料特性值，或材料之特性，温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数（5）. 定义DP材料：首先要定义EX和泊松比：MP，EX，MAT，……MP，NUXY，MAT，……定义DP材料单元表（这里不考虑温度）：TB，DP，MAT进入单元表并编辑添加单元表：TBDATA，1，CTBDATA，2，ψTBDATA，3，……如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：MP，EX，1，1E8MP，NUXY，1，0.3TB，DP，1TBDATA，1，27TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,deg（6）. 根据需要耦合某些节点自由度cp, nset, lab,,node1,node2,……node17nset: 耦合组编号lab: ux,uy,uz,rotx,roty,rotz ,allnode1-node17: 待耦合的节点号。