第2章-应变分析修改课件

2023/9/4周书敬第二章应变分析1在静力学理论中，通常假定物体是刚性的，即在力的作用下，构成该物体质点之间的距离保持不变。前章建立平衡条件时，就忽略了固体变形，即假定固体是刚体。实际上刚体是不存在，所有物体在某种程度上都是可以变形的，也即是说，在力的作用下，实际物体质点之间的距离总是要发生变化的。一个物体是否可以被假定为刚体，关键在于刚体假定的有效范围。本章从几何学的观点出发分析研究物体的变形。反映物体变形规律的数学方程也有两类，即几何方程和变形协调方程。由于这两类方程都是基于物体连续性的假定从几何学出发得到的，并不涉及产生变形的原因和物体的材料性质，所以它们均属于“普适方程”。2023/8/3周书敬第二章应变分析1在2023/9/4周书敬第二章应变分析2前面讨论了受力物体的应力，现在开始讨论物体的变形。在外力作用下，物体各点的位置要发生改变，即发生位移。如果物体各点发生位移后仍保持各点间初始状态的相对位置，则物体实际上只产生了刚体移动和转动，将这种位移称为刚体位移。如果物体各点发生位移后改变了各点间初始状态的相对位置，则物体就同时产生了形状的变化，统称该物体产生了变形。（书图2－1）2023/8/3周书敬第二章应变分析2前2023/9/4周书敬第二章应变分析3第一节

一点的应变状态应变与位移的关系

为了确定正应变的定义，在一受拉杆上有线段AB，在变形后，变为（见右图）。若线段AB的长度为，变形后的A点的物体不论是发生空间的刚体运动或实形状的变化，终归体现为物体内部每一点产生位移；因而，只要确定了物体内各点的位移，物体的变形状态也就确定了。因此研究物体内一点的变形是很重要的。2023/8/3周书敬第二章应变分析3第一节一点的应2023/9/4周书敬第二章应变分析4下面我们讨论一般情况，给出应变的概念。设在直角坐标系中，变形前A点的坐标是（x，y，z），变形后的坐标是（x+u，y+v，z+w），这里u，v，w是A点的位移在x，y，z三轴上的投影，它们都是坐标x，y，z的连续函数，而且位移的导数也是连续的。定义：正应变（2－1）显然，如果变形的分布是均匀的，则有:即：材料力学的拉伸应变。（2－2）位移是u，而B点的位移是u+

u，则线段增加了

u。2023/8/3周书敬第二章应变分析4下面我们讨2023/9/4周书敬第二章应变分析5

设由变形体中取出一个微小六面体（见书中图2－3变形体的投影），在研究微小六面体的变形时，采用的分析方法是将六面体的各面投影到直角坐标系的各个坐标平面上，研究这些平面投影的变形，并根据这些投影的变形规律来判断整个平行六面体的变形。由于变形很微小，所以可以认为两个平行面在坐标面上的投影只相差高阶的微量，因而，两个平行面的投影可以合并为一个投影面。2023/8/3周书敬第二章应变分析52023/9/4周书敬第二章应变分析6首先，研究平行六面体在xoz面上的投影ABCD（见书中图2－4）。在变形前六面体A点的坐标为（x，y，z），在六面体变形时，投影上的A点移到了点，同时而整个ABCD移到。设A点的位移是u，w，它们是坐标的函数，因此有：

（2－3）

2023/8/3周书敬第二章应变分析6首先，研究2023/9/4周书敬第二章应变分析7而B点的坐标为（x+dx，y，z），因此B点在x方向的位移为：

根据泰勒级数展开式，可得：

略去高阶项后得到：（2－4）

由于则AB在x轴上的投影的伸长量为，则有：

2023/8/3周书敬第二章应变分析7而B点的坐标为（x2023/9/4周书敬第二章应变分析8同理可得平行于y轴和z的边长的正应变，因此有：（2－5）

下面研究六面体的剪应变，即各直角的改变。（见图）

取变形前的直角BAC或，变形时，棱边转动一个角度，棱边转动一个角度，在xoz平面内，角应变用表示，其值为和之和，即：

（2－6）

若A点在z轴方向的位移为，

