Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Near-Infrared Tunable Laser Absorption Spectroscopic Acetylene Sensor System Using a Novel Three Mirror-Based, Dense Pattern Gas Cell
Previous Article in Journal
Low Energy Consumption Compressed Spectrum Sensing Based on Channel Energy Reconstruction in Cognitive Radio Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sector Piezoelectric Sensor Array Transmitter Beamforming MUSIC Algorithm Based Structure Damage Imaging Method

1
Institute of Chemical Materials, China Academy of Engineering Physics, Mianyang 621900, China
2
Research Center of Structural Health Monitoring and Prognosis, State Key Lab of Mechanics and Control of Mechanical Structures, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
*
Authors to whom correspondence should be addressed.
Sensors 2020, 20(5), 1265; https://doi.org/10.3390/s20051265
Submission received: 18 December 2019 / Revised: 19 February 2020 / Accepted: 25 February 2020 / Published: 26 February 2020
(This article belongs to the Section Physical Sensors)

Abstract

:
Elastic-wave-based structural health monitoring technology has a broad application potential for its sensitivity and ability to achieve regional monitoring. For structures with large damping and specific shapes, the traditional damage monitoring method is limited by the sensor arrangement area and affected by low signal-to-noise ratios, so it is difficult to accurately locate the damage in a structure. To solve this problem, this paper proposed a damage monitoring method based on a sector piezoelectric sensor array for multiple signal classification algorithm. By arranging two sector piezoelectric sensor arrays that are suitable for a specific structure, the damage scattering array signal under the multi-excitation source was obtained and synthesized, the signal-to-noise ratios were improved, and the damage location accuracy was thus improved. The effectiveness of the method was verified by monitoring the damage in a circular bonded structure with a metal ring. Compared with the damage localization methods based on the traditional single excitation source multiple signal classification algorithm, path imaging and delay-sum imaging, this method can achieve better damage location and has a higher localization accuracy.

1. Introduction

High-reliability, high-performance structural safety services are of great importance to major or important mechanical equipment and infrastructure in fields such as aerospace, civil engineering, and rail transit. These services also challenge the traditional structural safety assessment and maintenance guarantee [1,2,3]. The application of structural health monitoring (SHM) technology, which has gradually become a popular research area, can effectively prevent major accidents, achieve the maintenance of conditions, and reduce maintenance costs [4,5,6,7]. Among the SHM methods, the elastic-wave-based SHM method is regarded as a promising damage monitoring method because it is sensitive to metal cracks and the impact damage of composite materials, and can be used to monitor structures in real time. It is considered to be a promising damage monitoring method [8,9,10,11,12,13,14,15]. Piezoelectric sensors are commonly used to excite and receive elastic waves in structures, and thus, damage monitoring in structures can be combined with the damage imaging method [16]. According to different sensor arrangement methods, damage imaging methods can be divided into sparse sensor array imaging and dense sensor array imaging.
In the research on imaging methods based on sparse sensor arrays, Michaels et al. [17] used a multisensing method to acquire structural response signals and image fusion technology to realize hole damage imaging in aluminum plates. Yuan et al. [18] combined the delay-sun imaging method with the probabilistic hybrid model and proposed a 4D damage imaging method. The damage imaging of a carbon fiber reinforced structure under the influence of temperature is realized. Chang et al. [1] proposed an energy-factor-based damage factor calculation method for each excitation-sensing path. The delamination damage of a reinforced composite plate and the crack damage of a metal aluminum plate were identified. Qing et al. [19] proposed a damage factor calculation method based on signal cross-correlation and introduced intelligent sandwich technology to monitor the corrosion of metal pipelines. According to the phenomenon that damage has the greatest influence on the direct wave packet, Su et al. [20] used image fusion technology to image all possible points in a structure and introduced the cumulative distribution function to achieve damage prediction. Yuan et al. [21,22] proposed a time reversal method based on the dual PZT arrangement, which could effectively eliminate the influence of the environment and boundary reflection and realize nonreference damage imaging. The team then studied a type of high spatial resolution imaging based on virtual time reversal, which could clearly display both the single damage and the dual adjacent damage in the plate. Shen et al. [23] proposed an enhanced Lamb wave virtual time reversal technique for damage detection with transducer transfer function compensation. The method can achieve more accurate and robust damage imaging results. Fan et al. [24] used a self-designed tomographic imaging system based on a piezoelectric sensor array to achieve an accurate reconstruction of corrosion damage in the plate structures. Shen et al. [25] used tomographic imaging technology to monitor the cracks in the holes of structures and compared the tomographic imaging performance under different sensor array configurations.
However, the sparse array layout will be affected by region constraints and structural complexity. As a result, the sparse sensor arrangement is no longer applicable, and the localization accuracy of the algorithm is reduced. With the deepening of research, elastic wave damage monitoring methods based on dense sensor arrays and related array signal processing methods have attracted more and more attention [26,27,28]. As a new array signal processing method, the multiple signal classification (MUSIC) algorithm combined with the elastic wave monitoring method has been introduced for structural health monitoring due to its directional scanning capability and the potential of multisource monitoring [29,30]. Early MUSIC algorithms were mostly used for the angular localization of impact sources [31,32,33,34,35]. In recent years, Zhong and Yuan [36] proposed a near-field MUSIC algorithm based on a linear sensor array, which solved the near-field dead zone problem of the far-field MUSIC algorithm and verified the damage location on the composite plate structure.
For some high-damping materials or high attenuation structures with special structural forms, the elastic wave attenuates rapidly in the process of propagation, which means the signal-to-noise ratios (SNRs) of the damage scattering signals will be low, and thus, the damage location accuracy will be affected. For circular pipeline structures, Golato et al. [37] proposed a sparse array pipeline surface damage monitoring method based on two circular sensor rings. In this method, the linear signal model is reconstructed by grouping, and the damage location in a circular pipeline structure is realized. Chattopadhyay et al. [38] introduced elastic waves in their study and used space–time analysis without a reference method to achieve the damage localization of an X-COR sandwich high attenuation composite structure. Qu et al. [39] proposed a new damage detection technique for use in highly damped composite structures using the time reversal method with different excitation amplitudes. Qing et al. [40] proposed a damage monitoring method for torsional elastic waves, which realized the monitoring of high attenuation honeycomb structures. Bao et al. [41,42] proposed the transmitter beamforming MUSIC algorithm based on a bilinear array and achieved low SNRs corrosion damage monitoring in a reinforced aviation structure. The array error was also corrected, and thus, the localization accuracy of the algorithm was improved. The algorithm was verified on an anisotropic high attenuation composite reinforced structure.
In summary, the damage imaging method based on a sparse sensor array is susceptible to structural limitations, while the single excitation source MUSIC damage location method based on a dense sensor array will lead to localization errors due to the problem of low SNRs. In this paper, a transmitter beamforming MUSIC damage imaging method based on a sector sensor array is proposed to perform damage monitoring of a special circular bonded structure. In this method, according to the specific form of the circular bonded structure, the sector-shaped double sensor array is arranged, and the damage scattering array signal is superimposed at the search point according to the time delay rule to improve the SNRs. Through the eigenvalue decomposition of the covariance matrix of the superimposed sensor array signal, direction scanning is carried out on the structure, and the spatial spectrum of the monitoring area is constructed to locate the damage.
The paper consists of the following parts. Section 2 introduces the basic principle and damage monitoring process of the sector piezoelectric sensor array transmitter beamforming MUSIC (SPATB-MUSIC) algorithm. Section 3 studies the propagation properties of elastic waves in a circular bonded structure with a metal ring. In Section 4, the damage location experiment was performed on a circular bonded structure, and the localization details and results are given. Section 5 presents the conclusion.

