高级检索

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

留言板

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

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

激光超声检测表面裂纹深度的数值模拟

王明宇 周跃进 郭冲

引用本文:
Citation:

激光超声检测表面裂纹深度的数值模拟

    作者简介: 王明宇(1993-), 男, 硕士研究生, 现主要从事锚杆支护无损检测的研究.
    通讯作者: 周跃进, yjin_zhou@163.com
  • 基金项目:

    江苏省普通高校研究生科研创新计划资助项目 KYLX15_1407

    国家自然科学基金资助项目 51374201

  • 中图分类号: O426;TN247

Numerical simulation of laser ultrasonic detection of surface micro-crack depth

    Corresponding author: ZHOU Yuejin, yjin_zhou@163.com ;
  • CLC number: O426;TN247

  • 摘要: 为了研究激光冲击所引起的瑞利表面波在材料表面裂纹检测方面的作用,采用ABAQUS建立了激光超声无损检测表面裂纹的模型,并模拟激光在材料表面激发瑞利波及瑞利波沿材料表面传播的过程。从激光的峰值压力、空间分布、时间分布3个方面将激光冲击作用等效成力的作用,体现了激光参量在作用过程中所起到的作用;同时,使用无限单元消除了边界的反射回波。通过对接收点处瑞利表面波信号的分析,进行了不同深度表面裂纹对瑞利波回波峰值影响的定量比较。结果表明,瑞利反射波的正向峰值与下一次最大正向位移所对应的时间的比值和表面裂纹深度之间存在着线性关系,进而为实验测量表面裂纹深度提供了一种测量方法。
  • Figure 1.  Schematic diagram of model

    Figure 2.  Displacement-time curves of nodes 1, 2, 3

    Figure 3.  Displacement-time curves of nodes 4, 5, 6

    Figure 4.  Rayleigh wave signals with different surface crack depths

    Figure 5.  Peak values variation of Rayleigh wave signals with different surface crack depth

    Figure 6.  Displacement-time curve of node

    Figure 7.  Time ratio Δ versus surface crack depth

  • [1]

    YU Q, HE C X. Application of nondestructive testing in composite materials[J]. Engineering and Test, 2009, 49(2):24-29(in Chin-ese). 
    [2]

    HOU R F, WANG W H, MO R Y, et al. Laser ultrasound detection of surface crack in porcelain insulators[J]. Laser Technology, 2014, 38(1):35-38(in Chinese). 
    [3]

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

    LIU J S, XU Z H, GU G Q. Numerical study on improvement of signal-to-noise ratio of surface acoustic waves based on laser array[J]. Laser Technology, 2011, 35(3):403-406(in Chinese). 
    [5]

    XU B Q, LIU H K, XU G D, et al. Mixed stress-displacement finite element method for laser-generated ultrasound[J]. Laser Technology, 2014, 38(2):230-235(in Chinese). 
    [6]

    HUANG Y, JIANG Y F, JIN H, et al. Propagation of shock wave induced by ring laser and its effect on spalling[J]. Laser Technology, 2013, 37(3):301-305(in Chinese). 
    [7]

    ZHANG Z D, LU C, WEI Y F, et al. Numerical simulation to size surface-breaking tilt cracks by ultrasonic spectral method[J]. Journal of Nanchang Hangkong University(Natural Science Edition), 2009, 23(1):87-92(in Chinese). 
    [8]

    WANG G X, YAN G, GUAN J F. Effect of thickness of metal hollow cylinders on laser thermo-elastic generated Rayleigh wave[J]. Laser Technology, 2014, 38(1):58-64(in Chinese). 
    [9]

    SHEN Z H, XU B Q, NI X W, et al. Numerical simulation of pulsed laser induced ultrasound in monolayer and double layer materials[J].Chinese Journal of Lasers, 2004, 31(10):1275-1280(in Chin-ese). 
    [10]

    ZHU W H, LI Z Y, ZHOU G Q, et al. Measurement of laser-induced shock waves by PVDF[J]. Journal of Experimental Mecha-nics, 1997, 12(2):216-220(in Chinese). 
    [11]

    XU B Q, NI X W, SHEN Z H, et al. Numerical simulation of laser-generated ultrasonic by finite element method in the plate material[J]. Chinese Journal of Lasers, 2004, 31(5):621-625(in Chin-ese). 
    [12]

    GUAN J F, SHEN Z H, NI X W, et al. Numerical study on depth evaluation of micro-surface crack by laser generated ultrasonic waves[J]. Journal of Test and Measurement Technology, 2010, 24(1):15-21(in Chinese). 
    [13]

    DEVAUX D, FABBRO R, TOLLIER L, et al. Generation of shock waves by laser induced plasma in confined geometry[J]. Journal of Applied Physics, 1993, 74(4):2268-2273. doi: 10.1063/1.354710
    [14]

    WALED H, WILLIAM V. Finite elements analysis of Rayleigh wave interaction with finite-size, surface-breaking cracks[J]. Ultraso-nics, 2003, 41(1):41-52. 
    [15]

    GUAN J F, SHEN Z H, XU B Q, et al. Numerical simulation of scattering of Rayleigh wave by a surface——breaking crack[J]. Applied Acoustics, 2006, 25(3):138-144(in Chinese). 
  • [1] 徐志祥黄建华王铮恭黄义敏王雨 . 激光超声检测带涂层金属表面裂纹的数值研究. 激光技术, 2018, 42(6): 801-805. doi: 10.7510/jgjs.issn.1001-3806.2018.06.014
    [2] 周航张斌冯其波崔建英梁晨黄悦朗 . 环形光源激发超声进行缺陷检测的数值研究. 激光技术, 2021, 45(2): 168-173. doi: 10.7510/jgjs.issn.1001-3806.2021.02.007
    [3] 刘津升徐振华顾国庆 . 激光阵列用于改善声表面波信噪比的数值研究. 激光技术, 2011, 35(3): 403-406. doi: 10.3969/j.issn.1001-3806.2011.03.030
    [4] 马琪钧单奇谢婧石一飞沈中华 . 移动扫描激光线源检测圆柱表面缺陷实验研究. 激光技术, 2009, 33(1): 83-86.
    [5] 王贵鑫严刚关建飞 . 金属圆管管壁厚度对瑞利波影响的研究. 激光技术, 2014, 38(1): 58-64. doi: 10.7510/jgjs.issn.1001-3806.2014.01.013
    [6] 陈碧芳刘天夫 . 基于SPIDER法测量飞秒脉冲的数值计算与模拟. 激光技术, 2006, 30(1): 13-15.
    [7] 代保江陈烽张东石杜广庆孟祥卫 . 飞秒激光制备波导型光合波器的数值模拟. 激光技术, 2012, 36(2): 251-254,264. doi: 10.3969/j.issn.1001-3806.2012.02.029
    [8] 赵毕强魏旭立杨振刚刘劲松王可嘉 . 隔热板的太赫兹无损检测. 激光技术, 2015, 39(2): 185-189. doi: 10.7510/jgjs.issn.1001-3806.2015.02.008
    [9] 周宇王天一李文军游承武杨振刚王可嘉刘劲松 . 太赫兹无损检测可视化系统的设计. 激光技术, 2017, 41(3): 351-355. doi: 10.7510/jgjs.issn.1001-3806.2017.03.010
    [10] 严刚徐晓东陆建沈中华倪晓武 . 激光激发源对声表面波的影响. 激光技术, 2006, 30(3): 317-319.
    [11] 谢蒙萌陶宝祺姜锦虎 . 基于光栅大错位数字散斑干涉的无损检测技术. 激光技术, 2001, 25(4): 279-281.
    [12] 孟庆龙张艳尚静 . 光纤光谱结合模式识别无损检测苹果表面疤痕. 激光技术, 2019, 43(5): 676-680. doi: 10.7510/jgjs.issn.1001-3806.2019.05.017
    [13] 赵扬杨平华王铭振曹逸飞 . 钛合金增材制造孔隙缺陷的无损检测研究进展. 激光技术, 2024, 48(3): 432-437. doi: 10.7510/jgjs.issn.1001-3806.2024.03.020
    [14] 梅丽芳秦建红严东兵 . 活性激光焊接304不锈钢温度场的数值与试验研究. 激光技术, 2020, 44(4): 492-496. doi: 10.7510/jgjs.issn.1001-3806.2020.04.016
    [15] 汝连志刘凤德刘双宇张宏白頔 . 保护气体对接头形貌及熔滴过渡的影响与模拟. 激光技术, 2018, 42(3): 390-394. doi: 10.7510/jgjs.issn.1001-3806.2018.03.019
    [16] 汪靓杨宇黄敏朱启兵 . 基于偏振成像技术的油桃机械损伤检测. 激光技术, 2022, 46(6): 841-849. doi: 10.7510/jgjs.issn.1001-3806.2022.06.021
    [17] 郑彩英郭中华金灵 . 高光谱成像技术检测冷却羊肉表面细菌总数. 激光技术, 2015, 39(2): 284-288. doi: 10.7510/jgjs.issn.1001-3806.2015.02.029
    [18] 王聪毅高向东马女杰张艳喜游德勇 . 激光焊接缺陷多向磁场激励下磁光成像检测. 激光技术, 2020, 44(5): 592-599. doi: 10.7510/jgjs.issn.1001-3806.2020.05.011
    [19] 刘子烨刘建军洪治 . 基于光子混频连续太赫兹波的厚度检测. 激光技术, 2016, 40(4): 496-499. doi: 10.7510/jgjs.issn.1001-3806.2016.04.009
    [20] 胡林亭史德民李佩军任成才 . 激光监测系统测量精度的检测方法. 激光技术, 2008, 32(6): 670-672.
  • 加载中
