一种基于同步迭代重建算法合金结构腐蚀损伤监测方法

申请号 CN201611183452.3 申请日 2016-12-20 公开(公告)号 CN107045003A 公开(公告)日 2017-08-15
申请人 北京航空航天大学; 发明人 张卫方; 高晓岩; 王翔宇; 王邓江; 何晶靖; 张宁; 赵炎; 谢宇宽;
摘要 一种基于同步 迭代 重建 算法 的 铝 合金 结构 腐蚀 损伤监测方法,步骤如下:一:选用材料 铝合金 板和压电 传感器 ,并采用正方形阵列布局;二:采用50kHz的中心 频率 的 正弦波 激励 信号 ,接收到的兰姆波是一种多种模式同时存在的频散波,即反对称波;三:基于SIRT的 断层 扫描算法原理是将板的待 检测区域 离散化,分为很多小网格,对损伤进行成像;四:采用 氢氟酸 溶液腐蚀铝合金板来制造腐蚀损伤;五:均值滤波处理;通过以上步骤,针对铝合金板结构,通过多个 压电传感器 的布局优化及信号监测,采集不同压电传感器并通过 信号处理 提取特征参数,能定量表征出铝合金构件的腐蚀损伤程度,解决了能对铝合金结构的腐蚀损伤的实时监测的功能。
权利要求

1.一种基于同步迭代重建算法合金结构腐蚀损伤监测方法,其特征在于:具体包括如下步骤来实现:
步骤一:
选用材料铝合金板和压电传感器,并采用正方形阵列布局;其中,在每边布置了四个压电传感器,合计16个压电传感器共组成了96条路径;采用了一发一收的压电传感器收发方式;因为其激励和接收压电传感器分布在测试区域的两侧,会有更少的边界反射信号的干扰,所以一发一收布局会对远距离损伤更为敏感;
步骤二:
采用了50kHz的中心频率正弦波激励信号,接收到的兰姆波即Lamb波是一种多种模式同时存在的频散波,即所分析的反对称波;
步骤三:
针对压电传感器网络布局和选取的信号特征量,运用基于同步重建迭代即SIRT的断层扫描算法;基于同步重建迭代即SIRT的断层扫描算法原理是将板的待检测区域离散化,分为很多小网格;压电传感器采用一发一收的方法,每个压电传感器都分别作为发射和接收点;压电传感器发射的Lamb波从各个网格中经过,信号穿过不连续部位衰减后被接收信号的压电传感器记录,根据比尔定律,接收信号压电传感器记录的信息等效为射线所穿过路径上物体线性衰减系数之和,然后根据各网格中的衰减系数值成像;当板材无缺陷时,各个网格中的衰减系数基本相同,当存在缺陷时,缺陷处网格中的衰减系数不同于无缺陷处的衰减系数,基于以上原理,对损伤进行成像;对于同步重建迭代算法即SIRT算法,需要根据腐蚀损伤监测的特点,对迭代算法进行优化,包括迭代次数的选择,基于半收敛判断的终止条件;
步骤四:
采用氢氟酸溶液腐蚀铝合金板来制造腐蚀损伤,当腐蚀损伤在Lamb波传播的路径上时,信号直达波的相位发生了改变;
通过分析Lamb波的信号相位变化进行损伤定位成像;Lamb波的幅值为5V,放大倍数为
10,ART算法迭代次数选为10;
步骤五:
为了使得成像结果更为平滑,试验采用了均值滤波处理,对于每个待处理的当前网格,都会选择一个模板,该模板为其邻近的预定个网格所组成,用模板的均值来替代原像素的值,其公式为
式中,M为该网格与其周围的网格总数,f(x,y)为在第x行第y列的像素值,s为所有的网格,g(x,y)为求得的最终像素值;
试验成像损伤定位结果与实际损伤位置相一致,试验结果证明本损伤定位方法能准确的定位出两个损伤;
通过以上步骤,针对铝合金板结构,通过多个压电传感器的布局优化及信号监测,采集不同压电传感器并通过信号处理提取特征参数,能定量表征出铝合金构件的腐蚀损伤程度,解决了能对铝合金结构的腐蚀损伤的实时监测的功能。
2.根据权利要求1所述的一种基于同步迭代重建算法的铝合金结构腐蚀损伤监测方法,其特征在于:在步骤二中所述的“激励信号”,它是指针对激励信号的选取,压电传感器的激励信号由信号发生器产生,信号由下面的表达式描述:
式中,fc为信号的中心频率,N为信号的波峰数,A为信号的幅度,H(t)为Heaviside阶梯函数。

