ansys求解非线性问题

《ansys求解非线性问题》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、ANSYS求解非线性问题  牛顿一拉森方法ANSYS程序的方程求解器计算一系列的联立线性方程来预测工程系统的响应。然而，非线性结构的行为不能直接用这样一系列的线性方程表示。需要一系列的带校正的线性近似来求解非线性问题。逐步递增载荷和平衡迭代一种近似的非线性救求解是将载荷分成一系列的载荷增量。可以在几个载荷步内或者在一个载步的几个子步内施加载荷增量。在每一个增量的求解完成后，继续进行下一个载荷增量之前程序调整刚度矩阵以反映结构刚度的非线性变化。遗憾的是，纯粹的增量近似不可避免地随着每一个载荷增量积累误差，导种结

2、果最终失去平衡，如图1所示所示。                            （a） 纯粹增量式解            （b）全牛顿－拉普森迭代求解                                   图1   纯粹增量近似与牛顿－拉普森近似的关系ANSYS程序通过使用牛顿－拉普森平衡迭代克服了这种困难，它迫使在每一个载荷增量的末端解达到平衡收敛（在某个容限范围内）。图1（b）描述了在单自由度非线性分析中牛顿－拉普森平衡迭代的使用。在每次求解前，NR方法估算出残差矢量，这个矢量是回

3、复力（对应于单元应力的载荷）和所加载荷的差值。程序然后使用非平衡载荷进行线性求解，且核查收敛性。如果不满足收敛准则，重新估算非平衡载荷，修改刚度矩阵，获得新解。持续这种迭代过程直到问题收敛。ANSYS程序提供了一系列命令来增强问题的收敛性，如自适应下降，线性搜索，自动载荷步，及二分等，可被激活来加强问题的收敛性，如果不能得到收敛，那么程序或者继续计算下一个载荷前或者终止（依据你的指示）。对某些物理意义上不稳定系统的非线性静态分析，如果你仅仅使用NR方法，正切刚度矩阵可能变为降秩短阵，导致严重的收敛问题。这样的

4、情况包括独立实体从固定表面分离的静态接触分析，结构或者完全崩溃或者“突然变成”另一个稳定形状的非线性弯曲问题。对这样的情况，你可以激活另外一种迭代方法，弧长方法，来帮助稳定求解。弧长方法导致NR平衡迭代沿一段弧收敛，从而即使当正切刚度矩阵的倾斜为零或负值时，也往往阻止发散。这种迭代方法以图形表示在图2中。                   非线性问题（一）材料特性可是是线性或非线性、各向同性或正交各向异性，常数的或温度相关的。但必须遵循如下规则：（1）必须定义材料刚度（如弹性模量EX、超弹性系数等）（2）对

5、于惯性载荷（如重力、加速度等），必须定义材料的质量（如密度DENS）。对于温度载荷，必须定义热膨胀系数APLX。（3）对于网格密度。原则上网格密度越大，单元越多，结果就更符合实际，但是由于网格密度的提高会引起结构自由度的增大，直接影响求解的耗时，因此建议载应力或应变急剧变化的区域（通常也是用户感兴趣的区域），细分网格。考虑非线性因素时，要有足够的网格来得到非线性效应，例如在高塑性变形梯度区设置较密的网格来保证塑性分析所需的足够积分点密度。非线性有限元问题与线性有限元问题有很大的不同，主要有以下几个方面：l  

6、       非线性问题的方程是非线性的，因此一般需要进行迭代求解l         非线性问题不能采用叠加原理l         非线性问题不一定总有一致解。尽管问题的定义都是正确的，但可能会没有解。引起非线性的原因很多，大致可以分为3类：l         几何非线性l         材料非线性l         状态非线性1.       几何非线性：是由位移之间存在的非线性关系引起的。如果结构经受大变形，它变化的几何形状可能会引起结构的非线性响应。2.       材料非线性：是由应力应变关系引起的

7、。这些非线性关系不能单靠数学得到，而要基于试验数据。非线性材料特性有时可用数学模型模拟。3.       状态变化（包括接触）：由边界条件或载荷所引起的，许多普通结构标表现出一种与状态相关的非线性行为。接触是一种很普遍的非线性行为，也是状态变化中非线性类型中一个特殊而重要的子集。非线性问题（二）  非线性问题的求解ANSYS程序的方程求解器通过计算一系列的联立线性方程来预测工程系统的响应。然而，非线性结构的行为不能直接用这样一系列的线性方程表示，需要一系列的带校正的线性来近似求解非线性问题。1）  增量法：增

8、量方法是一种近似的非线性求解，将载荷分成一系列的载荷增量。可以在几个载荷步内或在一个载荷步的几个子步内施加载荷增量。在每一个增量求解完成后，继续进行下一个载荷增量之前，程序调整刚度矩阵以反映结构刚度的非线性变化。但是纯粹的增量近似不可避免的随着每一个载荷增量积累误差，从而导致解的漂移，而且随着增量数目的增加，这种漂移现象将会越来越严重，最终导致结果失去平衡。2）  New－Raphson方法：ANS