胡伟平 詹志新 孟庆春

(北京航空航天大学航空科学与工程学院固体力学研究所，北京 100191)

逐段分析求和法[1](或称为逐段刚化法[2]或逐段软化法[3])是梁变形计算中的一种重要方法。现有材料力学教材关于该方法的介绍主要有三种方式，第一种是只明确提出叠加法而不提逐段分析求和法[4-7]，仅把逐段分析求和法作为叠加法的一种实现方法。第二种是两种方法都明确提出来，但是对于这种方法的来源以及理论基础没有详细说明[2]。第三种是将该方法统一归于梁变形计算的叠加法一节，并明确指出逐段分析求和法与叠加法的理论基础其实并不相同，叠加法的理论基础是力的独立作用原理，而逐段分析求和法的理论基础是梁局部变形与总体位移之间的几何关系[1,3]，但未说明这个几何关系到底是什么，以及这种几何关系与具体应用时的方法是怎样联系起来的。一般来说，同学们对于力的独立作用原理都很清楚，但是对于梁局部变形与总体位移之间的几何关系，理解起来就比较模糊了。另外，教师在介绍逐段分析求和法时，往往偏重于该方法的具体应用[8]，忽略对这种方法的力学本质的分析，以至于同学们很容易认为这是一种全新的方法，并只注重采用这种方法做习题，而不去思考该方法背后的思想和原理。本文试图从钓鱼竿这一生活中的实例入手，通过对钓鱼竿的变形分析自然引入逐段分析求和法，并由变形分析导出具体计算时采用的方法，即逐段刚化法或逐段软化法。进而讨论应用逐段分析求和法的前提条件，以及其与叠加法的本质区别。最后，通过从离散逐步逼近连续的数学思想，阐明逐段分析求和法与积分法的内在联系。

1 钓鱼竿挠度计算的基本力学模型

钓鱼竿是大家非常熟悉的一种娱乐或生活用具，如图1(a) 所示。尽管一般不会去思考钓鱼竿的变形问题，但都会注意到钓鱼竿的结构形式。用日常生活用语描述就是一节一节的，用材料力学术语描述就是典型的变截面梁。因此，钓鱼竿的力学模型可以表示为如图1(b) 所示结构。

图1 钓鱼竿及其力学模型

当钓鱼竿钓上鱼的时候，鱼竿发生明显的弯曲变形，如图2(a)所示。此时，钓鱼竿的力学模型可以简化为变截面梁在自由端受铅垂集中力F，如图2(b)所示。观察钓鱼竿的变形，能够直观认识到，钓鱼竿端头的位移与鱼竿每一段的弯曲变形都相关，这便是教材上所说的“梁的总体位移与局部变形之间的关系”。但是，到底总体位移与局部变形之间是什么关系，以及如何来描述这种关系，这就需要进行定量的分析与推导了。

图2 钓上鱼的钓鱼竿及其对应的变形分析力学模型

为了便于进行问题的分析与说明，首先研究两段变截面梁的变形问题，其力学模型如图3(a)所示。按照一般的教学安排，在梁变形计算部分，首先介绍的是积分法，然后是叠加法。介绍这两种方法时都是以等截面均质梁为研究对象的。因此，变截面梁的变形计算对于同学们来说是一个新的问题。在材料力学的学习过程中，从变形观察入手来进行力学问题的分析是一个重要的方法。通过对变形规律的定性分析和描述，能够很好地启发定量建立相关方程。因此，在并不明确采用何种方法进行变截面梁的变形分析时，可以直观地画出两段变截面梁的变形图，如图3 (b) 所示。图中为了方便起见，定义A1A0为第1 段，A2A1为第2 段。

图3 钓鱼竿变形分析的基本力学模型与变形示意图

从图中可以看出，A0截面的挠度可以分成三部分，即wa，wb和wc，并且有

式中，wc为第2 段弯曲变形时A1截面挠度所引起的A0截面挠度，wb为第2 段弯曲变形产生的A1截面转角所引起的A0截面挠度，wa则为第1 段自身弯曲变形所引起的A0截面挠度。即，wc和wb为第2 段的变形效果，wa为第1 段的变形效果。于是，式(1)可以描述为：变截面梁A0截面的挠度等于各段变形效果的总和。将这一描述推广到一般情况，即为：梁总体位移等于各局部变形效果的总和。这就是梁总体位移与局部变形之间的几何关系。

2 钓鱼竿基本力学模型的计算方法

上一节建立了钓鱼竿基本力学模型，即两段变截面梁模型，并通过变形分析得到了变截面梁总体挠度与各段变形之间的关系，即由式(1) 给出的总体位移等于各段变形效果总和。按照这一思路，只需求出各段变形效果，问题就得以解决。但是在梁的真实变形中，两段的变形总是同时发生的，如何才能将两段变形效果分开进行计算。或者说，在什么情况下，变截面梁的挠度仅与其中一段的变形效果相关。按照这一思路，可以得到如图4(a)和图4(b)所示力学模型，即对梁每一段的变形效果进行分解。

图4 变截面梁每一段变形效果的分解

图4(a) 中，A2A1段视为刚体，在这种情况下A0截面的挠度仅包含A1A0段的变形效果。图4(b)中，A1A0段视为刚体，此时，A0截面的挠度仅包含A2A1段的变形效果。这一方法被称为逐段刚化法[2]，也有教材称之为逐段软化法[3]。这两种名称都比较直观形象，易于理解，只是各自的参考构型不同而已。逐段刚化法是以弹性梁为参考构型，通过刚化其他段来分析某一段的变形效果，而逐段软化法是以刚性梁为参考构型，通过软化某一段来求该段的变形效果。

