高压管汇活动弯头接触特性分析

摘要：为探究高压管汇活动弯头活动节在现场油气田使用过程中滚道易出现点蚀现象的原因，选取某厂型号为?130.2mm—105 MPa（51／8in—15 000 psi）的大通径活动弯头作为研究对象，采用有限元分析软件研究不同滚道轴向间距、径向间距、滚道半径对滚道接触应力的影响规律。研究结果表明：随着活动弯头滚道轴向间距的增大，最大接触应力呈现增长趋势；随着径向间距的增大，最大接触应力总体呈现上升趋势；随着滚道半径的增大，最大接触应力呈现先降低后增大的趋势，在滚道半径为9.66mm时达到最低。针对活动弯头活动节第一滚道相较于其余滚道接触应力较大的现象，通过对弯头结构分块解析，发现弯头偏载原因主要为弯头不同部位的结构刚度存在差异，且管壁内压是导致弯头偏载的载荷来源。指出管壁内压直接作用于外接头，导致外接头滚道最大接触应力比内接头大。研究结果可为活动弯头活动节进一步优化设计提供参考。

关键词：高压管汇；活动弯头；接触应力；滚道；偏载；结构优化

中图分类号：TE973 文献标识码：A DOI：10.16082／j．cnki．issn．1001-4578.2022.12.017

Contact Characteristics of Active Elbow in High-Pressure Manifold

Zhang Guoyou

（Sinopec SJ Petroleum Machinery Co.）

Abstract： In order to investigate the raceway pitting corrosion of the active elbow in high-pressure manifold in oil and gas fields， a large-diameter active elbow of Model 130.2 mm-105 MPa （5/8 in-15 000 psi） was taken as an example to investigate the influences of different raceway axial spacing， radial spacing and raceway radius on the raceway contact stress by using the finite element analysis software.The results show that，as the raceway axial spacing of active elbow increases， the maximum contact stress increases; as the raceway radial spacing increases， the maximum contact stress rises generally; as the raceway radius increases，the maximum contact stress decreases first and then increases， which reaches the minimum value when the raceway radius is 9.66 mm. In view of the phenomenon that the contact stress of the first raceway of the active elbow is higher than that of other raceways，the elbow structure was divided into 6 regions for analysis. It is found that the eccentric load of the elbow is mainly re- sulted from the difference of structural rigidity in different parts of the elbow， and the internal pressure of pipe wall is the source of the eccentric load. The intemal pressure of pipe wall directly acts on the outer joint; as a result， the maximum contact stress of the raceway of the outer joint is larger than that of the inner joint.The results provide references for the optimization design of the active elbow.

Keywords：high-pressure manifold;active elbow; contact stress;raceway;eccentric load1;structural opti- mization

0引言

压裂技术在油气田中使用日益广泛，在非常规油气开发中［2］，通常需要面对复杂的井身结构、复杂的岩性和复杂的流体性质，这对压裂设备提出了更高的要求，需要有良好的耐高压、抗腐蚀、耐冲蚀特性。但是目前高压管汇在实际工作中的寿命远远短于其他的压裂设备，严重制约了油气田整体开发效益，所以需要对高压管汇易损件进行研究［3-7］。

为了延长活动弯头的使用寿命，许多学者进行了大量研究并取得了一定进展。徐晓东等［8］通过对活动弯头结构和加工工艺改进来延长活动弯头使用寿命。敬佳佳等［9］利用CFD软件分析放喷管角度、位置、放喷量、含砂率、颗粒形状对壁面冲蚀率的影响，指出含砂率对最大冲蚀率的影响较大，弯管位置对最大冲蚀率的影响较小，放喷管最大冲蚀率随含砂率、颗粒形状系数的增大而增大，随弯管角度和直管段长度增大而减小，随着放喷量的增加最大冲蚀率先增大后减小。郭登明等［10］针对现有活动弯头出现明显压痕和点蚀的情况，以20CrNiMo活动弯头为研究对象，将活动弯头进行渗碳＋等温淬火＋回火热处理，观察活动弯头金相组织的变化，发现工件表面为下贝氏体金相组织，该组织具有良好的耐磨性和抗腐蚀性，内部为下贝氏体和马氏体的混合金相组织，该组织可以有效防止裂纹扩散。

以上研究多以弯头加工工艺、弯头最大冲蚀率、弯头热处理工艺为主，对大通径四滚道活动弯头接触特性规律研究有待完善。笔者以大通径四滚道活动弯头为研究对象，采用增广拉格朗日乘子法，分析不同滚道轴向间距、径向间距、滚道半径对活动弯头最大接触应力的影响规律。笔者对滚道接触应力进行分析总结，以期为活动弯头的设计、优化提供指导。

