Home Technology Nonlinear poro thermal vibration and parametric excitation in a magneto-elastic embedded nanobeam using homotopy perturbation technique
Article Open Access

Nonlinear poro thermal vibration and parametric excitation in a magneto-elastic embedded nanobeam using homotopy perturbation technique

  • Lakshmanan Anitha , Loganathan Vadivukarasi , Rajendran Selvamani EMAIL logo , Rossana Dimitri and Francesco Tornabene
Published/Copyright: October 19, 2024

Abstract

The primary focus of this study is to analyze the nonlinear vibration patterns and parametric excitation of embedded Euler–Bernoulli nanobeams subjected to thermo-magneto-mechanical loads. The Euler–Bernoulli nanobeam is developed with external parametric excitation. By utilizing nonlocal continuum theory and nonlinear von Karman beam theory, the governing equation of motion is derived. Subsequently, the homotopy perturbation technique is employed to determine the vibration frequencies. Finally, the modulation equation of Euler–Bernoulli nanobeams is derived for simply supported boundary conditions. The impacts of magnetic potential, temperature, damping coefficient, Winkler coefficient, and nonlocal parameters are tested numerically on nonlinear frequency–amplitude and parametric excitation–amplitude responses. Results demonstrate that physical variables significantly influence both nonlinear frequency behavior and parametric excitation.

1 Introduction

The captivating electromechanical properties of carbon nanostructures, such as carbon nanotubes and nanobeams, have drawn significant interest from researchers and scholars in the fields of advanced materials and engineering design. These nanostructures have found applications in various electromechanical devices, including light translucency [1,2], vibratory systems [3,4,5], gas atom diagnosis [6], storage units [7], and composite materials [8]. However, despite the recognized importance of small-scale effects on the characteristics and properties of nanostructures, classical plate theory is insufficient for assessing size effects in these structures [9]. To address this limitation, nonlocal elasticity theory, introduced by Eringen [10], has been widely adopted to examine size effects in nanostructures. The implementation of nonlocal elasticity theory has led to numerous theoretical investigations and significant advancements [11,12,13,14,15,16,17]. Much of the current research on micro/nanobeams focuses on their nonlinear properties. Nonlinear or high-amplitude vibration of beams – whether nano or micro – subjected to significant displacements occupies a crucial position in the engineering literature. For instance, Şimşek [18,19] explored the nonlinear vibration of nanobeams using nonlocal elasticity and strain gradient theories, highlighting the impact of small-scale effects on the nonlinear frequency response. Nazemnezhad and Hosseini-Hashemi [20] examined the nonlinear vibration behavior of functionally graded (FG) nanobeams under different boundary conditions, emphasizing the influence of the gradient index on nonlinear vibration characteristics. Additionally, Nourbakhsh et al. [21] employed the von Karman theory to analyze the effect of nonlinearity on the nonlinear frequency response of microbeams. Further contributions include Oskouie et al. [22], who presented the nonlinear frequency response of viscoelastic Euler–Bernoulli nanobeams, emphasizing the impact of viscoelastic properties on nonlinear vibration characteristics. Ghadiri et al. [23] utilized the multiple time scales method to study the nonlinear forced vibration of nanobeams experiencing a moving concentrated load supported by a viscoelastic foundation. He [24] investigated a coupling method that utilizes homotopy perturbation techniques (HPTs) for analyzing nonlinear problems. Barati [25] delved into nonlocal-strain gradient-forced vibration analysis of metal foam nanoplates, examining both uniform and graded porosities. Additional studies have explored various aspects of nanostructures. For example, Kováčik et al. [26] discussed the Poisson’s ratio of closed-cell aluminum foams, while Pourjabari et al. [27] explored the influence of porosity on free and forced vibration characteristics of graphene platelet reinforcement composite nanostructures. Chen et al. [28] contributed by discussing the free and forced vibrations of shear deformable FG porous beams, while Mirjavadi et al. [29] analyzed the nonlinear free and forced vibrations of graphene nanoplatelet-reinforced microbeams, considering geometrical imperfections. He [30] provided a new interpretation of the HPT. Meanwhile, Eltaher et al. [31] examined the coupling effects of nonlocal and surface energy on the vibration analysis of nanobeams. Reddy [32] discussed nonlocal theories related to the bending, buckling, and vibration of beams, and Aydogdu [33] investigated a comprehensive nonlocal beam theory to analyze these behaviors in nanobeams. Researchers have also recognized the importance of parametric excitation in electromechanical systems, particularly in energy harvesting systems [34,35,36], where a Duffing oscillator simulates the performance of energy harvesters. Parametric excitation has also been used to study nonlinear vibration and stability in various structures [37]. Darabi and Ganesan [37] and Wang [38] analyzed the effect of van der Waals interaction on the instability of double-walled nanobeams under parametric excitation. Similarly, Krylov et al. [39] explored pull-in instability in microdevices under parametric excitation using Mathieu and Hill’s equations. Yan et al. [40] provided insights into the behavior of Timoshenko beams subjected to parametric and external excitations. Moreover, Eringen [41,42] examined the theory of nonlocal polar elastic continua and related differential equations. Further studies by Reddy [43,44] have contributed to continuum mechanics and nonlocal nonlinear formulations for beam and plate bending, while Emam [45,46] explored static and dynamic analysis of post-buckling in geometrically imperfect composite beams. Murmu et al. [47] investigated the influence of in-plane magnetic fields on transverse vibration of graphene sheets, and Kitipornchai et al. [48] modeled vibration characteristics of multilayered graphene sheets. Additionally, researchers like Nayfeh and Mook [49], Nazemnezhad and Hosseini-Hashemi [20], and Azrar et al. [50,51] have contributed to the understanding of nonlinear oscillations and dynamic responses in beams. Characteristics and behaviors of various nanobeam configurations under different environmental and boundary conditions were investigated in previous studies [5257].

