CN101923027B - 一种结构损伤检测系统、设备以及结构损伤检测方法 - Google Patents

Abstract

本发明涉及一种用于检测结构的损伤的结构损伤检测方法，包括步骤：使用敲击装置扫描结构的表面，并且在扫描的同时以具有预定频率的敲击力进行敲击，该预定频率接近所述结构的一个或者多个自然频率之一；利用传感器感测由所述待检测结构表面传递到敲击装置的响应信号；对所述传感器感测到的信号进行变换处理以获取在所述结构表面的每个位置处的信号谱图；从所述信号谱图中截取与所述敲击力的预定频段相对应的谱图包络线；计算所述结构表面的每个位置处的损伤指示值，该损伤指示值反映了该位置处的谱图包络线与其他位置处的谱图包络线的相似程度；以及将其损伤指示值发生突变的位置确定为存在结构损伤的位置。本发明还涉及利用该结构损伤检测方法的结构损伤检测系统以及其中使用的结构损伤检测设备。

Description

一种结构损伤检测系统、设备以及结构损伤检测方法

技术领域

本发明涉及结构损伤检测系统、设备及结构损伤检测方法，尤其涉及一种无损地对工程结构的损伤进行检测的结构损伤检测系统、设备以及结构损伤检测方法。

背景技术

工程结构在制造和使用过程中由于工艺原因或使用载荷的作用，往往会产生一些局部损伤，比如裂纹、凹坑、脱粘等等，这些损伤积累到一定程度就会影响结构的正常使用。为确保结构能够正常运行，在其交付使用前或服役一段时间后都需要对之进行损伤检测。比如桥梁在通车前需要进行成桥试验以确保达到设计要求，在使用期内需要进行定期的常规检测或不定期的特殊检测，以便及时发现结构中出现的损伤。再如机翼蒙皮等航空构件，只有其内部的损伤小于规定的程度才能交付使用，而使用一段时间后，又需要对这些构件重新进行损伤检测以确保今后的飞行安全。

通常对这些工程结构进行检测不能另外造成这些工程结构的损伤，因此，这种检测方法也称为结构无损检测技术。现有的结构无损检测技术大致可分为局部检测和整体检测两个类别：

(1)局部检测技术通过人的肉眼观察或借助于X射线、红外线、超声波、雷达和磁涡流等方式进行损伤探测。这类技术虽然对局部区域具有较高的损伤检测精度，但如果事先不知道损伤的大致区域，要用它们对工程结构的每个区域进行依次检测是非常费时的。此外，有些工程结构还存在利用局部检测技术无法检测到的死角。

(2)整体检测技术通过在工程结构上布设传感器来测量工程结构在外加激励下的静位移、速度、加速度等响应，并以此反演出工程结构中的损伤。如果事先在工程结构中预埋传感器，还可以实现对损伤情况的实时监测，但此时只能利用工程结构工作时的环境激励所激发出来的结构响应进行损伤反演。相对于局部检测技术，整体检测技术的检测效率较高。但是整体检测技术具有如下缺点。首先，如果外加激励施加不当，工程结构的响应可能无法充分反映出损伤情况，从而影响检测精度。比如，如果静载荷较小，工程结构中的裂纹不能充分张开，那么静位移就对这种损伤不敏感。另外，由于损伤是工程结构的局部特征，它对结构动态特性的影响往往只体现在较高阶频段上，因此如果外加激励不能激发出工程结构的高阶响应，那么损伤反演的精度就较差，因此导致整体检测所获得的结果精度较低。此外，整体法还存在着传感器的合理布设、传感器的耐久性、环境噪声的干扰以及信号的采集和传输等问题。

基于Cawley建立的敲击力学模型的敲击法可以检测复合材料层合板中的损伤。该敲击力学模型认为当力锤敲击到结构损伤部位时，由于局部刚度的降低使得该处的敲击力时程信号较敲击到完整部位时产生的敲击力信号，峰值较小，宽度较大。据此，计算敲击力频谱曲线和频率坐标轴所围的面积，可以得到一个特征参数R，通过比较完好和损伤部位的R值就可以识别出损伤。另外也可以直接比较损伤前后的敲击力时程信号的宽度以识别损伤。敲击法也是一种局部检测技术，它不需要把传感器固定在结构上，而且由于能够采集到结构局部刚度的变化，所以对损伤比较敏感。不过这种方法使用了整个频段的敲击力频谱信息，这就在一定程度上降低了损伤识别的精度，这是因为：(1)损伤往往体现在敲击力的较高阶频段上，使用整个频段信息会掩盖对损伤敏感频段的信息；(2)整个频段的信息中包含了环境噪声的成份。另外敲击力谱线所覆盖的频率范围是由结构刚度以及锤头的刚度决定的。为了获得高阶频段谱线，虽然可以通过使用硬度较高的锤头实现，但这并不能精确地控制频谱的范围。即便能够激发出高阶频谱，这些高频成份在整个频段上所占的比例往往较小，这也会影响损伤检测精度。最后，该方法要求事先知道完整结构的特征参数，这也对其使用带来了不便。现有的敲击式损伤检测仪器只能可靠地检测出结构2mm厚度内直径约大于10mm的损伤。

在桥梁损伤检测领域，杨永斌等人提出了一种利用行驶的车辆作为激励、通过车辆的加速度信号来提取桥梁基频的方法。这种方法的优点是不需要在桥梁上安装传感器，而且测量过程不会影响桥梁的正常运行。但是该方法没有在车辆上安装敲击装置，因此很难激发出桥梁的较高频响应，所得的加速度信号对局部损伤不敏感，而仅根据桥梁的基频是很难检测出结构的损伤的。

因此，期望提供一种不需要事先知道工程结构的完整特征信息、就可以方便、快速且准确地检测工程结构中的损伤的结构损伤检测方法及结构损伤检测设备。

发明内容

鉴于上述问题，提出了本发明以便提供一种克服上述问题或者至少部分地解决上述问题的结构损伤检测方法及结构损伤检测设备。

受上述行进车辆检测方法和敲击检测方法的启示，本发明提供了一种通过安装有敲击设备的可以行进的检测设备来识别待检测结构的损伤的敲击-扫描方法。通过为待检测结构设计特定的敲击力和检测设备行进方式，即使在充满噪声的环境中，检测设备也可以相当灵敏地检测出损伤，并且不需要待检测结构的完整特征信息，因此，非常便于工程实施。

根据本发明的一个方面，提供了一种利用结构损伤检测设备来检测结构损伤的方法，其中该结构损伤检测设备包括敲击装置和传感器，所述方法包括步骤：使用所述结构损伤检测设备沿着所述结构的表面进行扫描，其中在进行扫描的同时，所述结构损伤检测设备中的敲击装置以包括预定频率的敲击力进行敲击，该预定频率接近所述结构的自然频率；利用所述结构损伤检测设备中的传感器感测在所述结构表面的每个位置处、由所述结构表面传递到所述敲击装置的响应信号，其中所述传感器感测的信号可以为位移、速度和加速度中的至少之一；对所述传感器感测到的信号进行变换处理以获取在所述结构表面的每个位置处的信号谱图；从所述信号谱图中截取与所述敲击力的预定频率相对应的谱图包络线；计算所述结构表面的每个位置处的损伤指示值，该损伤指示值反映了该位置处的谱图包络线与其他位置处的谱图包络线的相似程度；以及将损伤指示值发生突变的位置确定为存在结构损伤的位置。

根据本发明的另一个方面，提供了一种结构损伤检测设备，包括：敲击装置，用于在待检测结构表面上以具有预定频率的敲击力进行敲击，其中所述预定频率接近所述待检测结构的自然频率；以及传感器，感测所述由待检测结构表面传递到敲击装置的响应信号，其中所述传感器感测的响应信号可以为位移、速度和加速度中的至少之一；其中所述敲击装置和所述传感器集成在一起，以便所述敲击装置在所述待检测结构表面上进行敲击的同时，所述传感器感测由所述待检测结构表面传递到所述敲击装置的响应信号。

根据本发明的还有一个方面，提供了一种结构损伤检测系统，包括：根据本发明的结构损伤检测设备；以及信号处理部件，用于对所述结构损伤检测设备中的传感器所感测到的信号进行处理以确定待检测结构中的损伤位置，所述信号处理部件包括：谱图获取装置，用于对所述传感器感测到的信号进行变换处理以获取在所述结构表面的每个位置处的信号谱图；谱图包络线截取装置，用于从所述信号谱图中截取与所述敲击力的预定频率相对应的谱图包络线；损伤指示值计算装置，用于计算所述结构表面的每个位置处的损伤指示值，该损伤指示值反映了该位置处的谱图包络线与其他位置处的谱图包络线的相似程度；以及损伤位置确定装置，用于基于所述结构表面的每个位置处的损伤指示值来确定所述结构中的损伤位置。

根据本发明的结构损伤检测方法、设备和系统可以设定敲击力的大小和频段，这样既有利于激发出损伤信息，又可以避开环境噪声的干扰，此外，本发明仅仅对由待检测结构表面传递到敲击装置的响应中与敲击力的频段相对应的频段的响应信号进行进一步处理，这可以提高检测精度。

