APDL学习笔记

1、学习好资料欢迎下载APDL参数化语言（脚本语言）1.1简介程序的输入可根据指定的函数、变量以及选用的分析类型来确定，是进行优化设计和自适应 网格的最主要基础。APDL具有下列功能：（1）标量参数；（2）数组参数；（3）表达式和函数；（4）分支和循环；（5）重复功能和缩 写；（6）宏；（7）用户程序优点：在重复计算时，可以利用第一次（前一次）分析时生成的LOG文件，对这个文件作适当的修改即可得到自己的命令流文件，再添加些APDL控制命令，就得到 APDL命令的文件。在使用APDL命令流方式中，当 选择实体 模型时，建议先选择，再施加操作；如果网 格发生变化，最好根据几何特征、材料属性来选。操作命

2、令按使用范围分类:r特定环境下使用（在任何处理器中使用2.2参数化变量使用参数即APDL的变量，不必定义参数类型，默认为双精度型。参数类型有两种：标量 和数组，其中数组参数包括：数值型、字符型、字符串型和表格类型等。2.2.1参数名的命名规则和格式1参数的命名参数名的命名约定为：（1）必须以字母开头（2）只能包含字母、数字和下划线；（3）长度不能超过32个字符。参数名不能和ANSYS已有的标签名相同。名称ARG1至ARG9和AR10至AR99被保存为局部参数，一般只能在宏中使用，因此 不要将它们作为“通用”的参数使用。不要用下划线“ _”作为参数名的开头。2注意：1）2）3）2参数的定义（1）

3、在执行过程中指定参数值，可用“*SET ”命令对指定的参数赋值，格式为：*SET ,n ame,valve也可以用“=”调用“ *SET”。格式：Name=Value Name参数名，Value参数值（数值型或字符型，字符型时用括起来且长度不能超过8个字符）（2）从ANSYS系统中获取参数的值的两种方法：*GET命令和在线“取值函数”*GET命令格式：*GET，Par, Entity，ENTNUM，Item1，IT1NUM，Item2，IT2NUM （6 个部分，空的用， 代替）其中：Par：贮存返回值。Entity :将要返回项的关键词，女0 NODE,ELEM,KP,LINE,VOLU等。

4、ENTNUM :实体编号，若为 0则代表所有实体。Item1 :对于一个特定的实体，其项的名称。“取值函数”通过取值函数能够取出某项的值，并将它直接投入到当前软件的运行中。（取值函数见ANSYS参数化编程与命令手册P18）2.2.4参数表达式与函数1. 参数表达参数表达由参数、数字和加减乘除等运算符组成。运算符：+、-、* （乘）、/、* （乘方）运算顺序：1）先算圆括号；2）指数运算；3）乘法或除法；4）符号运算；5）加法或减法；6）逻辑运算 注意：圆括号能够嵌套 4层，每对圆括号内可以完成 9次运算。在表达式得算子之间要 避免 出现空格，特别在前不要有空格。注释符号最好采用“！”开头。2参

5、数函数包括：正弦余弦、开根号、取绝对值等等（见 ANSYS参数化编程与命令手册P29）*AFUN,RAD角度单位（弧度和度的转化）：改为度：*AFUN,DEG ;改为弧度：2.3参数化数组1. 数组类型四维（行、列、页、书）、按大小可分为：一维（行）、二维（行、列）、三维（行、列、页）、五维（行、列、页、书、书架）按数组元素类型分为：1开始的连续整数,数值型数组（ARRAY ）:默认数组类型，其行标、列标和页标必须是以数组元素的内容可以是整数或实数。字符型数组（CHAR ）:每个元素都由不超过8个字符的字符值组成，其行标、列标和页标必须是以1开始的连续整数。列、页标可以用实数定义，数组元素的内

6、容可以是整数或实数。特点见ANSYS参数化编程与命令手册命令：*DIM ，Par，Type, IMAX,JMAX,KMAX,Var1, Var2,表格型数组（TABLE ）:行、P31。字符串型数组（STRING ）:表格数组的生成方法：Var3,CSYSID2. 数组的定义（确定数组类型和大小）命令：*DIM, Par,Ty pe,IMAX,JMAX,KMAX,Var1, Var2, Var3,CSYSIDPar:指定的参数名Type：参数类型（ARR4指4维的数值型数组）IMAX,JMAX,KMAX:行列页标号的范围，默认值为1Var1, Var2, Var3 :标题（对于表格型的默认值为

