Dynamic Analysis of a Uniform Microbeam Resting on a Nonlinear Foundation Considering Its Curvature Subjected to a Mechanical Impact and Electromagnetic Actuation

This study proposes an investigation into the nonlinear vibration of a simply supported, flexible, uniform microbeam associated with its curvature considering the mechanical impact, the electromagnetic actuation, the nonlinear Winkler–Pasternak foundation, and the longitudinal magnetic field. The governing differential equations and the boundary conditions are modeled within the framework of a Euler–Bernoulli beam considering an element of the length of the beam at rest and using the second-order approximation of the deflected beam and the Galerkin–Bubnov procedure. In this work, we present a novel characterization of the microbeam and a novel method to solve the nonlinear vibration of the microactuator. The resulting equation of this complex problem is studied using the Optimal Homotopy Asymptotic Method, employing some auxiliary functions derived from the terms that appear in the equation of motion. An explicit closed-form analytical solution is proposed, proving that our procedure is a powerful tool for solving a nonlinear problem without the presence of small or large parameters. The presence of some convergence-control parameters assures the rapid convergence of the solutions. These parameters are evaluated using some rigorous mathematical procedures. The present approach is very accurate and easy to implement, even for complicated nonlinear problems. The local stability near the primary resonance is studied.


Introduction
The flexural vibration of a microbeam under mechanical impact, electromagnetic actuation, a moving body, or a nonlinear foundation is interesting for basic research on vibration problems that has in view their practical interest.The main component of an electromagnetic actuator is the solenoid actuator.This is comprised of housing, a moving rod, coils, and a spring [1].The moving rod is made of a ferroelectric material that is surrounded by coils.When a voltage is applied through the coils, a current flow appears through the coils, which results in an electromagnetic force exerted on the moving rod, causing the rod to move.An inherent nonlinear relationship between displacement and current input exists for solenoids.Some previous studies have been the works of Abiala [2], which used the finite element method and Newmark's integration to obtain the dynamic response of a beam under uniformly distributed moving loads.The effect of electromagnetic actuation on the dynamics of a single-sided Hertzian contact-forced oscillator is explored near primary resonance and secondary resonance by Bichri et al. [3].
Muscia [4] obtained an engineering solution to manufacture electromagnetic linear actuators for moving rudders and fin stabilizers on military ships.The finite element method analysis helped them to find the most convenient solution for joining the pole pieces to the respective bars.Bichri et al. [5] considered the effects of combined DC and fast AC electromagnetic actuations on the dynamic behavior of an excited cantilever beam considering a single-mode approximation.Also, the influence of the fast AC actuation and the air gap on the nonlinear behavior of the system was analyzed.The time-stepping solution of the electromagnetic field model resulting in the magnetic force densities, computed from the Lorentz force and based on the space-time variation of the induced eddy current, was employed by Abba and Rachek [6].In this way, this framework was dedicated to predicting the wave magnetic force density for the mechanical deformations occurring in electromagnetic actuators.
A novel nonlinear lumped parameter model of the electromagnetic actuator is proposed by Zuo et al. [7] for the engine's idling or low-speed condition.The LuGre friction model is applied to describe the friction during the mover's reciprocating motion.A specific experimental scheme is proposed for validation in which a suspended ring-shaped iron is added as the actuator's load.Another novel electromagnetic actuator that has been designed to produce energy via vibrations while suppressing said vibrations, which can be installed in many kinds of mechanisms to improve machining accuracy and reduce energy loss, is employed by Wei et al. [8].Based on this prototype, the characteristics of the electromagnetic actuator are employed by means of measuring experiments.On the other hand, a novel, bistable, nonlinear electromagnetic actuator with an elastic boundary to enhance the actuation performance when controlled by a harmonic input signal is attempted by Zhang et al. [9].This model has an inclined spring, one end of which is supported by an elastic boundary, which produces bistable nonlinearity to realize the large-amplitude inter-well actuation responses.The governing equations are controlled by a harmonic input signal.
A self-sensing electromagnetic actuator for vibration control of flexible structures is studied by Verma et al. [10].This scheme does not require any modifications in the design of the electromagnetic actuator and has the advantage that the voltage drop across the coil can be measured in relation to motion.Seebacher et al. [11] introduced a pseudodensity topology optimization method in nonlinear electromagnetism that follows the Solid Isotropic Material with Penalization procedure.The adjoint technique is used to determine the sensitivity of the objective functions.An applied linear characteristic of the ferromagnetic material results in a small or wrongly actuated electromagnetic dimension.The mathematical model of lateral vibration of a cantilever pipe discharging fluid using an electromagnetic actuator of the transformer type in the form of a partial differential equation coupled with two ordinary differential equations is presented by Szmidt et al. [12].The system is discretized in a set of ordinary differential equations and is solved numerically.The critical velocity and the amplitude of post-critical vibrations are determined and validated experimentally.
Al Bakri et al. [13] considered a sliding mode control method to design the closed-loop coil voltage such that the nonlinear system of an electromagnetic actuator slides along a certain sliding surface.The reliability of the system was approved under a wide range of steady-state actuator position dispersions and coil voltage disturbances using the Monte Carlo simulation method.
A novel miniaturized and modular dual-axis electromagnetic actuator for several applications is proposed by Mansour et al. [14].The modularity of this actuator gives an extra privilege to connect several modules together, providing more potential applications, and can be used to generate a quadrupled robot.The robot shows a high level of manipulability and has a light weight and a small size.A new concept of two-way electromagnetic actuation with the help of magnetic damping, which can extend its arm on both sides to facilitate active suspension mechanisms in both humps and potholes, has been derived by Prajwal et al. [15].Catia V5 software was used to design and simulate the actuator model.Repinaldo et al. [16] investigated the parameter identification and application of the active vibration control technique within a structure with two degrees of freedom using an electromagnetic actuator.The neuro-fuzzy control theory was applied to reduce the amplitude of vibration in the structure when subjected to an impulsive force.
Konig et al. [17] explored the design of a hysteresis-compensated self-sensing algorithm with low computational effort and integrator-based direct inductance.The measurement technique was applied for the resource-efficient estimation of the incremental inductance.Also, a modified Prandtl-Ishlinski approach was used for modeling this hysteretic behavior.Yang et al. [18] examined the broadband active vibration of a bistable, nonlinear electromagnetic actuator with an elastic boundary.A new feedback control law was proposed to control the actuator's input curves to significantly attenuate the broadband vibration transmissibility from the base excitation to the actuator's movements.The input-to-state stability of the control law for any non-negative control weights was proven.
Zhang et al. [19] proposed the physical model of the actuator via Ampere's loop law, the continuity principle of magnetic flux, and Ohm's law of magnetic circuits.To enhance hysteresis model accuracy, the magnetic flux leakage coefficients of the actuator with different positions were obtained.The experimental results demonstrate that the trajectory and tracking errors compared with the feedback control and the unoptimized hysteresis method feedforward control of the optimized hysteresis were substantially reduced.Lim et al. [20] investigated the stochastic vibration responses of a bistable electromagnetic actuator with an elastic boundary controlled by Gaussian white noise and low-pass-filtered stochastic noise.This actuator used an oblique spring and another linear spring to realize the bistability and the elastic boundary, respectively.
Jorshari et al. [21] investigated the vibration of single-walled boron nitride nanotube induced by a moving nanoparticle based on the Rayleigh beam model and used the Pasternak substrate to model the elastic medium.Baroudi et al. [22] predicted quantitatively the motion of a soft-beam clamped Euler's Elastica by adapting the Hencky technique in a regime of large displacement and deformations.A nonlinear planar beam with stretch and shear deformation is described by Ren et al. [23].The static equations of equilibria are discretized by the finite difference method.The nonlinear vibration of an elastic string with large amplitude and large curvature has been investigated by Zhao et al. [24].The modified compeer novel form method and the finite difference scheme are used to calculate the critical parameters, the time history diagram, configuration, total length, or fundamental frequency.
The present study is devoted to the nonlinear forced vibration of simply supported, flexible, uniform microbeams.The nonlinearities of the equations are caused by the curvature of the beam, mechanical impact, the electromagnetic actuator, the nonlinear Winkler-Pasternak elastic foundation, and the longitudinal magnetic field.The governing nonlinear differential equations and the corresponding boundary conditions are established for a Euler-Bernoulli beam considering an element of the length of the beam at rest and applying the second-order approximation of the deflected beam and the Galerkin-Bubnov procedure.The two differential equations are reduced to only one nonlinear differential equation, and this equation is solved by means of the Optimal Homotopy Asymptotic Method (OHAM) with the help of so-called auxiliary functions obtained from the terms that appear in the governing equation.A very accurate solution is obtained using a moderate number of convergence-control parameters.
Our original approach is a very powerful tool for solving a nonlinear problem without the presence of any small parameters in the governing equation or in the boundary conditions.The presence of some convergence-control parameters assures the rapid convergence of the approximate solutions after the first iteration.The convergence-control parameters are evaluated using some rigorous mathematical procedures.In general, we have great freedom to select both auxiliary functions and the number of convergence-control parameters.This approach has proven to be very accurate, simple, and easy to implement for any complicated nonlinear problems.The local stability near the primary resonance is studied using the Homotopy Perturbation Method, appealing to the Jacobian and eigenvalues of the system.