说明书全文

一种基于同步迭代重建算法合金结构腐蚀损伤监测方法

技术领域

[0001] 本发明提供一种基于同步迭代重建算法的铝合金结构腐蚀损伤监测方法,它具体涉及一种基于同步迭代重建算法(SIRT)的铝合金结构的腐蚀损伤监测方法,属于结构健康监测技术领域。

背景技术

[0002] 结构健康监测技术实现了对传感器技术、信号处理和信号分析技术等的综合应用,可对结构进行实时监测,为大型关键结构的安全性、可靠性和耐久性提供了重要参考。兰姆波是可以在薄板结构上传播很长一段距离且振幅损失很少的弹性扰动。在兰姆波的监测中,监控结构的传感器数量可以显著减少。基于兰姆波的无损检测因其经济性有着广泛的应用,并且越来越多的被应用在结构健康监测中。迭代重建算法是在解决离散求逆问题中可计算稳定正则解的一种通用方法,广泛应用于医学成像,地球物理学,材料科学,和许多其他学科。同步迭代重建算法(Simultaneous Iterative Reconstruction Techniques,SIRT)是迭代正则法的一个分支,SIRT算法与代数迭代重建算法相比具有更好的稳定性,对测量误差及噪声不敏感,在得到相同质量的图像条件下,迭代重建算法需要更少的投影数据,故其对减少传感器的个数有很大的意义。本专利把其引用于含孔铝合金板的腐蚀损伤成像中。当铝合金板存在腐蚀损伤时,获取所需的波形,并运用适当地数据处理分析方法对其进行处理分析,从而提取损伤特征,并通过同步迭代重建法得到损伤的位置信息,实现了基于SIRT算法的铝合金结构的腐蚀损伤监测研究。

发明内容

