A governing equation for web tension in a span considering thermal and viscoelastic effects is developed in this paper. The governing equation includes thermal strain induced by web temperature change and assumes viscoelastic material behavior. A closed-form expression for temperature distribution in the moving web is derived, which is utilized to obtain thermal strain. A model for web tension in a multispan roll-to-roll system can be developed using this governing equation. To evaluate the governing equation, measured data from an industrial web process line are compared with data from model simulations. Since the viscoelastic behavior of web materials is affected by the web temperature change, elevated temperature creep, and stress-relaxation experiments are conducted to determine the temperature-dependent viscoelastic parameters of the utilized viscoelastic model. Comparisons of the measured data with model simulation data are presented and discussed. An analysis of the web tension disturbance propagation behavior is also provided to compare transport behavior of elastic and viscoelastic materials.

## Introduction

The increased use of continuous processing of flexible materials by transporting them on rollers has resulted in considerable focus on investigating the encountered technical problems and providing suitable solutions to maintain productivity and good quality of products. The term web is used to refer to thin materials, which are manufactured and processed in continuous, flexible strip form, such as plastics, textiles, papers, and composites. Web tension and velocity are the two key transport variables that must be maintained at reference values to ensure the quality of web products. Tension variations in a web due to web treatment processes (such as heating, cooling, and printing) are encountered in many web handling industries; these processes create disturbances to web tension control.

Temperature in the moving web is one of the most important factors affecting the transport behavior as it changes both mechanical and physical properties [1]. In various web processes such as printing, coating, or other treatments, heat transfer between the web and environment or between the web and the contact area of a heated or cooled roller is encountered. Web tension is affected by the temperature change in the web span as it introduces thermal strain in the web.

Web tension behavior in the transport direction was investigated in Ref. [2]. Thermal effects from a heated roller and ambient air were included in web tension models. A dynamic model for web strain in a span that includes thermal and hygral effects was developed in Ref. [3]. A lumped capacitance model is considered to determine the heat transfer from the web surface in the free web span. A dynamic model for web tension behavior in a furnace of a continuous annealing line was developed in Ref. [4]. A new control scheme was proposed that considers thermal strain and the change in the Young's modulus. Research in utilizing models to design and develop roll-to-roll machines and design of model-based control systems has been active in the last several decades [5–9]. The model derived in this work can be used to design advanced model-based controllers to efficiently control web transport behavior.

Many web materials exhibit viscoelastic behavior when undergoing deformation during transport, especially in the presence of heating and cooling sources. In such situations, one needs to consider viscoelastic material properties in the derivation of the tension governing equation. Flexible and composite materials exhibit viscoelastic behavior that is directly attributed to the molecular structure of the underlying materials [10]. Viscoelastic behavior of web materials is considered in several papers [11–13]. The existing work did not consider methods to determine governing equations for web tension behavior that includes the combined effects of both temperature distribution in the moving web and the viscoelastic material behavior.

In this paper, thermal and viscoelastic effects are included in the developed web tension governing equation. The induced thermal strain in a web span due to temperature change is determined by considering heat transfer in the region of wrap due to conduction from the roller and in the free web span due to convection to ambient air. Temperature-dependent viscoelastic behavior is investigated using the standard linear solid (SLS) model. The viscoelastic parameters of the utilized viscoelastic model are determined by conducting elevated creep and stress-relaxation experiments on several web composite materials. A section of an industrial web process line containing multiple spans is considered for validation of the proposed model. Model simulations are conducted with an existing control strategy implemented on the production web process line. The contributions of the paper lie in simultaneously considering the thermal and viscoelastic effects on web tension, determining the properties of viscoelastic model parameters, conducting model simulations by considering a section of an industrial web processing line, and comparison of the model simulation results with measurements obtained from an industrial processing line to validate the proposed model.

The remainder of the paper is organized as follows: Web tension governing equation is developed in Sec. 2 by considering thermal strain and a viscoelastic constitutive relation between mechanical strain and web tension. Viscoelastic model parameter identification for some composite materials transported in the industrial web processing line, model simulations, and comparison of simulation results with experimental results are described in Sec. 3. Disturbance propagation in a multispan web system is discussed in Sec. 4. Conclusions and potential topics for future work are provided in Sec. 5.

