Processing math: 100%
  • ISSN 0258-2724
  • CN 51-1277/U
  • EI Compendex
  • Scopus 收录
  • 全国中文核心期刊
  • 中国科技论文统计源期刊
  • 中国科学引文数据库来源期刊

基于EWT-KLD的机械密封金刚石涂层磨损声发射降噪

林志斌 高宏力 吴昱东 谭咏文

孙凤, 邢大壮, 周冉, 金俊杰, 徐方超. 考虑能耗的电磁主动悬架LQR控制策略[J]. 西南交通大学学报, 2023, 58(4): 754-760, 798. doi: 10.3969/j.issn.0258-2724.20220815
引用本文: 林志斌, 高宏力, 吴昱东, 谭咏文. 基于EWT-KLD的机械密封金刚石涂层磨损声发射降噪[J]. 西南交通大学学报, 2024, 59(1): 177-184. doi: 10.3969/j.issn.0258-2724.20210599
SUN Feng, XING Dazhuang, ZHOU Ran, JIN Junjie, XU Fangchao. LQR Control Strategy for Electromagnetic Active Suspension Considering Energy Consumption[J]. Journal of Southwest Jiaotong University, 2023, 58(4): 754-760, 798. doi: 10.3969/j.issn.0258-2724.20220815
Citation: LIN Zhibin, GAO Hongli, WU Yudong, TAN Yongwen. Denoising of Acoustic Emission of Diamond-Coated Mechanical Seals Wear Based on Empirical Wavelet Transform and Kullback-Leibler Divergence[J]. Journal of Southwest Jiaotong University, 2024, 59(1): 177-184. doi: 10.3969/j.issn.0258-2724.20210599

基于EWT-KLD的机械密封金刚石涂层磨损声发射降噪

doi: 10.3969/j.issn.0258-2724.20210599
基金项目: 国家自然科学基金(51775452)
详细信息
    作者简介:

    林志斌(1990—),男,工程师,博士,研究方向为智能化状态监测与故障诊断技术,E-mail:freezinglin@163.com

    通讯作者:

    高宏力(1971—),男,教授,研究方向为设备智能化状态监测与故障诊断技术,E-mail:hongli_gao@swjtu.edu.cn

  • 中图分类号: TP277;TH117.1

Denoising of Acoustic Emission of Diamond-Coated Mechanical Seals Wear Based on Empirical Wavelet Transform and Kullback-Leibler Divergence

  • 摘要:

    为了准确获得机械密封金刚石涂层在磨损过程的声发射信号,在分析机械密封设备的噪声特性基础上,提出了基于经验小波变换(EWT)和相对熵(KLD)的声发射降噪方法;通过对磨损声发射信号进行经验小波变换得到划分其频带的滤波器组,对磨损声发射信号和背景噪声发射信号用相同的滤波器组划分频带;计算相应频带2种信号的相对熵,用累计和算法在升序排列的相对熵中找到首个大于3σ的值作为阈值,保留相对熵值大于阈值的频带重构信号,完成降噪. 研究结果表明:本文所提的EWT-KLD方法可以有效抑制不同工况、不同磨损状态的声发射信号的噪声,有效改善了磨损声发射信号的信噪比,尤其是微弱磨损信号的信噪比,提高了密封端面磨损声发射检测的精度和灵敏度;通过与传统降噪方法的对比发现,本文方法能够对不同工况下的密封磨损声发射信号降噪表现出更强的适应性和稳定性,对于及时检测早期密封磨损和准确监测磨损累积变化过程具有重要意义.

     

  • 高速磁悬浮列车因具有速度快、爬坡强、噪声小、平稳性好等优点,成为一种新型轨道交通工具. 在磁浮列车高速运行时,由于车体与轨道不接触,所以列车的受流需采用无接触受流方式. 因此,高速磁悬浮列车的辅助供电是关键技术之一,无接触式直线谐波发电机成为重要研究方向[1-3].

    高速磁悬浮列车用直线谐波发电机主要应用于常导电磁悬浮系统(electromagnetic suspension, EMS)和超导电动悬浮系统(electrodynamic suspension, EDS)[4-7]. 电磁悬浮系统直线发电机的研究主要包括:磁悬浮列车低速运行时,向长定子中注入高频电流和在列车的磁极上加绕一组直线发电机线圈,实现磁悬浮列车在低速运行时车载蓄电池的充足供电[8]. 文献[9]利用有限元法构建发电机数学模型,在发电机输出特性基础上,定义了目标函数,并针对目标函数优化发电机的配置,得到输出特性的最优解. 文献[10]利用有限元仿真求解出直线发电机感应电动势. 文献[11]通过解析计算求得发电机的磁通密度以及感应电动势. 文献[12]通过解析计算,推导出直线发电机电阻、电感等主要参数.

    在超导电动悬浮系统领域,文献[13]针对超导电动悬浮系统直线发电机,利用ANSYS有限元软件建立其三维模型,分析其在多运行工况下感应电动势、输出功率等特性. 文献[14]在前者基础之上,补充了超导线圈杜瓦外壳的涡流对直线发电机输出特性的影响,并对杜瓦外壳结构进行优化,优化后的杜瓦外壳能增大集电功率. 文献[15-16]提出直线发电机的电流控制方法,利用瞬时电流补偿无功功率,实现直线发电机等效电路功率因数为1. 文献[17-18]提出用于EDS系统的直线发电机,不但具备发电功能,而且具备调节电动悬浮列车转向架阻尼的特性,磁悬浮列车的悬浮力随集电线圈感应电流的相位变化而变化,通过改变相位获得电磁力,从而提高磁浮列车转向架阻尼. 文献[19]针对改进的集电线圈的尺寸和布置,提出一种采用零序电流控制的直线发电机系统,向集电线圈注入零序电流,在垂直方向产生电磁力,通过试验验证,发现该系统能有效抑制转向架振动. 文献[20]提出一种无传感器的垂向阻尼器,通过向直线发电机集电线圈中注入零序电流以产生垂直方向电磁力,提升转向架阻尼. 目前,对电动悬浮系统直线发电机的研究局限于有限元仿真以及直线发电机的阻尼特性分析,对超导电动悬浮系统直线谐波发电机发电特性的研究不足.

    综上,本文针对超导电动悬浮系统用直线谐波发电机的发电特性(超导线圈磁感应强度、悬浮线圈电流以及谐波磁场、集电线圈感应电动势)展开研究. 以日本山梨线MLX01型磁浮列车作为工程研究对象,提出超导线圈的磁动势分布模型,利用解析计算的方法推导出集电线圈感应电动势,通过与有限元、实测数据的比较,验证解析模型的有效性,该解析模型求解速度快,在电动悬浮系统直线谐波发电机的发电特性分析上具备较大优势,为电动悬浮系统直线发电机的设计和控制提供理论基础.

    超导电动悬浮系统直线发电机的拓扑结构如图1所示. 直线发电机主要由超导线圈、悬浮线圈、集电线圈组成,v为磁浮列车运行速度. 超导线圈安装在磁悬浮列车转向架两侧,集电线圈固定于超导线圈外侧,悬浮线圈位于轨道侧壁,沿线路铺设,用于悬浮和导向功能. 转向架单侧有4个超导线圈,相邻超导线圈的极性相反,一个超导线圈极距下对应一组悬浮线圈和一组集电线圈. 其中,悬浮线圈3个为一组,集电线圈4个为一组.

    图  1  直线谐波发电机结构与安装位置
    Figure  1.  Linear harmonic generator structure and installation location

    电动悬浮系统直线谐波发电机的线圈排布如图2所示. 超导线圈在x轴方向极距为τ,悬浮线圈极距为τ1. 超导线圈长度为a0,悬浮线圈长度为a1. 一组集电线圈由2个R相、1个S相和1个T相组成. R相、S相和T相集电线圈长度分别为a2Ra2Sa2T. S、T两相极距为4τ/15,R相之间极距为x2R.

    图  2  直线谐波发电机的线圈排布
    Figure  2.  Coil arrangement of linear harmonic generator

    直线谐波发电机的运行原理为:高速磁浮列车运行时,通电的超导线圈掠过悬浮线圈,在悬浮线圈内产生感应电动势,从而产生感应电流和磁场 [6]. 同时,集电线圈与悬浮线圈存在相对位移,集电线圈切割悬浮线圈感应磁场,产生感应电动势.

    利用空间谐波法,在xz方向分别构建超导线圈磁动势周期性模型. 其中:x方向为实际存在的超导线圈,z方向为假设存在无限多的超导线圈,如图3所示. 图中:Fx(x)、Fz(z)分别为超导线圈在xz轴方向的磁动势;Fz0(x) 为超导线圈在中心z = z0的磁动势;b0为超导线圈高度;τzz方向极距;Ns为超导线圈匝数;Is为单匝线圈电流.

    图  3  超导线圈磁动势分布模型
    Figure  3.  Magnetomotive force distribution model of superconducting coil

    由于相邻超导线圈极性相反,其磁动势沿x方向分布如图3(a)所示. 超导线圈中心z = z0,将超导线圈磁动势展开为傅里叶级数,如式(1)所示.

    Fz0(x)=m=1axmcos(mπ τx),
    (1)
    axm=1τττFx(x)cos(mπ τx)dx,
    (2)

    式中:axmx轴方向傅里叶级数系数.

    超导线圈磁动势沿z方向分布如图3(b)所示,其傅里叶级数系数aznn = 1,2,3,…)为

    azn=1τzτzτzFz(z)cos(nπ τzz)dz.
    (3)

    由此,可推出超导线圈磁动势在xOz平面的分布为

    F(x,z)=n=1azncos[nπτz(zz0)] = 16NsIsπ2×m=1n=1mncxmczncos(kxmx)cos[kzn(zz0)],
    (4)

    式中:kxm=mπ /τkzn=nπ /τz cxm=sin(kxma0/2)czn=sin(kznb0/2).

    由边界条件和拉普拉斯方程求出超导线圈的标量磁位ψx,y,z),如式(5)所示.

    {ψ(x,0,z)=F(x,z)/2,ψ(x,,z)=0,2ψ(x,y,z)=0.
    (5)

    利用磁位与磁感应强度的关系式,得到超导线圈在y方向的磁感应强度为

    By=μ0ψ(x,y,z)y,
    (6)

    式中:μ0为真空磁导率.

    求解式(6)得

    By=8μ0NsIsπ2m=1n=1mnλmncxmczneλmn(yy0)×cos(kxmx)cos[kzn(zz0)],
    (7)

    式中:λmn=k2xm+k2zny0为超导线圈中心纵坐标.

    超导线圈、集电线圈和悬浮线圈在yOz平面的位置关系如图4所示. 原点为列车转向架中心. 超导线圈中心位置为(x0y0z0),集电线圈中心坐标为(x2y2,0),悬浮线圈中心坐标为(x1y1,0). z1Uz1B分别为悬浮线圈上线圈中心、下线圈中心的z轴坐标;z2Uz2B分别为集电线圈上线圈中心、下线圈中心的z轴坐标.

    图  4  yOz平面内线圈间位置关系
    Figure  4.  Positional relationship between coils in the yOz plane

    悬浮线圈固定于轨道侧壁,与地面保持静止,而超导线圈固定于车体,并随列车一同移动. 设时间t = 0时超导线圈中心正对的悬浮线圈编号为0号,沿x正方向依次为0,1,2,…,第V个悬浮线圈中心的x轴坐标为

    xV=Vτ1vt
    (8)

    式中:τ1为磁浮线圈的极距.

    根据超导线圈的磁场与悬浮线圈存在交链关系,对超导线圈磁场与集电线圈所在区域进行积分,求解得到超导线圈磁感应强度与悬浮线圈上线圈的交链磁通ϕ1U以及超导线圈磁感应强度与悬浮线圈下线圈的交链磁通ϕ1B,分别如式(9)、(10)所示.

    ϕ1U = x1+a1/2x1a1/2z1U+b1/2z1Ub1/2Bydxdz,
    (9)
    ϕ1B = x1+a1/2x1a1/2z1B+b1/2z1Bb1/2Bydxdz,
    (10)

    式中:b1为悬浮线圈上线圈与下线圈的高度.

    超导线圈磁感应强度与悬浮线圈交链的总磁通量为

    ϕ1=ϕ1Bϕ1U.
    (11)

    悬浮线圈电流的回路方程为

    R1iV+L1di1dt=N1dϕ1dt
    (12)

    式中:R1为悬浮线圈电阻; L1为悬浮线圈电感; N1为悬浮线圈的匝数;iV为第V个悬浮线圈感应电流,如式(13)所示.

    iV=N1g1(ω)m=1m(PBPU)cos(mωtmπV/3)
    (13)
    ω=π τv,
    (14)
    g1(ω)=jωR1+jωL1,
    (15)

    式中:ω为悬浮线圈电流的基波频率;g1(ω)为能够反映悬浮线圈的导纳.

    Pi为悬浮线圈与超导线圈间的位置函数,反映了超导和悬浮线圈间的电磁耦合紧密程度,i∈{B,U},U代表上线圈,B代表下线圈.

    Pi=32μ0NsIsττz×n=1λmnSscSleveλmn(y1y2)cos[kzn(z1iz0)],
    (16)

    式中:SscSlev为系数,见式(17).

    {Ssc=sin(kxma0/2)sin(kznb0/2)kxmkzn,Slev=sin(kxma1/2)sin(kznb1/2)kxmkzn.
    (17)

    悬浮线圈电流的谐波分量随时间的变化曲线如图5所示. 由图可知:悬浮线圈电流主要由基波分量和三次谐波分量构成. 本文只考虑悬浮线圈电流的基波分量以及三次谐波分量产生的磁感应强度. 悬浮线圈电流的基波分量幅值IV1、三次谐波分量幅值IV3分别如式(18)、(19)所示.

    图  5  悬浮线圈电流谐波分析
    Figure  5.  Current harmonic analysis of suspension coil
    IV1=N1g1(ω)(PBPU),
    (18)
    IV3=3N1g1(3ω)(PBPU).
    (19)

    悬浮线圈电流的磁感应强度求解步骤和推导超导线圈磁感应强度的方法类似,此处不再赘述. 求解得到基波电流的磁感应强度B11i和三次谐波电流的磁感应强度B13i分别如式(20)、(21)所示.

    B11i={12N1μ0ττzm=1n=1IV1λmnSlevcos[kzn(zz1i)]eλmn|yy1|cos[kxmx(m1)ωt],m = 1,7,13,12N1μ0ττzm=1n=1IV1λmnSlevcos[kzn(zz1i)]eλmn|yy1|cos[kxmx + (m + 1)ωt],m = 5,11,17,
    (20)
    B13i=24N1μ0ττzm=1n=1IV3λ3m,nSlevcos(3mωt+3kxmx)cos(3ωt)cos[kzn(zz1i)]eλmn|yy1|,m=3915.
    (21)

    悬浮线圈电流产生的谐波磁场如图6所示,分析得知:悬浮线圈电流的感应磁场主要由基波磁场、五次谐波磁场构成. 由于基波磁场与随列车一同运动的集电线圈间不存在相对位移,无法在集电线圈中产生感应电动势. 而五次谐波磁场与集电线圈存在相对位移,五次谐波磁场切割集电线圈,继而产生感应电动势. 由于五次以上的高次谐波磁场含量很小,为简化计算,可忽略五次以上谐波磁场,故本文将五次谐波磁场作为发电磁场.

    图  6  悬浮线圈电流产生的谐波磁场
    Figure  6.  Harmonic magnetic field generated by suspension coil current

    令式(20)中m = 5,则五次谐波磁场的磁感应强度如式(22)所示.

    B15i=12N1μ0ττzm=1n=1IV1λmnSleveλmn|yy1|×cos[kzn(zz1i)]cos[kxmx + (m + 1)ωt].
    (22)

    悬浮线圈电流产生的五次谐波磁场在集电线圈所在区域进行积分,求解得到悬浮线圈上、下线圈与集电线圈上、下线圈的交链磁通,如式(23)、(24)所示.

    ϕiU=x2x2R2+a2R2x2x2R2a2R2z2U+b2U2z2Ub2U2B15idxdz,
    (23)
    ϕiB=x2x2R2+a2R2x2x2R2a2R2z2B+b2B2z2Bb2B2B15idxdz,
    (24)

    式中:ϕUU为集电线圈上线圈与悬浮线圈上线圈磁感应强度的交链磁通;ϕBU为集电线圈上线圈与悬浮线圈下线圈磁感应强度的交链磁通;ϕUB为集电线圈下线圈与悬浮线圈上线圈磁感应强度的交链磁通;ϕBB为集电线圈下线圈与悬浮线圈下线圈磁感应强度的交链磁通;b2B为集电线圈下线圈的高度;b2U为集电线圈上线圈的高度.

    R相集电线圈与悬浮线圈五次谐波磁场交链磁通ϕ2如式(25)所示.

    ϕ2=ϕUUϕUB+ϕBBϕBU=48N1μ0ττzm=1n=1IV1λmnSlevKeλmn|y2y1|cos[kxm(x2x2R/2)+6ωt],
    (25)
    K=SUcos[kzn(z1Uz2U)]SUcos[kzn(z1Bz2U)] +SBcos[kzn(z1Bz2B)]SBcos[kzn(z1Uz2B)],
    (26)

    式中:SBSU为系数,见式(27).

    {SB=sin(kxma2R/2)sin(kznb2B/2)kxmkzn,SU=sin(kxma2R/2)sin(kznb2U/2)kxmkzn.
    (27)

    R相集电线圈感应电动势eR与交链磁通ϕ2的关系如式(28)所示.

    eR=dϕ2dt.
    (28)

    推导得出

    eR=288ωN1μ0ττzm=1n=1IV1λmnSlevKeλmn|y2y1|×sin[6ωt+kxm(x2x2u/2)].
    (29)

    由感应电动势表达式可知:感应电动势幅值与列车运行速度相关,其频率是悬浮线圈电流基波频率的6倍. 其余相集电线圈感应电动势求解与上述求解过程类似,此处不再赘述.

    根据日本山梨线电动悬浮系统直线谐波发电机的基本参数[18,21](如表1所示),通过ANSYS有限元软件建立直线谐波发电机的三维有限元仿真模型(如图7所示),且利用日本山梨线实测数据[19,22]进行对比,进一步验证磁动势分布模型、解析公式的正确性. 图8为试验装置图[23].

    表  1  直线谐波发电机基本参数
    Table  1.  Parameters of linear harmonic generator
    线圈参数数值
    集电线圈a2R/mm200
    a2S/mm245
    a2T/mm245
    b2U/mm435
    b2B/mm245
    与超导线圈间距/mm75
    悬浮线圈a1/mm350
    b1/mm340
    N1/匝24
    τ1/mm450
    与超导线圈间距/mm185
    超导线圈a0/mm1070
    b0/mm500
    Ns/匝1400
    τ/mm1350
    Is/A500
    下载: 导出CSV 
    | 显示表格
    图  7  利用ANSYS 2019 R2建立的三维模型
    Figure  7.  3D model created with ANSYS 2019 R2
    图  8  试验装置
    Figure  8.  Test device

    超导线圈磁感应强度By沿着x轴方向分布如图9所示. 可知:两条曲线均呈马鞍形,解析计算峰值为885.7 mT,峰值的有限元结果为875.8 mT,相对误差为1.1%,验证了磁动势分布模型的有效性.

    图  9  超导线圈磁感应强度
    Figure  9.  Magnetic induction intensity of superconducting coil

    列车速度500 km/h,悬浮线圈电流的解析值与有限元仿真的对比如图10所示. 解析计算的峰值为918 A,有限元结果为888 A,相对误差为3.2%,验证了解析模型的正确性.

    图  10  悬浮线圈电流
    Figure  10.  Suspension coil current

    列车速度400 km/h,集电线圈感应电动势的解析值与有限元仿真的对比如图11所示. 感应电动势的解析值为133.0 V,而感应电动势的有限元仿真值为129.5 V,误差为2.6%,且两者相位几乎相等.

    图  11  感应电动势解析计算与有限元仿真对比
    Figure  11.  Comparison between analytical calculation and finite element simulation of induced electromotive force

    集电线圈感应电动势随速度的变化曲线如图12所示. 可知:感应电动势与列车速度呈线性关系,是由于速度大于100 km/h,悬浮线圈谐波磁场趋于饱和,此时感应电动势几乎只与速度相关. 感应电动势的解析值与测试数据相对误差约为10%,验证了本文推导的集电线圈感应电动势公式的有效性.

    图  12  感应电动势解析计算与试验数据对比
    Figure  12.  Comparison between analytical calculation and experimental data of induced electromotive force

    集电功率的解析值、有限元仿真和测试数据的对比如图13所示. 列车速度500 km/h,解析计算得到输出功率为43.3 kW,有限元结果为41.8 kW,试验数据为39.5 kW. 解析结果与试验数据之间平均误差为8.0%,有限元与试验数据平均误差为5.5%. 列车速度380 km/h,直线谐波发电机达到25.0 kW的目标集电功率,此时能够满足车载供电需求,证明了直线谐波发电机作为高速磁悬浮系统辅助供电的可靠性.

    图  13  集电功率解析计算与有限元、试验对比
    Figure  13.  Comparison of collector powers between analytical calculation, finite element simulation and experimental data

    图14(a)描述了列车速度500 km/h,悬浮线圈电流及其磁感应强度随速度的变化曲线. 由图可知:悬浮线圈电流先随速度增加,在速度约大于100 km/h后进入饱和区段,其数值基本保持恒定的910 A. 而悬浮线圈电流产生的磁感应强度类似,磁感应强度先随速度增大,在速度大于100 km/h后趋于饱和,其数值基本恒定在52.7 mT.

    图  14  悬浮线圈电流及其磁感应强度
    Figure  14.  Suspension coil current and its magnetic induction

    图14(b)描述了列车速度500 km/h,悬浮线圈电流谐波磁场的磁感应强度随速度的变化曲线. 基波、五次、七次谐波磁场的磁感应强度变化趋势类似,均在速度大于100 km/h后进入饱和区段.

    这是由于随着列车速度提升,悬浮线圈中的磁通变化率增大,使得悬浮线圈电流以及磁感应强度随之增大. 然而,列车速度不断提升,悬浮线圈的趋肤效应的影响随之变大,增加了悬浮线圈的有效电阻,使悬浮线圈电流的缓慢变化,同时也使其磁感应强度在高速区段近乎饱和[7].

    图15(a)描述了不同运行速度下,集电线圈感应电动势随时间的变化曲线. 随着列车速度增大,感应电动势和频率随之增大. 这是由于感应电动势频率为6ω,由ω=kx1v得知,感应电动势频率与列车速度成正比关系. 图15(b)描述了列车速度500 km/h,集电线圈各相感应电动势随时间的变化. 图15(c)描述了三相集电线圈感应电动随速度变化关系. 集电线圈的三相电压相位互差120°,R相感应电动势幅值略小于S/T相,是由集电线圈尺寸大小所决定. 三相感应电动势不平衡,连接时会存在零序电流,将零序电流注入集电线圈能够产生电磁力,提升转向架阻尼,进而衰减磁悬浮列车运行时转向架的振动.

    图  15  集电线圈感应电动势
    Figure  15.  Induced electromotive force of collector coil

    以日本山梨线MLX01为研究对象,通过与日本山梨线实测数据和有限元结果的比较验证了模型的有效性. 结论如下:

    1) 基于空间谐波法构建的超导线圈磁动势分布模型,解析求解的超导线圈磁感应强度和有限元结果的曲线高度吻合,验证了磁动势模型的有效性.

    2) 超导电动悬浮系统用直线谐波发电机感应发电主要利用悬浮线圈电流的五次谐波磁场,且感应电动势的频率是悬浮线圈电流基波频率的6倍. 磁浮列车速度大于100 km/h,由于趋肤效应的影响,悬浮线圈电流以及其磁感应强度趋于饱和.

    3) 集电线圈感应电动势与速度呈线性关系,感应电动势大小主要与列车运行速度、超导线圈极距相关和磁动势相关. 随着列车速度增加,速度约为380 km/h时,集电功率达到25.0 kW,此时便达到车载供电需求. 通过将感应电动势、功率的解析值与有限元、试验数据进行对比,验证了本文解析模型的正确性.

  • 图 1  基于经验小波变换和相对熵的机械密封金刚石涂层磨损声发射信号降噪流程

    Figure 1.  Denoising flowchart for AE signal of diamond-coated mechanical seal wear based on EWT and KLD

    图 2  傅里叶频谱的分割

    Figure 2.  Spectrum segmentation

    图 3  机械密封环的摩擦状态监测系统

    Figure 3.  Diagram of friction condition monitoring system for mechanical seals

    图 4  3种工况下实验样本的功率谱汇总

    Figure 4.  Power spectrum samples under three working conditions

    图 5  降噪前后信号原始波形与功率谱

    Figure 5.  Waveform and power spectrum of orignal signal and denoised signal

    图 6  传统EWT相关系数降噪方法

    Figure 6.  Traditional denoising method based on EWT and correlation coefficient

    图 7  本文提出的EWT-KLD降噪方法

    Figure 7.  Proposed denoising method method based on EWT-KLD

    图 8  不同降噪方法得到的去噪AE信号ASL

    Figure 8.  Denoised AE signal ASL via two approaches

    表  1  被测试机械密封环的基本参数及形貌

    Table  1.   Basic parameters and morphology of testing mechanical seals

    密封环 材料 端面内
    径/mm
    端面外
    径/mm
    端面涂
    层/μm
    表面维氏
    硬度/GPa
    外观 动环横截面
    SEM
    (3000 倍)
    金刚石涂层表面
    SEM
    (20000 倍)
    静环 掺杂石墨
    的碳化硅
    52 58 24.5
    动环 碳化硅基
    底金刚石
    涂层
    52 63 6 98
    下载: 导出CSV

    表  2  被测试机械密封环的设计参数

    Table  2.   Design parameters of testing mechanical seals

    设计参数
    表面粗糙度/μm0.2
    接触表面载荷/N80
    接触表面积/mm2518.1
    热传导系数/(W·(m·k)−142
    线速度/(m·s−12.26~5.28
    摩擦系数0.15
    下载: 导出CSV

    表  3  实验工况条件

    Table  3.   Experimental working conditions

    工况转速/(r·min−1载荷/N样本序号
    11780801~100
    21780100101~200
    3280080201~300
    下载: 导出CSV
  • [1] 矫健. 高温油泵机械密封失效分析与对策[J]. 内燃机与配件,2020,321(21): 75-76. doi: 10.3969/j.issn.1674-957X.2020.21.034

    JIAO Jian. Failure analysis and countermeasure of mechanical seal of high temperature oil pump[J]. Internal Combustion Engine and Parts, 2020, 321(21): 75-76. doi: 10.3969/j.issn.1674-957X.2020.21.034
    [2] 郝丽霞. 化工离心泵机械密封失效原因分析[J]. 中国设备工程,2020(3): 143-145.

    HAO Lixia. Failure analysis of chemical centrifugal pump mechanical seal[J]. China Equipment Engineering, 2020(3): 143-145.
    [3] 王超. 高温油泵机械密封失效分析及改进[J]. 石油化工设备,2019,48(6): 62-67.

    WANG Chao. Failure analysis and improvement of mechanical seal of high temperature oil pump[J]. Petro-demical Equipment, 2019, 48(6): 62-67.
    [4] 黄文涛. 机械密封在化工上的应用及泄漏的应对[J]. 化学工程与装备,2018(3): 185-186.

    HUANG Wentao. Application of mechanical seal in chemical industry and response to leakage[J]. Chemical Engineering and Equipment, 2018(3): 185-186.
    [5] 袁明蛟. 石油化工生产常用机械密封形式及失效原因探讨[J]. 中国化工贸易,2017,9(10): 204-206. doi: 10.3969/j.issn.1674-5167.2017.10.182

    YUAN Mingjiao. Discussion on commonly used mechanical seal forms and failure reasons in petrochemical production[J]. China Chemical Trade, 2017, 9(10): 204-206. doi: 10.3969/j.issn.1674-5167.2017.10.182
    [6] 陈志高. 某重水堆核电站主泵机械密封失效原因分析[J]. 设备管理与维修,2019(13): 161-164. doi: 10.16621/j.cnki.issn1001-0599.2019.07.79

    CHEN Zhigao. Failure analysis of mechanical seal of main pump in a heavy water reaction nuclear power plant[J]. Equipment Management and Maintenance, 2019(13): 161-164. doi: 10.16621/j.cnki.issn1001-0599.2019.07.79
    [7] 郭洪涛. 反应器冷却水泵机械密封失效原因及对策[J]. 齐鲁石油化工,2019,47(1): 51-52,56. doi: 10.3969/j.issn.1009-9859.2019.01.014

    GUO Hongtao. Failure causes and countermeasures of mechanical seal in cooling water pump for reactor[J]. Qilu Petrochemical Technology, 2019, 47(1): 51-52,56. doi: 10.3969/j.issn.1009-9859.2019.01.014
    [8] 赵远. 单端面机械密封故障的分析与处理分析[J]. 装备维修技术,2020,176(2): 119.

    ZHAO Yuan. Fault analysis and treatment analysis of single end mechanical seal[J]. Equipment Maintenance Technology, 2020, 176(2): 119.
    [9] 尹源,黄伟峰,刘向锋,等. 机械密封智能化的技术基础和发展趋势[J]. 机械工程学报,2021,57(3): 116-128. doi: 10.3901/JME.2021.03.116

    YIN Yuan, HUANG Weifeng, LIU Xiangfeng, et al. Technical Basis and Development Trend of Mechanical Seal Intelligence[J]. Journal of Mechanical Engineering, 2021, 57(3): 116-128. doi: 10.3901/JME.2021.03.116
    [10] TOWSYFYAN H, GU F S, BALL A D, et al. Modelling acoustic emissions generated by tribological behaviour of mechanical seals for condition monitoring and fault detection[J]. Tribology International, 2018, 125: 46-58. doi: 10.1016/j.triboint.2018.04.021
    [11] 李晓晖. 流体动压型机械密封端面状态监测及健康评估[D]. 成都: 西南交通大学, 2016.
    [12] 张尔卿. 机械密封端面状态监测及寿命预测关键技术研究[D]. 成都: 西南交通大学, 2015.
    [13] FAN Y E, GU F S, BALL A. A review of the condition monitoring of mechanical seals[C]//Proceedings of ASME 7th Biennial Conference on Engineering Systems Design and Analysis. Manchester: [s.n.], 2008: 179-184.
    [14] LI R Y, HE D. Rotational machine health monitoring and fault detection using EMD-based acoustic emission feature quantification[J]. IEEE Transactions on Instrumentation and Measurement, 2012, 61(4): 990-1001. doi: 10.1109/TIM.2011.2179819
    [15] HUANG W F, LIN Y B, LIU Y, et al. Face rub-impact monitoring of a dry gas seal using acoustic emission[J]. Tribology Letters, 2013, 52(2): 253-259. doi: 10.1007/s11249-013-0210-2
    [16] 刘忠,宋嘉城,邹淑云,等. 基于EMD的水轮机空化声发射信号阈值降噪方法[J]. 动力工程学报,2018,38(6): 501-507. doi: 10.3969/j.issn.1674-7607.2018.06.012

    LIU Zhong, SONG Jiacheng, ZOU Shuyun, et al. Threshold denoising method for acoustic emission signals from hydraulic turbine cavitation based on EMD[J]. Journal of Chinese Society of Power Engineering, 2018, 38(6): 501-507. doi: 10.3969/j.issn.1674-7607.2018.06.012
    [17] BACCAR D, SÖFFKER D. Wear detection by means of wavelet-based acoustic emission analysis[J]. Mechanical Systems and Signal Processing, 2015, 60/61: 198-207. doi: 10.1016/j.ymssp.2015.02.012
    [18] BIANCHI D, MAYRHOFER E, GRÖSCHL M, et al. Wavelet packet transform for detection of single events in acoustic emission signals[J]. Mechanical Systems and Signal Processing, 2015, 64/65: 441-451. doi: 10.1016/j.ymssp.2015.04.014
    [19] 关山,彭昶. 刀具磨损声发射信号的混沌特性分析[J]. 农业工程学报,2015,31(11): 60-65. doi: 10.11975/j.issn.1002-6819.2015.11.009

    GUAN Shan, PENG Chang. Chaotic characteristic analysis of tool wear acoustic emission signal[J]. Transactions of the Chinese Society of Agricultural Engineering, 2015, 31(11): 60-65. doi: 10.11975/j.issn.1002-6819.2015.11.009
    [20] 律方成,张瑜,董蒙,等. 基于复合评价指标的金属微粒声发射信号最优小波去噪及其特征提取[J]. 高压电器,2017,53(11): 1-8. doi: 10.13296/j.1001-1609.hva.2017.11.001

    LV Fangcheng, ZHANG Yu, DONG Meng, et al. Optimal wavelet denoising and feature extraction of metal particulate acoustic emission signal based on composite evaluation index[J]. High Voltage Appliances, 2017, 53(11): 1-8. doi: 10.13296/j.1001-1609.hva.2017.11.001
    [21] WEI N S, CHEN Z, GU F S, et al. Wavelet package denoising of acoustic emission signals for lubrication oil monitoring in engine systems[C]//Proceedings of IncoME-V & CEPE Net-2020. [S.l.]: Springer, 2021: 69-80.
    [22] GILLES J. Empirical wavelet transform[J]. IEEE Transactions on Signal Processing, 2013, 61(16): 3999-4010. doi: 10.1109/TSP.2013.2265222
    [23] JIANG W L, LI Z B, ZHANG S, et al. Hydraulic pump fault diagnosis method based on EWT decomposition denoising and deep learning on cloud platform[J]. Shock and Vibration, 2021, 2021(8): 6674351.1-6674351.18.
    [24] 张超,张晨,郭晓苏,等. 基于MCKD-EWT的滚动轴承故障诊断研究[J]. 轴承,2020(5): 43-48.

    ZHANG Chao, ZHANG Chen, GUO Xiaosu, et al. Study on fault diagnosis for rolling bearings based on MCKD-EWT[J]. Bearing, 2020(5): 43-48.
    [25] CHEGINI S N, BAGHERI A, NAJAFI F. Application of a new EWT-based denoising technique in bearing fault diagnosis[J]. Measurement, 2019, 144: 275-297. doi: 10.1016/j.measurement.2019.05.049
    [26] 王友仁,陈伟,孙灿飞,等. 基于能量聚集度经验小波变换的齿轮箱早期微弱故障诊断[J]. 中国机械工程,2017,28(12): 1484-1490. doi: 10.3969/j.issn.1004-132X.2017.12.015

    WANG Youren, CHEN Wei, SUN Canfei, et al. Early weak fault diagnosis of gearboxes based on energy aggregation and EWT[J]. China Mechanical Engineering, 2017, 28(12): 1484-1490. doi: 10.3969/j.issn.1004-132X.2017.12.015
    [27] 张瑞,邓艾东,司晓东,等. 一种新的声发射信号消噪及故障诊断方法[J]. 振动与冲击,2018,37(4): 75-81.

    ZHANG Rui, DENG Aidong, SI Xiaodong, et al. A new method for acoustic emission signal de-noised and fault diagnosis[J]. Journal of Vibration and Shock, 2018, 37(4): 75-81.
    [28] DUBITZKY W, WOLKENHAUER O, CHO K, et al. Encyclopedia of systems biology[M]. New York: Springer, 2013.
  • 加载中
图(8) / 表(3)
计量
  • 文章访问数:  302
  • HTML全文浏览量:  165
  • PDF下载量:  22
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-07-27
  • 修回日期:  2021-10-29
  • 网络出版日期:  2022-11-02
  • 刊出日期:  2021-11-03

目录

/

返回文章
返回