弹性力学-能量原理与变分法课件

1、 弹性理论问题需要解一系列偏微分方程组，并满足边界条件，这在数学上往往遇到困难。因此需要寻求近似的解法。变分法的近似解法是常用的一种方法。在数学上，变分问题是求泛函的极限问题。在弹性力学里，泛函就是弹性问题中的能量（功），变分法是求能量（功）的极值，在求极值时得到弹性问题的解，变分问题的直接法使我们比较方便地得到近似解。 本章首先给出计算形变势能的表达式。利用功与能的关系，主要介绍了位移变分法和应力变分法。1能量原理与变分法第十章 能量原理与变分法10-1 弹性体的变形比能与形变势能10-3 位移变分法10-4 应力变分方程与应力变分方法110-2 位移变分方程与极小势能原理能量原理与变分法1

2、0-1 弹性体的变形比能与形变势能一 变形比能 在复杂应力状态下，设弹性体受有全部六个应力分量 。根据能量守恒定理，形变势能的多少与弹性体受力的次序无关，而完全确定于应力及形变的最终大小。从而有弹性体的形变势能密度或比能：比能用应力分量表示2能量原理与变分法比能用应变分量表示其中因此，我们有比能对应力分量的偏导3能量原理与变分法比能对应变分量的偏导二 形变势能 由于应力分量和形变分量，进而比能 都是位置坐标的函数，所以整个弹性体的形变势能 为： 将比能的三种表达形式代入，得形变势能的三种积分形式4能量原理与变分法将几何方程代入，形变势能还可用位移分量来表示5能量原理与变分法一 变分及其性质 高

3、等数学我们学过微分的概念，微分是变量的增量。那么什么是变分呢？变分是函数的增量，通常用表示。变分具有以下的性质： 10-2 位移变分方程与极小势能原理6能量原理与变分法二 位移变分方程 设弹性体在一定外力作用下，处于平衡状态，发生的真实位移为 u,v,w，它们满足位移分量表示的平衡方程，并满足位移边界条件和用位移表示的应力边界条件。现在假设位移分量发生了位移边界条件所容许的微小改变（虚位移）u 、v、w，这时外力在虚位移上作虚功，虚功应和变形能泛函的增加相等，即这个方程就是所谓位移变分方程。其中X,Y,Z为体力分量， 为面力分量。7能量原理与变分法三 极小势能原理 由于虚位移是微小的，因此在虚

4、位移的过程中，外力的大小和方向可以当做保持不便，只是作用点有了改变。利用变分的性质，位移变分方程可改写为：设外力势能为则 该式的意义是：在给定的外力作用下，在满足位移边界条件的各组位移中，实际存在的一组位移应使总势能为最小。如果考虑二阶变分，进一步的分析证明，对于稳定平衡状态，这个极值是极小值。因此，该式又称为极小势能原理。8能量原理与变分法 显然，实际存在的位移，除了满足位移边界条件以外，还应当满足位移表示的平衡方程和应力边界条件；现在又看到，实际存在的位移，除了满足位移边界条件外，还满足位移变分方程。而且，通过运算，还可以从位移变分方程导出用位移表示的平衡微分方程和应力边界条件。于是可见：

5、位移变分方程可以代替平衡微分方程和应力边界条件。9能量原理与变分法 10-3 位移变分法其中u0,v0,w0 为设定的函数，它们的边界值等于边界上的已知位移；um 、vm、wm 为边界值等于零的设定函数，Am、Bm、Cm为待定的系数，位移的变分由它们的变分来实现。 先设定满足位移边界条件的位移分量的表达式，其中包含若干个待定的系数，再根据极小势能原理，决定这些系数。取位移分量的表达式如下：一 瑞次法10能量原理与变分法应变能的变分为外力势能的变分为位移分量的变分是11能量原理与变分法中，得到上面是个数为3m的线性代数方程组，求解后，代回位移分量的表达式，得到位移分量的近似解。这种方法称为瑞次法

6、。代入12能量原理与变分法二 伽辽金法 将变分看做形变分量的函数，则由于所以应用奥高公式，对上式中的第一项，我们有13能量原理与变分法对于其余各项也进行同样的处理，则将上式代入位移变分方程，并归项得14能量原理与变分法如果应力边界条件得到满足，则上式简化为这就是位移分量满足位移边界条件及应力边界条件时，位移变分所应满足的方程，称为伽辽金变分方程。15能量原理与变分法 若取位移分量的表达式如下：使得位移边界条件和应力边界条件都得到满足，则将位移变分代入伽辽金方程，就得到16能量原理与变分法由于的任意性，它们的系数应当分别为零于是得将上列三方程中的应力分量通过物理方程用形变分量表示，再通过几何方程

