Lamb Wave Phased Array Focus-Imaging Method Based on Frequency Response Function

Abstract

The disclosure discloses a Lamb wave phased array focus-imaging method based on a frequency response function. In the method, a piezoelectric sensor array is arranged on a surface of a tested structure, the frequency response function of an excitation and acquisition pair formed by an excitation array element and an acquisition array element is calculated according to a full-band response signal, and a dispersion pre-compensation signal is constructed; the dispersion pre-compensation signal and the frequency response function are multiplied in a frequency domain to obtain a frequency domain pre-compensation response signal; and according to a distance from the acquisition array element to a focal point at the coordinates, the dispersion of the frequency domain pre-compensation response signal is post-compensated, so as to obtain a frequency domain dispersion post-compensation signal until all sensor excitation and acquisition pairs are traversed.

Claims

1. A Lamb wave phased array focus-imaging method based on a frequency response function, the method comprising the following steps: S100, arranging a piezoelectric sensor array on a surface of a tested structure, sequentially exciting a full-band signal f.sub.m(t) by array elements P.sub.m (m=1, . . . N) of the piezoelectric sensor array, and collecting a response signal u.sub.mn(t) of the tested structure by other array elements P.sub.n (n=1, . . . N, n≠m), wherein the coordinates of the excitation array element P.sub.m are set as (x.sub.m, y.sub.m), the coordinates of the acquisition array element P.sub.n are set as (x.sub.n, y.sub.n), and subscripts m and n are the number identifiers of the excitation array element and the acquisition array element respectively; S200, according to the response signal u.sub.mn(t) calculating the frequency response function H.sub.mn(ω) of an excitation and acquisition pair composed of the excitation array element P.sub.m at coordinates (x.sub.m, y.sub.m) and the acquisition array element P.sub.n at coordinates (x.sub.n, y.sub.n), wherein ω is angular frequency and t is time, the excitation array element P.sub.m and the acquisition array element P.sub.n form an excitation and acquisition pair P.sub.m-P.sub.n, when one sensor in the sensor array acts as the excitation array element, there are (N−1) excitation and acquisition pairs in total, and when all sensors take turns to act as the excitation array element, there are N(N−1) excitation and acquisition pairs in total; S300, according to a distance L.sub.m from the excitation array element P.sub.m at coordinates (x.sub.m, y.sub.m) to a focal point F.sub.xy at coordinates (x, y) and a narrow-band signal S(ω) set to focus at t.sub.F, constructing a dispersion pre-compensation signal Pcom_S.sub.m(ω); S400, multiplying the dispersion pre-compensation signal Pcom_S.sub.m(ω) by the frequency response function H.sub.mn(ω) in the frequency domain to obtain a frequency domain pre-compensation response signal G.sub.mn(ω); S500, post-compensating for the dispersion of the frequency domain pre-compensation response signal G.sub.mn(ω) according to a distance L.sub.n from the acquisition array element P.sub.n at coordinates (x.sub.n, y.sub.n) to a focal point F.sub.xy at coordinates (x, y), so as to obtain a frequency domain dispersion post-compensation signal Comp_G.sub.mn(ω), and repeating S300 to S500 until all the sensor excitation and acquisition pairs P.sub.m-P.sub.n are traversed; S600, superposing time domain dispersion post-compensation signals comp_g.sub.mn(t) of all the sensor excitation and acquisition pairs P.sub.m-P.sub.n, taking the Hilbert envelope amplitude I.sub.x,y of a superposed signal at a focusing moment t.sub.F as an indicator of damage degree of the focal F.sub.xy at coordinates (x, y), and repeating S300 to S600 until all the focal points are detected; and S700, assigning all pixels with values for damage imaging, and taking the value of each pixel as a damage degree indication value I.sub.x,y corresponding to a focal position.

2. The method according to claim 1, wherein in S100, the excitation full-band signal f.sub.m(t) is a step signal and a chirp signal.

3. The method according to claim 1, wherein in S200, the frequency response function H.sub.mn(ω) is obtained according to the following formula: $H_{\mathrm{mn}}\left(\omega \right)=\frac{U_{\mathrm{mn}}\left(\omega \right)}{F_m\left(\omega \right)}$ wherein ω is an angular frequency, and U.sub.mn(ω) and F.sub.m(ω) are the frequency domain expression of the response signal u.sub.mn(t) and the excitation signal f.sub.m(t) respectively.

4. The method according to claim 1, wherein in S300, the surface of the tested structure is divided into imaging rectangular cells, each cell represents a pixel, and a center of each cell serves as a focal point of phased array focusing, with coordinates (x, y).

5. The method according to claim 1, wherein in S300, based on the focused narrow-band signal S(ω), the dispersion pre-compensation signal Pcom_S.sub.m(ω) is:
Pcom_S.sub.m(ω)=S(ω).Math.e.sup.ik(ω)L.sup.m wherein e is a natural logarithm, i is an imaginary unit, the wave number k(ω) is a function of angular frequency ω, the wave number k(ω) changes with an angular frequency ω, and the unit is 1/m.

6. The method according to claim 4, wherein the wave number k(ω) is obtained by numerical solution or experiment.

7. The method according to claim 4, wherein the focused narrow-band signal S(ω) is:
S(ω)=∫.sub.−∞.sup.+∞(H(t−1/f.sub.m)−H(t−1/f.sub.m−t.sub.F))(1−cos(2πf.sub.mt))sin(2πf.sub.ct)e.sup.−iωtdt wherein t is time, dt is the first-order differential of time, H(t) is a unit step function, e is a natural logarithm, i is an imaginary unit, f.sub.m and f.sub.c are a modulation frequency and a center frequency respectively, and t.sub.F is a set focusing moment and is set according to the size of a detection area.

8. The method according to claim 1, wherein in S500, the dispersion post-compensation signals Comp_G.sub.mn(ω) are calculated by the following formula:
Comp_G.sub.mn(ω)=G.sub.mn(ω).Math.e.sup.ik(ω)L.sup.n wherein e is a natural logarithm, i is an imaginary unit, the wave number k(ω) is a function of angular frequency ω, the wave number k(ω) changes with an angular frequency ω, and the unit is 1/m.

9. The method according to claim 1, wherein in S600, the time domain dispersion post-compensation signals comp_g.sub.mn(t) are calculated by the following formula:
comp_g.sub.mn(t)=∫.sub.−∞.sup.+∞Comp_G.sub.mn(ω)e.sup.iωtdω wherein e is a natural logarithm, i is an imaginary unit, ω is an angular frequency, and dω represents the first-order differential of the angular frequency.

10. The method according to claim 1, wherein in S600, the damage degree indication value I.sub.x,y of a focal point at coordinates (x, y) is calculated by the following formula: $I_{x,y}=\mathrm{Hilbert}\left(\underset{m=1}{\overset{N}{.Math.}}\underset{n=1,n\ne m}{\overset{N}{.Math.}}{\mathrm{comp\_g}}_{\mathrm{mn}}\left(t\right)\right){|}_{t=t_F}$ wherein Hilbert(⋅) is Hilbert transform, m represents an excitation signal of an m<th> array element, n represents an acquisition signal of an n<th> array element, and t.sub.F is a focusing moment.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0032] Various other advantages and benefits of the disclosure will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings of the specification are only for the purpose of illustrating preferred embodiments, and are not to be considered as limiting the disclosure. Obviously, the drawings in the following description are only some embodiments. For those skilled in the art, other drawings can be obtained according to these drawings without creative labor. Furthermore, throughout the drawings, the same parts are denoted by the same reference numerals.

[0033] In the drawings:

[0034] FIG. 1 is a diagram of the steps of a Lamb wave phased array focus-imaging method based on a frequency response function according to an embodiment of the disclosure;

[0035] FIG. 2 is a diagram of one-dimensional linear sensor array arrangement, focus division, signal excitation and acquisition in an aluminum plate structure in a Lamb wave phased array focus-imaging method based on a frequency response function according to an embodiment of the disclosure;

[0036] FIGS. 3(a) to 3(b) are waveform diagrams of a pre-compensation excitation signal designed according to a narrow-band focusing signal in a Lamb wave phased array focus-imaging method based on a frequency response function according to an embodiment of the disclosure;