7、行列页）3数组元素的赋值1）指定单个数组元素的值“=”进行赋值，每个“=”可以给10个数组元素给数组的某列赋值，可用 *SET命令或者赋值。4数组参数的运算1）向量运算（列向量）叉积等）矩阵x,y第一列相加，结果放在向量 z中*VOPER:对两个数组向量进行计算（求和，例：*VOPER , z（1）,x（1,1）,ADD,y（1,1）*VFUN :对一个数组向量进行函数运算对向量的第2列的每个元素进行平方运算，结果例：*VFUN,a3（1）,PWR,x（1,2）,2放在a3中2）矩阵运算数值型数组参数之间的数学运算，相乘，转置，解方程组等2.5.3 DO 循环1. DO-LOOPS格式：*DO

8、, Par,IVAL,FVAL,INCENDDOPar：循环变量名称（不能使用字符参数）IVAL,FVAL,INC :循环变量的初值、终值，增量，INC默认值为1,也可以是负值或实数。2.DO-WHILE 循环当循环次数不确定，满足某个条件循环才结束。*DOWHILE ,P arm*ENDDO3隐含的“：”循环格式：（X:Y:Z ）其中Z的默认值为1例如：n,（1:6）,（2:12:2）n,6,12执行6个命令：n,1,2n,2,42.5.4有条件分支（*IF）格式：*IF,VAL1，Oper,VAL2,BaseVAL1:比较过程中第1个数值或者数值型参数。Oper：比较运算：EQ（等于），N

9、E （不等于）、（LT）小于、（GT）大于、LE（小于或等于）、GE （大于或等于）、（ABLT ）绝对值小于、（ABGT ）绝对值大于VAL2:比较过程中第2个数值或者数值型参数。Base :当比较的结果成立（即为真）时，将要发生的动作1）如果将THEN赋给Base,那么*IF命令将成为一个IF-THEN-ELSE结构*IF 条件 THEN*ELSEIF*ELSE*ENDIFANSYS。2）如果将STOP赋给Base那么程序可根据某特定的条件退出前处理器4.1建立实体模型4.1.1生成关键点1“ K ”命令 功能：按坐标生成一个关键点格式：K,N PT,X,Y,Z0，自动指定当前可用关键点编

10、号的最小值；如果指定一个有相注意：只有当关键点没有依附在高级图元和划分网NPT:关键点编号，如果为同编号的关键点，将覆盖前一次的定义， 格之前，其坐标值才能重新指定。2.“ KL”命令功能：在已存在的线上按指定位置生成一个关键点。格式：KL,NL1，RATI0,NK1来说，线的方向将方向。0.01.0，默认为 0.5。NL1:线的编号。如果是一个负数，则对变量RATIORATIO :生成关键点的位置与线长之间的比率，其值NK1 :关键点的编号4.1.2生成线LSTR命令功能：由两点生成一条直线格式：LSTR，P1, P2XV2 , YV2 , ZV2 ,“L”命令 由两点生成一条线（直线或弧形

11、） 格式：L,P1，P2, NDIV,SPACE,XV1 ，YV1，ZV1，P1，P2:线开始、结尾的关键点编号LESIZENDIV :该线上将要划分单元的等分数，通常不用。划分线段用：SPACE:A间隔比率（通常不用）XV1，YV1，ZV1 :与直线P1端点斜率相关的斜率矢量位置（激活坐标系下）XV2，YV2，ZV2 :与直线P2端点斜率相关的斜率矢量位置LARC命令功能：三点生成一段圆弧格式：LARC ,P 1, P2,P C,RADP1, P2:圆弧线始端、末端关键点的编号PC:定义弧平面或线曲率中心的关键点，不能在P1,P2直线上，可以不在曲率中心。弧形的曲率半径，如果为空，系统自动计

12、算。CIRCLE命令 功能：生成圆弧线 格式：CIRCLE, PCENT,RAD, PAXIS, PZERO,ARC,NSEGPCENT:圆中心的关键点RAD :圆弧半径PAXIS :定义圆轴线的关键点，如果默认，轴线与工作面正交 PZERO :定义与圆正交平面的关键点ARC :圆弧的长度（度） NSEG :沿圆周生成的线段数4.1.3生成面1.“ A ”功能：由关键点生成任意形状的面格式：A,P1 , P2, P3, P4,.P18（至少三个关键点，最多 18），关键点必须按照顺时针或逆时针顺序输入。“ AL ”命令功能：通过已定义的边界线生成一个面格式：AL,L1 , L2 , L3 ,

