高级检索

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

留言板

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

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

激光扫描混凝土裂缝的超声检测

金谏 周若阳 赵金峰 潘永东 刘学增 桑运龙 顾盛

引用本文:
Citation:

激光扫描混凝土裂缝的超声检测

    作者简介: 金谏(1993-), 男, 硕士研究生, 现主要从事超声无损检测的研究.
    通讯作者: 潘永东, ypan@tongji.edu.cn
  • 基金项目:

    国家九七三重点基础研究发展计划资助项目 2015CB057806

    国家自然科学基金面上资助项目 11872282

    国家自然科学基金面上资助项目 51878497

  • 中图分类号: TN247

Ultrasonic detection of concrete cracks by laser scanning

    Corresponding author: PAN Yongdong, ypan@tongji.edu.cn ;
  • CLC number: TN247

  • 摘要: 为了探究激光扫描混凝土裂缝进行超声检测的可行性, 采用激光扫描混凝土表面激发声场、表面波探头接收超声信号的方法, 进行了相应的实验验证。通过对实验信号峰峰值和双极性特征分析, 实现了混凝土表面裂缝的定位。同时采用有限元方法模拟了扫描激光在材料表面激发声场的过程, 对超声信号进行了时域和频域上的分析。结果表明, 当激光与裂缝边沿距离在1mm时, 接收信号具有典型的双极性, 且峰峰值达到极大值; 当激光与裂缝边沿距离继续减小时, 峰峰值迅速降低。通过扫描激光激发超声, 在打点范围上形成由信号峰峰值表征的2维图像, 由峰峰值突变特征和信号的双极性特征, 可以实现混凝土表面缺陷边缘的定位。
  • Figure 1.  Schematic diagram of working procedure

    Figure 2.  Specimen size and scanning schematic diagram

    Figure 3.  Time domain signals with different xi coordinates at yi=12.5mm

    Figure 4.  Peak-to-peak value of time domain signals with different xi coordinates at yi=12.5mm

    Figure 5.  Peak-to-peak values of time domain signals with different xi coordinates at yi=10mm and 12.5mm, 15m

    Figure 6.  Peak-to-peak value imaging

    Figure 7.  Schematic diagram of finite element model

    Figure 8.  Comparison between damping and undamped signal when d=12mm, s=11mm

    Figure 9.  Time frequency graph of damping signal with d=12mm, s=11mm

    Figure 10.  Time-frequency graph of undamped signal with d=12mm, s=11mm

    Figure 11.  Comparisons between numerical and experimental signals

    a—typical test signal b—numerical signal

    Figure 12.  Acoustic signal at d=12mm and s=7mm

    Figure 13.  Time frequency diagram of acoustic signal at d=12mm and s=7mm

    Figure 14.  Time domain diagram when d=12mm, s=7mm~11m

    Figure 15.  Time domain diagram when d=12mm, s=9mm~11.8mm

    Figure 16.  Signal characteristics with 5mm defect depth

    Figure 17.  Signal characteristics with 7mm defect depth

    Figure 18.  Signal contrast when defect depth is 3mm, 5mm, 7mm and s=11mm

    Figure 19.  The change of signals' peak-to-peak value with s

    Table 1.  Material parameters of concrete for the finite element model


    Young modulus/
    GPa
    Poisson’s ratio density/(kg·m-3) thermal expansion
    coefficient/(10-6K-1)
    thermal conductivity/
    (W·m-1·K-1)
    specific heat capacity/
    (J·kg-1·K-1)
    4.5 0.2 2420 1 1.8 800
    下载: 导出CSV
  • [1]

    CHONGQING TRANSPORT COMMISSION. JTC H12-2003 Technical specification of maintenance for highway tunnel[S]. Beijing: China Communications Press, 2003: 7-17(in Chinese).
    [2]

    WANG R, QI T Y, WANG H X. Summary of crack detection methods for tunnel lining[J]. Sichuan Architecture, 2013, 33(3): 135-136(in Chinese).
    [3]

    YUAN W Q, XUE D. Review of tunnel lining crack detection algorithm based on machine vision[J]. Chinese Journal of Scientific Instrument, 2017, 38(12): 3100-3111(in Chinese). 
    [4]

    JIANG W B, LUO Q R, ZHANG X H. A review of concrete roads crack detection methods based on digital image[J]. Journal of Xihua University(Natural Science Edition), 2018, 37(1) :75-84(in Ch-inese). 
    [5]

    HUANG Y J, HE X J, WANG L L, et al. Research status of non destructive detection methon in concrete damage and cracking[J]. Architecture Technology, 2018, 49(1): 54-57(in Chinese).
    [6]

    LI Y L. Metal material internal defect detection with laser ultrasonic shear wave[D]. Nanjing: Nanjing University of Science and Techno-logy, 2014: 17-39(in Chinese).
    [7]

    FENG W W, JIN L, ZHAO J F, et al. Mode conversion of laser-excited shear waves interaction with the side of vertical cracks[J]. Laser Technology, 2018, 42(4): 487-493(in Chinese). 
    [8]

    LIU Y, XIANG Zh Q, TANG Zh F. Application of laser-induced ultrasound on rail flaw inspection[J]. Machinery Design & Manufacture, 2009(10): 60-61(in Chinese). 
    [9]

    SU K, REN D H, LI J, et al. Research of contacting micro-cracks detection technology[J]. Optical Technique, 2002, 28(6): 518-519 (in Chinese). 
    [10]

    ZHU H L, LIU Ch, ZHANG B, et al. Research on laser ultrasonic visual image processing[J]. Chinese Journal of Lasers, 2018, 45(1) : 0104004(in Chinese). doi: 10.3788/CJL
    [11]

    YASHIRO S, TAKATSUBO J, MIYAUCHI H, et al. A novel technique for visualizing ultrasonic waves in general solid media by pulsed laser scan[J]. NDT & E International, 2008, 41(2): 137-144. 
    [12]

    HALL J S, MCKEON P, SATYANARAYAN L, et al. Minimum variance guided wave imaging in a quasi-isotropic composite plate[J]. Smart Materials & Structures, 2011, 20(2): 025013. 
    [13]

    HALL J, MICHAELS J E. Minimum variance ultrasonic imaging app-lied to an in situ sparse guided wave array[J]. IEEE Transactions on Ultrasonics Ferroelectrics & Frequency Control, 2010, 57(10): 2311-2323. 
    [14]

    FAN W, QIAO P. A 2-D continuous wavelet transform of mode shape data for damage detection of plate structures[J]. International Journal of Solids & Structures, 2009, 46(25): 4379-4395. 
    [15]

    ZENG W, WANG H T, TIAN G Y, et al. Research on laser ultrasonic defect signal detection technologu based on energy analysis [J]. Chinese Journal of Scientific Instrument, 2014, 35(3): 650-655 (in Chinese).
    [16]

    KE G J, GUO Ch Q, HU Sh Q, et al. Study on the damping ratio of concrete[J]. Journal of Building Materials, 2004, 7(1): 35-40(in Chinese).
    [17]

    YAN G. Laser-generated surface acoustic waves technology for nondestructive detection of surface-breaking defects on metal[D]. Nanjing: Nanjing University of Science & Technology, 2004: 56-66(in Chinese).
    [18]

    GUAN J F. Mechanisms of laser generated surface acoustic wave and its application in the surface defect inspection[D]. Nanjing: Nanjing University of Science & Technology, 2006: 60-78(in Chinese).
  • [1] 宋艳马世榜张开飞 . 热弹激光超声激励及缺陷检测的有限元分析. 激光技术, 2021, 45(2): 246-251. doi: 10.7510/jgjs.issn.1001-3806.2021.02.021
    [2] 王玉庆王云霞马世榜 . 基于衍射横波的裂纹激光超声检测方法. 激光技术, 2019, 43(4): 546-550. doi: 10.7510/jgjs.issn.1001-3806.2019.04.019
    [3] 徐志祥黄建华王铮恭黄义敏王雨 . 激光超声检测带涂层金属表面裂纹的数值研究. 激光技术, 2018, 42(6): 801-805. doi: 10.7510/jgjs.issn.1001-3806.2018.06.014
    [4] 马健赵扬周凤艳孙继华 . 离焦量对激光超声测厚的影响. 激光技术, 2015, 39(3): 349-352. doi: 10.7510/jgjs.issn.1001-3806.2015.03.015
    [5] 许伯强刘洪凯徐桂东徐晨光李俊敏 . 基于应力-位移混合有限元法的激光超声数值模拟. 激光技术, 2014, 38(2): 230-235. doi: 10.7510/jgjs.issn.1001-3806.2014.02.018
    [6] 王贵鑫严刚关建飞 . 铝材管道管壁厚度对激光超声信号影响的探究. 激光技术, 2014, 38(2): 260-265. doi: 10.7510/jgjs.issn.1001-3806.2014.02.024
    [7] 颜鑫应恺宁戴鹭楠谭钧夫沈中华倪辰荫 . 基于物理信息神经网络的激光超声波场研究. 激光技术, 2024, 48(1): 105-113. doi: 10.7510/jgjs.issn.1001-3806.2024.01.017
    [8] 何宁谭智诚吕杏利蒋红艳 . 基于语音识别控制的激光超声水下遥感研究. 激光技术, 2017, 41(3): 391-395. doi: 10.7510/jgjs.issn.1001-3806.2017.03.017
    [9] 陈笑徐荣青沈中华陆建倪晓武 . 反射激光式超声测厚传感器及其应用. 激光技术, 2003, 27(6): 557-559.
    [10] 陈春明战宇刘常升方奔 . 激光超声测量T800/QY9511复合材料工程常数的有限元模拟. 激光技术, 2019, 43(1): 88-92. doi: 10.7510/jgjs.issn.1001-3806.2019.01.018
    [11] 马健赵扬郭锐宋江峰贾中青刘帅孙继华 . 激光辐照材料表层温升规律的数值模拟. 激光技术, 2013, 37(4): 455-459. doi: 10.7510/jgjs.issn.1001-3806.2013.04.009
    [12] 曹豆豆王开圣杨雁南 . 环状激光作用于薄管产生温度场的有限元模拟. 激光技术, 2010, 34(6): 753-756. doi: 10.3969/j.issn.1001-3806.2010.06.010
    [13] 周航张斌冯其波崔建英梁晨黄悦朗 . 环形光源激发超声进行缺陷检测的数值研究. 激光技术, 2021, 45(2): 168-173. doi: 10.7510/jgjs.issn.1001-3806.2021.02.007
    [14] 许伯强沈中华倪晓武关建飞陆建 . 涂层-基底系统激光激发表面波的时间-频率分析. 激光技术, 2004, 28(6): 609-612.
    [15] 蒋志伟龚时华王启行 . 双光束激光双路焊接的跟踪控制技术研究. 激光技术, 2013, 37(1): 1-5. doi: 10.7510/jgjs.issn.1001-3806.2013.01.001
    [16] 尹国福鲁建存刘彦义贺爱峰曹椿强 . 激光起爆系统光路完整性检测技术研究. 激光技术, 2011, 35(4): 554-558. doi: 10.3969/j.issn.1001-3806.2011.04.028
    [17] 赵兴海高杨程永生 . 激光点火技术综述. 激光技术, 2007, 31(3): 306-310,313.
    [18] 胡增荣童国权陈长军郭华锋周亮徐家乐 . 激光纳米表面工程技术. 激光技术, 2014, 38(6): 764-770. doi: 10.7510/jgjs.issn.1001-3806.2014.06.009
    [19] 易亚星李忠科邓方林姚森 . 牙颌模型激光扫描双轴系统设计. 激光技术, 2001, 25(2): 134-136.
    [20] 孙锴英赵圣之李桂秋杨克建 . 电光-GaAs主被动双调Q激光特性研究. 激光技术, 2008, 32(2): 204-206.
  • 加载中
图(19) / 表(1)
计量
  • 文章访问数:  7133
  • HTML全文浏览量:  5774
  • PDF下载量:  34
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-10-18
  • 录用日期:  2018-12-25
  • 刊出日期:  2019-07-25

激光扫描混凝土裂缝的超声检测

    通讯作者: 潘永东, ypan@tongji.edu.cn
    作者简介: 金谏(1993-), 男, 硕士研究生, 现主要从事超声无损检测的研究
  • 1. 同济大学 航空航天与力学学院, 上海 200092
  • 2. 中国建筑股份有限公司 北京 100029
  • 3. 同济大学 土木工程学院, 上海 200092
  • 4. 上海同岩土木工程科技股份有限公司, 上海 200092
  • 5. 上海地下基础设施安全检测与养护装备工程技术研究中心, 上海 200092
  • 6. 昆山市建设工程质量检测中心, 昆山 215337
基金项目:  国家九七三重点基础研究发展计划资助项目 2015CB057806国家自然科学基金面上资助项目 11872282国家自然科学基金面上资助项目 51878497

摘要: 为了探究激光扫描混凝土裂缝进行超声检测的可行性, 采用激光扫描混凝土表面激发声场、表面波探头接收超声信号的方法, 进行了相应的实验验证。通过对实验信号峰峰值和双极性特征分析, 实现了混凝土表面裂缝的定位。同时采用有限元方法模拟了扫描激光在材料表面激发声场的过程, 对超声信号进行了时域和频域上的分析。结果表明, 当激光与裂缝边沿距离在1mm时, 接收信号具有典型的双极性, 且峰峰值达到极大值; 当激光与裂缝边沿距离继续减小时, 峰峰值迅速降低。通过扫描激光激发超声, 在打点范围上形成由信号峰峰值表征的2维图像, 由峰峰值突变特征和信号的双极性特征, 可以实现混凝土表面缺陷边缘的定位。

