弹塑性力学-第8章 能量原理及其应用.doc

1、第八章 能量原理及其应用193第八章 能量原理及其应用弹塑性力学问题实质上是边值问题，即求解满足一定边界条件的偏微分方程组。然而只有对一些特殊的结构在特定加载条件下才能找到精确解，而对于一般的力学问题，如空间问题，泛定方程为含有 15 个未知量的 6 个偏微分方程，在给定边界条件时求解是极其困难的，而且往往足小对能的。因此，为了解决具体的工程结构力学问题，目前都广泛应用数值方法，如有限元法、无限元法、边界元法、无网格化法及样条元法等等。这些解法的依据都是能量原理。本章将讨论利用能量原理和极值原理求解弹塑性力学问题的近似解法。本章共讨论五个能量原理。首先是虚位移原理，由虚位移原理推导出最小势能原

2、理，其次介绍虚应力原理，和由虚应力原理推导出最小余能原理。另外，还简单介绍最大耗散能原理。本章还讲述了根据上述的能量原理建立的有关弹性力学问题的数值解法。8.1 基本概念1.1 物体变形的热力学过程由第四章知，物体在外界因素影响下的变形过程，严格来说都是一个热力学过程。因此研究物体的状态，不仅要知道物体的变形状态，而且还要知道物体中每一点的温度。如果物体在变形过程中，各点的温度与其周围介质的温度保持平衡，则称这一过程为等温过程；若在变形过程中，物体的温度没有改变，即既没有热量损失也没有热量增加，则称这一过程为绝热过程。物体的瞬态高频振动，高速变形过程都可视为绝热过程。令物体在变形过程中的动能为

3、 E，应变能为 U，则在微小的 时间间隔内，t物体从一种状态过渡到另一种状态时，根据热力学第一定律，总能量的变化为(a)QW其中， 为作用于物体上的体力和面力所完成的功； 是物体由其周围介质所W吸收(或向外发散)的热量，并以等量的功度量。假定弹性变形过程是绝热的，则对于静力平衡问题有(b)00上E将式(b)代入式(a)，则有(8.1-1)WU第八章 能量原理及其应用1941.2 应变能由第四章的式(4.1-5b)知，在线弹性情况下，单位体积的应变能为(8.1-2)ikjijijdU210对于一维应力状态，在 平面内，则x实际上就是应力应变曲线与 轴和 所0 x围成的面积(图 8.1)，即(8.

4、1-3)0Xxd其中 是物体变形过程某一指定时刻的应变，应 图 8.1 应变能与应变余能x变能 表示物体在变形过程中所储存的能量。 0U1.3 应变余能在图 8.1 中, 如果令 表示应力应变曲线与 轴和 所围成的面积，0Uxx即(8.1-4)0xd式中 是物体变形过程某一指定时刻的应力。称 为单位体积的应变余能，简x 0U称余能，有时又称其为应力能。由于 和 是物体变形过程中同一时刻的应力和应变，因此有x0x由此可见， 与 互补或互余对方为 矩形的面积。显然，在线弹性情况下0U 有 ，即余能与应变能在数值上相等。0尽管应变余能不像应变能那样具有明确的物理意义，但引入应变余能这一概念后，使讨论