1 有限元接触算法

1.1 罚函数法

罚函数法相当于在接触面与目标面之间引入弹簧系统。在未接触时，弹簧系统不发生作用；在接触发生时，弹簧系统就会阻止接触物体相互嵌入［12］。为方便系统计算，在不含接触总位能基础上引入一个惩罚势能：

式中：IIp（U）为惩罚势能，J；Ep为惩罚因子，N／m；P为嵌入深度，m；U为节点位移，m。

这样，接触问题就等价于无约束优化问题：

式中：II＊（U）为接触总位能，J；II（U）为不含接触总位能，J。

罚函数法受接触刚度的影响较大。接触刚度越大，穿透量越小；接触刚度越小，穿透量就越大。理论上而言，当我们将接触刚度设定到无穷大时，其穿透量即可以为0，但是刚度矩阵不可能无穷大，否则刚度矩阵呈现病态。

1.2 拉格朗日乘子法

拉格朗日乘子法［］是通过类比拉格朗日法的数学意义，以g（x）为约束条件，f（x）作为目标函数，引入拉格朗日乘子λ来构造公式：

在接触非线性问题上，拉格朗日乘子法以接触力的拉格朗日乘子λ乘以无侵彻条件，构造接触约束条件的附加泛函：

式中：IIc（U，λ）为含接触约束条件的附加泛函，J；λ为拉格朗日乘子；g为接触约束条件。

最终将约束条件极值问题转化为无条件极值问题：

拉格朗日法不再通过定义刚度来计算接触力，而是将接触力作为一个自由度，扩大了刚度矩阵。该算法优点是精度高；缺点是失去系数矩阵正定性，不易收敛。

1.3 增广拉格朗日乘子法

在增广拉格朗日乘子法［4］下，系统的总位能由不含接触约束条件的总位能II（U）和惩罚势能IIp（U）以及含接触约束条件的附加泛函IIc构成，构造修正的势能泛函如下：

增广拉格朗日是先将程序按照罚函数法进行开始计算，通过赋予允许穿透的最大值。如果在计算迭代过程中出现穿透量大于许用穿透量，则对刚度矩阵进行更新，将各个接触单元的接触刚度加上接触力乘以拉格朗日乘子的数值，形成不断更新刚度矩阵的罚函数法，也就是拉格朗日法。随着刚度矩阵的不断更新，最终穿透量将小于允许值，计算则停止。极大地克服了罚函数法的病态矩阵，以及拉格朗日法的刚度矩阵零对角元问题。

2 活动弯头有限元接触分析

2.1 物理模型及网格划分

本文选取某厂型号为?130.2mm—105 MPa（518 in—15000 psi）的大通径活动弯头为研究对象，采用Solidworks 对其进行三维建模。该活动弯头为四滚道活动弯头，4个滚道安装钢球数量分别为32、33、34和35颗，其中第一滚道钢球数量为32颗。因为活动弯头活动节整体结构比较复杂，无法对其直接进行接触分析。查阅文献可知，活动弯头第一滚道接触应力较大，所以模型以第一滚道32颗钢球为参考，建立活动弯头活动节1／32分析模型，如图1所示。

在活动弯头使用期间，活动弯头滚道表面容易出现压痕和点蚀。其中第一滚道最严重。因此，需要对活动弯头进行弹塑性接触分析，其中在活动弯头活动节塑性阶段采用双线性模型。材料的密度p＝7 850 kg／m3，弹性模量E＝210GPa，泊松比 v＝0.275，屈服强度σ，＝896 MPa，切线模量E＝21 GPa。

考虑到弹塑性接触分析对网格要求较高，本文采用如下网格划分方式：首先对球体进行切分，切分后可以自动生成规则的六面体网格；为保证滚道网格基本与钢球网格一致，采用球体加密网格方式，以钢球球心为坐标原点，对滚道一圈进行网格加密。对于非接触区域不存在塑性变形，顺延滚道六面体划分方式。网格大小设置为5mm，划分结果如图2所示。

2.2 边界条件

