Network wavefield imaging methods for quantification of complex discontinuity in plate-like structures

Abstract

Network wavefield imaging methods are able to image significantly complex discontinuities or shapes in plate-like structures for superior ultrasonic structural health monitoring (SHM)/nondestructive evaluation (NDE). The imaging provides high-resolution location, shape and/or size images of a structure, and for discontinuities with more complicated profiles. Guided wave (GW) network wavefield imaging methods combine tomography and wavefield/wavenumber imaging algorithms. Metallic plate damage detection uses guided ultrasonic waves and non-contact laser vibrometry. Guided waves are generated by piezoelectric transducers (PZT). A non-contact scanning laser Doppler vibrometer (SLDV) measures the full velocity plate guided wave wavefields. Developed network wavefield imaging algorithms account for multiple-actuator excitations from different angles enclosing the discontinuity, with algorithms using intrinsic wave characteristics such as wavefield, wavenumber, or reconstructed wave energy. Determined locations, sizes and shapes of highlighted areas in wavefield, wavenumber and/or filter reconstructed energy-based images correlate with location, size and shape of damage in metallic plates.

Claims

1. Methodology for structural health monitoring (SHM)/nondestructive evaluation (NDE) testing of a plate-like structure to detect discontinuities therein, comprising: providing a plurality of guided ultrasonic wave actuators; arranging said plurality of actuators to enclose a designated scanning area of a target plate structure to be checked for a potential discontinuity therein; selectively exciting said actuators to produce guided waves in said structure to be tested; measuring a wavefield v(t,x) of guided waves in said structure generated upon excitation of the respective plurality of actuators; and creating a synthetic image based on such wavefield v(t,x), to detect a discontinuity in said plate structure within said scanning area, wherein said synthetic image comprises wavefield and/or wavenumber information of damage induced waves for damage detection and quantification, and wherein creating said synthetic image includes at least one of: A: determining the maximum absolute value of the waveform at each point x, by:
v.sub.mag(x)=max(abs(v(t,x))); exciting each actuator in turn so that a corresponding plurality of wavefield data v.sup.i(t,x) are obtained for each actuator; and creating data for said synthetic wavefield image through performing one of a summation, multiplication, and weighting factor synthesis of the wavefield data; or B: determining at least in part from said wavefield data a spatial wavenumber image at a selected frequency range as: $k^{*}\left(\stackrel{\_}{x}\right)=\frac{1}{N}\underset{i=1}{\overset{N}{.Math.}}.Math.k^{*}\left(\stackrel{\_}{x},f_i\right).Math.;$ where f.sub.i (i=1,2,3 ... N) is the selected frequency band centered at the desired frequency f.sub.0, which corresponds to the particular guided wave mode, for N number of actuators; and creating data for said synthetic wavenumber image through performing one of a summation, multiplication, and weighting factor synthesis of the wavenumber image; or C: determining at least in part from said wavefield data a frequency wavenumber representation as: $v\left(f,k\right)={\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }v\left(t,x\right)e^{-j\left(2\pi \mathrm{ft}-k.Math.x\right)}\mathrm{dtdx};$ conducting a filtering process on said frequency wavenumber representation to obtain a filtered spectrum as:
V.sub.F(f,k)=V(f,k)F(f,k); determining waves corresponding to the filtered spectrum Vf.sub.F(f,k) as: $v_F\left(t,x\right)=\frac{1}{{\left(2\pi \right)}^2}{\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }v\left(f,k\right)e^{j\left(2\pi \mathrm{ft}-k.Math.x\right)}\mathrm{dfdk};$ creating an energy map to show energy change due to a discontinuity within said designated scanning area as: $E_F\left(x\right)={\int }_0^t\frac{1}{2}{v}_F^2\left(\tau ,x\right)d\tau ;$ and creating data for said synthetic wavefield image through performing one of a summation, multiplication, and weighting factor synthesis of the energy map E.sub.F.sup.i(x).

2. Methodology as in claim 1, wherein said respective actuators are arranged at respectively different angles from a reference point relative to the scanning area.

3. Methodology as in claim 1, wherein creating said synthetic image includes providing high-resolution images of a detected discontinuity in a structure being tested, to indicate its location, shape, and/or size in said structure.

4. Methodology as in claim 1, wherein said guided waves are generated by piezoelectric transducers, and a non-contact scanning laser Doppler vibrometer (SLDV) is used to measure the wavefield v(t,x) of guided waves in said structure, based on the Doppler effect.

5. Methodology as in claim 1, wherein said plate-like structure comprises a metal plate for industrial uses requiring detection of any potential discontinuity therein.

6. Methodology as in claim 5, wherein said plate-like structure comprises a metal plate of steel, aluminum, or metal alloys.

7. Methodology as in claim 1, wherein said guided waves are generated by at least one of lead zirconate titanate (PZT) transducers, pulsed lasers, electromagnetic acoustic transducers (EMATs), and air-coupled transducer (ACT).

8. Methodology as in claim 7, wherein sensors comprising at least one of high spatial resolution sensors SLDV, motion controlled scanning sensors EMAT or ACT are used to measure the wavefield v(t,x) of guided waves in said structure.

9. Methodology as in claim 1, wherein said actuators comprise at least one of lead zirconate titanate (PZT) transducers, pulsed lasers, electromagnetic acoustic transducers (EMATs), and air-coupled transducer (ACT); and at least one of high spatial resolution sensors SLDV, motion controlled scanning sensors EMAT or ACT are used to measure the wavefield v(t,x) of guided waves in said structure.

10. Methodology as in claim 1, wherein creating said synthetic image includes: determining the maximum absolute value of the waveform at each point x, by:
v.sub.mag(x)=max(abs(v(t,x))); exciting each actuator in turn so that a corresponding plurality of wavefield data v.sup.i(t,x) are obtained for each actuator; and creating data for said synthetic wavefield image through performing one of a summation, multiplication, and weighting factor synthesis of the wavefield data.

11. Methodology as in claim 10, wherein said creating data step includes creating data for said synthetic wavefield image through performing the summation: $v_{\mathrm{sum}}\left(x\right)=\underset{i=1}{\overset{N}{.Math.}}{v}_{\mathrm{mag}}^i\left(x\right).$

12. Methodology as in claim 1, wherein creating said synthetic image includes: determining at least in part from said wavefield data a spatial wavenumber image at a selected frequency range as: $k^{*}\left(\stackrel{\_}{x}\right)=\frac{1}{N}\underset{i=1}{\overset{N}{.Math.}}.Math.k^{*}\left(\stackrel{\_}{x},f_i\right).Math.;$ where f.sub.i (i=1,2,3 ... N) is the selected frequency band centered at the desired frequency f.sub.0, which corresponds to the particular guided wave mode, for N number of actuators; and creating data for said synthetic wavenumber image through performing one of a summation, multiplication, and weighting factor synthesis of the wavenumber image.

13. Methodology as in claim 12, wherein said creating data step includes creating data for said synthetic wavenumber image through performing the summation of the wavenumber image as: $k_{\mathrm{sum}}^{*}\left(\stackrel{\_}{x}\right)=\underset{i=1}{\overset{N}{.Math.}}{k}_i^{*}\left(\stackrel{\_}{x}\right).$

14. Methodology as in claim 1, wherein creating said synthetic image includes: determining at least in part from said wavefield data a frequency wavenumber representation as: $v\left(f,k\right)={\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }v\left(t,x\right)e^{-j\left(2\pi ft-k.Math.x\right)}\mathrm{dt}\mathrm{dx};$ conducting a filtering process on said frequency wavenumber representation to obtain a filtered spectrum as:
V.sub.F(f,k)=V(f,k)F(f,k); determining waves corresponding to the filtered spectrum V.sub.F(f,k) as: $v_F\left(t,x\right)=\frac{1}{{\left(2\pi \right)}^2}{\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }v\left(f,k\right)e^{j\left(2\pi ft-k.Math.x\right)}\mathrm{dfdk};$ creating an energy map to show energy change due to a discontinuity within said designated scanning area as: $E_F\left(x\right)={\int }_0^t\frac{1}{2}{v}_F^2\left(\tau ,x\right)d\tau ;$ and creating data for said synthetic wavefield image through performing one of a summation, multiplication, and weighting factor synthesis of the energy map E.sub.F.sup.i(x).

15. Methodology as in claim 14, wherein said creating data step includes creating data for said synthetic wavefield image through performing the summation of the energy map E.sub.F.sup.i(x) as: $E_F^{\mathrm{sum}}\left(x\right)=\underset{i=1}{\overset{N}{.Math.}}{E}_F^i\left(x\right).$

16. Methodology as in claim 1, wherein said plurality of actuators are arranged in configurations which are at least one of circle, square, rectangle, and ellipse shapes.

17. Methodology for enhanced imaging for detection and quantification of complex discontinuities in metallic plate-like structures, comprising: providing a plurality of guided ultrasonic wave actuators comprising respective transducers; arranging said plurality of guided ultrasonic wave actuators to enclose a designated scanning area of a target plate structure to be checked for a potential damage discontinuity therein, with said actuators each arranged at respectively different angles from a reference point relative to the scanning area; selectively exciting in turn said guided ultrasonic wave actuators to produce guided waves in said structure to be tested; using a non-contact scanning laser Doppler vibrometer (SLDV) to measure a wavefield v(t,x) of guided waves in said structure generated upon excitation of the respective plurality of actuators; and using intrinsic wave characteristics of such measured wavefield to create a synthetic image comprising wavefield and/or wavenumber information of determined locations and sizes of high intensity areas in correlation with location, size, and shape of damage discontinuities in a target plate structure.

18. Methodology as in claim 17, wherein said plurality of actuators are arranged in configurations which are at least one of circle, square, rectangle, and ellipse shapes.

19. Methodology as in claim 17, wherein said transducers comprise at least one of lead zirconate titanate (PZT) transducers, pulsed lasers, electromagnetic acoustic transducers (EMATs), and air-coupled transducer (ACT).

20. Methodology as in claim 17, wherein said using intrinsic wave characteristics step includes using algorithms for analysis based on at least one of wavefield, wavenumber, and reconstructed wave energy characteristics.

21. Methodology as in claim 20, wherein creating said synthetic image includes: determining the maximum absolute value of the waveform at each point x, by:
v.sub.mag(x)=max(abs(v(t,x))); exciting each actuator in turn so that a corresponding plurality of wavefield data v.sup.i(t,x) are obtained for each actuator; and creating data for said synthetic wavefield image through performing one of a summation, multiplication, and weighting factor synthesis of the wavefield data.

22. Methodology as in claim 21, wherein said creating data step includes creating data for said synthetic wavefield image through performing the summation: $V_{\mathrm{sum}}\left(x\right)=\underset{i=1}{\overset{N}{.Math.}}{v}_{\mathrm{mag}}^i\left(x\right).$

23. Methodology as in claim 20, wherein creating said synthetic image includes: determining at least in part from said wavefield data a spatial wavenumber image at a selected frequency range as: $k^{*}\left(\stackrel{\_}{x}\right)=\frac{1}{N}\underset{i=1}{\overset{N}{.Math.}}.Math.k^{*}\left(\stackrel{\_}{x},f_i\right).Math.;$ where f.sub.i (i=1,2,3 ... N) is the selected frequency band centered at the desired frequency f.sub.0, which corresponds to the particular guided wave mode, for N number of actuators; and creating data for said synthetic wavenumber image through performing one of a summation, multiplication, and weighting factor synthesis of the wavenumber image.

24. Methodology as in claim 23, wherein said creating data step includes creating data for said synthetic wavenumber image through performing the summation of the wavenumber image as: $k_{\mathrm{sum}}^{*}\left(\stackrel{\_}{x}\right)=\underset{i=1}{\overset{N}{.Math.}}{k}_i^{*}\left(\stackrel{\_}{x}\right).$

25. Methodology as in claim 20, wherein creating said synthetic image includes: determining at least in part from said wavefield data a frequency wavenumber representation as: $v\left(f,k\right)={\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }v\left(t,x\right)e^{-j\left(2\pi ft-k.Math.x\right)}\mathrm{dtdx};$ conducting a filtering process on said frequency wavenumber representation to obtain a filtered spectrum as:
V.sub.F(f,k)=V(f,k)F(f,k); determining waves corresponding to the filtered spectrum V.sub.F(f,k) as: $v_F\left(t,x\right)=\frac{1}{{\left(2\pi \right)}^2}{\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }v\left(f,k\right)e^{j\left(2\pi ft-k.Math.x\right)}\mathrm{dfdk};$ creating an energy map to show energy change due to a discontinuity within said designated scanning area as: $E_F\left(x\right)={\int }_0^t\frac{1}{2}{v}_F^2\left(\tau ,x\right)d\tau ;$ and creating data for said synthetic wavefield image through performing one of a summation, multiplication, and weighting factor synthesis of the energy map E.sub.F.sup.i(x).

26. Methodology as in claim 25, wherein said creating data step includes creating data for said synthetic wavefield image through performing the summation of the energy map E.sub.F.sup.i(x) as: $E_F^{\mathrm{sum}}\left(x\right)=\underset{i=1}{\overset{N}{.Math.}}{E}_F^i\left(x\right).$

27. Apparatus for structural health monitoring (SHM)/nondestructive evaluation (NDE) testing of a plate-like structure to detect discontinuities therein, comprising: a plurality of guided ultrasonic wave actuators, arranged on a target plate structure to be checked for a potential discontinuity therein, so as to enclose a designated scanning area of such structure; means for selective operation of said guided ultrasonic wave actuators to excite in turn guided waves in said structure to be tested; means for measuring a wavefield v(t,x) of guided waves in said structure generated upon excitation of the respective plurality of actuators; and processing means for creating a synthetic image based on such wavefield v(t,x), to detect a discontinuity in said plate structure within said scanning area.

28. Apparatus as in claim 27, wherein said plurality of actuators are arranged in configurations which are at least one of circle, square, rectangle, ellipse shapes and/or any other shapes which can enclose the discontinuity.

29. Apparatus as in claim 27, wherein said respective actuators are arranged at respectively different angles from a reference point relative to the scanning area.

30. Apparatus as in claim 27, wherein said plate-like structure comprises a metal plate.

31. Apparatus as in claim 27, wherein: said actuators comprise piezoelectric transducers; and said means for measuring includes a non-contact scanning laser Doppler vibrometer (SLDV) for measuring the wavefield v(t,x) of guided waves in said structure, based on the Doppler effect.

32. Apparatus as in claim 27, wherein said actuators comprise at least one of lead zirconate titanate (PZT) transducers, pulsed lasers, electromagnetic acoustic transducers (EMATs), and air-coupled transducer (ACT).

33. Apparatus as in claim 27, wherein said synthetic image comprises wavefield and/or wavenumber information of damage induced waves for damage detection and quantification.

34. Apparatus as in claim 33, wherein: said function generator is further operative for exciting each actuator in turn so that a corresponding plurality of wavefield data v.sup.i(t,x) are obtained for each actuator; and said processing means includes a processor further operative for: determining the maximum absolute value of the waveform at each point x, by:
v.sub.mag(x)=max(abs(v(t,x))); and creating data for said synthetic wavefield image through performing one of a summation, multiplication, and weighting factor synthesis of the wavefield data.

35. Apparatus as in claim 34, wherein said processing means includes a processor further operative for creating data for said synthetic wavefield image through performing the summation: $v_{\mathrm{sum}}\left(x\right)=\underset{i=1}{\overset{N}{.Math.}}{v}_{\mathrm{mag}}^i\left(x\right).$

36. Apparatus as in claim 33, wherein said processing means includes a processor further operative for: determining at least in part from said wavefield data a spatial wavenumber image at a selected frequency range as: $k^{*}\stackrel{\_}{\left(x\right)}=\frac{1}{N}\underset{i=1}{\overset{N}{.Math.}}.Math.k^{*}\left(\stackrel{\_}{x},f_i\right).Math.;$ where f.sub.i (i=1,2,3 ... N) is the selected frequency band centered at the desired frequency f.sub.0, which corresponds to the particular guided wave mode, for N number of actuators; and creating data for said synthetic wavenumber image through performing one of a summation, multiplication, and weighting factor synthesis of the wavenumber image.

37. Apparatus as in claim 36, wherein said processing means includes a processor further operative for creating data for said synthetic wavenumber image through performing the summation of the wavenumber image as: $k_{\mathrm{sum}}^{*}\left(\stackrel{\_}{x}\right)=\underset{i=1}{\overset{N}{.Math.}}{k}_i^{*}\left(\stackrel{\_}{x}\right).$