[0037] FIG. 4 is a waveform diagram of post-compensating for the dispersion of a dispersion signal in a Lamb wave phased array focus-imaging method based on a frequency response function according to an embodiment of the disclosure;

[0038] FIG. 5 is a diagram of superimposing dispersion post-compensation signals of all sensor pairs in FIG. 4 in a Lamb wave phased array focus-imaging method based on a frequency response function according to an embodiment of the disclosure;

[0039] FIG. 6 is a waveform diagram of a superimposed time domain dispersion post-compensation signal of an undamaged focal point (i.e., focal point D2 shown in FIG. 2) in a Lamb wave phased array focus-imaging method based on a frequency response function according to an embodiment of the disclosure;

[0040] FIG. 7 is a waveform diagram of a superimposed time domain dispersion post-compensation signal of a damaged focal point (i.e., focal point D1 shown in FIG. 2) in a Lamb wave phased array focus-imaging method based on a frequency response function according to an embodiment of the disclosure;

[0041] FIG. 8 is a non-reference damage detection imaging diagram in FIG. 2 of a Lamb wave phased array focus-imaging method based on a frequency response function according to an embodiment of the disclosure;

[0042] FIG. 9 is a reference damage detection imaging diagram in FIG. 2 of a Lamb wave phased array focus-imaging method based on a frequency response function according to an embodiment of the disclosure;

[0043] FIG. 10 is a diagram of two-dimensional annular sensor array arrangement and in-plate damage in an aluminum plate structure in a Lamb wave phased array focus-imaging method based on a frequency response function according to an embodiment of the disclosure;

[0044] FIG. 11 is a non-reference D1 single damage detection imaging diagram in FIG. 10 of a Lamb wave phased array focus-imaging method based on a frequency response function according to an embodiment of the disclosure; and

[0045] FIG. 12 is a non-reference D1 and D2 multi-damage detection imaging diagram in FIG. 10 of a Lamb wave phased array focus-imaging method based on a frequency response function according to an embodiment of the disclosure.

[0046] The disclosure will be further explained with reference to the drawings and embodiments.

DETAILED DESCRIPTION

[0047] Specific embodiments of the disclosure will be described in more detail with reference to FIGS. 1 to 12. Although specific embodiments of the disclosure are shown in the drawings, it should be understood that the disclosure may be embodied in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided in order to understand the disclosure more thoroughly and to fully convey the scope of the disclosure to those skilled in the art.

[0048] It should be noted that certain words are used in the specification and claims to refer to specific components. It should be understood by those skilled in the art that the same component may be referred to by different nouns. In this specification and claims, components are distinguished not according to nouns, but according to functions. The term “include” or “comprise” mentioned in the whole specification and claims is an open term, and should be interpreted as “including but not limited to”. The rest of the specification describes preferred embodiments for implementing the disclosure, but the description is for the purpose of the general principles of the specification, and is not intended to limit the scope of the disclosure. The scope of protection of the disclosure shall be as defined by the appended claims.

[0049] In order to facilitate the understanding of the embodiments of the disclosure, specific embodiments will be further explained with reference to the accompanying drawings, and the drawings do not constitute a limitation of the embodiments of the disclosure.

[0050] FIG. 1 is a diagram of the steps of a Lamb wave phased array focus-imaging method based on a frequency response function. As shown in FIG. 1, a Lamb wave phased array focus-imaging method based on a frequency response function includes the following steps:

[0051] S100, a piezoelectric sensor array is arranged on a surface of a tested structure, a full-band signal f.sub.m(t) is sequentially excited by array elements P.sub.m (m=1, . . . N) of the piezoelectric sensor array, and a response signal u.sub.mn(t) of the tested structure is collected by other array elements P.sub.n (n=1, . . . N, n≠m), wherein the coordinates of the excitation array element P.sub.m are set as (x.sub.m, y.sub.m), the coordinates of the acquisition array element P.sub.n are set as (x.sub.n, y.sub.n), and subscripts m and n are the number identifiers of the excitation array element and the acquisition array element respectively;

