ANSYS结构非线性分析指南.doc

第一章? 结构非线性分析概述 1.1? 什么是结构非线性 ??? 在日常生活中，经常会遇到结构非线性。例如，当用钉书针钉书时，金属钉书钉将永久地弯曲成一个不同的形状( 图1-1a )。如果你在一个木架上放置重物，随着时间的推移木架将越来越下垂( 图1-1b )。当在汽车或卡车上装载货物时，它的轮胎和下面路面间接触面将随货物重量而变化( 图1-1c )。如果将上述例子的载荷变形曲线画出来，用户将发现它们都显示了非线性结构的基本特征—结构刚度改变。 图1-1? 结构非线性行为的常见例子 ??? 引起结构非线性的原因很多，它可以被分成三种主要类型：状态改变、几何非线性、材料非线性。 1.1.1? 状态变化(包括接触) ??? 许多普通结构表现出一种与状态相关的非线性行为。例如，一根只能拉伸的电缆可能是松的，也可能是绷紧的。轴承套可能是接触的，也可能是不接触的。冻土可能是冻结的，也可能是融化的。这些系统的刚度由于系统状态的改变而变化。状态改变也许和载荷直接有关(如在电缆情况中)， 也可能由某种外部原因引起(如在冻土中的紊乱热力学条件)。 ??? 接触是一种很普遍的非线性行为。接触是状态变化非线性中一个特殊而重要的子集。参见第五章。 1.1.2? 几何非线性 ??? 如果结构经受大变形，它几何形状的变化可能会引起结构的非线性响应。一个例子是 图1-2 所示的钓鱼杆。随着垂向载荷的增加，杆不断弯曲以致于力臂明显地减少，导致杆端显示出在较高载荷下不断增长的刚性。几何非线性的特点是大位移、大转动。 图1-2? 钓鱼杆体现的几何非线性 1.1.3? 材料非线性 ??? 非线性的应力─应变关系是结构非线性行为的常见原因。许多因素可以影响材料的应力─应变性质，包括加载历史(如在弹─塑性响应情况下)、环境状况(如温度)、加载的时间总量(如在蠕变响应情况下)。 1.2? 非线性分析的基本知识 1.2.1? 方程求解 ??? ANSYS程序的方程求解器计算一系列的联立线性方程来预测工程系统的响应。然而，非线性结构的行为不能直接用这样一系列的线性方程表示。需要一系列的带校正的线性近似来求解非线性问题。 一种近似的非线性求解是将载荷分成一系列的载荷增量。可以在几个载荷步内或者在一个载荷步的几个子步内施加载荷增量。在每一个增量的求解完成后，继续进行下一个载荷增量之前程序调整刚度矩阵以反映结构刚度的非线性变化。但是，纯粹的增量近似不可避免地要随着每一个载荷增量积累误差，导种结果最终失去平衡，如 图1 - 3 （a） 所示。 . （b）纯粹增量式解???????? ?????（b）全牛顿－拉普森迭代求解（2个载荷增量） 图1-3??? 纯粹增量近似与牛顿－拉普森近似。 ANSYS程序通过使用牛顿－拉普森平衡迭代克服了这种困难，它迫使在每一个载荷增量的末端解达到平衡收敛（在某个容限范围内）。 图1 - 3 （b） 描述了在单自由度非线性分析中牛顿－拉普森平衡迭代的使用。在每次求解前，NR方法估算出残差矢量，这个矢量是回复力（对应于单元应力的载荷）和所加载荷的差值，然后使用非平衡载荷进行线性求解，且核查收敛性。如果不满足收敛准则，重新估算非平衡载荷，修改刚度矩阵，获得新解。持续这种迭代过程直到问题收敛。仅在非线性分析中使用这个选项。这个选项指定在求解期间每隔多久修改一次正切矩阵。你可以指定这些值中的一个。?程序选择（NROPT，ANTO）：程序基于你模型中存在的非线性种类选择用这些选项中的一个。在需要时牛顿－拉普森方法将自动激活自适应下降。?全〔NROPT，FNLL〕；程序使用完全的牛顿－拉普森处理方法，在这种处理方法中每进行一次平衡迭代修改刚度矩阵一次。如果自适应下降是关闭的，程序每一次平衡迭代都使用正切刚度矩阵。（我们一般不建议关闭自适应下降，但是你或许发现这样做可能更有效。）如果自适应下降是打开的（缺省），只要迭代保持稳定（也就是，只要残余项减小，且没有负主对角线出现）程序将仅使用正切刚度阵。如果在一次迭代中探测到发散倾向，程序抛弃发散的迭代且重新开始求解，应用正切和正割刚度矩阵的加权组合。当迭代回到收敛模式时，程序将重新开始使用正切刚度矩阵。对复杂的非线性问题自适应下降通常将提高程序获得收敛的能力。?修正的（NROPT，MODI）：程序使用修正的牛顿－拉普森方法，在这种方法中正切刚度矩阵在每一子步中都被修正。在一个子步的平衡迭 代期间矩阵不被改变。这个选项不适用于大变形分析。自适应下降是不可用的?初始刚度（NROPT，INIT）：程序在每一次平衡迭代中都使用初始刚度矩阵这一选项比完全选项似乎较不易发散，但它经常要求更多次的迭代来得到收敛。它不适用于大变形分析。自适应下降是不可用的。 ??? ANSYS程序提供了一系列命令来增强问题的收敛性，如自适应下