当大于零时，表示线段伸长，反之表示缩短。2023/8/3周书敬第二章应变分析8同理可得平行于y2023/9/4周书敬第二章应变分析9图：位移矢量在xoz平面上的投影返回2023/8/3周书敬第二章应变分析9图：位移矢量在xo2023/9/4周书敬第二章应变分析10则B点在Z轴方向的位移为，

B点与A点沿Z轴方向的位移之差为:

在直角三角形中，可得：

在分母中（）与1相比是一个微量，故可以略去，因而得出，

2023/8/3周书敬第二章应变分析10则B点在Z轴方2023/9/4周书敬第二章应变分析11同理可得：

所以有剪应变：同理可得另外两个剪应变。即有剪应变的表达式（2－7）

（2－7）

说明：剪应变的正负号2023/8/3周书敬第二章应变分析11同理可得：所以2023/9/4周书敬第二章应变分析12所以，正应变和剪应变的表达式为（2－8）：（2－8）

式（2－8）称为柯西（Cauchy）几何关系。[式(2－8)的提出者：法国工业学院的数学教授柯西（Cauchy）(1789－1857)，于1822年发表的论文提出的]注意：书中P48给出了帮助记忆的图形（图2－5）。可知：如果已知位移分量可以很简单的求出应变分量；反之，则问题比较复杂。2023/8/3周书敬第二章应变分析12所以，正应变和剪2023/9/4周书敬第二章应变分析13利用类似的方法，可以导出柱坐标表示的几何方程为式（2－9）：

（2－9）2023/8/3周书敬第二章应变分析13利用类似2023/9/4周书敬第二章应变分析14其中，分别表示一点位移在径向（r方向），环向（方向）以及轴向（z方向）的分量。

对于平面问题，柱坐标变为极坐标，则平面极坐标表示的几何方程为：

（2－10）

下面给出式（2－10）的推导过程。

2023/8/3周书敬第二章应变分析14其中，2023/9/4周书敬第二章应变分析15首先假定只有径向位移而没有环向位移：

如图（2－6）所示，在P点沿径向和环向取两个微段PA和PB，设PA移到了，位移为u；PB移到了，则P，A，B三点的位移分别为：径向位移图2023/8/3周书敬第二章应变分析15首先假定只有径向2023/9/4周书敬第二章应变分析16则PA的正应变为：

PB的正应变为：

径向线段PA的转角为：

环向线段PB的转角为：

所以有：

2023/8/3周书敬第二章应变分析16则PA的正应变为2023/9/4周书敬第二章应变分析17其次，假定只有环向位移而没有径向位移:

见图2－7，由于P点的环向位移v，径向线段PA移段到了，环向线段PB移到了，则P，A，B三点的位移分别为：

可见：径向线段PA的正应变为：图2-7环向位移图2023/8/3周书敬第二章应变分析17其次，假定只有环2023/9/4周书敬第二章应变分析18环向线段PB的正应变为：

径向线段PA的转角为：

环向线段PB的转角为：

2023/8/3周书敬第二章应变分析18环向线段PB的正2023/9/4周书敬第二章应变分析19所以剪应变为：

因此，如沿径向和环向都有位移，则根据叠加原理可得式（2－10）。

对于轴对称问题：，，则式（2－10）的平面极坐标几何方程为（2－11）

（2－11）

对于球对称问题：变形的几何方程为式（2－12）

（2－12）

2023/8/3周书敬第二章应变分析19所以剪应变为：2023/9/4周书敬第二章应变分析20注意：书中P47对方程（2－10）的相关项进行了解释，自己看一下。2023/8/3周书敬第二章应变分析20注意：2023/9/4周书敬第二章应变分析21第二节

应变状态分析

现在已知物体内任一点P的六个应变分量，试求经过该点（P点）的沿N方向的任一微小线段PN＝dr的正应变，以及经过P点的微小线段PN和的夹角的改变。

令PN的方向余弦为l、m、n，则PN在坐标轴上的投影为：

2023/8/3周书敬第二章应变分析21第二节应变2023/9/4周书敬第二章应变分析22（2－13）

设P点的位移分量为u，v，w，则N点的位移分量为：略去高阶项（小量）得：

同理可得：即有式（2－14）2023/8/3周书敬第二章应变分析22（2－13）设2023/9/4周书敬第二章应变分析23（2－14）

在变形后，线段PN在坐标轴上的投影为（2－15）式：即

（2－15）

2023/8/3周书敬第二章应变分析23（2－14）在2023/9/4周书敬第二章应变分析24令线段PN的正应变为，则该线段变形后的长度为：而且有

（2－16）

上式两边同除以，并利用（2－13）式得：

2023/8/3周书敬第二章应变分析24令线段P2023/9/4周书敬第二章应变分析25因为和位移分量的导数都是微小的，它们的平方和乘积可以不计，可得：

利用，上式可得：

再利用几何方程可得：

（2－17）

2023/8/3周书敬第二章应变分析25因为2023/9/4周书敬第二章应变分析26下面来求PN和的夹角的改变

设PN在变形后的方向余弦为，则由式（2－13）和式（2－15）可以得到：

注意到，都是微小量，在展开上式后，略去二阶以上的微小量得：

2023/8/3周书敬第二章应变分析26下面来求PN和2023/9/4周书敬第二章应变分析27同理可得出，即得出式（2－18）

（2－18）

与此类似，设线段在变形之前的方向余弦为，则其在变形后的方向余弦为：2023/8/3周书敬第二章应变分析27同理可得出2023/9/4周书敬第二章应变分析28（2－19）

（2－20）

其中，是的正应变。

令PN和在变形之前的夹角为，变形之后的夹角为，则有：

2023/8/3周书敬第二章应变分析28（2－19）（2023/9/4周书敬第二章应变分析29将式（2－18）和（2－19）代入，并略去高阶微量可得：

利用几何方程，并注意到，则有：（2－21）

由此可求出，进而可求得。

2023/8/3周书敬第二章应变分析29将式（2－18）2023/9/4周书敬第二章应变分析30由此可见：在物体内的任一点，如果已知六个应变分量，就可以求出经过该点的任一线段的正应变，也可以求得经过该点的任意两线段之间的夹角的改变。这就是说，六个应变分量完全决定了这一点的应变状态。

2023/8/3周书敬第二章应变分析30由此2023/9/4周书敬第二章应变分析31第三节主应变在研究一点的应力状态时，可以找到三个相互垂直的没有剪应力作用的平面，将这些面称为主平面，而这些平面的法线方向称为主方向。在研究应变问题时，同样可以找到三个相互垂直的平面，在这些平面上没有剪应变，将这些面称为应变主平面，而这些平面的法线方向称为应变主方向。对应于该主方向的正应变称为主应变。2023/8/3周书敬第二章应变分析31第三节主应变2023/9/4周书敬第二章应变分析32一点的应变状态也可以用张量表示，这时引进符号

（2－22）（书：2－13）

则应变张量为：

（2－23）（书：2－14）

通常称为“工程剪应变”2023/8/3周书敬第二章应变分析32一点的应2023/9/4周书敬第二章应变分析33应变张量还可以写为：

式中的不同符号可以交换使用，这就要看在某些特定用途中哪个用起来更方便。2023/8/3周书敬第二章应变分析33应变张量还可以写2023/9/4周书敬第二章应变分析34下面分析如何确定主应变：

在直角坐标系空间中取一微小线段，设A点在x方向的位移为u，则有B点在x方向的位移为：

图2-8略去高阶微量得：

显然（或由全微分概念）有：

2023/8/3周书敬第二章应变分析34下面分析如何确定2023/9/4周书敬第二章应变分析35进一步可写成式（2－24）（书：2－15）

（2－24）（书：2－15）

这里要注意的是：当一个物体从一个位置变形到另一个空间位置（图2－9）时，其中可能包括一部分刚体位移（平动或转动），而这部分位移不引起形变，其实式（2－24）中的和恰恰表示物体的微小刚性转动。（下页图）

图2-92023/8/3周书敬第二章应变分析35进一步可写成式（2023/9/4周书敬第二章应变分析36B点的三部分位移一般来说，对于可变形固体而言，与物体内任一点A无限临近的一点B的位移有三个部分组成：1、随同A点的一个平动位移，如图中的所示；2、绕A点的刚性转动在B点所产生的位移，如图中的所示；3、由于A点临近微元体的形状变化在B点引起的位移，如图所示，这部分位移与应变张量分量有关。2023/8/3周书敬第二章应变分析36B点的三部分位移2023/9/4周书敬第二章应变分析37因此，当考虑纯变形时有：