5、问题的范围扩大了。8.2 虚位移原理与最小势能原理2.1 虚位移原理设有变形体在外力作用下处于平衡状态。此处，外为包括体力分量 ， ， 及XYZ一部分表面的面力分量 ， ， 。假如有一组位移分量 ， 既能满足用位XYZu上vw移表示的平衡方程，又能满足位移边界条件以及用位移分量表示的应力边界条件。现在设想在变形体几何约束所允许的条件下，给它一个任意的微小变化，即所谓虚位移或位移变分 ， ， ，得到一组新的位移uvw第八章 能量原理及其应用195(a)wvu ,下面考察能量发生了什么变化。这时，外力在虚位移上所做的功(称为虚功)为(8.2-1a)dSZvYuXdVwZvYXWSV )()(或(8

6、.2-1b)FfiSiV式中， 为变形体的全部体积， 为变形体的全部表面积，其中给定外力的表面S记为 ，给定位移的表面记为露 。但面积分仅对给定面力的那一部分表面进Su行，对于给定位移的那一部分表面，因无虚位移，故不必考虑。应该指出，这里所说的虚位移般并不是由实际外力所引起的，而是由其他因素所引起，或者是为了分析问题而假想的。虚位移发生时，约束反力是不作功的，这是因为在约束力方向不可能产生位移。物体产生虚位移的过程中，物体必然产生微小的虚变形，因此在变形体中就产生虚应变能，即(8.2-2a) V zxyzxyzyx dVU )( 或写为(8.2-2b)VijdU假定变形体在虚位移的过程中，并没

7、有温度和速度的改变，因而也就没有热能和动能的改变。则按照能量守恒定律或热力学第一定律，应变能在虚位移上的增量 ，应当等于外力在虚位移所做的虚功，于是有U V zxyzxyzyx dV)( (8.2-3)SwZvYuXdVwZvYuXS()式(8.2-3)即为虚位移原理的位移变分方程，也称为拉格朗日(Lagrange)变分方程，有时也称为虚功方程。因此，虚位移原理可叙述为：在外力作用下处于平衡状态的可变形体，当给予物体微小虚位移时，外力在虚位移上所做虚功等于物体的虚应变能。现详细证明如下：若在虚位移原理的变分方程(8.2-3)中，因给定位移的部分表面 上 ，uS0i而在给定面力的部分表面 上，边

8、界条件 成立。则式(8.2-3)中右边对SijiFn的积分可以写为对整个物体表面 的积分，即有S第八章 能量原理及其应用196(8.2-4) SijiiViSiV dundufFdufW运用高斯散度定理，将上式对体积积分化为面积分，有(8.2-5)SwvmlzwyvxSzyxV )()()()(其中 为外边界法线方向单位矢量的方向余弦，即nml, ， ， ,cos(xl),cos(ym),cos(zn)注意到 SxVxVVxx udlduuddu )(以及 SxyVxyxyVxy dumvlduvduv )(其余的类似，因此由以上两式可得(8.2-6) SijijVij n,)(式中 。将式(

9、8.2-6)代入式(8.2-4)，有nmlj,(b)VjijiVji jijijiijVijiV duduffW,),( ),()当物体处于平衡状态时，因为 0,ijif所以式(b)中笫一项积分为零。又因 ),(21, ijjiijjiij u所以有 ijjiju,于是由式(b)得 VijdW将上式与式(8.2-2b)比较可知，有 U以上证明说明，当给予系统微小虚位移时，外力所作虚功与物体的虚应变能第八章 能量原理及其应用197相等是物体处于平衡状态的必要条件。另外，由应变位移关系以及先变分后微分与先微分后变分等价可知(c) ),()()(,uyvxyuxy 将式(c)代入(8.2-2a)，经

10、分部积分和利用格林公式，可得到如同下列形式的三个关系式(d) VxSx VxxVVx ududl dud)()()(以及下列形式的三个关系式(e)dVuyvxdSumvlyxVVxySyVyxy )()()( 将式(d)、(e) 所表示构六个关系式代人(8.2-2)式，则得(f) dVwzyx vzyxu dSwnmlnlnmlUzz yxzxyVS zzyxyzyzxV zxzy)( )( )()( ) 将式(f)代入式(8.2-3)，并加以整理，得 0)()( )()( dSwZnmlvYnml uXldVwZzyz vYzyxuXx zzyxyx SxzxyzV yxzxy 因为虚位移

11、各自独立，而且是完全任意的，因此上列积分式中括弧内的vu,系数均等于零，这样我们得到三个平衡方程和三个静力边界条件。因而证明是物体处于平衡状态的充分条件。UW第八章 能量原理及其应用198从以上讨沦可知，虚位移原理变分方程(8.2-3)式等价于平衡方程与应力边界条件。因此，满足变分方得(8.2-3)式的解就一定满足平衡方程和应力边界条件。所以，虚位移原理也可表述为：变形连续体平衡的必要与充分条件是，对于任意微小虚位移，外力所做虚功等于变形体所产生的虚变形能。应当指出，式(8.2-3)等号左边表示，由于产生虚位移 而引起物体内产生了虚iu应变能。这种虚位移实际上应理解为真实位移的变分，而不是其他

12、随便种位移函数。也就是说。式(8.2-3)中的虚应变 不是别的什么虚应变，而 xyx,是由 引起的，即它们之间满足下列条件wvu,(8.2-7),(21ijjiij u此外，位移 在己知位移边界 上还应满足 ，因此在己知位移边界 上i uSiiuS虚位移应为零，即(8.2-8)0i式(8.2-7)和(8.2-8)为方程(8.2-3)式的附加条件。因此，在应用虚位移方程式(8.2-3)时，所选取的解不必预先满足平衡力程和应力边界条件，但要求所给虚位移 能满足附加条件(8.2-7)式和(8.2-8)式，即应满足变形协调条件和wvu,几何边界条件。应当指出，由于虚位移原理的成立与材料本构关系无关，因

13、此，虚位移原理既适用于线弹性体、非线性弹性体，也适用于弹塑性体和理想塑料体等固体材料。例 8.1 如图 8.2 所示跨长为 ，抗弯刚度为 ，受分布荷重 作用的简支lEI)(xq梁，试用虚位移原理写出梁的挠曲线微分方程和边界条件。 解: 梁在平衡状态时，如果产生一虚 位移 ，由虚位移原理 w(1)WU此处 (2) dxAlx0由材料力学知，有 图 8.2 受均布荷重简支梁(3)“,Ezwzwxx 根据变分法则知(4)“)()(x将式(3)、(4)代入式(2)并整理后，得 llA dxEIdxzEU0“0“ )()(第八章 能量原理及其应用199对上式进行两次分布积分后，可化为(5)ll wdxw

14、EIU0)4()3(0“(外力所做虚功为(6)lxqW0)(将式(5)、(6)代入式(1)，则得(7)0)()( )3(0“0)4( lll wEIIwdxqEIw由于在支座处的虚位移应满足简支条件要求，所以边界条件为 0:0“lx考虑 的任意灶，于是要使式(7)成立，必有w(8)()4(xqEIw式(8)即为该梁的挠曲线微分方程。例 8.2 如图 8.3 所示受均布荷重 的简支梁，抗弯刚度为 ，跨中由弹性EI支座支承，试写出梁的边界条件和挠曲线微分方程。解 由例 8.1 知，变形能经两次分部积分后为ll wdxwEIU0)4()3(0“(2图 8.3 受均布荷重简支梁 令弹簧内的反力为 ，则

15、外力功为R(9)l cdxqW0式中 为梁在弹性支座 处的挠度。因 ，因此可由以上两式得cwC(10)02)(2)(2 0)3(0 0“4 cll l wEIwIdxqEI 于是，由式(10)可知，根据简支条件和对称条件，边界条件应为 ,)(,)( 0“0 xxlx同时，考虑到 除弹性支座处外，均为任意性，要使式(10)成立，必有w REIwlxclx)3(2,挠曲函数必须满足 0)4(qI第八章 能量原理及其应用200及。0,0,lxxw2.2 最小势能原理从位移变分方程(8.2-3)出发，可以导出虚功方程。假定物体从平衡位置有微小虚位移，物体的几何尺寸的变化略去不计，则原来作用在物体上的体

16、力和面力 的大小与方向都保持不变。于是，按照变分原理，式(8.2-3)中ZYX,ZY,变分的运算与积分的运算可以交换次序，故有(g)VSzxyzxyzyx dSwZvYuXdVZwYvXuV)()( 由第四章的(4.1-5a)知，有 ijijU0在式(g)中左边积分项中引入各向同性弹性体的广义虎克定律，并注意到 )(2122 zxyxzyx zzyzxyzyzxij 则有(h)20ijGU由第四章知，式中 。izyx在式(h)中代入应变位移关系，得(8.2-9) 222 22220 )()()( )(12),( xwzuvywxvyuGzwyvEwvuU当存在应变能 时，式(g)可写为0iU0

17、)()()( VSVi dSZvYXdVZYvXd于是有(8.2-10a)0)(WP也可写为(8.2-10b)VU其中第八章 能量原理及其应用201(8.2-10c)dSuFVufUiiVip )(0附加条件为(在 上) (8.2-10d)ii u式中 称为总势能， (=- )称为外力势能， 称为弹性变形体的应变势能。pWU当物体在不受外力作用的自然状态下，应变势能与外力的势能均为零。式(8.2-10)说明，在给走的外力作用下，实际的位移应使总势能的一阶变分为零，即使总势能取驻值。下面进一步让明有真实的位移总是使物体的总势能取最小值。对于稳定的平衡状态，物体偏离平衡状态而有虚位移时，其总势能的

18、增量恒为正。实际上可以让明，总势能量 的二阶变分为正。为此，令 为变形许可Piu的位移场， 为真实解的位移场，与之相应的应变张量分别为 和 ，于是当物iu jj体有虚位移时，有 ijijijiiu ,将 进行泰勒级数展开，并忽略二阶以上的高阶微量，可得)(0ijU(k)220000 )(1)()( ijijijijijij UU于是，变形许可状态的总势能与真实变形状态总势能之差为(l)dSuFdVufiSiVijijp )()(00因为(m)ppp2 1由式(8.2-9)知(n)00 dSuFdVufUiSiVp 比较式(k)，(m)可得 (q)VVklijklijp 020221当 足够小时

19、，式(q)必为正，因为如果令 ，则 ，则式(k)可化为ij ij 0ij2200)(1)(ijIJijU从而得 dVVijpijpijp )()()(2第八章 能量原理及其应用202由式(h)知， 为正定，所以)(0ijU(8.2-11)pp2,0上式表明如下一个原理：在给定外力作用下而保持平衡的弹性体，在满足位移边界条件的位移场中，真实的位移场使总势能取最小值。该原理称为最小势能原理。物体在外力作用下所产生的位移场，除了满足位移边界条件外，还必须满足以位移表示的平衡力程以及应力边界条件。最小总势能原理说明，真实的位移除满足几何边界条件外，还要满足最小势能原理的变分方程。实际上以上已经证明变分

20、方程(8.2-3)完全等价于平衡力程与应力边界条件。同样的结论也适用于式(8.2-9)。用最小势能原理和用泛定方程求解边值问题，只是形式上不同。以后将看到，这种解题手段的变更，在不少情况下将带来很大的力便，同时也扩大解题的范围。由最小总势能原理可导出熟知的卡氏(Castigliano)第一定理：当应变能 0U用广义位移表示为 对，则广义力 = 。 )(0iUiFiiU/)(0例 8.3 试由最小势能原理，弄略去剪应力影响，导出图 8.2 所示梁的挠曲线方程。解:根据应变能密度和虎克定律，梁的变形能为有(1)VxVdyzEd201由材料力学知，其中(2)Ax zdyI上dxwIMIy22,将式(

21、2)代入式(1)，经整理后得(3)lxEIU02)(1外力功为 lqwdW0根据最小势能原理 的变分量为 ，并注意到)(P )(x2)(2)( “ EIIEI因而第八章 能量原理及其应用203(4)l lp wdxqEI00“)(将上式等号右边第一项分部积分两次，可得(5) lll l dxEIIIwdxEI 0)4()3(0 0“)()( 对于简支端，有边界条件为(6)(0“0“ llEII将式(5)代入式(4)，并注意到(6)，得 )(04wdxqIl由于 任意性，因此有w 0)4(EI上式即为梁的挠曲线方程。8.3 位移变分法的应用基于虚位移原理的位移变分力程，提供了以位移作为基本未知数

22、的弹塑性力学问题的近似解法。瑞利里兹(Rayleigh-Ritz)和伽辽金(Galerkin)提出了各自解法，现分别介绍如下：3.1 瑞利 里兹法当给定面力和几何约束条件时，可以利用位移变分方程求解。因此时应力边界条件和位移边界条件为已知，由虚位移原理或最小势能原理所导出的变分方程(8.2-3)和(8.2-10a)均等价于平衡方程和应力边界条件，所以采用式(8.2-3)和(8.2-10a)求解肘，所选取的位移函数不需要先满足应力边界条件，只需满足位移边界条件。设位移函数为(8.3-1)nkknkzyxwcvbvzyxuau1010),(,),(式中， 为未知的待定的常数， 满足边界条件，即在已

23、知位移边kcba, 0,vu第八章 能量原理及其应用204界 上，应有uS wvu000,而 为坐标线性独立的识定函数，且在)21)(,),(),( nkzyxwzvzyxkkk 己知位移边界 上满足u ),21(, nkzyxvkkk 这样，无论 如何取值，位移函数总是能满足位移边界条件。n由于 是设定的已知函数，因此对位移进行一阶),(),(),(zyxzkkk变分时，只需对系数 取一阶变分，即cba(a)nk nkkknkk czyxwbzyxv上u1 11 ),(, 将式(8.3-1)代入虚位移原理变分方程(8.2-3)或最小势能原理变分方程(8.2-9a)，由 的任意性，可得确定全部

24、系数 的线性代数方程组。例如，cba, kca,将式(8.3-1)代入式(8.2-10a)，可得 0)(1 S kkkVkkknkk cZwbYvaXu dVcZwbYvXuUU 式中系数 的变分是完全任意的，彼此无关。于是，将上式整理后得kcba,(8.3-2)SkkVk kkkSkkVk dwZcUvYdbduXa由应变能函数表达式(8.2-9)及位移分量表达式(8.3-1)可知，应变能 应是待定U系数 的二次函数，因而式(8.3-2)将是各个待定系数的线性方程组，共有kcba,3 个方程。从方程组(8.3-2)可解全部系数 后，即可由式(8.3-1)求得位n kcba,移分量。式(8.3

25、-2)也可写为(8.3-3)00kpkkpcb),21(n第八章 能量原理及其应用205这一方法称为瑞利-里兹法，也称为里兹法。选择合适的 和 ，以及项数 ，可以获得精确度较高的位移0,wvukkv,n解。将求得的位移代入用位移表示的应力表达式，在计算应力分量时需对位移求导，通常近似解的精度往往会因求导而降低，因此应力近似解的精度一般都较差。这是因为应力分量并不精确地满足平衡方程，只是满足平衡方程与一个加权函数乘秋的积分为零的条件，即iu 0)(),( VSijiijiP dSuFnduf要提高精度，只有增加式(83.3-1)中位移函数的项数 ，当项数 时，则其n解将收敛于精确解。3.2 伽辽

26、金法如果选择的位移因数表达式(8.3-1)，不仅能满足位移边界条件，还能满足应力边界条件，那么变分方程(8.2-3)或(8.2-10a)为 WU或(8.3-4) SVijiji jiVVSVjiiijij dudndFufdU,0注意到所取位移函数满足应力边界条件，因此由上式得 0),(fiVji即(b)0 dVwZzyx vYzyxuXzz yVxzxy 将式(b)展开为三个方程，则每个方程均含有 个积分，并注意到变分关系式(a)n和 为任意值，所以耍使式(b)成立，则只能每个积分式均等于零，于kkcba,是可得第八章 能量原理及其应用206(8.3-5) Vkzzyzx kyzyxVkxz

27、xydVwZvYdVuX00 ),21(nk对于各向同性弹性变形体，将以上三个方程中的应力分量，通过广义虎克定律方程(4.2-14b)、几何方程(3.2-9)转换用位移分量表示，可得(8.3-6)nkdVwZzEvYyduXxV kkV k ,21021)( 021)(222 式中 。由式(8.3-1)可知，位移分量 是系数 的线性wyvxu vu,kcba,函数，所以式(8.3-5)和式(8.3-6)将是这些系数的线性方程组。求解此方程组，可解得 个系数，从而由式(8.3-1)求得位移分量。这个方法称为伽辽金法。n3比较以上两种基于虚位移原理的近似计算方法可知，在位移函数的选择上，伽辽金法比

28、瑞利-里兹法更为严格，它不仅要满足位移边界条件，还必须满足应力边界条件；但在应用上伽辽金法比较方便，因为可不必导出泛函。仅根据熟知的平衡方程就能列凼伽辽金方程。例 8.4 如图 8.4 所示受均布荷重悬臂梁，梁跨长为 ，抗弯刚度为 ，试根据lEI初等理论，不计体力和用瑞利-里兹法和伽辽金法分别求梁的挠度和固端弯矩。解：(1) 瑞利- 里兹法求解当不计剪切对挠度的影响，现设挠曲线方程为图 8.4 受均布荷重悬臂梁 (1)2cos1()lxaxw上式满足固定端条件第八章 能量原理及其应用2070,0)(xdw由初等理论知，梁的弯矩为 lalEIdxIM2cos)()(4梁的变形能为(2)23420

29、022 6cos41)(1 alEIdxlalIxEIUl l对于本问题，应用瑞利-里兹法，位移函数表达式(8.3-1)仅保留第三式，又因所选位移函数式(1)满足边界条件，所以相应的线性方程组(8.3-2)也仅保留第三式。在不计体力的条件下，该式成为(3)dxlqaUl0)2cos1(将式(2)代入式(3)后，并求解可得(4)(34EIl将式(4)代入式(1)，得梁的挠曲线表达式为(5)2cos1)(32)(4lxIqlxw最大挠度发生在 处，即lx EIqlEIql 44max19.0)(32这个结果与材料力学解解 误差仅为 4.5%。IlIlw4ax5.8利用 lxqldxEIM2cos)

30、1()(22可得 220max 945.0)(8qllx第八章 能量原理及其应用208与精确解 相比，误差达 41%。2205.qllMX如果取挠曲线函数为 )23cos1()2cos1() lxalxaxw此式也能满足前述边界条件。 llxlEIxMcscs)( 21 21230 212 73ocos alEIdlalalIUl 根据瑞利-里兹法 ll dxlqalEIall0232131 )23cos(7解以上联立方程，可得参数为 )31(27)1(4241 EIla上EIql故挠曲线为 lxlxIqlxw 2cos1327cos123)(4则最大挠度和最大弯矩分别为 2max4318.0

31、)(9qlMEIlw由上式可见，无论挠度还是弯矩，其精度均有所提高，但弯矩提高不大明显。如果将挠曲函数的项数增加，则随项数的增大，无论是挠度值还是应力值均会逐步接近精确值。(2)伽辽金法解依据伽辽金法位移函数不仅要满足位移边界条件，而且还要满足静力边界条件，因此设想取挠曲线函数为第八章 能量原理及其应用209(6)(234BxAxaw由梁固定端位移边界条件 ，可知满足位移边界条件。同时满足0,)0xd静力力边界条件。再由自由端的静力边界条件， ，可解得0)(“lwl26,4BlA因此，梁的挠曲线函数为(7)()234xlxaw由伽辽金方程(8.3-5)知，因在 方向无外载荷，又由初等理论知，第

32、一、y,二式自动满足，因此(8.3-5)式化为 0dVZxkVz即 l kzSdxwqxF00)(式中 系梁中的剪力, 。再由材料力学知zSF)64(23lwk，所以上式化为)()3(xEIw(8)0)( 23404( dxlxqEIl将式(6)代入式(7)，并积分后得Ia2将它代入式(6)，最后得挠曲线表达式、是大挠度和最大弯矩分别为 )64()( 23xlxEIqwIl8max2)0(2axqldwM这个结果与材料力学结果完全相同。第八章 能量原理及其应用210第八章 能量原理及其应用2118.4 虚应力原理和最小余能原理虚位移原理是从位移变分出发可直接求出位移分量。但在工程实际问题中，往

33、往感兴趣的是直接得到表征结构强度的应力分量。而位移变分法得到的位移分量，必须通过几何方程和本构方程求出烹形体内的应力分量，在计算过程中因需经多次微分往往会产生较大的误差。为此，真接以应力分量为未知数求解变形体问题的应力法便具有重要的价值。同时，对一些特殊的问题，例如，平面问题，柱体扭转等，可以引进应力而数，此时应力法更具有极大的方便。4.1 虚应力原理对于处于平衡状态的变形体应用变分原理时，取虚位移 ，即对位移分量iu进行变分，这些位移的变分必须是几何上可能的。为此，引入虚应力概念。所谓虚应力是满足平衡方程及指定的力边界料的、任意的、微小的应力。虚应力记为，即对应力分量进行变分，这些应力的变分

34、必须是静力上可能的，即经变分ij后，新的应力分量必须满足平衡方程和应力边界条件。设 为实际存在于变形体内的应力分量，则它应该满足平衡方程，应力边ij界条件和应力协调条件。现让这些应力分量发生静力许可的微小变化，得到新的应力分量为(a)ijijij它们必须满足平衡方程和应力边界条件，于是新的平衡方程为(b) 0)()()( 0)()()( Zzyz Yx Xzyzz yy xxx 式中 为给定的体力，设想没有改变。将式(B)与发生应力变分前的平衡方ZYX,程相减，得(8.4-1)0)()()( 0)()()(zzyzx yzyyx xzxyx第八章 能量原理及其应用212变形体的表面分为两帝分，

35、即给定面力的部分表面 和给定位移部分的表面 。SuS在表面力没有给定的边界 上，由于应力分量的变化，表面力也发生相应的变化，uS即 。新的表面力变成ZYX, ZYX,因此，新的应力分量在此边界面上应满足边界条件，即(c)Znml Xl zzzyzzxz yy xzzxyx )()()( )()()(将式(c)与原边界条件式相减，得(8.4-2)Znml YXlzzyzxyyxzxyx此外，对于表面力已给定的边界 上，表面力不能变化，即S0YX故应力变分在该边界上应满足的条件为(8.4-3)0zzyzxyyxzxyxnmll因此，为了使应力变分是静力许可的，它必须满足式(8.4-1)、(8.4-

36、2)和式(8.4-3)。于是，应力余能的变分，参照(8.1-4)式或(8.2-2a)可写为(d) V zxyzxyzyx dVU )( 将几何方程代入上式，得(e)V xyxvuu)( 将上式内各项进行分部积分并利用格林公式，得如同下列形式的三个关系式(f)dVxudSulxdxuSVxVV )()()( 和如同下列形式的另外三个关系式第八章 能量原理及其应用213(g) VxyxysxyxyV dVuvdSumvldVuv )()()()( 将式(f)、(g)代入式(d)，得(g)dVzyxwyxv udSnlwnlU zzyz xzxyxVxzyzxS zyyxz )()()()()()(

37、 )()( 由式(8.4-1)知，上式中对体积的积分项为零；并注意在已知表面力的边界 上，S应满足式(8.4-3)，而在已知位移的边界 上，因 ，且在该边界uSwv,因虚应力产生的附加面力应满足式(8.4-2)。因此，由式(8.3-2)可知，上式可简化为(f)dSZwYvXUuS)( 该变分方程等式的右边部分表示表面外力的增量 与实际位移YX,所作的功 ,同时注意式(d)，因而式(f)可写为wvu,(8.4-4)dSZwYvXudVSV zxyzxyzyx )( ) 式(8.4-4)表示虚应力原理，又称虚功原理，表述为：当物体在已知体力和面力作用下处于平衡状态时，微小虚面力在实际位移所作的虚功

38、，等于虚应力在真实应变所产生的虚应变余能。显然，式(8.4-4)成立的附加条件为式(8.4-1)和式(8.4-3)。该两式可简写为)(0),(上SFVij 由以上讨论可以看出，虚应力原理和虚位移原理在形式上是互补的。和虚位移原理一样，虚应力原理的成立也与材料的本构关系无关。应当指出，在虚位移原理中包含了实际的外力和内力，因而可理解为虚位移原理是对系统平衡的要求；而虚应力原理则包含有实际的位移和应变，所以可把虚应力原理看作是对物体变形协调的要求。实际上，由虚应力原理的变分方程(8.4-4)不难导出变形协调方程，这就是说式(8.4-4)等价于应变协调条件。于是，按式(8.4-4)解题时，对于所设解

39、答，不必预先满足变形协调条件，只须使虚应力 满足物体的平衡方程和应力边界条件。ij第八章 能量原理及其应用2144.2 最小余能原理由虚应力原理可直接导出最小总余能原理。为避免混乱，今后把用应变表示的弹性应变能函数 称为应变能函数，或应变能；而把用应力表示的应变余)(ijU能函数称为余应变能函数，或余变余能(或应力能)，记为 。如在虚应力原)(ijU理中引进广义虎克定律，并认为应变状态是有势的，应变分量可由余应变能函数导出，即 ijij)(0而 ijU0于是由上式可知，总的应变余能的变分为 VVijij dd)(0因此，式(8.4-4)可化为(8.4-5a)0)()(0 SZwYvXudUSV

40、ij 如果存在虚应力时，在边界 上，位移分量应保持不变。于是可将上式中的变分符号置于积分号外，即(8.4-5b)0)()(0 VSiijudFd其中 是在已知位移边界由虚应力引起伪附加表面面力。显然，在此情况下附加iF条件为(8.4-6)(0,上SFnVfijiiji 如果令变形体的余能为 uSiijVPdFU)(0则有 P上式说明，在所有满足平衡方程和应力边界条件的静力许可的应力场中，真实的应力场使余能取极值。进一步还可证明 02P因此，最小余能原理可表述为：在所有满足平衡方程和应力边界条件的静力许可的应力场中，真实的应力场使余能取最小值。第八章 能量原理及其应用215最小余能原理和最小势能

41、原理均适用于线性和非线性弹性体。由于真实的应力场既满足平衡方程，应力边界条件，又满足变形协调条件。可见，最小余能原理因真实的应力场满足平衡方程，应力边界条件，以及使余能取最小值的条件，所以最小余能原理与变形协调条件等价。下面指出最小余能原理的特殊情形。当物体全部表面力给定，则面力的变分为零，由式(8.4-5)得(8.4-7)0UP式(8.4-7)称为最小功原理。该原理可表述为：若物体的面力给定，则在所有满足平衡方程和边界条件的应力场中，真实的应力场必使余应变能取最小值。寸于线弹性体，因余应变能与应变能相等，因此又称(8.4-7)式为最小应变能原理。当最小余应变能原理用于线弹性力学问题则可导出熟知的卡氏第二定理。8.5 应力变分法与应用基于与位移变分法类似的思想，通过应力变分方程，以应力分量作为基本未知数，取得变形物体的近似解答。5.1 应力变分法应力变分法是设定某一个应力分量表达式，其中包含了若干待定常数，使其满足平衡方程和应力边界条件，然后通过应力变分方程决定这些常数。帕普考维奇(，.)建议取应力分量为(8.5-1)nk nkzxzxkzzk kyzyzkyynk nkxyxykxx AA1 1001 100,其中， 是选定的满足平衡方程和应力边条件的设定函数，而 ,0xy是选定的满足体力为零的