CN104034801B - 基于合成时反的结构损伤迭代聚焦成像监测方法 - Google Patents

Abstract

本发明公开了一种基于合成时反的结构损伤迭代聚焦成像监测方法，该方法包括：对需采用N－S阵列同步探测的实物结构进行N×S次独立的单通道激励-单通道接收的主动应力波探测；将S阵列中各单元的N路接收信号线性叠加得到S路合成信号；将S路合成信号进行时间反演放大，对实物结构进行N×S次独立单通道激励-单通道接收的时间反演探测；将N阵列中各单元的S路接收信号线性叠加，得到N路合成时反信号；采用合成时反继续对实物结构进行迭代探测，最终得到N路合成迭代信号；采用基于时域聚焦的时间反演成像方法，对合成迭代信号进行时间反演放大处理，并各自在成像模型N阵列单元激励回传，实现损伤聚焦成像。

Description

基于合成时反的结构损伤迭代聚焦成像监测方法

技术领域

本发明涉及测试测量技术及结构健康监测领域，尤其涉及一种基于合成时反的结构损伤迭代聚焦成像监测方法。

背景技术

时间反演是指对测量所得的时域信号的一种逆序处理，它将所测信号按照到达接收点的顺序进行前后反转。根据互易原理，时间反演具有时空聚焦特性，具体体现为：空间中多点测得的来自损伤目标的接收信号经过时间反演处理后，所得的多路时间反演信号能够无需先验知识、自适应地穿过复杂介质在同一时间在空间损伤目标处叠加出最大能量，并形成类似于原始损伤散射信号的逆时域波形的聚焦信号。时间反演技术利用时间反演信号的聚焦特性有效降低了频散效应和非均匀介质的影响。目前人们已经利用时间反演聚焦信号的峰值、能量、时空聚焦等特征，应用时间反演技术处理各种损伤测方法检测到的损伤信号后对多种结构进行了损伤检测。传统时间反演镜成像方法正是利用了时间反演技术的时空聚焦特性，对空间中多点测得的接收信号进行时间反演处理，然后把时间反演信号同步回传到多接收点上重新激励，记录成像空间中每一时刻的应力分布，只要找出成像空间中所有时间段最大信号出现的位置就能找到损伤目标。然而，传统时间演镜成像方法的目标成像依据只有在接收阵列单元数量较多时且传播路径较短、传播损耗较小、多损伤干扰小时，损伤目标点才可能成为空间最大信号所在的点。比如当接收阵列单元个数为1时，那么在时间反演过程中，到达损伤目标处的时间反演聚焦信号虽然能够在时域上恢复原始波形，但是该聚焦信号由于受到传播路径的损耗，其信号大小肯定要比用于激励时间反演信号的接收单元处的信号小，成像结果就会把目标错误定位在接收单元处。

根据以上分析，实际检测中若采用传统时间反演镜成像方法，需要针对多个发射阵列单元和接收阵列单元设计同步激励和同步采集装置，这对硬件设备的性能要求极高；此外，传统时间反演镜成像算法的准确度也受限于接收阵列单元个数、传播路径损耗、多损伤目标等因素限制。因此，建立一种满足测试测量需求并能有效实现最大目标检测的结构损伤成像监测方法具有重要的研究意义和前瞻性。

发明内容

为解决上述技术问题，本发明的目的是提供一种基于合成时反的结构损伤迭代聚焦成像监测方法。所述技术方案如下：

一种基于合成时反的结构损伤迭代聚焦成像监测方法，包括：

对需采用N－S阵列同步探测的实物结构进行N×S次独立的单通道激励-单通道接收的主动应力波探测；

将S阵列中各单元的N路接收信号线性叠加得到S路合成信号；

将S路合成信号进行时间反演放大，对实物结构进行N×S次独立单通道激励-单通道接收的时间反演探测；

将N阵列中各单元的S路接收信号线性叠加，得到N路合成时反信号；

