第三章 应力分析 本章引进了应力张量、导出了平衡方程,并且讨论了作为平衡方程通解的应力函数。 ? ?§1 应力张量 ? 1.1? 外力   弹性体所受的外力可以分为体力和面力两种。作用于弹性体上的重力、电磁力等超距力称为体力。单位体积上的体力记作 ,也可以按极限定义为 ???????????? ????????????????????????(1.1) 其中点 总在体积为 的微元之中, 是该微元上体力的合力。 弹性体与其它物体接触的面上,受有外界给它的力,称为面力,例如流体的压力、固体间的压力和摩檫力等。 ? 1.2? 内力??? ??? 在外力的作用下,弹性体内部的分子的初始状态发生变化,产生了分子之间的附加力,这种力称为内力。分子之间的内力作用距离很小,这种性质称为“短程性”。为显示内力,在弹性体内部过某点P作一小面元 , 面元两侧分别记作A和B(图3.1a,b)。在图3.1a中,向量 表示B部分通过面元 对A部分的作用力,在图3.1b中,向量 则表示A对B的作用力。内力仅通过面来作用是由于它的“短程性”所致。按Newton第三定律, 和 的大小相等、方向相反、作用在不同的部分上。按Cauchy的说法,将称 或 为应力向量。当 收缩至P点时,也可以用形如(1.1)式的极限来定义应力向量,我们仍记作 。显然, 不仅与P的位置有关也与面元 的方向有关。 ? (a)?????? ????????(b) 图3.1 1.3 ?六面体上的应力 为显示应力与方向有关,在弹性体内某点P的邻域内作一小六面体元,它的六个表面分别与坐标面平行,其中三个表面的外法向与坐标法向 ( )分别相同,其余三个表面的外法向则分别与坐标方向相反(图3.2)。 六面体外部关于外法向为 和 面上的应力向量,分别记为 和 ( =1,2,3)。将 在标架 上进行分解(图3.2a),得 ? ? 图3.2   ??? ???????????????????? ????(1.2) 式(1.2)中含有9个分量,可以排成一个矩阵 , ?= ??????????????????????(1.3) 其中 、 、 称为正应力, 、 、 、 、 、 称为剪应力。 ? ??? 引入记号 ?? ??????????????????????????????(1.4) 下面的1.5段中将证明 为张量,它是弹性力学中的一个重要的物理量,称为应力张量,??? 所对应的矩阵如(1.3)所示。式(1.2)的指标形式为 ???????????? ??(i=1,2,3)?????????????????? (1.5) 这里 的第一个脚标 与 所在的面的外法向 相对应,第二个脚标表示 在 方向上投影。对 则在标架 中进行分解(图3.2b), ???????????? ???(i=1,2,3)?????????????? (1.6) (1.5)和(1.6)表明,对 其投影的正方向为 , 其投影的正方向为 ,这种规定虽属人为,却与通常的拉伸为正、受压为负的习惯一致,对今后的应用将带来方便。 ? 1.4? 斜面上的应力? ?? ?在弹性体内某点P附近作一微四面体元PABC,其中PBC、PAC、PAB三个表面分别平行于相应的坐标面。表面ABC的外法向为 。表面ABC上所受的平均应力为 (图3.3)。四面体PABC上所有外力的合力为零,故有 ???????? ?????????????(1.7) 其中 、 、 ; 、 、 ; 、 、 分别为面PBC、PAC、PAB上的平均应力, 、 、 是 的分量, 、 、 为四面体PABC内的平均体力, 、 、 、 分别为面PBC、PAC、PAB、ABC的面积, 为四面体PABC的体积。  图3.3 ??? ? ??? 按解析几何可知 ? ????????? ?(1.8) 其中( ,h为点P至斜面ABC的高。将(1.8)代入(1.7),约去 ,并令 ,可以得到斜面应力公式 ???????????????????????????????????????? ?????? (1.9)??? 在(1.9)中的 和 ( )诸量都是先在各相应面上取平均,当 时,它们都与所取四面体微元无关。坐标形式的斜面应力公式(1.9)的指标形式和整体形式为 ????(i=1,2,3)????????????????? (1.10) ??? ?????????????????????????????????? (1.11)???? 式(1.9)-(1.11)表示出在点P截面上应力 与点P和法向 的关系,它表明应力张量 足以表征一点的应力状态。   ??? 斜面应力公式另一用途是表示弹性力学边值问题的应力边条件(见第五章)。 ?? ? 1.5? 应力张量 除标架 外,考虑一个新标架 ,新旧标架的关系为 , ( =1,2,3)??????????? ??????????(1.12) 过点P作法向为 的截面,其上的应力向量为 ,将它投影到截面的法向 和切向 ,分别记为 和 ,有 ????????????????????? (1.13) 类似的有 的表示式,可以将这些公式统一写成 ??????? ,(i,j=1,2,3)???????????? (1.14) 量 是法向为 截面上的应力向量在 上的投影, 的这个力学解释与量 的意义一致,仅坐标系不同。也就是说,具有明确力学意义的量 在不同坐标系下服从关系式(1.14),而它恰是关于变换(1.12)系数的二次齐次式,因此由(1.4)所定义的 为张量。 §2 平衡方程 ? 2.1 力的平衡 设坐标为 的点P位于六面体微元的中心(图3.4),微元的边长设为 。考虑 方向上外力的平衡,得  (2.1)  图3.4 上式中的前两项分别为外法向是 和 前后两个截面上的应力,其中略去了Taylor展开的高次项,式中第三、四项,第五、六项分别为右左和上下各个截面上的应力,最后一项为六面体所受的体力。在(2.1)中,先消去相同的项,再约去因子 ,最后令 趋于零,即六面体收缩至点P,得到下面(2.2)式中的第一式,  (2.2) 其中第二、三式可按导出第一式同样的过程得到, 和 仅与点P有关,为 的函数。方程(2.2)称为平衡方程,它是弹性理论的第二组方程。第一组方程是第一章中已导出的几何方程。方程(2.2)的指标形式和整体形式为  ( ) (2.3)  (2.4) §3 主应力 偏应力 ? 3.1? 主应力 设 和 分别为应力张量 在新旧坐标系下所对应的矩阵,按(1.14),有  (3.1) 其中 为坐标转换矩阵,上标“ ”表示转置。(3.1)表明对称矩阵 和 服从矩阵的合同变换。按矩阵理论,在适当的坐标变换下,应力张量 所对应的矩阵为对角形,即有 ??  (3.2) 与(3.2)相应的坐标系称为应力张量 的主坐标系。显然,对于不同的点,应力张量 不同,那么主坐标系也不同。主坐标系的三个方向称为主方向,(3.2)中的 称为主应力。我们有如下的结论 1° 主方向互相垂直, 2° 主方向间的剪应力为零。 应力张量 的三个不变量为 ???? ??????? ??????? (3.3) ???????  如果用主应力来表示,三个不变量为  (3.4) §4 应力函数 ?无体力的平衡方程为 ? (4.1) 或  (4.2) 众多学者研究过(4.1)或(4.2)的解。1863年,Airy首先给出(4.1)的一种解为  (4.3) 将(4.3)代入(4.1),不难验证它满足。(4.3)中的 称为Airy应力函数。1868年,Maxwell得到了较Airy广泛的解, , ,  , , ??????? (4.4) 其中 均为 的任意函数(当然假定连续可微到所需要的阶数),它们被称为Maxwell应力函数。1892年,Morera得到了平衡方程(4.1)的另一组解 , , , ,  (4.5) 其中 被称为Morera应力函数。1892年,Beltrami把Maxwell解和Morera解叠加,得到了平衡方程(4.2)相当广泛的解  (4.6) 其中 为对称张量, 称为Beltrami应力函数, 所对应的矩阵 与Maxwell应力函数 和Morera应力函数 的关系为, ??????????????  Schaefer(1953)在Beltrami解的基础上又增加了新的项,其解为  (4.7) 其中 为对称张量, 为调和函数向量。不难验证解(4.7)满足方程(4.2)。(4.7)称为Beltrami-Schaefer解。 ??? 长期以来,人们误认为Beltrami解(4.5)是平衡方程(4.2)的一般解。但是,Carlson(1966)指出仅仅对自平衡场(即任意封闭曲面上合力合力矩为零的场)解(4.5)才完备。由于存在非平衡场,由此(4.5)一般说来不完备。关于解(4.7),Gurtin(1972)证明了如下重要命题, 定理 (Beltrami-Scharfer解的完备性)设 是区域 上的应力场,它满足平衡方程  (4.8) 则在 上存在对称张量场 和调和向量场 ,使 可表成如下形式  (4.9) 证明 设 为 的Newton位势,即  (4.10) 由于 是对称张量, 也是对称张量,且按位势理论有  (4.11) 由于 的对称性,第一章的恒等式(4.15)可写成  (4.12) 将(4.12)代入(4.11),并令 , (4.13) 即得欲证的(4.9)式。对(4.13b)取Laplace算子,考虑到Newton位势(4.10)和平衡方程(4.8),可知 为调和向量。证毕。