高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

激光超声检测带涂层金属表面裂纹的数值研究

徐志祥 黄建华 王铮恭 黄义敏 王雨

引用本文:
Citation:

激光超声检测带涂层金属表面裂纹的数值研究

    作者简介: 徐志祥(1965-), 男, 教授, 现主要从事测试与控制技术研究。E-mail:zxxu@dlut.edu.cn.
  • 基金项目:

    中央高校基本科研业务费资助项目 DUT15ZD110

  • 中图分类号: TN247

Numerical study on coated metal surface crack by laser ultrasonic detection

  • CLC number: TN247

  • 摘要: 为了研究线源脉冲激光激发的超声波在带涂层金属板表面裂纹检测方面的应用,采用有限元模拟的方法,分别建立了含有裂纹的带镍涂层和不带镍涂层金属板模型,并模拟出激光激发出的瑞利波以及瑞利波的传播过程。通过对接收点处的波形进行理论分析,得出了涂层厚度、裂纹深度与瑞利波时频域信号的关系。结果表明,瑞利波波速随着涂层厚度h的不同而不断变化;当表面存在裂纹时,不带涂层模型的反射瑞利波与剪切瑞利波的到达时间差Δt与裂纹深度hc成线性关系,带涂层模型的Δthc以涂层厚度为分界点成分段线性关系。此研究结果为实际测量带涂层金属板的表面裂纹深度提供了参考。
  • Figure 1.  Model diagram

    Figure 2.  Schematic diagram of finite element grid division

    Figure 3.  Full-field displacement diagram with 0μm thickness coating model at different moments

    a—t=1.53μs b—t=2.58μs c—t=3.63μs d—t=4.52μs

    Figure 4.  Full-field displacement diagram with 100μm thickness coating model at different moments

    a—t=1.53μs b—t=2.58μs c—t=3.63μs d—t=4.52μs

    Figure 5.  The influence of different coating thickness models on time domain signals of the receiving point A

    Figure 6.  The results of short-time Fourier transform of time domain signal of different coating thickness models at the receiving point A

    Figure 7.  Relationship between Rayleigh wave velocity and coating thickness h

    Figure 8.  Rayleigh wave time domain waveform of uncoated 0μm crack and 300μm crack

    Figure 9.  Rayleigh wave frequency domain signal of uncoated 0μm crack and 300μm crack

    Figure 10.  Relationship between Δt and hc at different coating thickness

    表 1  Arrival time of Rayleigh wave at different receiving points and coating thickness h

    coating
    thickness h/μm
    0 50 100 150 200 250 300 350 400 450 500 550
    D/μs 0.652 0.716 0.776 0.800 0.796 0.788 0.776 0.772 0.752 0.746 0.736 0.728
    C/μs 1.344 1.444 1.588 1.636 1.636 1.620 1.608 1.588 1.572 1.548 1.532 1.524
    下载: 导出CSV
  • [1]

    MA J, ZHAO Y, ZHOU F Y. Effect of defocusing amount on thickness measurement based on laser ultrasound[J]. Laser Technology, 2015, 39(3):349-352(in Chinese). 
    [2]

    LIU P P, NAZIRAH A W, SOHN H. Numerical simulation of damage detection using laser-generated ultrasound[J]. Ultrasonics, 2016, 69:248-258. doi: 10.1016/j.ultras.2016.03.013
    [3]

    ZHAN Y, LIU Ch Sh, ZHANG F P, et al. Experimental study and finite element analysis based on equivalent load method for laser ultrasonic measurement of elastic constants[J]. Ultrasonics, 2016, 69:243-247. doi: 10.1016/j.ultras.2016.03.014
    [4]

    GUAN J F. Numerical analysis of inspection of micro-surface crack by laser generated surface acoustic waves[J]. Chinese Journal of Quantum Electronics, 2011, 28(3):362-368(in Chinese). 
    [5]

    SONG Y X, WANG J. The influence of laser parameters and laser ultrasonic detection methods on ultrasonic signals[J]. Infrared and Laser Engineering, 2014, 43(5):1433-1437(in Chinese). 
    [6]

    ZENG W, WANG H T, TIAN G Y, et al. Study on the oscillation effects of laser-induced SAW and near-surface defects[J]. Chinese Journal of Physics, 2015, 64(13):134302(in Chinese).
    [7]

    GUO H L, ZHENG B, LIU H. Numerical simulation and experimental research on interaction of micro-defects and laser ultrasonic signal[J]. Optics and Laser Technology, 2017, 96:58-64. doi: 10.1016/j.optlastec.2017.04.004
    [8]

    XU B Q, SHEN Z H, LU J. Numerical simulation of laser-induced transient temperature field in film-substrate system by finite element method[J]. International Journal of Heat and Mass Transfer, 2003, 46(25):4963-4968. doi: 10.1016/S0017-9310(03)00345-4
    [9]

    WANG M Y, ZHOU Y J, GUO Ch. Numerical simulation of laser ultrasonic detection of surface micro-crack depth[J]. Laser Technology, 2017, 41(2):178-181(in Chinese). 
    [10]

    LIU H, ZHEN B, WANG Zh B, et al.Time dependence of laser Rayleigh wave on detecting surface defect depth[J]. Laser & Infrared, 2017, 47(6):669-673(in Chinese). 
    [11]

    MA J, ZHAO Y, GUO R, et al. Numerical simulation of temperature rise of material surface irradiated by the laser[J]. Laser Technology, 2013, 37(4):455-459(in Chinese). 
    [12]

    XU B Q, SHEN Zh H, NI X W, et al. Thermal and mechanical finite element modeling of laser-generated ultrasound in coating-substrate system[J]. Optics and Laser Technology, 2006, 38(3):138-145. doi: 10.1016/j.optlastec.2004.12.002
    [13]

    WANG J J, SHEN Zh H, NI X W, et al. Numerical simulation of laser-generated surface acoustic waves in the transparent coating on a substrate by the finite element method[J]. Optics and Laser Technology, 2007, 39(1):21-28. doi: 10.1016/j.optlastec.2005.05.015
    [14]

    SHEN Zh H, YUAN L, ZHANG H Ch, et al. Laser ultrasound in solids[M]. Beijing:Post & Telecom Press, 2015:57-68(in Chinese).
    [15]

    XU B Q, SHEN Zh H, NI X W, et al. Time-frequency analysis of laser-generated ultrasonic wave in coating-substrate systems[J]. Laser Technology, 2004, 28(6):609-612(in Chinese). 
  • [1] 王明宇周跃进郭冲 . 激光超声检测表面裂纹深度的数值模拟. 激光技术, 2017, 41(2): 178-181. doi: 10.7510/jgjs.issn.1001-3806.2017.02.006
    [2] 王玉庆王云霞马世榜 . 基于衍射横波的裂纹激光超声检测方法. 激光技术, 2019, 43(4): 546-550. doi: 10.7510/jgjs.issn.1001-3806.2019.04.019
    [3] 周航张斌冯其波崔建英梁晨黄悦朗 . 环形光源激发超声进行缺陷检测的数值研究. 激光技术, 2021, 45(2): 168-173. doi: 10.7510/jgjs.issn.1001-3806.2021.02.007
    [4] 陈春明战宇刘常升方奔 . 激光超声测量T800/QY9511复合材料工程常数的有限元模拟. 激光技术, 2019, 43(1): 88-92. doi: 10.7510/jgjs.issn.1001-3806.2019.01.018
    [5] 许伯强沈中华倪晓武关建飞陆建 . 涂层-基底系统激光激发表面波的时间-频率分析. 激光技术, 2004, 28(6): 609-612.
    [6] 宋艳马世榜张开飞 . 热弹激光超声激励及缺陷检测的有限元分析. 激光技术, 2021, 45(2): 246-251. doi: 10.7510/jgjs.issn.1001-3806.2021.02.021
    [7] 金谏周若阳赵金峰潘永东刘学增桑运龙顾盛 . 激光扫描混凝土裂缝的超声检测. 激光技术, 2019, 43(4): 453-459. doi: 10.7510/jgjs.issn.1001-3806.2019.04.004
    [8] 马琪钧单奇谢婧石一飞沈中华 . 移动扫描激光线源检测圆柱表面缺陷实验研究. 激光技术, 2009, 33(1): 83-86.
    [9] 关建飞沈中华许伯强陆建倪晓武 . 激光超声在缺陷材料中散射波形的声谱分析. 激光技术, 2005, 29(3): 287-290,314.
    [10] 王贵鑫严刚关建飞 . 金属圆管管壁厚度对瑞利波影响的研究. 激光技术, 2014, 38(1): 58-64. doi: 10.7510/jgjs.issn.1001-3806.2014.01.013
    [11] 胡林亭史德民李佩军任成才 . 激光监测系统测量精度的检测方法. 激光技术, 2008, 32(6): 670-672.
    [12] 邓勇马响 . 激光回馈双折射测量系统波片光轴的自动定位. 激光技术, 2019, 43(2): 217-221. doi: 10.7510/jgjs.issn.1001-3806.2019.02.013
    [13] 徐达何凯平熊伟李华高源 . 线激光散斑检测弹幕武器炮口振动测量方法. 激光技术, 2017, 41(6): 876-880. doi: 10.7510/jgjs.issn.1001-3806.2017.06.022
    [14] 刘鹏高立民乐开端 . 应用于铁轨平直度检测的激光测量系统的研究. 激光技术, 2009, 33(6): 575-578. doi: 10.3969/j.issn.1001-3806.2009.06.005
    [15] 何宁谭智诚吕杏利蒋红艳 . 基于语音识别控制的激光超声水下遥感研究. 激光技术, 2017, 41(3): 391-395. doi: 10.7510/jgjs.issn.1001-3806.2017.03.017
    [16] 马健赵扬周凤艳孙继华 . 离焦量对激光超声测厚的影响. 激光技术, 2015, 39(3): 349-352. doi: 10.7510/jgjs.issn.1001-3806.2015.03.015
    [17] 陈笑徐荣青沈中华陆建倪晓武 . 反射激光式超声测厚传感器及其应用. 激光技术, 2003, 27(6): 557-559.
    [18] 许伯强刘洪凯徐桂东徐晨光李俊敏 . 基于应力-位移混合有限元法的激光超声数值模拟. 激光技术, 2014, 38(2): 230-235. doi: 10.7510/jgjs.issn.1001-3806.2014.02.018
    [19] 王贵鑫严刚关建飞 . 铝材管道管壁厚度对激光超声信号影响的探究. 激光技术, 2014, 38(2): 260-265. doi: 10.7510/jgjs.issn.1001-3806.2014.02.024
    [20] 颜鑫应恺宁戴鹭楠谭钧夫沈中华倪辰荫 . 基于物理信息神经网络的激光超声波场研究. 激光技术, 2024, 48(1): 105-113. doi: 10.7510/jgjs.issn.1001-3806.2024.01.017
  • 加载中
