您的位置:网站首页 > Ansys教程

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

时间:2010-11-14 11:34:33 来源:未知

3 /solu
u /solu 进入求解器
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 以同样工作名
#p#分页标题#e#(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:程序执行时间(秒)限制,缺省为无穷
#p#分页标题#e#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 在指定节点加集中荷载
#p#分页标题#e#node:节点号
lab: Fx,Fy,Fz,Mx,My,Mz
value: 力大小
value2: 力的第二个大小(如果有复数荷载)
nend,ninc:在从node到nend的节点(增量为ninc)上施加同样的力
注意:(1)节点力在节点坐标系中定义,其正负与节点坐标轴正向一致
u sfa, area, lkey, 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 显示变形前和变形后的结构
#p#分页标题#e#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、BEAM189截面结果
说明:只有刚计算完还未退出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 不显示未变形的结构
#p#分页标题#e#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 不显示符号
#p#分页标题#e#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) 等效应力
等效应力或若则有 (弹性状态下)
#p#分页标题#e#