《ANSYS18.0 有限元分析基础与实例教程》课件——第14章 非线性分析.pptxVIP

第14章非线性分析;14.1非线性分析概论;14.1.1非线性行为的原因;（3）材料非线性。非线性的应力-应变关系是造成结构非线性的常见原因。许多因素可以影响材料的应力-应变性质，包括加载历史（如在弹-塑性响应状况下），环境状况（如温度），加载的时间总量（如在蠕变响应状况下）。

ANSYS程序的方程求解器计算一系列的联立线性方程来预测工程系统的响应。然而，非线性结构的行为不能直接用这样一系列的线性方程表示。需要一系列的带校正的线性近似来求解非线性问题。;1．非线性求解方法

一种近似的非线性求解是将载荷分成一系列的载荷增量。可以在几个载荷步内或者在一个载荷步的几个子步内施加载荷增量。在每一个增量的求解完成后，继续进行下一个载荷增量之前程序调整刚度矩阵以反映结构刚度的非线性变化。遗憾的是，纯粹的增量近似不可避免地随着每一个载荷增量积累误差，导至结果最终失去平衡，如图14-3a)所示。

ANSYS程序通过使用牛顿－拉普森平衡迭代克服了这种困难，它迫使在每一个载荷增量的末端解达到平衡收敛（在某个容限范围内）。图14-3b)描述了在单自由度非线性分析中牛顿－拉普森平衡迭代的使用。在每次求解前，NR方法估算出残差矢量，这个矢量是回复力（对应于单元应力的载荷）和所加载荷的差值。程序然后使用非平衡载荷进行线性求解，且核查收敛性。如果不满足收敛准则，重新估算非平衡载荷，修改刚度矩阵，获得新解。持续这种迭代过程直到问题收敛。

对某些物理意义上不稳定系统的非线性静态分析，如果仅仅使用NR方法，正切刚度矩阵可能变为降秩矩阵，导致严重的收敛问题。这样的情况包括独立实体从固定表面分离的静态接触分析，结构或者完全崩溃或者“突然变成”另一个稳定形状的非线性弯曲问题。对这样的;;2．非线性求解级别

非线性求解被分成3个操作级别：

（1）“顶层”级别由在一定“时间”范围内你明确定义的载荷步组成。假定载荷在载荷步内是线性地变化的。

（2）在每一个载荷子步内，为了逐步加载可以控制程序来执行多次求解（子步或时间步）。

（3）在每一个子步内，程序将进行一系列的平衡迭代以获得收敛的解。

图14-5说明了一段用于非线性分析的典型的载荷历史。

3．载荷和位移的方向改变

当结构经历大变形时应该考虑到载荷将发生了什么变化。在许多情况中，无论结构如何变形施加在系统中的载荷保持恒定的方向。而在另一些情况中，力将改变方向，随着单元方向的改变而变化。

ANSYS程序对这两种情况都可以建模，依赖于所施加的载荷类型。加速度和集中力将不管单元方向的改变而保持它们最初的方向，表面载荷作用在变形单元表面的法向，且可被用来模拟“跟随”力。图14-6说明了恒力和跟随力。

4．非线性瞬态过程分析

非线性瞬态过程的分析与线性静态或准静态分析类似：以步进增量加载，;程序在每一步中进行平衡迭代。静态和瞬态处理的主要不同是在瞬态过程分析中要激活时间积分效应。因此，在瞬态过程分析中“时间”总是表示实际的时序。自动时间分步和二等分特点同样也适用于瞬态过程分析。;小转动（小挠度）和小应变通常假定变形足够小以至于可以不考虑由变形导致的刚度阵变化，但是大变形分析中，必须考虑由于单元形状或者方向导致的刚度阵变化。使用命令NLGEOM,ON（GUI：MainMenuSolutionAnalysisTypeSolnControl(:BasicTab)或者MainMenuSolutionUnabridgedMenuAnalysisTypeAnalysisOptions)，可以激活大变形效应（针对支持大变形的单元）。对于大多数实体单元（包括所有大变形单元和超弹单元）和大多数梁单元和壳单元都支持大变形。

大变形过程在理论上并没有限制单元的变形或者转动（实际的单元还是要受到经验变形的约束，即不能无限大），但求解过程必须保证应变增量满足精度要求，即总体载荷要被划分为很多小步来加载。

1．大应变大挠度（大转动）

所有梁单元和大多数壳单元以及其他的非线性单元都有大挠度（大转动）效应，可以通过命令NLGEOM,ON（GUI：MainMenuSolutionAnalysisTypeSolnControl(:BasicTab)或者MainMenuSolutionUnabridgedMenuAnalysisTypeAnalysisOptions）来激活该选项。;2．应力刚化

结构的面外刚度有时候会受到面内应力的明显影响，这种面内应力与面外刚度的耦合，即所谓的应力刚化，在面内应力很大的薄结构（例如缆索、隔膜）中非常明显。

因为应力刚化理论通