图(10) / 表(1)
计量
  • 文章访问数:  6295
  • HTML全文浏览量:  4011
  • PDF下载量:  161
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-12-25
  • 录用日期:  2018-01-23
  • 刊出日期:  2018-11-25

激光超声检测带涂层金属表面裂纹的数值研究

    作者简介: 徐志祥(1965-), 男, 教授, 现主要从事测试与控制技术研究。E-mail:zxxu@dlut.edu.cn
  • 大连理工大学 机械工程学院, 大连 116024
基金项目:  中央高校基本科研业务费资助项目 DUT15ZD110

摘要: 为了研究线源脉冲激光激发的超声波在带涂层金属板表面裂纹检测方面的应用,采用有限元模拟的方法,分别建立了含有裂纹的带镍涂层和不带镍涂层金属板模型,并模拟出激光激发出的瑞利波以及瑞利波的传播过程。通过对接收点处的波形进行理论分析,得出了涂层厚度、裂纹深度与瑞利波时频域信号的关系。结果表明,瑞利波波速随着涂层厚度h的不同而不断变化;当表面存在裂纹时,不带涂层模型的反射瑞利波与剪切瑞利波的到达时间差Δt与裂纹深度hc成线性关系,带涂层模型的Δthc以涂层厚度为分界点成分段线性关系。此研究结果为实际测量带涂层金属板的表面裂纹深度提供了参考。

