# Machine Directional Register System Modeling for Shaft-Less Drive Gravure Printing Machines

Authors: Shanhui Liu; Xuesong Mei; Jian Li; Li'e Ma
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

In the latest type of gravure printing machines referred to as the shaft-less drive system, each gravure printing roller is driven by an individual servo motor, and all motors are electrically synchronized. The register error is regulated by a speed difference between the adjacent printing rollers. In order to improve the control accuracy of register system, an accurate mathematical model of the register system should be investigated for the latest machines. Therefore, the mathematical model of the machine directional register (MDR) system is studied for the multicolor gravure printing machines in this paper. According to the definition of the MDR error, the model is derived, and then it is validated by the numerical simulation and experiments carried out in the experimental setup of the four-color gravure printing machines. The results show that the established MDR system model is accurate and reliable.


---

## Body

## 1. Introduction

Recently, with the development of high-speed and high-precision printing, shaft-less drive technology has been widely used in the latest gravure printing machines. In this technology, the printing rollers are driven by independent servo motors, and the register error is controlled through directly adjusting the speed of the pertinent servo motor instead of changing the position of the compensation roller, so the traditional mechanical shaft is eliminated to make printing faster and more accurate. This register mode has the advantages such as higher response frequency, better control accuracy, and lower additional cost. However, the register error in the shaft-less drive mode has a cumulative effect; namely, the adjustment of register error in upstream will affect all register errors in downstream. Hence, the mathematical model of the register system which can accurately reflect the variation of the register errors is needed to be studied to control the register errors accurately.It is difficult to get an exact register system model because the register errors are affected by various factors, such as the web tension and speed, the precision of parts and assembling, and the web material and production environment. Register errors are defined as two-dimensional errors: machine directional register (MDR) errors and cross-directional register (CDR) errors. Brandenburg [1] derived a linear mathematical model of the MDR errors of a moving web in a first-order differential equation by using an equilibrium equation of mass which is transported by printing rollers, but this model was based on the mechanical shaft drive mode and could not adapt to the shaft-less drive mode. Yoshida et al. [2–4] developed the MDR system model for the mechanical shaft and shaft-less drive mode, but the model still ignores many influencing factors and is imperfect. Hua et al. [5] got the MDR errors model of gravure printing machines via the system identification method, but this method has poor universality and bad practicability. Although the MDR errors model of the two-color system was established in [6], the MDR system model of the multicolor system was not derived, some key deducing steps were ignored, and some wrong symbols appeared in some formulae, which makes some processes of derivation difficult to be understood for readers. The modeling and control of CDR system were studied in [7–9]. The CDR errors are controlled by the lateral motion of printing rollers, which is greatly different from the control of MDR errors. In addition, initial works of modeling and control in correlation between MDR and CDR errors were researched in [10, 11]. However, the in-depth research on the MDR system model based on shaft-less drive mode is still not enough in gravure printing machines.Therefore, in this study, a mathematical model of the MDR system in shaft-less drive mode is derived for gravure printing machines by using an equilibrium equation of mass and considering the fluctuation of tension, speed, and other factors. This paper is structured as follows. Section2 introduces the schematic diagram of shaft-less drive gravure printing machines. Section 3 establishes the MDR errors model of the two-color system and then obtains the MDR system model of the multicolor gravure printing machines. Section 4 verifies the derived model by simulation and experiments. The results show that the model is accurate enough to describe the MDR errors in shaft-less drive gravure printing machines.


## 2. Shaft-Less Drive Gravure Printing Machines

The schematic diagram of then-color gravure printing machines is presented in Figure 1, and n stands for the total number of printing units. This system is composed of an unwinding unit, an infeeding unit, an n-color printing system, an outfeeding unit, and a rewinding unit. Shaft-less drive mode has been used in the machines; that is, all of the driving shafts are driven by independent servo motors. The n-color printing system consists primarily of n printing rollers, and c1 is the first printing roller whose main task is to maintain the system speed constant. In the printing subsystem, tension is same with that in infeeding zone at steady state. Load cells are installed at idle rollers in the middle of continuous process in both the infeeding and outfeeding units for tension pickup. Two passive dancer rollers are used to reduce tension fluctuations and measure the tension signals simultaneously in the unwinding and rewinding units. Photoelectric sensors are installed in the machine for registering errors pickup.Figure 1
Schematic diagram of then-color gravure printing machines.Control of the gravure printing machines is basically composed of tension control and register control. Register control is required for high printing resolution in the multistage printing systems. Hence, an accurate MDR system model is essential to ensure adequate printing performance for the gravure printing machines.


## 3. Mathematical Modeling

### 3.1. Basic Assumption

In the gravure printing machines, the MDR error of a moving web is defined in two adjacent printing rollers as the relative difference of the distance between the previous printed pattern by the upstream printing roller and the later printed one in the downstream printing roller. To remove unnecessary interference and simplify the mathematical model, the following assumptions are made in the derivation of model.(1)
There is no slippage between web material and rollers and no mass transfer between the web material and the environment.(2)
The web material is isotropic, and its properties do not change with temperature and humidity.(3)
The web cross-section in the unstretched state does not vary along machine direction in a web span.(4)
The density and modulus of elasticity of the web in the unstretched state are constant in the cross-section.(5)
The strain in web is small (much less than unity) and uniform within the web span.(6)
The printing system is in a stable state, in which there are no disturbances from tension and other factors in timet≤0.


