第八章外压容器8.1外压薄壁壳体的稳定性分析薄壁壳体承受外压作用时,在壳壁内会产生压缩薄膜应力,壳体的失效形式可能有两种:一种是由于强度不足而发生压缩屈服破坏;另一种是由于刚度不足而发生失稳破坏。所谓失稳,就是当外压达到一定数值时,壳体失去其原有形状,直到壳体被压瘪的现象。8.1.1概述图示是外压薄壁圆筒失稳后被压瘪的实例图示表示了外压薄壁圆筒失稳时在圆筒的横截面上可能出现的几何形状,即在周向呈现波纹状。外压壳体失稳时的相应压力称为临界压力,以pcr表示,此时壳壁中产生的压缩应力称为临界应力,以σcr表示。外压薄壁圆筒的临界压力与圆筒的几何尺寸(主要为筒体的厚度δ、筒体外直径Do和计算长度L,其中计算长度L是指筒体上两相邻支撑线之间的距离)及其材料性能(主要为材料的弹性模量E和泊松比μ)有关,此外,载荷的均匀性和对称性、边界条件等因素也对筒体的临界压力有一定影响。8.1.2外压圆筒的稳定性计算8.1.2.1外压圆筒的分类外压圆筒按其失效情况,可分为长圆筒、短圆筒和刚性圆筒三种。其特点如下:(1)长圆筒所谓长圆筒是指筒体的长径比L/D0值较大,筒体的临界压力pcr与圆筒的长度无关,可以忽略筒体两端边界对筒体的支撑作用,只与筒体材料的机械性能(E、μ)和筒体的厚径比δe/D0有关。(2)短圆筒筒体两端边界对筒体的影响较大,不能忽略。筒体失稳时的临界压力pcr不仅与筒体材料的机械性能(E、μ)和筒体的厚径比δe/D0有关,而且也与筒体的长径比L/D0有关。(3)刚性圆筒所谓刚性圆筒是指筒体的长径比L/D0较小而厚径比δe/D0较大的圆筒。由于筒体的相对厚度较大,故筒体的刚性较好,不会发生失稳,而其破坏原因是由于在外压作用下,在筒体壁内产生的压应力超过材料的屈服极限所致。对于这种圆筒,在设计时只需满足强度要求即可。8.1.2.2外压圆筒的临界压力计算32)(12oecrDEp对于钢质圆筒,μ=0.3,故上式可简写成:3)(2.2oecrDEp(1)长圆筒上式常称为Bresse公式。在临界压力作用下,圆筒壳的周向应力为:值得注意的是,以上计算公式仅适用于弹性失稳的情况。)(1.120DEDpeeocrcr(2)短圆筒eoecrDLDEp0259.2ooeeocrcrDLDEDp5.1)(3.12同样,以上计算公式也仅适用于弹性失稳的情况。上式也称为B.M.Pamm(拉姆)公式。在临界压力作用下,圆筒壳的周向应力为:(3)刚性圆筒刚性筒在外压作用下一般不存在失稳问题。即在外压作用下不会被压瘪,只是由外压引起的最大周向应力超过筒体材料的屈服极限时而产生强度破坏。因此,只需校验其强度是否足够就可以了,其强度计算公式与内压圆筒的应力计算公式完全一样。(4)外压圆筒的临界长度对于已知直径和壁厚的圆筒,临界长度就是区分长、短圆筒的界限,用Lcr表示。将以上长、短圆筒的临界压力计算公式若L=Lcr,则按两式求得的临界压力应该相同,便得到了临界长度为:当圆筒计算长度L>Lcr时,按长圆筒公式计算。反之,按短圆筒公式计算。ecrDDL0017.1[例8-1]有一圆筒,材料为Q245R,弹性模量E=2.1×105MPa,内直径Di=1000mm,厚度δ=9.7mm,计算长度L=20m,圆筒内的介质无腐蚀性、常温操作,试求其临界压力pcr。解(1)计算圆筒的外直径Do=Di+2δ=1000+2×9.7=1019.4mm(2)计算该圆筒的临界长度由于L=20m>Lcr,故此圆筒受外压时为长圆筒。(3)计算该圆筒的临界压力mmmDDLecr23.12122277.94.10194.101917.117.100MPaDEpocr40.0)4.10197.9(101.22.2)(2.23538.2.1外压圆筒的设计(1)图算法的由来根据前述外压圆筒的稳定性计算,对于仅受均匀外压的圆筒临界压力,无论是长圆筒或短圆筒,均可用如下公式表示:式中δe——外压圆筒的有效厚度,mm;K——外压圆筒的几何特征系数,对于长圆筒,K=2.2;对于短圆筒,K值则与L/Do和Do/δe有关。3)(oecrDKEp8.2外压容器的工程设计圆筒在临界压力作用下,筒壁产生相应的周向应力σcr和周向应变εcr分别为从上式可看出,周向应变εcr只是L/D0和D0/δe的函数,与材料性质无关。在GB150中,εcr采用字母A表示,称为外压应变系数;以A为横坐标,L/D0为纵坐标,D0/δe为参变量,得到适用于所有材料的外压应变系数A的曲线图,如图14-6所示。2)(22oeeocrcrDEKDp2)(2oecrcrDKE在图14-6的曲线中,与纵坐标平行的直线簇表示长圆筒情况,失稳时的应变量与L/D0无关;图下方的斜线簇属短圆筒情况,失稳时的应变量与L/D0和D0/δe都有关。圆筒的许用外压力令,称为外压应力计算系数,GB150取圆筒的稳定性安全系数m=3,将B和m代入上式可得3)(1oecrDKEmmppcroeoeeomDKEmDmKEDp2)(22)(][22eoDpB][crB32由于A=εcr,而,故B与A的关系即为与εcr的关系。若利用材料单向拉伸应力—应变关系,将纵坐标乘以,就可作出B与A的关系曲线,即外压应力计算系数B的计算图。若圆筒失稳时发生了塑性变形,工程上通常采用正切弹性模量,即应力—应变曲线上任一点的斜率,因此,图算法对于非弹性失稳也同样适用。图14-7至图14-9为几种常用钢材的外压应力计算系数B的计算图。crB32cr3232ddEt(2)工程设计方法对于D0/δe≥20的圆筒和管子,其步骤如下:①假设一个名义厚度δn,则有效厚度δe=δn-C1-C2,计算出L/Do和Do/δe。②在图14-6的左方找到L/Do值,将此点沿水平方向右移与Do/δe线相交(遇中间值用内插法),过此交点沿垂直方向下移,在图的下方得到系数A。若L/Do值大于50,则用L/Do=50查图;若L/Do值小于0.05,则用L/Do=0.05查图。③确定外压应力计算系数Bⅰ.根据所用材料选用相应的外压应力计算系数B曲线图(图14-7~图14-9),由A值查取B值。ⅱ.若A值超出设计温度曲线的最大值,则取对应温度曲线右端点之纵坐标值为B值;ⅲ.若A值小于设计温度曲线的最小值,即所得A值落在设计温度曲线的左方,则按下式计算B值:④根据B值,按下式计算许用外压力[p]:比较计算外压力pc和许用外压力[p],若pc≤[p]且较接近,则所假设的名义厚度δn合理;若p>[p],则需再假设一较大的名义厚度δn,重复上述步骤直到满足设计要求为止。AEB32)/(][eoDBp对于D0/δe<20的圆筒和管子,圆筒和管子,求取B值的计算步骤大致相同,这里不作介绍,必要时可参阅GB150。(3)外压容器有关设计参数①设计压力p外压容器(例如真空容器、液下容器和埋地容器)的设计压力确定时应考虑在正常工作情况下可能出现的最大内外压力差。ⅰ.确定真空容器的壳体厚度时,设计压力按承受外压考虑。当装有安全控制装置(如真空泄放阀)时,设计压力取1.25倍最大内外压力差或0.1MPa两者中的低值;当无安全控制装置时,取0.1MPa。ⅱ.由两个或两个以上压力室组成的容器,如夹套容器,应分别确定各压力室的设计压力。确定公用元件的设计压力时,应考虑相邻室之间的最大压力差。②稳定性安全系数m上节外压圆筒的稳定性计算中已经得到了外压圆筒的临界压力计算公式,但均是在一定假设条件下按理想状态推导得到的,与实际情况有一定的差别,而制造技术(如机加工和焊缝结构等)所能保证的质量使得不能加工出绝对形状的圆筒。这些因素都直接影响圆筒的临界压力,因此,在计算许用设计外压力时,应考虑一定的安全裕度,即:式中[p]——许用外压力,MPa;m——稳定性安全系数,GB150规定:对于圆筒,取m=3。mppcr][③计算长度L外压圆筒的计算长度L是指筒体上两相邻支撑线之间的距离,其中支撑线就是刚性构件所处的环向线。通常封头、法兰、加强圈等均可视作刚性支撑构件。对于凸形封头,应计入其直边高度以及封头曲面深度的1/3。不同结构的外压圆筒计算长度取法见图14-10。其他设计参数的选取与内压容器相同。8.2.2外压球壳设计推导过程:钢制球形壳体弹性失稳的临界压力为:工程上:图算法。2)/(21.1oecrREp取稳定性安全系数m=14.52,得球壳许用外压力:2)/(0833.052.14][eocrREppeoR]P[B令eoRPEAB][32根据)/(32][eoREAp得将[p]代入式得)/R(125.0Aeo由B和[p]的关系式得半球形封头的许用外压力为:)/R(B]p[eoa.假定名义厚度δn,令δe=δn-C,用式计算出A,根据所用材料选用厚度计算图,由A查取B,再按式计算许用外压力[p]。图算步骤:b.若A值落在设计温度下材料线左方,用式计算[p]。若[p]≥pc且较接近,则该封头厚度合理;否则重新假设δn,重复上述步骤,直到满足要求为止。不用几何算图外压容器封头结构形式与内压容器的封头相同,主要包括有半球形、椭圆形、碟形、球冠形和锥形封头等。外压容器的封头设计除了满足强度条件外还须满足稳定性要求。原则:凸形封头按球壳计算、锥形封头按筒体计算。(1)半球形封头受外压的半球形封头厚度计算,与受外压球壳的计算步骤相同。8.2.3外压封头设计(2)椭圆形封头受外压的椭圆形封头的厚度计算,与球壳的计算步骤相同,唯一不同的是R0=KD0,对于标准半椭圆封头来说R0=0.9D0。(3)碟形封头受外压的碟形封头的厚度计算与受外压球壳的计算步骤相同,其中R0为碟形封头的球面部分外半径。(4)其他封头对于受外压的球冠形封头、锥形封头的设计详见GB150。[例8-2]某圆筒形容器,其内直径为2400mm,长14000mm,两端为标准椭圆形封头,直边高度为40mm,材料为Q345R,最高操作温度为200°C,其弹性模量E=1.86×105MPa,真空下操作,无安全控制装置,介质为腐蚀性较小的气体,可取腐蚀裕量C2=2mm。试用图算法求筒体和封头的厚度。解:(1)筒体厚度①设筒体名义厚度δn=22mm,题中已知C1=0.30mm,C2=2mm,则δe=δn-C1-C2=22-0.30-2=19.70mmDo=2400+2×22=2444mmL——计算长度,等于圆筒长度加上两封头的直边高度及曲面深度的1/3,即:L=14000+2×40+2×1/3×600=14480mmpc——计算外压力,真空下操作且无安全控制装置,取pc=0.1MPa。②计算L/Do、Do/δe③在图14-6的左方找出L/Do=5.92的点,将此点沿水平方向右移,与Do/δe=124.06线相交于一点,过此交点垂直下移,在图的下方得到系数A=0.00014。④材料为Q345R,选用算图14-8,A=0.00014落在设计温度曲线的左方,故按式(14-12)计算系数B:92.5244414480oDL06.12470.192444eoD36.17101.8600014.032325AEB⑤根据B值,计算许用外压力[p]:由于pc<[p]且较接近,则所选筒体名义厚度合适,即δn=22mm。(2)封头厚度①设封头名义厚度δn=8mm,则δe=δn-C1-C2=8-0.30-2=5.70mm由表14-6查得K1=0.90,于是Ro=K1Do=0.90×2444=2199.6mm②利用式(14-15)计算系数AMPaDBpeo14.006.12436.17)/(][89.38570.56.2199eoR00032.089.385125.0)/(125.0eoRA③材料为Q345R,选用算图14-8,