工程科学学报，第37卷，第7期：942949,2015年7月 Chinese Journal of Engineering,Vol.37,No.7:942-949,July 2015 DOI:10.13374/j.issn2095-9389.2015.07.017:http://journals.ustb.edu.cn 金属材料内部非金属夹杂超声检测的数值模拟 陈 丹”，肖会芳)区，黎敏”，王善超”，徐金梧” 1)北京科技大学机械工程学院，北京1000832)北京科技大学国家板带生产先进装备工程技术研究中心，北京100083 ☒通信作者，E-mail:huifangxiao(@usth.edu.cm 摘要在金属材料内部夹杂物的超声检测中，如何通过检测获得的回波信号辨识夹杂物的属性和位置，一直是其重点和难 点问题.通过建立包含夹杂物缺陷的二维金属板模型，采用有限元数值模拟的方法，对材料内部超声波场进行计算，获得了 两种最典型的夹杂物A山，0，和TN,以及二者在材料内部不同深度时的超声回波信号.研究了夹杂物类型和夹杂物深度对超 声回波时域波形以及对界面波、夹杂物缺陷回波和底面回波频谱分布的影响规律 关键词金属材料：非金属夹杂物：超声检测：数值模拟：频谱分析 分类号TG115.285:G142.1 Numerical simulation of ultrasonic testing for non-metallic inclusions in metallic materials CHEN Dan,XIAO Hui-fang?,LI Min,WANG Shan-chao),XU Jin-wu) 1)School of Mechanical Engineering,University of Science and Technology Beijing,Beijing 100083,China 2)National Engineering Research Center of Flat Rolling Equipment,University of Science and Technology Beijing,Beijing 100083 Corresponding author,E-mail:huifangxiao@ustb.edu.cn ABSTRACT In the area of internal inclusion detection using ultrasonic,identifying the location and type of internal inclusions using the collected back echo signals has always been a difficult problem.A two-dimensional metallic plate containing inclusions is estab- lished.Ultrasonic wave fields inside this material are calculated by the finite element method.Echo waves for two typical inclusions of Al,O,and TiN and for these inclusions with different depths are obtained.The influences of the type and depth of inclusions on the time-domain and the spectral distribution of interfacial waves,echo waves of inclusions and bottom echo waves are analyzed especially. KEY WORDS metallic materials:non-metallic inclusions:ultrasonic testing:numerical simulation:spectral analysis 钢中非金属夹杂物是指钢中不具有金属性质的氧 夹杂物上形核，并进一步聚合、长大和扩展，最终导致 化物、硫化物、硅酸盐或氮化物.它们是钢在治炼过程 工程结构材料断裂.夹杂物数量越多，尺寸越大，裂纹 中因加入脱氧剂而形成的氧化物、硅酸盐和钢在凝固 萌生的几率越高，钢材的韧性值就越低.大尺寸夹杂 过程中由于某些元素溶解度下降而形成的硫化物和氮 物对钢疲劳强度的有害性更强，而小于临界尺寸(4~ 化物，以及炉渣或耐火材料来不及排除而留在钢中形 7.5μm)的夹杂物往往无害甚至有益四.随着洁净钢 成的夹杂物四.非金属夹杂物对钢的使用性能和工艺 的发展，控制钢中夹杂物的存在已成为钢铁生产企业 性能均产生严重影响，包括强度、塑性、断裂韧性、切 技术发展战略的关键.为了判断工艺条件对钢中夹杂 削、疲劳、热脆、耐腐蚀性等.由于非金属夹杂物的存 物类型和位置的影响，需要在炼钢过程中建立一套评 在破坏了金属的连续性，造成应力集中，引起微裂纹在 价钢中夹杂物水平的在线检测方法，找出钢水中夹杂 收稿日期：2014-04-20 基金项目：“十二五”国家科技支撑计划资助项目(2012BAF04B02):北京高等学校青年英才计划资助项目(YET0373)
工程科学学报，第 37 卷，第 7 期: 942--949，2015 年 7 月 Chinese Journal of Engineering，Vol． 37，No． 7: 942--949，July 2015 DOI: 10． 13374 /j． issn2095--9389． 2015． 07． 017; http: / /journals． ustb． edu． cn 金属材料内部非金属夹杂超声检测的数值模拟 陈 丹1) ，肖会芳2) ，黎 敏1) ，王善超1) ，徐金梧1) 1) 北京科技大学机械工程学院，北京 100083 2) 北京科技大学国家板带生产先进装备工程技术研究中心，北京 100083 通信作者，E-mail: huifangxiao@ ustb． edu． cn 摘 要 在金属材料内部夹杂物的超声检测中，如何通过检测获得的回波信号辨识夹杂物的属性和位置，一直是其重点和难 点问题． 通过建立包含夹杂物缺陷的二维金属板模型，采用有限元数值模拟的方法，对材料内部超声波场进行计算，获得了 两种最典型的夹杂物 Al2O3和 TiN，以及二者在材料内部不同深度时的超声回波信号． 研究了夹杂物类型和夹杂物深度对超 声回波时域波形以及对界面波、夹杂物缺陷回波和底面回波频谱分布的影响规律． 关键词 金属材料; 非金属夹杂物; 超声检测; 数值模拟; 频谱分析 分类号 TG115. 285; G142. 1 Numerical simulation of ultrasonic testing for non-metallic inclusions in metallic materials CHEN Dan1) ，XIAO Hui-fang2) ，LI Min1) ，WANG Shan-chao1) ，XU Jin-wu1) 1) School of Mechanical Engineering，University of Science and Technology Beijing，Beijing 100083，China 2) National Engineering Ｒesearch Center of Flat Ｒolling Equipment，University of Science and Technology Beijing，Beijing 100083 Corresponding author，E-mail: huifangxiao@ ustb． edu． cn ABSTＲACT In the area of internal inclusion detection using ultrasonic，identifying the location and type of internal inclusions using the collected back echo signals has always been a difficult problem． A two-dimensional metallic plate containing inclusions is established． Ultrasonic wave fields inside this material are calculated by the finite element method． Echo waves for two typical inclusions of Al2O3 and TiN and for these inclusions with different depths are obtained． The influences of the type and depth of inclusions on the time-domain and the spectral distribution of interfacial waves，echo waves of inclusions and bottom echo waves are analyzed especially． KEY WOＲDS metallic materials; non-metallic inclusions; ultrasonic testing; numerical simulation; spectral analysis 收稿日期: 2014--04--20 基金项目: “十二五”国家科技支撑计划资助项目( 2012BAF04B02) ; 北京高等学校青年英才计划资助项目( YETP0373) 钢中非金属夹杂物是指钢中不具有金属性质的氧 化物、硫化物、硅酸盐或氮化物． 它们是钢在冶炼过程 中因加入脱氧剂而形成的氧化物、硅酸盐和钢在凝固 过程中由于某些元素溶解度下降而形成的硫化物和氮 化物，以及炉渣或耐火材料来不及排除而留在钢中形 成的夹杂物［1］． 非金属夹杂物对钢的使用性能和工艺 性能均产生严重影响，包括强度、塑性、断裂韧性、切 削、疲劳、热脆、耐腐蚀性等． 由于非金属夹杂物的存 在破坏了金属的连续性，造成应力集中，引起微裂纹在 夹杂物上形核，并进一步聚合、长大和扩展，最终导致 工程结构材料断裂． 夹杂物数量越多，尺寸越大，裂纹 萌生的几率越高，钢材的韧性值就越低． 大尺寸夹杂 物对钢疲劳强度的有害性更强，而小于临界尺寸( 4 ～ 7. 5 μm) 的夹杂物往往无害甚至有益［2］． 随着洁净钢 的发展，控制钢中夹杂物的存在已成为钢铁生产企业 技术发展战略的关键． 为了判断工艺条件对钢中夹杂 物类型和位置的影响，需要在炼钢过程中建立一套评 价钢中夹杂物水平的在线检测方法，找出钢水中夹杂
陈丹等：金属材料内部非金属夹杂超声检测的数值模拟 ·943· 物指标和产品等级的关系，确定最优的炼钢工艺) 本文针对钢板中夹杂缺陷的超声检测问题，建立 常见的夹杂物主要分为内生夹杂物和外生夹杂物，如 包含夹杂物缺陷的二维金属板模型，采用有限元数值 最典型的内生夹杂物Al,0,和TN.不同夹杂物的尺寸 模拟的方法，在设置了符合收敛性要求的网格尺寸和 范围一般从几微米到几百微米，有些显微夹杂物甚至 网格划分情况下对材料内部的超声波场进行计算，获 小于1m.与一般的物理、化学分析方法相比，采用超 得了两种典型夹杂物A,0,和TN,以及二者在材料内 声波检测的方法不仅可以实现钢材内部夹杂物的定性 不同深度位置的超声回波信号，研究了夹杂物属性和 和定量检测，而且具有不损伤材料和快速便捷的优点， 夹杂物深度对超声回波时域波形以及对界面波、夹杂 有望实现夹杂物的工业在线检测 物缺陷回波和底面回波频谱分布的影响规律，为材料 然而，对于带钢内部夹杂物的超声检测，如何通过 内部夹杂物缺陷的超声检测提供指导. 检测获得的超声回波信号辨识夹杂物的属性和位置， 1模型描述 一直是夹杂物超声检测的重点和难点问题.与实验测 量相比，数值模拟可以任意改变模型参数以实现不同 1.1几何模型 的检测目标，且仿真结果不受环境噪声、材料表面形貌 由于实际带钢的长度尺寸远大于其宽度和厚度尺 等不确定因素的影响，具有便捷、经济、高效等优点，因 寸，因而对实际的三维钢板模型进行二维的平面应力 此在夹杂物的超声检测方面具有重要的理论和实际意 简化处理，建立宽度和厚度截面的二维金属板模型. 义.Ogilvy通过建立三维有限元模型，给出了夹杂物 研究材料内部夹杂物缺陷对超声波传播特性影响规律 团簇对超声波的透射衰减与背散射信号强度之间的经 的模型示意图见图1.二维矩形金属板的宽度为L,厚 验关系，指出夹杂物的存在会引起检测信号敏感度的 度为H:宽度为l,厚度为h的矩形夹杂物缺陷位于金 下降.Darmon等因模拟并预测了钢板中铝夹杂物的 属板的半宽位置，且夹杂物方向平行于金属板表面，夹 超声响应问题，采用修正后的波恩近似方法可以提高 杂物在金属板内部的位置仅由其深度d决定，定义为 实际检测的精度，并适用于不同形状的夹杂物 夹杂物上表面到金属板上表面的距离.由于声波反射 Agg©lism采用数值模拟的方法研究了不均匀弹性体中 法检测采用同端发射同端接收的方式，因此模型中的 夹杂物的尺寸效应对表面波传播特性的影响，指出夹 加载区域与接收区域均位于模型上表面中心区域附 杂物形状对材料中纵波和瑞利波速度的影响高达 近，二维模型的加载区域长度与普通聚焦探头的焦柱 15%以上. 直径相当. 加载区域 声吸收边界 夹杂物 声吸收边界 金属基体 L 图1二维模型示意图 Fig.I Schematic illustration of a two-dimensional model 计算中，钢板的尺寸为H=5.0mm,L=5.0mm 表1仿真模型参数 夹杂物的尺寸为h=0.2mm,l=1.5mm,吸收边界的宽 Table 1 Parameters of the simulation model 度S=2L=10.0mm.钢板基体及两种夹杂物Al,03和 材料类型 密度/(kgm3)杨氏模量/GPa 泊松比 TN的材料参数如表1所示. 45*钢板（基体） 7890 209 0.269 为了便于分析夹杂物的材料属性和夹杂深度对反 TiN(夹杂物) 5440 450 0.220 射回波波形在时域和频域上的影响规律，选取尺寸较大 A山203（夹杂物） 2700 12 0.300 的矩形夹杂物进行数值模拟，其几何参数如表2所示
陈 丹等: 金属材料内部非金属夹杂超声检测的数值模拟 物指标和产品等级的关系，确定最优的炼钢工艺［3］． 常见的夹杂物主要分为内生夹杂物和外生夹杂物，如 最典型的内生夹杂物 Al2O3和 TiN． 不同夹杂物的尺寸 范围一般从几微米到几百微米，有些显微夹杂物甚至 小于 1 μm． 与一般的物理、化学分析方法相比，采用超 声波检测的方法不仅可以实现钢材内部夹杂物的定性 和定量检测，而且具有不损伤材料和快速便捷的优点， 有望实现夹杂物的工业在线检测［4］． 然而，对于带钢内部夹杂物的超声检测，如何通过 检测获得的超声回波信号辨识夹杂物的属性和位置， 一直是夹杂物超声检测的重点和难点问题． 与实验测 量相比，数值模拟可以任意改变模型参数以实现不同 的检测目标，且仿真结果不受环境噪声、材料表面形貌 等不确定因素的影响，具有便捷、经济、高效等优点，因 此在夹杂物的超声检测方面具有重要的理论和实际意 义． Ogilvy［5］通过建立三维有限元模型，给出了夹杂物 团簇对超声波的透射衰减与背散射信号强度之间的经 验关系，指出夹杂物的存在会引起检测信号敏感度的 下降． Darmon 等［6］模拟并预测了钢板中铝夹杂物的 超声响应问题，采用修正后的波恩近似方法可以提高 实际 检 测 的 精 度，并 适 用 于 不 同 形 状 的 夹 杂 物． Aggelis［7］采用数值模拟的方法研究了不均匀弹性体中 夹杂物的尺寸效应对表面波传播特性的影响，指出夹 杂物形状对材料中纵波和瑞利波速度的影响高达 15% 以上． 本文针对钢板中夹杂缺陷的超声检测问题，建立 包含夹杂物缺陷的二维金属板模型，采用有限元数值 模拟的方法，在设置了符合收敛性要求的网格尺寸和 网格划分情况下对材料内部的超声波场进行计算，获 得了两种典型夹杂物 Al2O3和 TiN，以及二者在材料内 不同深度位置的超声回波信号，研究了夹杂物属性和 夹杂物深度对超声回波时域波形以及对界面波、夹杂 物缺陷回波和底面回波频谱分布的影响规律，为材料 内部夹杂物缺陷的超声检测提供指导． 1 模型描述 1. 1 几何模型 由于实际带钢的长度尺寸远大于其宽度和厚度尺 寸，因而对实际的三维钢板模型进行二维的平面应力 简化处理，建立宽度和厚度截面的二维金属板模型． 研究材料内部夹杂物缺陷对超声波传播特性影响规律 的模型示意图见图 1． 二维矩形金属板的宽度为 L，厚 度为 H; 宽度为 l，厚度为 h 的矩形夹杂物缺陷位于金 属板的半宽位置，且夹杂物方向平行于金属板表面，夹 杂物在金属板内部的位置仅由其深度 d 决定，定义为 夹杂物上表面到金属板上表面的距离． 由于声波反射 法检测采用同端发射同端接收的方式，因此模型中的 加载区域与接收区域均位于模型上表面中心区域附 近，二维模型的加载区域长度与普通聚焦探头的焦柱 直径相当． 图 1 二维模型示意图 Fig． 1 Schematic illustration of a two-dimensional model 计算中，钢板的尺寸为 H = 5. 0 mm，L = 5. 0 mm． 夹杂物的尺寸为 h = 0. 2 mm，l = 1. 5 mm，吸收边界的宽 度 S = 2L = 10. 0 mm． 钢板基体及两种夹杂物 Al2O3和 TiN 的材料参数如表 1 所示． 为了便于分析夹杂物的材料属性和夹杂深度对反 射回波波形在时域和频域上的影响规律，选取尺寸较大 的矩形夹杂物进行数值模拟，其几何参数如表 2 所示． 表 1 仿真模型参数 Table 1 Parameters of the simulation model 材料类型 密度/( kg·m － 3 ) 杨氏模量/GPa 泊松比 45# 钢板( 基体) 7890 209 0. 269 TiN( 夹杂物) 5440 450 0. 220 Al2O3 ( 夹杂物) 2700 72 0. 300 · 349 ·
·944· 工程科学学报，第37卷，第7期 表2夹杂物几何参数列表 型的计算量和计算精度，通常二维仿真模型的宽度尺 Table 2 Geometrical parameters of inclusions 寸远小于实际检测样品的宽度，而通过在宽度方向设 深度，d/mm 置吸收边界，从而抑制两侧声波的反射对求解结果的 类型 尺寸/m2 1 2 影响.通过边界处理的方式抑制两侧的反射横波使得 TiN 1.5×0.2 0.3 2.5 4.5 表面节点位移为纯纵波分量，从而保证与实际测试的 AL203 1.5×0.2 0.3 2.5 4.5 一致性.针对上述尺寸的模型，考虑以无限元网格作 为两侧的边界条件并不能完全抑制声波的反射，本文 1.2有限元计算模型 采用尺寸逐渐增加的有限单元作为声吸收边界.该边 采用商用有限元软件ABAQUS求解超声波在金 界区域是测试区域宽度的2倍，这样既可以有效地抑 属材料内的波动方程.对几何模型进行网格划分，使 制声波在两侧边界的反射，同时又可以降低模型的计 整个求解区域离散化，然后通过定义载荷类型和边 算量，提高计算机的运算效率m.该声吸收边界的控 界条件，在给定的分析类型和输出要求下进行数值 制方程具有如下的形式网 求解.在分析步长的定义中，采用隐式动力学求解算 u=0,x1=±a,i=1,2. 法，通过增加单位时间内的增量步数提高结果的输 {[emay是-c 出精度，时间增量为1ns,对应实际测试中1GHz的 (1) 采样频率。 式(1)建立在空间域2={(x12):x1>0}上.式中 在上述仿真模型的建立过程中要重点满足以下三 的正方向指向2域内，x,=±a表示两侧人工边界所 点要求. 在的位置，u表示x,方向的位移，c为纵波波速，a表示 (1)网格划分.网格分主要考虑网格的单元类 波的第j个传播方向与x,方向的夹角，m为所考虑的 型及相应的网格尺寸.由于模型为规则的矩形板，因 总阶数.此外，底面边界采取固定Y方向上的位移分 此采用四节点的矩形应力单元，既可以提高网格划分 量进行约束限制，从而保证底面回波的完全反射 的质量，又可提高计算的求解精度.其次，对于上述 (3)加载方式.加载方式是压电晶片产生的超声 45钢板，当采用中心频率为10MHz的超声波进行数 波在材料表面作用的等效形式.高阻尼的聚焦探头可 值模拟时，材料中的纵波声速为5748m·s·,对应的声 以提高检测的轴向分辨率四，因此本文采用空间上高 波波长为574.8μm.因此网格尺寸定为10um,小于声 斯分布，时间上汉宁窗调制（图2(a)所示）的三周期正 弦波瞬时动态压力载荷，近似模拟换能器声场分布，其 波波长的1/60，即使对大于16MHz的高频成分，上述 尺寸的网格也足以满足计算的收敛性要求s-0 整体表达式为阳： (2)边界条件.实际夹杂物超声检测中，一般采 y=sin(2mfi)e4)' (2) 用收发一体的纵波直探头，受探头结构及耦合剂的影 式中：a代表载荷区域的半宽，等于探头焦柱直径的 响，接收信号为纯纵波.因此，即使检测对象为小尺寸 12:f代表声波的中心频率：1为载荷脉冲的持续时 样品，检测信号也不会受到两侧反射横波影响.对超 间.当a=250μm,l=300ns,f=10MHz时的声波时域 声传播的有限元仿真计算而言，仿真结果是以表面节 波形和频谱分布，分别如图2(a)和2(b)所示.从图中 点的纵向位移来表征实测A扫波形，考虑到有限元模 可知，声波的频带宽度为6MHz,中心频率为10MHz 10 (aj 0.8 (b) 0.6 0.4 0.2 0 02 -04 0.6 -0.8 50 100150200250300 002 68101214161820 时间/ns 率MHx 图2激励载荷特性时域(a)及频域(b) Fig.2 Excitation load in time (a)and frequency domains (b)
工程科学学报，第 37 卷，第 7 期 表 2 夹杂物几何参数列表 Table 2 Geometrical parameters of inclusions 类型 尺寸/m2 深度，d /mm 1 2 3 TiN 1. 5 × 0. 2 0. 3 2. 5 4. 5 Al2O3 1. 5 × 0. 2 0. 3 2. 5 4. 5 1. 2 有限元计算模型 采用商用有限元软件 ABAQUS 求解超声波在金 属材料内的波动方程． 对几何模型进行网格划分，使 整个求解区域离散化，然后通过定义载荷类型和边 界条件，在给定的分析类型和输出要求下进行数值 求解． 在分析步长的定义中，采用隐式动力学求解算 法，通过增加单位时间内的增量步数提高结果的输 出精度，时间增量为 1 ns，对应实际测试中 1 GHz 的 采样频率． 在上述仿真模型的建立过程中要重点满足以下三 点要求． 图 2 激励载荷特性时域( a) 及频域( b) Fig． 2 Excitation load in time ( a) and frequency domains ( b) ( 1) 网格划分． 网格划分主要考虑网格的单元类 型及相应的网格尺寸． 由于模型为规则的矩形板，因 此采用四节点的矩形应力单元，既可以提高网格划分 的质量，又可提高计算的求解精度． 其次，对于上述 45# 钢板，当采用中心频率为 10 MHz 的超声波进行数 值模拟时，材料中的纵波声速为 5748 m·s － 1，对应的声 波波长为574. 8 μm． 因此网格尺寸定为10 μm，小于声 波波长的 1 /60，即使对大于 16 MHz 的高频成分，上述 尺寸的网格也足以满足计算的收敛性要求［8--10］． ( 2) 边界条件． 实际夹杂物超声检测中，一般采 用收发一体的纵波直探头，受探头结构及耦合剂的影 响，接收信号为纯纵波． 因此，即使检测对象为小尺寸 样品，检测信号也不会受到两侧反射横波影响． 对超 声传播的有限元仿真计算而言，仿真结果是以表面节 点的纵向位移来表征实测 A 扫波形，考虑到有限元模 型的计算量和计算精度，通常二维仿真模型的宽度尺 寸远小于实际检测样品的宽度，而通过在宽度方向设 置吸收边界，从而抑制两侧声波的反射对求解结果的 影响． 通过边界处理的方式抑制两侧的反射横波使得 表面节点位移为纯纵波分量，从而保证与实际测试的 一致性． 针对上述尺寸的模型，考虑以无限元网格作 为两侧的边界条件并不能完全抑制声波的反射，本文 采用尺寸逐渐增加的有限单元作为声吸收边界． 该边 界区域是测试区域宽度的 2 倍，这样既可以有效地抑 制声波在两侧边界的反射，同时又可以降低模型的计 算量，提高计算机的运算效率［11］． 该声吸收边界的控 制方程具有如下的形式［12］ { : ∏ m j = [ 1 cos( aij) t － c x ] } 1 ui = 0，x1 = ± a，i = 1，2． ( 1) 式( 1) 建立在空间域 Ω = { ( x1，x2 ) ∶ x1 ＞ 0} 上． 式中 x1 的正方向指向 Ω 域内，x1 = ± a 表示两侧人工边界所 在的位置，ui表示 xi方向的位移，c 为纵波波速，αij表示 波的第 j 个传播方向与 xi方向的夹角，m 为所考虑的 总阶数． 此外，底面边界采取固定 Y 方向上的位移分 量进行约束限制，从而保证底面回波的完全反射． ( 3) 加载方式． 加载方式是压电晶片产生的超声 波在材料表面作用的等效形式． 高阻尼的聚焦探头可 以提高检测的轴向分辨率［13］，因此本文采用空间上高 斯分布，时间上汉宁窗调制( 图 2( a) 所示) 的三周期正 弦波瞬时动态压力载荷，近似模拟换能器声场分布，其 整体表达式为［14］: y = sin( 2πft) e － 4( x / a) 2 ． ( 2) 式中: a 代表载荷区域的半宽，等于探头焦柱直径的 1 /2; f 代表声波的中心频率; t 为载荷脉冲的持续时 间． 当 a = 250 μm，t = 300 ns，f = 10 MHz 时的声波时域 波形和频谱分布，分别如图 2( a) 和 2( b) 所示． 从图中 可知，声波的频带宽度为 6 MHz，中心频率为 10 MHz． · 449 ·
陈丹等：金属材料内部非金属夹杂超声检测的数值模拟 ·945· 1.3有限元计算方法有效性验证 角度向两侧传播，瑞利波在材料表面一定深度向两 分别从材料内部的声场分布和时域波形两方面验 侧传播。由于纵波速度近似为横波波速的2倍，因此 证有限元模型的有效性.对不含夹杂物的金属材料进 图中任意时刻横波的传播距离接近纵波传播距离的 行仿真计算，将得到的声场分布与动态光弹实验得到 12.图3(b)拍摄的是贴在玻璃界面上的方形压电 的声场分布进行对比，其结果如图3所示.图3()中， 晶片所构成的纵波换能器辐射的全脉冲声场.场中 为了与动态光弹实验结果相对应，施加的压力载荷的 为主的是波前为平面的体纵波，光亮最强.由于动态 脉宽为100s.但是，与水浸超声实验相对应，夹杂物 光弹实验采用等幅加载方式，产生的声波将分别以 仿真计算采用的压力载荷脉宽仍然为300s.在整个 加载区域的两侧端点为中心向外传播，因此边缘横 材料内部自下而上依次是体纵波、体横波和瑞利波. 波及边缘纵波较为明显，而对于二维平面模型采用 体纵波沿着加载方向的能量比较集中，表明体纵波主 瞬时的高斯分布载荷，使得声波能量集中于载荷中 要沿该方向传播，体横波在此方向上的能量较低，而在 心，边缘效应不明显.此外在玻璃表面同样存在明显 两侧能量较高，表明体横波以偏离纵波方向一定的 的瑞利波. al (b) 瑞利波 体横设 图3声场分布有限元法(a)及动态光弹法(b)u的 Fig.3 Ultrasonie field distribution of the finite element method (a)and the dynamic photoelastic method (b) 将仿真波形和实测波形进行对比，如图4所示 图4(a)是对厚度H=2.0mm的模型进行仿真计算得 CL= E 1-σ Vp(1+o)(1-2o) (3) 到的节点Y方向速度随时间的变化曲线.图4(b)是 得到45钢板中纵波的理论声速大小为5748ms,对 对相同厚度的钢板样品进行水浸超声实验得到的实测 应的界面波与底面回波之间的时间间隔为695s.式中 A扫波形.图4(a)和图(b)显示，二者的波形基本一 Ep和σ分别为材料的杨氏模量、密度和泊松比.仿 致，界面波与底面回波之间不存在附加波形，表明两侧 真、实测及理论计算结果在时间间隔上的一致性，验证 声吸收边界有效地抑制了边界声波的反射.图4(a) 了本文有限元模型的准确性。但是，由于仿真计算忽 中界面波与底面回波之间的时间间隔为680s,而实 略了实际材料内部晶界多重背散射引起的衰减和实验 验测试结果显示该时间间隔为686s.同时，由无限大 测试的背景噪声，因此仿真计算结果的信噪比要高于 介质中纵波声速的计算公式a 实测波形 25 150 (b 20 界而波 仿真波形 界面波 实测波形 15 100 10 50 底面问波 底面回波 680s 686ns 汤 100 -25 1002003004005006007008009001000 -150 01002003004005006007008009001000 时间/ns 时国ns 图4仿真(a)与实测(b)波形对比 Fig.4 Waveform contrast of simulation (a)and experimental test (b)
陈 丹等: 金属材料内部非金属夹杂超声检测的数值模拟 1. 3 有限元计算方法有效性验证 分别从材料内部的声场分布和时域波形两方面验 证有限元模型的有效性． 对不含夹杂物的金属材料进 行仿真计算，将得到的声场分布与动态光弹实验得到 的声场分布进行对比，其结果如图3 所示． 图3( a) 中， 为了与动态光弹实验结果相对应，施加的压力载荷的 脉宽为 100 ns． 但是，与水浸超声实验相对应，夹杂物 仿真计算采用的压力载荷脉宽仍然为 300 ns． 在整个 材料内部自下而上依次是体纵波、体横波和瑞利波． 体纵波沿着加载方向的能量比较集中，表明体纵波主 要沿该方向传播，体横波在此方向上的能量较低，而在 两侧能量较高，表明体横波以偏离纵波方向一定的 角度向两侧传播，瑞利波在材料表面一定深度向两 侧传播． 由于纵波速度近似为横波波速的 2 倍，因此 图中任意时刻横波的传播距离接近纵波传播距离的 1 /2． 图 3( b) 拍摄的是贴在玻璃界面上的方形压电 晶片所构成的纵波换能器辐射的全脉冲声场． 场中 为主的是波前为平面的体纵波，光亮最强． 由于动态 光弹实验采用等幅加载方式，产生的声波将分别以 加载区域的两侧端点为中心向外传播，因此边缘横 波及边缘纵波较为明显，而对于二维平面模型采用 瞬时的高斯分布载荷，使得声波能量集中于载荷中 心，边缘效应不明显． 此外在玻璃表面同样存在明显 的瑞利波． 图 3 声场分布有限元法( a) 及动态光弹法( b) ［15］ Fig． 3 Ultrasonic field distribution of the finite element method ( a) and the dynamic photoelastic method ( b) ［15］ 将仿真波形和实测波形进行对比，如图 4 所示． 图 4( a) 是对厚度 H = 2. 0 mm 的模型进行仿真计算得 到的节点 Y 方向速度随时间的变化曲线． 图 4( b) 是 对相同厚度的钢板样品进行水浸超声实验得到的实测 A 扫波形． 图 4( a) 和图( b) 显示，二者的波形基本一 致，界面波与底面回波之间不存在附加波形，表明两侧 声吸收边界有效地抑制了边界声波的反射． 图 4( a) 中界面波与底面回波之间的时间间隔为 680 ns，而实 验测试结果显示该时间间隔为 686 ns． 同时，由无限大 介质中纵波声速的计算公式［16］ CL = E ρ · 1 － σ 槡 ( 1 + σ) ( 1 － 2σ) ( 3) 得到 45# 钢板中纵波的理论声速大小为 5748 m·s － 1，对 应的界面波与底面回波之间的时间间隔为695 ns． 式中 E、ρ 和 σ 分别为材料的杨氏模量、密度和泊松比． 仿 真、实测及理论计算结果在时间间隔上的一致性，验证 了本文有限元模型的准确性． 但是，由于仿真计算忽 略了实际材料内部晶界多重背散射引起的衰减和实验 测试的背景噪声，因此仿真计算结果的信噪比要高于 实测波形． 图 4 仿真( a) 与实测( b) 波形对比 Fig． 4 Waveform contrast of simulation ( a) and experimental test ( b) · 549 ·
·946· 工程科学学报，第37卷，第7期 2 结果与讨论 底面回波之间的A扫波形，如图5所示.其中，夹杂物 距离样品上表面的距离均为d=2.5mm.图中横坐标 采用上述有限元建模方法，建立具有不同类型夹 表示声波在材料中沿厚度方向的传播时间，纵坐标表 杂物和不同深度夹杂物的仿真模型，计算获得回波信 示载荷区域节点的y方向速度分量.图5(a)为两种夹 号，从时域和频域分析夹杂物对超声回波信号的影响 杂物波形的对比结果，其中蓝色实线表示AL,0,夹杂 规律. 物，红色虚线表示TN夹杂物：图5(b)为不含夹杂物 2.1不同类型的夹杂物 的A扫波形.图5(a)显示，在界面波和底面回波之间 对于AL,O,和TN夹杂物，计算获得了从界面波到 可以明显地观察到由夹杂物引起的缺陷回波波形 30 30 (a -ALO h 20 界面波 一TiN 20 界面波 无夹杂物 夹杂物回波 底面回波 底面回波 0 0 0.85us -10 18s 1.7u8 20 0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 2.5 时间/以s 时间/μs 图5不同夹杂物的仿真波形.(a)A,O3和TN夹杂物的对比：(b)无夹杂物 Fig.5 Waveforms of different inclusions:(a)comparison of Al,0 and TiN:(b)no inclusion 图5(a)和(b)的对比结果显示：(1)两种夹杂物 的反射率和透射率分别为ru=0.558、rm=0.294和 回波均出现在0.85μs时刻，与理论计算结果0.87us 【u=0.829、lm=0.956.AL0,夹杂物的声波反射率是 基本一致；(2)含有夹杂物时，底面回波的到达时刻约 TN夹杂物的1.9倍，从而引起A山，0，夹杂物回波波幅 为1.7μ$，而图5(b)所示的不含夹杂物的底面回波约 大于TiN夹杂物. 为1.8μ$，即夹杂物的存在使相同条件下底面回波提 同时，考虑声压的透射率及夹杂物对声波的衰减 前0.1us,这是由于夹杂物的尺寸效应对材料中的声 作用，二者底面回波幅值的差异较小.由于仿真计算 速产生影响”：(3)夹杂物回波和底面回波相对于界 模型没有考虑晶粒尺寸和晶界等材料微结构，因而无 面波的展宽增加，这主要是由于夹杂物薄层的上下表 散射衰减，模型中的衰减主要来自扩散衰减和吸收衰 面对声波的多次反射和多次透射叠加，但薄层厚度仅 减叨.在声波频率∫相同的条件下，吸收衰减系数α。 为声波波长的1/3，因此不会引起明显的干涉现象a, =C∫中常数C,取决于材料的密度、杨氏模量等.由于 图5(a)中不同夹杂物的回波波形显示，l,0,夹 夹杂物的密度和杨氏模量不同，因此对声波的吸收衰 杂物的反射回波较TN夹杂物更为明显.这是由于夹 减也不同.仿真计算结果表明，山0，夹杂物对超声波 杂物的厚度仅为0.2mm,小于声波波长，其声压反射 的吸收衰减小于TN夹杂物. 和透射满足异质薄层的声压反射率和透射率公式： 进一步对不含夹杂物、含有AL,0,夹杂物和TN夹 4(m-）2 杂物模型的A扫波形进行快速傅里叶变换(FFT)得到 m 波形频谱图，如图6所示.其中，图6(a)和(b)分别对 (4) ）sin22m4 比了夹杂物回波和底面回波的频谱特性.图6(a)显 入， 示：A山O,夹杂物的频率幅值均大于TN夹杂物：两种 夹杂物均有各自的多个敏感频率，其中L,0,夹杂物 (5) 的敏感频率是8MHz和13MHz,而TiN夹杂物则是 8MHz和12.5MHz,由于夹杂物对超声波存在敏感频 式中，m是基体与夹杂物的声阻抗之比，山2为夹杂物的 率，从而使得该频率成分的超声波衰减较大.由 厚度，入，为夹杂物中的声波波长.根据表1所示材料 图6(b)可知：不含夹杂物的底面回波频谱分布呈高斯 参数可计算得到45钢、Al,0,夹杂物和TN夹杂物的 状，存在唯一的峰值：而含有夹杂物的底面回波频谱分 声阻抗分别为Zm=4.535×103kg·m2·s、Zu= 布变化较大，出现多个峰值.这是由于材料中部的夹 1.728×102kgm-2·s和Z=3.210×107kgm2. 杂物会相应地对其底面回波频率的衰减产生影响.同 s,将其代入式(4)和式(5)计算得到的夹杂物薄层 时图6(b)显示，与没有夹杂物时相比，夹杂物的存在
工程科学学报，第 37 卷，第 7 期 2 结果与讨论 采用上述有限元建模方法，建立具有不同类型夹 杂物和不同深度夹杂物的仿真模型，计算获得回波信 号，从时域和频域分析夹杂物对超声回波信号的影响 规律． 2. 1 不同类型的夹杂物 对于 Al2O3和 TiN 夹杂物，计算获得了从界面波到 底面回波之间的 A 扫波形，如图 5 所示． 其中，夹杂物 距离样品上表面的距离均为 d = 2. 5 mm． 图中横坐标 表示声波在材料中沿厚度方向的传播时间，纵坐标表 示载荷区域节点的 y 方向速度分量． 图 5( a) 为两种夹 杂物波形的对比结果，其中蓝色实线表示 Al2 O3 夹杂 物，红色虚线表示 TiN 夹杂物; 图 5( b) 为不含夹杂物 的 A 扫波形． 图 5( a) 显示，在界面波和底面回波之间 可以明显地观察到由夹杂物引起的缺陷回波波形． 图 5 不同夹杂物的仿真波形． ( a) Al2O3和 TiN 夹杂物的对比; ( b) 无夹杂物 Fig． 5 Waveforms of different inclusions: ( a) comparison of Al2O3 and TiN; ( b) no inclusion 图 5( a) 和( b) 的对比结果显示: ( 1) 两种夹杂物 回波均出现在 0. 85 μs 时刻，与理论计算结果 0. 87 μs 基本一致; ( 2) 含有夹杂物时，底面回波的到达时刻约 为 1. 7 μs，而图 5( b) 所示的不含夹杂物的底面回波约 为 1. 8 μs，即夹杂物的存在使相同条件下底面回波提 前 0. 1 μs，这是由于夹杂物的尺寸效应对材料中的声 速产生影响［7］; ( 3) 夹杂物回波和底面回波相对于界 面波的展宽增加，这主要是由于夹杂物薄层的上下表 面对声波的多次反射和多次透射叠加，但薄层厚度仅 为声波波长的 1 /3，因此不会引起明显的干涉现象［16］． 图 5( a) 中不同夹杂物的回波波形显示，Al2O3 夹 杂物的反射回波较 TiN 夹杂物更为明显． 这是由于夹 杂物的厚度仅为 0. 2 mm，小于声波波长，其声压反射 和透射满足异质薄层的声压反射率和透射率公式［16］: r = ( 1 4 m － 1 ) m 2 sin2 2πd2 λ2 1 + ( 1 4 m － 1 ) m 2 sin2 2πd2 槡 λ2 ， ( 4) t = 1 1 + ( 1 4 m － 1 ) m 2 sin ( 2 2πd2 λ 槡 ) 2 ． ( 5) 式中，m 是基体与夹杂物的声阻抗之比，d2为夹杂物的 厚度，λ2为夹杂物中的声波波长． 根据表 1 所示材料 参数可计算得到 45# 钢、Al2O3夹杂物和 TiN 夹杂物的 声阻抗分别为 Zsteel = 4. 535 × 107 kg·m － 2·s － 1、ZAl = 1. 728 × 107 kg·m － 2·s － 1 和 ZTiN = 3. 210 × 107 kg·m － 2· s － 1，将其代入式( 4) 和式( 5) 计算得到的夹杂物薄层 的反射率和透射率分别为 rAl = 0. 558、rTiN = 0. 294 和 tAl = 0. 829、tTiN = 0. 956． Al2O3夹杂物的声波反射率是 TiN 夹杂物的 1. 9 倍，从而引起 Al2O3夹杂物回波波幅 大于 TiN 夹杂物． 同时，考虑声压的透射率及夹杂物对声波的衰减 作用，二者底面回波幅值的差异较小． 由于仿真计算 模型没有考虑晶粒尺寸和晶界等材料微结构，因而无 散射衰减，模型中的衰减主要来自扩散衰减和吸收衰 减［17］． 在声波频率 f 相同的条件下，吸收衰减系数 αa = C1 f 中常数 C1取决于材料的密度、杨氏模量等． 由于 夹杂物的密度和杨氏模量不同，因此对声波的吸收衰 减也不同． 仿真计算结果表明，Al2O3夹杂物对超声波 的吸收衰减小于 TiN 夹杂物． 进一步对不含夹杂物、含有 Al2O3夹杂物和 TiN 夹 杂物模型的 A 扫波形进行快速傅里叶变换( FFT) 得到 波形频谱图，如图 6 所示． 其中，图 6( a) 和( b) 分别对 比了夹杂物回波和底面回波的频谱特性． 图 6( a) 显 示: Al2O3夹杂物的频率幅值均大于 TiN 夹杂物; 两种 夹杂物均有各自的多个敏感频率，其中 Al2 O3 夹杂物 的敏感频率是 8 MHz 和 13 MHz，而 TiN 夹杂物则是 8 MHz和 12. 5 MHz，由于夹杂物对超声波存在敏感频 率，从 而 使 得 该 频 率 成 分 的 超 声 波 衰 减 较 大． 由 图 6( b) 可知: 不含夹杂物的底面回波频谱分布呈高斯 状，存在唯一的峰值; 而含有夹杂物的底面回波频谱分 布变化较大，出现多个峰值． 这是由于材料中部的夹 杂物会相应地对其底面回波频率的衰减产生影响． 同 时图 6( b) 显示，与没有夹杂物时相比，夹杂物的存在 · 649 ·
陈丹等：金属材料内部非金属夹杂超声检测的数值模拟 ·947 1.2 一无夹杂物 0.6 一TiN 1.0 05 一TiN 0.8 0.4 0.6 0.3 0.4 0.2 0.1 0.2 0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 额率/MHz 顿率/MHz 图6夹杂物回波(a)和底面回波(b)的声波频谱 Fig.6 Spectra of inclusions'echo waves (a)and bottom echo waves (b) 使得底面回波各频率成分的幅值均有所降低 与底面回波相互混叠，观察不到完整的底面回波波形， 图5所示的时域波形结果表明：夹杂物与样品的 并且整个回波的幅值减小.近底面含有夹杂物和无夹 声阻抗差异越小，即夹杂物属性与本体材料属性越接 杂物的波形对比结果，如图7(b)所示.由于夹杂物的 近，其反射回波越弱，在有各类噪声和衰减存在的情况 存在，引起底面回波的提前和展宽增加，由此可以判断 下检测难度越大.图6所示的不同夹杂物的缺陷回波 夹杂物距底面的具体位置.近表面含有夹杂物和无夹 及底面回波的频谱对比结果表明：不同夹杂物对超声 杂物的波形对比结果，如图7(c)所示.与图7(b)对 波的敏感频率不同，AL,O,夹杂物对8MHz和13MHz 比可知，二者波形存在较大差异.界面波的展宽增 的超声波较为敏感，并由此引起较大衰减，而TN夹杂 加，而底面回波起始时刻的完全重合充分说明了夹 物对频率为8MHz和12.5MHz的超声波衰减较大. 杂物存在于样品上表面.上述结果显示，在实际检测 2.2不同深度的夹杂物 中，将样品翻转后进行测量可以得到近表面夹杂物 图5结果显示，当夹杂物位于材料中部，使得夹杂 的具体位置. 物缺陷回波信号与界面波和底面回波明显分离时，通 对AL,O,夹杂物不同夹杂位置的A扫波形作快速 过缺陷回波出现的时间可以得到夹杂物缺陷的深度信 傅里叶变换得到相应的频谱图，如图8所示.图8(a) 息.但是，对于材料内部的近表面或近底面夹杂物检 是夹杂物位于近表面时的界面波与底面回波的频谱对 测，夹杂物回波信号淹没在表面波或底面回波信号 比结果，其中界面波的频谱分布包含了完整的夹杂物 中而引起辨识困难.为了更好地观察夹杂物在材料 缺陷回波.从图可知，对于中心频率为10MHz的超声 内近表面不同位置引起的回波波形变化，选择声压 波，夹杂物引起高频部分的衰减要大于低频部分的衰 反射率较高的山0，夹杂物进行分析.其中，夹杂物 减，使得包含夹杂物回波的界面波中心频率向低频部 上表面分别距离材料上表面0.3mm和4.5mm(考 分偏移.其次，底面回波的频谱存在两个明显的峰值， 虑夹杂物的厚度为0.2mm,因而夹杂物下表面距材 峰值频率分别为8MHz和13MHz.需要注意的是，该 料底面的距离为0.3mm).计算获得材料中的A扫 峰值频率不同于图6()得到的敏感频率，前者是底面 波形，如图7所示. 回波的频谱成分，而后者是缺陷回波的频谱成分.图8 图7(a)显示，夹杂物距离样品上表面0.3mm时， (b)是夹杂物位于近底面时的界面波与底面回波的频 夹杂物回波与界面波发生混叠，引起界面波的展宽增 谱对比结果，其中底面回波包含了完整的夹杂物回波. 加.根据超声检测盲区的计算公式可知，当薄板中 界面波的频谱成分保持不变，而底面回波的频谱成分 的纵波声速为5748m·s,界面波的宽度为300ns时 发生劈裂，劈裂幅度相对较小，二者峰值处的频率分别 (界面波的宽度是由探头频率及加载时间共同决定 为10.5MHz和12.5MHz.上述劈裂产生的原因均是 的)，检测盲区为(300×10-9s×5748ms)12=0.86 超声波在敏感频率处产生较强的衰减而使单一的频谱 mm.图7(a)所示的夹杂物位于检测盲区内，其引起的 峰转变为多频谱峰.综上所述，夹杂物位置的不同，超 反射回波有部分淹没在界面波中，夹杂物回波在二者 声波敏感频率的衰减程度不同，相比较而言，夹杂物距 界面及夹杂物内部界面间发生多次反射和多次透射， 离材料上表面对敏感频率的衰减程度较大.通过超声 并且对底面回波也产生类似的影响，这种现象与多层 波的频谱分布可以定性地判断夹杂物的位置.关于峰 复合材料中声波的传播特性相近.当夹杂物位于材料 值劈裂的位置及大小是否与夹杂物的性质及尺寸有 近底面0.3mm时，界面波保持不变，夹杂物回波 关，还有待于进一步研究
陈 丹等: 金属材料内部非金属夹杂超声检测的数值模拟 图 6 夹杂物回波( a) 和底面回波( b) 的声波频谱 Fig． 6 Spectra of inclusions' echo waves ( a) and bottom echo waves ( b) 使得底面回波各频率成分的幅值均有所降低． 图 5 所示的时域波形结果表明: 夹杂物与样品的 声阻抗差异越小，即夹杂物属性与本体材料属性越接 近，其反射回波越弱，在有各类噪声和衰减存在的情况 下检测难度越大． 图 6 所示的不同夹杂物的缺陷回波 及底面回波的频谱对比结果表明: 不同夹杂物对超声 波的敏感频率不同，Al2 O3 夹杂物对 8 MHz 和 13 MHz 的超声波较为敏感，并由此引起较大衰减，而 TiN 夹杂 物对频率为 8 MHz 和 12. 5 MHz 的超声波衰减较大． 2. 2 不同深度的夹杂物 图 5 结果显示，当夹杂物位于材料中部，使得夹杂 物缺陷回波信号与界面波和底面回波明显分离时，通 过缺陷回波出现的时间可以得到夹杂物缺陷的深度信 息． 但是，对于材料内部的近表面或近底面夹杂物检 测，夹杂物回波信号淹没在表面波或底面回波信号 中而引起辨识困难． 为了更好地观察夹杂物在材料 内近表面不同位置引起的回波波形变化，选择声压 反射率较高的 Al2O3夹杂物进行分析． 其中，夹杂物 上表面分别距离材料上表面 0. 3 mm 和 4. 5 mm ( 考 虑夹杂物的厚度为 0. 2 mm，因而夹杂物下表面距材 料底面的距离为 0. 3 mm) ． 计算获得材料中的 A 扫 波形，如图 7 所示． 图 7( a) 显示，夹杂物距离样品上表面 0. 3 mm 时， 夹杂物回波与界面波发生混叠，引起界面波的展宽增 加［18］． 根据超声检测盲区的计算公式可知，当薄板中 的纵波声速为 5748 m·s － 1，界面波的宽度为300 ns时 ( 界面波的宽度是由探头频率及加载时间共同决定 的) ，检测盲区为( 300 × 10 － 9 s × 5748 m·s － 1 ) /2 = 0. 86 mm． 图 7( a) 所示的夹杂物位于检测盲区内，其引起的 反射回波有部分淹没在界面波中，夹杂物回波在二者 界面及夹杂物内部界面间发生多次反射和多次透射， 并且对底面回波也产生类似的影响，这种现象与多层 复合材料中声波的传播特性相近． 当夹杂物位于材料 近底面 0. 3 mm 时，界面波保持不变，夹杂物回波 与底面回波相互混叠，观察不到完整的底面回波波形， 并且整个回波的幅值减小． 近底面含有夹杂物和无夹 杂物的波形对比结果，如图 7( b) 所示． 由于夹杂物的 存在，引起底面回波的提前和展宽增加，由此可以判断 夹杂物距底面的具体位置． 近表面含有夹杂物和无夹 杂物的波形对比结果，如图 7( c) 所示． 与图 7( b) 对 比可知，二者波形存 在 较 大 差 异． 界面波的展宽增 加，而底面回波起始时刻的完全重合充分说明了夹 杂物存在于样品上表面． 上述结果显示，在实际检测 中，将样品翻转后进行测量可以得到近表面夹杂物 的具体位置． 对 Al2O3夹杂物不同夹杂位置的 A 扫波形作快速 傅里叶变换得到相应的频谱图，如图 8 所示． 图 8( a) 是夹杂物位于近表面时的界面波与底面回波的频谱对 比结果，其中界面波的频谱分布包含了完整的夹杂物 缺陷回波． 从图可知，对于中心频率为 10 MHz 的超声 波，夹杂物引起高频部分的衰减要大于低频部分的衰 减，使得包含夹杂物回波的界面波中心频率向低频部 分偏移． 其次，底面回波的频谱存在两个明显的峰值， 峰值频率分别为 8 MHz 和 13 MHz． 需要注意的是，该 峰值频率不同于图 6( a) 得到的敏感频率，前者是底面 回波的频谱成分，而后者是缺陷回波的频谱成分． 图 8 ( b) 是夹杂物位于近底面时的界面波与底面回波的频 谱对比结果，其中底面回波包含了完整的夹杂物回波． 界面波的频谱成分保持不变，而底面回波的频谱成分 发生劈裂，劈裂幅度相对较小，二者峰值处的频率分别 为 10. 5 MHz 和 12. 5 MHz． 上述劈裂产生的原因均是 超声波在敏感频率处产生较强的衰减而使单一的频谱 峰转变为多频谱峰． 综上所述，夹杂物位置的不同，超 声波敏感频率的衰减程度不同，相比较而言，夹杂物距 离材料上表面对敏感频率的衰减程度较大． 通过超声 波的频谱分布可以定性地判断夹杂物的位置． 关于峰 值劈裂的位置及大小是否与夹杂物的性质及尺寸有 关，还有待于进一步研究． · 749 ·
·948· 工程科学学报，第37卷，第7期 一距表面0.3mm 20 一距底面0.3m 10 0 AAMAA -20 306 0.5 1.0 1.5 2.0 2.5 时间s 30 30 无夹杂物 无夹杂物 20 一距底而0.3mm 20 手底而03m 10 10 -10 10 20 20 0% 05 1.0 300 1.5 2.0 2.5 1.0 15 2.0 2.5 时间/s 时间/μs 图7A山20不同夹杂深度的对比.(a)近表面与近底面的对比：(b)近底面与无夹杂对比：()近表面与无夹杂对比 Fig.7 Comparison of different depths of Al,0:(a)near-surface and near-bottom:(b)near-bottom and no-inclusion:(c)near-surface and no-in- clusion 2.0 25 1.8 一夹杂物回波与界面波交叠 界面波 夹杂物回波与张面回波交叠 1.6 底面回波 2.0 1.4 1.2 1.5 1.0 0.8 1.0 0.6 0.4 05 0.2 0.5 .0 15 2.D 2.5 0.5 1.0 15 2.0 2.5 频率/MHz 频率/MHz 图8不同夹杂深度的频谱图.(a)近表面0.3mm处：(b)近底面0.3mm处 Fig.8 Spectral distribution of inclusion:(a)0.3 mm near the surface:(b)0.3 mm near the bottom 图9(a)是无缺陷处的A扫波形，图9(b)是缺陷处的 3实验验证 A扫波形.实验结果显示，当样品近底面存在缺陷时， 由于含有较大尺寸夹杂物的样品难以获得，并且 缺陷回波与底面回波相耦合，其回波较无缺陷处的底 目前尚无其他方法能够确定材料内部夹杂物的具体位 面回波提前约200s,对应的缺陷深度约为590μm,该 置.因此，本文采用与仿真模型具有相同厚度H= 深度与夹杂物距底面距离相当.同时，实测结果与仿 5.0mm,但底面具有与夹杂物尺寸相同的微孔样品，验 真结果中图7(b)的底面回波波形存在差异，这是由于 证上述方法对材料内部微缺陷检测的有效性 微孔与近底面夹杂物属于不同的缺陷类型，底面微孔 实验样品为45钢，尺寸为10mm×10mm×5mm, 只包含固-一液界面，而近底面夹杂物存在上下两个固一 采用PVA SAM3O0超声显微镜设备，对样品进行检测， 固界面，因此实验中的底面回波不存在展宽拖尾现象. 其中聚焦探头的频率为10MHz.分别提取缺陷处和非 上述实验结果表明，采用超声方法可以实现材料内部 缺陷处的A扫波形进行对比，如图9所示.其中， 一定深度的大尺寸缺陷或夹杂物的检测与定位
工程科学学报，第 37 卷，第 7 期 图 7 Al2O3不同夹杂深度的对比． ( a) 近表面与近底面的对比; ( b) 近底面与无夹杂对比; ( c) 近表面与无夹杂对比 Fig． 7 Comparison of different depths of Al2O3 : ( a) near-surface and near-bottom; ( b) near-bottom and no-inclusion; ( c) near-surface and no-inclusion 图 8 不同夹杂深度的频谱图． ( a) 近表面 0. 3 mm 处; ( b) 近底面 0. 3 mm 处 Fig． 8 Spectral distribution of inclusion: ( a) 0. 3 mm near the surface; ( b) 0. 3 mm near the bottom 3 实验验证 由于含有较大尺寸夹杂物的样品难以获得，并且 目前尚无其他方法能够确定材料内部夹杂物的具体位 置． 因此，本文采用与仿真模型具有相同厚度 H = 5. 0 mm，但底面具有与夹杂物尺寸相同的微孔样品，验 证上述方法对材料内部微缺陷检测的有效性． 实验样品为 45# 钢，尺寸为 10 mm × 10 mm × 5 mm， 采用 PVA SAM300 超声显微镜设备，对样品进行检测， 其中聚焦探头的频率为 10 MHz． 分别提取缺陷处和非 缺 陷 处 的 A 扫 波 形 进 行 对 比，如 图 9 所 示． 其中， 图 9( a) 是无缺陷处的 A 扫波形，图 9( b) 是缺陷处的 A 扫波形． 实验结果显示，当样品近底面存在缺陷时， 缺陷回波与底面回波相耦合，其回波较无缺陷处的底 面回波提前约 200 ns，对应的缺陷深度约为 590 μm，该 深度与夹杂物距底面距离相当． 同时，实测结果与仿 真结果中图 7( b) 的底面回波波形存在差异，这是由于 微孔与近底面夹杂物属于不同的缺陷类型，底面微孔 只包含固--液界面，而近底面夹杂物存在上下两个固-- 固界面，因此实验中的底面回波不存在展宽拖尾现象． 上述实验结果表明，采用超声方法可以实现材料内部 一定深度的大尺寸缺陷或夹杂物的检测与定位． · 849 ·
陈丹等：金属材料内部非金属夹杂超声检测的数值模拟 ·949· 40 30 (a) 尤缺陷处 缺陷处 20 2 10 0 -10 -10 -20 -20 -30 -30 -40 -40 50 -50 -60 -60 50 1000 15002000 2500 500 1000 1500 20002500 时间/s 时s 图9实验波形对比.(a)无缺陷：(b)有缺陷 Fig.9 Experimental waveform comparison:(a)without flaw:(b)with flaw [5]Ogilvy J A.A model for the effects of inclusions on ultrasonie in- 4 结论 spection.Ultrasonics,1993,31(4):219 [6]Darmon M,Calmon P,Bele B.Modelling of the ultrasonic re- 通过建立包含夹杂物缺陷的二维钢板模型，采用 sponse of inclusions in steels /AlP Conference Proceedings.Bell- 有限元数值模拟的方法，对具有不同属性夹杂物和不 ingham,2003:101 同深度夹杂物的材料内部超声波场、超声回波时域波 7]Aggelis D G.Numerical simulation of surface wave propagation in 形、界面波、夹杂物缺陷回波和底面回波频谱分布进行 material with inhomogeneity:inclusion size effect.NDT E Int, 2009,42(6):558 了研究.对于尺寸为0.2mm×1.5mm的大型夹杂物 8]Dattal D,Kishoret NN.Features of ultrasonic wave propagation to 而言，采用中心频率为10MHz的超声波探头，可以获 identify defects in composite materials modeled by finite element 得有效的夹杂物回波信号，实现夹杂物的定位，从而为 method.NDT&E/m,1996,29(4):213 工业上采用中低频探头检测材料内部大尺寸夹杂物及 9]Liu D S,Chiou DY.Modeling of inclusions with interphases in 其成像提供可靠依据和指导.具体结论如下： heterogeneous material using the infinite element method.Comput Mater Sci,2004,31(3-4):405 (1)对不同属性夹杂物的超声检测，在时域波形 0] Autrique J C,Magoules F.Studies of an infinite element method 和频域波形上都会引起超声波形的显著变化.夹杂物 for acoustical radiation.Appl Math Modell,2006,30(7):641 与基体材料的声阻抗差异越小，通过时域和频域波形 01] Wang X D,Huang G L.Identification of embedded cracks using 进行检测的难度越大，声阻抗差异的大小是影响夹杂 back-propagating elastic waves.Inperse Pro,2004,20(5):1393 D2] 物检测的重要因素 Higdon R L.Absorbing boundary conditions for difference ap- proximations to the multi-dimensional wave equation.Math Com- (2)不同属性夹杂物对超声波的敏感频率不同. put,1986,47(176):437 与没有夹杂物时相比，夹杂物的存在使得底面回波各 03] Rose J L.Ultrasonic Wares in Solid Media.Cambridge:Cam- 频率成分的幅值均有所降低 bridge University Press,2004 (3)对不同深度夹杂物的检测，在界面波和底面 [14]Sun J G,Li S G,Li S B,et al.The transducer study of ultrason- ic in solid by finite element simulation.Sens World,2008,14 回波之间会有明显的夹杂物回波，在频谱分布上会引 (3):10 起底面回波的峰值劈裂，近表面夹杂物还会引起界面 (孙继华，李书光，李树榜，等.有限元模拟固体中超声传播 波频谱的峰值劈裂. 的换能器研究.传感器世界，2008,14(3)：10) 05] Ying C F.Zhang S Y,Shen J Z.Ultrasound Scattering in Sol- 参考文献 ids.Beijing:Defense Industry Press,1994 Zhang Y,Wu JF,Miu L D.Research progress in analysis of non- (应崇福，张守玉，沈建中超声在固体中的散射.北京：国 metallic inclusions in steel.Baosteel Technol,2008(2):35 防工业出版社，1994) (张毅，郭君飞，缪乐德.钢中非金属夹杂物的分析研究进 [16]Shi Y W.Ultrasonic Testing.Beijing:Machinery Industry Press, 展.宝钢技术，2008(2)：35) 2005 Jiang X S.Nonferrous Inclusion in Steel.Beijing:Metallurgy In- (史亦伟.超声检测.北京：机械工业出版社，2005) dustry Press,2011 [17]Liu D,Tumer JA.Ultrasonic backscatter in two-phase media (姜锡山.钢中非金属夹杂物.北京：治金工业出版社， and its dependency on the correlation function /Proceedings of 2011) Meetings on Acoustics.Hong Kong,2012:045004 B]Zhang L,Thomas B G.State of the art in evaluation and control of 18] QiF,Chen JJ,Wang R T.The opposite phase superposition steel cleanliness.IS//Int,2003,43(3):271 method of decreasing blind zone in ultrasonic detection.Tech 4]Atkinson H V,Shi G.Characterization of inclusions in clean Acoust,2011,30(4):325 steels:a review including the statistics of extremes methods.Prog (齐飞，陈品品，王润田.反相叠加减小超声检测盲区的方 Mater Sci,2003,48(5):457 法.声学技术，2011,30(4)：325)
陈 丹等: 金属材料内部非金属夹杂超声检测的数值模拟 图 9 实验波形对比． ( a) 无缺陷; ( b) 有缺陷 Fig． 9 Experimental waveform comparison: ( a) without flaw; ( b) with flaw 4 结论 通过建立包含夹杂物缺陷的二维钢板模型，采用 有限元数值模拟的方法，对具有不同属性夹杂物和不 同深度夹杂物的材料内部超声波场、超声回波时域波 形、界面波、夹杂物缺陷回波和底面回波频谱分布进行 了研究． 对于尺寸为 0. 2 mm × 1. 5 mm 的大型夹杂物 而言，采用中心频率为 10 MHz 的超声波探头，可以获 得有效的夹杂物回波信号，实现夹杂物的定位，从而为 工业上采用中低频探头检测材料内部大尺寸夹杂物及 其成像提供可靠依据和指导． 具体结论如下: ( 1) 对不同属性夹杂物的超声检测，在时域波形 和频域波形上都会引起超声波形的显著变化． 夹杂物 与基体材料的声阻抗差异越小，通过时域和频域波形 进行检测的难度越大，声阻抗差异的大小是影响夹杂 物检测的重要因素． ( 2) 不同属性夹杂物对超声波的敏感频率不同． 与没有夹杂物时相比，夹杂物的存在使得底面回波各 频率成分的幅值均有所降低． ( 3) 对不同深度夹杂物的检测，在界面波和底面 回波之间会有明显的夹杂物回波，在频谱分布上会引 起底面回波的峰值劈裂，近表面夹杂物还会引起界面 波频谱的峰值劈裂． 参 考 文 献 ［1］ Zhang Y，Wu J F，Miu L D． Ｒesearch progress in analysis of nonmetallic inclusions in steel． Baosteel Technol，2008( 2) : 35 ( 张毅，邬君飞，缪乐德． 钢中非金属夹杂物的分析研究进 展． 宝钢技术，2008( 2) : 35) ［2］ Jiang X S． Nonferrous Inclusion in Steel． Beijing: Metallurgy Industry Press，2011 ( 姜锡山． 钢 中 非 金 属 夹 杂 物． 北 京: 冶 金 工 业 出 版 社， 2011) ［3］ Zhang L，Thomas B G． State of the art in evaluation and control of steel cleanliness． ISIJ Int，2003，43( 3) : 271 ［4］ Atkinson H V，Shi G． Characterization of inclusions in clean steels: a review including the statistics of extremes methods． Prog Mater Sci，2003，48( 5) : 457 ［5］ Ogilvy J A． A model for the effects of inclusions on ultrasonic inspection． Ultrasonics，1993，31( 4) : 219 ［6］ Darmon M，Calmon P，Bele B． Modelling of the ultrasonic response of inclusions in steels / / AIP Conference Proceedings． Bellingham，2003: 101 ［7］ Aggelis D G． Numerical simulation of surface wave propagation in material with inhomogeneity: inclusion size effect． NDT ＆ E Int， 2009，42( 6) : 558 ［8］ Dattal D，Kishoret N N． Features of ultrasonic wave propagation to identify defects in composite materials modeled by finite element method． NDT ＆ E Int，1996，29( 4) : 213 ［9］ Liu D S，Chiou D Y． Modeling of inclusions with interphases in heterogeneous material using the infinite element method． Comput Mater Sci，2004，31( 3--4) : 405 ［10］ Autrique J C，Magoulès F． Studies of an infinite element method for acoustical radiation． Appl Math Modell，2006，30( 7) : 641 ［11］ Wang X D，Huang G L． Identification of embedded cracks using back-propagating elastic waves． Inverse Probl，2004，20( 5) : 1393 ［12］ Higdon Ｒ L． Absorbing boundary conditions for difference approximations to the multi-dimensional wave equation． Math Comput，1986，47( 176) : 437 ［13］ Ｒose J L． Ultrasonic Waves in Solid Media． Cambridge: Cambridge University Press，2004 ［14］ Sun J G，Li S G，Li S B，et al． The transducer study of ultrasonic in solid by finite element simulation． Sens World，2008，14 ( 3) : 10 ( 孙继华，李书光，李树榜，等． 有限元模拟固体中超声传播 的换能器研究． 传感器世界，2008，14( 3) : 10) ［15］ Ying C F，Zhang S Y，Shen J Z． Ultrasound Scattering in Solids． Beijing: Defense Industry Press，1994 ( 应崇福，张守玉，沈建中． 超声在固体中的散射． 北京: 国 防工业出版社，1994) ［16］ Shi Y W． Ultrasonic Testing． Beijing: Machinery Industry Press， 2005 ( 史亦伟． 超声检测． 北京: 机械工业出版社，2005) ［17］ Liu D，Turner J A． Ultrasonic backscatter in two-phase media and its dependency on the correlation function / / Proceedings of Meetings on Acoustics． Hong Kong，2012: 045004 ［18］ Qi F，Chen J J，Wang Ｒ T． The opposite phase superposition method of decreasing blind zone in ultrasonic detection． Tech Acoust，2011，30( 4) : 325 ( 齐飞，陈晶晶，王润田． 反相叠加减小超声检测盲区的方 法． 声学技术，2011，30( 4) : 325) · 949 ·