另外，根据本发明的结构损伤检测方法、设备和系统不需要事先知道无损结构的完整特征信息，更便于实施；并且适用于在航空航天、土木和机械等任何可能发生结构材料的各种形式损伤的领域中应用。

附图说明

通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中：

图1是用于说明本发明原理的结构损伤检测设备的简化视图；

图2是根据本发明一个实施例的结构损伤检测系统100的示意图；

图3是根据本发明另一个实施例的结构损伤检测系统200的俯视示意图；

图4是根据本发明另一个实施例的结构损伤检测系统200的沿图3中的A向的正视示意图；

图5是MAC矩阵元素的分布图；以及

图6示意性地示出了根据本发明一个实施例的结构损伤检测方法的流程图。

具体实施例

下面结合附图和具体的实施方式对本发明作进一步的描述。

在具体描述本发明的实施例之前，首先结合附图阐述一下本发明的原理。

图1是用于说明本发明原理的结构损伤检测设备的简化模型。如图1所示，其中采用了一个简化模型来说明本发明所提出的敲击-扫描方案的理论基础。在这个模型中，将结构损伤检测设备建模为在具有k_V刚度的弹簧上的质量块M，并且该质量块上安装有两组相对的不同轴转轮

这些转轮中的每个从初始角θ_0i开始以恒定角速度(i＝1，2，…，N)旋转，以模拟敲击力的发生器。此外，在该模型中，还将待检测结构模拟为长度为L的二维简支梁，该简支梁被建模为具有弯曲刚度EI以及每单位长度质量为m的平面伯努利-欧拉(Bernoulli-Euler)梁。结构损伤检测设备以恒定速度v移动通过该梁。为了简化并仍然保持该问题的实质，这里暂时省略了其他的实际因素，如阻尼和表面粗糙度等。然而，如下面的数值模拟结果中所示的那样，即使在将这些实际因素考虑在内时，本发明所提出的敲击-扫描方案的仍然是有效的。

为了获得对结构损伤灵敏的检测设备动态响应，必须解决考虑了敲击力的检测设备与结构的交互作用问题。根据图1描述的简化模型，检测设备和梁(即待检测结构)的动力学方程可以描述为：

$M{\stackrel{..}{y}}_V+k(y_V-y_B{|}_{x=\mathrm{vt}})=F\left(t\right)---\left(1\right)$

$m{\stackrel{..}{y}}_B+\mathrm{EI}{y}_B^{\left(\mathrm{IV}\right)}=f\left(t\right)\mathrm{\δ}(x-\mathrm{vt})---\left(2\right)$

其中y_V和y_B分别是根据静态平衡位置而测量得到的检测设备和梁的位移。这里我们忽略由于检测设备和梁的重力而引起的梁变形。因此，当检测设备处于整个梁上的不同位置时，静态平衡位置总是处于相同水平。此外，符号上的双点表示相对于时间t的偏导数；

k表示在检测设备的弹簧刚度k_V和梁刚度k_B之间的串联等效刚度：

$k=\frac{k_V{k}_B}{k_V+k_B}---\left(3\right)$

可以将在接触点处的梁刚度或者弹性常数k_B计算为在接触点处作用的载荷除以在该点处产生的扰度的比值。使用材料力学的基本理论，可以容易得知在单位载荷下、在位置x处的无损伤的梁的挠度为

因此，在接触点x＝vt处的梁刚度k_B为：

$k_B=\frac{3\mathrm{EIL}}{{(L-x)}^2{x}^2}---\left(4\right)$

因为k_B沿着梁变化，所以等效刚度k也随着位置x而改变。

在方程(1)中的F(t)是施加到检测设备上的敲击力：

$F\left(t\right)=-2\stackrel{\‾}{m}r\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{\stackrel{\‾}{\mathrm{\ω}}}_i^2\cos ({\stackrel{\‾}{\mathrm{\ω}}}_{i}t+{\mathrm{\θ}}_{0i})---\left(5\right)$

以及f(t)δ(x-vt)是在检测设备和梁之间的接触力。这里，δ是克罗内克(Kronecker)delta函数，其指示接触力的移动。函数f(t)可以表述为：

$\left(t\right)=k(y_V-y_B{|}_{x=\mathrm{vt}})-2N\stackrel{\‾}{m}g-\mathrm{Mg}---\left(6\right)$

对于梁结构，可以将力学损伤定义为刚度或者阻尼的突然变化(这可以由局部阻抗来表示，或者可以有效地由瞬时刚度Z(x)来表示。根据方程(1)、(2)和(6)，可以得到：

$Z\left(x\right)=\frac{f\left(t\right)}{y_B\left(t\right)}\mathrm{\δ}(x-\mathrm{vt})=\frac{-M{\stackrel{..}{y}}_V\left(\frac{x}{v}\right)+F\left(\frac{x}{v}\right)-(2N\stackrel{\‾}{m}+M)g}{y_B\left(\frac{x}{v}\right)}---\left(7\right)$

上式可以被重写为：

${\stackrel{..}{y}}_V\left(x\right)=-\frac{y_B\left(\frac{x}{v}\right)}{M}Z\left(x\right)+\frac{F\left(\frac{x}{v}\right)}{M}-\frac{2N\stackrel{\‾}{m}+M}{M}g---\left(8\right)$

由于敲击力F和检测设备的质量

为常数，因此

相对于损伤的灵敏度系数中，y_B的值起着关键的作用。使用模态叠加方法，梁的位移可以表示为：

其中

是第j个模态以及q_Bj(t)是对应的模态坐标。

将方程(9)代入方程(2)中，在方程(2)两端乘以

并且在整个梁长度上进行积分，可以得到：

其中ω_Bj是梁的第j个自然角频率：

${\mathrm{\ω}}_{\mathrm{Bj}}={\left(\frac{\mathrm{j\π}}{L}\right)}^2\sqrt{\frac{\mathrm{EI}}{m}}---\left(11\right)$

将方程(1)、(6)以及

代入到方程(10)中，可以得出：

${\stackrel{..}{q}}_{\mathrm{Bj}}+{\mathrm{\ω}}_{\mathrm{Bj}}^2{q}_{\mathrm{Bj}}=-\frac{2}{\mathrm{mL}}[M{\stackrel{..}{y}}_V-F\left(t\right)+2N\stackrel{\‾}{m}g+\mathrm{Mg}]\sin \frac{\mathrm{j\πvt}}{L}---\left(12\right)$

通常

否则就不能确保在检测设备和梁之间的稳固接触。因此，方程(12)可以近似为：

${\stackrel{..}{q}}_{\mathrm{Bj}}+{\mathrm{\ω}}_{\mathrm{Bj}}^2{q}_{\mathrm{Bj}}\≈-\frac{2}{\mathrm{mL}}[-F\left(t\right)+(2N\stackrel{\‾}{m}+M)g]\sin \frac{\mathrm{j\πvt}}{L}---\left(13\right)$

利用方程(5)，可以通过Duhamel积分来求解方程(13)：

$q_{\mathrm{Bj}}=-\frac{2}{{\mathrm{\ω}}_{\mathrm{Bj}}\mathrm{mL}}{\∫}_0^t[2\stackrel{\‾}{m}r\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{\stackrel{\‾}{\mathrm{\ω}}}_i^2\cos ({\stackrel{\‾}{\mathrm{\ω}}}_i\mathrm{\τ}+{\mathrm{\θ}}_{0i})+(2N\stackrel{\‾}{m}+M)g]\sin \left(\frac{\mathrm{j\πv\τ}}{L}\right)\sin \left[{\mathrm{\ω}}_{\mathrm{Bj}}(t-\mathrm{\τ})\right]\mathrm{d\τ}$

$={\mathrm{\Δ}}_W\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{\left(R_j^i\right)}^2\left\{S_j\right[\frac{\sin ({\mathrm{\ω}}_{\mathrm{Bj}}t+{\mathrm{\θ}}_{0i})}{{(R_j^i-1)}^2-S_j^2}+\frac{\sin ({\mathrm{\ω}}_{\mathrm{Bj}}t-{\mathrm{\θ}}_{0i})}{{(R_j^i+1)}^2-S_j^2}]+\frac{\sin [({\stackrel{\‾}{\mathrm{\ω}}}_i+\frac{\mathrm{j\πv}}{L})t+{\mathrm{\θ}}_{0i}]}{{(R_j^i+S_j)}^2-1}$

$-\frac{\sin [({\stackrel{\‾}{\mathrm{\ω}}}_i-\frac{\mathrm{j\πv}}{L})t+{\mathrm{\θ}}_{0i}]}{{(R_j^i-S_j)}^2-1}\}+\frac{{\mathrm{\Δ}}_{\mathrm{Bj}}}{1-S_j^2}[\sin \left(\frac{\mathrm{j\πv}}{L}t\right)-S_j\sin \left({\mathrm{\ω}}_{\mathrm{Bj}}t\right)]$

(14)

其中

${\mathrm{\Δ}}_W=\frac{2\stackrel{\‾}{m}r}{\mathrm{mL}}---\left(15\right)$

${\mathrm{\Δ}}_{\mathrm{Bj}}=-\frac{2(2N\stackrel{\‾}{m}+M)g{L}^3}{{\left(\mathrm{j\π}\right)}^4\mathrm{EI}}---\left(16\right)$

$S_j=\frac{\mathrm{j\πv}}{{\mathrm{L\ω}}_{\mathrm{Bj}}}---\left(17\right)$

$R_j^i=\frac{{\stackrel{\‾}{\mathrm{\ω}}}_i}{{\mathrm{\ω}}_{\mathrm{Bj}}}---\left(18\right)$

在方程(14)中，由于初始速度和位移因为阻尼而快速衰减，因此可以被忽略。将方程(14)代入到方程(9)中，可以获得梁的位移：

$y_B(x,t)={\mathrm{\Δ}}_W\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}\sin \frac{\mathrm{j\πx}}{L}{\left(R_j^i\right)}^2\left\{S_j\right[\frac{\sin ({\mathrm{\ω}}_{\mathrm{Bj}}t+{\mathrm{\θ}}_{0i})}{{(R_j^i-1)}^2-S_j^2}+\frac{\sin ({\mathrm{\ω}}_{\mathrm{Bj}}t-{\mathrm{\θ}}_{0i})}{{(R_j^i+1)}^2-S_j^2}]$

$+\frac{\sin [({\stackrel{\‾}{\mathrm{\ω}}}_i+\frac{\mathrm{j\πv}}{L})t+{\mathrm{\θ}}_{0i}]}{{(R_j^i+S_j)}^2-1}-\frac{\sin [({\stackrel{\‾}{\mathrm{\ω}}}_i-\frac{\mathrm{j\πv}}{L})t+{\mathrm{\θ}}_{0i}]}{{(R_j^i-S_j)}^2-1}\}$

