弹性力学-平面问题的极坐标解答课件.pptx

1、第四章第四章 平面问题的极坐标解答平面问题的极坐标解答4-1 4-1 极坐标中的平衡微分方程极坐标中的平衡微分方程4-9 4-9 圆孔的孔边应力集中圆孔的孔边应力集中4-4 4-4 应力分量的坐标变换式应力分量的坐标变换式4-3 4-3 极坐标中的应力函数与相容方程极坐标中的应力函数与相容方程4-2 4-2 极坐标中的几何方程及物理方程极坐标中的几何方程及物理方程4-5 4-5 轴对称应力和相应的位移轴对称应力和相应的位移4-6 4-6 圆环或圆筒受均布压力。压力隧洞圆环或圆筒受均布压力。压力隧洞4-7 4-7 曲梁的纯弯曲曲梁的纯弯曲4-8 4-8 圆盘在匀速转动中的应力及位移圆盘在匀速转动

2、中的应力及位移4-10 4-10 楔形体在楔顶或楔面受力楔形体在楔顶或楔面受力4-11 4-11 半平面体在边界上受法向集中力半平面体在边界上受法向集中力习题课习题课4-1 4-1 极坐标中的平衡微分方程极坐标中的平衡微分方程 在处理弹性力学问题时，选择什么形式的坐标系统，虽不会影响对问题本质的描绘，但却直接关系到解决问题的难易程度。如坐标选得合适，可使问题大为简化。例如对于圆形、楔形、扇形等物体，采用极坐标求解比用直角坐标方便的多。图41 考虑平面上的一个微分体 ，沿 方向的正应力称为径向正应力，用 表示，沿 方向的正应力称为切向正应力，用 表示，剪应力用 表示，各应力分量的正负号的规定和直

