资源预览内容
第1页 / 共9页
第2页 / 共9页
第3页 / 共9页
第4页 / 共9页
第5页 / 共9页
第6页 / 共9页
第7页 / 共9页
第8页 / 共9页
第9页 / 共9页
亲,该文档总共9页全部预览完了,如果喜欢就下载吧!
资源描述
非线性有限元方程组的解法不论材料非线性问题还是几何非线性问题,其有限元 方程都是非线性的:R-外部载荷的等效节点力矢量; P-内力的等效节点力矢量。 在全量方法的位移有限元解法中,u是结构的位移矢量 ,在增量方法的位移有限元解法中, u是结构的位移 增量矢量。非线性有限元方程组的解法 对于线弹性小变形问题,其有限元方程组是线性的 其解答利用直接方法很容易得到 但是对于非线性有限元方程组则不能利用直接方法 得到其解答。 一般地说,不能期望得到非线性方程组的精确界精确界。 通常利用各种数值方法,用一系列一系列的线性方程组去 逼近非线性方程组的解。非线性有限元方程组的解法 直接迭代法 牛顿法 增量法非线性有限元方程组的解法(直接迭代法)固体力学中非线性有限元方程通常可以写成:其中设初始近似解为 ,那么可得一个近似矩阵 于是由(1)可得到一个改进的近似解:据此容易写出直接迭代法的迭代公式:按照这种迭代公式可以得到一个解数列 ,当这个数列收敛时停止计 算,其数列收敛值就是方程(1)的解。非线性有限元方程组的解法(直接迭代法)关于数列收敛的判据,可以采用各种各样的范数定义和收敛判据。若设第 n 次和第 n+1 的解分别为 、 ,则“偏差”为:范数的定义可取 或于是收敛判据可取为: (位移收敛判据位移收敛判据) 在这里注意到,对于非线性方程(1),将 代入一般不是严格满足的,即该值可作为对偏离平衡偏离平衡的一种度量(称为失衡力),收敛判据可相应地取为 :(失衡力收敛判据)收敛判据)非线性有限元方程组的解法(牛顿法)把非线性有限元方程记为:现在设 是方程(1)的第 n 次近似解。一般地,这时该值可作为对偏离平衡的一种度量(称为失衡力失衡力)。设修正值为 , 此时新的近似解为:将(3)代入(1)中并在 附近将 泰勒(Taylor)展开:若记可得 从而可解出修正量 为:非线性有限元方程组的解法(牛顿法)这样,牛顿法的迭代公式:其收敛判据与直接迭代法的直接迭代法的收敛判据雷同。非线性有限元方程组的解法(增量法)求解非线性方程组的另一类方法是增量方法。使用增量方法的一个优点是 可以得到整个载荷变化过程的一些中间的数值结果。当问题的性质与加载的加载的 历史历史有关时,例如弹塑性问题,则必须采用增量方法。 在增量方法中通常引入载荷因子,用 表示载荷, 于是非线性有限 元方程可写成:用载荷因子系列:相应于不同的载荷。 若相应于载荷因子 的解已经求得,记为 ,则现在来求相应于载荷因子为 时的解。 设 为其解, 于是有将 在 处泰勒展开得非线性有限元方程组的解法(增量法)若记作:考虑到 ,于是方程(5)可近似为或若考虑到相应于载荷因子 的解 并不是精确解,亦即:于是方程的解为
收藏 下载该资源
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号