Next Article in Journal
Advanced Chemical Looping Materials for CO2 Utilization: A Review
Next Article in Special Issue
Properties of Love Waves in Functional Graded Saturated Material
Previous Article in Journal
A Review of the Characteristics, Synthesis, and Thermodynamics of Type-II Weyl Semimetal WTe2
Previous Article in Special Issue
Neuro-Fuzzy Modelling of the Metallic Surface Characterization on Linear Dry Contact between Plastic Material Reinforced with SGF and Alloyed Steel
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Computation of Complex Dispersion and Properties of Evanescent Lamb Wave in Functionally Graded Piezoelectric-Piezomagnetic Plates

School of Mechanical and Power Engineering, Henan Polytechnic University, Jiaozuo 454000, China
*
Author to whom correspondence should be addressed.
Materials 2018, 11(7), 1186; https://doi.org/10.3390/ma11071186
Submission received: 24 May 2018 / Revised: 2 July 2018 / Accepted: 7 July 2018 / Published: 10 July 2018

Abstract

:
Functionally graded piezoelectric-piezomagnetic material (FGPPM), with a gradual variation of the material properties in the desired direction(s), can improve the conversion of energy among mechanical, electric, and magnetic fields. Full dispersion relations and wave mode shapes are vital to understanding dynamic behaviors of structures made of FGPPM. In this paper, an analytic method based on polynomial expansions is proposed to investigate the complex-valued dispersion and the evanescent Lamb wave in FGPPM plates. Comparisons with other related studies are conducted to validate the correctness of the presented method. Characteristics of the guided wave, including propagating modes and evanescent modes, in various FGPPM plates are studied, and three-dimensional full dispersion and attenuation curves are plotted to gain a deeper insight into the nature of the evanescent wave. The influences of the gradient variation on the dispersion and the magneto-electromechanical coupling factor are illustrated. The displacement amplitude and electric potential and magnetic potential distributions are also discussed in detail. The obtained numerical results could be useful to design and optimize different sensors and transducers made of smart piezoelectric and piezomagnetic materials with high performance by adjusting the gradient property.

1. Introduction

Due to the excellent coupling behavior among mechanical, electric, and magnetic fields, piezoelectric-piezomagnetic composites (or magneto-electro-elastic material) composed of piezoelectric and piezomagnetic phases have been increasingly applied to different engineering structures, especially to the smart or intelligent systems as intelligent sensors, damage detectors, etc. [1]. It is found that the smart structures made of functionally graded materials (FGM) possess a better structural performance than traditional composite materials. The concept of FGM has been extended to the development of new piezoelectric-piezomagnetic materials appointed functionally graded piezoelectric-piezomagnetic materials, which can realize the smooth transition of the physical constitutive parameters of the piezoelectric and piezomagnetic materials. FGPPMs have been used in some devices to improve their efficiency and other features. Many applications are closely connected with the vibration and wave propagation of FGM and FGPPM [2,3,4,5]. Dispersion relations and wave mode shapes are very important for understanding dynamic behaviors of structures. Wave propagation features in FGPPM plates could be also useful in designing and optimizing the high-accuracy sensors and transducers [6,7].
With the remarkable achievements in fabrication of FGPPM during the decades, many investigators have turned attention to the study of wave propagation in such materials. Wang and Rokhlin [8] presented the differential equations governing the transfer and stiffness matrices for a functionally graded generally anisotropic magneto-electro-elastic medium and calculated the surface wave velocity dispersion. Pan et al. [9] derived an exact solution for the multilayered plate made of functionally graded, anisotropic and linear magneto-electro-elastic materials based on Pseudo-Stroh formalism. Bhangale et al. [10] carried out the free vibration studies on the simply supported functionally graded magneto-electro-elastic plate by semi-analytical finite element method. Wu et al. [11] investigated the wave propagating characteristics in the non-homogeneous magneto-electro-elastic plates by using orthogonal polynomial approach. By employing the power series technique, Cao [12] investigated the Lamb wave propagation in FGPPM plates. Singh and Rokne [13] investigated the SH wave propagating in FGPPM structures. Xiao et al. [14] investigated the dispersion properties of wave propagation in the functionally graded magneto-electro-elastic plate by the Chebyshev spectral element method.
As is reviewed above, so far, studies on the guided wave in FGPPM structures are limited to the propagating waves, but the evanescent waves have not been investigated. Recently, some studies on pseudo surface acoustic waves (PSAW) in piezoelectric half-spaces find that the PSAW modes have higher velocities and lower attenuations, compared to the classical surface acoustic waves [15,16]. Such modes make the piezoelectric device possess higher resolution. Evanescent wave modes also have the similar features. According to the classification of Auld [17], the complete wave modes consist of propagating modes with real wave number and evanescent modes with complex or purely imaginary wave number. Note that evanescent modes represent local modes that would exist at discontinuities and decay with propagating distance (so referred to as evanescent or non-propagating wave). As early as 1955, Lyon [18] obtained the purely imaginary roots of the dispersion equation for an elastic plate. Remarkable is the work done by Mindlin who demonstrated the presence of complex roots of the Rayleigh-Lamb equation [19]. Freedman [20] studied the imaginary valued Lamb mode spectra covering virtually the full range of the Poisson ratio. Quintanilla et al. [21] calculated the full spectrum for guided wave problems in plates and layered cylinders using a spectral collocation method. More recently, Yan and Yuan [22,23] discussed the potential application of evanescent waves in structural health monitoring and investigated the conversion of evanescent SH and Lamb waves into propagating waves using a semi-analytical approach. Chen et al. [24] studied theoretically the real-valued and imaginary-valued SH waves in a piezoelectric plate of cubic crystals. These researches focused on the simple material and purely imaginary modes. In fact, the search of complex roots corresponding to evanescent waves is a difficult task for FGM with material properties of variable coefficients. To the best of the authors’ knowledge, the evanescent waves in FGM or piezoelectric-piezomagnetic composite have not been studied before, which is the motivation of this study.
In this paper, an analytic method based on polynomial expansions is proposed to calculate guided waves in FGPPM plates. The presented method can replace the problem of computing a transcendental dispersion equation by a general eigenvalue problem in wave number. The complete solutions of the dispersion equation, including the purely real, purely imaginary and complex solutions, can be obtained. We plot the full dispersion curves in three dimensional (3D) frequency-complex wave number space to gain a better and deeper insight into the characteristics of evanescent waves. Two known cases are given to validate this approach. The characteristics of evanescent guided waves in various FGPPM plates are illustrated. The effects of different graded fields on the dispersion curves and the coupled electromechanical factor are investigated. The displacement amplitude and electric potential and magnetic potential distributions are also discussed in detail.