English Abstract

    • 在公路、桥梁及隧道等交通基础设施的建设中, 混凝土发挥着重要作用。同时在混凝土结构表面的维护中, 裂缝是最常见病害之一。对裂缝的早期发现与检测具有现实意义[1-2]

      现阶段混凝土表面裂缝检测多采用视觉检测法和一些较为原始的无损检测方法[3-5],如超声脉冲法、脉冲回波法、雷达检测法、红外热谱法、声发射法。这些方法依赖工人的熟练程度,不仅大量耗费人力物力,在精度和效率方面存在很大劣势,而且无法发现微裂缝。根据国家对混凝土允许裂缝宽度限值的规定,微裂缝定义为宽度小于0.2mm的裂缝,是裂缝最根源的形态。

      激光超声是一种非接触式的无损检测方法,具有非接触、宽带、高灵敏度、高的空间分辨率等优点,广泛地被用来检测表面和埋藏的缺陷[6-8]。同时,激光超声技术作为一种宽带检测技术,可以通过精确测定超声位移幅度,检测到毫米量级的缺陷,甚至是微米量级的缺陷[9]。通过高速扫查的检测手段,结合图像处理方法以及信号分类方法,可以实现超声波的可视化,突显缺陷信息[10-14]

      本文中将有限元模拟与实验相结合,对激光扫描混凝土裂缝的超声检测技术进行了探究。在实验室条件下对混凝土试样尝试了2维成像方法;然后在有限元模拟中引入瑞利阻尼模拟了混凝土中的高频衰减特性,同时对存在缺陷时的信号特征进行了分析。

    • 本实验中采用实验室自主研发的激光超声场检测仪,包含激光激发和高速信号采集两个模块,超声波由输出波长为1064nm和脉宽为10ns的Nd:YAG激光器激发,激光脉冲通过焦距为100mm的扫描振镜聚焦镜头实现在试样表面的聚焦与扫点,超声波信号由压电传感器接收,再通过检测仪中的高速信号采集模块采样并存入计算机,采样频率为100MHz。检测流程如图 1所示。

      Figure 1.  Schematic diagram of working procedure

      混凝土试件按照C50标准制作,尺寸大小为150mm×150mm×150mm,用特定尺寸的铁片在表面预制宽度为0.2mm的裂缝,裂缝深度10mm。通过激光扫描控制程序调节打点范围和间隔。如图 2所示,激光沿y轴共计打点50排,每排沿x轴扫描打点50次。打点间隔标定为0.508mm±0.005mm。

      Figure 2.  Specimen size and scanning schematic diagram

      关于接收探头的选择,金属常用的探头工作频率通常在1MHz以上,不适合混凝土材料的检测,且普通的低频直探头半径过大,对试样尺寸要求高,因此从信号特征和实验条件考虑,尝试使用最高工作频率为500kHz的表面波探头在距离裂缝左边界17.5mm处接收信号。

      yi一定时,可得到一组不同xi坐标下的信号波形及对应的传播时间ti。取yi=12.5mm时,不同xi坐标下的时域信号作图。当xi坐标大于12.5mm(见图 3中长虚线以上),此时激光源处于缺陷与接收探头之间时,激发方式为近场激发,信号具有明显的双极性特征。由信号的时间差和传播距离,可以求得该信号传播速率为2.11km/s,与理论表面波波速2.22km/s的误差率为4.96%。

      Figure 3.  Time domain signals with different xi coordinates at yi=12.5mm

      取时域信号的峰峰值作图 4。可以观察到在xi坐标从12.5mm~13.5mm时,峰峰值由0.461V增加到1.169V,当xi坐标小于12.5mm时,信号特征以背景噪声为主,整体能量较小;当xi=13.5mm时,信号呈双极性特征,与xi坐标小于12.5mm的信号相比信号峰峰值发生了突增。

      Figure 4.  Peak-to-peak value of time domain signals with different xi coordinates at yi=12.5mm

      yi为10mm, 12.5mm, 15mm时,不同xi坐标的时域信号峰峰值如图 5所示。由图可知,信号峰峰值在xi=13.5mm~14mm范围内取得最大值,且峰峰值在xi=12.5mm后突增,该位置为裂缝的边缘。此时,激光处于裂缝与接收点之间,激光光源近场激发声场。

      Figure 5.  Peak-to-peak values of time domain signals with different xi coordinates at yi=10mm and 12.5mm, 15m

      对打点矩阵得到的信号做滤波等预处理后,可以得到打点矩阵上每一点信号的峰峰值,将峰峰值大小由色阶表示,可以在打点范围上形成由信号峰峰值表征的2维图像,如图 6所示。图中xi坐标为12.5mm与14mm之间存在明显的能量反差,与实际裂缝边缘位置一致。这与参考文献中提到的能量分布技术得到的结论一致[15]。当xi坐标继续增大,信号以直达表面波信号为主导,信号能量继续增加,如果进一步调整阈值,可以反映出更真实、明显的裂缝形貌。

      Figure 6.  Peak-to-peak value imaging

      2维图像中信号表征的裂缝并不占满整个打点区域,这可能是由于表面波接收探头在接收角度上对上下区域的信号不敏感造成的,另外由于接收探头的接收特性,可以发现与远离探头中心线时接收到的信号能量相比,在探头中心所对应的中线yi=12.5mm±2mm范围内,信号能量更大。

    • 为了进一步研究激光近场激发混凝土信号特征,采用有限元法对激光在材料表面激发超声波进行研究。考虑到混凝土材料的复杂性以及有限元模拟的效率,现将混凝土当作各向同性均匀材料进行研究。

      声表面波的激光激发机理主要有两种,在目前的理论中,主要是针对热弹模型,理论相对完备,融蚀模型的完备解必须建立在融蚀模型的基础之上。在实际情况中,两种机制往往同时起作用。本文中有限元采用热弹模型。在线弹性范围内,材料的热弹场是由热弹耦合方程所控制的,具体方程如下:

      $ \kappa {▽^2}T = \rho {c_V}\dot T + {T_0}\beta ▽ \cdot \mathit{\boldsymbol{\dot u}} - q $

      (1)

      $ \begin{array}{*{20}{c}} {{{▽}^2}\mathit{\boldsymbol{u}}(x, t) + (\lambda + \mu )▽[▽ \cdot \mathit{\boldsymbol{u}}(x, t)] = }\\ {p\mathit{\boldsymbol{\ddot u}} + \beta ▽T} \end{array} $

      (2)

      式中,TT0分别是材料的实际温度和环境温度,u(x, t)是材料内部不同时刻的位移向量,κ是热传导系数,cV是材料的比定容热容,β是热弹耦合系数, 且β=(3λ+2μ)αα是线性热膨胀系数,q表示单位体积在单位时间内吸收的热量,ρ表示材料密度,λμ表示材料的拉梅常数。▽为向量微分算子,其形式化定义为$▽ = \frac{{\rm{d}}}{{{\rm{d}}\mathit{\boldsymbol{r}}}}$,在n维空间中,分母dr为含有n个分量的向量。$\mathit{\boldsymbol{\dot u}}, \mathit{\boldsymbol{\ddot u}}$分别为u(x, t)对时间变量t的1阶求导与2阶求导。当激光辐照表面产生超声波时,${T_0}\beta ▽ \cdot \mathit{\boldsymbol{\dot u}}$超声场对温度场的影响项,而βT则是温度场对超声场的影响项。本文中采用耦合求解方式,考虑超声场和温度场之间的相互影响。

    • 采用ABAQUS有限元软件中的Explicit显式求解器求解瞬态的波传播问题。图 7为2维半无限平面中激光激发超声波的有限元模型,其中平面左、右、下侧均为无限边界。在模型中采用规则的矩形凹痕来代替实际材料中的表面缺陷,为与实验进行对照,另外根据微裂缝定义和国家对混凝土裂缝宽度的允许限值,将凹痕的宽度取为0.2mm。为叙述方便,d表示接收点至缺陷的距离,s表示接收点至激光源的距离。激光源处于传感器和缺陷之间。

      Figure 7.  Schematic diagram of finite element model

      有限元模型中采用4节点平面矩形单元。MOSER等人推荐每个波长范围内应包含20个单元,这样才能有效捕捉波动效应。在本次研究中,网格大小取为10μm。考虑到计算的稳定性,时间步长取为0.2ns。根据实验中使用的激光检测仪实际参量,有限元脉冲激光的上升时间为5ns,光斑半径取为100μm。计算中,混凝土材料的热物理参量见表 1

      Table 1.  Material parameters of concrete for the finite element model


      Young modulus/
      GPa
      Poisson’s ratio density/(kg·m-3) thermal expansion
      coefficient/(10-6K-1)
      thermal conductivity/
      (W·m-1·K-1)
      specific heat capacity/
      (J·kg-1·K-1)
      4.5 0.2 2420 1 1.8 800
    • 阻尼比是用以描述材料或结构在振动过程中某种能量耗散方式的术语,是材料或结构的动力特性参量之一,阻尼比的大小对动态响应有非常明显的影响。

      一般情况下,阻尼矩阵是依赖于频率的,但在实际分析中,要精确地确定阻尼矩阵是相当困难的。实际计算过程中通常采用瑞利(Rayleigh)阻尼。瑞利阻尼假设结构的阻尼矩阵C是质量矩阵[M]和刚度矩阵[K]的组合,即:

      $ \mathit{\boldsymbol{C}} = {\mathit{\boldsymbol{\alpha }}_{\rm{R}}}[\mathit{\boldsymbol{M}}] + {\beta _{\rm{R}}}[\mathit{\boldsymbol{K}}] $

      (3)

      式中, αRβR是两个比例系数,分别具有s-1和s的量纲,不依赖于频率。其中第1项是针对低频的粘性阻尼,计算时可以取0;第2项是针对高频的阻尼。瑞利阻尼系数αRβR与阻尼比ξn之间的关系为:

      $ {\xi _n} = \frac{{{\alpha _{\rm{R}}}}}{{2{\omega _n}}} + \frac{{{\beta _{\rm{R}}}{\omega _n}}}{2} $

      (4)

      式中,ωn为材料的固有频率。

      在ABAQUS材料属性中试取瑞利阻尼系数αR=0s-1βR=1×10-8s。模型中凹痕宽度为0.2mm,深度为3mm,接收点至缺陷的距离d=12mm,接收点至激光源的距离s=11mm。有限元脉冲激光的上升时间为5ns,光斑半径取为100μm。其余模型参量如上面所述。

      图 8~图 10显示了添加瑞利阻尼与无阻尼时有限元结果,同时该信号将在下面详细展开。由图 8可知,位移信号中的高频信号在添加瑞利阻尼的条件下得到有效抑制,且信号峰值减小。另外,由小波变换得到的时域信息中可知,存在瑞利阻尼条件下,信号频率范围在0.1MHz~2MHz,中心频率为0.65MHz,由(4)式可反求材料的阻尼系数为1%。普通混凝土的阻尼比在0.9%~2.0%之间。有限元计算与文献中的结论一致[16],说明瑞利阻尼试取合适。

      Figure 8.  Comparison between damping and undamped signal when d=12mm, s=11mm

      Figure 9.  Time frequency graph of damping signal with d=12mm, s=11mm

      Figure 10.  Time-frequency graph of undamped signal with d=12mm, s=11mm

      为了验证有限元模拟的正确性,根据实验接收条件,对有限元信号进行低通滤波,截止频率为500kHz。信号对比如图 11所示。图 11a为实验中扫描激光源与接收探头在裂缝同侧时的信号,矩形窗内信号为激光近场激发时接收到的首波信号,时长10μs,与有限元信号时长一致;图 11b展示了图 9中数值模拟信号经过低通滤波后的时域波形,观察可知,数值模拟信号波形特征与实验一致,验证了有限元模拟的正确性。

      Figure 11.  Comparisons between numerical and experimental signals

    • 在有限元模型中添加瑞利阻尼,利用有限元模拟移动激光源扫描混凝土进行超声检测的过程。激发激光源从s=7mm位置开始向缺陷位置移动。当d=12mm,s=7mm时的声信号如图 12所示。图中,R表示直达表面波,Rr表示反射回波。RRr峰之间的振荡由缺陷产生。对表面波频率进行分析,由小波变化得到信号的时频图如图 13所示。由图 13可知,表面波中心频率为626kHz,中心频率信号能量最大的时间在3.308μs,与时域图上直达表面波信号R峰值出现时间一致。可以观察到,Rr信号由于传播过程中的高频衰减,中心频率减小为427kHz,且峰值只有直达波信号R的35.79%。由此可见,能量衰减和高频衰减特性是混凝土材料中超声信号与金属材料的主要区别。

      Figure 12.  Acoustic signal at d=12mm and s=7mm

      Figure 13.  Time frequency diagram of acoustic signal at d=12mm and s=7mm

      设置近场范围内激光作用在不同的激发位置模拟扫描过程,分别建立平面应变有限元模型。模型中凹痕宽度0.2mm,深度3mm,当d=12mm,s分别为7mm(见曲线1), 8mm(见曲线2), 9mm(见曲线3), 10mm(见曲线4), 11mm(见曲线5)时的时域图如图 14所示。图中可以清晰地读取直达波信号R和反射信号Rr的到达时,由信号传播距离和时间差可知,RRr信号的波速为2.24km/s,符合材料中的表面波波速。由混凝土材料参量计算可得表面波理论波速为2.22km/s,有限元误差率为0.90%。

      Figure 14.  Time domain diagram when d=12mm, s=7mm~11m

      在得到表面波的波速与频率的基础上,可以求得该混凝土有限元模型中,当激发条件为脉冲激光的上升时间5ns,光斑半径取为100μm时,直达表面波波长为3.58mm。该波长可指导后续的研究。

      另外,在激光源逼近缺陷边沿过程中,在RRr峰之间的振荡也越来越明显,这是由于表面缺陷自由边界的反射增强,同时,激光激发多种声波模式受缺陷的影响发生了模式转换,在近场区域叠加在一起,从而引起信号幅值的显著增大[17-18]

      读取信号的峰峰值,并在s为9mm~12mm的范围内以0.2mm为单位,增加研究数量,呈现更丰富的细节。部分时域图如图 15所示。从时域图上看,当s < 11mm,仍可以清晰地分辨出直达波信号和反射波信号,但当信号峰峰值达到最大值后,当激光源继续向缺陷移动,s为11mm,11.4mm,11.8mm对应信号中,原直达波信号时间分别为4.96μs, 5.06μs, 5.14μs,原反射波时间分别为6.176μs, 6.196μs, 6.248μs,由信号传播距离和时间差可得,RRr信号的波速大于5km/s,与表面波波速2.22km/s不一致,也与纵波波速3.78km/s、横波波速2.44km/s不同,是在激光近场激发这种特殊条件下,叠加之后形成的声信号。

      Figure 15.  Time domain diagram when d=12mm, s=9mm~11.8mm

      同时,由图 15可以发现,s为11mm, 11.4mm, 11.8mm时,信号中的正极性峰值到达时不改变。

      当缺陷深度为5mm, 7mm时,保持接收与激发条件与上面一致时,接收点至激光源的距离分别为10.5mm, 11mm, 11.5mm的时域图如图 16图 17所示。当缺陷深度为5mm、峰峰值对应s为10.5mm, 11mm, 11.5mm时, 分别为0.3160×10-12m, 0.3128×10-12m, 0.2233×10-12m;当缺陷深度为7mm、峰峰值对应s为10.5mm, 11mm, 11.5mm时, 分别为0.3187×10-12m, 0.3135×10-12m, 0.2236×10-12m。当缺陷深度改变时,信号仍然在s=11mm处拥有典型的双极性特征,随后在s=11.5mm时峰峰值降低,正极性信号的峰值到达时不改变。另外,当缺陷深度大于表面波波长时,表面波反射能量增强,使得在s=10.5mm处峰峰值大于s=11mm处的峰峰值,但该现象不会影响根据峰峰值判断缺陷位置。

      Figure 16.  Signal characteristics with 5mm defect depth

      Figure 17.  Signal characteristics with 7mm defect depth

      当缺陷深度改变时,固定接收点和激发点,信号对比如图 18所示。峰峰值的改变可以忽略不记,即峰峰值不随缺陷深度的改变而改变。而双极性信号之后的振荡信号包含了缺陷深度的信息。

      Figure 18.  Signal contrast when defect depth is 3mm, 5mm, 7mm and s=11mm

      下面对峰峰值-接收点至激光源的距离关系图(见图 19)进行详细的说明,凹痕宽度为0.2mm、深度为3mm、d=12mm、s为4mm~9.8mm时,由于混凝土中的衰减特性,首波幅值减少,使得峰峰值降低;当激光源近一步接近缺陷边缘时,波在缺陷处的反射增强,同时叠加上缺陷处的声波模式转化,形成一个正向的波形,此时,激光源与缺陷的距离为2.2mm,在一个波长以内,由路程和波速可知,直达表面波与叠加后的波形到达时间差在2μs以内,以上因素最终使得波形呈双极性的特征,并在s=11mm处取得最大值;当激光源继续接近缺陷边缘,由于激光源空间上的高斯分布特性,激发条件在近缺陷处发生改变,入射能量减少,峰峰值减低,并在缺陷边沿处取得最小值。

      Figure 19.  The change of signals' peak-to-peak value with s

      综上所述,当激光扫描带有裂缝的混凝土表面,表面波探头接收超声信号,由于缺陷的存在使得信号会有峰峰值的改变,可以根据峰峰值-接收点至激光源的距离图判断缺陷位置。

    • 本文中探究了激光扫描混凝土裂缝的超声检测技术,实现了混凝土裂缝的快速定位。

      (1) 在实验中,通过激光超声场检测仪扫描含有裂缝的混凝土表面,使用最高工作频率为500kHz的表面波探头接收超声信号。通过激光矩阵激发的表面波信号的峰峰值成像,结合信号的双极性特征,可以检测到混凝土表面的裂缝边界。

      (2) 采用有限元方法模拟并分析了激光在混凝土中激发的声场,发现当固定接收点,随着接收点至激光源的距离的改变,信号峰峰值有3个变化区间。另外,双极性特征只与激发位置有关,而与裂缝深度没有直接关系。

参考文献 (18)

目录

    /

    返回文章
    返回