### 3.2. MDR Error Model of the Two-Color System

The mathematical MDR error model of the two-color system in shaft-less drive gravure printing machines is illustrated in Figure2 by the coupling between the ith and (i+1)th units, and the ith unit is register errors reference. As shown in Figure 2, a coordinate axis along the machine direction is set up and the position of the ith unit is set as the origin of the coordinate axis. Figures 2(b) and 2(c) show that the color marks on printing roller ci and ci+1 are printed on the web at the printing instants t1 and t2, respectively.Figure 2
MDR error model of the two-color system.The MDR error in the moving coordinate system is to calculate the difference between the printing point on the roller and the printed point on the web representing the same point of the pattern at the printing instants. Therefore, according to the moving coordinate system, the MDR error is just the deviation from the coordinate of color mark on the printing rollerXFi+1 and that on the web XBi. At instant t2, the MDR error is defined by
(1)ei(t2)=XBi(t2)-XFi+1(t2).
It is to be noted that the movement of the printing roller can be viewed as a combination of two movements: one is the rotation of the roller and the other is the virtual translational movement of the roller along the machine direction which is caused by the variation of the web length. This variation can be affected by many factors such as dryer system, tension sensors, or some other disturbances.By using assumptions (1) and (2), the following conclusion can be obtained:
(2)ddt∫0Liρi(x,t)Ai(x,t)dx=ρi-1(x,t)Ai-1(x,t)ViR(t)-ρi(x,t)Ai(x,t)Vi+1R(t).
Hence, the mass of the web material entering the printing unit i in the time interval from t1 tot2 can be derived as follows:
(3)ddt∫XBi(t1)XBi(t2)ρi(x,t)Ai(x,t)dx=ρi-1(x,t)Ai-1(x,t)ViR(t).
According to assumptions (3) and (4) and the definition of strain, the relationship between the stretched and un-stretched webs can be written as follows:
(4)ρi(x,t)Ai(x,t)ρuAu=11+εi(x,t),
where u means the unstretched state of the web.Combining (3) with (4), we can obtain
(5)∫XBi(t1)XBi(t2)11+εi(x,t)dx=∫t1t2ViR(t)1+εi-1(x,t)dt.
Under assumption (5), we can obtain
(6)11+εi(x,t)=11+εi(t)≈1-εi(t).
Combining (5) with (6), we can get
(7)(XBi(t2)-XBi(t1))(1-εi(t))=∫t1t2(1-εi-1(t))ViR(t)dt.
As shown in Figure 2(a), the initial phase angle of the color mark on the printing roller ci is negative and can be calculated by
(8)αi(t)=αi(0)+∫0tωi(t)dt.
The coordinate of color mark on the printing roller ciXFi can be calculated by (9) including two parts, XFiR and XFiT, which are caused by rotation and virtual translational movement of the printing roller ci, respectively:
(9)XFi(t)=XFiR(t)+XFiT(t)=XFiR(0)+XFiT(0)+∫0t[ViR(t)+ViT(t)]dt.
Under assumption (6), the following equation could be obtained
(10)XFiR(0)=Riαi(0)≤0,XFiT(0)=0.
So, (9) can be rewritten as follows:
(11)XFi(t)=Riαi(0)+∫0t[ViR(t)+ViT(t)]dt.
At t = t1, as shown in Figure 2(b), the color mark on the printing roller ci is printed on the web, the following equation can be gotten
(12)αi(t1)=αi(0)+∫0t1ωi(t)dt=0.
We can rewrite (12) as follows:
(13)∫0t1ωi(t)dt=-αi(0)⟹∫0t1ViR(t)dt=-Riαi(0).
The coordinate of the color mark on the printing roller ci coincides with that of the color mark on the web at t = t1, so combining (11) with (13), we can get the following:
(14)XFi(t1)=XBi(t1)=∫0t1ViT(t)dt.
Because the ith unit is MDR error reference, the printing roller ci has no virtual translational movement, so ViT(t)=0, and we can obtain
(15)XBi(t1)=∫0t1ViT(t)dt=0.
Hence, (7) can be rewritten as follows:
(16)XBi(t2)=11-εi(t)∫t1t2(1-εi-1(t))ViR(t)dt.
Similarly, as shown in Figure 2, the coordinate of the color mark on the printing roller ci+1 can be derived by.(17)XFi+1(t)=XFi+1R(t)+XFi+1T(t)=XFi+1R(0)+XFi+1T(0)+∫0t[Vi+1R(t)+Vi+1T(t)dt].
At t = 0, (18) can be obtained as follows:
(18)XFi+1T(0)=Lui.
Derived with the same procedure as XFi, the following conclusion can be obtained:
(19)∫0t2Vi+1R(t)dt=-XFi+1R(0).
Combining (17), (18) with (19), we can get
(20)XFi+1(t2)=XBi+1(t2)=Lui+∫0t2Vi+1T(t)dt.
From the literature [12], we can obtain
(21)[1-εi(t)]dLi(t)dt-Li(t)dεi(t)dt=ViR(t)-Vi+1R(t)-εi-1(t)ViR(t)+εi(t)Vi+1R(t).
Since dεi(t)/dt≈0, (21) can be rewritten as follows:
(22)dLi(t)dt=1-εi-1(t)1-εi(t)ViR(t)-Vi+1R(t).
Because the change of the web length mapped to the coordinate is equivalent to the speed difference between the virtual translational movements of the two printing rollers, it can be expressed as follows:
(23)dLi(t)dt=Vi+1T(t)-ViT(t).
Due to ViT(t)=0, combining (22) with (23), (24) can be obtained as follows:
(24)Vi+1T(t)=1-εi-1(t)1-εi(t)ViR(t)-Vi+1R(t).
Combining (20) with (24), (20) can be rewritten as follows:
(25)XFi+1(t2)=Lui+∫0t2[1-εi-1(t)1-εi(t)ViR(t)-Vi+1R(t)]dt.
By using assumptions (5), we can get the following:
(26)11-εi(t)≈1+εi(t),εi-1(t)εi(t)≈0.
Therefore, (25) can be simplified as follows:
(27)XFi+1(t2)=Lui+∫0t2[ViR(t)-Vi+1R(t)mmmmmmmmmmm-(εi-1(t)-εi(t))ViR(t)]dt.
According to the definition of MRD error, and combining (1), (16), and (27), the MDR error between the ith and unit (i+1)th printing units at the instant t2 can be calculated as follows:
(28)ei(t2)=11-εi(t)∫t1t2(1-εi-1(t))ViR(t)dt-Lui-∫0t2[ViR(t)-Vi+1R(t)-(εi-1(t)-εi(t))ViR(t)]dt.
The running time tTi = t2-t1 represents the time that a color mark on the web takes to move from printing roller ci to printing roller ci+1, so t1 = t2-tTi, (28) can be simplified as follows:
(29)ei(t2)=∫0t2(Vi+1R(t)-εi(t)ViR(t))dt-∫0t2-tTi(1-εi-1(t))ViR(t)dt-Lui.
Extending (28) to any time interval, then we can get
(30)ei(t)=∫0t[Vi+1R(t)-εi(t)ViR(t)]dt-∫0t-tTi(1-εi-1(t))ViR(t)dt-Lui.
The differential form of MDR error between the ith and (i+1)th units can be derived as follows:
(31)dei(t)dt=Vi+1R(t)-ViR(t-tTi)+εi-1(t-tTi)ViR(t-tTi)-εi(t)ViR(t).
In the shaft-less drive printing system, there is no change of the web length caused by compensation rollers, so ΔLi(t) is mainly caused by load cells, dryer systems, and other disturbances. In practice, ΔLi(t) ≪Lui, so we can obtain
(32)Li(t)=Lui+ΔLi(t)≈Lui.
In (21), we substitute Lui for Li(t), and then (21) can be rewritten as follows:
(33)Luidεi(t)dt=[1-εi(t)]Vi+1R(t)-[1-εi-1(t)]ViR(t).
Combining (31) with (33), the MDR error model of the two-color system can be obtained as follows:(34)dei(t)dt=Vi+1R(t)-ViR(t-tTi)+εi-1(t-tTi)ViR(t-tTi)-εi(t)ViR(t),Luidεi(t)dt=[1-εi(t)]Vi+1R(t)-[1-εi-1(t)]ViR(t).
Considering ViR(t)=Riωi(t) and εi(t)=Ti(t)/AE, we can obtain
(35)AEdei(t)dt=AE(Ri+1ωi+1(t)-Riωi(t-tTi))+Ri(Ti-1(t-tTi)ωi(t-tTi)-Ti(t)ωi(t)),LuidTi(t)dt=AE(Ri+1ωi+1(t)-Riωi(t))+Ti-1(t)Riωi(t)-Ti(t)Ri+1ωi+1(t).
In the practical printing process of gravure printing machines, the change of the web cross-section in different web spans is very small and negligible, so we substitute average A for Ai(t) in (35). Equation (35) shows that the MDR mathematical model of the two-color gravure printing machines is expressed by two equations representing the dynamics of tension fluctuation between units and the dynamics of the register error fluctuation, respectively.