图(7)
计量
  • 文章访问数:  9339
  • HTML全文浏览量:  7069
  • PDF下载量:  478
  • 被引次数: 0
出版历程
  • 收稿日期:  2016-01-13
  • 录用日期:  2016-03-21
  • 刊出日期:  2017-03-25

激光超声检测表面裂纹深度的数值模拟

    通讯作者: 周跃进, yjin_zhou@163.com
    作者简介: 王明宇(1993-), 男, 硕士研究生, 现主要从事锚杆支护无损检测的研究
  • 1. 中国矿业大学 力学与建筑工程学院, 徐州 221116
  • 2. 中国矿业大学 深部岩土力学与地下工程国家重点实验室, 徐州 221116
基金项目:  江苏省普通高校研究生科研创新计划资助项目 KYLX15_1407国家自然科学基金资助项目 51374201

摘要: 为了研究激光冲击所引起的瑞利表面波在材料表面裂纹检测方面的作用,采用ABAQUS建立了激光超声无损检测表面裂纹的模型,并模拟激光在材料表面激发瑞利波及瑞利波沿材料表面传播的过程。从激光的峰值压力、空间分布、时间分布3个方面将激光冲击作用等效成力的作用,体现了激光参量在作用过程中所起到的作用;同时,使用无限单元消除了边界的反射回波。通过对接收点处瑞利表面波信号的分析,进行了不同深度表面裂纹对瑞利波回波峰值影响的定量比较。结果表明,瑞利反射波的正向峰值与下一次最大正向位移所对应的时间的比值和表面裂纹深度之间存在着线性关系,进而为实验测量表面裂纹深度提供了一种测量方法。