13、L4,L10（至少 3 条）注意：线号的顺序可以任意，但要收尾相接。“ BLC5 ”命令功能：通过中心点和角点生成矩形面或块体格式：BLC5 , XCENTER,YCENTER,WIDTH , HEIGHTQEPTHXCENTER,YCENTER :矩形面或块体中心的坐标值WIDTH :定义矩形面和块体的宽度，与X轴平行HEIGHT :定义矩形面和块体的高度，与丫轴平行DEPTH :块体深度，与 Z轴平行“ BLC4 ”命令通过两个角点或 Z方向的深度生成一个矩形或块体格式：BLC4 , XCORNER,YCORNER , WIDTH , HEIGHT,DEPTH矩形面或块体第一个角点在工作面

14、的X,Y坐标值。“ RECTNG ”命令功能：在工作平面上生成平面格式：RECTNG,X1 , X2, Y1 , Y2X1 , X2 :矩形在工作平面 X方向坐标值变化范围Y1 , Y2 :矩形在工作平面 丫方向坐标值变化范围CYL4 ”命令功能：在工作平面生成一个圆面或圆柱体格式：CYL4 , XCENTER,YCENTER , RAD1 , THETA1 , RAD2 , THETA2 , DEPTH XCENTER,YCENTER:圆面或圆柱体中心 X,Y值RAD1 , RAD2 :圆面或圆柱体的内、外半径（任意一个为0或空，或相等即为圆或圆柱）THETA1 , THETA2，圆面或圆柱

15、体开始、结束的角度。CYL5 ”命令：功能：通过端点生成一个圆形和圆柱体格式：CYL5 , XEDGE1 , YEDGE1, XEDGE2 , YEDGE2,DEPTHXEDGE1 , YEDGE1:圆面或圆柱体直径一个端点X,Y坐标值XEDGE2 , YEDGE2 :圆面或圆柱体直径另一个端点X,Y坐标值DEPTH : Z方向的深度 4.3实体模型的修改1. ”KM0DIF ”命令功能：修改已存在的关键点格式：KMODIF,N PT,X,Y,ZNPT:关键点的编号，也可以为ALLX,Y,Z :新的坐标2. “ KDIST ”命令计算并输出两关键点间的距离格式：KDISTKP1 , KP2“

16、NDELE ”命令功能：删除所选的节点格式：NDELE,NODE1 ，NODE2，NINCNINC 从 N0DE1 至U NODE2。NODE1，NODE2，NINC:确定要删除节点的范围，按增量“ EDELE ”命令 功能：删除所选的单元格式：EDELE，IEL1，IEL2，NINC“ UPGEOM ”命令并更新有限元模型的几何模型（不改功能：将分析所得到的位移加到有限元模型的节点上， 变实体模型几何体）格式：UP GEOM,FACTOR,LSTE P, SBSTE P,Fn ame,Ext,-TP)MP”命令 功能：指定一个与温度相关的线性材料性能或常数（非线性材料用 格式：MP, Lab

17、,MAT,CO,C1,C2,C3,C4Lab,MAT（参考ANSYS参数化编程与命令手册P178）C0, C1,:多项式的常数项和系数 注意：一般单载荷步求解用线性材料性能，多载荷步用非线性材料属性。（载荷步就是为获取解而设置的加载条件）若为2次或更高次的多项式中的常数项和系数可以通过对离散的温度点进行线性插值求得。 性能指定适宜的温度步时，必须使用命令“MPTEMP ”，或MPTGEN “ MPTEMP ”命令功能：为材料属性定一个温度表格式：MPTEMP,STLOC,T1 , T2, T3, T4, T5, T6STLOC:确定输入温度的起始位置，默认值为最后填充值加1T1 , T2, T

18、3, T4, T5, T6 :指定从STLOC开始将要赋给 6个位置的温度值。442设置实常数1.”R”命令功能：定义单元的实常数格式：R,NSET,R1，R2，R3, R4，R5，R6NSET:设置实常数的编号R1，R2，R3，R4，R5，R6 :实常数的值或为表格边界条件输入的表格名称（tabname%），如果超过6个，可使用命令“ RMORE ”4.5单元设置与网格划分ANSYS软件平台提供了网格映射划分和自由适应划分的策略。映射划分用于曲线、曲面、实体的网格划 分方法，可使用三角形、四边形、四面体、五面体和六面体，通过指定单元边长、网格数量等参数对网格 进行严格控制，映射划分只用于规则