## Governing Equation for Web Tension Considering Thermal and Viscoelastic Effects

To develop a governing equation for web tension by considering both thermal and viscoelastic effects, the following procedure is employed. First, a governing equation for longitudinal strain in the web is developed by utilizing the law of conservation of mass in a control volume encompassing a span (web between two adjacent rollers); the longitudinal strain is a combination of both mechanical strain and thermal strain. Second, to determine thermal strain due to heating/cooling of the web, one requires the temperature distribution in the moving web. Heat transfer due to conduction is assumed in the region of web wrap on the roller and due to convection in the free span between rollers to determine the temperature distribution. Third, to determine a constitutive law that relates web mechanical strain and tension, the web is assumed to be linearly viscoelastic. These three aspects are combined to develop a governing equation for web tension.

### Governing Equation for Web Strain.

where the following notation is used: *x,* transport direction coordinate, *t*: time, $\rho (x,t)$, material density; *A,* web cross section area; *L _{i}*, length of the

*i*th web span;

*ε*, web strain in the

_{i}*i*th web span;

*t*web tension in the

_{i,}*i*th web span;

*v*web transport velocity on the

_{i,}*i*th roller;

*ω*angular velocity of the

_{i,}*i*th roller. Note that

*t*without a subscript is used for time and with an

*i*is used to denote web tension in the

*i*th span. Also, note that small strain assumption is utilized in subsequent derivations, which is true for many flexible materials (such as metals, polymers, and paper) processed using roll-to-roll methods.

*x*in

*ε*is used to indicate that this is longitudinal (or transport) direction strain and depends on the transport direction coordinate

_{x}*x*. Substitution of Eq. (4) into Eq. (1) for the

*i*th span results in

*i*th span is a function of the web velocities on the adjacent rollers and the strains in the

*i*th and $(i\u22121)$th spans. The total strain in the

*i*th span is assumed to be given by strain due to mechanical stress and strain due to temperature distribution in the moving web as follows:

where *α* is the coefficient of thermal expansion, $\u03d1i(x,t)$ is the temperature distribution in the moving web, and $\u03d10$ is the reference temperature. In the following, a closed-form expression for temperature distribution in a moving web will be derived.

### Temperature Distribution in the Moving Web.

where $K\u03d1$ is the thermal conductivity constant and *y* is the transverse direction coordinate (along web thickness). The web is treated as a slender rectangular bar within the region of wrap. Since the wrap region is small compared to the length of the web span, the temperature along the transport direction within the wrap region is assumed to be the same. If $\varphi $ is the angle of wrap, *R* is the roller radius, and $v\xaf$ is the average web transport velocity on the roller, then the resident time of a representative volume element is given by $R\varphi /v\xaf$; in general, this resident time is small because the average transport speed is much larger than the product of roller radius and wrap angle.

where *h* is the web thickness, $\theta w,i\u22121$ is the temperature of the web prior to engaging the (*i* − 1)th roller, $\theta R,i\u22121$ is the temperature of the roller, and $\theta U,i\u22121$ is the ambient temperature in the region of wrap. Note that ϑ is used to denote temperature that is a function on both the spatial coordinate and time, whereas *θ* is used for a boundary condition or a temperature of the web for a fixed spatial location.

*ρ*is the web density,

*c*is the specific heat of the web material, and $\gamma *$ is the coefficient of heat transmission, which represents the heat flux for a unit difference in temperature between the web surface and the surrounding ambient air. To simplify the solution of Eq. (15), it is assumed that the web transport velocity is equal to the average web velocity, i.e., $v(x,t)=vr$ and the ambient temperature $\u03d1U,i$ is a function of time only. The temperature in the web with these assumptions is given by

### Mechanical Strain Governing Equation.

*v*, and reference strain,

_{r}*ε*, by setting $vi(t)=vr+Vi(t)$ and $\epsilon t,i(t)=\epsilon r+Si(t)$ (where $Vi(t)$ and $Si(t)$ are the velocity and strain variations, respectively), and ignoring the second-order terms, results in the following governing equation for variation in strain:

_{r}### Linear Viscoelastic Model Relating Mechanical Strain and Web Tension.