通过比较图4 与图3(b)，还可以得到应用逐段分析求和法的前提条件。也就是说，在何种情况下，图3(b)所示梁的变形与图4 所示两种变形效果的总和是相等的。很显然，只有在原始构型下进行受力和变形分析时，这种变形分解才成立。因此，小变形条件是逐段分析求和法的前提条件。需要说明的是，逐段分析求和法并没有线弹性的限制条件。从图4 可以看出，这个分解过程完全是基于变形几何的，与材料的物理特性无关。这是逐段分析求和法与叠加法的本质区别。关于这一点，文献[9] 从能量法的角度进行了相关讨论。此处结合一个具体例子进行简单说明。

图5 为一非线性悬臂梁，梁材料的应力σ与应变ε满足关系式σ=kε1/2，k为刚度，求A0截面挠度。

图5 非线性悬臂梁的位移分析

根据式(6)，采用逐段分析求和法分别求解图5(b)刚化A1A0段时A0截面的挠度，和图5(c) 刚化A2A1段时A0截面的挠度。图5(b) 梁A0截面的挠度为(2M0)2a2/(2Λ)+(2M0)2a/Λ·a，图5(c) 梁A0截面的挠度为M20a2/(2Λ)。相加得到A0截面的总挠度为13M20a2/(2Λ)，与直接根据式(4) 和式(5) 计算所得结果一致。

若采用叠加法，计算仅在A1截面作用力偶矩M0时A0截面挠度为M20a2/(2Λ)+M20a/Λ·a，以及仅在A0截面作用力偶矩M0时A0截面挠度为M20(2a)2/(2Λ)，相加得总挠度为7M20a2/(2Λ)，结果错误。

3 钓鱼竿总位移的计算

前面建立了钓鱼竿基本力学模型以及具体的计算方法，现需将基本模型得到的结果应用于钓鱼竿整体模型中。从整体模型中任意截取两段，如图6(a)所示，其变形形式可以用如图6(b) 所示的基本模型来描述，即相当于图3 所示基本模型中固支端存在挠度和转角。

图6 钓鱼竿基本模型与整体模型的关系

对于图6(b) 所示基本模型，可以写出Ak−1截面挠度和转角的表达式

式中，右边第一项就是各段自身的变形效果，第二项就是各段之间的转角连续条件引起的附加效果。

4 逐段分析求和法与积分法的内在联系

通过前面的分析，对逐段分析求和法的方法来源和物理本质已经有了比较清楚的认识，下面来探讨逐段分析求和法与积分法的关系。

首先，在线弹性情况下，积分法是基于挠曲轴的二阶近似微分方程，通过积分得到梁截面挠度和转角的方法，是进行梁变形分析的最基本方法。其基本方程为

式中的积分常数C和D需要根据位移边界和连续条件确定。

为了研究逐段分析求和法与积分法之间的关系，需要对图3 所示基本力学模型进行简单的数学处理。考虑到积分法的研究对象是等截面梁问题，而阐述逐段分析求和法时研究对象是分段的变截面梁问题，可以通过将离散问题连续化的方法建立这两种模型之间的关系。假设图3(a) 中的变截面梁每一段长度为dx，则可以得到图7(a) 所示力学模型，该模型的变形图如图7(b) 所示。

图7 连续变化的变截面梁的力学模型

由图7(b)，可以写出

将式(12) 移项，变为

进一步，在平面假设的前提下，挠曲轴切线转角与梁的截面转角相等。因此，式(15) 正是积分法中挠度和转角的关系方程。因此，在线弹性小变形条件以及平面假设成立的情况下可以将逐段分析求和法理解为积分法在分段的变截面梁或者分段的变刚度梁时的具体体现。

将这种内在关系推广到非线性情况。若梁为非线性材料，在小变形和平面假设成立的前提下有

式中，ρ为梁中性层曲率半径。此时梁中性层曲率半径就不能再表示成M(x)/(EI) 的形式了，需要根据材料的物理方程进行具体推导。但是，在小变形和平面假设成立的前提下，截面转角与挠度之间的导数关系仍然成立，而基于图7 推导出来的式(12)～式(15)亦成立，其并不需要利用梁变形与内力之间的物理关系。因此，在小变形和平面假设成立的前提下，对于非线性情况，逐段分析求和法亦等价于挠度和转角的积分关系在分段的变截面梁或者分段的变刚度梁时的具体体现。

5 结论

本文以钓鱼竿为例，通过建立基本力学模型，从变形分析入手，自然引入逐段分析求和法，进而得到在计算每一段变形效果时的逐段刚化法或逐段软化法。并讨论了应用逐段分析求和法的前提条件，以及与积分法的内在关系。可以总结如下。

(1)逐段分析求和法是基于梁总体位移等于各局部变形效果的总和这一几何关系建立的。

(2) 应用逐段分析求和法的前提条件是小变形，并不要求结构为线弹性。这是与叠加法最本质的区别。

(3)对于线弹性小变形情况，在平面假设成立的前提下，逐段分析求和法与积分法是等价的，逐段分析求和法是积分法在分段变刚度梁时的具体体现。对于非线性材料，当满足小变形条件和平面假设时，逐段分析求和法亦等价于挠度和转角的积分关系在分段的变截面梁或者分段的变刚度梁时的具体体现。