2. Mathematics and Formulation of the Problem

Consider a FGPPM plate with varying material properties with regard to thickness (the z-axis). The plate described in Cartesian coordinate system (x, y, z), is infinite horizontally but finite in the z direction with a thickness h, occupies the region 0 ≤ zh, as shown in Figure 1. The wave propagates along the x direction, and the upper and bottom surface of the plate are traction free.
For the piezoelectric-piezomagnetic medium, the governing field equations (Equation (1)) and the generalized constitutive relations (Equation (2)) can be expressed as [12]:
σ i j , j = ρ u ¨ i , D i , i = 0 , B i , i = 0
σ i j = C i j k l ε k l e k i j E k q i j k H k , D i = e i k l ε k l + i k E k + g i k H k , B i = q i j k ε j k + g i k E k + μ i k H k
ε i j = 1 2 ( u i , j + u j , i ) , E i = φ , i , H i = ψ , i
Generalized geometric equations under a rectangular coordinates system are in the above Equations (1)–(3), σ i j is the stress tensor, D i is the electric displacement and B i is the magnetic induction C i j k l , e k i j , q i j k , i k , g i k and μ i k are the elastic, piezoelectric, piezomagnetic, dielectric, magnetic, and magnetoelectric parameters of the FGPPM, respectively, while all of them, including the mass density ρ, are functions of z. According to Einstein summation convention, where i, j, k and l = 1, 2, 3 corresponding to x, y, z directions, respectively. ε k l , E k and H k are the strain tensor, the electric field, and magnetic field, respectively. ui (i = x, y, z) denotes the mechanical displacement component in the ith direction. φ and ψ are the electric potential and magnetic potential. Comma in subscripts and superposed dot denote spatial and time derivatives, respectively.
For Lamb waves propagating along the x direction, the displacement components, electric potential and magnetic potential can be expressed as
u i ( x , z , t ) = exp ( i k x i ω t ) U i ( z ) , φ ( x , z , t ) = exp ( i k x i ω t ) X ( z ) , ψ ( x , z , t ) = exp ( i k x i ω t ) Y ( z )
where Ui (i = x, z) represents the amplitude of the displacements in the ith directions, X and Y represent the amplitude of electric potential and magnetic potential, respectively. k is the wave number, ω is the angular frequency, and i is the imaginary number.
Since the material properties change gradually with thickness and are the functions of z, they can be fitted into the following form:
f ( z ) = f ( l ) ( z / h ) l , l = 0 , 1 , 2 , L
where f (f = ρ, C, e, , q, g and μ) denotes material parameters, l is the order number, f(l) is the coefficient. For homogeneous material, f ( z ) = C ( 0 ) , and when l > 0, f(l) is zero.
The following boundary and continuous conditions should be satisfied as follows. For the traction-free boundary condition, it requires that σ z z | z = 0 , h = 0 , σ x z | z = 0 , h = 0 , σ y z | z = 0 , h = 0 . For electric and magnetic open circuit, D z | z = 0 , h = 0 , B z | z = 0 , h = 0 , and for electric and magnetic shorted circuit, φ | z = 0 , h = 0 , ψ | z = 0 , h = 0 .
Then take the traction-free and electrical and magnetic open-circuit boundary conditions as an example. Considering the boundary of material, the position-dependent material parameters are given by:
f ( z ) = f ( z ) π ( z )
where π(z) is a rectangular window function defined by π ( z ) = { 1 , 0 z h 0 , e l s e w h e r e , whose derivative is a Dirac’s delta function, δ ( z h ) δ ( z ) . Then the boundary conditions can be automatically incorporated in the constitutive relations [25].
To reduce the number of resolving equations, we substitute Equations (3)–(6) into Equation (2) with following substitution into Equation (1). Consequently, the governing differential equations in terms of the displacement, electric potential and magnetic potential components can be obtained. Here, the case of an orthotropic FGPPM plate with the z direction polarization is given:
( z h ) l [ C 55 ( l ) U + l z 1 C 55 ( l ) U + i k ( C 13 ( l ) + C 55 ( l ) ) W + i k ( e 15 ( l ) + e 31 ( l ) ) X + i k ( q 15 ( l ) + q 31 ( l ) ) Y k 2 C 11 ( l ) U + l i k z 1 ( C 55 ( l ) W + e 15 ( l ) X + q 15 ( l ) Y ) ] π ( z ) + ( δ ( z 0 ) δ ( z h ) ) ( z h ) l ( C 55 ( l ) U + i k C 55 ( l ) W + i k e 15 ( l ) X + i k q 15 ( l ) Y ) = ρ ( l ) z l ω 2 h l U π ( z )
( z h ) l [ C 33 ( l ) W + e 33 ( l ) X + q 33 ( l ) Y + i k ( C 13 ( l ) + C 55 ( l ) ) U + l z 1 ( C 33 ( l ) W + e 33 ( l ) X + q 33 ( l ) Y ) + l i k z 1 C 13 ( l ) U k 2 ( C 55 ( l ) W + e 15 ( l ) X + q 15 ( l ) Y ) ] π ( z ) + ( δ ( z 0 ) δ ( z h ) ) ( z h ) l ( C 33 ( l ) W + e 33 ( l ) X + q 33 ( l ) Y + i k C 13 ( l ) U ) = ρ ( l ) z l ω 2 h l W π ( z )
( z h ) l [ e 33 ( l ) W 33 ( l ) X g 33 ( l ) Y + i k ( e 15 ( l ) + e 31 ( l ) ) U + l z 1 ( e 33 ( l ) W 33 ( l ) X g 33 ( l ) Y ) + l i k z 1 e 31 ( l ) U k 2 e 15 ( l ) W + k 2 11 ( l ) X + k 2 g 11 ( l ) Y ] π ( z ) + ( δ ( z 0 ) δ ( z h ) ) ( z h ) l ( e 33 ( l ) W 33 ( l ) X g 33 ( l ) Y + i k e 31 ( l ) U ) = 0
( z h ) l [ q 33 ( l ) W g 33 ( l ) X μ 33 ( l ) Y + i k ( q 15 ( l ) + q 31 ( l ) ) U + l i k z 1 q 31 ( l ) U + l z 1 ( q 33 ( l ) W g 33 ( l ) X μ 33 ( l ) Y ) k 2 q 15 ( l ) W + k 2 g 11 ( l ) X + k 2 μ 11 ( l ) Y ] π ( z ) + ( δ ( z 0 ) δ ( z h ) ) ( z h ) l ( q 33 ( l ) W g 33 ( l ) X μ 33 ( l ) Y + i k q 31 ( l ) U ) = 0
where U and W respectively represent the amplitude of vibration in the x and z directions. The superscript (’) is the derivative with respect to z.
The four amplitudes can be expanded into Legendre orthogonal polynomial series as:
U ( z ) = m = 0 p m 1 Q m ( z ) ,   W ( z ) = m = 0 p m 2 Q m ( z ) ,   X ( z ) = m = 0 p m 3 Q m ( z ) ,   Y ( z ) = m = 0 p m 4 Q m ( z )
where p m α ( α = 1 , 2 , 3 , 4 ) are the expansion coefficients, Q m ( r ) are an orthonormal set of polynomials in the interval [0,h].
Q m ( z ) = 2 m + 1 h P m ( 2 z h h )
where Pm is the Legendre polynomial of order m.
Substituting Equations (8) and (9) into Equation (7), then multiplying both sides of the modified Equation (7) by the complex conjugate Q j * ( z ) with j running from 0 to M, integrating over z from 0 to h, taking advantage of the orthonormality of the polynomial, yields:
k 2 [ A l 11 j , m A l 12 j , m A l 13 j , m A l 14 j , m A l 21 j , m A l 22 j , m A l 23 j , m A l 24 j , m A l 31 j , m A l 32 j , m A l 33 j , m A l 34 j , m A l 41 j , m A l 42 j , m A l 43 j , m A l 44 j , m ] { p m 1 p m 2 p m 3 p m 4 } + k [ B l 11 j , m B l 12 j , m B l 13 j , m B l 14 j , m B l 21 j , m B l 22 j , m B l 23 j , m B l 24 j , m B l 31 j , m B l 32 j , m B l 33 j , m B l 34 j , m B l 41 j , m B l 42 j , m B l 43 j , m B l 44 j , m ] { p m 1 p m 2 p m 3 p m 4 } + [ C l 11 j , m C l 12 j , m C l 13 j , m C l 14 j , m C l 21 j , m C l 22 j , m C l 23 j , m C l 24 j , m C l 31 j , m C l 32 j , m C l 33 j , m C l 34 j , m C l 41 j , m C l 42 j , m C l 43 j , m C l 44 j , m ] { p m 1 p m 2 p m 3 p m 4 } = ω 2 [ M l m j 0 0 0 0 M l m j 0 0 0 0 0 0 0 0 0 0 ] { p m 1 p m 2 p m 3 p m 4 }
or is abbreviated as
k 2 A p + k 1 B p + C p = ω 2 M p
where A, B, C and M are matrices of order 4(M + 1)·(M + 1), p = [ p m 1 p m 2 p m 3 p m 4 ] T , the elements of the matrices are as following,
A l 11 j , m = 1 h l C 11 ( l ) β ( m , l , 0 , j ) A l 22 j , m = 1 h l C 55 ( l ) β ( m , l , 0 , j ) A l 23 j , m = 1 h l e 15 ( l ) β ( m , l , 0 , j ) A l 24 j , m = 1 h l q 15 ( l ) β ( m , l , 0 , j ) A l 32 j , m = 1 h l e 15 ( l ) β ( m , l , 0 , j ) A l 33 j , m = 1 h l 11 ( l ) β ( m , l , 0 , j ) A l 34 j , m = 1 h l g 11 ( l ) β ( m , l , 0 , j ) A l 42 j , m = 1 h l q 15 ( l ) β ( m , l , 0 , j ) A l 43 j , m = 1 h l g 11 ( l ) β ( m , l , 0 , j ) A l 44 j , m = 1 h l μ 11 ( l ) β ( m , l , 0 , j ) A l 12 j , m = A l 21 j , m = 0 A l 13 j , m = A l 31 j , m = 0 A l 14 j , m = A l 41 j , m = 0 ; B l 12 j , m = 1 h l { i ( C 13 ( l ) + C 55 ( l ) ) β ( m , l , 1 , j ) + l i C 55 ( l ) β ( m , l 1 , 0 , j ) + i C 55 ( l ) γ ( m , l , 0 , j ) } , B l 13 j , m = 1 h l { i ( e 15 ( l ) + e 31 ( l ) ) β ( m , l , 1 , j ) + l i e 15 ( l ) β ( m , l 1 , 0 , j ) + i e 15 ( l ) γ ( m , l , 0 , j ) } B l 14 j , m = 1 h l { i ( q 15 ( l ) + q 31 ( l ) ) β ( m , l , 1 , j ) + l i q 15 ( l ) β ( m , l 1 , 0 , j ) + i q 15 ( l ) γ ( m , l , 0 , j ) } B l 21 j , m = 1 h l { i ( C 13 ( l ) + C 55 ( l ) ) β ( m , l , 1 , j ) + l i C 13 ( l ) β ( m , l 1 , 0 , j ) + i C 13 ( l ) γ ( m , l , 0 , j ) } B l 31 j , m = 1 h l { i ( e 15 ( l ) + e 31 ( l ) ) β ( m , l , 1 , j ) + l i e 31 ( l ) β ( m , l 1 , 0 , j ) + i e 31 ( l ) γ ( m , l , 0 , j ) } B l 41 j , m = 1 h l { i ( q 15 ( l ) + q 31 ( l ) ) β ( m , l , 1 , j ) + l i q 31 ( l ) β ( m , l 1 , 0 , j ) + i q 31 ( l ) γ ( m , l , 0 , j ) } , B l 11 j , m = B l 22 j , m = B l 33 j , m = B l 44 j , m = 0 B l 23 j , m = B l 32 j , m = 0 B l 24 j , m = B l 42 j , m = 0 B l 34 j , m = B l 43 j , m = 0 ; C l 11 j , m = 1 h l { C 55 ( l ) β ( m , l , 2 , j ) + l C 55 ( l ) β ( m , l 1 , 1 , j ) + C 55 ( l ) γ ( m , l , 1 , j ) } C l 22 j , m = 1 h l { C 33 ( l ) β ( m , l , 2 , j ) + l C 33 ( l ) β ( m , l 1 , 1 , j ) + C 33 ( l ) γ ( m , l , 1 , j ) } C l 23 j , m = 1 h l { e 33 ( l ) β ( m , l , 2 , j ) + l e 33 ( l ) β ( m , l 1 , 1 , j ) + e 33 ( l ) γ ( m , l , 1 , j ) } C l 24 j , m = 1 h l { q 33 ( l ) β ( m , l , 2 , j ) + l q 33 ( l ) β ( m , l 1 , 1 , j ) + q 33 ( l ) γ ( m , l , 1 , j ) } , C l 32 j , m = 1 h l { e 33 ( l ) β ( m , l , 2 , j ) + l e 33 ( l ) β ( m , l 1 , 1 , j ) + e 33 ( l ) γ ( m , l , 1 , j ) } C l 33 j , m = 1 h l { 33 ( l ) β ( m , l , 2 , j ) l 33 ( l ) β ( m , l 1 , 1 , j ) 33 ( l ) γ ( m , l , 1 , j ) } C l 34 j , m = 1 h l { g 33 ( l ) β ( m , l , 2 , j ) l g 33 ( l ) β ( m , l 1 , 1 , j ) g 33 ( l ) γ ( m , l , 1 , j ) } C l 42 j , m = 1 h l { q 33 ( l ) β ( m , l , 2 , j ) + l q 33 ( l ) β ( m , l 1 , 1 , j ) + q 33 ( l ) γ ( m , l , 1 , j ) } C l 43 j , m = 1 h l { g 33 ( l ) β ( m , l , 2 , j ) l g 33 ( l ) β ( m , l 1 , 1 , j ) g 33 ( l ) γ ( m , l , 1 , j ) } C l 44 j , m = 1 h l { μ 33 ( l ) β ( m , l , 2 , j ) l μ 33 ( l ) β ( m , l 1 , 1 , j ) μ 33 ( l ) γ ( m , l , 1 , j ) } , C l 12 j , m = C l 21 j , m = 0 C l 13 j , m = C l 31 j , m = 0 C l 14 j , m = C l 41 j , m = 0 M l m j = 1 h l ρ ( l ) β ( m , l , 0 , j ) ; with   β ( m , l , n , j ) = 0 h Q j * ( z ) z l n Q m ( z ) z n d z , γ ( m , l , n , j ) = 0 h Q j * ( z ) z l π ( z ) z n Q m ( z ) z n d z .
The objective is to find wave numbers k that satisfy the Equation (11). It is simple and useful for propagating wave, by specifying real k and then solving for ω. But if interest is the evanescent wave, the approach is useless because k is complex and the solving of Equation (11) involves a multivariable search. In order to overcome this difficulty, we develop a new solution procedure as shown below.
We introduce two new vectors:
q = k p ,   N = ω 2 M .
Substitution Equation (12) into Equation (11), and then multiplying both sides of the modified Equation (11) by inverse matrix A 1 , yields
A 1 ( N C ) p ( A 1 B ) q = k q .
Combining Equation (13) and the above vector q = k p , we obtain
[ Z I 4 ( M + 1 ) A 1 ( N C ) A 1 B ] [ p q ] = k [ p q ] .
where I is the identity matrix and Z is a zero matrix.
If we define R = [ p q ] T , then Equation (14) can be written as
[ Z I 4 ( M + 1 ) A 1 ( N C ) A 1 B ] R = k R .
Up to this stage, the problem is reduced to a typical eigenvalue problem, which can be easily solved using an eigensolver routine that yields the complex eigenvalues k. All the developments performed in this paper were implemented in Mathematica software (version 8.0, Wolfram company, Champaign, IL, USA). The calculation technique in the short-circuit case is similar to that which is used in the open-circuit case. The deduction process is not shown to save space.

3. Numerical Results and Discussion

Based on the previous formulations, the computer program in terms of the presented method has been written using Mathematica software to calculate the dispersion and phase velocity curves for the FGPPM plate composed of CoFe2O4 (top) and Ba2TiO3 (bottom), h = 1 mm. The material parameters are from literature [26] and are listed in Table 1.
We use the Voigt-type model, as described in the literature [27], to calculate the effective material property of the FGPPM plate:
F ( z ) = F B V B ( z ) + F C V C ( z ) , V B ( z ) + V C ( z ) = 1
where FB and FC respectively represent the material property of the Ba2TiO3 and CoFe2O4 materials, and VB and VC are volume fraction.
Equation (16) can be rewritten as
F ( z ) = F B + ( F C F B ) V C ( z )
Similar to Equation (5), VC(z) can be expressed as a power expansion, Here we consider four different gradient fields, V C ( z ) = ( z / h ) n , n = 1, 2 and 3, namely linear, quadratic and cubic graded fields, and sinusoidal graded field V C ( z ) = sin ( 0.5 π z / h ) .

3.1. Approach Validation and Convergence of the Problem

To check the validity and the efficiency of our approach, we make a comparison between our results and the literature results. Because there is no investigation on the evanescent waves in FGPPM so far, we compute the full spectrum of Lamb wave in a steel plate and make a comparison with the available results in literature [17] from a spectral collocation method. The calculating parameters are ρ = 7932 kg/m3, C11 = 281.757 GPa, C12 = 113.161 GPa, C44 = 84.298 GPa, and h = 10 mm. The non-dimensional frequency and wave number are defined as Ω = ( ω h ρ / C 44 ) / π , Ψ = k h / π , respectively. The resulting dispersion curves are given in Figure 2. It clearly shows that the numerical results obtained by the present polynomial approach agree well with those obtained by the spectral collocation method, which validates our approach and program.
The material in the above verification example is isotropic. We calculate the dispersion curves of Lamb wave in an orthotropic plate and make a comparison with the available results in literature [28] from the reverberation-ray matrix method, which serves as a further validation of our approach. The material is PZT-4, and the material parameters are listed in Table 2. Figure 3 shows the obtained frequency spectra. Here again, the agreement is quite good between our results and those from the reverberation ray matrix method.
Then we discuss the convergence of the present polynomial approach. We present dispersion curves of propagating Lamb-like wave in a linear FGPPM plate with electric and magnetic open circuit and h = 1 mm, when the truncation order M takes 7, 8, 9 and 15, respectively, as shown in Figure 4. It can be seen that more and more order modes converge as M increases. When M = 7, the first three modes are convergent. The first four when M = 8, and the first seven when M = 9. So, we can think that at least the first (M − 1)/2 modes are convergent. Similarly, this can be concluded for the purely imaginary modes, and we don’t present the dispersion curves of purely imaginary branches for saving space. For evanescent Lamb-like waves, we tabulate the results in Table 3 since graph is not convenient for comparison. These numerical results also show that the complex solutions are convergent as M increases. When M = 10 and M = 11, the first three modes are convergent. The first four when M = 12, the first five when M = 13, and the first six when M = 14. Obviously, the real solution is easier to converge than the complex one. From these results, good convergence of the present approach can be observed. We take M = 30 in this paper.

3.2. Full Dispersion Curves of Lamb Wave

Propagating waves have received a lot of attention, and here we put the emphasis on evanescent waves. We plot the full dispersion curves in 3D frequency-complex wave number space for a clearer visualization of the solutions and a better understanding of the nature of the modes, when necessary, with a different color for clarity. Figure 5a plots the full dispersion curves of Lamb wave for a linear FGPPM plate with electric and magnetic open circuit. Since the eigenvalues are computed for one ω at a time, the dispersion curves are constructed of unconnected dots and the points near the cut-off frequencies become sparse. We can observe that purely real and purely imaginary solutions appear in pairs of opposite signs and the complex ones appear in quadruples of complex conjugates and opposite signs. Purely real wave numbers correspond to the propagating wave, and purely imaginary and complex wave numbers correspond to the evanescent wave. For a given frequency, a certain small number of real branches exist together with an infinite number of complex and purely imaginary branches (mostly imaginary with few complexes in the given range). For clarity, Figure 5b shows one quadrant dispersion curves in a small range. For complex branches, most of them start from 0 frequency and end at the minima of the purely real branches. Occasionally, one connecting two purely imaginary branches appears. The real part of the complex branches is usually small. For purely imaginary branches, most of them start from 0 frequency and end at cut-off frequencies with increasing frequency, and some with small wave numbers start from one cut-off frequency and terminate the other one.
Figure 6 shows the phase velocity dispersion and attenuation curves of the first three propagating and complex branches. The dimensionless phase velocity and frequency and attenuation are defined by Vp = ω / ( Re ( k ) C 55 / ρ ) , fh = ω h / ( 2 π C 55 / ρ ) and Im ( kh ) . We can find from these curves that the phase velocity of a propagating mode is decreased and gradually tends to a steady value with increasing frequency, but the velocity of an evanescent mode becomes bigger as well as the attenuation decreases. At high frequency, the evanescent mode has a very small attenuation, and its phase velocity is noticeably bigger than that of a propagating mode. For example, at fh = 2–3, the phase velocity of the second evanescent mode is about 8, but that of the propagating mode is below 2. Also, the wave dispersion is quite weak in this frequency range.

3.3. Influences of Graded Field on Dispersion Curves and the Electromechanical Coupling Factor

Considering two graded shapes, cubic and sinusoidal graded fields. Figure 7 shows their dispersion curves of Lamb wave. The results show that the effect of the graded field on dispersion characteristics of Lamb wave is significant, including the propagating modes and evanescent modes. Comparison between Figure 7 and Figure 5b, we can notice that the imaginary part of the complex branches for the sinusoidal graded case, at Ω = 0 plane, is bigger than that for the linear and cubic cases. Interestingly, for the cubic graded cases, the complex branch connecting two purely imaginary branches disappears and turns into a different one connecting a purely imaginary branch and a real branch. For clarity, Figure 8 shows the frequency spectra and phase velocity spectra of Lamb propagating wave for the three graded fields. Obviously, the effect of the graded field is little on the low mode, but becomes significant with increasing the mode order and wave number. The phase velocity for the sinusoidal graded field is bigger than that for the linear graded field, while the linear bigger than the cubic. The reason lies in that the different graded fields result in different material volume distributions, and the wave velocity depends on the material properties. Figure 9 gives the variation curves of the three gradient fields in the z direction. The CoFe2O4 content for the sinusoidal graded field is the highest, and the wave velocity of Ba2TiO3 is slower than that of CoFe2O4.
The magneto-electromechanical coupling factor K2 is an important parameter for designing acoustic wave devices. A high magneto-electromechanical coupling factor is expected in engineering applications. It is defined as [29]
K 2 = 2 | V o c V s c | V o c
where Voc and Vsc are the phase velocities for the electric and magnetic open circuit and short circuit, respectively.
To illustrate the effect of graded field on the K2, we calculate the K2 for S0 modes of four different FGPPM plates, as shown in Figure 10. We can find that the K2 reaches a maximum at a certain wave number and tend to the same little value with increasing wave number, which implies the influence of the graded field on the energy propagation of Lamb wave in high-frequency zone is insignificant. It reaches a maximum from 4.4% for the sinusoidal graded field to 9.5% for the cubic graded field. They are located near kh = 2 and kh = 1.5 respectively. The K2 for the cubic graded field is always bigger than that of the other three graded cases. Also the maximum of K2 shifts to the smaller wave number when the graded power exponent is increasing.

3.4. Wave Structure Analysis

The distributions of displacement and electric potential and magnetic potential fields can be obtained according to Equations (4) and (8). Considering a special position where the complex branch firstly collapses onto the real branch at about Ω = 1.0, as marked with a circle in Figure 5b. Figure 11 and Figure 12 present the distributions of the physical quantities in the z and x directions when Ω = 1.01115, Ψ = 0.23216 − 0.04612i, and Ω = 1.01911, Ψ = 0.17938, respectively. As seen in these figures, the real branch propagates without any attenuation, and the complex branch exhibits an oscillatory distribution and propagates a very long distance, about a few tens of thicknesses of the plate. The displacement uz and electric potential and magnetic potential distributions change along the z direction in a nearly anti-symmetric manner. The displacement ux exhibits a nearly symmetric manner. The distribution of displacement uz of the complex branch is very similar to that of the real branch, implying the evanescent wave mode converts into the propagating wave mode.

3.5. Merits of the Presented Method

Based on the above calculation of wave propagation in a FGPPM plate, we can summarize the following advantages of the presented method, which makes the method attractive.
(1)
The complex mathematical issue is reduced to solve an eigenvalue problem, which is capable of accurately determining all the real, imaginary and complex solutions of a transcendental dispersion equation.
(2)
The conventional approaches (root-finding routines or finite element simulations) require an iterative search procedure or a far greater coding effort, to find complex roots. The present method can avoid tedious iterative two-variable search and is simple to program. It needs to take a larger polynomial order to obtain solutions of the higher modes, which will cause more computer memory and long time.
(3)
The method is easy to implement and can be extended to complex structures such as multilayered or curved structures.

4. Conclusions

This paper presents an analytic method based on polynomial expansions for the determination of the full dispersion spectrum of the guided waves in the FGPPM plate. The correctness of the present method is verified via numerical comparison with available reference results. For the first time, the complete 3D dispersion curves of Lamb wave in a FGPPM plate are illustrated in a wide frequency range. The characteristics of the Lamb waves including the propagating and evanescent modes in various FGPPM plates are investigated. The emphasis on evanescent waves makes this work relevant for applications in the nondestructive evaluation of material or structural properties. Based on the above numerical results, some interesting conclusions can be drawn:
(1)
Superior to the conventional methods that necessitate an iterative search procedure to solve the complex roots of a dispersion equation, the presented analytic method can transform the set of differential equations for the acoustic waves into an eigenvalue problem in the form AX = kX to find the complex solutions.
(2)
Complex branches of the Lamb wave usually collapse onto the extremum of the real branches. They exhibit both local vibration and local propagation, and some can propagate a quite long distance (more than ten times of the plate thickness). They will turn into the propagating modes with increasing frequency.
(3)
Some evanescent modes have a noticeably higher phase velocity than the propagating modes. The phase velocity of the low order evanescent modes is more than four times larger than that of the propagating modes. Also, the wave dispersion of the evanescent mode is quite weak in a certain frequency range.
(4)
The magneto-electromechanical coupling factor of the guided wave in a FGPPM plate may be improved by adjusting the graded field. The coupling factor reaches a maximum from 4.4% for the sinusoidal graded field to 9.5% for the cubic graded field. The maximum of the magneto-electromechanical coupling factor for the S0 mode shifts to lower frequencies with increasing the gradient index.

Author Contributions

X.Z. and J.Y. proposed the studied problem and the corresponding solving method; X.Z. and Z.L. conducted the theoretical derivation and the computation; X.Z. and J.Y. analyzed and discussed the results; X.Z. and J.Y. wrote the paper.
Funding: This research was funded by the National Natural Science Foundation of China (No. U1504106), the fundamental research funds for the national outstanding youth project of Henan Polytechnic University (No. NSFRF140301), and the Program for Innovative Research Team of Henan Polytechnic University (T2017-3).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, J.G.; Chen, L.F.; Fang, S.S. State vector approach to analysis of multilayered magneto-electro-elastic plates. Int. J. Solids Struct. 2003, 40, 1669–1680. [Google Scholar] [CrossRef]
  2. Tornabene, F. Free vibration analysis of functionally graded conical, cylindrical shell and annular plate structures with a four-parameter power-law distribution. Comput. Methods Appl. Mech. Eng. 2009, 198, 2911–2935. [Google Scholar] [CrossRef]
  3. Kandasamy, R.; Dimitri, R.; Tornabene, F. Numerical study on the free vibration and thermal buckling behavior of moderately thick functionally graded structures in thermal environments. Compos. Struct. 2016, 157, 207–221. [Google Scholar] [CrossRef]
  4. Čanađija, M.; Barretta, R.; Sciarra, F.M.D. On functionally graded timoshenko nonisothermal nanobeams. Compos. Struct. 2016, 135, 286–296. [Google Scholar] [CrossRef]
  5. Barretta, R.; Čanađija, M.; Feo, L.; Lucianod, R.; de Sciarraa, F.M.; Pennac, R. Exact solutions of inflected functionally graded nano-beams in integral elasticity. Compos. Part B 2018, 142, 273–286. [Google Scholar] [CrossRef]
  6. Rostami, J.; Pwt, T.; Fang, Z. Sparse and dispersion-based matching pursuit for minimizing the dispersion effect occurring when using guided wave for pipe inspection. Materials 2017, 10, 622. [Google Scholar] [CrossRef] [PubMed]
  7. Chen, J.Y.; Guo, J.H.; Pan, E.N. Wave propagation in magneto-electro-elastic multilayered plates with nonlocal effect. J. Sound Vib. 2017, 400, 550–563. [Google Scholar] [CrossRef]
  8. Wang, L.G.; Rokhlin, S.I. Recursive geometric integrators for wave propagation in a functionally graded multilayered elastic medium. J. Mech. Phys. Solids 2004, 52, 2473–2506. [Google Scholar] [CrossRef]
  9. Pan, E.; Han, F. Exact solution for functionally graded and layered magneto-electro-elastic plates. Int. J. Eng. Sci. 2005, 43, 321–339. [Google Scholar] [CrossRef]
  10. Bhangale, R.K.; Ganesan, N. Free vibration of simply supported functionally graded and layered magneto-electro-elastic plates by finite element method. J. Sound Vib. 2006, 294, 1016–1038. [Google Scholar] [CrossRef]
  11. Wu, B.; Yu, J.G.; He, C.F. Wave propagation in non-homogeneous magneto-electro-elastic plates. J. Sound Vib. 2008, 317, 250–264. [Google Scholar]
  12. Cao, X.S.; Shi, J.P.; Jin, F. Lamb wave propagation in the functionally graded piezoelectric-piezomagnetic material plate. Acta Mech. 2012, 223, 1081–1091. [Google Scholar] [CrossRef]
  13. Singh, B.M.; Rokne, J. Propagation of SH waves in layered functionally gradient piezoelectric-piezomagnetic structures. Philos. Mag. 2013, 93, 1690–1700. [Google Scholar] [CrossRef]
  14. Xiao, D.L.; Han, Q.; Liu, Y.J.; Li, C. Guided wave propagation in an infinite functionally graded magneto-electro-elastic plate by the Chebyshev spectral element method. Compos. Struct. 2016, 153, 704–711. [Google Scholar] [CrossRef]
  15. Roshchupkin, D.; Ortega, L.; Snigirev, A.; Snigirev, I. X-ray imaging of the surface acoustic wave propagation in La3Ga5SiO14 crystal. Appl. Phys. Lett. 2013, 103, 154101. [Google Scholar] [CrossRef]
  16. Glushkov, E.; Glushkova, N.; Zhang, C.Z. Surface and pseudo-surface acoustic waves piezoelectrically excited in diamond-based structures. J. Appl. Phys. 2012, 112, 064911–064920. [Google Scholar] [CrossRef]
  17. Auld, B.A. Acoustic Fields and Waves in Solids, 2nd ed.; Krieger Publishing Company: Malabar, FL, USA, 1990. [Google Scholar]
  18. Lyon, R.H. Response of an elastic plate to localized driving forces. J. Acoust. Soc. Am. 1955, 27, 259–265. [Google Scholar] [CrossRef]
  19. Mindlin, R.D.; Medick, M.A. Extensional vibrations of elastic plates. J. Appl. Mech. 1959, 26, 561–569. [Google Scholar]
  20. Freedman, A. The vibration, with the Poisson ratio, of Lamb modes in a free plate, I: General spectra. J. Sound Vib. 1990, 137, 209–230. [Google Scholar] [CrossRef]
  21. Quintanilla, F.H.; Lowe, M.J.S.; Craster, R.V. Full 3D dispersion curve solutions for guided waves in generally anisotropic media. J. Sound Vib. 2016, 363, 545–559. [Google Scholar] [CrossRef] [Green Version]
  22. Yan, X.; Yuan, F.G. Conversion of evanescent Lamb waves into propagating waves via a narrow aperture edge. J. Acoust. Soc. Am. 2015, 137, 3523–3533. [Google Scholar] [CrossRef] [PubMed]
  23. Yan, X.; Yuan, F.G. A semi-analytical approach for SH guided wave mode conversion from evanescent into propagating. Ultrasonics 2018, 84, 430–437. [Google Scholar] [CrossRef] [PubMed]
  24. Chen, H.; Wang, J.; Du, J.K.; Yang, J. Propagation of shear-horizontal waves in piezoelectric plates of cubic crystals. Arch. Appl. Mec. 2016, 86, 517–528. [Google Scholar] [CrossRef]
  25. Datta, S.; Hunsinger, B.J. Analysis of surface waves using orthogonal functions. J. Appl. Phys. 1978, 49, 475–479. [Google Scholar] [CrossRef]
  26. Li, J.Y.; Dunn, M.L. Micromechanics of Magnetoelectroelastic Composite Materials: Average Fields and Effective Behavior. J. Intell. Mater. Syst. Struct. 1998, 9, 404–416. [Google Scholar] [CrossRef]
  27. Han, X.; Liu, G.R. Elastic waves in a functionally graded piezoelectric cylinder. Smart Mater. Struct. 2003, 12, 962–971. [Google Scholar] [CrossRef]
  28. Guo, Y.Q.; Chen, W.Q.; Zhang, Y.L. Guided wave propagation in multilayered piezoelectric structures. Sci. China Ser. G Phys. Mech. Astron. 2009, 52, 1094–1104. [Google Scholar] [CrossRef]
  29. Ezzin, H.; Amor, M.B.; Ghozlen, M.H.B. Propagation behavior of SH waves in layered piezoelectric/piezomagnetic plates. Acta Mech. 2017, 228, 1071–1081. [Google Scholar] [CrossRef]
Figure 1. Geometry of the problem.
Figure 1. Geometry of the problem.
Materials 11 01186 g001
Figure 2. Dispersion curves of Lamb wave in a steel plate; hollow dots-our results, solid lines-literature results from the spectral collocation method; real branch in blue, purely imaginary branch in black, complex branch in red.
Figure 2. Dispersion curves of Lamb wave in a steel plate; hollow dots-our results, solid lines-literature results from the spectral collocation method; real branch in blue, purely imaginary branch in black, complex branch in red.
Materials 11 01186 g002
Figure 3. Dispersion curves of Lamb wave in a PZT-4 plate; red dots—our results, black dotted lines—literature results from the reverberation-ray matrix method.
Figure 3. Dispersion curves of Lamb wave in a PZT-4 plate; red dots—our results, black dotted lines—literature results from the reverberation-ray matrix method.
Materials 11 01186 g003
Figure 4. Dispersion curves of propagating Lamb-like wave in a linear functionally graded piezoelectric-piezomagnetic material (FGPPM) plate with various “M”.
Figure 4. Dispersion curves of propagating Lamb-like wave in a linear functionally graded piezoelectric-piezomagnetic material (FGPPM) plate with various “M”.
Materials 11 01186 g004
Figure 5. 3D dispersion curves of Lamb wave: (a) four quadrants, (b) one quadrant; blue—real solutions, green—imaginary solutions, red—complex solutions.
Figure 5. 3D dispersion curves of Lamb wave: (a) four quadrants, (b) one quadrant; blue—real solutions, green—imaginary solutions, red—complex solutions.
Materials 11 01186 g005
Figure 6. The phase velocity dispersion and attenuation curves; propagating wave in blue, evanescent wave in red.
Figure 6. The phase velocity dispersion and attenuation curves; propagating wave in blue, evanescent wave in red.
Materials 11 01186 g006
Figure 7. 3D dispersion curves of Lamb wave: (a) cubic graded field, (b) sinusoidal graded field. blue—real branches, green—purely imaginary branches, red—complex branches.
Figure 7. 3D dispersion curves of Lamb wave: (a) cubic graded field, (b) sinusoidal graded field. blue—real branches, green—purely imaginary branches, red—complex branches.
Materials 11 01186 g007
Figure 8. Dispersion curves of propagating Lamb wave for FGPPM plates with different graded fields; (a) Frequency spectra; (b) Phase velocity spectra.
Figure 8. Dispersion curves of propagating Lamb wave for FGPPM plates with different graded fields; (a) Frequency spectra; (b) Phase velocity spectra.
Materials 11 01186 g008
Figure 9. Variation curves of the three graded functions.
Figure 9. Variation curves of the three graded functions.
Materials 11 01186 g009
Figure 10. Effect of graded field on magneto-electromechanical coupling factor for S0 mode.
Figure 10. Effect of graded field on magneto-electromechanical coupling factor for S0 mode.
Materials 11 01186 g010
Figure 11. Distributions of the physical quantities when Ω = 1.01115, Ψ = 0.23216 − 0.04612i. (a) displacement distribution, (b) electric potential and magnetic potential distribution.
Figure 11. Distributions of the physical quantities when Ω = 1.01115, Ψ = 0.23216 − 0.04612i. (a) displacement distribution, (b) electric potential and magnetic potential distribution.
Materials 11 01186 g011aMaterials 11 01186 g011b
Figure 12. Distributions of the physical quantities when Ω = 1.01911, Ψ = 0.17938. (a) displacement distribution, (b) electric potential and magnetic potential distribution.
Figure 12. Distributions of the physical quantities when Ω = 1.01911, Ψ = 0.17938. (a) displacement distribution, (b) electric potential and magnetic potential distribution.
Materials 11 01186 g012
Table 1. Material parameters of two piezoelectric-piezomagnetic materials.
Table 1. Material parameters of two piezoelectric-piezomagnetic materials.
MaterialsProperty
C11C12C13C22C23C33C44C55C66
Ba2TiO3166777816678162434344.6
CoFe2O428617317028617026945.345.346.5
Property
e15e24e31e32e33112233ρ
Ba2TiO311.611.6−4.4−4.318.6196201285.8
CoFe2O4000000.80.80.935.3
Property
q15q24q31q32q33μ11μ22μ33
Ba2TiO3000005510
CoFe2O4550550580.3580.3699.7−590−590157
Units: Cij (109 N/m2), ij (10−10 F/m), eij (C/m2), qij (N/Am), μij (10−6 Ns2/C2), ρ (103 kg/m3).
Table 2. Material parameters.
Table 2. Material parameters.
MaterialProperty
C11C12C13C22C23C33C44C55C66
PZT-413978741397411525.625.630.5
Property
e15e24e31e32e33112233ρ
PZT-412.712.7−5.2−5.215.16565567.5
Units: Cij (109 N/m2), ij (10−10 F/m), eij (C/m2), ρ (103 kg/m3).
Table 3. Convergence of complex wave numbers of the first six modes (Ω = 0.1).
Table 3. Convergence of complex wave numbers of the first six modes (Ω = 0.1).
MMode
123456
100.31308
+0.67653i
0.32055
+1.23591i ™
0.16593
+1.87713i ™
0.33292
+2.58928i ™
0.34536
+3.14395i ™
0.24222
+3.74226i ™
110.31308
+0.67653i
0.32055
+1.23591i ™
0.16593
+1.87713i ™
0.33184
+2.59041i
0.34722
+3.14407i ™
0.24481
+3.73460i ™
120.31308
+0.67653i
0.32055
+1.23591i ™
0.16593
+1.87713i ™
0.33063
+2.59061i
0.35203
+3.14418i ™
0.25047
+3.72524i
130.31308
+0.67653i
0.32055
+1.23591i ™
0.16593
+1.87713i ™
0.33063
+2.59061i
0.35211
+3.14436i ™
0.25504
+3.72496i ™
140.31308
+0.67653i
0.32055
+1.23591i ™
0.16593
+1.87713i ™
0.33063
+2.59061i
0.35211
+3.14436i ™
0.25561
+3.72471i ™
200.31308
+0.67653i
0.32055
+1.23591i ™
0.16593
+1.87713i ™
0.33063
+2.59061i
0.35211
+3.14436i ™
0.25561
+3.72471i ™

Share and Cite

MDPI and ACS Style

Zhang, X.; Li, Z.; Yu, J. The Computation of Complex Dispersion and Properties of Evanescent Lamb Wave in Functionally Graded Piezoelectric-Piezomagnetic Plates. Materials 2018, 11, 1186. https://doi.org/10.3390/ma11071186

AMA Style

Zhang X, Li Z, Yu J. The Computation of Complex Dispersion and Properties of Evanescent Lamb Wave in Functionally Graded Piezoelectric-Piezomagnetic Plates. Materials. 2018; 11(7):1186. https://doi.org/10.3390/ma11071186

Chicago/Turabian Style

Zhang, Xiaoming, Zhi Li, and Jiangong Yu. 2018. "The Computation of Complex Dispersion and Properties of Evanescent Lamb Wave in Functionally Graded Piezoelectric-Piezomagnetic Plates" Materials 11, no. 7: 1186. https://doi.org/10.3390/ma11071186

APA Style

Zhang, X., Li, Z., & Yu, J. (2018). The Computation of Complex Dispersion and Properties of Evanescent Lamb Wave in Functionally Graded Piezoelectric-Piezomagnetic Plates. Materials, 11(7), 1186. https://doi.org/10.3390/ma11071186

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop