高级检索

  • ISSN 1006-3080
  • CN 31-1691/TQ
引用本文:
Citation:

基于微观接触力学的旋塞阀密封面疲劳寿命数值研究

    作者简介: 陆俊杰(1990-),男,江苏无锡人,博士生,研究方向为工程摩擦学。E-mail:lujunjie1990@126.com;
    通讯作者: 安琦, anqi@ecust.edu.cn
  • 中图分类号: TH 137.52+3

Numerical Study on the Fatigue Life of Sealing Surfaces of Plug Valve Based on Micro-Contact Mechanics

    Corresponding author: Qi AN, anqi@ecust.edu.cn
  • CLC number: TH 137.52+3

  • 摘要: 以旋塞阀为研究对象,基于旋塞阀紧定螺钉和旋塞阀阀芯的变形协调关系模型,建立接触面受力计算模型。基于微观接触力学,采用自相关函数构建微观粗糙表面,建立了考虑旋塞阀密封面微凸体疲劳损伤及密封接触面疲劳寿命的计算方法。以F-2 ISO-STANDARD型旋塞阀为研究对象,研究了旋塞阀接触应力随预紧圈数变化时旋塞阀疲劳寿命与拧紧圈数的数值计算关系。研究表明:旋塞阀密封面的寿命随螺钉拧紧圈数增加而减小,且在相同拧紧圈数下随旋塞阀密封面表面粗糙度增加而增加;旋塞阀密封面的寿命随旋塞压紧力增加而减小,且在相同旋塞压紧力下随旋塞阀密封面表面粗糙度增加而减小;旋塞阀密封面疲劳寿命随管内流体压力幅值比增加而增加,随管内流体极限压力幅值增加而减小,随管内流体平均压力增加先增加后减小,并得出相关定量变化曲线。
  • 图 1  旋塞阀结构及工作原理图

    Figure 1.  Scheme of the plug valve structure and working principle

    图 3  旋塞阀密封面接触模型示意图

    Figure 3.  Diagram of plug valve sealing surface contact model

    图 2  工作过程中旋塞阀阀芯受力示意图

    Figure 2.  Force diagram of the cock during working

    图 4  旋塞阀密封面微单元接触模型

    Figure 4.  Plug valve sealing surface contact model

    图 5  单个微凸体接触简化模型

    Figure 5.  Deformation of a single pair of asperities contact

    图 6  管道内流体压力波动随时间变化曲线

    Figure 6.  Curve of fluid pressure fluctuation over time in a pipe

    图 7  微凸体发生疲劳失效后密封面微观接触示意图

    Figure 7.  Schematic diagram of micro-contact on sealing surface after fatigue failure of asperities

    图 8  旋塞阀阀芯密封性能计算流程图

    Figure 8.  Calculation flow chart of sealing performance of plug valve

    图 9  旋塞阀密封面泄漏前疲劳寿命与螺钉拧紧圈数之间的关系

    Figure 9.  Relationship between fatigue life before sealing surface leakage of plug valve and number of screw tightening cycles

    图 10  旋塞与阀座疲劳失效前的旋塞阀工作寿命与螺钉拧紧圈数之间的关系

    Figure 10.  Relationship between the working life of plug valve and the number of screw tightening cycles before the fatigue failure of cock and valve seat

    图 11  旋塞阀密封面泄漏前疲劳寿命与对应的旋塞压紧力之间的关系

    Figure 11.  Relationship between the fatigue life of the sealing surface of the plug valve before leakage and the corresponding pressure of the cock

    图 12  旋塞阀密封面疲劳寿命与管内波动流体压力幅值之间的关系

    Figure 12.  Relationship between the fatigue life of the sealing surface of the plug valve and the amplitude of fluid pressure fluctuation in the pipe

    图 13  旋塞阀密封面疲劳寿命与管内波动流体压力幅值比之间的关系

    Figure 13.  Relationship between the fatigue life of the sealing surface of the plug valve and the amplitude ratio of fluid pressure fluctuation in the pipe

    图 14  旋塞阀密封面疲劳寿命与管内波动流体压力幅值之间的关系

    Figure 14.  Relationship between the sealing surface fatigue life of plug valve and the amplitude ratio of fluid pressure fluctuating in the pipe

    表 1  F-2 ISO-STANDARD(1'')旋塞阀的基本参数

    Table 1.  Parameters of F-2 ISO-STANDARD (1'')

    NameModel
    PlugDN 25
    SeatDN 25
    Prestressed bolts A4
    Connecting boltsA4
    下载: 导出CSV
  • [1] 陆培文. 旋塞阀的功能与应用[J]. 化工设备与管道, 2001, 38(2): 52-55. doi: 10.3969/j.issn.1009-3281.2001.02.013
    [2] 陈浩. 方钻杆旋塞阀的失效分析[J]. 天然气工业, 2007(12): 80-82. doi: 10.3321/j.issn:1000-0976.2007.12.024
    [3] YUZAWA S, HASHIZUME T, OKUTSU R, et al. Lifetime of a Control Valve Plug Related With Cavitation Erosion at Differential Pressure of 20 MPA[C]//ASME 2002 Joint US-European Fluids Engineering Division Conference. Montreal, Quebec, Canada : American Society of Mechanical Engineers, 2002: 255-264.
    [4] 王杭州, 秦华伟, 杨灿军, 等. 深海热液采样阀PEEK阀芯疲劳寿命研究[J]. 浙江大学学报(工学版), 2010, 44(1): 14-18. doi: 10.3785/j.issn.1008-973X.2010.01.003
    [5] 李玉坤, 尤泽广, 陈晓红, 等. RMG530减压阀随机振动疲劳分析[J]. 中国石油大学学报(自然科学版), 2016, 40(5): 142-150.
    [6] 高扬. 基于修正古德曼曲线的阀芯疲劳分析[J]. 流体传动与控制, 2015(4): 4-8. doi: 10.3969/j.issn.1672-8904.2015.04.002
    [7] 王海涛, 吴莘馨. 高温气冷堆主蒸汽隔离阀高温蠕变疲劳特性研究[J]. 核动力工程, 2011, 32(S1): 18-21.
    [8] RICHINS W D, SNOW S D, MILLER G K, et al. Effects on fatigue life of gate valves due to higher torque switch settings during operability testing[R]. Lockheed Idaho Technologies Co. Idaho Falls, ID (United States): [s.n.], 1995.
    [9] 陈浩, 宋周成, 王晓萍, 等. 钻柱内防喷系统研究方向探讨[J]. 石油机械, 2007, 35(7): 73-75.
    [10] 金静静, 李玲. 某油田钻杆旋塞阀断裂原因分析[J]. 理化检验(物理分册), 2016, 52(05): 355-359.
    [11] 冯文荣, 曹佳, 朱海燕, 等. 方钻杆旋塞阀流场的三维数值模拟研究[J]. 石油机械, 2012, 40(05): 7-10.
    [12] 蔄靖宇, 李正美, 安琦. 油管接头台阶密封面的接触压力及密封性能[J]. 华东理工大学学报(自然科学版), 2015, 41(3): 417-423. doi: 10.3969/j.issn.1006-3080.2015.03.022
    [13] 蕳靖宇, 黄龙文, 安琦. 油套管螺纹端面密封微观接触性能数值模拟[J]. 华东理工大学学报(自然科学版), 2016, 42(5): 722-729.
    [14] THOMAS T R. Rough Surfaces[M]. 2nd, London :Imperial College Press, 1999.
    [15] YANG S Z, WU Y (Eds.). Engineering application of time series analysis[M]. Wu Han: China University of science and Technology Press, 1990.
    [16] HARRIS T A, KOTALAS M N. Rolling bearing analysis[M]. 5th, Boca Raton: CRC Press, 2006.
    [17] BORODICH F M. The hertz frictional contact between nonlinear elastic anisotopic bodies[J]. International Journal of Solids and Structures, 1993, 30(11): 1513-1526. doi: 10.1016/0020-7683(93)90075-I
    [18] 焦中良, 李志文, 李志勇, 等. 油气管道波动压力的雨流计数与分析[J]. 油气储运, 2011, 30(8): 624-628.
    [19] MINER M A. Cumulative damage in fatigue[J]. Journal of Applied Mechanics, 1945(67): AI59-AI64.
    [20] 李海梅, 宋刚, 刘永志. 金属材料疲劳极限的估算[J]. 郑州大学学报(工学版), 2002(4): 26-29+39.
    [21] MAN J, ZHOU Q, TAO Z, et al. Micro-scale numerical simulation on metal contact seal[J]. Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, 2014, 228(12): 2168-2177. doi: 10.1177/0954406213515644
    [22] 刘俭辉. 基于损伤力学-临界面法的多轴疲劳寿命预估方法研究[D]. 西安: 西北工业大学, 2015.
  • [1] 徐海涛缪宏钢翁涛安琦 . 汽车发动机张紧轮主轴力学分析及其疲劳寿命计算方法. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20191107001
    [2] 顾烨晨杨家鹏安琦 . 空调贯流风扇转子振动对支承轴承受力及寿命影响. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180604001
    [3] 杨家鹏熊文莉安琦 . 滚子直径误差对双列调心轴承疲劳寿命的影响. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180411001
    [4] 粟爽格安琦 . 内外圈装配过盈量对圆柱滚子轴承力学性能的影响. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190121002
    [5] 高源安琦 . 轴承座同心度误差对深沟球轴承-转子系统振动性能的影响. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180415001
    [6] 许叶龙刘迎圆惠虎於晔鸿 . 固液混合过程的数值模拟及实验研究. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180709001
    [7] 罗安王汉奎王建文 . 基于小冲杆试验数据的力学性能的数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180609002
    [8] 许浩杰高磊陆俊杰安琦 . 考虑波箔变形的波箔型气体箔片轴承润滑性能的数值研究. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190108001
    [9] 魏梓轩周家乐 . 基于VAE的编码DNA载体阻断事件聚类分析与研究. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190424001
    [10] 潘斌翁涛杨家鹏安琦 . 汽车发动机皮带系统臂式张紧轮力学性能研究. 华东理工大学学报(自然科学版), doi: 10.14135/j.cmki.1006-3080.20190114006
    [11] 魏江平林家骏陈宁 . 多特征非接触式测谎技术. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190619002
    [12] 沈冰洁黄婕曹银平张福行 . 低温氧化过程中煤的宏观特性与微观结构变化. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20191028001
    [13] 胡贵华隆建杜文莉 . 考虑汽化影响的乙烯裂解炉对流段数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180913002
    [14] 蔡冬莹王跃林段先建胡丹李玉冰胡彦杰 . 预混合高速射流燃烧器内温度场的数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190221004
    [15] 于苗钱锋胡贵华隆建李天越 . 工业级MIP提升管反应器气固两相流动特性的数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180913003
    [16] 贾翔飞钱嘉澍吴幼青陈剑佩 . 鼓泡塔反应器中两相流动CFD-PBM耦合数值模拟. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20191121001
    [17] 姚永康唐俊马杨琴傅杨彭以元郭志前 . 硅罗丹明荧光探针的研究进展. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180711001
    [18] 赵剑沈阳曹旭妮 . 基于铁蛋白的溶栓蛋白纳米粒子的构建及活性分析. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20180601001
    [19] 柏德鸿宗原赵玲 . 扭曲片强化裂解炉管多尺度混合特性分析. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190315003
    [20] 陈兰萍牛玉刚 . 基于多代理的微电网分区分布式最优潮流分析. 华东理工大学学报(自然科学版), doi: 10.14135/j.cnki.1006-3080.20190523004
  • 加载中
