1、第 十 四 章 杆件的应变能及其应用一、教学目标和教学内容1教学目标让学生掌握杆件弹性应变能的有关概念。理解和掌握在工程力学有广泛应用的能量方法。掌握功能原理、功的互等定理、位移互等定理、卡氏定理。能够熟练地计算基本变形杆件和常见的组合变形杆件的应变能。对于简单结构应变能，也能够完成应变能的计算。能够较为熟练地应用卡氏第二定理，完成杆件的位移计算，并可以求解简单超静定问题。为进一步在结构力学等后续课程中，学习和应用能量方法奠定基础。2教学内容介绍能量法的有关概念。例如，外力的功、应变能、比能等等。介绍基本变形杆件应变能计算和组合变形杆件应变能计算。讲解功能原理、功的互等定理和位移互等定理。讲解
2、余能概念和卡氏定理。二、重点难点重点：建立应变能等有关概念。基本变形杆件和常见的组合变形杆件的应变能的计算。卡氏第二定理及其应用。难点：杆件应变能计算中的可否叠加问题。对于广义力和相应广义位移的正确理解和认识。应用卡氏第二定理求位移时，如何正确地选取或设定与位移相应的广义力。能否正确写出内力方程，灵活地进行先求偏导数再积分的运算。三、教学方式采用启发式教学，通过提问，引导学生思考，让学生回答问题。四、建议学时6 学时五、讲课提纲1、弹性应变能与功能原理弹性体在荷载作用下将发生变形，外力作用点要产生位移因此，在弹性体的变形过程中，外力沿其作用方向做了功，称为外力功。对于弹性体，因为变形是可逆的，
3、外力功将以一种能量形式积蓄在弹性体内部。当将荷载逐渐卸除时，该能量又将重新释放出来作功，使弹性体恢复到变形前的形状。例如钟表里的发条在被拧紧的过程中，发生了弹性变形而积蓄了能量，在它放松的过程中可带动指针转动，从而发条就作了功。弹性体伴随弹性变形积蓄了能量，从而具有对外界作功的潜在能力，通常把这种形式的能量称为弹性应变能（Dlastic strain energy）或弹性变形能（Dlastic deformation energy） ，用 表示。 V根据物理学中的功能原理，积蓄在弹性体内的应变能 及能量损耗V在数值上应等于荷载所作的功，即EWEV如果在加载过程中动能及其它形式的能量损耗不计，应
4、有（14.1）利用上述的这种功能概念解决固体力学问题的方法统称为能量法，相应的基本原理统称为功能原理（Principle for work and energy） 。弹性体的功能原理的应用非常广泛，它是目前在工程中得到广泛应用的有限单元法的重要理论基础。2 、杆件的应变能计算如前所述，若外力在加载过程中所作的功全部以应变能的形式积蓄在弹性体内，即在加载和卸载的过程中能量没有任何损失，则只要得到加载过程中外力功的数值，弹性体应变能的数值也就可以计算出来，所以说外力功是应变能的一种度量。2.1 外力功的计算外力作功分为以下两种情况。一种情况为常力作功。 这里所谓常力，是指工程动力学中，作用在不变形
5、的刚体上使刚体产生运动的力。当外力在作功过程中保持不变时，它所作的功等于外力与其相应位移的乘积。例如，在沿外力 方向线上有线位移 ，则FFW另一种情况为静荷载作功。所谓静荷载，是指构件所承受的荷载从零开始缓慢地增加到最终值，然后不随时间改变。所以静荷载的施加过程均为变力。静荷载作功，可以解释为在其施加过程中的一种变力作功。例如图 14.1 所示的简单受拉杆，拉力由零逐渐增加到定值 ，由 产生的伸长变形由零逐渐增加F到 ，这就是拉力 的作用点的位移。如果材料服从胡克定律，则外力 与位l F移 成线性关系（图 14.2 ） 。设 表示加载过程中拉力的一个值，相应的位a1移为 ，此时将拉力增加一微量
6、 ，使其产生相应的位移增量 ，这时，1l d)(1ld已经作用在杆上的拉力 将在该位移增量上作全功，其值为1F（14.2）dW)(1l图 14.1在上式中略去了 在 上作的功，这部分功为二阶微量。 在图 14.21dF)(1ldW中以阴影面积来表示。拉力从零增加到 的整个加载过程中所作的总功则为a F这种单元面积的总和，也就是 的面积，即OABF01)(1ldl2上述积分是与静荷载施加过程有关的积分，可以称为静荷载作功的过程积分。积分结果的系数 1/2，既是已经完成过程积分的标志，又表示构件材料为线性弹性材料。将以上的分析推广到其它的受力情况，因而静荷载下外力功的计算式可写为（14.3）FW2
7、1图 14.2式中的 是广义力，它可以是集中力或集中力偶； 是与广义力 相对应FF的位移，称为广义位移，它可以是线位移或角位移。上式表明，当外力是由零逐渐增加的变力时，在符合胡克定律的范围内，外力在其相应位移上所作的功，等于外力最终值与相应位移最终值乘积的一半。2.2 杆件的应变能计算2.2.1 应变能的有关概念按照功能原理，应变能可以由计算外力的功得到，这是应变能的一种计算方法。同时，也表明线弹性材料杆件的应变能，在完成了过程积分，也始终具有1/2 系数。应变能和外力的功，它们在杆件受力变形过程中的积累，也可以由荷载伸长图和应力应变图(见图 14.2)考察到。2.2.2 杆件的应变能计算2.
8、2.2.1 杆件在各种基本变形时应变能的计算如前所述，应变能是根据能量守衡原理通过外力功来计算的。以下我们讨论的均为静荷载问题，动能和其他能量的损耗不计。1. 轴向拉伸或压缩杆的应变能及比能当拉（压）杆的变形处于线弹性范围内时，外力所作的功为lFW21则杆内的应变能为lV由图 14.1 知，杆件任一横截面上的轴力lFldVWl 21)(01 lFldWVl 21)(0121 01dlFldWVl21)(01FN考虑到胡克定律有EAlN所以，拉（压）杆的应变能为（14.4 ）EAlFVN2 a或 （14.4 ）l)(2 b若外力较复杂，轴力沿杆轴线为变量 ，可以先计算长度为)(xFN微段内的应变
9、能，再按积分的方法计算整个杆件的应变能，即dx=dVEAxFN2)(（14.5）lNdx)(2为了对构件的弹性变形能有更全面的了解，我们不但要知道整个构件所能积蓄的应变能，而且要知道杆的单位体积内所能积蓄的应变能。对于承受均匀拉力的杆（图 14.1） ，杆内各部分的受力和变形情况相同，所以每单位体积内积蓄的应变能相等，可用杆的应变能 除以杆的体积 来计算。这种单位体积内的应VV变能，称为应变比能（Density of strain energy） ，简称比能，并用 表示，v于是=vVAlFN21可见应变比能 的数值也可以用 图中 的面积来表示（图 14.2 ） 。vOabb根据胡克定律 ，比能
10、又可以写成下列形式 E（14.6）v221E2剪切变形时的应变能及比能为了分析的方便，从受剪切杆中截取如图 14.3 所示的单元体，该单元体a处于纯剪切应力状态，假想其在一个面（如左侧面）上被固定起来，则在剪应力由零逐渐增加到 值的过程中，单元体将发生如图所示的变形，与此对应的剪应变由零增加到 值，其右侧面向下的位移为 。当材料在线弹性范围 dx内工作时，其 与 成正比（图 14.3 ） ，与图 14.2 、 中所示受拉杆的相应bab图形类似。所以，单元体各表面上的剪力在单元体变形过程中所作的功为)(21dyzW)(21dxyz上式中，作功的力是单元体右侧面上的剪力。由于剪应变 很小，其余各面
11、上的剪力，在其作用方向上没有位移，都没有在其作用方向上作功。故单元体内积蓄的应变能为图 14.3dVdVW21单元体内积蓄的应变比能则为vdV21这表明， 等于 直线下的面积。由剪切胡克定律 ，比能又可以 G写成下列形式 v（14.7）221G3. 圆轴扭转时的应变能及比能如图 14.4 所示的受扭圆轴，若扭转力偶矩由零开始缓慢增加到最终值 ，a T则在线弹性范围内，相对扭转角 与扭转力偶矩 间的关系是一条直线（图T14.4 ） 。与轴向拉伸杆件相似，扭转圆轴的应变能应为b图 14.4VTW21由于圆轴横截面上的扭矩 ，且MxPxGIl所以，受扭圆轴的应变能为（14.8）PxIlV2l2实际上
12、，受扭圆轴中各点的应力状态均为纯剪切应力状态，因而可以直接采用公式（14.7） ，求积分即得杆件的应变能。因为剪应力 ，PxIM所以VlAdxGv2APxPxGIlMdIGl222当扭矩 沿轴线为变量时，上式变为xlPxGIdV2)(（14.9）可见利用比能计算全杆内积蓄的应变能应用范围更广，该方法适用于杆各横截面上内力变化（相应横截面上各点处的应力也不同）的情况。 4. 弯曲变形时的应变能及比能（1）纯弯曲梁设如图 14.5 所示的简支梁在两端的纵向对称平面内受到外力偶a作用而发生纯0M弯曲，在加载过程中，梁的各横截面上的弯矩均有 = ，故梁在线弹性M0范围内工作时，其轴线弯曲成为一段圆弧（
13、图 14.5 ） ，两端横截面有相对的转a动，其夹角为图 14.5l且EIM01故EIl0与前面的情况相似，在线弹性范围内，当弯曲外力偶矩由零逐渐增加到时，梁两端截面 上 相 对 转 动 产 生 的 夹 角 也 从 零 逐 渐 增 加 到 ， 与 的 关 系0M0M也 是 斜 直 线 （ 图 14.5 ） ，所以杆件纯弯曲变形时的应变能为bWV（14.10） EIlM202（2）横力弯曲梁 在工程实际中，最常遇到的是受横力弯曲的梁（如图 14.6 所示） 。这时，a梁横截面上同时有剪力和弯矩，所以梁的应变能应包括两部分：弯曲应变能和剪切应变能。由于剪力和弯矩通常均随着截面位置的不同而变化，都是
14、 的函x数，因此，计算梁的应变能应从分析梁上长为 的微段开始（图 14.6 ） 。dxb在弯矩的作用下，微段产生弯曲变形，两端横截面有相对的转动（图 14.6） ；在剪力的作用下，微段产生剪切变形，两端横截面有相对的错动（图 14.6c） 。由于在小变形的情况下，弯曲正应力不会引起剪应变，剪应力也不会引起d线应变，或者说，由弯矩产生的位移与由剪力产生的位移互相垂直，因此，可以先分别计算出弯矩和剪力在各自相应的变形位移上所作的功，然后根据叠加原理将它们叠加起来。但由于在工程中常用的梁往往为细长梁，与剪应力对应的剪切应变能，比与弯矩对应的弯曲应变能小得多，可以不计，所以只需要计算弯曲应变能。图 1
15、4.6微段梁左右两端横截面上的弯矩应分别为 和 。在计)(xM)()(xd算其应变能时，弯矩增量 所作的功为二阶微量，可忽略不计，因此可将)(xdM该微段看作是纯弯曲的情况。应用式（14.10）可求得微段的弯曲应变能EIdxdV2)(全梁的弯曲应变能则可积分上式得到（14.11）lIdxM2)(如 果 梁 中 各 段 内 的 弯 矩 由 不 同 的 函 数 表 示 ， 上 列 积 分 应 分 段 进 行 ， 然)(后 再 求 其 总 和。由以上各种变形形式下应变能的计算式可以看出，应变能是力的二次函数，也是变形的二次函数。当构件同时受几个力（或力偶）作用时，能否用叠加原理求应变能？2.2.2.
16、2 复杂受力情况下应变能的计算1有关应变能的两个重要概念。（1）是否可以应用叠加原理计算应变能。下面以图 14.7 所示的拉杆为例加以说明。拉杆在 、 同时作a 1F2用下的应变能为（14.12）EAlFlEAFV22)(11 而当 、 单独作用时（图 14.7 、 ） ，杆的应变能分别为1F2 bc图 14.7和 EAlFV21 EAlFV2显然21可见对图 14.7 所示的情况不能用叠加原理计算应变能。其原因是这些荷a载所作的功是互相影响的，即荷载除在其自身引起的位移上作功外，在其它荷载引起的位移上也要作功，所以不能将各荷载单独分析再进行叠加。这样的荷载称为属于同类型荷载。例如若先将 作用
17、在拉杆上，杆件有伸长 ，则 所1F1lF作的功为112lW在 不卸除的情况下，再施加 ，杆件又伸长了 ，故变力 、常力1FF2l2F所作的功为和 221l213lW则整个加载过程外力所作的功为32121212lFllF将上式转化为应变能则同样得到式（14.12） 。其中 就是两力所作功互相3W影响的结果。（2）应变能是否与加载次序及过程有关。对于上述的拉杆，若先施加 再施加 ，通过类似的计算可以证2F1明，杆件内积蓄的应变能与上述分析结果一样，当然也与 、 同时作用时一F2样。可见，积蓄在弹性体内的弹性应变能只决定于弹性体变形的最终状态，或者说只决定于作用在弹性体上的荷载和位移的最终值，与加载
18、的先后次序无关。2组合变形时的应变能如果作用在杆件上的某一荷载作用方向上，其它荷载均不在该荷载方向上引起位移，则前一荷载与其它荷载将属于不同荷载类型，则仍可应用叠加原理计算应变能，即是说，可以单独计算前一荷载作用下杆件的应变能，单独计算其它荷载作用下杆件的应变能，然后叠加得出杆件的总的应变能。组合变形时的应变能就属于这种情况。图 14.8 所示的微段杆是从处于拉、弯、扭组合变形下的圆杆中取出的，其长度为 ，横截面上的轴力 、弯矩 和扭矩 均只在各自引起的位dx)(xFN)(xMx移 、 和 上作功，各类荷载所作的功互相没有影响，故微段杆内的应)(l变能可用叠加原理计算，即图 14.8dWV )
19、(21ldxFNdxM)(21x)(21PxNGIdMEIAF)()()(2整个圆杆的应变能则为+lNEAdxFV2)(（14.13）lPxl GIdMEId2)()(22.5 应变能的普遍表达式以上讨论了杆件在基本变形和简单组合变形下应变能的计算，现在研究更普遍的情况。设有 个广义力 、 、 作用在如图 14.9 所示的物体上，n1F2n且设物体的约束条件足以使它只会发生由于变形引起的位移，不会发生刚体位移。 、 、 表示荷载沿各自作用方位上的广义位移（图 14.9） 。由前12n面的分析我们已经知道，弹性体在变形过程中积蓄的应变能，只决定于作用在弹性体上的荷载和位移的最终值，与加载的先后次
20、序无关。于是，不管实际加载的情况如何，在计算应变能时，为计算方便起见，可以假设这些荷载按同一比例从零开始逐渐增加到最终值，则弹性体的应变能等于各广义力在加载过程中所作功的总和，即图 14.9（14.14）inidFVi10当作用于弹性体上的荷载与其相应位移之间的关系是线性时，即物体为线弹性体，则应变能的计算式为（14.14）niiFV12这表示线弹性体的应变能等于各荷载与其相应位移乘积的二分之一的总和。这一结论称为克拉贝依隆原理。3、 功的互等定理和位移互等定理由前面的讨论可知，对线弹性体结构，积蓄在弹性体内的弹性应变能只决定于作用在弹性体上的荷载的最终值，与加载的先后次序无关。由此可以导出功
21、的互等定理和位移互等定理。它们在结构分析中有重要应用。功的互等定理又称互等功定理，是意大利的 E.贝蒂（E.Betti）1872 年和英国的瑞利（Rayleigh）1873 年分别独立提出的，所以又称贝蒂.瑞利互等功定理。位移互等定理又称互等位移定理， 是英国的 J.C.麦克斯韦（J.C.Maxwell）于 1864 年提出的，又称麦克斯韦位移互等定理。 下面以一处于线弹性阶段的简支梁为例进行说明。图 14.10 、 代表梁的ab两种受力状态，1、2 截面为其上任意两截面。如图 14.10 所示， 使梁在截1F面 1、2 上的位移分别为 和 ；在图 14.10 中，当 作用时，在截面12b21
22、、2 上的产生的位移则分别为 和 。在位移符号的角标中，第一个表示截2面位置，第二个指是由哪个力引起的。图 14.10现在用两种办法在梁上加载，来计算 、 共同作用时外力的功。1F2先施加 再施加 时（图 14.10 ） ，外力的功1F2a12211 FW而当先施加 再施加 时（图 14.10 ） ，外力的功2 b211221由于杆件的应变能等于外力的功，与加载次序无关，即 ，所以21WV有（14.16）2112F这表明，第一个力在第二个力引起的位移上所作的功，等于第二个力在第一个力引起的位移上所作的功。这就是功的互等定理（Reciprocal Theorem of Work） 。图 14.1
23、0当 时，由（14.16）式可推出一个重要的推论，即21F（14.17）21这表明，作用在方位 1 上的荷载使杆件在方位 2 上产生的位移 ，等于将21此荷载作用在方位 2 上而在方位 1 上产生的位移 。这就是位移互等定理1（Reciprocal Theorem of Displaceme.nt） 。若令 =1（即为单位力） ，且此时用 表示位移，则有21F21由于 1、2 两截面是任意的，故上述关系可写为以下一般形式jiij即 处作用的单位力在 处产生的位移，等于 处作用的单位力在 处产生ji ij的位移。这是位移互等定理的特殊表达形式，在结构分析中十分有用。以上分析对弹性体上作用的集中力
24、偶显然也是适用的，不过相应的位移是角位移，所以上述互等定理中的力和位移泛指广义力和广义位移。4、余能概念及卡氏第一定理4.1 余能 以上推出的公式均只在线弹性范围内成立。下面，我们进一步讨论非线性弹性体的应变能表达式，并介绍非线性弹性体的应变余能（简称余能）概念及表达式。我们仍以图（14.1）所示的拉杆为例，但材料非线性弹性的，这时力 与F相应的位移 的关系就是非线性的（图 14.11 ） 。对比图 14.2 ，不难看出仍aa可用下式计算外力作的功（14.18）FdW图 14.11 表示，外力功的大小与位移从 0 到 之间一段 曲线下的面积相aF当。式（14.18）是以位移作为积分变量的，若以
25、力作为积分变量，则有（14.19）FCd称为余功（Complementary work） 。余功没有具体的物理意义，但有明确的几何意义，那就是外力从 0 到 之间一段 . 曲线与纵坐标轴间的面积。从图14.11 不难看出，功和余功互补为常力功。a图 14.11由于材料是弹性的，如果将加载和卸载过程中的能量损耗略去不计，则同样有与线弹性体类似的结论，即积蓄在弹性体内的应变能 在数值上应等于V外力所作的功（14.20）FdW同样地，余功 与余应变能 在数值上也相等，即CCV（14.21）F此即为由外力余功来计算余应变能（Complementary strain energy）的表达式。如果从拉杆中
26、取出一个边长为 1 的单元体，该单元体处于单轴应力状态，上、下表面的力为 ，对于单元体而言它们是外力。与 相应的 F伸长量为 ，于是在对拉杆加载过程中，作用在单元体上的外力功1l为dw该外力功在数值上等于积蓄在单元体内的应变能，即比能 ，于是v（14.22）v同样地，若以应力作为积分变量，则有（14.23）dC式中 称为余应变比能（density of complementary strain energy） ，其大Cv小就代表 曲线与纵坐标轴间的面积（14.17 ） 。 b例如，材料的应力应变关系为 时，物体的应变比能和余应E变比能分别为032Edv23E0231EdvC材料在拉伸压缩时的应
27、力应变关系可写成 ，显然，这时应变能和余应变能在数值上是相等的，对于线弹性体，当变形在线弹性范围内时，应变能和余应变能在数值上也是相等的。但应该注意，余功、余应变能如前述都没有明确的物理意义，只是因为它们具有与外力功一样的量纲，才把它们作为一种能量参数，而在求解非线性弹性问题时，它们非常有用。4.2 卡氏第一定理卡氏定理（Castiglianos Theorem） ，是意大利工程师卡斯蒂利亚诺（A.Castigliano ）于 1873 年提出的，故得其名。我们以图 14.9 所示的弹性体为例来说明。设弹性体上作用有 个广义力 ，与niF这些力对应的广义位移为 ，其中 、2、 ，如果将弹性体的
28、应变能i1i表示为位移的函数 、 、 ，则应变能函数对某个广义位移的V1(V2n)偏导数，等于与该位移相应的广义力，即（ 、2、 ） （14.24）iiFi这就是卡氏第一定理（First Castigliano Theorem） 。下面进行证明。现假设沿第 个作用力方向的位移有一微小增量 ，则弹性体的应变i id能 有相应的增量为V（ ）dVi a这时弹性体内的应变能为（ ）i b式 中 ： 可 由 应 变 能 的 普 遍 表 达 式 （ 14.13） 计 算 ， 代 表 应 变 能 对 于 位 移 的V iV i变 化 率。此外，由于只有沿第 个作用力方向的位移有一微小增量，沿其余作i用力方
29、向无位移变化，故外力功的增量为（ ）idFWc前面我们已经知道，外力功在数值上等于应变能，它们的变化量也应相等，即dWV将（ ） 、 （ ）两式代入上式中，则有ab（ 、2、 ） iiVF1in由上述证明过程可见，卡氏第一定理同时适用于线性弹性体和非线性弹性体。 5、 卡氏第二定理及其应用5.1 卡氏第二定理下面仍以图 14.9 所示的弹性体为例来说明。根据余能的计算公式（14.21） ，并仿照应变能的普遍表达式（14.14） ，弹性体的余能 可写成下列CV形式：（14.25）niFiCidV10它是外力的函数 、 、 。(2n)现假设第 个广义力有一微小增量 ，则弹性体的余能有相应的增量i
30、idF为（ ）iiCdFV a此外，由于除 外其余外力均维持常量不变，故外力余功的增量为iF（ ）iCdWb则由CV可得到（14.26）iCiF利用上式可计算非线性弹性体在广义力 作用方位上与 相应的广义位移 。i iFi对于线性弹性体，正如前节所述，此时的应变能与余能相等，即，则式（14.26）可改写为CV（ 、2、 ） （14.27）iiFV1in这就是卡氏第二定理（Second Castigliano Ttheorem） 。应用卡氏第二定理计算线性弹性体的位移很简便。5.2 卡氏第二定理的应用在第二节中我们已经得到了各种基本变形及组合变形情况下应变能的计算式，这些式子中的内力均为外力的函
31、数，分别代入式（14.27） ，便可得到各种基本变形及组合变形情况下计算位移的卡氏第二定理的应用式如下：（1）轴向拉伸或压缩杆iiFVi lNEAdx2)(dxFiNl)(（2）扭转圆轴iiFVi lpxGIM2)(dxFMGIxilP)(（3）平面弯曲梁=iiFVi lGIdxM2)(lidxFMEI)(（4）组合变形杆件+ +iiVdxFEAiNl)(lixPxdFMGI)(lidxFMEI)(（5）简单桁架结构。由于桁架的每根杆件均受均匀拉伸或压缩，若桁架共有根杆件，故niiFVniiNiAEl1根据卡氏定理的表达式，我们知道在计算结构某处的位移时，该处应有与所求位移相应的外力作用，如果这种外力不存在，可在该处附加虚设的外力 ，从而仍然可以采用卡氏定理求解.F
