压力容器在化工、核能、石油、轻工及医疗等领域被广泛运用,随着科学技术的进步,压力容器朝着高压大型化的方向发展[1-3],其工作条件苛刻[4-5],受力情况的复杂性和结构的多样性[6]对容器的开盖和密封结构设计提出更严格的要求[7]。压力容器设计中,封头与筒体的连接结构有多种形式,如螺栓法兰连接系统[8-9],传统的螺栓法兰系统在高压条件下,内压引起的轴向力较大,一旦垫片的回弹性不能满足要求,甚至可能发生垫片与法兰分离的情况,导致密封失效[10-11]。此外,在高压容器向大型化、高参数化方向发展中,传统螺栓法兰系统受主螺栓截面积限制,筒体端盖和法兰端盖的径向截面尺寸受限。卡箍作为一种无主螺栓连接的紧固结构,尺寸更为紧凑[12]。采用螺栓连接的分瓣式卡箍高压容器,其主要零件包括筒体、封盖和分瓣式卡箍,分瓣式卡箍是此类压力容器的关键结构,通过自紧密封结构配合使用,达到容器内部密封效果[13-14],具有结构紧凑、便于制造、操作方便等优点[15-18],其性能的好坏直接影响产品质量与生产效率[19]。四瓣卡箍形式如图1所示。
<G:\武汉工程大学\2024\第4期\孔德洲-1.tif>[(b)][(a)][卡箍][筒体][封盖]
图1 高压容器四瓣式卡箍示意图:
(a) 法兰、卡箍剖视图, (b) 卡箍侧视图
Fig. 1 High-pressure vessel four-lobe clamp structure:
(a) flange, clamp section view, (b) clamp side view
在内压作用下,高压容器卡箍将筒体与封盖的端面压紧,以保证筒体和封盖间的密封性。通常将卡箍等分为多瓣结构,如两瓣式、三瓣式和四瓣式,瓣间通过螺栓或其他结构拉紧,卡箍内锥面被压在封盖和筒体的外锥面上,从而使筒体与封盖的端面紧密贴合。若卡箍间连接螺栓受力(即卡箍间分瓣力)过大,则封盖和筒体分离导致压力泄漏。
为了保证高压大型容器卡箍连接机构刚度及强度,提出了一种基于刚体向量变换法的卡箍间分瓣力解析方法。同时在实例基础上建模,利用ANSYS Workbench软件进行有限元静力学分析,采用提取支反力的方法提取卡箍分瓣面间的分瓣力与理论分析结果作对比。该分瓣力的解析方法能准确计算分瓣式卡箍间的分瓣力,为高压大型容器卡箍的螺栓结构或其他类似夹紧结构的设计提供理论依据。
1 卡箍受力理论分析计算
容器内压强为[p],容器承受内压面积为[A],则力[F=pA]作用在筒体和封盖间使筒体和封盖分离,该力通过锥面由卡箍抵消。
以卡箍为受力体进行分析,将作用在卡箍内锥面上的力均布到锥面中间的圆环曲线(称为受力圆环)上,圆环直径为[D]。卡箍单侧锥面受力如图2所示,其中[α]为锥面角(锥面与底平面夹角),[f]为受力圆环上线比压(单位长受力),与锥面正交,[fH]为受力圆环上的线比压沿轴向分量。根据力平衡有:
[fH=pA/(πD)] (1)
而
[f=fH/cosα] (2)
其中,[f]在受力圆环处构成与卡箍锥面正交的力圆锥。
<G:\武汉工程大学\2024\第4期\孔德洲-2.tif>[α]
图2 分瓣式卡箍锥面受力图
Fig. 2 Cone force diagram of split clamp
建立右手坐标系{o-xyz},其中z向为卡箍轴向,筒体端面处于xoy平面内,如图3所示。
由图3可知,受力圆环在xoy平面内。平面yoz与受力圆环相交点处线比压为[f(0)],其中线比压的模[f=f],[f]绕z轴的转动即为受力圆环上其他点处的线比压。在yoz平面内,[f(0)=f=0fsinαfcosαT],绕z轴转动[θ]角,可得受力圆环任意位置处的线比压,如式(3)所示。
[f(θ)=cosθ-sinθ0sinθ cosθ00 01]
[f(0)=-fsinαsinθ fsinαcosθfcosα] (3)
在[f]绕z轴转动[θ]角处取受力圆环弧长l的微分[dl=D2dθ],[f(θ)]的轴向分力[fcosα]沿整个受力圆环积分,得[f(θ)]轴向分力之和:
[02πfcosαdl=02πfcosαD2dθ=fH?πD=F] (4)
即[f(θ)]的轴向分力之和与力[F]平衡。
同样,[f(θ)]的径向分力沿受力圆环积分,有:
[02π-fsinαsinθ fsinαcosθdl=00] (5)
即整圈卡箍所受径向合力自平衡。
卡箍分为多瓣后,将积分区间限定为卡箍等分角(称为分瓣角),积分可得每瓣卡箍的受力。定义分瓣面方向与分瓣面正交且指向体外。在定义的坐标系下,分瓣面方向总在xoy面内,记为[s=cosκ,sinκ,0T],其中[κ]为分瓣面方向与ox轴夹角。将受力圆环上某点线比压[f(θ)]分解到分瓣面方向,在分瓣角内积分得到卡箍分瓣面处的力。则[f(θ)]沿[s]方向的分力为:
[fs(θ)=fsinαsin(κ-θ)cosκ, sinκ, 0T] (6)
设ox轴与2个卡箍分瓣面中1个面正交(如图3所示),记卡箍分瓣数为[n],则积分区间为[0,κ0=2π/n],2个分瓣面的角度分别为[κ1=0]和[κ2=π+2πn]。2个分瓣面的受力分别为:
[F1=0κ0fsinαsin(0-θ)1, 0, 0TD2dθ] (7)
[F2=0κ0fsinαsin(π+κ0-θ)cos(π+κ0)sin(π+κ0)0D2dθ]
(8)
积分得
[F1=D2fsinα(cosκ0-1)1, 0, 0T] (9)
[F2=D2fsinα(cosκ0-1)-cosκ0-sinκ00] (10)
由式(9)和式(10)得分瓣面受力的幅值为:
[F1=F2=D2fsinα(1-cosκ0)] (11)
式(11)所反映的对称性与受力圆环的对称性及等分的卡箍分瓣相对应。将式(1)和式(2)代入式(11),则沿分瓣面方向的分瓣力[Fs]为:
[Fs=F1=F2=pA2πtanα(1-cosκ0)] (12)
即螺栓拉紧卡箍分瓣时,螺栓应在式(12)所计算分瓣力的作用下正常工作。
2 卡箍结构的有限元分析
2.1 卡箍结构模型建立
以筒体直径为1 200 mm的四瓣式卡箍为例,采用三维建模软件Solidworks建立卡箍连接结构模型及部分筒体模型,为便于分析,选取每瓣卡箍模型的1/2进行建模,将卡箍与两法兰接触面配合,所建模型如图4所示,四瓣式卡箍将封头法兰与筒体法兰轴向夹紧闭合。
<G:\武汉工程大学\2024\第4期\孔德洲-4.tif>[封头法兰][卡箍][连接螺栓][筒体法兰]
图4 四瓣式卡箍模型示意图
Fig. 4 Schematic diagram of the four-lobe clamp model
2.2 接触面设置
将三维模型导入ANSYS Workbench中进行静力学分析,为了使分析结果更加符合实际情况,在分析之前需要根据实际的接触面移动情况对接触对的类型进行设置:封头法兰和筒体法兰接触面设置为不分离类型;封头法兰齿边与卡箍接触面以及筒体法兰齿边与卡箍接触面由于在卡箍闭合夹紧与开启放松时法线方向无分离,切线方向发生轻微的无摩擦运动,因此设置为不分离接触类型;连接螺栓与卡箍之间法线方向可分离,切线方向存在相对滑动及摩擦力,因此设为摩擦接触类型,摩擦系数为0.5。
2.3 模型网格划分
在添加载荷和边界条件求解之前需要对模型划分网格,网格尺寸设为5 mm,采用六面体网格划分方法,模型的网格划分如图5所示。
<G:\武汉工程大学\2024\第4期\孔德洲-5.tif>
图5 模型网格划分示意图
Fig. 5 Schematic diagram of model grid division
2.4 边界条件和载荷设定
对模型的边界条件和载荷进行设定,筒体法兰远端做轴向位移约束,封头法兰截面处施加平衡载荷,对称表面设置为无摩擦支撑约束,螺栓端面施加无摩擦支撑约束,在封头法兰和筒体法兰的内表面添加100 MPa内压。
2.5 结果提取
在分析结果中插入支反力选项,提取螺栓端面处支反力,并与理论计算值对比。
2.6 实例计算
将筒体内径作为筒体的受压直径(受压直径约为密封圈圈径,稍大于筒体内径),锥面角取为5°,分别考虑筒体直径800、1 200、1 500 mm与内压50、100、150 MPa工况下的两瓣式、三瓣式与四瓣式卡箍分瓣力计算,理论计算值与有限元计算值如表1所示。
表1表明了分瓣式卡箍分瓣面拉力理论解与有限元解基本相符合,同时展示了筒内压力与筒体直径对分瓣面拉力的影响。
筒体内径及分瓣数相同的条件下,分瓣面拉力随着筒内压力变大而变大;在两瓣式卡箍、筒体直径1 200 mm和内压100 MPa条件下分瓣面拉力为2.96 MN,在四瓣式卡箍、筒体直径1 200 mm和内压100 MPa条件下分瓣面拉力为1.46 MN,约为两瓣式卡箍的2倍。
筒内压力及分瓣数相同的条件下,分瓣面拉力随着筒体直径变大而变大;在两瓣式卡箍、筒体直径1 500 mm和内压100 MPa条件下分瓣面拉力为4.88 MN,在两瓣式卡箍、筒体直径1 200 mm和内压100 MPa条件下分瓣面拉力为2.96 MN,是直径1 500 mm两瓣式卡箍的1.65倍。
3 结 论
为了解决分瓣式卡箍的螺栓连接结构在连接高压容器筒体和封盖的过程中出现刚度过大及强度破坏问题,提出了一种基于刚体向量变换法的卡箍间分瓣力解析方法。首先将筒体和封盖的轴向水平力均布到卡箍内锥面中间的受力圆环上,再计算出受力圆环上任意一点的线比压,将其径向分力分解到分瓣面方向,再沿受力圆环在分瓣角内积分计算出分瓣力。根据计算结果可得分瓣面拉力随着分瓣数的增加而减小。在相同内径下,分瓣面拉力与筒体内压成正比,且内压越大,分瓣面拉力随分瓣数变化越快;在相同内压下,分瓣面拉力与筒体内径成正比,且内径越大,分瓣面拉力随分瓣数变化越快。最后对设计条件下的装配体结构进行仿真实验,理论分析与有限元仿真结果吻合良好,验证了卡箍分瓣面拉力的理论的准确性,为分瓣式卡箍连接结构设计提供理论依据和参考。