采用合成时反继续对实物结构进行迭代探测，最终得到N路合成迭代信号；

采用基于时域聚焦的时间反演成像方法，对合成迭代信号进行时间反演放大处理，并各自在成像模型N阵列单元激励回传，实现损伤聚焦成像。

与现有技术相比，本发明的一个或多个实施例可以具有如下优点：

通过采用主动应力波方式进行结构损失进行检测，实现了对大范围结构损伤的长距离、实时主动监测；

通过合成时反技术，采用单通道激励－单通道接收系统即可等效合成多通道同步激励－多通道同步接收信号，使实物探测试验能够不受多通道同步激励－多通道同步接收系统的极高性能要求限制，实现探测信号增强，提高损伤目标散射信号分辨率；

通过应用迭代探测技术，实现了对结构最大损伤目标的聚焦增强检测，提高最大损伤目标检测准确性、最大损伤目标散射信号分辨率；

通过基于时域聚焦的时间反演镜成像方法，采用时域有限差分数值模拟实现了对损伤目标的聚焦成像，成像效果不受限于阵列单元个数和传播路径损耗。

附图说明

图1是基于合成时反的结构损伤迭代聚焦成像方法流程图；

图2是基于合成时反的结构损伤迭代检测模型图；

图3是二维仿真实施案例模型图；

图4a₁、4a₂、4b₁、4b₂、4c₁、4c₂、4d₁和4d₂是实施案例模型中S阵列各单元的接收信号及合成信号图；

图5是实施案例模型中S阵列合成信号的时间反演放大信号图；

图6a₁、6a₂、6b₁、6b₂、6c₁、6c₂、6d₁和6d₂是实施案例模型中N阵列各单元接收到的时反信号及合成时反信号图；

图7是案例模型中一次迭代探测的合成迭代信号图；

图8是案例模型中一次合成迭代信号的损伤聚焦成像图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述：

如图1所示，本实施例提供了一种基于合成时反的结构损伤迭代聚焦成像检测方法，该方法包括以下步骤：

步骤10对需采用N－S阵列同步探测的实物结构进行N×S次独立的单通道激励-单通道接收的主动应力波探测；

采用主动应力波探测技术初始激励N阵列各单元，N路初始激励信号均选用一致的宽带高斯脉冲信号。N－S探测阵列的阵列单元均采用基于正逆压电效应的PZT换能器，N阵列PZT换能器受初始激励电信号激励时可在被测结构上产生相应应力波。

步骤20将S阵列中各单元的N路接收信号线性叠加得到S路合成信号；

S阵列PZT换能器感应到结构中传播的应力波时可将应力波转化为电信号，将S阵列中各单元接收到的N路信号线性叠加，可得到S路合成信号，该S路合成信号与N阵列同步激励-S阵列同步接收时采集到的S路信号是等价的。

步骤30将S路合成信号进行时间反演放大，对实物结构进行N×S次独立单通道激励-单通道接收的时间反演探测；

对S路合成信号进行时间反演处理，考虑到S路合成信号相对N路初始激励信号的传播衰减，将各时间反演信号按同等倍数放大到激励水平，并各自回传到S阵列各单元中进行重新激励，从而实现对实物结构再次进行S×N次独立单通道激励－单通道接收的主动应力波时间反演探测。

步骤40将N阵列中各单元的S路接收信号线性叠加，得到N路合成时反信号；

根据合成时反思想，该N路合成时反信号与S阵列同步时反激励-N阵列同步接收时采集到的N路时反信号是等价的。合成时反处理实现了与传统同步时间反演一致的探测信号增强、损伤目标散射信号分辨率提高的效果。

步骤50采用合成时反继续对实物结构进行迭代探测，最终得到N路合成迭代信号；

