第2章 能量原理与变分法2.ppt

里滋方法的步骤 上述方法称为里滋法，其要点是要找到满足全部边界条件的位移函数，而这种函数一般仍然难以找到，尤其在边界不规整的情况下。所以里滋方法的应用在这一点上受到极大的限制。五十年代以来，在这基础上发展起来的有限元方法，采用了单元离散，分片插值的方法，这就避免了这一困难，虽然带来了计算工作量大的问题，但由于电子计算机的产生和发展，计算工作量大的问题得到解决，有限元方法得到迅速的发展。 我们学习变分方法的意义，主要因为它是学习有限元等数值方法的基础。 伽辽金方法 里滋方法要求位移函数满足位移边界条件，如果进一步要求根据位移函数求得的应力还满足应力边界条件，公式还可以简化，这种方法称为伽辽金方法。 由前面我们得到的最小势能原理 如果我们所取的位移不仅满足位移边界条件，而且根据它们求得的应力还满足应力边界条件(不要求满足平衡方程），则上式简化为 总势能为 取位移为 则上式为 由各待定常数变分的独立性，各系数为零 以位移表示 在平面应变问题中，方程为 在平面应力问题中，方程为 这样就得到位移函数待定常数的线性方程组，解之，再代回，就得到位移的近似解，根据应力的位移表示式，就可求得应力。通常位移在所取项不多时，就可得到较精确的解，但应力解需要较多的项，这也是通常位移法所遇到的问题。 伽辽金方法的计算工作量较小，但对位移函数的要求较高，要求应力应满足应力边界条件。在特殊情况，如仅有位移边界，而无应力边界，这也表示着应力边界条件得到满足，这时用伽辽金方法十分方便。 作为例子，可见徐芝纶“弹性力学”p358-p361，解法参见MCAD. 解：首先用瑞利—里茨法 位移试函数 例1： 两端简支的等截面梁，受均匀分布载荷q作用如图所示。试求解梁的挠度w（x）。 满足梁的位移边界条件 在x=0，l处，w=0 总势能 根据 则 所以 回代 挠曲线表达式是无穷级数——精确解 这个级数收敛很快，只要取少数几项就可以得到足够的精度。 如果取一项 这一结果与精确值十分接近 解：应用Гапёркин法 位移试函数 同时满足面力边 界条件 根据Гапёркин法 分析可得 解： 位移试函数 例2：矩形薄板，四边固定，受有平行于板面的体力作用。设坐标轴如图所示，试用Rayleigh—Ritz法求解。 m和n为正整数 在边界x=0，a，和y=0，b上，u=v=0，所以试函数满足位移边界条件。 平面应力问题 因此 将位移试函数代入求导数后再积分 因此 如果体力已知，积分可求待定系数Amn和Bmn 应力变分方法 设有任一弹性体，在外力的作用下处于平衡。命σij为实际存在的应力分量，它们满足平衡微分方程和应力边界条件，也满足相容方程。现在，假想体力不变,而应力分量发生了微小的变化δσij，即所谓虚应力或应力的变分，使应力分量成为σij+δσij，设它们只满足平衡微分方程和应力边界条件。 既然两组应力分量都满足同样体力作用下的平衡微分方程，应力分量的变分必然满足无体力时的平衡方程，即 (a) 同时，在位移给定的边界上（面力不可能给定），应力分量的变分必然伴随着面力分量的变分δpx、δpy、δpz 根据应力边界条件的要求，应力分量的变分在边界上必须满足 由于应力分量的变分，形变势能必有相应的变分。把形变势能看做应力分量的函数，则形变势能的变分应为 (b) 将下式代入 再将几何方程代入，得 根据分步积分和奥－高公式，对上式右边的各项进行处理，例如 * * 第二章 能量原理和变分法 真实的位移除了满足位移边界条件外，根据它们求得的应力还应满足应力边界条件和平衡微分方程。求解微分方程的边值问题，只有在简单的情况下，才能得到解析解。多数情况下，只能采用数值计算的方法。 基于能量原理的变分法为数值计算提供了理论基础。其中基于最小势能原理的里滋方法等可用于数值计算。 应变能 虚位移原理 最小势能原理 位移变分法 里滋方法 伽辽金方法 应力变分法 设弹性体在一定外力作用下，处于平衡状态，发生的真实位移为u,v,w，它们满足位移分量表示的平衡方程，并满足位移边界条件和用位移表示的应力边界条件。弹性体受力后，发生变形，外力作功，外力功转化为应变能，储存在弹性体内，单元体内的应变能为 第一节 应变能 或 整个弹性体内的应变能为 应变余能的概念 以一维应力状态为例，U0实际是应力应变曲线下的面积（不限于线弹性） 定义 为单位体积的应变余能，在一维情况下为 εx σx dεx O dσx 应变余能没有明显的物理意义，在一维情况下，表示应力应变曲线在应力一侧下的面积。 1 应变余能与应变能互补 2 应变余能的积分式中，积分变量为应力分量 3 在线