Thickness Control for a Steel Beam

This example shows how to design a MIMO LQG regulator to control the horizontal and vertical thickness of a steel beam in a hot steel rolling mill.

Rolling Stand Model

Figures 1 and 2 depict the process of shaping a beam of hot steel by compressing it with rolling cylinders.

Figure 1: Beam Shaping by Rolling Cylinders.

Figure 2: Rolling Mill Stand.

The desired H shape is impressed by two pairs of rolling cylinders (one per axis) positioned by hydraulic actuators. The gap between the two cylinders is called the roll gap. The goal is to maintain the x and y thickness within specified tolerances. Thickness variations arise primarily from variations in thickness and hardness of the incoming beam (input disturbance) and eccentricities of the rolling cylinders.

An open-loop model for the x or y axes is shown in Figure 3. The eccentricity disturbance is modeled as white noise w_e driving a band-pass filter Fe. The input thickness disturbance is modeled as white noise w_i driving a low-pass filter Fi. Feedback control is necessary to counter such disturbances. Because the roll gap delta cannot be measured close to the stand, the rolling force f is used for feedback.

Figure 3: Open-Loop Model.

Building the Open-Loop Model

Empirical models for the filters Fe and Fi for the x axis are

and the actuator and gap-to-force gain are modeled as

To construct the open-loop model in Figure 3, start by specifying each block:

Next construct the transfer function from u,we,wi to f1,f2 using concatenation and append as follows. To improve numerical accuracy, switch to the state-space representation before you connect models:

T = append([ss(Hx) Fex],Fix);

Finally, apply the transformation mapping f1,f2 to delta,f:

Px = [-gx gx;1 1] * T;
Px.OutputName = {'x-gap' , 'x-force'};

Plot the frequency response magnitude from the normalized disturbances w_e and w_i to the outputs:

bodemag(Px(: , [2 3]),{1e-2,1e2}), grid

Note the peak at 6 rad/sec corresponding to the (periodic) eccentricity disturbance.

LQG Regulator Design for the X Axis

First design an LQG regulator to attenuate the thickness variations due to the eccentricity and input thickness disturbances w_e and w_i. LQG regulators generate actuator commands u = -K x_e where x_e is an estimate of the plant states. This estimate is derived from available measurements of the rolling force f using an observer called "Kalman filter."

Figure 4: LQG Control Structure.

Use lqry to calculate a suitable state-feedback gain K. The gain K is chosen to minimize a cost function of the form

where the parameter beta is used to trade off performance and control effort. For beta = 1e-4, you can compute the optimal gain by typing

Assuming the x- and y-axis are decoupled, you can use these two regulators independently to control the two-axis rolling mill.

Cross-Coupling Effects

Treating each axis separately is valid as long as they are fairly decoupled. Unfortunately, rolling mills have some amount of cross-coupling between axes because an increase in force along x compresses the material and causes a relative decrease in force along the y axis.

Cross-coupling effects are modeled as shown in Figure 5 with gxy=0.1 and gyx=0.4.

Figure 5: Cross-Coupling Model.

To study the effect of cross-coupling on decoupled SISO loops, construct the two-axis model in Figure 5 and close the x- and y-axis loops using the previously designed LQG regulators:

Note the high thickness variations along the x axis. Treating each axis separately is inadequate and you need to use a joint-axis, MIMO design to correctly handle cross-coupling effects.

MIMO Design

The MIMO design consists of a single regulator that uses both force measurements fx and fy to compute the actuator commands, u_x and u_y. This control architecture is depicted in Figure 6.

Figure 6: MIMO Control Structure.

You can design a MIMO LQG regulator for the two-axis model using the exact same steps as for earlier SISO designs. First, compute the state feedback gain, then compute the state estimator, and finally assemble these two components using lqgreg. Use the following commands to perform these steps:

The MIMO design shows no performance loss in the x axis and the disturbance attenuation levels now match those obtained for each individual axis. The improvement is also evident when comparing the principal gains of the closed-loop responses from input disturbances to thickness gaps x-gap, y-gap:

Note how the MIMO regulator does a better job at keeping the gain equally low in all directions.

Simulink® Model

If you are a Simulink® user, click on the link below to open a companion Simulink® model that implements both multi-loop SISO and MIMO control architectures. You can use this model to compare both designs by switching between designs during simulation.