38. Apparatus as in claim 33, wherein said processing means includes a processor further operative for: determining at least in part from said wavefield data a frequency wavenumber representation as: $v\left(f,k\right)={\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }v\left(t,x\right)e^{-j\left(2\pi \mathrm{ft}-k.Math.x\right)}\mathrm{dtdx};$ conducting a filtering process on said frequency wavenumber representation to obtain a filtered spectrum as:
V.sub.F(f,k)=V(f,k)F(f,k); determining waves corresponding to the filtered spectrum V.sub.F(f,k) as: $v_F\left(t,x\right)=\frac{1}{{\left(2\pi \right)}^2}{\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }v\left(f,k\right)e^{j\left(2\pi \mathrm{ft}-k.Math.x\right)}\mathrm{dfdk};$ creating an energy map to show energy change due to a discontinuity within said designated scanning area as: $E_F\left(x\right)={\int }_0^t\frac{1}{2}{v}_F^2\left(\tau ,x\right)d\tau ;$ and creating data for said synthetic wavefield image through performing one of a summation, multiplication, and weighting factor synthesis of the energy map E.sub.F.sup.i(x).

39. Apparatus as in claim 38, wherein said processing means includes a processor further operative for creating data for said synthetic wavefield image through performing the summation of the energy map E.sub.F.sup.i(x) as: $E_F^{\mathrm{sum}}\left(x\right)=\underset{i=1}{\overset{N}{.Math.}}{E}_F^i\left(x\right).$

Description

BRIEF DESCRIPTION OF THE FIGURES

(1) A full and enabling disclosure of the presently disclosed subject matter, including the best mode thereof, to one of ordinary skill in the art, is set forth more particularly including the specification, and including reference to the accompanying figures in which:

(2) FIG. 1 is a generally top and side perspective view of a representative first discontinuity which may be identified in accordance with the presently disclosed subject matter;

(3) FIG. 2 is a generally top plan view of a representative second discontinuity which may be identified in accordance with the presently disclosed subject matter;

(4) FIGS. 3A and 3B are an illustration and schematic diagram, respectively, of exemplary arrangements for practicing presently disclosed imaging methodology;

(5) FIG. 4 is a schematic diagram of multiple actuators placed relative to the representative first discontinuity of subject FIG. 1 herewith, for being imaged in accordance with presently disclosed subject matter;

(6) FIGS. 5A, 5B, and 5C are respective images of the representative first discontinuity of subject FIG. 1 as obtained with various embodiments of the presently disclosed subject matter;

(7) FIGS. 6A, 6B, and 6C are respective images of the representative second discontinuity of subject FIG. 2 as obtained with various embodiments of the presently disclosed subject matter; and

(8) FIGS. 7A, 7B, 7C, and 7D are schematic diagrams of respective exemplary embodiments of multiple actuator placement arrangements for use in presently disclosed subject matter.

(9) Repeat use of reference characters in the present specification and drawings is intended to represent the same or analogous features or elements or steps of the presently disclosed subject matter.

DETAILED DESCRIPTION OF THE PRESENTLY DISCLOSED SUBJECT MATTER

(10) Reference will now be made in detail to various embodiments of the presently disclosed subject matter, one or more examples of which are set forth below. Each embodiment is provided by way of explanation of the subject matter, not limitation thereof. In fact, it will be apparent to those skilled in the art that various modifications and variations may be made in the presently disclosed subject matter without departing from the scope or spirit of the subject matter. For instance, features illustrated or described as part of one embodiment, may be used in another embodiment to yield a still further embodiment, and corresponding and/or associated methodologies may be practiced relative to apparatus disclosed and/or suggested herewith, all of which comprise various embodiments of the presently disclosed subject matter.

(11) As discussed, relatively improved guided wave (GW) structural health monitoring (SHM)/nondestructive evaluation (NDE) results from practice of the network wavefield imaging methods disclosed herewith which variously combine tomography and wavefield/wavenumber imaging algorithms. Such improved methodologies result in useful imaging of relatively significantly more complex discontinuities or shapes in plate-like structures. The presently disclosed imaging techniques can provide users high-resolution images of one or more discontinuities in a structure including its location, size and shape. Further, the presently disclosed methods can provide high quality images of discontinuities with more complicated profiles.

(12) Presently disclosed algorithms advantageously account for multiple actuator-excitations from different angles enclosing a discontinuity for imaging, with exemplary such algorithms advantageously using intrinsic wave characteristics such as (but not limited to) wavefield, wavenumber, or reconstructed wave energy. Results include simultaneous improvement in both resolution (compared with prior GW imaging methods) and in functioning with circumstances where a target discontinuity is relatively more complicated.

(13) For some presently disclosed embodiments, advantages reside in detection and high-resolution quantification of a complex discontinuity in plate-like structures, as well as the ability to so detect multiple of such discontinuities since the enclosure setup of actuators can provide views from different angles for the final imaging synthesis.

(14) As referenced, various of presently disclosed embodiments for GW SHM/NDE network wavefield imaging methodologies make use of combining certain concepts of tomography with wavefield/wavenumber imaging, to obtain reliable inspection and quantification of complex discontinuities in plate-like structures. Further advantageously, in accordance with presently disclosed subject matter, subject algorithms can be implemented per different exemplary embodiments with different types of sensors and sensing systems to perform damage inspection in plate-like structures.

(15) Different exemplary algorithm embodiments are disclosed herewith for presently disclosed synthetic guided wave SHM/NDE network wavefield imaging methods. A first example relates to a presently disclosed guided wave network wavefield imaging method. An initial aspect is based on the use of wavefield data. In this instance, velocity wavefield data v(t,x) can be obtained using a scanning sensing system, where t is time and x is the spatial vector. Various scanning approaches and/or specific devices may be practiced in accordance with presently disclosed subject matter, as further discussed herein. The following portion of this disclosure initially addresses specific algorithmic aspects of the subject network wavefield imaging exemplary method as presently disclosed.