### 3.3. MDR System Model of then-Color Gravure Printing Machines

The schematic diagram of then-color MDR system is presented in Figure 3. T0(t) is the tension of the infeeding zone at steady state and can be measured by load cells.Figure 3
Schematic diagram of then-color gravure printing machine.According to (35) and Figure 3, the MDR system model of the n-color gravure printing machines can be expressed as follows:
(36)AEde1(t)dt=AE(R2ω2(t)-R1ω1(t-tT1))+R1(T0(t-tT1)ω1(t-tT1)-T1(t)ω1(t)),Lu1dT1(t)dt=AE(R2ω2(t)-R1ω1(t))+T0(t)R1ω1(t)-T1(t)R2ω2(t),⋮AEden-1(t)dt=AE(Rnωn(t)-Rn-1ωn-1(t-tT(n-1)))+Rn-1(Tn-2(t-tT(n-1))ωn-1(t-tT(n-1))+Rn-1-Tn-1(t)ωn-1(t)),Lu(n-1)dTn-1(t)dt=AE(Rnωn(t)-Rn-1ωn-1(t))+Tn-2(t)Rn-1ωn-1(t)-Tn-1(t)Rnωn(t).
The ωi(t) is composed of a stationary component ω* which corresponds to web feed speed and a variable component Δωi(t) determined by the MDR control input. In particular, ω1(t) is used as a speed reference, so we can get (37) as follows:
(37)ω1(t)=ω*,ω2(t)=ω*+Δω2(t),⋮ωn(t)=ω*+Δωn(t).
In (36), tTi is also known as time delay and is defined as follows:
(38)tTi=LuiRiωi≈LuiRiω*.
Equations (36) and (37) show that then-color MDR system in the gravure printing machine is a multi-input multioutput, nonlinear, strong-coupling, and time-delay system. In the system, ω2(t)~ωn(t) are the input signals, and e1(t)~en-1(t) are the output signals.