（2－25）（书：2－16）

如果用张量表示，则为

其中，j称为“哑标”（表示求和）。

现在取一微小四面体O123（图2－10），为法线方向，设斜面123上只有正应变（即主平面），则有：2023/8/3周书敬第二章应变分析37因此，当考虑纯变2023/9/4周书敬第二章应变分析38并且一定为要求的主应变。

（成比例是因为与方向一致）

（书：2－17）

代入式（2－25）（书：2－16）得出：（书：2－18）

（2－26）（书：2－18）

2023/8/3周书敬第二章应变分析38并且一定为2023/9/4周书敬第二章应变分析39若上式有非零解，必须有“系数行列式为零”，可得：

（2－27）（书：2－19）

其中，为应变第一、二、三不变量，且有：

（2－28）（书：2－20）

2023/8/3周书敬第二章应变分析39若上式有非零解，2023/9/4周书敬第二章应变分析40若方程式（2－27）可以因式分解，则应有：

式中，为主应变。用主应变表示的应变不变量将为：

（书：2－20）

在主应变平面上，剪应变为零。则由方程（2－27）可以求出三个主应变。

2023/8/3周书敬第二章应变分析40若方程式（2－22023/9/4周书敬第二章应变分析41增例1：已知物体中任意一点的位移分量如下式表示，试比较点A（1，2，3）与点B（0.5，-1，0）的最大伸长值（绝对值）。解：利用几何方程求得应变分量为：2023/8/3周书敬第二章应变分析41增例2023/9/4周书敬第二章应变分析42点A的应变分量值为：应变不变量为：该点的主应变值可由下式确定，即为计算方便，令代入上式，得2023/8/3周书敬第二章应变分析42点A的应变分量2023/9/4周书敬第二章应变分析43以代入上式，消去二项式，得此方程的解为：由此得A点的主应变为：故点A的最大伸长的绝对值为可以用获得的三个主应变之和是否等于第一应变不变量的值，检验所得结果是否正确。2023/8/3周书敬第二章应变分析43以2023/9/4周书敬第二章应变分析44用同样的方法可以求得点B的主应变为：故点B的最大伸长的绝对值为由以上计算可知，点A最大伸长值大于点B的最大伸长的绝对值。2023/8/3周书敬第二章应变分析44用同2023/9/4周书敬第二章应变分析45增例2：已知物体中某点的应变分量为：试求该点的主应变方向。解：首先计算应变不变量，并解三次方程，求得主应变值为为求解主应变方向，利用下列方程组：2023/8/3周书敬第二章应变分析452023/9/4周书敬第二章应变分析46将代入上式，第一式自然满足，其余两个方程式为以上两式的唯一解为。为满足，则有。即的方向余弦为（1，0，0）。2023/8/3周书敬第二章应变分析46将2023/9/4周书敬第二章应变分析47将代入方程组，得由第一式得。由二、三式可得。再由得，由该式求得，而。即的方向余弦为（0，0.585，0.811）。同样可求得的方向余弦为（0，-0.811，0.585，）2023/8/3周书敬第二章应变分析47将2023/9/4周书敬第二章应变分析48第四节

应变张量和应变偏量

仿照应力张量分解，应变张量可以分解为与体积变化有关的“球形应变张量”和与物体形状变化有关的“应变偏量”。利用书中（2－14）式可以分解为：

其中球形应变张量为：

（2－30）（书：2－22）

一、应变张量的分解2023/8/3周书敬第二章应变分析48第四节应变2023/9/4周书敬第二章应变分析49应变偏量可写为：

式中，为平均正应变。

其中，，，称为“应变偏量分量”。可写为：

2023/8/3周书敬第二章应变分析49应变偏量2023/9/4周书敬第二章应变分析50（2－32）（书：2－23）2023/8/3周书敬第二章应变分析50（2－32）（书2023/9/4周书敬第二章应变分析51若用主应变表示应变偏量，则有式（2－33）（书：2－24）（2－33）（书：2－24）