English Abstract

    • 表面裂纹大量存在于长期处于交变载荷下的金属构件中,裂纹的存在会导致机器致命性的损坏,并且引发安全问题,因此如何快速、便捷、准确地识别并定位出表面裂纹是当前无损检测的一个热点。由于传统的无损检测技术(如磁粉检测、电涡流检测)不适用于恶劣环境,而激光超声检测技术以其时空分辨率高、信号带宽宽、无损坏、非接触等优点,在无损检测领域得到广泛的应用[1-3]。为了研究激光超声的产生、传播过程以及表面裂纹的检测机理,近十几年来,国内外很多学者进行了大量的研究,特别是利用激光激发出的瑞利波检测表面裂纹的特征,并取得了重大的研究成果[4-11]。但是工程领域中, 为了提高金属构件表面的硬度以及耐磨性等性能,通常会在构件表面镀上不同材料及厚度的涂层,这些涂层的存在使得瑞利波的产生机理与传播更为复杂。虽然有学者研究了涂层/基底系统中激光激发的瑞利波[12-15],但是很少有人研究激光激发的瑞利波与带涂层金属表面裂纹相互作用下的规律。

      作者采用数值模拟的方法,通过对模拟得到的瑞利波时域信号进行短时傅里叶变换与傅里叶变换,分别研究了在镍涂层下的涂层厚度、裂纹深度与瑞利波时频域信号的关系, 为通过瑞利波的时域特征定量表征裂纹的深度提供了参考。

    • 为了使计算出的数据更接近真实情况,将脉冲激光的作用等效成热力源的形式。由于脉冲激光的能量符合高斯分布,材料吸收的能量Q(r, z, t)表示为:

      $ Q\left( r, z, t \right)={{I}_{0}}Af\left( r \right)g\left( t \right) $

      (1)

      对于高斯热源,有:

      $ f(r)=\text{exp}(-\frac{{{r}^{2}}}{r_{0}^{2}})~ $

      (2)

      $ g\left( t \right)=\frac{t}{\tau }~\text{exp}(-\frac{t}{\tau }) $

      (3)

      式中, A为材料对激光的吸收系数,I0为激光峰值功率密度,f(r)为脉冲激光的空间分布,g(t)为脉冲激光的时间分布, r是空间半径, zz方向长度, t是时间。计算中取激光波长λ =1064nm,脉冲上升时间τ=10ns,光斑半径r0=300μm,单脉冲能量E0=13.5mJ。

    • 当线源激光沿宽度方向均匀照射到板上时,该问题可以简化为一个2维的平面问题,又由于入射激光束光斑具有轴对称分布的特性,因此所有数值计算都是基于2维轴对称模型。数值计算过程中,取带涂层铝板的长度L=12mm,总厚度ha=10mm(含涂层),裂纹距离激光束作用中心dc=8mm。信号接收点距离缺陷左边沿d的值分别为0mm,2mm,4mm,6mm,对应点简称为ABCD。模型如图 1所示。

      Figure 1.  Model diagram

      为了精确计算出温度场的分布以及模拟激光激发出的超声波在模型内的传播过程,本文中采用变网格方法,在激光作用区域和缺陷区域采用10μm的密网格,在远离激光作用区域采用90μm的松网格,网格如图 2所示。

      Figure 2.  Schematic diagram of finite element grid division

    • 分别建立涂层厚度h为0μm(无涂层)和100μm、裂纹深度hc均为500μm、缺陷距离激光辐射中心dc均为8mm的模型,并进行数值模拟,得到激光照射后不同时刻的全场位移,如图 3图 4所示。从图 3a图 4a可以清晰地区分激光激发出的头波H、纵波P、横波S和瑞利波R,并看到其能量的差别,其中瑞利波能量较大;从图 3b图 4b可以看到界面反射纵波Pr; 图 3c图 4c是透射瑞利波Rt图 3d图 4d是瑞利波与缺陷作用后的反射瑞利波Rr和模式转换得到的剪切波Rs。从图 3图 4可以看出, 带涂层与无涂层模型的瑞利波在相同时刻到达位置不同,即波速不同;由于涂层对能量的衰减,带涂层的模型中所有波的幅值都有所降低;带涂层的瑞利波后面紧跟着更多频率的波,出现了色散;带涂层模型的负向位移被缺陷反射的掠面纵波正向位移抵消了一部分,因此带涂层的Rr波负向位移值较小。以上表明, 该数值模型可以有效模拟出激光激发的声波以及瑞利波与缺陷作用的过程,为进一步分析瑞利波与裂纹的相互作用奠定基础。

      Figure 3.  Full-field displacement diagram with 0μm thickness coating model at different moments

      Figure 4.  Full-field displacement diagram with 100μm thickness coating model at different moments

    • 为研究镍涂层对激光激发出的瑞利波信号的影响,在裂纹深度hc=0μm的前提下取不同厚度的涂层进行数值模拟。当h分别为0μm和100μm、接收点为A时,时域信号如图 5所示。

      Figure 5.  The influence of different coating thickness models on time domain signals of the receiving point A

      图 5中看到,无涂层模型瑞利波到达时间为2.67μs,带100μm涂层模型的瑞利波到达时间为2.95μs,即镍涂层的存在影响了瑞利波的传播速度。这是因为瑞利波的高频成分因波长短透入深度浅,低频成分因波长长透入深度深,因此大部分高频成分在镍涂层内传播,低频成分在基底铝内传播;又由于铝的瑞利波速度比镍的快,因此高频成分相对于低频成分传播得慢,导致在带100μm涂层模型上的瑞利波在3.7μs时刻出现一定的色散,如图 5b所示。对上述信号进行短时傅里叶变换,得到图 6。从图 6a中可以看到, 对于涂层厚度为0μm的模型,无论高频还是低频的分布都集中在一个时间段内;从图 6b中可以看到, 当涂层厚度为100μm时,低频成分比高频成分出现得更早,因此图 6同样说明高频成分传播速度比低频成分慢,存在着色散。

      Figure 6.  The results of short-time Fourier transform of time domain signal of different coating thickness models at the receiving point A

    • 表 1中给出了在不同涂层厚度、接收点CD的瑞利波到达时间。通过表中数据算出对应波速,并取其平均值,采用B样条插值方法拟合出波速和涂层厚度的曲线图,如图 7所示。

      Figure 7.  Relationship between Rayleigh wave velocity and coating thickness h

      表 1  Arrival time of Rayleigh wave at different receiving points and coating thickness h

      coating
      thickness h/μm
      0 50 100 150 200 250 300 350 400 450 500 550
      D/μs 0.652 0.716 0.776 0.800 0.796 0.788 0.776 0.772 0.752 0.746 0.736 0.728
      C/μs 1.344 1.444 1.588 1.636 1.636 1.620 1.608 1.588 1.572 1.548 1.532 1.524

      图 7中可以看出,瑞利波的波速从接近基底铝的瑞利波波速开始,随着涂层厚度的增大,减小到最小值,然后逐渐增大为趋近涂层镍的瑞利波波速。为分析带涂层金属表面裂纹深度hcRrRs到达时间差Δt的关系提供了参考。

    • 图 8中给出了无涂层、信号接收点B处0μm裂纹和300μm裂纹的声表面波时域波形。可以看到,当表面存在裂纹时,除了R波峰外,会多出两个明显的特征波峰Rr波峰和Rs波峰。

      Figure 8.  Rayleigh wave time domain waveform of uncoated 0μm crack and 300μm crack

      图 9为上述时域信号经过快速傅里叶变换得到的频域信号。为解决频域分辨率较低的问题,这里采用了局部的傅里叶变换的方法,加窗位置为1.6μs~5.6μs。观察频域波形发现,有裂纹时的频谱幅度有微量的增强,这是由于缺陷的存在,使传播到缺陷位置处的声表面波发生反射,致使接收到的声表面波的总体能量变大所造成的。

      Figure 9.  Rayleigh wave frequency domain signal of uncoated 0μm crack and 300μm crack

    • 图 10中分别给出了0μm, 200μm, 400μm的涂层厚度模型接收点B处的Rs波与Rr波到达时间差Δt与裂纹深度hc的关系。从图 10中可以看到, 对于0μm涂层厚度模型,Δthc呈现出线性关系;对于200μm和400μm涂层厚度模型,Δthc在整体上是非线性的关系,但是它们都以各自涂层厚度为分界点的分段线性关系。这是因为当缺陷深度hc小于涂层厚度h时,Rs波大部分都在涂层内传播;当hc大于涂层厚度时,Rs波首先在基底铝内传播,然后在镍涂层内传播,两者属于不同的介质,因此导致不同的线性关系。

      Figure 10.  Relationship between Δt and hc at different coating thickness

    • 作者建立了带涂层金属表面裂纹的有限元数值模型,研究了涂层厚度对瑞利波波速的影响,以及瑞利波时域特征与表面裂纹深度之间的关系。

      (1) 当涂层厚度从0μm增加到500μm时,瑞利波波速逐渐减小到一个极小值,再逐渐增大为接近涂层的瑞利波波速。

      (2) 无涂层模型的Rr波与Rs波的到达时间差与裂纹深度成线性关系,有涂层模型的Rr波与Rs波的到达时间差与裂纹深度以涂层厚度为分界点成分段线性关系。

      (3) 涂层的存在降低了激光激发出的波的幅值,且Rr波的负极性被缺陷反射回来的掠面纵波正极性抵消,幅值被极大地衰弱。

      以上结论为激光超声无损检测的实际测量提供了一种参考。

参考文献 (15)

目录

    /

    返回文章
    返回