## 3.1. Basic Assumption

In the gravure printing machines, the MDR error of a moving web is defined in two adjacent printing rollers as the relative difference of the distance between the previous printed pattern by the upstream printing roller and the later printed one in the downstream printing roller. To remove unnecessary interference and simplify the mathematical model, the following assumptions are made in the derivation of model.(1)
There is no slippage between web material and rollers and no mass transfer between the web material and the environment.(2)
The web material is isotropic, and its properties do not change with temperature and humidity.(3)
The web cross-section in the unstretched state does not vary along machine direction in a web span.(4)
The density and modulus of elasticity of the web in the unstretched state are constant in the cross-section.(5)
The strain in web is small (much less than unity) and uniform within the web span.(6)
The printing system is in a stable state, in which there are no disturbances from tension and other factors in timet≤0.


## 3.2. MDR Error Model of the Two-Color System

The mathematical MDR error model of the two-color system in shaft-less drive gravure printing machines is illustrated in Figure2 by the coupling between the ith and (i+1)th units, and the ith unit is register errors reference. As shown in Figure 2, a coordinate axis along the machine direction is set up and the position of the ith unit is set as the origin of the coordinate axis. Figures 2(b) and 2(c) show that the color marks on printing roller ci and ci+1 are printed on the web at the printing instants t1 and t2, respectively.Figure 2
MDR error model of the two-color system.The MDR error in the moving coordinate system is to calculate the difference between the printing point on the roller and the printed point on the web representing the same point of the pattern at the printing instants. Therefore, according to the moving coordinate system, the MDR error is just the deviation from the coordinate of color mark on the printing rollerXFi+1 and that on the web XBi. At instant t2, the MDR error is defined by
(1)ei(t2)=XBi(t2)-XFi+1(t2).
It is to be noted that the movement of the printing roller can be viewed as a combination of two movements: one is the rotation of the roller and the other is the virtual translational movement of the roller along the machine direction which is caused by the variation of the web length. This variation can be affected by many factors such as dryer system, tension sensors, or some other disturbances.By using assumptions (1) and (2), the following conclusion can be obtained:
(2)ddt∫0Liρi(x,t)Ai(x,t)dx=ρi-1(x,t)Ai-1(x,t)ViR(t)-ρi(x,t)Ai(x,t)Vi+1R(t).
Hence, the mass of the web material entering the printing unit i in the time interval from t1 tot2 can be derived as follows:
(3)ddt∫XBi(t1)XBi(t2)ρi(x,t)Ai(x,t)dx=ρi-1(x,t)Ai-1(x,t)ViR(t).
According to assumptions (3) and (4) and the definition of strain, the relationship between the stretched and un-stretched webs can be written as follows:
(4)ρi(x,t)Ai(x,t)ρuAu=11+εi(x,t),
where u means the unstretched state of the web.Combining (3) with (4), we can obtain
(5)∫XBi(t1)XBi(t2)11+εi(x,t)dx=∫t1t2ViR(t)1+εi-1(x,t)dt.
Under assumption (5), we can obtain
(6)11+εi(x,t)=11+εi(t)≈1-εi(t).
Combining (5) with (6), we can get
(7)(XBi(t2)-XBi(t1))(1-εi(t))=∫t1t2(1-εi-1(t))ViR(t)dt.
As shown in Figure 2(a), the initial phase angle of the color mark on the printing roller ci is negative and can be calculated by
(8)αi(t)=αi(0)+∫0tωi(t)dt.
The coordinate of color mark on the printing roller ciXFi can be calculated by (9) including two parts, XFiR and XFiT, which are caused by rotation and virtual translational movement of the printing roller ci, respectively:
(9)XFi(t)=XFiR(t)+XFiT(t)=XFiR(0)+XFiT(0)+∫0t[ViR(t)+ViT(t)]dt.
Under assumption (6), the following equation could be obtained
(10)XFiR(0)=Riαi(0)≤0,XFiT(0)=0.
So, (9) can be rewritten as follows:
(11)XFi(t)=Riαi(0)+∫0t[ViR(t)+ViT(t)]dt.
At t = t1, as shown in Figure 2(b), the color mark on the printing roller ci is printed on the web, the following equation can be gotten
(12)αi(t1)=αi(0)+∫0t1ωi(t)dt=0.
We can rewrite (12) as follows:
(13)∫0t1ωi(t)dt=-αi(0)⟹∫0t1ViR(t)dt=-Riαi(0).
The coordinate of the color mark on the printing roller ci coincides with that of the color mark on the web at t = t1, so combining (11) with (13), we can get the following:
(14)XFi(t1)=XBi(t1)=∫0t1ViT(t)dt.
Because the ith unit is MDR error reference, the printing roller ci has no virtual translational movement, so ViT(t)=0, and we can obtain
(15)XBi(t1)=∫0t1ViT(t)dt=0.
Hence, (7) can be rewritten as follows:
(16)XBi(t2)=11-εi(t)∫t1t2(1-εi-1(t))ViR(t)dt.
Similarly, as shown in Figure 2, the coordinate of the color mark on the printing roller ci+1 can be derived by.(17)XFi+1(t)=XFi+1R(t)+XFi+1T(t)=XFi+1R(0)+XFi+1T(0)+∫0t[Vi+1R(t)+Vi+1T(t)dt].
At t = 0, (18) can be obtained as follows:
(18)XFi+1T(0)=Lui.
Derived with the same procedure as XFi, the following conclusion can be obtained:
(19)∫0t2Vi+1R(t)dt=-XFi+1R(0).
Combining (17), (18) with (19), we can get
(20)XFi+1(t2)=XBi+1(t2)=Lui+∫0t2Vi+1T(t)dt.
From the literature [12], we can obtain
(21)[1-εi(t)]dLi(t)dt-Li(t)dεi(t)dt=ViR(t)-Vi+1R(t)-εi-1(t)ViR(t)+εi(t)Vi+1R(t).
Since dεi(t)/dt≈0, (21) can be rewritten as follows:
(22)dLi(t)dt=1-εi-1(t)1-εi(t)ViR(t)-Vi+1R(t).
Because the change of the web length mapped to the coordinate is equivalent to the speed difference between the virtual translational movements of the two printing rollers, it can be expressed as follows:
(23)dLi(t)dt=Vi+1T(t)-ViT(t).
Due to ViT(t)=0, combining (22) with (23), (24) can be obtained as follows:
(24)Vi+1T(t)=1-εi-1(t)1-εi(t)ViR(t)-Vi+1R(t).
Combining (20) with (24), (20) can be rewritten as follows:
(25)XFi+1(t2)=Lui+∫0t2[1-εi-1(t)1-εi(t)ViR(t)-Vi+1R(t)]dt.
By using assumptions (5), we can get the following:
(26)11-εi(t)≈1+εi(t),εi-1(t)εi(t)≈0.
Therefore, (25) can be simplified as follows:
(27)XFi+1(t2)=Lui+∫0t2[ViR(t)-Vi+1R(t)mmmmmmmmmmm-(εi-1(t)-εi(t))ViR(t)]dt.
According to the definition of MRD error, and combining (1), (16), and (27), the MDR error between the ith and unit (i+1)th printing units at the instant t2 can be calculated as follows:
(28)ei(t2)=11-εi(t)∫t1t2(1-εi-1(t))ViR(t)dt-Lui-∫0t2[ViR(t)-Vi+1R(t)-(εi-1(t)-εi(t))ViR(t)]dt.
The running time tTi = t2-t1 represents the time that a color mark on the web takes to move from printing roller ci to printing roller ci+1, so t1 = t2-tTi, (28) can be simplified as follows:
(29)ei(t2)=∫0t2(Vi+1R(t)-εi(t)ViR(t))dt-∫0t2-tTi(1-εi-1(t))ViR(t)dt-Lui.
Extending (28) to any time interval, then we can get
(30)ei(t)=∫0t[Vi+1R(t)-εi(t)ViR(t)]dt-∫0t-tTi(1-εi-1(t))ViR(t)dt-Lui.
The differential form of MDR error between the ith and (i+1)th units can be derived as follows:
(31)dei(t)dt=Vi+1R(t)-ViR(t-tTi)+εi-1(t-tTi)ViR(t-tTi)-εi(t)ViR(t).
In the shaft-less drive printing system, there is no change of the web length caused by compensation rollers, so ΔLi(t) is mainly caused by load cells, dryer systems, and other disturbances. In practice, ΔLi(t) ≪Lui, so we can obtain
(32)Li(t)=Lui+ΔLi(t)≈Lui.
In (21), we substitute Lui for Li(t), and then (21) can be rewritten as follows:
(33)Luidεi(t)dt=[1-εi(t)]Vi+1R(t)-[1-εi-1(t)]ViR(t).
Combining (31) with (33), the MDR error model of the two-color system can be obtained as follows:(34)dei(t)dt=Vi+1R(t)-ViR(t-tTi)+εi-1(t-tTi)ViR(t-tTi)-εi(t)ViR(t),Luidεi(t)dt=[1-εi(t)]Vi+1R(t)-[1-εi-1(t)]ViR(t).
Considering ViR(t)=Riωi(t) and εi(t)=Ti(t)/AE, we can obtain
(35)AEdei(t)dt=AE(Ri+1ωi+1(t)-Riωi(t-tTi))+Ri(Ti-1(t-tTi)ωi(t-tTi)-Ti(t)ωi(t)),LuidTi(t)dt=AE(Ri+1ωi+1(t)-Riωi(t))+Ti-1(t)Riωi(t)-Ti(t)Ri+1ωi+1(t).
In the practical printing process of gravure printing machines, the change of the web cross-section in different web spans is very small and negligible, so we substitute average A for Ai(t) in (35). Equation (35) shows that the MDR mathematical model of the two-color gravure printing machines is expressed by two equations representing the dynamics of tension fluctuation between units and the dynamics of the register error fluctuation, respectively.