$+\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}\sin \frac{\mathrm{j\πx}}{L}\frac{{\mathrm{\Δ}}_{\mathrm{Bj}}}{1-S_j^2}[\sin \left(\frac{\mathrm{j\πv}}{L}t\right)-S_j\sin \left({\mathrm{\ω}}_{\mathrm{Bj}}t\right)]$

(19)

此外，将方程(5)和(19)代入到方程(1)中，得到：

${\stackrel{..}{y}}_V+{\mathrm{\ω}}_V^2{y}_V=G\left(t\right)---\left(20\right)$

其中ω_V是检测设备的自然角频率：

${\mathrm{\ω}}_V=\sqrt{\frac{k}{M}}---\left(21\right)$

以及

$G\left(t\right)={\mathrm{\ω}}_V^2{\mathrm{\Δ}}_W\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}\sin \left(\frac{\mathrm{j\πv}}{L}t\right){\left(R_j^i\right)}^2\left\{S_j\right[\frac{\sin ({\mathrm{\ω}}_{\mathrm{Bj}}t+{\mathrm{\θ}}_0)}{{(R_j^i-1)}^2-S_j^2}+\frac{\sin ({\mathrm{\ω}}_{\mathrm{Bj}}t-{\mathrm{\θ}}_0)}{{(R_j^i+1)}^2-S_j^2}]$

$+\frac{\sin [({\stackrel{\‾}{\mathrm{\ω}}}_i+\frac{\mathrm{j\πv}}{L})t+{\mathrm{\θ}}_0]}{{(R_j^i+S_j)}^2-1}-\frac{\sin [({\stackrel{\‾}{\mathrm{\ω}}}_i-\frac{\mathrm{j\πv}}{L})t+{\mathrm{\θ}}_0]}{{(R_j^i-S_j)}^2-1}\}$

$+{\mathrm{\ω}}_V^2\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}\sin \left(\frac{\mathrm{j\πv}}{L}t\right)\frac{{\mathrm{\Δ}}_{\mathrm{Bj}}}{1-S_j^2}[\sin \left(\frac{\mathrm{j\πv}}{L}t\right)-S_j\sin \left({\mathrm{\ω}}_{\mathrm{Bj}}t\right)]-\frac{2\stackrel{\‾}{m}r}{M}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{{\stackrel{\‾}{\mathrm{\ω}}}_i}^2\cos ({\stackrel{\‾}{\mathrm{\ω}}}_{i}t+{\mathrm{\θ}}_0)$

(22)

类似地，通过忽略检测设备的初始位移和速度，可以通过Duhamel积分来从方程(20)中求解检测设备的振动位移：

$y_V=\frac{1}{{\mathrm{\ω}}_V}{\∫}_0^{t}G\left(\mathrm{\τ}\right)\sin \left[{\mathrm{\ω}}_V(t-\mathrm{\τ})\right]\mathrm{d\τ}$

$={\mathrm{\ω}}_V\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}[\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{\Δ}}_W{\left(R_j^i\right)}^2{P}_j^i\left(t\right)+\frac{{\mathrm{\Δ}}_{\mathrm{Bj}}}{1-S_j^2}{Q}_j\left(t\right)]$

$-\frac{\stackrel{\‾}{m}r}{M{\mathrm{\ω}}_V}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{{\stackrel{\‾}{\mathrm{\ω}}}_i}^2[\frac{\cos ({\mathrm{\ω}}_{V}t+{\mathrm{\θ}}_{0i})}{{\stackrel{\‾}{\mathrm{\ω}}}_i-{\mathrm{\ω}}_V}-\frac{\cos ({\mathrm{\ω}}_{V}t-{\mathrm{\θ}}_{0i})}{{\stackrel{\‾}{\mathrm{\ω}}}_i+{\mathrm{\ω}}_V}-\frac{2{\mathrm{\ω}}_V\cos ({\stackrel{\‾}{\mathrm{\ω}}}_{i}t+{\mathrm{\θ}}_{0i})}{{{\stackrel{\‾}{\mathrm{\ω}}}_i}^2-{\mathrm{\ω}}_v^2}]---\left(23\right)$

其中

$P_j^i\left(t\right)=\frac{1}{2[{(R_j^i+S_j)}^2-1]}\{-\frac{\frac{\mathrm{j\πv}}{L}\cos ({\mathrm{\ω}}_{V}t+{\mathrm{\θ}}_{0i})}{{\left(\frac{\mathrm{j\πv}}{L}\right)}^2-{({\stackrel{\‾}{\mathrm{\ω}}}_i+\frac{\mathrm{j\πv}}{L}-{\mathrm{\ω}}_V)}^2}+\frac{\frac{\mathrm{j\πv}}{L}\cos ({\mathrm{\ω}}_{V}t-{\mathrm{\θ}}_{0i})}{{\left(\frac{\mathrm{j\πv}}{L}\right)}^2-{({\stackrel{\‾}{\mathrm{\ω}}}_i+\frac{\mathrm{j\πV}}{L}+{\mathrm{\ω}}_V)}^2}$

$+{\mathrm{\ω}}_V\frac{\cos [(\frac{j2\mathrm{\πv}}{L}+{\stackrel{\‾}{\mathrm{\ω}}}_i)t+{\mathrm{\θ}}_{0i}]}{{(\frac{j2\mathrm{\πv}}{L}+{\stackrel{\‾}{\mathrm{\ω}}}_i)}^2-{\mathrm{\ω}}_v^2}-{\mathrm{\ω}}_V\frac{\cos ({\stackrel{\‾}{\mathrm{\ω}}}_{i}t+{\mathrm{\θ}}_{0i})}{{{\stackrel{\‾}{\mathrm{\ω}}}_i}^2-{\mathrm{\ω}}_v^2}\}$

$-\frac{1}{2[{(R_j^i-S_j)}^2-1]}\{-\frac{\frac{\mathrm{j\πv}}{L}\cos ({\mathrm{\ω}}_{V}t+{\mathrm{\θ}}_{0i})}{{\left(\frac{\mathrm{j\πv}}{L}\right)}^2-{({\stackrel{\‾}{\mathrm{\ω}}}_i-\frac{\mathrm{j\πv}}{L}-{\mathrm{\ω}}_V)}^2}+\frac{\frac{\mathrm{j\πv}}{L}\cos ({\mathrm{\ω}}_{V}t-{\mathrm{\θ}}_{0i})}{{\left(\frac{\mathrm{j\πv}}{L}\right)}^2-{({\stackrel{\‾}{\mathrm{\ω}}}_i-\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_V)}^2}$

$+{\mathrm{\ω}}_V\frac{\cos ({\stackrel{\‾}{\mathrm{\ω}}}_{i}t+{\mathrm{\θ}}_{0i})}{{{\stackrel{\‾}{\mathrm{\ω}}}_i}^2-{\mathrm{\ω}}_v^2}-{\mathrm{\ω}}_V\frac{\cos [(\frac{j2\mathrm{\πv}}{L}-{\stackrel{\‾}{\mathrm{\ω}}}_i)t-{\mathrm{\θ}}_{0i}]}{{(\frac{j2\mathrm{\πv}}{L}-{\stackrel{\‾}{\mathrm{\ω}}}_i)}^2-{\mathrm{\ω}}_V^2}\}$