三个坐标平面为应变主平面在主应变为坐标的应变空间中有：由应变偏量张量的定义式（书2-23）可见，它是一个实对称二阶张量，因此，存在三个主值及其相应的主方向。可以证明，应变偏量张量的主方向与应变张量的主方向一致，而且它的主值e1，e2，e3与应变张量的主应变存在如左的关系。2023/8/3周书敬第二章应变分析51若用主应变表示应2023/9/4周书敬第二章应变分析52注意：纯剪应变状态的条件与纯剪应力状态的条件相同，即纯剪变形的必要且充分条件是，因此，为纯剪状态且与有相同的主轴。同样，应变偏量增量也存在三个不变量，它们分别表示为：当用张量给出一点的应变状态时，需注意2023/8/3周书敬第二章应变分析52注意：纯剪应变状2023/9/4周书敬第二章应变分析53其三次方程为：二、体积应变在考虑塑性变形时，经常采用“体积不变”假设，这时球形应变张量为零，应变偏量等于应变张量，即“应变分量与应变偏量的分量相等”，这一假设，对于简化计算来了方便。现在我们来研究每单位体积的体积改变，即体积应变。

2023/8/3周书敬第二章应变分析53其三次方程为：二2023/9/4周书敬第二章应变分析54设有微小的正平行六面体，它的棱边长度是：

变形前它的体积为：

变形后它的体积称为：因此，它的体积应变为：

对于小应变（忽略高阶微量）有：

验证体积不变假设的成立2023/8/3周书敬第二章应变分析54设有微小2023/9/4周书敬第二章应变分析55（2－34）

由此则有：显然，若体积不变，则必有球形应变张量为零成立，且有。在主应变空间：

对于小应变有：

2023/8/3周书敬第二章应变分析55（2－34）由2023/9/4周书敬第二章应变分析561、主剪应变（工程主剪应变）

（2－35）（书：2－25）

三、相关结论

与应力分析类似。在应变分析中也有一些相应的公式，下面给出有关结论：如果，则最大剪应变为：（2－36）（书：2－26）

2023/8/3周书敬第二章应变分析561、主剪应变（工2023/9/4周书敬第二章应变分析57（1）等倾面（或称八面体面）的剪应变为，则有：

（2－37）（书：2－27）

2、八面体应变（正应变、剪应变）对任意一组坐标轴x，y，z的应变分量的八面体剪应变可写为：2023/8/3周书敬第二章应变分析57（1）等倾面（或2023/9/4周书敬第二章应变分析58单向拉伸情况：可得

此时的应变张量为：

平均应变为:3、单向拉抻时的应变

（2）等倾面（或称八面体面）的正应变为，则有：

（三个主应变的平均值）2023/8/3周书敬第二章应变分析58单向拉伸情况：可2023/9/4周书敬第二章应变分析59应变偏量的分量为：

（书：2－28）

球形应变张量为：

2023/8/3周书敬第二章应变分析59应变偏量的分量为2023/9/4周书敬第二章应变分析60应变偏量为：

在以主应变为坐标轴的主应变空间内讨论。4、应变强度（等效应变）

（2－39）（书：2－30）当体积不可压缩时，令，称为应变强度或等效应变。

这里之所以不称为应变强度，而又引进符号，是因为要与应力分析中的情况相一致。

2023/8/3周书敬第二章应变分析60应变偏量为：2023/9/4周书敬第二章应变分析615、应变率应变率：在变形过程中，单位时间中应变值的增量称为“应变率”。即：

（2－41）（书：2－31）

根据小变形的几何关系，可得应变率分量：即：应变率分量等于位移率分量对相应坐标的偏导数，也等于应变分量对时间的偏导数。见书中解释。

2023/8/3周书敬第二章应变分析615、应变率2023/9/4周书敬第二章应变分析62（2－42）（书：2－34）

增加：关于应变率的推导

