第九部分变分法课程.pptVIP

* * 即 可得 由 即 解得 * * 里滋方法的步骤: 上述方法称为里滋法，其要点是要找到满足全部边界条件的位移函数，而这种函数一般仍然难以找到，尤其在边界不规整的情况下。所以里滋方法的应用在这一点上受到极大的限制。五十年代以来，在这基础上发展起来的有限元方法，采用了单元离散，分片插值的方法，这就避免了这一困难，虽然带来了计算工作量大的问题，但由于电子计算机的产生和发展，计算工作量大的问题得到解决，有限元方法得到迅速的发展。 * * 我们学习变分方法的意义，主要因为它是学习有限元等数值方法的基础。 * * 第三节 伽辽金方法 * * 里滋方法要求位移函数满足位移边界条件，如果进一步要求根据位移函数求得的应力还满足应力边界条件，公式还可以简化，这种方法称为伽辽金方法。 由前面我们得到的最小势能原理 * * 如果我们所取的位移不仅满足位移边界条件，而且根据它们求得的应力还满足应力边界条件(不要求满足平衡方程），则上式简化为 总势能为 * * 取位移为 则上式为 * * 由各待定常数变分的独立性，各系数为零 以位移表示 * * 在平面应变问题中，方程为 在平面应力问题中，方程为 这样就得到位移函数待定常数的线性方程组，解之，再代回，就得到位移的近似解，根据应力的位移表示式，就可求得应力。通常位移在所取项不多时，就可得到较精确的解，但应力解需要较多的项，这也是通常位移法所遇到的问题。 * * 伽辽金方法的计算工作量较小，但对位移函数的要求较高，要求应力应满足应力边界条件。在特殊情况，如仅有位移边界，而无应力边界，这也表示着应力边界条件得到满足，这时用伽辽金方法十分方便。 作为例子，可见徐芝纶“弹性力学”p358-p361. * * 第四节 应力变分方法 * * 设有任一弹性体，在外力的作用下处于平衡。命σij为实际存在的应力分量，它们满足平衡微分方程和应力边界条件，也满足相容方程。现在，假想体力不变,而应力分量发生了微小的变化δσij，即所谓虚应力或应力的变分，使应力分量成为σij+δσij，设它们只满足平衡微分方程和应力边界条件。 * * 既然两组应力分量都满足同样体力作用下的平衡微分方程，应力分量的变分必然满足无体力时的平衡方程，即 (a) * * 同时，在位移给定的边界上（面力不可能给定），应力分量的变分必然伴随着面力分量的变分δpx、δpy、δpz 根据应力边界条件的要求，应力分量的变分在边界上必须满足 由于应力分量的变分，形变势能必有相应的变分。把形变势能看做应力分量的函数，则形变势能的变分应为 (b) * * 将下式代入 再将几何方程代入，得 根据分步积分和奥－高公式，对上式右边的各项进行处理，例如 * * 最后可得 再将式（a）和（b）代入，即得 * * 这就是所谓应力变分方程(卡斯提安诺变分方程)。这方程的右边代表面力的变分在实际位移上所做的功。由此可见，由于应力发生的变分，形变势能的变分等于面力的变分在实际位移上所做的功。上式可改写为 括号中的表达式称为弹性体的余能，因此，在满足平衡方程和应力边界条件的各组应力中，实际存在的应力应使弹性体的余能成为极值，称为最小余能原理。通过推导，可以从余能原理可以得到相容条件。 * * 如果在某一部分边界上，面力是给定的，则该部分边界上的面力不能有变分，应力变分方程中右边该部分的积分为零，如果位移边界上的位移为零，方程中右边该部分的积分也为零。实际积分只对位移边界上的位移不为零的那部分进行，面力的变分按下式进行： 如果在全部边界上，面力是给定的，或位移边界上的位移为零，方程中右边该部分的积分为零，应力变分方程简化为 * * 注:在线弹性时,弹性体的应变余能等于应变能,在非线弹性时,上面的弹性应变能应为应变余能. 近似解法： 设定应力分量的表达式，使其满足平衡方程和应力边界条件，但其中包含若干待定系数，然后根据应力变分方程决定这些系数。应力分量一般可设为 其中Am为互不依赖的m个系数。σ0给出的应力满足实际的应力边界条件，σm给出的应力满足无面力时的应力边界条件。 * * 由于应力分量的数量较多，确定起来有困难，通常用应力函数方法。 在平面应力问题中，如果体力分量为常数，则存在应力函数。将应力函数设定为 其中Am为互不依赖的m个系数。φ0给出的应力满足实际的应力边界条件，φm给出的应力满足无面力时的应力边界条件。 * * 在平面应力状态，用应力分量表示的变形势能为 如果考虑单连体，而且是应力边界问题，应力分量应与ν无关，可设其为零，则有 用应力函数表示为 * * 在应力边界问题中，变分方程简