2. SPATB-MUSIC Algorithm

In this section, based on the transmitter beamforming and MUSIC algorithm, a SPATB -MUSIC algorithm is introduced in detail. The method uses a dual array consisting of two sector sensor arrays. The method then produces an elastic wave by driving one of the sensor array rows and receives elastic waves with the other sensor array row. Then, in the search for the monitoring area, the focus signal is obtained by making the elastic wave of each element in the sensor array reach the damage at the same time. The focus array signal is used as the input of the MUSIC algorithm to obtain the noise subspace and the steering vector space, the spatial spectrum is used to describe the orthogonality between them, and the damage location is realized by jointly searching the entire monitoring region.

2.1. Elastic Wave Signal Propagation Model under Sector Piezoelectric Array

In the sound source estimation method, it can be divided into far field and near field, according to the distance of the sound source to the sensor array. In the far field model, the wavefront can be approximated as the linear wavefront, and the location of damage can be determined only by searching the angle. However, for the structure monitored in this paper, the distance between the damage and the sensor array is short. If the approximation is made in the far field, a large positioning error will occur. Therefore, the near field model is introduced. Figure 1 shows the near field elastic wave propagation model of a single sector piezoelectric array. There are a total of 2M + 1 elements, which are located on the circumference of the radius R, and the distance between the adjacent elements is the same. The coordinate system is established by taking the central array element as the coordinate origin and the reference array element. The coordinates of the ith array element are (Ri, wi), where Ri is the distance between the ith array element and the reference array element, and wi is the clip with the x-axis angle. In the near-field case, the wavefront cannot be approximated as the linear wavefront but as the circular wavefront, and the distance between the signal source and the array elements of the sensor array is different, while the angles of the array elements with respect to the x-axis are not consistent.
If the center frequency of the signal is ω0, the response signal x0(t) of the reference element PZT0 can be expressed as
x 0 ( t )   =   u ( t ) e j ( ω 0 t k r )
where u(t) is the damage scattering signal; r is the distance from the damage to the reference array element; k is the wavenumber (k = ω0/c); and c is the elastic wave propagation velocity. Then, the response signal xi(t) of the array element PZTi in the array can be represented as
x i ( t ) = x 0 ( t ) e j ω 0 τ i + n i ( t ) , i = M , , M
where ni(t) is the background noise, and τi is the time difference caused by different arrival times of the damage scattered signal due to the different distances between PZTi and PZT0 from the damage. The method of calculation is as follows:
τ i = Δ r i c = r r i c
The path propagating from the signal source to PZT0 and PZTi and the connection with PZT0 and PZTi can form a triangle. According to the cosine theorem, the distance ri from the signal source to each element of the piezoelectric sensor array can be expressed as
r i = r 2 + R i 2 2 r R i cos ( θ w i ) , i = M , M + 1 , , M
and ri is not parallel to r. In this equation, θ is the angle between the damage source and the reference element and the x-axis, and in Equation (4) instead of Equation (3), the time delay can be rewritten as
τ i = r r 2 + R i 2 2 r R i cos ( θ w i ) c , i = M , M + 1 , , M
From this equation, it can be found that the time delay is related not only to the sound source angle, but also to the distance. The guidance vector of the signal source in the near field case is defined as
a i r , θ = e j ω 0 τ i
Writing the array response signal in the form of a vector,
X t = x M t x M + 1 t x M t T A r , θ = a M r , θ a M + 1 r , θ a M r , θ T X t = A r , θ s 0 t + N t ,   N t = n M t n M + 1 t n M t T

2.2. Transmitter Beamforming under Dual Sector Piezoelectric Array

