第5章弹性静力学小位移变形理论的变分原理(16K).doc

第5章 弹性静力学小位移变形理论的变分原理 对连续体来说，其数学上的处理方法是利用给定的边界条件下的微分方程（或偏微分方程），并在一定的边界条件下求得其解，这种解析方法，实际做起来往往遇到很大的困难，使许多工程实际问题的计算模型很难建立，满足不了实际需要。自从五十年代直刚法问世以来，利用离散化的方法，将一个连续体划分为有限数量及具有一定几何形状的单元体，即有限单元，再按照一定的过程进行计算，这就使得过去许多工程计算感到困难的问题得到解决，这种方法不受结构特殊几何形状的限制，因此，它的适应范围是相当广泛的。有限元素法的提出和应用，是工程分析方法上的一次重大的变革，随着理论探讨上的深入及计算机性能的不断提高，使得解的精确性不断地得到改进，以至使得有限元素法成为当前计算领域方面的一个强有力的工具，无论对结构问题（如静力学、动力学）、非结构问题（如流体力学、光学、电磁学）以及许多边缘学科等都得到广泛的应用。 有限元素法的解题过程和步骤在一般的有关有限元法教课书和著作中均有详细讨论，本章不再赘述。 变分原理是有限元素法的基础，要很好地理解有限元素法，则应该对能量变分原理有一个较系统地了解。本章的目的是尽可能地对这些能量变分原理作系统性的介绍，从一般常用的最小位能原理和最小余能原理，引深到引用拉格朗日乘子法（Lagrange Multiple Method）的完全及不完全广义变分原理和为分区集合体的分区（Sub-region）广义变分原理，这将涉及到以混合（Mixed）模型和杂交（Hybrid）模型为基础的变分原理。在此基础上，针对不同变分原理，进一步说明了有限元素法中的元素的刚度特性和推导元素刚度矩阵的一般过程及表达显式，以及变分原理在结构分析中的若干应用实例，使读者能比较清晰地了解各类变分原理与建立有限元模型之间的关系。 §5.1 小位移弹性理论的最小位能原理与最小余能原理 设在卡氏直角坐标系中，坐标参数为，体积为的弹性体中任意一点的位移参数为、应力分量为以及应变分量为。由线弹性力学理论，我们可以得到如下的用于描述一个弹性静力学小位移变形问题的基本方程式。 （1）力的平衡方程 （在内） （5-1） 式中表示体力，表示应力分量对坐标分量的偏导数（以下相同）。 （2）应变位移关系式（几何关系） （在内） （5-2） （3）应力应变关系式（物理关系） （5-3） （5-3’） 式中为弹性模量系数，为劲度系数，和都具有对称性。 （4）在弹性体的边界上，表面S可划分为两部分：外力已知的边界及位移为已知的边界，前者称为力的边界，后者称为位移边界，即 （5-4） 在力的边界上， （5-5） 式中为已知边界力，为的边界外法线向量与坐标轴夹角的方向余弦。 在位移边界上， （5-6） 式中为已知边界位移。（5-5）式和（5-6）式统称为“边界条件”。 上述的诸方程共有15个，即3个平衡方程，6个应变位移关系方程，6个物理关系方程。而未知变量也共计15个：6个应力分量，6个应变分量和3个位移分量。因此该问题是可以求解的。 小位移变形弹性体的应变能泛函（或应变能密度）和余应变能泛函（余应变能密度）可表示为 （5-9） （5-10） 不难看出，和有以下关系， （5-11） 并且容易证明 （5-12） （5-13） （一）虚功原理与总位能原理 这里用和分别表示应变变分和位移变分，在虚功原理中可视为虚应变和虚位移。则由虚功原理可写出虚功方程为 （5-14） (5-14)式成立是有条件的，要求和在弹性体内部满足应变位移关系和在位移边界上满足给定位移边界条件，即 （在V内） （5-15a） （在上） （5-15b） 虚功原理表明，如果弹性体在给定的体力和边界力作用下处于平衡状态，则对于为位移边界条件所容许的任意虚位移，（5-14）式成立。反过来，如果（5-14）式对于为位移边界条件所容许的任意虚位移成立，则弹性体处于平衡状态。值得提出的是，不管材料的应力应变关系是线性还是非线性，虚功原理都成立。 如果用下