Research Article Volume 7 Issue 1
School of Metallurgy and Material Engineering, Iran University of Science and Technology, Iran
Correspondence: Saeed Shahsavari, School of Metallurgy and Material Engineering, Iran University of Science and Technology, Tehran, Iran, Tel 00989376081167
Received: June 19, 2023 | Published: June 28, 2023
Citation: Shahsavari S, Boutorabi SMA. Nonlinear nonlocal damped free and forced vibrations of piezoelectric SWCNTs under longitudinal magnetic field due to surface effects using a two steps perturbation method. MOJ App Bio Biomech. 2023;7(1):88-99. DOI: 10.15406/mojabb.2023.07.00180
In the present work, damped free and forced vibrations of single-walled piezoelectric carbon nanotubes under longitudinal magnetic field due to surface effects surrounded on a non-linear viscoelastic medium using the nonlocal Euler-Bernoulli beam theory and multiple time scales method are investigated. Lorentz force equation is used to obtain the vertical force due to the applied voltage to the system. The surface effects as well as a combinational non-linear viscoelastic foundation are considered, and finally, the dynamic equilibrium equations are used, and non-linear equations of motion are extracted. In the following, the Galerkin and multiple time scales methods are used, and finally, analytical solutions are extracted as the non-linear free and forced vibrational responses of the system. The relevant coefficients of the extracted analytical solutions are discovered for two both simple support (S-S) and clamped (C-C) boundary conditions. In the following, , and the effects of the different parameters such as non-local parameter as well as electric-magnetic fields, effect of hardness-linear damping parameters of nonlinear considered viscoelastic foundation, applied magnetic field, base modes for different forms considering surface effects, and etc. will be studied. As one the results of this study, the presence of a non-local parameter has increased the curvature deviation to the right and the stiffening effect. In other words, the non-local parameter is a factor to increase the nonlinear effect of the system. Also, it is predictable that as the load affect position moves away from the center of the single-walled piezoelectric carbon nanotube toward the supports, the amplitude of the dynamic response decreases significantly, and this relative reduction is greater for the C-C boundary condition than for the S-S boundary condition. It is also important to note that the location of the load has no effect on the rate of deviation of the curve peak, and the degree of nonlinearity of the vibrational response of the system.
Keywords: nonlocal euler-bernoulli beam theory, non-linear non-local damped forced vibration, viscoelastic piezoelectric carbon nanotubes, magnetic field, non-linear viscoelastic foundation, surface effects
Having unique physical and chemical properties such as extremely high tensile strength and rigidity, combined with very low density, makes CNTs an essential part of all industrial fields. Nanotubes are critical members of the design and fabrication of nanoelectromechanical systems.1 Microelectromechanical and nanoelectromechanical systems are the essential applications of piezoelectric generators.2 Surface effects are other phenomena that occur due to the static balance of atoms on the surface.3 These effects on elastic bodies are divided into three categories: surface tension, residual surface stress, and surface density.4 Indeed, surface effects will play a serious role when the surface-to-volume ratio increases.5–6 Finally, surface effects may have an important role in the vibration analysis of nanotubes.7–8
The damped forced vibration of SWCNTs was analyzed using a new shear deformation beam theory.9 In this new shear deformation beam theory, there was no need to use any shear correction factor, and also the number of unknown variables was the only one that was similar to the Euler-Bernoulli beam hypothesis. The torque effect of an axial magnetic field on a functionally graded (FG) nano-rod has been studied to capture size effects under magnetic field using Maxwell’s relation.10 As an important result of this study could be mentioned that an FG nano-rod model based on the nonlocal elasticity theory behaves softer and has a smaller natural frequency. Investigation of free vibration of viscoelastic nanotube under longitudinal magnetic field also was studied in recent years.11 According to a deep investigation on the natural frequencies and effect of different parameters such as the nonlocal parameter, structural damping coefficient, material length scale parameter, and the longitudinal magnetic field, the results of this research may be helpful for understanding the potential applications of nanotubes in Nano-Electromechanical System. Many studies in the literature are focused on the investigation of the vibration analysis of nanostructures in an elastic medium. The vibrational response of a SWCNT member considered in an elastic medium to transport a viscous fluid,12 Investigating the vibration of a DWCNT embedded in an elastic medium despite initial axial forces,13,14 Investigating the vibration of nanotubes in conditions where the surrounding medium is considered elastic,15 investigating the nonlinear free vibration of a DWCNT member considering the von Karman assumption in order to apply nonlinear effects and nonlinear,16 nonlinear vibration of a DWCNT by C-C boundary condition in an elastic medium using the nonlinear van der Waals forces, and also the von Kármán geometric to consider nonlinearity,17 Investigating the forced vibration of a DWCNT with the potential to carry a moving nanoparticle,18 investigation of the nonlinear vibration of an MWCNT in thermal environments,19 vibration analysis of an MWCNT by thermal effects, and also considering the size effects on the large amplitude,20 studying free vibration of SWCNT with elastic effects for different boundary conditions,21 vibration of SWCNT with an elastic medium due thermal conditions,22 thermal-mechanical effects consideration to study vibration and buckling instability of a SWCNT that carry out fluid, and also is rested on an elastic foundation,23 electro-thermo-mechanical effects consideration to study vibration of boron nitride nanorod that is in elastic medium by non-uniform and non-homogeneous properties,24 buckling analysis of SWCNT on a viscoelastic foundation for different boundary conditions,25 buckling investigation of SWCNT under thermal effects that is in an elastic medium with one elasticity parameter,26 and buckling behavior of single-walled CNT considering thermal effect in an elastic medium27 were studied using non-local theory of elasticity. another theory that can be considered as size-dependent continuum theories, and also can be used to study the electro-thermal transverse vibration of CNT in an elastic medium was presented based on the non-local shell theory as well as piezoelasticity properties of system.28 Continuum shell theory by considering nonlocal effects was used to study free vibration of single-walled carbon nanotubes by non-homogenous elastic medium.29 Using multiple elastic models of beams as well as continuum mechanics theory, non-linear free vibration of multi-walled carbon nanotubes was studied.30 Nonlinear vibration as well as thermal stability and of pre/post buckling due to temperature effects, and also size-dependent FG beams designed on an elastic foundation was studied based on the modified couple stress theory.31 The multiple time scales method, as a perturbation method, is known as an efficient technique in nonlinear differential equations analysis. Multiple scale method was used to study free and forced vibration of beams as well as DWCNT on an elastic medium by considering geometric nonlinearity and S-S boundary condition.32–34 Non-local elasticity theory was used to study tensioned nanobeam nonlinear vibration with considering various boundary conditions.35,36 The non-local Euler–Bernoulli theory of elasticity37 can be considered as a theory to modify nanoscale as well as nanostructers models that greatly has been used in modelling carbon nanotubes due to the size-effects in nanoscale behavior. Several references in the literatures can be found on the basis of the nonlocal theory of elasticity.38,39 In some of the recent relevant papers, the nonlinear frequency response of single-walled carbon nanotubes to primary resonance has been studied based on the nonlocal Euler-Bernoulli beam theory.40–48 The nonlocal Euler–Bernoulli elasticity theory also has been greatly used to study Nonlocal instability of cantilever piezoelectric carbon nanotubes by considering surface effects subjected to axial flow,49 Linear free vibration analysis of piezoelectric SWCNTs that the results are extracts considering linear part of a new general nonlinear viscoelastic foundation without considering external force,50 Vibration analysis of piezoelectric nanowires with surface and small scale effects,51 Flexoelectric and surface effects on size-dependent flow-induced vibration and instability analysis of fluid-conveying nanotubes based on flexoelectricity beam model,52 and Size-dependent nonlinear vibration of functionally graded composite micro-beams reinforced by carbon nanotubes with piezoelectric layers in thermal environments.53 Consequently, due to lots of the presented works in the literature, the objective of the present work is to study damped free and forced vibration of single-walled piezoelectric carbon nanotubes under longitudinal magnetic field considering surface effects resting on a non-linear viscoelastic foundation based on the nonlocal Euler-Bernoulli beam theory as well as multiple time scales and Galerkin methods. A theoretical and numerical study on non-linear nonlocal free and forced vibration responses are performed. This study, considering surface effects, develops a theoretical response for the non-linear nonlocal vibration of piezoelectric carbon nanotubes located on a non-linear viscoelastic foundation, and attenuation factors due to the external magnetic field and piezoelectric voltage Using Galerkin and multiple time scales methods. Surface effects, which harden the surface of nanotubes, are created by the formation of homogeneous masses by Van forces in the waltz. In the following, considering different boundary conditions of the system (S-S and C-C), the relevant coefficients of the developed analytical response are determined. Finally, the vibrational response parameters such as Amplitude-frequency response curves of nonlinear forced vibration, and etc. will be generally investigated for different boundary conditions, and the effects of the different parameters such as non-local parameter and electric-magnetic fields, Effect of hardness-linear damping parameters of viscoelastic foundation, applied magnetic field, base modes for different forms considering surface effects, and etc. will be studied.
Scheme of the considered system
Figure 1 shows structural schematics of a piezoelectric single-walled carbon nanotube with an inner diameter , outer diameter , thickness , length , mass per unit length , and elastic modulus , subjected to voltage and a magnetic field of magnitude . This nanotube, subjected to a harmonic external point load, resting on a nonlinear viscoelastic foundation. Figure 2 shows discrete model of the single walled carbon nanotube by coordinate system.
Figure 1 Scheme of the single-walled piezoelectric carbon nanotubes under longitudinal magnetic field subjected to voltage and resting on nonlinear viscoelastic foundation with surface layer.
Non-local effects
Generally, nanoscale or microstructure analysis methods can be classified into three categories: molecular simulation, non-classical continuous environment, and multi-scale.54 The molecular simulation methods, which simulated atoms and their bonds, investigate nanostructures' behavior more profound than other methods. The high computational volume of this method makes it possible to simulate only small nanostructured dimensions using this method.55,56 Molecular methods are divided into quantum mechanics, molecular dynamics, and molecular mechanics.57 Complex energy functions are used to explain the motion of elementary particles in quantum mechanics. However, this method is associated with computational limitations due to the complexity of the equations used. In the molecular dynamics method, the motion of atoms is studied concerning their adjacent bonds. Atoms are known as rigid particles that are affected by the potential field of neighbouring atoms. The last one, the molecular mechanical method, is based on the displacement of interatomic bonds with beam and spring elements. Potential energy can express the field of interatomic forces in atomic structures. Equalizing the beam and spring elements' potential energy and strain energy, the mechanical properties are obtained equally. In classical theory, the stress state at a particular point can be calculated if the strain state is known at that point. Whereas, in non-classical theory, such as non-local tension theory, the strain state must be known throughout the range to obtain the stress state at a particular point. All other proposed theories based on non-classical continuum mechanics, including strain gradients and coupled stress theories, relate to two-dimensional elasticity.58–60 Therefore, nonlocal stress tensor can written as follows:61,62
σ(X)=∫Vk(|X′−X|, τ)T(X′)dVσ(X)=∫Vk(|X′−X|, τ)T(X′)dV (1)
T(X)=C(X):ε(X)T(X)=C(X):ε(X) (2)
That σ(X)σ(X) is nonlocal stress tensor at point XX , (|X′−X|, τ)(|X′−X|, τ) is nonlocal modulus, is a material constant, T(X)T(X) is classical stress tensor, is ε(X)ε(X) strain tensor and C(X)C(X) is elasticity tensor. Because of that solving of the integral constitutive Equation (1) is complicated, a simplified differential for is greatly used as follows:
T=(1−μ∇2)σ, μ=τ2l2T=(1−μ∇2)σ, μ=τ2l2 (3)
Where μμ is non-local coefficient and is ∇2∇2 Laplacian operator. Also, ττ and ll are material constant.61 Finally, nonlocal stress tensor for a beam must satisfy equation (4):
σ(X)−(τl)2∂2σ(X)∂X2=Eε(X)σ(X)−(τl)2∂2σ(X)∂X2=Eε(X) (4)
Where is a constant that known as the classical elasticity modulus.
Now, using the existing stress relations in the non-local theory, the relationship between the local bending torque Ml(x,t)Ml(x,t) and the non-local torque Mnl(x,t)Mnl(x,t) at the cross section of the nanotube is written as follows63:
(1−(τl)2∂2∂x2)Mnl(x,t)=Ml(x,t)(1−(τl)2∂2∂x2)Mnl(x,t)=Ml(x,t) (5)
In equation (5), μ μ is the parameter of the nanotube size scale, which is also known as the non-local parameter.
External force related to Lorentz force and the effect of the magnetic fieldAssume →U=(u,v,w)→U=(u,v,w) as the displacement field of the nanostructure; the Lorentz force flfl is equal to64:
→fl=η (→J×→H)→fl=η (→J×→H) (6)
Where J stands for the current density and H is the magnetic field strength vector in the environment. Also, in equation (6):
→J=∇×→h→J=∇×→h (7)
That:
→h=∇×(→U×→H)→h=∇×(→U×→H) (8)
Therefore, by considering longitudinal magnetic field:
→h=−Hx(∂v∂y+∂w∂z)→i+Hx∂v∂x→j+Hx∂w∂x→k→h=−Hx(∂v∂y+∂w∂z)→i+Hx∂v∂x→j+Hx∂w∂x→k (9)
→J=Hx(−∂2v∂x∂z+∂2w∂x∂y)→i−Hx(∂2v∂y∂z+∂2w∂x2+∂2w∂z2)→j+Hx(∂2v∂x2+∂2v∂y2+∂2w∂y∂z)→J=Hx(−∂2v∂x∂z+∂2w∂x∂y)→i−Hx(∂2v∂y∂z+∂2w∂x2+∂2w∂z2)→j+Hx(∂2v∂x2+∂2v∂y2+∂2w∂y∂z) (10)
Finally, the Lorentz force is:
→fl=fx→i+fy→j+fz→k=η[0→i+Hx2(∂2v∂x2+∂2v∂y2+∂2w∂y∂z)→j+Hx2(∂2w∂x2+∂2w∂y2+∂2v∂y∂z)]→fl=fx→i+fy→j+fz→k=η[0→i+Hx2(∂2v∂x2+∂2v∂y2+∂2w∂y∂z)→j+Hx2(∂2w∂x2+∂2w∂y2+∂2v∂y∂z)] (11)
Therefore, along the x,y,zx,y,z directions:
fx=0fx=0
fy=Hx2(∂2v∂x2+∂2v∂y2+∂2w∂y∂z)fy=Hx2(∂2v∂x2+∂2v∂y2+∂2w∂y∂z)
fz=Hx2(∂2w∂x2+∂2w∂y2+∂2v∂y∂z),fz=Hx2(∂2w∂x2+∂2w∂y2+∂2v∂y∂z), (12)
Therefore, considering w(x,y,z,t)=w(x,t)w(x,y,z,t)=w(x,t) , the vertical force affecting the lateral vibration per unit length of the Euler-Bernoulli nanotube is according to the following equation:
Flz=∫AfzdA=(ηA Hx2)∂2w∂x2Flz=∫AfzdA=(ηA Hx2)∂2w∂x2 (13)
Applied force from nonlinear viscoelastic foundationIn this study, vibration analysis of SWCNTs embedded in viscoelastic medium is presented. In this case, it is assumed that the chemical bonds of SWCNT to be formed generally between the external surface of the carbon nanotube and the viscoelastic medium (see Figures 1 and 2). In this study, the present viscoelastic foundation has linear Winkler stiffness, nonlinear stiffness, linear viscosity damping, and the attenuation is nonlinear. Therefore, in the present work, the force coming from the nonlinear viscoelastic medium could be considered in the general form as follows:
Fmedium=Flinear+Fnonlinear=Fmedium=Flinear+Fnonlinear=
(Flineark+Fnonlineark)+(Flinearc+Fnonlinearc)=(Flineark+Fnonlineark)+(Flinearc+Fnonlinearc)=
[(k1w(x,t)+(k3w3(x,t))]+(c0+c2w2(x,t))∂w(x,t)∂t),[(k1w(x,t)+(k3w3(x,t))]+(c0+c2w2(x,t))∂w(x,t)∂t), (14)
That includes spring element k1k1 , nonlinear spring k3k3 , linear damper c0c0 , and nonlinear damper c2c2 .
Piezoelectric equations considering surface effectsOne of the characteristics that distinguish nanomaterials is their surface effects. There are two essential and mechanically distinct surface effects in nanostructures: surface stress and surface elasticity. Atoms on the surface have different bonds than atoms in the balk, and surface stresses lead to these atoms being at the minimum energy level. Also, in nanostructures, surface elasticity occurs due to the lack of bonded neighbors in surface atoms.65,66
According to Euler-Bernoulli beam theory, the axial and transverse displacement fields can be considered as follows:
u(x,y,z,t)=U0−z∂w∂xu(x,y,z,t)=U0−z∂w∂x
w(x,y,z,t)=w(x,t),w(x,y,z,t)=w(x,t), (15)
Therefore, the only nonzero (axial) strain of this beam at any point is expressed by the following equation:
ϵxx=∂U0∂x−z∂2w(x,t)∂x2ϵxx=∂U0∂x−z∂2w(x,t)∂x2 (16)
Where ϵ0 is the axial strain applied to the structure and w (x, t) is also the transverse displacement (bending) at the desired location x on the center plate of the nanotube. Let's suppose Ez is the electric field created by the electric potential Φ assuming that this field exists only in the direction of the z-axis where Z indicates the distance from the center of the arrow, then the relationship between the electric current created and the existing electric potential is:67
Ez=−∂Φ∂z (17)
For the surface piezoelectric model, the surface-related structural relationships are different from the bulk volume equations. If the direction of polarization of the piezoelectric medium is assumed in the direction of the z-axis, the stress-strain structural equations related to the surface of the piezoelectric carbon nanotube-based on the relationships presented are written as follows:3–5
σsx=σ0x+cs11∈x−es31Ez (18)
Dsx=D0x (19)
In equations (18) and (19), σsx and Dsx are axial surface stress and electric surface displacement, respectively. The parameters σ0x , D0x , es31 and cs11 indicate the residual surface stress, the residual electrical displacement, the surface piezoelectric constant, and the elastic surface constant, respectively. On the other hand, in piezoelectric carbon nanotubes, the structural relationship between bulk stress and bulk strain would be obtained as follows:68
σx=C11∈x− e31 Ez
σx= C11∈x− e31 Ez (20)
Dz=e31ϵx+k33Ez (21)
That σx σx and Dz represent the classical stress tensor component and the electrical displacement, respectively. Also, C11 , e31 , and k33 are the elastic constant, piezoelectric constant, and dielectric constant of nanotubes, respectively.
In the absence of free electric charge ∂Dz∂z=0 , and the potential energy boundary conditions related to the voltage applied to the external environment of the nanotube is obtained from the following equation:
Φ (D2)=V, Φ (−D/2) = 0 (22)
By placing the relation (16) in equation (18), and then placing them in the constitutive equations (20) and (21) as well as using Gauss's law, the axial stress equations (in the local area) on the surface and volume of bulk nanotubes in the local space will be written as follows:
σx=C11ϵ0−e31VD−(C11+ e312k33)z ∂2w(x,t)∂x2 (23)
σsx=σ0x+cs11∈0− es31VD−(Cs11+es31 e31k33)z ∂2w(x,t)∂x2 (24)
Nonlocal governing differential equation of motion for equivalent continuum structure of embedded SWCNTThe equations of motion for the considered SWCNT, as shown in Figure 1, can be extracted as follows:
[(V(x,t))−(V(x,t)+∂V(x,t)∂xdx)+Fmediumdx+Fextdx=(ρAdx)∂2w(x,t)∂t2
−∂V(x,t)∂x+Fmedium+Fext=(ρA)∂2w(x,t)∂t2, (25)
And also for the momentum components:
Fvaeff∂w(x,t)∂x+V(x,t)−∂M(x,t)∂x=0 (26)
Where Fvaeff=Peff+Fmedium ; that Peff is the equivalent axial force which includes the force induced in the system due to the application of axial strain, surface electric charge and the electric field created in the whole nanotube structure, V (x, t) is the shear force at the cross section and M (x, t) is the bending moment at the cross section. Also, Fmedium is the force applied to the system by the viscoelastic foundation and FLz is the force applied in the nanotube by the magnetic field to the nanotube per unit length of the nanotube. m is also the mass per unit nanotube length. Based on volume and surface stresses of SWCN, the resultant bending moment can be written as follows for a beam model:
M=∫AσxzdA+∫SσxzdS=−EIeff∂2w(x,t)∂x2=
[πD464(E+e312k33)+πD38(cs11+es31e31k33)]∂2w(x,t)∂x2, (27)
where A,S are the circumference and surface area of the circular point of the nanotube, respectively, and EIeff is the equivalent bending stiffness of the whole piezoelectric nanotube. The equivalent axial force to the nanotube due to the application of strain and electric charge on the volume and surface of the nanotube can be considered as follows:3–5
Peff=πD24(E∈0+e31VD)+2(D(σ0x+cs11∈x)+es31V) (28)
So, the equation of motion in local theory is developed as follows:
∂2M(x,t)∂x2−Peff∂2w(x,t)∂x2−Fmedium−(ηA Hx2)∂2w∂x2=−(ρA)∂2w(x,t)∂t2 (29)
To apply the non-local effects, the equation (5) is used. First, equation (5) is rewritten as follows:
∂2M(x,t)∂x2=∂2∂x2(1−(τl)2∂2∂x2)Mnl(x,t) (30)
So:
∂2M(x,t)∂x2=−(τl)2∂4Mnl(x,t)∂x4+∂2Mnl(x,t)∂x2 (31)
Therefore, equation (29) can be rewritten as follows:
−(τl)2∂4Mnl(x,t)∂x4+∂2Mnl(x,t)∂x2−Peff∂2w(x,t)∂x2−{[(k1w(x,t)+(k3w3(x,t))]+(c0+c2w2(x,t))∂w(x,t)∂t)}−(ηA Hx2)∂2w∂x2=−(ρA)∂2w(x,t)∂t2, (32)
Finally, based on the equations (25) and (26), and for a flexural moment and cross-sectional shear force in non-local theory space, the differential equation of motion of the piezoelectric single-walled carbon nanotube considering non-local effects is developed as follows:
EIeff∂4w∂x4+[(k1w+k3w3)+(c0+c2w2)∂w∂t]−(ηA Hx2)∂2w∂x2−Peff∂2w∂x2+(ρA)∂2w∂t2−Fext−(τl)2
{[(k1∂2w∂x2+k3∂2∂x2w3)+(c0∂3w∂x2∂t+c2∂2∂x2(∂w∂tw2))]−
(ηA Hx2)∂4w∂x4−Peff∂4w∂x4+(ρA)∂4w∂x2∂t2−∂2∂x2Fext}=0,
(33)
The boundary conditions of nanotubes are also expressed according to the following relations, as shown in figures 3:
For both side simple support:
x=0,L; w=0,∂2w∂x2=0 (34)
For both side clamp support:
x=0,L; w=0,∂w∂x=0 (35)
Next, the equation (33) is investigated by nonlinear equation solving methods, and an analytical solution for nonlinear free and force vibration responses are extracted extracted.
Reducing the order of the equations of motion using Galerkin method
The Galerkin method is used to reduce the differential equations obtained. Therefore, the general form of vibrational response is considered as follows:
w(x,t)=∑i=Ni=1φi(x)qi(t) (36)
Where φi(x) is a mode shape function that satisfies the geometric boundary conditions, and qi(t) is the generalized time coordinates of the system.
By placing equation (36) in equation (33), and applying the Galerkin method, equation (37) is extracted:
[ˉmd2qidt2+β0dqidt+α1qi]+[β2qi2dqidt+α3qi3]=ˉfeiΩt (37)
Where:
α1=EIeffL∫0φiφi(4)dx+k1L∫0φi2dx−(Fverticaltotal)L∫0φiφi(2)dx−μk1L∫0φiφi(2)dx+μ(Fverticaltotal)L∫0φiφi(4)dx
α3=k3L∫0φi4dx−μk3(6L∫0φi'2φi2dx+3L∫0φi3φi(2)dx)
β0=c0L∫0φi2dx−μc0L∫0φiφi(2)dx
β2=c2L∫0φi3dx−μc2(6L∫0φi'2φi2dx+3L∫0φi3φi(2)dx)
ˉm=(ρA)L∫0φi2dx−μ(ρA)L∫0φiφi(2)dx , (38)
Solving nonlinear differential equation (the multiple time scales method)The multiple time scales method, as a perturbation method, is known as an efficient technique in nonlinear differential equations analysis. In this part, the multiple time scales method is used to develop the nonlinear vibrational response of the system.
To be able to apply this method, equation (37) is rewritten in the following form:
[d2qidt2+ε(β0+β2qi2)dqidt+ω2nqi]+[ε(α3)qi3]=0 (39)
In equation (39), ε ε is a small and dimensionless parameter. In the multiple time scales method, we first define the following parameters and equations:
Tn= ∈nt , ∂∂Tn=Dn , n=0,1,2,…
ddt=(D0+εD1+ε2D2+…)
d2dt2=(D02+2εD0D1+ε2D12+…) (40)
The general form of the response is also expressed as follows:
q(t)=q0(T0,T1,T2,…)+εq1(T0,T1,T2,…)+ε2q2(T0,T1,T2,…)+… = q0(T0,T1,…)+εq1(T0,T1,…)+Ο(ε2) (41)
Due to the attenuation in the equations of motion, the amplitude is expanded as the following equation:
a(ε)=a0+εa1++ε2a2+… (42)
Therefore, the nonlinear equation of the system is rewritten as follows:
[d2q1dt2+ε(β0+β2a2(ε)q12)dq1dt+ω2nq1]+ε[(α3a2(ε))q13]=0 (43)
Now, by placing equation (41) in the equation (43) and separating by different powers on both sides of the resultant equation, a set of differential equations will be obtained as follows:
Ο(ε0):[∂2∂T02+ω2n]q0(T0,T1,T2) (44)
Therefore:
q0=cos(ωnT0+γ(T1,T2)) (45)
Also:
Ο(ε1):[∂2∂T02+ω2n]q1(T0,T1,T2)=
[−β0∂∂T0−2∂2∂T0∂T1]q0(T0,T1,T2)−(a02q02(T0,T1,T2))[−α3q0(T0,T1,T2)+β2∂∂T0q0(T0,T1,T2)], (46)
By placing equation (46) in equation (43):
(D02+ω2n)q1=−14a02α3cos(3ωnT0+3γ(T1,T2))+14ωna02β2sin(3ωnT0+3γ(T1,T2)) +84(∂∂T1γ(T1,T2))ωncos(ωnT0+γ(T1,T2))−34a02α3cos(ωnT0+γ(T1,T2)) +(β0+14a02β2)ωnsin(ωnT0+γ(T1,T2)) (47)
The condition for the solvability of the high differential equation is to prevent the formation of very extensive terms in the time response. Therefore, the coefficients of the terms sin(ωnT0+γ(T1,T2)) and cos(ωnT0+γ(T1,T2)) in equation (47) should be considered equal to zero. This prevents the following terms from appearing in the system response, and does not take the vibrational response of the system to infinity over a long period time.
Therefore, the equations of the system are obtained as follows:
8(∂∂T1γ(T1,T2))ωn−3a02α3=0 (48)
4β0+a02β2=0 (49)
By solving the equations (48) and (49) together, the unknown variables are obtained in the following relation:
γ(T1,T2)=12α3a0232ωnT1+γ1(T2) (50)
So, the second-order estimation of the response will be obtained as follows:
q1=a02α332ω2ncos(3ωnT0+3γ(T1,T2))−a02β232ωnsin(3ωnT0+3γ(T1,T2)) (51)
Higher estimates can be obtained similarly.
To investigate forced vibrational response of the system, equation (37) is rewritten as follows:
d2qdt2+ε(β0+β2q2)dqdt+ω2nq+εα3q3=εˉfeiΩt (52)
In equation (52), the excitation frequency for the initial resonance is according to the relation Ω=εσ+ωn , where σ is the parameter of deviation from resonance.
So:
(D02+ω2n)q0=0 →q0=A(T1)eiT0ωn (53)
(D02+ω2n)q1=−β0D0q0−2D0D1q0−β2q02D0q0−α3q03+ˉfei(σT1+ωnT0) (54)
By placing equation (53) in equation (54), and to prevent significant time response terms, the coefficient of a term eiT0ωn must be set equal to zero, resulting is the following equation:
−2iωn∂A∂T1−[ˉA(α3+iβ2ωn)A+iβ0ωn]A+ˉfeiT1σ=0 (55)
Now, assuming the polar form A=a2eiγ , and placing it in the equation (55), the following set of equations will be obtained:
β0ωna+14β2ωna3−2ˉfsin(σT1−γ)=0 (56)
18α3a3−ˉfcos(σT1−γ)=0 (57)
The frequency equation response of the system with the definition of θ=σT1−γ , and considering that in the case of uniform response, the changes of a and θ will be zero. Finally, the nonlinear vibrational response equation is presented to the following final form in the form of a closed-form equation:
4ωn2[( 116β2a3+14β0a)2+(116α3a3ωn)2]−ˉf2=0 (58)
Dimensionless form of input and output parameters are considered by the following equations:
ˉx=xL , ˉw=wL , ˉfexc=fexcEI
¯EIeff=EIeffEI , ˉPeff=PeffL2EI , ˉFlz=FlzL2EI
K1=k1L4EI , K3=k3L6EI (59)
C0=c0L2√EIm , C2=c2L4√EIm
μ=(e0a)2L2 , T=tL2√EIm , Ω=ωL2√mEI
Validation for the results of the free vibration analysisIn this case, to validate the results with the references,69 the cross-sectional area of the nanowire is considered to be a rectangle with a small thickness and height h which its length is L. On the other hand, in this paper, using the classical local theory, the results of free vibration of a nanowire coated with piezoelectric crystals of PZT-5H for two boundary conditions, S-S and C-C, are obtained. Those properties are:
e31=−6.5;k33=1.3×10−8;
cs11=7.56;es31=−3×10−8;
E=126×109;L=20h;
In the following figures, the dimensionless frequency of the first mode of nanowire vibration in terms of its cross-sectional height for different values of an external voltage applied in local and non-local theory is compared with the results obtained in local formulation references.70 In the S-S boundary condition:
Figure 4 Validation for local70 and non-local (present) theory in S-S boundary conditions And also, for the C-C boundary condition:
Figure 5 Validation for local70 and non-local (present) theory in C-C boundary conditions
Furthermore, for the other forms of the surface effects:
Form 1:σ0x≠0 , cs11≠0,es31≠0
Form 2:σ0x=0 , cs11=0,es31=0
Form 3:σ0x≠0 , cs11=0,es31=0
Form 4:σ0x=0 , cs11≠0,es31=0
In the below curves, the changes of dimensionless frequency of the free oscillations at the first mode with two different boundary conditions in terms of the applied voltage and different values of dimensionless non-local parameter, with input parameters σ0x=1,Hx=0,C0=0,K1=10 , are shown.
Figure 8 The dimensionless frequency change curve of the free oscillations of the first mode of the system in terms of the applied voltage for different values of the dimensionless non-local parameter of the nanotube, in SS boundary conditions.
Figure 9 The dimensionless frequency change curve of the free oscillations of the first mode of the system in terms of the applied voltage for different values of the dimensionless non-local parameter of the nanotube, in CC boundary conditions.
In the below curves, the frequency changes of dimensionless free oscillations at the first mode of free vibration of the linear system for two different boundary conditions according to the magnitude of the applied magnetic field in different values of dimensionless non-local parameter, with input parameters σ0x=1,V=0.5,C0=0,K1=10 , are shown.
Figure 10 The dimensionless frequency change curve of the free oscillations of the first mode of the system in terms of the magnitude of the magnetic field for different values of the dimensionless non-local parameter of the nanotube, in SS boundary conditions.
Figure 11 The dimensionless frequency change curve of the free oscillations of the first mode of the system in terms of the magnitude of the magnetic field for different values of the dimensionless non-local parameter of the nanotube, in CC boundary conditions.
In the below curves, the frequency changes of the dimensionless free oscillations at the first mode of the free vibration of the linear system for two different boundary conditions according to the magnitude of the linear stiffness of the foundation in different forms considering surface effects, with input parameters σ0x=1,Hx=1×107,μ=0.1,C0=0,V=1 , are shown.
Figure 12 The dimensionless frequency change curve of free oscillations of the first mode of the system in terms of dimensionless linear stiffness of the substrate for different forms considering surface effects, in SS boundary conditions.
Figure 13 The dimensionless frequency change curve of free oscillations of the first mode of the system in terms of dimensionless linear stiffness of the substrate for different forms considering surface effects, in SS boundary conditions.
In the below curves, the frequency changes of the dimensionless free oscillations at the first mode of the free vibration of the linear system for two different boundary conditions according to the linear damping coefficient of the foundation in different forms considering surface effects, with input parameters σ0x=1,Hx=1×107,μ=0.1,K1=100,V=0.5 , are shown.
Figure 14 The dimensionless frequency change curve of real and imaginary parts of the special values of the first mode of the system in terms of dimensionless linear damping of the substrate for different forms considering surface effects, in SS boundary conditions.
Figure 15 The dimensionless frequency change curve of real and imaginary parts of the special values of the first mode of the system in terms of dimensionless linear damping of the substrate for different forms considering surface effects, in CC boundary conditions.
In the below curves, the results of the nonlinear vibration for the non-damped system are presented.
For Hx=0,μ=0.1,K3=100,V=0.5 :
Figure 16 Graph of Nonlinear Frequency to Linear Frequency Ratio of the Immortal System Base Mode for Different Values of the foundation Linear Hardness Parameter, Left: SS boundary condition, right: CC boundary condition.
For Hx=0,μ=0.1,K1=100,K3=100 :
Figure 17 Graph of Nonlinear to Linear Frequency Ratio Changes of the Immortal System Base Mode for Different Values of Applied Voltage, Left: SS Boundary Conditions, Right: CC Boundary Conditions.
For μ=0.1,K1=100,K3=100,V=0.5 :
Figure 18 Graph of Nonlinear Frequency to Linear Frequency Ratio of the Immortal System Base Mode for Different Values of Applied Magnetic Field, Left: SS Boundary Conditions, Right: CC Boundary Conditions.
For Hx=1×107,K1=100,K3=100,V=0.5
Figure 19 Graph of changes in the ratio of nonlinear to linear frequency ratio of the base mode of the immortal system for different values of the non-local parameter of the nanotube, left: SS boundary condition, right: CC boundary condition.
In the following two diagrams for zero voltage and negative voltage, nonlinear to linear frequency ratio diagrams are drawn for SS boundary conditions:
Figure 20 Graph of Nonlinear Frequency to Linear Frequency Ratio of Immortal System Base Mode for Different Forms Considering Surface Effects without Voltage Applied to SS Boundary Conditions.
Figure 21 Graph of changes in the ratio of nonlinear to linear frequency ratio of the base mode of the immortal system for different forms Considering surface effects with applied voltage V = -0.5, for SS boundary conditions.
For V=0,Hx=0,K1=0,K3=50,C0=0.15,C2=0,¯x0=0.5 :
Figure 22 Amplitude-frequency vibration curves for different values of μ (μ = 0 simple line, μ = 0.2 bold fold line); Figure above: SS boundary conditions, bottom figure: CC boundary conditions.
For V=0,Hx=0,K1=0,K3=50,C0=0.15,μ=0.1,¯x0=0.5 :
V=0, Hx=0, K1=0,K3=50, C0=0.15, μ=0.1, ˉx0=0.5:
Figure 23 Amplitude-frequency vibration curves for different values of C_2, (C_2 = 1 simple line, C_2 = 3 bold fold lines), Figure above: SS boundary conditions, bottom figure: CC boundary conditions.
For V=0,Hx=0,K1=0,K3=50,C2=1,μ=0.1,¯x0=0.5 :
V=0, Hx=0, K1=0,K3=50, C2=1, μ=0.1, ˉx0=0.5:
Figure 24 Amplitude-frequency vibration curves for different values of C_0, (C_0 = 0.1 simple line, C_0 = 0.3 bold fold line), high: SS boundary conditions, low: CC boundary conditions.
For Hx=0,K1=100,K3=50,C0=0.1,C2=1, μ=0.1,¯x0=0.5 :
Hx=0, K1=100,K3=50, C0=0.1, C2=1, μ=0.1, ˉx0=0.5:
Figure 25 Amplitude-frequency vibration curves for different values of V, (V = -0.5 purple line, V = 0 blue fold line, V = 0.5 bold black), high: SS boundary conditions, low: CC boundary conditions.
According to the presented diagrams, the presence of a non-local parameter has increased the curvature deviation to the right and the stiffening effect. In other words, the non-local parameter is a factor to increase the nonlinear effect of the system. On the other hand, the presence of this factor reduces the maximum range of fluctuations.
By comparing the graphs related to two different boundary conditions, the following results are inferred:
None.
None.
The authors declare that they have no competing interests.
©2023 Shahsavari, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.