In practical applications, for the bonding structure monitored in this paper, the thickness and damping are large, and the SNRs of the damage scattering signal is small, which leads to the error of the final localization result. To solve this problem, a double sector piezoelectric array was used, and a transmitter beamforming MUSIC algorithm based on the sector piezoelectric array was proposed. By arranging double sector array sensors, the distance damage delay of the excitation array elements can be calculated. The damage scattering signal is superimposed to enhance the signal and effectively improve the SNRs of the damage scattering signal.
As shown in Figure 2, the elastic wave propagation models for arranging the double sector array sensors are named as array A and array S, respectively, where array A is the excitation source array, and array S is the sensor array. The number of elements of each sector array is 2M + 1, each element is named from left to right, the excitation source array is AM to AM, and the sensor array is S−M to SM. The coordinates of element p in the excitation array are R p A , w p A , and r p A is the distance from the pth array element of the excitation source array to the damage.
Therefore, according to Equation (8), the distance from each element of the excitation source array to the damage location and the time delay of each element signal arriving at the damage relative to the reference element signal can be calculated.
r p A = R p A cos w p A r cos θ 2 + R p A sin w p A r sin θ 2 t p = Δ r p A c = r p A r 0 A c , p = M , , M
The focused damage scattering signal can be represented by Equation (9). Xp(t) is the damage scatter signal when the pth array element is excited in the array. At this time, X(t) is an enhanced damage scattering array signal.
X t = X M t e j ω 0 t M + + X p t e j ω 0 t p + + X M t e j ω 0 t M

2.3. MUSIC Algorithm

The eigenvalue decomposition is performed by using the enhanced array signal as an input of the MUSIC algorithm as shown in Equation (10):
R = E X X H = A E S S H A H + E N N H = A R S A H + R N
RS and RN are the signal covariance matrix and the noise covariance matrix, respectively. Assuming that the noise power is σ2, Equation (10) is transformed into
R = A R S A H + R N = A R S A H + σ 2 I
The eigen decomposition of the array signal covariance matrix is as follows:
R = U Σ U H
where U is an eigenvector matrix, and Σ is a diagonal matrix composed of eigenvalues, as shown in Equation (13).
Σ = λ 1 λ 2 λ 2 M + 1 , λ 1 λ 2 λ K λ K + 1 = = λ 2 M + 1 = σ 2
The diagonal matrix composed of large eigenvalues is defined as ΣS, and the diagonal matrix composed of small eigenvalues is defined as ΣN.
Σ S = λ 1 λ 2 λ K Σ N = λ K + 1 λ K + 2 λ 2 M + 1
The eigenvector corresponding to ΣS is then tensioned into a signal subspace US = [e1 e2 ek], the eigenvector corresponding to ΣN is tensioned into a noise subspace UN = [ek+1 ek+2 e2M+1], and Equation (12) can be further written as Equation (15).
R = U S Σ S U S H + U N Σ N U N H
As the signal and noise are independent of each other, the noise subspace and the signal subspace are orthogonal, and the guidance vector of the signal subspace is also orthogonal to the noise subspace.
A H r , θ U N = 0
Due to the existence of the background noise, A(r, θ) and UN are not completely orthogonal. Therefore, Equation (17) is used to realize the minimum optimal search location.
r ^ , θ ^ = argmina A H r , θ U N U N H A r , θ
Therefore, the spectral estimation equation of the MUSIC algorithm is defined as
P MUSIC r , θ = 1 A H r , θ U N U N H A r , θ
The entire monitoring area is searched to obtain the spatial spectrum on the monitoring area. When the searched location is the damage location, a spectral peak will appear in the spatial spectrum because the strongest array scattering signal is obtained. The damage location can be obtained by Equation (19).
r d , θ d = argmax P MUSIC r , θ
In summary, the basic procedure of the damage monitoring method of the transmitter beamforming MUSIC algorithm based on the sector piezoelectric array can be divided into the following four parts, as shown in Figure 3:
(1)
The signal of the structural damage scattering array is obtained.
(2)
By using the transmitter beamforming equation, the damage scattering array signal is obtained.
(3)
The synthesized damage scattering array signal is decomposed by eigenvalue decomposition, and the orthogonal signal subspace and noise subspace are obtained.
(4)
The spatial spectrum of the corresponding monitoring region is obtained by spectral peak searching, and the corresponding coordinates of the spectral peak are the results of the damage location.

3. Propagation Characteristics of Elastic Waves in a Circular Bonded Structure

In this section, the circular bonded structure is introduced in detail. The elastic wave propagation characteristics experiment in this structure is described in detail. The amplitudes and wave velocities of the signals in a typical sensing channel are counted.

3.1. Circular Bonded Structure and Experimental Setup

The experimental setup includes an integrated elastic wave signal scanning system and a circular bonded structure with a smart piezoelectric sensor array layer attached, as shown in Figure 4. The integrated scanning system is used to generate an excitation signal and to obtain a response signal from the PZT sensor output [43]. In the experiment, a sinusoidal modulation five-peak signal was used as the excitation signal. The amplitude was ±70 V, the central frequency was set to 50 kHz, and the sampling rate was set to 10 MSPS. A schematic view of the dimensions of the structure, which is formed by bonding plexiglass and steel with a two-component epoxy adhesive, is shown in Figure 5. The plexiglass part was cylindrical and had a thickness of 35 mm and a radius of 45 mm. There were cylindrical grooves with a radius of 15 mm, and a hemispherical groove with a radius of 10 mm in the center. The inner diameter and outer diameter of the steel ring were 10 mm and 15 mm, respectively, and the thickness was 2 mm.
A group of piezoelectric sensors numbered from 1 to 18 were attached to the surface of the test piece. The diameter of the sensor was 8 mm, and the thickness was 0.48 mm. The center of the sensor is located on a circumference with a radius of 40 mm. During the monitoring process, sensors numbered from 1 to 7 and from 12 to 18 are selected to form a double sector array, named array A and array S, respectively. Each array contained seven piezoelectric sheets, and array A was used to excite the elastic waves. Array S was used as the sensing receiving elastic wave signal. The distance between adjacent sensors in the array was 9 mm, and a total of 49 sensing channel signals were obtained.

3.2. Propagation Characteristics of Elastic Waves in Circular Bonded Structures

This section investigates the propagation properties of elastic waves in a circular bonded structure. To obtain the propagation characteristics of the elastic wave in the structure, the elastic wave signals in the form of three channels were selected for analysis, as shown in Figure 6. The signals in the PZT A1-PZT S18 channel did not pass through the bonding and variable thickness region, while the signals in the PZT A1-PZT S15 channel passed through the bonding area, but did not pass through the variable thickness region. The signals in the PZT A1-PZT S12 channel passed through the bonded and variable thickness regions. Figure 7 shows the waveform in the form of three channels. The first wave packet of the signal was not analyzed due to the crosstalk caused by the line, followed by the signal direct wave. When the path length of the PZT A1-PZT S18 channel was only 62.4 mm, the direct band was more obvious, and the crosstalk was lightly aliased. When the path length of the PZT A1-PZT S15 channel was 75.5 mm, the signal direct wave was clear and no aliasing occurred. As the path length of the PZT A1-PZT S12 channel was 80.0 mm, the boundary forms were very complex. When the elastic wave propagates in the structure, it encounters the structure boundary, which results in the reflection of the signal and generates the boundary reflection wave. When the time difference between the direct wave and the boundary reflection wave is small, the waveform will be aliased, as shown in Figure 7c. However, the signals used in this paper were damage scattering signals; most crosstalk and boundary reflections were suppressed; and aliasing had little effect on the calculation of wave velocity. Therefore, the effect of aliasing on the final imaging results was small.
To further study the propagation law of elastic wave signals in the structure, the propagation velocity of the signal direct wave was calculated, and the Shannon plural wavelet transform method was used to obtain the excitation signal and response signal envelope [9]. As shown in Figure 8, the envelope of the PZT A1-PZT S18 channel signal was extracted, the flight time of the signal was obtained by the peak corresponding to time subtraction, and the signal propagation speed was calculated in combination with the direct path length of the signal.
The calculation results of the velocity are given in Table 1, and the amplitude of the direct wave crest of the signal was statistically analyzed. The amplitude of the signals in the PZT A1-PZT S15 channel, which passed through the bonded area without passing through the region of variable thickness, was 0.814 V. The amplitudes of the other two channels were only approximately half of their amplitudes at 0.486 V and 0.383 V, and the wave velocity was 699 m/s. The wave velocities of the other two channels were 722 m/s and 755 m/s, and the velocity difference was small. It can be seen that when the elastic wave signal propagates in this kind of bonding structure, the amplitude is small, the difference of the signal propagation speed is small, and the average velocity of signal propagation can be used for damage location in the process of monitoring. According to the above velocity calculation method, the average velocity of the 49 channels was 736 m/s, which was used for damage location in the next section.

4. Damage Location in Circular Bonded Structures

In this section, the SPATB-MUSIC algorithm is used to evaluate the damage in the circular bonded structure, and this method is compared with the traditional single excitation source MUSIC algorithm, the path imaging method, and the delay-sum imaging method.

4.1. Damage in the Structure and Monitoring

As shown in Figure 9, the part indicated by the red line is the location of the damage, and four damages were constructed in the structure by peeling the metal bonded to the plexiglass. As shown in Figure 10a, in undamaged conditions, the elastic wave array signal in the structure is obtained as a reference signal, and then, as shown in Figure 10b–e, the damage is artificially constructed in the structure in turn. The elastic wave array signal is obtained after the damage is generated at each place of the structure, and the signal in the former damage state is used as the reference signal in the latter damage state. The signal excitation acquisition system and parameter setting were the same as those in Section 3.1. The elastic waves were excited by the sensors in array A, and the sensors in array S were used to obtain the response signals under the 49 sensing channels.
Table 2 shows the damage information. A total of four damages were constructed, which were numbered D1, D2, D3, and D4. To verify the monitoring ability of the proposed method for different sizes of damage, dimensions of 10 × 5 mm2 and 5 × 5 mm2 were set. The dimensions of D1, D3, and D4 were 10 × 5 mm2, and the dimensions of D2 were only 5 × 5 mm2. The specific location of the damage is shown in the table.

4.2. Effect of Damage on the Signal

To study the influence of damage on the signal in the structure, Figure 11 shows the reference signals and damage signals for the PZT A4-PZT S18 and PZT A4-PZT S14 channels with damage sizes of 10 × 5 mm2 and 5 × 5 mm2, respectively. As the damage is located on the direct path of the signal, the signal given here is the direct part of the signal. It can be seen that when there was damage with a size of 10 × 5 mm2 in the structure, the damage signal changed with respect to the reference signal, and the amplitude increased. Although the damage scattering signal had clear wave packets, the SNRs were low. The amplitude only changed by approximately 0.070 V. For damage with a size of 5 × 5 mm2, the signal change was smaller, the damage scattering signals were aliased, and the wave packet was not obvious.

4.3. Transmitter Beamforming of Damage Scattering Array Signal