19、的几何图素，要求对应的线或面的单元数要相同。1、网格疏密网格疏密是指在结构不同部位采用大小不同的网格，这是为了适应计算数据的分布特点。在计算数据变化梯度较大的部位（如应力集中处），为了较好地反映数据变化规律，需要采用比较密集的网格。划分疏密不同的网格主要用于应力分析（包括静应力和动应力）。在结构温度场计算中也趋于采用均匀网格。2、单元阶次许多单元都具有线性、二次和三次等形式，其中二次和三次形式的单元称为高阶单元。选用高阶单元可提 高计算精度，因为高阶单元的曲线或曲面边界能够更好地逼近结构的曲线和曲面边界，且高次插值函数可 更高精度地逼近复杂场函数，所以当结构形状不规则、应力分布或变形很复杂时可

20、以选用高阶单元。但高 阶单元的节点数较多，在网格数量相同的情况下由高阶单元组成的模型规模要大得多，因此在使用时应权 衡考虑计算精度和时间。在有应力集中和刚度突变的地方，应该采用高阶单元来对其进行网格划分。增加网格数量和单元阶次都可以提高计算精度。因此在精度一定的情况下，用高阶单元离散结构时应选择适当的网格数量，太多的网 格并不能明显提高计算精度，反而会使计算时间大大增加。为了兼顾计算精度和计算量，同一结构可以采 用不同阶次的单元，即精度要求高的重要部位用高阶单元，精度要求低的次要部位用低阶单元。不同阶次 单元之间或采用特殊的过渡单元连接，或采用多点约束等式连接。平面单元类型：Plane42 :

21、 2维实体。该元素即可用于平面单元（平面应力或平面应变）也可用于轴对称单元。该元素由 个节点定义，每个节点 2个自由度：x,y方向。具有塑性，徐变，膨胀，应力强化，大变形，大应变能力。 PLANE82 :二维8节点实体。该元素是 plane42的高次形式。它为混合（四边形三角形）自动网格划分 提供了更精确的求解结果，并能承受不规则形状而不会产生任何精度上的损失。8节点元素具有位移协调形状，适用于模拟弯曲边界。该元素由8个节点定义，每个节点 2个自由度，x,y方向。可用于平面单元也可用于轴对称单元。具有塑性，徐变，膨胀，应力强化，大变形，大应变能力。并提供不同的输出选项。Plane83 :二维8

22、节点实体。用于承受非轴对称荷载的2维轴对称结构。如弯曲，剪切或扭转。该元素每个节点3个自由度：x,y,z方向。对于非扭转节点，这 3个方向分别代表半径，轴向和切线方向。该元素是 plane25的高次形式。它为混合（四边形-三角形）自动网格划分提供了更精确的求解结果，并能承受不 规则形状而不会产生任何精度上的损失。该元素也是plane82的一般轴向形式，其荷载不需要对称。Plane183 : 2维8节点实体。具有二次位移，适用于模拟不规则网格。该元素由8个节点定义，每个节点2个自由度，x,y方向。可用于平面单元也可用于轴对称单元。具有塑性，超弹性，应力强化，大变形，大 应变能力。可用来模拟几乎不

23、能压缩的次弹性材料和完全不能压缩的超弹性材料的变形。支持初始应力。并提供不同的输出选项。加载与求解5.1指定分析类型5.1.1指定分析类型与重启动格式：ANTYPE,Antype,Status,LDSTEP,SUBSTEP,Action Antype:默认值为：静态（STATIC）STATICH或0:静态分析，适用于所有自由度。BUCKLE或1:稳定性分析，仅对结构自由度有效MODAL或2:模态分析Status:指定分析状态：NEW:指定一次新的分析（默认）5.2.2施加载荷对结构施加载荷时，可以将载荷施加在实体模型或有限元模型上，但在系统求解时，载荷将被自动转换到 有限元模型上。载荷施加在实

24、体上的优缺点:1）施加在实体模型上的载荷与有限元网格无关，当网格被改变时，不会影响所施加的载荷。2）不能显示所有实体模型载荷。载荷施加在有限元模型上的优缺点：1)2）可直接应用于缩减分析中可以直接将约束施加在所选的节点上3）命令:单元的修改会使载荷失效，要求用户删除在修改前的载荷，并在修改后重新施加载荷功能：DA ”:在所选的面上施加自由度约束格式：DA,AREA,Lab,Value1, Value2AREA :要施加约束的面号Lab:有效的标签名，自由度标签:Value1:自由度值（实数部分）UX、UY、UZ、ROTX、ROTY、ROTZ、TEMPValue2:为FLOTRAN自由度，当值为

