PAPER No. 8
Estimating Web Strain, Modulus of Elasticity and Loss Torque Using an Extended Kalman Filter
1999 International Web Handling Conference
Brian Thomas Boulter
Applied Industrial Control Solutions
231 Skyview Drive
Seven Hills, OH, USA 44131
© ApICS ® LLC 2000
Abstract: A set of suboptimal Extended Kalman Filters (EKF's) are described. The filters are used to estimate web tension, modulus of elasticity, and losstorque for all sections in a small web processing line located at Rockwell Automation's Euclid Research Facility. Conditions necessary for observability, and estimate convergence are discussed
I. Introduction
A. Nomenclature
I 
State estimate vector at the k_{th} epoch given the predicted state vector from the previous epoch. 

Noisefree predicted state vector 

Predicted state vector (with residuals) 

Discrete sample time [sec] 

Updated probability density matrix from the previous epoch. 

Linearized Predicted probability density matrix, given the update from the previous epoch. 

Updated probability density matrix 

State Transition Matrix at the k_{th}epoch 

Updated Kalman filter gain matrix 

Sensor noise covariance matrix at the k_{th} epoch 

Process noise covariance matrix at the k_{th} epoch 

Measured states at the k_{th} epoch 

Predicted outputs, given the noisefree predicted states. 
I 
The identity matrix 
e 
Web strain [engineering units of choice] 
t 
Motor torque [engineering units of choice] 
t_{L} 
losstorque [same units as t] 
A_{i} 
i^{th} web cross sectional area [engineering units of choice] 
T_{i} 
i^{th} tenion [units of force] 

i^{th} normalized tension 
E_{i} 
i^{th} modulus of elasticity [units of force per units of A] 

i^{th} normalized modulus of elasticity [1=E_{REF}] 
z_{i} 
i^{th} strain factor 1/(1+ e) 
w_{i} 
i^{th} roll rotational speed [rad/sec] 

i^{th} normalized rotational speed [1=w_{b,i}] 
w_{b,i} 
i^{th} motor base speed [rad/sec] 
L_{i} 
i^{th} length [engineering units of choice] 
R_{i} 
i^{th} roll radius [engineering units of choice] 
GR_{i} 
i^{th }gear ratio 
w_{j,i} 
j^{th} state in the i^{th} zone, zeromean white noise terms. 
B. Problem Definition
The EKF is often used in the aerospace industry to estimate states and/or important system parameters in noisy systems. Why do we need to concern ourselves with this kind of technology in the drive systems business?
There are many sources of uncertainty in a web processing system, for example, tension measurements can be corrupted by incorrect calibration of the load cell or strain gauge. The same can be said for the uncertainty associated with calculation of roll diameter as a function of an estimated surface speed and its ratio to the winder motor speed in ratio detector algorithms. Time varying frictional losses add uncertainty to friction compensation algorithms as do slight changes in material density to inertia compensation algorithms. The estimation of motor torque from a current reference in DC or torque reference in AC as compared to the actual motor shaft torque introduces another source of uncertainty. In the paper and film industries, timevarying modulus of elasticity adds another source of uncertainty. A Kalman filter that is designed to estimate states based on a known model of the uncertainty in the measured states is desirable.
A Luenberger observer, can be implemented [5,6] to estimate system states, but such an implementation is not robust in the presence of noise. This is especially true for metals applications. Because of the high gain in the relationship between stress and strain (Young’s modulus) in metals applications, noise in the speed feedback can result in large errors in the strain estimation. A Kalman filter that is designed to provide reliable estimates of important system parameters in the presence of well modeled system noise would be most useful. An additional benefit provided by such a filter, is that it could estimate unmeasured parameters, such as web strain, and losstorque.
The Kalman filter can provide the user with a measure of how "good" the estimate is. This measure, the covariance of the estimate, is most accurate in systems where the uncertainty is Gaussian, and all process and sensor noise is bandlimited and white. In systems where there are cyclical disturbances (e.g. machine resonances) and the resulting noise is pink at certain frequencies, a Kalman filter can be designed to filter the cyclic noise and calculate unbiased parameter and state estimates [3], however, the design of the Kalman filter in these instances is significantly more complex, for this reason pink noise sources will not be considered in this presentation.
The usual drawback to the implementation of any optimal EKF is the need to heuristically derive optimum noise covariance matrices Q and R. Another is the tradeoff between the need for computing power, and the need for filter bandwidth (i.e. update time). Most implementations usually result in such long update times that the Filter estimates are of little practical use in real time systems.
Figure 7. Shows the optimal EKF implementation for the Euclid labline system. In this system there are 22 system states, and 12 system inputs. To execute the EKF algorithm with this many states on a 7010 Automax ® processor would take approximately 20 [sec] per update. Clearly, with such an update time, the tension estimates could not be used to close a sufficiently responsive tension loop.
However, if a set of suboptimal filters are designed, and linked with estimated parameters, as shown in Figure 8, the resulting suboptimal filter computational requirements are significantly reduced. In this case the most number of estimated states is 3 and the most number of inputs is 2. The resulting matrices are at most 3x3, and all matrix computations can be accomplished with simple routines. The EKF algorithm can now be executed in about 15 [msec] using a 7010 Automax processor. Given that the scantime of a standard Rockwell Automation tension loop is 22 [msec], tension estimates, for example, can be used to close a tension loop.
Having overcome the implementation hurdle, the process of exercising the EKF on the Euclid labline revealed that the estimated web strain, losstorque, and web modulus of elasticity states, in those sections where the tension feedback was not made available to the filter, would not converge. A dual state analysis performed by Angus Andrews at the Rockwell Science Center [8] provided a proof that tension feedback was required if these states were to be estimated. The conclusions drawn in the proof were extensible to Luenberger observers, implying that any observer or state estimator of web strain, losstorque, or modulus of elasticity was not feasable if it is implemented without tension feedback. This was an important conclusion given that the impetus for the study was to investigate the possibility of developing tension sensorless web process lines.
C. Paper Overview
The first section contains a brief description of the EKF algorithm. The second section describes the nonlinear process model used in the final suboptimal Kalman filter implemented on the Euclid labline. The third section describes the coding of the EKF for use with the Automax distributed control system, and results from the integration of the filter onto the Euclid labline. Conditions necessary for observability and estimate convergence, along with conclusions are presented in the fourth section.
II. The Extended Kalman FIlter Algorithm
The EKF algorithm described by Figures 4 and 5, and derived in numerous texts on the subject ([1,2,3,4]) is summarized below. The noise statistics are modeled as follows: Let Q and R be the covariances of zeromean Gaussian process noise ‘w’ and the measurement noise ‘v’ respectively: Then:
State prediction equations:
Covariance prediction equations:
Kalman gain equations for each suboptimal section
Covariance update equations:
State update equations:
III. System Modeling and Filter Implementation
The Euclid labline (Figure 1.) is a four section web processing line designed for the evaluation of web handling control strategies. For this study the web material used in the process was Kraft paper. A oneline diagram of the labline is shown in Figure 2. A simplified oneline diagram is shown in Figure 3.
Figure 1. Euclid LabLine Facility
Figure 2. LabLine Physical Model
Figure 3. Simplified LabLine Model
A. State Variable Designations
The EKF algorithm described above (Figures 4, and 5), requires the apriori designation of all system state variables, inputs, outputs, measured state variables, estimated state variables, process noise sources, and sensor noise sources. These designations are made based on a satisfactorily complete model of the process.
Figure 4. Algorithm Structure
Figure 5. Extended Kalman Filter Algorithm
Unwinder + Downstream Web (With Dancer Disabled as a Torque Regulator)
Measured states:
Roll (angular) speed
Roll diameter (from ratio detector)
Roll torque: t = f(current & field strength)
Estimated states:
Web strain factor: z = 1/(1+e ), where e = strain
Roll losstorque: t _{L}
Web modulus of elasticity: E
Tensions:
Input (woundon): assumed zero
Output (downstream): calculate, T = AEe = AE(1/z –1)
SWrap Bridle + Downstream Web
Measured states
Roll (angular) speed
Roll torque: t = f(current & field strength)
Estimated states
Web strain factor: z = 1/(1+e ), where e = strain
Roll losstorque: t _{L}
Web modulus of elasticity: E
Tensions
Input: from upstream subsystem
Output (downstream): calculate, T = AEe = AE(1/z –1)
PullRoll + Downstream Web
Measured states:
Roll (angular) speed
Roll torque: t = f(current & field strength)
Web tension (from load cell)
Estimated states:
Web strain factor: z = 1/(1+e ), where e = strain
Roll losstorque: t _{L}
Web modulus of elasticity: E
Tensions:
Input: from upstream subsystem
Output (downstream): calculate, T = AEe = AE(1/z –1)
Rewinder
Measured states:
Roll (angular) speed
Roll diameter (from ratio detector)
Roll torque: t = f(current & field strength)
Estimated state:
Roll losstorque: t _{L}
Tensions:
Input: from upstream subsystem
Output (downstream, woundin): assumed zero
B. First Principal Equations
The following first principal equations were used to describe each section of the simplified system (Figures 3 and 6) in each of the suboptimal EKF’s (Figure 8). Note: for the sake of brevity, the Unwinder was configured as a torque regulator, and the estimator for the dancer position was not implemented. Note: the choice of units must be made appropriately so that the equation dimensions are correct.
State Equations
Upstream: 

Downstream: 

Motor speed: 

Motor torque: 

Losstorque: 

Strain factor: 

Modulus of elasticity: 

Where: 

Note: The w_{i} terms represent the zeromean noise terms (covariance = Q_{i} ).
The subsystem model equations for the Unwinder and Rewinder are slightly different than those shown above for constant diameter rolls. First, an additional state equation is necessary in both the Unwinder and Rewinder to estimate the varying roll diameter; a measurement of the diameter is available in each case (from a ratio detector). And second, the Rewinder does not have a downstream web, so the state equation for strain factor is eliminated. The derivation of these equations is left up to the reader.
C. LabLine Implementation
Figure 6. Subsystem Model Structure
Figure 7. Optimal Kalman Filter Design
Figure 8. SubOptimal Kalman Filter Design
Each subsystem model was implemented as one roll plus downstream web, as illustrated in Figure 6. The subscript 0 represents upstream and the subscript 1 represents downstream. The implementation was structured with one AutoMax 7010 processor card per suboptimal filter. Programming was done using AutoMax enhanced BASIC equations. Software switches were included for correctly reassigning variables when reversing line speed direction. The various measured state variables from each subsystem were observed using the Sigma ® data system during line operation in the forward and reverse directions. Data was extracted and plotted using the Mathworks Corp. MATLAB ® data presentation algorithms.
IV. Observations
Estimated states for the line between the Pullroll and the Rewinder were the most wellbehaved. The estimated modulus of elasticity for the web section between the two rolls tended to settle at a steadystate value that was reasonably close to the expected value. The reason for these wellbehaved results was later proven to be due to the presence of a tension measurement in this zone.
Estimated states for other regions of the line tended to deviate from expected values, similar to the deviations noted during MATLAB simulations of the Labline. In particular, the modulus of elasticity values for the two other web sections – i.e., the Swrap to Pullroll section and the Unwinder to Swrap section, appeared to diverge either above or below their expected values, depending on the covariance values selected for the Qmatrix and Rmatrix elements. This indicated that the modulus of elasticity and losstorque are not observable, unless there is a measurement of tension available.
The results presented below for the SWrap section losstorque estimation illustrate these points. The measured states are angular speed, w _{0} , and motor torque, t _{0}. The estimated states also included: strain factor, z_{1} , and a calculated modulus of elasticity, E_{1}. For the diagonal Q matrix covariances given in Table 1, the results are shown in Figures 9 through 11:
Table 1.
Figure 9. Actual tension (green) versus Kalman filter estimated tension (yellow)
Figure 10. Calculated modulus of elasticity (green) versus Kalman filter estimate (yellow)
Figure 11. Losstorque estimated by Kalman filter
Figure 12. Estimated losstorque for cross coupling terms Q(1,3) = Q(3,1) = (0.006)^2.
Figure 13. Estimated losstorque for cross coupling terms Q(1,3) = Q(3,1) = (0.01)^2
Increasing the process covariance for the modulus of elasticity from (.003)^{2} to (.01)^{2} caused the estimated value of E_{1} to diverge above the actual. Decreasing the process covariance for E_{1} from (.003)^{2} to (.001)^{2} caused the estimated modulus of elasticity to diverge below the actual value. Similar heuristic strategies were employed for optimally setting the covariances of all other state variables. The effect of offdiagonal values in the Q matrix were then investigated.
The most likely offdiagonal terms to be present in the above example were assumed to be the ones linking losstorque and roll angular speed; i.e., most likely, t _{L} » B_{v} w . If roll angular speed is state 1 and losstorque is state 3, then the crosscoupled process covariance becomes the matrix elements Q(1,3) = Q(3,1). Estimation of losstorque appeared to approach a steadystate value only when these crosscoupling terms were at or near a critical value. Much effort was expended in investigating and identifying optimal settings for all diagonal and crosscoupling covariances. For example, decreasing the crosscoupled process covariance, losstorque to speed, Q(1,3) = Q(3,1) from (.008)^{2} to (.006)^{2} resulted in the losstorque estimate shown in Figure 12. Whereas increasing the crosscoupled process covariance, losstorque to speed, Q(1,3) = Q(3,1) from (.008)^{2} to (.01)^{2} resulted in the unbounded losstorque estimate shown in Figure 13. Similar observations were made while setting other estimated state process and sensor noise covariance values.
It is believed that the estimation of Young’s modulus would be significantly improved if a model relating stress to strain for the material in the web process line was implemented in the above state equations. The model used in the Kalman filter on the labline was composed of only a noise model. Several models were investigated, they were obtained heuristically with measurements of web stretch as a function of tension, however, slight changes in humidity rendered them inaccurate. The stress strain relationship of metals appears to be much better understood and modeled. It is believed that the Kalman filter will perform significantly better in metal processes where these models can be implemented. The author noticed that, in the area of paper, film, and plastics, there is a lack of published material on this subject.
V. Conclusions
The author corresponded with Angus Andrews at the Rockwell Science Center for assistance when it was observed that the state estimates of tension, strain factor, web modulus of elasticity, and losstorque would not converge in those sections where a tension measurement was not fed back into the EKF. A dual state analysis of the system was performed [8] and a comparison of optimal and suboptimal filters revealed that the above states were only observable if a measurement of at least one of the above mentioned states was made available to each the suboptimal filters. Tension is the only state that can be practically measured. This conclusion implied that the implementation of a tension sensorless system is not feasible. However, important states, such as web strain, or losstorque, can be estimated if the tension feedback is made available, and an EKF is used to estimate these states, even in the presence of significant process and sensor noise.
The structure and results of the dual state analysis will be the subject of another paper on this topic. However, by inspection of the equations describing strain factor, tension, losstorque (i.e. speed), it is clear that the above states exhibit dependence. For example, if the estimate of web tension is lower than the actual, the estimate of losstorque must be higher than the actual, this is because the estimate of shaft torque is accompanied by a measurement. It can be considered known with a degree of certainty. At a steady state speed, the torque on the shaft minus losstorque, will produce strain in the material. What is lost in one estimate must be made up in the other estimate.
VI. Acknowledgements
The author would like to express his gratitude to the following people who, provided assistance in the completion of this study. Dr. Alan S. Wilsky, and Dr. Robert R. Tenney of the Massachusetts Institute of Technology for their "Fundamentals of Detection, Parameter Estimation, and Kalman Filtering" summer course [2]. Dr. Angus Andrews of the Rockwell Science Center for his insights into Kalman Filtering and dual state analysis techniques [7], [8]. Harold W. Fox of ApICS LLC for his help in analyzing the theory, his work in coding the nonoptimal Kalman Filters into Automax Extended Basic, and his help in undertaking the testing. We would also like to thank Rockwell Automation for their permission to use the Euclid Lab Line for testing the algorithms.
References
[1] Kalman, R. E., "A New Approach to Linear Filtering and Prediction Problems", Transactions of the ASME, J. Basic Engineering, Vol. 82 pp 3445, 1960.
[2] Willsky, Alan S., Tenney, Robert R., "Fundamentals of Detection, Parameter Estimation and Kalman Filtering", M.I.T. Professional Institute Course Materials, 1997.
[3] Gardner, William A., "Introduction to Random Processes With Applications to Signals and Systems", McGraw Hill ISBN 0070228558 1989.
[4] Kailith, T., "Linear Systems", Englewood Cliffs, N.J.: PrenticeHall. 1980.
[5] Wolferman, W., "Tension Control of Webs  A Review of the Problems and Solutions in the Present and Future", Tab 15, Proceedings of the 3^{rd} IWHC International Web Handling Conference, Oklahoma State University June 1995.
[6] Wolferman, W., "Sensorless Tension Control of Webs", Tab 23, Proceedings of the 4^{th} IWHC International Web Handling Conference, Oklahoma State University June 1997.
[7] Andrews, Angus, "Appendices on Reducing Tension Sensor Requirements in Web Processing", Unpublished, Internal Rockwell Automation engineering document, 1997.
[8] Andrews, Angus, "Dual State Analysis of Web Processing Line Models", Unpublished, Internal Rockwell Automation engineering document, 1998.