(16) Once such velocity wavefield data v(t,x) is obtained using a scanning sensing system as disclosed herein, a resulting wavefield image can be represented by amplitude, magnitude, and by an energy map. Per the presently disclosed subject matter, the terminology “magnitude” is intended to be defined as the maximum absolute value of the waveform at each point x, which can be expressed by the following Eq. (1) (with the term “equation” used interchangeably with “algorithm” herewith):
v.sub.mag(x)=max(abs(v(t,x)))  (1)

(17) When N number of actuators are employed for enclosing a possible target discontinuity (or an area or region which will be checked for the existence of any such discontinuity or discontinuities), and guided wave (GW) excitation is used from different directions in the system, N wavefield data v.sup.i(t,x) are obtained for each actuator sent excitations in round-robin. Meanwhile, using Eq. (1) herewith, N magnitude images v.sup.i.sub.mag (x) are obtained, where i=1 . . . N. The synthesis aspect of this exemplary embodiment of the presently disclosed subject matter is performed as a summation represented in this instance by Eq. (2) herewith, as follows:

(18) $\begin{array}{cc}{v}_{\mathrm{sum}}\left(x\right)=\underset{i=1}{\overset{N}{.Math.}}{v}_{\mathrm{mag}}^i\left(x\right)& \left(2\right)\end{array}$

(19) It should be understood that various alternative embodiments of the presently disclosed subject matter may make use of a synthesis process which has a final synthesis aspect based on mathematical function other than summation. For example, mathematical forms such as the use of multiplication or weighting factors may be used in place of direct summation, and all such variations are intended to come within the spirit and scope of the presently disclosed subject matter. It should be further understood that various implementations of general purpose computers and/or specific application computers, programmed with algorithms/equations disclosed herewith, may comprise various embodiments of processing means (or processors or other hardwired or programmable devices) disclosed herewith for creating a synthetic image based on such wavefield v(t,x), to detect a discontinuity in said plate structure within said scanning area.

(20) A second exemplary embodiment of guided wave SHM/NDE network wavefield imaging methodology disclosed herewith relates to a presently disclosed synthetic guided wave wavenumber imaging method. With obtained wavefield data v(t,x) comprising data in a time-space domain, space-frequency-wavenumber representation of such data in a frequency and wavenumber domain together with the spatial information can be determined per Eq. (3), as follows:

(21) $\begin{array}{cc}S\left(\stackrel{\_}{x},f,k\right)={\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }v\left(t,x\right)W\left(t,x-\stackrel{\_}{x}\right)e^{-j\left(2\pi \mathrm{ft}-k.Math.x\right)}\mathrm{dtdx}& \left(3\right)\end{array}$
where x is the retained spatial vector (x, y), f is the frequency, k is the wavenumber vector (k.sub.x,k.sub.y). W(x) is constructed using a Hanning function as,