采用上述合成时反思想继续对实物结构进行迭代探测，实现探测信号增强和最大损伤目标聚焦，最终得到N路合成迭代信号。所述迭代探测是指将N阵列上一次接收的合成时反信号时间反演放大处理后重新激励N阵列产生探测应力波，探测应力波经过结构传播后被S阵列感应接收，将S阵列的合成信号进行时间反演后放大到激励水平重新激励S阵列产生时反探测应力波，时反探测应力波经过结构传播后又被N阵列感应接收形成新的合成时反信号，以上就是一次迭代探测过程。

步骤60采用基于时域聚焦的时间反演成像方法，对合成迭代信号进行时间反演放大处理，并各自在成像模型N阵列单元激励回传，实现损伤聚焦成像。

所述成像模型是与无损实物结构参数一致的数值模拟模型，采用时域有限差分方法构建，该模型可利用已知空间介质网格尺寸、介质密度、刚度张量、传播速度以及发射阵列、接收阵列坐标等信息计算出空间中所有坐标点在时域上的应力信息和速度信息。采用基于时域聚焦的时间反演镜成像方法，对合成迭代信号再次进行时间反演放大处理，并各自在成像模型N阵列单元中激励回传，记录各合成迭代时反信号在模型中每一时刻的空间应力分布并进行归一化，把模型空间中同坐标点的同一时刻归一化信号进行相乘，并把所有时刻的结果相加得到损伤聚焦成像图形。所述损伤聚焦成像图形中坐标点的像素值表示所有合成迭代时反信号在该坐标点的归一化合成应力，对于多损伤结构，最大损伤目标处的合成应力最大，则损伤聚焦成像图形中像素值最大的区域即为最大损伤目标。

参见图2，为基于合成时反的结构损伤迭代检测模型。出于简化的目的，我们忽略阵列单元的时域非线性响应特性。设空间中有N个阵列单元P_n(1≤n≤N)，有S个阵列单元P_s(1≤s≤S)，有D个损伤目标P_d(1≤d≤D)，把各个损伤目标的散射过程简化为线性响应λ_d。N阵列单元与损伤目标之间的频域传递函数为H(P_n，P_d，ω)，损伤目标与S阵列单元之间的频域传递函数为H(P_d，P_s，ω)。

假设N个阵列单元的初始激励信号均为E(ω)。N阵列同步激励初始信号后，损伤目标P_d接收到的来自N阵列的播信号为

$S(P_d,\mathrm{\ω})=\underset{n=1}{\overset{N}{\mathrm{\Σ}}}E\left(\mathrm{\ω}\right)H(P_n,P_d,\mathrm{\ω})---\left(1\right)$

损伤目标P_d的散射信号为：

$S_e(P_d,\mathrm{\ω})={\mathrm{\λ}}_d\underset{n=1}{\overset{N}{\mathrm{\Σ}}}E\left(\mathrm{\ω}\right)H(P_n,P_d,\mathrm{\ω})---\left(2\right)$

当损伤目标的散射信号传播到S阵列时，P_s单元接收到的来自全体散射目标的信号为

$\begin{array}{c}S(P_s,\mathrm{\ω})=\underset{d=1}{\overset{D}{\mathrm{\Σ}}}{S}_e(P_d,\mathrm{\ω})H(P_d,P_s,\mathrm{\ω})\\ =\underset{d=1}{\overset{D}{\mathrm{\Σ}}}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{\λ}}_{d}E\left(\mathrm{\ω}\right)H(P_n,P_d,\mathrm{\ω})H(P_d,P_s,\mathrm{\ω})\end{array}---\left(3\right)$

S阵列接收到的来自全体散射目标的散射信号为

$\begin{array}{c}S(P_{s=1\→S},\mathrm{\ω})=\underset{s=1}{\overset{S}{\mathrm{\Σ}}}S(P_s,\mathrm{\ω})\\ =\underset{s=1}{\overset{S}{\mathrm{\Σ}}}\underset{d=1}{\overset{D}{\mathrm{\Σ}}}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{\λ}}_{d}E\left(\mathrm{\ω}\right)H(P_n,P_d,\mathrm{\ω})H(P_d,P_s,\mathrm{\ω})\end{array}---\left(4\right)$

