1. Introduction
1.1. State of the Art
The structural health monitoring (SHM) technology is increasingly used for the detection of progressive defects in the critical aerospace structures [
1,
2,
3,
4]. SHM has been introduced for both active and passive damage identification for the real-time monitoring of the structures. Monitoring of the acoustic emission (AE) from the progressive fatigue damage is categorized as the passive online monitoring [
5,
6,
7]. The AE signals from the fatigue crack have always been an interest for the researchers [
8,
9,
10,
11]. Lee et al. experimentally showed that the AE waves from fatigue crack growth propagate as guided waves [
8]. They used the commercially available resonant type (250 kHz) and wideband (100 kHz–600 kHz) AE sensors to capture the AE signals. These sensors were designed to effectively measure the out-of-plane wave motion. However, the AE signals from the fatigue crack are usually low-amplitude signals and therefore challenging to detect using the conventional AE sensors [
12,
13].
Piezoelectric AE sensors such as PICO, S9225, WSA had been used to detect the high-amplitude simulated AE signals [
14,
15]. They were also used for detecting the AE signals from the pencil lead break (PLB), impact damage [
16,
17]. However, the AE signals from these sources were generally high-amplitude signals. They had higher-amplitude out-of-plane components and these AE sensors were also well-constructed to sense the out-of-plane wave motion. Thus, the AE sensors could effectively sense the high-amplitude AE signals from these sources. However, the AE signals due to the fatigue crack in the thin plate-like aerospace structures are much more complex than the simple PLB and impact damage [
18,
19,
20,
21].
AE hits from the crack in a notched thick specimen were captured by the piezoelectric AE sensors [
12]. However, the surface waved were dominant in these AE signals. The statistical Bayesian approach was developed to analyze the AE hits by Agletdinov et al. [
22]. Nam and Mal [
23] identified three types of elastic during fatigue crack in the aluminum specimen. However, no load information was mentioned about at what load level the AE signals had happened. Andreykiv et al. analytically showed that the formation of a penny-shaped crack in the aluminum alloys emits acoustic waves [
24]. They simulated different modes of crack formation by a sudden drop of the stress on the crack surface.
The experimental analyses supported by numerical simulations provide deeper insight of the elastic wave propagation phenomena [
25,
26,
27]. Willberg et al. [
28] reviewed the research on guided wave simulations and compared various numerical methods. They provided a comprehensive list of references for the simulation of guided waves. They showed that depending on the application certain methods could be more appropriate than others. The higher order finite element method (FEM) for Lamb wave had been performed [
29]. Gresil et al. [
30] investigated the general guidelines for Lamb wave modeling as referred in the review paper as well. These research provided very important guidelines that have been followed in the present research. However, the fatigue crack related AE source modeling had not been explicitly investigated in these studies.
Many theoretical works were developed to correlate the AE waveforms to their sources [
31,
32,
33]. In seismology, the moment tensor approach was used to describe the movement of a fault during an earthquake [
34,
35,
36]. The moment tensor comprises with nine generalized couples including three dipoles. Buried monopole and dipole with various step functions and fracture mechanics-based methods were used to model the AE source [
37,
38,
39]. However, the interaction of the AE waveform to the crack boundary in a thin plate is yet to be analyzed. The finite element (FE) methods are the most effective tool to estimate the acoustic wave field [
28,
40].
Piezoelectric wafer active sensors (PWAS) were commonly used for SHM applications [
41,
42]. These sensors successfully captured the guided waves in active detection technique [
43,
44]. The active detection method required an user-defined excitation signal where the high-amplitude signals can be used to excite the structure.
1.2. Principles of the PWAS Transducer
PWAS transducers are small, lightweight, inexpensive, and can be produced in different geometries [
45]. They can be permanently bonded on host structures in large quantities and achieve real-time monitoring of the structural health status. They can be conveniently used for transmitting and receiving guided waves. A PWAS mounted on the structure is shown in
Figure 1a. The sensing modes of a PWAS are illustrated in
Figure 1b. It can measure both symmetric (wave motions are predominantly in-plane) and antisymmetric Lamb wave modes (wave motions are predominantly out-of-plane).
PWAS transducer couples the electrical and mechanical effects (mechanical strain,
, mechanical stress,
, electrical field,
, and electrical displacement,
). The piezoelectric constitutive equations in tensor notations can be written as
where
is the mechanical compliance of the material measured at zero electric field (
),
is the dielectric permittivity measured at zero mechanical stress (
), and
represents the piezoelectric coupling effect. PWAS utilizes coupling between in-plane strains
and transverse electric field
. The
coupling is the main effect but other couplings maybe present due to manufacturing limitations. However, the influences of the other couplings on the voltage signals are relatively smaller.
The AE signals have much lower amplitude than the typical guided waves generated by the PWAS transducer for active sensing applications. The suit of the PWAS transducers in application to low-amplitude fatigue-crack related AE signals has yet to be done. Assessment of the PWAS transducer to capture the AE signals using multiphysics FEM simulation and experiment would be an interest in the SHM and NDE community.
1.3. Scope of the Article
The novelty of the article consists in using PWAS transducers to capture AE signals generated by an advancing crack during cyclic fatigue loading. In the past, PWAS transducers had been used for active SHM sensing where a PWAS transmitter is excited with a relatively large electric tone burst (10–20 V amplitude) and the resulting elastic waves are captured by a PWAS receiver. In the present work, we have used only PWAS receivers during passive SHM sensing to capture the elastic waves generated by AE events during fatigue cyclic loading of a cracked specimen. The elastic waves generated by AE events are of relatively much lower amplitude than those generated during active SHM sensing. We found that the AE signals are much weaker than the active SHM signals; the AE signals were also contaminated with environmental and electromagnetic noise. We had to implement denoising techniques to clean the AE signals and reveal their frequency content related to fatigue cracking as contained in the AE waveform [
46]. A combination of multiphysics simulations and experiments were performed. Some important aspects of the AE signals have been obtained that could potentially be used for the fatigue crack length estimation. To validate the simulation results, in-situ AE-fatigue experiments were performed. The experimental results supported the simulation results.
2. Description of the 3D Multiphysics FEM Simulation
The multiphysics simulations were performed on the three-dimensional finite element (FE) model by combining the physics of piezoelectric and structural materials. The wafer type piezoelectric sensors also known as PWAS were simultaneously modeled with the host structure to accomplish the multiphysics simulation. The transient FE analysis was performed to obtain the time-domain signals and the acoustic wave propagation in the structure. Aircraft grade aluminum 2024-T3 material properties were used for the host structure. The material properties of the plate were: modulus of elasticity of the material, E = 73.1 GPa; Poisson’s ratio, ν = 0.33, the density, ρ = 2780 kg/m3; and the shear wave speed, cs = 3140 m/s.
2.1. Practical Aspects of the Model
The geometry of the model is shown in
Figure 2. A thin (thickness = 1 mm) square plate with 120 mm side was modeled using ANSYS 15.0 (ANSYS, Inc., Canonsburg, PA, USA). Since it has been an interest for the fatigue crack related AE signal simulation, a fatigue crack was introduced in the plate. In practice, the fatigue crack always initiates from some kind of discontinuities such as holes, notch, joints etc. By keeping this in mind, a small 1-mm diameter hole was modeled in the center of the crack. The hole was modeled with a 9.5-mm fatigue crack on each side. The total length (2a) of the crack was 20-mm (tip to tip). The host structure was meshed with eight nodes SOLID45 structural element. This geometry was chosen by keeping in mind that we could run an AE-fatigue experiment to validate the simulation results.
2.2. Modeling of the Fatigue Crack
There were numerous ways to model an actual fatigue crack in FEM [
47,
48]. We adopted a popular and simpler method to model the crack. The crack was modeled using the discontinuity at the adjacent pair of nodes along the two crack surfaces. There were two sets of nodes along the crack faces and each set was representing the nodes on each face. The nodes are discontinuous along the crack faces and there was no interaction between these nodes. Two sets of nodes were adjacent to each other and the solid elements were disbanded along the crack faces. This type of modeling can represent a crack opening situation while crack grows near the peak load of the fatigue loading in practice. The modeling of the cracks using the above approach is fair enough to model the actual cracks with a small hole in the plate-like structure in the finite element methods [
47,
49]. This crack modeling approach had also been implemented in active SHM sensing [
50]. The symmetric crack was modeled about the center small hole as shown in
Figure 2a.
2.3. Modeling of PWAS Transducers
Two PWAS transducers were modeled together with the host aluminum structure. The size of each PWAS was 7-mm. One PWAS was modeled in the near-field of crack and the center of the PWAS was at 5-mm from the crack. The other PWAS was modeled in the far-field of crack and the center of the PWAS was at 25-mm from the crack. The PWAS was meshed with eight node SOLID5 coupled field element. This element with properly selected coupled field options can handle the physics of piezoelectric and structural materials. The properties of the PWAS were chosen based on the experimentally used PWAS transducers supplied from the manufacturer. The in-plane polarization was chosen for the modeled PWAS so that it could sense the in-plane strain of the host structure. The following material properties for the PWAS transducer had been used [
51].
where
is the elastic stiffness matrix,
is the dielectric matrix, and
is the piezoelectric stress matrix.
The piezoelectric transduction used through the PWAS material properties directly gave the voltage signal as an output. Any mismatch due to meshing should be avoided for better FEM results. Thus the PWAS and the underneath host structure was meshed using the same topology of meshing. This ensured that each node of the PWAS element coincided with the corresponding node of the host structure element. The nodes between the PWAS elements and the structural elements are then merged together to simulate a perfect bonding between PWAS and host structure. However, in practice, there is a very thin adhesive layer in between the PWAS and host structure. The thickness of the adhesive layer depends on the operator and this variability is not considered in the present work. The effect of adhesive layer had been investigated by Santoni-Bottai [
52], Lin et al. [
53] in our LAMSS group. It is recommended to use a very thin layer of adhesive of better agreement between simulation and experimental signal.
2.4. Modeling of the Acoustic Emission Event at the Crack Tip
The acoustic emission source at the crack tip was modeled using the dipole concept as suggested by Hamstad et al. [
37] and Prosser et al. [
26]. The dipoles basically represent the components of the moment tensor as discussed by Aki and Richards [
34]. The moment tensor components were used to simulate the seismological acoustic events. Two dipoles were used to simulate the AE event due to fatigue crack growth. The dipoles were applied at the tip of the crack as shown in
Figure 2b. The dipoles were placed at one element apart from each other that represented a unit growth of the fatigue-crack. Since the crack in the thin plate usually happens along the entire thickness, the dipoles were placed along the thickness as line dipole source as shown in
Figure 2c.
Each dipole consisted of two self-equilibrating point forces of equal strengths. The magnitude of the point force follows a sharp rising step function. This type of step function represented a wideband AE source at the crack tip. The rise time was the characteristic parameter of the step function. The half-cosine bell rising step function with 1.5 μs rise time was used in the modeling following the Ref. [
37]. The time variation and the frequency spectrum of the half-cosine bell step function are shown in
Figure 3. This type of the smoothing step up function was used for the modeling since it provided a better simulation of the AE event as mentioned in Ref. [
37]. The transient integration parameters were: γ = 0.0, α = 0.25, δ = 0.5, θ = 0.5, oscillation limit criterion = 0.5, time step = 0.12 µs. These values were carefully chosen to obtain numerical stability.
The element size of the model was selected such that the simulation result could predict a wide range of frequencies (40 kHz to 700 kHz). In general, 10–15 nodes per wavelength were considered throughout the model. The element sizes were varied in the model. The smaller elements were used in the high-stress gradient region especially near the crack, interface of the sensors and structure. The relatively larger element size was used outside the crack region. This resulted in more than one million (1.092 M) elements in the 3D space. The multiple FEM simulations were performed (not shown here) to obtain the convergence of the results. Four elements were used in the thickness direction. The meshing was used following the criteria mentioned in Ref. [
37] to achieve better FEM results. The non-reflective boundary (NRB) was used in the FEM model to avoid the reflections from the plate edges. The NRB was modeled using the criteria mentioned in Ref. [
54].
The aspect ratio of the element was maintained close to one which allowed a better accuracy in the results. The full integration method was used in the solution to avoid the hourglass effect as recommended by ANSYS. Also, with the proper refinement of meshing and maintaining reasonable aspect ratio, the shear locking effect has been diminished in the FEM results, but, may still be present in smaller amounts, since the full integration method was used. To completely avoid these artefacts, higher order elements or different finite element schemes as presented in Ref. [
29] are suggested for further investigation.
3. Multiphysics FEM Simulation Results
The multiphysics simulation results showed that the acoustic event at the crack tip generated acoustic waves in the structures. These acoustic waves were captured by the both PWAS transducers as AE signals. The simulated results can be discussed in two significant aspects as follows:
3.1. Acoustic Emission Propagation and Secondary Source of AE Signals
The effect of the AE event at the crack tip propagated as the guided acoustic waves in the plate. The animation snapshots at various time of the simulation are shown in
Figure 4. An AE event occurred at the one crack-tip, referred to as “primary AE event” and the generated AE waves travel in the plate. The acoustic waves hit the near-field PWAS 1 and then passed it as the time had progressed. Subsequently, the guided waves hit the far-field PWAS 2 and the other crack-tip. When the waves hit the other crack-tip it generated another AE event referred to as “secondary AE event”.
Both primary and secondary AE event may generate the AE signals that were captured by the two PWAS transducers. The local dynamics of the crack surface and the acoustic waves may generate a crack resonance as we verified by the experimental measurement in our recently published paper [
55]. This local crack resonance phenomenon is highly related to the fatigue crack length. Thus it could potentially be used for fatigue crack length estimation.
It was observed that the acoustic waves propagate with Rayleigh surface wavespeed along the crack surface. The near-field PWAS may sense complex guided wave signals as it was placed near to the source. The far-field PWAS may sense the plate guided waves and some feature of the secondary AE event.
3.2. Simulated AE Signals Captured by Two PWAS Transducers
From the simulation, it was shown that both PWAS transducers captured the AE signals. The time domain signals of the two PWAS transducers are shown in
Figure 5a,b. The fast Fourier transform (FFT) was performed to show the frequency spectra of the two signals.
Figure 5c shows that the highest peak of the PWAS 1 signal is at 50-kHz frequency. The next consecutive peaks are located at 100, 200, 320-kHz.
Figure 5d shows that the highest peak of the PWAS 2 signal is at 200-kHz frequency. The next consecutive peaks are located at 50, 100, 320, 450-kHz. The higher frequency contents such as 100, 200, 320, 450-kHz peaks gained amplitude as compared to the 50-kHz frequency peak. It may be because these high-frequency signals may become fully-developed as they travel away from the source.
3.3. Simulated AE Signals with Two PWAS Transducers on the Two Sides of the Crack
A new FEM model has been set up to study any dynamic characteristics, mode conversion or shadowing effect on the far PWAS due to the obstruction of the near PWAS. In this case, the far PWAS was kept same as the previous case but the near PWAS was modeled on the other side of the crack. The relative position of the two PWAS transducers with respect to the fatigue crack is illustrated in
Figure 6. In this new configuration of the sensors, it was expected that the far-field PWAS signal would not be affected by the near-field PWAS.
The animation snapshot of the acoustic wave propagation at different time of the simulation is illustrated in
Figure 6. The wave propagation phenomenon is almost same as the previous situation (
Figure 4). The three types of guided waves were observed: crack surface waves, the plate guided Lamb waves and SH waves. The crack resonance feature was also observed in this case. The simulated AE signals were recorded by the two PWAS transducers.
3.4. Comparison between the Two Simulated AE Signals Corresponding to the Two Sensor Configurations
Two sensor configurations were considered. Configuration 1: the two sensors were on the same side of the crack as illustrated in
Figure 2a; configuration 2: the two sensors were placed on two sides of the crack (
Figure 6). These configurations are referred to as “1” and “2”. The AE waveforms and their frequency spectra of the near-field PWAS signals are shown in
Figure 7. It indicates that the both signals are almost same as expected. Hence, the simulated AE signals were consistent with each other.
The AE waveforms and their frequency spectra of the far-field PWAS signals are shown in
Figure 8. The waveforms and frequency spectra were very similar to each other. Hence, we could conclude that the near-field PWAS had very minimal effect on the dynamic characteristic of the far-field PWAS.
3.5. Simulated AE Signals with Individual PWAS Transducer
We also conducted FEM simulations by placing one sensor at a time (not shown here). The similar results as discussed in
Section 3.3 and
Section 3.4 were obtained. As far as the two sensors were sufficiently far from each other, the one sensor dynamics has no effect on the dynamics of another sensor.
4. Description of the In-Situ AE-Fatigue Experimental Setup
To validate the multiphysics simulation results, an in-situ AE-fatigue experiment was designed using two PWAS transducers. The AE signals emanated from the fatigue crack was measured with simultaneous measurement of the fatigue loading. Aircraft grade aluminum Al-2024 T3 material was used to make the specimen. The dimension of the specimen was 305 mm in length, 100 mm in width, and 1 mm in thickness. Like the non-reflective boundary in the simulation, the wave-absorbing clay boundary was used in the experiment. The schematic diagram of the experimental setup is shown
Figure 9a. Two PWAS transducers were bonded at 5-mm and 25-mm from the crack (
Figure 9b) by mimicking the simulation configuration.
For the fatigue crack initiation, a small 1-mm hole is drilled at the center of the specimen and the cyclic fatigue loading was applied by using the MTS machine (MTS Systems Corporation, Eden Prairie, MN, USA). An axial tensile cyclic fatigue loading was varied sinusoidally between 2.3 kN to 23 kN. These load levels gave a stress level of 6.5% and 65% of the yield limit (345 MPa) of the material which is commonly used for practical aircraft testing for structural integrity [
56]. A 20-mm fatigue crack was created to mimic the simulation condition. During the initial crack experiment, no AE instrumentation was employed as shown in
Figure 10a.
The specimen with 20-mm fatigue crack was then equipped with two PWAS transducers and wave absorbing clay boundary. The diameter of each PWAS transducer was 7 mm. The instrumented specimen was then subjected to fatigue loading as shown in
Figure 10b. The use of wave absorbing boundary would provide cleaner AE signals without having any plate edge reflections.
Three measurement systems were used simultaneously during the latter stage of the experiment. They were: (a) AE signal measurement by the AE system; (b) MTS load application by the MTS system; and (c) fatigue crack growth measurement by the high-resolution video recording. The load level was reduced to the 60% of the previous load level and the cyclic loading was applied slowly (0.05 Hz). This loading frequency had no interference with the captured AE signals. Also, a band-pass filter (30 kHz–700 kHz) was used to avoid any interference from the low-frequency noises such as hydraulic loading, MTS grips, and mechanical vibrations. A 40 dB preamplifier was also used with the band-pass filter which was recommended by the manufacturer of the AE system. Under the axial cyclic loading, the fatigue crack grew from 20-mm to 25-mm. For the triggering of the AE signal measurement, a threshold had been chosen 2 dB above the environmental noise level.
5. AE-Fatigue Experimental Results
From the experimental results, we observed that the PWAS transducers successfully captured the AE signals coming from the fatigue crack. The experimental results can be discussed in two aspects as follows:
5.1. Fatigue Crack Related AE Hit Captured by the Two PWAS Transducers
The fatigue crack released multiple AE hits in each cycle as the crack grew. The AE hits captured by the PWAS transducers are plotted with the cyclic loadings on
Figure 11. This shows that AE hits happened at almost every cycle of the fatigue loading. Most of the AE hits happened near the peak load. Because, when the load reaches to the peak load, the stress level at the tip of the crack reaches a critical level causing the failure of material-bonding and may cause the crack growth. A portion of the total cycles was shown in
Figure 11 for clear illustration of the AE hits per loading cycle. It was observed that as the fatigue cycle continued more AE hits per cycle were captured by the PWAS transducers (not shown here).
Figure 11 also shows that near-field PWAS captured more AE-hits than the far-field PWAS, as expected. The both sensors may capture the propagating AE signals. The near-field PWAS may capture the evanescent AE signals that would diminish as they travel further distance and could not meet the threshold of AE hit detection. The propagating AE signals showed geometric spreading and structural attenuation as they traveled far away from the crack.
Figure 12 shows the zoomed-in view of two fatigue cycles (100 s to 145 s) of
Figure 11. This clearly shows the one-to-one correspondence between the AE hits captured by the two PWAS transducers. In the first cycle of
Figure 12, the two AE-hits happened during loading period and both of them are captured by the two sensors. The similar observation can be made for the unloading period. However, there were two weaker AE hits that were captured by only near-field PWAS transducer. The similar observation holds true for the second cycle of
Figure 12. Note that each AE hit corresponds to an AE signal.
5.2. Fatigue Crack Generated AE Signals Captured by the Two PWAS Transducers
The AE signals from the fatigue crack were captured by the two PWAS transducers. The time domain signals and their frequency spectra are illustrated in
Figure 13. The same AE event captured by the two sensors is shown here. The same AE event was confirmed by the one-to-one correspondence of the AE-hits as depicted earlier in
Figure 12. The amplitude of the near-field PWAS signal was higher than that of the captured the far-field PWAS signal, as expected. This is because of the geometric spreading of the acoustic wave. The frequency spectra of the two signals (
Figure 13) show some considerable differences.
The highest peak of PWAS 1 signal is located at 50-kHz frequency. There are also other peaks at 100-kHz and 170-kHz frequencies. On the PWAS 2 frequency spectrum, the highest peak is located at 100-kHz. There are also some other peaks at 50, 170, 220, 320-kHz frequencies.
5.3. The Effect of AE Sensor Location on the AE Signals
The AE sensor position has a significant effect on the captured AE signals. The geometric spreading of the acoustic wave signal has been observed. Overall, as the distance increases from the source, the AE signals become weaker. Sometimes, the AE signal may not be detectable if the sensor is too far from the source and the hit-amplitude falls under the threshold. Since the near-field PWAS was placed very close to the fatigue crack source, the local dynamics of the crack surface may affect the near-field AE signals. In the simulation, the crack resonance (local dynamics of the crack surface) phenomenon resulted in a low-frequency peak (50-kHz) in the AE signal of the near-field PWAS. In the experimental AE signal, the near-field PWAS may capture this local crack resonance phenomenon that resulted in a similar low-frequency peak (50-kHz). Because of the dominant low-frequency peak, the higher frequency peaks may look smaller in the plot. Thus the multiphysics simulation helped us to understand the experimental AE signals.
6. Comparison between the Multiphysics Simulation and Experimental Results
The multiphysics simulation and the experimental study suggested that the PWAS transducer were capable of the sensing the AE signals. The experimental results show the agreement with the simulation results in the following aspects
6.1. The Lower Amplitude Signals as the AE Signals Travel Far Away from the Crack
The multiphysics simulation showed that the amplitude of the far-field PWAS AE signal was lower than the near-field AE signal. The same observation was found true for the experimentally measured AE signals. This is intuitive and expected for any guided wave generated from a point source or a thickness-wise distributed line-source in a plate (this is also true for bulk waves). Thus we may conclude that the AE signals travel as guided waves in the structure. This finding was well-supported by the other researchers [
8]. Furthermore, we may construe that the experimentally measured AE signals were generated from the tip of the fatigue crack.
6.2. The Effect of the Local Dynamics of the Crack on the AE Signals
The first two peaks of the AE signals were observed near 50-kHz and 100-kHz. Both multiphysics simulation and the experiment had supported these results. Note that the 50-kHz peak is dominant in the near-field sensor. Since the local crack dynamics was observed from the simulation (
Figure 14), the similar phenomenon may happen during the experiment. This local dynamics may result in the 50-kHz peak in the experimental AE signals.
Also, this 50-kHz peak was becoming weaker as the AE signal travel far away from the crack but did not disappear as it reached to the far-field sensor. The effect of local dynamics (especially the secondary source from the other crack tip as observed from the simulation,
Figure 14b) may propagate to the far-field sensor. However, it became weaker as observed from both experiment and simulation. On the other hand, the higher frequencies such as 170, 220, 320-kHz may be related to the plate guided waves. These frequencies may be related to the primary effect of the fatigue AE event at the crack tip.
6.3. Discussion of the Differences in Simulation and Experimental Results
In addition to the similarities between the simulation and experimental AE signal we have also observed some differences. The simulated near-field PWAS signal showed some frequency peaks at 200, 320-kHz. These frequency peaks were not so dominant in the experimental near-field PWAS signal. Also, on the far-field AE signal, the 450-kHz peak was observed in the simulation but not in the experimental signal. The several reasons could be responsible. For example, the exact features of the real fatigue crack were not the same as the simulation. These features include the zigzag nature of the actual fatigue crack propagation, through-thickness waviness, the slope of the fracture path that dictates the direction of the dipole force, and the actual amount of crack propagation in each AE event. Also, one component of moment tensor may not be enough to simulate the fatigue AE event as assumed in the simulation which would be considered for future investigation.
Despite the differences discussed above, the multiphysics simulation results were in very good agreement with the experimental results. The main objective that was to assess the capability of PWAS transducer to capture the low-amplitude fatigue crack related AE-signals was well-supported by both multiphysics simulation and the experimental measurement.
7. Conclusions
Piezoelectric wafer active sensors (PWAS) successfully captured the fatigue crack related acoustic emission signals. Both multiphysics simulations and experiments supported this result. The advancement of the fatigue crack generates the AE signals at every loading cycle. This indicates that the AE events happen at every cycle as the fatigue crack grows. The simulation results suggested that the dynamics of near-field PWAS had a minimal effect on the far-field PWAS with the prescribed distance mentioned in the paper. The effect of the AE event at the crack tip travels as guided acoustic waves. The AE signals show the geometric spreading of the amplitudes which is similar to the Lamb waves originating from a point source. The distance of the PWAS transducer has a significant effect on the sensing AE signals. The near-field PWAS captures the evanescent AE signals which resulted in higher number of AE hits in the near-field PWAS than that in the far-field PWAS. The local dynamics of the AE signals and the fatigue crack surface has been identified that could potentially be used for the fatigue crack length estimation. This has been supported by both FEM simulation and experimental measurement.
8. Future Work
The non-linear crack analysis would involve self-contact algorithms for crack modeling and could be investigated. The simulated AE signals and the actual experimental fatigue AE signals cannot be exactly matched or matched within some tolerance unless the exact features of the fatigue crack could be modeled. The exact features include the zigzag nature of the actual fatigue crack propagation, through-thickness waviness, the slope of the fracture path that dictates the direction of the dipole force, and the actual amount of crack propagation in each AE event. The AE-fatigue experiment would be performed using the PWAS transducers and the commercial AE sensors. The AE signals from the two types of sensors would be compared to each other. The AE signals corresponding to the larger fatigue crack could be recorded to see if there is any change in the AE signals. The AE hits could be further analyzed to find any possible groups of AE signals. The multiphysics simulation could be performed for other components of the moment tensors. The simulated AE signals could be compared for different moment tensor components.
Acknowledgments
This research work is supported from Office of Naval Research, the U.S. Navy, grant # N000141410655, and we would like to thank the program manager Ignacio Perez. The authors would like to acknowledge the partial funding support from the experiment.com, New York, NY, USA and the Office of the Vice President for Research, University of South Carolina, Columbia, SC, USA.
Author Contributions
Md Yeasin Bhuiyan developed the idea of the multiphysics simulation, performed the simulation, conducted the experiments, analyzed the data and wrote the manuscript. Victor Giurgiutiu advised throughout the processes, give feedback to improve the experiments and revised the manuscript.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Wang, D.; Zhang, W.; Wang, X.; Sun, B. Lamb-wave-based tomographic imaging techniques for hole-edge corrosion monitoring in plate structures. Materials 2016, 9, 916. [Google Scholar] [CrossRef] [PubMed]
- Cai, J.; Yuan, S.; Wang, T. Signal construction-based dispersion compensation of lamb waves considering signal waveform and amplitude spectrum preservation. Materials 2017, 10, 4. [Google Scholar] [CrossRef] [PubMed]
- Ihn, J.-B.; Chang, F.-K. Pitch-catch Active Sensing Methods in Structural Health Monitoring for Aircraft Structures. Struct. Health Monit. Int. J. 2008, 7, 5–19. [Google Scholar] [CrossRef]
- Yuan, S.; Liang, D.; Shi, L.; Zhao, X.; Wu, J.; Li, G.; Qiu, L. Recent Progress on Distributed Structural Health Monitoring Research at NUAA. J. Intell. Mater. Syst. Struct. 2008, 19, 373–386. [Google Scholar] [CrossRef]
- Mo, C.; Wisner, B.; Cabal, M.; Hazeli, K.; Ramesh, K.T.; El Kadiri, H.; Al-Samman, T.; Molodov, K.D.; Molodov, D.A.; Kontsos, A. Acoustic emission of deformation twinning in magnesium. Materials 2016, 9, 662. [Google Scholar] [CrossRef] [PubMed]
- Rekatsinas, C.S. Time Domain Spectral Finite Element Methods for the Simulation of Wave Propagation in Composite Layered Structures with Active Piezoelectric Sensors. Ph.D. Thesis, University of Patras, Patras, Greece, September 2016. [Google Scholar]
- Bhuiyan, M.Y.; Shen, Y.; Giurgiutiu, V. Guided Wave Based Crack Detection in the Rivet Hole Using Global Analytical with Local FEM Approach. Materials 2016, 9, 602. [Google Scholar] [CrossRef] [PubMed]
- Lee, C.K.; Wilcox, P.D.; Drinkwater, B.W.; Scholey, J.J.; Wisnom, M.R.; Friswell, M.I. Acoustic Emission during Fatigue Crack Growth in Aluminium Plates. In Proceedings of the ECNDT, Berlin, Germany, 25–29 September 2006; pp. 1–8. [Google Scholar]
- Lindley, T.C.; Palmer, I.G.; Richards, C.E. Acoustic emission monitoring of fatigue crack growth. Mater. Sci. Eng. 1978, 32, 1–15. [Google Scholar] [CrossRef]
- Scruby, C.B.; Baldwin, G.R.; Stacey, K.A. Characterization of fatigue crack extension by quantitative acoustic emission. Int. J. Fract. 1985, 28, 201–222. [Google Scholar]
- Mitchell, J.R.; Egle, D.M.; Appl, F.J. Detecting Fatigue Cracks with Acoustic Emission. Proc. Okla. Acad. Sci. 1973, 53, 121–126. [Google Scholar]
- Chang, H.; Han, E.H.; Wang, J.Q.; Ke, W. Acoustic Emission Study of Fatigue Crack Closure of Physical Short and Long Cracks for Aluminum Alloy LY12CZ. Int. J. Fatigue 2009, 31, 403–407. [Google Scholar] [CrossRef]
- Abdelrahman, M.; ElBatanouny, M.K.; Ziehl, P.H. Acoustic emission based damage assessment method for prestressed concrete structures: Modified index of damage. Eng. Struct. 2014, 60, 258–264. [Google Scholar] [CrossRef]
- Ernst, R.; Zwimpfer, F.; Dual, J. One sensor acoustic emission localization in plates. Ultrasonics 2016, 64, 139–150. [Google Scholar] [CrossRef] [PubMed]
- Perez, I.M.; Cui, H.; Udd, E. Acoustic emission detection using fiber Bragg gratings. In Proceedings of the SPIE’s 8th Annual International Symposium on Smart Structures and Materials, Newport Beach, CA, USA, 4–8 March 2001; Volume 4328, pp. 209–215. [Google Scholar]
- Gorman, M.R.; Prosser, W.H. Application of Normal Mode Expansion to Acoustic Emission Waves in Finite Plates. J. Appl. Mech. 1996, 63, 555–557. [Google Scholar] [CrossRef]
- Guo, D.; Mal, A.K.; Hamstad, M.A. AE Wavefield Calculations in a Plate. In Proceedings of the International AE Conference, Big Island, HI, USA, 9–14 August 1998; pp. 19–29. [Google Scholar]
- Buttle, D.J.; Scruby, C.B. Acoustic Emission Source Location in Fiber Reinforced Plastic Composites. J. Acoust. Emiss. 1988, 7, 211–223. [Google Scholar]
- Morton, T.M.; Harrington, R.M.; Bjeletich, J.G. Acoustic emissions of fatigue crack growth. Eng. Fract. Mech. 1973, 5, 691–697. [Google Scholar] [CrossRef]
- Scruby, C.B.; Wadley, H.N.G.; Sinclair, J.E. The origin of acoustic emission during deformation of aluminium and an aluminium-magnesium alloy. Philos. Mag. A 1981, 44, 249–274. [Google Scholar] [CrossRef]
- Hamel, F.; Bailon, J.P.; Bassim, M.N. Acoustic emission mechanisms during high-cycle fatigue. Eng. Fract. Mech. 1981, 14, 853–860. [Google Scholar] [CrossRef]
- Agletdinov, E.; Pomponi, E.; Merson, D.; Vinogradov, A. A novel Bayesian approach to acoustic emission data analysis. Ultrasonics 2016, 72, 89–94. [Google Scholar] [CrossRef] [PubMed]
- Nam, K.; Mal, A. Characteristics of elastic waves generated by crack initiation in aluminum alloys under fatigue loading. J. Mater. Res. 2001, 16, 1745–1751. [Google Scholar] [CrossRef]
- Andreykiv, O.; Skalsky, V.; Serhiyenko, O.; Rudavskyy, D. Acoustic Emission Estimation of Crack Formation in Aluminium Alloys. Eng. Fract. Mech. 2010, 77, 759–767. [Google Scholar] [CrossRef]
- Castellano, A.; Foti, P.; Fraddosio, A.; Marzano, S.; Piccioni, M.D. Mechanical characterization of CFRP composites by ultrasonic immersion tests: Experimental and numerical approaches. Compos. Part B Eng. 2014, 66, 299–310. [Google Scholar] [CrossRef]
- Prosser, W.H.; Hamstad, M.A.; Gary, J.; O’Gallagher, A. Finite Element and Plate Theory Modeling of Acoustic Emission Waveforms. J. Nondestruct. Eval. 1999, 18, 83–90. [Google Scholar] [CrossRef]
- Zelenyak, A.; Hamstad, M.; Sause, M. Modeling of Acoustic Emission Signal Propagation in Waveguides. Sensors 2015, 15, 11805–11822. [Google Scholar] [CrossRef] [PubMed]
- Willberg, C.; Duczek, S.; Vivar-Perez, J.M.; Ahmad, Z.A.B. Simulation Methods for Guided Wave-Based Structural Health Monitoring: A Review. Appl. Mech. Rev. 2015, 67, 10803–10820. [Google Scholar] [CrossRef]
- Willberg, C.; Duczek, S.; Vivar Perez, J.M.; Schmicker, D.; Gabbert, U. Comparison of different higher order finite element schemes for the simulation of Lamb waves. Comput. Methods Appl. Mech. Eng. 2012, 241, 246–261. [Google Scholar] [CrossRef]
- Gresil, M.; Giurgiutiu, V.; Shen, Y.; Poddar, B. Guidelines for Using the Finite Element Method for Modeling of Guided Lamb Wave Propagation in SHM Processes. In Proceedings of the 6th European Workshop on Structural Health Monitoring, Dresden, Germany, 3–6 July 2012. [Google Scholar]
- Ohtsu, M.; Ono, K. A generalized theory and source representations of acoustic emission. J. Acoust. Emiss. 1986, 5, 124–133. [Google Scholar]
- Pollock, A.A. Acoustic emission inspections. In Metals Handbook Ninth edition: Volume 17; ASM International: Materials Park, OH, USA, 1989; pp. 278–294. [Google Scholar]
- Prosser, W.H. Applications of advanced, waveform based AE techniques for testing composite materials. Proc. SPIE 1996, 2944, 146–153. [Google Scholar]
- Aki, K.; Richards, P.G. Quantitative Seismology; University Science Books: Sausalito, CA, USA, 2002. [Google Scholar]
- Hudson, J.A.; Pearce, R.G.; Rogers, R.M. Source Type Plot for Inversion of the Moment Tensor. J. Geophys. Res. 1989, 94, 765–774. [Google Scholar] [CrossRef]
- Frohlich, C.; Apperson, K.D. Earthquake focal mechanisms, moment tensors, and the consistency of seismic activity near plate boundaries. Tectonics 1992, 11, 279–296. [Google Scholar] [CrossRef]
- Hamstad, M.A.; O’Gallagher, A.; Gary, J. Modeling of Buried Monopole and Dipole Sources of Acoustic Emission with a Finite Element Technique. J. Acoust. Emiss. 1999, 17, 97–110. [Google Scholar]
- Lysak, M.V. Development of the Theory of Acoustic Emission by Propagating Cracks in Terms of Fracture Mechanics. Eng. Fract. Mech. 1996, 55, 443–452. [Google Scholar] [CrossRef]
- Sause, M.G.R.; Richler, S. Finite element modelling of cracks as acoustic emission sources. J. Nondestruct. Eval. 2015, 34, 4. [Google Scholar] [CrossRef]
- Joseph, R.; Bhuiyan, M.Y.; Giurgiutiu, V. Acoustic emission source modeling in a plate using buried moment tensors. In Proceedings of the SPIE (Health Monitoring of Structural and Biological Systems), Portland, OR, USA, 25–29 March 2017; p. 1017028. [Google Scholar]
- Giurgiutiu, V.; Zagrai, A.; Bao, J. Piezoelectric Wafer Embedded Active Sensors for Aging Aircraft Structural Health Monitoring. Struct. Health Monit. 2002, 1, 41–61. [Google Scholar] [CrossRef]
- Qiang, W.; Shenfang, Y. Baseline-free Imaging Method based on New PZT Sensor Arrangements. J. Intell. Mater. Syst. Struct. 2009, 20, 1663–1673. [Google Scholar] [CrossRef]
- Köhler, B.; Gaul, T.; Lieske, U.; Schubert, F. Shear horizontal piezoelectric fiber patch transducers (SH-PFP) for guided elastic wave applications. NDT E Int. 2016, 82, 1–12. [Google Scholar] [CrossRef]
- Yu, L.; Santoni-Bottai, G.; Xu, B.; Liu, W.; Giurgiutiu, V. Piezoelectric wafer active sensors for in situ ultrasonic-guided wave SHM. Fatigue Fract. Eng. Mater. Struct. 2008, 31, 611–628. [Google Scholar] [CrossRef]
- Giurgiutiu, V. Structural Health Monitoring with Piezoelectric Wafer Active Sensors, 2nd ed.; Academic Press: Oxford, UK, 2014. [Google Scholar]
- Kuang, K.S.C.; Li, D.; Koh, C.G. Acoustic emission source location and noise cancellation for crack detection in rail head. Smart Struct. Syst. 2016, 18, 1063–1085. [Google Scholar] [CrossRef]
- Friswell, M.I.; Penny, J.E.T. Crack Modeling for Structural Health Monitoring. Struct. Health Monit. 2002, 1, 139–148. [Google Scholar] [CrossRef]
- Jiang, Y.; Feng, M. Modeling of Fatigue Crack Propagation. J. Eng. Mater. Technol. 2004, 126, 77–86. [Google Scholar] [CrossRef]
- Trendafoilova, I.; Gorman, D.G.; Manoach, E. An Investigation on Vibration-based Damage Detection in Circular Plates. Struct. Health Monit. 2009, 8, 291–302. [Google Scholar] [CrossRef] [Green Version]
- Bhuiyan, M.Y.; Shen, Y.; Giurgiutiu, V. Interaction of Lamb waves with rivet hole cracks from multiple directions. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2017. [Google Scholar] [CrossRef]
- Shen, Y.; Giurgiutiu, V. Predictive modeling of nonlinear wave propagation for structural health monitoring with piezoelectric wafer active sensors. J. Intell. Mater. Syst. Struct. 2014, 25, 506–520. [Google Scholar] [CrossRef]
- Santoni-Bottai, G. Fundamental Studies in the Lamb-Wave Interaction Between Piezoelectric Wafer Active Sensor and Host Structure during Structural Health Monitoring. Ph.D. Thesis, University of South Carolina, Columbia, SC, USA, January 2010. [Google Scholar]
- Lin, B.; Gresil, M.; Cuc, A.; Giurgiutiu, V. Predictive Modeling of Piezoelectric Wafer Active Sensors for Structural Health Monitoring. Ferroelectrics 2014, 470, 168–182. [Google Scholar] [CrossRef]
- Shen, Y.; Giurgiutiu, V. Effective non-reflective boundary for Lamb waves: Theory, finite element implementation, and applications. Wave Motion 2015, 58, 22–41. [Google Scholar] [CrossRef]
- Bhuiyan, M.Y.; Bao, J.; Poddar, B.; Giurgiutiu, V. Toward identifying crack-length-related resonances in acoustic emission waveforms for structural health monitoring applications. Struct. Health Monit. Int. J. 2017, 1–9. [Google Scholar] [CrossRef]
- Bhuiyan, M.Y.; Bao, J.; Poddar, B.; Giurgiutiu, V. Analysis of acoustic emission waveforms from fatigue cracks. In Proceedings of the SPIE (Health Monitoring of Structural and Biological Systems), Portland, OR, USA, 25–29 March 2017; Voulme 10170, pp. 101702A-1–101702A-8. [Google Scholar]
Figure 1.
(
a) Piezoelectric wafer active sensor (PWAS) [
45]; (
b) acoustic waves sensing mode of the PWAS (in-plane and out-of-plane wave motion sensing capability of the PWAS).
Figure 1.
(
a) Piezoelectric wafer active sensor (PWAS) [
45]; (
b) acoustic waves sensing mode of the PWAS (in-plane and out-of-plane wave motion sensing capability of the PWAS).
Figure 2.
The meshed FEM model with in-plane dipole excitation: (a) overall view; (b) top view zoomed into crack tip area showing dipole loading; (c) Side view showing thickness wise assignment of dipole components.
Figure 2.
The meshed FEM model with in-plane dipole excitation: (a) overall view; (b) top view zoomed into crack tip area showing dipole loading; (c) Side view showing thickness wise assignment of dipole components.
Figure 3.
The dipole loading used for simulating the acoustic emission (AE) event at the crack tip. (a) A typical half-cosine bell step function and its frequency response (Rise time = 1.5 µs); (b) The half-cosine bell step-up from 0 to 1 in the amplitude scale.
Figure 3.
The dipole loading used for simulating the acoustic emission (AE) event at the crack tip. (a) A typical half-cosine bell step function and its frequency response (Rise time = 1.5 µs); (b) The half-cosine bell step-up from 0 to 1 in the amplitude scale.
Figure 4.
The animation snapshots of the acoustic wave propagation at different time of simulation: (a) the acoustic wave traveling toward the PWAS 1; (b) the waves reached to the PWAS 2 and the other tip of the crack; (c) the other crack-tip generated secondary AE event; (d) the reflected AE waves traveling back toward the first crack tip; (e) the reflected AE waves hit the PWAS 1 again; (f) after back and forth AE traveling along the crack and another secondary AE event at the other crack-tip.
Figure 4.
The animation snapshots of the acoustic wave propagation at different time of simulation: (a) the acoustic wave traveling toward the PWAS 1; (b) the waves reached to the PWAS 2 and the other tip of the crack; (c) the other crack-tip generated secondary AE event; (d) the reflected AE waves traveling back toward the first crack tip; (e) the reflected AE waves hit the PWAS 1 again; (f) after back and forth AE traveling along the crack and another secondary AE event at the other crack-tip.
Figure 5.
Multiphysics simulated AE signals from the two PWAS transducers. The typical AE signals coming from the same AE event captured by the (a) near-field PWAS 1; (b) far-field PWAS 2; (c,d) Some high-frequency contents (200, 320, 450-kHz) have relatively higher amplitude than the low-frequency component as the acoustic wave traveled away from the crack.
Figure 5.
Multiphysics simulated AE signals from the two PWAS transducers. The typical AE signals coming from the same AE event captured by the (a) near-field PWAS 1; (b) far-field PWAS 2; (c,d) Some high-frequency contents (200, 320, 450-kHz) have relatively higher amplitude than the low-frequency component as the acoustic wave traveled away from the crack.
Figure 6.
The animation snapshots of the acoustic emission signal propagation when the two PWAS transducers were on the two sides of the crack. The near PWAS 1 was at 5-mm from the crack and the far PWAS 2 was at 25-mm from the crack but on the other side of the crack. The traveling fatigue-crack-generated AE signals at different time of simulation are illustrated.
Figure 6.
The animation snapshots of the acoustic emission signal propagation when the two PWAS transducers were on the two sides of the crack. The near PWAS 1 was at 5-mm from the crack and the far PWAS 2 was at 25-mm from the crack but on the other side of the crack. The traveling fatigue-crack-generated AE signals at different time of simulation are illustrated.
Figure 7.
Simulated AE signals captured by the near PWAS 1 for two different sensor configurations: 1 → both sensors are on the same side of the crack, 2 → Sensors are on the opposite side of the crack.
Figure 7.
Simulated AE signals captured by the near PWAS 1 for two different sensor configurations: 1 → both sensors are on the same side of the crack, 2 → Sensors are on the opposite side of the crack.
Figure 8.
Simulated AE signals captured by the near PWAS 1 for two different sensor configurations: 1 → both sensors are on the same side of the crack, 2 → Sensors are on the opposite side of the crack.
Figure 8.
Simulated AE signals captured by the near PWAS 1 for two different sensor configurations: 1 → both sensors are on the same side of the crack, 2 → Sensors are on the opposite side of the crack.
Figure 9.
(a) Schematic diagram of the specimen with fatigue crack and a PWAS transducer bonded at 5-mm and 25-mm from the crack; (b) Actual specimen with two PWAS transducers bonded at 5-mm and 25-mm from the crack.
Figure 9.
(a) Schematic diagram of the specimen with fatigue crack and a PWAS transducer bonded at 5-mm and 25-mm from the crack; (b) Actual specimen with two PWAS transducers bonded at 5-mm and 25-mm from the crack.
Figure 10.
(a) Actual test coupon in the MTS grip without any AE instrumentation for the initial 20-mm crack; (b) Instrumented Test coupon in the MTS grip for capturing fatigue-crack related AE signals. Two PWAS transducers were used to capture the AE signals.
Figure 10.
(a) Actual test coupon in the MTS grip without any AE instrumentation for the initial 20-mm crack; (b) Instrumented Test coupon in the MTS grip for capturing fatigue-crack related AE signals. Two PWAS transducers were used to capture the AE signals.
Figure 11.
The AE-hits captured by the near-field and far-field PWAS transducers. The AE hits are plotted with the cyclic fatigue loading. The AE hits from the crack happened when the fatigue loading reached near the peak load (A portion of the total cycles was shown here for clear illustration of the AE hits per loading cycle).
Figure 11.
The AE-hits captured by the near-field and far-field PWAS transducers. The AE hits are plotted with the cyclic fatigue loading. The AE hits from the crack happened when the fatigue loading reached near the peak load (A portion of the total cycles was shown here for clear illustration of the AE hits per loading cycle).
Figure 12.
Zoomed-in view of the AE hits with the cyclic fatigue loading. The AE hits from the crack happened near the peak load level during both loading and unloading period of the cycle (PWAS 1 and PWAS 2 AE hits were distinguished by blue-square and red-rhombus legend, respectively).
Figure 12.
Zoomed-in view of the AE hits with the cyclic fatigue loading. The AE hits from the crack happened near the peak load level during both loading and unloading period of the cycle (PWAS 1 and PWAS 2 AE hits were distinguished by blue-square and red-rhombus legend, respectively).
Figure 13.
Experimentally measured AE signals by the two PWAS transducers. The typical AE signals coming from the same AE event captured by the (a) near-field PWAS; (b) far-field PWAS. Some high-frequency contents gained amplitude as the acoustic wave had traveled away from the crack.
Figure 13.
Experimentally measured AE signals by the two PWAS transducers. The typical AE signals coming from the same AE event captured by the (a) near-field PWAS; (b) far-field PWAS. Some high-frequency contents gained amplitude as the acoustic wave had traveled away from the crack.
Figure 14.
FEM simulation results of the (a) local dynamics of the crack; (b) reflection of AE waves from the crack tip.
Figure 14.
FEM simulation results of the (a) local dynamics of the crack; (b) reflection of AE waves from the crack tip.
© 2017 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).