ANSYS非线性行为分析

　本文分析介绍了ANSYS的非线性行为。

　　了解程序的运作方式和结构的表现行为

　　如果以前没有使用过某一种特别的非线性特性，在将它用于大的，复杂的模型前，构造一个非常简单的模型(也就是，仅包含少量单元)，以及确保你理解了如何处理这种特性。

　　%26#8226;通过首先分析一个简化模型，以便使你对结构的特性有一个初步了解。对于非线性静态模型，一个初步的线性静态分析可以使你知道模型的哪一个区域将首先经历非线性响应，以及在什么载荷范围这些非线性将开始起作用。对于非线性瞬态分析，一个对梁，质量块及弹簧的初步模拟可以使你用最小的代价对结构的动态有一个深入了解。在你着手最终的非线性瞬时动态分析前，初步非线性静态，线性瞬时动态，和/或模态分析同样地可以有助于你理解你结构的非线性动态响应的不同的方面。

　　%26#8226;阅读和理解程序的输出信息和警告。至少，在你尝试后处理你的结果前，确保你的问题收敛。对于与路程相关的问题，打印输出的平衡迭代记录在帮助你确定你的结果是有效还是无效方面是特别重的。

　　简化

　　%26#8226;尽可能简化最终模型。如果可以将3─D结构表示为2─D平面应力，平面应变或轴对称模型，那么这样做，如果可以通过对称或反对称表面的使用缩减你的模型尺寸，那么这样做。(然而，如果你的模型非对称加载，通常你不可以利用反对称来缩减非线性模型的大小。由于大位移，反对称变成不可用的。)如果你可以忽略某个非线性细节而不影响你模型的关键区域的结果，那么这样做。

　　%26#8226;只要有可能就依照静态等效载荷模拟瞬时动态加载。

　　%26#8226;考虑对模型的线性部分建立子结构以降低中间载荷或时间增量及平衡迭代所 需要的计算时间。

　　采用足够的网格密度

　　%26#8226;考虑到经受塑性变形的区域要求一个合理的积分点密度。每个低阶单元将提供和高阶单元所能提供的一样多积分点数，因此经常优先用于塑性分析。在重要塑性区域网格密度变得特别地重要，因为大挠度要求对于一个精确的解，个单元的变形(弯曲)不能超过30度。

　　%26#8226;在接触表面上提供足够的网格密度以允许接触应力以一种平滑方式分布。

　　%26#8226;提供足够用于分析应力的网格密度。 那些应力或应变关心的面与那些需要对位移或非线性解析处的面相比要求相对好的网格。

　　%26#8226;使用足够表征最高的重要模态形式的网格密度。所需单元数目依赖于单元的假定位移形状函数，以及模态形状本身。

　　%26#8226;使用足够可以用来分析通过结构的任何瞬时动态波传播的网格密度。如果波传播是重要的，那么至少提供20个单元来分析一个波长。

　　逐步加载

　　%26#8226;对于非保守的，与路径相关的系统，你需要以足够小的增量施加载荷以确保你的分析紧紧地跟随结构的载荷响应曲线。

　　%26#8226;有时你可以通过逐渐地施加载荷提高保守系统的收敛特性，从而使所要求的Newton_Raphson平衡迭代次数最小。

　　合理地使用平衡迭代

　　%26#8226;务必允许程序使用足够多的平衡迭代〔NEQIT〕。在缓慢收敛，路径无关的分析中这会是特别重要的。

　　%26#8226;相反地，在与路径严重相关的情况下，可能不应该增加平衡迭代的最大次数超过程序的缺省值(25)。如果路径相关问题在一个给定的子步内不能快速收敛，那么你的解可能偏离理论载荷响应路径太多。这个问题当你的时间步长太大时出现。通过强迫你的分析在一个较小的迭代次数后终止，你可以从最后成功地收敛的时间步重起动〔ANTYPE〕，建立一个较小的时间步长，然后继续求解。打开二分法%26sup2;AUTOTS，ON〕会自动地用一个较小的时间步长重起动求解。

　　克服收敛性问题

　　如果问题中出现负的主对角元，计算出过度大的位移，或者仅仅没能在给定的最大平衡迭代次数内达到收敛，则收敛失败发生。收敛失败可能表明出结构物物理上的不稳定性，或者也可能仅是有限无模型中某些数值问题的结果。ANSYS程序提供几种可以用来在分析中克服数值不稳性的工具。如果正在模拟一个实际物理意义上不稳定的系统(也就是，具有零或者负的刚度)，那么将拥有更多的棘手问题。有时你可以应用一个或更多的模拟技巧来获得这种情况下的一个解。

　　让我们来探讨一下某些你可以用来尝试提高你的分析的收敛性能的技术。

　　打开自动时间步长

