Spatially dependent transfer functions for web span lateral dynamics which provide web lateral position and slope as outputs at any location in the web span are derived in this paper. The proposed approach overcomes one of the key limitations of the existing methods which provide web lateral position only on the rollers. The approach relies on taking the Laplace transform with respect to the temporal variable of both the web span lateral governing equation and the boundary conditions on the rollers, and solving the resulting equations. A general web span lateral transfer function, which is an explicit function of the spatial position along the span, is obtained first followed by its application to common guide configurations. The approach also significantly simplifies the consideration of shear (relevant to short spans), in addition to bending, which has been found to be difficult to handle in past studies. We first develop spatially dependent lateral transfer functions by considering only bending which is relevant to most web handling situations, and then add shear to the formulation and develop spatially dependent lateral transfer functions that include both bending and shear. Results from model simulations and pertinent discussions are provided. The spatially dependent transfer functions derived in this paper are a significant improvement over existing lateral transfer functions and provide mechanisms to analyze web lateral behavior within spans, study propagation of lateral disturbances, and aid in the development of closed-loop lateral control systems in emerging applications that require precise lateral positioning of the web.

## Introduction

In Roll-to-Roll (R2R) manufacturing, flexible materials called webs are transported on rollers through processes (such as printing, coating, heat treatment, lamination, etc.). Studies in the literature have mostly focused on modeling and control of moving webs in the longitudinal or transport direction. In many applications, control of lateral web motion (motion perpendicular to the transport direction and in the plane of the web) has been mostly relegated to just keeping the web on rollers during transport. Increased use of R2R manufacturing in recent years on a variety of polymer materials under different processing conditions, in both conventional products and emerging products in flexible electronics, has led to additional requirements on the control of lateral motion of the web. For example, in R2R printing where multiple print cylinders are employed to sequentially print and register patterns on the web [1], there have been stringent requirements on minimizing print registration in both longitudinal and lateral directions.

There have been several studies on modeling the lateral behavior of moving webs. The first seminal work on the topic was reported by Shelton in his Ph.D. thesis in 1968 [2] and subsequently published in this journal [3,4]. Subsequent work in modeling and control of web lateral dynamics based on this treatment was reported in Refs. [5–9]. For the purpose of deriving the governing equations of the web lateral position on rollers, the moving web between two rollers is treated as a tensioned Euler–Bernoulli beam. For most web materials, the web mass between two rollers is negligible, i.e., the force due to acceleration of web mass is negligible when compared to web tension. Thus, the lateral motion of the web between two rollers is treated as the motion of a static beam; and the web between two rollers is treated as a tensioned beam. Four boundary conditions (web lateral position and slope on each roller) are utilized to solve the fourth-order partial differential equation describing the lateral motion of the web. A key observation/principle is utilized to setup two of the boundary conditions—a web approaching a roller aligns itself normal to the axis of rotation of the roller. This is also well known in the belt transport literature. This principle is utilized to setup two normal entry conditions: web lateral velocity and acceleration in terms of roller lateral velocity and acceleration, web entry angle at the roller, and angle of the roller. Based on this approach, transfer functions from the guide roller lateral position (input variable) to the web lateral position on the roller (controlled or output variable) are determined for various guide roller mechanisms, such as the end-pivoted guide, center-pivoted guide, offset pivot guide (OPG), remotely pivoted guide (RPG), etc. Figure 1 provides an illustration of a web span with upstream (entry) and downstream (exit) rollers and the definition of web wrap angle on a roller.

There are several limitations to the existing approach: (1) it provides a governing equation only for the lateral position behavior on the roller and not for any position within the span and (2) it does not the slope of the web which may contribute to the creation of web lateral oscillations and their propagation to downstream spans with web transport. Further, in the existing approach, the solution (lateral position) of the lateral governing equation is obtained by assuming constant downstream boundary conditions (on the downstream roller of the span); subsequently, the time derivative of the solution is used in the normal entry conditions to determine the transfer function for the lateral position on the roller. The purpose of the guide roller is to modify the downstream boundary conditions, and therefore, the assumption on the downstream boundary conditions is counter to the notion that the axis of rotation and translation of the guide roller is utilized to change the web boundary conditions on the roller. In this paper, we address these issues and derive spatially dependent lateral transfer functions that provide the web lateral position and slope at any point in the web span and not just on the roller.