$+\frac{S_j}{2[{(R_j^i-1)}^2-S_j^2]}\{-\frac{\frac{\mathrm{j\πv}}{L}\cos ({\mathrm{\ω}}_{V}t+{\mathrm{\θ}}_{0i})}{{\left(\frac{\mathrm{j\πv}}{L}\right)}^2-{({\mathrm{\ω}}_{\mathrm{Bj}}-{\mathrm{\ω}}_V)}^2}+\frac{\frac{\mathrm{j\πv}}{L}\cos ({\mathrm{\ω}}_{V}t-{\mathrm{\θ}}_{0i})}{{\left(\frac{\mathrm{j\πv}}{L}\right)}^2-{({\mathrm{\ω}}_{\mathrm{Bj}}+{\mathrm{\ω}}_V)}^2}$

$+{\mathrm{\ω}}_V\frac{\cos [(\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})t+{\mathrm{\θ}}_{0i}]}{{(\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})}^2-{\mathrm{\ω}}_V^2}-{\mathrm{\ω}}_V\frac{\cos [(-\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})t+{\mathrm{\θ}}_{0i}]}{{(\frac{\mathrm{j\πv}}{L}-{\mathrm{\ω}}_{\mathrm{Bj}})}^2-{\mathrm{\ω}}_V^2}\}$

$-\frac{S_j}{2[{(R_j^i+1)}^2-S_j^2]}\{-\frac{\frac{\mathrm{j\πv}}{L}\cos ({\mathrm{\ω}}_{V}t-{\mathrm{\θ}}_{0i})}{{\left(\frac{\mathrm{j\πv}}{L}\right)}^2-{({\mathrm{\ω}}_{\mathrm{Bj}}-{\mathrm{\ω}}_V)}^2}+\frac{\frac{\mathrm{j\πv}}{L}\cos ({\mathrm{\ω}}_{V}t+{\mathrm{\θ}}_{0i})}{{\left(\frac{\mathrm{j\πv}}{L}\right)}^2-{({\mathrm{\ω}}_{\mathrm{Bj}}+{\mathrm{\ω}}_V)}^2}$

$+{\mathrm{\ω}}_V\frac{\cos [(\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})t-{\mathrm{\θ}}_{0i}]}{{(\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})}^2-{\mathrm{\ω}}_V^2}-{\mathrm{\ω}}_V\frac{\cos [(-\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})t-{\mathrm{\θ}}_{0i}]}{{(\frac{\mathrm{j\πv}}{L}-{\mathrm{\ω}}_{\mathrm{Bj}})}^2-{\mathrm{\ω}}_V^2}\}$

(24)

$Q_j\left(t\right)=\frac{{\mathrm{\ω}}_V\cos \left(\frac{j2\mathrm{\πv}}{L}t\right)}{2[{\left(\frac{j2\mathrm{\πv}}{L}\right)}^2-{\mathrm{\ω}}_V^2]}+\frac{1}{2{\mathrm{\ω}}_V}-\frac{2{\left(\frac{\mathrm{j\πv}}{L}\right)}^2\cos \left({\mathrm{\ω}}_{V}t\right)}{{\mathrm{\ω}}_V[{\left(\frac{j2\mathrm{\πv}}{L}\right)}^2-{\mathrm{\ω}}_V^2]}$

$-S_j{\mathrm{\ω}}_V\{\frac{\cos \left[(\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})t\right]}{2[{(\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})}^2-{\mathrm{\ω}}_V^2]}-\frac{\cos \left[(\frac{\mathrm{j\πv}}{L}-{\mathrm{\ω}}_{\mathrm{Bj}})t\right]}{2[{(\frac{\mathrm{j\πv}}{L}-{\mathrm{\ω}}_{\mathrm{Bj}})}^2-{\mathrm{\ω}}_V^2]}$

$+\frac{2\frac{\mathrm{j\πv}}{L}{\mathrm{\ω}}_{\mathrm{Bj}}\cos \left({\mathrm{\ω}}_{V}t\right)}{[{(\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})}^2-{\mathrm{\ω}}_V^2][{(\frac{\mathrm{j\πv}}{L}-{\mathrm{\ω}}_{\mathrm{Bj}})}^2-{\mathrm{\ω}}_V^2]}\}---\left(25\right)$

因此，检测设备的加速度为：

${\stackrel{..}{y}}_V={\mathrm{\ω}}_V\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}[\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{\Δ}}_W{\left(R_j^i\right)}^2{\stackrel{..}{P}}_j^i\left(t\right)+\frac{{\mathrm{\Δ}}_{\mathrm{Bj}}}{1-S_j^2}{\stackrel{..}{Q}}_j\left(t\right)]$

$+\frac{\stackrel{\‾}{m}r}{M}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{{\stackrel{\‾}{\mathrm{\ω}}}_i}^2[\frac{{\mathrm{\ω}}_V\cos ({\mathrm{\ω}}_{V}t+{\mathrm{\θ}}_{0i})}{{\stackrel{\‾}{\mathrm{\ω}}}_i-{\mathrm{\ω}}_V}-\frac{{\mathrm{\ω}}_V\cos ({\mathrm{\ω}}_{V}t-{\mathrm{\θ}}_{0i})}{{\stackrel{\‾}{\mathrm{\ω}}}_i+{\mathrm{\ω}}_V}-\frac{2{{\stackrel{\‾}{\mathrm{\ω}}}_i}^2\cos ({\stackrel{\‾}{\mathrm{\ω}}}_{i}t+{\mathrm{\θ}}_{0i})}{{{\stackrel{\‾}{\mathrm{\ω}}}_i}^2-{\mathrm{\ω}}_V^2}]---\left(26\right)$

其中

和

可以容易地从方程(24)和(25)获得：

${\stackrel{..}{P}}_j^i\left(t\right)=\frac{1}{2[{(R_j^i+S_j)}^2-1]}\{\frac{\frac{\mathrm{j\πv}}{L}{{\mathrm{\ω}}_V}^2\cos ({\mathrm{\ω}}_{V}t+{\mathrm{\θ}}_{0i})}{{\left(\frac{\mathrm{j\πv}}{L}\right)}^2-{({\stackrel{\‾}{\mathrm{\ω}}}_i+\frac{\mathrm{j\πv}}{L}-{\mathrm{\ω}}_V)}^2}-\frac{\frac{\mathrm{j\πv}}{L}{{\mathrm{\ω}}_V}^2\cos ({\mathrm{\ω}}_{V}t-{\mathrm{\θ}}_{0i})}{{\left(\frac{\mathrm{j\πv}}{L}\right)}^2-{({\stackrel{\‾}{\mathrm{\ω}}}_i+\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_V)}^2}$

$-{\mathrm{\ω}}_V{(2\frac{\mathrm{j\πv}}{L}+{\stackrel{\‾}{\mathrm{\ω}}}_i)}^2\frac{\cos [(2\frac{\mathrm{j\πv}}{L}+{\stackrel{\‾}{\mathrm{\ω}}}_i)t+{\mathrm{\θ}}_{0i}]}{{(2\frac{\mathrm{j\πv}}{L}+{\stackrel{\‾}{\mathrm{\ω}}}_i)}^2-{\mathrm{\ω}}_V^2}+{\mathrm{\ω}}_V{{\stackrel{\‾}{\mathrm{\ω}}}_i}^2\frac{\cos ({\stackrel{\‾}{\mathrm{\ω}}}_{i}t+{\mathrm{\θ}}_{0i})}{{{\stackrel{\‾}{\mathrm{\ω}}}_i}^2-{\mathrm{\ω}}_V^2}\}$

$-\frac{1}{2[{(R_j^i-S_j)}^2-1]}\{\frac{\frac{\mathrm{j\πv}}{L}{{\mathrm{\ω}}_V}^2\cos ({\mathrm{\ω}}_{V}t+{\mathrm{\θ}}_{0i})}{{\left(\frac{\mathrm{j\πv}}{L}\right)}^2-{({\stackrel{\‾}{\mathrm{\ω}}}_i-\frac{\mathrm{j\πv}}{L}-{\mathrm{\ω}}_V)}^2}-\frac{\frac{\mathrm{j\πv}}{L}{{\mathrm{\ω}}_V}^2\cos ({\mathrm{\ω}}_{V}t-{\mathrm{\θ}}_{0i})}{{\left(\frac{\mathrm{j\πv}}{L}\right)}^2-{({\stackrel{\‾}{\mathrm{\ω}}}_i-\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_V)}^2}$

$+{\mathrm{\ω}}_V{(2\frac{\mathrm{j\πv}}{L}-{\stackrel{\‾}{\mathrm{\ω}}}_i)}^2\frac{\cos [(2\frac{\mathrm{j\πv}}{L}-{\stackrel{\‾}{\mathrm{\ω}}}_i)t-{\mathrm{\θ}}_{0i}]}{{(\frac{\mathrm{j\π}2v}{L}-{\stackrel{\‾}{\mathrm{\ω}}}_i)}^2-{\mathrm{\ω}}_V^2}-{\mathrm{\ω}}_V{{\stackrel{\‾}{\mathrm{\ω}}}_i}^2\frac{\cos ({\stackrel{\‾}{\mathrm{\ω}}}_{i}t+{\mathrm{\θ}}_{0i})}{{{\stackrel{\‾}{\mathrm{\ω}}}_i}^2-{\mathrm{\ω}}_V^2}\}$

