Stability of Axially Moving Piezolaminated Viscoelastic Plate Subjected to Follower Force

The stability of the moving viscoelastic plate with the piezoelectric layer subjected to uniformly distributed tangential follower force is investigated. The force excited by the piezoelectric layer due to external voltage is modeled as the follower tensile force. The differential equation of the axially moving viscoelastic rectangular plate with piezoelectric layer subjected to uniformly distributed tangential follower force is formulated on the basis of the Kirchhoff thin plate theory and the two-dimensional viscoelastic differential constitutive relation. The complex eigenvalue equations are established by the differential quadrature method. Via numerical calculation, the curves of real parts and imaginary parts of dimensionless complex frequencies versus uniformly distributed tangential follower force and dimensionless moving speed are obtained. The effects of nonconservative force, dimensionless axially moving speed, and dimensionless applied voltages on the stability of axially moving nonconservative viscoelastic plate with piezoelectric layer are analyzed.


Introduction
Piezoelectric materials are used in the smart structure as sensors and/or actuators in the form of patches or layers.A plate with distributed, surface-bonded piezoelectric sheets has been widely used in active vibration, in structural health monitoring, in acoustic control, in buckling control, and among other applications.The behavior of laminated piezoelectric plates has gained much attention by researchers [1][2][3].For example, numerous scholars have studied the piezoelectric effects on the dynamic behavior of composite structures.The influence of using smart materials on the free vibration and natural frequencies of laminated piezoelectric plates has also been addressed [4][5][6][7][8].Structural systems lose their stability due to divergence, and such systems under partial follower compressive loading were investigated both qualitatively and quantitatively [9].Chandrashekhara and Bhatia [10] investigated the active buckling control of smart composite plates based on finite element analysis.In [11], both cross force and voltage, together with the bending of the plate, were analyzed, and the bending of the plate was controlled through the application of different voltages to the actuators.In accordance with the inverse piezoelectric effect, the equivalent action of the actuators was obtained.Using Hamilton's principle, the finite element formula of the bending deformation of the piezoelectric structure was derived.The bending of the annular plate made of piezoelectric material with four boundary conditions was also calculated.Chase and Bhashayam [12] used the piezoelectric sheet to increase the critical load of laminated beams and performed optimal stabilization on plate buckling.Wang [13] and Wang and Quek [14] studied the buckling of column structures with a pair of piezoelectric layers and improve the flutter and buckling capacity of the column using the piezoelectric layers.Ha et al. [15] examined a composite structure containing distributed piezoceramic sensors and actuators using finite element analysis.
Stability and transverse vibrations of axially moving systems are present in various industrial applications, such as the paper webs and plastic sheets, steel strip in a thin 2 Advances in Materials Science and Engineering steel sheet production line, band saw blade, conveyor belts and chain in power transmission lines, and aerial cableways.The axially moving considerations have been studied widely, and dynamics and stability have been reviewed by Wickert and Mote [16].The research on two-dimensional axially moving plate can be dated back to Ulsoy and Mote Jr. [17], who analyzed the coupled transverse and torsional vibration of band saw blade.The stability and vibration characteristics of an axially moving plate with two simple supported edges and two free edges under homogeneous tension have been investigated by Lin [18,19].A few researches on transverse vibrations and stability of axially moving viscoelastic plates have also been done.Zhou et al. [20,21] studied transverse vibration characteristics of axially moving viscoelastic rectangular plates and parabolically varying thickness plate.Marynowski [22] compared a threeparameter Zener model and a two-parameter Kelvin-Voigt model for the viscoelasticity.Saksa et al. [23] investigated stability and dynamic behaviour of axially moving viscoelastic panels with the help of the classical modal analysis.Tang and Chen [24] studied stability in parametric resonance of moving viscoelastic plates with time-dependent travelling speed.
It is well known that instability may be induced in structures systems because of nonconservative forces.Indeed, much of the research in this field has focused on the stability of the structures subjected to different types of nonconservative forces.However, few papers have focused on the stability of axially moving viscoelastic plates with piezoelectric layers under nonconservative forces.
Here, a relatively novel solution technique, the differential quadrature method (DQM), is used to analyze dynamic stability behavior of axially moving viscoelastic plates with piezoelectric layers subjected to nonconservative force.The current work focuses on the stability improvement of the axially moving nonconservative viscoelastic plate that is surface-bonded by piezoelectric layers.An analytical model that fully embodies the piezoelectric effects is first obtained based on the behavior of the follower force imposed by the external voltage on the piezoelectric layer.This paper aims to present the differential equation of the axially moving viscoelastic plate with piezoelectric layers subjected to uniformly distributed tangential follower force.Based on the 2D viscoelastic differential constitutive relation and Kirchhoff thin plate theory, the differential equation of the axially moving viscoelastic rectangular plate with piezoelectric layer subjected to uniformly distributed tangential follower force in the Laplace domain is formulated.And carrying out the Laplace inverse transformation, the differential equation of motion of the axially moving viscoelastic rectangular plate with piezoelectric layer subjected to uniformly distributed tangential follower force in time domain is yielded.The complex eigenvalue equations are established by the differential quadrature method.Via numerical calculation, the complex frequency and the instability type of the axially moving nonconservative viscoelastic plate constituted by elastic behavior in dilatation and the Kelvin-Voigt laws in distortion are obtained.The effects of nonconservative force, dimensionless axially moving speed, and dimensionless applied electric potential difference on the stability of axially moving nonconservative viscoelastic plate with piezoelectric layer are analyzed.

