第四章 材科非线性限元法.docVIP

第四章 材科非线性有限元法 本章所研究的非线性有限元方法是材料非线性有限元。。，，，，。 4.1.1 非线性弹性有限元法的基本公式 对于小变形的非线性弹性问题，: 平衡微分方程 (3-8) 几何方程 (3-2) 物理方程 (3-19) 或 (3-20) 式中，的表达式见（3-37）~（3-40），的函数。 位移边界条件： (3-4) 应力边界条件： (3-9) 虚功方程是 (3-11) 可见，，，3-11）不涉及材料性质， 元和整体平衡方程完全适用于非线性弹住和弹塑性问题， (4-1) (4-2) (4-3) (4-4) 将式（4-1）代入（3-2）， (4-5) 其中单元割线劲度矩阵为 (4-6) 同样可由（4-4）得整体平衡方程 (4-7) 整体割线劲度矩阵 (4-8) 由于与位移有关，4-7）是一个非线性方程组。4-7）式，4-7）要用直接迭代法，， 在求解非线性方程组时，，。， 由式（4-4） (4-9) 由上式及Newton法迭代公式（2-7） (4-10) 4.1.2 求解的选代过程 不同的非线性方程组求解方法，， 处理一般都分级按增量方法计算。 1. 荷载分级 对于实际的工程问题来说，，， 施工加载过程，仿真的数值分析，， 载过程将得到不同的位移和应力计算成果。其次，， 荷载的分级首先要考虑荷载的性质， 加。为了更精细地模拟施工过程并使迭代过程收敛，2.2介绍的荷载系数 法分成更小的增量。图4-1为三峡永久船闸的闸室横剖面，，，，，，，，。4-1 2. 增量迭代法 增量迭代法实际上就是非线性方程组求解的Euler-Newton法， 增量，，m级荷载， (2-46) 式中 一一第m级荷载增量施加后的总荷载， ——第m级荷载增量， ——第m级荷载第i次迭代结束时的结点力，按下式计算 (4-11) ——第m级荷载第i次迭代的不平衡力 切线劲度矩阵根据式（4-9 (4-12) 如果已知第m级荷载增量时第i次迭代的近似解，相应的应变为 (4-13) 于是切线弹性矩阵可以确定，4-12）和（4-11）分别计算和，最后利用式（2-46）求出，这个迭代过程由i=0开始直至足够小，i=I， (4-14) 按（2-46）后， (4-15) 第m级荷载增量的最终应力为 (4-16) 一维问题的迭代过程见图4-2。 可以看出，， 2．3初应力法 如果在弹性材料内确实存在初应力， (4-17) 由上式及虚功原理可导出单元的结点力为 (4-18) 集合单元得出以下的有限元支配方程 (4-19) 式中，引起的等效结点荷载 图4-2 (4-20) 对于非线性问题，，（4-19，4-17）表示，并非真正的初应力，和弹性应力的差，4-3可以看出，，（）看作是变化的，， 前面所谈的增量迭代法，并完整地求解一次线性方程组。如果将式（2-46）中的改用初始的切线劲度矩阵或，Euler修正Newton法，，可以认为是第m级荷载增量开始施加时的“弹性”劲度矩阵 (4-21) 式中的“弹性”矩阵实际上是。式（2-46） (4-22) 根据式（4-21）可得出下面的恒等式