$+\frac{S_j}{2[{(R_j^i-1)}^2-S_j^2]}\{\frac{\frac{\mathrm{j\πv}}{L}{{\mathrm{\ω}}_V}^2\cos ({\mathrm{\ω}}_{V}t+{\mathrm{\θ}}_{0i})}{{\left(\frac{\mathrm{j\πv}}{L}\right)}^2-{({\mathrm{\ω}}_{\mathrm{Bj}}-{\mathrm{\ω}}_V)}^2}-\frac{\frac{\mathrm{j\πv}}{L}{{\mathrm{\ω}}_V}^2\cos ({\mathrm{\ω}}_{V}t-{\mathrm{\θ}}_{0i})}{{\left(\frac{\mathrm{j\πv}}{L}\right)}^2-{({\mathrm{\ω}}_{\mathrm{Bj}}+{\mathrm{\ω}}_V)}^2}$

$-{\mathrm{\ω}}_V{(\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})}^2\frac{\cos [(\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})t+{\mathrm{\θ}}_{0i}]}{{(\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})}^2-{\mathrm{\ω}}_V^2}+{\mathrm{\ω}}_V{(\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})}^2\frac{\cos [(-\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})t+{\mathrm{\θ}}_{0i}]}{{(\frac{\mathrm{j\πv}}{L}-{\mathrm{\ω}}_{\mathrm{Bj}})}^2-{\mathrm{\ω}}_V^2}\}$

$-\frac{S_j}{2[{(R_j^i+1)}^2-S_j^2]}\{\frac{\frac{\mathrm{j\πv}}{L}{{\mathrm{\ω}}_V}^2\cos ({\mathrm{\ω}}_{V}t-{\mathrm{\θ}}_{0i})}{{\left(\frac{\mathrm{j\πv}}{L}\right)}^2-{({\mathrm{\ω}}_{\mathrm{Bj}}-{\mathrm{\ω}}_V)}^2}-\frac{\frac{\mathrm{j\πv}}{L}{{\mathrm{\ω}}_V}^2\cos ({\mathrm{\ω}}_{V}t+{\mathrm{\θ}}_{0i})}{{\left(\frac{\mathrm{j\πv}}{L}\right)}^2-{({\mathrm{\ω}}_{\mathrm{Bj}}+{\mathrm{\ω}}_V)}^2}$

$-{\mathrm{\ω}}_V{(\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})}^2\frac{\cos [(\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})t-{\mathrm{\θ}}_{0i}]}{{(\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})}^2-{\mathrm{\ω}}_V^2}+{\mathrm{\ω}}_V{(-\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})}^2\frac{\cos [(-\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})t-{\mathrm{\θ}}_{0i}]}{{(\frac{\mathrm{j\πv}}{L}-{\mathrm{\ω}}_{\mathrm{Bj}})}^2-{\mathrm{\ω}}_V^2}\}$

(26A)

${\stackrel{..}{Q}}_j\left(t\right)=\frac{2{\mathrm{\ω}}_V{\left(\frac{\mathrm{j\πv}}{L}\right)}^2[\cos \left({\mathrm{\ω}}_{V}t\right)-\cos \left(\frac{j2\mathrm{\πv}}{L}t\right)]}{{\left(\frac{j2\mathrm{\πv}}{L}\right)}^2-{\mathrm{\ω}}_V^2}$

$-S_j{\mathrm{\ω}}_V\{-{(\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})}^2\frac{\cos \left[(\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})t\right]}{2[{(\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})}^2-{\mathrm{\ω}}_V^2]}+{(\frac{\mathrm{j\πv}}{L}-{\mathrm{\ω}}_{\mathrm{Bj}})}^2\frac{\cos \left[(\frac{\mathrm{j\πv}}{L}-{\mathrm{\ω}}_{\mathrm{Bj}})t\right]}{2[{(\frac{\mathrm{j\πv}}{L}-{\mathrm{\ω}}_{\mathrm{Bj}})}^2-{\mathrm{\ω}}_V^2]}$

$-\frac{2\frac{\mathrm{j\πv}}{L}{\mathrm{\ω}}_{\mathrm{Bj}}{{\mathrm{\ω}}_V}^2\cos \left({\mathrm{\ω}}_{V}t\right)}{[{(\frac{\mathrm{j\πv}}{L}+{\mathrm{\ω}}_{\mathrm{Bj}})}^2-{\mathrm{\ω}}_V^2][{(\frac{\mathrm{j\πv}}{L}-{\mathrm{\ω}}_{\mathrm{Bj}})}^2-{\mathrm{\ω}}_V^2]}\}$

(26B)

当待检测结构出现局部损伤时，梁的全局属性不会改变太多，也就是说，当梁存在局部小损伤时，梁的弯曲刚度EI和每单位长度的质量m几乎不会改变。因此，方程(26)中的参数ω_Bj，Δ_W，Δ_Bj，S_j和

几乎不会因为小损伤而改变。然而，梁刚度k_B会沿着梁的长度方向发生改变。局部损伤对在损伤附近的k_B有很大的影响。考虑方程(3)和(21)，可以看出检测设备的自然角频率ω_V对损伤很灵敏。因此，检测设备的加速度对损伤的灵敏度可以计算为：

$\frac{\∂{\stackrel{..}{y}}_V}{\∂{\mathrm{\ω}}_V}=\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}[\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{\Δ}}_W{\left(R_j^i\right)}^2{\stackrel{..}{P}}_j^i+\frac{{\mathrm{\Δ}}_{\mathrm{Bj}}}{1-S_j^2}{\stackrel{..}{Q}}_j]+{\mathrm{\ω}}_V\underset{j=1}{\overset{\∞}{\mathrm{\Σ}}}[\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{\Δ}}_W{\left(R_j^i\right)}^2\frac{\∂{\stackrel{..}{P}}_j^i}{\∂{\mathrm{\ω}}_V}+\frac{{\mathrm{\Δ}}_{\mathrm{Bj}}}{1-S_j^2}\frac{\∂{\stackrel{..}{Q}}_j}{\∂{\mathrm{\ω}}_V}]$

$+\frac{\stackrel{\‾}{m}r}{M}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{{\stackrel{\‾}{\mathrm{\ω}}}_i}^2[\frac{\cos ({\mathrm{\ω}}_{V}t+{\mathrm{\θ}}_{0i})}{{\stackrel{\‾}{\mathrm{\ω}}}_i-{\mathrm{\ω}}_V}-\frac{{\mathrm{\ω}}_V\sin ({\mathrm{\ω}}_{V}t+{\mathrm{\θ}}_{0i})t}{{\stackrel{\‾}{\mathrm{\ω}}}_i-{\mathrm{\ω}}_V}+\frac{{\mathrm{\ω}}_V\cos ({\mathrm{\ω}}_{V}t+{\mathrm{\θ}}_{0i})}{{({\stackrel{\‾}{\mathrm{\ω}}}_i-{\mathrm{\ω}}_V)}^2}$

$-\frac{\cos ({\mathrm{\ω}}_{V}t-{\mathrm{\θ}}_{0i})}{{\stackrel{\‾}{\mathrm{\ω}}}_i+{\mathrm{\ω}}_V}+\frac{{\mathrm{\ω}}_V\sin ({\mathrm{\ω}}_{V}t-{\mathrm{\θ}}_{0i})t}{{\stackrel{\‾}{\mathrm{\ω}}}_i+{\mathrm{\ω}}_V}+\frac{{\mathrm{\ω}}_V\cos ({\mathrm{\ω}}_{V}t-{\mathrm{\θ}}_{0i})}{{({\stackrel{\‾}{\mathrm{\ω}}}_i+{\mathrm{\ω}}_V)}^2}$

$-\frac{4{\mathrm{\ω}}_V{{\stackrel{\‾}{\mathrm{\ω}}}_i}^2\cos ({\stackrel{\‾}{\mathrm{\ω}}}_{i}t+{\mathrm{\θ}}_{0i})}{{({{\stackrel{\‾}{\mathrm{\ω}}}_i}^2-{\mathrm{\ω}}_V^2)}^2}]$

(27)

其中和

可以分别容易地从方程(26A)和(26B)中获得。根据方程(26)、(27)和(15-18)，可以发现，当

即，

时可以增强检测设备的加速度对损伤的灵敏度。这意味着最好使得敲击频率接近待检测结构的自然频率，然而，根据方程(19)，还必须注意避免待检测结构的发生共振现象。

此外，通过将上述方程(26)转换到频率域中，可以发现，该方程中涉及到以下7个频率，即ω_V、

和

由于

的值通常非常小，因此，其中仅仅ω_V、

和ω_Bj是主要的频率。如上所述，结构损伤会导致ω_V的改变，而ω_V的改变又会改变在方程(26)中的相应频率相关项的系数。也就是说，在将上述方程(26)转换到频率域中之后，各个相应频率处的幅值中包含了结构损伤信息，通过分析这些幅值，就可以获得有关损伤的信息。本发明基于上述分析而做出。

图2示出了根据本发明一个实施例的结构损伤检测系统200。如图2所示，结构损伤检测系统200包括结构损伤检测设备210和信号处理部件220。结构损伤检测设备210用于在待检测结构表面上进行扫描，在扫描的同时施加敲击，并感测由待检测结构表面传递到结构损伤检测设备的响应信号。信号处理部件220对结构损伤检测设备210所感测到的响应信号进行处理以确定待检测结构中有无损伤以及损伤的位置。