*b*denotes the damping coefficient of the damper; the subscript

_{i}*i*denotes that these parameters correspond to the

*i*th span. The sum $E1,i+E2,i$ is called the initial modulus of the SLS model and $E1,i$ is called the equilibrium modulus. By taking Laplace transform of Eq. (20), the following relationship between the stress and strain in the frequency domain can be obtained:

*t*as $Ti(t)=ti(t)\u2212tr$. Let $Ti(s)$ be the Laplace transform of $Ti(t)$. Then, since the stress is given by $\sigma i(s)=Ti(s)/A$, Eq. (21) can be written as

_{r}Equation (23) provides the tension behavior of a moving web in the *i*th span considering the effect of temperature change and viscoelastic material behavior. The following key assumptions were utilized to derive this governing equation for web span tension: (1) strain is assumed to be small, (2) only heat conduction in the thickness direction is assumed between the web and the roller surface in the region of wrap, (3) only convection between the web surface and ambient air is assumed in the free span, and (4) the web material is assumed to be linearly viscoelastic.

It is evident that there is transport of tension from upstream spans ($Ti\u22121$) to the downstream spans. The first term of the equation includes the viscoelastic effects of the upstream and downstream spans. Web tension is affected by temperature difference between adjacent web spans ($\u03d1i\u22121\u2212\u03d1i$), and the temperature of the web at the beginning of the web span ($\u03d1i(0,t)$); these are obtained from Eqs. (13), (14), and (16). This span governing equation allows (1) analysis of web transported through heating/cooling sections and (2) the control system designed to test the robustness of controllers applied to web handling processes by addressing viscoelastic and thermal effects on web tension.

## Model Validation

To evaluate the developed model, experiments and model simulations were conducted for a section of an industrial web process line used for manufacturing flooring materials in web form. To conduct model simulations, one requires temperature-dependent viscoelastic parameters of the SLS model. In the following, we first provide a discussion on the stress-relaxation and creep experiments that were conducted to determine the SLS model parameters, followed by a description of the experimental process line and a discussion of the model validation results.

### Determination of Viscoelastic Model Parameters.

The viscoelastic parameters of the SLS model were determined by conducting elevated temperature creep and stress-relaxation experiments. These experiments were conducted on a tensile testing machine equipped with an environmental chamber to allow for testing at elevated temperatures; test specimens were prepared according to the ASTM D638-08 [18]. Experiments were conducted for two composite web materials, Royelle and Stratamax Better, that are manufactured by Armstrong World Industries; we will refer to these materials as web 1 (Royelle) and web 2 (Stratamax Better). Both materials are made of felt and gel (with different felt combinations). Felt is made of cellulose, limestone, and styrene-butadiene-rubber latex, and gel is made of polyvinyl chloride, di-iso phenyl phthalate, limestone, titanium dioxide, and other process aids. Experiments were conducted at different levels of temperature (i.e., 18 °C, 66 °C, 110 °C, and 154 °C). These temperature levels were chosen as they are close to the operating temperatures in the industrial web processing line. Since the procedure is the same, we will show the experimental findings for only one web 1.

Stress-relaxation tests were carried out by deforming the specimens to a given amount and the decrease in stress was recorded for a duration of 1000 s. Specimens were deformed gradually until they reached the desired strain of 1%. After that, strain was maintained constant for 1000 s and the decrease in stress with respect to time was recorded. Data from these tests are provided in Fig. 4.

*E*) for the web materials and the linear spring constant (

*E*

_{1}). The initial modulus of elasticity is the sum of the moduli of the two linear springs in the SLS model, which can be determined from the relaxation data as the modulus at the initial time using the following equation:

The modulus *E*_{1} represents the steady-state response of the stress-relaxation function as shown in Eq. (24). Therefore, *E*_{1} was determined from the stress-relaxation curves after large time. *E*_{2} was determined as the difference between the initial modulus *E* and *E*_{1}.

*b*. The creep function is given by

where $\tau c=\tau r((E1+E2)/E1)$ is the creep or retardation time and *τ _{r}* is the relaxation time. The viscosity parameter

*b*was evaluated from the regression analysis of the creep data using Eq. (25).