Formulation of the Problem
In Figure 1 is represented the physical model of a simply supported microbeam of length L subjected to a mechanical impact by the transversal constant force, an electromagnetic load V DC , the nonlinear Winkler-Pasternak elastic foundation F WP , and longitudinal magnetic field F m .parameters are evaluated using some rigorous mathematical procedures.In general, we have great freedom to select both auxiliary functions and the number of convergence-control parameters.This approach has proven to be very accurate, simple, and easy to implement for any complicated nonlinear problems.The local stability near the primary resonance is studied using the Homotopy Perturbation Method, appealing to the Jacobian and eigenvalues of the system.

Formulation of the Problem
In Figure 1 is represented the physical model of a simply supported microbeam of length L subjected to a mechanical impact by the transversal constant force, an electromagnetic load VDC, the nonlinear Winkler-Pasternak elastic foundation FWP, and longitudinal magnetic field Fm.The transversal and longitudinal displacements are W(x,t) and u(x,t), respectively.The Young's modulus E, the mass density ρ, and the cross-sectional area A of the microbeam are supposed to be constant.
Considering an element of length dx of the microbeam at rest, the corresponding deflection can be written as The unit vector parallel to the deflected element is of the form The length of the microbeam is variable during the motion, such that the tension becomes The governing equations are The transversal and longitudinal displacements are W(x,t) and u(x,t), respectively.The Young's modulus E, the mass density ρ, and the cross-sectional area A of the microbeam are supposed to be constant.
Considering an element of length dx of the microbeam at rest, the corresponding deflection can be written as The unit vector parallel to the deflected element is of the form The length of the microbeam is variable during the motion, such that the tension becomes The governing equations are where I is the moment of inertia of the microbeam cross-section, F mi is the impact force, and F em is the electromagnetic actuation.
Neglecting the terms of the forms ∂x 2 , and ∂ 4 u ∂x 4 , from Equation (1), one can obtain Micromachines 2024, 15, 969 then, from Equations ( 6) and ( 7), it follows that e = ∂u ∂x Equation ( 4) can be rewritten, taking into account Equations ( 1)-( 3) and ( 7), as Usually, the longitudinal inertial term ∂ 2 u ∂x 2 can be neglected, such that, from Equation ( 9), it is clear that e is a constant, e = C, and, therefore, from Equation ( 8), it holds that By integrating the last equation with the boundary conditions the constant C is given by Substituting Equation (12) into Equation (5), we obtain The term F mi represents the mechanical impact of load F 0 : where v is the speed of the load, and δ is the Dirac function.The vibration of the microbeam depends on the action of the DC voltage source, such that electromagnetic actuation can be considered as in which C 0 is the capacitance of the actuator, d is the gap width, and V DC is the voltage.Expression ( 15) can be simplified in the following form: The nonlinear elastic Winkler-Pasternak and shear foundation of the beam is represented by where K 1 and K 3 are linear and nonlinear coefficients, respectively, and K s is the shear coefficient of the elastic foundation.
The longitudinal magnetic field is given by magnetic permeability, area A, and the component H x in the x-direction.
According to the governing Equation ( 13), the term that defines the curvature of the microbeam can be rewritten in the following form: Inserting Equations ( 14)-( 19) into Equation ( 13), one retrieves The above equation is obtained for a simply supported microbeam such that the boundary conditions are To express the governing Equations ( 20) and ( 21) in nondimensional form, the following parameters are defined: Omitting the bar, Equation ( 20) can be rewritten in nondimensional form as The solution of Equation ( 23) can be assumed to be of the following form: where, taking account of the boundary conditions (21), the eigenfunction X(x) can be expressed as which satisfies the orthonormality condition 1 0 Substituting Equation ( 24) into (23), multiplying with X(x), integrating on the domain [0, 1], using the expression and then using the Galerkin-Bubnov procedure, Equation (23) will be written as follows: ..
T + ω 2 T + aT 3 + bT 5 = f sin πvt (28) where the dot denotes the derivative with respect to time and dX dx The initial conditions for Equation (28) are Equation ( 28) with the initial conditions (30) is a second-order nonlinear forced differential equation of the fifth order and thus is very difficult to analytically solve by an exact solution.In what follows, for Equations ( 28) and (29), the OHAM is applied to study the nonlinear forced vibration.
It should be mentioned that the coefficients ω, a, b, and f include the effect of the curvature, mechanical impact, electromagnetic actuation, and nonlinear elastic foundation.

