In Guided Wave Structural Health Monitoring (GW-SHM), a strong need for reliable and fast simulation tools has been expressed throughout the literature in order to optimize SHM systems or demonstrate performance. Even though guided wave simulations can be conducted with most finite elements software packages, computational and hardware costs are always prohibitive for large simulation campaigns. A novel SHM module has been recently added to the CIVA software and relies on unassembled high order finite elements to overcome these limitations. This paper focuses on the thorough validation of CIVA for SHM to identify the limits of the models. After introducing the key elements of the CIVA SHM solution, a first validation is presented on a stainless steel pipe representative of the oil and gas industry. Second, validation is conducted on a composite panel with and without stiffener representative of some structures in the aerospace industry. Results show an excellent match between the experimental and simulated datasets, but only if the input parameters are fully determined prior to the simulations.

Adverse environmental conditions result in corrosion during the life cycle of marine structures such as pipelines, offshore oil platforms, and ships. Generalized corrosion leading to the loss of wall thickness can cause the degradation of the integrity, strength, and load bearing capacity of the structure. Nondestructive detection and monitoring of corrosion damage in difficult to access areas can be achieved using high-frequency guided waves propagating along the structure. Using standard ultrasonic wedge transducers with single-sided access to the structure, specific high-frequency guided wave modes (overlap of both fundamental Lamb wave modes) were generated that penetrate through the complete thickness of the structure. The wave propagation and interference of the guided wave modes depend on the thickness of the structure and were measured using a noncontact laser interferometer. Numerical simulations using a two-dimensional finite element model were performed to visualize and predict the guided wave propagation and energy transfer across the plate thickness. During laboratory experiments, the wall thickness was reduced uniformly by milling of one steel plate specimen. In a second step, wall thickness reduction was induced using accelerated corrosion for two mild steel plates. The corrosion damage was monitored based on the effect on the wave propagation and interference (beating effect) of the Lamb wave modes in the frequency domain. Good agreement of the measured beatlengths with theoretical predictions was achieved, and the sensitivity of the methodology was ascertained, showing that high-frequency guided waves have the potential for corrosion damage monitoring at critical and difficult to access locations.

Timber poles are widely used in electricity transmission and telecommunication sectors throughout the world. The stress wave propagation for the condition assessment of timber poles is identified as a promising non-destructive testing (NDT) technique due to its simplicity and cost-effectiveness compared to other traditional methods. In this paper, a novel damage severity evaluation criterion for timber poles is proposed on the basis of short-time wavelet entropy of the reflected stress waves. The stress waves are generated by transverse impacts close to the ground level of the pole. The reflected stress waves are recorded and processed in the time frequency domain using the discrete wavelet transform. The decomposed signal components using discrete wavelet analysis are used to determine the wavelet entropy. The wavelet entropies of intact and damaged poles are compared to obtain the relative wavelet entropy (RWE) for damage severity estimation. Further, a numerical model for an in situ pole system is developed to simulate the transverse stress wave propagation and to evaluate the capability of the proposed defect severity estimation method. The developed numerical model is validated with experimental data from controlled testing and the data from field tests. The validated numerical model is then used to simulate different defect scenarios. The wavelet entropy is sensitive to the damage severity in timber poles and can be used as an effective tool to evaluate the severity of damages.

This paper presents a modeling and simulation method for studying ultrasonic guided wave propagation in hybrid metal-composites, also known as fiber-metal laminates. The objective is to develop an efficient and versatile modeling tool to aid in the design of cost-effective nondestructive evaluation technologies. The global–local method, which combines finite element discretization and Lamb wave modal expansion is used. An extension to the traditional global–local method is made to couple the source problem with the scattering problem to deal with a surface source generating Lamb waves that interact with defects in multilayered structures. This framework is used to study the sensitivity of different excitation frequencies to ply gap defects of various sizes. The coupled model considers the transducer contact conditions and the ultrasonic system response in the Lamb wave excitation, along with the scattering phenomenon caused by the defects. This combined result is used to define the optimal excitation frequency for the strongest transmission or reflection for a given defect size that can be observed in a physical experiment. Such results can be applied to the design of a damage detection scheme in realistic aerospace structures.

