Elsevier

Ultrasonics

Volume 54, Issue 6, August 2014, Pages 1553-1558
Ultrasonics

Nonlinear guided waves in plates: A numerical perspective

https://doi.org/10.1016/j.ultras.2014.04.009Get rights and content

Highlights

  • Harmonic generation from fundamental modes in plate is numerically investigated.

  • Role of material and geometric nonlinearities in harmonic generation is assessed.

  • It was found that group velocity matching is not necessary for harmonic generation.

  • Phenomenon of mode-interaction is demonstrated.

Abstract

Harmonic generation from non-cumulative fundamental symmetric (S0) and antisymmetric (A0) modes in plate is studied from a numerical standpoint. The contribution to harmonic generation from material nonlinearity is shown to be larger than that from geometric nonlinearity. Also, increasing the magnitude of the higher order elastic constants increases the amplitude of second harmonics. Second harmonic generation from non-phase-matched modes illustrates that group velocity matching is not a necessary condition for harmonic generation. Additionally, harmonic generation from primary mode is continuous and once generated, higher harmonics propagate independently. Lastly, the phenomenon of mode-interaction to generate sum and difference frequencies is demonstrated.

Introduction

Use of nonlinear ultrasound for characterizing microstructure of structural materials, especially metals, has been a topic of interest for several decades. Initial investigations by Breazeale and Thompson [1] and Hikata et al. [2] put forth elastic material nonlinearity and dislocations as predominant causes for nonlinear ultrasonic behavior i.e., higher harmonic generation. These results motivated the use of higher harmonic generation for investigating nonlinear behavior and hence material degradation in structures. Cantrell and Yost [3] investigated fatigued microstructures using nonlinear ultrasound. Cantrell [4] presented a comprehensive approach to relate the nonlinearity parameter (β) to the dislocation substructures in metals. For polycrystalline nickel, a monotonic increase in β with the fatigue cycles was predicted. The above theoretical/experimental investigations employed bulk waves (which travel in unbounded media) to characterize nonlinearity. On the other hand, guided waves (which travel in bounded structures) offer several advantages from an inspection standpoint in that long-range inspection can be carried out from a single location. Hence these are more amenable for structural health monitoring applications. Nonlinear guided waves combine the penetration power of guided waves with the early damage detection capabilities of nonlinear ultrasound. Hence they have emerged as an attractive alternative for detecting microstructural changes preceding macro-scale damage in the structures.

Deng [5], [6] analyzed second harmonic generation from guided waves in plates. De Lima and Hamilton [7] presented a perturbation based approach to analyze second harmonic generation using the normal mode expansion technique [8] and arrived at two conditions necessary for cumulative second harmonic generation, namely phase matching and non-zero power flux. While bulk waves satisfy the above criterion for all frequencies of excitation, only specific guided wave modes satisfy them. These were identified [9], [10] from the dispersion relations governing Rayleigh Lamb modes in the plate and experimental investigations [11], [12], [13], [14] corroborated theoretical predictions. While the above investigations dealt with second harmonic generation, theoretical investigations were also carried out for predicting the nature of higher harmonic generation from guided waves in plates. Srivastava and Lanza di Scalea [15] concluded that cumulative even harmonics exist only as symmetric modes while odd harmonics can exist either as symmetric or antisymmetric modes. Chillara and Lissenden [16] presented a generalized theory to study the nonlinear interaction of guided wave modes. They concluded that the interaction of guided wave modes of the same nature generate symmetric modes while those of opposite nature generate antisymmetric modes. They also proposed a procedure to predict the nature of higher harmonics from the theory of mode interaction developed.

While a comprehensive theoretical framework is now available to study higher harmonic generation from guided waves in plates, some issues still need to be addressed. These stem from the following issues:

  • 1.

    The theoretical analysis is carried out for time-harmonic (single-frequency continuous wave) excitations while the experiments employ transducers with finite band-width.

  • 2.

    Group velocity matching in addition to the phase velocity matching between the primary mode and other higher harmonics was assumed to be necessary for cumulative harmonic generation [17]. However, Deng et al. [18] presented an analysis and an experimental investigation confirming that it is not a necessary condition. In general, all the primary Rayleigh-Lamb modes capable of cumulative second harmonic generation in plates satisfy group velocity matching in addition to phase velocity matching [9]. However, with finite wave packets, it would be interesting to examine if the higher harmonic modes can propagate without group velocity matching as the non-zero-power-flux criterion which requires power flow from primary to higher harmonic modes may not be satisfied once the primary and higher modes separate due to difference in group velocity.

  • 3.

    The analysis is generally carried out using a perturbation approach as no closed form solution is available.

  • 4.

    The capability of the existing constitutive model to capture the effect of different kinds of nonlinearities has not been established.

In this article, we address some of the above issues using results from numerical simulations carried out in COMSOL Multiphysics 4.3, a commercial finite element software. Specifically, we consider the harmonic generation from the primary Symmetric (S0) and Antisymmetric (A0) modes in a plate. While neither of the above modes are exactly phase matched, they are deliberately chosen to highlight some of the aspects to be discussed later. As part of the study, we also comment on some aspects of the constitutive model like the contribution of material and geometric nonlinearities to higher harmonic generation.

The content of the article is organized as follows. Section 2 presents some continuum mechanics preliminaries. Section 3 then presents the results obtained from numerical simulations. Finally, conclusions are drawn in Section 4.

Section snippets

Preliminaries

In this section, we present some preliminaries intended to enhance understanding of the results to be presented. The discussion is brief so we refer the reader to [19] for more details.

We denote the deformation gradient by F and the Lagrangian strain by E. Also, the displacement gradient is denoted by H and the following relations exist between the above quantities.F=I+HE=12(FTF-I)=12(H+HT+HTH)We denote the linearized strain by El, which is related to the displacement gradient (H) byEl=12(H+HT).

Numerical simulations

In this section, we present results obtained from numerical simulations performed in COMSOL 4.3. All the simulations were carried out using the Murnaghan model (Eq. (5)) for Aluminum, the elastic constants of which are tabulated in Table 1. As mentioned earlier, fundamental symmetric and antisymmetric modes are used in the simulation. The schematic of the 2D model used for the simulation is shown in Fig. 1. The thickness of the plate is chosen to be 1 mm and the length of the plate is assumed to

Conclusions

In this article we investigated the behavior of nonlinear guided waves in plates from a numerical perspective. Simulations of nonlinear wave propagation were carried out in COMSOL using the Murnaghan model of hyper-elasticity. Harmonic generation from the fundamental guided wave modes S0 and A0 is studied. Solutions from ‘linear’, ‘nonlinear’ and ‘geometrically nonlinear’ variants of the Murnaghan model are compared for S0 mode at 0.5 MHz. While time domain signals show very little difference,

Acknowledgements

This material is based upon work supported by the Nuclear Energy Universities Program under Award number 00102946 and the National Science Foundation under Award number 1300562.

References (24)

  • B.A. Auld

    Acoustic Fields and Waves in Solids

    (1990)
  • N. Matsuda et al.

    Phase and group velocity matching for cumulative harmonic generation in Lamb waves

    J. Appl. Phys.

    (2011)
  • Cited by (0)

    View full text