The Optimal Homotopy Asymptotic Method (OHAM)
The OHAM was proposed for the first time in 2008 by Marinca and Herisanu [25][26][27][28].To apply the OHAM, we consider the following nonlinear differential equation: where L[T(t)] is a linear operator, t denotes independent variable, T(t) is an unknown function, N[T(t)] is a nonlinear operator, D is the domain of interest, and B is a boundary operator.By means of the OHAM, one first constructs a family of equations: where p ∈ [0, 1] is an embedding parameter, H is a nonzero auxiliary function for p ̸ = 0 and H(t,0,C i ) = 0, T is the unknown approximate solution, and C i , i = 1,2,. ..,n are n unknown parameters which will be determined later.Obviously, when p = 0, it holds that from which can be determined the initial approximation T 0 (t).If T(t) = T 0 (t)+ pT 1 (t, C 1 , . . . ,C n ), then, from Equation (32), the first approximation T 1 (t, C 1 , C 2 , . . . C n ) can be determined for p = 1: where H(t, C 1 , C 2 , . .., C n ) is an arbitrary auxiliary function which can be determined such that the products HN[T 0 (t)] and N[T 0 (t)] are of the same shape.The parameters C 1 , C 2 ,. ..., C n which appear in the first-order approximate solution T 1 can be determined in many ways, using the least square method, the Galerkin method, the collocation method, the Ritz method, or by minimizing the square residual error as where R is given by Equation (31): The unknown parameter C i can be identified optimally from the condition With these parameters known (namely, the convergence-control parameters), the approximate analytical solution T(t) is well determined.
It should be emphasized that the auxiliary function and the number n of the convergencecontrol parameters are not unique.
The OHAM is a very efficient procedure since any nonlinear differential equation is reduced to only two linear differential equations, and, on the other hand, the convergencecontrol parameters are optimally determined.

Application of Optimal Homotopy Asymptotic Method to Study the Nonlinear Vibration of the Microbeam
In what follows, we apply the OHAM to obtain an approximate solution for Equations ( 28) and (30).Making the transformations where Ω is the unknown natural frequency of the beam, Equation (28) becomes and corresponding initial conditions are where the prime denotes the derivative with respect to τ.The linear and nonlinear operators are, respectively, The initial approximation is determined from Equation (33), which is rewritten as The first approximation θ 1 (τ) is obtained from Equation (34) as where N[θ 0 (τ)] is determined by inserting Equation (43) into the nonlinear operator (41): The coefficients from the above equations are The auxiliary function H from Equation ( 44) is chosen such that the product between HN[θ 0 (τ)] and N[θ 0 (τ)] is of the same form.For example, we can choose these auxiliary functions in the following expressions: Using only expression (47), Equation ( 44) can be rewritten after some manipulations as follows: No secular term in Equation (51) requires that the coefficient of cosτ be zero.From this condition, one can determine the natural frequency: 4  16 The solution of linear Equation ( 51) is (cos τ − cos 3τ)+ The approximate solution T is obtained from Equations (38), (41) and (55):
The OHAM can be applied to any nonlinear differential equations without a supplementary hypothesis, such as, for example, the existence of a small parameter in the differential equation or in the boundary/initial conditions.

Numerical Results
The efficiency of the OHAM can be proved through numerical investigations for the following case considering α = 0.007, β = 0.001, γ 1 = 0. To validate the obtained analytical solution, the analytical results are compared with numerical integration results obtained by means of a fourth-order Runge-Kutta approach.The explicit analytical solution benefits from the optimal values identified for the convergence-control parameters and is plotted in Figure 2, demonstrating the high accuracy of the employed method.In Figure 2, we present the approximate solution (54) in comparison with the numerical integration results.It can be observed that the error between the approximate analytical solution obtained by the OHAM and the solution obtained by numerical integration is very small, these solutions being in an excellent agreement.In this way, the results obtained by the OHAM are validated through this comparison, and an excellent accuracy of the analytical solution is proved.
The OHAM can be applied to any nonlinear differential equations without a supplementary hypothesis, such as, for example, the existence of a small parameter in the differential equation or in the boundary/initial conditions.

Numerical Results
The efficiency of the OHAM can be proved through numerical investigations for the following case considering α = 0.007, β = 0.001, γ1 = 0. To validate the obtained analytical solution, the analytical results are compared with numerical integration results obtained by means of a fourth-order Runge-Kutta approach.The explicit analytical solution benefits from the optimal values identified for the convergence-control parameters and is plotted in Figure 2, demonstrating the high accuracy of the employed method.In Figure 2, we present the approximate solution (54) in comparison with the numerical integration results.It can be observed that the error between the approximate analytical solution obtained by the OHAM and the solution obtained by numerical integration is very small, these solutions being in an excellent agreement.In this way, the results obtained by the OHAM are validated through this comparison, and an excellent accuracy of the analytical solution is proved.

Analysis of the Stability of the Steady-State Motion near the Primary Resonance
The procedure applied in this section proposes to distinguish between two time-scales by associating a separate independent variable with each one and using the Homotopy Perturbation Method [29].We reconsider the primary resonance for Equation (28): with λ being a detuning parameter and p ∈ [0, 1] the embedding parameter.As in Section 3, In accordance with the Homotopy Perturbation Method, for p = 0, any nonlinear differential equation reduces to a linear differential equation, and, for p = 1, one obtains the original equation, such that the original Equation ( 28) is rewritten as ..
The variable T can be written as a power series of p: where S i are the solutions of the i-th order homotopy equations.
From Equations ( 57), ( 58) and ( 61), the homotopy equations of order p i , i = 0,1 are, respectively, The solution of the linear Equation (62) with partial derivatives is where A and B are functions depending on η.Now, substituting Equation (64) into Equation ( 63) and avoiding the secular terms yields The equilibrium points correspond to the periodic solution of Equation (28).To determine them, set the following: In this way, the equilibrium points can be obtained from the algebraic equations It is possible only in the case where              It can be seen that there is a single equilibrium point B e in every case, except in situations where there are two equilibrium points B e in a small neighborhood of f = 0 for a = 0.5 and b = 1 (Figure 4).
The stability of the steady-state motion is determined by the eigenvalues of the Jacobian matrix obtained from Equations (66) and (67) (the Routh-Hurwitz criteria): [J] = The eigenvalues of the Jacobian matrix are obtained from the characteristic equation det( where I is the unity matrix of the second order, and λ is the eigenvalue of the Jacobian matrix.The characteristic equation can be rewritten as We have the following possible cases: Case 6.1 If ab ≥ 0 and a 2 + b 2 ̸ = 0, then the eigenvalues are purely imaginary with opposite signs of their imaginary parts.This is a center, such that there is no net motion towards or away from the equilibrium point.Case 6.2 If ab < 0 and B 2 e = − 18a 25b or B 2 e = − 6a 5b , then λ 1 = λ 2 = 0.The situation is similar to the above.Case 6.3 If ab < 0 and − 18a 25b < B 2 e < − 6a 5b , then λ 1 and λ 2 are real, and λ 1 = −λ 2 .This case corresponds to a saddle point.Case 6.4 If ab < 0 and B 2 e < − 18a 25b or B 2 e > − 6a 5b , then the eigenvalues are purely imaginary, as in the situation of case 6.1.

Results and Discussion
In Figures 6-9, the effects of different parameters on the vertical displacements of the Euler-Bernoulli beam are shown.
where The eigenvalues of the Jacobian matrix are obtained from the characteristic equation where I is the unity matrix of the second order, and λ is the eigenvalue of the Jacobian matrix.The characteristic equation can be rewritten as 0 256 We have the following possible cases: Case 6.1 If ab ≥ 0 and a 2 + b 2 ≠ 0, then the eigenvalues are purely imaginary with opposite signs of their imaginary parts.This is a center, such that there is no net motion towards or away from the equilibrium point.

Results and Discussion
In Figures 6-9, the effects of different parameters on the vertical displacements of the Euler-Bernoulli beam are shown.From Figure 6, we deduce that, if the curvature parameter β increases, then the abscissa of the intersection point of the vertical displacement and the horizontal axis decreases.Different values of β are taken as 0.01, 0.02, and 0.03.The effect of the nonlinear Winkler-Pasternak parameter K 3 on the vertical displacement is shown in Figure 7.This figure shows that, if the parameter increases from K 3 = 0.4 to K 3 = 0.6, then the abscissa of the intersection point of the vertical displacement and the horizontal axis increases.The effect of the gap parameter α is illustrated in Figure 8; the abscissa of the intersection point of the vertical displacement and the horizontal axis decreases with the growth of the gap parameter.From Figure 9, it follows that, if the electromagnetic parameter γ 1 increases, then the above-mentioned abscissa increases even more substantially for the time t > 14.From Figure 6, we deduce that, if the curvature parameter β increases, then the abscissa of the intersection point of the vertical displacement and the horizontal axis decreases.Different values of β are taken as 0.01, 0.02, and 0.03.The effect of the nonlinear Winkler-Pasternak parameter K3 on the vertical displacement is shown in Figure 7.This figure shows that, if the parameter increases from K3 = 0.4 to K3 = 0.6, then the abscissa of the intersection point of the vertical displacement and the horizontal axis increases.The effect of the gap parameter α is illustrated in Figure 8; the abscissa of the intersection point   From Figure 6, we deduce that, if the curvature parameter β increases, then the abscissa of the intersection point of the vertical displacement and the horizontal axis decreases.Different values of β are taken as 0.01, 0.02, and 0.03.The effect of the nonlinear Winkler-Pasternak parameter K3 on the vertical displacement is shown in Figure 7.This figure shows that, if the parameter increases from K3 = 0.4 to K3 = 0.6, then the abscissa of the intersection point of the vertical displacement and the horizontal axis increases.The effect of the gap parameter α is illustrated in Figure 8; the abscissa of the intersection point   From Figure 6, we deduce that, if the curvature parameter β increases, then the abscissa of the intersection point of the vertical displacement and the horizontal axis decreases.Different values of β are taken as 0.01, 0.02, and 0.03.The effect of the nonlinear Winkler-Pasternak parameter K3 on the vertical displacement is shown in Figure 7.This figure shows that, if the parameter increases from K3 = 0.4 to K3 = 0.6, then the abscissa of the intersection point of the vertical displacement and the horizontal axis increases.The effect of the gap parameter α is illustrated in Figure 8; the abscissa of the intersection point It can be observed that the longitudinal inertia term u(x,t) can be neglected only if E/ρ satisfies the condition (83).

Conclusions
According to the present results, the oscillatory behavior of a simply supported flexible uniform microbeam is studied.The Euler-Bernoulli microbeam is subjected to a mechanical impact using the Dirac delta function, an electromagnetic actuation, a nonlinear Winkler-Pasternak foundation, and a longitudinal magnetic field.In the governing differential equations, the curvature of the beam, the second-order approximation of the deflected beam, and the boundary conditions are taken into consideration.The OHAM is applied to solve the complex nonlinear differential equation, introducing so-called auxiliary functions and some convergence-control parameters without supplementary hypotheses.These parameters are optimally determined by rigorous mathematical procedures.The proper procedure leads to a very accurate solution using only the first iteration.It should be emphasized that the nonlinear dynamical system is reduced to only two linear differential equations.
The main aspects of novelty presented in this work, which are not present in previous related research, are summarized as follows: (a) The combination of the effects of the curvatures, mechanical impact, electromagnetic actuation, nonlinear elastic Winkler-Pasternak foundation, and longitudinal magnetic field.This combination is ensured by the new dynamical model which is considered in the study.(b) The Optimal Homotopy Asymptotic Method is used to solve a very complex differential equation by means of some auxiliary functions, providing novel solutions.(c) The values of the convergence-control parameters are optimally determined, ensuring a high degree of accuracy.(d) It is proved that the longitudinal inertia term can be neglected, but only if E/ρ satisfies the condition (83).(e) Local stability for the primary resonance is studied by means of the Homotopy Perturbation Method.
The effect of some physical parameters is highlighted.The local stability depends on the nature of the solutions to the characteristic equation, and different cases are discussed.The local stability near the primary resonance is achieved by using the Homotopy Perturbation Method, the variable expansion method, equilibrium points, and the Jacobian matrix.The effects of different parameters on the vertical displacements of the beam are considered.

Figure 1 .
Figure 1.Schematic diagram of the microbeam and electromagnetic microactuator.

Figure 1 .
Figure 1.Schematic diagram of the microbeam and electromagnetic microactuator.

5bB 5 e + 6aB 3 e
− 8 f = 0 (70) In Figures 3-5 is depicted the equilibrium point B e in respect to the parameter f for different values of the parameters a and b: a = 0.25, b = 0.3 (Figure 3); a = 0.5, b = 1 (Figure 4); a = −1, b = 0.25 (Figure 5).It is possible only in the case where depicted the equilibrium point Be in respect to the parameter f for different values of the parameters a and b: a = 0.25, b = 0.3 (Figure 3); a = 0.5, b = 1 (Figure 4); a = −1, b = 0.25 (Figure 5).
depicted the equilibrium point Be in respect to the parameter f for different values of the parameters a and b: a = 0.25, b = 0.3 (Figure 3); a = 0.5, b = 1 (Figure 4); a = −1, b = 0.25 (Figure 5).

Case 6 . 2
If ab < 0 and  = − or  = − , then λ 1= λ 2 = 0.The situation is similar to the above.Case 6.3 If ab < 0 and − <  < − , then λ1 and λ2 are real, and λ1 = −λ2.This case corresponds to a saddle point.Case 6.4 If ab < 0 and  < − or  > − , then the eigenvalues are purely imaginary, as in the situation of case 6.1.