牛顿一拉森方法之ANSYS选项

2017-03-02  by:CAE仿真在线  来源:互联网

1、进入ANSYS求解器

命令:/Solution

GUI:Main Menu>Solution

2、定义分析类型及分析选项。

分析类型和分析选项在第一个载荷步后(也就是,在你发出你的第一个SOLVL命令之后)不能被改变。ANSYS提供这些选项用于静态分析。

表1─1 分析类型和分析选项

这些选项中的每一个都将在下面详细地解释。

选项:新的分析〔ANTYPE〕

一般情况下会使用New Analysis(新的分析)。

选项:分析类型:静态〔ANTYPE〕

选择Static(静态)。

选项:大变形或大应变选项(GEOM)

并不是所有的非线性分析都将产生大变形。

选项:应力刚化效应〔SSTIF〕

如果存在应力刚化效应选择ON。

选项:牛顿-拉普森选项〔NROPT〕

仅在非线性分析中使用这个选项。这个选项指定在求解期间每隔多久修改一次正切矩阵。你可以指定这些值中的一个。

程序选择(NROPT,ANTO):程序基于你模型中存在的非线性种类选择用这些选项中的一个。在需要时牛顿-拉普森方法将自动激活自适应下降。

完全〔NROPT,FNLL〕;程序使用完全的牛顿-拉普森处理方法,在这种处理方法中每进行一次平衡迭代修改刚度矩阵一次。如果自适应下降是关闭的,程序每一次平衡迭代都使用正切刚度矩阵。(我们一般不建议关闭自适应下降,但是你或许发现这样做可能更有效。)如果自适应下降是打开的(缺省),只要迭代保持稳定(也就是,只要残余项减小,且没有负主对角线出现)程序将仅使用正切刚度阵。如果在一次迭代中探测到发散倾向,程序抛弃发散的迭代且重新开始求解,应用正切和正割刚度矩阵的加权组合。当迭代回到收敛模式时,程序将重新开始使用正切刚度矩阵。对复杂的非线性问题自适应下降通常将提高程序获得收敛的能力。

修正的(NROPT,MODI):程序使用修正的牛顿-拉普森方法,在这种方法中正切刚度矩阵在每一子步中都被修正。在一个子步的平衡迭 代期间矩阵不被改变。这个选项不适用于大变形分析。自适应下降是不可用的

初始刚度(NROPT,INIT):程序在每一次平衡迭代中都使用初始刚度矩阵这一选项比完全选项似乎较不易发散,但它经常要求更多次的迭代来得到收敛。它不适用于大变形分析。自适应下降是不可用的。

选项:方程求解器

对于非线性分析,使用前面的求解器(缺省选项)。

3、在模型上加载,记住在大变型分析中惯性力和点载荷将保持恒定的方向,但表面力将“跟随”结构而变化。

4、指定载荷步选项。这些选项可以在任何载荷步中改变。下列选项对非线性静态分析是可用的:

普通选项

普通选项包括下列:

Time(TIME)

ANSYS程序借助在每一个载荷步末端给定的TIME参数识别出载荷步和子步。使用TIME命令来定义受某些实际物理量(如先后时间,所施加的压力,等等。)限制的TIME值。程序通过这个选项来指定载荷步的末端时间。注意──在没有指定TIME值时,程序将依据缺省自动地对每一个载荷步按1.0 增加TIME(在第一个载荷步的末端以TIME=1.0开始)。

子步的数目〔NSUBST〕

时间步长〔DELTIM〕

非线性分析要求在每一个载荷步内有多个子步(或时间步;这两个术语是等效的)从而ANSYS可以逐渐施加所给定的载荷,得到精确的解。NSUBST和DELTIM命令都获得同样的效果(给定载荷步的起始,最小,及最大步长)。NSNBST 定义在一个载荷步内将被使用的子步的数目,而DELTIM明确地定义时间步长。如果自动时间步长是关闭的,那么起始子步长用于整个载荷步。缺省时是每个载荷步有一个子步。

渐进式或阶跃式的加载

在与应变率无关的材料行为的非线性静态分析中通常不需要指定这个选项,因为依据缺省,载荷将为渐进式的阶跃式的载荷〔KBC,1〕除了在率─相关材料行为情状下(蠕变或粘塑性),在静态分析中通常没有意义。

自动时间步长〔AUTOTS〕

这一选项允许程序确定子步间载荷增量的大小和决定在求解期间是增加还是减 小时间步(子步)长。缺省时是OFF(关闭)。你可以用AUTOTS命令打开自动时间步长和二分法。通过激活自动时间步长,可以让程序决定在每一个载荷步内使用多少个时间步。在一个时间步的求解完成后,下一个时间步长的大小基于四种因素预计:

在最近过去的时间步中使用的平衡迭代的数目(更多次的迭代成为时间步长减小的原因)

对非线性单元状态改变预测(当状态改变临近时减小时间步长)

塑性应变增加的大小

蠕变增加的大小

非线性选项

程序将连续进行平衡迭代直到满足收敛准则(或者直到达到允许的平衡迭代的 最大数〔NEQIT〕。我们可以用缺省的收敛准则,也可以自己定义收敛准则。

收敛准则〔CNVTOL〕

缺省的收敛准则

依据缺省,程序将以VALUE?TOLER的值对力(或者力矩)进行收敛检查。VALUE的缺省值是在所加载荷(或所加位移,Netwton-Raphson回复力)的SRSS,和MINREF(其缺省为1.0)中,取值较大者。TOLER的缺省值是0.001。你应当几乎总是使用力收敛检查。可以添加位移(或者转动)收敛检查。对于位移,程序将收敛检查建立在当前(i)和前面(i─1)次迭代之间的位移改变上。注意──如果你明确地定义了任何收敛准则(CNVTOL〕,缺省准则将“失效”。因此,如果你定义了位移收敛检查,你将不得不再定义力收敛检查(使用多个CNVTOL命令来定义多个收敛准则)。

用户收敛准则

你可以定义用户收敛准则,替代缺省的值。使用严格的收敛准则将提高你的结果的精度,但以多更次的平衡迭代为代价。如果你想严格(加放松)你的准则,你应当改变TOLER两个数量级。一般地,你应当继续使用VALUE的缺省值;也就是,通过调整TOLER,而不是VALUL 改变收敛准则。你应当确保MINREF=1.0的缺省值在你的分析范围内有意义。

在单一和多DOF系统中检查收敛

要在单自由度(DOF)系统中检查收敛,你对这一个DOF计算出不平衡力,然后对照给定的收敛准则(VALUE*TOLER)参看这个值(同样也可以对的单一DOF的位移(和旋度)收敛进行类似的检查。)然而,在多DOF系统中,你也许想使用不同的比较方法。

ANSYS程序提供三种不同的矢量规范用于收敛核查。

无限规范在你模型中的每一个DOF处重复单-DOF核查。

LI规范将收敛准则同所有DOFS的不平衡力(力矩)的绝对值的总和相对照。

L2规范使用所有DOFS不平衡力(或力矩)的平方总和的平方根进行收敛检查。

实例

对于下面例子,如果不平衡力(在每一个DOF处单独检查)小于或等于5000?0.0005(也就是2.5),且如果位移的改变(以平方和的平方根检查)小于或等于10?0.001(也就是0.01),子步将认为是收敛的。

CNVTOL,F,5000,0.005,0

CNVTOL,U,10,0.001,2

平衡迭代的最大次数〔NEQIT〕

使用这个选项来对在每一个子步中进行的最大平衡迭代次数实行限制(缺省=25)。如果在这个平衡迭代次数之内不能满足收敛准则,且如果自动步长是打开的〔AUTOTS〕,分析将尝试使用二分法。如果二分法是不可能的,那么,分析将或者终止,或者进行下一个载荷步,依据你在NCNV命令中发出的指示。

求解终止选项〔NCNV〕

这个选项处理五种不同的终止准则:

如果位移“太大”它建立一个用于终止分析和程序执行的准则。

它对累积迭代次数设置限制。

它对整个时间设置限制。

它对整个CPU时间设置限制。

弧长选项〔ARCLEN〕

如果你预料结构在它的载荷历史内在某些点将变得物理意义上不稳定(也就是, 结构的载荷—位移曲线的斜度将为0或负值),你可以使用弧长方法来帮助稳定数值求解。注意──当合适时,你可以和弧长方法一起使用许多其它的分析和载荷步选项。然而,你不应和弧长方法一起使用下列选项:不要使用线搜索〔LNSRCH〕,时间步长预测〔PRED〕,自适应下降〔NROPT,,,ON〕,自动时间步长〔AUTOTS,TIME,DELTIM〕,或打开时间-积分效应(TIMINT)。

时间步长预测──纠正选项〔PRED〕

对于每一个子步的第一次平衡迭代你可以激活和DOF求解有关的预测。这个特点加速收敛且如果非线性响应是相对平滑的,它特别的有用。在包含大转动或粘弹 的分析中它并不是非常有用。

线搜索选项〔LNSRCH〕

这个选项是对自适应下降的替代。当被激活时,无论何时发现硬化响应。这个收敛提高工具用程序计算出的比例因子(具有0和1之间的值)乘以计算出的位移增量。因为线搜索算法是用来对自适应下降选项〔NROPT〕进行的替代,如果线搜索选项是开,自适应下降不被自动激活。不建议你同时激活线搜索和自适应下降。

当存在强迫位移时,直到迭代中至少有一次具有一个的线搜索值运算才会收敛。ANSYS调节整个DU矢量,包括强迫位移值;否则,除了强迫DOF处一个小的位移值将随处发生。直到适代中的某一次具有1的线搜索值,ANSYS才施加全部位移值。

蠕变准则〔CRPLIM,CRCR〕

如果结构表现出蠕变行为,可以指定蠕变准则用于自动时间步调整。(如果自动时间步长〔AUTOTS〕不是打开的,这个蠕变准则将无效。)程序将对所有单元计算蠕应变增量(在最近时间步中蠕变的变化)对弹性应变的比值。如果最大比值比判据大,程序将减小下一个时间步长;如果小,程序或许增加下一个时间步长。(同样地程序将把自动时间步长建立在平衡迭代次数,即将发生的单元状态改变,以及塑性应变增量的基础上。时间步长将被调整到对应这些项目中的任何一个所计算出的最小值。)如果比值高于0.25的稳定界限,且如果时间增量不能被减小,解可能发散且分析将由于错误信息而终止。这个问题可以通过使最小时间步长足够小避免〔DELTIM,NSUBST)。

激活和杀死选项

在ANSYS/Mechanical和ANSYS/LS-DYNA产品中,你可以去杀死和激活单元来模拟材料的消去和添加。

程序通过用一个非常小的数(它由ESTIF命令设置)乘以它的刚度从总质量矩阵消去它的质量“ 杀死”一个单元。对无活性单元的单元载荷(压力,热通量,热应变,等等)同样地设置为零。你需要在前处理中定义所有可能的单元;你不可能SOLUTION中产生新的单元。要在你的分析的后面阶段中“出生”的那些单元,在第一个载荷步前应当被杀死,然后去在适当的载荷步的开始被重激活,当单元被重激活时,它们具有零应变状态,且(如果NLGEOM, ON) 它们的几何(开头长度,面积等等) 被修改来与它们的的现偏移位置相适应。

杀死(EKILL)

激活(EALIVE)

改变材料性质参考号〔MPCHG〕

另一种在求解期间影响单元行为的办法是来改变它的材料性质参考号。这个选项允许你在载荷步间改变一个单元的材料性质。

EKILL适用于大多数单元类型。MPCHG适用于所于单元类型。

输出控制选项

输出控制选项包括下列:

打印输出(OUTPR)

使用这个选项来在输出文件(Jobname.out)中包括进便所想要的结果数据。

结果文件输出〔OUTRES〕

这个选项控制结果文件中的数据(Jobname.rst)。

OUTPR和OUTRES用来控制结果被写入这些文件的频率。

结果外推〔ERESX〕

这个选项,依据缺省,拷贝一个单元的积分点应力和弹性应变结果到结点而替代外推它们,如果在单元中存在非线性(塑性,蠕变,膨胀)的话。积分点非线性变化总是被拷贝到结点。

注意:对输出行使下列警告:

恰当使用多个OUTRES或OUTPR命令有时可能有一点小的技巧。

依据缺省,在非线性分析中只有最后一个子步被写入结果文件。要写入所有子步,设置OUTRES中的FREQ域为ALL。

依据缺者,只有1000个结果集(子步)可以被写入结果文件。如果超过了这个数目(基于你的OUTRES指定),程序将由于错误而终止。使用命令/CONFIG,NRES来增加这个界限。

5、存储基本数据的备份副本于另一文件。

命令:SAVE

GUI:Utility Menu>File>Save As

6、开始求解计算。

命令:SOLVE

GUI:Main Menn>Solution>-Solve-Current LS

7、如果你需要定义多个载荷步,对每一个其余的载荷步重复步骤3至6。

8、离开SOLUTION处理器

命令:FINISH

GUI:关闭Solution菜单。

相关标签搜索：牛顿一拉森方法之ANSYS选项 Ansys有限元培训 Ansys workbench培训 ansys视频教程 ansys workbench教程 ansys APDL经典教程 ansys资料下载 ansys技术咨询 ansys基础知识 ansys代做 Fluent、CFX流体分析 HFSS电磁分析 Abaqus培训