《非线性分析》由会员分享，可在线阅读，更多相关《非线性分析（17页珍藏版）》请在金锄头文库上搜索。
1、ANSYS 非线性分析命令解析ANSYS 2009-02-27 19:46:43 阅读 1056 评论 2 字号：大中小 ANSYS 应用基于问题物理特性的自动求解控制方法，把各种非线性分析控制参数设置到合适的值。如果用户对这些设置不满意，还可以手工设置。下列命令的缺省设置已进行了优化处理：AUTOTS PRED MONITOR DELTIM NROPT NEQIT NSUBST TINTP SSTIF CNVTOL CUTCONTROL KBC LNSRCH OPNCONTROL EQSLV ARCLEN CDWRITE LSWRITE 这些命令及其设置在将在后面讨论。参见ANSYS Com
2、mands Reference。如果用户选择自己的设置而不是 ANSYS 的缺省设置，或希望用以前版本的 ANSYS 的输入列表，则可用/ SOLU 模块的 SOLCONTROL ，OFF 命令，或在 / BATCH 命令 后用/ CONFIG ，NLCONTROL ，OFF 命令。参见 SOLCONTROL 命令的详细描述。ANSYS 对下面的分析激活自动求解控制单场的非线性或瞬态结构以及固体力学分析，在求解自由度为UX、 UY、UZ 、ROTX、ROTY、 ROTZ 的结合时；单场的非线性或瞬态热分析，在求解自由度为 TEMP 时；注意 - 本章后面讨论的求解控制对话框，不能对热分析做设置
3、。用户必须应用标准的 ANSYS求解命令或 GUI 来设置。 2.2 非线性静态分析步骤尽管非线性分析比线性分析变得更加复杂，但处理基本相同。只是在非线形分析的过程中，添加了需要的非线形特性。 非线性静态分析是静态分析的一种特殊形式。如同任何静态分析，处理流程主要由以下主要步骤组成：建模；设置求解控制；设置附加求解控制；加载；求解；考察结果。2.2.1 建模 这一步对线性和非线性分析基本上是一样的，尽管非线性分析在这一步中可能包括特殊的单元或非线性材料性质，参考4材料非线性分析，和6.1单元非线性。如果模型中包含大应变效应，应力 应变数据必须依据真实应力和真实(或对数)应变表示。参见ANSYS
4、 Modeling and Meshing Guide。在 ANSYS 中建立了模型后，应该设置求解控制(分析类型、分析选项、荷载步等)选项，施加荷载，最后求解。非线性分析与线性分析的不同之处是，前者需要许多荷载增量，并且总是需要平衡迭代。下面讨论一般过程。参见本章的例子。2.2.2 设置求解控制设置求解控制包括定义分析类型、设置分析的常用选项和指定荷载步选项。在做结构非线性静态分析时，可以应用求解控制对话框来设置。该对话框对许多非线性静态分析提供了缺省设置。这样，用户需要的设置降低到最少。求解控制框的缺省设置，基本上与2.1 所述的自动求解控制的设置相同。由于求解控制对话框是非线性静态分析的
5、推荐工具，我们在下面将详细论述，如用户不想用这个对话框(GUI：Main MenuSolution-Analysis Type-Soln Control)，可以应用标准的 ANSYS 求解命令集或相应的菜单(GUI ：Main Menu SolutionUnabridged Menuoption)。求解控制对话框的概况，见ANSYS Basic Analysis Guide3.11 。注意 - 对于非线性结构完全瞬态分析，建议应用求解控制对话框，但并不是必须如此，见 2.3。 2.2.2.1 求解控制对话框 进入选择(GUI：Main MenuSolution-Analysis Type-So
6、ln Control)进入求解控制对话框。下面几节将论述这个求解对话框中的内容。对于其详细说明，可以在相应标签下，按 HELP 按钮进入帮助系统。2.2.2.2 求解控制对话框 -Basic 标签求解控制对话框共有五个标签，其中最基本的选项位于第一个标签上，其他标签依此提供更高级的控制。进入对话框后，缺省的标签就是 Basic 标签。Basic 标签中的内容，提供了 ANSYS 分析所需要的最少设置。如果用户对 Basic 标签中的设置满意，就不必调整其他标签中的更高级的设置。在按 OK 按钮以后，设置才作用于 ANSYS 数据库，并关闭对话框。可用的 Basic 标签选项见 表 2-1 。按
7、 HELP 可得更多的说明。表 2-1选项 参见ANSYS Basic Analysis Guide指定分析类型 ANTYPE , NLGEOM 1.2.6.13.16控制时间设置，包括:荷载步末的时间 TIME , 自动时间步 AUTOTS ,一个荷载步中的子步数 NSUBST 或 DELTIM 2.42.7.1指定写入数据库中的结果数据 OUTRES 2.7.4在非线性静态分析中的一些特殊考虑如下：1、在设置 ANTYPE 和 NLGEOM 时，如果是执行新的分析，选择Large Displacement Static，但要记住并不是所有的非线性分析都产生大变形，见3。如果想重启动一个已失
8、败的非线性分析，选择Restart Gurrent Analysis。在第 1 荷载步以后( 即在首次运行 SOLVE 命令后) ，用户不能改变这个设置。通常用户要作一个新的分析，而不是重启动分析。重启动分析的讨论见ANSYS Basic Analysis Guide。2、在进行时间设置时，记住这些选项可在任何荷载步改变。参见ANSYS Basic Analysis Guide 2。高级的时间/频率选项，参见 2.2.2.8。非线性分析要求在一个时间步上有多个子步，以使ANSYS 能够逐渐地施加荷载，并取得精确解。 NSUBST 和 DELTIM 命令产生相同的效果(建立荷载步的开始、最小和最
9、大时间步)，但互为倒数。 NSUBST 定义一个荷载步上的子步数，而 DECTIM 显式地定义时间步大小。如果自动时间步 AUTOTS 关闭，则起始子步大小用于整个荷载步。3、 OUTRES 控制结果文件(Jobname.RST)中的数据。缺省时，在非线性分析中把最后一个子步的结果写入此文件。结果文件只能写入 1000 个结果集( 子步)，但用户可以用 / CONFIG ，NRES 命令来增大这一限值，参见ANSYS Basic Analysis Guide。2.2.2.3 求解控制对话框 -Transient 标签这个标签的内容是瞬态分析控制，只有在 Basic 标签中选择了瞬态分析时这个标
10、签才能应用，否则呈灰色。所以在这里不论述，参见2.3。2.2.2.4 求解控制对话框 -Soln Options 标签这个标签设置的选项见 表 2-2 。按本标签的 HELP 可得到更多的说明。表 2-2选项 参见指定方程求解器 EQSLV 2.2.2.7.1 ANSYS Basic Analysis Guide 3.2-3.10对多重启动指定参数 RESCONTROL ANSYS Basic Analysis Guide 3.16.22.2.2.5 求解控制对话框 -Nonlinear 标签用 Nonlinear 标签设置的选项见 表 2-3 。按 HELP 按钮进入帮助系统可得到更多的说明
11、。表 2-3选项 参见激活线性搜索 LNSRCH 2.2.2.8.52.4激活自由度求解预测 PRED 2.2.2.8.4指定一个荷载步中的最大子步数 NEQIT 2.2.2.8.3指定是否需要包括蠕变计算 RATE 4.42.2.3.2.1设置收敛准则 CNVTOL 2.2.2.8.2控制二分 CUTCONTROL 2.2.2.8.62.2.2.6 求解控制对话框 -Advanced NL 标签 用 Advanced NL 标签设置的选项见 表 2-4 。按 HELP 按钮进入帮助系统可得到更多的说明。表 2-4选项 参见 指定分析终止准则 NCNV 2.2.2.8.3激活和终止弧长法的控制
12、 ARCLEN , ARCTRM 2.4ANSYS Basic Analysis Guide22.2.2.7 求解控制对话框 -设置其他高级分析选项2.2.2.7.1 方程求解器ANSYS 的自动求解控制在大多数情况下，激活稀疏矩阵直接求解器( EQSLV ，SPARSE)。这是缺省的求解器，除了在子结构分析的生成步骤外(这时用波前直接求解器) 。其他选项包括波前直接求解器和 PCG 求解器。对于实体单元 (如 SOLID92 或 SOLID45)，使用 PCG 求解器可能更快，尤其是在三维模型中。如果用户采用 PCG 求解器，可以考虑用 MSAVE 命令降低内存应用。 MSAVE 命令对于线
13、性材料特性的 SOLID92 单元，触发单元方法。为了应用这一命令，必须是小应变( NLGEOM ,OFF)静力或完全瞬态分析。模型中不符合上述条件的其他部分，应用总体集成刚度矩阵来求解。对于符合上述条件的模型部分，用 MSAVE ,ON 可能可节省 70%的内存，但求解时间可能增加，这与计算机的配置和 CPU速度有关。与 ANSYS 中的迭代求解器不同，稀疏矩阵求解器是一个强大的求解器。虽然 PCG 求解器能够求解不定矩阵方程，但在它碰到一个病态矩阵时，如果不能收敛，求解器将迭代至指定的迭代次数后停止迭代。在发生这种问题时，它触发二分。在完成二分后，如果矩阵是良态的，求解器继续求解。最后整个
14、非线性荷载步可以得到求解。在结构非线性分析中，选择稀疏矩阵求解器，还是选择 PCG 求解器，可参照下面的建议：1、如果是梁、壳或者梁、壳、实体结构，选择稀疏矩阵求解器；2、如果是三维结构，而且自由度数相对较大(200000 个自由度或以上)，选择 PCG 求解器；3、如果问题是病态(由不良单元形状引起)，或在模型的不同区域材料特性相差巨大，或者位移边界条件不足，选择稀疏矩阵求解器。2.2.2.8 求解控制对话框 -设置其他高级荷载步选项2.2.2.8.1 自动时间步ANSYS 的自动求解控制打开自动时间步长 AUTOTS ，ON。这一选项允许程序确定子步间载荷增量的大小和决定在求解期间是增加还
15、是减小时间步(子步) 长。在一个时间步的求解完成后，下一个时间步长的大小基于四种因素预计：在最近过去的时间步中使用的平衡迭代的数目(更多次的迭代成为时间步长减小的原因) ；对非线性单元状态改变预测(当状态改变临近时减小时间步长 )；塑性应变增加的大小；蠕变增加的大小。2.2.2.8.2 收敛准则 程序将连续进行平衡迭代直到满足收敛准则 CNVTOL (或者直到达到允许的平衡迭代的最大次数 NEQIT 。如果缺省的收敛准则不满意，可以自己定义收敛准则。ANSYS 的自动求解控制应用等于 0.5%的力( 或力矩)的 L2-范数容限 (TOLER)，这对于大部分情况合适。在大多数情况下，除了进行力范数的检查外，还进行 TOLER 等于 5%的位移 L2-范数的检查。缺省时，程序将通过比较不平衡力的平方和的平方根（SRSS）与 VALUETOLER 的值来检查力(在包括转动自由度时，还有力矩 )的收敛。 VALUE 的缺省值是所加载荷 (或在施加位移时，Netwton-Raphson 回复力)的 SRSS，或 MINREF (其缺省为 0.001)，取较大者。如果 SOLCONTROL ，OFF，则对于力的收敛， TOLER 的