Our approach relies on taking the Laplace transform of the web span lateral governing equation and the boundary conditions with respect to the temporal variable. We consider the web lateral position and slope on the upstream roller to the span as two boundary conditions and the two normal entry conditions (for lateral velocity and lateral acceleration) on the downstream roller as the other two boundary conditions. The idea of incorporating the normal entry conditions as boundary conditions has been considered in Ref. [6] where the normal entry conditions are applied to a system modeled using a dynamic beam equation; a two-dimensional Laplace transform (for both spatial and temporal variables) was applied to the dynamic beam equation which provides a solution to the beam equation in the frequency domain. Due to the complexity of determining the inverse Laplace transform of the resulting solution, the frequency domain solution was ignored and the spatial derivatives in the beam model were discretized using a finite difference method to obtain a set of ordinary differential equations for web lateral response on the roller; a two-span system example was used to illustrate the procedure. In our approach, we employ the Laplace transform in the temporal variable not only for the beam governing equation, but also for the two boundary conditions on the upstream roller and the two normal entry conditions on the downstream roller. This allows us to solve the resulting equations and obtain spatially dependent lateral transfer functions. The application of the Laplace transforms in the temporal variable for distributed parameter systems is discussed in Ref. [10].

The contributions of the paper are summarized in the following. First, we derive spatially dependent transfer functions for web lateral dynamics where the output can be any point spatially along the span. Further, we also derive transfer functions for web slope as the output at any point in the span. This is important as one can derive controllers to regulate lateral position at any point in the span, and investigate strategies that minimize propagation of lateral position and slope errors into the downstream spans. Second, the derived spatially dependent transfer functions can be applied to any guide configuration and specific transfer functions for those guiding situations can be determined. Third, we incorporate shear into the formulation by modifying the lateral governing equation and normal entry conditions appropriately. We use the same approach as before to derive the spatially dependent lateral transfer functions that include both bending and shear. The pure bending case is relevant to long spans, i.e., spans with large web span length to web width ratio, whereas the effect of shear is more prevalent in short spans. Further, since the proposed approach allows us to obtain not only position but also slope at any location in the span, one need not use multiple sensors to determine the position and slope of the web within the span as is done in some existing studies.

The rest of the paper is organized as follows: In Sec. 2, we provide a discussion of the existing approach and its assumptions and limitations. The spatially dependent transfer functions for the pure bending case are derived, and application of the spatially dependent lateral transfer function to specific guiding configurations is given in Sec. 3. Inclusion of shear into the lateral governing equation and modification of the normal entry conditions are discussed in Sec. 4; spatially dependent lateral transfer functions for lateral position and slope for combined bending and shear are also given in this section. Numerical simulations for different guiding configurations are discussed in Sec. 5. Concluding remarks and future work are provided in Sec. 6.

## Existing Approach and Assumptions

*y*) and slope ($\u2202y/\u2202x$) at the two ends of the web span, which are assumed to be known and given by

where the subscripts 0 and *L* are utilized to denote variables at the entry and exit rollers, respectively, of the web span under consideration. Note that these boundary conditions imply that both ends of the tensioned beam (web span) are free.

Note that the normal entry rule is used as a mechanism by which a guide roller can control the lateral position of the web via the rotation and translation of the guide roller by an actuating mechanism.

*x*is substituted into Eq. (5). The resulting equation contains the web angle,

*θ*, which is replaced by the slope term from the entry rule given by Eq. (4). Taking the Laplace transform of the resulting equation with respect to time results in the following lateral response at the downstream roller due to various inputs [7]:

_{wL}*Remark 2.1.* This existing approach has several drawbacks: (1) it provides the evolution of the lateral position only on the roller; it does not provide lateral web position within the span. (2) Although lateral position at the downstream roller is of interest, this is assumed as a known boundary condition in the development. (3) The normal entry conditions are used in an indirect manner in the sense that the knowledge of the lateral position and slope on the downstream roller are assumed to be known (resulting in the response at a specific location) and then these are applied to the resulting solution to fit the normal entry conditions.

## Spatially Dependent Transfer Functions

In the existing approach (summarized in Sec. 2), the effect of the boundary conditions on web lateral behavior within the free span is not clear. The “normal entry rule” is applied to introduce a dynamical behavior on the roller and to obtain a relationship between the lateral web position on the guide roller and the control action (guide roller motion). In this work, instead of assuming the downstream boundary conditions for lateral position and slope (at *x* = *L*) (whose evolution is of interest), the two normal entry conditions are employed as downstream boundary conditions. This allows for directly incorporating the effect of the roller into the solution of the governing equation. Further, this also allows for directly coupling the dynamic effects of the rollers with the span lateral dynamics.

We define a free span as the length of material between two rollers that is not wrapped on the rollers. For the upstream roller, we establish the boundary conditions for the span at the exit of the region of wrap of the upstream roller. Due to the application of the normal entry rule to the upstream roller, which stipulates that the web aligns perpendicular to the roller at contact, the web leaves this roller perpendicularly. Then, the upstream roller angle becomes a boundary condition for web slope. In the region of wrap of the upstream roller, the lateral displacement remains the same throughout, which is taken as the second boundary condition at the beginning of the span.

*x*= 0 represent web lateral position and web slope at the entry of the span (exit of the region of wrap of the upstream roller); the slope is the same as the upstream roller angle $\theta 0(t)$, due to the interpretation of the entry rule for the web on that roller. The other two conditions at

*x*=

*L*are the normal entry rules at the entry of the region of wrap for the downstream roller. We will apply the following Laplace transform for the time variable:

*C*are functions of the variable

_{i}*s*. Substituting the boundary conditions into Eq. (11) results in the following solution:

where $glx(x)$ is the first partial derivative of $gl(x)$ with respect to *x* ($l=1,2)$ and $glxx(x)$ is the second partial derivative. Note that $g1x(L)=1,\u2009g2x(L)=0,\u2009g1xx(L)=0,\u2009g2xx(L)=\u22121$.

*x*=

*L*, Eq. (12) can be simplified to

Equation (17) provides the spatially dependent Laplace transform of the lateral position in terms of various inputs and their associated transfer functions. Further, one can obtain transfer functions for slope and moment at any location along the web span by substituting Eq. (16) in Eqs. (14) and (15).

### Comparison With Existing Transfer Functions.

*x*=

*L*. Note that at

*x*=

*L*, we can write the following relations:

*x*=

*L*

However, in our approach, the free span does not include the region of wrap in either the upstream or downstream rollers; the interpretation of the normal entry rule is that the web will acquire the roller angle and keep it for the entire region of wrap; due to this the first term in the numerator of $y0(s)$ and $z0(s)$ cancel each other. The following remarks provide some observations based on the results of this section.

*Remark 3.1.* By employing normal entry conditions on the downstream roller of the span as the boundary conditions, we incorporated the effect of web/roller contact into the solution of the governing equation. The proposed method allows us to directly obtain higher order spatial partial derivatives of the lateral response and thus can be used for obtaining web slope, moment, shear force, etc. The method also further opens up the opportunity to develop controllers for processes which require control of lateral position within the span.

*Remark 3.2.* The proposed method can be extended to include shear by establishing appropriate boundary conditions. The inclusion of shear modifies the boundary condition for the lateral acceleration and introduces an additional pole and zero in the lateral transfer functions. This is discussed later in Sec. 4.

Equation (17) is the general expression for spatially dependent web lateral position and further simplification of this equation can be achieved by considering the specific roller configuration corresponding to a given situation. This is shown in Table 1 for the most common situations.

Setup | Conditions | Equation |
---|---|---|

Downstream fixed roller | $\theta \u0302L(s)=0,\u2009z\u0302L(s)=0$ | $y\u0302(x,s)=P1(x,s)Db(s)\theta \u03020(s)+P2(x,s)Db(s)y\u03020(s)$ |