　　%26#8226;当打开自动时间步长时，往往需要一个小的最小的时间步长(或者大的最大的步长数)。

　　%26#8226;当有接触单元(如CONTACT48，CONTACT12，等等)时使用自动时间分步，程序可能趋向于重复地进行二分法直到它达到最小时间步长。然后程序将在整个求解期间使用最小时间步长，这样通常产生一个稳定但花费时间的解。接触单元具有一个控制程序在它的时间步选择中将是多么保守的选项设置(KEYOPT(7))，这样，允许你加速在这些情况下的运行时间。#p#分页标题#e#

　　%26#8226;对于其它的非线性单元，你需要仔细地选择你的最小时间步。如果你选择一个太小的最小时间步，自动时间分步算法可能使你的运行时间太长。相反地，使你的最小时间步长太大可能导致不收敛。

　　%26#8226;务必对时间步长设置一个最大限度(〔DELTIM〕或者〔NSUBST〕)，特别别是对于复杂的模型。这确保所有重要的模态和特性将被精确地包含进。这在下列情况下可能是重要的。

　　%26#8226;具有局部动态行为特性的问题(例如，涡轮叶片和轮毂部件)，在这些问题中系统的低频能量含量以优势压倒高频范围。

　　%26#8226;具有很短的渐进加载时间问题。如果时间步长允许变得太大，载荷历程的渐进部分可能不能被精确地表示出来。

　　%26#8226;包含在一个频率范围内被连续地激励的结构的问题(例如，地震问题)。

　　%26#8226;当模拟运动结构(具有刚体运动的系统)时注意。分析输入或系统驱动频率所要求的时间步通常比分析结构的频率所要求的大几个数量级。采用这样粗略的一个时间步会将相当大的数值干扰引入解中;求解甚至可能变得不稳定。

　　下面这些准则通常可以帮助你获得一个好的解：

　　%26#8226;如果实际可行，采用一个至少可以分析系统的第一阶非零频率的时间步长。

　　%26#8226;把重要的数值阻尼(在TINTP命令中0.05〈P〈1)加到求解中以过滤出高频噪音，特别是如果采用了一个精略的时间步长时，由于阻尼(质量矩阵乘子，ALPHAD命令)会阻碍系统的刚体运动(零频率模态)，在一个动态运动分析中不要使用它。

　　%26#8226;避免强加的位移历程说明，因为强加的位移输入具有(理论上)加速度上的无限突跃，对于Newmark时间积分算法其导致稳定性问题。

　　使用二分法

　　无论何时你打开自动时间步长〔AUTOTS，ON〕，二分法被自动激活。 这个特性通常会使你能够从由于采用一个太大的时间步导致的收敛失败中恢复。它受最小时间步长限制(〔NSUBST，DELTIM〕)。二分法对于任何对加载步长敏感的分析一般是有益的。对于发现一个非线性系统的屈曲临界负载它同样是有用的。

　　使用Newton-Raphson选项和自适应下降因子

　　Newton-Raphson选项的最佳选择将依据存在于你模型中的非线性种类变化。尽管通过让程序选择Newton-Raphson选项〔NROPT，AUTO〕通常你会获得最佳的收敛特性，但也可能偶尔遇到使用一些其它选择会更有效的情况。例如，如果非线性材料的行为发生在你模型的一个相对小的区域中，采用修正的Newton-Raphson或者初始刚度选项可以降低分析的总体CPU代价。自适应下降因子〔NROPT〕和塑性以及某些非线性单元，包括接触单元同时使用。 在几乎没有载荷重新分配的情况下，通过关闭这个特性你可以获得更快的收敛性。自适应下降在仅有大挠度的非线性的问题中几乎没有效果。

　　使用线性搜索

　　线性搜索〔LNSRCH〕作为一个对自适应下降〔NROPT〕的替代会是有用的。(一般地，你不应同时既激活线性搜索又激活自适应下降。)线性搜索方法通常导致收敛，但在时间上它可能是缓慢的和昂贵的(特别是具有塑性时)，在下列情况下你可以设置线搜索为打开状态：

　　%26#8226;当你的结构是力加载的(其与位移控制的相反)时。

　　%26#8226;如果你正在分析一个刚度增长的%26ldquo;薄膜%26rdquo;结构(如一根钓鱼杆)。

　　%26#8226;如果你注意到(从程序的输出信息)你的分析正导致自适应下降频频被激活。

　　应用预测

　　预测〔PRED〕基于基于前一个时间 步的求解预估在这个时间步中的求解情况，因此可能减少所需的平衡迭代次数。如果非线性响应相对地平滑这个特性会是有益的。在大转动和粘弹性分析中它一般不是有益的。

　　应用弧长方法

