Vibration Characteristics for Moving Printing Membrane with Variable Density along the Lateral Direction

The vibration model of moving membrane with variable density distribution is established, and the density distribution of the moving membrane varies along the lateral direction. The transverse vibration differential equations of moving membrane are established based on D’Alembert’s principle and discretized by using the differential quadrature method (DQM).The relationships of the first three dimensionless complex frequencies between dimensionless speed, density coefficient, and tension ratio of the membrane are analyzed by numerical calculation.The effects of the density coefficient and the tension ratio on transverse vibration characteristics of the membrane are investigated. The relationship between density coefficient and critical speed is obtained. The numerical results show that the density coefficient and the tension ratio have important influence on the stability of moving membrane. So the study provides a theoretical basis for improving the working stability of themembrane in the high-speed printing process.


Introduction
The membrane including plastic film, paper web, cloth, metal foil, and other types of film-like material is widely used to make printing and packaging products.The processing schematic of axially moving membrane is shown in Figure 1.The membrane deformation, the membrane folding, the surface scratches, and other processing defects can be caused by the transverse vibrations in printing process, which will affect the printing accuracy and printing quality seriously [1].Therefore, the systematic study of transverse vibration characteristics of the membrane is important, which has caught the attention of many scholars.
Kulachenko and his coworkers [2,3] studied nonlinear dynamics problem; the stability of transverse vibration of the web was studied by finite element method.The elementfree Galerkin method was used to analyze free vibration of thin plates resting on Pasternak elastic foundations with all possible types of classical boundary conditions by Bahmyari et al. [4].Nguyen et al. [5,6] analyzed the stability and the control of transverse vibration of web by changing axial velocity and axial tension, respectively.Banichuk and his coinvestigators [7] considered the dynamic characteristics of the moving web with nonhomogeneous tension, and the analytical approaches method was used.Vedrines et al. [8] measured the vibrations of the printing web with a laser sensor and determined the vibration in plane.Wang et al. [9] applied finite difference method to acquire natural vibration frequency of the annular membrane with wrinkle.The results showed that the frequency increased with the increase of wrinkling level.Wu et al. [10,11] established the vibration equations and analyzed the stability of the membrane with an intermediate elastic support and variable speed.
In these literatures above, the membrane chosen as sample was all supposed with constant density.But the membrane density is changing in many situations.In printing process, the plate is wetted and inked based on the image distribution, and then the ink is transferred to the substrate through rolling of the plate cylinder and impression cylinder.The ink and fountain solution which are absorbed into the paper web or membrane cause the change of density.Besides, the nonuniform thickness of membrane also can change the surface density [12].Nonuniform density will affect the vibration characteristics and stability of membrane.Therefore, the study of the transverse vibration characteristics of the moving membrane with variable density has a great significance on controlling the membrane vibration during the printing process and improving overprint precision of printed products.
Transverse natural vibration of the annular membrane with variable surface density and without axial velocity was studied by using modified perturbation method in [13].A numerical analysis of axisymmetric transverse-free vibration of the circular membrane with nonhomogeneous changing density was analyzed by Li [14,15], and the first-order natural frequencies were given.Buchanan [16] examined the circular membrane dynamic properties, and the density distribution of the circular membrane was a linear variation along the diameter direction.Gupta and Khanna [17] and Zhou and Wang [18] studied the vibration characteristics of the viscoelastic plate with parabolic variable thickness.The suboptimal control method was applied to control the variable density of the moving web by Ma and her coworkers [19].Sheet flutter and the interaction between sheet and air were analyzed by Pramila [20]; the results showed that surrounding air had an impact on the critical velocities and the eigenfrequencies.He [21] also studied the natural frequencies of an axially moving band vibrating in an ideal fluid.And the influence of the surrounding fluid was taken into account by using hydrodynamic added mass.Niemi and Pramila [22] researched transverse vibrations of an axially moving membrane submerged in ideal fluid by using the FEM.The effect of the density of the element mesh, the truncation distance, and the various lumping techniques on the accuracy of the results was also analyzed.Gutierrez et al. [23] employed a series of numerical experiments to deal with the transverse vibration of annular membranes where the density varies with the radial variable linearly, quadratically, and cubically.A general quasi-analytical model based on the Frobenius power series expansion method was described so as to handle vibrations of solid circular and annular membranes with continuously varying density by Willatzen [24].Bala Subrahmanyam and Sujith [25] also studied the traverse vibration of annular membranes with continuously varying densities.In this study, the stability of the moving membrane with variable density in the lateral direction is studied.The transverse vibration differential equation of the moving membrane with variable density is established and discretized by using DQM.The relationships between the dimensionless vibration frequency and the dimensionless velocity, the density coefficient, and the tension ratio of the printing membrane are studied.The influence of the density coefficient and the tension ratio on vibration characteristics of the membrane is discussed.

