第3章 压力容器安全设计的理论与基础知识3(1)

1、第三章 压力容器安全设计的理论与基础知识 3-4壳体的边缘应力 3-4壳体的边缘应力壳体的边缘应力 薄膜应力是假设壳壁很薄，根本不能随弯矩，没有弯曲应力。实际上压力容器的壳体总有一定的刚度，受压时半径增大，壳壁曲率发生变化，壳体总是存在一些弯曲应力。另外，几何形状不连续，也会产生“不连续应力”或“边缘应力”(变形不同，互相牵制”。壳体不连续应力的影响范围很小，即它只存在于联接处两边附近的很窄的一个区域内，而且它也不直接影响到壳体的破坏强度。但在一些不合理的结构中，不连续应力可以达到很高的数值，而高的局部应力对受反复载荷的容器的疲劳寿命是有很大影响的。 边缘应力的概念 (1)圆筒受内压直径增大时

2、，筒壁金属的环向“纤维”不但被拉长了，而且它的曲率半径由原来的R变到R+R，如图所示。根据力学可知，：有曲率变化就有弯曲应力。所以在内压圆筒壁的纵向截面上，除作用有环向拉应力2外，还存在着弯曲应力2b。但由于这一应力数值相对很小，可以忽略不计。边缘应力的概念 (2)圆筒与封头、圆筒与法兰、不同厚度或不同材料的筒节、裙式支座与直立壳体相联接处的平行圆等。此外，当壳体经线曲率有突变或载荷沿轴向有突变的接界平行圆，亦应视作联接边缘，以上各种情况参见图325。 边缘应力的概念 (3)圆筒形容器受内压之后，由于封头刚性大，不易变形而筒体刚性小，容易变形，连接处二者变形大小不同，即圆筒半径的增长值大于封头

3、半径的增长值，如图326a左侧虚线所示。如果让其自由变形，必因两部分的位移不同而出现边界分离现象，显然。这与实际情况不符。实际上由于边缘联接并非自由，必然发生如图326a右侧虚线所示的边缘弯曲现象，伴随这种弯曲变形。也要产生弯曲应力，因此，联接边缘附近的横截面内，除作用有轴(经)向拉伸应力1外，还存在着轴(经)向弯曲应力1b，这就势必改变了无力矩应力状态，用无力矩理论就无法求解。 边缘应力的概念 (3)分析这种边缘弯曲的应力状态，可以将边缘弯曲现象看作是附加边缘力和弯矩作用的结果，如图326b所示。意思就是在壳体两部分受薄膜力之后出现了边界分离，若再加上边缘力和弯曲使之协调，才能满足边缘联接的

4、连续性。因此联接边缘处的应力就特别大。如果确定这种有力矩的应力状态就可以简单地将薄膜应力与边缘弯曲应力叠加。边缘应力的特点： 今有一内径为Di=1000毫米，壁厚S=10毫米的钢制内压圆筒，其一端为平板封头，且封头厚度远远大于筒体壁厚。内压为P=1MPa，经理论计算和实测其内、外壁轴向应力(薄膜应力与边缘弯曲应力的叠加值)分布情况如图327所示。 边缘应力的特点： 其一，局部性。不同性质的联接边缘产生不同的边缘应力，但它们都有一个明显的衰减波特性。以圆筒壳为例，其沿轴向的衰减经过一个周期之后，即离开边缘距离为2.5 (其中r与s分别为圆筒的半径与壁厚)之处边缘应力已经基本衰减完了。 rs边缘应

5、力的特点： 其二，自限性。从根本上说，发生边缘弯曲的原因是由于薄膜变形不连续。自然，这是指弹性变形。当边缘两侧的弹性变形相互受到约束，则必然产生边缘力和边缘弯矩，从而产生边缘应力。但是当边缘处的局部材料发生屈服进入塑性变形阶段时，上述这种弹性约束就开始缓解，因而原来不同的薄膜变形便趋于协调，结果边缘应力就自动限制。这就是边缘应力的自限性。 边缘应力的特点： 边缘应力与薄膜应力不同，薄膜应力是由介质压力直接引起的，而边缘应力则是由联接边缘两部分变形协调所引起的附加应力，它具有局部性和自限性，通常把薄膜应力称为一次应力，把边缘应力称为二次应力。根据强度设计准则，具有自限性的应力，一般使容器直接发生

6、破坏的危险性较小。 对边缘应力的处理 ： (1) 在边缘区作局部处理。由于边缘应力具有局部性，在设计中可以在结构上只作局部处理。例如，改变连接边缘的结构，如图3-28所示；边缘应力区局部加强；保证边缘区内焊缝的质量；降低边缘区的残余应力（如进行消除应力热处理）；避免边缘区附加局部应力或应力集中，如不在连接边缘区开孔等。对边缘应力的处理 ： (2) 只要是塑性材料，即使边缘局部某些点的应力达到或超过材料的屈服点，邻近尚未屈服的弹性区能够抑制塑性变形的发展，使塑性区不再扩展，故大多数塑性较好的材料制成的容器，例如低碳钢、奥氏体不锈钢、铜、铝等压力容器，当承受静载荷时，除结构上作某些处理外，一般并不

7、对边缘应力作特殊考虑。 但是，某些情况则不然。例如，塑性较差的高强度钢制的重要压力容器，低温下铁素体钢制的重要压力容器，受疲劳载荷作用的压力容器等。对于这些压力容器，如果不注意控制边缘应力，则在边缘高应力区有可能导致脆性破坏或疲劳破坏。因此必须正确计算边缘应力。对边缘应力的处理 ： (3)由于边缘应力具有自限性，它的危害性没有薄膜应力的大。薄膜应力随着外力的增大而增大，是非自限性的。如前所述，具有自限性的应力属二次应力。当分清应力性质以后，在设计中考虑边缘应力可以不同于薄膜应力。实际上，无论设计中是否计算边缘应力， 在边缘结构上作妥善处理显然都是必要的。 弹性基础梁的弯曲关系 弹性基础梁搁置在

8、能连续支承而且具有弹性的基础上，它在集中载荷下即产生弯曲。这种梁上任一点的挠度y与作用于该点梁上的的基础反力q成正比k（基础反力与挠度的比例常数）。即q=ky弹性基础梁的弯曲关系 现从梁上该任意点处割取长度为dx的一个微体，微体两端的受力情况如图所示。由于微体处于平衡状态，则作用于微体上垂直方向上的力的总和应为零，即： FQ-(Q+dQ)+kydx0 dQ/dx=ky 同样，取距载荷作用点为x+dx处的弯矩总和为零，则得：MQdx+M-(M+dM)+kydxdx/2=0略去高阶微量(dx)2, 则得：dM/dxQ则 kydxdQdxMd22弹性基础梁的弯曲关系设距载荷作用点为x处梁的中性层的曲

9、率半径为，则其倒数、1/便是梁轴的曲率。由材料力学得知弯曲梁的弹性曲线和弯矩M成正比，而与它的惯性矩J及材料的纵性弹性模量E成反比。关系式为：其中乘积EJ成为梁的弯曲刚度。EJM1弹性基础梁的弯曲关系而由微积分中得到曲率的表达式为：对于很小的挠度y，dy/dx的数值要比1小得多，略去分母中的(dy/dx)2即得：2322211dxdydxyd221dxyd弹性基础梁的弯曲关系设因此弯曲梁的弹性曲线的一般方程式是：对M两次求导得：即令 ，即： 得到这个四阶微分方程。EJMdxyd224422dxydEJdxMd044EJkydxyd44EJk44EJk则04444ydxyd弹性基础梁的弯曲关系这

10、个四阶微分方程的特征方程为：其根为： 。 故微分方程的解为：在远离载荷处，即x=，y0,则C1,C2都为零。即得弹性基础梁的弯曲关系式为: 常数c3和c4可以由梁的各种特殊条件求出。 0444)1)1 (ii（和)sincos()sincos(4321xCxCexCxCeyxx)sincos(43xCxCeyx圆筒形壳体的弯曲设有一圆筒形壳体受到相邻壳体的牵制而产生半径缩小(或增大)的变形，从壳体中割取一个纵向微条，则壳体发生半径缩小的变形就是这纵向微条产生挠度y=的弯曲变形。在壳体中与联接处的距离不同的各个截面上半径缩小是不同的，也就是这一微条沿着它的长度方向具有不同的挠度y。设微条的宽度为

11、一单位长度，y=在距离联接处为x的壳体半径的形变=，s=壳体厚度。M0、Q0=在联接处单位圆周长度的弯矩和剪力。因为壳体的纵向微条的弯曲与弹性基础梁的弯曲相似。则弯曲曲线方程仍为： )sincos(43xCxCeyx圆筒形壳体的弯曲它的常数C3、C4则可以由微条的端点条件(联接处)求出：以平板式壳壁的弯曲刚度D代替梁的弯曲刚度EJ，则 式中 为平板弯曲刚度。03300220)()(xxdxydDQdxydDM)1 (1223EsD圆筒形壳体的弯曲将曲线方程求导带入得：代入弯曲曲线方程，求出挠度y、 DMCMQDC20400332);(213322dxydDQdxydDMdxdy、剪力、弯矩斜率

12、圆筒形壳体的弯曲下面是具体表达式：令：A1=e-x(cosx+sinx)；A2=e-xsinx A3=e-x(cosx-sinx)；A4=e-xcosx20302010102403401203024032043022422222AMAQQAQAMMAkQAkMADMADQAkMAkQADMADQy四个公式都具有振幅迅速递减的衰减波特征，它的变化周期为正弦或余弦函数的1/，即等于2/，所以梁的挠度距离加载点半个周期/处，就已经非常小了，所以只要圆筒体长度大于2/（1个周期）就可以应用上面公式来计算它的挠度，斜率，弯矩和剪力。圆筒形壳体的弯曲挠度y及斜率在端点(即壳体联接处)x=0处，取得最大值，

13、kQkMkMkQyyxx02030max0200max24)(22)(圆筒形壳体的弯曲挠度y及斜率在端点(即壳体联接处)x=0处，取得最大值，xA1A2A3A4xA1A2A3A401.0000 0.0000 1.0000 1.0000 2.20.0244 0.0896 -0.1548 -0.0652 0.10.9907 0.0903 0.8100 0.9003 2.30.0080 0.0748 -0.1416 -0.0668 0.20.9651 0.1627 0.6398 0.8024 3/40.0000 0.0670 -0.1340 -0.0670 0.30.9267 0.2189 0.48

14、88 0.7077 2.4-0.0056 0.0613 -0.1282 -0.0669 0.40.8784 0.2610 0.3564 0.6174 2.5-0.0166 0.0491 -0.1149 -0.0658 0.50.8231 0.2908 0.2415 0.5323 2.6-0.0254 0.0383 -0.1019 -0.0636 0.60.7628 0.3099 0.1431 0.4530 2.7-0.0320 0.0287 -0.0895 -0.0608 0.70.6997 0.3199 0.0599 0.3798 2.8-0.0369 0.0204 -0.0777 -0.0

15、573 /40.6448 0.3224 0.0000 0.3224 2.9-0.0403 0.0132 -0.0666 -0.0534 0.80.6354 0.3223 -0.0093 0.3131 3.0-0.0423 0.0070 -0.0563 -0.0493 0.90.5712 0.3185 -0.0657 0.2527 3.1-0.0431 0.0019 -0.0469 -0.0450 1.00.5083 0.3096 -0.1108 0.1988 -0.0432 0.0000 -0.0432 -0.0432 1.10.4476 0.2967 -0.1457 0.1510 3.2-0

16、.0431 -0.0024 -0.0383 -0.0407 1.20.3899 0.2807 -0.1716 0.1091 3.4-0.0408 -0.0085 -0.0237 -0.0323 1.30.3355 0.2626 -0.1897 0.0729 3.6-0.0366 -0.0121 -0.0124 -0.0245 1.40.2849 0.2430 -0.2011 0.0419 3.8-0.0314 -0.0137 -0.0040 -0.0177 1.50.2384 0.2226 -0.2068 0.0158 5/4-0.0279 -0.0139 0.0000 -0.0139 /20

17、.2079 0.2079 -0.2079 0.0000 4.2-0.0204 -0.0131 0.0057 -0.0074 1.60.1959 0.2018 -0.2077 -0.0059 4.4-0.0155 -0.0117 0.0079 -0.0038 1.70.1576 0.1812 -0.2047 -0.0235 4.6-0.0111 -0.0100 0.0089 -0.0011 1.80.1234 0.1610 -0.1985 -0.0376 3/2-0.0090 -0.0090 0.0090 0.0000 1.90.0932 0.1415 -0.1899 -0.0484 7/40.

18、0000 -0.0029 0.0058 0.0029 2.00.0667 0.1231 -0.1794 -0.0563 20.0019 0.0000 0.0019 0.0019 2.10.0439 0.1057 -0.1675 -0.0618 9/40.0012 0.0006 0.0000 0.0006 现在求式中的值：我们曾经令 ，现以平板或壳壁的弯曲刚度 来代替梁的弯曲刚度EJ则得： ，其中k未知，求k：单位宽度的微条在距壳体联接处为x处的半径缩小量为，则该截面圆周的相对压缩变形为2=/r，因此在此处产生的环向应力便为 则单位长度上的纵向微条的环向力为 ，由于纵向微条所对的圆周角=1/r，

19、所以环向力的径向分力为 。这个微条的径向分力P反抗微条产生挠度，而且它与挠度成比例地沿着微条的长度分布，它相当于弹性基础梁上的基础反力q，故得 44EJk)1 (1223EsD44DkryErEE2SryEN yrESrNP21yrESqky22rESk rsrsDk285. 13 . 0)1 (34424圆筒体与封头联接组成的容器的应力 M为x处的单位圆周长度的弯矩，每单位圆周长度的抗弯模量w=s2/6（上述公式也适用于封头在联接处附近的窄小区域。）211116)(sMbT曲应力由于弯曲产生的纵向弯222b226),()()()(sMrEycT增大取正缩小取负周向弯曲应力由于纵向弯曲而产生的

20、力而产生的压缩或拉伸应壳体直径被缩小或增大圆筒体与封头联接组成的容器应力计算的步骤第一步：把圆筒体与封头看作是在联接处分离开的独立壳体；第二步：把圆筒体与封头的有关系数如：、D以及由它们组成的公式中Q0及M0项的系数先计算出来。第三步：分别计算出圆筒体与封头由于压力的作用而产生的径向形变及斜率(角转移)。第四步：令在联接处圆筒体的径向总形变(包括由于压力产生的和由弯矩M0及径向剪力Q0引起的径向形变)与封头总形变相等，并令在联接处圆筒体的总弯曲变形曲线的斜率与封头的总弯曲斜率相等，并联解这两等式而求得M0及Q0;第五步：根据M0，Q0写出距联接处为X的壳体截面上的弯矩方程式第六步：写出距联接处

21、X的壳体截面上合成应力的方程式。第七步：求出最大的径向合成应力和周向合成应力以及所在的位置，必要时应绘出应力变化曲线。2010AQAMM一般取圆筒体与半球形封头具有相同的厚度s。圆筒体在内压P作用下,经向应力: 周向应力: 周向相对形变亦即半径的相对形变为:圆筒体的半径增量为:(一)圆筒体与半球形封头联接sPR21sPR2)2(2122sEPREE)2(222sEPRR一般取圆筒体与半球形封头具有相同的厚度s。半球形封头经向和周向应力均为:半径相对形变为:半球形封头的半径增量为:(一)圆筒体与半球形封头联接)1 (2122sEPREEsPR221)1 (222sEPRR在联接处，圆筒体由于弯矩

22、M0及径向剪力Q0而产生的挠度、斜率由公式得:半球封头由于弯矩M0及径向剪力Q0而产生的挠度、斜率，因为这种弯曲变形只局限在很短的一段筒节中，所以也可按求得，因为R，S与圆筒体相同，则也相同。则它的挠度和斜率应分别为：(一)圆筒体与半球形封头联接0203024002222)(MkQkAkMAkQyx020310240302424)(QkMkAkQAkMx020022)(MkQkyx0203024)(QkMkx在联接处，筒体与封头的总弯曲变形曲线的斜率必相等得： 所以 M0=0在联接处圆筒体与封头的径向总形变必相等 得：(正值，说明假设方向是正确的)(一)圆筒体与半球形封头联接00)()(xx0

23、2022)1 (22)2(2QksEPRQksEPR8/8220PREsksEkPRQ020302032424QkMkQkMk00)()()(xxyy则在与联接处距离为x的任意截面的挠度及弯矩按求得：求出圆筒体的经向合成应力: 经向合成应力的第二项弯曲应力与函数值A2=e-xsinx有关，A2 在x=/4取得最大值;周向合成应力与函数A2，A4之值有关，约在x=/2处有最大值。半球形封头应力计算方法相同，但周向合成应力第一项应为球体薄膜应力，第二项取正值，因为径向剪力使它的半径增大。(一)圆筒体与半球形封头联接4230240422AEsPRAkMAkQy2220108APAQAMM222211

24、4326AsPsPRsMT22242224346AsPAsPRsPRsMyRET(一)圆筒体与半球形封头联接2222114326AsPsPRsMT22242224346AsPAsPRsPRsMyRETrsrsDk285. 13 . 0)1 (34424 A1=e-x(cosx+sinx)； A2=e-xsinx A3=e-x(cosx-sinx)； A4=e-xcosx例题 半球形封头的钢制圆筒型容器，半径R1m，设计压力p2MPa，封头和筒体的壁厚S0.02m，计算圆筒体的合成应力并求出最大合成应力及所在位置。 例题 086. 902. 01285. 1285. 1RS22222221142.455056.8202. 042302. 02124326AAAsPsPRsMTxeAxsin24xmx086. 0086. 944解：， 2=82.56圆筒体的经向合成应力: 式中，由计算壳体不连续应力所用的函数表可查得，当时，A2取得最大值0.3224，故得最大经向合成应力为：（1T）max5