由式(3)可知，对于P_s单元，其接收信号可以等效为N个项的线性叠加，每一项表示的意思是某一个N阵列单元单独激励信号时，P_s单元接收到的信号。由式(4)亦可知，整个S阵列的接收信号为N×S个类似项的线性叠加，而每一项表

$S_{\mathrm{TR}}(P_s,\mathrm{\ω})=\underset{d=1}{\overset{D}{\mathrm{\Σ}}}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{\λ}}_d{E}^{*}\left(\mathrm{\ω}\right)H^{*}(P_n,P_d,\mathrm{\ω})H^{*}(P_d,P_s,\mathrm{\ω})---\left(5\right)$

进行反演传播时，将时间反演信号重新加载在P_s单元上，则模型中任意点P_k接收到的来自P_s单元的时间反演传播信号为

$\begin{array}{c}{S}_{\mathrm{TR}}^{\left(s\right)}(P_s,\mathrm{\ω})=S_{\mathrm{TR}}(P_s,\mathrm{\ω})H(P_s,P_k,\mathrm{\ω})\\ =\underset{d=1}{\overset{D}{\mathrm{\Σ}}}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{\λ}}_d{E}^{*}\left(\mathrm{\ω}\right)H^{*}(P_n,P_d,\mathrm{\ω})H^{*}(P_d,P_s,\mathrm{\ω})H(P_k,P_s,\mathrm{\ω})\end{array}---\left(6\right)$

其中H(P_k，P_s，ω)为P_s单元与任意点P_k间的频域传递函数。则任意点P_k接收到的来自S阵列的时间反演传播信号为

$\begin{array}{c}{S}_{\mathrm{TR}}(P_k,\mathrm{\ω})=\underset{s=1}{\overset{S}{\mathrm{\Σ}}}{S}_{\mathrm{TR}}^{\left(s\right)}(P_k,\mathrm{\ω})\\ =\underset{s=1}{\overset{S}{\mathrm{\Σ}}}\underset{d=1}{\overset{D}{\mathrm{\Σ}}}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{\λ}}_d{E}^{*}\left(\mathrm{\ω}\right)H^{*}(P_n,P_d,\mathrm{\ω})H^{*}(P_d,P_s,\mathrm{\ω})H(P_k,P_s,\mathrm{\ω})\end{array}---\left(7\right)$

当P_k位于损伤目标P_d所处位置时，则此时式(7)可以改写为

$\begin{array}{c}{S}_{\mathrm{TR}}(P_k,\mathrm{\ω})=S_{\mathrm{TR}}(P_d,\mathrm{\ω})\\ =\underset{s=1}{\overset{S}{\mathrm{\Σ}}}\underset{d=1}{\overset{D}{\mathrm{\Σ}}}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{\λ}}_d{E}^{*}\left(\mathrm{\ω}\right)H^{*}(P_n,P_d,\mathrm{\ω})H^{*}(P_d,P_s,\mathrm{\ω})\\ =\underset{s=1}{\overset{S}{\mathrm{\Σ}}}\underset{d=1}{\overset{D}{\mathrm{\Σ}}}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{\λ}}_d{E}^{*}\left(\mathrm{\ω}\right)H^{*}(P_n,P_d,\mathrm{\ω}){\left|H(P_d,P_s,\mathrm{\ω})\right|}^2\end{array}---\left(8\right)$