在小变形的条件下，设物体内任一点速度在坐标轴上的投影为：2023/8/3周书敬第二章应变分析62（2－42）增加2023/9/4周书敬第二章应变分析63其中,此处用到“小变形”假设：即在小变形条件下，（a）物体内各点的位置坐标因变形而有的改变可以忽略不计，即初始位置与瞬时位置坐标可以不加区别；（b）此处，还可以略去物体内各点的位移梯度分量的影响，用对时间的偏导数代替全导数。应变率分量用符号表示为：则有：类推得（书中2-34）式，与应变张量相似，可得应变率张量。2023/8/3周书敬第二章应变分析63其中,此2023/9/4周书敬第二章应变分析64在塑性力学中经常使用应变增量的概念。实验证明，静力学中塑性变形规律和时间因素是没有关系的，此处dt并不代表时间，因此，用应变增量来代替应变率往往更能表示塑性静力学应变不受时间参数影响的特点。即：通常使用的不是应变率张量，而是在时间步长或dt内的应变增量。应变增量：

6、应变增量2023/8/3周书敬第二章应变分析642023/9/4周书敬第二章应变分析65有了应变增量的概念，则可描述应变成比例变化或不成比例变化时的规律。7、应变强度增量

在比例变形的条件下，在应变空间是一条直线，而当各应变分量的变化不成比例时，在三个坐标轴互成120°角的坐标中，将是一条折线或曲线，如图2－17（书P63）。2023/8/3周书敬第二章应变分析65有了应变2023/9/4周书敬第二章应变分析66如果变形过程是两条折线组成的，则这一情况表示两个比例变形的过程。当变形由曲线表示时，则可认为变形过程的总应变强度等于曲线的总长度，即应变强度不仅与初始应变状态和最终应变状态有关，而且还与应变历史即变形过程有关。

由图2－17可见，应变强度增量与应变增量分量，和有关。（2－43）（书：2－36）

[注：的表达式中，只有简单加载条件下才有]2023/8/3周书敬第二章应变分析662023/9/4周书敬第二章应变分析67此式即为应变强度增量的表达式，它是各应变分量增量的函数。在塑性力学中，当应变较大时，需采用另外一种表示应变的方法。

8、工程应变有一截面为而长度为的受拉构件，在某一时刻其长度达到而截面积为，且杆件伸长量为，则应变增量及应变的表达方法如下：

工程应变：假设两质点相距，变形后为，则有工程应变表达式（2－44）：

2023/8/3周书敬第二章应变分析672023/9/4周书敬第二章应变分析68（2－44）

对数应变：设某瞬时的应变增量为，积分后得到对数应变的表达式（2－45）：

（2－45）（书：2－37）

（2－46）（书：2－38）

显然有对数应变和工程应变之间的关系为：2023/8/3周书敬第二章应变分析68（2－44）2023/9/4周书敬第二章应变分析69截面收缩率：

（2－47）

其中A0为初始时截面面积，A为某一时刻的截面面积。

若材料为不可压缩，则有：

（

）

或

（2－48）

不同应变指数之间的关系见书中表2-2（P65）2023/8/3周书敬第二章应变分析69截面收缩率2023/9/4周书敬第二章应变分析70增例3：给定一点的应变张量计算：（a）主应变、和；（b）最大剪应变；（c）八面体应变和。解：（a）计算应变不变量，求主应变。2023/8/3周书敬第二章应变分析70增例3：给定一点2023/9/4周书敬第二章应变分析71特征方程变为或求得三个主应变为2023/8/3周书敬第二章应变分析71特征方程变为或求2023/9/4周书敬第二章应变分析72校核：用、和的值代入三个不变量的表达式，以校核所得结果。（b）计算最大剪应变。（c）八面体应变和。2023/8/3周书敬第二章应变分析72校核：用2023/9/4周书敬第二章应变分析73增例4：一点的应变状态由给定的应变张量表示确定：（a）应变偏量张量；（b）应变偏量不变量和；（c）单位体积的体积变化（膨胀）。2023/8/3周书敬第二章应变分析73增例4：一点的应2023/9/4周书敬第二章应变分析74解：（a）计算平均应变。所以有应变偏量张量为：（b）计算不变量。2023/8/3周书敬第二章应变分析74解：（a）计算平2023/9/4周书敬第二章应变分析75（c）单位体积的体积变化（膨胀）。即在该应变张量表示的应变状态下，该点附近体元的体积减小。2023/8/3周书敬第二章应变分析75（c）单位体积的2023/9/4周书敬第二章应变分析76第五节应变协调方程（连续性方程、