## 3.3. MDR System Model of then-Color Gravure Printing Machines

The schematic diagram of then-color MDR system is presented in Figure 3. T0(t) is the tension of the infeeding zone at steady state and can be measured by load cells.Figure 3
Schematic diagram of then-color gravure printing machine.According to (35) and Figure 3, the MDR system model of the n-color gravure printing machines can be expressed as follows:
(36)AEde1(t)dt=AE(R2ω2(t)-R1ω1(t-tT1))+R1(T0(t-tT1)ω1(t-tT1)-T1(t)ω1(t)),Lu1dT1(t)dt=AE(R2ω2(t)-R1ω1(t))+T0(t)R1ω1(t)-T1(t)R2ω2(t),⋮AEden-1(t)dt=AE(Rnωn(t)-Rn-1ωn-1(t-tT(n-1)))+Rn-1(Tn-2(t-tT(n-1))ωn-1(t-tT(n-1))+Rn-1-Tn-1(t)ωn-1(t)),Lu(n-1)dTn-1(t)dt=AE(Rnωn(t)-Rn-1ωn-1(t))+Tn-2(t)Rn-1ωn-1(t)-Tn-1(t)Rnωn(t).
The ωi(t) is composed of a stationary component ω* which corresponds to web feed speed and a variable component Δωi(t) determined by the MDR control input. In particular, ω1(t) is used as a speed reference, so we can get (37) as follows:
(37)ω1(t)=ω*,ω2(t)=ω*+Δω2(t),⋮ωn(t)=ω*+Δωn(t).
In (36), tTi is also known as time delay and is defined as follows:
(38)tTi=LuiRiωi≈LuiRiω*.
Equations (36) and (37) show that then-color MDR system in the gravure printing machine is a multi-input multioutput, nonlinear, strong-coupling, and time-delay system. In the system, ω2(t)~ωn(t) are the input signals, and e1(t)~en-1(t) are the output signals.


## 4. Simulation and Experiments

### 4.1. Experimental Setup and MDR Simulation Model of the Four-Color System

To investigate the reliability of the derived MDR mathematical model, numerical simulation and experiments under open-loop control are carried out by using MATLAB and experimental setup for the four-color gravure printing machines. Figures4 and 5 show the setup configuration of the shaft-less drive gravure printing machine.Figure 4
Experimental setup of the four-color gravure printing machines.Figure 5
Control configuration of the experimental setup.The system includes an unwinding unit, an infeeding unit, an outfeeding unit, a rewinding unit, and a four-color printing subsystem, and all of the driving shafts are driven by independent servo motors. The printing rollers are driven by servo motors (YASKAWA SGMGH-44), and printing speed is obtained by the encoder on the servo motor. The register controller (KESAI ST3000) works in shaft-less drive mode, and the MDR errors are measured using photoelectric sensors (XINMAO XM08). Tension signals are tested by load cells (MITSUBISHI LX-030TD) and tension measuring amplifiers (MITSUBISHI LM-10TA) mounted on the shaft of idle rollers. The shaft-less drive system is controlled by multiaxis controller (GOOGOLTECH T8VME). The experimental web material is biaxially oriented polypropylene (BOPP). As shown in Figure5, tension of the printing subsystem is same with that in the infeeding zone at steady state and is not controlled; in other words, T1(t), T2(t), and T3(t) are not controlled. In the infeeding unit, T0(t) is adjusted by changing the speed of the roller 2. After preregistration, e1(t), e2(t), and e3(t) work in the open-loop control mode and are sampled once per revolution of the printing rollers. The parameters of the experimental setup are summarized in Table 1.Table 1
System parameters.
Parameters
Values
Units
Entry tension to 1th unit (T0)
100
N
Web length (Lui)
9.1
m
Radius of the printing roller (Ri)
0.2
m
Web cross-section area (A)
2.0
×
10
-
5
m2
Modulus of elasticity of the web (E)
160
MPaAccording to Table1, we can get
(39)R1=R2=R3=R4=R,Lu1=Lu2=Lu3=Lu,tT1=tT2=tT3=tT.
Hence, combining (36) with (39), the MDR system model of the four-color gravure printing machines in the simulation can be expressed as follows:
(40)AEde1(t)dt=AER(ω2(t)-ω1(t-tT))+R(T0(t-tT)ω1(t-tT)-T1(t)ω1(t)),AEde2(t)dt=AER(ω3(t)-ω2(t-tT))+R(T1(t-tT)ω2(t-tT)-T2(t)ω2(t)),AEde3(t)dt=AER(ω4(t)-ω3(t-tT))+R(T2(t-tT)ω3(t-tT)-T3(t)ω3(t)),LudT1(t)dt=AER(ω2(t)-ω1(t))+R(T0(t)ω1(t)-T1(t)ω2(t)),LudT2(t)dt=AER(ω3(t)-ω2(t))+R(T1(t)ω2(t)-T2(t)ω3(t)),LudT3(t)dt=AER(ω4(t)-ω3(t))+R(T2(t)ω3(t)-T3(t)ω4(t)).


### 4.2. Results and Analysis of the Simulation and Experiments

In order to ensure the consistency of the simulation and experiments, the parameters used in the numerical simulation are taken from Table1. Because the simulation and experiments under open-loop control are carried out, the speed of all servo motors is the same in the four-color printing system. In order to verify the developed MDR system model of the n-color gravure printing machines, T0(t) has a step change from 100 N to 110 N at t=5 s when ω1 is equal to 100 r/min and 200 r/min. The responses of the MDR errors and tension signals in simulation and experiments are shown in Figures 6 and 7, respectively.Open-loop dynamics of the MDR errors in simulation.
(a)
Tension dynamics atω1=100r/min
(b)
Tension dynamics atω1=200r/min
(c)
MDR errors dynamics atω1=100r/min
(d)
MDR errors dynamics atω1=200r/minOpen-loop dynamics of the MDR errors in experiments.
(a)
Tension dynamics atω1=100r/min
(b)
Tension dynamics atω1=200r/min
(c)
MDR errors dynamics atω1=100r/min
(d)
MDR errors dynamics atω1=200r/minAs shown in Figure6, the step change of T0(t) causes the step change of T1(t), T2(t), and T3(t), and with an increase in ω1(t), the time thatT1(t), T2(t), and T3(t) take to reach the steady state value is gradually shortened. For example, when ω1(t) is equal to 100 r/min and 200 r/min, T1(t) reaches 110 N at 35 s and 20 s, respectively. The step change of T0(t) also directly causes the change of e1(t) and then the changes of e2(t) and e3(t) by coupling characteristics, and as the speed increases, the interval from the peak of ei(t) to the peak of ei+1 (t) is gradually shortened; for example, when ω1(t) is equal to 100 r/min and 200 r/min, the interval from the peak of e1(t) to the peak of e2(t) is 3.11 s and 1.56 s, respectively.Figures6 and 7 illustrate that although there are some disturbances in the experiments, the experimental results are almost coincident with the simulation curves. For example, with an increase in ω1(t), the time that tension signals take to reach the steady state value is gradually shortened and so on. Comparing with Figures 6(c), 6(d), 7(c), and 7(d), respectively, the peak of ei(t) in the experiments is slightly more than that in the simulation, such as when ω1(t) is equal to 200 r/min, the peaks of e1(t) in the experiments and simulation are −0.82 mm and −0.78 mm, respectively. There are two primary reasons: one is that there are external disturbances from environment in the experiments, and the other is that there are errors between the simulation parameters and real experimental parameters.The results of the simulation and experiments indicate that the established MDR system model is accurate and reliable and can be used to design a MDR controller. At the same time, the results also offer a powerful proof that the MDR system of the multicolor gravure printing machines is a multi-input multioutput, nonlinear, strong-coupling, and time-delay system.


## 4.1. Experimental Setup and MDR Simulation Model of the Four-Color System