7、用位移分量表示，简化后即得：17能量原理与变分法 这样就得到位移函数待定常数的线性方程组，求解后，代回位移分量的表达式，得到位移分量的近似解。这种方法称为伽辽金法。 要注意的是：用位移变分法求位移分量，只须取几项就可达到较高的精度，然而由此求出的应力却很不精确。为了求得的应力充分精确，必须取更多的项。18能量原理与变分法三 应用与举例 将位移变分法应用于平面问题，瑞次法和伽辽金法都将得到简化。由于两种平面问题都不必考虑 z 方向的位移 w ，且u 和 v 都不随坐标 z 而变，所以位移分量的表达式可设为 在采用瑞次法时，为了决定系数Am及 Bm ，在z方向取一个单位长度，只须应用如下二式来求解

8、线性方程组19能量原理与变分法其中形变势能用位移分量表示形式简化如下，平面应变问题平面应力问题 在采用伽辽金法时，对于平面应变问题，要应用如下二式来求解线性方程组20能量原理与变分法对于平面应力问题，要求解的方程组如以下二式形式 伽辽金方法的计算工作量较小，但对位移函数的要求较高，除了要求满足位移边界条件外，还要求根据位移函数求得的应力应满足应力边界条件。在特殊情况，如仅有位移边界，而无应力边界，这也表示着应力边界条件得到满足，这时用伽辽金方法十分方便。 瑞次法的要点是要找到满足全部边界条件的位移函数，而这种函数一般仍然难以找到，尤其在边界不规整的情况下。所以瑞次方法的应用在这一点上受到极大的

9、限制。21能量原理与变分法例1：如图所示的薄板，不计体力，求薄板的位移解： 设位移它们是满足位移边界（左边和下边）的边界条件的。 在平面应力状态下可得xyq1q2bao22能量原理与变分法即可得由即解得23能量原理与变分法例2：如图所示，宽为2a而高度为b的矩形薄板，左右两边及下边均被固定，而上边的位移给定为不计体力，试求薄板的位移和应力。 解： 取坐标轴如图所示。设位移分量为aabboxy24能量原理与变分法可以满足位移边界条件，即 在该问题中，并没有应力边界条件，因此可以认为所设位移既然满足了位移边界条件，也就满足了全部边界条件，这就可以应用伽辽金法求解，使数学运算比较简单一些。 注意体力

10、X=Y=0而m = 1，伽辽金方程成为25能量原理与变分法将位移分量的各二阶导数以及代入伽辽金方程，进行积分并求解得26能量原理与变分法 为简单起见，取b=a而=0.2，将A1和B1代入所设位移函数得应用几何方程及物理方程，可有上式求得应力分量27能量原理与变分法 10-4 应力变分方程应力变分方法 设有任一弹性体，在外力的作用下处于平衡。命ij为实际存在的应力分量，它们满足平衡微分方程和应力边界条件，也满足相容方程。现在，假想体力不变,而应力分量发生了微小的变化ij，即所谓虚应力或应力的变分，使应力分量成为ij +ij ，设它们只满足平衡微分方程和应力边界条件。一 应力变分方程 既然两组应力

11、分量都满足同样体力作用下的平衡微分方程，应力分量的变分必然满足无体力时的平衡方程，即28能量原理与变分法(a)(b)同时，在位移给定的边界上（面力不可能给定），应力分量的变分必然伴随着面力分量的变分 。根据应力边界条件的要求，应力分量的变分在边界上必须满足29能量原理与变分法 由于应力分量的变分，形变势能必有相应的变分。把形变势能看做应力分量的函数，则形变势能的变分应为将下式代入30能量原理与变分法再将几何方程代入，得 根据分步积分和奥高公式，对上式右边的各项进行处理，例如最后可得31能量原理与变分法再将式（a）和（b）代入，即得这就是所谓应力变分方程(卡斯提安诺变分方程)。方程的右边代表面力

12、的变分在实际位移上所做的功。由此可见，由于应力发生的变分，形变势能的变分等于面力的变分在实际位移上所做的功。32能量原理与变分法 如果在某一部分边界上，面力是给定的，则该部分边界上的面力不能有变分，于是 ，而应力变分方程右边的相应积分项成为零；如果在某一部分边界上，给定的位移等于零，则应力变分方程右边的相应积分项也成为零。因此，应力变分方程右边的积分，只须在这样的边界上进行：面力没有给定，而给定的位移又不等于零。二 极小余能原理 将应力变分方程改写为由于在需要积分的边界上，位移是给定的，在变分过程中保持不变，所以上式可以改写为33能量原理与变分法 中括号内的表达式称为弹性体的余能。因此，在满足

13、平衡方程和应力边界条件的各组应力中间，实际存在的一组应力应使弹性体的余能成为极值。如果考虑二阶变分，可以证明这个极值是极小值，所以上述结论称为极小余能原理。 以前看到，实际存在的应力，除了满足平衡微分方程以外，还应当满足相容方程。现在又看到，实际存在的应力，除了满足平衡微分方程和应力边界条件以外，还满足应力变分方程。而且，通过运算，还可以从应力变分方程导出相容条件。于是可见，应力变分方程可以代替相容条件。34能量原理与变分法三 应力变分法 设定应力分量的表达式，使其满足平衡方程和应力边界条件，但其中包含若干待定系数，然后根据应力变分方程决定这些系数。应力分量一般可设为 其中Am为互不依赖的m个

14、系数。ij0是满足平衡微分方程和应力边界条件的设定函数，ijm 是满足“无体力和面力作用时的平衡微分方程和应力边界条件”的设定函数。这样，不论系数Am如何取值，ij总能满足平衡微分方程和应力边界条件。 注意：应力的变分只是由系数Am的变分来实现。(c)35能量原理与变分法 如果在弹性体的每一部分边界上，不是面力被给定，便是位移等于零，则应力变分方程简化为显然，形变势能U是Am的二次函数，因而(d)式将是Am的一次方程。这样的方程共有 m 个，恰好可以用来求解系数Am，回代(c)式，求得应力分量。即(d) 如果在某一部分边界上，位移是给定的，但并不等于零，则在这一部分边界上须直接应用变分方程，即

15、36能量原理与变分法在这里，u、v、w是已知的，积分只包括该部分边界。将面力的变分与应力的变分两者之间的关系，即代入方程的右边积分后，将得出如下的结果：其中Bm是常数。另一方面，方程的左边37能量原理与变分法因而得上式仍然是Am的一次方程，总共有m个，且各个Am是互不相关的，因而可以求出所有的Am ，回代(c)式，求得应力分量。 在应用应力变分法时，要使设定的应力分量既满足应力边界条件，又满足平衡微分方程，这往往是很困难的。但是，在某些类型的问题中存在着应力函数，而且用应力函数表示的应力分量又能满足平衡微分方程。这时，我们就只须设定应力函数的表达式，使它给出的应力分量能满足应力边界条件，困难就

16、大大减少了。38能量原理与变分法 其中Am为互不依赖的m个系数。0给出的应力满足实际的应力边界条件，m 给出的应力满足无面力时的应力边界条件。 由于应力分量的数量较多，确定起来有困难，通常用应力函数方法。 在平面应力问题中，如果体力分量为常数，则存在应力函数。将应力函数设定为四 应力函数方法 在平面应力状态，用应力分量表示的形变势能为39能量原理与变分法 如果考虑单连体，且是应力边界问题，应力分量应与 无关，可设其为零。则两类平面问题皆简化为用应力函数表示为 对于平面应变问题40能量原理与变分法 在应力边界问题中，因为面力不能有变分，所以变分方程简化为U=0,因此系数应满足 上式为线性方程组，

17、求解Am后，得到应力函数的近似解，最后得到各应力分量。将应力函数的表达式代入，得 41能量原理与变分法 由于是近似解，应力分量不能精确满足相容条件，由应力分量求得的应变分量也不能精确满足变形协调条件，不能根据几何方程求得位移分量。 由上式即可解得系数Am。从而确定应力函数，再由应力函数求得各应力分量。 应力函数法的要点是要找到满足全部边界条件的应力函数，而这种函数一般仍然难以找到，尤其在边界不规整的情况下。所以应力方法的应用在这一点上受到极大的限制。42能量原理与变分法例题3：设有矩形薄板，体力不计，在两对边上受有按抛物线分布的拉力，其最大集度为 q ，图示。试求板的应力。 解：在图示坐标下，

18、边界条件是取xyaabbq43能量原理与变分法则由给出的应力分量满足边界条件。而由所对应的应力能满足无面力时的边界条件。 将以及体力X=Y=0代入下式中44能量原理与变分法进行积分，并简化以后得为简单起见，令a=b，则代入中，并令a=b，再求应力分量得45能量原理与变分法能量原理与变分法习题10.1 铅直平面内的正方形薄板，边长为2a，四边固定，只受重力作用，设 ，试取位移分量的表达式为用瑞次法求解位移分量(取A1项及B1项)及应力分量。aaaaxy解： 当只取A1项及B1项时：能量原理与变分法形变势能现计算 和 能量原理与变分法在用瑞次法时，要求由解之得能量原理与变分法所以能量原理与变分法习题10.2 简支梁受均布荷