Analytical and higher order finite element hybrid approach for an efficient simulation of ultrasonic guided waves I: 2D-analysis

In recent years the interest in online monitoring of lightweight structures with ultrasonic guided waves is steadily growing. Especially the aircraft industry is a driving force in the development of structural health monitoring (SHM) systems. In order to optimally design SHM systems powerful and efficient numerical simulation tools to predict the behaviour of ultrasonic elastic waves in thin-walled structures are required. It has been shown that in real industrial applications, such as airplane wings or fuselages, conventional linear and quadratic pure displacement finite elements commonly used to model ultrasonic elastic waves quickly reach their limits. The required mesh density, to obtain good quality solutions, results in enormous computational costs when solving the wave propagation problem in the time domain. To resolve this problem different possibilities are available. Analytical methods and higher order finite element method approaches (HO-FEM), like p-FEM, spectral elements, spectral analysis and isogeometric analysis, are among them. Although analytical approaches offer fast and accurate results, they are limited to rather simple geometries. On the other hand, the application of higher order finite element schemes is a computationally demanding task. The drawbacks of both methods can be circumvented if regions of complex geometry are modelled using a HO-FEM approach while the response of the remaining structure is computed utilizing an analytical approach. The objective of the paper is to present an efficient method to couple different HO-FEM schemes with an analytical description of an undisturbed region. Using this hybrid formulation the numerical effort can be drastically reduced. The functionality of the proposed scheme is demonstrated by studying the propagation of ultrasonic guided waves in plates, excited by a piezoelectric patch actuator. The actuator is modelled utilizing higher order coupled field finite elements, whereas the homogenous, isotropic plate is described analytically. The results of this "semi-analytical" approach highlight the opportunities to reduce the numerical effort if closed-form solutions are partially available.