To investigate the reliability of the derived MDR mathematical model, numerical simulation and experiments under open-loop control are carried out by using MATLAB and experimental setup for the four-color gravure printing machines. Figures4 and 5 show the setup configuration of the shaft-less drive gravure printing machine.Figure 4
Experimental setup of the four-color gravure printing machines.Figure 5
Control configuration of the experimental setup.The system includes an unwinding unit, an infeeding unit, an outfeeding unit, a rewinding unit, and a four-color printing subsystem, and all of the driving shafts are driven by independent servo motors. The printing rollers are driven by servo motors (YASKAWA SGMGH-44), and printing speed is obtained by the encoder on the servo motor. The register controller (KESAI ST3000) works in shaft-less drive mode, and the MDR errors are measured using photoelectric sensors (XINMAO XM08). Tension signals are tested by load cells (MITSUBISHI LX-030TD) and tension measuring amplifiers (MITSUBISHI LM-10TA) mounted on the shaft of idle rollers. The shaft-less drive system is controlled by multiaxis controller (GOOGOLTECH T8VME). The experimental web material is biaxially oriented polypropylene (BOPP). As shown in Figure5, tension of the printing subsystem is same with that in the infeeding zone at steady state and is not controlled; in other words, T1(t), T2(t), and T3(t) are not controlled. In the infeeding unit, T0(t) is adjusted by changing the speed of the roller 2. After preregistration, e1(t), e2(t), and e3(t) work in the open-loop control mode and are sampled once per revolution of the printing rollers. The parameters of the experimental setup are summarized in Table 1.Table 1
System parameters.
Parameters
Values
Units
Entry tension to 1th unit (T0)
100
N
Web length (Lui)
9.1
m
Radius of the printing roller (Ri)
0.2
m
Web cross-section area (A)
2.0
×
10
-
5
m2
Modulus of elasticity of the web (E)
160
MPaAccording to Table1, we can get
(39)R1=R2=R3=R4=R,Lu1=Lu2=Lu3=Lu,tT1=tT2=tT3=tT.
Hence, combining (36) with (39), the MDR system model of the four-color gravure printing machines in the simulation can be expressed as follows:
(40)AEde1(t)dt=AER(ω2(t)-ω1(t-tT))+R(T0(t-tT)ω1(t-tT)-T1(t)ω1(t)),AEde2(t)dt=AER(ω3(t)-ω2(t-tT))+R(T1(t-tT)ω2(t-tT)-T2(t)ω2(t)),AEde3(t)dt=AER(ω4(t)-ω3(t-tT))+R(T2(t-tT)ω3(t-tT)-T3(t)ω3(t)),LudT1(t)dt=AER(ω2(t)-ω1(t))+R(T0(t)ω1(t)-T1(t)ω2(t)),LudT2(t)dt=AER(ω3(t)-ω2(t))+R(T1(t)ω2(t)-T2(t)ω3(t)),LudT3(t)dt=AER(ω4(t)-ω3(t))+R(T2(t)ω3(t)-T3(t)ω4(t)).


## 4.2. Results and Analysis of the Simulation and Experiments

In order to ensure the consistency of the simulation and experiments, the parameters used in the numerical simulation are taken from Table1. Because the simulation and experiments under open-loop control are carried out, the speed of all servo motors is the same in the four-color printing system. In order to verify the developed MDR system model of the n-color gravure printing machines, T0(t) has a step change from 100 N to 110 N at t=5 s when ω1 is equal to 100 r/min and 200 r/min. The responses of the MDR errors and tension signals in simulation and experiments are shown in Figures 6 and 7, respectively.Open-loop dynamics of the MDR errors in simulation.
(a)
Tension dynamics atω1=100r/min
(b)
Tension dynamics atω1=200r/min
(c)
MDR errors dynamics atω1=100r/min
(d)
MDR errors dynamics atω1=200r/minOpen-loop dynamics of the MDR errors in experiments.
(a)
Tension dynamics atω1=100r/min
(b)
Tension dynamics atω1=200r/min
(c)
MDR errors dynamics atω1=100r/min
(d)
MDR errors dynamics atω1=200r/minAs shown in Figure6, the step change of T0(t) causes the step change of T1(t), T2(t), and T3(t), and with an increase in ω1(t), the time thatT1(t), T2(t), and T3(t) take to reach the steady state value is gradually shortened. For example, when ω1(t) is equal to 100 r/min and 200 r/min, T1(t) reaches 110 N at 35 s and 20 s, respectively. The step change of T0(t) also directly causes the change of e1(t) and then the changes of e2(t) and e3(t) by coupling characteristics, and as the speed increases, the interval from the peak of ei(t) to the peak of ei+1 (t) is gradually shortened; for example, when ω1(t) is equal to 100 r/min and 200 r/min, the interval from the peak of e1(t) to the peak of e2(t) is 3.11 s and 1.56 s, respectively.Figures6 and 7 illustrate that although there are some disturbances in the experiments, the experimental results are almost coincident with the simulation curves. For example, with an increase in ω1(t), the time that tension signals take to reach the steady state value is gradually shortened and so on. Comparing with Figures 6(c), 6(d), 7(c), and 7(d), respectively, the peak of ei(t) in the experiments is slightly more than that in the simulation, such as when ω1(t) is equal to 200 r/min, the peaks of e1(t) in the experiments and simulation are −0.82 mm and −0.78 mm, respectively. There are two primary reasons: one is that there are external disturbances from environment in the experiments, and the other is that there are errors between the simulation parameters and real experimental parameters.The results of the simulation and experiments indicate that the established MDR system model is accurate and reliable and can be used to design a MDR controller. At the same time, the results also offer a powerful proof that the MDR system of the multicolor gravure printing machines is a multi-input multioutput, nonlinear, strong-coupling, and time-delay system.


## 5. Conclusions

In the gravure printing machines, the register accuracy is the most important in performance. Therefore, in order to improve the control accuracy of MDR system, a nonlinear mathematical model of the MDR system in shaft-less drive mode is presented and derived in detail for the multicolor gravure printing machines in this paper. First, we introduce the schematic diagram of the shaft-less drive gravure printing machines. Secondly, we make a deep analysis of the causes and influencing factors of the MDR errors, then establish the MDR error of the two-color system, and further obtain the MDR system model of then-color gravure printing machines. At last, the proposed model is verified by the numerical simulations and experimental studies. Therefore, it can be concluded that the MDR system model is accurate enough to describe dynamics of the MDR errors in multicolor printing system and is worthwhile to be developed as an effective MDR controller in future study.



---

*Source: 251636-2013-11-04.xml*