Hence, a literature review suggests that the impacts of nonlinear vibration and parametric excitation of magneto-thermo elastic embedded nanobeams have not been extensively investigated in existing studies. This research conducts a thorough examination of the nonlinear vibration behaviors and parametric excitation effects on embedded Euler–Bernoulli nanobeams under thermo-magneto-mechanical loads, along with external parametric excitation. To begin, a succinct model of the nanobeam is developed, followed by the application of an external axial force to induce parametric excitation. Subsequently, employing the nonlocal continuum theory and nonlinear von Karman beam theory, the governing nonlinear differential equation of motion is derived. The HPT is then employed to solve this equation. Finally, the modulation equation and the dynamic instability of the Euler–Bernoulli nanobeam are derived, leading to an examination of both trivial and nontrivial steady-state solutions. Results demonstrate that the magnetic potential, temperature, damping coefficient, Winkler coefficient, and nonlocal parameters have a significant impact on both nonlinear frequency behavior and parametric excitation.

2 Modeling of porous metal nanobeam

The metal’s material traits are contingent upon the distribution of voids or pores. These voids can be distributed uniformly or in non-uniform patterns. In cases of non-uniform distribution, it can be further categorized as symmetric (non-uniform 1) or asymmetric (non-uniform 2). Subsequently, the forthcoming section will introduce the expressions for the material properties, specifically the elastic modulus (E) and mass density (ρ), pertaining to metal foam [25].

(2.1) E = E 2 ( 1 e 0 ϒ ) , ρ = ρ 2 ( 1 e 0 ϒ )   ( Uniform ) ,

where ϒ = 1 e 0 1 e 0 2 π 1 e 0 2 π + 1 2 .

The above equation defines ϒ as a function of the coefficient e 0 , which relates to the pore amount in the metal foam. This variable ϒ is used in the expression for E and ρ for uniformly distributed voids in metal foam.

(2.2) E ( z ) = E 2 1 e 0 cos π z h , ρ ( z ) = ρ 2 1 e m cos π z h ( Non uniform  1 ) ,

E ( z ) = E 2 1 e 0 cos π z 2 h + π 4 ,

(2.3) ρ ( z ) = ρ 2 1 e m cos π z h + π 4 ( Non uniform  2 ) ,

where z is the spatial coordinate and h is the characteristic thickness. In the above definitions, the index 2 refers to a material property at its highest value. Also, there are two coefficients e 0 and e m elated to pore amount and mass distribution as

(2.4)   e 0 = 1 E 2 E 1 = 1 G 2 G 1 , e m = 1 1 e 0 ,

where E 1 and G 1 are the elastic modulus and shear modulus of the metal without voids. E 2 and G 2 are the highest values of the elastic modulus and shear modulus in the presence of voids.

3 Formulating the problem

Figure 1 illustrates the schematic of a porous nanobeam embedded in a visco-Pasternak foundation, subjected to an axial force along the x-axis with a height h and length L. The axial force is represented as a function undergoing harmonic excitation with frequency (Ω̄). Moreover, the vertical displacement of the nanobeam is indicated by w along the z-axis.

Figure 1 
               (a) Geometry of the beam. (b) Porosity distribution. (i) Uniform porosity. (ii) Non-uniform porosity.
Figure 1

(a) Geometry of the beam. (b) Porosity distribution. (i) Uniform porosity. (ii) Non-uniform porosity.

3.1 Governing equations

In accordance with Eringen’s nonlocal elasticity theory [9,10,41,42], the stress experienced at a reference point X is postulated to be contingent upon the strain field throughout the body at each point X′. The nonlocal stress tensor σ at point X is formulated as follows:

(3.1) σ = V K ( X X , τ ) σ ( X ' ) d X .

Here σ denotes the classical stress tensor and K ( X X ) represents the Kernel function, which signifies the nonlocal modulus. Eringen [10,42] illustrates that it is feasible to express the integral constitutive relation in an equivalent differential form as:

(3.2) ( 1 ( e 0 a ) 2 2 ) σ = σ ' ,

where 2 = 2 x 2 + 2 y 2 represents the Laplacian operator, and ( e 0 a ) introduces the nonlocal parameter, where e 0 is a material-specific constant and a is the internal characteristic length. The determination of the value of e 0 typically involves experimental methods or matching the dispersion relation of plane waves with those of atomic lattice dynamics. Subsequently, the nonlocal constitutive relation for the Euler–Bernoulli nanobeam can be expressed as:

(3.3) σ x x ( e 0 a ) 2 2 σ x x x 2 = E ε x x ,

where σ x x and ε x x denote the normal stress and strain, respectively, while E represents Young’s modulus. Following the Euler–Bernoulli beam model, the axial force and the resultant bending moment can be formulated as

(3.4) { N , M } = A σ x ( 1 , z ) d A ,

where z represents the transverse coordinate in the deflection direction, and A denotes the area of the cross-section of the nanobeam. Utilizing the classical beam theory as outlined by Reddy [43,44], the displacements can be expressed as follows:

(3.5) u 1 ( x , z , t ) = u ( x , t ) z w x , u 2 = 0 ,   u 3 ( x , z , t ) = w ( x , t ) .

In Eq. (3.5), u and w represent the axial and transverse displacements of the nanobeam along x and z directions, respectively. Now, considering the nonlinear von Karman strain, we can express it as:

(3.6) ε = ε 0 + z k ,

where ε is the strain vector, and ε 0 and k are the nonlinear strain vector and the change in the curvature vector, respectively, defined as follows:

(3.7) ε 0 = u 0 x + 1 2 w x 2 , k = 2 w x 2 .

In this context, u 0 represents the initial axial displacement in the strain expression, which captures the axial deformation before accounting for the additional effects due to transverse displacements and curvature. From Eqs. (3.3)–(3.7) the axial load and the bending moment are as follows:

(3.8) ( 1 ( e 0 a ) 2 2 ) N = E A ε 0 ( 1 ( e 0 a ) 2 2 ) M = E I k ,

where I = A z 2 d A represents the moment of inertia. Therefore, the equation of motion can be expressed as [45,46]:

(3.9) E I 4 w x 4 x N w x + ( e 0 a ) 2 3 x 3 N w x + ρ A 2 t 2 w ( e 0 a ) 2 2 w x 2 = f ( e 0 a ) 2 2 f x 2 .

The axial normal force N can be determined as follows:

(3.10) N = M x + T x + N x + F cosΩ t E A 2 L 0 L w x 2 d x .

In Eq. (3.10), M x , T x , and N x represent a uniaxial magnetic field, thermal load caused by temperature change, and in-plane load caused by initial stress, respectively. Additionally, the term of F cosΩ t is also the axial force capable of inducing parametric excitation, we can define the parameters of the axial normal force as follows:

(3.11) M x = η H x 2 2 w x 2 ,    T x = α E A T   ,    N x = ξ σ 0 .

Here H x and η represent the in-plane uniaxial magnetic field and the magnetic field permeability, respectively. Specifically, M x explains the Lorentz force along the x-axis [47]. T x , α , A, and T denote the coefficient of thermal expansion, the cross-sectional area, the difference between the temperature, and its initial reference temperature, respectively. Moreover, ξ and σ 0 are the compression ratio and the initial stress, respectively. In this study, it is assumed that ξ = 1 and initial stress is along the x -axis direction. Additionally, in Eq. (3.9), f is defined as follows:

(3.12) f = K w w + c d w t ,

where k w and c d represent the linear coefficient of Winkler and damper modulus parameter, respectively. The Winkler-type foundation can be characterized based on the model from the study of Kitipornchai et al. [48]. Finally, to obtain the equation of motion, we substitute Eqs. (3.10)–(3.12) into Eq. (3.9) as follows:

(3.13) E I 4 w x 4   M x + T x + N x + F cosΩ t EA 2 L 0 L w x 2 d x   2 w x 2 + c d w t + K w w T x + N x + F cosΩ t EA 2 L 0 L w x 2 d x ( e 0 a ) 2 K w 2 w x 2 ( e 0 a ) 2 c d 3 w t x 2 = ρ A 2 w t 2 ( e 0 a ) 2 4 w t 2 x 2 .

To facilitate a good comparison between results, indirect parameters can be articulated as follows:

(3.14) X = x L ,   w 0 = W L ,   γ = e 0 a L ,   K w = K w L 4 EI , T x = T x L 2 EI , N x = N x L 2 EI , M x = M x L 2 EI ,   c d = c d L 4 EI , Ω = Ω ρ A L 4 EI .

By incorporating these indirect parameters and substituting them into Eq. (3.9), the governing equation of nanobeam can be derived as follows:

(3.15) 4 w 0 X 4 M x + T x + N x + F cosΩ t E A 2 L 0 L w 0 x 2 d x   2 w 0 X 2 + c d w 0 t + K w w 0 +   ( γ ) 2 4 w 0 X 4 M x + T x + N x + F c o s Ω t E A 2 L 0 L w 0 x 2 d x   ( γ ) 2 K w 2 w 0 X 2   ( γ ) 2 c d 3 w 0 t X 2 = 2 w 0 t 2 ( γ ) 2 4 w 0 t 2 X 2

4 HPT

HPT offers an analytical approximate solution for problems that exhibit continuity within the solution domain. This technique involves considering a differential equation.

(4.1) L y + N y = f ( x ) ,   x ϵ Ω .

Under the boundary condition, B y , y x = 0 , x ∈ Γ . Here L is the linear operator, N is the nonlinear operator, B is the boundary operator, Γ represents the boundary of the domain Ω, and (x) is the known analytic function. HPT defines a homotopy as υ ( x , p ) = Ω × [ 0 , 1 ] R that satisfies the following inequalities:

(4.2) H ( υ , p ) = ( 1 p ) [ L ( v ) L ( y 0 ) ] + p [ L ( v ) + N ( v ) f ( x ) ] = 0 ,

or

(4.3) H ( υ , p ) = L ( v ) L ( y 0 ) + p L ( y 0 ) + p [ N ( v ) f ( x ) ] = 0 ,

where x ϵ Ω and p ϵ [ 0 , 1 ] , and y 0 is an initial approximation that satisfies the boundary condition. Now from Eqs. (4.2) to (4.3) one can obtain

H ( υ , 0 ) = L ( v ) L ( y 0 ) = 0 ,

(4.4) H ( υ , 1 ) = L ( v ) + N ( v ) f ( x ) = 0 .

In topology, L ( v ) L ( y 0 ) and L ( v ) + N ( v ) f ( x ) are homotopic functions, indicating that there exists a continuous deformation between them through the homotopy parameter p .

Consider the power series solution of (4.2)–(4.3) as follows:

(4.5) υ = υ 0 + p υ 1 + p 2 υ 2 + p 3 υ 3 + .

Hence, the approximate solution of (4.2) can be obtained

(4.6) y = lim p 1 v = v + υ 0 + υ 1 + υ 2 + .

4.1 Implementation of boundary conditions in HPT

Simply Supported–Simply Supported (S–S)

w 0 = d 2 w 0 d x 2 = d 4 w 0 d x 4 = 0  at x = 0 ,

(4.7) w 0 = d 2 w 0 d x 2 = d 4 w 0 d x 4 = 0  at x = 1 .

So after using the above boundary conditions in the nth order approximate solution the system of homogeneous equation can be written as

(4.8) M ( λ ) [ A   B   C   D ] T = 0 .

For a nontrivial solution, determinant of coefficient matrix must be zero. The determinant of coefficient matrix yields a characteristic equation in terms of λ . Positive real roots of this equation are the normalized free vibration frequencies for the case considered.

4.2 HPT formulation for present problem

Consider the nondimensional differential Eq. (3.15), this equation can be reformulated as

(4.9) 4 w 0 X 4 = N γ 2 N 1 2 w 0 X 2 + 1 γ 2 N 1 2 w 0 t 2 γ 2 γ 2 N 1 4 w 0 t 2 X 2 c d w 0 t K w w 0 .

The homotopy can be applied as [35]

(4.10) 4 w 0 X 4 = p N γ 2 N 1 2 w 0 X 2 + 1 γ 2 N 1 2 w 0 t 2 γ 2 γ 2 N 1 4 w 0 t 2 X 2 c d w 0 t K w w 0 ,

where p is the homotopy parameter, p ∈ [0, 1]. It is obvious that when p = 0, the equation becomes homogeneous, i.e.,

(4.11) 4 w X 4 = 0 .

The initial approximation W 0 is obtained by solving the homogeneous Eq. (4.11), hence

(4.12) w 0 = A x 4 + B x 3 + C x 2 + D .

The basic assumption of the HPT is that the solution of Eq. (4.10) can be written as a power series in p

(4.13) w aprprox = w 0 + p w 1 + p 2 w 2 + .

By substituting the value of w aprprox from Eq. (4.13) in Eq. (4.10)

(4.14) 4 X 4 [ w 0 + p w 1 + p 2 w 2 + ] = p N γ 2 N 1 2 X 2 [ w 0 + p w 1 + p 2 w 2 + ] + 1 γ 2 N 1 2 t 2 [ w 0 + p w 1 +   p 2 w 2 + ] γ 2 γ 2 N 1 4 t 2 X 2 w 0 + p w 1 + p 2 w 2 + c d t w 0 + p w 1 + p 2 w 2 + ]   k w [ w 0 + p w 1 + p 2 w 2 + .

Now comparing the coefficients of p in Eq. (4.14), the recurrence relation can be obtained as

(4.15) 4 w i X 4 = N γ 2 N 1 2 w i 1 X 2 + 1 γ 2 N 1 2 w i 1 t 2 γ 2 γ 2 N 1 4 w i 1 t 2 X 2 c d w i 1 t K w w i 1 ,

where i ≥ 1 and initial guess w 0 is given in Eq. (4.12). The solution to Eq. (4.9) can be approximated as per the study of He [24].

(4.16) w 0 = lim p 1 w aprrox = w 0 + w 1 + w 2 + .

The convergence of the series in Eq. (4.16) is proved in previous studies [24,30].

5 Numerical results

In this section, numerical results are examined based on the application of both thermo-magneto-mechanical loading and external parametric excitation. The focus lies in understanding the impact of parametric excitation through the examination of instability regions and bifurcation points. To aid comprehension, key parameters are defined across various regions of the graph. This helps to elucidate the concepts and enhance clarity. The system’s material properties, including those of the nanobeam and elastic matrix, consist of the following parameters: Temperature T = 300 K, Poisson’s ratio of the beam material ν = 0.3 , Winkler coefficient k w = 0.5 GPa, and a viscoelastic damping coefficient c d = 3 × 10 7  Pa s . Additionally, the nanobeam diameter d = 3 nm and the small-scale parameter are considered to be smaller than 2 nm according to the study of Reddy [32]. One can explore its dependence on several factors, including the scaling parameter, porosity pattern, and nonlocal effects. To do this, a set of material constants are taken in Table 1.

Table 1

Beam dimensions and their material properties [26]

Parameter Value (unit) Description
a 0.5 m Length of the beam
b 0.05 m Width of the beam
h 0.05 m Height of the beam
E 1,100 GPa Young’s modulus
ρ 1.3 g/cm3 Mass density

First, to verify the accuracy of the formulation, Table 2 is presented. The numerical results of the present study reported in the table are compared with other available research studies and literature [20,50] so that, they are partly similar and close to our research. Table 2 shows the nonlinear frequency ratio ( ω NL , ω L ) for amplitude-to-radius ( ω max / r ) ratio of isotropic beam with simply supported boundary conditions. The nonlinear frequency ratio is tabled for different amplitude ratios [1,2,3]. The results presented in the study of Nazemnezhad and Hosseini-Hashemi [20] are consistent with the numerical results of the present work as a similar analytical approach was used in both.

Table 2

Frequency ratio ( ω NL , ω L ) at different maximum amplitude-to-radius ( ω max / r ) ratios of isotropic beam with simply supported boundary conditions

Amplitude ratio Ref. [20] Ref. [50] Present work
1 1.0937 1.0892 1.0758
2 1.3750 1.3178 1.3690
3 1.8438 1.6257 1.8401

Figure 2 illustrates that the relationship between the nonlinear frequency and the amplitude of the parametric excitation varies depending on the nondimensional damping coefficient c d . When c d is low, indicating less damping in the system, higher oscillation amplitudes are achievable. Consequently, the nonlinear frequency tends to be higher in such cases. This implies that the nonlinear frequency increases together with a rise in the parametric excitation amplitude, producing a more intense nonlinear response. In summary, lower values of c d cause the nonlinear frequency to decrease as the amplitude increases, whereas lower values of the nondimensional damping coefficient cause the nonlinear frequency to increase with increasing amplitude of parametric stimulation. Figure 3 displays the relationship between the force amplitude and the nonlinear frequency amplitude concerning the nondimensional Winkler coefficient k w . The graph illustrates that as the Winkler coefficient increases and the force amplitude decreases. However, it is noteworthy that the Winkler coefficient does not influence the occurrence of bifurcation points. In other words, while changes in k w affect the force amplitude; they do not impact the system’s bifurcation behavior. Figure 4 shows how, for a range of uniaxial magnetic field values, the nonlinear frequency varies with the magnitude of parametric excitation. The force nonlinear frequency of a system subjected to a uniaxial magnetic field denotes the frequency at which the system’s response turns nonlinear as a result of the magnetic field’s effect. The system’s response to the magnetic field may exhibit nonlinear behavior as the amplitude of parametric stimulation increases, resulting in changes in the nonlinear frequency. Figure 5 illustrates the nonlocal parameter’s effect on the nonlinear frequency–response curves. It is clear from the graph that increasing the nonlocal parameter reduces the hardening behavior while enhancing the bending stiffness. The relationship between the parametric excitation amplitude and the force amplitude of the nonlinear frequency is shown for different temperatures in Figure 6. Additionally, it is observed that an increase in the temperature gradient results in a higher amplitude at the lower limit point bifurcation. Conversely, elevating the initial imperfection amplitude leads to greater response amplitude at the higher limit point bifurcation. Figure 7 illustrates the impact of varying the parameter excitation on the amplitude concerning the Winkler coefficient k w . Notably, as the parameter excitation declines in the amplitude of the Winkler coefficient and damps at some point (0.1–0.15). Nevertheless, amplitude 1 shows a growing trend while the parametric excitation keeps increasing. Moreover, it is worth noting that the smallest amplitudes occur at a parametric excitation level of k w = 1.5. In Figure 8, the amplitude for a system under parametric excitation depends on both the force amplitude k and the nondimensional damping coefficient c d . Lower damping generally leads to larger amplitudes near resonance, while higher damping suppresses oscillations and keeps the amplitude relatively low. Critical damping minimizes the amplitude of the system’s response. In Figure 9, the effect of changing the parametric excitation on the amplitude is demonstrated concerning the uniaxial magnetic field H x . Initially, the parameter excitation declines in the amplitude and damps at some point (0.1–0.15). However, as the parametric excitation continues to increase, the amplitude exhibits a rising trend. Additionally, it is noteworthy that the smallest amplitudes are observed at a parametric excitation level of H x = 0.5. The relationship between the parametric excitation and amplitude for various temperature levels is depicted in Figure 10. It is found that a gradual decrease in temperature starting around 1.0 and damping at some point (0.1–0.15). On the other hand, the amplitude shows a rising trend while the parametric excitation keeps increasing. Additionally, it is noteworthy that the smallest amplitudes are observed at a parametric excitation level of T = 300. In Figure 11, the relationship between parametric excitation and amplitude is illustrated across various values of the nondimensional nonlocal parameter γ . As the γ value gradually decreases, starting from approximately 1.0 and tapering off between 0.1 and 0.15, a resonant behavior is evident. This resonance is attributed to the combined effects of parametric excitation and amplitude on the system. Figures 12 and 13 show the stress contours of the magneto-thermo-elastic nanobeam for different porosity distributions. We can find that the stress change zone is restricted in a finite area and the stress does not change out of this area. The blue color in these figures refers to a stress variation of zero in this region. We can observe that the region with changes in stress becomes larger with lower nonlocal values.

Figure 2 
               The effect of nonlinear frequency on the amplitude of parametric excitation for different values of damping coefficient 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       c
                                    
                                    
                                       d
                                    
                                 
                              
                              )
                           
                        
                        ({c}_{\text{d}})
                     
                  .
Figure 2

The effect of nonlinear frequency on the amplitude of parametric excitation for different values of damping coefficient ( c d ) .

Figure 3 
               The effect of nonlinear frequency on the amplitude of parametric excitation for different values of Winkler coefficient 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       k
                                    
                                    
                                       w
                                    
                                 
                              
                              )
                           
                        
                        ({k}_{\text{w}})
                     
                  .
Figure 3

The effect of nonlinear frequency on the amplitude of parametric excitation for different values of Winkler coefficient ( k w ) .

Figure 4 
               The effect of nonlinear frequency on the amplitude of parametric excitation for different values of uniaxial magnetic field (
                     
                        
                        
                           
                              
                                 H
                              
                              
                                 x
                              
                           
                           )
                        
                        {H}_{x})
                     
                  .
Figure 4

The effect of nonlinear frequency on the amplitude of parametric excitation for different values of uniaxial magnetic field ( H x ) .

Figure 5 
               The effect of nonlinear frequency on the amplitude of parametric excitation for different values of nonlocal parameter (
                     
                        
                        
                           γ
                           )
                        
                        \gamma )
                     
                  .
Figure 5

The effect of nonlinear frequency on the amplitude of parametric excitation for different values of nonlocal parameter ( γ ) .

Figure 6 
               The effect of nonlinear frequency on the amplitude of parametric excitation for different values of temperature (T).
Figure 6

The effect of nonlinear frequency on the amplitude of parametric excitation for different values of temperature (T).

Figure 7 
               The effect of parametric excitation on the amplitude of parametric excitation for different values of Winkler coefficient 
                     
                        
                        
                           (
                           
                              
                                 k
                              
                              
                                 w
                              
                           
                           )
                        
                        ({k}_{\text{w}})
                     
                  .
Figure 7

The effect of parametric excitation on the amplitude of parametric excitation for different values of Winkler coefficient ( k w ) .

Figure 8 
               The effect of parametric excitation on the amplitude of parametric excitation for different values of damping coefficient 
                     
                        
                        
                           (
                           
                              
                                 c
                              
                              
                                 d
                              
                           
                           )
                        
                        ({c}_{\text{d}})
                     
                  .
Figure 8

The effect of parametric excitation on the amplitude of parametric excitation for different values of damping coefficient ( c d ) .

Figure 9 
               The effect of parametric excitation on the amplitude of parametric excitation for different values of uniaxial magnetic field (
                     
                        
                        
                           
                              
                                 H
                              
                              
                                 x
                              
                           
                           )
                        
                        {H}_{x})
                     
                  .
Figure 9

The effect of parametric excitation on the amplitude of parametric excitation for different values of uniaxial magnetic field ( H x ) .

Figure 10 
               The effect of parametric excitation on the amplitude of parametric excitation for different values of temperature (T).
Figure 10

The effect of parametric excitation on the amplitude of parametric excitation for different values of temperature (T).

Figure 11 
               The effect of parametric excitation on the amplitude of parametric excitation for different values of nondimensional nonlocal parameter 
                     
                        
                        
                           
                              (
                              
                                 γ
                              
                              )
                           
                        
                        (\gamma )
                     
                  .
Figure 11

The effect of parametric excitation on the amplitude of parametric excitation for different values of nondimensional nonlocal parameter ( γ ) .

Figure 12 
               The contour plot of the stress distribution for uniform porosity with 
                     
                        
                        
                           
                              
                                 e
                              
                              
                                 0
                              
                           
                           =
                           0
                        
                        {e}_{0}=0
                     
                  .
Figure 12

The contour plot of the stress distribution for uniform porosity with e 0 = 0 .

Figure 13 
               The contour plot of the stress distribution for Nonuniform porosity with 
                     
                        
                        
                           
                              
                                 e
                              
                              
                                 0
                              
                           
                           =
                           0.5
                           .
                        
                        {e}_{0}=0.5.
Figure 13

The contour plot of the stress distribution for Nonuniform porosity with e 0 = 0.5 .

6 Conclusions

The aim of this research is to examine the dynamic parametric excitation and nonlinear vibration behavior of Euler–Bernoulli porous nanobeams under thermo-magneto-mechanical loading. Initially, a concise model of the Euler–Bernoulli nanobeam is developed and subjected to parametric external excitation. Utilizing the nonlocal continuum theory and nonlinear von Karman beam theory, the governing nonlinear differential equation of motion is derived. The partial differential equation is then converted into an ordinary differential equation using the HPT. Next, the Euler–Bernoulli nanobeam modulation equation is found. Special attention is given to the influence of parametric excitation, and bifurcation points are scrutinized to delineate instability regions. Notably, it is observed that the damping coefficient, along with parametric excitation, significantly affects the system stability and frequency responsiveness. Thermo-magneto-mechanical loads are found to induce either growth or decay in the amplitude. The following is a list of the study’s other main results:

  • The influence of parametric excitation induced by an external axial force on system stability is substantial.

  • The damping coefficient significantly influences system stability, while factors such as the nonlocal parameter and Winkler coefficient are of less importance.

  • Amplitude response is observed to vary as a function of the excitation frequency. For initial amplitudes of significant magnitude, the response decays until reaching a steady-state solution.

  • An increase in force amplitude leads to a notable separation between stable and unstable curves, creating a gap between them.

  • Nano-size beam having nonuniform pores 2 results in greater vibration frequency.

  • Results demonstrate that physical variables significantly influence both nonlinear frequency behavior and parametric excitation. The numerical results serve as reference points for conducting further analyses of nanobeams, which serve as fundamental components in nano-electromechanical systems.

  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 consented to its submission to the journal, reviewed all the results, and approved the final version of the manuscript. RS and LA conceived of the presented idea, developed the theory, and performed the computations. LV verified the analytical methods. RD and FT supervised the finding of this entire work. All authors discussed the results and contributed to the final manuscript.

  3. Conflict of interest: Authors F.T. and R.D., 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, and it was handled entirely by other editors of the journal. The authors declare no other conflict of interest.

References

[1] Eda G, Fanchini G, Chhowalla M. Large-area ultrathin films of reduced graphene oxide as a transparent and flexible electronic material. Nat Nanotechnol. 2008;3(5):270–4. 10.1038/nnano.2008.83.Search in Google Scholar PubMed

[2] Li D, Müller MB, Gilje S, Kaner RB, Wallace GG. Processable aqueous dispersions of graphene nanosheets. Nat Nanotechnol. 2008;3(2):101–5. 10.1038/nnano.2007.451.Search in Google Scholar PubMed

[3] Potekin R, Kim S, McFarland DM, Bergman LA, Cho H, Vakakis AF. A micromechanical mass sensing method based on amplitude tracking within an ultra-wide broadband resonance. Nonlinear Dyn. 2018;92:287–304. 10.1007/s11071-018-4055-y.Search in Google Scholar

[4] Mahmoud MA. Validity and accuracy of resonance shift prediction formulas for microcantilevers: a review and comparative study. Crit Rev Solid State Mater Sci. 2016;41(5):386–429. 10.1080/10408436.2016.1142858.Search in Google Scholar

[5] Ji Y, Choe M, Cho B, Song S, Yoon J, Ko HC, et al. Organic nonvolatile memory devices with charge trapping multilayer graphene film. Nanotechnology. 2012;23(10):105202. 10.1088/0957-4484/23/10/105202.Search in Google Scholar PubMed

[6] Arash B, Wang Q. Detection of gas atoms with carbon nanotubes. Sci Rep. 2013;3(1):1782. 10.1038/srep01782.Search in Google Scholar

[7] Bunch JS, Van Der Zande AM, Verbridge SS, Frank IW, Tanenbaum DM, Parpia JM, et al. Electromechanical resonators from graphene sheets. Science. 2007;315(5811):490–3. 10.1126/science.1136836.Search in Google Scholar PubMed

[8] Kuilla T, Bhadra S, Yao D, Kim NH, Bose S, Lee JH. Recent advances in graphene based polymer composites. Prog Polym Sci. 2010;35(11):1350–75. 10.1016/j.progpolymsci.2010.07.005.Search in Google Scholar

[9] Eringen AC, Edelen D. On nonlocal elasticity. Int J Eng Sci. 1972;10(3):233–48.10.1016/0020-7225(72)90039-0Search in Google Scholar

[10] Eringen AC. Theories of nonlocal plasticity. Int J Eng Sci. 1983;21(7):741–51. 10.1016/0020-7225(83)90058-7.Search in Google Scholar

[11] Ghadiri M, Shafiei N, Akbarshahi A. Influence of thermal and surface effects on vibration behavior of nonlocal rotating Timoshenko nanobeam. Appl Phys A. 2016;122:1–9. 10.1007/s00339-016-0196-3.Search in Google Scholar

[12] Sudak LJ. Column buckling of multiwalled carbon nanotubes using nonlocal continuum mechanics. J Appl Phys. 2003;94(11):7281–7. 10.1063/1.1625437.Search in Google Scholar

[13] Zhang YQ, Liu GR, Wang JS. Small-scale effects on buckling of multiwalled carbon nanotubes under axial compression. Phys Rev B. 2004;70(20):205430. 10.1103/PhysRevB.70.205430.Search in Google Scholar

[14] Barretta R, Feo L, Luciano R, de Sciarra FM. Variational formulations for functionally graded nonlocal Bernoulli–Euler nanobeams. Compos Struct. 2015;129:80–9. 10.1016/j.compstruct.2015.03.033.Search in Google Scholar

[15] Ghadiri M, Safi M. Nonlinear vibration analysis of functionally graded nanobeam using homotopy perturbation method. Adv Appl Math Mech. 2017;9(1):144–56. 10.4208/aamm.2015.m899.Search in Google Scholar

[16] Ehyaei J, Akbarshahi A, Shafiei N. Influence of porosity and axial preload on vibration behavior of rotating FG nanobeam. Adv Nano Res. 2017;5(2):141. 10.12989/anr.2017.5.2.141.Search in Google Scholar

[17] Ebrahimi F, Hosseini SH. Thermal effects on nonlinear vibration behavior of viscoelastic nanosize plates. J Therm Stresses. 2016;39(5):606–25. 10.1080/01495739.2016.1160684.Search in Google Scholar

[18] Şimşek M. Large amplitude free vibration of nanobeams with various boundary conditions based on the nonlocal elasticity theory. Compos Part B: Eng. 2014;56:621–8. 10.1016/j.compositesb.2013.08.082.Search in Google Scholar

[19] Şimşek M. Nonlinear free vibration of a functionally graded nanobeam using nonlocal strain gradient theory and a novel Hamiltonian approach. Int J Eng Sci. 2016;105:12–27. 10.1016/j.ijengsci.2016.04.013.Search in Google Scholar

[20] Nazemnezhad R, Hosseini-Hashemi S. Nonlocal nonlinear free vibration of functionally graded nanobeams. Compos Struct. 2014;110:192–9. 10.1016/j.compstruct.2013.12.006.Search in Google Scholar

[21] Nourbakhsh H, Mohammadzadeh R, Rafiee M, Rafiee R. Nonlinear effects on resonance behaviour of beams in micro scale. Appl Mech Mater. 2012;110:4178–86. 10.4028/www.scientific.net/AMM.110-116.4178.Search in Google Scholar

[22] Oskouie MF, Ansari R, Sadeghi F. Nonlinear vibration analysis of fractional viscoelastic Euler–Bernoulli nanobeams based on the surface stress theory. Acta Mech Solida Sin. 2017;30(4):416–24. 10.1016/j.camss.2017.07.003.Search in Google Scholar

[23] Ghadiri M, Rajabpour A, Akbarshahi A. Non-linear forced vibration analysis of nanobeams subjected to moving concentrated load resting on a viscoelastic foundation considering thermal and surface effects. Appl Math Model. 2017;50:676–94. 10.1016/j.apm.2017.06.019.Search in Google Scholar

[24] He JH. A coupling method of a homotopy technique and a perturbation technique for non-linear problems. Int J Non-Linear Mech. 2000;35(1):37–43. 10.1016/S0020-7462(98)00085-7.Search in Google Scholar

[25] Barati MR. Nonlocal-strain gradient forced vibration analysis of metal foam nanoplates with uniform and graded porosities. Adv Nano Res. 2017;5(4):393. 10.12989/anr.2017.5.4.393.Search in Google Scholar

[26] Kováčik J, Marsavina L, Linul E. Poisson’s ratio of closed-cell aluminium foams. Materials. 2018;11(10):1904. 10.3390/ma11101904.Search in Google Scholar PubMed PubMed Central

[27] Pourjabari A, Hajilak ZE, Mohammadi A, Habibi M, Safarpour H. Effect of porosity on free and forced vibration characteristics of the GPL reinforcement composite nanostructures. Comput Math Appl. 2019;77(10):2608–26. 10.1016/j.camwa.2018.12.041.Search in Google Scholar

[28] Chen D, Yang J, Kitipornchai S. Free and forced vibrations of shear deformable functionally graded porous beams. Int J Mech Sci. 2016;108:14–22. 10.1016/j.ijmecsci.2016.01.025.Search in Google Scholar

[29] Mirjavadi SS, Afshari BM, Barati MR, Hamouda AM. Nonlinear free and forced vibrations of graphene nanoplatelet reinforced microbeams with geometrical imperfection. Microsyst Technol. 2019;25:3137–50. 10.1007/s00542-018-4277-4.Search in Google Scholar

[30] He JH. Addendum: new interpretation of homotopy perturbation method. Int J Mod Phys B. 2006;20(18):2561–8. 10.1142/S0217979206034819.Search in Google Scholar

[31] Eltaher MA, Mahmoud FF, Assie AE, Meletis E. Coupling effects of nonlocal and surface energy on vibration analysis of nanobeams. Appl Math Comput. 2013;224:760–74. 10.1016/j.amc.2013.09.002.Search in Google Scholar

[32] Reddy J. Nonlocal theories for bending, buckling and vibration of beams. Int J Eng Sci. 2007;45(2–8):288–307. 10.1016/j.ijengsci.2007.04.004.Search in Google Scholar

[33] Aydogdu M. A general nonlocal beam theory: its application to nanobeam bending, buckling and vibration. Phys E: Low-Dimension Syst Nanostruct. 2009;41(9):1651–5. 10.1016/j.physe.2009.05.014.Search in Google Scholar

[34] Alevras P, Theodossiades S, Rahnejat H. Broadband energy harvesting from parametric vibrations of a class of nonlinear Mathieu systems. Appl Phys Lett. 2017;110(23):233901. 10.1063/1.4984059.Search in Google Scholar

[35] Amer YA, El-Sayed AT, Kotb AA. Nonlinear vibration and of the Duffing oscillator to parametric excitation with time delay feedback. Nonlinear Dyn. 2016;85:2497–505. 10.1007/s11071-016-2840-z.Search in Google Scholar

[36] Bobryk RV, Yurchenko D. On enhancement of vibration-based energy harvesting by a random parametric excitation. J Sound Vib. 2016;366:407–17. 10.1016/j.jsv.2015.11.033.Search in Google Scholar

[37] Darabi M, Ganesan R. Non-linear vibration and dynamic instability of internally-thickness-tapered composite plates under parametric excitation. Compos Struct. 2017;176:82–104. 10.1016/j.compstruct.2017.04.059.Search in Google Scholar

[38] Wang YZ. Nonlinear internal resonance of double-walled nanobeams under parametric excitation by nonlocal continuum theory. Appl Math Model. 2017;48:621–34. 10.1016/j.apm.2017.04.018.Search in Google Scholar

[39] Krylov S, Harari I, Cohen Y. Stabilization of electrostatically actuated microstructures using parametric excitation. J Micromech Microeng. 2005;15(6):1188. 10.1088/0960-1317/15/6/009.Search in Google Scholar

[40] Yan Q, Ding H, Chen L. Nonlinear dynamics of axially moving viscoelastic Timoshenko beam under parametric and external excitations. Appl Math Mech. 2015;36(8):971–84. 10.1007/s10483-015-1966-7.Search in Google Scholar

[41] Eringen AC. Nonlocal polar elastic continua. Int J Eng Sci. 1972;10(1):1–6. 10.1016/0020-7225(72)90070-5.Search in Google Scholar

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

[43] Reddy JN. An introduction to continuum mechanics. United Kingdom: Cambridge University Press; 2013.Search in Google Scholar

[44] Reddy J. Nonlocal nonlinear formulations for bending of classical and shear deformation theories of beams and plates. Int J Eng Sci. 2010;48(11):1507–18. 10.1016/j.ijengsci.2010.09.020.Search in Google Scholar

[45] Emam SA. A static and dynamic analysis of the postbuckling of geometrically imperfect composite beams. Compos Struct. 2009;90(2):247–53. 10.1016/j.compstruct.2009.03.020.Search in Google Scholar

[46] Emam SA, Nayfeh AH. Postbuckling and free vibrations of composite beams. Compos Struct. 2009;88(4):636–42. 10.1016/j.compstruct.2008.06.006.Search in Google Scholar

[47] Murmu T, McCarthy MA, Adhikari S. In-plane magnetic field affected transverse vibration of embedded single-layer graphene sheets using equivalent nonlocal elasticity approach. Compos Struct. 2013;96:57–63. 10.1016/j.compstruct.2012.09.005.Search in Google Scholar

[48] Kitipornchai S, He XQ, Liew KM. Continuum model for the vibration of multilayered graphene sheets. Phys Rev B. 2005;72(7):075443. 10.1103/PhysRevB.72.075443.Search in Google Scholar

[49] Nayfeh AH, Mook DT. Nonlinear oscillations. USA: John Wiley & Sons; 2008.Search in Google Scholar

[50] Azrar L, Benamar R, White RG. Semi-analytical approach to the non-linear dynamic response problem of S–S and C–C beams at large vibration amplitudes part I: general theory and application to the single mode approach to free and forced vibration analysis. J Sound Vib. 1999;224(2):183–207. 10.1006/jsvi.1998.1893.Search in Google Scholar

[51] Azrar L, Benamar R, White RG. A semi-analytical approach to the non-linear dynamic response problem of beams at large vibration amplitudes, part ii: multimode approach to the steady state forced periodic response. J Sound Vib. 2002;255(1):141. 10.1006/jsvi.2000.3595.Search in Google Scholar

[52] Ebrahimi F, Kokaba M, Shaghaghi G, Selvamani R. Dynamic characteristics of hygro-magneto-thermo-electrical nanobeam with non-ideal boundary conditions. Adv Nano Res. 2020;8(2):169–82. 10.12989/anr.2020.8.2.169.Search in Google Scholar

[53] Selvamani R, Loganathan R, Dimitri R, Tornabene F. Nonlocal state-space strain gradient wave propagation of magneto thermo piezoelectric functionally graded nanobeam. Curved Layer Struct. 2023;10(1):20220192. 10.1515/cls-2022-0192.Search in Google Scholar

[54] Selvamani R, Rexy J, Ebrahimi F. Two phase local/non local waves in a magneto thermo electrical composite nano beam reinforced with graphene oxide powder. Waves Random Complex Media. 2021;34(5):1–26. 10.1080/17455030.2021.1985745.Search in Google Scholar

[55] Anitha L, Rajalakshmi L, Selvamani R, Ebrahimi F. Forced nonlinear vibrations in a smart magneto-viscoelastic multiscale composite nanobeam in a humid thermal environment. Eng Trans. 2023;71(4):617–44. 10.24423/EngTrans.3114.20231121.Search in Google Scholar

[56] Selvamani R, Rexy J, Ebrahami F. Vibration in an electrically affected hygro-magneto-thermo-flexo electric nanobeam embedded in winkler-pasternak foundation. Mech Adv Compos Struct. 2021;8(2):401–14. 10.22075/macs.2021.22068.1311.Search in Google Scholar

[57] Selvamani R, Rubine L, Rexy J, Ebrahimi F. Dispersion analysis of electrically actuated hygro-magneto-thermo-flexo electric nanobeam embedded on silica aerogel foundation. Mater Phys Mech. 2022;50(1):1–19. 10.18149/MPM.5012022_1.Search in Google Scholar

Received: 2024-06-01
Revised: 2024-08-29
Accepted: 2024-09-09
Published Online: 2024-10-19

© 2024 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. Flutter investigation and deep learning prediction of FG composite wing reinforced with carbon nanotube
  3. Experimental and numerical investigation of nanomaterial-based structural composite
  4. Optimisation of material composition in functionally graded plates for thermal stress relaxation using statistical design support system
  5. Tensile assessment of woven CFRP using finite element method: A benchmarking and preliminary study for thin-walled structure application
  6. Reliability and sensitivity assessment of laminated composite plates with high-dimensional uncertainty variables using active learning-based ensemble metamodels
  7. Performances of the sandwich panel structures under fire accident due to hydrogen leaks: Consideration of structural design and environment factor using FE analysis
  8. Recycling harmful plastic waste to produce a fiber equivalent to carbon fiber reinforced polymer for reinforcement and rehabilitation of structural members
  9. Effect of seed husk waste powder on the PLA medical thread properties fabricated via 3D printer
  10. Finite element analysis of the thermal and thermo-mechanical coupling problems in the dry friction clutches using functionally graded material
  11. Strength assessment of fiberglass layer configurations in FRP ship materials from yard practices using a statistical approach
  12. An enhanced analytical and numerical thermal model of frictional clutch system using functionally graded materials
  13. Using collocation with radial basis functions in a pseudospectral framework to the analysis of laminated plates by the Reissner’s mixed variational theorem
  14. A new finite element formulation for the lateral torsional buckling analyses of orthotropic FRP-externally bonded steel beams
  15. Effect of random variation in input parameter on cracked orthotropic plate using extended isogeometric analysis (XIGA) under thermomechanical loading
  16. Assessment of a new higher-order shear and normal deformation theory for the static response of functionally graded shallow shells
  17. Nonlinear poro thermal vibration and parametric excitation in a magneto-elastic embedded nanobeam using homotopy perturbation technique
  18. Finite-element investigations on the influence of material selection and geometrical parameters on dental implant performance
  19. Study on resistance performance of hexagonal hull form with variation of angle of attack, deadrise, and stern for flat-sided catamaran vessel
  20. Evaluation of double-bottom structure performance under fire accident using nonlinear finite element approach
  21. Behavior of TE and TM propagation modes in nanomaterial graphene using asymmetric slab waveguide
  22. FEM for improvement of damage prediction of airfield flexible pavements on soft and stiff subgrade under various heavy load configurations of landing gear of new generation aircraft
  23. Review Article
  24. Deterioration and imperfection of the ship structural components and its effects on the structural integrity: A review
  25. Erratum
  26. Erratum to “Performances of the sandwich panel structures under fire accident due to hydrogen leaks: Consideration of structural design and environment factor using FE analysis”
  27. Special Issue: The 2nd Thematic Symposium - Integrity of Mechanical Structure and Material - Part II
  28. Structural assessment of 40 ft mini LNG ISO tank: Effect of structural frame design on the strength performance
  29. Experimental and numerical investigations of multi-layered ship engine room bulkhead insulation thermal performance under fire conditions
  30. Investigating the influence of plate geometry and detonation variations on structural responses under explosion loading: A nonlinear finite-element analysis with sensitivity analysis
Downloaded on 20.12.2025 from https://www.degruyterbrill.com/document/doi/10.1515/cls-2024-0013/html
Scroll to top button