25、0时,工作平面和坐标系转换“ CSWPLA ”命令功能：在工作平面的原点处定义一个局部坐标系格式：CSWPLA,KCN,KCS,PAR1 ，PAR2KCN:赋给该坐标系的参考编号，必须大于10,KCS :坐标系统的类型。若为 0或CART,为直角坐标系；若为 系；3或TORO为环形坐标系1或CYLIN为柱坐标系；2或SPHE球坐标PAR1 :适应于椭圆或类似球体或环形系统，如果KCS=1或2,值，默认为1 ；PAR1是椭圆的Y轴半径与X轴半径的比PAR2 :仅适应于类似球体系统。如果KCS=2，PAR2是椭圆的Z轴半径与X轴半径的比值,默认为 1 ；2 :球坐标系;11或更大的数：由用户已定义

26、的局部坐标系CSYS”功能：激活一个预先已经定义的坐标系 格式：CSYSKCNKCN:将要激活的坐标系的参考号。若为 0:直角坐标系；1 :用Z轴为旋转轴的柱坐标系; 4或WP :工作平面；5:用Y轴为旋转轴的柱坐标系；若为ANSY諏值函数有关实体状态的取值函数NSEL(N)ESEL(E)KSEL(K)LSEL(L)ASEL(A)VSEL(V)表示某个实体状态，其返回值-1，没有选中，0，没有定义，1，被选中有关下一个被选实体的取值函数NDNEXT(N)ELNEXT(E)KP NEXT(K)LSNEXT(L)ARNEXT(A)VLNEXT(V)表示编号大于N,E,K,L,A,V的下一个被选实体

27、有关实体位置的取值函数CENTRX(E)CENTRY(E)CENTRZ(E)单元E在中心位置的X,Y,Z的坐标系（直角坐标系），有所选的节点决定NX(N)NY(N)NZ(N)KX(K)KY(K)KZ(K)节点N或关键点K在激活坐标系中X,Y,Z的坐标值LX(L,LFRAC)LY(L,LFRAC)LZ(L,LFRAC)线段L在长度比率为LFRAC(01)时的X,Y,Z的坐标值有关最靠近某位置的节点或关键点编号的取值函数NODE(X,Y,Z)KP (X,Y,Z)被选择的节点嘴靠近 X,Y,Z位置的节点或关键点编号(在激活的坐标系下，如果存在多个节点或关键点,那么取其最小值)有关距离的取值函数DIS

28、TND(N1,N2)DISTK P(K1,K2)节点或关键点两点之间的距离DISTEN(E,N)单元E的中心点与节点N之间的距离，中心点将由单元上被选择的节点确定有关角度的取值函数ANGLEN(N1,N2,N3)ANGLEK(K1,K2,K3)3个节点中，N1或K1是顶点节点或关键点两条边之间的夹角，缺省时单位为弧度，其中所选择的有关最靠近实体的节点，关键点和单元的取值函数NNEAR(N)最靠近节点N的被选节点KNEAR(K)最靠近关键点K的被选关键点ENEARN(N)最靠近节点N的被选单元，单元的位置将由被选节点确定有关面积的取值函数AREAND(N1,N2,N3)AREAK P(K1,K2

29、,K3)由节点或关键点3点为顶点围成的三角形的面积ARNODE(N)从选择且与节点N的相关的单元分配给节点 N的面积，对于二维平面，返回与节点相关的边面积，对于轴对称体，返回与节点相关的边缘表面积，对于三维实体，返回与节点相关的端面积有关节点和关键点的方向余弦NORMNX(N1,N2,N3)NORMNY(N1,N2,N3)NORMNZ(N1,N2,N3)NORMKX(K1,K2,K3)NORMKY(K1,K2,K3)NORMKZ(K1,K2,K3)由节点或关键点三点确定的平面法线在X或Y或Z方向的方向余弦有关节点与单元连接性的取值函数ENEXTN(N,LOC)与节点N相连接的单元。当多个单元共享同一个节点时，在结果列表中,LOC表示为位置，并在列表的结束处返回一个零NELEM(E, NPOS)单元E中，NPOS位置处的节点编号与单元面相关的取值函数ELADJ(E,FA