根据活动弯头的实际工况，对活动弯头有限元网格模型施加载荷和边界条件。活动弯头主要承受2种载荷：内压和轴向载荷。内压作用在外接头和内接头内壁面，外接头在管壁内压作用下产生径向位移，钢球在外接头和内接头共同作用下产生挤压变形；轴向载荷是在管内流体压力作用下，将管道看作密闭容器，外接头和内接头都将受到管内流体对管道的横向作用力。在施加载荷时，根据施加工况直接在外接头和内接头管道内表面施加内压105MPa的压力；轴向载荷可以在内接头端面施加，外接头端面采用位移约束，约束外接头端面轴向位移。内接头端面载荷采用均布载荷方式施加，均布载荷p可以通过下式求得：

式中：p.为管壁内压，MPa；S，为内接头外端面内圆的面积，m㎡；S。为内接头外端面外圆的面积，m㎡；d，为内接头外端面的内径，mm；d。为内接头外端面的外径，mm。

活动弯头钢球多，钢球与滚道接触较多，不利于分析。上文对模型进行简化，取1／2进行分析，所以对外接头和内接头侧面需要进行合理约束，本文采用法向无位移约束。施加载荷和边界条件后的模型如图3所示。

2.3 网格无关性验证

通过对数值模型的试算，钢球的法向接触刚度因子与网格数量对弹塑性接触非线性仿真分析结果有较大影响，为确保计算结果的正确性，进行无关性验证。首先计算最大接触应力与网格大小的关系，计算结果如图4所示。

通过对网格无关性验证可以得到，随着网格不断减小，网格量不断增大，结果趋于收敛。网格尺寸在0.65～1.05mm时接触应力波动较为平稳，分布较为均匀。若将网格继续加密，接触应力开始出现不规则，分布不均，开叉等不正常现象。网格尺寸在0.65～1.05mm范围内的最大接触应力为2 264.5 MPa，最小接触应力为2145.2MPa，误差为5.26％，精度满足要求。综合上述分析，本文选取网格大小为1.05mm，在此基础上选择合适的法向接触刚度因子，计算最大穿透量、最大接触应力与接触面法向接触刚度因子的影响，如图5所示。

查阅文献可知［15］法向刚度控制接触面和目标面之间的穿透量。较高的法向刚度会降低贯入量，但会导致整体刚度矩阵的病态和收敛困难。较低的法向刚度可能导致一定数量的穿透，并产生不准确的解决方案。理想情况下，需要一个足够大的法向刚度，使渗透足够小。通过对不同法向刚度分析计算，可以看出随着接触面法向刚度的增加，接触面与目标面的最大穿透量不断减少趋于收敛；最大接触应力逐渐增大直至稳定，综合考虑，选取法向接触刚度因子为10。

2.4 仿真计算结果

采用 Workbench 软件对活动弯头进行弹塑性接触非线性分析，得到滚道接触应力如图6和图7所示。由图6和图7可知，接触应力云图呈现椭圆形，与赫兹理论假设一致，分析结果合理。外接头和内接头第一滚道接触应力都明显大于其余滚道，呈现逐渐减小趋势，与现场第一滚道压痕最为严重的现象一致。对比外接头与内接头滚道最大接触应力可以发现，外接头滚道最大接触应力比内接头大。

2图8为活动弯头结构示意图。本文将结构划分为6个区域。由上文可知，活动弯头所受载荷主要有2种。活动弯头内壁在承受载荷时，活动弯头所有区域在径向上都会外扩，区域1壁厚较大，在区域5共同作用下，不易变形；区域4壁厚较大，在区域6共同作用下，不易变形；相较于区域1、区域4，区域2、区域3径向变形较为容易。在内压作用下，区域3推动区域1，区域4推动区域2，区域3可以过滤内压能量不如区域4，剩余能量都由区域1承担；区域4过滤较多的内压能量，剩余的能量到区域2，又因为区域2径向易变形，而区域1不易变形，所以造成滚道严重偏载。活动弯头在承受轴向载荷时，所有区域在径向上都会有远离钢球的趋势，在第一滚道和第二滚道都有径向易变形区和不易变形区，所以轴向载荷并不是造成滚道偏载的主要原因。

以外接头滚道最大接触应力为判定依据，通过对活动弯头内压和轴向载荷单独分析，再与载荷共同作用进行对比，得出不同载荷下滚道所受最大接触应力，分析结果如表1所示。

由表1可以看出，在轴向载荷单独作用下，滚道最大接触应力相差不大；在管壁内压单独作用下，第一滚道相对第四滚道最大接触应力降低了44.68％；在轴向载荷和管壁内压共同作用下，最大接触应力降低了31.70％。由此可见，轴向载荷对管壁内压导致的滚道接触应力偏载有一定的缓解作用。

