1. Introduction
Structural Health Monitoring (SHM) systems [
1], i.e., inspections systems based on permanently installed sensors, have been intensively investigated in recent years, due to the high potential that they have in reducing maintenance costs in a wide range of industrial application fields. In particular, an increasing demand to develop SHM systems comes from the automotive and aerospace industries [
2]. In these fields, open challenges are: (i) the reduction of the monitoring system weight (including the associated cabling and circuitry), the minimization of both (ii) the power consumption, and (iii) the amount of data that are to be transferred and processed
Among the various sensor technologies proposed for SHM, those that are based on the usage of ultrasonic guided waves (GWs), are considered to be particularly attractive because of the GW property to travel across long distances with low attenuation. In particular, passive GW methods can be used to detect and locate acoustic sources (AS) in plates or shells, such as those that are generated by impacts, or acoustic emissions generated by cracks and corrosion, by using piezoelectric sensors [
3,
4].
The approaches presented in literature for AS localization with GW can be divided in three main groups: (i) inverse methods [
5,
6], (ii) hyperbolic location methods [
7], and (iii) methods that use Direction of Arrival (DoA) estimations [
8,
9]. However, if the first two methods are well validated at laboratory scale, these approaches are rarely adopted for field deployment, because of the stringent hypotheses that they rely on. In particular, the first group implies onerous calculations because it requires a very accurate modeling of the guided ultrasonic propagation, whereas the second group requires a large number of sensors in order to reduce the estimation uncertainty. Moreover, this latter approach is unsuitable when the wave propagation characteristics vary in the monitored area. For these reasons, the methods belonging to the third group, which estimate DoAs by using clusters of closely located sensors that do not rely on the knowledge of the characteristics of propagation, are often the only viable strategy to locate AS with embedded systems.
Kundu et al. in [
9,
10] proposed a simple cluster of 3 circular sensors. In this work, this solution will be referred to as “conventional cluster”. The DoA is performed by estimating the Difference in Time of Arrival (DToA) of wave-fronts between sensors. The DToAs, i.e., the time delays, can be estimated with conventional cross-correlation procedures (i.e., the optimal method for measurements that are affected by additive white Gaussian noise
with zero mean [
11]). This cluster typology allows for estimating the DoA with minimum amount of data to acquire (only three signals) and reduced computational cost (only two time delays estimations). Such an approach was proven to be suitable for real-field applications to perform real-time impact detection and localization [
10]. It is worth noting that Park, Packo, Kundu, and Sen in [
12,
13] showed that it is possible to estimate the AS localization even in the case of heavily anisotropic structures (i.e., characterized by elliptical or rhomboidal wave fronts).
In this work, we further investigate the concept of DoA detection via DToAs estimations, by analyzing how beneficial the adoption of clusters constituted by irregularly shaped piezoelectric transducers can be. In particular, in the proposed approach, the sensors are shaped with a procedure based on the Radon Transform (RT) and its Inverse (IRT). This procedure allows to define the DToA between two sensors as an arbitrary DoA function, divided by the wave velocity
v, as better detailed in
Section 3. The RT approach was firstly presented in [
14] in order to estimate the DoA by using two shaped sensors when the actual wave velocity is assumed to be known. In another previous work [
15], the concept has been extended to the case of unknown velocity of propagation, designing a cluster with three shaped sensors. It is worth noting that the relationship between DToA and DoA affects the uncertainty in DoA estimation. In this paper, we define a new cluster of three shaped sensors in order to further reduce the uncertainty with respect to [
15]. The conducted numerical tests show that the DoA estimation performed by the new cluster is more accurate than the one achieved with the conventional cluster or with the cluster proposed in [
15], and it is robust, even when the measurements are affected by noise or acquired at low sampling frequencies. It is also worth noting that the computational cost that is associated to the processing of the shaped cluster signals for DToA estimation is fully equivalent to the one of conventional clusters, while the DoA estimation is even simpler because it purely linear w.r.t the DToAs ratio, while the conventional approach is based on trigonometric functions.
The experimental validation of the shaping concept is beyond the scope of this paper. Nevertheless, it is worth mentioning that the proposed cluster of shaped sensors can be realized by relying on different piezoelectric materials and manufacturing techniques. e.g., metallized PVDF (polyvinylidene fluoride) sheets can be used by shaping the electrodes on the upper surface with a laser cut as in [
16]. Alternatively, the shaping strategy can be based on printing metallic electrodes on PVDF films in order to obtain the desired shape sensors, as proposed in [
17], or by using litographic procedures as in [
18]. In the practical adoption of these devices, a fundamental step is the definition of a well controlled bonding procedure, because it may heavily affect the sensor response [
19].
The paper is organized, as follows. In
Section 2, by exploiting the Propagation of Uncertainty theory [
20], the functions that bind the DoA to the Differences in Time of Arrival (DToAs) of the wavefront at piezo-patches P1 and P2, and at P1 and P3, are defined to reduce the uncertainty in DoA estimation, when the wave mode velocity is not known.
Section 3 shows how to implement the desired functions by properly shaping the sensors in the context of guided waves, in particular of Lamb waves. To this aim, the Radon Transform and its inverse are introduced.
Section 4 validates the theoretically expected performances with numerical simulations. The estimation of the DoA is simulated in challenging working conditions, i.e., with noisy measurements or very low sampling frequencies, to emulate real field acquisitions. Moreover, a numerical comparison between the shaped sensors cluster and the conventional cluster is provided. Finally,
Section 5 depicts the conclusions.
2. Minimization of Uncertainty in Doa Estimation
The simplest sensor configuration to perform the wave DoA estimation is based on a couple of isotropic transducers, i.e., the piezopatches P1 and P2 illustrated in
Figure 1. Assuming a planar wavefront, the Difference in Time of Arrival (DToA,
) is related to the wave DoA (
) by this formula:
where
d is the distance among the sensors and
v is wave velocity. Consequently, the estimation of the DoA can be performed from the measured DToA (
) by applying the inverse formula:
The Propagation of Uncertainty theory [
20] that is applied to this case reveals that the error in the estimation can be very large for low values of
. This is due to the non-linear characteristics of the
arccos relation. To counteract this problem, in [
14], a sensor reshaping procedure was proposed in order to achieve a linear dependence between the DToA and the DoA, so that:
where
and
can be chosen by the designer. Thanks to new linear relation, the worst-case error is constant over the considered range of angles (
), and it is considerably lower than the maximum error that is achieved using the conventional sensors. It is worth noting that the worst case error is inversely proportional to the selected value of
. For this reason, the larger
the lower is the error in the DoA estimation; however, there is a trade off involved in the selection of
, because, the larger is
, the larger are the physical dimensions of the piezoelectric patches and the maximum value of the DToA, which defines the time-window length of the acquired signal (i.e., the information to be stored and processed).
It is also worth noting that a cluster that is composed by two sensors can be used to estimate
just if the actual wave velocity
v is assumed to be known. Therefore, if the velocity is not known a priori, a third sensor P3 is needed in order to estimate the DoA. In the proposed approach, the sensor P3 is shaped, so that it establishes a constant reference to estimate the wave velocity, i.e.,
(where
is a constant), the relation between the DToAs ratio and DoA is:
In the last equation, it can be observed that the DoA
is linked to the ratio of DToA by a simple linear relation. Moreover, in order to estimate the DoA, the knowledge of the wave propagation velocity is not necessary, because such quantity does not appear in (
5). However, the third sensor P3 introduces an additional uncertainty in the DoA estimation due to
.
In particular, the worst case error can be very large when the DoA is around 90 [deg], if
is set to 0. The
Figure 2a shows the upper and lower bounds of worst case error for the shaped sensors cluster proposed in [
15]. It was designed for a work-range of [0, 90] [deg] and it was defined by the following parameter values:
,
,
,
. The value of
d defines the maximum absolute value of DToA, equal to
(in this case, obtained when
[deg], for piezo P2). Therefore, the value of
d is redefined as the Maximum Difference in Distance of Propagation (MDDP). It is worth noting that the MDDP,
d, in general, does not correspond to the physical distance between the shaped sensors. The error bounds were computed and illustrated for all of the clusters here considered by setting the value of the wave velocity
and DToAs Standard Deviation equal to 2 [
s].
The worst case error can be reduced by selecting an appropriate value for
(i.e.,
so that
), or by working on a shifted range of angles, centered on the zero angle, i.e.,
. In this case, the parameters
and
can be used in order to minimize the uncertainty in DoA estimation. In particular,
should be equal to 0, while
and
are defined by selecting the maximum absolute values of
, i.e.,
, equal to the maximum absolute values of DToAs for the conventional cluster. Limiting
implies the limitation of the maximal dimensions of the transducers, and also the relaxation of memory requirements for the storage and subsequent cross-correlation of the acquired time-waveforms. To summarize, the values of the parameters
and the functions
are:
In particular, the new shaped sensors cluster was designed, so that the MDDP is again
.
Figure 2b illustrates the worst-case error for this configuration, in terms of upper and lower bounds.
Figure 3b shows the upper and lower bounds of the conventional cluster with three circular sensors [
9,
10] with the MDDP equal to
(this distance corresponds, for the circular sensors, to the physical distances between sensors in
Figure 3a). The cluster of shaped sensors that implements the relations (
3) and (
4) with the parameters being detailed in Equation (
6) allows for a theoretical average reduction of the worst case error of 11%. It is important to stress the fact that the comparison is provided with the same maximum absolute value for the DToAs (i.e.,
). This means that the two clusters have the same memory requirements and same computational cost to estimate the time delays. It is also worth highlighting that the two clusters differ considerably in terms of DoA estimation function: the conventional cluster requires the evaluation of the
atan function of the DToAs ratio, while the estimation function of the new cluster is simply linear w.r.t. the DToAs ratio. It is given by inverting (
5) with respect to
. Therefore it is compatible with embedded systems that are equipped with micro-controller units.
Unfortunately, implementing a cluster of piezo-patches that are characterized by the DToA-DoA relationships in (
5) and (
6) is not possible with simple piezo-disks. Shaped sensors have to be used, but this introduces an additional complication: in fact, shaped sensors are generally characterized by an angle-dependent wavenumber tuning effect [
19], i.e., the frequency response of the transducer depends on the direction of propagation of the incident wave. In turn, this effect may hamper the possibility to extract the time difference of arrival with conventional cross-correlation procedures, in which the response among different transducers is assumed to be the same. In the next section, the method that was proposed in [
14] to tackle this problem is reviewed and applied to the considered three-sensor cluster case.
3. The Radon Transform as Sensor Shape Design Tool
The goal of having transducers with the desired relations (
3) and (
4), is achieved if the sensors are shaped so that the angle-dependent frequency responses
and
(of P2 and P3, respectively) are linked to
(i.e., the isotropic one of the reference circular sensor P1), so that:
i.e., the frequency response of the three piezo patches should be equal but for a phase shift which is proportional to the
and
quantities which we have introduced in the previous section. The term
is the wave vector-frequency relationship of a guide wave mode, i.e., the so called dispersion curve of the Lamb wave that is captured by the sensor. Replacing
with
where
is the phase wave velocity, under the hypothesis of narrow-band waves, the variations of
can be neglected and, given
the center frequency, a constant velocity
can be considered. The relationships (
7) then create the relationships (
3) and (
4) over time.
To impose the relationships (
7) three key elements can be exploited:
- (a)
The model of the frequency response of a sensor in the presence of a Lamb wave mode.
- (b)
The Radon Transform (RT) and its inverse (IRT).
- (c)
The Projection Slice Theorem.
For the element (a), the model that was proposed by Senesi and Ruzzene [
21] can be adopted:
where
denotes the amplitude and the polarization of the wave component relevant to the piezo properties of the patch at the considered frequency,
is the wave vector that characterizes the propagation,
is a quantity that is related to the material properties of the piezo-structure system, and finally
is the sensor Directivity function, which is the only function that is dependent on the shape of the sensor. It can be computed by the following integral:
where
is referred to as shape function and it describes the shape of sensor as a step function of form:
with
the area of the piezoelectric path. Without lack of generality, here we consider the case of a single piezoelectric patch, with a single polarization and costant piezoelectric properties. From Equation (
9), the Directivity function can be computed as the bidimensional Fourier Transform (FT) at the angle
of the shape function
. Furthermore, the output variable is evaluated in
, the wave vector of the Lamb wave mode considered (e.g.,
,
, etc. modes).
The element (b), which is the Radon Transform [
22], is defined, as follows:
and it consists of multiple line-integrals. It admits inverse.
Finally (Point (c)), the Projection Slice Theorem [
23] states that the bidimensional FT at the angle
of a given function, is equal to the mono-dimensional FT of the Radon Transform of that function. This implies that the Directivity function, and therefore the frequency response, can simply be calculated as the monodimensional FT of the RT of the shape function and by evaluating into
. As a consequence, the relations (
3) and (
4), i.e., (
9) in frequency, with
and
specified in (
6) are obtained by imposing two suitable relations between the RT of sensor P1 and the RT of sensors P2 and P3:
In conclusion, the synthesis procedure of the shape functions for sensors P2 and P3, which realizes the cluster of the desired re-shaped sensors, can be summarized, as follows: consider an isotropic circular reference sensor 1, then calculate its RT (constant over
) and impose the RT of sensors P2 and P3 by using (
12) with
and
given by (
6) (see
Figure 4a,b). In particular, the shaped sensors cluster was designed with
(the Maximum Difference in Distance of Propagation).
By inverting these two RTs, we obtain two shape functions that are continuously modulated (
Figure 4c,d). In order to obtain two step functions, as in (
10), a quantization procedure is imposed by setting all values greater than a certain threshold to 1 and the others to 0. In this way, the desired shape (step) functions can be obtained (as illustrated in
Figure 4e,f).
Figure 4g,h shows the actual RTs (i.e., the RTs computed by the post quantization shape functions). It is worth noting that they differ from the imposed RTs (
Figure 4c,d), due to the binary quantization procedure. To achieve a better matching between the imposed RT and the actual (post-quantization), it can be beneficial to set
constant, even beyond the considered work-range of angles, as shown in
Figure 4b.
Figure 5 illustrates the designed cluster.
4. Numerical Results and Discussion
In order to validate the design procedure and the proposed shaped sensors cluster, as represented in
Figure 5, multiple impacts occurring on an aluminium plate (3 mm thick, Young’s modulus 70 GPa, Poisson’s coefficient
and material density
) were simulated. The response of the shaped piezo-patches to impacts generating the fundamental antisymmetric mode (
) was computed using the Greens function formalism adopted in [
14]. An impulse response of a band pass Batterworh filter (10th order) with narrow bandwidth (i.e., 7.5 [KHz]) was used in order to emulate the actuated pulse. Subsequently, the estimated DToAs between P1 and P2, and the one between the P1 and P3 were computed by the locating the peak of the cross-correlation of acquired signals (i.e., the Maximum Likelihood (ML) estimator of DToA for measurements affected by additive white Gaussian noise (AWGN) with zero mean [
11]), and their ratio
was used to estimate the DoA (
).
The ideal estimation function, linearly dependent with respect to the ratio
, could be obtained by inverting (
5) with respect to
. However, due to binary quantization procedure, the actual values of parameters
and
might be slightly different with respect to the imposed ones (
6). In
Figure 6, it can be seen that, when a narrow-band Lamb wave mode (e.g.,
) is simulated, the relationship between the DToA ratio and the DoA (
) is still linear. The parameters of the best-fit line, i.e., y-intercept
m and slope
q, can be computed by applying conventional linear regression procedures.
The estimated DoA (
) as a function of its actual value (i.e.,
) for different working conditions is shown in
Figure 7 and
Figure 8, when the
mode propagation is simulated with r = 0.4 [m] as impacts distance. In particular, the two figures report the results at different values of sampling frequency (Fs) and at different values of Peak Signal to Noise Ratio (PSNR) when the measurements are affected by additive white Gaussian noise (AWGN) with zero mean, respectively. As can be seen, the DoA estimation is very accurate, with a standard deviation (SD) of DoA below three degrees, even when the measurements are affected by Gaussian noise (up to the PSNR = 18 dB) and the simulated "physical" sampling frequency is low (down to Fs = 100 KHz).
It is worth noting that:
In order to generate the
Figure 7, the acquired signals were resampled (upsampling factor equal to 10) to improve the estimation accuracy of DToAs and so of DoA. The sampling rate conversion is the optimal method to reduce the uncertainty on DToAs (i.e., using
functions to interpolate), but the computational cost can be significant. In embedded systems, such onerous processing may be unfeasible. Possible alternatives consist in: (i) transmitting the acquired signals to a central processing unit (CPU) or (ii) using different interpolation techniques with reduced computational complexity, such as those proposed in [
24,
25].
The upper- and lower- bound errors that are displayed in
Figure 7 are calculated by applying the propagation of uncertainty theory [
26].
In the
Figure 8, to numerically assess the SD of DoA, 1000 simulations were performed with measurements being affected by AWGN with zero mean, for each value of PSNR.
A numerical comparison of the performance in terms of SD between the conventional cluster of three circular sensors and the designed cluster is given in
Table 1, assuming that measurements are affected by AWGN with zero mean. The simulation of the A0 mode propagation was run when considering r = 0.4 [m] as impacts distance and Fs = 1 [MHz]. In order to assess the SD of DoA numerically, 1000 simulations were performed. The comparison is provided with the same maximum absolute value for the DToAs (i.e.,
), so with the same memory requirement and same computational cost to estimate the time delays. As can be seen from
Table 1, the SD value that is obtained with the cluster of shaped sensors is considerably reduced for all PSNR values.