English Abstract

    • 表面裂纹大量地存在于工程构件当中,如何快速准确地识别它们已成为无损检测的任务。普通的检测方法不适用于高温、高压、有毒[1]等恶劣环境,而激光超声无损检测凭借其非破坏、非接触[2-3]等优点克服了这种局限性,并且具有较高的精度,成为了无损检测的一个重要研究方向[4]。近些年来,随着计算机技术的飞速发展,计算机数值仿真有了越来越广泛的应用。关于激光超声无损检测的数值模拟已成为无损检测领域的一个热点[5-6],而瑞利波进行无损检测的应用也越来越广[7-8]。通过模拟激光激发超声波的物理过程与超声波的传播过程,可以实现对表面裂纹的检测,从而对实验研究提供一种判定方法。

      本文中将激光冲击作用进行简化,将激光的简化模型与激光的参量联系起来,利用ABAQUS软件模拟了激光激发出超声波的物理过程和超声波在材料表面进行传播的过程,并借此模型建立了表面裂纹深度与表面波信号特征量之间的联系,为激光超声检测表面裂纹深度提供了一种有效的方法。

    • 建立一块各向同性且均匀的铝制薄板,铝板的长度为100mm,宽度为30mm,厚度为2mm,裂纹左边界距离板左端75.93mm。现有一束脉冲线源激光沿厚度方向照射到试样表面上,该问题可以简化成平面应变问题来分析, 如图 1所示。

      Figure 1.  Schematic diagram of model

      为了精确模拟激光激发超声波的物理过程和超声波在表面的传播过程,网格的尺寸应不超过波长的1/4[9]。又考虑到计算效率,将网格尺寸取为0.040mm×0.013mm。在模型的左边界添加全约束。有限单元取平面应变单元CPE4R。用有限元来模拟超声波,模型必须要有足够大的几何尺寸才能避免来自边界的反射波对计算的影响,但这样会消耗大量的计算资源和计算时间。为了提高计算效率,同时也为了消除边界回波对表面波信号的干扰,在模型的左边界、右边界以及下边界使用无限单元,无限单元为CINPE4。超声波若在物体中能够完整的传播,需要有足够的时间,因此总时间长度取为5×10-5s,同时保证激光加载的精度,取时间步长为1×10-9s。

    • 为了使激光冲击作用的简化更接近于实际情况,体现激光在激发过程中的作用,需要将激光冲击作用与力作用相互联系。本文中考虑了激光的主要参量,即脉冲宽度、波长、能量、功率密度,从峰值功率密度、时间分布、空间分布3个角度将激光作用简化成力的作用。

      激光波长λ=1064nm,脉冲宽度τ=10ns,激光作用面宽度d=0.6mm,单次脉冲能量E0=13.5mJ,峰值功率密度${I_0} = \frac{{{E_0}}}{{\tau A}} = 7.5 \times {10^6}{\text{W/c}}{{\text{m}}^{\text{2}}} $。其中,A为面积,线源脉冲激光在作用面上产生的峰值压力为[10]

      $ {p_{\max }} = k{\tau ^{ - \frac{1}{8}}}{\lambda ^{ - \frac{1}{4}}}{I_0}^{\frac{3}{4}} = 2.43{\text{MPa}} $

      (1)

      式中,k为与材料有关的量,取k=2.3×10-7

      激光脉冲的空间分布函数[11]为:

      $ \mathit{f}\left( r \right) = \exp \left( { - \frac{{{r^2}}}{{{r_0}^2}}} \right) $

      (2)

      式中, r0是激光辐照的半径, r为距激光作用面中心线的距离。

      激光脉冲的时间分布函数[12]为:

      $ g\left( t \right) = \alpha \frac{{8{t^3}}}{{{\tau ^4}}}\exp \left( { - \frac{{ - 2{t^2}}}{{{\tau ^2}}}} \right) $

      (3)

      式中,α为可由峰值压力确定的修正系数,取α=21,且载荷作用时间t取为脉冲宽度的3倍[13], 即载荷作用时间为30ns。

    • 由于之后的研究是建立在瑞利表面波基础上的,所以有必要验证在模型上加载的脉冲是否能够产生瑞利表面波。一般瑞利表面波的验证是根据其性质来判断的。

      首先在模型表面上离激励处较远的区域取3个点,得到如图 2所示的位移时间曲线。

      Figure 2.  Displacement-time curves of nodes 1, 2, 3

      图 2中1, 2, 3这3条曲线分别代表了离激励处越来越远的3个接收点的瑞利波信号。从图 2可以看出,在离激励作用位置较远的区域上,各节点的相对幅值的差距变小了,波形在传播的过程中将基本不变,此时产生了瑞利波,此点符合瑞利表面波的性质[14]

      然后,在模型不同深度的位置取3个点,得到如图 3所示的位移时间曲线。

      Figure 3.  Displacement-time curves of nodes 4, 5, 6

      图 3中4, 5, 6这3条曲线分别代表了深度越来越大的3个接收点的表面波信号。从图 3可以看出,随着深度的增加,节点振动的位移幅值逐渐衰减,此点符合瑞利表面波的衰减效应。

      由此看出,简化后的激光冲击作用的确可以产生瑞利表面波。

    • 在表面上距裂纹左边界10mm处取信号的接收点,同时裂纹深度分别为0μm,230μm,360μm时,接收点处瑞利波信号的位移时间曲线如图 4所示。

      Figure 4.  Rayleigh wave signals with different surface crack depths

      图 4得知,当材料表面没有裂纹时,接收点处信号并没有出现瑞利反射回波。当材料表面存在裂纹时,接收点处信号出现了较强的瑞利反射回波。并且随着裂纹深度的增加,瑞利波信号的峰值呈显著增加的状态。因此这可以成为实验判断表面裂纹是否存在的主要依据。

      图 5中给出了瑞利反射波的幅值与裂纹深度之间的关系。可以看出,裂纹深度在0μm~260μm时,瑞利波的峰值随着深度的增加而显著增加;裂纹深度在260μm~420μm时,瑞利波的峰值随着深度的增加而缓慢增加;裂纹深度在420μm以上时,瑞利波的峰值随着深度的增加而基本保持不变。这一变化趋势[15]的原因是瑞利波的能量分布主要集中在表面以下、一个波长的深度范围内,随着裂纹深度的增加,瑞利波被裂纹阻挡,发生反射的能量成分增长,回波的幅值也随深度的增加而增加。当裂纹深度超过一个波长量级时,深度的增加将仅会对反射的能量成分产生微小的作用,此时回波的幅值会随深度增加而产生微小的浮动。这一趋势进一步验证了有限元模型的正确性,但同时也表明了不能简单地根据瑞利波的峰值来判断裂纹深度。

      Figure 5.  Peak values variation of Rayleigh wave signals with different surface crack depth

      为了使瑞利波的位移时间曲线能够反映与裂纹深度的关系,取反射波的正向峰值与下一次最大正向位移所对应的时间(t1t2图 6所示)的比值Δ=t2/t1为指标。图 7显示了该指标与裂纹深度之间的关系。

      Figure 6.  Displacement-time curve of node

      Figure 7.  Time ratio Δ versus surface crack depth

      图 7中的黑点为通过上述模型计算出来的比值,直线为拟合直线。从图 7可以清晰地看到,时间比值Δ与裂纹深度之间存在着近似线性的关系,并且线性相关系数为98.9%。通过这种关系,可以通过实验测量进而得到表面裂纹的深度。

    • 本文中建立了激光超声无损检测表面裂纹深度的有限元数值模型,考虑了激光工作参量将激光冲击作用等效成脉冲载荷作用,得到了性能良好的瑞利波。与传统的有限元模拟超声波的方法相比,引入无限单元能缩小模型尺寸,提高计算效率,并且也能有效消除边界反射波对计算的影响。模拟结果表明, 材料存在表面裂纹时, 接收点处会有明显的反射回波,借此可以判断材料表面是否存在裂纹。并且瑞利反射波的正向峰值与下一次最大正向位移所对应的时间比值Δ和裂纹深度之间存在线性关系。利用这一关系,可以为激光超声无损检测的实验测量提供一种有效方法。

参考文献 (15)

目录

    /

    返回文章
    返回