结合自身经验，谈ANSYS中的APDL命令(一)
2014-12-21 11:07阅读：
结合自身经验，谈ANSYS中的APDL命令(一)
发表时间：2009-4-7
作者: 倪欣
来源:
e-works
关键字:
ansys APDL 命令流
在ANSYS中，命令流是由一条条ANSYS的命令组成的一个命令组合，这些命令按照一定顺序排布，能够完成一定的ANSYS功能，本文是作者结合自身经验所总结的一些命令。
在ANSYS中，命令流是由一条条ANSYS的命令组成的一个命令组合，这些命令按照一定顺序排布，能够完成一定的ANSYS功能，这些功能一般来说通过菜单操作也能够实现（而那些命令流能够实现，菜单操作实现不了的单个命令比较少见）。
以下命令是结合我自身经验，和前辈们的一些经验而总结出来的，希望对大家有帮助。
（1）．Lsel, type, item, comp, vmin, vmax, vinc, kswp 选择线
type: s
从全部线中选一组线
r
从当前选中线中选一组线
a
再选一部线附加给当前选中组
au
none
u(unselect)
inve: 反向选择
item: line 线号
loc 坐标
length 线长
comp: x,y,z
kswp: 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，C
TBDATA，2，ψ
TBDATA，3，……
如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：
MP，EX，1，1E8
MP，NUXY，1，0.3
TB，DP，1
TBDATA，1，27
TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,deg
（6）. 根据需要耦合某些节点自由度
cp, nset, lab,,node1,node2,……node17
nset: 耦合组编号
lab: ux,uy,uz,rotx,roty,rotz ,all
node1-node17:
待耦合的节点号。如果某一节点号为负，则此节点从该耦合组中删去。如果node1=all,则所有选中节点加入该耦合组。
注意：1，不同自由度类型将生成不同编号
2，不可将同一自由度用于多套耦合组
cpintf, lab, toler 将相邻节点的指定自由度定义为耦合自由度
lab: ux,uy,uz,rotx,roty,rotz ,all
toler: 公差，缺省为0.0001
说明：先选中欲耦合节点，再执行此命令
（7）. D, node, lab, value, value2, nend, ninc, lab2, lab3, ……lab6
定义节点位移约束
Node : 预加位移约束的节点号，如果为all,则所有选中节点全加约束，此时忽略nend和ninc.
Lab: ux,uy,uz,rotx,roty,rotz,all
Value,value2: 自由度的数值（缺省为0）
Nend, ninc: 节点范围为：node-nend，编号间隔为ninc
Lab2-lab6: 将lab2-lab6以同样数值施加给所选节点。
（8）. 设置求解选项
antype, status, ldstep, substep, action
antype: static or 1 静力分析
buckle or 2 屈曲分析
modal or 3 模态分析
trans or 4 瞬态分析
status: new 重新分析（缺省），以后各项将忽略
rest 再分析，仅对static,full transion 有效
ldstep: 指定从哪个荷载步开始继续分析，缺省为最大的，runn数（指分析点的最后一步）
substep: 指定从哪个子步开始继续分析。缺省为本目录中，runn文件中最高的子步数
action, continue: 继续分析指定的ldstep,substep
（9）. pred,sskey, --,lskey….. 在非线性分析中是否打开预测器
sskey: off 不作预测（当有旋转自由度时或使用solid65时缺省为off）
on
第一个子步后作预测（除非有旋转自由度时或使用solid65时缺省为on）
-- ： 未使用变量区
lskey: off 跨越荷载步时不作预测（缺省）
on 跨越荷载步时作预测（此时sskey必须同时on）
注意：此命令的缺省值假定solcontrol为on分页
（10）. autots, key 是否使用自动时间步长
key:on: 当solcontrol为on时缺省为on
off: 当solcontrol为off时缺省为off
1: 由程序选择（当solcontrol为on且不发生autots命令时在 .log文件中纪录“1”
注意：当使用自动时间步长时，也会使用步长预测器和二分步长
（11）. NROPT, option,--,adptky 指定牛顿拉夫逊法求解的选项
OPTION: AUT程序选择
FULL:完全牛顿拉夫逊法
MODI:修正的牛顿拉夫逊法
INIT：使用初始刚阵
UNSYM：完全牛顿拉夫逊法，且允许非对称刚阵
ADPTKY:ON: 使用自适应下降因子
OFF：不使用自适应下降因子
（12）. NLGEOM，KEY
KEY: OFF:不包括几何非线性（缺省）
ON：包括几何非线性
（13）. ncnv, kstop, dlim, itlim, etlim, cplim 终止分析选项
kstop: 0 如果求解不收敛，也不终止分析
1 如果求解不收敛，终止分析和程序（缺省）
2 如果求解不收敛，终止分析,但不终止程序
dlim：最大位移限制，缺省为1.0e6
itlim: 累积迭代次数限制，缺省为无穷多
etlim：程序执行时间（秒）限制，缺省为无穷
cplim：cpu时间（秒）限制，缺省为无穷
（14）. solcontrol ,key1, key2,key3,vtol 指定是否使用一些非线性求解缺省值
key1: on 激活一些优化缺省值（缺省）
CNVTOL
Toler=0.5%Minref=0.01(对力和弯矩)
NEQIT
最大迭代次数根据模型设定在15~26之间
ARCLEN
如用弧长法则用较ansys5.3更先进的方法
PRED
除非有rotx,y,z或solid65，否则打开
LNSRCH
当有接触时自动打开
CUTCONTROL
Plslimit=15%, npoint=13
SSTIF
当NLGEOM,on时则打开
NROPT,adaptkey
关闭（除非：摩擦接触存在；单元12,26,48,49,52存在；当塑性存在且有单元20,23,24,60存在）
AUTOS 由程序选择
off 不使用这些缺省值
key2: on 检查接触状态（此时key1为on）
此时时间步会以单元的接触状态（据keyopt(7)的假定）为基础
当keyopt(2)=on 时，保证时间步足够小
key3: 应力荷载刚化控制，尽量使用缺省值
空：缺省，对某些单元包括应力荷载刚化，对某些不包括
nopl:对任何单元不包括应力刚化
incp:对某些单元包括应力荷载刚化
（15）. outres, item, freq, cname 规定写入数据库的求解信息
item: all 所有求解项
basic 只写nsol, rsol, nload, strs
nsol 节点自由度
rsol 节点作用荷载
nload 节点荷载和输入的应变荷载（？）
strs 节点应力
freq: 如果为n,则每n步（包括最后一步）写入一次
none: 则在此荷载步中不写次项
all: 每一步都写
last: 只写最后一步（静力或瞬态时为缺省）
（16）. nsubst, nsbstp, nsbmx, nsbmn, carry 指定此荷载步的子步数
nsbstp: 此荷载步的子步数
如果自动时间步长使用autots,则此数定义第一子步的长度；如果solcontrol打开，且3D面-面接触单元使用，则缺省为1-20步；如果solcontrol打开，并无3D接触单元，则缺省为1子步；如果solcontrol关闭，则缺省为以前指定值；如以前未指定，则缺省为1）
nsbmx, nsbmn：最多，最少子步数（如果自动时间步长打开）
（17）. time, time 指定荷载步结束时间
注意：第一步结束时间不可为“0”
（18）. f, node, lab, value, value2, nend, ninc 在指定节点加集中荷载
node:节点号
lab: Fx,Fy,Fz,Mx,My,Mz
value: 力大小
value2: 力的第二个大小（如果有复数荷载）
nend,ninc：在从node到nend的节点（增量为ninc）上施加同样的力
注意：（1）节点力在节点坐标系中定义，其正负与节点坐标轴正向一致
（19）. sfa, area, lkey, lab, value, value2 在指定面上加荷载
area: n 面号
all 所有选中号
lkey: 如果是体的面，忽略此项
lab: pres
value: 压力值
（20）. SFBEAM, ELEM, LKEY, LAB, VALI, VALJ, VAL2I, VAL2J, IOFFST,
JOFFST
对梁单元施加线荷载
ELEM: 单元号，可以为ALL,即选中单元
LKEY: 面载类型号，见单元介绍。对于BEAM188，1为竖向；2为横向；3为切向
VALI,VALJ: I, J节点处压力值
VAL2I,VAL2J: 暂时无用
IOFFST, JOFFST: 线载距离I, J 节点距离
u lswrite, lsnum 将荷载与荷载选项写入荷载文件中
lsnum ：荷载步文件名的后缀，即荷载步数
当 stat 列示当前步数
init 重设为“1”
缺省为当前步数加“1”
注意：
1. 尽量加面载，不加集中力，以免奇异点
2. 面的切向荷载必须借助面单元
（21）. lssolve, lsmin, lsmax, lsinc 读入并求解多个荷载步
lsmin, lsmax, lsinc ：荷载步文件范围
