Home Technology Nonlocal state-space strain gradient wave propagation of magneto thermo piezoelectric functionally graded nanobeam
Article Open Access

Nonlocal state-space strain gradient wave propagation of magneto thermo piezoelectric functionally graded nanobeam

  • Rajendran Selvamani EMAIL logo , Rubine Loganathan , Rossana Dimitri and Francesco Tornabene
Published/Copyright: April 3, 2023

Abstract

In this work, the state -space nonlocal strain gradient theory is used for the vibration analysis of magneto thermo piezoelectric functionally graded material (FGM) nanobeam. An analysis of FGM constituent properties is stated by using the power law relations. The refined higher order beam theory and Hamilton’s principle have been used to obtain the motion equations. Besides, the governing equations of the magneto thermo piezoelectric nanobeam are extracted by developed nonlocal state-space theory. And to solve the wave propagation problems, the analytical wave dispersion method is used. The effect of magnetic potential, temperature gradient, and electric voltage in variant parameters are presented in graph.

1 Introduction

Functionally graded materials (FGMs) are a type of composites initiated by a group of Japanese scientists to control the volume fraction of the mixture of two or more materials. The nonlinear vibration of the piezoelectric nanobeams based on the nonlocal and Timoshenko theory, the influence of the nonlocal parameter, temperature change, and external electric voltage on the size dependent nonlinear vibration characteristics of the piezo electric nanobeam are exposed [1]. Researchers [2] studied the natural frequencies along with mechanical and thermo electric vibration of piezoelectric nanobeams based on the nonlocal theory. Ebrahimi [3] reported the scattering of waves of FG nanobeam of viscoelastic nature. In the framework of third-order shear deformation theory [4], the vibration characteristics of functionally graded (METE-FG) nanobeams were analyzed. And the free vibrations of FG nano plates resting on elastic foundation via Hamilton principle was dealt in detail [5]. Alibeigi et al. [6] introduced the buckling retaliation of nanobeams on the basis of the Euler–Bernoulli beam model with the von Kármán geometrical nonlinearity. Bending of flexo electric magneto-electro-elastic (MEE) nanobeams lying over Winkler–Pasternak according to nonlocal elasticity theory has been studied [7]. Several studies were conducted on [8,9,10] hygro-thermal loading, the bending analysis of magneto-electro piezoelectric nanobeams system, dynamic analysis of smart nanostructures, and frequency analysis of thermally post buckled FGM thin beams. Stress-driven vs strain-driven elastic nanobeams have been discussed via integral elasticity [11,12]. Using the kinematic model, Kiani and Eslami [13] reported the buckling of beams made of FG under different types of thermal loading. The propagation of wave of infinite functionally graded plate in thermal environment was reported by Sun and Luo [14]. A consistent refined HSDT is designed to probe the free vibration of GF plates on elastic foundation and the influence of boundary condition on the natural frequency [15]. The different working conditions of the nano sized elements were studied by Thai et al. [16]. By considering the nonlocal elasticity [17], a prediction has been made that the essential behaviors of the nanostructures cannot be same as the macro scale structures. The Euler–Bernoulli beam theory was used to study the bending analysis of microtubules (MT) by Eringen [18]. Based on Euler–Bernoulli beam theory, the bending analysis of MT was studied using the method of Differential Quadrature (DQ) by Civalek and Demir [19]. Using finite element method (FEM) [20], the nonlinear bending in nanobeams were discussed. Reddy and El-Borgi [21] investigated the dispersion of waves with the effects of surface stresses in smart piezoelectric nanoplates. Bi-directional FGM nanobeams with the characteristics of bending, buckling, and vibrational nonlocal elements were concentrated in some of the previous studies [22,23,24]. The natural frequency variation located on a viscoelastic sheet was surveyed by using the nonlocal theory [25]. The size-dependent elements of beam were analyzed by Ebrahimi and Barati [26]. Nonlocal elasticity and its running conditions are discussed in details in the literature [27]. Based on the nonlocal strain gradient theory (NSGT) [28], the thermo-mechanical buckling problem of graphene sheets was proposed. Stiffness, softening, and hardening effect of FG beam were studied by Li et al. [29]. Solving the wave dispersion problem of nanoplates was accomplished by Ebrahimi and Dabbagh [30] with the application of infusing NSGT and surface-related elasticity for responsive piezoelectric materials. With the small-scale effect, the free vibration of 3D FGM Euler–Bernoulli nanobeam was studied by Hadi et al. [31]. Alibeigi et al. [32] exposed the buckling response of a nanobeam on the basis of the Euler–Bernoulli beam model using a couple stress theory under various types of thermal loading and an electrical and magnetic field. Timoshenko beam theory was investigated by Ke and Wang [33] with the rise in uniform temperature, magnetic potential, and external electric potential via nonlocal form to MEE vibrations. Bending of MEE nanobeam was studied in detail by Ebrahimi et al. [34]. Along with that, Ebrahimi et al. [35] investigated the bending of MEE nanobeams relating the nonlocal elasticity theory under hygro-thermal loading embedded in Winkler–Pasternak foundation. The size dependent problems using nonlocal elasticity theory, nonlocal couple stress theory, and shear deformation theory were reported [36,37]. Ebrahimi et al. [38] discussed the effects of various parameters on the wave dispersion characteristics of size-dependent nanoplates. The thermal effects on the buckling and free vibration of the FG nanobeams is documented well in the literature [39]. Ebrahimi and Barati [40] discussed the damping vibration characteristics of the hygro-thermally affected FG viscoelastic nanobeams. The thermal effect on buckling and free vibration characteristics of size-dependent Timoshenko nanobeams, and the free vibration of curved FG nano size beam in thermal environment been discussed in the literature [41,42]. The buckling and vibration properties of sandwich FG beams were studied by Vo et al. [43]. Jalaei et al. [44] studied the thermal and magnetic effects on the FG Timoshenko nanobeam. Studies over the hygro-thermal wave characteristic of nanobeam of an inhomogeneous material with porosity under magnetic field is notable [45].