Pure displacement | $\theta \u0302L(s)=0$ | $y\u0302(x,s)=P1(x,s)Db(s)\theta \u03020(s)+P2(x,s)Db(s)y\u03020(s)+P4(x,s)Db(s)z\u0302L(s)$ |

Pure rotation | $z\u0302L(s)=0$ | $y\u0302(x,s)=P1(x,s)Db(s)\theta \u03020(s)+P2(x,s)Db(s)y\u03020(s)+P3(x,s)Db(s)\theta \u0302L(s)$ |

RPG (see Fig. 3) | $z\u0302L(s)=X1\theta \u0302L(s)$ | $y\u0302(x,s)=P1(x,s)Db(s)\theta \u03020(s)+P2(x,s)Db(s)y\u03020(s)+P5(x,s)Db(s)\theta \u0302L(s)$ where $P5(x,s)=P3(x,s)+X1P4(x,s)$ |

OPG (see Fig. 4) | $\theta \u03020(s)=\theta \u0302L(s)=\theta \u0302(s)$ | $y\u0302(x,s)=P6(x,s)Db(s)\theta \u0302(s)+P2(x,s)Db(s)y\u03020(s)$ where $P6(x,s)=P1(x,s)+P3(x,s)+X1P4(x,s)$ |

Setup | Conditions | Equation |
---|---|---|

Downstream fixed roller | $\theta \u0302L(s)=0,\u2009z\u0302L(s)=0$ | $y\u0302(x,s)=P1(x,s)Db(s)\theta \u03020(s)+P2(x,s)Db(s)y\u03020(s)$ |

Pure displacement | $\theta \u0302L(s)=0$ | $y\u0302(x,s)=P1(x,s)Db(s)\theta \u03020(s)+P2(x,s)Db(s)y\u03020(s)+P4(x,s)Db(s)z\u0302L(s)$ |

Pure rotation | $z\u0302L(s)=0$ | $y\u0302(x,s)=P1(x,s)Db(s)\theta \u03020(s)+P2(x,s)Db(s)y\u03020(s)+P3(x,s)Db(s)\theta \u0302L(s)$ |

RPG (see Fig. 3) | $z\u0302L(s)=X1\theta \u0302L(s)$ | $y\u0302(x,s)=P1(x,s)Db(s)\theta \u03020(s)+P2(x,s)Db(s)y\u03020(s)+P5(x,s)Db(s)\theta \u0302L(s)$ where $P5(x,s)=P3(x,s)+X1P4(x,s)$ |

OPG (see Fig. 4) | $\theta \u03020(s)=\theta \u0302L(s)=\theta \u0302(s)$ | $y\u0302(x,s)=P6(x,s)Db(s)\theta \u0302(s)+P2(x,s)Db(s)y\u03020(s)$ where $P6(x,s)=P1(x,s)+P3(x,s)+X1P4(x,s)$ |

The following remarks provide a perspective for each item in Table 1. With a fixed downstream roller, the response depends only on the perturbations ($y\u03020(s),\u2009\theta \u03020(s)$) at the entry of the span. To regulate web lateral position, one has to control both position and slope; much of the existing work has focused only on regulating the lateral position at the exit of the guide without accounting for the roller angle. Note that in most guide installations, the plane of the guiding span and the span downstream of it are perpendicular to each other; if this is not the case, then the roller angle will affect downstream lateral position. We will illustrate this effect in subsequent model simulations. Note that the second and third items in Table 1 correspond to a pure displacement guide roller (with $z\u0302L(s)$ as the input) and pure rotation guide roller (with $\theta \u0302L(s)$ as the input); these could be used for terminal guiding on unwind and rewind rollers. The last two items correspond to two commonly used intermediate guides.

## Consideration of Shear

*n*is the correction factor and is equal to 1.2 for webs with rectangular cross-sectional area,

*A*is the web cross section area, and

*G*is the shear modulus. We will use subscript

*t*to refer to the total deflection due to combined bending and shear and subscripts

*b*and

*s*, respectively, for pure bending and pure shear. The tensioned beam governing equation that includes shear may be obtained by using the force equilibrium in the lateral direction. A free-body diagram of the web is provided in Fig. 5(a). The shear force may be expressed as

