钛合金扩散焊水浸超声相控阵相位迁移成像方法*
2022-01-17杨宸旭李秋锋沈佳卉王海涛胡强陈
杨宸旭李秋锋沈佳卉王海涛胡 强陈 尧*
(1.无损检测技术教育部重点实验室(南昌航空大学),江西 南昌 330063;2.中国电建集团江西省电力设计院有限公司,江西南昌 330096;3.中建一局集团第二建筑有限公司,北京 100161)
近年来,对于扩散焊的研究相关检测技术引起了许多国内外学者的关注。扩散焊优点众多,不需要钎料以及保护气体;焊后,不需要机械加工;扩散焊由于其原理为高温高压下的原子扩散,不会改变母材特性[1-3]。扩散连接技术是钛合金构件连接的一种主要形式,该类焊接件在航空航天、工业制造等领域得到广泛应用[4-6]。
由于钛合金扩散焊常见缺陷比较微小,需要高精度和高灵敏度的检测方法,一般检测方法无法满足需求[7]。超声水浸检测是无损检测领域中广泛使用的一种检测方法,此检测方法操作简单方便,检测精度较高、易实现自动化检测。传统单探头超声水浸检测成像需要根据步进值移动探头位置,增加了检测时间和实验误差;且相同频率下单探头对于微小缺陷的检出率较低。
随着检测需求的日益增加,对超声水浸检测系统的检测效率及成像分辨率提出了更高的要求[8-10]。超声相控阵可以独立控制探头内的每个阵元使其声束可以进行偏转和聚焦,可以在不移动探头的情况下进行较大范围的成像。水浸超声相控阵检测结合了水浸检测和阵列检测的优点,相比传统单探头水浸超声检测可成倍的提升检测效率[11-12]。
水浸超声相控阵检测方法采用高频探头可获得较高的分辨率,同时相控阵探头可以灵活调整聚焦深度,利用这一点,检测过程中可以实现平行聚焦扫描同一深度缺陷,适用于钛合金板扩散焊检测。本文首次提出采用水浸超声相控阵检测钛合金扩散焊常见缺陷,相比传统超声相控阵楔块检测,该方法以水作为耦合剂实现超声检测;利用有限元仿真软件模拟相控阵平行聚焦扫描检测过程并提取各个阵元数据,将一个阵元组的各个阵元数据进行去延时叠加处理后进行双介质频域合成孔径成像。
1 基本理论
1.1 超声相控阵平行聚焦扫查原理
超声相控阵检测时,为提高实际检测缺陷时的检出率和准确性,通常采用声束聚焦检测微小缺陷[13]。图1中,第n个阵元距离中心距离为x n,聚焦深度为F,声波在介质中的传播速度为C,则根据阵元在声场区域的几何特性,第n个阵元的延时为:
图1 声束聚焦在P点
即:
式中:Δs为两个相邻阵元之间的声程差;d为阵元间距;t0为避免t n为负值的时间常数。
声束聚焦平行扫描为探头阵元组的各阵元按照设定好的聚焦法则,形成聚焦声束对缺陷进行实时扫查。在进行聚焦平行扫描前可选择一定数量的阵元作为阵元组中的子阵元,形成的多个阵元组分别激励并接收信号实现相控阵聚焦平行扫描。以4阵元为一阵元组为例,相控阵声束聚焦平行扫描示意图如图2所示。阵元组中的各个阵元延迟时间根据式(2)进行计算,利用MATLAB编程可以实现对各阵元延时的快速计算。
图2 相控阵声束平行聚焦扫描
1.2 双介质频域合成孔径聚焦成像
由于超声相控阵多采用高频探头检测缺陷,其声场声束波长短,声束扩散角较小,普通时域算法成像无法对偏离扩散角位置较远的点进行叠加成像。频域相控阵根据其爆炸反射原理,对时域信号的傅里叶变换并通过相位迁移(Phase Shift Migration,PSM)算法对检测区域进行成像,可以用于超声相控阵聚焦缺陷信号的成像。
频域合成孔径成像算法利用傅里叶变换所得的精确解进行成像,最早应用于地震波成像[14-15]。爆炸模型不仅适用于单介质声场,同样也适用于多介质声场中,频域合成孔径算法主要有两种成像方法,ω-k算法和PSM算法[16-17],其中PSM算法是直接根据波动方程推导,以得到声场精确解,是本文采用的成像算法。
根据波动方程,平面上任意一点声压可以表示为[18-19]:
式中:c为声波在介质中的传播速度;∇2为拉普拉斯算子。从获取数据P(x,z=0,t)傅里叶频谱中可以得到P(x,z,t=0)的成像数据:
式中:k x为x轴方向波矢。
根据式(3)和式(4)结合傅里叶变换唯一性和边界条件P(k x,z=0,ω),将缺陷处当做声源可得到:
假定在z>Z范围内,根据角谱法[20],若将所得到的频域波场Z外推至z,则有:
式中:k z为z轴方向波矢。
结合上式和逆傅里叶变换得到任意点(x,z)的二维图像:
本次实验中水层界面为规则平面,界面两侧均为均匀介质,即声速在某一深度为定值,根据Snell定理,将C表示为Z轴方向的函数,利用声压相位行进方向的连续性和双层介质界面的声透射系数,通过改变每一层z轴方向的波矢即可完成双介质PSM成像。本文将各阵元组聚焦声束的回波数据进行去延时叠加形成(探头阵元数-子阵元组所含阵元数+1)条A扫数据,将这些数据进行双介质PSM成像并与传统相控阵B扫成像对比。
2 水层厚度分析和模型建立
2.1 水层厚度分析
水浸检测的水层厚度是检测过程中重要参数的之一,根据多元高斯叠加定理可以获得声束在双层介质中的声场分布,高斯叠加系数[21]如表1所示。
表1 高斯叠加系数
利用瑞利积分形式表示声场中任意一点的声压;通过叠加高斯声束即可获得点声压幅值;计算两层介质中的质点振动速度并根据声束透射系数即可得到多元高斯叠加声场模型,模拟水-钛合金声场分布情况。本次实验设定探头频率为5 MHz,一阵元组孔径大小8.5 mm,水和钛合金传播速度分别为1 482 m/s、6 332 m/s,水层厚度5 mm、10 mm时水-钛合金声场分布如图3所示。
图3 水-钛合金声场分布
根据图3(a),水层厚度为5 mm时,钛合金声场部分声束带有旁瓣,旁瓣会分散主瓣能量甚至产生伪像,从而降低成像分辨率。将水层厚度增至10 mm时,如图3(b)所示,旁瓣问题得到明显改善,若继续增加水层厚度,高频率声波会随着声程的增加而使其能量受损严重,这将导致缺陷回波幅值过低从而无法成像的问题。
2.2 模型建立
根据实际待检测对象在COMSOL中建立有限元模型,模型如图4所示。水层厚度为10 mm,宽为50 mm,钛合金厚10 mm,宽度50 mm。扩散连接为原子间的扩散,所以除缺陷区域外,两块钛合金薄板可看做一个整体[22]。
图4 测试模型结构图
在扩散连接过程中,受焊接工艺与其他因素影响,常产生中部高两端方向逐渐变窄的大面积未复合缺陷;类圆球状微小气孔缺陷以及方向、长短不同的弥散分布型条状缺陷[23]。根据钛合金扩散焊常见缺陷及其基本形貌特点,设置宽6 mm、厚0.2 mm椭圆形缺陷为未焊合缺陷;直径0.2 mm点状缺陷为微孔缺陷;宽1.5 mm、厚0.1 mm矩形缺陷为紧贴型缺陷,三种缺陷中心距上表面15 mm。
设置阵元中心距为0.6 mm,阵元宽度为0.5 mm,阵元数目为64,入射角度设定为0°,又以水层作为耦合剂,故忽略横波干扰。在64个阵元中,以16阵元为一阵元组,超声相控阵平行聚焦扫描的每组阵元延时法则均相同,组中16个阵元延时时间如图5所示,通过这种激励方式可以将每组阵元发射的声束聚焦在水层上表面15 mm深的位置。
图5 阵元组中1-16阵元延时激励方式
通过计算缺陷位置在远场区内,不会受到声压急剧起伏的影响,可以良好的进行检测。在本模型中激励信号选择高斯脉冲,其具有能量集中、余波能量小、旁瓣对主波的干扰小、指向性好等优点,数学表达式为:
在有限元中,网格划分得越小,计算的准确性越高,但需要较高计算机配置和较长计算时间。为保证计算精度,有限元的最大网格尺寸应介于波长的0.05倍到0.1倍之间。考虑到本模型大小和实际条件,选择每个波长包含10个单元间隔即可满足仿真需求。由于波束在水和钛合金中的波长不同,所以在两种介质的网格划分需要分开计算,其各自最大网格单元大小为:
相似的,时间步长的选取对计算结果的准确性同样影响较大,为保证声压模拟的稳定性,采样时间应不大于声波在单元网格中的传播时间,即:
式中:Cmax为超声波在两种介质中的最大声速。另外,模型中的材料根据检测方案选择COMSOL仿真软件内置材料:H2O(water)[liquid]和Ti-6Al-4V[solid,polished]。图6为第一阵元组声束聚焦在缺陷附近,时间为7.735μs。
图6 相控阵声束聚焦
3 钛合金扩散焊缺陷成像研究
3.1 回波信号分析
通过有限元仿真软件获得所有阵元组中的各个阵元所接收的信号,将组中全部阵元进行去延时叠加得到一条合成信号,最终可得到49条A扫数据。以第一阵元组为例,图7为组中一部分阵元信号与组中全部阵元的合成信号。
根据图7(a)与图7(c)可以看出,一组阵元中的单阵元能量有限或基于位置限制,回波信号较差,若直接使用单阵元数据成像必定会造成图像缺陷不清晰。图7(b)所表示的5号阵元回波信号相比1号阵元和12号阵元的回波信号更“干净”,但缺陷的回波信号不明显,且无缺陷区域信号不平稳,幅值逼近缺陷信号幅值,若直接成像会出现伪像。将阵元组中1-16号阵元分别去掉对应的延时并叠加得到图7(d),从时域信号可以看出,非缺陷区域信号平稳且幅值明显低于缺陷信号幅值;从频域信号可以看出,从左至右的高亮云图分别对应了声场中的始波、水-钛合金界面波、缺陷波、底波,信号效果得到显著提升。
图7 A扫信号时频图
3.2 成像对比及分析
将49条合成信号分别进行超声相控阵B扫成像以及相位迁移算法成像,由于本次实验以16阵元为一阵元组,聚焦声束主声轴为组中中心位置阵元间的轴线,所以成像区域的横向范围为30 mm,纵向成像区域为25 mm,如图8所示。对图像进行-55 dB处理,由8(a)可以看出,超声相控阵传统B扫成像可以较清晰的呈现未焊合缺陷和紧贴型缺陷,但未焊合缺陷成像出现伪像情况,会影响缺陷位置判断。对于微孔缺陷,其缺陷特征无法呈现,在钛合金扩散焊实际应用中,微孔缺陷会影响工件在役服务寿命或造成更加严重的后果。
由图8(b),相比超声相控阵传统B扫成像,PSM算法成像的最大特点在于可以检出直径为0.2 mm的微孔缺陷;其次由于PSM算法根据波动方程求解的精确解并逐层递推成像,使得图像分辨率显著提高,水-钛合金界面与三种缺陷成像更加平滑,可以轻易分辨三种缺陷的类型。
图8 成像结果对比
3.3 缺陷图像量化对比
为直观对比两种情况下的相控阵聚焦声束成像效果,三种缺陷的对比率(Contrast Ratio,CR)和定位误差如表2、表3所示。
表2 三种缺陷CR对比 单位:dB
表3 三种缺陷位置误差 单位:%
根据表2、表3,PSM算法下的相控阵聚焦声束成像所对应的未焊合、紧贴型缺陷CR分别提高了6.93 dB、5.84 dB,未焊合缺陷与紧贴型缺陷位置误差率降低了5.37%、4.66%。
4 结论
钛合金扩散焊过程中产生的未焊合、微孔、紧贴型缺陷,采用水浸超声相控阵平行聚焦扫查检测缺陷,制定49个以16阵元为一组的阵元组。根据49条合成A扫信号进行超声相控阵传统B扫成像,在-55 dB处理下未检测到0.2 mm直径微孔缺陷,且缺陷位置误差较大。
基于相控阵合成信号进行PSM算法成像,不仅检测出微孔缺陷,且图像分辨率缺陷位置精度都有显著提升,CR平均提升6.39 dB,未焊合、紧贴型两种缺陷位置误差率偏平均降低5.02%,本次实验为扩散焊水浸超声相控阵检测提供理论依据。