[0003] (1)目的
[0004] 本专利针对航空铝合金腐蚀损伤难以实时有效监测的问题,提出了一种基于同步迭代重建算法的铝合金结构腐蚀损伤监测方法,即基于同步迭代重建算法的铝合金结构腐蚀损伤监测方法。本专利利用压电传感器阵列进行信号激励与采集,并利用先进的信号处理方法和信号分析方法提取损伤特征量。
[0005] (2)技术方案
[0006] 本发明针对铝合金结构腐蚀损伤的监测,主要通过同步迭代重建算法实现对铝合金结构腐蚀损伤的实时监测。
[0007] 本发明一种基于同步迭代重建算法的铝合金结构腐蚀损伤监测方法,具体包括如下步骤来实现:
[0008] 步骤一:
[0009] 本专利选用材料铝合金板和压电传感器,并采用正方形阵列布局;其中,在每边布置了四个压电传感器,合计16个压电传感器共组成了96条路径,如图2所示;本专利采用了一发一收的压电传感器收发方式;因为其激励和接收压电传感器分布在测试区域的两侧,会有更少的边界反射信号的干扰,所以一发一收布局会对远距离损伤更为敏感;
[0010] 步骤二:
[0011] 本发明采用了50kHz的中心频率正弦波激励信号,接收到的兰姆波即Lamb波是一种多种模式同时存在的频散波,在本发明中分析的为反对称波;
[0012] 步骤三:
[0013] 针对压电传感器网络布局和选取的信号特征量,运用基于同步重建迭代(SIRT)的断层扫描算法;基于同步重建迭代(SIRT)的断层扫描算法原理是将板的待检测区域离散化,分为很多小网格;压电传感器采用一发一收的方法,每个压电传感器都分别作为发射和接收点;压电传感器发射的Lamb波从各个网格中经过,信号穿过不连续部位衰减后被接收信号的压电传感器记录,根据比尔定律,接收信号压电传感器记录的信息可以等效为射线所穿过路径上物体线性衰减系数之和,然后根据各网格中的衰减系数值成像;当板材无缺陷时,各个网格中的衰减系数基本相同,当存在缺陷时,缺陷处网格中的衰减系数不同于无缺陷处的衰减系数,基于以上原理,对损伤进行成像;对于同步重建迭代算法即SIRT算法,需要根据腐蚀损伤监测的特点,对迭代算法进行优化,包括迭代次数的选择,基于半收敛判断的终止条件;
[0014] 步骤四:
[0015] 本发明采用氢氟酸溶液腐蚀铝合金板来制造腐蚀损伤,当腐蚀损伤在Lamb波传播的路径上时,信号直达波的相位发生了改变;
[0016] 本发明通过分析Lamb波的信号相位变化进行损伤定位成像;Lamb波的幅值为5V,放大倍数为10,ART算法迭代次数选为10;
[0017] 步骤五:
[0018] 为了使得成像结果更为平滑,试验采用了均值滤波处理,对于每个待处理的当前网格,都会选择一个模板,该模板为其邻近的若干个网格所组成,用模板的均值来替代原像素的值,其公式为
[0019]
[0020] 式中,M为该网格与其周围的网格总数,f(x,y)为在第x行第y列的像素值,s为所有的网格,g(x,y)为求得的最终像素值。
[0021] 试验成像损伤定位结果与实际损伤位置相一致,试验结果证明本损伤定位方法能准确的定位出两个损伤。
[0022] 其中,步骤一中所述的“压电传感器”是美国Steminc公司的压电传感器,它是SMD07T05R412WL压电陶瓷圆盘,其具体参数如表1所示。
[0023] 表1
[0024]
[0025] 其中,步骤二中所述的“激励信号”,它是指针对激励信号的选取,压电传感器的激励信号由任意信号发生器产生,信号可由下面的表达式描述:
[0026]
[0027] 式中,fc为信号的中心频率,N为信号的波峰数,A为信号的幅度,H(t)为Heaviside阶梯函数。
[0028] 其中,步骤三中所述的“SIRT算法”是指同步重建迭代算法,当信号穿过不连续部位衰减后被接收信号的压电传感器记录,根据比尔定律,接收信号压电传感器记录的信息可以等效为射线所穿过路径上物体线性衰减系数之和,然后根据各网格中的衰减系数值成像;
[0029] 先将图像离散化,设图像大小为N×N,每个像素的值为xi(i=1,2,...,N×N),压电传感器单元记录的投影数据为M个,设bj(j=1,2,...,M);由此可得,[0030]
[0031] 式中ai,j为第i个像素对第j个路径的贡献因子。
[0032] 整个成像方法即为解一个如下公式所示的静不定方程,
[0033] AX=b
[0034] 式中,A表示需要检测区域的真实形态,x为检测成像结果,b为A的投影;
[0035] 设bi是第i条射线经过被测区域的衰减值,则有公式:
[0036]
[0037] 式中,X(s)为衰减系数;设X(s)在第j个网格上为常数,Lij为第i条射线在第j个网格中的实际长度;
[0038] Lamb波的传播路径实际就是每个网格中各条射线实际走的长度Lij;该Lij求解的具体思路:
[0039] (1)以网格左上为原点,建立坐标轴,各网格点坐标和压电传感器坐标已知,压电传感器数为m;
[0040] (2)根据发射和接收压电传感器的坐标建立该射线的直线方程,求出直线与各个网格的交点;
[0041] (3)求解后一点与前一个点的距离,即可得Lij;
[0042] 则有公式:
[0043] bi=∑jLijxj
[0044] 式中,Lij组成的矩阵成为A矩阵;最终求解的方程:
[0045]
[0046] 式中m为射线的条数,n为网格的数量。
[0047] 可通过各种迭代算法求解该线性方程组,X矩阵即待重建图像的像素值。
[0048] SIRT的迭代公式如下:
[0049]
[0050] 式中h代表了第j条投影所穿过的像素。对于每条信号该算法都会求其与实际投影值的差,再加上权重并对所有投影线进行修正和求和,最终反投影得到对像素 的修正量。
[0051] 采用SIRT同步重建迭代算法求解以上方程,其基本思想是对图像每个像素的更新量都为对投影线的修正,修正通过贡献因子取加权平均,然后反投影得到。因为SIRT算法对所有的投影线的贡献进行了综合,所以可以避免其中一条投影线上的误差对整个图像造成较大的影响,稳定性更好。
[0052] 通过以上步骤,针对铝合金板结构,通过多个压电传感器的布局优化及信号监测,采集不同压电传感器并通过信号处理提取特征参数,能定量表征出铝合金构件的腐蚀损伤程度,解决了能对铝合金结构的腐蚀损伤的实时监测的功能。
[0053] (3)本发明的优点
[0054] 本发明提供了一种同步迭代重建算法的铝合金结构腐蚀损伤的监测方法。本方法通过选择实验材料,选择压电传感器网络布局,选取激励信号,最终开展试验后通过同步迭代重建算法来进行成像,以此来对腐蚀损伤进行监测,对腐蚀区域进行定位和定量,定量结果较为准确,腐蚀损伤面积相对误差不超过20%。并通过与实验原始图像对比证实了此方法的有效性。本专利具有准确性高,分析速度快,算法稳定等特点。本专利将同步迭代重建的算法用到航空铝合金结构腐蚀损伤监测为本专利的创新点。附图说明
[0055] 图1(a):压电传感器左视图示意图。
[0056] 图1(b):压电传感器主视图示意图。
[0057] 图2:压电传感器信号路径图。
[0058] 图3:监测方法图。
[0059] 图4:SIRT算法示意图。
[0060] 图5:基于氢氟酸溶液的腐蚀健康信号与损伤信号。
[0061] 图6:压电传感器布局图。
[0062] 图7:本发明所述方法流程图
[0063] 图中序号、符号、代号说明如下:
[0064] 图3中,Lij是指每个网格中各条射线实际走的长度;
[0065] 图4中,xk是指该像素的第k次迭代;L是指平面。

