ansys命令流----前后处理和求解常用命令之求解与后处理

3.1 加边界条件

u 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以同样数值施加给所选节点。

注意：在节点坐标系中讨论

3.2 设置求解选项

u 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

说明：继续以前的分析（因某种原因中断）有两种类型

singleframe restart: 从停止点继续

需要文件：jobname.db 必须在初始求解后马上存盘

jobname.emat 单元矩阵

jobname.esav 或 .osav : 如果.esav坏了，将.osav改为.esav

results file: 不必要，但如果有，后继分析的结果也将很好地附加到它后面

注意：如果初始分析生成了.rdb, .ldhi, 或rnnn 文件。必须删除再做后继分析

步骤： （1）进入anasys 以同样工作名

（2）进入求解器，并恢复数据库

（3）antype, rest

（4）指定附加的荷载

（5）指定是否使用现有的矩阵（jobname.trl）（缺省重新生成）

kuse: 1 用现有矩阵

（6）求解

multiframe restart:从以有结果的任一步继续（用不着）

u pred,sskey, --,lskey….. 在非线性分析中是否打开预测器

sskey: off 不作预测（当有旋转自由度时或使用solid65时缺省为off）

on 第一个子步后作预测（除非有旋转自由度时或使用solid65时缺省为on）

-- ： 未使用变量区

lskey: off 跨越荷载步时不作预测（缺省）

on 跨越荷载步时作预测（此时sskey必须同时on）

注意：此命令的缺省值假定solcontrol为on

u autots, key 是否使用自动时间步长

key:on: 当solcontrol为on时缺省为on

off: 当solcontrol为off时缺省为off

1: 由程序选择（当solcontrol为on且不发生autots命令时在 .log文件中纪录“1”

注意：当使用自动时间步长时，也会使用步长预测器和二分步长

u NROPT, option,--,adptky 指定牛顿拉夫逊法求解的选项

OPTION: AUTO:程序选择

FULL:完全牛顿拉夫逊法

MODI:修正的牛顿拉夫逊法

INIT：使用初始刚阵

UNSYM：完全牛顿拉夫逊法，且允许非对称刚阵

ADPTKY:ON: 使用自适应下降因子

OFF：不使用自适应下降因子

u NLGEOM，KEY

KEY: OFF:不包

括几何非线性（缺省）

ON：包括几何非线性

u ncnv, kstop, dlim, itlim, etlim, cplim 终止分析选项

kstop: 0 如果求解不收敛，也不终止分析

1 如果求解不收敛，终止分析和程序（缺省）

2如果求解不收敛，终止分析,但不终止程序

dlim：最大位移，缺省为1.0e6

itlim: 累积迭代次数，缺省为无穷多

etlim：程序执行时间（秒），缺省为无穷

cplim：cpu时间（秒），缺省为无穷

u 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:对某些单元包括应力荷载刚化（查）

vtol：

u outres, item, freq, cname 规定写入数据库的求解信息

item: all 所有求解项

basic 只写nsol, rsol, nload, strs

nsol 节点自由度

rsol 节点作用荷载

nload 节点荷载和输入的应变荷载（？）

strs 节点应力

freq: 如果为n,则每n步（包括最后一步）写入一次

none: 则在此荷载步中不写次项

all: 每一步都写

last: 只写最后一步（静力或瞬态时为缺省）

3.3 定义载荷步

u nsubst, nsbstp, nsbmx, nsbmn, carry 指定此荷载步的子步数

nsbstp: 此荷载步的子步数

如果自动时间步长使用autots,则此数定义第一子步的长度；如果solcontrol打开，且3D面-面接触单元使用，则缺省为1-20步；如果solcontrol打开，并无3D接触单元，则缺省为1子步；如果solcontrol关闭，则缺省为以前指定值；如以前未指定，则缺省为1）

nsbmx, nsbmn：最多，最少子步数（如果自动时间步长打开）？

u time, time 指定荷载步结束时间

注意：第一步结束时间不可为“0”