Figure 5 shows the temperature-dependent modulus for web 1 and web 2. It can be seen that the Young's modulus decreases with increasing temperature for both web materials. The cross marks in the graph represent the experimental data, whereas the solid lines are for curve fitting. The temperature-dependent moduli determined are $362\u2009exp(\u22120.0025\theta )$ for web 1 and $463\u2009exp(\u22120.0029\theta )$ for web 2.

Figure 6 shows the data corresponding to the temperature-dependent modulus *E*_{1} of the SLS model, which can be computed at a given temperature using: ($312\u22121.633\u2009\theta 0.854$) for web 1 and ($368\u22120.0064\u2009\theta 1.97$) for web 2. The modulus *E*_{2} is obtained by subtracting *E*_{1} from *E* for each web material.

The viscosity parameter *b* was determined using the creep data and based on the methods provided in Refs. [19] and [20]. The resulting expressions are: $7.58(\theta r/\theta )0.34$ for web 1 and $5.27(\theta r/\theta )0.37$ for web 2, where *θ _{r}* and

*θ*are the room temperature and operating temperature, respectively.

### Web Process Line Description and Model Validation Results.

A section of the coating and fusion line at Armstrong World Industries is considered for model simulations; a line schematic of the section is provided in Fig. 7. In this section, the web is coated by a clear coat and transported through heating and cooling rolls to modify web properties. Heating the web melts and cures the coating to get a uniform coat thickness on the surface of the web material. The cooling process is necessary after heating to relieve internal stresses that are induced by the heating process. This section was selected as it involves thermal process as the web is transported over heated/chilled rolls. Measured data from this section are compared with data from the model simulations. In model simulations, we employed the same control strategies that are used in the industrial process line. Simulations were conducted for the two web flooring products, Royelle (web 1) and Stratamax Better (web 2).

It is common practice in web handling to divide the web process line into tension zones, where a tension zone is defined as the web section between two independent driven rollers [21]. Since a typical web process line contains a significantly more number of idle rollers when compared to the number of driven rollers, it considerably simplifies both the analysis of the web process line and the design of control systems. This simplification means that the dynamics of the idle rollers are treated as perturbations or disturbances to the tension zone control systems. In the case of Fig. 7, there is a tension zone from the coating roller to the heating roller and a tension zone from the heating roller to the set of four cooling rollers. Note that the cooling rollers are either mechanically slaved with a single motor shaft driving the four roller shafts using a transmission or electronically slaved with a single controller driving all four motors, and so these are treated as a single drive unit. With the definition of the tension zones, the coating and fusion line web section shown in Fig. 7 can be simplified to the figure shown in Fig. 8. In Fig. 8, *M*_{0} is the coating roller, *M*_{1} is the heating roller, and *M*_{2} denotes the integration of the four cooling rollers into one unit. We will refer to the tension zone between *M*_{0} and *M*_{1} as zone 1 and the one between *M*_{1} and *M*_{2} as zone 2. Note that the span tension equation is applied to the length of the web in each tension zone; in this regard, one can consider the tension in each of these zones as an average value. We will denote tension in zone 1 as *t*_{1} and in zone 2 as *t*_{2}. It is noted that both the span lengths between the driven rolls and the wrap angle of the web around the rolls will be maintained even though this is not apparent from the simplified line sketch shown in Fig. 8.

The transport velocity of the web on the *i*th roller is denoted by $\u2009vi$ and the input torque from the *i*th motor is denoted by $\u2009ui$. Coating, heating, and cooling rollers are all driven and coating and heating rollers are under pure speed control, whereas cooling roller is under speed and tension control; the speed-based tension control strategy that is utilized is provided in Fig. 9.

*u*is the input torque applied on the roller

_{i}*i*,

*n*is the gear ratio of the motor to roller

_{i}*i*, and

*b*is the coefficient of the friction in the roller shaft bearings. The web velocity is related to the roller angular velocity by the relation, $vi=Ri\omega i$, assuming there is no slip between the web and the roller. The governing equation for the dynamics of web velocity variation is given by

_{fi}*T*and

_{i}*U*are the variations in tension and control input, respectively,

_{i}*u*is the control input that maintains the forced equilibrium at the reference values [21], and

_{oi}*t*is the reference tension. Thus, the dynamics for web velocity variations in the time domain can be written as