相容方程）

（Equationsofcompatibility）

在研究物体变形时，一般都取一个平行六面体进行分析，物体在变形时，各相邻的小单元不能是互相无关的，必然是相互有联系的，因此应该认为是物体在变形前是连续的，变形后仍然是连续的，连续物体应变之间关系的数学表达式即为“应变协调方程”。

2023/8/3周书敬第二章应变分析76第五节应变协2023/9/4周书敬第二章应变分析77（2－8）

方程组(2-8)表示的几何方程表明，六个应变分量是通过三个位移分量表示的，这六个应变分量不是互不相关的，它们之间必然存在着一定的联系。这一事实很重要，因为如果我们知道了位移分量，则容易通过(2-8)式获得应变分量；但是反过来，如果纯粹从数学角度任意给出一组“应变分量”,2023/8/3周书敬第二章应变分析77（2－8）2023/9/4周书敬第二章应变分析78则几何方程给出了包含六个方程而只有三个未知函数的偏微分方程组，由于方程的个数超过了未知函数的个数，方程组可能是矛盾的。要使这方程组不矛盾，则六个应变分量必须满足一定的条件。下面的任务就是建立这个条件。为此，我们要设法从方程组（2-8）中消去所有的位移分量。设物体中的某一点的坐标是（x，y，z），其位移是u、v、w，应变为，，若已知u、v、w，则应变便可用位移表示；如果在表达式中消去位移u、v、w，则可得到应变之间的关系。

2023/8/3周书敬第二章应变分析78则几何方程给出了2023/9/4周书敬第二章应变分析79处理方式：现对正应变分别对y、x

取两次偏微分，则有：

将以上两式相加，可得：

这里,我们利用了位移分量具有三阶的连续偏导数的性质。因为，

所以有：

2023/8/3周书敬第二章应变分析79处理方式2023/9/4周书敬第二章应变分析80同理可得另外两个类似的方程，故有式（2－49）（书：2－39）

（2－49）（书：2－39）

这是一组相容方程。

2023/8/3周书敬第二章应变分析80同理可得另外两个2023/9/4周书敬第二章应变分析81若取剪应变的表达式：

将上式的分别对求一阶偏导数，可得：

2023/8/3周书敬第二章应变分析81若取剪应变的表达2023/9/4周书敬第二章应变分析82将上式中的第一式与第三式相加，然后减去第二式，则可得：再对求导得出：

同理可得另外两式，即有式（2－50）（书：2－40）：

（2－50）（书：2－40）

这是又一组相容方程。

2023/8/3周书敬第二章应变分析82将上式中的第一式2023/9/4周书敬第二章应变分析83综合以上（2-49、50）[书2-39、40]两式，有：该式称为“变形协调方程式”或“变形的协调方程”，又称为圣维南（Saint-Venant）方程。是圣维南首次导出的。

(2-51)2023/8/3周书敬第二章应变分析83综合以上（2-42023/9/4周书敬第二章应变分析84其实，通过上述相似的变化，可以导出无穷多组相容方程，但是可以证明，如果满足了上式[(2－49)和(2－50)两组相容方程]，就可以保证位移的连续性。

上式表示要使以位移分量为未知函数的六个几何方程不相矛盾，则六个应变分量必须满足应变协调方程。

方程意义的几何解释：如将物体分割成无数个微分平行六面体，并使每一个微元体发生变形。这时如果表示微元体变形的六个应变分量不满足一定的关系，则在物体变形后，微元体之间就会出现“撕裂”或“套叠”等现象，从而破坏了变形后物体的整体性和连续性。为使变形后的微元体能重新拼2023/8/3周书敬第二章应变分析84其实，通2023/9/4周书敬第二章应变分析85合成连续体，则应变分量就要满足一定的关系，这个关系就是应变协调方程。因此说，应变分量满足应变协调方程，是保证物体连续的一个必要条件。

需要说明的几点：

1、可以证明：如果物体是单联通的，则应变分量满足应变协调方程还是物体连续的充分条件。从数学的观点来看，也就是说，如果应变分量满足应变协调方程，则对于单联通物体，就一定能通过几何方程的积分求得单值连续的位移分量。

