第36卷第6期 北京科技大学学报 Vol.36 No.6 2014年6月 Journal of University of Science and Technology Beijing Jun.2014 金属内部裂纹的超声衰减特性与近表面缺陷识别 肖会芳)四,徐金梧”,何飞”,黎敏 1)北京科技大学国家板带生产先进装备工程技术研究中心,北京1000832)北京科技大学机械工程学院,北京100083 ☒通信作者,E-mail:huifangxiao@usth.edu.cn 摘要建立了包含裂纹缺陷的二维金属板模型。采用有限元方法,对具有不同深度裂纹的材料内部超声波场进行计算,获 得了不同裂纹深度金属材料上表面的回波信号,分析了裂纹深度对超声波传播特性的影响规律.进一步分离提取不同阶次底 面回波的频谱特征,获得了由裂纹缺陷引起的超声衰减系数随频率的变化关系.最后提出了通过底面回波频谱图辨识近表面 裂纹缺陷的方法. 关键词近表面:裂纹:无损检测;超声;衰减:辨识 分类号TG115.28*5 Ultrasonic attenuation characteristics of metal materials with cracks and identifi- cation of subsurface flaws XA0 Hui fang》,XU Jin+-,HE Fei'》,山Min2 1)National Engineering Research Center of Flat Rolling Equipment,University of Science and Technology Beijing.Beijing 100083,China 2)School of Mechanical Engineering,University of Science and Technology Beijing,Beijing 100083,China Corresponding author,E-mail:huifangxiao@ustb.edu.cn ABSTRACT A two-dimensional isotropic plate model with an embedded horizontal crack was established.Finite element analysis with absorbing boundary condition was employed to simulate the experimental pulse-echo mode to obtain ultrasonic data.The effect of crack depth on the wave propagation was studied and the spectra of different backwall echoes were further extracted.A relationship was determined between ultrasonic attenuation due to the crack and frequency.Finally,a method was proposed for identifying subsurface cracks based on the spectra of backwall echoes KEY WORDS subsurface;cracks;non-destructive testing:ultrasonic:attenuation:identification 超声无损检测技术(ultrasonic non-destructive 识近表面裂纹缺陷,对超声无损检测带钢缺陷具有 testing,UNDT),是一种利用超声波在材料中的传播 重要的理论和实际意义. 特性来判断被测材料的缺陷和异常的无损检测方 本文通过建立包含裂纹缺陷的二维金属板模 法n-).其中,超声脉冲回波检测法(pulse-echo 型,采用有限元数值仿真的方法对具有不同深度裂 mode)以测量精度高、检测速度快而广泛应用于金 纹的材料内部超声波场进行计算,获得了不同裂纹 属材料内部裂纹缺陷无损检测. 深度金属上表面的回波信号,研究了裂纹深度对超 然而,对金属内部的近表面裂纹缺陷,由于裂纹 声波传播特性的影响规律,并进一步分离提取底面 回波信号会淹没在表面回波信号中而无法辨识缺 回波的频域特征,获得了由裂纹缺陷引起的超声衰 陷,导致存在检测盲区,成为超声回波检测法最主要 减系数随频率的变化关系,提出了通过底面回波频 的缺点和难点之一9.如何实现从回波信号中辨 谱图辨识近表面裂纹缺陷的方法,为检测近表面裂 收稿日期:201306-26 基金项目:“十二五”国家科技支撑计划资助项目(2012BAF04B02):国家自然科学基金青年科学基金资助项目(51304019) DOI:10.13374/j.issn1001-053x.2014.06.019:http://journals.ustb.edu.cn
第 36 卷 第 6 期 2014 年 6 月 北京科技大学学报 Journal of University of Science and Technology Beijing Vol. 36 No. 6 Jun. 2014 金属内部裂纹的超声衰减特性与近表面缺陷识别 肖会芳1) ,徐金梧1) ,何 飞1) ,黎 敏2) 1) 北京科技大学国家板带生产先进装备工程技术研究中心,北京 100083 2) 北京科技大学机械工程学院,北京 100083 通信作者,E-mail: huifangxiao@ ustb. edu. cn 摘 要 建立了包含裂纹缺陷的二维金属板模型. 采用有限元方法,对具有不同深度裂纹的材料内部超声波场进行计算,获 得了不同裂纹深度金属材料上表面的回波信号,分析了裂纹深度对超声波传播特性的影响规律. 进一步分离提取不同阶次底 面回波的频谱特征,获得了由裂纹缺陷引起的超声衰减系数随频率的变化关系. 最后提出了通过底面回波频谱图辨识近表面 裂纹缺陷的方法. 关键词 近表面; 裂纹; 无损检测; 超声; 衰减; 辨识 分类号 TG115. 28 + 5 Ultrasonic attenuation characteristics of metal materials with cracks and identification of subsurface flaws XIAO Hui-fang1) ,XU Jin-wu1) ,HE Fei 1) ,LI Min2) 1) National Engineering Research Center of Flat Rolling Equipment,University of Science and Technology Beijing,Beijing 100083,China 2) School of Mechanical Engineering,University of Science and Technology Beijing,Beijing 100083,China Corresponding author,E-mail: huifangxiao@ ustb. edu. cn ABSTRACT A two-dimensional isotropic plate model with an embedded horizontal crack was established. Finite element analysis with absorbing boundary condition was employed to simulate the experimental pulse-echo mode to obtain ultrasonic data. The effect of crack depth on the wave propagation was studied and the spectra of different backwall echoes were further extracted. A relationship was determined between ultrasonic attenuation due to the crack and frequency. Finally,a method was proposed for identifying subsurface cracks based on the spectra of backwall echoes. KEY WORDS subsurface; cracks; non-destructive testing; ultrasonic; attenuation; identification 收稿日期: 2013--06--26 基金项目: “十二五”国家科技支撑计划资助项目( 2012BAF04B02) ; 国家自然科学基金青年科学基金资助项目( 51304019) DOI: 10. 13374 /j. issn1001--053x. 2014. 06. 019; http: / /journals. ustb. edu. cn 超声无损检测技术( ultrasonic non-destructive testing,UNDT) ,是一种利用超声波在材料中的传播 特性来判断被测材料的缺陷和异常的无损检测方 法[1 - 3]. 其 中,超声脉冲回波检测法 ( pulse-echo mode) 以测量精度高、检测速度快而广泛应用于金 属材料内部裂纹缺陷无损检测. 然而,对金属内部的近表面裂纹缺陷,由于裂纹 回波信号会淹没在表面回波信号中而无法辨识缺 陷,导致存在检测盲区,成为超声回波检测法最主要 的缺点和难点之一[4 - 5]. 如何实现从回波信号中辨 识近表面裂纹缺陷,对超声无损检测带钢缺陷具有 重要的理论和实际意义. 本文通过建立包含裂纹缺陷的二维金属板模 型,采用有限元数值仿真的方法对具有不同深度裂 纹的材料内部超声波场进行计算,获得了不同裂纹 深度金属上表面的回波信号,研究了裂纹深度对超 声波传播特性的影响规律,并进一步分离提取底面 回波的频域特征,获得了由裂纹缺陷引起的超声衰 减系数随频率的变化关系,提出了通过底面回波频 谱图辨识近表面裂纹缺陷的方法,为检测近表面裂
第6期 肖会芳等:金属内部裂纹的超声衰减特性与近表面缺陷识别 ·839· 纹缺陷提供指导 纹表面相互作用并发生反射,反射回波被检测探头 接收.由于裂纹回波的传播距离小于底面回波,因 1二维金属板模型 此其到达接收探头的时间较早.假设从激发到接收 1.1模型描述 裂纹回波信号的时间间隔为t。,纵波波速为C。,当 研究材料内部裂纹缺陷对超声波传播及衰减特 时间t大于始波或底面回波的时间宽度时,裂纹 性影响规律的模型示意图,如图1所示.二维矩形 的深度位置可估算为h=C,t。2,通过裂纹回波可以 金属板的宽为B,厚度为H.宽度为l、高度为d的矩 提取裂纹的深度信息,如图2(a)所示. 形裂纹位于金属板的半宽位置,且裂纹方向平行于 但是,当裂纹位于近表面位置,使得从激发到接 金属板表面,因而裂纹在金属板内部的位置仅由其 收裂纹回波信号的时间t。小于始波的时间宽度 深度h决定.沿宽度方向等幅分布的法向激励载荷 (如图2(b)所示),或缺陷回波宽度1小于底面回波 作用在金属板的上表面,引起表面的振动并向金属 的时间宽度t(如图2(c)所示),则裂纹回波信号淹 板内部传播,产生平行于y方向的纵波L和平行于 没在始波或底面回波信号中,从而无法直接从回波 x方向的横波S,并与金属板的上、下表面以及裂纹 信号中辨识缺陷的存在,即存在缺陷检测的盲区. 相互作用.由于横波$几乎不与裂纹相互作用,因 同时,裂纹会引起超声波的散射和衰减,对底面 此主要表现为纵波L与裂纹的相互作用,纵波在裂 回波也会产生影响,因而底面回波也包含了丰富的 纹表面反射的同时,在裂纹尖端发生散射,并产生波 裂纹特征信息.文献6]研究了裂纹宽度对反射回 型转换为横波6-刀 波的影响.结果表明:当裂纹宽度较小时(13),反射系数随着裂纹宽度递减.本文 保持裂纹尺寸不变,通过建立不同裂纹深度h的金 裂纹 属板模型,采用有限元方法,获得不同裂纹深度的反 射回波曲线 1.2有限元计算模型 采用商用有限元软件ABAQUS求解超声波在 板内的传播.钢板的尺寸为5mm×5mm,对应的材 料参数为弹性模量E=200GPa,泊松比v=0.28,密 度p=7850kgm-3.根据Christoffel方程阅,金属板 图1具有内部裂纹缺路的二维金属板模型示意图 内超声波的波速分别为纵波C。=5707m·s,横波 Fig.1 Sketch diagram of a two-dimensional isotropic panel with an C.=3155ms1.为了模拟超声换能器声场,在金属 embedded crack 板上表面施加法向的瞬时压力载荷P(t).激励信 通过检测上表面的回波信号,采用同端激励、同 号P(t)为Hanning窗调制,中心频率为f=l0MHz, 端接收的方式,模拟实际超声检测中的脉冲回波法 时间长度为△1=200s的单周期正弦脉冲,其时域波 当内部存在裂纹时,超声波在到达底面之前,先与裂 形和频谱图,如图3所示.中心频率为10MH的探头 A) 底面回波 底面回波 4c) 底面回波 始波 始波 始波 缺陷回波 缺陷回波 缺箔向波 图2裂纹回波时间t,与表面回波的时间宽度h之间的关系.()t。>h:(b)。<t,裂纹近上表面:(c)t<,裂纹近下表面 Fig.2 Relationship between time t and:(a)t (b)t<,the crack is close to the frontwall;(b)t<,the crack is close to the backwall
第 6 期 肖会芳等: 金属内部裂纹的超声衰减特性与近表面缺陷识别 纹缺陷提供指导. 1 二维金属板模型 1. 1 模型描述 研究材料内部裂纹缺陷对超声波传播及衰减特 性影响规律的模型示意图,如图 1 所示. 二维矩形 金属板的宽为 B,厚度为 H. 宽度为 l、高度为 d 的矩 形裂纹位于金属板的半宽位置,且裂纹方向平行于 金属板表面,因而裂纹在金属板内部的位置仅由其 深度 h 决定. 沿宽度方向等幅分布的法向激励载荷 作用在金属板的上表面,引起表面的振动并向金属 板内部传播,产生平行于 y 方向的纵波 L 和平行于 x 方向的横波 S,并与金属板的上、下表面以及裂纹 相互作用. 由于横波 S 几乎不与裂纹相互作用,因 此主要表现为纵波 L 与裂纹的相互作用,纵波在裂 纹表面反射的同时,在裂纹尖端发生散射,并产生波 型转换为横波[6 - 7]. 图 2 裂纹回波时间 tc与表面回波的时间宽度 tb之间的关系. ( a) tc > tb ; ( b) tc < tb,裂纹近上表面; ( c) tc < tb,裂纹近下表面 Fig. 2 Relationship between time tc and tb : ( a) tc > tb ; ( b) tc < tb,the crack is close to the frontwall; ( b) tc < tb,the crack is close to the backwall 图 1 具有内部裂纹缺陷的二维金属板模型示意图 Fig. 1 Sketch diagram of a two-dimensional isotropic panel with an embedded crack 通过检测上表面的回波信号,采用同端激励、同 端接收的方式,模拟实际超声检测中的脉冲回波法. 当内部存在裂纹时,超声波在到达底面之前,先与裂 纹表面相互作用并发生反射,反射回波被检测探头 接收. 由于裂纹回波的传播距离小于底面回波,因 此其到达接收探头的时间较早. 假设从激发到接收 裂纹回波信号的时间间隔为 tc,纵波波速为 Cp,当 时间 tc大于始波或底面回波的时间宽度 tb时,裂纹 的深度位置可估算为 h = Cp tc /2,通过裂纹回波可以 提取裂纹的深度信息,如图 2( a) 所示. 但是,当裂纹位于近表面位置,使得从激发到接 收裂纹回波信号的时间 tc 小于始波的时间宽度 tb ( 如图 2( b) 所示) ,或缺陷回波宽度 te小于底面回波 的时间宽度 td ( 如图 2( c) 所示) ,则裂纹回波信号淹 没在始波或底面回波信号中,从而无法直接从回波 信号中辨识缺陷的存在,即存在缺陷检测的盲区. 同时,裂纹会引起超声波的散射和衰减,对底面 回波也会产生影响,因而底面回波也包含了丰富的 裂纹特征信息. 文献[6]研究了裂纹宽度对反射回 波的影响. 结果表明: 当裂纹宽度较小时( l < 3λ,λ 为纵波波长) ,反射系数( 定义为反射波幅值与入射 波幅值的比值) 随着裂纹宽度递增; 当裂纹宽度较 大时( l > 3λ) ,反射系数随着裂纹宽度递减. 本文 保持裂纹尺寸不变,通过建立不同裂纹深度 h 的金 属板模型,采用有限元方法,获得不同裂纹深度的反 射回波曲线. 1. 2 有限元计算模型 采用商用有限元软件 ABAQUS 求解超声波在 板内的传播. 钢板的尺寸为 5 mm × 5 mm,对应的材 料参数为弹性模量 E = 200 GPa,泊松比 v = 0. 28,密 度 ρ = 7850 kg·m - 3 . 根据 Christoffel 方程[8],金属板 内超声波的波速分别为纵波 Cp = 5707 m·s - 1 ,横波 Cs = 3155 m·s - 1 . 为了模拟超声换能器声场,在金属 板上表面施加法向的瞬时压力载荷 P( t) . 激励信 号 P( t) 为 Hanning 窗调制,中心频率为 fc = 10 MHz, 时间长度为 Δt = 200 ns 的单周期正弦脉冲,其时域波 形和频谱图,如图 3 所示. 中心频率为 10 MHz 的探头 ·839·
·840· 北京科技大学学报 第36卷 是中厚板缺陷检测使用最为广泛的探头之一4,网 型的宽度远小于实际带钢宽度,为了消除计算模型 对中心频率为I0MHz的激励载荷,纵波的波长 宽度尺寸有限而引入的边界反射效应对波场的影 为入=0.57mm.采用四节点CPS4R单元离散,为了 响,在金属板的左右边界采用了吸收边界0-山.求 精确模拟超声波的传播,网格尺寸的大小选为入/ 解脉冲激励载荷下的波场,提取金属板上表面所有 60,即每个纵波波长内有60个单元.时间步长为 节点在y方向的位移,并对所有节点位移求平均值. 0.1ns,以保证在单个时间步内,即使波速最大的纵 这与实际超声测量中换能器作用面积上所接收到的 波的传播距离也小于一个网格尺寸.由于有限元模 回波声压能量相对应 1.0 01 0.6 0.5 0.4 0.3 0.2 0.1 100 150 200 15 30 t/ns f/MHz 图3瞬时激励载荷曲线.(a)时域:(b)颍域 Fig.3 Transient excitation load in time (a)and frequency (b)domains 1.3有限元计算方法有效性验证 起实验测试结果和仿真结果差异的原因 采用有限元方法模拟无缺陷金属板上表面反射 0.08 仿真 回波,并与实验测量结果进行对比以检验本文采用 0.06 实验 有限元方法模拟金属内超声波传播的有效性.金属 0.04 板上表面所有节点y方向位移的平均值,如图4所 0.02 示.图4显示,第1次底面回波到达上表面的时间t =1.80μs,相邻两次底面回波之间的时间间隔也为 -0.02 t=1.80us,与理论计算值1=1.75us基本一致.由 -0.04 于有限元计算采用显示中心差分算法,在数值计算 过程中存在舍入误差和截断误差,从而引起计算中 -0.06 的数值噪声☒ -0.08 为了验证有限元计算结果的准确性,采用IOS t/us AIR-1550TWM激光超声检测设备测量了无缺陷钢 图4厚度为5m无缺陷金属板y方向波形的仿真结果和实测 板的超声波信号.激光发生器的波长为532nm,发 结果对比 Fig.4 Comparison of the numerical and experimental reflected ydis- 射的激光垂直作用在被测钢板表面,产生的超声信 placement in a defect-free plate with the thickness of 5 mm 号被带宽为125MHz的激光超声接收仪接收.测量 钢板为矩形断面,厚度H=5mm,宽度B=20mm. 2 计算结果与分析 采样点数为500,采样率为0.02us,信号总长度为 10μs.测量获得的1~8μs时间内钢板上表面的y 2.1内部无裂纹 方向位移响应,如图4所示.仿真结果和实验结果 对图4所示的无裂纹金属板上表面仿真曲线的 的各底面回波到达时间基本完全一致,验证了本文 四次底面回波进行快速傅里叶变换后的频域曲线如 有限元方法的准确性。但是,由于仿真计算忽略了 图5所示.图5显示,对入射中心频率为10MHz的 实际材料内部品界多重背散射引起的衰减,因此仿 载荷,随着超声波在金属板内传播,底面回波的频率 真计算结果的回波幅值衰减要小于实测结果.同 范围扩大,不同底面回波的频率范围均在0~20 时,实验测试的噪声和测试样品的有限宽度也是引 MHz.图5同时显示,随着回波次数的增加,不同频
北 京 科 技 大 学 学 报 第 36 卷 是中厚板缺陷检测使用最为广泛的探头之一[4,9]. 对中心频率为 10 MHz 的激励载荷,纵波的波长 为 λ = 0. 57 mm. 采用四节点 CPS4R 单元离散,为了 精确模拟超声波的传播,网格尺寸的大小选为 λ / 60,即每个纵波波长内有 60 个单元. 时间步长为 0. 1 ns,以保证在单个时间步内,即使波速最大的纵 波的传播距离也小于一个网格尺寸. 由于有限元模 型的宽度远小于实际带钢宽度,为了消除计算模型 宽度尺寸有限而引入的边界反射效应对波场的影 响,在金属板的左右边界采用了吸收边界[10 - 11]. 求 解脉冲激励载荷下的波场,提取金属板上表面所有 节点在 y 方向的位移,并对所有节点位移求平均值. 这与实际超声测量中换能器作用面积上所接收到的 回波声压能量相对应. 图 3 瞬时激励载荷曲线. ( a) 时域; ( b) 频域 Fig. 3 Transient excitation load in time ( a) and frequency ( b) domains 1. 3 有限元计算方法有效性验证 采用有限元方法模拟无缺陷金属板上表面反射 回波,并与实验测量结果进行对比以检验本文采用 有限元方法模拟金属内超声波传播的有效性. 金属 板上表面所有节点 y 方向位移的平均值,如图 4 所 示. 图 4 显示,第 1 次底面回波到达上表面的时间 t = 1. 80 μs,相邻两次底面回波之间的时间间隔也为 t = 1. 80 μs,与理论计算值 t = 1. 75 μs 基本一致. 由 于有限元计算采用显示中心差分算法,在数值计算 过程中存在舍入误差和截断误差,从而引起计算中 的数值噪声[12]. 为了验证有限元计算结果的准确性,采用 IOS AIR--1550 TWM 激光超声检测设备测量了无缺陷钢 板的超声波信号. 激光发生器的波长为 532 nm,发 射的激光垂直作用在被测钢板表面,产生的超声信 号被带宽为 125 MHz 的激光超声接收仪接收. 测量 钢板为矩形断面,厚度 H = 5 mm,宽度 B = 20 mm. 采样点数为 500,采样率为 0. 02 μs,信号总长度为 10 μs. 测量获得的 1 ~ 8 μs 时间内钢板上表面的 y 方向位移响应,如图 4 所示. 仿真结果和实验结果 的各底面回波到达时间基本完全一致,验证了本文 有限元方法的准确性. 但是,由于仿真计算忽略了 实际材料内部晶界多重背散射引起的衰减,因此仿 真计算结果的回波幅值衰减要小于实测结果. 同 时,实验测试的噪声和测试样品的有限宽度也是引 起实验测试结果和仿真结果差异的原因. 图 4 厚度为 5 mm 无缺陷金属板 y 方向波形的仿真结果和实测 结果对比 Fig. 4 Comparison of the numerical and experimental reflected y-displacement in a defect-free plate with the thickness of 5 mm 2 计算结果与分析 2. 1 内部无裂纹 对图 4 所示的无裂纹金属板上表面仿真曲线的 四次底面回波进行快速傅里叶变换后的频域曲线如 图 5 所示. 图 5 显示,对入射中心频率为 10 MHz 的 载荷,随着超声波在金属板内传播,底面回波的频率 范围扩大,不同底面回波的频率范围均在 0 ~ 20 MHz. 图 5 同时显示,随着回波次数的增加,不同频 ·840·
第6期 肖会芳等:金属内部裂纹的超声衰减特性与近表面缺陷识别 ·841· 率成分的幅值均有所减小,即存在衰减.由于金属 2.2内部有裂纹 板为各向同性且无缺陷,该衰减主要由于几何扩散 为了获得裂纹深度对超声波的反射、散射以及 引起团 衰减的影响规律,通过建立不同裂纹深度h的有限 0.08 元仿真模型,计算获得反射回波信号,并提取其衰减 1次回波 0.07 2次回波 系数随频率的变化关系.其中,裂纹的尺寸为!= 3次回波 0.06 4次同波 0.5mm,d=0.1mm. 0.05 2.2.1时域波形 0.04 裂纹深度h=0.25,1,2,4.75mm时,金属板 0.03 上表面接收的反射回波的时域波形,如图7所示. 0.02 图7(b)和(c)显示,对缺陷深度h=1mm和h=2 mm,在始波Fn和第1次底面回波(14B)之间, 0.01 以及相邻的底面回波1Bo2Bho3Bh和4 10 15 20 25 30 f/MHz B之间,可以观察到明显的多次缺陷回波信号I Cn2C3Cn和4Co·对低阶底面回波, 图5内部无裂纹金属板的不同阶次底面回波的频谱图 回波信号更为集中且幅值较大.随着回波次数增 Fig.5 Spectra of different backwall echoes for the defect-free plate 加,底面回波信号的幅值逐渐减小而时间跨度逐渐 超声波在金属板内传播的衰减特性随频率的变 增大,即第4次底面回波时间跨度T,最大,第1次 化关系,可以采用下式所示的衰减系数描述网: 底面回波时间跨度T最小,表明经过多次反射后, 1 1Fa+1(w) 超声波发生了衰减和散射. a=-2hn F.(o) (1) 多次裂纹回波信号具有与底面回波信号类似的 式中:α为衰减系数,定义为单位距离上振幅比值的 特征,即较早的裂纹回波信号(1C)更为集中, 自然对数,Np·cm;H为金属板厚;Fn+1(o)和 随着回波次数增加(2Co、3Ceh.和4曲Cn),可 F,(w)分别为相邻的底面回波在频域内的幅值,n= 以观察到明显的衰减和散射.图7(b)和(c)同时显 1,2,3.对四次底面回波,可以计算获得三组衰 示,随着裂纹深度h增加,入射波传播到裂纹表面的 减系数,对三组衰减系数进行平均处理获得的平 时间增加,相邻底面回波之间的缺陷回波次数减少, 均衰减系数随频率的变化关系曲线,如图6所 然而,对裂纹深度h=0.25mm,裂纹缺陷非常 示.图6显示,由几何扩散引起的衰减系数较小, 靠近上表面,始波和缺陷回波之间的时间间隔小于 最大值仅为0.075Np·cm-1,对应的频率为f= 始波的时间宽度,缺陷回波与始波发生混叠,导致无 18.6MHz.同时可以看到,扩散衰减系数随频率 法识别裂纹回波,如图7(a)所示.对裂纹深度h= 的增大逐渐增加.在计算由裂纹引起的衰减时, 4.75mm,裂纹缺陷非常靠近底面,由于底面回波和 需要将由几何扩散引起的衰减量从计算获得的 缺陷回波之间的时间间隔小于底面回波的时间宽 总衰减量中减除 度,缺陷回波与底面回波发生混叠,如图7(d)所示. 0.08 对采用同端激励、同端接收的脉冲回波法,这些靠近 0.07 检测物体表面的区域为检测盲区. 2.2.2波场分布 E006 当裂纹深度为h=1mm时,不同时刻t=0.30, 0.45,1.15,2.00μs,材料内部的波场分布,如图8 0.04 所示.从图8可以清楚地观察波在传播过程中,在 裂纹表面的反射、散射和波型转换特征.图8(a)显 0.02 示,产生的纵波沿着厚度方向传播,向底面传播的波 0.015 4 和从裂纹表面反射的波在t=0.30μs相遇.当波到 68101214161820 f/MHz 达裂纹表面后,以裂纹的尖端和表面作为新的波源, 图6内部无裂纹金属板的衰减系数随频率的变化关系曲线 产生不同的散射波,如图8(b)所示.一部分透射波 Fig.6 Plot of attenuation coefficient vs.frequency for the defect-free W。,继续向底面传播,而另一部分反射波W,则向 plate 材料上表面传播,在上表面和裂纹表面之间发生多
第 6 期 肖会芳等: 金属内部裂纹的超声衰减特性与近表面缺陷识别 率成分的幅值均有所减小,即存在衰减. 由于金属 板为各向同性且无缺陷,该衰减主要由于几何扩散 引起[13]. 图 5 内部无裂纹金属板的不同阶次底面回波的频谱图 Fig. 5 Spectra of different backwall echoes for the defect-free plate 超声波在金属板内传播的衰减特性随频率的变 化关系,可以采用下式所示的衰减系数描述[14]: α = - 1 2Hln Fn + 1 ( ω) Fn ( ω) . ( 1) 图 6 内部无裂纹金属板的衰减系数随频率的变化关系曲线 Fig. 6 Plot of attenuation coefficient vs. frequency for the defect-free plate 式中: α 为衰减系数,定义为单位距离上振幅比值的 自然 对 数,Np·cm - 1 ; H 为 金 属 板 厚; Fn + 1 ( ω) 和 Fn ( ω) 分别为相邻的底面回波在频域内的幅值,n = 1,2,3. 对四次底面回波,可以计算获得三组衰 减系数,对三组衰减系数进行平均处理获得的平 均衰减系数随频率的变化关系曲线,如 图 6 所 示. 图 6 显示,由几何扩散引起的衰减系数较小, 最大值 仅 为 0. 075 Np·cm - 1 ,对 应 的 频 率 为 f = 18. 6 MHz. 同时可以看到,扩散衰减系数随频率 的增大逐渐增加. 在计算由裂纹引起的衰减时, 需要将由几何扩散引起的衰减量从计算获得 的 总衰减量中减除. 2. 2 内部有裂纹 为了获得裂纹深度对超声波的反射、散射以及 衰减的影响规律,通过建立不同裂纹深度 h 的有限 元仿真模型,计算获得反射回波信号,并提取其衰减 系数随频率的变化关系. 其中,裂纹的尺寸为 l = 0. 5 mm,d = 0. 1 mm. 2. 2. 1 时域波形 裂纹深度 h = 0. 25,1,2,4. 75 mm 时,金属板 上表面接收的反射回波的时域波形,如图 7 所示. 图 7( b) 和( c) 显示,对缺陷深度 h = 1 mm 和 h = 2 mm,在始波 Fecho和第 1 次底面回波( 1st Becho ) 之间, 以及相邻的底面回波 1st Becho、2nd Becho、3rd Becho和 4th Becho之间,可以观察到明显的多次缺陷回波信号 1st Cecho、2nd Cecho、3rd Cecho和 4th Cecho . 对低阶底面回波, 回波信号更为集中且幅值较大. 随着回波次数增 加,底面回波信号的幅值逐渐减小而时间跨度逐渐 增大,即第 4 次底面回波时间跨度 T4 最大,第 1 次 底面回波时间跨度 T1 最小,表明经过多次反射后, 超声波发生了衰减和散射. 多次裂纹回波信号具有与底面回波信号类似的 特征,即较早的裂纹回波信号( 1st Cecho ) 更为集中, 随着回波次数增加( 2nd Cecho、3rd Cecho和 4th Cecho ) ,可 以观察到明显的衰减和散射. 图 7( b) 和( c) 同时显 示,随着裂纹深度 h 增加,入射波传播到裂纹表面的 时间增加,相邻底面回波之间的缺陷回波次数减少. 然而,对裂纹深度 h = 0. 25 mm,裂纹缺陷非常 靠近上表面,始波和缺陷回波之间的时间间隔小于 始波的时间宽度,缺陷回波与始波发生混叠,导致无 法识别裂纹回波,如图 7( a) 所示. 对裂纹深度 h = 4. 75 mm,裂纹缺陷非常靠近底面,由于底面回波和 缺陷回波之间的时间间隔小于底面回波的时间宽 度,缺陷回波与底面回波发生混叠,如图 7( d) 所示. 对采用同端激励、同端接收的脉冲回波法,这些靠近 检测物体表面的区域为检测盲区. 2. 2. 2 波场分布 当裂纹深度为 h = 1 mm 时,不同时刻 t = 0. 30, 0. 45,1. 15,2. 00 μs,材料内部的波场分布,如图 8 所示. 从图 8 可以清楚地观察波在传播过程中,在 裂纹表面的反射、散射和波型转换特征. 图 8( a) 显 示,产生的纵波沿着厚度方向传播,向底面传播的波 和从裂纹表面反射的波在 t = 0. 30 μs 相遇. 当波到 达裂纹表面后,以裂纹的尖端和表面作为新的波源, 产生不同的散射波,如图 8( b) 所示. 一部分透射波 Wp,继续向底面传播,而另一部分反射波 Wr,则向 材料上表面传播,在上表面和裂纹表面之间发生多 ·841·
·842· 北京科技大学学报 第36卷 0.04 0.04r (a) LT 0.02 0.02 4C n i -0.02 -02 2. -0.04 0.04 0.060 0060 3 0.04 0.06 0.02 0.04 0.02 -0.02 -0.02 0.04 0.04 0.06 0.06 0.0801234567 0.08 2345678 图7不同裂纹深度时金属板上表面接收的反射回波波形.(a)h=0.25mm:(b)h=1mm:(c)h=2mm:(d)h=4.75mm Fig.7 Reflected y-displacement for four different crack depths:(a)h=0.25 mm;(b)h=I mm;(e)h=2 mm;(d)h=4.75 mm d 图8裂纹深度为h=1mm时不同时刻的波场图.(a)t=0,30us:(b)t=0.455:(c)1=1.15s:(d)=2.00μ5 Fig.8 Propagation of ultrasonic waves for the crack depth of I mm at different time steps:(a)t=0.30 ps;(b)t=0.45 us;(e)=1.15 ps:(d) 1=2.00us
北 京 科 技 大 学 学 报 第 36 卷 图 7 不同裂纹深度时金属板上表面接收的反射回波波形. ( a) h = 0. 25 mm; ( b) h = 1 mm; ( c) h = 2mm; ( d) h = 4. 75 mm Fig. 7 Reflected y-displacement for four different crack depths: ( a) h = 0. 25 mm; ( b) h = 1 mm; ( c) h = 2 mm; ( d) h = 4. 75 mm 图 8 裂纹深度为 h = 1 mm 时不同时刻的波场图. ( a) t = 0. 30 μs; ( b) t = 0. 45 μs; ( c) t = 1. 15 μs; ( d) t = 2. 00 μs Fig. 8 Propagation of ultrasonic waves for the crack depth of 1 mm at different time steps: ( a) t = 0. 30 μs; ( b) t = 0. 45 μs; ( c) t = 1. 15 μs; ( d) t = 2. 00 μs ·842·
第6期 肖会芳等:金属内部裂纹的超声衰减特性与近表面缺陷识别 ·843· 次反射,被探头接收,如图7(b)所示的多次裂纹回 幅值均减小,且衰减量远大于图5所示的几何散射 波信号.图8(d)显示,第1次底面回波在t=2.00 衰减量. μs到达上表面.之后,底面回波从上表面反射,并 图9(a)和(d)显示,由于裂纹回波完全淹没在 以第2次入射波再次进入材料内部,继续向底面传 始波或底面回波中,高次回波的衰减量较小,尤其是 播.由于第2次入射波、第1次传播波的波场和裂 在中心频率附近.同时,裂纹对回波信号的调制作 纹散射波的相互作用,波场更为复杂 用较弱,回波信号的频谱图与无缺陷回波信号相似 2.2.3衰减系数 (图5),不同之处在于中心频率附近的带宽变窄. 对图7所示的不同裂纹深度金属板上表面接收 与图9(a)和(d)所对应的裂纹回波完全淹没在始 的反射回波中的四次底面回波进行快速傅里叶变换 波或底面回波中的情况相比,图9(b)和(c)显示, 后的频域曲线,如图9所示.与无裂纹类似,由于超 由于裂纹回波与表面回波分离,缺陷引起的波场与 声波的传播,底面回波的频率成分增加,不同底面回 底面回波发生混叠,对底面回波信号产生明显的调 波的频率范围均在0~20MHz.同时,由于存在内部 制作用,高次回波的频率衰减现象严重,且衰减频率 裂纹,随着回波次数的增加,不同频率成分超声波的 发生在低频部分 0.06 0.06 (ay 1次回波 1次回波 0.05 2次回波 0.05 一2次回波 3次回波 3次回波 0.04 0.04 4次回波 4次回波 号0.03 0.02 0.02 0.01 00 10 15 20 25 30 10 15 20 25 30 fMMHz f/MHz 0.06 0.08m (e) 1次回波 d 1次回波 0.05 一一2次回波 0.07 一2次回波 3次回波 0.06 3次回波 1.04 4次回波 0.05 4次回波 0.03 0.04 0.03 00 0.02 0.01 0.01 10 15 20 30 00 10 15 20 25 30 f/MHz MHz 图9不同裂纹深度时底面回波的频谱图.(a)h=0.25mm:(b)h=1mm:(c)h=2mm:(d)h=4.75mm Fig.9 Spectra of the four backwall echoes for different crack depths:(a)h=0.25 mm:(b)h=1 mm;(c)h=2 mm:(d)h=4.75 mm 图9(a)和(d)所示的频域特征可用于辨识材 化关系曲线,如图10所示.其中,图10(a)的裂纹 料内部是否存在近表面缺陷:首先,通过脉冲回波法 深度为h=0.25mm和h=4.75mm,对应于裂纹回 测试无缺陷样品的表面回波信号,并通过快速傅里 波淹没在表面回波中;图10(b)的裂纹深度为h= 叶变换获得其底面回波的频谱曲线,作为参考曲线 1mm和h=2mm,对应于裂纹回波与表面回波分 (如图5);对被检测样品,同样获得其底面回波的频 离.图10显示,由裂纹缺陷导致的衰减量远大于 谱曲线,与参考曲线进行对比.若其中心频率附近 由几何扩散引起的衰减量(如图6所示),但在整 的带宽变窄,则其近表面可能存在缺陷 个频域范围内,其高频衰减量的相对幅值却小于 采用式(1)计算不同裂纹深度时的衰减系数, 几何扩散 并将几何扩散引起的衰减量从计算获得的总衰减量 图10(a)显示,当裂纹回波淹没在表面回波中 中减除,获得不同裂纹深度的衰减系数随频率的变 时,对不同的裂纹深度,衰减系数随频率的变化关系
第 6 期 肖会芳等: 金属内部裂纹的超声衰减特性与近表面缺陷识别 次反射,被探头接收,如图 7( b) 所示的多次裂纹回 波信号. 图 8( d) 显示,第 1 次底面回波在 t = 2. 00 μs 到达上表面. 之后,底面回波从上表面反射,并 以第 2 次入射波再次进入材料内部,继续向底面传 播. 由于第 2 次入射波、第 1 次传播波的波场和裂 纹散射波的相互作用,波场更为复杂. 2. 2. 3 衰减系数 对图 7 所示的不同裂纹深度金属板上表面接收 的反射回波中的四次底面回波进行快速傅里叶变换 后的频域曲线,如图 9 所示. 与无裂纹类似,由于超 声波的传播,底面回波的频率成分增加,不同底面回 波的频率范围均在 0 ~ 20 MHz. 同时,由于存在内部 裂纹,随着回波次数的增加,不同频率成分超声波的 幅值均减小,且衰减量远大于图 5 所示的几何散射 衰减量. 图 9( a) 和( d) 显示,由于裂纹回波完全淹没在 始波或底面回波中,高次回波的衰减量较小,尤其是 在中心频率附近. 同时,裂纹对回波信号的调制作 用较弱,回波信号的频谱图与无缺陷回波信号相似 ( 图 5) ,不同之处在于中心频率附近的带宽变窄. 与图 9( a) 和( d) 所对应的裂纹回波完全淹没在始 波或底面回波中的情况相比,图 9( b) 和( c) 显示, 由于裂纹回波与表面回波分离,缺陷引起的波场与 底面回波发生混叠,对底面回波信号产生明显的调 制作用,高次回波的频率衰减现象严重,且衰减频率 发生在低频部分. 图 9 不同裂纹深度时底面回波的频谱图. ( a) h = 0. 25 mm; ( b) h = 1 mm; ( c) h = 2 mm; ( d) h = 4. 75 mm Fig. 9 Spectra of the four backwall echoes for different crack depths: ( a) h = 0. 25 mm; ( b) h = 1 mm; ( c) h = 2 mm; ( d) h = 4. 75 mm 图 9( a) 和( d) 所示的频域特征可用于辨识材 料内部是否存在近表面缺陷: 首先,通过脉冲回波法 测试无缺陷样品的表面回波信号,并通过快速傅里 叶变换获得其底面回波的频谱曲线,作为参考曲线 ( 如图 5) ; 对被检测样品,同样获得其底面回波的频 谱曲线,与参考曲线进行对比. 若其中心频率附近 的带宽变窄,则其近表面可能存在缺陷. 采用式( 1) 计算不同裂纹深度时的衰减系数, 并将几何扩散引起的衰减量从计算获得的总衰减量 中减除,获得不同裂纹深度的衰减系数随频率的变 化关系曲线,如图 10 所示. 其中,图 10( a) 的裂纹 深度为 h = 0. 25 mm 和 h = 4. 75 mm,对应于裂纹回 波淹没在表面回波中; 图 10( b) 的裂纹深度为 h = 1 mm 和 h = 2 mm,对应于裂纹回波与表面回波分 离. 图 10 显示,由裂纹缺陷导致的衰减量远大于 由几何扩散引起的衰减量( 如图 6 所示) ,但在整 个频域范围内,其高频衰减量的相对幅值却小于 几何扩散. 图 10( a) 显示,当裂纹回波淹没在表面回波中 时,对不同的裂纹深度,衰减系数随频率的变化关系 ·843·
·844· 北京科技大学学报 第36卷 基本完全一致:最小衰减系数发生在中心频率附近, 由于裂纹回波信号对底面回波信号的频谱调制作用 对应于最小衰减量;随着频率成分偏离中心频率,衰 较严重,衰减系数随频率的变化关系曲线存在多个 减量增大,最大衰减量发生在高频处,对应的频率约 峰值与谷值,最大的衰减系数仍出现在高频处,对不 为f=16MHz. 同的裂纹深度h=1mm和h=2mm,最大衰减量对 图10(b)显示,当裂纹回波与表面回波分离时, 应的频率约为f=15MHz和f=16.5MHz. 0.8 0.45 b 0.7 --0.25mm 0.40 一h=lmm -h=4.75mm -h=2 mm 0.6 0.35 0.5 0.30 0.4 03 0.25 0.2 0.20 01 0.15 0.10 6810121416 6810121416 f/MHz f/MHz 图10不同裂纹深度时衰减系数随频率的变化关系曲线 Fig.10 Plots of attenuation coefficient vs.frequency for different crack depths for interal flaws in as-hot-rolled steel strip using ultrasonic probe 3结论 array.Mater Trans,2011,52:531 (1)不同深度裂纹对底面回波信号的调制作用 5]Chi DZ,Gang T,Yao Y X,et al.Mode for near surface defect 不同,随着深度增大,调制作用增强,但是近表面裂 detection based on ultrasonic TOFD.Trans China Weld Inst, 2011,32(2):25 纹对底面回波信号的调制作用较弱. (迟大钊,刚铁,姚英学,等。一种基于超声TOD法的近表 (2)对近表面缺陷,与无缺陷回波信号相比,其 面缺陷检测模式.焊接学报,2011,32(2):25) 底面回波信号的频谱图中心频率附近的带宽变窄. 6] Rajagopal P,Lowe M JS.Scattering of the fundamental shear hor- 该频域特征可用于辨识材料内部是否存在近表面 izontal guided wave by a part-thickness crack in an isotropic plate. 缺陷. J Acoust Soc Am,2008,124(5):2895 (3)对近表面缺陷,最小衰减量发生在中心频 ] Tsaur D H.Exact scattering and diffraction of antiplane shear waves by a vertical edge crack.Geophys nt,2010.181:1655 率附近,随着频率成分偏离中心频率,衰减量增大, 8] Kline R A.Nondestructive Characterization of Composite Media. 最大衰减量发生在高频处 New York:Technomic Publishing Company,1992 (4)当裂纹回波与表面回波分离时,衰减系数 9) Yeh H Y,Cheng J H.NDE of metal damage:ultrasonics with a 随频率的变化关系曲线存在多个峰值与谷值,最大 damage mechanics model.Int J Solids Struct.2003,40:7285 的衰减系数仍出现在高频处 [Dattal D,Kishoret N N.Features of ultrasonic wave propagation to identify defects in composite materials modelled by finite ele- 参考文献 ment method.NDT E Int,1996,29:213 [Achenbach J D.Quantitative nondestructive evaluation.Int J Sol- [11]Wang X D,Huang G L.Identification of embedded cracks using ids Struct,2000,37:13 backpropagating elastic waves.Inrerse Probl,2004,20:1393 Wagle S,Kato H.Ultrasonic wave intensity reflected from fretting [12]Ghoshal G,Tumner J A.Numerical model of longitudinal wave fatigue cracks at bolt joints of aluminum alloy plates.NDT E scattering in polyerystals.IEEE Trans Ultrason Ferroelectr Freq mt,2009,42:690 Control,2009,56(7):1419 B3]Wang X G,Chang J J,Shan YC,et al.Testing of characteristic [13]Liu D L,Tumer J A.Ultrasonic backscatter in two-phase media of rubber thin layer by ultrasonic echo signal.Chin Mech Eng, and its dependency on the correlation function.Proc Meet Acoust, 2008,44(10):114 2012:045004 (王兴国,常俊杰,单英春,等.超声回波信号检测橡胶薄层 [14]Islam M N,Arai Y.Ultrasonic back reflection evaluation of crack 的特性.机械工程学报,2008,44(10):114) growth from PSBs in loweycle fatigue of stainless steel under con 4]Takada H,Tomura Y,Aratani M,et al.On-ine detection system stant load amplitude.Mater Sci Eng A,2009,520:49
北 京 科 技 大 学 学 报 第 36 卷 基本完全一致: 最小衰减系数发生在中心频率附近, 对应于最小衰减量; 随着频率成分偏离中心频率,衰 减量增大,最大衰减量发生在高频处,对应的频率约 为 f = 16 MHz. 图 10( b) 显示,当裂纹回波与表面回波分离时, 由于裂纹回波信号对底面回波信号的频谱调制作用 较严重,衰减系数随频率的变化关系曲线存在多个 峰值与谷值,最大的衰减系数仍出现在高频处,对不 同的裂纹深度 h = 1 mm 和 h = 2 mm,最大衰减量对 应的频率约为 f = 15 MHz 和 f = 16. 5 MHz. 图 10 不同裂纹深度时衰减系数随频率的变化关系曲线 Fig. 10 Plots of attenuation coefficient vs. frequency for different crack depths 3 结论 ( 1) 不同深度裂纹对底面回波信号的调制作用 不同,随着深度增大,调制作用增强,但是近表面裂 纹对底面回波信号的调制作用较弱. ( 2) 对近表面缺陷,与无缺陷回波信号相比,其 底面回波信号的频谱图中心频率附近的带宽变窄. 该频域特征可用于辨识材料内部是否存在近表面 缺陷. ( 3) 对近表面缺陷,最小衰减量发生在中心频 率附近,随着频率成分偏离中心频率,衰减量增大, 最大衰减量发生在高频处. ( 4) 当裂纹回波与表面回波分离时,衰减系数 随频率的变化关系曲线存在多个峰值与谷值,最大 的衰减系数仍出现在高频处. 参 考 文 献 [1] Achenbach J D. Quantitative nondestructive evaluation. Int J Solids Struct,2000,37: 13 [2] Wagle S,Kato H. Ultrasonic wave intensity reflected from fretting fatigue cracks at bolt joints of aluminum alloy plates. NDT & E Int,2009,42: 690 [3] Wang X G,Chang J J,Shan Y C,et al. Testing of characteristic of rubber thin layer by ultrasonic echo signal. Chin J Mech Eng, 2008,44( 10) : 114 ( 王兴国,常俊杰,单英春,等. 超声回波信号检测橡胶薄层 的特性. 机械工程学报,2008,44( 10) : 114) [4] Takada H,Tomura Y,Aratani M,et al. On-line detection system for internal flaws in as-hot-rolled steel strip using ultrasonic probe array. Mater Trans,2011,52: 531 [5] Chi D Z,Gang T,Yao Y X,et al. Mode for near surface defect detection based on ultrasonic TOFD. Trans China Weld Inst, 2011,32( 2) : 25 ( 迟大钊,刚铁,姚英学,等. 一种基于超声 TOFD 法的近表 面缺陷检测模式. 焊接学报,2011,32( 2) : 25) [6] Rajagopal P,Lowe M J S. Scattering of the fundamental shear horizontal guided wave by a part-thickness crack in an isotropic plate. J Acoust Soc Am,2008,124( 5) : 2895 [7] Tsaur D H. Exact scattering and diffraction of antiplane shear waves by a vertical edge crack. Geophys J Int,2010,181: 1655 [8] Kline R A. Nondestructive Characterization of Composite Media. New York: Technomic Publishing Company,1992 [9] Yeh H Y,Cheng J H. NDE of metal damage: ultrasonics with a damage mechanics model. Int J Solids Struct,2003,40: 7285 [10] Dattal D,Kishoret N N. Features of ultrasonic wave propagation to identify defects in composite materials modelled by finite element method. NDT & E Int,1996,29: 213 [11] Wang X D,Huang G L. Identification of embedded cracks using back-propagating elastic waves. Inverse Probl,2004,20: 1393 [12] Ghoshal G,Turner J A. Numerical model of longitudinal wave scattering in polycrystals. IEEE Trans Ultrason Ferroelectr Freq Control,2009,56( 7) : 1419 [13] Liu D L,Turner J A. Ultrasonic backscatter in two-phase media and its dependency on the correlation function. Proc Meet Acoust, 2012: 045004 [14] Islam M N,Arai Y. Ultrasonic back reflection evaluation of crack growth from PSBs in low-cycle fatigue of stainless steel under constant load amplitude. Mater Sci Eng A,2009,520: 49 ·844·