3、角坐标中一样。径向及环向的体力分量分别用 及 表示。如图4-1。PACBrrrrKKrrrrdrrrrrdddrrdrrdrKrKyxoPABC 考虑图示单元体的平衡，有三个平衡方程：0, 0, 0MFFr由 ，可以得出剪应力互等关系： 0Mrr 0rF由 ，有：0)(22)()(drrdKdrdrdddrddrdrdddrrdrrrrrrrrr 0F由 ，有：022)()()(drrdKddrddrdrdddrrdrrdrdrdrrrrrr因为 很微小，所以取 ， ，并用 代替 ，整理以上两式，得：d22sindd12cosdrr02101KrrrKrrrrrrrrr这就是极坐标的平衡微分方

4、程。 两个平衡微分方程中包含三个未知函数 、 和 ，所以问题是静不定的。因此必须考虑变形条件和物理关系。rrr 上述方程和直角坐标系下的平衡方程有所不同，直角坐标系中，应力分量仅以偏导数的形式出现，在极坐标系中，由于微元体垂直于半径的两面面积不等，而且半径愈小差值愈大，这些反映在方程里带下划线的项中。一、几何方程一、几何方程位移与形变间的微分关系4-2 4-2 极坐标中的几何方程及物理方程极坐标中的几何方程及物理方程在极坐标中规定:rrruu -径向正应变-环向正应变-剪应变(径向与环向两线段之间的直角的改变)-径向位移-环向位移用叠加法讨论极坐标中的形变与位移间的微分关系。图4-2drdrr

5、uo（1）假定只有径向位移，而无环向位移。如图4-2所示。径向线段 的正应变为：PArudrudrruurrrrr)(环向线段 的正应变为：PBrurdrddurrr)(径向线段 的转角为：PA0环向线段 的转角为：PBrrrrurrduduu1)(可见剪应变为：rrur1drPP BB A Adruo（2）假定只有环向位移，而无径向位移。如图4-3所示。图4-3径向线段 的正应变为：PA0r环向线段 的正应变为：PBurrduduu1)(径向线段 的转角为：PArudrudrruu)(环向线段 的转角为：PBru可见剪应变为：rurur 如果同时存在径向和环向位移，则由叠加法得：ruruur

6、urrururrrrr11这就是极坐标中的几何方程。二、物理方程二、物理方程（1）平面应力情况：rrrrrrEGEE)1 (21)(1)(1（2）平面应变情况：rrrrrEEE)1 (2)1(1)1(122 将上式中的 换为 ， 换为 。E21E14-3 4-3 极坐标中的应力函数与相容方程极坐标中的应力函数与相容方程 为了得到极坐标中用应力函数表示的应力和相容方程，利用极坐标和直角坐标的关系：sin,cosarctan,222ryrxxyyxr得到：rrxyrryxryyrrxxrcos,sin,sin,cos22rryyrryrrxxrrxcossinsincos2222222222222

7、22222222222coscossin2coscossin2sinsincossin2sincossin2cosrrrrrrryrrrrrrrx222222222222cossinsincoscossinsincoscossinrrrrrrryx 在=0时，极坐标的各分量和直角坐标各分量相同。将上面各式代入应力分量的表达式（常体力）：yxxyxyyx22222（a）（b）（c）得到：)1()()()()(11)()(0202202202220220rryxrxrrryyxryxr可以证明，当体力为零时，这些应力分量确能满足平衡微分方程。由（a）+（b），得：22222222211rrrryx

8、于是由直角坐标的相容方程：0)(22222yx得到极坐标中的相容方程：0)11(222222rrrr 用极坐标求解平面问题时（体力不计），就只须从相容方程求解应力函数 ，然后求出应力分量，再考察应力分量是否满足边界条件，多连体还要满足位移单值条件。),(r4-4 4-4 应力分量的坐标变换式应力分量的坐标变换式 在一定的应力状态下，如果已知极坐标中的应力分量，就可以利用简单的关系式求得直角坐标中的应力分量。反之，如果已知直角坐标中的应力分量，也可以利用简单的关系式求得极坐标中的应力分量。 设已知极坐标中的应力分量 、 、 。试求直角坐标中的应力分量 、 、 。rrxyxyrryyxrrrrxy

9、xcaboyxAB图4-4 如图4-4，在弹性体中取微小三角板 ，各边上的应力如图所示。三角板的厚度取为一个单位。令 边的长度为 ，则 边及 边的长度分别为 及 。 Abcdsabacsindscosds 根据三角板 的平衡条件 ，可得平衡方程：A 0 xF0cossinsincossincos22dsdsdsdsdsrrrx用 代替 ，得：rrcossin2sincos22rrx同理，由平衡条件 ，可得： 0yF)sin(coscossin)(22rrxy另取微小三角板 ，如图4-4，根据平衡条件 ，得到： 0yFBcossin2cossin22rry 综合以上结果，得出应力分量由极坐标向直

10、角坐标的变换式为：)sin(coscossin)(cossin2cossincossin2sincos222222rrxyrryrrx利用简单的三角公式，上式可改写为：2cos2sin22sin2cos222sin2cos22rrxyrrryrrrx4-5 4-5 轴对称应力和相应的位移轴对称应力和相应的位移 如果应力分量仅是半径的函数，如受内外压的圆环，称为轴对称问题。 采用逆解法，假定应力函数 仅是径向坐标 的函数:r)(r相容方程简化为:0dd1dd222rrr这是一个四阶常微分方程，它的通解为:DCrrBrrA22lnln 这时，应力的表达式为:02)ln23(2)ln21 (22rr

11、rCrBrACrBrA 正应力分量仅是 的函数，与 无关，并且剪应力为零，应力分量对称于通过z轴的任一平面，称为轴对称应力。r 将上述应力的表达式代入应力应变关系式中，可以得到应变的表达式,再代入位移与应变积分后的几何方程,得到轴对称应力状态下的位移分量:cossin4sincos)1 (2)31 () 1(ln)1 (2)1 (1KIHrEBruKICrBrrBrrAEur 对于平面应变问题，须将上面公式 换为 ， 换为 。E21E14-6 4-6 圆环或圆筒受均布压力。压力隧洞圆环或圆筒受均布压力。压力隧洞 如图4-5，圆环的内半径为a，外半径为b，受内压力qa，外压力qb。为轴对称问题。

12、根据上节有解为:02)ln23(2)ln21 (22rrrCrBrACrBrA图4-5边界条件为:bbrraarrbrrarrqq)(,)(0)(, 0)(一、圆环或圆筒受均布压力一、圆环或圆筒受均布压力 在这里只有两个方程，而有三个待定常数，需要从多连体的位移单值条件补充一个方程。在环向位移表达式：cossin4KIHrEBru中，第一项是多值的，在同一r处， =1和=1+2时，环向位移成为多值，这是不可能的，因此，从位移单值条件必须有B=0。 这样从上面两个方程中可解出A和C，代入应力分量表达式，得到拉密解答：baqCbAqCaA2222于是:baqCbBbAqCaBaA2)ln1 (2)

13、ln1 (22由边界条件得到:babarqbaraqabrbqbaraqabrb222222222222222211111111 下面分别讨论内压力和外压力单独作用的情况。（1）只作用均匀内压时（ ），例如液压缸，上面解答化为：0bqaarqabrbqabrb11,1122222222r图4-6应力分布大致如图4-6所示。当 时，得到具有圆孔的无限大薄板，或具有圆形孔道的无限大弹性体，这时上面的解答成为：baarqraqra2222,（2）只有外压时（ ），例如液压柱塞，上面解答化为：0aqbbrqbaraqbara2222222211,11应力分布大致如图4-7所示。r图4-7二、压力隧洞二

14、、压力隧洞qo,E,Errr图4-8 如图4-8所示，受均匀内压力 作用的圆筒埋在无限大弹性体中，圆筒和无限大弹性体的材料不同。试分别讨论两者的应力和位移情况。q 两者都属于轴对称应力问题，采用半逆解法。设圆筒的应力表达式为：CrACrAr2,222设无限大弹性体的应力表达式为：CrACrAr2,222由应力边界条件求待定常数 、 、 、 。ACAC（1）在圆筒的内表面：qarr)(由此得：qCaA22（2）在无限大弹性体内距离圆筒很远处几乎没有应力。0)( , 0)(rrr由此得：02C（3）在圆筒和无限大弹性体的接触面上，应当有：brrbrr)()(（1）（2）由此得：CbACbA2222

15、三个方程不足以确定四个常数，下面来考虑位移。 由于圆筒和无限大弹性体都是多连体，并属于平面应变问题，可以写出两者的径向位移的表达式。圆筒：sincos)11 (2)11 (12KICrrAEur无限大弹性体：sincos)11 (2)11 (12KIrCrAEur将以上两式简化后得：sincos)21 (21KIrACrEursincos)21 (21KIrArCEur（3）在接触面上，两者应具有相同的位移，即：brrbrruu)()(因此有：sincos)21 (21sincos)21 (21KIbAbCEKIbACbE因为这一方程在接触面上的任意一点都应当成立，也就是在 取任何数值时都应当

16、成立，所以方程两边的自由项必须相等。于是有：)21 (21)21 (21bAbCEbACbE简化后，得：0)21 (222bAbACn其中：)1 ()1 (EEn（4）联立方程（1）、（2）、（3）、（4）求出 、 、 、 ，代入应力分量的表达式，得： ACAC)1 ()21 (1 )1 (2)1 ()21 (1 )1 ()21 (1 )1 ()21 (1 )1 ()21 (1 222222222222nabnrbnqnabnnrbnqnabnnrbnqrr 当 时，应力分布大致如图4-8所示。1n4-7 4-7 曲梁的纯弯曲曲梁的纯弯曲r 内半径为a，外半径为b的狭矩形截面的圆轴曲梁，在两端

17、受大小相等、方向相反的弯矩，为轴对称问题。有：02)ln23(2)ln21 (22rrrCrBrACrBrA边界剪应力都为零：0)( , 0)(0)( , 0)(0rrbrrarr图4-9 在梁的内外两面，正应力要求:0)(,0)(brrarr02)ln21 (02)ln21 (22CbBbACaBaA从而可得：在梁端的边界条件要求:Mrrrbabad0d0dddd22arrbrrbarbababaabrrrrr Marabrrrrrrrrrrrrbarbrrbabarbababababa22222)()(dddddddddd则：将 的表达式：DCrrBrrA22lnln代入，并由边界条件得：

18、MabCaabbBabA)()lnln(ln2222 在这里有三个方程和三个待定常数，解出A、B和C，代入应力分量表达式，得到郭洛文解答:0lnlnln14lnlnln4222222222222rrrabrbarrbababNaMabrbarrbabNaM其中：222222ln41abababN 4-8 4-8 圆盘在匀速转动中的应力及位移圆盘在匀速转动中的应力及位移一、等厚度圆盘在匀速转动中的应力及位移一、等厚度圆盘在匀速转动中的应力及位移 设有等厚度圆盘，绕其回转轴以匀角速度 旋转。圆盘可以认为是在下面的体力作用下处于平衡状态：0,2KrKr 由于这里是轴对称的物体受轴对称的体力，所以应力

19、分布也是轴对称的。 即：应力分量 及 都只是 的函数，而 rr0rr 。所以有平衡微分方程：02102Krrrrrdrdrrrr令：22,rdrdrr(1) 在这里，由于圆盘只受回转轴的约束，而这种约束是轴对称的，所以它的弹性位移也是轴对称的。即：径向位移 ，而环向位移 。于是几何方程简化为：)(ruurr0u0,rrrrruudrd消去 ，得到相容方程：ru)(rdrdr32222)3(rdrdrdrdr解方程得到：将物理方程代入，再联立式(1)，得到由应力函数 表示的相容方程：rBArr28332联立式（1），得：2222222831283rBArrBArr(2)其中 和 是任意常数。AB

20、 盘边的边界条件：0)(arr其中 是圆盘的半径。代入式（2），得：a2243aA取 ，代入式（2）得应力分量的表达式为：0B)3311 (83)1 (8322222222araarar最大应力在圆盘的中心：2200max83)()()(arrr径向位移：)1 ()3(8)1 ()(3332ararEaErrurr在圆盘的中心（ ）， 。最大弹性位移发生在圆盘的边缘（ ）：0r0ruar Eaur4)1 ()(32max二、变厚度圆盘在匀速转动中的应力及位移二、变厚度圆盘在匀速转动中的应力及位移 假定圆盘的厚度为 ，而应力不沿厚度变化，则等厚度圆盘的微分方程可以近似地应用于每单位厚度的圆盘。于

21、是可得全厚度内的平衡微分方程为：)(rtt 0)(2rtrtttdrdrr令：22,trdrdtrtr可得：32222)3()1 ()1 (trdrdttrdrdrdrdttrdrdr 取厚度的变化规律为： Crt其中 是常数， 为任意正数。则上式成为：C32222)3()1 ()1 (Crdrdrdrdr解方程，得：32)3(83rCBrArnm其中 和 是任意常数，而：AB)1 ()2(22nm由此可得出应力分量：2211222211)3(8311)3(83rnrCBmrCArdrdtrrCBrCAtrnmnmrmaCA32)3(83由边界条件 ，求得： 0)(arr为了应力在圆盘的中心（

22、 ）处不成为无限大，取 。0r0B从而得应力分量为：)(331)()3(83)()()3(8321222122ararmaararammr)(1 ()()(3()3(8)(3232ararmEaErrumrr且，有：4-9 4-9 圆孔的孔边应力集中圆孔的孔边应力集中 板中开有小孔，孔边的应力远大于无孔时的应力，也远大于距孔稍远处的应力，称为孔边应力集中。 应力集中的程度与孔的形状有关。一般说来，圆孔孔边的集中程度最低。这里简略讨论圆孔孔边应力集中问题，较为复杂的孔边应力集中问题一般用复变函数方法，在第五章中进行讨论。rrAb图4-10一、一、 矩形板左右两边受集度为矩形板左右两边受集度为q

23、q的均布拉力的均布拉力 设有矩形薄板，在离开边界较远处有半径为 的小圆孔，在左右两边受均布拉力，其集度为 ，如图4-10。aq2sin2)(,2cos22)(qqqbrrbrr 以远大于 的某一长度 为半径，以小孔中心为圆心作圆，根据直角坐标与极坐标的变换公式，得到大圆的边界条件：ba上述面力可以分解成两部分，其中第一部分是：0)(,2)(brrbrrq第二部分是：2sin2)(,2cos2)(qqbrrbrr求面力（a）所引起的应力。令： 。得： 2qqb0,112,11222222222rrbaraqbaraq(a)(b)由于 ，所以可近似地取 ，从而得到解答：ab 0ba0),1 (2)

24、,1 (22222rrraqraq求面力（b）所引起的应力。采用半逆解法：假设 为 的某一函数乘以 ，而 为 的另一函数乘以 。即：rr2cosrr2sin2sin)(,2cos)(21rrrr又应力函数和应力分量之间的关系为：)1(,11222rrrrrrr因此可以假设：2cos)(rf代入相容方程，得：0)(9)(9)(2)(2cos32223344drrdfrdrrfdrdrrfdrdrrfd删去 ，求解常微分方程，得：2cos224)(rDCBrArrf从而得应力函数：2cos)(224rDCBrAr从而得应力分量：2sin)6226(2cos)6212(2cos)642(422424

25、2rDrCBArrDBArrDrCBrr 对上式应用边界条件（b），并由边界条件：0)(, 0)(arrarr得到方程：264242qbDbCB0622606422622642242422aDaCBAaaDaCBqbDbCBAb求解 、 、 、 ，令 ，得：ABCD0ba4,4, 042qaDqaCqBA将各已知量代入应力分量表达式，即得齐尔西的解答：2sin)31)(1 (22cos)31 (2)1 (22cos)31)(1 (2)1 (222224422222222raraqraqraqraraqraqrrr二、矩形板四边受二、矩形板四边受q的均布拉力的均布拉力 如果矩形薄板在左右两边受有

26、均布拉力 ，并在上下两边受有均布拉力 ，如图4-11，也可由前面解答得出应力分量。首先命该解答中的 等于 ，然后命该解答中的 等于 ，将 用 代替，最后将两个结果相叠加。得到：1q2qq1qq2q0902sin)31)(1 (22cos)31 (2)1 (22cos)31)(1 (2)1 (2222221442122212222212221raraqqraqqraqqraraqqraqqrr图4-111q1q2q2q4-10 4-10 楔形体在楔顶或楔面受力楔形体在楔顶或楔面受力rrrP图4-12 楔形体的中心角为 ，下端为无限长。1. 顶部受集中力P 设楔形体在楔顶受有集中力，与楔形体的中心

27、线成角 。取单位宽度的部分来考虑，并令单位宽度上所受的力为 。 楔形体内一点的应力分量决定于、P、r、，因此，应力分量的表达式中只包含这几个量。其中、是无量纲的量，因此根据应力分量的量纲，应力分量的表达式应取PN/r的形式，其中N是、组成的无量纲的量。由应力函数的表达式可以看出应力函数中r的幂次应当比各应力分量的幂次高出两次，因此可设： P)(rf代入相容方程后得：0)(d)(d2d)(d122443fffr求解这一微分方程，得：)sincos(sincos)sincos(sincos)(DCrBrArDCBAfByAxBrArsincos不影响应力，取：其中)sincos(DCr于是得：0)

28、1(0)sincos(21122222rrrCDrrrrrrr楔形体左右两面的边界条件：0)(, 0)(2/2/ara 上述应力分量满足该边界条件。集中力P按圣维南原理处理，取出任一圆柱面ab，则该截面上的应力和P合成平衡力系：0sinsind:00coscosd:02/2/2/2/PrFPrFryrx将 的表达式代入，可求出C、D，最后得到密切尔解答：r00)sincossinsincoscos(2rrrrP2.顶部受有力偶M作用rrr图4-13 设楔形体在楔顶受有力偶，而每单位宽度内的力偶矩为M ，如图4-13。 根据和前面相似的分析，应力分量应为MN/r2的形式，而应力函数应与r无关。)

29、(代入相容方程后，得：0dd4dd122444r求解这一微分方程，得：DCBA2sin2cos 力偶可看成反对称力，正应力和应力函数应当是 的奇函数，从而A=D=0,于是：22222222cos2)1(02sin411rCBrrrrBrrrrrr楔形体左右两面边界条件：0)( ,0)(2/2/ara上述应力分量自动满足第一式，根据第二式，可得：cos2BC于是：22)cos2(cos202sin4rBrBrrr 集中力偶M按圣维南原理处理，取出任一圆柱面ab，则该截面上的应力M成平衡力系：cossin20d:022/2/MBMrMro最后得到英格立斯的解答：22)cos(sin)cos2(co

30、s0)cos(sin2sin2rMrMrrr3.一面受均布压力qr图4-14 设楔形体在一面受有均布压力 ,如图4-14。q 应力分量应为qN的形式，而应力函数应为qNr2的形式：)(2fr代入相容方程后，得：0d)(d4d)(d122442ffr求解这一微分方程，得:)2sin2cos(2DCBArCBADCBADCBArrr2cos22sin2222sin22cos2222sin22cos2边界条件为:0)(, 0)(0)(,)(00rrq求解常数，得应力分量的李维解答：qqqqqrrr)(tg22sintg)2cos1 ()(tg2)2sin2()2cos1 (tg)(tg2)2sin2

31、()2cos1 (tg4-11 4-11 半平面体在边界上受法向集中力半平面体在边界上受法向集中力xyPro00cos2rrrrP利用坐标变换可得到直角坐标中的应力分量式（2）：rPrPrPxyyx223cossin2cossin2cos2（1）（2） 命楔形体的中心角等于一个平角，这楔形体的两个侧边就连成一个直边，而楔形体就成为一个半平面体，如图4-15。一、应力分量一、应力分量0P 当平面体在边界上受有垂直于边界的力 时，在密切尔解答中令 、 。于是得式（1）：图4-15或将其中的极坐标改为直角坐标而得：222222222223)(2)(2)(2yxyxPyxxyPyxxPxyyx二、位移

32、分量二、位移分量 假设是平面应力情况。将应力分量代入物理方程，得形变分量：0,cos2,cos2rrrEPrEP再将形变分量代入几何方程，得：01cos21cos2ruruurrEPurrurEPrurrr于是可以得出位移分量：cossincos)1 (sin)1 (lnsin2sincossin)1 (lncos2KIHrEPEPrEPuKIEPrEPur其中 、 、 都是任意常数。HIK 由对称条件 ，得：0)(0u0, 0KH代入式（3），得：（3）sinsin)1 (cos)1 (lnsin2cossin)1 (lncos2IEPEPrEPuIEPrEPur 如果半平面体不受铅直方向的

33、约束，则常数 不能确定。如果半平面体受有铅直方向的约束，就可以根据这个约束条件来确定常数 。II（4） 边界上任意一点 向下的铅直位移，即所谓沉陷。由式（4）中的第二式可得 点的沉陷为：MMIEPrEPu)1 (ln2)(2 如果常数 没有确定，则沉陷也不能确定。这时只能求出相对沉陷。I 在边界上取定一个基点 ，它距载荷作用点的水平距离为 。则边界上一点 对于基点 的相对沉陷，等于 点的BsMBM沉陷减去 点的沉陷，如图4-16：B)1 (ln2)1 (ln2IEPsEPIEPrEP简化以后，得：rsEPln2POsyxBMr图4-16 半平面体在边界上受法向分布力作用时的应力和沉陷,可以由半

34、平面体在边界上受法向集中力用叠加法得出。0a0c0b练习练习1 如图1所示，由内外筒组成的组合筒（长度有限，两端自由），装配前内筒的外半径比外筒的内半径大 ，求接触压力 ，并导出环向预应力的表达式。Ap解：解：1.设装配后接合处的公共半径为 ，接触压力 使内筒的外半径减小了 ，而使外筒的内半径增大了 ，按位移协调条件有：0cp21212.将 代入只受内压力作用圆环的位移公式，得：pqbbcacra,000)(120202020101cbcbEpc（1）（2）平面问题的极坐标解答平面问题的极坐标解答习题课习题课图1)(220202020202acacEpc（3）将式（2）、（3）代入式（1），得

35、：)()(2202020202012020202010acacEpccbcbEpc3.若内、外筒为同一种材料，则 ，从上式可解得：2121,EEE)(2)(20203020202020abcaccbEp4.内、外筒的环向应力为：)1 ()(),1()(220202020220202020rbcbpcraacpc内外将 代入只受外压力作用圆环的位移公式，得：pqcbaacrb,000解：解：oyxPrr 练习练习22楔形体顶端受集中力 作用， 与 轴的夹角为 ，如图2所示。取单位厚度考虑，试确定楔形体内的应力分量。Px1.由于描述楔形体几何特征的角度 是无量纲的，故可由量纲分析法得知，应力函数中

36、 只能以一次幂形式出现，即：r)(rf2.由调和方程求出 后，即可求得应力函数为：)(f)sincos(sincosDCBAr 由于 不影响应力分量，故可删去，因此有：ByAxBrArsincos0),sincos(2)sincos(rrCDrDCr（1）（2）（3）（4）图23.楔形体两侧面的边界条件能自然满足：0)(, 0)(ara考虑半径为 的楔形体上部的静力平衡条件：raaroaaraaryaaaarrxdrMPrdrdFPrdrdF0:00sincossin:00cossincos:02由前两式可解出 和 ，从而求出应力分量（密切尔解）：CD2sin2cos,2sin2sinPDPC

37、0,)2sin2(sinsin2)2sin2(coscos2rrrPrP 练习练习33求图3所示问题的截面m-n上的应力 。xmnxPPaM oxymnxyPoxac2解：解：（a）（b）将图（a）所示力系向 点简化，便得图（b）所示与原力系静力等效的力系，其中 。根据圣维南原理，此类代换对远离楔顶之处的应力的影响可不计。将楔顶受集中力作用与受力矩作用下的应力解答叠加，得原问题的应力：oPaM 图32cos22sin2cos2cos, 02cos22sin2sin22sin2sin222rParParPrr由应力分量的坐标变换式可得：)(2cos)()(32cos22sin2)(2sin22c

38、ossin2sincos22232233223222222yxxyyxxyyxyxPayxyxPrrx) 1 (0)11)(11(2222222222rrrrrrrr 练习练习4 4 试将以 表示的相容方程式 展开。011222222rrrr)1()1(1)1(1)1()(1)(1)(222222222222222222222rrrrrrrrrrrrrrrrrrr)2(0)1(1)1(1222222222rrrrr分项求偏导数，最后相加，得：解：解：22422222233224422221)(1,1)(1,)(rrrrrrrrrrrr233222322222233322221)1(1,11)1

39、(11211)1(rrrrrrrrrrrrrrrrrrrrrrrr4442222222332242222242233233224222221)1(112)1(11226)1(rrrrrrrrrrrrrrrrrr)3(01422112444224233224232223344rrrrrrrrrrrrr 练习练习55 等厚度圆环的内外半径分别为a和b，以等角速度 旋转，见图4a），试求其应力和位移。yzxabOxOaba)b)1、本题为位移轴对称问题（平面应力情况），其特征是 只是r的函数，根据基本方程，求解本题的平衡微分方程为式（1）,其中 是材料的比重，ruu,0g是重力加速度，几何方程与物理

40、方程分别为：)2(0,) 1 (02rrrrrrrudrdurgrrdrd解：解：图4把式（3）代入式（1），得：)4(8)1 (011222222rEgrrBArurgrErudrdurdrdurrrr式（4）代入式（3）得应力表达式为：)5(8)31 (8)3(2222122221rgrCCrgrCCr) 3(0),(1),(122rrrrrrdrduruErudrduE式中， 。1,121EBCEAC2、由内外圈的边界条件确定常数，进而求出位移和应力：31)1 ()1(8)3(8)31 (8)3(8)3(8)3(,8)3(0)(, 0)(2222222222222222222222222

41、2221rrbarbagurgrbabagrbarbagbagCbagCrrbrrarr 练习练习6 6 楔形体右侧面受均布荷载q作用，见图5。试求其应力分量。yqxrO) 1 ()(2fr2、由双调和方程确定f（），由平衡微分方程在一般情况下的通解得应力表达式：DCBAfdfddfdr2sin2cos)(, 0)(4)(1224421、楔形体内任一点的应力分量决定于 其中，q的量纲为力长度 ，与应力的量纲相同。因此，各应力分量的表达式只可能取Nq的形式，而N是以 表示的无量纲函数，也就是应力表达式中不可能出现 ，再由平衡微分方程在一般情况下的通解知，应力函数 应是 的某一函数乘以 ，即可设：,rq2,r2r图5解：解：)3(2sin22cos2)2sin2cos(2)2sin2cos(2)2()2sin2cos(2CBADCBADCBADCBArrr3、应力边界条件为：0)(, 0)(, 0)(,)(00arraq由此可求出常数，代回式（3）得应力分量（李维解答）为：)4()(tan22sintan)2cos1()(tan2)2sin2()2cos1(tan)(tan2)2sin2()2cos1(tanqqqqrr