(22) $\begin{array}{cc}W\left(x\right)=\{\begin{array}{cc}0.5\left[1+\cos \left(2\pi \frac{.Math.x.Math.}{D_x}\right)\right]& \mathrm{if}.Math.x.Math.\le D_x/2\\ 0& \mathrm{otherwise}\end{array}& \left(4\right)\end{array}$
For each spatial location, a wavenumber vector k*(x,f) is determined by corresponding to the largest value of |S(x, f, k)| as,

(23) $\begin{array}{cc}{k}^{*}\left(\stackrel{\_}{x},f\right)=\arg \underset{k}{\max }.Math.S\left(\stackrel{\_}{x},f,k\right).Math.& \left(5\right)\end{array}$

(24) At a selected frequency range, the spatial wavenumber image k*(x) is determined as follows, per Eq. (6):

(25) $\begin{array}{cc}{k}^{*}\left(\stackrel{\_}{x}\right)=\frac{1}{N}\underset{i=1}{\overset{N}{.Math.}}.Math.k^{*}\left(\stackrel{\_}{x},f_i\right).Math.& \left(6\right)\end{array}$
where f.sub.i (i=1, 2, 3 . . . N) is the selected frequency band centered at the desired frequency f.sub.0, which corresponds to the particular guided wave mode.

(26) When N number of actuators are employed for enclosing a possible target discontinuity (or an area or region which will be checked for the existence of any such discontinuity or discontinuities), and guided wave (GW) excitation is used from different directions in the system, N wavenumber images k.sub.i*(x) are obtained using the space-frequency-wavenumber representation method based on Eq. (6) herewith, where i=1 . . . N. The synthesis aspect of this exemplary embodiment of the presently disclosed subject matter is performed as a summation for the wavenumber image, represented in this instance by Eq. (7) herewith, as follows:

(27) $\begin{array}{cc}{k}_{\mathrm{sum}}^{*}\left(\stackrel{\_}{x}\right)=\underset{i=1}{\overset{N}{.Math.}}{k}_i^{*}\left(\stackrel{\_}{x}\right)& \left(7\right)\end{array}$

(28) As with the first exemplary embodiment, other forms of mathematical synthesis such as the use of multiplication or weighting factors may be used in place of the direct summation approach of Eq. (7) herewith.

(29) A third exemplary embodiment of guided wave SHM/NDE network wavefield imaging methodology disclosed herewith relates to a presently disclosed synthetic guided wave filtering reconstructed energy-based imaging method. Using such filter reconstruction imaging method, based on the wavefield data v(t,x) the frequency-wavenumber representation is determined by Eq. (8) herewith, as follows:

(30) $\begin{array}{cc}v\left(f,k\right)={\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }v\left(t,x\right)e^{-j\left(2\pi \mathrm{ft}-k.Math.x\right)}\mathrm{dtdx}& \left(8\right)\end{array}$

(31) In this instance, per the exemplary embodiment disclosed herewith, a filter F(f, k) is designed to retain the additional wavenumber components in order to obtain the resulting filtered spectrum V.sub.F(f,k). Such filtering process may be performed, for example, as the algorithm or Eq. (9) herewith, as follows:
V.sub.F(f,k)=V(f,k)F(f,k)  (9)

(32) Waves corresponding to the filtered spectrum V.sub.F(f,k) are determined by Eq. (10) herewith, as follows:

(33) $\begin{array}{cc}{v}_F\left(t,x\right)=\frac{1}{{\left(2\pi \right)}^2}{\int }_{-\infty }^{\infty }{\int }_{-\infty }^{\infty }v\left(f,k\right)e^{j\left(2\pi \mathrm{ft}-k.Math.x\right)}\mathrm{dfdk}& \left(10\right)\end{array}$

(34) An energy map to show energy change due to a discontinuity (within the enclosed, targeted area being scanned) may be constructed through use of Eq. (11) herewith, as follows:

(35) $\begin{array}{cc}{E}_F\left(x\right)={\int }_0^t\frac{1}{2}{v}_F^2\left(\tau ,x\right)d\tau & \left(11\right)\end{array}$

(36) A synthetic energy map in accordance with the presently disclosed exemplary embodiment (using a summation approach) may be obtained from the summation of the energy map E.sub.F.sup.i (x), as provided by Eq. (12) herewith, as follows:

(37) 0$\begin{array}{cc}{E}_F^{\mathrm{sum}}\left(x\right)=\underset{i=1}{\overset{N}{.Math.}}{E}_F^i\left(x\right)& \left(12\right)\end{array}$

(38) As with the first and second exemplary embodiments, other forms of mathematical synthesis such as the use of multiplication or weighting factors may be used in place of the direct summation approach of Eq. (12) herewith.

(39) Proof-of-concept experiments have been successfully performed in the context of achieving detection of complex discontinuity subject matter. Specifically, an exemplary embodiment of the presently disclosed subject guided wave SHM/NDE network wavefield imaging method shows the ability to detect and quantify a simulated complex discontinuity in a 1-mm aluminum 2024-T3 plate. In this instance, use of the term “simulated” is with reference to the fact that the subject or target being detected has been deliberately created and placed relative to the subject target plate, as opposed to being a discontinuity that inadvertently formed during fabrication of the plate. Other than that context, the “simulated” discontinuity is certainly a literal and real discontinuity relative to the otherwise existent plate.

(40) Different simulated discontinuities with different shapes such as cube, square tube, triangle bar, quartz rod and plastic letters, etc. may be bonded on a target or demonstration surface.

(41) FIG. 1 is a generally top and side perspective view of a representative first discontinuity generally 100 which may be identified in accordance with the presently disclosed subject matter. More specifically, the exemplary structure (“simulated” discontinuity) comprises an 18-mm thick equilateral triangle bar with 21-mm side lengths, glued on a target plate surface 102.

(42) FIG. 2 is a generally top plan view of a representative second discontinuity which may be identified in accordance with the presently disclosed subject matter. More specifically, such exemplary structure (“simulated” discontinuity) comprises in fact three separate structures each with 20-mm height, 10-mm thickness and about 18-mm width. As illustrated, they comprise the three letters “U”, “S”, and “C” marked respectively as structures 104, 106, and 108, each glued on plate surface 110.

(43) FIGS. 3A, 3B, and 4 represent demonstration (or experimental setups) for presently disclosed subject matter. For example, FIGS. 3A and 3B are an illustration and schematic diagram, respectively, of exemplary arrangements for practicing presently disclosed imaging methodology. FIG. 4 is a schematic diagram of multiple actuators placed relative to the representative first discontinuity of subject FIG. 1 herewith, for being imaged in accordance with presently disclosed subject matter.

(44) FIGS. 3A and 3B represent an exemplary test specimen (plate) generally 112, having a target scan area generally 114 (denoted in FIG. 3A with reflective tape). As represented particularly by FIG. 3B, a grid coordinate system x, y may be used for indicating and/or determining horizontal and vertical positions on the test specimen relative to a reference point.

(45) A non-contact scanning laser Doppler vibrometer (SLDV) generally 116 with laser head 118, as shown variously in both FIGS. 3A and 3B, may be used as a guided wave (GW) receiver over an area of 70 mm×70 mm with 1-mm spatial resolution. A function generator generally 120 may provide a trigger (such as an excitation chirp signal) generally 122 to actuate the SLDV, as well known. General operational details of using an SLDV are well known to one of ordinary skill in the art and can be found in readily available literature on such subject matter, such as follows. See, for example: Sohn, H., Dutta, D., Yang, J. Y., Park, H. J., DeSimio, M., Olson, S. and Swenson, E., “Delamination detection in composites through guided wave field image processing,” Composites science and technology, vol. 71, No. 9, pp. 1250-1256, 2011. Yu, L., Tian, Z. and Leckey, C. A., “Crack imaging and quantification in aluminum plates with guided wave wavenumber analysis methods,” Ultrasonics, vol. 62, pp. 203-212, 2015. Tian, Z., Yu, L., Leckey, C. and Seebo, J., “Guided wave imaging for detection and evaluation of impact-induced delamination in composites,” Smart Materials and Structures, vol. 24, No. 10, pp. 105019, 2015. Yu, L., Tian, Z., Li, X., Zhu, R. and Huang, G., “Core-skin debonding detection in honeycomb sandwich structures through guided wave wavefield analysis,” Journal of Intelligent Material Systems and Structures, pp. 1045389X18758180, 2018.

(46) FIG. 3B also represents placement of an exemplary actuator generally 124 relative to scan area generally 114. Thus, such FIG. 3B provides a schematic of an exemplary sensing system for acquiring wavefields of guided waves. In this instance, such exemplary actuator may serve as the representative reference point for the mentioned x, y coordinate system. Also, it may comprise in this exemplary embodiment a lead zirconate titanate (PZT) transducer. The inspection system uses a PZT wafer as the actuator to generate interrogation guided waves and an SLDV as the receiver to receive diagnosis guided waves in the scan area. The SLDV is a rapid, high spatial resolution and non-contact means for guided wave sensing and wave propagation visualization. Thus, an SLDV is one example of means for measuring a wavefield v(t,x) of guided waves in said structure generated upon excitation of the respective plurality of actuators, in accordance with presently disclosed subject matter. In particular, such means for measuring may include a non-contact scanning laser Doppler vibrometer (SLDV) for measuring the wavefield v(t,x) of guided waves in said structure, based on the Doppler effect. As known and understood by those of ordinary skill in the art, displacement and velocity components of guided waves at scan points of the SLDV can be acquired based on the Doppler effect on light waves. A laser Doppler vibrometer, even if having only one laser beam, can direct the laser beam to various locations for measurements, by employing two moving scanning mirrors driven by galvanometric actuators.

(47) The spatial resolution of scanning is related to the standoff distance (distance between laser head and test object) and the angular resolutions of the scanning mirrors. For example, for a Polytec PSV-400-M2 device, angular resolutions of the two scanning mirrors are less than 0.002°, which allow achieving high spatial resolutions. Through point-by-point measurement of guided waves at multiple scanning points on a predefined scan grid, the SLDV acquires a velocity wavefield v(t, x) of guided wave propagation as a function of both time t and space x.

(48) The presently disclosed subject matter provides the ability to perform reliable guided SHM/NDE inspections using different actuators including, but not limited to, lead zirconate titanate (PZT) transducers, pulsed lasers, electromagnetic acoustic transducers (EMAT's), or air-coupled transducer (ACT's) and the like. Likewise, while an SLDV arrangement is illustrated, presently disclosed methodologies may be practiced using different sensors including but not limited to high spatial resolution sensors SLDV, or motion controlled scanning sensors EMAT and/or ACT or the like.

(49) As noted, FIG. 4 is a schematic diagram of multiple actuators placed relative to the representative first discontinuity generally 100 of subject FIG. 1 herewith, for being imaged in accordance with presently disclosed subject matter, relative to scanning area 126 of a test specimen 128. In other words, FIG. 4 represents an experimental test setup to detect a triangle shape discontinuity generally 100 in an aluminum plate 128.

(50) As illustrated in the exemplary arrangement of present FIG. 4, in accordance with presently disclosed subject matter, there are eight (8) respective PZT transducers PZT-1 through PZT-8 (generally 130 through 144, respectively) which are bonded onto the surface of specimen 128. As shown, they are preferably in this exemplary embodiment arranged in a circle generally 146 enclosing the simulated discontinuity 100, with the specimen center (and also the circle center) generally 148 set as the coordinate origin, or (0, 0) reference. As will be understood by those of ordinary skill in the art, the actuators are used to generate guided waves. As an example only, a PZT excitation signal of a 3-cycle toneburst at 120 kHz may be used.

(51) Once the i.sup.th PZT is actuated, wavefield data v.sub.i(t,x) is obtained over the SLDV scanning area generally 126. Also, the wavefield image v.sub.mag.sup.t(x), wavenumber image k.sub.i*(x), and filter reconstructed energy-based image E.sub.F.sup.i(x) are obtained using algorithms as otherwise described herein.

(52) All eight PZT actuators 130 through 144 are, in turn, actuated for the SLDV generally 116 to collect the wavefield data and for the algorithms to, in turn, process such data to generate the images. With image data for each respective PZT in hand, the wavefield image v.sub.mag.sup.sum(x), the wavenumber image k.sub.sum*(x), and the filter reconstructed energy-based image E.sub.p.sup.sum(x) are obtained. As understood from the complete disclosure herewith, such information allows the user to quantify the overall features of a detected discontinuity, including its location, size, and shape.

(53) FIGS. 5A, 5B, and 5C are respective synthetic images of the representative first discontinuity generally 100 of subject FIG. 1 as obtained with various embodiments of the presently disclosed subject matter. More particularly, FIGS. 5A, 5B, and 5C synthetic images for the representative triangle-shaped discontinuity 100 are formed, respectively, through (a) wavefield imaging; (b) wavenumber imaging; and (c) filter reconstructed energy-based imaging methodologies of presently disclosed subject matter.

(54) The physical positions of triangular discontinuity generally 100 of subject FIG. 1 are shown as dotted-line overlays (for demonstration and disclosure purposes only) for the images of FIGS. 5A, 5B, and 5C. Specifically, dotted-line overlay representations 150, 152, and 154, respectively, represent a high correlation of the highlighted portion of the generated synthetic image with the real discontinuity generally 100, including its location, size, and shape.

(55) FIGS. 6A, 6B, and 6C are respective synthetic images of the representative second discontinuity (collectively, the respective multiple lettering structures “U” 104, “S” 106, and “C” 108) of subject FIG. 2 as obtained with various embodiments of the presently disclosed subject matter. More particularly, FIGS. 6A, 6B, and 6C synthetic images for the representative multi-component discontinuity (104/106/108) are formed, respectively, through (a) wavefield imaging; (b) wavenumber imaging; and (c) filter reconstructed energy-based imaging methodologies of presently disclosed subject matter.

(56) The highlighted portions of the resulting wavefield image (FIG. 6A) and wavenumber image (FIG. 6B) all show the subject three letters 104/106/108 very clearly and match their actual/real size, location, and shape. In this instance, the filter reconstructed energy-based image (FIG. 6C) shows the overall shape and location of the letters 104/106/108 while not as clearly in this instance indicating their exact shape. Such relatively less good imaging performance for the format of FIG. 6C is potentially caused by relatively weaker interaction between detection waves and the subject discontinuity items 104/106/108 due to the relatively small line width of the letters (at only 4 mm). In accordance with the presently disclosed subject matter, relatively more accurate imaging based on filter reconstructed energy-based imaging methodology may be obtained by arranging the subject system for generating relatively higher resolution images, such as by the number and alignment of the subject presently disclosed multiple actuator arrangement.

(57) FIGS. 7A, 7B, 7C, and 7D are schematic diagrams of respective exemplary embodiments of multiple actuator network arrangements for use in presently disclosed subject matter. In particular, FIGS. 7A, 7B, 7C, and 7D represent multiple actuator configurations which are, respectively, (a) circle, (b) square, (c) rectangle, and (d) ellipse shapes. Note the actuator arrangements are not limited to the presented configurations given here. Any network arrangement can be used as long as the discontinuity to be inspected is enclosed within.

(58) As understood from the complete disclosure herewith, the presently disclosed subject matter makes use of specified algorithms in conjunction with synthetic guided wave imaging methodologies. The presently disclosed and referenced algorithms can be implemented in SHM/NDE systems with different market actuators including, but not limited, to PZT, pulsed laser, electromagnetic acoustic transducer (EMAT), is a transducer [22], air-coupled transducer (ACT), or similar, and different high-spatial-resolution sensors such as SLDV, or motion-controlled sensors such as ACT and EMAT, or similar.

(59) Further, as represented in FIGS. 7A through 7D, actuators can be arranged in different configurations, such as circle (FIG. 7A), square (FIG. 7B), rectangle (FIG. 7C), ellipse (FIG. 7D), or polygon or other shapes. As shown in representative FIG. 7A, the plural represented actuators 156 are not only equi-distant from their adjacent actuators, but all the same distance from the reference or center point generally 158. With respect to each of FIGS. 7B, 7C, and 7D, however, the plural represented actuators are equi-distant from their adjacent actuators but not necessarily from the reference or center point. For example, in FIG. 7B, the distance between actuators 160 and 162 is the same as the distance between 162 and 164, while the distance between actuator 160 and center point 166 is different from the distance between actuator 162 and center point 166. Similar characteristics apply for FIGS. 7C and 7D. Compare, for example, equi-distances between 168 and 170 versus 170 and 172, but not equal distances between 168 to 174 compared with 170 to 174. Also, compare equi-distances between 176 to 178 versus 178 to 180, but not equal distances between 176 to 182 compared with 180 to 182.

(60) Those of ordinary skill in the art from the complete disclosure herewith will understand that various actuator arrangements may be practiced with various sensors to form particular embodiments of the presently disclosed subject matter. Further, various of the network wavefield imaging methodologies disclosed herewith may be alternatively practiced with each of such embodiments to comprise yet further alternative embodiments herewith. In particular, experience of various users may show particular embodiments relatively more useful than others for detection of discontinuities typical to the context in which they practice the presently disclosed subject matter. All such variations and selective applications are intended to come with the spirit and scope of the presently disclosed subject matter.

(61) Thus, presently disclosed subject matter provides for guided wave SHM/NDE network wavefield imaging with relatively superior resolution for complex structural discontinuities using a combined tomography and wavefield/wavenumber imaging algorithm in plate-like structures. In some instances, and in exemplary embodiments, such results are obtained through the combination of traditional and wavefield/wavenumber imaging through synthesizing algorithms and/or the use of multiple-actuators arrangements for multi-path inspection enclosing the wavefield scanning area with the discontinuity (using an approach that relatively more information can be captured using more actuators). Still further, various present exemplary embodiments make use of guided wave synthesis imaging using wavefield inspection data from different views, and others provide for complex discontinuity quantification from different angles enclosing the discontinuity. Thus, the presently disclosed subject matter in some instances helps capture discontinuity feature or features from different angles of view, for improving high-resolution synthesis images for highlighting overall discontinuity information including location, size, and shape.

(62) Further, potential industrial applications of the presently disclosed subject matter are applicable to aviation, aerospace, and automotive engineering, where plate-like materials are heavily used. Presently disclosed guided wave network wavefield imaging methodology provides a way to quantify more complex discontinuity (singular or plural) in plate-like structures. Moreover, the presently disclosed subject matter for network wavefield imaging methodologies provides a more interactive way of online structural health mentoring of plate-like structures. By such application types, the presently disclosed network wavefield imaging methodologies can be used for reliable discontinuity inspection for plate-like structures, online structural health monitoring of plate-like structures, and quality control of plate-like structures in product design and manufacturing processes. That makes the presently disclosed subject matter potentially useful in many circumstances, such as in federal and industrial laboratories, plate-like structure manufactures, aerospace, automotive, energy generation, and related industries that are required to assure the safety of their products by structural health monitoring and nondestructive evaluation. All of the foregoing variations and potential applications are intended to come within the spirit and scope of the presently disclosed subject matter.