The Transverse Vibration Model of Moving Membrane
The membrane between two guide rollers can be simplified to a moving membrane model shown in Figure 2. The rollers and the brush mechanisms (in Figure 1) which support the membrane can be regarded as the boundary condition with the four edges fixed.Here, the membrane is soft, homogeneous, and inextensible.The thickness of the membrane is neglected.The translating direction of membrane is  direction.The membrane lateral or width direction is  direction.The displacement direction of transverse vibration is the  direction.The membrane length is , the width is , and V is axially translating velocity.The function (, , ) denotes the transverse displacement of the printing membrane, where  is time.  and   are uniform tension per unit length in the ,  direction.The density function of membrane is (), as is shown in Figure 3.The membrane surface density is varied along  direction that can be expressed as follows: where  denotes the density coefficient.
The speed V() in the transverse displacement direction can be given as follows: Then the transverse acceleration  can be expressed as follows: ( As shown in Figure 4, consider an elemental area of the membrane  and analyze the force on them.The composition force in  direction caused by the membrane tension on  and  can be derived as follows: It is assumed that the membrane is subjected to a transverse load (, , ) in  direction; here,  represents time.The motion equation of vibration system can be given based on D' Alembert's principle [26]: Irrespective of the transverse load, that is, (, , ) = 0, then the transverse vibration differential equation of membrane can be expressed as follows: The piecewise functions () can be expressed as follows: Substituting ( 4) into (3) yields Introduce the dimensionless quantities as follows: Substituting ( 9) into (8) yields Supposes the solution to (10) is Substituting ( 11) into (10) obtains the transverse vibration differential equations of moving membrane with variable density: where  1 = √ −1,  is dimensionless time,  is dimensionless natural frequency,  is dimensionless speed,  is tension ratio, and  is aspect ratio.
Figure 5 shows the classical boundary conditions for an axially moving membrane, which are four sides fixed and three sides fixed and one free and two sides fixed and two free, and we only consider the boundary condition with four sides fixed.
The boundary conditions of moving membrane with four sides fixed are expressed as follows:

Establishment of Complex Characteristic Equation
The principle of DQM [27] is using the weighted sum of the function value of all nodes on the whole domain to replace the value of the function and its derivative at a given node, so the differential equations are turned into algebraic equations which regard the function value at the node as unknown.
Introduce  ×  grid points according to DQM.They are The value of each of the order partial derivatives of the unknown function can be described as [2]     −  2  ∑ =1  [2] [2] [2] The weight coefficient of each order can be acquired by ) The complex characteristic equation is established by using DQM.Then ( 15) is turned into differential quadrature form: [2]     −  2  ∑ =1  [2] [2] [2] [2] [1] The boundary conditions are Converting ( 18) and ( 19) to matrix form obtains Matrix R, G, K contains tension ratio , dimensionless velocity , density coefficient , aspect ratio , and other parameters.