3 活动弯头接触应力影响因素分析

由上述仿真可知，活动弯头存在接触应力过大和偏载现象，笔者对接触应力偏载做出了解释，并在此基础上分析活动弯头结构参数对滚道接触应力的影响。弯头结构参数主要包括滚道轴向间距、滚道径向间距和滚道半径。

3.1 滚道轴向间距

在滚道结构参数不变的情况下，选取不同大小的滚道轴向间距分别为27.4、31.4、35.4、39.4和43.4mm。研究滚道轴向间距对活动弯头滚道接触应力的影响规律，整理数值仿真数据可得活动弯头滚道最大接触应力随滚道轴向间距的变化规律，如图9所示。

图9表明，随着活动弯头滚道轴向间距的增大，第一滚道、第二滚道最大接触应力呈现上升趋势，第三滚道最大接触应力上下波动，第四滚道最大接触应力呈现下降趋势。因为随着滚道轴向间距的增大，外接头内壁面承载内压面积增大，单个钢球所需承载的载荷增加，但是区域3过滤内压能力没有区域4强，所以导致活动弯头滚道最大接触应力的增加。

3.2 滚道径向间距

在滚道结构参数不变的情况下，选取不同大小的滚道径向间距分别为0、1.1、2.1、3.1、4.1、5.1和6.1mm。研究滚道径向间距对活动弯头滚道接触应力的影响规律，整理数值仿真数据可得活动弯头滚道最大接触应力随滚道径向间距的变化规律，如图10所示。

图10表明，随着活动弯头滚道径向间距的增大，第一滚道、第二滚道最大接触应力呈现上升趋势，第三滚道最大接触应力上下波动，第四滚道最大接触应力呈现下降趋势。这是因为随着弯头滚道径向距离变大，区域1和区域4的壁厚变大，区域2和区域3的壁厚没有改变，导致区域4承载的内压比例增加，分配给区域2的内压能减少；区域3承载的内压比例基本不变，但是区域1壁厚增加，导致结构刚度增加，更加不易发生形变，致使钢球形变增大，所以接触应力增大。

3.3 滚道半径对滚道接触应力的影响

在钢球直径、外接头与内接头滚道中心距一定的情况下，在9.60～9.72mm范围内均匀选取7组不同大小的滚道半径，研究滚道半径对活动弯头滚道接触应力的影响规律。整理数值仿真数据可得活动弯头滚道最大接触应力随滚道半径的变化规律，如图11所示。

图11表明，随着滚道半径的增大，滚道最大接触应力先减小后增大，在滚道半径为9.66mm时达到最低。这是由于在钢球半径不变情况下，改变滚道半径，相当于改变了钢球与滚道的接触角，导致滚道最大接触应力产生变化。

4 结论

（1）活动弯头各滚道接触应力云图均呈现椭圆形，符合赫兹理论假设；外接头和内接头第一滚道接触应力都明显大于其余滚道，呈现逐渐减小趋势，与现场第一滚道压痕最为严重的现象一致；因管壁内压直接作用于外接头，造成外接头滚道最大接触应力比内接头大。

（2）活动弯头滚道最大接触应力随着滚道轴向间距的增大而增大，随着滚道径向间距的增大大致呈现增长趋势，随着滚道半径的增大先减小后增大，在滚道半径为9.66mm时达到最低。

（3）因为活动弯头各个部分弯头的壁厚不一致，所以其径向刚度有所差别，最终造成滚道接触产生偏载。活动弯头产生偏载的动力原因主要来自管壁内压，活动弯头在轴向载荷的作用下，有利于缓解内压产生的偏载。

参考文献

[1]王洪伟.工程采油技术的现状及展望框架构建［J］.装备维修技术，2020（1）：22．

WANG H W.Present situation and prospect framework construction of engineering oil recovery technology [J]. Equipment Technology，2020（1）：22.

[2]郑华林，张益维，刘少胡.水力压裂冲蚀磨损对连续管剩余寿命影响研究［J］.中国安全生产科学技术，2016，12（7）：110—115．

ZHENG H L，ZHANG Y W， LIU S H. Study on effect of erosion wear to residual life of coiled tubing for hy- draulic fracturing [J].Joumal of Safety Science and Technology，2016，12（7）：110-115.

[3]金雪梅，张祥来，廖浩，等.加砂压裂过程中高压管汇失效爆裂分析［J］．安全，2017，38（1）：17-18.

JIN X M， ZHANG X L，LIAO H，et al.Failure and burst amnalysis of high-pressure tube during sand fractu- ring [J].Safety&Security，2017，38（1）：17-18.