2、如果能正确地求出物体各点的位移函数u，v，w，并根据几何方程求出各应变分量，则应变协调方程自然满足。2023/8/3周书敬第二章应变分析85合成连续体，则应2023/9/4周书敬第二章应变分析86

3、从物理意义来看，如果位移函数是连续的，变形自然也就是可以协调。

4、计算时，采用位移法求解，应变协调方程可以自然满足；而采用应力法求解，则需要同时考虑应变协调方程。

5、对于多联通物体，我们总可以作适当的截面使它变成单联通物体，如此则上述的结论完全适用。具体的说，如果应变分量满足应变协调方程，则在此被割开后的区域里，一定能求得单值连续的函数u，v，w。但是对求得的u，v，w，他们在截面两侧趋向于截面上某一点的值一般是不相同的，为了使考察的多联通物体在变形后仍2023/8/3周书敬第二章应变分析862023/9/4周书敬第二章应变分析87保持为连续体，必须加上下列的补充条件：式中：分别为与截面同一点无限临近的两侧点的位移。因此，对于多联通物体，应变分量满足应变协调方程，只是物体连续的必要条件，只有加上补充条件，条件才是充分的。

6、对于平面应变问题，有：

则相容方程只有（2－49）中的第一式。

2023/8/3周书敬第二章应变分析87保持为连续体，必2023/9/4周书敬第二章应变分析88柱坐标中的相容方程：用相同的方法可以导出柱坐标中的变形协调条件为式（2－52）（书：2－41），即：

已知柱坐标系中物体内任意一点的六个应变分量所满足的几何方程的形式为：（2－9）2023/8/3周书敬第二章应变分析88柱坐标中的2023/9/4周书敬第二章应变分析89（2－52）（书：2－41）

2023/8/3周书敬第二章应变分析89（2－52）2023/9/4周书敬第二章应变分析90极坐标中相容方程（平面应变问题）我们知道：对平面问题，柱坐标变为极坐标()，几何方程为（2-10）：由于，变形协调条件只剩下（2－52）中的第三式，即：

（2－10）

2023/8/3周书敬第二章应变分析90极坐标中相容方程2023/9/4周书敬第二章应变分析91（2－53）（书：2－42）

轴对称问题的相容方程：

对于轴对称平面应变问题，应变分量与无关，变形协调条件简化为式（2－54），即：（2－54）（书：2－43）

式（2－54）的左边项可由式（2－53）的第二项获得：即：

2023/8/3周书敬第二章应变分析91（2－53）轴对2023/9/4周书敬第二章应变分析92

应变协调方程的物理意义：如果将变形体分解为许多微元体，每个微元体的变形都用六个应变分量描述。若应变分量不满足应变协调方程，则这些微元体将不能构成一个连续体，因为这时可能会出现裂纹或发生重叠。满足应变协调方程便能保证变形前后物体的连续性，因此，连续介质的应变状态是否可能，需要利用应变协调方程来检验。

球坐标系下的相容方程：几何方程和应变协调方程见相关书籍。

2023/8/3周书敬第二章应变分析92应变协调2023/9/4周书敬第二章应变分析93举例：例1（书中P68）已知下列的应变分量是物体变形时产生的，试求系数之间应满足的关系式。平面应变问题解：该应变状态属于平面应变状态，这些应变分量应满足变形协调条件，即书中（2-39）式的第一式。本题的目的是应变协调方程的应用。2023/8/3周书敬第二章应变分析932023/9/4周书敬第二章应变分析94由应变分量可得：将以上各式代入应变协调条件可得：在物体内任一点上，即x、y为任意值时，上式皆应成立，因此得上式即为系数应满足的条件，而系数可为任意常数。2023/8/3周书敬第二章应变分析94由应变分量可得：2023/9/4周书敬第二章应变分析95例2（书P69）：在平面轴对称情况下，轴向应变为常数，试确定其余两个应变分量和的表达式（材料是不可压缩的）。该问题是轴对称平面应变问题。解释：

轴对称平面应变问题的相容方程（书中式2-43）为：

积分该式可得出上式。

令：2023/8/3周书敬