u f, node, lab, value, value2, nend, ninc 在指定节点加集中荷载

node:节点号

lab: Fx,Fy,Fz,Mx,My,Mz

value: 力大小

value2: 力的第二个大小（如果有复数荷载）

nend,ninc：在从node到nend的节点（增量为ninc）上施加同样的力

注意：（1）节点力在节点坐标系中定义，其正负与节点坐标轴正向一致

u sfa, area, lke

y, lab, value, value2 在指定面上加荷载

area: n 面号

all 所有选中号

lkey: 如果是体的面，忽略此项

lab: pres

value: 压力值

u 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”

3.3.1 注意

1. 尽量加面载，不加集中力，以免奇异点

2. 面的切向荷载必须借助面单元

3.4 求解载荷步

u lssolve, lsmin, lsmax, lsinc 读入并求解多个荷载步

lsmin, lsmax, lsinc ：荷载步文件范围

4 /post1（通用后处理）

u set, lstep, sbstep, fact, king, time, angle, nset 设定从结果文件读入的数据

lstep :荷载步数

sbstep：子步数，缺省为最后一步

time：时间点（如果弧长法则不用）

nset： data set number

u dscale, wn, dmult 显示变形比例

wn: 窗口号（或all）,缺省为1

dmult, 0或auto : 自动将最大变形图画为构件长的5%

u pldisp, kund 显示变形的结构

kund： 0 仅显示变形后的结构

1 显示变形前和变形后的结构

2 显示变形结构和未变形结构的边缘

u *get, par, node, n, u, x(y,z) 获得节点n的x(y,z)位移给参数par

等价于函数 ux(n),uy(n),uz(z)

node(x,y,z): 获得(x,y,z)节点号

arnode(x,y,z)：获得和节点n相连的面

注意：此命令也可用于/solu模块

u fsum, lab, item 对单元之节点力和力矩求和

lab: 空 在整体迪卡尔坐标系下求和

rsys 在当前激活的rsys坐标系下求和

item: 空 对所有选中单元（不包括接触元）求和

cont: 仅对接触节点求和

u PRSSOL, ITEM, COMP 打印BEAM188、BEAM1截面结果

说明：只有刚计算完还未退出ANSYS时可用，重新进入ANSYS时不可用

item comp 截面数据及分量标志

S COMP X,XZ,YZ应力分量

PRIN S1,S2,S3主应力SINT应力强度，SEQV等效应力

EPTO COMP 总应变

PRIN 总主应变，应变强度，等效应变

EPPL COMP 塑性应变分量

PRIN 主塑性应变，塑性应变强度，等效塑性应变

u plnsol, item, comp, kund, fact 画节点结果为连续的轮廓线

item: 项目（见下表）

comp: 分量

kund: 0 不显示未变形的结构

1 变形和未变形重叠

2 变形轮廓和未变形边缘

fact: 对于接触的2D显示的比例系数，缺省为1

item comp discription

u x,y,z,sum 位移

rot x,y,z,sum 转角

s x,y,z,xy,yz,xz 应力分量

1，2，3 主应力

Int,eqv 应力intensity,等效应力

epeo x,y,z,xy,yz,xz 总位移分量

1,2,3 主应变

Int,eqv 应变intensity,等效应变

epel x,y,z,xy,yz,xz 弹性应变分量

1，2，3 弹性主应变

Int,eqv 弹性intensity,弹性等效应变

eppl x,y,z,xy,yz,xz 塑性应变分量

u PRNSOL, item, comp 打印选中节点结果

item: 项目（见上表）

comp: 分量

u PRETAB, LAB1, LAB2, ……LAB9 沿线单元长度方向绘单元表数据

LABn : 空： 所有ETABLE命令指定的列名

列名：任何ETABLE命令指定的列名

u PLLS, LABI, LABJ, FACT, KUND 沿线单元长度方向绘单元表数据

LABI:节点I的单元表列名

LABJ:节点J的单元表列名

FACT: 显示比例，缺省为1