_{r}where $\tau vi=Ji/bfi$ and $Ui=ui\u2212uoi$. Equation (31) represents the deviation (variation) in velocity from its reference value ($Vi=vi\u2212vr$, where *v _{i}* is the actual velocity and

*v*is the reference velocity).

_{r}Model simulations were conducted with the same controllers (proportional–integral controllers) and gains that were used in the industrial process line; proportional and integral gain values of 6.3 and 1 for the velocity controller, and 7.8 and 0.1 for the tension controller are used. The velocity reference used in the simulations is 0.5 m/s for both web 1 and web 2 to match that of the industrial web processing line. The web tension reference in the industrial web processing line was 978 N (220 lb) for both web 1 and web 2, which was also used in the model simulations. It is assumed that the reference velocity is identical for all tension zones in the simulations. The other web line parameters used in the model simulations are shown in Table 1. The web material parameters used in the model simulations are provided in Table 2. The average web temperature in span *L*_{1} is 26 °C, whereas the average temperature in web span *L*_{2} is 125 °C. These temperatures were measured in the industrial web processing line using a laser temperature sensing device. This average temperature was used only for the determination of the temperature-dependent viscoelastic parameters discussed in Sec. 3.1. The model simulations use Eqs. (27) and (31) for each span of the given web section, and as such consider the evolution of temperature distribution as given in Sec. 2.2. An idle roller mounted on load cells is used to measure web tension in web span *L*_{2,} which is fed back to the tension controller that gives a speed correction to the inner speed loop. The inner speed-loop controls the speed of the cooling rollers (see Fig. 9).

Data from model simulations were compared with measured data from the industrial process line and are shown in Fig. 10 for web 1. Figure 11 shows the fast Fourier transform for the model tension output results and the measured data for web 1. The frequency content of both data is also correlated well. Similar trends are obtained for web 2.

Tension zone 1 between the coating and heating rollers is not equipped with a load cell roller for measuring tension in the industrial process line. Thus, web tension was determined using only model simulations and the data are shown in Fig. 12 for web 1. This tension zone does not have any cooling/heating areas as the coating process is implemented at room temperature, but its behavior is affected by its adjacent web spans and transport of disturbances from the upstream spans. Similar trends were obtained for web 2. In Sec. 4, we provide a discussion on how to study propagation of tension disturbances in a multispan system using the governing equations derived in this paper.

## Disturbance Propagation in Multispan System

*V*

_{0}and

*T*

_{0}, are considered as disturbances to the multispan system. In the figure,

*T*

_{1},

*T*

_{2}, and

*T*

_{3}are the tension variations in web spans

*L*

_{1},

*L*

_{2}, and

*L*

_{3}, respectively, and

*V*

_{0},

*V*

_{1},

*V*

_{2}, and

*V*

_{3}are the velocity variations of the rollers

*R*,

_{o}*R*

_{1},

*R*

_{2}, and

*R*

_{3}, respectively. Define the following variables: $Kve,i=AE1,i/vr,\u2009\tau c,i=bi(E1,i+E2,i)/(E1,iE2,i),\u2009\tau r,i=bi/E2,i,\u2009\tau i=Li/vr$ and $Kth=\alpha vr$,

*i*= 1, 2, and 3. Also, define the following transfer functions:

*n*spans provides the following generalized transfer function from

*V*

_{0}to tension variation in the

*n*th zone,

*T*as:

_{n}We provide a brief discussion on these disturbance propagation effects by considering a web section with three tension zones as shown in Fig. 13, where with web lengths in the three tension zones, *L*_{1}, *L*_{2}, and *L*_{3}, are at 26 °C, 70 °C, and 125 °C, respectively. The temperature of the upstream span to the roller *R*_{0} is considered the same as the downstream roller (tension zone *L*_{1}). The viscoelastic parameters of the SLS model at these temperatures are determined using the equations derived in Sec. 3. For this study, we considered all web span lengths to be equal to 3 m and the web reference velocity to be 0.5 m/s.