Differential Equation of Motion.
The axially moving viscoelastic thin plate with surface-bonded piezoelectric layers subjected to uniformly distributed tangential follower force  0 is shown in Figure 1.The piezoelectric layers are symmetrical to the midplane of the structure.The thin plate has a thickness ℎ  in the -direction, a material density   , and Poisson's ratio   .The piezoelectric layers have a thickness ℎ  , a layer density   , and Poisson's ratio   .
Using Kirchhoff 's hypothesis of classical thin plates, the strain and displacement relations for the viscoelastic plate can be written as where  * is the transverse displacement of the plate or deflection and  is the distance of the arbitrary point of the plate from the neutral plane.The strain of the piezoelectric bonding layers pertains to the strains of both surfaces of the viscoelastic plate, expressed as Considering mechanical and electrical behavior of piezoelectric materials.Proper care is taken when constitutive equations for piezoelectric materials are used to model induced strain actuation, as they are poling direction dependent.The planar isotropy of poled ceramics is expressed by their piezoelectric strain constants, such that  31 =  32 Advances in Materials Science and Engineering 3 [26,27].The constitutive equations of the piezoelectric layers for plane stress are where  denotes applied electric potential difference.Substituting (2) into (3) yields The constitutive relations of linear viscoelastic materials in the Laplace domain are In the Laplace domain, the constitutive relations of linear viscoelastic materials using Laplace transformation  of deflection  * are expressed as where the differential operators , and the polynomials  0 ,  0 , and  1 about Laplace variable  are independent of spatial coordinates.
Assume that the piezoelectric sheets are of infinitesimal thickness and disregard its influence on the bending stiffness of the structure.The following terms    and    result in uniform stress in the piezoelectric layer.The force excited by the piezoelectric layer due to external electric potential difference is modeled as a follower tensile force, and the follower tensile force is a kind of nonconservative force.The tension force caused by the piezoelectric layer can be written as In the Laplace domain, the relations between internal torque and Laplace transformation w of deflection  * are given by The transverse speed of moving plate is The differential operator is given by / = (/) + V(/).
The equilibrium equation of axially moving viscoelastic plate with piezoelectric layer subjected to uniformly distributed tangential follower force is given by where Carrying out the Laplace transformation of ( 10) and ( 11) and multiplying it by  0   .When the partial derivative is continuous, in the Laplace domain the differential equation of motion of the axially moving viscoelastic rectangular plate with piezoelectric layers subjected to uniformly distributed tangential follower force is obtained: where Equation ( 11) is suitable for various viscoelastic differential models.
Here, we assume that the material of the plate obeys elastic behavior in dilatation and the Kelvin-Voigt laws in distortion.Substituting the polynomials   = 1,   = 2 + 2,   = 1, and   = 3 into (11) and carrying out the Laplace inverse transformation yield the differential equation of motion of axially moving nonconservative viscoelastic plate with elastic dilatation and Kelvin-Voigt distortion with piezoelectric layers in the time domain where Substituting ( 15) into (13) yields the following dimensionless differential equation: where  is the dimensionless time,  denotes the dimensionless time delay,  is the dimensionless moving speed, and  represents the dimensionless electric potential difference parameter.One has ∇ 4  = ( 4 / 4 ) + 2 2 ( 4 / 2  2 ) +  4 ( 2 / 4 ).Suppose that the solution to ( 16) takes the form where  = √ −1 and  is a dimensionless complex frequency.Substituting ( 17) into (16) gives where Boundary conditions for four-edge simply supported (SSSS) rectangular plates are Boundary conditions for two opposite edges simply supported and other edges clamped (CSCS) are as follows:

Complex Eigenvalue Equation.
The complex eigenvalue equation of the axially moving viscoelastic rectangular plate with piezoelectric layers subjected to uniformly distributed tangential follower force is derived by the DQM [28][29][30][31].
The DQM is used to approximate the partial derivatives of a function with respect to a spatial variable at any discrete point as the weighted linear sum of the function values at all the discrete points chosen in the solution domain of spatial variable.
Postulating smooth function (, ) in the region 0 ≤  ≤ , 0 ≤  ≤ , the th order partial derivative of (, ) with respect to , the th order partial derivative of (, ) with respect to , and the mixed partial derivative of the th order with respect to  and the th order with respect to  are, respectively, approximated as where  and  are the number of grid points in the  and  direction, respectively,  ()  and  ()  are the weighted coefficients, and they are defined by In the case of  = 2, 3, . . .,  − 1,  = 2, 3, . . .,  − 1,  ( In this paper,  = .The distribution forms of the grid points are nonuniform; the distribution forms of the grid points are  According to the DQM procedures, ( 18) can be discretized into the following forms:
After eliminating the boundary degrees of freedom from (22) using the boundary conditions, the equation can be written in the matrix form as The complex eigenvalue equation of the axially moving viscoelastic rectangular plate with piezoelectric layers subjected to uniformly distributed tangential follower force is that coefficient determinant equal to zero; that is, In (30),  is a generalized complex eigenvalue, where the matrices Q, R, G, and K involve such parameters as dimensionless moving speed , the dimensionless follower force , the dimensionless time delay , the dimensionless time , the dimensionless electric potential difference parameter , the dimensionless electric potential difference geometrical parameter , and aspect ratio  of the plate.

Results and Discussion
In the case of  = 0,  = 0,  = 0, and  = 0, (18) is reduced to the differential equation of motion of the elastic plate.Equation ( 30) is reduced to the complex eigenvalue equation of the elastic plate.In order to verify the DQM program, the first three-order natural frequencies of the transverse free vibration of elastic plate with four-edge simply supported boundary conditions are calculated firstly, and the results are in good agreement with those exhibited in [31], which can be seen from Table 1.
In the numerical examples, the number of grid points is  =  = 9.Numerical studies have been conducted to investigate the effects of several key parameters on the dynamics stabilities of axially moving viscoelastic plates with the piezoelectric layer subjected to uniformly distributed tangential follower force.The critical moving speed of axially moving plates for  = 10 −5 ,  = 0,  = 1,  = 0, and  = 0 with different boundary conditions is seen from Table 2.

The Influence of the Dimensionless Electric Potential
Difference. Figure 2 displays the variation of the curve for the cases  = 10 −5 ,  = 50,  = 1,  = 0, and  = 0.05.It can be seen that, with increase of axially moving speed , the real part Re() of complex frequencies in the first mode becomes zero and the imaginary part Im() of complex frequencies has two branches.This shows that the first mode behaves unstably by the divergence instability when the axially moving speed becomes equal to or larger than the critical speed  = 3.13.When the moving speed increases further, the piezolaminated plate regains stability in the first-order mode.By maintaining an increase in the moving velocity, the plate undergoes divergence instability in the second mode.When the speed of motion is  = 8.78, the first mode behaves unstably by the divergence instability again.When  = 10.85, the second mode behaves divergence instability again.
In Figure 3, the case of  = 10 −5 ,  = 50,  = 1,  = 10, and  = 0.05 is shown.The first-order mode exhibits divergence instability when the moving speed  = 5.77; at  = 7.41 the first-order divergence instability is over.When the moving speed increases, the plate undergoes divergence instability in the second mode.By maintaining an increase in the moving velocity, the first mode behaves unstably by the divergence instability again.Then, the second mode behaves divergence instability again.In comparison with Figure 2, with the dimensionless electric potential difference increase, the critical speed increases, while the types of instability of the moving piezolaminated plate are in agreement.
Figure 4 shows the variation of the first three-order dimensionless complex frequencies of the plate with dimensionless axially moving speed  for  = 10 −5 ,  = 50,  = 1,  = 20, and  = 0.05.It can be seen that the influence of the dimensionless electric potential difference on the critical speed is obvious.With the dimensionless electric potential difference increase, the critical speed increases, while the types of instability of the moving piezolaminated plate subjected to uniformly distributed tangential follower force have not been changed.

The Influence of the Dimensionless Moving Speed on
Stability.Figures 5 and 6 display the variation of the curve of the first-order dimensionless complex frequencies versus dimensionless follower force  for the cases  = 10 and  = 20.It can be seen that the viscoelastic plate exhibits divergence instability in the first-order mode under different moving speed .This result indicates that increases in the dimensionless speed of moving piezolaminated viscoelastic plate subjected to uniformly distributed tangential follower force reduce the critical load but will not change the type of instability experienced by the moving plate.

The Influence of the Dimensionless Uniformly Distributed
Tangential Follower Force on Stability.Figures 7 and 8 indicate the variation of the first-order dimensionless complex frequencies with the dimensionless moving speed for the different values of , respectively.In Figure 7, for dimensionless follower force  = 0, it shows that the first-order mode behaves divergence instability; then the plate undergoes single-mode flutter.In the case of  = 50 and  = 100, it shows that the first-order mode behaves divergence instability first, when the dimensionless moving speed increases, the plate regains stability in the first-order mode.By maintaining an increase in the moving speed, the first mode behaves divergence instability again.It can be seen that the real part of the dimensionless complex frequencies decreases with the increase of the dimensionless follower force .The critical speed of the first mode decreases when the dimensionless follower force  increases.
In comparison with Figures 7 and 8, respectively, it is observed that the same dimensionless follower force , the real parts of the dimensionless complex frequencies under  = 20, is greater than that in the case of  = 10.The critical speed of the first mode increases with the increase of electric potential difference.

Conclusions
The formulation and results for the stability of the moving viscoelastic plate with the piezoelectric layer subjected to uniformly distributed tangential follower force are presented.The general complex frequencies are calculated.This paper analyzes stabilities of the moving viscoelastic plate with the piezoelectric layer subjected to follower force.The results of the analysis of the present study are summarized below.
(1) With the dimensionless electric potential difference increase, the critical speed increases, while the types  Imag() q = 0 q = 50 q = 100  of instability of the moving piezolaminated plate subjected to uniformly distributed tangential follower force have not been changed.The moving piezolaminated plate subjected to follower force behaves unstably by the divergence instability.(2) By increasing in the dimensionless speed of moving piezolaminated viscoelastic plate subjected to follower force, the critical load (i.e., critical follower force) decreases but will not change the type of divergence instability experienced by the moving plate.
(3) The influence of the dimensionless uniformly distributed tangential follower force on stability of moving piezolaminated viscoelastic plate is obvious.By maintaining an increase in the dimensionless follower force, the critical speed of the first mode decreases.
This research demonstrates the potential of using piezoelectric materials in enhancing the stability of axially moving viscoelastic plates subjected to uniformly distributed tangential follower force.The conclusions provide a theoretical basis and effective approach for improving structure designing and working stability.

Figure 1 :
Figure 1: Configuration of axially moving viscoelastic plate with piezoelectric layers subjected to uniformly distributed tangential follower force.

Table 1 :
The first three-order natural frequencies of the uniform thickness elastic plate with different boundary conditions.

Table 2 :
Critical moving speed of axially moving plates with different boundary conditions.