结构损伤检测设备210包括敲击装置212和传感器214。敲击装置212用于在当结构损伤检测设备210在待检测结构的表面上进行扫描的同时以具有一定频率和大小的力敲击该待检测结构。如上理论分析所述，敲击装置212所施加的敲击力的频率应当接近待检测结构的一个或者多个自然频率之一，以便高灵敏度地检测出结构损伤。根据本发明的一个实施例，离散频率接近待检测结构的某个自然频率意指将离散频率设置为待检测结构的该自然频率的±20％范围之内的某个值，更可取的是±10％范围之内的某个值，以及±5％范围之内的某个值。可选地，敲击力可以具有多个离散频率，而多个离散频率中的至少一个频率接近待检测结构的多个自然频率中的至少一个。进一步可选地是，敲击力可以具有一个或者多个频段，而所述一个或者多个频段中的至少一个涵盖了接近待检测结构的多个自然频率的至少一个的值。同时，为了防止出现待检测结构的共振，可以不将敲击力的频率设置为刚好等于待检测结构的自然频率，或者敲击力的频率设置为在某个范围内快速变化，从而使得待检测结构发生共振的时间很短而不会导致待检测结构的损伤。

传感器214感测由待检测结构表面传递到敲击装置的响应。一般而言，待检测结构表面在敲击力的作用下会产生位移、速度和加速度响应，由于待检测结构表面和敲击装置的相互作用，这些响应会传递到敲击装置上，其瞬时值的大小反映待检测结构在敲击位置处的结构特征，因此，当待检测结构中出现损伤时，该位置处由待检测结构表面传递到敲击装置的位移、速度或加速度响应会与没有损伤位置的响应明显不同。传感器214可以是能够感测由待检测结构表面传递到敲击装置的响应的任何传感器。例如，传感器214可以是能够感测由待检测结构表面传递到敲击装置的位移、速度和加速度响应之一的任何传感器。根据本发明的一个实施例，传感器214为用于感测由待检测结构表面传递到敲击装置的加速度响应的加速度传感器，这是因为，在传感器114所感测的加速度a和结构位移响应b、结构局部阻抗Z以及敲击力F之间的关系存在如下关系：

$a\left(x\right)\≈-\frac{b\left(x\right)}{M}Z\left(x\right)+\frac{F\left(x\right)}{M}-g$

其中，x代表敲击的位置，M是敲击装置的质量，g是重力加速度。

为了精确地感测出在敲击装置212的敲击位置处被检测结构的响应，可以将传感器214和敲击装置212集成在一起。

另外，结构损伤检测设备210还包括承载装置216，用于承载传感器214和敲击装置212，并且在被检测结构的表面上进行扫描。

信号处理部件220包括谱图获取装置222、谱图包络线截取装置224、损伤指示值计算装置226和损伤位置确定装置228。

谱图获取装置222获取传感器214所感测到的传感器信号。根据检测精度的要求，可以将待检测结构的表面划分为具有一定大小的多个部分，而结构损伤检测设备210扫描每个表面部分都需要一定时间，传感器214所获取的信号为在该段时间上的传感器信号分布。随后，谱图获取装置222对该传感器信号在时间上的分布进行变换处理以获取该信号在频率域或者尺度域上的表示。这种变换处理可以利用现有技术领域的任何变换处理进行，例如该变换处理可以是短时傅立叶变换、小波变换或者Hilbert-Huang变换等。因此，谱图获取装置122获取了在所述结构表面的每个位置部分处的信号频谱或尺度谱。

谱图包络线截取装置224从所述信号谱图中截取与敲击装置的敲击力频率相对应的谱图包络线。由于敲击力的至少一个频率接近待检测结构的至少一个自然频率，因此，谱图包络线截取装置224所截取的谱图包络线应当涵盖待检测结构的至少一个自然频率。在敲击力具有某个频段的情况下，所截取的谱图包络线具有与该敲击力的频段相对应的频率范围。上面所给出的、有关所截取的谱图包络线的频率范围的示例仅仅为示范性的。本领域技术人员应当理解，只要所截取的谱图包络线的频段与敲击装置的敲击频率相对应，并且适于在后续数据处理中使用，则所有的频率范围都在本发明的保护范围之内。由于仅仅对与敲击频段相对应的谱图进行分析，可以减少由于处理全部频段谱图所需要的处理开销，同时不影响检测精度。另外由于环境噪声的频段一般不同于被检测结构的自然频率，因此不同于敲击装置的敲击频率，所以还可以显著减少环境噪声对检测结果的影响，从而提高检测的精度。

损伤指示值计算装置226基于所截取的谱图包络线来计算待检测结构表面的每个位置处的损伤指示值，该损伤指示值反映了该位置处的谱图包络线与其他位置处的谱图包络线的相似程度。

可以有多种方式来计算损伤指示值，根据本发明的一个实施例，可以如下所述来计算损伤指示值，首先将谱图包络线转换为谱图向量，这可以通过获取谱图包络线中的多个频率处的幅值作为谱图向量的分量来进行转换。谱图向量的分量数量可以根据检测精度以及系统的处理性能等确定。分量对应的频率可以均匀分布在谱图包络线的频段中，也可以不均匀分布。但是，对于待检测结构而言，每个位置处的频率选择方式应当是相同的。

在生成了每个位置处的谱图向量之后，损伤指示值计算装置226计算反映当前位置处的谱图向量和其他位置处的谱图向量的相似程度的损伤指示值。可以有多种方式来计算两个谱图向量之间的相似程度。根据本发明的一个实施例，可以将损伤指示值计算为：

$\left(\underset{j=1...n}{\mathrm{\Σ}}\frac{Y_i\·Y_j}{\left|Y_i\right|\×\left|Y_j\right|}\right)/n,$

其中，Y_i和Y_j分别表示在所述结构上的当前位置和其他位置处的谱图向量，·表示谱图向量的内积运算，|Y_i+|和|Y_j|表示向量Y_i和Y_j的长度，n为待检测结构所划分的单位检测部分的数量。

在Y_i和Y_j越相似，即包络线i和包络线j的形状越相似时，就越接近1。因此，针对某个谱图向量Y_i的损伤指示值可以是该谱图向量Y_i和所有其他谱图向量Y_j相似度的总和的平均值。

根据本发明的另一个实施例，可以利用MAC系数来获取损伤指示值。MAC系数矩阵的计算公式为：

$\mathrm{MAC}(i,j)=\frac{Y_i\·Y_j}{\left|Y_i\right|\×\left|Y_j\right|},$

其中，Y_i和Y_j分别表示在所述结构上第i个和第j个部分的谱图向量，·表示谱图向量的内积运算，|Y_i|和|Y_j|表示向量Y_i和Y_j的长度。MAC系数矩阵中的每个元素都表示了两个谱图向量之间的相似程度，其中主对角线上的元素一定等于1。MAC系数矩阵的第i行或第i列元素的大小就能反映出对应于某个谱图向量Y_i的损伤情况，因而也就是该谱图向量Y_i所对应位置的损伤指示值。

损伤位置确定装置228基于待检测结构表面的每个位置处的损伤指示值来确定待检测结构中的损伤位置。例如，损伤位置确定装置128可以将损伤指示值发生突然下降的位置确定为存在损伤的位置。因为如果结构没有损伤，那么损伤指示值在空间的分布是比较光滑的；如果某处的损伤指示值发生了突变，那么就意味着该处局部阻抗和其他地方明显不同，很可能出现了损伤，而突变的大小也就反映了损伤的严重程度。比如，取MAC系数矩阵第一行元素的值为纵轴元素，并以所对应的结构位置为横轴画一条曲线，如果结构没有损伤，那么该曲线应该是比较光滑的；如果曲线在某处突然下降，就说明该处存在损伤，下降的幅度越大，损伤也就越严重。根据本发明的一个实施例，可以将曲线的一阶导数值的绝对值超过预定阈值的位置确定为存在损伤的位置。

图3和4示出了根据本发明另一个实施例的、适于对类似于桥梁的结构进行结构损伤检测的损伤检测系统300的示意图。在该实施例中，待检测结构为轴向长度远大于垂直于轴向的截面的板结构，该板结构沿轴向的横截面都非常相似。这种结构在工程实践中较为常见，例如公路桥梁，铁路桥梁等都具有这种结构。而且在工程实践中，这类桥梁也经常需要进行检测以确保桥梁没有损伤，或者在出现损伤后及时加以修复以防止由于桥梁断裂而出现事故。

图3和4所示的损伤检测系统300非常适合于上述类似桥梁结构的损伤检测。如图3所示，待检测结构为层合板5，其固定在地面上的三根钢棒4上，并且离地具有一定间隙以保证结构损伤检测设备310在层合板5上移动时层合板5不接触地面，从而层合板5可以很好地模拟现实中的桥梁结构。另外，为了测试检测结构，在层合板5中制造了损伤6，该损伤表现为锯缝。结构损伤检测设备310的承载装置为钢制小车1，其中的敲击装置3可以在来自计算机11的控制信号的控制下产生具有预定频段的敲击力。根据本发明的一个实施方式，该敲击力为在20-100Hz之间快速扫描的等幅正弦控制信号所产生的、频段在20-100Hz之间的敲击力。层合板5的前两阶自然频率通常在20-100Hz之间，因此，该敲击力涵盖了层合板5的自然频率，从而如上所述，可以灵敏地检测出在层合板5中的局部损伤。