[4]傅建斌，靳彦欣，李勇，等.页岩气压裂高压管汇失效风险分析及防控措施［J］.安全、健康和环境，2021，21（2）：39-44.

FU J B，JIN Y X，LI Y，et al. Failure risk analysis and prevention measures of high-pressure tube in shale gas fracturing construction [J].Safety Health & Envi- ronment，2021，21（2）：39-44.

[5]王森，朱丽云，王振波，等.弯头间连接管长对π形管冲蚀影响的数值模拟［J］．油气储运，2021，40（11）：1285-1292.

WANG S， ZHU L Y， WANG Z B， et al. Numerical simulation on effect of connection pipe length between elbows on erosion of π-shaped pipeline [J].Oil& Gas Storage and Transportation，2021， 40（11）：1285- 1292.

[6]张晓琳，青松铸，周秀兰，等.管道内检测数据坐标化技术［J］．油气储运，2021，40（6）：692—698.

ZHANG X L，QING S Z，ZHOU XL，et al. Data co- ordinate technology of pipeline in-line inspection [J]. Oil &Gas Storage and Transportation，2021，40（6）： 692-698.

[7]徐军，李世勇，陈世波.川南页岩气地面流程测试管线刺漏原因及机理分析［J］．钻采工艺，2021，44（6）：83-87.

XU J， LI S Y， CHEN S B. Cause and mechanism anal- ysis of spillage of shale gas ground flow testing pipelines in Southern Sichuan [J].Drilling& Production Tech- nology，2021，44（6）：83-87.

[8]徐晓东，黄志国.提高高压活动弯头体使用寿命的工艺途径［J］．机械设计与制造，2005（8）：118—119.

XU X D，HUANG Z G. A technological way to the lm- provement of the extension of service life of high pres- sure movable bend [J].Machinery Design &Manufac- ture，2005（8）：118-119.

[9]敬佳佳，唐曦，陈文斌，等.高产天然气井放喷管汇弯管冲蚀磨损特性研究［J］．表面技术，2021，50（12）：329-339，355.

JING JJ，TANG X， CHEN W B，et al. Study on ero- sion characteristics of elbow erosion of manifold for relief pressure of high-yield natural gas well [J].Surface Technology，2021，50 （12）：329-339，355.

[10]郭登明，周靖力，肖信武，等.延长高压活动弯头寿命的工艺研究［J］．石油机械，2017，45（5）：86-88，94.

GUO D M，ZHOU JL，XIAO X W，et al.Study on the technology for extending service life of high pres- sure movable elbow [J].China Petroleum Machinery， 2017，45（5）：86-88，94.

[11]孙林松，王德信，谢能刚.接触问题有限元分析方法综述［J］．水利水电科技进展，2001，21（3）：18-20.

SUN LS，WANG D X，XIE N G. A summary of fi- nite element analysis for contact problems [J]. Ad- vances in Science and Technology of Water Resources， 2001，21（3）：18-20.

[12]KONYUKHOV A， SCHWEIZERHOF K. Geometrical- ly exact theory for contact interactions of 1D manifolds. Algorithmic implementation with various finite element models [J].Computer Methods in Applied Mechanics and Engineering，2012，205/206/207/208：130- 138.

[13]PENG Z，DU C B，ZHAO W H， et al.Dynamic crack face contact and propagation simulation based on the scaled boundary finite element method [J].Com- puter Methods in Applied Mechanics and Engineering， 2021，385：114044.

[14]姜育松，苏超.工程接触问题数值分析方法［J］.水力发电，2010，36（4）：75—78，94．

JIANG Y S， SU C. Numerical analysis methods for solving contact problems in engineering [J].Water Power，2010， 36 （4）： 75-78，94.

[15]安劭侠，简晓辉，马利，等.基于有限元法的活动弯头接触特性分析［J］.组合机床与自动化加工技术，2021（10）：33—37．

AN S X， JIAN X H， MA L，et al. Analysis of con- tact characteristics of the swivel joints based on finite elementmethod [J]. Modular Machine Tool & Auto- matic Manufacturing Technique，2021（10）：33- 37.

作者简介：张国友，高级工程师，生于1970年，2009年毕业于长江大学，获硕士学位，现从事科研及管理工作。地址（434024）湖北省荆州市。电话：（0716）8429057。E-mail：zhanggy.oset@sinopec.com。

收稿日期：2022—08—22

（本文编辑 刘锋）