kund: 0 不显示未变形的结构

1 变形和未变形重叠

2 变形轮廓和未变形边缘

5 /post26 （时间历程后处理）

u nsol, nvar, node, item, comp,name

在时间历程后处理器中定义节点变量的序号

nvar：变量号（从2到nv（根据numvar定义））

node: 节点号

item comp

u x, y,z

rot x, y,z

u ESOL, NVAR, ELEM, NODE, ITEM, COMP, NAME 将结果存入变量

NVAR: 变量号，2以上

ELEM: 单元号

NODE: 该单元的节点号，决定存储该单元的哪个量，如果空，则给出平均值

ITEM:

COMP:

NAME: 8字符的变量名， 缺省为ITEM加COMP

u rforce, nvar, node, item, comp, name 指定待存储的节点力数据

nvar: 变量号

node: 节点号

item comp

F x, y.z

M x, y,z

name: 给此变量一个名称，8个字符

u add, ir, ia,ib,ic,name,--,--,facta, factb, factc

将ia,ib,ic变量相加赋给ir变量

ir, ia,ib,ic：变量号

name: 变量的名称

u /grid, key

key: “0” 或“off” 无网络

“1”或“on” xy网络

“2”或“x” 只有x线

“3”或“y” 只有y线

u xvar, n

n: “0”或“1” 将x轴作为时间轴

“n” 将x轴表示变量“n”

“-1” ？

u /axlab, axis, lab 定义轴线的标志

axis: “x”或“y”

lab: 标志，可长达30个字符

u plvar, nvar, nvar2, ……,nvar10 画出要显示的变量（作为纵坐标）

u prvar, nvar1, ……,nvar6 列出要显示的变量

6 PLOTCONTROL菜单命令

u pbc, ilem, ……,key, min, max, abs 在显示屏上显示符号及数值

item: u 所加的位移约束

rot 所加的转角约束

key: 0 不显示符号

1 显示符号

2 显示符号及数值

u /SHOW, FNAME, EXT, VECT, NCPL 确定图形显示的设备及其他参数

FNAME: X11:屏幕

文件名：各图形将生成一系列图形文件

JPEG: 各图形将生成一系列JPEG图形文件

说明：没必要用此命令，需要的图形文件可计算后再输出

7 参数化设计语言

u *do, par, ival, fval, inc 定义一个do循环的开始

par: 循环控制变量

ival, fval, inc：起始值，终值，步长（正，负）

u *enddo 定义一个do循环的结束

u *if,val1, oper, val2, base: 条件语句

val1, val2: 待比较的值（也可是字符，用引号括起来）

oper: 逻辑操作（当实数比较时，误差为1e-10）

eq, ne, lt, gt, le, ge, ablt, abgt

base: 当oper结果为逻辑真时的行为

lable: 用户定义的行标志

stop: 将跳出anasys

exit

: 跳出当前的do循环

cycle: 跳至当前do循环的末尾

then: 构成if-then-else结构

注意：不允许跳出、跳进一个do,if循环至label句 ？

8 理论手册

1.方程组解法：(1)直接解法；(2)迭代解法

(1) 直接解法：a.稀疏矩阵法；b. 波前解法

a. 稀疏矩阵法：占内存大，但运算次数少；通过变换刚度矩阵的顺序使得非零元素最少

b. 波前解法： 占内存小

波前是指在还没有一个单元被解完的时候激活的方程数？

(2) 迭代解法：JCG法；PCG法；ICCG法

JCG法：可解实数、对称、非对称矩阵

PCG法：高效求解各种矩阵（包括病态），但仅解实、对称矩阵

ICCG法：类似JCG,但更强

2. 应变密度，等效应变，应力密度，等效应力

（1）应变密度（strain intensity）

应变密度

是三个主应变

（2）等效应变

有效泊松比 ：用户由avprin 命令设定；0（如果不设定）

（3）应力密度(stress intensity)

应力密度

(4) 等效应力

等效应力或若则有 （弹性状态下）