ANSYS非线性分析命令解释
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
ANSYS应用基于问题物理特性的自动求解控制方法,把各种非线性分析控制参数设置到合适的值。
如果用户对这些设置不满意,还可以手工设置。
下列命令的缺省设置已进行了优化处理:
AUTOTS PRED MONITOR
DELTIM NROPT NEQIT
NSUBST TINTP SSTIF
CNVTOL CUTCONTROL KBC
LNSRCH OPNCONTROL EQSLV
ARCLEN CDWRITE LSWRITE
这些命令及其设置在将在后面讨论。
参见《ANSYS Commands Reference》。
如果用户选择自己的设置而不是ANSYS的缺省设置,或希望用以前版本的ANSYS的输入列表,则可用/ SOLU 模块的SOLCONTROL ,OFF命令,或在/ BATCH 命令后用/ CONFIG ,NLCONTROL,OFF命令。
参见SOLCONTROL 命令的详细描述。
ANSYS对下面的分析激活自动求解控制
单场的非线性或瞬态结构以及固体力学分析,在求解自由度为UX、UY、UZ、ROTX、ROTY、ROTZ的结合
时;
单场的非线性或瞬态热分析,在求解自由度为TEMP时;
注意--本章后面讨论的求解控制对话框,不能对热分析做设置。
用户必须应用标准的ANSYS求解命令或GUI来设置。
2.2 非线性静态分析步骤
尽管非线性分析比线性分析变得更加复杂,但处理基本相同。
只是在非线形分析的过程中,添加了需要的非线形特性。
非线性静态分析是静态分析的一种特殊形式。
如同任何静态分析,处理流程主要由以下主要步骤组成:
建模;
设置求解控制;
设置附加求解控制;
加载;
求解;
考察结果。
2.2.1 建模
这一步对线性和非线性分析基本上是一样的,尽管非线性分析在这一步中可能包括特殊的单元或非线性材料性
质,参考§4《材料非线性分析》,和§6.1《单元非线性》。
如果模型中包含大应变效应,应力─应变数据必
须依据真实应力和真实(或对数)应变表示。
参见《ANSYS Modeling and Meshing Guide》。
在ANSYS中建立了模型后,应该设置求解控制(分析类型、分析选项、荷载步等)选项,施加荷载,最后求解。
非线性分析与线性分析的不同之处是,前者需要许多荷载增量,并且总是需要平衡迭代。
下面讨论一般过程。
参见本章的例子。
2.2.2 设置求解控制
设置求解控制包括定义分析类型、设置分析的常用选项和指定荷载步选项。
在做结构非线性静态分析时,可以应用求解控制对话框来设置。
该对话框对许多非线性静态分析提供了缺省设置。
这样,用户需要的设置降低到最少。
求解控制框的缺省设置,基本上与§2.1所述的自动求解控制的设置相同。
由于求解控制对话框是非线性静态分析的推荐工具,我们在下面将详细论述,如用户不想用这个对话框(GUI:Main Menu>Solution>-Analysis Type-Sol"n Control),可以应用标准的ANSYS求解命令集或相应的菜单(GUI:Main Menu> Solution>Unabridged Menu>option)。
求解控制对话框的概况,见《ANSYS Basic Analysis Guide》§3.11。
注意--对于非线性结构完全瞬态分析,建议应用求解控制对话框,但并不是必须如此,见§2.3。
2.2.2.1 求解控制对话框—进入
2.2.2.8.1 自动时间步
ANSYS的自动求解控制打开自动时间步长[ AUTOTS ,ON]。
这一选项允许程序确定子步间载荷增量的大小和决定在求解期间是增加还是减小时间步(子步)长。
在一个时间步的求解完成后,下一个时间步长的大小基于四种因素预计:
在最近过去的时间步中使用的平衡迭代的数目(更多次的迭代成为时间步长减小的原因);
对非线性单元状态改变预测(当状态改变临近时减小时间步长);
塑性应变增加的大小;
蠕变增加的大小。
2.2.2.8.2 收敛准则
程序将连续进行平衡迭代直到满足收敛准则[ CNVTOL ](或者直到达到允许的平衡迭代的最大次数〔NEQIT 〕。
如果缺省的收敛准则不满意,可以自己定义收敛准则。
ANSYS的自动求解控制应用等于0.5%的力(或力矩)的L2-范数容限(TOLER),这对于大部分情况合适。
在大多数情况下,除了进行力范数的检查外,还进行TOLER等于5%的位移L2-范数的检查。
缺省时,程序将通过比较不平衡力的平方和的平方根(SRSS)与VALUE×TOLER的值来检查力(在包括转动自由度时,还有力矩)的收敛。
VALUE的缺省值是所加载荷(或在施加位移时,Netwton-Raphson回复力)的SRSS,或MINREF (其缺省为0.001),取较大者。
如果SOLCONTROL ,OFF,则对于力的收敛,TOLER的缺省值是0.001,而MINREF的缺省为1.0。
用户应当几乎总是使用力收敛检查。
可以添加位移(或者转动)收敛检查。
对于位移,程序将收敛检查建立在当前(i)和前面(i-1)次迭代之间的位移改变(Δu)上,Δu =u i -u i-1 。
注意─如果用户明确地定义了任何收敛准则[ CNVTOL ],缺省准则将失效。
因此,如果用户定义了位移收敛检查,用户将不得不再定义力收敛检查(使用多个CNVTOL 命令来定义多个收敛准则)。
使用严格的收敛准则将提高用户的结果的精度,但以更多次的平衡迭代为代价。
如果用户想紧缩(或放松-但不推荐)收敛准则,用户应当改变TOLER一到两个数量级。
一般地,用户应当继续使用VALUE的缺省值;也就是,通过调整TOLER,而不是VALUE,来改变收敛准则。
用户应当确保MINREF=0.001的缺省值在用户的分析范围内有意义。
如果应用某一单位系统,使荷载变得十分小,可能需要指定较小的MINREF值。
在非线性分析中,不推荐把两个或多个不相连的结构放在一起分析,因为收敛检查试图把这些彼此不相连的结构联系起来,通常会产生不希望的残余力。
在单一和多自由度系统中检查收敛
要在单自由度系统中检查收敛,用户对这一个自由度计算出不平衡力,然后将这个值与给定的收敛准则(VALUE×TOLER)比较(同样也可以对单自由度的位移或旋转收敛进行类似的检查)。
然而,在多自由度系统中,用户也许想使用不同的比较方法。
ANSYS程序提供三种不同的矢量范数用于收敛检查:
无穷范数在用户模型中的每一个自由度处重复单-自由度检查;
L1范数将收敛准则同所有自由度的不平衡力(或力矩)的绝对值的总和相比较;
L2范数使用所有自由度不平衡力(或力矩)的SRSS进行收敛检查。
当然,对于位移收敛检查,可以执行附加的L1、L2检查。
实例
对于下面例子,如果不平衡力(在每一个自由度处单独检查)小于或等于5000×0.0005(也就是2.5),且如果位移的改变(以SRSS检查)小于或等于10×0.001(也就是0.01),子步将认为是收敛的。
CNVTOL,F,5000,0.0005,0
CNVTOL,U,10,0.001,2
2.2.2.8.3 平衡迭代的最大次数
ANSYS的自动求解控制把NEQIT的值,根据问题的物理特性,设置为15到26次平衡迭代。
应用小时间步,可减少二次收敛迭代次数。
这个选项限制了一个子步中进行的最大平衡迭代次数(如关闭求解控制,缺省=25)。
如果在这个平衡迭代次数之内不能满足收敛准则,且如果自动步长是打开的[ AUTOTS ],分析将尝试使用二分法。
如果二分法是不可能的,那么,分析将或者终止,或者进行下一个载荷步,依据用户在NCNV 命令中发出的指示。
2.2.2.8.4 预测─修正选项
如不存在梁或壳单元,ANSYS的自动求解控制设置PRED ,ON。
如果当前子步的步长大大减小,PRED将关闭。
对于瞬态分析,将关闭预测选项。
对于每一个子步的第一次平衡迭代,用户可以激活自由度求解的预测。
这个特点将加速收敛,且如果非线性响应是相对平滑的,它特别的有用。
在包含大转动或粘弹的分析中它并不是非常有用。
在大转动分析中,预测可能引起发散,因而不推荐使用。
2.2.2.8.5 线性搜索选项
ANSYS的自动求解控制,将根据需要关闭或打开线性搜索。
对大多数接触问题,LNSRCH 打开。
对大多数非接触问题,LNSRCH 关闭。
这个收敛增强工具用程序计算出的比例因子(具有0和1之间的值)乘以计算出的位移增量。
因为线性搜索算法是用来对自适应下降选项[ NROPT ]进行的替代,如果线性搜索选项是开,自适应下降不被自动激活。
不建议用户同时激活线性搜索和自适应下降。
当存在强制位移时,只有至少有一次迭代的线性搜索值为1,计算才可以收敛。
ANSYS调节整个ΔU矢量,包括强制位移值,否则,除了强制自由度处以外,一个小的位移值将随处发生。
直到迭代中的某一次具有1的线性搜索值,ANSYS才施加全部位移值。
2.2.2.8.6 步长缩减准则
为了更好地控制时间步长上的二分和缩减,应用[ CUTCONTROL ,Lab ,VALUE,Option]。
缺省时,对于Lab=PLSLIMIT(最大塑性应变增量极限),VALUE设置为15%。
设这么大的值,是为避免由高塑性应变引起的不必要的二分,因为高塑性应变可能是由用户并不感兴趣的局部奇异引起。
对于显式蠕变(Option=0),Lab=CRPLIM(蠕变增量极限),VALUE设置为10%。
这对蠕变分析是一个合理的极限。
对于隐式蠕变(Option=1),缺省为无最大蠕变准则。
但是用户可以指定蠕变率控制。
对于二阶动力方程,每个周期的点数(Lab=NPOINT),缺省为VALUE=13,这样可以很小的代价获得有效精度。
2.2.3 设置附加求解选项
本节论述的选项,不出现在求解对话框中。
这些选项的缺省值,一般很少需要改变。
2.2.
3.1 求解控制对话框不能设置的高级分析选项
2.2.
3.1.1 应力刚化效应
为了考虑屈曲、分叉行为,ANSYS在所有几何非线性分析中,包括了应力刚化。
如果用户有信心放弃这种效应,则可以关闭应力刚化效应( SSTIF ,OFF)。
在一些单元中,这个命令无作用,见《ANSYS Elements Reference》。
命令:SSTIF
GUI:Main Menu>Solution>Unabridged Menu>Analysis Options
2.2.
3.1.2 牛顿-拉普森选项
在存在非线性时,ANSYS的自动求解控制将应用自适应下降关闭的完全牛顿-拉普森选项。
但在应用点-点,点-面接触单元的摩擦接触分析中,自适应下降功能是自动打开的(如CONTAC12、CONTAC48、CONTAC49、CONTAC52单元)。
下伏接触单元需要自适应下降才能收敛。
命令:NROPT
GUI:Main Menu>Solution>Unabridged Menu>Analysis Options
仅在非线性分析中使用这个选项。
这个选项指定在求解期间每隔多久修改一次正切矩阵。
如果用户不想采用缺省值,可以指定这些值中的一个:
·程序选择( NROPT ,ANTO):程序基于用户模型中存在的非线性种类选用这些选项中的一个。
需要时牛顿-拉普森方法将自动激活自适应下降。
·完全牛顿-拉普森法( NROPT ,FULL);程序使用完全的牛顿-拉普森方法。
在这种处理方法中,每进行一次平衡迭代,就修改刚度矩阵一次。
如果自适应下降是打开(可选),只要迭代保持稳定(也就是只要残余项减小,且没有负主对角线出现),程序将仅使用正切刚度阵。
如果在一次迭代中探测到发散倾向,程序抛弃发散的迭代且重新开始求解,应用正切和正割刚度矩阵的加权组合。
当迭代回到收敛模式时,程序将重新开始使用正切刚度矩阵。
对复杂的非线性问题自适应下降通常将提高程序获得收敛的能力,但它只支持《ANSYS Element Reference》中由单元输入汇总表中的―Special Features‖指明的单元(见《ANSYS Element Reference》表4.n.1,其中n为单元编号)。
·修正的牛顿-拉普森法( NROPT ,MODI):使用修正的牛顿-拉普森方法。
在这种方法中,正切刚度矩阵在每一子步中都被修正。
在一个子步的平衡迭代期间矩阵不被改变。
这个选项不适用于大变形分析。
自适应下降不可用。
·初始刚度牛顿-拉普森法( NROPT ,INIT):在每一次平衡迭代中都使用初始刚度矩阵。
这一选项比完全选项似乎较不易发散,
但它经常要求更多次的迭代来得到收敛。
它不适用于大变形分析。
自适应下降不可用。
·不对称矩阵完全牛顿-拉普森方法( NROPT ,UNSYM):应用完全牛顿-拉普森方法,刚度矩阵在每一次平衡迭代中都修正。
此外,它生成并使用在下面任何一种情况中可以应用的不对称矩阵:
如用户在运行压力产生的破坏分析,不对称的压力荷载刚度可能有助于取得收敛。
可应用SOLCONTROL ,INCP 命令来包括荷载刚度。
如果应用TB ,USER 命令定义不对称材料模型,则需要用NROPT ,UNSYM 命令来充分应用所定义的特性。
如进行接触分析,不对称接触刚度矩阵可以完全地耦合滑动和法向刚度。
见§5.4。
用户应首先试验NROPT ,FULL 命令;然后如果收敛困难的话,再试验NROPT ,UNSYM 命令。
注意,应用不对称求解器需要
比对称求解器更多的计算机时间。
·如果模型有多态单元,则将在状态改变时进行叠代修正,而不管牛顿-拉普森选项设置如何。
2.2.
3.2 求解控制对话框不能设置的高级荷载步选项
2.2.
3.2.1 蠕变准则
如果结构表现出蠕变行为,可以指定蠕变准则用于自动时间步调整[ CRPLIM , CRCR , Option ](如果自动时间步长[ AUTOTS ]关闭,蠕变准则无效)。
程序将对所有单元计算蠕应变增量(在最近时间步中蠕变的变化Δε cr )对弹性应变ε el 的比值。
如果最大比值比判
据CRCR大,程序将减小下一个时间步长;如果小,程序或许增加下一个时间步长(同样,程序将把自动时间步长建立在平衡迭代次数、即将发生的单元状态改变以及塑性应变增量的基础上。
时间步长将被调整到对应这些项目中的任何一个所计算出的最小值)。
对于显式蠕变( OPTION=0),如果比值Δε cr/ ε el 高于0.25的稳定界限,且如果时间增量不能被减小,解可能发散且分析将由于错误信息而终止。
这个问题可以通过使最小时间步长足够小来避免[ DELTIM 和NSUBST ]。
对于隐式蠕变( OPTION=1),缺省无最大蠕变极限,但用户可以指定任意的蠕变率控制。
命令:CRPLIM
GUI:Main Menu>Solution>Unabridged Menu>-Load Step Opts-
Nonlinear>Creep Criterion
注意--如果在分析中不需要包括蠕变效应,则应用RATE 命令及Option=OFF,或把时间步设置成比前一个时间步长些,但不
大于1.0e-6。
2.2.
3.2.2 时间步开放控制
这个选项可用于热分析(记住用户不能通过求解控制对话框来设置热分析选项,必须用ANSYS标准命令集或相应菜单来设置)。
这个选项的主要应用是最终温度达到稳态的非稳态热分析。
在这种情况下,时间步可很快开放。
其缺省值是,如果TEMP增量在三个连
续子步中小于0.1(NUMSTEP=3),则时间步大小可以为―开放‖(缺省值=0.1)。
然后时间步被连续增加以加快求解效率,。
命令:OPNCONTROL
GUI:Main Menu>Solution>Unabridged Menu>-Load Step Opts-
Nonlinear>Open Control
2.2.
3.2.3 求解监视
这个选项为监视指定节点上的指定自由度的求解值提供了方便。
这个命令为用户快速观察求解收敛效率提供了可能,而不必通过冗长的输出文件来取得这些信息。
例如,在一个子步上尝试次数过大,这个文件包含的信息将提供指示:要么降低初始时间步,要么增加最小的子步数,这可通过NSUBST 命令来避免二分次数过多。
命令:MONITOR
GUI:Main Menu>Solution>Unabridged Menu>-Load Step Opts-Nonlinear>Monitor
2.2.
3.2.4 激活和杀死选项
根据需要指定―生‖、―死‖选项。
对选定的单元,可以―杀死‖[ EKILL ]和―激活‖[ EALIVE ],以模拟在结构中移走或添加材料。
作为标准的―生‖、―死‖方法以外的另一个方法,用户可以对所选择的单元在荷载步之间改变材料特性[ MPCHG ]。
命令:EKILL
EALIVE
GUI:Main Menu>Solution>-Load Step Opts-Other>Kill Elements
Main Menu>Solution>-Load Step Opts-Other>Activate Elem
程序通过用一个非常小的数(它由ESTIF 命令设置)乘以它的刚度并从总质量矩阵消去它的质量来―杀死‖一个单元。
对杀死单元的单元载荷(压力、热通量、热应变等等)同样地设置为零。
用户需要在前处理中定义所有可能的单元,用户不可能在SOLUTION 中产生新的单元。
要在用户的分析的后面阶段中―激活‖的那些单元,在第一个载荷步前应当被―杀死‖,然后在适当的载荷步的开始被重新―激活‖。
当单元被重新―激活‖时,它们具有零应变状态,且(如果NLGEOM ,ON )它们的几何构形(长度、面积等等) 被修改来与它们现在变形后的位置相适应。
参见《ANSYS Advanced Analysis Techniques Guide》。
另一个在求解过程中影响单元行为的方法是修改选定单元的材料特性:
命令:MPCHG
GUI:Main Menu>Solution>-Load Step Opts-Other>Change Mat Props>Change Mat Num
注意--应用[ MPCHG ]是要注意。
在求解期间改变它的材料性质参考号,可能产生不希望的结果,特别是如果用户改变材料非线性特性[ TB ]。
2.2.
3.2.5 输出控制选项
除了可以通过求解控制对话框可以设置的OUTRES 外,用户还可以设置其他输出选项。
命令:OUTPR
ERESX
GUI:Main Menu>Solution>Unabridged Menu>-Load Step Opts-Output
Ctrls>Solu Printout
Main Menu>Solution>Unabridged Menu>-Load Step Opts-Output
Ctrls>Integration Pt
打印输出选项[ OUTPR ]可在输出文件( Jobname . OUT )中包括所想要的任何结果数据。
结果外推[ ERESX ]拷贝一个单元的积分点应力和弹性应变结果到结点来替代外推,如果在单元中存在非线性(塑性、蠕变、膨胀)的话。
积分点非线性应变总是被拷贝到结点。
参见《ANSYS Basic Analysis Guide》§2。
2.2.4 施加荷载
在这一步把荷载施加到模型中,参见《ANSYS Basic Analysis Guide》§2。
记住惯性荷载和点荷载将保持方向不变,但表面荷载在大变形分析中将跟随结构的变形而变化。
用户可以定义一维数据表(TABLE类型的数组参数)来施加复杂边界条件。
2.2.5 求解
1、把数据库保存为一个文件。
命令:SAVE
GUI:Utility Menu>File>Save as
2、求解
命令:SOLVE
GUI:Main Menu>Solution>-Solve-Current LS
3、如用户定义了多个荷载步,则必须指定时间设置、荷载步选项等,然后保存和求解每个附加的荷载步。
参见《ANSYS Basic Analysis Guide》。
4、退出求解器
命令:FINISH
GUI:关闭求解菜单
2.2.6 考察结果
非线性静态分析的结果,主要由位移、应力、应变以及反作用力组成。
可以用通用后处理器POST1,或者时间历程后处理器POST26,来考察这些结果。
记住用POST1一次仅可以读取一个子步,且来自那个子步的结果应当已被写入Jobname.RST 。
(载荷步选项命令OUTRES 控制哪一个子步的结果被存储入Jobname.RST )。
典型的POST1后处理顺序将在下面描述。
2.2.6.1 要记住的要点
用POST1考察结果,数据库中的模型必须与用于求解计算的模型相同。
结果文件( Jobname.RST )必须是可用的。
2.2.6.2 用POST1考察结果
1、检查用户的输出文件( Jobname.OUT )是否在所有的子步分析都收敛。
如果不收敛,用户可能不想进行后处理,而是想确定为什么收敛失败。
如果用户的解收敛,那么继续进行后处理。
2、进入POST1。
如果用于求解的模型现在不在数据库中,发出RESUME 命令。
命令:/ POST1
GUI:Main Menu>General Postproc
3、读取需要的载荷步和子步结果,这可以依据载荷步和子步号或者时间来识别,然而不能依据时间来识别出弧长法结果。
命令:SET
GUI:Main Menn>General Postproc>Read Results-Load step
同样地用户可以使用SUBSET 或者APPEND 命令来只对选出的部分模型读取或者合并结果数据。
这些命令中的任何一个中
的LIST参数列出结果文件中可用的解。
用户同样地可以通过INRES 命令限制从结果文件到基本数据被写的数据总量。
另外可以用ETABLE 命令对选出的单元进行后处理,见《ANSYS Commands Reference》
警告:如果用户指定了一个没有结果可用的TIME值,ANSYS程序将进行线性内插来计算出那Time处的结果。
认识到在非线分析中这种线性内插通常将导致某些精度损失(参看图2-1 )。
因此,对于非线性分析,通常用户应当在一个精确地对应于要求子步的TIME 处进行后处理。
图2-1 非线性结果的线性内插可能引起某些误差
4、使用下列任意选项显示结果
1)显示已变形的形状
命令:PLDISP
GUI:Main Menu>General Postproc>Plot Results>Deformed Shapes
在大变形分析中,一般优先使用真实比例显示[ DSCALE ,,1]。
2)等值线显示
命令:PLNSOL 或PLESOL
GUI:Main Menu>General Postproc>Plot Results>-Contour Plot-
Nodal Solu 或Element Solu
使用这些选项来显示应力、应变或者任何其它可用项目的等值线。
如果邻接的单元具有不同材料行为(可能由于塑性或多线性弹性的材料性质,由于不同的材料类型,或者由于邻近的单元的死活属性不同而产生),用户应当注意避免结果中的结点应力平均错误。
PLNSOL 和PLESOL 命令的KUND域使用户可以在原始图形上叠加变形图。
同样地用户可以绘制单元表数据和线单元数据的等值线:
命令:PLETAB ,PLLS
GUIS:Main Menu>General Postproc>Element Table>Plot Element Table
Main Menu>General Postproc>Plot Results>-Contour Plot-Line Elem Res
使用PLETAB 命令来绘制单元表数据的等值线,用PLLS 命令来绘制线单元数据的等值线。
3)列表
命令:PRNSOL (结点结果),
PRESOL (结果),
PRRSOL (反作用力数据)
PRETAB
PRITER (子步总计数据)等等。
NSORT
ESORT
GUIS:Main Menu>General Postproc>List Results>Nodal Solution
Main Menu>General Postproc>List Results>Element Solution
Main Menu>General Postproc>List Results>Reaction Solution
使用NSORT 和ESORT 命令在将数据列表前对它们进行排序。
4)其它的性能
在POST1中还可用许多其它的后处理功能(如在路径上映射结果,记录参量列表,等等),见《ANSYS Basic Analysis Guide》。
对于非线性分析,载荷工况组合通常是无效的。
2.2.6.3 用POST26考察结果
用户可以使用时间─历程后处理器POST26来考察非线性结构的载荷─历程响应。
使用POST26比较一个ANSYS变量对另一个变量的关系。
例如,用户可以用图形表示某一结点处的位移与对应的所加载荷的关系,或者用户可以列出某一结点处的塑性应变和对应的TIME值之间的关系。
典型的POST26后处理顺序可以遵循这些步骤:
1、根据用户的输出文件( Jobname.OUT )检查是否在所有要求的载荷步内分析都收敛。
用户不应当将设计决策建立在不收敛结果的基础上。
2、如果用户的解收敛,进入POST26,如果用户的模型不在数据库内,发出RESUME 命令。
命令:/ POST26
GUI:Main Menu>TimeHist Postpro
3、定义在后处理期间使用的变量
命令:NSOL,ESOL,RFORCE
GUI:Main Menu>Time Hist Postproc>Define Variables
4、图形或者列表显示变量
命令:PLVAR (图形表示变量),PRVAR ,EXTREM (列表变量)
GUIS:Main Menu>Time Hist Postprac>Graph Variable S
Main Menu>Time Hist Postproc>List Variables
Main Menu>Time Hist Postproc>List Extremes
5、其它的性能
许多其它的后处理函数可用于POST26,参考《ANSYS Basic Analysis Guide》§6。
此外还可参见
NLGEOM,SSTIF,NROPT,TIME,NSUBST,AUTOTS,KBC, CNVTOL,NEQIT, NCNV,PRED,OUTES 和SOLU 命令的说明。
2.2.7 终止正在运行的工作,重起动
用户可以通过产生一个―abort‖文件(Jobname.ABT)停止一个非线性分析,见《ANSYS Basic Analysis Guide》§3。
一旦求解成功地完成,或者收敛失败发生,程序也将停止分析。
如果一个分析在终止前已成功地完成了一次或多次迭代,用户可以屡次重启动它。
见《ANSYS Basic Analysis Guide》§3.16。
2.3 非线性瞬态分析步骤
许多需要进行非线性瞬态分析的任务,与非线性静力分析(参见§2.2)和线性完全瞬态分析相同或相似。
本节论述非线瞬态分析的一些附加考虑。
请记住§2.2论述的求解控制对话框,不能应用于热分析的求解控制,只能应用标准的ANSYS命令集或菜单来进行热分析的设置。
2.3.1 建模
这一步骤与非线性静力分析相同,参见§2.2。
但是,如果分析中包含时间的积分效应,则必须输入质量密度[ MP ,DENS]。
如果需要,还可以定义与材料相关的结构阻尼[ MP ,DAMP]。
2.3.2 施加荷载和求解
1、指定瞬态分析类型,定义分析选项,与非线性静力分析相同:
新的分析或重启动[ ANTYPE ]
分析类型:瞬态[ ANTYPE ]
大变形效应[ NLGEOM ]
大位移瞬态(如果用求解控制对话框设置分析类型)。
2、施加荷载,并指定荷载步选项,这与线性完全瞬态动力分析中相同。
瞬态时间历程通常需要多个荷载步,其中第1荷载步典型地用于建立初始条件,见《ANSYS Basic Analysis Guide》。
此外,非线性静力分析中所用的一般的非线性、生和死、输出控制等,在非线性瞬态分析中也可应用。
在非线性瞬态分析中,时间必须大于0。
对于非线性瞬态分析,用户必须说明是阶梯荷载还是斜坡荷载[ KBC ]。
见《ANSYS Basic Analysis Guide》对此的进一步论述。
命令:ALPHAD
BETAD
TIMINT
TINTP
GUI:Main Menu>Solution>-Analysis Type-Sol"n Control:Transient Tab
Main Menu>Solution>Unabridged Menu>-Load Step Opts-Time/
Frequenc>Damping
Main Menu>Solution>Unabridged Menu>-Load Step Opts-Time/
Frequenc>Time Integration
动力选项解释:
⑴阻尼--Rayleigh阻尼常数用常数质量[ ALPHAD ]和刚度[ BETAD ]矩阵乘子定义。
在非线性分析中,刚度可能激烈改变--除特殊情况外,不要应用BETAD 。
⑵时间积分效应[ TIMINT ]。
只在瞬态分析中,时间积分效应才缺省打开。
对于蠕变、粘弹性、粘塑性、膨胀,应当关闭时间积分效应(也就是说明进行静力分析)。
这些时间相关效应通常不包括在动力分析中,因为瞬态动力时间步,对于任何明显的长期变形来说,时间太短。
除了在运动学(刚体运动)分析中,用户应当很少需要调整瞬态积分参数[ TINTP ]--它对Newmark方程提供数值阻尼,参见《ANSYS Theory Reference》。
ANSYS的自动求解控制,把缺省设为一个新的时间积分方案,对于应用一阶瞬态方程。
这通常用于不稳定状态热问题(θ=1)(由SOLCONTROL ,ON 设置),这是反向EULER 方案。
它是无条件稳定的。
对于象相变这样的高度非线性热问题,这种方案更有效。
振荡极限容限缺省为0.0,以使响应的一阶特征值可用于更精确地决定一个新的时间步值。
注意--如果用求解控制对话框设置求解控制,用户可在Transient标签中进入所有这些选项。
3、把各个荷载步的荷载数据写到荷载步文件中。
命令:LSWRITE
GUI:Main Menu>Solution>Write LS File
4、把数据库备份到一个命名文件中。
命令:SAVE
GUI:Utility Menu>File>Save As
5、开始求解。
对于多荷载步的求解参见《ANSYS Basic Analysis Guide》§1。
命令:LSSOLVE
GUI:Main Menu>Solution>-Solve-From LS Files
6、在求解完所有荷载步后,退出求解。
命令:FINISH
GUI:关闭Solution菜单。