To study the beamforming effect, a comparison chart of the damage scattering array signals for damage sizes of 10 × 5 mm2 and 5 × 5 mm2 before and after enhancement is presented. As shown in Figure 12a, when the damage size was 10 × 5 mm2, PZT A4 was used as the excitation sensor, and array S was used to receive the array signals. Although there were obvious damage scattering segments and the wavefront was formed, the amplitude of the signal was less than that of the signal after the scatter segment. The damage scattering array signal shown in Figure 12b is the result of beamforming. It can be seen that the damage scattering segment was enhanced, and the PZT S12 receiving signal was removed. The remaining sensor receiving array signals had the largest amplitude of the damage scattering segment. As shown in Figure 13a, for a damage size of 5 × 5 mm2, the PZT A4 was used as the excitation sensor, and array S was used as the received damage scattering array signal. It can be seen that the amplitude of the damage scattering section was small at this time, and no obvious wavefront could be observed. The signal wave packets were mixed. Conversely, the damage scattering array signal obtained by transmitter beamforming enhanced the damage scattering section and formed a more obvious wavefront. However, time of flight (ToF) changed after beamforming, for which there are two reasons. (1) In the beamforming algorithm, the average speed of signal propagation in all channels is used to calculate the ToF when the signal reaches the damage. However, there are differences in the signal propagation speed in each channel, resulting in changes in ToF after beamforming. (2) Before beamforming, the damage sensitive segment of the damage scattering signal has a small amplitude and is aliased with the boundary reflections, while the array signal after beamforming is synthesized by seven groups of array signals, so that the damage sensitive segment is enhanced and becomes obvious, resulting in the visual changes of ToF before and after beamforming.

4.4. Damage Localization Result

Using PZT A4 as the excitation source, the damage scattering array signal obtained by the difference between the reference signal received by array S and the damage signal was used as the input of the MUSIC algorithm. The angle search range was from 0° to 180°, and the distance was from 0 cm to 8 cm. The lengths were set to 1° and 0.1 cm, respectively, and the results of the four damage locations are shown in Figure 14. The yellow circle is the actual damage location, and the maximum pixel value in the cloud image is the actual localization result. The localization errors were (3°, 0.8 cm), (1°, 1.0 cm), (11°, 0.1 cm), and (6°, 0.1 cm). It can be seen that there were errors in the location of the four damages. The angle error was up to 11°, the distance error was up to 1.0 cm, and the background noise in the 5 × 5 mm2 damage location cloud map was relatively large.
Figure 15 shows the four damage localization cloud maps using the SPATB-MUSIC algorithm. Table 3 shows the localization results and errors of the four damages. In this condition, sensor array A was used as the excitation source, and array S was taken as the reception. The four localization errors were (0°, 0.3 cm), (1°, 0.5 cm), (1°, 0.2 cm), and (0°, 0.2 cm). It can be seen that the localization accuracy for the damage was higher using the method, and the localization accuracy was greatly improved relative to the single excitation source. The angle error was only 1° at the maximum, the distance error was only 0.5 cm, and the background noise of the damage localization cloud image of 5 × 5 mm2 was effectively suppressed.
The results show that this method could also accurately locate the damage of the structure in the case of a low SNRs.
The localization results of two of the more commonly used damage localization methods for the damage size of 5 × 5 mm2 are also given. One is the path imaging method combined with the damage factor, which assigns the probability of the points in the ellipse monitoring area formed by each path and then superimposes the attenuation factor extraction results of each sensing channel. The damage location imaging result is shown in Figure 16a. The localization result was (92°, 3.7 cm). This was 1.5 cm away from the actual location, and the error was 18.8% compared with the monitoring area diameter of 8 cm. The calculation method is as shown in Equation (20),
E r r o r σ = E r r o r d R
where Errord is the distance from the damage localization result and the actual damage location; R is the diameter of the monitoring area; and Errorσ is the relative error. In the other method, the delay is calculated by using the propagation velocity of the elastic wave in the structure, and the wave packet amplitude of the damage scattering signal was stacked according to the time delay rule to obtain a pixel map reflecting the location of the damage. When the corresponding point of the actual damage location is searched, the pixel value increases abnormally, and thus, the damage location is realized. The localization result is shown in Figure 16b. The localization result was (90°, 6.3 cm) and the relative error was 13.8%. For comparison, the localization result of the STAPB-MUSIC method was also presented as shown in Figure 16c. It can clearly be seen that in the imaging diagram of the proposed method in this study, the pixel focus location was closer to the actual damage location. Table 4 shows the damage localization errors between the methods proposed in this paper and the above two imaging methods. It can be seen that the SPATB-MUSIC algorithm has a higher damage location accuracy than the two traditional imaging methods, with a relative error of only 6.3%.

5. Conclusions

The small signal variation due to damage and structural limitations causes difficulties in the traditional structural damage imaging method. In this paper, a structural damage monitoring method based on the SPATB-MUSIC algorithm was proposed for circular structures. In this method, a double sector piezoelectric sensor array is arranged on the surface of the structure. The elastic wave array signal under a multi-excitation source is obtained. The damage scattering array signal is obtained by subtracting the reference signal from the damage signal, and the transmitter beamforming method is then used to enhance the damage scattering signal and improve the signal SNRs. The enhanced damage scattering array signal is used as the input of the MUISC algorithm to realize the accurate location of the damage. The effectiveness of the method was verified by a damage monitoring experiment on a large circular damping structure formed by bonding plexiglass and metal materials. The experimental results showed that the proposed method has a higher localization accuracy than the traditional single transmitter source MUSIC algorithm, the path imaging method, and the delay-sum imaging method. The angular localization error was less than 1°, and the distance error was less than 0.5 cm. Although the STAPB-MUSIC method proposed in this paper uses the direct wave to locate the damage of the structure with high accuracy, there is still a lack of further research on the propagation characteristics of the signals and wave modes in the structure. In the future, we will carry out more detailed research on this aspect.

Author Contributions

