第三章 非线性有限元理论.ppt

G 为如下9×3m的矩阵 第五节 几何非线性问题的有限元法 式中：(A)L为如下6×9的矩阵 第五节 几何非线性问题的有限元法 由(A)L可见，格林应变-位移关系是非线性的。 非线性部分 所以 为用虚位移原理建立单元特性方程，还得建立应变增量和 位移增量间的关系。对线性部分 综上所述，格林应变增量为 如果记 则 第五节 几何非线性问题的有限元法 对弹性问题，在物质坐标系描述下本构关系为 在上述基础上，由虚位移原理可得 由于应变、应力以矩阵表示，因此弹性矩阵应按下式并考 虑应变矩阵定义来建立 式中 为单元节点力矩阵， 为单元等效节点荷载矩阵 第五节 几何非线性问题的有限元法 或 Pd和PE分别为直接和等效节点荷载矩阵，P为综合等效节点荷载矩阵。上式也可写为 将克希霍夫应力表达式代入，可得 按集成规则组装后可得 第五节 几何非线性问题的有限元法 根据本构方程，则有 式中K(U)是非对称的 对非线性弹性问题 、 和 都是位移的函数。 根据非线性切线刚度矩阵的定义，可得 又因 B、G为已知矩阵 第五节 几何非线性问题的有限元法 式中 所以 因 所以 第五节 几何非线性问题的有限元法 由此可得 式中 基于上述说明，可得 第五节 几何非线性问题的有限元法 如果引入如下记号 则单元切线刚度矩阵为（三个部分）： 初应力或几何刚度矩阵 线弹性刚度矩阵 大位移刚度矩阵 “结构”切线刚度矩阵为 建立了切线 刚度矩阵， 用牛顿法等 即可求解。 第五节 几何非线性问题的有限元法 1.上述讨论没涉及具体单元，因此具有普遍性。具体单元分 析时，因形函数一般是自然坐标的函数，故需作坐标变换后代入 相关公式，从而建立具体单元的切线刚度矩阵。 2.本节讨论了弹塑性大变形问题有限元法，而为了保证收 敛，应该用增量迭代法进行求解。 3.对于本节讨论的几何大变形问题，如果用增量形式来解表 达，刚度矩阵只是在弹性基础上增加了两项，牛顿－拉夫森算法 完全可以求解。 注： 第五节 几何非线性问题的有限元法 （3）结构稳定性问题 对分支点稳定问题，关键是建立几何刚度矩阵。此时，以失稳前的平衡位置作初始位形，以失稳形态作现时位形。 分支点失稳问题和有限几何大变形问题的区别有： （2）大位移刚度矩阵可忽略。 （3）应变仅包含失稳位移的非线性项。 （4）分支点处相应失稳位移的综合荷载为零。 （1）失稳前变形是微小的。 第五节 几何非线性问题的有限元法 下面以二维薄板单元为例，讨论稳定性问题。 有限元里的二维问题单元位移为 但失稳时有平面外位移w，应变需考虑w及其非线性项 非线性项 第五节 几何非线性问题的有限元法 设平面外位移为 面内弹性应力为 对平面外位移，其G矩阵为 对应的M 矩阵 第五节 几何非线性问题的有限元法 对各种具体单元，将Nw的具体形函数代入G矩阵的表达式，即可积分得到具体单元的几何刚度矩阵。 基于上述结果，单元几何刚度矩阵为 对于分支点稳定问题，结构力学已经指出，在比例加载下，最终归结为一个特征值问题 解得特征值后，即可得到临界荷载，一般只关心最小临界荷载。 第五节 几何非线性问题的有限元法 第五节 几何非线性问题的有限元法 总 结： 在这一节里，我们分别讨论了三种几何非线性问题，小变形几何非线性问题、有限变形几何非线性问题以及结构稳定性问题，后两个分析将是我们以后学习的重点，应该搞清楚它们之间的关系和应用领域的区别。 对于非线性有限大变形问题，应该搞清楚几个应力、应变的概念和区别。 参考文献 徐兴、郭乙木等编著，非线性有限元及程序设计，浙江大学出版社，1993 2. 郭乙木等编著，线性与非线性有限元及其应用，机械工业出版社，2003 3. 孙柄楠等编著，工程弹塑性力学，浙江大学出版社， 1999 4. 杨伯源等编著，工程弹塑性力学，机械工业出版社，2003 5. 丁浩江等编，弹性和塑性力学中的有限元法，机械工业出版社，2004 思考题 1.固体力学中有几种“非线性”？试叙述它们的异同。 2.目前关于塑性力学的本构关系存在哪两种理论？ 试分析其异同。 3.试叙述非线性弹性和弹塑性之间的区别。 4.求解非线性有限元方程的方法有几种？试分析它们 的优缺点？ 5.几何非线性问题可以分为几种？试叙述它们之间的 区别和各自的应用领域。 6.在有限元分析中，应力和应变有哪几种？试叙述它 们之间的区别。 思考题 7.弹塑性问题进行有限元分析过程中，必然会出现结构 某一部分进入塑性，而另外一部分仍属于弹性阶段，