根据传统时间反演镜理论，式(8)就是S阵列同步时间反演后聚焦于损伤目标P_d的聚焦信号，很明显，这个信号是S个类似于项的线性叠加，该项表示的意思是损伤目标接收到的来自P_s阵列单元的时间反演聚焦信号。式(8)等效于多个P_s阵列单元单独进行时间反演后聚焦到损伤目标的损伤聚焦信号的叠加结果，即可以通过单通道时间反演激励的线性叠加来达到多通道同步时间反演激励的探测效果。根据这个特性，在时间反演过程，可以不必像传统时间反演理论中介绍的对多路时间反演信号进行同步回传激励，取而代之可以采用多个阵列单元单独时间反演激励下的目标聚焦信号的线性叠加等效合成同步时间反演的目标聚焦信号。这就是合成时反思想。

接下来，我们在合成时反的基础上继续进行迭代探测。在式(8)基础上，P_n接收到的来自损伤目标P_d的合成时反信号为

$\begin{array}{c}{S}_{\mathrm{TR}}(P_n,\mathrm{\ω})={\underset{d=1}{\overset{D}{\mathrm{\Σ}}}{\mathrm{\λ}}_{d}S}_{\mathrm{TR}}(P_d,\mathrm{\ω})H(P_n,P_d,\mathrm{\ω})\\ =\underset{s=1}{\overset{S}{\mathrm{\Σ}}}\underset{d=1}{\overset{D}{\mathrm{\Σ}}}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}{{{\mathrm{\λ}}_d}^{2}E}^{*}\left(\mathrm{\ω}\right)H^{*}(P_n,P_d,\mathrm{\ω})H(P_n,P_d,\mathrm{\ω}){\left|H(P_d,P_s,\mathrm{\ω})\right|}^2\\ =\underset{s=1}{\overset{S}{\mathrm{\Σ}}}\underset{d=1}{\overset{D}{\mathrm{\Σ}}}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}{{\mathrm{\λ}}_d}^2{E}^{*}\left(\mathrm{\ω}\right){{\left|H(P_{n,}{P}_d,\mathrm{\ω})\right|}^2\left|H(P_d,P_s,\mathrm{\ω})\right|}^2\end{array}---\left(9\right)$

把如式(9)所示的N阵列所有合成时反信号进行时间反演处理，并各自重新激励在N阵列单元上进行第一次迭代探测。根据合成时反思想，则N阵列中P_n阵列单元接收回来的第一次合成迭代信号为

$S_{\mathrm{TR}1}(P_n,\mathrm{\ω})=\underset{s=1}{\overset{S}{\mathrm{\Σ}}}\underset{d=1}{\overset{D}{\mathrm{\Σ}}}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}{{\mathrm{\λ}}_d}^4{E}^{*}\left(\mathrm{\ω}\right){\left|H(P_n,P_d,\mathrm{\ω})\right|}^4{\left|H(P_d,P_s,\mathrm{\ω})\right|}^4---\left(10\right)$

则N阵列中P_n阵列单元接收回来的第(2i-1)(i＝1,2,···)次合成迭代信号为

$S_{\mathrm{TR}(2i-1)}(P_n,\mathrm{\ω})=E^{*}\left(\mathrm{\ω}\right)\underset{s=1}{\overset{S}{\mathrm{\Σ}}}\underset{d=1}{\overset{D}{\mathrm{\Σ}}}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}{{\mathrm{\λ}}_d}^{4i}{\left|H(P_n,P_d,\mathrm{\ω})\right|}^{4i}{\left|H(P_d,P_s,\mathrm{\ω})\right|}^{4i}---\left(11\right)$

N阵列中P_n阵列单元接收回来的第(2i)(i＝1,2,···)次合成迭代信号为

$S_{\mathrm{TR}\left(2i\right)}(P_n,\mathrm{\ω})=E^{*}\left(\mathrm{\ω}\right)\underset{s=1}{\overset{S}{\mathrm{\Σ}}}\underset{d=1}{\overset{D}{\mathrm{\Σ}}}\underset{n=1}{\overset{N}{\mathrm{\Σ}}}{{\mathrm{\λ}}_d}^{4i+2}{\left|H(P_n,P_d,\mathrm{\ω})\right|}^{4i+2}{\left|H(P_d,P_s,\mathrm{\ω})\right|}^{4i+2}---\left(12\right)$