The work was finished by the cooperation of all authors. All authors gave detailed suggestions to improve the work. L.Q. proposed the idea and the principle of the SPATB-MUSIC damage localization method. Y.W. and T.F. programmed the algorithm. Y.W. and X.T. performed the relevant experimental validation. This paper was written by T.F. and Y.W. All the authors read and approved the final paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 51921003, Grant No. 11604309), the Program for Distinguished Talents of Six Domains in Jiangsu Province of China (Grant No. GDZB-035), the Fok Ying Tung Education Foundation of China (Grant No. 161048), and the Priority Academic Program Development of Jiangsu Higher Education Institutions of China.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ihn, J.B.; Chang, F.K. Pitch-catch active sensing methods in structural health monitoring for aircraft structures. Struct. Health Monit. 2008, 7, 5–19. [Google Scholar] [CrossRef]
  2. Brownjohn, J. Structural health monitoring of civil infrastructure. Philos. Trans. Math. Phys. Eng. Sci. 2007, 365, 589–622. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Diamanti, K.; Soutis, C. Structural health monitoring techniques for aircraft composite structures. Prog. Aerosp. Sci. 2010, 46, 342–352. [Google Scholar] [CrossRef]
  4. Qiu, L.; Yuan, S.; Mei, H.; Fang, F. An improved Gaussian mixture model for damage propagation monitoring of an aircraft wing spar under changing structural boundary conditions. Sensors 2016, 16, 291. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Yuan, S.; Lai, X.; Zhao, X.; Xu, X.; Zhang, L. Distributed structural health monitoring system based on smart wireless sensor and multi-agent technology. Smart Mater. Struct. 2006, 15, 1–8. [Google Scholar] [CrossRef]
  6. Wu, J.; Yuan, S.; Ji, S.; Zhou, G.; Wang, Y.; Wang, Z. Multi-agent system design and evaluation for collaborative wireless sensor network in large structure health monitoring. Expert Syst. Appl. 2010, 37, 2028–2036. [Google Scholar] [CrossRef]
  7. Giurgiutiu, V.; Zagrai, A.; Bao, J.J. Piezoelectric wafer embedded active sensors for aging aircraft structural health monitoring. Struct. Health Monit. 2002, 1, 41–61. [Google Scholar] [CrossRef]
  8. Su, Z.; Ye, L.; Lu, Y. Guided Lamb waves for identification of damage in composite structures: A review. J. Sound Vib. 2006, 295, 753–780. [Google Scholar] [CrossRef]
  9. Qiu, L.; Yuan, S.; Zhang, X.; Wang, Y. A time reversal focusing based impact imaging method and its evaluation on complex composite structures. Smart Mater. Struct. 2011, 20, 105014. [Google Scholar] [CrossRef]
  10. Raghavan, A. Guided-wave structural health monitoring. Diss. Abstr. Int. 2007, 975, 91–114. [Google Scholar] [CrossRef]
  11. Chen, J.; Yuan, S.; Jin, X. On-line prognosis of fatigue cracking via a regularized particle filter and guided wave monitoring. Mech. Syst. Signal. Process. 2019, 131, 1–17. [Google Scholar] [CrossRef]
  12. Chen, X.; Michaels, J.E.; Michaels, T.E. A methodology for estimating guided wave scattering patterns from sparse transducer array measurements. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2015, 62, 208–2019. [Google Scholar] [CrossRef] [PubMed]
  13. Yuan, S.; Ren, Y.; Qiu, L.; Mei, H. A multi-response-based wireless impact monitoring network for aircraft composite structures. IEEE Trans. Ind. Electron. 2016, 63, 7712–7722. [Google Scholar] [CrossRef]
  14. Qing, X.; Li, W.; Wang, Y.; Sun, H. Piezoelectric transducer-based structural health monitoring for aircraft applications. Sensors 2019, 19, 545. [Google Scholar] [CrossRef]
  15. Croxford, A.J.; Wilcox, P.D.; Drinkwater, B.W.; Konstantinidis, G. Strategies for guided-wave structural health monitoring. Proc. R. Soc. A Math. Phys. Eng. Sci. 2007, 463, 2961–2981. [Google Scholar] [CrossRef]
  16. Mei, H.; Mohammad, F.H.; Joseph, R.; Migot, A.; Giurgiutiu, V. Recent advances in piezoelectric wafer active sensors for structural health monitoring applications. Sensors 2019, 19, 383. [Google Scholar] [CrossRef] [Green Version]
  17. Michaels, J.E.; Michaels, T.E. Guided wave signal processing and image fusion for in situ damage localization in plates. Wave Motion. 2007, 44, 482–492. [Google Scholar] [CrossRef]
  18. Ren, Y.; Qiu, L.; Yuan, S.; Fang, F. Gaussian mixture model and delay-and-sum based 4D imaging of damage in aircraft composite structures under time-varying conditions. Mech. Syst. Signal. Process. 2020, 135, 106390. [Google Scholar] [CrossRef]
  19. Qing, X.; Beard, S.; Shen, S.B.; Banerjee, S.; Bradley, I.; Salama, M.M.; Chang, F.K. Development of a real-time active pipeline integrity detection system. Smart Mater. Struct. 2009, 18, 115010. [Google Scholar] [CrossRef]
  20. Su, Z.; Cheng, L.; Wang, X.M.; Yu, L.; Zhou, C. Predicting delamination of composite laminates using an imaging approach. Smart Mater. Struct. 2009, 18, 74002. [Google Scholar] [CrossRef] [Green Version]
  21. Wang, Q.; Yuan, S. Baseline-free imaging method based on new PZT sensor arrangements. J. Intell. Mater. Syst. Struct. 2009, 20, 1663–1673. [Google Scholar]
  22. Cai, J.; Shi, L.; Yuan, S.; Shao, Z. High spatial resolution imaging for structural health monitoring based on virtual time reversal. Smart Mater. Struct. 2011, 20, 055018. [Google Scholar] [CrossRef]
  23. Wang, J.; Shen, Y. An enhanced Lamb wave virtual time reversal technique for damage detection with transducer transfer function compensation. Smart Mater. Struct. 2019, 28, 085017. [Google Scholar] [CrossRef]
  24. Rao, J.; Ratassepp, M.; Lisevych, D.; Caffoor, M.H.; Fan, Z. On-line corrosion monitoring of plate structures based on guided wave tomography using piezoelectric sensors. Sensors 2017, 17, 2882. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Wang, S.; Wu, W.; Shen, Y.; Liu, Y.; Jiang, S. Influence of the PZT sensor array configuration on Lamb wave tomography imaging with the RAPID algorithm for hole and crack detection. Sensors 2020, 20, 860. [Google Scholar] [CrossRef] [Green Version]
  26. Wilcox, P. Omni-directional guided wave transducer arrays for the rapid inspection of large areas of plate structures. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2003, 50, 699–709. [Google Scholar] [CrossRef]
  27. Malinowski, P.; Wandowski, T.; Trendafilova, I.; Ostachowicz, W. A phased array-based method for damage detection and localization in thin plates. Struct. Health Monit. 2009, 8, 5–15. [Google Scholar] [CrossRef] [Green Version]
  28. Yoo, B.; Purekar, A.S.; Zhang, Y.; Pines, D.J. Piezoelectric-paint-based two-dimensional phased sensor arrays for structural health monitoring of thin panels. Smart Mater. Struct. 2010, 19, 075017. [Google Scholar] [CrossRef]
  29. Bencheikh, M.L.; Wang, Y. Joint DOD-DOA estimation using combined ESPRIT-MUSIC approach in MIMO radar. Electron. Lett. 2010, 46, 1081–1083. [Google Scholar] [CrossRef]
  30. Santosh, S.; Sharma, K. A review on multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 2003, 34, 276–280. [Google Scholar]
  31. Engholm, M.; Stepinski, T. Direction of arrival estimation of Lamb waves using circular arrays. Struct. Health Monit. 2011, 9, 467–480. [Google Scholar] [CrossRef]
  32. Yang, H.J.; Lee, Y.J.; Lee, S.K. Impact source localization in plate utilizing multiple signal classification. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2013, 227, 703–713. [Google Scholar] [CrossRef]
  33. Yang, H.; Shin, T.J.; Lee, S. Source location in plates based on the multiple sensors array method and wavelet analysis. J. Mech. Sci. Technol. 2014, 28, 1–8. [Google Scholar] [CrossRef]
  34. Yuan, S.; Zhong, Y.; Qiu, L.; Wang, Z. Two-dimensional near-field multiple signalclassification algorithm-based impact localization. J. Intell. Mat. Syst. Struct. 2015, 26, 400. [Google Scholar] [CrossRef]
  35. Yuan, S.; Bao, Q.; Qiu, L.; Zhong, Y. A single frequency component-based re-estimated MUSIC algorithm for impact localization on complex composite structures. Smart Mater. Struct. 2015, 24, 105021. [Google Scholar] [CrossRef]
  36. Zhong, Y.; Yuan, S.; Qiu, L. Multiple damage detection on aircraft composite structures using near-field MUSIC algorithm. Sens. Actuat. A Phys. 2014, 214, 234–244. [Google Scholar] [CrossRef]
  37. Golato, A.; Ahmad, F.; Santhanam, S.; Amin, M.G. Multi-helical path exploitation in sparsity-based guided-Wave imaging of defects in pipes. J. Nondestruct. Eval. 2018, 37, 27. [Google Scholar] [CrossRef]
  38. Li, G.; Chattopadhyay, A. Reference-free damage localization in time-space domain for structural health monitoring of X-COR sandwich composites. J. Intell. Mater. Syst. Struct. 2019, 30, 371–385. [Google Scholar] [CrossRef]
  39. Xiao, H.; Shen, Y.; Xiao, L.; Qu, W.; Lu, Y. Damage detection in composite structures with high-damping materials using time reversal method. Nondestruct. Test. Eval. 2018, 33, 329–345. [Google Scholar] [CrossRef]
  40. Song, F.; Huang, G.L.; Hu, G.K. Online guided wave-based detection in honeycomb sandwich structures. Aiaa J. 2012, 50, 284–293. [Google Scholar] [CrossRef] [Green Version]
  41. Bao, Q.; Yuan, S.; Guo, F.; Qiu, L. Transmitter beamforming and weighted image fusion–based multiple signal classification algorithm for corrosion monitoring. Struct. Health Monit. 2019, 18, 621–634. [Google Scholar] [CrossRef]
  42. Bao, Q.; Yuan, S.; Wang, Y.; Qiu, L. Anisotropy compensated MUSIC algorithm based composite structure damage imaging method. Compos. Struct. 2019, 214, 293–303. [Google Scholar] [CrossRef]
  43. Qiu, L.; Yuan, S. On development of a multi-channel PZT array scanning system and its evaluating application on UAV wing box. Sens. Actuators A Phys. 2009, 151, 220–230. [Google Scholar] [CrossRef]