Numerical Results and Analysis
When the density coefficient  = 0, the vibration equation of moving membrane degenerates the transverse vibration equation in which density distribution is uniform.To illustrate the effectiveness of the differential quadrature method, assuming density coefficient  = 0, dimensionless speed  = 0, tension ratio  = 1, and aspect ratio  = 1,  = 2, respectively, the results comparison between this study and [28] is shown in Table 1.Then assuming density coefficient  = 0, dimensionless velocity  = 0, aspect ratio  = 2, and tension ratio  = 0.5,  = 0.8, respectively, the results  comparison between this study and [28] is shown in Table 2.It can be seen from Tables 1 and 2 that the degenerated results in this study are in good consistency with the results in [28].Therefore, it proves that the method used in this paper is feasible and believable.
Figures 6 and 7 show the relationship between first three dimensionless natural frequencies and the dimensionless speed, when the density coefficient  = 0, the aspect ratio  = 2, and the tension ratio  = 1 and  = 0.5, respectively.As illustrated in Figure 6, the dimensionless speed increases, the real part Re() of the first three order dimensionless complex frequencies tends to decrease gradually.When 0 <  < 1, the first three dimensionless complex frequencies diminish with the increase of dimensionless speed, and the imaginary part of the complex frequencies is zero consistently, so the membrane works in steady state.When the dimensionless speed  = 1, the real part of the complex frequencies becomes zero simultaneously; here, the dimensionless speed  = 1 is called the critical speed   .When  > 1, the imaginary part of the complex frequencies is not zero anymore, and the membrane begins to be in a divergent instability state.Compared with Figure 7, for any same dimensionless speed in stable work region, the dimensionless complex frequency decreases gradually with the decrease of tension ratio.
Figure 8 shows the relationship between density coefficient  and critical speed   , when  = 1,  = 2 and  = 0.5, and  = 2, respectively.It is found from the figure that the critical speed diminishes with the increase of density coefficient.The critical speed is   = 1 for the different tension ratio when  = 0.When  > 0, the density coefficient is fixed, and the smaller the tension ratio, the smaller the critical speed.
Figure 9 shows the relationship between the first three dimensionless natural frequencies and the density coefficient when  = 1,  = 2 and  = 0,  = 0.5, respectively.The results show that the real part of the dimensionless complex frequency is decreasing with the increasing density coefficient and the real part of the dimensionless complex frequency is decreasing overall when the dimensionless velocity increases.
Figures 10(a)-10(d) show the relationship between the first three dimensionless natural frequencies and the tension ratio when ( = 2,  = 0,  = 0), ( = 2,  = 0.5,  = 0), ( = 2,  = 0,  = 0.8), and ( = 2,  = 0.5,  = 0.8), respectively.Comparative analysis shows that, for any order dimensionless complex frequency when dimensionless speed  = 0, the vibration frequency of the moving printing membrane increases with the increase of the tension ratio; the vibration frequency is downward with increasing dimensionless velocity; the real part of dimensionless complex frequency decreases with the increasing of density coefficient when the dimensionless velocity is the same.In this process, the imaginary part of complex frequency Im() = 0, so the motion system is in a stable state.

Conclusions
The DQM is used to analyze the transverse vibration characteristics of moving membrane with variable density distribution along the lateral direction.The results show that the tension ratio and the density coefficient have important impacts on the stability of moving membrane.The conclusions are as follows: (1) When  = 0, for the same aspect ratio, with the increase of tension ratio, the motion system is in a stable state in region 0 <  < 1.When the actual speed reaches the critical speed (  = 1), the motion system comes into a divergent instability state.The tension ratio of the membrane has no influence on the critical speed.
(2) When density coefficient  > 0, for the same aspect ratio, with the decrease of tension ratio, the critical speed is becoming smaller, the stable work region of motion system is smaller too.For any same dimensionless speed of stable work region, with the decrease of tension ratio, the dimensionless complex frequencies tend to be smaller and smaller.The tension ratio has important influence on the stability of moving membrane.
(3) The density coefficient  has an influence on the critical speed and the dimensionless natural frequencies.The critical speed and the natural frequencies diminish with the increasing density coefficient.Tension ratio does not affect the change trends of critical speed with density coefficient, but it affects the critical speed value.
The study provides a theoretical basis for improving the working stability of the membrane in the high-speed printing process.51305341) and Natural Science Foundation of Shaanxi (Grant no.2016JM5023).

Figure 2 :Figure 3 :
Figure 2: The model of printing moving membrane under the four edges fixed.

Figure 4 :
Figure 4: The force on an elemental area of the membrane.

Figure 5 :
Figure 5: The boundary conditions of the printing moving membrane.