由式(11)和式(12)可以看出，P_n阵列单元接收到的合成迭代信号主要是E^*(ω)、λ_d、H(P_n，P_d，ω)、H(P_d，P_s，ω)三个因子的乘积项叠加，该乘积项的实际意义是探测信号在存在损伤目标P_d的传播路径上的传播过程。实际测试中，由于N阵列和S阵列往往平行布置，因此处在该平行狭窄空间中的不同n、d、s取值下的传递函数H(P_n，P_d，ω)、H(P_d，P_s，ω)不会有巨大差异，而散射过程线性响应λ_d则因为损伤目标大小的不同会产生较大差异，损伤目标越大，其散射过程线性响应λ_d也越大。当迭代探测的次数越多时，不同损伤目标的散射过程线性响应差异将以2i指数级增长。因此，迭代探测次数的增加将使得最大损伤目标因子在合成迭代信号中占的比重越大，即最大损伤目标的分辨率越高，从而实现对最大损伤目标的聚焦探测。这就是合成迭代的思想。

结合本发明方法的内容提供仿真实施案例，参见图3，为二维仿真模型。该仿真模型由时域有限差分法构建，该模型网格步长为0.04mm，模型大小为800×800，模型中有两处损伤，大损伤目标中心坐标为(400,400)，尺寸为20×20，小损伤中心坐标为(300,500)，尺寸为10×10。仿真模型选用混凝土作为传播介质，损伤介质为水。N阵列由N1、N2、N3、N4单元组成，S阵列由S1、S2、S3、S4单元组成。初始激励采用中心频率为1MHz的一阶高斯调制信号，如步骤10所述，分别以该信号对N1、N2、N3、N4单元进行激励，并单独记录S1、S2、S3、S4单元的接收信号，共需进行16次主动应力波探测试验。

参见图4a₁、4a₂、4b₁、4b₂、4c₁、4c₂、4d₁和4d₂，为S阵列中各单元的接收信号及合成信号图。如步骤20所述，将S阵列中各单元的4路接收信号线性叠加得到合成信号。图4a展示了S1单元单独接收到的来自N1、N2、N3、N4阵列单元的4路接收信号，并将该4路接收信号合成一路可以等效N1、N2、N3、N4阵列单元同步发射激励信号时的合成信号。同理，图4a₂、4b₂、4c₂和4d₂展示了S2单元、S3单元、S4单元各自的4路独立接收信号和相应的合成信号。图4a₂、4b₂、4c₂和4d₂的合成信号中包含了大损伤信号、小损伤信号和基波信号，其中损伤信号相对基波信号的分辨率仍不是很高。

参见图5，为S阵列合成信号的时间反演放大信号图。如步骤30所述，将图4a₂、4b₂、4c₂和4d₂中S阵列的合成信号进行时间反演处理，并按同等倍数放大到激励水平，各自回传到S阵列各单元中进行重新激励，从而实现对仿真模型再次进行16次独立的单通道激励－单通道接收的主动应力波时间反演探测。

参见图6a₁、6a₂、6b₁、6b₂、6c₁、6c₂、6d₁和6d₂，为N阵列各单元接收到的时反信号及合成时反信号图。如步骤40所述，图6a₂、6b₂、6c₂和6d₂分别展示了N1单元、N2单元、N3单元和N4单元各自独立接收到的4路时反信号和相应的合成时反信号。由图6a₂、6b₂、6c₂和6d₂所示的合成时反信号明显看出，经过时间反演处理后大损伤目标分辨率相较基波信号和小损伤信号提高很多。