Elastodynamic Green's function for anisotropic solids is required for wave propagation modeling in composites. Such modeling is needed for the interpretation of experimental results generated by ultrasonic excitation or mechanical vibration-based nondestructive evaluation tests of composite structures. For isotropic materials, the elastodynamic Green’s function can be obtained analytically. However, for anisotropic solids, numerical integration is required for the elastodynamic Green's function computation. It can be expressed as a summation of two integrals—a singular integral and a nonsingular (or regular) integral. The regular integral over the surface of a unit hemisphere needs to be evaluated numerically and is responsible for the majority of the computational time for the elastodynamic Green's function calculation. In this paper, it is shown that for transversely isotropic solids, which form a major portion of anisotropic materials, the integration domain of the regular part of the elastodynamic time-harmonic Green's function can be reduced from a hemisphere to a quarter-sphere. The analysis is performed in the frequency domain by considering time-harmonic Green's function. This improvement is then applied to a numerical example where it is shown that it nearly halves the computational time. This reduction in computational effort is important for a boundary element method and a distributed point source method whose computational efficiencies heavily depend on Green's function computational time.

Ultrasonic guided waves (GWs) are being extensively investigated and applied to nondestructive evaluation and structural health monitoring. Guided waves are, under most circumstances, excited in a frequency range up to several hundred kilohertz or megahertz for detecting defect/damage effectively. In this regard, numerical simulation using finite element analysis (FEA) offers a powerful tool to study the interaction between wave and defect/damage. Nevertheless, the simulation, based on linear/quadratic interpolation, may be inaccurate to depict the complex wave mode shape. Moreover, the mass lumping technique used in FEA for diagonalizing mass matrix in the explicit time integration may also undermine the calculation accuracy. In recognition of this, a time domain spectral element method (SEM)—a high-order FEA with Gauss–Lobatto–Legendre (GLL) node distribution and Lobatto quadrature algorithm—is studied to accurately model wave propagation. To start with, a simplified two-dimensional (2D) plane strain model of Lamb wave propagation is developed using SEM. The group velocity of the fundamental antisymmetric mode ( A 0 ) is extracted as indicator of accuracy, where SEM exhibits a trend of quick convergence rate and high calculation accuracy (0.03% error). A benchmark study of calculation accuracy and efficiency using SEM is accomplished. To further extend SEM-based simulation to interpret wave propagation in structures of complex geometry, a three-dimensional (3D) SEM model with arbitrary in-plane geometry is developed. Three-dimensional numerical simulation is conducted in which the scattering of A 0 mode by a through hole is interrogated, showing a good match with experimental and analytical results.

This paper addresses the predictive simulation of acoustic emission (AE) guided waves that appear due to sudden energy release during incremental crack propagation. The Helmholtz decomposition approach is applied to the inhomogeneous elastodynamic Navier–Lame equations for both the displacement field and body forces. For the displacement field, we use the usual decomposition in terms of unknown scalar and vector potentials, Φ and H . For the body forces, we hypothesize that they can also be expressed in terms of excitation scalar and vector potentials, A * and B * . It is shown that these excitation potentials can be traced to the energy released during an incremental crack propagation. Thus, the inhomogeneous Navier–Lame equation has been transformed into a system of inhomogeneous wave equations in terms of known excitation potentials A * and B * and unknown potentials Φ and H . The solution is readily obtained through direct and inverse Fourier transforms and application of the residue theorem. A numerical study of the one-dimensional (1D) AE guided wave propagation in a 6 mm thick 304-stainless steel plate is conducted. A Gaussian pulse is used to model the growth of the excitation potentials during the AE event; as a result, the actual excitation potential follows the error function variation in the time domain. The numerical studies show that the peak amplitude of A0 signal is higher than the peak amplitude of S0 signal, and the peak amplitude of bulk wave is not significant compared to S0 and A0 peak amplitudes. In addition, the effects of the source depth, higher propagating modes, and propagating distance on guided waves are also investigated.

This paper presents the local interaction simulation approach (LISA) for efficient modeling of linear and nonlinear ultrasonic guided wave active sensing of complex structures. Three major modeling challenges are considered: material anisotropy with damping effects, nonlinear interactions between guided waves and structural damage, as well as geometric complexity of waveguides. To demonstrate LISA's prowess in addressing such challenges, carefully designed numerical case studies are presented. First, guided wave propagation and attenuation in carbon fiber composite panels are simulated. The numerical results are compared with experimental measurements obtained from scanning laser Doppler vibrometry (SLDV) to illustrate LISA's capability in modeling damped wave propagation in anisotropic medium. Second, nonlinear interactions between guided waves and structural damage are modeled by integrating contact dynamics into the LISA formulations. Comparison with commercial finite element software reveals that LISA can accurately simulate nonlinear ultrasonics but with much higher efficiency. Finally, guided wave propagation in geometrically complex waveguides is studied. The numerical example of multimodal guided wave propagation in a rail track structure with a fatigue crack is presented, demonstrating LISA's versatility to model complex waveguides and arbitrary damage profiles. This article serves as a comprehensive, systematic showcase of LISA's superb capability for efficient modeling of transient dynamic guided wave phenomena in structural health monitoring (SHM).