*x*and employing $\u2202N(x)/\u2202x=\u2212T\u22022yt(x,t)/\u2202x2$ (obtained from Eq. (23)), one can establish the following governing equation for the web that includes both bending and shear:

where $Ke2=T/(EI[1+nT/AG])$. Equation (27) is similar to the pure bending case, except that *K* is replaced by *K _{e}* which is related to

*K*as $Ke2=K2/[1+nT/AG]$. Note that if shear is not considered, we can set

*n*= 0 to obtain

*K*=

_{e}*K*.

*y*; these are given by

_{t}*y*would aid in determining the coefficients of the general solution of the governing equation. To determine such a relationship, the shear force expressed by Eq. (26) is used in Eq. (22). This results in

_{t}*K*with

*K*. Using the same approach as in the pure bending case, the total lateral response due to both bending and shear may be simplified to

_{e}In the previously-mentioned definitions for $P\u0303j,j=1:4$, terms that carry forward from the pure bending case are indicated. Inclusion of shear results in the additional terms in those definitions. Note that if shear is not considered, then *γ* = 0, $g3e(x)=0,\u2009g1e(x)=g1(x)$, and $g2e(x)=g2(x)$, and Eq. (39) reduces to Eq. (12)

## Numerical Simulations

Numerical simulations were conducted with the lateral transfer functions that include both bending and shear. A three-roller system shown in Fig. 6 is considered. The numerical values of the web and guide parameters used in the simulations are provided in Table 2 [7]. For the first set of numerical simulations, we consider the governing equation for the span between rollers R1 and R2 (entry span), and apply different disturbance scenarios provided in Table 3. Figure 7 provides the evolution of the lateral web position and slope in the entry span at four different locations. One can observe that the web slope changes direction as we move spatially toward roller R2. Figure 8 shows results when the roller R2 is an RPG guide roller which rotates in sinusoidal motion around a pivot point which is at a distance *X*_{1} from R2 in the entry span; note that this also induces a roller displacement $zL(t)=X1\theta L(t)$. In essence, this action mimics a combination of pure displacement and pure rotation at R2. In this case, the web slope is in the same direction as the lateral motion direction as R2 is directing motion in the entry span (note that R1 is fixed).

Definition | Symbol | Value | Units |
---|---|---|---|

Entry span length | L | 3.2808 (1) | ft (m) |

Exit span length | L | 3.806 (1.16) | ft (m) |

Integral gain | k_{i} | 10 | |

Pivoting distance for OPG | X_{1} | 3.52808 (1) | ft (m) |

Pivoting distance for RPG | X_{1} | 2.5833 (0.7874) | ft (m) |

Proportional gain | k_{p} | 90 | |

Reference tension | T | 10 (44.48) | lbf (N) |

Shear modulus | G | 0.154 (1.062 × 10^{9}) | Mpsi (Pa) |

Transport speed | v | 500 (2.54) | ft/min (m/s) |

Web thickness | h | 0.005 (0.127) | in (mm) |

Web width | W | 5.4 (137.16) | in (mm) |

Wrap angle | β | 1.553 (89) | rad (deg) |

Young's modulus | E | 0.40466 (2.76 × 10^{9}) | Mpsi (Pa) |

Definition | Symbol | Value | Units |
---|---|---|---|

Entry span length | L | 3.2808 (1) | ft (m) |

Exit span length | L | 3.806 (1.16) | ft (m) |

Integral gain | k_{i} | 10 | |

Pivoting distance for OPG | X_{1} | 3.52808 (1) | ft (m) |

Pivoting distance for RPG | X_{1} | 2.5833 (0.7874) | ft (m) |

Proportional gain | k_{p} | 90 | |

Reference tension | T | 10 (44.48) | lbf (N) |

Shear modulus | G | 0.154 (1.062 × 10^{9}) | Mpsi (Pa) |

Transport speed | v | 500 (2.54) | ft/min (m/s) |

Web thickness | h | 0.005 (0.127) | in (mm) |

Web width | W | 5.4 (137.16) | in (mm) |

Wrap angle | β | 1.553 (89) | rad (deg) |

Young's modulus | E | 0.40466 (2.76 × 10^{9}) | Mpsi (Pa) |