参见图7，为一次迭代探测的合成迭代信号图。如步骤50所述，在图6a₂、6b₂、6c₂和6d₂所示的合成时反信号基础上继续进行一次迭代探测，图7分别展示了N1单元、N2单元、N3单元和N4单元接收回来的一次合成迭代信号。由图7明显看出，大损伤目标信号分辨率得到增强。

参见图8，为一次合成迭代信号的损伤聚焦成像图。如步骤60所述，基于时域聚焦的时间反演镜成像方法，将图7所示的一次合成迭代信号进行时间反演放大处理，并各自在无损仿真模型中的N1、N2、N3、N4单元激励回传，记录各合成迭代时反信号在模型中每一时刻的空间应力分布并进行归一化，把同一时刻相同坐标点的各归一化信号进行相乘，并把所有时刻的结果相加得到损伤聚焦成像图形。由图8明显看出，像素值最大点位于原始大损伤目标位置中心，实现对大损伤目标的聚焦成像。

虽然本发明所揭露的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的专利保护范围，仍须以所附的权利要求书所界定的范围为准。

Claims (6)

1.基于合成时反的结构损伤迭代聚焦成像监测方法，其特征在于，所述方法包括：

对需采用N－S阵列同步探测的实物结构进行N×S次独立的单通道激励-单通道接收的主动应力波探测；

将S阵列中各单元的N路接收信号线性叠加得到S路合成信号；

将S路合成信号进行时间反演放大，对实物结构进行N×S次独立单通道激励-单通道接收的时间反演探测；

将N阵列中各单元的S路接收信号线性叠加，得到N路合成时反信号；

采用合成时反继续对实物结构进行迭代探测，最终得到N路合成迭代信号；

采用基于时域聚焦的时间反演成像方法，对合成迭代信号进行时间反演放大处理，并各自在成像模型N阵列单元激励回传，实现损伤聚焦成像；

所述迭代探测是指将N阵列上一次接收的合成时反信号时间反演放大处理后重新激励N阵列产生探测应力波，探测应力波经过结构传播后被S阵列感应接收，将S阵列的合成信号进行时间反演后放大到激励水平重新激励S阵列产生时反探测应力波，时反探测应力波经过结构传播后又被N阵列感应接收形成新的合成时反信号，完成一次迭代探测过程。

2.根据权利要求1所述的基于合成时反的结构损伤迭代聚焦成像监测方法，其特征在于，采用主动应力波探测技术初始激励N－S阵列中N阵列各单元时，N路初始激励信号均选用相同的宽带高斯脉冲信号。

3.根据权利要求1所述的基于合成时反的结构损伤迭代聚焦成像监测方法，其特征在于，所述N－S阵列的阵列单元均采用基于正逆压电效应的PZT换能器，PZT换能器受电信号激励时在被测结构上产生相应应力波，PZT换能器感应到结构中传播的应力波时，将应力波转化为电信号。

4.根据权利要求1所述的基于合成时反的结构损伤迭代聚焦成像监测方法，其特征在于，所述基于时域聚焦的时间反演成像方法包括：

对实物探测所得合成迭代信号再次进行时间反演放大处理，并各自在成像模型N阵列单元激励回传；

记录各合成迭代时反信号在模型中每一时刻的空间应力分布并进行归一化；

将同一时刻相同坐标点的各归一化信号进行相乘，并将所有时刻的结果相加得到损伤聚焦成像图形。

5.根据权利要求1或4所述的基于合成时反的结构损伤迭代聚焦成像监测方法，其特征在于，所述成像模型是与无损实物结构参数一致的数值模拟模型，采用时域有限差分方法构建，利用已知空间介质网格尺寸、介质密度、刚度张量、传播速度及发射阵列接收阵列坐标信息，计算出空间中所有坐标点在时域上的应力信息和速度信息。

6.根据权利要求4所述的基于合成时反的结构损伤迭代聚焦成像监测方法，其特征在于，所述损伤聚焦成像图形中坐标点的像素值表示所有合成迭代时反信号在该坐标点的归一化合成应力。