图(14)表(1)
计量
  • 文章访问数:  1489
  • HTML全文浏览量:  1539
  • PDF下载量:  5
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-01-03
  • 网络出版日期:  2019-12-28

基于微观接触力学的旋塞阀密封面疲劳寿命数值研究

    作者简介:陆俊杰(1990-),男,江苏无锡人,博士生,研究方向为工程摩擦学。E-mail:lujunjie1990@126.com
    通讯作者: 安琦, anqi@ecust.edu.cn
  • 1. 华东理工大学机械与动力工程学院 上海 200237
  • 2. 安策阀门(太仓)有限公司 江苏太仓 215400

摘要: 以旋塞阀为研究对象,基于旋塞阀紧定螺钉和旋塞阀阀芯的变形协调关系模型,建立接触面受力计算模型。基于微观接触力学,采用自相关函数构建微观粗糙表面,建立了考虑旋塞阀密封面微凸体疲劳损伤及密封接触面疲劳寿命的计算方法。以F-2 ISO-STANDARD型旋塞阀为研究对象,研究了旋塞阀接触应力随预紧圈数变化时旋塞阀疲劳寿命与拧紧圈数的数值计算关系。研究表明:旋塞阀密封面的寿命随螺钉拧紧圈数增加而减小,且在相同拧紧圈数下随旋塞阀密封面表面粗糙度增加而增加;旋塞阀密封面的寿命随旋塞压紧力增加而减小,且在相同旋塞压紧力下随旋塞阀密封面表面粗糙度增加而减小;旋塞阀密封面疲劳寿命随管内流体压力幅值比增加而增加,随管内流体极限压力幅值增加而减小,随管内流体平均压力增加先增加后减小,并得出相关定量变化曲线。