Hence, this work shows the wave propagation analysis of FG nanobeam with the help of a nonlocal state-space strain gradient viscoelasticity. The magneto thermo material properties of the nanobeam also graded and implemented via power law relations and the motion equations are deduced through the Hamilton’s principle. Furthermore, the dispersion for computed external electric voltage, magnetic effect, and the gradient of temperature are presented with graphical solution.

2 Mathematical formulations

Based on the state-space nonlocal strain gradient theory, the length L, width b, and thickness h of a viscoelastic FG nanobeam has been investigated (Figure 1). The two parts of constituent FGM are composed of ceramic part and metallic part. The components of FGM are considered to be temperature dependent to evolve a realistic viscoelastic study.

Figure 1 
               Vibration analysis of FG nanobeam.
Figure 1

Vibration analysis of FG nanobeam.

In this section, power-law relations have been used to compute the properties. To calculate the temperature variable towards the thickness direction, the volume fraction of each phase must be calculated by using the power-law model. Hence, the volume fraction of ceramic partis given as follows:

(2.1) V c = z h + 1 2 p ,

where the thickness h and the exponential power law p probe each phase in the material with its distribution of properties. By considering any desired material property at the local temperature

(2.2) P = P 0 ( P 1 T 1 + 1 + P 1 T + P 2 T 2 + P 3 T 3 ) ,

where P 0 , P 1 , P 1 , P 2 , P 3 are the coefficients of material phases. The volume fraction of the metallic phase gives the volume fraction of the ceramic phase by V m = 1 V c .

According to Eringen’s nonlocal theory, the stress state at a point inside a body is a function of the strains at all points in the neighboring regions. The basic equations with zero body force can be defined as follows:

(2.3) σ i j = v α ( y y , τ ) [ C i j k l ε k l ( y ) e m i j E m ( y ) Ω n ( y ) C i j k l α k l Δ T ] d V ( y ) D i = v α ( y y , τ ) [ e i k l ε k l ( y ) E m ( y ) + Ω n ( y ) Δ T ] d V ( y ) B i = v α ( y y , τ ) [ ε k l ( y ) + E m ( x ) + χ n i Ω n ( y ) λ i Δ T ] d V ( y ) ,

where σ i j , ε i j , D i , E i , respectively, represent the stress, strain, electric displacement, and electric field and B i , Ω i are the magnetic induction and magnetic field. α k l and Δ T stand for thermal expansion and temperature difference. C i j k l , e m i j , and χ i j are the elastic, piezoelectric, and magnetic constants, respectively, and τ = e 0 a / l defines the scale coefficient, e 0 is the material constant, and a and l are the characteristic length in the internal and external sides.

The governing equations of the nanobeams are obtained by an accurate kinematic theory. Higher order shear deformation theory also reveals stress-strain changes in solid bodies. From the previous study [8], we can take the refined shear deformable beam’s displacement as follows:

(2.4) x ( x , z , t ˜ ) = ( x , t ˜ ) w b ( x , t ˜ ) x z w s ( x , t ˜ ) x f ( z ) z ( x , z , t ˜ ) = w b ( x , t ˜ ) + w s ( x , t ˜ ) ,

where Π , w b , and w s are the longitudinal displacement, bending, and shear components of the transverse displacement. Furthermore, in order to distribute the shear strain, f ( z ) is the shape function which is designed as follows [4]:

(2.5) f ( z ) = h e z h 2 + π 2 π sin π z h + h cos π z h h 2 h 2 + π 2 ,

To capture the shear strain and stress, the deformed structural cross section is uncertain with this function. At free surfaces, it is required to satisfy the assumption of shear strain nonexistence. By continuum infinitesimal strain tensor, the nonzero strains can be measured as follows:

(2.6) ε x x = Π x z 2 w b x 2 f ( z ) 2 w s x 2 γ x z = g ( z ) w s x ,

where g ( z ) = 1 d f ( z ) d z .

2.1 Motion equations

In accordance with Hamilton’s principle, the extended Lagrangian can be given as follows:

(2.7) L = T + V ,

So, the Hamilton’s principle can be given as follows:

(2.8) δ t ˜ 1 t ˜ 2 ( ̇ ̇ T ̇ + V ̇ ) d t ˜ = 0 ,

In Eq. (2.8), variables , V , and T are the strain energy, work done, and kinetic energy, respectively. Hence, the virtual strain energy can be given as follows:

(2.9) δ ̇ = ( σ i j δ ε i j d D x E x D z E z ) d z d x ,

(2.10) ϕ ( x , z , t ˜ ) = cos ( β z ) ϕ ( x , t ˜ ) + 2 z V 0 h e i k t ˜ ,

where D x = e 15 γ x z + e 11 E x and D z = e 31 ε x x + e 33 E z . The variation in electric potential in the x direction is β = π / h ; ϕ ( x , t ˜ ) ; V 0 and Ω are the external electric voltage and natural frequency of the piezoelectric nanobeam, respectively.

By infusing Eq. (2.6) in Eq. (2.9),

(2.11) δ Π = 0 L N δ Π x M b 2 δ w b x 2 M s 2 δ w s x 2 + Q δ w s x D x E x D z E z d z d x ,

and the stress resultants can be obtained as follows:

(2.12) [ N , M b , M s ] = A [ 1 , z , f ( z ) ] σ x x d A ,

(2.13) Q = A g ( z ) σ x z d A ,

The kinetic energy of the system can be determined as follows:

(2.14) δ T = ρ ( z ) Π x t ˜ δ Π x t ˜ + Π z t ˜ δ Π z t ˜ d V ,

Infusion of Eq. (2.4) in Eq. (2.14) results in the following:

(2.15) δ T = 0 L I 0 Π t ˜ δ Π t ˜ + ( w b + w s ) t ˜ δ ( w b + w s ) t ˜ I 1 Π t ˜ 2 δ w b x t ˜ + 2 w b x t ˜ δ Π t ˜ J 1 Π t ˜ 2 δ w s x t ˜ + 2 w s x t ˜ δ Π t ˜ + I 2 2 w b x t ˜ 2 δ w b x t ˜ + K 2 2 w s x t ˜ 2 δ w s x t ˜ + J 2 2 w b x t ˜ 2 δ w s x t ˜ + 2 w s x t ˜ 2 δ w b x t ˜ ,

In accordance with the magnetic and temperature effect,

f B = η A Ω x 2 2 w x 2 N T = h / 2 h / 2 E ( z ) λ ( z ) Δ T d z ,

where f B , η , A, and Ω x stand for the magnetic force, magnetic field permeability, cross sectional area of the nanobeam, and the magnetic potential of the longitudinal magnetic field. For an FG nanobeam, it is assumed that the temperature can be distributed uniformly across its thickness and the temperature gradient at stress free state is Δ T .

In the above definition, the inertia of mass moments can be defined as follows:

[ I 0 , I 1 , I 2 , J 1 , J 2 , K 2 ] = A [ 1 , z , z 2 , f ( z ) , z f ( z ) , f 2 ( z ) ] ρ ( z ) d A

and the work done N x with a temperature gradient of thermal effect can be defined as follows:

(2.16) δ V = 1 2 0 l N x + N T w x δ w x + η Ω x 2 2 w x 2 d x ,

By inserting the Eqs. (2.11) and (2.15) in Eq. (2.8), the equation of the beam in Euler–Lagrange can be derived and the outcome can be coupled as follows:

(2.17) N x = I 0 2 Π t ˜ 2 I 1 3 w b x t ˜ 2 J 1 3 w s x t ˜ 2 ,

(2.18) 2 M b x 2 = I 0 2 ( w b + w s ) t ˜ 2 + I 1 3 Π x t ˜ 2 I 2 4 w b x 2 t ˜ 2 J 2 4 w s x 2 t ˜ 2 N x + N T 2 w x 2 + η Ω x 2 2 w x 2 ,

(2.19) 2 M b x 2 + Q x N x 2 w x 2 = I 0 2 ( w b + w s ) t ˜ 2 + J 1 3 Π x t ˜ 2 J 2 4 w b x 2 t ˜ 2 K 2 4 w s x 2 t ˜ 2 N x + N T 2 w x 2 + η Ω x 2 2 w x 2 ,

(2.20) δ ϕ = h 2 h 2 cos ( β z ) D x x + β sin ( β z ) D z .

3 Nonlocal state-space model

This section demonstrates the nonlocality stress and strain effects on the time–space domains, when the wave length or excitation frequency interferes with time and intrinsic characteristic length. Nonlocal time–space viscoelasticity problems are based on the combination of the Boltzmann superposition integral and the Eringen concept of nonlocal elasticity. Accordingly, integral stress in nonlinear state and strain equations are stated as follows:

(3.1) t v ( K σ ( t ˜ τ , | r r | ) σ i j ( r , τ ) d r d τ ) = t v K ε ( t ˜ τ , | r r | ) C i j k l ε k l ( r , τ ) d r d τ ,

where σ i j and ε k l are the stress and strain tensor arrays and the nonlocal kernel functions are K σ ( t ˜ τ , | r r | ) and K ε ( t ˜ τ , | r r | ) . Eq. (3.1) can be read in the following form via Fourier, inverse Fourier, and Taylor series,

(3.2) 1 l σ 2 2 + τ σ t ˜ σ i j = C i j k l 1 l ε 2 2 + τ ε t ˜ ε k l ,

To balance the absence of stiffness-hardening behavior, the nonlocal strain gradient elasticity must be incorporated in the equation. The following relation can be used to derive the nonlocal strain gradient viscoelasticity with fraction.

(3.3) 1 l σ 2 2 + τ σ t ˜ σ i j = C i j k l 1 l ε 2 2 + τ ε t ˜ ε k l ,

The Kelvin–Voigt relation of viscoelastic material with three parameters in a solid state is given as follows:

(3.4) ( 1 μ 2 2 ) σ i j = C i j k l 1 λ 2 2 + τ t ˜ ε k l ,

where λ = l ε and μ = l σ are length scale and nonlocal parameters, respectively. The relation between rheological character and spatial nonlocality of the system is designed as τ = μ τ c E c . From Eqs. (2.12) and (2.13), the displacement field having resultant stress over the area of cross section of the beam can be computed as follows:

(3.5) ( 1 μ 2 2 ) N = 1 λ 2 2 + τ t ˜ × A x x u x B x x 2 w b x 2 B x x s 2 w s x 2 e 31 E z ,

(3.6) ( 1 μ 2 2 ) M b = 1 λ 2 2 + τ t ˜ × B x x u x D x x 2 w b x 2 D x x s 2 w s x 2 ( N x + N T η Ω x 2 ) 2 w x 2 e 31 E z ,

(3.7) ( 1 μ 2 2 ) M s = 1 λ 2 2 + τ t ˜ × B x x s u x D x x s 2 w b x 2 H x x s 2 w s x 2 ( N x + N T η Ω x 2 ) 2 w x 2 e 31 E z ,

(3.8) ( 1 μ 2 2 ) Q x z = 1 λ 2 2 + τ t ˜ × A x x s w s x ( N x + N T η Ω x 2 ) 2 w x 2 e 15 E x ,

(3.9) ( 1 μ 2 2 ) D x = e 15 γ x z + e 11 E x ,

(3.10) ( 1 μ 2 2 ) D z = e 31 ε x x + e 33 E z ,

where cross sectional rigidities are

(3.11) [ A x x , B x x , B x x s , D x x , D x x s , H x x s ] = A [ 1 , z , f ( z ) , z 2 , z f ( z ) , f 2 ( z ) ] ,

(3.12) A s = A g 2 ( z ) G ( z ) d A ,

where e 31 E x = e 31 ( cos ( β z ) ) ϕ x and e 15 E z = e 15 β sin ( β z ) ϕ 2 V 0 e i k t h .

4 Governing equations

Eqs. (3.5)–(3.10) must be substituted in equations (2.17)–(2.20). Now the governing equations are as follows:

(4.1) 1 λ 2 2 + τ t ˜ A x x 2 Π x 2 B x x 3 w b x 3 B x x s 3 w s x 3 + ( 1 μ 2 2 ) I 0 Π ̈ + I 1 w ̈ b x + J 1 w ̈ s x N x + N T η Ω x 2 2 w x 2 e 31 E z = 0 ,

(4.2) 1 λ 2 2 + τ t ˜ B x x 3 Π x 3 D x x 4 w b x 4 D x x s 4 w s x 4 + ( 1 μ 2 2 ) I 0 ( w ̈ b + w ̈ s ) I 1 Π ̈ x + I 2 2 w ̈ b x 2 + J 2 2 w ̈ s x 2 N x + N T η Ω x 2 2 w x 2 e 31 E z = 0 ,

(4.3) 1 λ 2 2 + τ t ˜ B x x s 3 Π x 3 D x x s 4 w b x 4 H x x s 4 w s x 4 + A s 2 w s x 2 + ( 1 μ 2 2 ) I 0 ( w ̈ b + w ̈ s ) J 1 Π ̈ x + J 2 2 w ̈ b x 2 + K 2 2 w ̈ s x 2 N x + N T η Ω x 2 2 w x 2 e 31 E z e 15 E x = 0 .

5 Analytical solution

The governing equation obtained in Section 4 will be solved in this section. The analytical wave dispersion method is used in this case to solve the problems of wave propagation for various types of structures, including beams, plates, and shells. A higher order beam’s solution function can be in the form as follows:

(5.1) Π w b w s = Π exp [ i ( β x ω t ˜ ) ] w b exp [ i ( β x ω t ˜ ) ] w s exp [ i ( β x ω t ˜ ) ] ,

where Π , w b , and w s are the anonymous amplitudes of propagating waves. Here ω & β are frequency and wave number, respectively. By substituting the above expression in Eqs. (4.1)–(4.3), the achieved form will be as follows:

(5.2) { [ K ] + [ C ] ω + [ M ] ω 2 } Π w b w s = 0 ,

where damping and mass matrices are given by [ K ] , [ C ] , and [ M ] respectively. The components of these symmetric matrices are as follows:

(5.3) k 11 = ( 1 + λ 2 β 2 ) A x x β 2 k 12 = i ( N x + η Ω x 2 ) ( 1 + λ 2 β 2 ) B x x β 5 k 13 = i ( N x + η Ω x 2 ) ( 1 + λ 2 β 2 ) B x x s β 5 k 14 = e 31 ( 1 λ 2 β 2 ) β sin ( β z ) ϕ 2 V 0 e i k t h k 22 = ( N x + η Ω x 2 ) N T ( 1 + λ 2 β 2 ) D x x β 6 k 23 = ( N x + η Ω x 2 ) N T ( 1 + λ 2 β 2 ) D x x s β 6 k 24 = e 31 ( 1 λ 2 β 2 ) β sin ( β z ) ϕ 2 V 0 e i k t h k 33 = ( N x + η Ω x 2 ) N T ( 1 + λ 2 β 2 ) ( H x x s β 6 + A s β 4 ) k 34 = e 31 ( 1 λ 2 β 2 ) β 2 sin ( β z ) ϕ 2 V 0 e i k t h ,

(5.4) c 11 = A x x τ β 2 c 12 = i ( N x + η Ω x 2 ) B x x τ β 5 c 13 = i ( N x + η Ω x 2 ) B x x s τ β 5 c 22 = ( N x + η Ω x 2 ) D x x τ β 6 c 23 = ( N x + η Ω x 2 ) D x x s τ β 6 c 24 = τ e 31 β 2 cos ( β z ) d z d t ϕ 2 V 0 e i k t h ( i k ) c 33 = ( N x + η Ω x 2 ) τ β 4 ( H x x s β 2 + A s ) c 34 = τ e 31 β 2 cos ( β z ) d z d t ϕ 2 V 0 e i k t h ( i k ) ,

(5.5) m 11 = ( 1 + μ 2 β 2 ) I 0 m 12 = i ( N x + η Ω x 2 ) β 2 ( 1 + μ 2 β 2 ) I 1 m 13 = i ( N x + η Ω x 2 ) β 2 ( 1 + μ 2 β 2 ) J 1 m 22 = ( 1 + μ 2 β 2 ) ( I 0 + I 2 ( N x + η Ω x 2 ) β 4 ) m 23 = ( 1 + μ 2 β 2 ) ( I 0 + J 2 ( N x + η Ω x 2 ) β 4 ) m 33 = ( 1 + μ 2 β 2 ) ( I 0 + K 2 ( N x + η Ω x 2 ) β 2 ) m 34 = ( 1 + μ 2 β 2 ) e 31 β 2 cos ( β z ) d z d t ϕ 2 V 0 e i k t h ( i k ) .

6 Results and discussion

This section illustrates the magneto-thermo vibration of FG nanobeam with the numerical examples. The material properties are composed of BaTiO 3 and CoFe 2 O 4 , which are presented in Table 1. Table 2 shows the comparison of the buckling load for various power-law exponent values. Whereas, Table 3 gives the variation in the magnetic potential and electric voltages with the varying and increasing nonlocal parameter. Since an increasing nonlocal parameter shows a decrement with the magnetic potential ( Ω ) and the electric voltage ( V ).

Table 1

MEE coefficients of material properties

Material Properties
BaTiO 3 E (Pa) 166 e 33 (N/m2 K) (7.124 × 10−9)
ρ (kg/m3) 5,800 e 15 (c/m2) 14.1
ν (–) 1.1945
CoFe 2 O 4 E (Pa) 286 e 31 (c/m2) −(4.1)
ρ (kg/m3) 5,300 e 11 (c/V m) (5.841 × 10−9)
ν (–) 1.167
Table 2

Comparison of FGM beam non-dimensional buckling for various power-law exponents

L/h p = 0 p = 0.5 p = 1 p = 2 p = 5 p = 10
5 Nguyen et al. (2015) 48.8406 32.0013 24.6894 19.1577 15.7355 14.1448
Present 48.835 31.967 24.6870 19.1605 15.7401 14.13
10 Nguyen et al. (2015) 52.3083 34.0002 26.1707 20.3909 17.1091 15.5278
Present 52.3082 34.0087 26.1727 20.3936 17.1118 15.5291
Table 3

Dimensionless frequency of an FG nanobeam varies with nonlocal parameters, electric voltages, and magnetic potentials

µ p = 0.2 p = 1 p = 5
V = −5 V = 0 V = +5 V = −5 V = 0 V = +5 V = −5 V = 0 V = +5
0 Ω = −0.05 8.34927 7.68726 7.24086 8.3637 7.69636 7.24398 8.3781 7.70545 7.24711
Ω = 0 8.34708 7.68034 7.22898 8.36151 7.68946 7.23212 8.37592 7.69856 7.23525
Ω = −0.05 8.34489 7.67343 7.21709 8.35933 7.68255 7.22023 8.37373 7.69166 7.22337
1 Ω = −0.05 7.96429 7.33365 6.9088 7.97941 7.34319 6.91208 7.9945 7.35273 6.91536
Ω = 0 7.96199 7.32641 6.89636 7.97712 7.33596 6.89964 7.99222 7.3455 6.90293
Ω = −0.05 7.95969 7.31915 6.88389 7.97483 7.32872 6.88718 7.98993 7.33827 6.89047
2 Ω = −0.05 7.62789 7.02471 6.61874 7.64368 7.03468 6.62216 7.65944 7.04462 6.62558
Ω = 0 7.6255 7.01715 6.60575 7.64129 7.02712 6.60917 7.65705 7.03708 6.6126
Ω = −0.05 7.6231 7.00958 6.59273 7.6389 7.01956 6.59617 7.65466 7.02953 6.5996
3 Ω = −0.05 7.33065 6.75176 6.3625 7.34708 6.76213 6.36606 7.36347 6.77248 6.36962
Ω = 0 7.32816 6.74389 6.34898 7.34459 6.75427 6.35255 7.36099 6.76463 6.35612
Ω = −0.05 7.32566 6.73601 6.33544 7.3421 6.7464 6.33902 7.3585 6.75678 6.34259
0 Ω = −0.05 8.55703 8.31260 8.37748 9.43205 8.84645 8.55259 10.2325 9.34987 8.72418
Ω = 0 8.41640 7.88376 7.67745 9.30465 8.44476 7.86815 10.1152 8.97075 8.05433
Ω = −0.05 8.27337 7.43021 6.90683 9.17548 8.02299 7.11820 9.99651 8.57489 7.32348
1 Ω = −0.05 8.08984 7.91790 8.03872 9.01035 8.47664 8.22104 9.84516 9.00077 8.39941
Ω = 0 7.94094 7.46642 7.30630 8.87690 8.05654 7.50644 9.72317 8.60629 7.70137
Ω = −0.05 7.78919 6.98583 6.49177 8.74141 7.61329 6.71622 9.59964 8.19284 6.93341
2 Ω = −0.05 7.67792 7.57252 7.74446 8.64241 8.15496 7.93354 9.50958 8.69849 8.11823
Ω = 0 7.52087 7.09911 6.98123 8.50319 7.71737 7.19042 9.38323 8.28964 7.39369
Ω = −0.05 7.36046 6.59180 6.12362 8.36165 7.25342 6.36107 9.25516 7.85955 6.58997
3 Ω = −0.05 7.31058 7.26690 7.48595 8.31777 7.87198 7.68141 9.21554 8.43377 7.87201
Ω = 0 7.14545 6.77216 6.69332 8.17302 7.41771 6.91123 9.08510 8.01142 7.12247
Ω = −0.05 6.97642 6.23831 5.79324 8.02566 6.93375 6.04369 8.95277 7.56553 6.28416

Now, Figures 2 and 3 highlight the effect of external voltage (V) with the variation in the rising temperature ( Δ T ) and the gradient index (p) for μ = 1.0 & 1.5 . Whereas, with the increase in the external voltage (V = 0, 0.5, and 1), the temperature gradually decreases with respect to the gradient index (p) and nonlocal parameter. Figure 4 shows the decrement in the buckling load with an increase in the magnetic effect ( Ω ) and the nonlocal value ( μ ). Also, when there is no magnetic effect found, the buckling load at some point decreases with the increse in the value of the nonlocal ( μ ). Figure 5 presents the variation in the dimensionless buckling load with the nonlocal parameter ( μ ) in the effect of the external voltage (V). Since there is an increase in the voltage (V), the dimensionless buckling load declines and gets neutralized at some point of the nonlocal parameter. Figures 6 and 7 depict the effect of gradient index with respect to the presence of the dimensionless buckling load and magnetic potential for different nonlocal values. Therefore, with the surge effect of the gradient index (p), there is an increase in the magnetic potential with the decrease in the dimensionless bucklingload and ampilifies over the nonlocal parameter ( μ ). Figures 8 and 9 represent the impact of the electric voltage and dimensionless buckling load in the presence of the gradient index p = 0, 1, 5 over nonlocal parameter ( μ ). Hence, when the gradient index is null, there is no fluctuation in the buckling load and forms a linear effect and nonlocal values ampilify the buckling load. So, if the gradient index rises, then the electric voltage increases with the decrease in buckling load gradually. Figures 10 and 11 interpret that temperature ( Δ T ) stabilizes with an increase in the gradient index and magnetic potential. Thus, when the magnetic potential Ω = 0 , the temperature calms down and stabilizes at some point with respect to an increase in the gradient index (p) and nonlocal values. Figures 12 and 13 show the effect of an internal damping factor ( τ ) in the existence of wave number and wave frequency over nonlocal parameter ( μ ). In these figures, the damping factor with τ 1.5 shows a raising linear value and when τ 0.5 , there is an oscillation in the wave number and wave frequency.

Figure 2 
               External voltage with the presence of rising temperature to gradient index for μ = 1.0.
Figure 2

External voltage with the presence of rising temperature to gradient index for μ = 1.0.

Figure 3 
               External voltage in the presence of temperature with respect to gradient index for μ = 1.5.
Figure 3

External voltage in the presence of temperature with respect to gradient index for μ = 1.5.

Figure 4 
               Magnetic potential on the dimensionless buckling with respect to nonlocal parameter.
Figure 4

Magnetic potential on the dimensionless buckling with respect to nonlocal parameter.

Figure 5 
               External voltage on the dimensionless buckling load with respect to nonlocal parameter.
Figure 5

External voltage on the dimensionless buckling load with respect to nonlocal parameter.

Figure 6 
               Gradient index with the presence of rising dimensionless buckling load to magnetic potential for μ = 1.0.
Figure 6

Gradient index with the presence of rising dimensionless buckling load to magnetic potential for μ = 1.0.

Figure 7 
               Gradient index in the presence of dimensionless buckling load to magnetic potential μ = 1.5.
Figure 7

Gradient index in the presence of dimensionless buckling load to magnetic potential μ = 1.5.

Figure 8 
               Gradient index in the presence of dimensionless buckling load to electric voltage when μ = 1.0.
Figure 8

Gradient index in the presence of dimensionless buckling load to electric voltage when μ = 1.0.

Figure 9 
               Gradient index in the presence of dimensionless buckling load to electric voltage when μ = 1.5.
Figure 9

Gradient index in the presence of dimensionless buckling load to electric voltage when μ = 1.5.

Figure 10 
               Magnetic potential in the presence of temperature with respect to gradient index for μ = 1.0.
Figure 10

Magnetic potential in the presence of temperature with respect to gradient index for μ = 1.0.

Figure 11 
               Magnetic potential in the presence of temperature with respect to gradient index for μ = 1.5.
Figure 11

Magnetic potential in the presence of temperature with respect to gradient index for μ = 1.5.

Figure 12 
               Damping factor in the presence of wave number with respect to wave frequency for μ = 1.0.
Figure 12

Damping factor in the presence of wave number with respect to wave frequency for μ = 1.0.

Figure 13 
               Damping factor in the presence of wave number with respect to wave frequency for μ = 1.5.
Figure 13

Damping factor in the presence of wave number with respect to wave frequency for μ = 1.5.

7 Conclusion

The above study shows the wave propagation analysis of piezoelectric FGM nanobeam. Magneto thermo properties of the FG nanobeam are considered to be the function of thickness according to the power-law model. The governing equations are extracted by substituting the structure displacement field equations in the beam’s Euler–Lagrange equations and are framed as symmetric matrices components to arrive at required solutions. Hence, the upshots of the work are as follows:

  • The stability behaviors of FGM nanobeam are affected .by magneto thermo piezo electricity and nonlocal values.

  • Physical variants could be controlled via applying a suitable value of damping factor.

  • Natural frequency reduces, while the nonlocal parameter and gradient index of the FG nanobeam amplify.

  • The increase in power-law index softens the volume fraction.

  • The bending rigidity and phase velocities are high in amplified wave numbers and get reversed in low wave numbers.

Acknowledgments

The authors would like to thank the reviewers for their useful comments and appropriate revisions of the manuscript.

  1. Funding information: The authors state no funding involved.

  2. Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

  3. Conflict of interest: Rossana Dimitri and Francesco Tornabene, who are the co-authors of this article, are current Editorial Board members of Curved and Layered Structures. This fact did not affect the peer-review process. The authors declare no other conflict of interest.

  4. Data availability statement: The datasets analysed during the current study are available from the corresponding author on reasonable request.

References

[1] Ke LL, Wang YS, Wang ZD. Nonlinear vibration of piezoelectric based on the nonlocal theory. Compos Struct. 2012;94:2038–47.10.1016/j.compstruct.2012.01.023Search in Google Scholar

[2] Ke LL, Wang YS. Thermoelectric-mechanical vibration of piezoelectric nanobeams based on the nonlocal theory. Smart Mater Struct. 2012;21:025018.10.1088/0964-1726/21/2/025018Search in Google Scholar

[3] Ebrahimi F. Wave dispersion in viscoelastic FG nanobeam via a novel spatial – temporal nonlocal strain gradient framework. Waves Random Complex Media. 2021. https://doi.org/10.1080/17455030.2021.1970282.10.1080/17455030.2021.1970282Search in Google Scholar

[4] Ebrahimi F, Barati MR. Vibration analysis of smart piezoelectrically actuated nanobeams subjected to magneto-electrical field in thermal environment. J Vib Control. 2018;24:549–64.Search in Google Scholar

[5] Zaoui FZ, Ouinas D, Tounsi A. New 2D and quasi-3D shear deformation theories for free vibration of functionally graded plates on elastic foundations. Compos B: Eng. 2019;159:231–47.10.1016/j.compositesb.2018.09.051Search in Google Scholar

[6] Alibeigi B, Tadi Beni Y, Mehralian F. On the thermal buckling of magneto-electro-elastic piezoelectric nanobeams. Eur Phys J Plus. 2018;133:133.Search in Google Scholar

[7] Shariati A, Ebrahimi F, Karimiasl M, Selvamani R, Toghroli A. On bending characteristics of smart magneto-electro-piezoelectric nanobeams system. Adv Nano Res. 2020;9:183–91.Search in Google Scholar

[8] Ebrahimi F, Karimiasl M, Selvamani R. Bending analysis of magneto-electro piezoelectric nanobeams system under hygro-thermal loading. Adv Nano Res. 2020;8:203–14.Search in Google Scholar

[9] Ebrahimi F, Dabbagh A. Wave propagation analysis of smart nanostructures. 1st ed. Boca Raton (FL): CRC Press; 2020.10.1201/9780429279225-1Search in Google Scholar

[10] Li SR, Su HD, Cheng CJ. Free vibration of functionally graded material beams with surface-bonded piezoelectric layers in thermal environment. Appl Math Mech. 2019;30:969–82.10.1007/s10483-009-0803-7Search in Google Scholar

[11] Romano G, Barretta R. Stress-driven versus strain-driven nonlocal integral model for elastic nano-beams. Compos B: Eng. 2017;114:184–8.10.1016/j.compositesb.2017.01.008Search in Google Scholar

[12] Barretta R, Fabbrocino F, Luciano R, de Sciarra FM. Closed-form solutions in stress-driven two-phase integral elasticity for bending of functionally graded nano-beams. Physica E Low Dimens Syst Nanostruct. 2018;97:13–30.10.1016/j.physe.2017.09.026Search in Google Scholar

[13] Kiani Y, Eslami MR. Thermal buckling analysis of functionally graded material beams. Int J Mech Mater Des. 2010; 6:229–38.10.1007/s10999-010-9132-4Search in Google Scholar

[14] Sun D, Luo SN. Wave propagation of functionally graded material plates in thermal environments. Ultrasonics. 2011;51:940–52.10.1016/j.ultras.2011.05.009Search in Google Scholar PubMed

[15] Thai HT, Choi DH. A refined shear deformation theory for free vibration of functionally graded plates on elastic foundation. Compos B: Eng. 2012;43:2335–47.10.1016/j.compositesb.2011.11.062Search in Google Scholar

[16] Thai HT, Park T, Choi DH. An efficient shear deformation theory for vibration of functionally graded plates. Arch Appl Mech. 2013;83:137–49.10.1007/s00419-012-0642-4Search in Google Scholar

[17] Shahsavari D, Shahsavari M, Li L, Karami B. A novel quasi-3D hyperbolic theory for free vibration of FG plates with porosities resting on Winkler/Pasternak/Kerr foundation. Aerosp Sci Technol. 2018;72:134–49.10.1016/j.ast.2017.11.004Search in Google Scholar

[18] Eringen AC. On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves. J Appl Phys. 1983;54:4703–10.10.1063/1.332803Search in Google Scholar

[19] Civalek Ö, Demir Ç. Bending analysis of microtubules using nonlocal Euler–Bernoulli beam theory. Appl Math Model. 2011;35:2053–67.10.1016/j.apm.2010.11.004Search in Google Scholar

[20] Arani AG, Amir S, Shajari AR, Mozdianfard MR. Electro-thermo-mechanical buckling of DWBNNTs embedded in bundle of CNTs using nonlocal piezoelasticity cylindrical shell theory. Compos B: Eng. 2012;43:195–203.10.1016/j.compositesb.2011.10.012Search in Google Scholar

[21] Reddy JN, El-Borgi S. Eringen’s nonlocal theories of beams accounting for moderate rotations. Int J Eng Sci. 2014;82:159–77.10.1016/j.ijengsci.2014.05.006Search in Google Scholar

[22] Zhang LL, Liu JX, Fang XQ, Nie GQ. Effects of surface piezoelectricity and nonlocal scale on wave propagation in piezoelectric nanoplates. Eur J Mech A Solids. 2014;46:22–9.10.1016/j.euromechsol.2014.01.005Search in Google Scholar

[23] Nejad MZ, Hadi A. Eringen’s non-local elasticity theory for bending analysis of bi-directional functionally graded Euler–Bernoulli nano-beams. Int J Eng Sci. 2016;106:1–9.10.1016/j.ijengsci.2016.05.005Search in Google Scholar

[24] Nejad MZ, Hadi A. Non-local analysis of free vibration of bi-directional functionally graded Euler–Bernoulli nano-beams. Int J Eng Sci. 2016;105:1–11.10.1016/j.ijengsci.2016.04.011Search in Google Scholar

[25] Nejad MZ, Hadi A, Rastgoo A. Buckling analysis of arbitrary two-directional functionally graded Euler–Bernoulli nano-beams based on nonlocal elasticity theory. Int J Eng Sci. 2016;103:1–10.10.1016/j.ijengsci.2016.03.001Search in Google Scholar

[26] Ebrahimi F, Barati MR. Size-dependent vibration analysis of viscoelastic nanocrystalline silicon nanobeams with porosities based on a higher order refined beam theory. Compos Struct. 2017;166:256–67.10.1016/j.compstruct.2017.01.036Search in Google Scholar

[27] Lim CW, Zhang G, Reddy JN. A higher-order nonlocal elasticity and strain gradient theory and its applications in wave propagation. J Mech Phys Solids. 2015;78:298–313.10.1016/j.jmps.2015.02.001Search in Google Scholar

[28] Farajpour A, Haeri Yazdi MR, Rastgoo A, Mohammadi M. A higher-order nonlocal strain gradient plate model for buckling of orthotropic nanoplates in thermal environment. Acta Mech. 2016;227:1849–67.10.1007/s00707-016-1605-6Search in Google Scholar

[29] Li L, Li X, Hu Y. Free vibration analysis of nonlocal strain gradient beams made of functionally graded material. Int J Eng Sci. 2016;102:77–92.10.1016/j.ijengsci.2016.02.010Search in Google Scholar

[30] Ebrahimi F, Dabbagh A. Wave propagation analysis of embedded nanoplates based on a nonlocal strain gradient-based surface piezoelectricity theory. Eur Phys J Plus. 2017;132:449.10.1140/epjp/i2017-11694-2Search in Google Scholar

[31] Hadi A, Zamani Nejad M, Hosseini M. Vibrations of three-dimensionally graded nanobeams. Int J Eng Sci. 2018;128:12–23.10.1016/j.ijengsci.2018.03.004Search in Google Scholar

[32] Alibeigi B, Tadi Beni Y, Mehralian F. On the thermal buckling of magneto-electro-elastic piezoelectric nanobeams. Eur Phys J Plus. 2018;133:1–18.10.1140/epjp/i2018-11954-7Search in Google Scholar

[33] Ke LL, Wang YS. Free vibration of size-dependent magneto-electro-elastic nanobeams based on the nonlocal theory. Physica E: Low Dimens Syst Nanostruct. 2014;63:52–61.10.1016/j.physe.2014.05.002Search in Google Scholar

[34] Ebrahimi F, Karimiasl M, Singhal A. Magneto-electro-elastic analysis of piezoelectric–flexoelectric nanobeams rested on silica aerogel foundation. Eng Comput. 2021;37:1007–14.10.1007/s00366-019-00869-zSearch in Google Scholar

[35] Ebrahimi F, Karimiasl M, Selvamani R. Bending analysis of magneto-electro piezoelectric nanobeams system under hygro-thermal loading. Adv Nano Res. 2020;8:203–14.Search in Google Scholar

[36] Ramezani SR, Mojra A. Stability analysis of conveying-nanofluid CNT under magnetic field based on nonlocal couple stress theory and fluid-structure interaction. Mech Based Des Struct Mach. 2020;51:583–600.10.1080/15397734.2020.1851254Search in Google Scholar

[37] Ebrahimi F, Barati MR. Vibration analysis of smart piezoelectrically actuated nanobeams subjected to magneto-electrical field in thermal environment. J Vib Control. 2018;24:549–64.10.1177/1077546316646239Search in Google Scholar

[38] Ebrahimi F, Barati MR, Dabbagh A. A nonlocal strain gradient theory for wave propagation analysis in temperature-dependent inhomogeneous nanoplates. Int J Eng Sci. 2016;107:169–82.10.1016/j.ijengsci.2016.07.008Search in Google Scholar

[39] Ebrahimi F, Salari E. Effect of various thermal loadings on buckling and vibrational characteristics of nonlocal temperature-dependent functionally graded nanobeams. Mech Adv Mater Struct. 2016;23:1379–97.10.1080/15376494.2015.1091524Search in Google Scholar

[40] Ebrahimi F, Barati MR. Hygrothermal effects on vibration characteristics of viscoelastic FG nanobeams based on nonlocal strain gradient theory. Compos Struct. 2017;159:433–44.10.1016/j.compstruct.2016.09.092Search in Google Scholar

[41] Ebrahimi F, Salari E. Thermal buckling and free vibration analysis of size dependent Timoshenko FG nanobeams in thermal environments. Compos Struct. 2015;128:363–80.10.1016/j.compstruct.2015.03.023Search in Google Scholar

[42] Ebrahimi F, Daman M. Nonlocal thermo-electro mechanical vibration analysis of smart curved FG piezoelectric Timoshenko nanobeam. Smart Struct Syst. 2017;20(3):351–68.Search in Google Scholar

[43] Vo TP, Thai HT, Nguyen TK, Inam F, Lee J. A quasi-3D theory for vibration and buckling of functionally graded sandwich beams. Compos Struct. 2015;119:1–12.10.1016/j.compstruct.2014.08.006Search in Google Scholar

[44] Jalaei MH, Arani AG, Nguyen-Xuan H. Investigation of thermal and magnetic field effects on the dynamic instability of FG Timoshenko nanobeam employing nonlocal strain gradient theory. Int J Mech Sci. 2019;161–162:105043.10.1016/j.ijmecsci.2019.105043Search in Google Scholar

[45] Karami B, Shahsavari D, Karami M, Li L. Hygrothermal wave characteristic of nanobeam-type inhomogeneous materials with porosity under magnetic field. Proc Inst Mech Eng C: J Mech Eng Sci. 2019;233(6):2149–69.10.1177/0954406218781680Search in Google Scholar

Received: 2022-11-17
Accepted: 2023-03-12
Published Online: 2023-04-03

© 2023 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Research Articles
  2. Investigation of differential shrinkage stresses in a revolution shell structure due to the evolving parameters of concrete
  3. Multiphysics analysis for fluid–structure interaction of blood biological flow inside three-dimensional artery
  4. MD-based study on the deformation process of engineered Ni–Al core–shell nanowires: Toward an understanding underlying deformation mechanisms
  5. Experimental measurement and numerical predictions of thickness variation and transverse stresses in a concrete ring
  6. Studying the effect of embedded length strength of concrete and diameter of anchor on shear performance between old and new concrete
  7. Evaluation of static responses for layered composite arches
  8. Nonlocal state-space strain gradient wave propagation of magneto thermo piezoelectric functionally graded nanobeam
  9. Numerical study of the FRP-concrete bond behavior under thermal variations
  10. Parametric study of retrofitted reinforced concrete columns with steel cages and predicting load distribution and compressive stress in columns using machine learning algorithms
  11. Application of soft computing in estimating primary crack spacing of reinforced concrete structures
  12. Identification of crack location in metallic biomaterial cantilever beam subjected to moving load base on central difference approximation
  13. Numerical investigations of two vibrating cylinders in uniform flow using overset mesh
  14. Performance analysis on the structure of the bracket mounting for hybrid converter kit: Finite-element approach
  15. A new finite-element procedure for vibration analysis of FGP sandwich plates resting on Kerr foundation
  16. Strength analysis of marine biaxial warp-knitted glass fabrics as composite laminations for ship material
  17. Analysis of a thick cylindrical FGM pressure vessel with variable parameters using thermoelasticity
  18. Structural function analysis of shear walls in sustainable assembled buildings under finite element model
  19. In-plane nonlinear postbuckling and buckling analysis of Lee’s frame using absolute nodal coordinate formulation
  20. Optimization of structural parameters and numerical simulation of stress field of composite crucible based on the indirect coupling method
  21. Numerical study on crushing damage and energy absorption of multi-cell glass fibre-reinforced composite panel: Application to the crash absorber design of tsunami lifeboat
  22. Stripped and layered fabrication of minimal surface tectonics using parametric algorithms
  23. A methodological approach for detecting multiple faults in wind turbine blades based on vibration signals and machine learning
  24. Influence of the selection of different construction materials on the stress–strain state of the track
  25. A coupled hygro-elastic 3D model for steady-state analysis of functionally graded plates and shells
  26. Comparative study of shell element formulations as NLFE parameters to forecast structural crashworthiness
  27. A size-dependent 3D solution of functionally graded shallow nanoshells
  28. Special Issue: The 2nd Thematic Symposium - Integrity of Mechanical Structure and Material - Part I
  29. Correlation between lamina directions and the mechanical characteristics of laminated bamboo composite for ship structure
  30. Reliability-based assessment of ship hull girder ultimate strength
  31. Finite element method on topology optimization applied to laminate composite of fuselage structure
  32. Dynamic response of high-speed craft bottom panels subjected to slamming loadings
  33. Effect of pitting corrosion position to the strength of ship bottom plate in grounding incident
  34. Antiballistic material, testing, and procedures of curved-layered objects: A systematic review and current milestone
  35. Thin-walled cylindrical shells in engineering designs and critical infrastructures: A systematic review based on the loading response
  36. Laminar Rayleigh–Benard convection in a closed square field with meshless radial basis function method
  37. Determination of cryogenic temperature loads for finite-element model of LNG bunkering ship under LNG release accident
  38. Roundness and slenderness effects on the dynamic characteristics of spar-type floating offshore wind turbine
Downloaded on 20.12.2025 from https://www.degruyterbrill.com/document/doi/10.1515/cls-2022-0192/html
Scroll to top button