[0052] S200, according to the response signal u.sub.mn(t), the frequency response function H.sub.mn(ω) of an excitation and acquisition pair composed of the excitation array element P.sub.m at coordinates (x.sub.m, y.sub.m) and the acquisition array element P.sub.n at coordinates (x.sub.n, y.sub.n) is calculated, wherein ω is angular frequency and t is time, the excitation array element P.sub.m and the acquisition array element P.sub.n form an excitation and acquisition pair P.sub.m-P.sub.n, when one sensor in the sensor array acts as the excitation array element, there are (N−1) excitation and acquisition pairs in total, and when all sensors take turns to act as the excitation array element, there are N(N−1) excitation and acquisition pairs in total;

[0053] S300, according to a distance L.sub.m from the excitation array element P.sub.m at coordinates (x.sub.m, y.sub.m) to a focal point F.sub.xy at coordinates (x, y) and a narrow-band signal S(ω) set to focus at t.sub.F, a dispersion pre-compensation signal Pcom_S.sub.m(ω) is constructed;

[0054] S400, the dispersion pre-compensation signal Pcom_S.sub.m(ω) is multiplied by the frequency response function H.sub.mn(ω) in the frequency domain to obtain a frequency domain pre-compensation response signal G.sub.mn(ω);

[0055] S500, the dispersion of the frequency domain pre-compensation response signal G.sub.mn(ω) is post-compensated according to a distance L.sub.n from the acquisition array element P.sub.n at coordinates (x.sub.n, y.sub.n) to a focal point F.sub.xy at coordinates (x, y), so as to obtain a frequency domain dispersion post-compensation signal Comp_G.sub.mn(ω), and S300 to S500 are repeated until all the sensor excitation and acquisition pairs P.sub.m-P.sub.n are traversed;

[0056] S600, time domain dispersion post-compensation signals comp_g.sub.mn(t) of all the sensor excitation and acquisition pairs P.sub.m-P.sub.n are superposed, the Hilbert envelope amplitude I.sub.x,y of a superposed signal at the focusing moment t.sub.F is taken as an indicator of damage degree of the focal F.sub.xy at coordinates (x, y), and S300 to S600 are repeated until all the focal points are detected; and

[0057] S700, all pixels are assigned with values for damage imaging, and the value of each pixel is taken as a damage degree indication value I.sub.x,y corresponding to a focal position.

[0058] In a preferred embodiment of the method, in S100, the excitation signal f.sub.m(t) is a full-band signal, that is, the amplitude of each frequency point in the frequency domain expression F.sub.m(ω) of f.sub.m(t) is not 0. In order to obtain the frequency response function H.sub.mn(ω), the excitation signal f.sub.m(t) can be a step signal and a chirp signal, but is not limited to the two signals.

[0059] In a preferred embodiment of the method, in S200, the frequency response function H.sub.mn(ω) is obtained according to the following formula:

[00003]$H_{\mathrm{mn}}\left(\omega \right)=\frac{U_{\mathrm{mn}}\left(\omega \right)}{F_m\left(\omega \right)}$

[0060] wherein ω is an angular frequency, t is time, and U.sub.mn(ω) and F.sub.m(ω) are the frequency domain expression of the response signal u.sub.mn(t) and the excitation signal f.sub.m(t) respectively.

[0061] In a preferred embodiment of the method, in S300, the surface of the tested structure is divided into imaging rectangular cells, each cell represents a pixel, and a center of each cell serves as a focal point F.sub.xy of phased array focusing, with coordinates (x, y).

[0062] In a preferred embodiment of the method, in S300, based on the focused narrow-band signal S(ω), the dispersion pre-compensation signal Pcom_S.sub.m(ω) is:

Pcom_S.sub.m(ω)=S(ω).Math.e.sup.ik(ω)L.sup.m,

[0063] wherein e is a natural logarithm, i is an imaginary unit, the wave number k(ω) is a function of angular frequency ω, the wave number k(ω) changes with the angular frequency ω, and the unit is 1/m.

[0064] In a preferred embodiment of the method, the wave number k(ω) is obtained by numerical solution or experiment.

[0065] In a preferred embodiment of the method, the focused narrow-band signal S(ω) is:

S(ω)=∫.sub.−∞.sup.+∞(H(t−1/f.sub.m)−H(t−1/f.sub.m−t.sub.F))(1−cos(2πf.sub.mt))sin(2πf.sub.ct)e.sup.−iωtdt

[0066] wherein t is time, dt is the first-order differential of time, H(t) is a unit step function, e is a natural logarithm, i is an imaginary unit, and f.sub.m and f.sub.c are a modulation frequency and a center frequency respectively, which are determined according to the characteristics of the damage to be detected. In the case of single damage, the wavelength corresponding to the selected center frequency should be in the same order of magnitude as the size of the detected damage, and the modulation frequency is obtained by reducing the center frequency by 1-5 times. When there are adjacent multiple damages, the width of a wave packet corresponding to the modulation frequency should be smaller than a distance between adjacent damages, and the center frequency should be determined according to the damage size. t.sub.F is a set focusing moment and is set according to the size of a detection area.

[0067] In a preferred embodiment of the method, in S500, the dispersion post-compensation signals Comp_G.sub.mn(ω) are calculated by the following formula:

Comp_G.sub.mn(ω)=G.sub.mn(ω).Math.e.sup.ik(ω)L.sup.n

[0068] wherein e is a natural logarithm, i is an imaginary unit, the wave number k(ω) is a function of angular frequency ω, the wave number k(ω) changes with an angular frequency co, and the unit is 1/m.

[0069] In a preferred embodiment of the method, in S600, the time domain dispersion post-compensation signals comp_g.sub.mn(t) are calculated by the following formula:

comp_g.sub.mn(t)=∫.sub.−∞.sup.+∞Comp_G.sub.mn(ω)e.sup.iωtdω

[0070] wherein e is a natural logarithm, i is an imaginary unit, ω is an angular frequency, and dω represents the first-order differential of the angular frequency.

[0071] In a preferred embodiment of the method, in S600, the damage degree indication value I.sub.x,y of a focal point F.sub.xy at coordinates (x, y) is calculated by the following formula:

[00004]$I_{x,y}=\mathrm{Hilbert}\left(\underset{m=1}{\overset{N}{.Math.}}\underset{n=1,n\ne m}{\overset{N}{.Math.}}{\mathrm{comp\_g}}_{\mathrm{mn}}\left(t\right)\right){|}_{t=t_F}$

[0072] wherein Hilbert(⋅) is Hilbert transform, m represents an excitation signal of an m<th> array element, n represents an acquisition signal of an n<th> array element, and t.sub.F is a focusing moment.

[0073] In an embodiment, S300 further includes the following steps:

[0074] S301, a narrow-band signal S(ω) focused at a set moment t.sub.F is designed:

S(ω)=∫(H(t−1/f.sub.m)−H(t−1/f.sub.m−t.sub.F))(1−cos(2πf.sub.mt))sin(2πf.sub.ct)e.sup.−iωtdt

[0075] wherein t is time, dt is the first-order differential of time, H(t) is a unit step function, e is a natural logarithm, i is an imaginary unit, f.sub.m and f.sub.c are a modulation frequency and a center frequency respectively which are flexibly determined according to the characteristics of the damage to be detected, and t.sub.F is a set focusing moment and is set according to the size of a detection area.

[0076] S302, a wave number k(ω) of the tested structure is obtained through numerical calculation or experiment; and

[0077] S303, after obtaining the narrow-band signal S(ω) and the wave number k(ω), a dispersion pre-compensation response Pcom_S.sub.m(ω) is calculated by the following formula:

Pcom_S.sub.m(ω)=S(ω).Math.e.sup.ik(ω)L.sup.m

[0078] wherein e is a natural logarithm, i is an imaginary unit, the wave number k(ω) is a function of angular frequency ω, the wave number k(ω) changes with an angular frequency co, and the unit is 1/m.

[0079] A dispersion post-compensation response Comp_G.sub.mn(ω) in S500 is obtained by multiplying the dispersion pre-compensation response by a corresponding dispersion wave number term in the frequency domain to eliminate the dispersion of Lamb waves at a specific propagation distance:

Comp_G.sub.mn(ω)=G.sub.mn(ω).Math.e.sup.ik(ω)L.sup.n

[0080] wherein e is a natural logarithm, k(ω) is a wave number of the tested structure related to an angular frequency ω, and the unit is 1/m.

[0081] In another embodiment, the time domain dispersion post-compensation signal comp_g.sub.mn(t) in S600 is obtained by inverse Fourier transform, so as to use time domain energy to characterize the energy in a distance domain:

comp_g.sub.mn(t)=∫.sub.−∞.sup.+∞Comp_G.sub.mn(ω)e.sup.iωtdω

[0082] wherein e is a natural logarithm, i is an imaginary unit, ω is an angular frequency, and dω represents the first-order differential of the angular frequency.

[0083] In another embodiment, in S600, the damage degree indicator I.sub.x,y of a focal point F.sub.xy at coordinates (x, y) is calculated by the following formula:

[00005]$I_{x,y}=\mathrm{Hilbert}\left(\underset{m=1}{\overset{N}{.Math.}}\underset{n=1,n\ne m}{\overset{N}{.Math.}}{\mathrm{comp\_g}}_{\mathrm{mn}}\left(t\right)\right){|}_{t=t_F}$

[0084] wherein Hilbert(⋅) is Hilbert transform, m represents an excitation signal of an m<th> array element, n represents an acquisition signal of an n<th> array element, and t.sub.F is a focusing moment.

[0085] In this embodiment, for the same detection focal point F.sub.xy, it is necessary to constantly switch the sensor excitation and acquisition pairs, and because a hardware channel cannot excite and acquire signals at the same time, the excitation array element and the acquisition array element cannot be of the same label, that is, n≠m.

[0086] See FIG. 2, which is a diagram of one-dimensional linear sensor array arrangement and Lamb wave focusing of an aluminum plate structure in this embodiment. The sensor array is a linear phased array with 9 elements, the pixels are divided into rectangular frames with side lengths dx and dy, and the geometric center of the rectangular frame is the focal point. P.sub.1˜P.sub.9 sensors excite the full-band signal f.sub.mn(t) in turn, while other sensors acquire the response signal u.sub.mn(t). When P.sub.1 excites the signal f.sub.1n(t), n=2 . . . 9, the response acquired by P.sub.2˜P.sub.9 is u.sub.12(t), u.sub.13(t) . . . (t). 9×(9−1) groups of response signals are acquired when excitation is conducted by P.sub.1˜P.sub.9 in turn. For the two focal points D1 and D2 in FIG. 2, D1 is a focal point with damage, marked with a black triangle, and D2 is a focal point without damage.

[0087] FIGS. 3(a) to 3(b) are waveform diagrams of a pre-compensation excitation signal designed according to a narrow-band focusing signal in this embodiment, which are used to explain the relationship between the narrow-band focusing signal and the pre-compensation excitation signal. FIG. 3(a) shows a narrow-band focusing signal, and FIG. 3(b) shows a dispersion pre-compensation excitation signal. When the dispersion pre-compensation signal propagates to 500 μs, it will have the waveform of FIG. 3(a), that is, the focusing moment set in this embodiment is t.sub.F=500 μs.

[0088] FIG. 4 is a waveform diagram of post-compensating for the dispersion of a dispersion signal in this embodiment, which is used to explain the relationship between the dispersion signal and the dispersion post compensation signal. In the figure, the wide waveform is the dispersion response signal of the excitation and acquisition pair P.sub.m-P.sub.n, the excitation and acquisition pair P.sub.m-P.sub.n is noted as P.sub.mn, and the narrow waveform at t.sub.F=500 μs is the corresponding dispersion post-compensation signal. It can be seen that compared with the dispersion response signal, the dispersion post-compensation signal has a more compact waveform, so the dispersion is suppressed.

[0089] FIG. 5 is a diagram of superimposing the dispersion post-compensation signals of all the excitation and acquisition pairs in FIG. 4 in this embodiment. By superimposing the dispersion post-compensation signals of all the excitation and acquisition pairs, the amplitude will be more intense, and the Hilbert envelope amplitude at t.sub.F is the damage degree indication value I.sub.x,y of the detection focal point. This value I.sub.x,y is the pixel value of a pixel at the focal point during damage imaging

[0090] FIG. 6 is a waveform diagram of a superimposed time domain dispersion post-compensation signal of an undamaged focal point (i.e., focal point D2 shown in FIG. 2) in this embodiment. Since there is no damage at D2, the reflected signal will be weak, so the value of the superimposed dispersion post-compensation signal at the focusing moment t.sub.F is small, and the damage degree indication value of this point I.sub.x,y=0.088, meaning that the pixel value of the pixel here is 0.088.

[0091] FIG. 7 is a waveform diagram of a superimposed time domain dispersion post-compensation signal of a damaged focal point (i.e., focal point D1 shown in FIG. 2) in this embodiment. Because there is damage at D1, the reflected signal will be intense. The damage degree indication value of D1 is I.sub.x,y=3.076, meaning that the pixel value of the pixel here is 3.076, which is about 35 times of the pixel value at D2, so the damaged point and undamaged point can be easily distinguished.

[0092] FIG. 8 is a non-reference damage detection imaging diagram of the aluminum plate illustrated in FIG. 2, with pixel side lengths dx=2 mm and dy=2 mm. Non-reference means that there is no need to subtract the step response signal in the damaged state from the step response signal in the undamaged state, and the imaging result obtained in this way is shown in FIG. 8. Because there is no reference signal, original signals will include complex and intense edge reflection signals. Therefore, in order to avoid interference, the detection area will be reduced, and the reduced areas include the part inward from the boundary of the tested structure and the part outward from the center of the sensor. These areas are dead zones of damage detection, which are inevitable for traditional Lamb wave phased array focusing, so the final damage imaging result is a circle.

[0093] FIG. 9 is a reference damage detection imaging diagram of the aluminum plate illustrated in FIG. 2, with pixel side lengths dx=2 mm and dy=2 mm Because the method of the disclosure can quickly extract the dispersion post-compensation response signal of each focal point, damage detection with reference signals can be quickly implemented, the size of the dead zone can be reduced by subtracting the reference signals, and the influence of multimodality of Lamb waves can be weakened by subtracting operation. The imaging diagram of reference damage detection in this embodiment is shown in FIG. 9. Compared with FIG. 8, it can be easily seen that the external and internal dead zones of the tested structure are reduced, which not only increases the detection area, but also improves the ability to reflect the damage. It is worth noting that there is a symmetrical pseudo-damage point in FIG. 8 and FIG. 9, which is caused by the symmetrical focus problem of the linear phased array, and can be eliminated by arranging a two-dimensional phased array, see the explanation of FIG. 10 to FIG. 12.

[0094] FIG. 10 is a diagram of two-dimensional sensor array arrangement and in-plate damage of an aluminum plate structure in this embodiment, in which there is damage at D1 and D2.

[0095] FIG. 11 is an imaging result of single damage detection when there is damage at D1 in FIG. 10 in this embodiment. The “mirror pseudo-damage” in FIGS. 8 and 9 no longer exists, and the damage can be accurately positioned.

[0096] FIG. 12 is an imaging result of multi-damage detection when there is damage at both D1 and D2 in FIG. 10 in this embodiment. The “mirror pseudo-damage” in FIGS. 8 and 9 no longer exists, and the damage can be accurately and clearly positioned. The disclosure is not limited to the annular phased array form, but also includes other forms of two-dimensional phased arrays, such as rectangular, spiral and concentric circular phased arrays, all of which can be used to carry out damage detection research based on the method of the disclosure.

[0097] From the above description with reference to the drawings and embodiments, it can be seen that the disclosure can accurately locate and identify damage, and also overcomes the disadvantage that dead zones are inevitable during phased array imaging.

[0098] Although the embodiments of the disclosure have been described above with reference to the drawings, the disclosure is not limited to the above-mentioned specific embodiments and application fields, and the above-mentioned specific embodiments are only illustrative and instructive, not restrictive. Under the inspiration of this specification and without departing from the scope protected by the claims of the disclosure, those of ordinary skill in the art can also obtain various forms, which all fall within the protection scope of the disclosure.