For the second set of numerical simulations, we consider the following scenario which often comes up in practical situations and does not have a clear answer in the literature. Does a web guide completely eliminate propagation of lateral position oscillations into downstream spans? To address this issue, we consider the setup shown in Fig. 6 with a controlled remotely pivoted guide at R2; a PI controller is often used in industry based on measurement of web lateral position on the guide roller R2. Figure 9 provides the lateral position response for different locations in the entry span. It is clear that the web guide is able to regulate the web lateral position, $yL(t)$, to zero. This controller is able to regulate lateral position to zero at the guide roller. Further, there is no propagation of these oscillations into the exit span provided the exit span is perpendicular to the entry span, i.e., the angle of wrap on the guide roller is 90 deg. If the wrap angle is 90 deg, the rotation of the guide roller in the plane of the entry span simply twists the web in the exit span. If the wrap angle is not 90 deg, then the in-plane rotation of the guide roller is projected as an initial web slope for the exit span, resulting in twisting and bending of the web in the exit span.

Figure 10 provides the guide roller angle or control action ($\theta L(t)$) to regulate the web lateral position at zero as shown in Fig. 9. When the wrap angle on R2 is $\beta =89deg$, the projection of the guide roller rotation into the exit span plane is $\theta 0R2(t)=\theta L(t)cos(\beta )$ which is provided in Fig. 11; this acts as a disturbance to the exit span. Figure 12 provides the response of the lateral position and slope at different locations of the exit span due to this disturbance. Although the lateral position at the guide roller has been regulated at zero, both the web lateral position and slope at subsequent locations of the span are not zero. There is propagation of lateral oscillations into downstream spans. These typical numerical simulations illustrate the benefit of the developed spatially dependent transfer functions for understanding lateral behavior in ideal as well as nonideal situations.

## Conclusions

In this paper, we have derived spatially dependent transfer functions for web lateral dynamics. The obvious benefits of such governing equations are that one can obtain the evolution of lateral position response at any location in the span as well as all higher-order spatial partial derivatives, such as slope, moment, shear force, etc. Further, these transfer functions may be used to control the lateral position and slope at a prescribed location in the span other than on the roller. In addition, R2R manufacturing of flexible and printed electronics requires positioning the web precisely in both lateral and longitudinal directions. Traditional printing systems have relied solely on longitudinal registration for printing presses with multiple print units. With the goal of achieving print registration accuracy within a few microns in R2R printing of electronics, this work is expected to aid in a more precise analysis of lateral behavior and facilitate the design of model-based lateral control systems for achieving tight regulation of lateral print registration. In the future, we plan to conduct experiments to validate the proposed models and design and evaluate model-based controllers.

## Funding Data

Directorate for Engineering National Science Foundation under Grant No. (1635636).

## Nomenclature

*A*=cross-sectional area of web

*E*=modulus of elasticity of web material

*f*=_{k}defined functions, $k=1,..,3$

*G*=shear modulus

- $gl(x)$ =
defined functions, l = 1, 2, 3

- $glx(x)$ =
first-order spatial derivatives of defined functions

- $glxx(x)$ =
second-order spatial derivatives of defined functions

*I*=web span moment area of inertia

*K*=constant parameter, $K2=T/EI$

*K*=_{e}constant parameter K with shear effect, $Ke2=K2/[1+nT/AG]$

*L*=web span length

- m =
mass per unit length

*n*=correction factor

*N*=shear force

*P*=_{m}numerator polynomial functions in transfer functions, $m=1,..,5$

*s*=Laplace variable

*t*=time

- τ =
time constant, $=L/v$

*T*=web tension

- v =
web transport velocity

*v*=_{y}web lateral velocity

*x*=transport direction distance

- X
_{1}= pivoting distance from guide roller

*y*=web lateral displacement

*z*=_{i}roller lateral displacement

*β*=wrap angle

*θ*=_{i}roller angle

*θ*=_{wi}web angle

- $L{\u2022}$ =
Laplace transform

- $\u2022\u0303$ =
shear variable that can be reduced to bending

- $\u2022\u0302(x,s)$ =
transformed function