English Abstract

  • 旋塞阀的密封靠阀芯与阀座表面的压紧接触,其密封性能主要由阀芯和阀座接触表面的力学性能决定[1-2]。旋塞阀在正常工作过程中,由于所密封的管道内的流体压力始终处于变化状态,导致密封接触面易于发生疲劳破坏,如何对旋塞阀接触面的疲劳寿命进行计算,从而实现对其疲劳寿命的预测,是目前面临的问题之一。

    Yuzawa等[3]基于阀座泄漏流量不能超过阀可控的最小流量的准则,提出了一种估计阀寿命的方法。王杭州等[4]研究PEEK阀芯在峰值为70 MPa脉动工作压力下的疲劳寿命,运用传统名义应力法,由PEEK材料的应力-疲劳寿命曲线求得PEEK阀芯在工作过程中的应力-疲劳寿命曲线,结合疲劳损伤累积理论估算出阀芯的疲劳寿命。李玉坤等[5]利用有限元分析软件Workbench建立RMG530减压阀阀体有限元模型,根据现场测得的随机振动载荷对减压阀进行随机振动疲劳分析。高扬[6]基于修正古德曼曲线并结合有限元分析实际应用的阀芯,实现在长时间的疲劳试验前就预判出失效点。王海涛等[7]对隔离阀阀体的高温蠕变疲劳特性进行研究,计算了阀体关键部位的蠕变与疲劳损伤程度,并进一步探索阀体寿命对应力类型和水平的敏感性。Richins等[8]根据USNRC Generic Letter 89-10的测试要求在阀门关闭期间,利用阀芯和阀座之间的超调负载和打开阀门所需的阀杆推力,确定了峰值应力范围来估计由于较高的阀杆推力载荷而导致的疲劳使用系数的增加。陈浩等[9]通过有限元方法对方钻杆旋塞阀本体的疲劳强度进行了分析和试验。金静静等[10]通过宏观分析、化学成分分析、力学性能测试、金相分析及微观形貌分析,对某油田发生的一起钻杆旋塞阀横向断裂事件进行了研究,结果表明旋塞阀的断裂机理为疲劳断裂。冯文荣等[11]采用CFD软件对方钻杆旋塞阀流场进行了三维数值模拟,发现旋钮水平面处的作用力存在3 MPa差异,增加了该处的疲劳损坏的可能性。

    综上所述,目前虽然有一些关于阀体密封面疲劳寿命的研究,但这些研究仍然处于较为宏观的层面,鲜有从基于接触面微凸体真实接触的力学行为对旋塞阀接触性能进行的研究,且无法有效地对旋塞阀性能进行定量计算。因此,本文通过采用自相关函数构建真实的旋塞阀三维粗糙接触表面,且基于微凸体接触理论构建旋塞阀预紧及密封面接触的力学模型,在此基础上,也对旋塞阀疲劳失效过程进行了数值模拟。

    • 旋塞阀的结构及工作原理如图1所示,旋塞阀分别由螺钉、阀芯、压力盘、垫片、PTFE垫片、旋塞阀阀盖与阀体联接螺栓、阀座等组成。密封面接触为阀座金属表面和旋塞金属表面接触,其微观的接触是粗糙表面接触,接触面间的压力主要来自螺栓拧紧力以及管路中的流体压力。

      图  1  旋塞阀结构及工作原理图

      Figure 1.  Scheme of the plug valve structure and working principle

      旋塞阀截断关闭时,管道内流体由于外部因素干扰引起压力波动,进而引起旋塞与阀体密封面间接触压力波动,导致了旋塞和阀座密封面微观接触的疲劳积累损伤。随着工作时间增加旋塞与阀座密封面微凸体由于疲劳积累损伤而脱落,接触密封面微凸体脱落个数增加最终导致旋塞与阀座密封面发生泄漏失效。

    • 在安装过程中,旋塞阀螺钉拧紧$N$圈后实现阀芯与阀座的接触面压紧,此时螺钉、压力盘、垫片等发生弹性变形。螺钉拧紧$N$圈后,将有$n$对内外螺纹相互啮合,螺钉与阀盖第$i$圈内、外螺纹接触面沿齿向中点,根据文献[12],建立紧定螺钉拧紧圈数与压紧力变形协调关系,通过数值计算,可以求出螺钉拧紧$N$圈后对旋塞所产生的压紧力${F_N}$

      旋塞处于闭合状态时,进口端压力为${p_{{\rm{in}}}}$,出口端压力为${p_{{\rm{out}}}}$,旋塞左、右侧密封面接触应力分别为$\alpha_1 $$\alpha_r $,这时旋塞和阀座之间的接触压力使其产生密封效果,管内流体压力作用在旋塞上的合力(进出口压力积分得)${F_P}$为:

      式中:$p = {p_{{\rm{in}}}} - {p_{{\rm{out}}}}$r为旋塞曲率半径;α为旋塞相位角(如图3(a)所示);R2R3为管道流体压力z方向上作用界曲率半径。

      图  3  旋塞阀密封面接触模型示意图

      Figure 3.  Diagram of plug valve sealing surface contact model

      旋塞阀密封面接触模型示意图如图3所示,当管内通有流体时,旋塞与阀座接触面上的压力不呈均匀分布状态,由于受力的不平衡,导致阀座沿圆周方向上的弹性接触变形不均匀,表面微观接触面积远小于宏观尺度上接触面积,阀芯宏观尺度上虽然几乎不发生改变,但密封接触面的弹性压缩量的不均匀将会导致阀芯受压发生偏移。将旋塞接触面沿轴向展开,并进行网格化,如图2所示,每个网格上的接触压力都不相等。

      图  2  工作过程中旋塞阀阀芯受力示意图

      Figure 2.  Force diagram of the cock during working

      每个网格的弹性压缩量$\Delta {s_{i,j}}$可以表达为:

      式中:$\Delta {s_x}$为阀芯位移;${R_{vi}}$为接触微单元上阀芯平均半径;${R_{pi}}$为接触微单元上阀座内圈平均半径,${R_{pi}} = {R_{vi}} - $$ \tan 2^\circ \times \Delta s$$\Delta s$为阀芯$z$方向上位移。

      微单元接触面积${A_{i,j}}$为:

      式中:ΔR为旋塞密封面微单元曲率半径的增量。

      微单元接触高度$\Delta {s_i}$的数值计算基于赫兹接触力学,每个接触单元的接触力${F_{i,j}}$,接触力${F_{i,j}}$水平方向合力${\bar F_p}$

      接触力${F_i}$竖直方向合力${\bar F_N}$

      旋塞密封面接触力${F_{i,j}}$水平方向合力${\bar F_p}$与流体压力、摩擦力${F_f}$平衡;接触力${F_{i,j}}$竖直方向合力${\bar F_N}$与预紧螺栓传递到阀芯压力${F_N}$平衡。

    • 旋塞与阀座密封面实际的接触表面是粗糙的表面,在计算每个接触单元的接触力${F_{i,j}}$时,应考虑微观接触的受力变形。为此,本文采用自相关函数[13-15]构建旋塞阀座套筒的接触表面:

      其中,$\sigma $为表面均方根粗糙度;$x$为密封圈接触面圆周方向坐标,$y$为密封圈接触面轴向坐标,${\tau _x}$${\tau _y}$分别为$x$$y$方向任意两点之间的距离;${\beta _x}$${\beta _y}$分别为$x$$y$方向上的相关长度。

      旋塞阀密封面微单元接触模型如图4所示,数值计算微单元的三维随机粗糙表面(Ra=0.5),本文参与接触密性数值计算为抛光后加工的不锈钢钢材材料三维粗糙表面为密封表面研究对,并赋予随机及表面相应的力学性能参数。旋塞阀螺钉拧紧$N$圈后,管内流体压力波动值在${p_i}$时,旋塞与阀座两表面第(i, j)对密封微单元如图4(c)所示叠放接触,该微单元接触应力云图如图4(d)所示。

      图  4  旋塞阀密封面微单元接触模型

      Figure 4.  Plug valve sealing surface contact model

      微观上,两密封表面的接触实际上是两个表面上微凸体之间接触。对于单个微凸体间的接触,可以将接触面近似看作为球面接触,两微凸体的微观接触状态如图5所示。其中:${R_{i1}}$${R_{i2}}$分别为两微凸体接触面曲率半径;${\delta _i}$为单对微凸体的接触高度;${P_i}$为单对微凸体接触力平衡的外载荷;${a_i}$为单对微凸体的接触面积。

      图  5  单个微凸体接触简化模型

      Figure 5.  Deformation of a single pair of asperities contact

      本文假设接触过程中微凸体只发生弹性变形,采用赫兹接触理论来研究接触力、接触变形与外载荷之间关系。根据赫兹接触理论,两球形微凸体发生点接触,接触力${P_i}$、最大接触应力$\sigma _{\max }^i$与微凸体的接触高度${\delta _i}$关系式为[16-17]

      式中:${E_1}$${E_2}$分别为旋塞、阀座弹性模量;${\nu _1}$${\nu _2}$分别为旋塞、阀座泊松比。

      旋塞阀接触单元的接触力${F_{i,j}}$表示为

      结合旋塞阀密封面粗糙度的影响,利用上述模型可计算出在不同过盈量下旋塞阀阀芯与阀座在不同过盈量下密封面接触力分布。根据文献[18]中油气管道雨流计数的波动压力采集数据结果分析,得到旋塞阀管道流体压力波动分布规律为随机分布。按此分布规律使用Matlab内置函数randsample数值模拟管道内压力随时间波动,结果如图6所示,数值模拟结果得到旋塞阀管内流体随时间波动值,管内流体压力极限${P_{{\rm{max}}}}$、流体压力平均值${P_{{\rm{mean}}}}$、流体压力幅值${P_{\rm{a}}}$和管内流体压力比${r_{\rm{P}}}$

      图  6  管道内流体压力波动随时间变化曲线

      Figure 6.  Curve of fluid pressure fluctuation over time in a pipe

      管道内流体压力随机变化引起旋塞阀密封面接触应力变化,将导致接触面的疲劳破坏。根据Palmgren和Miner[19]提出的线性疲劳累积损伤计算方法,将各应力相对应的累积循环次数与各应力相对应的材料发生疲劳破坏时的极限循环次数商的总和定义为材料的总寿命损伤,理论上该值为1,即:

      其中:D为总寿命损伤;${\sigma _i}$为密封圈接触表面的应力幅;Ni为各应力相对应的累积循环次数;${N'_i}$为与各应力相应力相对应材料发生疲劳破坏时的极限循环次数。

      考虑工程应用数据难易度,本文选取丁氏的估算公式[20]

      其中:${\sigma _r}$为应力比为$r$的交变应力中最大应力的极限值;$f$为材料的疲劳比,($f = {{{\sigma _{ - 1}}} / {{\sigma _b}}}$),${\sigma _{ - 1}} $为对称循环下材料的疲劳极限,${\sigma _b}$为材料静抗拉强度极限;$n$为材料常数,与$f$有关,$n = {1 / {\left( {c \cdot f} \right)}}$$c$为待定系数)。

      其中:${\sigma _f}$为应力比为$r$的交变应力中应力幅的极限值;m是疲劳曲线指数;${N'_f}$为达到疲劳极限的循环次数。

      微凸体发生疲劳失效后密封面微观接触示意图如图7所示,如图7(a)中,计算微凸体发生疲劳失效后密封面微观接触为微单元的弹性压缩量$\Delta {s_{i,j}}$,得到该状态下旋塞阀接触单元的接触力${F_{i,j}}$。根据文献[22]中的密封面泄漏识别算法,判定密封面是否泄漏,如果密封面出口通道出现标记值−1,即旋塞阀密封面发生泄露,则判定旋塞阀与阀座密封失效,如图7(b)所示。

      图  7  微凸体发生疲劳失效后密封面微观接触示意图

      Figure 7.  Schematic diagram of micro-contact on sealing surface after fatigue failure of asperities

    • 旋塞阀密封性能数值计算方法步骤如下,流程图如图8所示:

      图  8  旋塞阀阀芯密封性能计算流程图

      Figure 8.  Calculation flow chart of sealing performance of plug valve

      (1)给定旋塞阀预紧圈数$N$,密封面粗糙度${R_{\rm{a}}}$,旋塞阀工作时间${T_i} = {T_{i - 1}} + \Delta T$,进、出口端流体波动压力值为${p_{{\rm{in}}}}$${p_{{\rm{out}}}}$

      (2)计算管内流体压力作用在旋塞上的合力${F_P}$,旋塞所产生的压紧力${F_N}$,微凸体发生点接触最大接触应力$\sigma _{\max }^i$

      (3)判断密封面每一个微凸体是否发生疲劳失效,如果微凸体发生失效,失效微凸体不再参与密封面接触力计算;否则,密封面保持不变;

      (4)根据文献[21]中的密封面泄漏识别算法,判定密封面是否泄漏,如果密封面发生泄露,那么旋塞阀密封面失效;否则,重复计算步骤(1)至步骤(4)。

      (5)输出旋塞阀密封面关闭状态的疲劳寿命T

    • 基于上述构建的计算模型,本文以金属接触密封圆锥形旋塞阀为研究对象,研究了旋塞阀密封面疲劳泄漏前预期寿命$T$随着紧定螺栓预紧紧圈数$N$、旋塞压紧力${F_N}$、密封面粗糙度${R_a}$以及管内入口处流体波动压力${p_{{\rm{in}}}}$等参数的变化规律。本文以安策阀门公司F-2 ISO-STANDARD型1''旋塞阀外形尺寸为数值计算研究对象,具体参数如表1所示。

      NameModel
      PlugDN 25
      SeatDN 25
      Prestressed bolts A4
      Connecting boltsA4

      表 1  F-2 ISO-STANDARD(1'')旋塞阀的基本参数

      Table 1.  Parameters of F-2 ISO-STANDARD (1'')

    • 根据文献[18]中油气管道雨流计数的波动压力采集数据结果分析得到管道流体压力波动分布规律,同时按此分布规律探讨了管道入口压力随时间波动的情况,结果得流体压力平均值为${P_{{\rm{in,\, mean}}}} = 2.5\,{\rm{Mpa}}$,流体压力幅${P_{{\rm{ina}}}} = 2.5\;{\rm{Mpa}}$,流体压力比${r_p} = 0$。旋塞阀密封机构基体疲劳失效寿命即密封面最大接触应力引起的疲劳失效疲劳寿命,根据文献[22]的实验数据,得到304不锈钢疲劳寿命实验常数$m = 8.849\;6$,304不锈钢的疲劳极限($N = {10^7}$)为$187.51\;{\rm{Mpa}}$

      旋塞阀密封面疲劳泄漏前预期寿命随旋塞压紧螺钉拧紧圈数变化规律如图9所示。旋塞阀密封面在相同粗糙度下,旋塞阀密封面疲劳泄漏前的预期疲劳寿命随着旋塞阀紧定螺钉拧紧圈数减小将会增加,这是由于拧紧圈数减少,其接触面的接触应力也会减少所致。在相同螺钉拧紧圈数下,接触面的粗糙度越大,其所预期疲劳寿命也会越大。研究发现旋塞阀密封面粗糙度${R_{\rm{a}}} = 1.1$时紧定螺钉拧紧圈数小于0.11时,旋塞阀密封面泄漏通道已经形成,导致旋塞阀密封面预期疲劳寿命最大值为$6.382 \times {10^6}\;{\rm{h}}$,将不会出现密封面无限疲劳寿命区。

      图  9  旋塞阀密封面泄漏前疲劳寿命与螺钉拧紧圈数之间的关系

      Figure 9.  Relationship between fatigue life before sealing surface leakage of plug valve and number of screw tightening cycles

      为了对图9的计算结果进行验证,在不考虑微观粗糙表面接触泄漏,仅考虑最大接触应力旋塞与阀座疲劳失效前的旋塞阀工作寿命与螺钉拧紧圈数之间的关系(如图10所示)。结果发现,本文建立旋塞阀密封面疲劳泄漏前预期寿命对应紧定螺钉拧紧圈数的变化规律与考虑最大接触应力旋塞与阀座疲劳失效前的旋塞阀工作寿命对应紧定螺钉拧紧圈数变的化规律基本一致。在数值上,在相同旋塞阀密封面粗糙,相同螺钉拧紧圈数下,根据本文建立基于微观接触微凸体疲劳泄漏模型数值计算的旋塞阀密封面泄漏前的疲劳寿命大于仅考虑最大接触应力旋塞与阀座疲劳失效前的旋塞阀工作寿命。表明本文建立旋塞阀密封面疲劳寿命算法是可靠的。

      图  10  旋塞与阀座疲劳失效前的旋塞阀工作寿命与螺钉拧紧圈数之间的关系

      Figure 10.  Relationship between the working life of plug valve and the number of screw tightening cycles before the fatigue failure of cock and valve seat

      旋塞阀数值计算旋塞阀密封面疲劳泄漏前预期寿命随旋塞压紧力变化规律如图11所示。旋塞阀密封面在相同粗糙度下,随着旋塞压紧力减小,旋塞阀密封面疲劳泄漏前的预期寿命将会增加。在相同旋塞压紧力下,接触面的粗糙度越大,其所预期疲劳寿命越小,这是由于相同旋塞压紧力下,粗糙度高的接触微观表面实际接触面积较小。研究发现旋塞阀密封面粗糙度${{R}_{\rm{a}}} = 1.1$旋塞压紧力小于$3.25 \times {10^4}\;{\rm{N}}$时,旋塞阀密封面泄漏通道已经形成,导致旋塞压紧力继续减小密封面预期疲劳泄漏寿命都为0,将不会出现密封面无限疲劳寿命区。

      图  11  旋塞阀密封面泄漏前疲劳寿命与对应的旋塞压紧力之间的关系

      Figure 11.  Relationship between the fatigue life of the sealing surface of the plug valve before leakage and the corresponding pressure of the cock

      据本文建立基于接触表面微观疲劳泄漏模型数值计算研究管内流体压力波动幅值对旋塞阀密封面疲劳寿命影响,当管内流体平均压力值为${P_{{\rm{inmean}}}} = 5\;{\rm{MPa}}$、螺钉拧紧圈数为$N = 0.05$时,结果如图12所示,旋塞阀密封面疲劳寿命随管内流体压力波动幅值增加而减小,增速随管内流体压力波动幅值增加而增加。这是由于管内平均压力不变时,管内流体压力波动幅增加,其密封面微凸体波动接触应力幅增加且微凸体接触应力比减小所致。

      图  12  旋塞阀密封面疲劳寿命与管内波动流体压力幅值之间的关系

      Figure 12.  Relationship between the fatigue life of the sealing surface of the plug valve and the amplitude of fluid pressure fluctuation in the pipe

      据本文建立基于接触表面微观疲劳泄漏模型数值计算研究管内流体波动压力比对旋塞阀密封面疲劳寿命影响,当管内流体压力极限值为${P_{{\rm{inmax}}}} = 10\;{\rm{MPa}}$,螺钉拧紧圈数为$N = 0.05$时,结果如图13所示,旋塞阀密封面疲劳寿命随管内流体压力波动压力比增加而增加,增速随管内流体波动压力比增加而放缓。这是由于管内流体压力波动比增加,管内流体压力极限值不变时,其密封面微凸体接触应力比增加且接触应力幅值减小所致。

      图  13  旋塞阀密封面疲劳寿命与管内波动流体压力幅值比之间的关系

      Figure 13.  Relationship between the fatigue life of the sealing surface of the plug valve and the amplitude ratio of fluid pressure fluctuation in the pipe

      据本文建立基于接触表面微观疲劳泄漏模型进一步数值计算研究管内波动流体平均压力对旋塞阀密封面疲劳寿命影响,当管内流体压力幅值为Pina=1 MPa、螺钉拧紧圈数为N=0.05时,结果如图14所示,旋塞阀密封面疲劳寿命随管内波动流体平均压力增加,先是增加且增速放缓,达到极值后随管内波动流体平均压力增加而减小,且减速增加。这是由于管内流体波动压力幅不变时,管内流体平均压力增加,其密封面微凸体接触平均应力增加且接触应力波动幅值减小,管内流体在相对低压处接触面微凸体疲劳寿命受接触应力幅值减小影响较大,而管内流体在在相对高压处接触面微凸体疲劳寿命受平均接触应力增加影响较大所致。

      图  14  旋塞阀密封面疲劳寿命与管内波动流体压力幅值之间的关系

      Figure 14.  Relationship between the sealing surface fatigue life of plug valve and the amplitude ratio of fluid pressure fluctuating in the pipe

    • (1)以旋塞阀密封性能为研究对象,对旋塞阀的结构进行分析,考虑流体压力波动作用,建立了旋塞阀紧定螺钉和旋塞阀阀芯的变形协调关系模型;同时采用自相关函数构建微观粗糙表面,对旋塞阀密封面接触受力分析,建立旋塞阀密封面泄漏识别算法;通过数值计算方法对旋塞阀密封面微凸体线性疲劳损伤进行求解,实现了对旋塞阀阀芯与阀座接触面密封疲劳寿命的计算。

      (2)以安策阀门公司F-2 ISO-STANDARD型$1''$旋塞阀作为数值计算研究对象,应用本文所建立的旋塞阀阀芯与阀座金属接触力学数值计算模型,对旋塞阀密封面疲劳泄漏前预期寿命随着紧定螺栓拧紧圈数、旋塞压紧力、密封面粗糙度变化规律进行数值计算,研究表明:随着螺钉拧紧圈数、旋塞压紧力增加旋塞阀密封面疲劳泄漏前预期寿命将减小;在相同螺钉拧紧圈数下,旋塞阀密封面疲劳泄漏前预期寿命随粗糙度增加而增加;在相同旋塞压紧力下,旋塞阀密封面疲劳泄漏前预期疲劳寿命随粗糙度增加而减小。应用本文建立基于微观接触微凸体线性疲劳积累损伤引起密封面泄漏模型的旋塞阀密封面疲劳寿命数值计算结果,由于考虑密封面微观泄漏失效,不会出现无限寿命区。

      (3)据本文建立的基于微观接触微凸体疲劳泄漏模型数值计算研究管内流体压力波动值对旋塞阀密封面疲劳寿命影响,结果表明:若管内流体平均压力值相等,旋塞阀密封面疲劳寿命随管内流体压力幅值增加而减小;若管内流体压力极限值保持不变,旋塞阀密封疲劳寿命随管内流体压力比增大而增大;若管内流体压力幅值保持不变,旋塞阀密封疲劳寿命随管内流体平均压力增加先是增加,随着平均压力增加而减小。

(14)  表(1) 参考文献 (22) 相关文章 (20)

目录

    /

    返回文章