结构损伤检测设备310中的传感器2为加速度传感器，其感测由待检测结构表面传递到敲击装置3的加速度信号，并且将其发送到计算机11中的信号处理部件220进行进一步处理。信号处理部件220的结构和参考图1所示的实施方式基本相同，这里就不再进行赘述。

如图3和4所示，在对层合板5进行检测时，将结构损伤检测设备310置于层合板5的中心处，并以大约V＝1.3cm/s的速度沿层合板5的长度方向(即轴向)进行扫描。在扫描的同时，用上述频段在20-100Hz之间的敲击力来敲击，并同时由加速度传感器2采集加速度信号以发送到信号处理部件220进行进一步处理。

在该实施例中，以1.3cm为单位来沿板长度方向划分层合板5，由于小车1的行进速度为1.3cm/s，而且加速度传感器采集信号的采样率为1024次/秒。因此，谱图获取装置222按照1024个加速度信号的数据长度获取感测信号，并采用短时傅里叶变换方法获得频谱图。谱图包络线截取装置224截取与敲击力的频段相对应的、即20-100Hz频段的频谱包络线。损伤指示值计算装置126选择大小为20-100Hz的频率作为频谱向量中的分量相对应的频率，即损伤指示值计算装置226将频谱包络线向量化为长度为101的向量。然后，损伤指示值计算装置126计算MAC系数矩阵：

$\mathrm{MAC}(i,j)=\frac{Y_i\·Y_j}{\left|Y_i\right|\×\left|Y_j\right|}.$

随后，损伤位置确定装置228可以根据MAC系数矩阵的值来确定损伤位置。

图5示出了MAC矩阵元素的分布图，这是对尺寸为120×15×3cm的层合板中存在人为制造的一条约20×8mm的锯缝6的层合板进行检测所获取的结果。从图5中可以明显看出，中心钢棒和锯缝位置的MAC系数明显低于其它位置的MAC系数，由此可以说明本发明方法能够有效地检测出结构刚度发生突然变化的位置，并可由此定位损伤。

图6示意性地示出了根据本发明一个实施例的结构损伤检测方法的流程图。如图6所示，该方法始于步骤S610，在该步骤S 610中，使用根据本发明的结构损伤检测设备沿着待检测结构的表面进行扫描，其中在进行扫描的同时，结构损伤检测设备中的敲击装置以具有预定频率的敲击力进行敲击。如上理论分析所述，敲击装置所施加的敲击力的频率应当接近待检测结构的一个或者多个自然频率之一，以便高灵敏度地检测出结构损伤。根据本发明的一个实施例，离散频率接近待检测结构的某个自然频率意指将离散频率设置为待检测结构的该自然频率的±20％范围之内的某个值，更可取的是±10％范围之内的某个值，以及±5％范围之内的某个值。可选地，敲击力可以具有多个离散频率，而多个离散频率中的至少一个频率接近待检测结构的多个自然频率中的至少一个。进一步可选的是，敲击力可以具有一个或者多个频段，而所述一个或者多个频段中的至少一个涵盖了接近待检测结构的多个自然频率的至少一个的值。同时，为了防止出现待检测结构的共振，可以不将敲击力的频率设置为刚好等于待检测结构的自然频率，或者敲击力的频率设置为在某个范围内快速变化，从而使得待检测结构发生共振的时间很短而不会导致待检测结构的损伤。

随后，在步骤S620中，利用根据本发明的结构损伤检测设备中的传感器来感测该敲击装置的响应信号。如上参考结构损伤检测设备所述，所感测的信号可以是能够反映由待检测结构表面传递到敲击装置的在敲击位置处的任何响应信号。例如传感器感测的信号可以为位移、速度或加速度。特别的，传感器所感测的信号为加速度信号更便于工程实施。在步骤S630中，对所述传感器感测到的信号进行变换处理以获取在待检测结构表面的每个位置处的信号谱。根据检测精度的要求，需要将待检测结构的表面分为多个部分。在每个部分上感测到的信号是结构损伤检测设备扫描该部分所需时间上的分布，因此可以对信号在时间域上分布进行变换处理以获取该信号在频率域或尺度域上的表示。这种变换处理可以利用现有技术领域的任何变换处理进行，例如该变换处理可以是短时傅立叶变换、小波变换或Hi lbe八-Huang变换等。

随后，在步骤S640中，从信号谱图中截取与敲击装置的敲击力的预定频率相对应的部分，并获取所截取谱图的包络线。由于敲击力的至少一个频率接近待检测结构的至少一个自然频率，因此，所截取的谱图包络线应当涵盖待检测结构的至少一个自然频率。在敲击力具有某个频段的情况下，所截取的谱图包络线具有与该敲击力的频段相对应的频率范围。上面所给出的、有关所截取的谱图包络线的频率范围的示例仅仅为示范性的。本领域技术人员应当理解，只要所截取的谱图包络线的频段与敲击装置的敲击频率相对应，并且适于在后续数据处理中使用，则所有的频率范围都在本发明的保护范围之内。。

在步骤S650中，根据所截取的谱图包络线来计算待检测结构表面的每个位置处的损伤指示值，该损伤指示值反映了该位置处的谱图包络线与其他位置处的谱图包络线的相似程度。可以有多种方式来计算损伤指示值。

根据本发明的一个实施例，可以如下所述来计算损伤指示值：

1、将谱图包络线转换为谱图向量，这可以通过获取谱图包络线中的多个频率处的幅值作为谱图向量的分量来进行转换。谱图向量的分量数量可以根据检测精度以及系统的处理性能等确定。分量对应的频率可以均匀分布在谱图包络线的频段中，也可以不均匀分布。但是，对于待检测结构而言，每个位置处的频率选择方式应当是相同的；

2、在生成了每个位置处的谱图向量之后，计算反映当前位置处的谱图向量和其他位置处的谱图向量的相似程度的损伤指示值。可以有多种方式来计算两个谱图向量之间的相似程度。根据本发明的一个实施例，可以将损伤指示值计算为：

$\left(\underset{j=1...n}{\mathrm{\Σ}}\frac{Y_i\·Y_j}{\left|Y_i\right|\×\left|Y_j\right|}\right)/n,$

其中，Y_i和Y_j分别表示在所述结构上的当前位置和其他位置处的谱图向量，·表示谱图向量的内积运算，|Y_i|和|Y_j|表示向量Y_i和Y_j的长度，n为待检测结构所划分的单位检测部分的数量。

在Y_i和Y_j越相似，即包络线i和包络线j的形状越相似时，就越接近1。因此，针对某个谱图向量Y_i的损伤指示值可以是该谱图向量Y_i和所有其他谱图向量Y_j相似度的总和的平均值。

根据本发明的另一个实施例，可以利用MAC系数来获取损伤指示值。MAC系数矩阵的计算公式为：

$\mathrm{MAC}(i,j)=\frac{Y_i\·Y_j}{\left|Y_i\right|\×\left|Y_j\right|},$

其中，Y_i和Y_j分别表示在所述结构上第i个和第j个部分的谱图向量，·表示谱图向量的内积运算，|Y_i|和|Y_j|表示向量Y_i和Y_j的长度。MAC系数矩阵中的每个元素都表示了两个谱图向量之间的相似程度，其中主对角线上的元素一定等于1。MAC系数矩阵的第i行或第i列元素的大小就能反映出对应于某个谱图向量Y_i的损伤情况，因而也就是该谱图向量Y_i所对应位置的损伤指示值。

在为待检测结构表面的每个位置计算了损伤指示值之后，在步骤S660中，基于待检测结构表面的每个位置处的损伤指示值来确定待检测结构中的损伤位置。例如，可以将损伤指示值发生突然下降的位置确定为存在损伤的位置。比如，取MAC系数矩阵第一行元素的值为纵轴元素，并以所对应的结构位置为横轴画一条曲线，如果结构没有损伤，那么该曲线应该是比较光滑的；如果曲线在某处突然下降，就说明该处存在损伤，下降的幅度越大，损伤也就越严重。根据本发明的一个实施例，可以将该曲线的一阶导数值的绝对值超过预定阈值的位置确定为存在损伤的位置。

应当注意的是，在本发明的结构损伤检测设备的信号处理部件中，根据其要实现的功能而对其中的部件进行了逻辑划分，但是，本发明不受限于此，可以根据需要对信号处理部件中的各个部件进行重新划分或者组合，例如，可以将一些部件组合为单个部件，或者可以将一些部件进一步分解为更多的子部件。

本发明的信号处理部件实施例可以以硬件实现，或者以在一个或者多个处理器上运行的软件模块实现，或者以它们的组合实现。本领域的技术人员应当理解，可以在实践中使用微处理器或者数字信号处理器(DSP)来实现根据本发明实施例的信号处理部件中的一些或者全部部件的一些或者全部功能。本发明还可以实现为用于执行这里所描述的方法的一部分或者全部的设备或者装置程序(例如，计算机程序和计算机程序产品)。这样的实现本发明的程序可以存储在计算机可读介质上，或者可以具有一个或者多个信号的形式。这样的信号可以从因特网网站上下载得到，或者在载体信号上提供，或者以任何其他形式提供。