Figure 14 shows the bode magnitude plot for tensions in three successive spans (with lengths *L*_{1}, *L*_{2}, and *L*_{3}) in response to the upstream tension disturbance *T*_{0}, when considering web 1 to be viscoelastic. It can be observed that the disturbance *T*_{0} affects web tension dynamics at low frequencies as the temperature increases in successive spans in the transport direction. The viscoelastic parameters $E1,i,\u2009E2,i$, and *b _{i}* decrease with increase in web temperature. As a result, the creep time $\tau c,i$ and the relaxation time $\tau r,i$ are decreased. For the viscoelastic case, the magnitude of the transfer function $Gtn(s)$ decreases with increase in web temperature in the successive spans in the transport direction. This behavior can be attributed to the fact that the viscoelastic parameters of web materials decrease with increasing temperature. The magnitude of these transfer functions is proportional to the temperature difference in the adjacent web spans. At frequencies greater than $1/\tau i$, attenuation of the disturbances

*V*

_{0}and

*T*

_{0}is observed as the web travels into subsequent spans.

Figure 15 shows the effect of propagation of disturbance *T*_{0} on web tensions, $T1,\u2009\u2009T2,\u2009\u2009and\u2009T3$, with increasing temperature in the same three tension zone section by considering elastic behavior for the web 1 material. The evolution of web tension the elastic web was obtained by setting the viscoelastic parameter *b* = 0; this means that $\tau ci=0$ and $\tau ri=0$. It can be seen that the magnitude of the transfer functions from the input *T*_{0} to *T _{n}* is constant at low frequencies, although there is a temperature change in successive web spans. A comparison of the elastic and viscoelastic cases indicates that the low frequency behavior is different for the parameters considered. The developed governing equations will assist in such disturbance propagation studies with various viscoelastic material parameters. We have also studied the effect of the upstream velocity disturbance

*V*

_{0}via the transfer function $Gvn(s)$ on downstream span tensions; similar disturbance propagation magnitudes were observed for both viscoelastic and elastic cases.

## Conclusions

In this paper, we have derived a governing equation for web span tension by considering heating/cooling of moving webs and viscoelastic behavior for the web material. A simple heat transfer model that provides a closed-form equation for temperature distribution and a first-order viscoelastic material model are employed in the derivation of the web span tension governing equation. Although model simulations and experimental results were provided for specific flooring materials that are processed and manufactured in the industrial web processing line, the models are applicable to other web materials if the viscoelastic parameters of those webs are determined by appropriate stress-relaxation and creep tests. The advantage of the developed governing equation for web tension lies in the development of advanced model-based tension control systems to improve transport behavior. Development of advanced tension control systems based on this tension equation will be considered as part of a future study.

## Funding Data

U.S. National Science Foundation (Grant Nos. 0854612 and 1246854).

## Nomenclature

*A*=cross-sectional area of web

*b*=viscosity of the web

*b*=_{fi}coefficient of friction of roller

*i**c*=specific heat of the web material

*E*=modulus of elasticity (Young's modulus)

*E*_{1},*E*_{2}=springs moduli of the standard linear solid model

*h*=web thickness

*J*=_{i}inertia of roller

*i*- $K\u03d1$ =
thermal conductivity constant

*L*=_{i}length of span

*i**m*=mass of the web

*n*=_{i}gear ratio of the motor to roller

*i*for driven rollers*R*=_{i}radius of roller

*i**S*=_{i}strain variation

*t*=time

*t*=_{i}web tension in the span

*i**T*=_{i}web tension variation in the span

- $Tw,i$ =
span time constant

- $T\u03d1$ =
thermal time constant

*u*=_{i}control input to roller

*i**v*=_{i}surface velocity of roller

*V*=_{i}velocity variation in the roller

*α*=coefficient of thermal expansion

- $\gamma *$ =
coefficient of heat transmission

*δ*=exponential function of the ratio of time constants

*ε*=_{i}strain in span

*i*- $\epsilon \u02d9i$ =
strain rate in span

*i* *θ*=temperature, function of

*t*only- ϑ =
temperature, function of

*x*and*t* - $v\xaf$ =
average web velocity

*ρ*=web density

*σ*=applied stress

- $\sigma \u02d9$ =
stress rate

*τ*=_{c}creep time

*τ*=_{i}time constant

*τ*=_{r}relaxation time

- $\varphi $ =
angle of wrap

*ω*=roller angular velocity