具体实施方式

[0066] 本发明一种基于同步迭代重建算法的铝合金结构腐蚀损伤监测方法,见图7所示,其具体实施步骤如下:
[0067] 步骤一:
[0068] 本专利选用材料Al-2024-T3,其密度杨氏模量和泊松比分别为2.78g/cm3,73.1GPa,0.33。其尺寸设计为500mm*500mm*2mm。
[0069] 本专利选用美国Steminc公司的压电传感器。压电传感器的尺寸和外观如图1、2所示。
[0070] 本专利根据待测的铝合金材料和压电传感器的几何特征,从线性阵列、矩形阵列、圆形阵列和平行阵列中选择矩形阵列。图4为压电传感器的正方形阵列布局。其中,在每边布置了四个压电传感器,16个压电传感器共组成了96条路径。
[0071] 步骤二:
[0072] 对于压电传感器的收发方式,目前可被分为两类:脉冲-回声(pulse-echo)布局和一发一收(pitch-catch)布局。对于脉冲-回声布局,激励压电传感器和接收压电传感器必须要在测试区域的同侧。因为边界反射信号的叠加的干扰,所以这种布局不适用于远距离检测损伤。而对于一发一收布局,因为其激励和接收压电传感器分布在测试区域的两侧,会有更少的边界反射信号的干扰,所以一发一收布局会对远距离损伤更为敏感。本专利采用一发一收的方式。
[0073] 接收到的Lamb波是一种多种模式同时存在的频散波。本专利分析的为第一个反对称波包。
[0074] 步骤三:
[0075] 针对压电传感器网络布局和选取的信号特征量,运用基于同步重建迭代(SIRT)的断层扫描算法。
[0076] 基于同步重建迭代算法(SIRT)的断层扫描算法原理是将板的待检测区域均匀的划分网络,如图3中的网格。压电传感器采用一发一收,每个压电传感器都分别作为发射和接收点。压电传感器发射的Lamb波从各个网格中经过,根据各网格中的衰减系数值成像。当板材无缺陷时,各个网格中的衰减系数基本相同,当存在缺陷时,缺陷处网格中的衰减系数不同于无缺陷处的衰减系数,基于以上原理,对损伤进行成像。
[0077] 对于SIRT算法,需要根据腐蚀损伤监测的特点,对迭代算法进行优化,包括迭代次数的选择,基于半收敛判断的终止条件。
[0078] 步骤四:
[0079] 本专利采用氢氟酸溶液腐蚀铝合金板来制造腐蚀损伤。图5给出了两组典型的氢氟酸腐蚀损伤对Lamb波传播信号影响的测试结果。从图5中健康和损伤信号的对比可以看出,当腐蚀损伤在Lamb波传播的路径上时,信号直达波的相位发生了改变。
[0080] 试验的压电传感器布局如图6所示,为正方形布局,共16个压电传感器。试样和压电传感器与本研究拟采用的材料一致。本专利通过分析Lamb波的信号相位变化进行损伤定位成像。Lamb为五周期的正弦调制波,中心频率分别为100kHz,幅值为5V,放大倍数为10,SIRT算法迭代次数选为10。
[0081] 步骤五:
[0082] 为了使得成像结果更为平滑,试验采用了均值滤波处理,对于每个待处理的当前网格,都会选择一个模板,该模板为其邻近的若干个网格所组成,用模板的均值来替代原像素的值。其公式为
[0083]
[0084] 式中,M为该网格与其周围的网格总数,f(x,y)为在第x行第y列的像素值,s为所有的网格,g(x,y)为求得的最终像素值。
[0085] 试验时把损伤置于所测区域的右下方。试验成像损伤定位结果与实际损伤位置相一致。
[0086] 对于两个损伤,试验结果证明本损伤定位方法能准确的定位出两个损伤。
[0087] 通过实验室分析结果,验证了所提出的基于断层扫描的铝合金监测方法的有效性。
[0088] 基于断层扫描的铝合金结构腐蚀损伤监测的应用验证是针对铝合金板结构,通过多个压电传感器的布局优化及信号监测,采集不同压电传感器并通过信号处理提取特征参数,定量表征出铝合金构件的腐蚀损伤程度,结合本专利腐蚀损伤的分析结果,验证所提出方法的有效性和适用性。
[0089] 其中,步骤一中所述的“美国Steminc公司的压电传感器”,它是指SMD07T05R412WL压电陶瓷圆盘,其具体参数如表1所示。
[0090] 表2
[0091]
[0092] 其中,步骤二中所述的“激励信号”,它是指针对激励信号的选取,压电传感器的激励信号由任意信号发生器产生,信号可由下面的表达式描述:
[0093]
[0094] 式中,fc为信号的中心频率,N为信号的波峰数,A为信号的幅度,H(t)为Heaviside阶梯函数。
[0095] 其中,步骤三中所述的“SIRT算法”是指同步重建迭代算法,整个算法过程如下:
[0096] 先将图像离散化,设图像大小为N×N,每个像素的值为xi(i=1,2,...,N×N),压电传感器单元记录的投影数据为M个,设bj(j=1,2,...,M)。由此可得,[0097]
[0098] 式中ai,j为第i个像素对第j个路径的贡献因子。
[0099] 由上式可得整个成像方法即为解一个如下公式所示的静不定方程,[0100] Ax=b
[0101] 式中,A表示需要检测区域的真实形态,x为检测成像结果,b为A的投影。
[0102] 设bi是第i条射线经过被测区域的衰减值,则有公式:
[0103]
[0104] 式中,X(s)为衰减系数。设X(s)在第j个网格上为常数,Lij为第i条射线在第j个网格中的实际长度。
[0105] Lamb波的传播路径实际就是每个网格中各条射线实际走的长度Lij。求解的具体思路:
[0106] (1)以网格左上角为原点,建立坐标轴,各网格点坐标和压电传感器坐标已知,压电传感器数为m。
[0107] (2)根据发射和接收压电传感器的坐标建立该射线的直线方程,求出直线与各个网格的交点。
[0108] (3)求解后一点与前一个点的距离,即可得Lij。
[0109] 则有公式:
[0110] bi=∑jLijxj
[0111] 式中,Lij组成的矩阵成为A矩阵。最终求解的方程:
[0112]
[0113] 式中m为射线的条数,n为网格的数量。
[0114] 可通过各种迭代算法求解该线性方程组,X矩阵即待重建图像的像素值。
[0115] SIRT的迭代公式如下:
[0116]
[0117] 式中h代表了第j条投影所穿过的像素。对于每条信号该算法都会求其与实际投影值的差,再加上权重并对所有投影线进行修正和求和,最终反投影得到对像素 的修正量。
[0118] 采用SIRT同步重建迭代算法求解以上方程,其基本思想是对图像每个像素的更新量都为对投影线的修正,修正通过贡献因子取加权平均,然后反投影得到。因为SIRT算法对所有的投影线的贡献进行了综合,所以可以避免其中一条投影线上的误差对整个图像造成较大的影响,稳定性更好。
[0119] 通过以上步骤,基于同步迭代算法的铝合金结构腐蚀损伤监测的应用验证是针对铝合金板结构,通过多个压电传感器的布局优化及信号监测,采集不同压电传感器并通过信号处理提取特征参数,定量表征出铝合金构件的腐蚀损伤程度,结合本专利腐蚀损伤的分析结果,验证所提出方法的有效性和适用性。
QQ群二维码
意见反馈