应该注意的是上述实施例对本发明进行说明而不是对本发明进行限制，并且本领域技术人员在不脱离所附权利要求的范围的情况下可设计出替换实施例。在权利要求中，不应将位于括号之间的任何参考符号构造成对权利要求的限制。单词“包含”不排除存在未列在权利要求中的元件或步骤。位于元件之前的单词“一”或“一个”不排除存在多个这样的元件。本发明可以借助于包括有若干不同元件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的单元权利要求中，这些装置中的若干个可以是通过同一个硬件项来具体体现。单词第一、第二、以及第三等的使用不表示任何顺序。可将这些单词解释为名称。

Claims (24)

1.一种利用结构损伤检测设备来检测结构的损伤的结构损伤检测方法，其中该结构损伤检测设备包括敲击装置和传感器，所述方法包括步骤：

使用所述结构损伤检测设备沿着所述结构的表面进行扫描，其中在进行扫描的同时，所述结构损伤检测设备中的敲击装置以包括预定频率的敲击力进行敲击，该预定频率接近所述结构的一个或者多个自然频率之一；

利用所述结构损伤检测设备中的传感器感测在所述结构表面的每个位置处、由所述结构表面传递到所述敲击装置的响应信号，其中所述传感器感测的信号可以为位移、速度和加速度中的至少之一；

对所述传感器感测到的信号进行变换处理以获取在所述结构表面的每个位置处的信号谱图；

从所述信号谱图中截取与所述敲击力的预定频率相对应的谱图包络线；

计算所述结构表面的每个位置处的损伤指示值，该损伤指示值反映了该位置处的谱图包络线与其他位置处的谱图包络线的相似程度；以及

将损伤指示值发生突变的位置确定为存在结构损伤的位置。

2.如权利要求1所述的结构损伤检测方法，其中所述预定频率的值为所述结构的所述一个或者多个自然频率之一的±20％范围之内的某个值。

3.如权利要求2所述的结构损伤检测方法，其中所述预定频率的值为所述结构的所述一个或者多个自然频率之一的±10％范围之内的某个值。

4.如权利要求3所述的结构损伤检测方法，其中所述预定频率的值为所述结构的所述一个或者多个自然频率之一的±5％范围之内的某个值。

5.如权利要求1-4中任一个所述的结构损伤检测方法，其中所述敲击力具有多个离散频率，且多个离散频率中的至少一个接近所述结构的至少一个自然频率。

6.如权利要求1-4中任一个所述的结构损伤检测方法，其中所述敲击力具有一个或者多个频段，且所述一个或者多个频段中的至少一个涵盖了接近所述结构的至少一个自然频率的值。

7.如权利要求1-4中的任一项所述的方法，其中还包括步骤：根据所述损伤指示值的突变程度来确定损伤大小。

8.如权利要求1-4中的任一个所述的结构损伤检测方法，其中计算所述结构表面的每个位置处的损伤指示值的步骤包括：

为所述结构表面中的所有位置，生成与所述位置处的谱图包络线相对应的谱图向量，其中所述谱图向量的每个分量为谱图包络线中某个频率所对应的谱值；以及

为所述结构表面中的每个位置，计算反映当前位置处的谱图向量和其他位置处的谱图向量的相似程度的损伤指示值。

9.如权利要求8所述的结构损伤检测方法，其中所述计算反映该位置处的谱图向量和其他位置处的谱图向量的相似程度的损伤指示值包括：

将损伤指示值计算为：

其中，Y_i和Y_j分别表示在所述结构上的当前位置和其他位置处的谱图向量，·表示谱图向量的内积运算，|Y_i|和|Y_j|表示向量Y_i和Y_j的长度，以及n为所述结构中的检测位置数量。

10.如权利要求8所述的结构损伤检测方法，其中所述计算反映该位置处的谱图向量和其他位置处的谱图向量的相似程度的损伤指示值包括：

计算MAC系数矩阵，其被计算为：

$\mathrm{MAC}(i,j)=\frac{Y_i\·Y_j}{\left|Y_i\right|\×\left|Y_j\right|},$

其中，Y_i和Y_j分别表示在所述结构上第i个和第j个部分的谱图向量，·表示谱图向量的内积运算，|Y_i|和|Y_j|表示向量Y_i和Y_j的长度；以及

根据所述MAC系数矩阵中与该位置相对应的行或者列元素的相对大小来确定该位置的损伤指示值。

11.如权利要求1-4中的任一个所述的结构损伤检测方法，其中所述结构的形状为细长型，沿所述结构轴向的每个截面相似，所述敲击装置在所述结构的表面上沿轴向进行扫描，所述传感器感测所述结构沿轴向方向每个位置处的响应信号。

12.如权利要求1-4中的任一个所述的结构损伤检测方法，其中所述传感器为感测由所述结构表面传递到敲击装置的加速度响应的传感器。

13.一种结构损伤检测设备，包括：

敲击装置，用于在待检测结构表面上以具有预定频率的敲击力进行敲击，其中所述预定频率接近所述待检测结构的一个或者多个自然频率之一；以及

传感器，感测所述由待检测结构表面传递到敲击装置的响应信号，其中所述传感器感测的响应信号可以为位移、速度和加速度中的至少之一；

其中所述敲击装置和所述传感器集成在一起，以便所述敲击装置在所述待检测结构表面上进行敲击的同时，所述传感器感测由所述待检测结构表面传递到所述敲击装置的响应信号。

14.如权利要求13所述的结构损伤检测设备，其中所述传感器可以为加速度传感器。

15.如权利要求13所述的结构损伤检测设备，其中所述预定频率的值为所述待检测结构的一个或者多个自然频率之一的±20％范围之内的某个值。

16.如权利要求15所述的结构损伤检测设备，其中所述预定频率的值为所述待检测结构的一个或者多个自然频率之一的±10％范围之内的某个值。

17.如权利要求16所述的结构损伤检测设备，其中所述预定频率的值为所述待检测结构的一个或者多个自然频率之一的±5％范围之内的某个值。

18.如权利要求13-17中的任一个所述的结构损伤检测设备，其中所述敲击力具有多个离散频率，且多个离散频率中的至少一个接近所述待检测结构的至少一个自然频率。

19.如权利要求13-17中的任一个所述的结构损伤检测方法，其中所述敲击力具有一个或者多个频段，且所述一个或者多个频段中的至少一个涵盖了接近所述待检测结构的至少一个自然频率的值。

20.一种结构损伤检测系统，包括：

如权利要求13-19中的任一项所述的结构损伤检测设备；以及

信号处理部件，用于对所述结构损伤检测设备中的传感器所感测到的信号进行处理以确定待检测结构中的损伤位置，所述信号处理部件包括：

谱图获取装置，用于对所述传感器感测到的信号进行变换处理以获取在所述结构表面的每个位置处的信号谱图；

谱图包络线截取装置，用于从所述信号谱图中截取与所述敲击力的预定频率相对应的谱图包络线；

损伤指示值计算装置，用于计算所述结构表面的每个位置处的损伤指示值，该损伤指示值反映了该位置处的谱图包络线与其他位置处的谱图包络线的相似程度；以及

损伤位置确定装置，用于基于所述结构表面的每个位置处的损伤指示值来确定所述结构中的损伤位置。

21.如权利要求20所述的结构损伤检测系统，其中所述损伤指示值计算装置为所述结构表面中的所有位置，生成与所述位置处的谱图包络线相对应的谱图向量，其中所述谱图向量的每个分量为谱图包络线中某个频率所对应的谱值；以及

所述损伤指示值计算装置为所述结构表面中的每个位置，计算反映当前位置处的谱图向量和其他位置处的谱图向量的相似程度的损伤指示值。

22.如权利要求21所述的结构损伤检测系统，其中所述损伤指示值计算装置将损伤指示值计算为：

$\left(\underset{j=1...n}{\mathrm{\Σ}}\frac{Y_i\·Y_j}{\left|Y_i\right|\×\left|Y_j\right|}\right)/n,$

其中，Y_i和Y_j分别表示在所述结构上的当前位置和其他位置处的谱图向量，·表示谱图向量的内积运算，|Y_i|和|Y_j|表示向量Y_i和Y_j的长度，以及n为所述结构中的检测位置数量。

23.如权利要求21所述的结构损伤检测系统，其中所述损伤指示值计算装置被配置为：

计算MAC系数矩阵，该MAC系数矩阵被计算为：

$\mathrm{MAC}(i,j)=\frac{Y_i\·Y_j}{\left|Y_i\right|\×\left|Y_j\right|},$

其中，Y_i和Y_j分别表示在所述结构上第i个和第j个部分的谱图向量，·表示谱图向量的内积运算，|Y_i|和|Y_j|表示向量Y_i和Y_j的长度；以及

根据所述MAC系数矩阵中与该位置相对应的行或者列元素的相对大小来确定该位置的损伤指示值。

24.如权利要求21-23中的任一个所述的结构损伤检测系统，其中所述待检测结构的形状为细长型，沿所述待检测结构轴向的每个截面相似，所述敲击装置在所述待检测结构的表面上沿轴向进行扫描，所述传感器感测由所述待检测结构表面传递到所述敲击装置在所述待检测结构沿轴向方向每个位置处的响应信号。

Images