非线性有限元 第三章 连续介质力学.ppt

5 守恒方程 4 能量守恒 在域内由体积力,和在表面上由面力做的功率为 在物体中总能量的变化率为 由热源s和热流q提供的功率为 其中热流一项的符号是负的，因为正的热流是向物体外面流出的 能量守恒 5 守恒方程 4 能量守恒 即物体内总能量的变化率（包括内能和动能）等于外力的功率和由热流及热能源提供的功率。这是已知的热力学第一定律。 内能的支配依赖于材料。在弹性材料中，它以内部弹性能的形式存储起来，并在卸载后完全恢复；在弹塑性材料中，部分内能转化为热，部分由于材料内部结构的变化而耗散了。 应用Reynold定理将求导数移入积分内，然后将面积分转换为域积分 5 守恒方程 4 能量守恒 将Cauchy定律和Gauss定理应用于面力边界积分，得到 代入能量守恒公式，对热流积分应用Gauss定理，并整理各项得到 动量方程,为0 5 守恒方程 4 能量守恒 由域的任意性,得到能量守恒的偏微分方程 当没有热流和热源时，即为一个纯机械过程，能量方程成为 这不再是一个偏微分方程，它以应力和应变率度量的形式，定义了给予物体单位体积的能量变化率；称为内能变化率或内部功率。由变形率和Cauchy应力的缩并给出内部功率。变形率和Cauchy应力在功率上是耦合的。 功率上的耦合有助于弱形式的建立：在功率上耦合的应力和应变率的度量可以用于构造虚功原理或虚功率原理，即动量方程的弱形式。在功率上耦合的变量也可以说在功或者能量上是耦合的，但是常常使用功率耦合的说法，因为它更加准确。 5 守恒方程 6 Lagrangian守恒方程 以应力和应变的Lagrangian度量形式，在参考构形中直接建立守恒方程是有益的。在连续介质力学的文献中，这些公式称为Lagrangian描述，而在有限元的文献中，这些公式称为完全的Lagrangian格式。 对于完全的Lagrangian格式，总是使用Lagrangian网格。在Lagrangian框架中的守恒方程与刚刚建立的守恒方程基本上是一致的；它们只是以不同的变量表示。实际上将看到，可以通过框3.2中的转换关系和链规则得到它们。 6 Lagrangian守恒方程 在完全的Lagrangian格式中，独立变量是材料坐标X和时间t。主要的相关变量是初始密度ρ0(X，t)，位移u(X，t)以及应力和应变的Lagrangian度量。 使用名义应力P(X，t)作为应力的度量。这导致动量方程与Eulerian描述的动量方程(3.5.33)惊人的相似，所以非常容易记忆。变形将通过变形梯度F(X，t)描述。 对于构造本构方程，使用成对的P和F不是特别有用的，因为F在刚体运动中不为零，而P是不对称的。因此，本构方程通常表示为PK2应力S和Green应变E的形式。然而，通过框3.2中的转换关系，S和E之间的关系可以很容易的转换为P和E之间的关系。 6 Lagrangian守恒方程 7 极分解和框架不变性 目的是探讨刚体转动的作用： 表述极分解定理，该定理能够从任何运动中得到刚体转动。 考虑刚体转动对于本构方程的影响。证明对于Cauchy应力，需要对时间导数进行修改建立率－本构方程。这就是框架不变性或者应力的客观率。 表述三种框架不变率：Jaumann率。Truesdell率和Green－Naghdi率。 展示了由于次弹性本构方程和这些不同变化率的错误应用，在结果中的惊人误差。 极分解定理 7 极分解和框架不变性 在大变形问题中，阐明转动作用的基本原理就是极分解定理。这个定理表述为，任何变形梯度张量F可以乘法分解为一个正交矩阵R和一个对称张量U的乘积，称U为右伸长张量 (先伸长再转动)。 物体的任何运动包括一个变形，由对称映射U表示，和一个刚体转动R；所有的正交变换都是转动。在这个方程中没有出现刚体平动，因为dx和dX分别是在当前和参考构形中的微分线段，而且微分线段的映射不受平动的影响。 如果将方程积分得到x＝Φ(X,t)的形式，那么刚体平动将作为一个积分常数出现。在刚体平动中，F＝I，和dx＝dX。 其中 有 7 极分解和框架不变性 极分解定理证明 得到 右边总是一个正矩阵，所以矩阵U的所有特征值总是正值，故U的逆矩阵存在 矩阵U与工程应变联系得非常紧密。它的主值是在矩阵U的主方向上线段的伸长。其吸引人之处在于建立本构方程。张量U－I 称为Biot应变张量。 一个运动也可以分解为一个左伸长张量和一个转动的形式 称V为左伸长张量（先转动再伸长）。 7 极分解和框架不变性 通过极分解定理分别求在t＝1.0 和 t＝0.5 时的刚体转动和伸长张量 考虑三角形单元的运动，其中节点坐标xI (