北京航空航天大学超声无损检测实验室2022年度工作进展.pdf

无损检测2023年第45卷第9 期85IDT院企风采北京航空航天大学超声无损检测实验室2022年度工作进展2022年，北京航空航天大学超声无损检测实验室在超声检测基础理论研究、超声仿真系统设计以及超声检测系统研发上取得了系列成果。在碳纤维增强树脂基复合材料纤维褶皱缺陷的超声检测方面，建立了利用复合材料层间反射信号对褶皱缺陷进行表征的理论模型，针对层间信号微弱、超声信号相位敏感等特性，提出了利用线性相位滤波器缓解噪声引起错位的方法。在碳纤维复合材料的阵列超声检测方面，提出了基于Viterbi搜索算法的声线示踪方法，建立了CFRP材料阵列超声检测的声束路径计算模型，改进了CFRP材料的阵列超声全聚焦检测成像算法。在异种粗晶组织材料扩散焊界面的超声检测方面，建立了粗晶母材、扩散焊界面、复杂形面声学检测模型，实现了粗晶组织材料扩散焊界面静态声场计算与缺陷响应回波信号仿真，为异种粗晶组织材料扩散焊检测提供了理论基础。在超声检测仿真系统设计方面，建立了阵列超声声线路径仿真、超声声场分布仿真以及超声检测动态声场仿真等多个核心仿真模块，结合超声组件三维建模功能，可实现超声无损检测的快速仿真成像分析。在超声检测系统研发方面，研制了“5+1”轴水浸超声C扫描检测系统和双通道空气耦合超声检测系统，可分别应用于复杂曲面试样和复合材料板材的自动化检测。1碳纤维增强树脂基复合材料纤维褶皱缺陷的超声检测研究纤维褶皱缺陷是碳纤维复合材料中常见的几何形状缺陷。在碳纤维复合材料零部件制造过程中，由于预浸料滑移引起铺层弯曲变形，这种内部铺层偏离预设方向的缺陷称为褶皱缺陷。相关研究表明，碳纤维复合材料中的褶皱缺陷通常是应力集中点,严重降低了材料的力学性能。与分层、脱黏、夹杂等引起复合材料内部材料不连续的缺陷不同，褶皱缺陷基本上不改变局部材料的连续性，直接利用幅值信息的常规超声无损检测方法难以有效评估。复合材料褶皱缺陷检测利用常规超声检测中视为结构噪声的层间反射信号反演褶皱缺陷的几何形态特征。褶皱缺陷是铺层连续变化的缺陷，褶皱的几何形态特征是影响力学性能的重要因素，褶皱检测需要利用超声信号反演出其几何形态特征。基于全矩阵数据的全聚焦方法（TFM-FMC）通过虚拟延时聚焦成像，可以达到很高的空间精度。由于层间反射信号微弱，全聚焦检测结果易受到各种噪声的干扰，需要设计离散时间滤波器来消除噪声干扰。由于全聚焦算法是对全矩阵数据中A型信号进行延时叠加处理，所以要求滤波器对A型信号无损检测2023年第45卷第9 期86IDT院企风采的时间信息没有影响或可以校正。由于褶皱缺陷的几何形态特征是空间信息，其对信号的相位信息要求较高，应尽量避免相位失真，需要选择线性相位滤波器或者频域滤波器，因此，提出了一种使用线性相位滤波的超声阵列数据进行皱褶检测的方法。为了尽量减小层间反射回波之间的干涉影响，选择利用中心频率为10 MHz的线阵换能器（10 L6 4-0.35-D 7 7 EJA 557）采集数据。利用短时傅里叶变换对全矩阵数据内的A型信号进行时频分析（见图1），确定层间回波信息主要集中在9MHz附近，据此设计线性相位滤波器，其幅频特性与相频特性如图2 所示。对含有随机孔洞的人工褶皱缺陷（见图3）进行检测，其滤波前、后全聚焦成像结果与全聚焦相位成像的余弦值如图 4 所示。对比图 4(c)与图 4(d)可知,滤波后成像结果中的错位现象明显得到改善。图5为滤波后全聚焦相位成像余弦值的峰值与褶皱横截面叠加对比图，可见，滤波后成像结果可以准确地反演出褶皱缺陷的几何形态特征。20上表面回波(1,1)5015层间川波(sr)/010-505-100底血回波01501618202224频率/(MHz)(a)分析图例一20(1.8)(H/8P)/5015(srl)/E0105010001501618202224频率/(MHz)(b)分析图例二20(1.16)5015(srl)/elF8010505-10001501618202224频率/(MHz)(c)分析图例三20(1.24)(ZHW/8P)/率影率5015(srl)/El0105010001501618202224频率/(MHz)(d)分析图例四图1全矩阵数据中部分A型信号时频分析图1.5O(e)/0.560102030405060频率/MHz图2滤波器的幅频特性与相频特性10000孔洞孔洞图3褶皱试样图0(ww)/z205-1071320 x/(mm)(a）滤波前全聚焦成像结果错位0(ww)/z2071071320 x/(mm)(b）滤波前全聚焦相位成像的余弦值无损检测2023年第45卷第9 期87余弦值IDT院企风采0(ww)/z20571071320 xl(mm)(c）滤波后全聚焦成像结果02(w)/z05-1071320 x/(mm)(d）滤波后全聚焦相位成像的余弦值图4滤波前后全聚焦相位成像结果与全聚焦相位的10000止炫值峰值图5滤波后全聚焦相位成像余弦值的峰值与褶皱横截面叠加图2石碳纤维增强树脂基复合材料缺陷的阵列超声全聚焦成像检测研究在先进树脂基复合材料（CFRP）的超声检测中，常规超声检测效率低，而阵列超声全聚焦成像检测技术则依赖准确的声传播延时。针对CFRP中的各向异性和多层折射界面而导致声波延时计算困难的问题，提出了一种使用Viterbi 搜索算法的声线示踪方法，用于计算阵列超声全聚焦成像检测的延时法则，从而实现对CFRP中典型缺陷的阵列超声全聚焦成像检测。首先，通过分析碳纤维铺叠角度和不同角度铺层的弹性系数矩阵，计算qP波在0、45和9 0 纤维铺层内不同传播角度下的传播速度，如图6 所示。0.50元0010000459080000.75元0.25600040002000F1.00元0.0020004000600080001.25元1.75.元100001.50元图6 在0、45和9 0 单向铺层中qP波群速度随传播角度的变换其次，通过费马原理和Viterbi搜索算法计算CFRP内部声传播路径和声传播时间，CFRP离散模型中，使用底面反射法（BRM）计算第1个阵元发射、不同位置阵元接收的声线路径示意图如图7 所示。第1个件元发射，第16 个牌元按收第1个降元发时，第3 2 个陈元接收第1个间元发射，第48 个阵元楼收一第1个体元发时，864个味元接收图7 第1个阵元发射、不同位置阵元接收的声线路径示意图然后，为了验证声线示踪方法的准确性，使用Viterbi声线示踪法、BRM试验以及固定声速法分别计算3 种方法的声传播时间差，如图8 所示。结果表明，使用固定声速法，只有前16 个阵元计算的声时接近BRM实验测量值，当接收阵元与发射阵元的距离更远时，两者差值迅速增加，说明CFRP试块中多层结构和层间各向异性对声波传播影响不可忽略，不能直接等效于在各向同性介质中传播；而声线示踪法计算的声时与BRM实验测量值接近，验证了声线示踪法的准确性。固定声速法(3055m/s)2.0声线示踪法BRM试验1.51.00.50.0010203040506070接收阵元图8 三种不同方法计算的qP波传播时间差曲线无损检测2023年第45卷第9 期88P+J院企风采最后，根据Viterbi声线示踪法计算的延时法则，使用阵列超声全聚焦算法实现了对CFRP构件中典型缺陷的成像检测，如图9 所示。结果表明，使用声线示踪法的阵列超声全聚焦图像，可以有效避免CFRP各向异性和层间反射的干扰,相干叠加的信号幅值增加，铺层显示连贯，可以有效检出复合材料中的缺陷。0-5-102-15(ww)/z-203褶镀缺陷-254-305-35-406-457502345678910111213141516171819201dBx/(mm)(a)成像一0-5-102-15(w)/z-203-25-305-35-406-455023567891011121314151617181920dBx/(mm)(b)成像二图9 Viterbi声线示踪法分别对CFRP构件上下表面的纤维褶皱缺陷进行全聚焦成像结果3异种粗晶材料扩散焊界面超声检测研究为满足航空航天极端环境的使用要求，异种粗晶组织材料常通过扩散焊连接在一起形成多功能复合材料。由于粗晶组织材料会引起强烈的超声散射衰减和波形畸变，使现有的用于细晶材料扩散焊界面检测的高频超声和非线性方法难以适用。为研究适用于粗晶组织材料扩散焊界面检测的新方法，建立了考虑晶粒特性的粗晶组织金属材料仿真模型，探究了超声波在各种类型粗晶组织材料中的传播规律，如图10 所示。扩散焊缺陷尺寸远小于超声波波长，难以使用常规幅值方法进行表征，建立了扩散焊微孔缺陷、考虑复杂形面特性的声学模型，为后续粗晶扩散焊微小缺陷的识别分析提供了理论基础，如图11 所示。0.90.80.70.60.50.40.30.20.1000.20.40.60.81(a)晶粒建模（b）陈列超声声场仿直时网-0.9 16 6 7 us表蛋：速准Y分童(m/s)mm￥10 1968472654-2210246810mm(c)动态声场仿真反10业口之吉图10 粗晶组织金属仞料声场仿具声压/(Pa)3x105位移t-2.075us/(mm)1410-8缺陷回波扩散焊微孔缺陷0(a）扩散焊微孔缺陷响应无损检测2023年第45卷第9 期89IDT院企风采换能器耦合层入射面（水层）钛合金(TC4)吸收层一扩散焊缺陷（弹性薄层简化）(b）复杂形面模型0.5始波一扩散焊缺陷回波0-0.5底面回波012345时间/（s）(c）提取的A波信号图11扩散焊界面声场仿真4超声检测仿真软件系统设计通过超声检测仿真，可明确试样的超声传播特性，指导并优化超声检测工艺设计，显著提高超声检测效果。在采用有限元等通用仿真工具进行超声检测仿真时，存在操作繁琐、技术门槛高等局限性。为此，开发了超声检测仿真专用软件，可显著降低超声检测工艺设计过程的经济成本和时间成本。首先，基于OpenCascade三维CAD建模库，开发了检测模型的三维建模模块，可实现多种类型换能器、多种复杂结构试样、楔块以及缺陷的建模、显示以及交互，允许用户自主定义组件的结构特征和声学特征，进而准确表征超声检测工况，检测模型如图12 所示。此外，基于射线追踪算法和组件三维模型，建立了阵列超声声线仿真功能模块。图12 阵列超声检测三维模型其次，基于Rayleigh-Sommerfeld方法,开发了基于半解析方法的超声声场分布计算模块，可快速计算超声换能器在试样中所激发声场的三维分布，并可据此分析超声检测参数与声场分布特征的映射规律，声场分布显示结果界面如图13所示。BUAA,UT_SIMU口X30934252259焦斑深度范围：16.0 0-57.0 0 mm焦斑长度：47.51m11焦斑宽度：4.6 3 mm29读联保存图13超声声场分布显示界面最后，基于时域有限差分仿真方法，建立了动态声场仿真模块。该仿真模块通过三维检测模型交互的方式实现有限差分网格及阵列超声聚焦检测方案的设置，并利用GPU并行计算，极大地提升了仿真计算效率。通过计算并记录仿真区域和换能器晶片处应力和速度的变化，实现声场的动态传播过程和超声回波信号的仿真计算，如图14所示。在此基础上，基于超声回波计算结果，仿真系统提供了阵列超声扇形扫查、线性扫查以及全聚焦成像等多种检测仿真成像功能，可供用户进行仿真成像分析，软件界面如图15所示。t=4.65us1=8.60st=10.69usT=11.39s图14阵列超声动态声场传播过程无损检测2023年第45卷第9 期902023年5院企风采图15超声检测仿真软件界面5超声检测系统研发针对曲面构件的无损检测需求，研制了“5+1”轴水浸超声C扫描检测系统，如图16 所示。检测系统采用五轴加回转台的机械结构，可实现回转体构件和复杂曲面构件的自动化检测。支持常规单探头和阵列超声两种检测模式，检测软件具有工件仿形、轨迹规划、运动控制、C扫成像、3D显示成像等功能。系统现已进入检测实验测试阶段。图16“5+1”轴水浸超声C扫描检测系统碳纤维增韧碳化硅陶瓷基（C/SiC）复合材料是一种极具发展潜力的耐高温结构材料。为实现C/SiC复合材料的质量无损评估且满足工业自动化检测的需求，研制了一套适用于C/SiC复合材料的双通道空气耦合超声自动检测系统，如图17所示。系统的检测范围为3.0 m2.0m，双通道模式下检测速度可达1.0 m/s，支持矩形、梯形、三角形工件的自动扫描。检测软件具有运动控制、A型信号显示、C扫描成像、缺陷统计分析、检测报告生成等功能。C/SiC复合材料空气耦合超声检测系统可实现厚度2 0 mm的C/SiC复合材料3 0 mm分层缺陷的有效检出（尺寸误差小于15%），现已用于航天系统某部门的实际生产检测中。图17复合材料双通道空气耦合超声检测系统6结语2022年北京航空航天大学超声无损检测实验室取得的成绩离不开国内同行的大力支持与帮助。未来，实验室期望继续保持与同行的深入交流与合作，攻克相关领域的无损检测难题，促进我国先进无损检测技术的发展！北京航空航天大学机械工程及自动化学院周正干，马腾飞，朱甜甜，杨功鹏，王俊，于敬涛供稿