非线性有限元方法及实例分析

《非线性有限元方法及实例分析》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、http://www.paper.edu.cn非线性有限元方法及实例分析梁军河海大学水利水电工程学院，南京（210098）摘要：对在地下工程稳定性分析中常用的非线性方程组的求解方法进行研究，讨论了非线性计算的迭代收敛准则，并利用非线性有限元方法分析了一个钢棒单轴拉伸的实例。关键词：非线性有限元，方程组求解，实例分析1引言有限单元法已成为一种强有力的数值解法来解决工程中遇到的大量问题，其应用范围从固体到流体，从静力到动力，从力学问题到非力学问题。有限元的线性分析已经设计工具被广泛采用。但对于绝大多数水利工程中遇到的实际问题如地下洞室等，将其作为非线性问题[1]加以考虑更符

2、合实际情况。根据产生非线性的原因，非线性问题主要有3种类型：1．材料非线性问题（简称材料非线性或物理非线性）2．几何非线性问题3．接触非线性问题（简称接触非线性或边界非线性）2非线性方程组的求解在非线性力学中，无论是哪一类非线性问题，经过有限元离散后，它们都归结为求解一个非线性代数方程组[2]：ψ()δδΛδ=0112nψ()δδΛδ=0212n……ψ()δδΛδ=0n12n（1.1）δ,δ,Λ,δψ,ψ,Λ,ψδ,δ,Λ,δ其中12n是未知量，12n是12n的非线性函数，引用矢量记号[]Tδ＝δδΛδ12n（1.2）[]Tψ=ψψΛψ12n（1.3）上述方程组（1.1

3、）可表示为ψ()δ=0（1.4）可以将它改写为ψ()δ≡F()δ−R≡K()δδ−R=0（1.5）其中K()δ是一个n×n的矩阵，其元素kij是矢量的函数，R为已知矢量。在位移有限F(δ)ψ()δ=0元中，δ代表未知的结点位移，是等效结点力，R为等效结点荷载，方程表示结点平衡方程。在线弹性有限元中，线性方程组-1-http://www.paper.edu.cnKδ－R＝0（1.6）ψ(δ)=0可以毫无困难地求解，但对线性方程组则不行。一般来说，难以求得其精确解，通常采用数值解法，把非线性问题转化为一系列线性问题。为了使这一系列线性解收敛于非线性解，曾经有过许多方法，但这

4、些解法都有一定的局限性。某解法对某一类非线性问题有效，但对另一类问题可能不合适。因而，根据问题性质正确选用求解方法成为非线性有限元的一个极重要的问题。常见的求解非线性方程组的数值方法有迭代法、增量法和混合法[3][4][5]。2.1.迭代法在每次迭代过程中都施加全部荷载，但逐步修改位移和应变，使之满足非线性的应力－应变关系。迭代法还分为直接迭代法、Newton-Raphson方法、修正的Newton-Raphson方法和拟Newton法。①直接迭代法对非线性方程组K()δδ－R＝0（1.7）0设其初始的近似解为δ＝δ，由此确定近似的K矩阵0(0)K＝Kδ（1.8）可得改

5、进的近似解1()0−1δ＝KR（1.9）重复这一过程，以第i次近似解求出第i+1次近似解的迭代公式为i(i)K＝Kδ（1.10）i+1()i−1δ＝KR（1.11）ii+1i直到Δδ=δ−δ变得充分小，即近似收敛时，终止迭代。对于单变量问题，这一迭代过程是收敛的，单对于多自由度情况，由于未知量通过矩阵K耦合，迭代过程可能不收敛。K(δ)δ－R＝0在迭代过程中，得到的近似解一般不会满足即(i)(i)iψδ≡Kδδ−R≠0（1.12）ψ()δ作为对平衡偏离的一种度量，称为失衡力。②Newton-Raphson方法Newton-Raphson方法是求解非线性方程组-2-htt

6、p://www.paper.edu.cnψ()δ≡F()δ−R≡K()δδ−R=0（1.13）的一个著名方法，简称Newton法。ψ()δδ＝δi设为具有一阶导数的连续函数，是方程（1.13）的第i次近似解。若i(i)()ψ=ψδ≡Fδi−R≠0（1.14）希望能找到一个更好的、方程（1.13）的近似解为i+1iiδ＝δ=δ+Δδ（1.14）ii将（1.15）代入（1.13），并在δ＝δ附近按一阶Tayor级数展开，则ψ()δ在δ处的线性近似公式为ii+1i⎛∂ψ⎞iψ=ψ+⎜⎟Δδ⎝∂δ⎠（1.16）iii⎛∂ψ⎞KT=KT()δ≡⎜⎟引入记号⎝∂δ⎠Newton法的

7、迭代公式可归纳为i()i−1i()(i−1iΔδ=−Kψ=KR−F)TT（1.17）iii⎛∂ψ⎞⎛∂F⎞KT=⎜⎟=⎜⎟⎝∂δ⎠⎝∂δ⎠（1.18）i+1iiδ=δ+Δδ（1.19）Newton法的收敛性是比较好的，但对于某些非线性问题，如理想塑性和塑性软化问题，KK在迭代过程中T可能事奇异或病态的，于是T的求逆就会出现困难。为此，可引入一个iiηK+ηI阻尼因子，使矩阵T或者成为非奇异的，或者使它的病态减弱。③修正的Newton-Raphson方法iK上述两种方法求解非线性方程组时，在迭代过程中的每一步都需要重新计算T。如将i00KK