　　对于许多物理意义上不稳定的结构你可以应用弧长方法〔ARCLEN〕，〔ARCTRM〕来获得数值上稳定的解，当应用弧长方法时，请记住下列考虑事项：

　　%26#8226;弧长方法限制于仅具有渐进加载方式的静态分析。

　　%26#8226;程序由第一个子步的第一次迭代的载荷(或位移)增量计算出参考弧长半径，采用下列公式：

　　参考弧长半径=总体载荷(或位移)÷NSBSTP

　　这里NSBSTP是NSUBST命令中指定的子步数。

　　当选择子步数时，考虑到更多的子步将导致很长的求解时间。理想地，你会选择一个最佳有效解所需的最小子步数。或许你不得不对所需的子步数进行%26ldquo;评诂%26rdquo;，按照需要调整后再重新求解。

　　%26#8226;当弧长方法是激活 的时，不要使用线搜索〔LNSRCH〕，预测〔PRED〕，自适应下降〔NROPT，，，ON〕自动时间分步〔AUTOTS，TIME，DELTIM〕，或时间积分效应〔TIMINT〕。

　　%26#8226;不要尝试将收敛建立在位移的基础上〔CNVTOL，U〕。使用力的收敛准则(CNVTOL，F〕

　　%26#8226;要用弧长方法来帮助使求解时间最小化，一个单一子步中的最大平衡迭代数应当小于或等于15。#p#分页标题#e#

　　%26#8226;如果一个弧长求解在规定的最大迭代次数内〔NEQIT〕没能收敛，程序将自动进行二分且继续分析。直到获得一个收敛的解，或者最小的弧长半径被采用(最小半径由NSUBST〔NSUBST〕和MINARC 〔ARCLEN〕定义)。

　　%26#8226;一般地，你不能应用这种方法来在一个确定的载荷或位移值处获得一个解因为这个值随获得的平衡态改变(沿球面弧)。

　　%26#8226;类似地，当在一个非线性屈曲分析中应用弧长方法来在某些已知的容限范围内确定一个极限载荷或位移的值可能是困难的。通常你不得不通过尝试─错误─再尝试调整参考弧长半径(使用NSUBST)来在极限点处获得一个解。应用带二分〔AUTOTS〕的标准NEWTON-RAPHSON迭代来确定非线性载荷屈曲临界负载的值可能会更方便。

　　%26#8226;通常你应当避免和弧长方法一起使用JCG或者PCG求解器〔EQSLV〕，因为弧长方法可能会产生一个负定刚度矩阵(负的主对角线)，用这些求解器其可能导致求解失败。

　　%26#8226;在任何载荷步的开始你可以从Newton-Raphson迭代方法到弧长方法自由转换。然而，要从弧长到Newton-Raphson迭代转换，你必须终止分析然后重起动，且在重起动的第一个载荷步中去杀死弧长方法〔ARCLEN，OFF〕。

　　一个弧长求解在这些情况下终止：

　　%26#8226;当由ARCTRM或NCNV命令定义的极限达到时。

　　%26#8226;当在所施加的载荷范围内求解收敛时。

　　%26#8226;当你使用一个放弃文件时(Jobname.ABT)。

　　%26#8226;使用载荷位一移曲线作为用于评价和调整你的分析以帮助你获得所需结果的准则。通常对于每一个分析都绘 制你的载荷一偏移曲线(采用POST26命令)是一种好的作法。

　　%26#8226;经常地，一个不成功的弧长分析可以归因于弧长半径或 者太大或 者太小。 沿载荷一偏移曲线原路返回的%26ldquo;回漂%26rdquo;是一种由于使用太大或 太小弧长半径导致的典型难点。研究载荷偏移曲线来理解这个问题。然后使用NSUBST和ARCLEN命令来调整弧长半径的大小和范围为合适的值。

　　%26#8226;总体弧长载荷因子(SOLU命令中的ALLF项)或者会是正的或者会是负的。类似地，TIME，其在弧长分析中相关于总体弧长载荷因数，同样会不是正的就是负的。ALLF或TIME的负值表示弧长特性正在以反方向加载，以便保持结构中的稳定性。负的ALLF或者TIME值一般会在各种突然转换分析中遇到。

　　%26#8226;当将弧长结果读入基本数据用于POSTI后处理时〔SET〕，你总是应当引用由它的载荷步和子步号〔LSTEP和SBSTEP〕或者进它的数据设置号所设定的所需结果数据。不要引用用TIME值的结果，因为TIME值在一个弧长分析中并不总是单调增加的。(单一的一个TIME值可能涉及多于一个的解。)此外，程序不能正确地解释负的TIME值(C其可能在一个突然转换分析中遇到。)

　　%26#8226;如果TIME为负的，记住在产生任何POST26图形前定义一个合适的变化范围(〔IXRANGE〕或者〔IYRANGE〕)。

　　在你的模型响应中人为地抑制发散

　　如果你不想使用弧长方法来分析一个在奇异(零刚度)形状时开始开，或者通过奇异形状的力加载的结构时，有时你可以使用其它的技术来人工地抑制模型响应中的发散。

　　%26#8226;在某些情况下，你可以使用强加的位移来替代所施加的力。这种方法可以用于在较靠近平衡位置处开始一个静态分析，或者用于控制整个不稳定响应期间(如突然转换或后翘曲)的位移。

　　%26#8226;其它在阻止由于初始不稳定性所造成的问题时有效的技术包括：使用带有强加的初始应变的应力刚化〔SSTIF〕，%26ldquo;致冷%26rdquo;(也就是，增加暂时的人工热应变)，或者将一个静态问题执行为一个%26ldquo;缓慢动态%26rdquo;分析(也就是，在任意一个载荷步尝试使用时间积分效应阻止解发散。

　　%26#8226;你也可以应用控制单元(如COMBIN37)，或者应用其它单元的出生和死亡选项对不稳定的DOFs施加暂时的人工刚度。这里的想法是在中期的载荷步期间人为地约束系统，以阻止不符合实际的大位移被计算出。随着系统变位到稳定的形态，人工刚度被移去。

　　应用雅各比共轭梯度求解器

　　这个求解器(通过EQSLV命令获得)在经历某一奇异划(零 (零刚度)状态的分析中会是有用的。叶%26Ocirc;JCG求解器来说相对大的求解容差有时会%26ldquo; 涂抹掉%26rdquo;这种奇异性，导致载荷一位移曲线的斜度具有某些假的非零值。(在EQSLV中这个求解器的容限不是非线性收敛容限。)

　　雅各比共轭梯度求解器仅是一种求解线性矩阵方程的替代方法。这种求解器的使用不能替代任何方式的非线性处理。

　　关闭特殊的单元形状

　　有时在非线性分析中使用无中节点单元的形状选项会产生收敛困难。

　　合理地使用出生和死亡

　　认识到结构的刚度矩阵的任何突然改变可能会导致收敛问题。当激活或杀死单元时，试着将变化分散在若干子步内。(如果需要，采用一个小的时间步长来完成这种变化。)也要注意到随着你激活或杀死单元可能会产生的奇异性(如尖的再生角)。像这样的奇异性可能产生收敛问题。#p#分页标题#e#

　　检验你的分析结果

　　好的有限无分析(FEA)过程总是要求你检验你的结果。你需要自己证明你理解了程序，你正在正确地使用它，以及你的分析结果正确地体现出你的结构的物理特性。在检验你的非线性分析时你可以使用若干标准验证技术。

　　标准分析

　　一个确保你了解如何恰当地施加程序的特殊特性的好的方法是通过进行一个或多个标准分析。在一个标准分析中，一般是你对一个有%26ldquo;理论%26rdquo;解存在的简单结构进行独立地分析。这里的想法是通过将你的FEA结果与已知结果相对照以验证你可以正确地运用程序的特性。当然，标准分析结构应当与要分析的完整结构非常相似。ANSYS Verification Manual是标准问题的一种较好的来源。

　　结果合理么

　　大多数工程师在他们职业的早期就认识到要对他们的数值结果的有效性提出疑问，无论这些结果是通过%26ldquo;手工%26rdquo;计算，计算机分析，还是一些其它方法得到的。在你开始任何分析前，你总是应当对你期望获得的结果至少具有一个粗略的概念(通过经验、试验、标准分析等等获得)。如果你最终的结果似乎不合理，也就是，如果它们不同于你的期望值，你应当确信你理解了这是为什么。好的工程实际要求你总是使你的分析结果和合理的期望值相一致。

　　理解你的输出

　　记住ANSYS程序将一个非线性分析作为一系列带修正的线性近似来完成。程序的打印输出给出你关于这些近似和修正发展的连续反馈。(打印输出或者直接出现在你的屏幕上，记录在Jobname.OUT中，或者被写入某些其它人文件〔OUTPUT〕。)你可以在POST中应用PRITER命令，或者在POST26中应用SOLU和PRVAR命令检查这种类似的信息。在你接受结果前，你应当确信你理解了你的分析的迭代历程。特别地，不要忽视任何还没有完全理解它们意思的程序错误和警告声明。

　　作载荷和响应历程的曲线图

　　这种检验技巧可以认为是两种其它技巧的图形结合：对合理性的检查和考察迭代历程。载荷和响应历程的POST26图形表示应当和你所知道的你结构特性的期望值相一致。重要的结果(位移，反作用力，应力，等等)应当显示出相对平滑的响应历程。任何非平滑性可能表示采用了一个太粗略的时间步。