有限元七材料非线性问题有限元.docVIP

第七章　材料非线性问题的有限元法 7.1　引　　言 前面各章所讲述的问题，都属于线性变形体系。所谓线性变形体系是指位移与载荷呈线性关系的体系，而且当载荷全部撤除后，体系将完全恢复原始状态。这种体系也称为线性弹性体系，它需满足下列条件： （1）材料的应力与应变关系满足虎克定律： （2）位移是微小的； （3）所有约束均为理想约束。 在分析线性弹性体系时，可以按照体系变形前的几何位置和形状建立平衡方程，并且可以应用叠加原理。根据这种理论建立起来的方程是线性的，对于小应变和小位移的情形这种分析是适用的。 实际结构的位移与载荷可以不呈线性关系，这样的体系称为非线性变形体系。如果体系的非线性是由于材料应力与应变关系的非线性引起的，则称为材料非线性，如材料的弹塑性性质、松驰、徐变等。如果结构的变位使体系的受力发生了显著的变化，以至不能采用线性体系的分析方法时就称为几何非线性，如结构的大变形、大挠度的问题等。还有一类非线性问题是边界条件非线性，或状态非线性，如各种接触问题等。但本书只讨论前两类非线性问题的有限元解法，即材料非线性和几何非线性问题的有限元解法，对接触问题的有限元解法，读者可参考其它书籍。 材料非线性问题的处理相对比较简单，通常不必修改整个问题的表达式，而只需将应力—应变关系线性化，求解一系列的线性问题，并通过某种校正方法，最终将材料特性调整到满足给定的木构关系，从而获得了问题的解。 对于几何非线性问题，那就需要对公式进行根本的修改，这个问题将在后面详细讨论，不过应该指出，用于求解材料非线性问题的基本迭代方法也同样适用于几何非线性问题的求解。事实上，有些工程结构问题同时具有这两类非线性性质，它们可以统一地加以处理。 本章将首先介绍用有限元方法处理非线问题的一般方法，然后讨论这些方法在非线性弹性、弹塑性和蠕变问题中的应用。在介绍弹塑性问题的处理方法前，为便于讨论，需扼要叙述一下Mises屈服准则和Prandtl-Reuss塑性流动理论，并据此写出弹塑性矩阵表达式。最后对平面刚架的极限分析做了简要介绍。 7.2　非线性问题的一般处理方法 非线性问题用有限元法离散化应得到如下形式的一组代数方法： 或写成 （7.1） 其中　。虽然线性方程组直接求解并无困难，但对于方程组（7.1），单元刚变矩阵是单元节点位移向量的函数，直接求解就行不通。然而，下面介绍的非线性方程组的各种解法，仍以反复地求解线性方程组去获得满足一定精度要求的非线性方程组的解答。 7.2.1　直接迭代法 对于方程（7.1） （7.2） 最简单的求解方法是直接迭代法。开始求解时先假定一组初始值代入上式的中，可求得改进了的一次近似值 式中 重复上述过程，将迭代格式写成 迭代一直进行到误差的某种范数小于预先规定的容许值er，即满足 则停止。 可以看出，该法的每一次迭代都需形成一次系数矩阵，并求解一次线性代数方程组。这里还隐含着一个假定，系数矩阵可以表示成的显函数，因此该法只适用于与变形历史无关的非线性问题，例如非线性弹性问题及可利用形变理论分析的弹塑性问题。而对于依赖于变形历史的非线性问题，直接迭代法是不适用的，例如加载路径不断变化或涉及卸载及反复加载等必须利用增量理论分析的弹塑性问题。 图7.1显示了单变量问题中这种迭代过程收敛和发散的可能性。通常，如曲线是凹的，则迭代发散。 图7.1 直接迭代法 7.2.2　Newton-Raphson方法（简称N－R方法） 若已获得方程（7.1）的第n次近似解，为了求得改进的近似解，可利用仅保留线性项的Taylor级数展开式 （7.3） 则有 （7.4） 前式中　为切线矩阵，即 （7.5） 于是从式（7.3）可以得到式（7.4）中的 （7.6） 重复上述迭代过程，直至达到所要求的精度。 Newton-Raphson方法的迭代过程示于图7.2中，通常迭代过程是收敛的。但当所选取的初始值偏离真实解较大时，正如图7.2（b）所表示的那样，发散也是可能的。 图7.2　Newton-Raphson迭代法 由式（7.3）看出，该法在每项迭代中必须重新形成一个系数阵并求解一次线性代数方程组。应该指出，如果原始的离散化方程组是通过变分原理导出的，则切线刚度矩阵总是对称的，而利用直接迭代法，系数矩阵的这种对称性不一定能保持。 7.2.3　修正的Newton-Raphson方法（简称修正的N－R法） 为克服Newt