Figure 1. Elastic wave propagation model in the near field.
Figure 1. Elastic wave propagation model in the near field.
Sensors 20 01265 g001
Figure 2. Transmitter beamforming based on the sector piezoelectric array.
Figure 2. Transmitter beamforming based on the sector piezoelectric array.
Sensors 20 01265 g002
Figure 3. Basic flow of the SPATB-MUSIC damage monitoring method.
Figure 3. Basic flow of the SPATB-MUSIC damage monitoring method.
Sensors 20 01265 g003
Figure 4. Experimental system and the circular bonded structure with a smart sensor layer.
Figure 4. Experimental system and the circular bonded structure with a smart sensor layer.
Sensors 20 01265 g004
Figure 5. Dimensions of the circular bonded structures and sector sensor arrays.
Figure 5. Dimensions of the circular bonded structures and sector sensor arrays.
Sensors 20 01265 g005
Figure 6. Three typical channels.
Figure 6. Three typical channels.
Sensors 20 01265 g006
Figure 7. Signal waveforms of typical channels. (a) The signal of the PZT A1-PZT S18 channel. (b) The signal of the PZT A1-PZT S15 channel. (c) The signal of the PZT A1-PZT S12 channel.
Figure 7. Signal waveforms of typical channels. (a) The signal of the PZT A1-PZT S18 channel. (b) The signal of the PZT A1-PZT S15 channel. (c) The signal of the PZT A1-PZT S12 channel.
Sensors 20 01265 g007
Figure 8. Calculation of flight time using the Shannon complex wavelet transform.
Figure 8. Calculation of flight time using the Shannon complex wavelet transform.
Sensors 20 01265 g008
Figure 9. The location of damage to the structure.
Figure 9. The location of damage to the structure.
Sensors 20 01265 g009
Figure 10. Details of damage. (a) Health. (b) Damage 1. (c) Damage 2. (d) Damage 3. (e) Damage 4.
Figure 10. Details of damage. (a) Health. (b) Damage 1. (c) Damage 2. (d) Damage 3. (e) Damage 4.
Sensors 20 01265 g010
Figure 11. Effects of different sizes of damage on the signals. (a) The size of the damage is 10 × 5 mm2. (b) The size of the damage is 5 × 5 mm2.
Figure 11. Effects of different sizes of damage on the signals. (a) The size of the damage is 10 × 5 mm2. (b) The size of the damage is 5 × 5 mm2.
Sensors 20 01265 g011
Figure 12. Damage scattering array signals before and after damage beamforming with a damage size of 10 × 5 mm2. (a) Before beamforming. (b) After beamforming.
Figure 12. Damage scattering array signals before and after damage beamforming with a damage size of 10 × 5 mm2. (a) Before beamforming. (b) After beamforming.
Sensors 20 01265 g012
Figure 13. Damage scattering array signals before and after beamforming with a damage size of 5 × 5 mm2. (a) Before beamforming. (b) After beamforming.
Figure 13. Damage scattering array signals before and after beamforming with a damage size of 5 × 5 mm2. (a) Before beamforming. (b) After beamforming.
Sensors 20 01265 g013
Figure 14. Damage localization results using the single transmitter MUSIC algorithm. (a) D1, (b) D2, (c) D3, (d) D4.
Figure 14. Damage localization results using the single transmitter MUSIC algorithm. (a) D1, (b) D2, (c) D3, (d) D4.
Sensors 20 01265 g014
Figure 15. Damage localization results using the SPATB-MUSIC algorithm. (a) D1, (b) D2, (c) D3, (d) D4.
Figure 15. Damage localization results using the SPATB-MUSIC algorithm. (a) D1, (b) D2, (c) D3, (d) D4.
Sensors 20 01265 g015
Figure 16. Results of damage localization for a damage size of 5 × 5 mm2. (a) Path imaging localization result. (b) Delay-sum imaging localization result. (c) STAPB-MUSIC localization result.
Figure 16. Results of damage localization for a damage size of 5 × 5 mm2. (a) Path imaging localization result. (b) Delay-sum imaging localization result. (c) STAPB-MUSIC localization result.
Sensors 20 01265 g016
Table 1. Amplitudes and velocities of the three propagation paths.
Table 1. Amplitudes and velocities of the three propagation paths.
Propagation PathsAmplitudeVelocity
PZT A1-PZT S180.486 V722 m/s
PZT A1-PZT S150.814 V699 m/s
PZT A1-PZT S120.383 V755 m/s
Table 2. Damage information.
Table 2. Damage information.
NumberSizeLocations
D110 × 5 mm2(105°, 4.6 cm)
D25 × 5 mm2(87°, 5.2 cm)
D310 × 5 mm2(95°, 5.1 cm)
D410 × 5 mm2(80°, 5.0 cm)
Table 3. Damage localization results.
Table 3. Damage localization results.
NumberLocationsSingle Transmitter MUSICSPATB-MUSIC
Localization ResultsErrorLocalization ResultsError
AngleDistanceAngleDistance
D1(105°, 4.6 cm)(103°, 3.8 cm)0.8 cm(105°, 4.3 cm)0.3 cm
D2(87°, 5.2 cm)(88°, 4.2 cm)1.0 cm(86°, 4.7 cm)0.5 cm
D3(95°, 5.1 cm)(84°, 5.0 cm)11°0.1 cm(94°, 5.3 cm)0.2 cm
D4(80°, 5.0 cm)(86°, 4.9 cm)0.1 cm(80°, 4.8 cm)0.2 cm
Table 4. Localization results and errors for a damage size of 5 × 5 mm2.
Table 4. Localization results and errors for a damage size of 5 × 5 mm2.
Imaging MethodsPath ImagingDelay-Sum ImagingSPATB-MUSIC
Errord1.5 cm1.1 cm0.5 cm
Errorσ18.8%13.8%6.3%

Share and Cite

MDPI and ACS Style

Fu, T.; Wang, Y.; Qiu, L.; Tian, X. Sector Piezoelectric Sensor Array Transmitter Beamforming MUSIC Algorithm Based Structure Damage Imaging Method. Sensors 2020, 20, 1265. https://doi.org/10.3390/s20051265

AMA Style

Fu T, Wang Y, Qiu L, Tian X. Sector Piezoelectric Sensor Array Transmitter Beamforming MUSIC Algorithm Based Structure Damage Imaging Method. Sensors. 2020; 20(5):1265. https://doi.org/10.3390/s20051265

Chicago/Turabian Style

Fu, Tao, Yi Wang, Lei Qiu, and Xin Tian. 2020. "Sector Piezoelectric Sensor Array Transmitter Beamforming MUSIC Algorithm Based Structure Damage Imaging Method" Sensors 20, no. 5: 1265. https://doi.org/10.3390/s20051265

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop