Modeling and Control of The Paper Machine Drying Section: Slätteke, Ola
Modeling and Control of The Paper Machine Drying Section: Slätteke, Ola
Slätteke, Ola
2006
Document Version:
Publisher's PDF, also known as Version of record
Link to publication
General rights
Copyright and moral rights for the publications made accessible in the public portal are retained by the authors
and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the
legal requirements associated with these rights.
• Users may download and print one copy of any publication from the public portal for the purpose of private study
or research.
• You may not further distribute the material or use it for any profit-making activity or commercial gain
• You may freely distribute the URL identifying the publication in the public portal
Take down policy
If you believe that this document breaches copyright please contact us providing details, and we will remove
access to the work immediately and investigate your claim.
LUNDUNI
VERSI
TY
PO Box117
22100L und
+4646-2220000
Modeling and Control of the
Paper Machine Drying Section
Ola Slätteke
Automatic Control
Modeling and Control of the
Paper Machine Drying Section
Modeling and Control of the
Paper Machine Drying Section
Ola Slätteke
ISSN 0280í5316
ISRN LUTFD2/TFRT--1075--SE
The topic of this thesis is modeling and control of the last part of the paper
machine – the drying section. Paper is dried by letting it pass through a
series of steam heated cylinders and the evaporation is thus powered by
the latent heat of vaporization of the steam. The moisture in the paper is
controlled by adjusting the set point of the steam pressure controllers.
There exist several commercial incentives to focus on the performance
of the moisture control. The time to perform a grade change is often
limited by the moisture and shorter grade change time is directly
correlated to economic profit. Studies have shown that the drying section
uses Ҁ of the total energy requirement in paper making. Reduced
variations in moisture gives opportunity for target shifts (changed set
point) which reduces the amount of raw material and steam requirement.
It also creates opportunity for increased production rate.
The thesis is divided in two parts. The first part deals with the control
of the steam pressure inside the cylinders. Both a black-box model and a
physical model are given for the steam pressure process. A tuning rule for
both PI and PID control is derived and various other controller structures
are investigated. Many of the results are verified by experiments on paper
machines at different paper mills.
The second part of the thesis treats the moisture controller. The
physical model from the first part is expanded with a model for the paper.
This gives a complete simulation model for the drying section that is
implemented in the object-oriented modeling language Modelica. Two
new approaches to control the moisture by feedback are evaluated. The
first utilizes the air around the paper in combination with the drying
cylinders to improve the controller performance. The second uses only the
last part of the drying section to control the moisture, while the first part is
put at an appropriate level. Finally, feedforward of a surface temperature
signal is examined.
Acknowledgements
There are a number of people who have contributed to this thesis. First of
all I would like to thank my advisors Björn Wittenmark, Tore Hägglund,
and Krister Forsman. Our regular meetings have been very constructive
and fruitful, and this thesis would not have been possible without their
outstanding support. At the same time, I have been given a large amount
of independence in my research which is something I have appreciated.
I would also like to acknowledge some of the people at ABB; Per
Sandström, Jonas Warnqvist, Jonas Berggren, and Alf Isaksson. It has
been a great experience working with all of you.
There are many people I have come in contact with at different paper
mills during my research. I would particularly like to mention all of my
old colleagues at Stora Enso Nymölla. It has also been a pleasure getting
acquainted with Stefan Snygg at Stora Enso Hylte, Stefan Ericsson and
Lars Jonhed at AssiDomän Frövi.
I have had the opportunity to work with a few people at the
Department of Chemical Engineering in Lund; Magnus Karlsson, Stig
Stenström, Bernt Nilsson, and Erik Baggerud. Magnus really deserves an
extra salute for the work we have done together; I have learnt a lot from
him.
Much of the work on physical modeling in the last chapter was carried
out on account of a large amount of inspiration by Karl Johan Åström. It
all started as a minor discussion and ended up as a major piece of work.
Working at the Department of Automatic Control in Lund is an honor
and it is a great atmosphere to operate in. I would like to thank all my
colleagues for the years we have had together. I will miss you.
During the last year of my PhD-studies I had the privilege to work for
a month at the Pulp and Paper Centre, University of British Columbia,
Vancouver, under direction of Prof. Guy Dumont. This was an instructive
and very interesting time for me.
Finally I would like to thank ABB and the Swedish Foundation for
Strategic Research (SSF) within the project CPDC for the financial
support of the project.
Preface
I first encountered process control in the summer of 1990. I was working
as a summer intern at a pulp and paper mill at one of their winders (a
machine that slits and winds the paper from the paper machine into the
roll widths ordered by the customer). A winder does not have much
process control but one night shift I was assigned to manage a pulper (a
unit for slushing paper into pulp). I got a two minute crash course in
control theory by one of the operators. For the first time in my life I heard
words like set point and control signal. I remember that I did not
understand much of it at that time. There were two important control
loops to keep an eye on, the level control and the consistency control.
Both were controlled by single-loop controllers, manufactured by Fisher
& Porter, if I remember it correctly. A dangerous operating point was if
the consistency was too high to physically empty the pulper at the same
time as the level was too high to dilute the pulp mix. I promised the
operator to not reach that point and hoped that I was right. Luckily I
managed to do fine through the night and I was placed there the following
nights too.
The next summer I was working at the same site but this year at the
instrument department. One day we were replacing a malfunctioning flow
gauge at the pulp dryer and I was watching a level controller at the
instrument panel, trying to understand how it worked. I noticed that the
level was too low but the controller only opened the valve by 40% and it
was increasing slowly. I asked the maintenance guy who was dismounting
the flow meter, why the valve was not fully opened. I thought that was the
appropriate thing for the controller to do if the level was low. He then
explained to me the concepts of dynamics, overshoot and stability, and
from that day on I was hooked on the exciting field of process control.
During my studies I continued to work at the instrument department
each summer. I learned a lot, things that are still useful for me today,
every thing from repairing old pneumatic controllers with liquid solvent,
programming the DCS-system and understanding different control
structures. After my degree I worked there for a few years more before I
went back to the university to become a PhD student.
Contents
1. Introduction ........................................................................................ 13
1.1 Introduction and motivation .................................................... 13
1.2 Outline and contribution of the thesis ..................................... 19
2. Fundamentals of the Paper Drying Process ..................................... 22
2.1 Cylinder configurations in the drying section ......................... 23
2.2 The steam and condensate system........................................... 25
2.3 The moisture control loop ....................................................... 29
2.4 Disturbances in the drying section .......................................... 38
2.5 A note on the choice of units................................................... 40
PART 1.
Modeling and Control of the Steam and Condensate System
3. Black-box Models and Controller Structures .................................. 45
3.1 A black-box model structure í the IPZ transfer function ....... 46
3.2 PID control of the steam pressure ........................................... 53
3.3 Improved set point response by feedforward .......................... 58
3.4 A state feedback controller...................................................... 63
3.5 A two-pole model of the steam pressure ................................. 69
3.6 The differential pressure loop ................................................. 74
3.7 Summary ................................................................................. 77
4. A Physical Model of a Steam Heated Cylinder................................ 79
4.1 The model................................................................................ 80
4.2 Time and frequency domain analysis...................................... 89
4.3 Comparisons with plant data ................................................... 91
4.4 A modified model ................................................................... 94
4.5 Summary ................................................................................. 97
5. A Tuning Method for IPZ Models .................................................... 98
5.1 A design method based on optimization ................................. 99
5.2 The IPZ tuning rule for PI control......................................... 103
5.3 The IPZ tuning rule for PID control...................................... 109
5.4 Stability regions..................................................................... 113
5.5 Industrial verification of the tuning rule................................ 115
5.6 Comparison between PI and PID control .............................. 118
5.7 Comparison to other design methods .................................... 121
5.8 Summary ............................................................................... 138
PART 2.
Modeling and Control of Paper Moisture in the Drying Section
6. Enhanced Moisture Control Using the Air System ....................... 143
6.1 A literature review of drying section models ........................ 144
6.2 The model.............................................................................. 144
6.3 A prestudy ............................................................................. 149
6.4 Mid-ranging........................................................................... 150
6.5 Moisture control by mid-ranging the air system ................... 155
6.6 Summary ............................................................................... 166
7. Feedforward from a Paper Surface Temperature Measurement 168
7.1 The peak position í the position of a dry surface ................. 169
7.2 Design of a feedforward controller ....................................... 175
7.3 Simulations............................................................................ 178
7.4 Summary ............................................................................... 182
8. Object-Oriented Modeling and Predictive Control of the Moisture
Content .............................................................................................. 183
8.1 The model.............................................................................. 184
8.2 Steady-state model validation ............................................... 194
8.3 Open loop simulations........................................................... 196
8.4 Control of moisture by mid-range MPC................................ 199
8.5 Summary ............................................................................... 207
9. Conclusions ....................................................................................... 209
9.1 Summary ............................................................................... 209
9.2 Future work ........................................................................... 211
1
Introduction
Paper is used for printing and writing, for wrapping and packaging, and
for a variety of other applications ranging from kitchen towels to the
manufacture of building materials. It simply comes in an enormous
variety of qualities. Some common types of paper qualities include the
following:
13
Chapter 1. Introduction
to the super market. An empty mailbox each day you come home from
work. No thesis to hold in your hand right now. The world simply became
a better place to live in with the advent of paper some 2000 years ago.
The pulp and paper industry is a highly competitive and capital-
intensive market that is under increasing price pressure. The price
pressure on the finished products implies that the margins are often small
and a producer can only be profitable by manufacturing high volumes
[Duncan, 2003]. In Europe, the total production of paper in 2003 was 95
million tonnes with a turnover of €72 billion [CEPI, 2004]. Customers are
demanding lower costs, better terms of delivery, and higher product
quality. In the last decade a large number of company acquisitions and
mergers has taken place in the forest industry all over the world as an
answer to the high competition [FFIF, 2004], see Table 1.1. Compared to
other industries such as food, chemical, and pharmaceutical, the paper
industry has delivered a relatively low return on capital employed
(ROCE). As a result the forest industry companies have grown by size and
the industry has become more consolidated. The main objectives behind
the mergers and acquisitions are lower production costs, less sensitivity to
economic fluctuations, reduced transportation costs, reduced labor costs,
and other positive synergy effects. Companies have realized that it might
be cheaper (and certainly quicker) to buy production capacity rather than
building it. At the same time there is a steady overcapacity in the world,
the industry is facing increasing environmental requirements and there is
an increased competition from other industries as alternatives to fiber
products appear [Dumont, 1988]. The plastic packaging demand is e.g.
expected to have a rapid growth in coming years. Therefore the
production of paper requires constant attention on process efficiency,
increasing productivity, and lower costs.
Table 1.1 Figures illustrating the consolidation trend in Europe with less number of
companies and paper machines, and yet a higher capacity [CEPI, 2004].
1991 2001 2002 2003
Number of companies 1 042 918 901 884
Number of paper machines 2 181 1 863 1 811 1 815
Employment 362 100 288 700 285 000 279 400
Turnover (million euros) 39 263 77 028 74 235 71 866
Capacity (1000 tonnes) 72 343 100 713 103 489 104 978
Consumption (1000 tonnes) 62 140 83 306 85 674 86 186
14
Chapter 1. Introduction
The function of a paper machine is to form the paper sheet and remove the
water from the sheet. A paper machine is divided into three main parts,
the wire section, the press section, and the drying section, see Figure 1.1.
When the stock enters the head box in the wire section, it contains roughly
1 % of fibers or less. This low viscous mix is dispensed through a long
slice onto the wire. As it travels on the wire, much of the water drains
away by gravitational forces or is pulled away by suction from
underneath. As the water disappears, the cellulose fibres start to adhere to
one another by hydrogen bonds and form a paper web. When the paper
web leaves the wire section and enters the press section, the dry solids
content is around 20 %. In the press section, the newly formed sheet is
pressed between rotating steel rolls and water is displaced into a press felt.
After a few press nips the web enters the drying section with a solid
content of approximately about 50 %. It now encounters the dryer
cylinders. These are large hollow metal cylinders, heated internally with
steam, which dry the paper as it passes them. Finally, the paper is wound
up on a big roll and removed from the paper machine. The moisture
content is now roughly 5í10 %.
Although the drying section is only responsible for removing less than
1 % of the water volume in the original stock to the head box, this is the
part of the paper machine that, by far, consumes most energy. Studies
have shown that the drying section uses around Ҁ of the total energy
requirement in paper making [Fellers and Norman, 1998]. This implies
that the drying section is the most expensive part of the paper machine in
terms of energy use per kg removed water. Moreover, the drying section
affects a lot of the important physical properties of the final product, such
as paper sheet elasticity, twist, and curl.
~ 140 m
Wire section Press section Drying
Drying section
section
99 % 80 % 50 % 5%
Figure 1.1 The principle of paper production is simple. The water is separated from the
original stock which is smoothened out to a thin and endless paper sheet. By adding
different types of fillers the paper surface obtains different properties. Typical values of
moisture content are indicated. By courtesy of Skogsindustrierna.
15
Chapter 1. Introduction
For a paper mill, and even for a group of companies, erecting a new paper
machine is a large investment. A high production rate and capacity is
therefore essential to achieve a high return on the investment. One of the
most important quality variables in paper manufacturing is moisture
content. Below are a few reasons of why a well tuned moisture control
system provides economic yield.
16
Chapter 1. Introduction
to model errors in the feedforward loop the moisture will
deviate from the set point when the feedback is turned on
again. Hence, the moisture control is important in the last part
of a grade change and a shorter grade change time is directly
correlated to economic profit.
x The moisture control loop is indirectly involved during a web
break by the steam pressure in the steam cylinders. A very
common problem is that the cylinders become overheated since
there is no longer any cooling paper around them. When the
paper web is put back, picking and new web breaks easily
occurs. By having an optimized steam control system during a
web break, the time it takes to get the paper sheet back on the
reel can be reduced.
x Many paper properties depend on moisture content, e.g. curl,
stretch, tear, strength and stiffness [Gavelin, 1972].
These are some of the reasons why the drying section plays a vital role in
paper manufacturing. As the title reflects, this thesis is focused on both
modeling and control of the drying section of a paper machine.
0.8
0.6 Tolerance
limit
Probability density
0.4
0.2
0.0
5 6 7 8 9 10 11 12 13 14 15
Moisture (%)
Figure 1.2 The reduction of moisture variation makes it possible to increase the set point
(target shift). The solid curve represents a condition where the standard deviation has been
reduced by 50 % compared to the dashed curve. Hence, it is possible to increase the set
point from 10 to 11 %. The difference between the tolerance limit and mean value is
sometimes called the give-away.
17
Chapter 1. Introduction
Figure 1.3 A drawing of the first paper machine from 1808, also known as the Fourdrinier
paper machine [Clapperton, 1967]. It was invented in 1798 by Nicholas-Louis Robert,
while working for the French paper mill owned by the Didot family. His machine used a
belt of wire screen to produce a continuous web of paper. He was backed in England by
the Fourdrinier brothers, who built and sold the first paper machines. By 1810, the
Fourdrinier brothers found themselves in bankruptcy and Bryan Donkin, their engineer,
continued to improve the basic design. Soon he was successfully manufacturing a machine
that mechanized the process of making paper. A water and pulp mixture flowed across a
moving, vibrating web of woven wire cloth, forming a wet mat of interlocking fibers.
From the wire, the newly formed paper transferred to a moving web of woolen cloth (the
felt), before being dried.
18
Chapter 1. Introduction
1.2 Outline and contribution of the thesis
19
Chapter 1. Introduction
x Introduction of a new signal which can be used in feedforward
to improve control performance. The signal is based on the
temperature profile in the machine direction.
x Design of a model library in Modelica® including components
for a drying section, with possibility to easily build a dynamic
simulation model of a whole drying section.
Publications
The thesis is based on the following publications:
20
Chapter 1. Introduction
x Slätteke, O.: “Object oriented modeling and predictive control of
the moisture content in paper production,” submitted to American
Control Conference 2006.
Outline
This thesis is divided into two main parts. The first treats modeling and
control of the steam and condensate system, and the second focuses on
modeling and control of paper moisture. The different chapters are
organized as follows.
Chapter 3 presents a black-box model for the steam pressure process and
a few different controller structures are investigated.
Chapter 4 derives a physical model for the steam pressure process that is
compared with the black box model. It is also validated against plant data.
Chapter 6 shows how the air system in the dryer hood can be combined
with the conventional steam pressure control to enhance the moisture
control performance.
21
Chapter 2. Fundamentals of the Paper Drying Process
2.
2
Fundamentals of the Paper
Drying Process
The most common way to evaporate water from the paper web is to use
the latent heat of vaporization in steam. A steam-filled dryer is a cost
effective method to transfer heat into the sheet. The energy in steam has
proven to cost less than a quarter of any other available method [Pauksta,
1998]. Other advantages are low toxicity, easy of transportability, and
high heat capacity. Since most of the heat content of steam is stored as
latent heat, large quantities of heat can be transferred efficiently at a
constant temperature, which is a useful attribute in paper drying and many
other heating applications. Also, the energy can be extracted as
mechanical work through a turbine which makes many mills more or less
self-supporting in terms of electricity. For chemical pulp mills, steam is
obtained simply as a by-product in the chemical recovery process line.
The moist paper can be led around a single large steam heated
cylinder, called Yankee cylinder (mainly used for the drying of tissue) or
a large number of steam heated cast iron cylinders in series (commonly
called cans), called multi-cylinder drying. In this thesis, attention is only
given to the multi-cylinder dryer but most of the theory can also be
applied to the Yankee cylinder, see Figure 3.3. Two thorough textbooks
about paper drying are [Karlsson, 2000], and [Gavelin, 1972]. A glossary
can be found in Appendix A.
22
Chapter 2. Fundamentals of the Paper Drying Process
2.1 Cylinder configurations in the drying section
When the steam enters the cylinder it releases its thermal energy to the
cast iron shell and condenses into water. This condensate is drawn off by
suction with a siphon and fed back to the boiler house. The steam is
typically fed to the cylinders on the backside of the machine (called the
drive side), and the condensate is evacuated on the front side (called the
operator side) or the backside. At some machines, especially on wide
machines, the condensate is removed on both sides. Effective condensate
removal is important for the heat transfer of the dryer cylinder. Therefore,
it is desirable to let some steam pass through the siphon together with the
condensate. This so-called blow-through steam ensures removal of
condensate, air, and other noncondensable gases from the cylinder.
Noncondensable gases reduce the partial pressure of the steam in the
cylinders and lower the condensation temperature at a given total
pressure. In addition, the air molecules tend to accumulate at the cylinder
surface as they can hardly diffuse fast enough against the direction of flow
of the steam and as a result the heat transfer between the steam and the
cylinder shell is reduced. The effect of air in a cylinder is therefore much
greater than would be expected from the average percentage of air in the
cylinder [Gavelin, 1972].
On slow machines (< 300í400 m/min) the condensate forms a pool at
the bottom of the cylinder. It is mainly old board machine running at these
speeds. As the speed increases, the condensate starts cascading
23
Chapter 2. Fundamentals of the Paper Drying Process
24
Chapter 2. Fundamentals of the Paper Drying Process
runnability problems at higher machine speeds, the single-tier
configuration was invented (in 1975 at Stora Enso Hylte mill [Carlberg,
1989]), see Figure 2.2. Using this technique, a single fabric is supporting
the web on both the top and the bottom cylinders, as well as in the passage
between them. Since the fabric is between the web and the cylinders in the
bottom row, no significant drying occurs there. In modern machines, the
bottom row of cylinders is therefore replaced by smaller vacuum rolls to
increase the runnability even more [Asensio and Seyed-Yagoobi, 1992].
PC Steam header
of the group
CYLINDER 1
CYLINDER 2
CYLINDER 3
CYLINDER 4
CYLINDER 5
CYLINDER 6
Condensate
pipe
Figure 2.3 Sketch of a dryer group consisting of 6 cylinders and one common pressure
gauge and controller.
25
Chapter 2. Fundamentals of the Paper Drying Process
Steam header
PC2 PC1
PDC
To the condenser
Tank B Tank A
LC1
LC2
Figure 2.4 Part of a drying section with a cascade system. Sometimes the PDC-valve is
installed in the pipe between group A and tank A.
26
Chapter 2. Fundamentals of the Paper Drying Process
Steam header
2
1 PC
Thermo
compressor
Dryer group
1
PDC
To the condenser
Tank B Tank
LC
Figure 2.5 A steam group with a thermo compressor unit. Both the PC and PDC controller
work in split-range and the numbers indicate in which order the actuators are manipulated.
Some care must be taken since noncondensibles might be recirculated and accumulated.
27
28
system.
PC FI TI PI 1
PC
TC 1 1 2
PC PC PC
2
Condenser TC PI 2
PC PC
2 PC
1
To steam box
PDC PDC
2 PDC
PDC PDC 1
1 PDC To steam box
To condenser
2
LC T1
3.4 m3
FI
Chapter 2. Fundamentals of the Paper Drying Process
From boiler
LC T2 LC T3 LC T4 LC T5
LC 2 m3 6 m3 6 m3 3 m3
T0
3.1 m3
FI
Figure 2.6 A typical piping and instrumentation diagram of the steam and condensate
Chapter 2. Fundamentals of the Paper Drying Process
controller primarily uses the recirculated steam from the thermo
compressor but have the possibility to also use make-up steam from the
header, if necessary. In the same manner, PDC leads the blow-through
steam to the compressor and dumps the steam to the condenser only in
exceptional cases e.g. web breaks or grades with low drying demand.
Sometimes the thermo compressor is used in cascading configurations too.
Naturally, there are many other ways to structure a steam- and
condensate system, and Figure 2.4 and Figure 2.5 are purely simplified
cases to illustrate some basic ideas. The steam and condensate system
almost always uses a case-by-case design. Figure 2.6 shows a typical, but
still somewhat simplified, P&ID of a drying section, indicating both the
complexity and the control engineering challenge.
Figure 2.7 The scanner platform moves the measuring sensor back and forth across the
sheet. By courtesy of ABB Ltd.
29
Chapter 2. Fundamentals of the Paper Drying Process
~1000 m
~10 m
Figure 2.8 The path of the scanning sensor. The large arrow points out the direction of
machine speed. Notice the different length scale in the machine direction and cross
direction.
[Kastanakis and Lizr, 1991]. In [Natarajan, et al, 1988] an algorithm is
developed, which uses least squares to estimate the CD component and
Kalman-filtering for the MD component. It is then further developed in
[Dumont, et al, 1991], and [Chen, 1992]. A similar decomposing
algorithm based on Karhunen-Loeve expansion is [Rigopoulos, et al,
1997], and [Chen and Subbarayan, 1999]. [Chang, et al, 2000] proposes
an elliptic sensor trajectory by variable scanning speed or the use of two
scanners traveling in opposite direction to improve the MD/CD-
estimations.
As stated above, the ultimate objective of these measurements is
control. The primary mechanism today for the control of the moisture MD
variations is the dryer steam pressure. Other methods have been proposed,
like infrared drying [Kuang, et al, 1995], and [Seyed-Yagoobi, et al,
2001], impulse drying [Orloff and Crouse, 1999], and [Martinez, et al,
2001], and Condebelt drying [Lehtinen, 1995], and [Retulainen, 2001].
Most of these methods have been tested in lab-scale for many years but
have not yet found acceptance in industry for various reasons [Crotogino,
2001]. The exception is Condebelt who has one installation in Finland,
which has been running since 1996.
The CD profile, on the other hand, is controlled either by
remoisturizing showers, steam boxes (a device that improves the
vaporization in the paper by adding superheated steam directly onto the
sheet), or by infrared heating boxes located at intervals across the
machine’s width [Dumont, et al, 1993]. An even moisture content in the
CD is easiest to achieve if it is low, therefore it occurs that the paper is
over-dried and then remoisturized. Of course, then the gain in higher
quality has to be weighted against the cost of higher energy use. This
30
Chapter 2. Fundamentals of the Paper Drying Process
4.40
Moisture (%)
4.35
4.30
4.25
0 500 1000 1500 2000 2500 3000 3500
80.2
Basis weight (g/m )
2
80.0
79.8
79.6
79.4
0 500 1000 1500 2000 2500 3000 3500
t (s)
Figure 2.9 Moisture content and basis weight measurements taken from a fine paper
machine. The set point for moisture in this case was 4.3% and the basis weight set point
was 80 g/m2. The 2-sigma values were 0.056% and 0.3 g/m2 respectively.
thesis focuses solely on the MD-control. More details about CD
estimation and control can be found in [Stewart, et al, 2003], [Heaven, et
al, 1994], and [Kjaer, et al, 1995].
The performance of the control system has, in the pulp and paper
industry, historically been described in “2-sigma” or two standard
deviations of the controlled variables. All produced reels of paper leaves
the paper machine together with a “reel-report” that include statistics like
“2-sigma MD”, “2-sigma CD”, and “2-sigma total” for both the moisture
Scanner
Machine direction
Figure 2.10 Structure for the moisture control loop with one scanner device and six steam
groups.
31
Chapter 2. Fundamentals of the Paper Drying Process
Steam pressure
ǻp
and basis weight [Sell, 1995]. These are the average values for the whole
reel, but it is the short term 2-sigma values that are used to make the
decision if the product meets the quality requirements. To focus on the
variability in this way makes sense since a consistent and uniform product
is an important objective, as pointed out in Chapter 1, and the set points of
the quality variables are constant during long periods and only altered at
grade changes. For grade change control, see [Murphy and Chen, 1999],
[Kuusisto, et al, 2002], and [Viitamäki, 2004]. An example of scanner
measurements, in machine direction, during a normal run are shown in
Figure 2.9, taken from a machine producing 80 g/m2 of high quality copy
paper. At 1500 s, there is a short period of time when the measurements
are not updated, most distinct in the basis weight. This is due to the
automatic calibration of the scanner, performed at constant intervals,
when the measuring head is positioned at one of the ends in the CD. Also,
see Section 2.5 for a comment on moisture units.
32
Chapter 2. Fundamentals of the Paper Drying Process
in terms of control. Traditionally, this has been solved by letting all steam
pressure controllers follow the same signal. The moisture controller then
manipulates the steam pressure set point of one dryer group and the others
follow that one, yielding a SISO system for the moisture controller.
Figure 2.10 shows how this can be arranged with one scanner device, also
called measuring frame, and the six dryer groups in Figure 2.6. Dryer
group 5 (called lead group) operates at the highest steam pressure and
receives the control signal from the moisture controller. The set points of
the other groups are then calculated from that value, either as a ratio or a
difference, see Figure 2.11 and Figure 2.12. The purpose of this is
twofold. Firstly, the constant relation between the pressure in the groups
gives good conditions for the function of the cascade system, and
secondly it is also important for both runnability and the quality of the
paper.
The functions f in Figure 2.12 are given by
0 k n d 1, m n 0
°
fn k n r mn where ® or , (2.1)
° k 1, m d 0
¯ n n
f1 ( r ) f 2 (r ) f 3 (r ) f 4 (r) f 6 (r)
Figure 2.12 The set point r from the moisture controller is distributed to the steam
pressure controllers by passing it through a ratio/bias-function except for the lead group,
in this case group 5.
33
Chapter 2. Fundamentals of the Paper Drying Process
–1
the scanner at the end of the machine controls the second part. The middle
scanner can of course also be used for feedforward control.
As indicated above, the moisture control loop is a cascade loop. Drawn
as a block diagram, it looks like in Figure 2.13. The inner loop controls
the steam pressure in the dryer groups. This is in general accomplished by
a PI- or PID controller. In the outer loop there is in general a model based
dead-time compensating controller, typically of the internal model control
(IMC) concept [Morari and Zafiriou, 1989] or based on the Dahlin type
[Dahlin, 1968] (which is a subset of IMC). The performance of these
controllers are evaluated in [Bialkowski, 1996] and [Makkonen, et al,
1995]. The IMC controls the moisture in the paper sheet, by giving set
point values to the PI-controllers in the inner loop. In Chapter 6, a mid-
ranging control structure by combining two IMC-controllers is investi-
Steam pressure
Condition of fabrics
Process air
Web tension
Leakage air
Condensate distribution
Moisture Tuning of controllers
Dry weight
Production speed Dryer section
Layer distribution
Formation
Bulk
Web temperature
Retention
Filler
Refining
Freeness
Blow through steam
Condensate flow
Exhaust air
Figure 2.14 A list of variables that affect the final moisture in the paper during the drying
process. A short description of some of the terms can be found in Appendix A.
34
Chapter 2. Fundamentals of the Paper Drying Process
gated. In this way, two single-loop controllers form a quasi-multivariable
controller. The same structure is also implemented by a true multi-
variable controller but with different manipulated variables, see Chapter 8.
Other non-conventional moisture control schemes can also be found in
[Åström, 1967], [Brown and Millard, 1993], [Xia, et al, 1993], [Rudd and
Schweiger, 1994], [Murphy, et al, 1996], [Wang, 1996] and [Wells,
1999].
Apart from the steam pressure in the cylinders, there are a large
number of variables that determine the moisture in the paper sheet. To
indicate the complexity of the problem some of them are listed below and
given in Figure 2.14.
35
Chapter 2. Fundamentals of the Paper Drying Process
x Ply loading: In paperboard, different layers consist of different
pulps, hence different physical properties. This influences the
drying.
x Blow through steam: In case of improper amount of blow through
steam, the cylinder may be flooded. This has a large influence on
the heat transfer to the paper (and the load on the drives).
Some of these disturbances are controlled variables and can therefore be
regarded as known. This opens for possibilities of feedforward, which
often is the case for production speed and dry weight (see below). Other
variations like inlet moisture, leakage air, or amount of condensate in the
cylinders can be very difficult to measure, and can only be reduced by
feedback. Variations in dryer performance due to conditions of dryer
fabrics can be considered as constant since this is a slowly degrading
process, unless it is unevenly distributed on the fabric.
13
Moisture (%)
12
11
10
142
Cond. weight (g/m )
2
141
140
139
138
137
0 200 400 600 800 1000 0 200 400 600 800
time (s)
Figure 2.15 A case study from a fluting machine. At a first inspection it was found that
there was a very large moisture variation with a period time of two minutes and the steam
control system was thoroughly examined to find the cause. Later it turned out that the
source of the disturbance was a large variation in dry weight and since the drying demand
is correlated with the amount of fibers, the moisture is also affected. The set points are
indicated with dotted lines.
36
Chapter 2. Fundamentals of the Paper Drying Process
1.2
1.0
2-V moisture (%)
0.8
0.6
0.4
0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0
2-V dry weight (g/m2)
Figure 2.16 Each ‘*’ corresponds to the mean value of a 30 minute sample from the
machine in Figure 2.15. It clearly shows the strong correlation between variations in weight
and moisture.
Steam
a)
pressure
Drying Moisture
process
Paper sheet
process Total weight
Stock
flow
Steam
b) pressure
Drying Moisture
process
Paper sheet
Stock process Dry weight
flow
Figure 2.17 There are cross-connections between both steam pressure í total weight and
stock flow í moisture (a), but one of these is easily eliminated by using dry weight as a
controlled variable (b).
37
Chapter 2. Fundamentals of the Paper Drying Process
Apart from moisture, basis weight is also measured at the reel and
controlled by using the stock flow as the manipulated variable (this is the
conventional configuration, there exist others where also the machine
speed is included). The total weight is naturally affected by the amount of
water in the web. However, this coupling is eliminated by instead
controlling the dry weight and leaving only the cross-coupling stock-
moisture behind, see Figure 2.17. Since the moisture is measured, the
amount of water is easy to deduct from the total weight measurement.
As a matter of curiosity, it can be mentioned that one of the first paper
companies to use digital computers to control one of their machines was
Billerud AB in Sweden [Åström, 1967] and [Åström, 2000a]. This was in
the middle of the 1960’s, and the system was an IBM 1710 with a CPU
running at 100 kHz and 80 kB of memory. The system had a special real
time operating system, written as a part of the installation project. All
control was done in a supervisory mode, the digital computer provided set
points to the analog system and it was based on stochastic control theory.
The history of process control in relation to the development of computers
is overviewed in [Balchen, 1999].
38
Chapter 2. Fundamentals of the Paper Drying Process
3
10
1
10
-1
10 -4 -3 -2 -1 0
10 10 10 10 10
0
10
Power spectrum moisture
-2
10
-4
10 -4 -3 -2 -1 0
10 10 10 10 10
Frequency (Hz)
Figure 2.18 Power spectrum for steam pressure (above) and sheet moisture (below), from a
fine paper machine. The dotted lines indicate the 95 % confidence interval.
215
Steam pressure (kPa)
210
205
200
195
0 100 200 300 400 500 600 700 800 900 1000
4.4
Sheet moisture (%)
4.3
4.2
0 100 200 300 400 500 600 700 800 900 1000
Time (s)
Figure 2.19 Time series for the steam pressure and moisture, used for spectrum estimations.
39
Chapter 2. Fundamentals of the Paper Drying Process
2.5 A note on the choice of units
Moisture
There exist two alternatives to express the amount of moisture in the
sheet. Which one is used depends on the context. Books and articles
treating control of the sheet moisture use moisture content, defined as
100mwc
w [%] (2.2)
mwc mds
where mwc is the mass of the water content in the sheet and mds is the mass
of dry solids. This is the quantity used by most control system vendors
and also by staff at the mills.
Alternatively, the amount of moisture can also be expressed as
moisture ratio, defined as
m wc
u [kg/kg] (2.3)
m ds
100u w
w , u (2.4)
1 u 100 w
The advantage of the moisture ratio is the linearity and it better reflects
variations of water content. A change in moisture content by 98 % í 99 %
corresponds to a change of 50 kg water, while a change in moisture
content by 8 % í 9 % corresponds to a change of 0.012 kg of water, see
also Figure 2.20.
Throughout this thesis, moisture content is used in all figures and also
in black-box models regarding sheet moisture. The reason for this choice
is that the, by the author, expected target group for thesis, is much more
familiar to that unit and different results will therefore be more prominent.
The exception for this choice is in Chapter 8 where instead moisture
ratio is used in a few physical relations, since this simplifies them slightly.
40
Chapter 2. Fundamentals of the Paper Drying Process
80
70
60
Moisture content (%)
50
40
30
20
10
0
0 0.5 1.0 1.5 2.0 2.5 3
Moisture ratio (kg/kg)
Pressure
All pressures are given in Pa. The value is always given in absolute
pressure apart from a few exceptions where it is given in gauge pressure.
It is then clearly notified. It is assumed that the absolute value is 101.325
kPa above the gauge pressure.
Temperature
Temperatures in formulas are given in K. In figures, temperatures are
shown in °C for simple interpretation.
41
Chapter 2. Fundamentals of the Paper Drying Process
42
Part Modeling and control
of the steam and
1 condensate system
3.
3
Black-box Models and
Controller Structures
1 sT1 sL
G IPZ ( s ) kv e T1 ! T2 . (3.1)
s 1 sT2
210
200
Pressure (kPa)
190
180
170
0 20 40 60 80 100 120 140 160 180 200
85
Valve position (%)
80
75
70
65
0 20 40 60 80 100 120 140 160 180 200
Time (s)
Figure 3.1 Open loop step response, taken from a liner machine. The response can be
approximated by (3.1) with the process parameters kv=0.0027, T1=58, T2=3.1, and L=2.0.
46
Chapter 3. Black-box Models and Controller Structures
Note that T1 is always larger than T2, typically by a factor of 5 to 50. The
transfer function can then be regarded as an integrator in series with a
lead-network. This gives a characteristic open loop step response. It is
different from most processes normally encountered in the process
industry, which are often approximated by first order systems or pure
integrators, possibly with a time delay. Figure 3.1 shows a step response
from a liner machine. It is performed on the second of totally nine
cylinder groups, which consisted of six dryer cylinders. The basis weight
is 246 g/m2 and the moisture content set point is 7.5 %.
Initially, the steam pressure makes a rapid decrease after the first step
in the valve position, and then there is a distinct breakpoint in the curve
with a significant change in pressure decrease rate. The same thing occurs
at the second step in the control signal. This is the characteristic
appearance of the IPZ-process.
[Stenström et al, 2002] indicates the physics behind this phenomenon.
When the steam valve position is increased it will increase the pressure
inside the cylinder, and consequently the saturation temperature of the
steam and the condensation rate. However, the increasing condensation
525
520
Pressure (kPa)
515
510
505
500
0 100 200 300 400 500 600
31
Valve position (%)
30
29
28
27
26
0 100 200 300 400 500 600
Time (s)
Figure 3.2 Open loop step response, taken from a board machine. The response can be
approximated by (3.1) with the process parameters kv=0.0020, T1=73, T2=21, and L=1.0.
47
Chapter 3. Black-box Models and Controller Structures
rate lags behind the increasing steam inlet flow as the condensate layer
heats up to the new steady state temperature. Therefore, there will be a
fast initial build-up in steam pressure, before the steam consumption has
reached its new value. This and many other physical properties in the
drying section will be further analyzed in Chapter 8. As will be shown
later, the integrator in (3.1) is an approximation of a real pole close to the
origin.
In Figure 3.2, the effect of the integrator in the model structure is
particularly apparent. The valve is first opened a small amount and then
equally closed. When the control signal is put back to the original
position, a new level of steady state pressure is reached. This is the same
behavior as for the level in a tank when the influent or effluent is
manipulated. In closed loop control, this implies that the controller output
always returns to the original level in steady state, in absence of
disturbances acting on the process.
To examine if the dimension of the cylinder has any impact on the
model structure, a step response has been done on a Yankee dryer. These
are mainly used for drying of tissue and have diameters up to 5.5 m
710
700
Pressure (kPa)
690
680
670
660
0 50 150 250 350 450 550
74
Valve position (%)
72
70
68
66
64
0 50 150 250 350 450 550
Time (s)
Figure 3.3 Open loop step response taken from a Yankee dryer. The model corresoinds
to (3.1) with the process parameters kv=0.0026, T1=269, T2=87, and L=2.0
48
Chapter 3. Black-box Models and Controller Structures
[Karlsson, 2000]. This should be compared with multi-cylinder dryers,
which normally are 1.5 or 1.8 m in diameter. The result is shown in
Figure 3.3. It is not a big surprise that the response has an IPZ structure
since the fundamental physical principles must be the same as in the case
of a smaller cylinder. In Chapter 4, we will investigate how physical
properties of the steam and dimensions of the cylinder will affect the
parameters of the linear process model.
q 2(t)
2.5
(t0 ,y0)
Output (kPa)
q 1(t)
1.5
0.5
0
-20 0 20 40 60 80 100 120
Time (s)
Figure 3.4 Simulated unit step response of an IPZ-process with kv=0.01, T1=220, T2=20,
and L=0.
49
Chapter 3. Black-box Models and Controller Structures
tangent to the step response at time t = 0 (called q1(t)) and the asymptote
as time tends to infinity (called q2(t)), both marked in the figure. Suppose
that the size of the step in the control signal is u0 and that the slope of
q1(t) is k1 and the slope of q2(t) is k2. Also, suppose that the two lines
intersect at the coordinates (t0,y0). Then we have
k2 y0
kv , T1 , T2 t0 . (3.2)
u0 k2
The time delay L is obtained in a standard fashion as the time that elapses
between the time when the controller output is changed and the time at
which the response of the process output begins.
To derive the expressions in (3.2), start by denoting the step response
by y(t). To get the final slope of the step response (which also is the
steady state value of the impulse response), we use the final value
theorem
dy (t ) u0 k v 1 sT1
k2 lim lim s 2 G ( s ) lim u kvu0 (3.3)
t of dt s o0 s s o0 1 sT2 0
By the initial value theorem, the initial derivative is
dy (t ) dy (t ) u0 kv 1 sT1 kvT1
k1 lim lim s 2G ( s ) lim u u0 .
dt t 0
t o0 dt s of s s of 1 sT2 0 T2
(3.4)
k v T1u 0
q1 ( t ) t, (3.5)
T2
since the response starts in the origin. By the inverse Laplace transform,
the step response can also be written as
y (t )
k v u 0 T1 T2 1 e t / T2 t , (3.6)
and to get the equation for the asymptote of the step response we observe
that
50
Chapter 3. Black-box Models and Controller Structures
lim y (t ) k v u 0 t k v u 0 T1 T2 , (3.7)
t of
q 2 (t ) k v u 0 t T1 T 2 , (3.8)
and the two lines q1(t) and q2(t) intersect at t=T2, since
q1 (t ) k1t , q2 (t ) k 2t y2 . (3.10)
k1 1
y0 y2 , t0 y2 . (3.12)
k1 k 2 k1 k 2
We also have
k2 k1 y2 y2
kv , T1 , T2 , (3.13)
u0 (k1 k 2 )k 2 k1 k 2
51
Chapter 3. Black-box Models and Controller Structures
and
k1 y2 kvT1 k1
E1 kvT1 , E2 , (3.14)
(k1 k 2 )u0 T2 u0
dT1 k1 y2 (k1 2k 2 )
, (3.15)
dk 2 (k1 k 2 ) 2 k 22
dT1 k1 2k 2 dk 2
(3.16)
T1 k1 k 2 k 2
dk v dk 2 dE1 k 2 dk 2 dE2
, , 0. (3.17)
kv k2 E1 k1 k 2 k 2 dk 2
§ kv · sL
G IPZ ( s ) T 0
¨ k vT1 ¸e (3.18)
2
© s ¹
52
Chapter 3. Black-box Models and Controller Structures
kv
kvT1
Figure 3.5 A unit step response of the IPZ process, if the low-pass filtering part is
removed or equivalently T2 Ł 0.
approximation done when the slow time constant is assumed to be an
integrator.
kc t d
uc (t ) kc E r (t ) y (t ) ³ r (W ) y(W )dW k T
c d J r (t ) y (t )
Ti 0 dt
(3.19)
where uc(t) is the control signal, y(t) is the process output, r(t) is the set
point, and kc, Ti, Td, E, and Ȗ are controller parameters. Often the set point
is not included in the derivative part (Ȗ = 0) and there is also a first order
filter on the derivative part. Written in the frequency domain using the
Laplace operator, (3.19) becomes
§ 1 · § 1 ·
C ff k c ¨¨ E J Td s ¸¸, C c k c ¨¨1 Td s ¸¸, (3.20a)
© Ti s ¹ © Ti s ¹
53
Chapter 3. Black-box Models and Controller Structures
r
Cff (s)
d n
uc y
Cc (s) Ȉ Ȉ P(s) Ȉ
1
Figure 3.6 Block diagram of the PID-controller used in this thesis, together with the
process. r is the set point, uc the control signal, y the measurement, d a load disturbance,
and n measurement noise.
where Cff and Cc are transfer functions, defined by the block diagram in
Figure 3.6. A low-pass filter added to the derivative part gives
§ · § ·
¨ ¸ ¨ ¸
1 J T s 1 T s
C ff ( s ) kc ¨ E d ¸, C ( s )
c k c ¨1 d ¸,
¨ Ti s Td ¸ ¨ Ti s Td ¸
¨ 1 s ¸ ¨ 1 s¸
© N ¹ © N ¹
(3.20b)
where N = 2í20 [Åström and Hägglund, 2005]. The purpose of the filter
is to attenuate the high frequency gain of the derivative part but it can
also be regarded as an inherent part of the control design, like the other
parameters in the controller [Isaksson and Graebe, 2002].
54
Chapter 3. Black-box Models and Controller Structures
745
740
Pressure (kPa)
735
730
725
0 50 100 150 200 250 300
56
Valve position (%)
54
52
50
48
46
0 50 100 150 200 250 300
Time (s)
Figure 3.7 Typical closed loop step response of the pressure loop when using a PID
controller, taken from a fluting machine.
1 ˆ
P( s) P( s ), (3.21)
s
where
and since
where Gry(s) is the transfer function from set point to process output and
Gre(s) is the transfer function from set point to control error, we can write
the transfer function from set point to control error as
55
Chapter 3. Black-box Models and Controller Structures
§ 1 ·1
kc ¨¨ E J Td s ¸¸ Pˆ ( s )
Ti s ¹s
Gre ( s ) 1 Gry ( s ) 1 ©
§ 1 ·1
1 k c ¨¨1 Td s ¸¸ Pˆ ( s )
© Ti s ¹s
§ 1 ·1 § 1 ·1
1 k c ¨¨1 Td s ¸¸ Pˆ ( s ) kc ¨¨ E J Td s ¸¸ Pˆ ( s )
© Ti s ¹s © Ti s ¹s
§ 1 ·1
1 k c ¨¨1 Td s ¸¸ Pˆ ( s )
© Ti s ¹s
Ti s 2 k c Ti s 1 TiTd s 2 Pˆ ( s ) k c E Ti s 1 J TiTd s 2 Pˆ ( s )
Ti s 2 k c (Ti s 1 TiTd s 2 ) Pˆ ( s )
Ti s 2 1 E (1 J )Td s Pˆ ( s ) kcTi s
.
Ti s 2 k c (Ti s 1 TiTd s 2 ) Pˆ ( s )
(3.24)
Now let us study the integral of the controller error as time goes to
infinity and the set point is a step.
f
1
³ e(t )dt
0
lim s
s o0 s2
Gre ( s )
(3.25)
T s 1 E (1 J )Td s Pˆ ( s )kcTi
lim i 2
s o0 T s k (T s 1 T T s 2 ) P ˆ ( s)
1 E Ti .
i c i i d
In the classical case, when E=1, the integral of the control error is zero.
This implies that there must be an overshoot in the set point step response
to compensate the initial positive area in the time plot. We can also
conclude that E>1 will give an overshoot in the step response since the
integral has a negative value. The range of E is normally defined as
0 d E d 1, (3.26)
56
Chapter 3. Black-box Models and Controller Structures
1.5
1
G ry
0.5
0 -4 -3 -2 -1 0 1
10 10 10 10 10 10
1.2
1.0
0.8
G ny
0.6
0.4
0.2
0 -4 -3 -2 -1 0 1
10 10 10 10 10 10
Frequency (1/s)
Figure 3.8 Frequency response from set point to output (above) and frequency response
from noise to output (below) for the closed loop system in Figure 3.7, when governed by
a PID-controller.
57
Chapter 3. Black-box Models and Controller Structures
3.3 Improved set point response by feedforward
P(CM y M u ) PM u M y
G ry My , (3.27)
1 PC 1 PC
Mu P 1 M y . (3.28)
Mu
r
uff
e ufb y
My Ȉ C Ȉ u P
c
í1
Figure 3.9 A two-degree-of-freedom controller (2DOF), where P is the process, C the
controller, and Mu and My are feedforward filters. The advantage of this configuration is
that the servo and regulation problems are separated.
58
Chapter 3. Black-box Models and Controller Structures
Let P be given by the IPZ transfer function
1 sT1 sL
P kv e , (3.29)
s (1 sT2 )
1
My e sL , (3.30)
1 sTcl
then
s (1 sT2 )
Mu
k v (1 sT1 )(1 sTcl )
(3.31)
T2 Tcl T2 1 T1 T2 1
,
k v Tcl T1 k v Tcl (T1 Tcl ) (1 sTcl ) k v T1 (T1 Tcl ) (1 sT1 )
and uff can be generated as the output of the sum of a constant gain and
two low pass filters. This is necessary if the DCS-systems have no
component for high order filters or no possibility to implement the pure
derivative. Observe that the calculations made above are independent of
the feedback controller C. Under ideal conditions, the control signal ufb
remains constant during a set point change and the purpose of C is simply
CONTR
y MV
OUTP uc
My SP
r
Mu FF
Figure 3.10 Implementation of the feedforward structure in Figure 3.9 in a DCS system.
To let the operator manipulate the correct set point, the objects should be grouped
together and presented as one controller element in the operator station, as indicated by
the dashed line.
59
Chapter 3. Black-box Models and Controller Structures
to give good disturbance rejection and robustness. In reality, there will be
a control error during a set point change due to modeling errors which the
feedback loop takes care of by driving e to zero. Figure 3.10 shows how
the feedforward part could be implemented in an industrial control
system to make use of its anti-windup, bumpless transfer, and other
functions. Observe that normally there is a pure time delay in My, which
also should be implementable in the control system.
How should Tcl be chosen? Looking at the shape of the control signal
gives a hint. The transfer function from set point to control signal is
s (1 sT2 )
G ru Mu , (3.32)
k v (1 sT1 ) (1 sTcl )
1.2
1.0
Process output
0.8
0.6
0.4
0.2
0
0 10 20 30 40 50 60 70 80 90 100
0.8
Control signal
0.6
0.4
0.2
0
0 10 20 30 40 50 60 70 80 90 100
Time (s)
Figure 3.11 Set point step response of the 2DOF-controller. Tcl = 5 (solid), Tcl = 10
(dashed), and Tcl = 20 (dash-dotted). The dotted line is the set point.
60
Chapter 3. Black-box Models and Controller Structures
If r is a unit step, the initial value theorem gives
T2
u c ( 0) . (3.33)
kvTclT1
1.2
1.0
Process output
0.8
0.6
0.4
0.2
0
0 50 100 150 200 250 300
0.8
Control signal
0.6
0.4
0.2
0
0 50 100 150 200 250 300
Time (s)
Figure 3.12 Dashed í 2DOF-controller Tcl = 10, solid í PI with aggressive tuning, and
dash-dotted í PI with robust tuning. The thin dashed curve is the set point. Note that the
2DOF-controller cancels process dynamics from r ĺ y but not from load disturbance d,
therefore there is a slight overshoot when the disturbance acts on the system.
61
Chapter 3. Black-box Models and Controller Structures
1.0
Process output
0.5
0
0 20 40 60 80 100 120 140 160 180
1.0
Control signal
0.5
0
0 20 40 60 80 100 120 140 160 180
Time (s)
Figure 3.13 Robustness analysis to modeling errors of the 2DOF-controller. Perfect
process model (solid), doubled time delay in process (dotted), and doubled process gain
(dashed). The dash-dotted line is the set point.
desired value. The physical explanation is that the energy flow from the
steam to the cylinder and paper is slowly increasing which implies that
the steam consumption is slowly increasing. To compensate for this the
controller must add extra steam to the cylinder long after the set point has
been reached. Therefore, even though the set point has been reached the
system has not reached steady-state, and the cylinder and paper
temperatures are still increasing.
Figure 3.12 shows how a well-tuned PI controller behaves without the
feedforward compared to the 2DOF-controller. What is meant by well
tuned is obviously relative. The tuning method used for the PI controller
is introduced in Chapter 5. The method has one user parameter that
determines the robustness of the loop and two different settings are used
in the figure, here denoted as aggressive and robust. The two settings are
chosen to give both a faster and slower response compared to the 2DOF-
controller. The 2DOF-controller gives a smoother performance in set
point response. Since the transfer function from r to y is a first order
system, the frequency response does not have a peak as in Figure 3.8,
which is a nice feature of the 2DOF-controller structure. It is also an
advantage that this structure is easy to implement into most commercial
62
Chapter 3. Black-box Models and Controller Structures
DCS-systems, which is not the case for the controller structure presented
in the following section.
Finally, the sensitivity to modeling errors is investigated. Figure 3.13
shows closed loop simulations where the time delay and gain in the
process is increased by a factor two. The response in process output is no
longer equal to the desired response in presence of modeling errors since
the feedback loop is active, and the second term in (3.27) is not canceled
by Mu. However, the control system proves to handle model errors well.
k v T1 k
s v
sT1 1 T2 T2 b1 s b2
G(s) k v . (3.34)
s ( sT2 1) § 1 · s ( s a1 )
s¨¨ s ¸¸
© T2 ¹
ª a1 1º ªb º
° x (t ) Ax(t ) Buc (t ) « » x(t ) « 1 » uc (t )
® ¬ 0 0¼ ¬b2 ¼ . (3.35)
°
¯ y (t ) Cx (t ) >1 0@ x(t )
63
Chapter 3. Black-box Models and Controller Structures
By definition this system is observable. It is also controllable unless T1 =
T2 which implies that we have a pure integrator instead of an IPZ process.
Let the state feedback be given by
ª a1 b1l1 1 b1l 2 º ªb º
°° x (t ) « b l » x(t ) l r « 1 » r (t )
® ¬ 2 1 b2 l 2 ¼ ¬b2 ¼ , (3.37)
°
°¯ y (t ) >1 0@ x(t )
(b1 s b2 )l r
G ry ( s ) . (3.38)
s a1 b1l1 b2 l 2 s a1b2 l 2 b2 l1
2
By choosing
lr l1 a1l 2 , (3.39)
we get unit steady-state gain. If we choose to place the poles in Ȝ1 and Ȝ2,
we get the equation system
ª b1 b2 º ª l1 º ª O1 O 2 a1 º
«b ». (3.40)
¬ 2 a1b2 »¼ «¬l 2 »¼ «
¬ O1O 2 ¼
Note that we have assumed that the two states of the system are
measurable, which practically is unlikely. What is also needed then is a
Kalman filter to estimate the states, which is feasible since the system is
observable. However, since an observer will not affect the dynamics of
the closed loop system from set point r to output y, except for initial
transients in the estimation, it is adequate to consider simply the state
feedback in this analysis.
64
Chapter 3. Black-box Models and Controller Structures
A Simulation Example
Assume that we have the system
The open loop poles are in p1 = 0 and p2 = í0.1, and the zero is in
z = í0.02. To examine a cancellation of the process zero, place one pole
at Ȝ1 = í0.2, and let the other one be at Ȝ2 = í0.04, í0.02, and í0.01 (in
three different simulations). In other words, we have one case where both
closed loop poles are faster than the process zero, one case where one
pole cancels the zero, and finally one case where one pole is slower than
the zero.
The controller parameters can then be calculated as
l1 1.25 12.5O 2 ,
l2 12.5 1125O 2 , (3.42)
lr 100O 2 .
Figure 3.14 shows the simulation results. We can see that when both
poles are faster than the process zero we get an over-shoot due to the lead
action of the zero (which is closer to the origin). When Ȝ2 = í0.02, we
cancel the zero and get a first order response. In this case, the
performance of the system is given only by the position of the other pole.
Finally, we let Ȝ2 = í0.01, and the closed loop system has one fast and
one slow pole. The slow pole is clearly visible in the step response.
65
Chapter 3. Black-box Models and Controller Structures
1.5
1.0
y, r
0.5
0
0 20 40 60 80 100 120 140 160 180 200
3
uc
0
0 20 40 60 80 100 120 140 160 180 200
Time (s)
Figure 3.14 Simulation of the three different cases of state feedback. Set point (solid), Ȝ2
= í0.04 (dashed), Ȝ2 = í0.02 (dotted), and Ȝ2 = í0.01 (dash-dotted).
ª x1 º ª a1 1 0º ª x1 º ª b1 º ª0 º
« x » « 0 0 0» « x 2 » «b2 » u c ««0»» r.
» « » « » (3.44)
« 2» «
«¬ x 3 »¼ «¬ 1 0 0»¼ «¬ x3 »¼ «¬ 0 »¼ «¬1»¼
uc Lx l 3 x3 , (3.45)
66
Chapter 3. Black-box Models and Controller Structures
r ki uc y
Ȉ Ȉ Process
s
íL
í1
Figure 3.15 Adding integral action to a state feedback system
l3 (b2 b1s )
Gry ( s ) . (3.47)
s (a1 b1l1 b2l2 ) s 2 (a1b2l2 b2l1 b1l3 ) s b2l3
3
The new state, x3, can be added by having an additional feedback loop as
shown in Figure 3.15. The integrator gain, ki, is then given by
ki l 3 . (3.48)
a1 b2 l 2 b1l1 0.42
°
®a1b2 l 2 b2 l1 b1l 3 0.048 , (3.50)
° b2 l 3 0.0008
¯
which has the solution l1 = 3, l2 = 10 and l3 = í0.4. Figure 3.16 shows the
simulation result. Since there is an extra pole in the closed loop system
67
Chapter 3. Black-box Models and Controller Structures
1.0
0.5
y, r
0
1.0
0.5
0
0 20 40 60 80 100 120 140 160 180 200
1.5
uc
0.5
0
0 20 40 60 80 100 120 140 160 180 200
Time (s)
Figure 3.16 A simulation where a disturbance acts on the controller output for the case
when there is a cancellation of the process zero. Without integral action in the controller
(dotted í same as the dotted curve in Figure 3.14), and with integral action (dash-dotted).
Observe that the dynamics are different in the two cases since the number of closed loop
poles is different.
due to the integrator state, x3, the dynamics of the set point step response
is different from the system in (3.42).
A Remark
In Section 3.3 the overshoot in set point step response, as the standard
PID-controller introduces, is eliminated by canceling process dynamics
by feedforward. Here we have shown how the process zero can be
cancelled by one of the closed loop poles introduced by feedback. By
letting the remaining two poles (in the case of integral action in the
controller) be real, we avoid both overshoot in the step response and
resonance in the frequency response.
The practical disadvantage with state feedback is that most
commercial DCS do not have any support for this controller type. Instead,
it must be implemented by the user together with features like windup
protection, bumpless transfer, and bumpless parameter change.
68
Chapter 3. Black-box Models and Controller Structures
3.5 A two-pole model of the steam pressure
Since the steam pressure in a cylinder group finally must reach some
steady state value, the integrator in (3.1) is more likely to be a long time
constant. This implies that the model does not capture the low frequency
behavior well and this has been one of the criticisms against it. However,
the purpose of the model is not simulation and the question is if its
complexity is sufficient for a control design usage. With an extra pole, we
get
1 sT1
G PPZ ( s ) kp e sL , T3 ! T1 ! T2 , (3.51)
1 sT2 1 sT3
instead of (3.1). Here the integrator is replaced by a new pole close to the
origin, and instead we call the model structure PPZ (pole-pole-zero).
Since it in most practical cases is not desirable to let the steam pressure in
an open loop step response reach steady state, because of the influence on
the paper qualities, this slow pole can be difficult to identify graphically
280
Pressure (kPa)
260
240
220
1600 1800 2000 2200 2400 2600 2800 3000 3200
45
40
Valve position (%)
'u
f
35
30
25
20
15
1600 1800 2000 2200 2400 2600 2800 3000 3200
Time (s)
Figure 3.17 Closed loop step response. Taken from a liner machine.
69
Chapter 3. Black-box Models and Controller Structures
from a single step. It can be identified from well excited closed or open
loop data with system identification techniques for linear parametric
models, such as the prediction-error methods, though.
In Figure 3.17 a closed loop step response from a steam group in a
liner machine, consisting of 13 cylinders, is shown. We will now see that
the fact that the control signal does not return to the same steady state
value as before the step in the set point, exposes the PPZ-structure of the
process.
Start by defining the change in control signal as time goes to infinity,
and assume a step in the set point with size 'r
§ 'r ·
'uf { lim uc (t ) u0 lim¨ s Gru ( s ) u0 ¸ u0 Gru 0 'r , (3.52)
t of s o0 © s ¹
where Gru(s) is the transfer function from the set point r to the control
signal uc. Assume that the process dynamics are described by the IPZ-
model given in (3.1). Also, assume that the controller is given by (3.20a).
Then we have
§ 1 ·
kc ¨¨ E sJTd ¸¸
C ff ( s ) © sTi ¹
Gru ( s )
1 Cc ( s )GIPZ ( s ) § 1 · sT1 1 sL
1 kv kc ¨¨1 Td s ¸¸ e (3.53)
© Ti s ¹ s sT2 1
kc s s 2JTiTd sE Ti 1 sT2 1
.
2
2
s Ti sT2 1 kv kc s TiTd sTi 1 sT1 1 e sL
and we can see that the final value of the control signal must be equal to
zero. This is obviously because of the integrator in the process and it
means that the control signal in Figure 3.17 should have returned to its
original steady state value (to the level it had before the step in the set
point). If we instead let the process be of PPZ-type, see (3.51), then the
steady state value of the control signal will become
70
Chapter 3. Black-box Models and Controller Structures
C ff 0
'r
1 Cc 0GPPZ 0
(3.55)
kc 1 sETi s 2JTiTd 1 sT2 1 sT3 'r 'r
sTi 1 sT2 1 sT3 kc k p 1 sTi s 2TiTd 1 sT1 e sL s 0
kp
and it depends of the process gain and the size of the step.
Normally, the difference between the PPZ- and the IPZ-model is only
exposed in long time series. Figure 3.18 shows an example from a paper
board machine (a group consisting of 12 cylinders) where almost two
430
420
Validation
410
400
Steam pressure (kPa)
390
380
370
360
350
Identification
340
330
60
valve position (%)
40
20
0
0 1000 2000 3000 4000 5000 6000
Time (s)
Figure 3.18 Model validation showing the difference between the IPZ (dashed) and the
PPZ (dotted) process. Measurements taken from a board machine.
71
Chapter 3. Black-box Models and Controller Structures
hours of logged data is used. Both an IPZ (dotted line) and a PPZ (dashed
line) model have been identified and their simulations are plotted together
with the measured value (solid line). The first part of the time series is
used for identification (marked in the figure) and the second part is a pure
simulation. The estimated PPZ model is
G1 ( s ) 15.5
1 291s e 2.83 s (3.56)
1 797 s 1 74.2 s
and the estimated IPZ model is
G2 ( s ) 0.00421
1 1896 s e 2.83 s (3.57)
s 1 138 s
The difference between the two models is apparent. The IPZ-model does
not capture the low frequency component of the process as well as the
PPZ-model does. By manipulating the parameters in the IPZ-model, the
graphical fit can be improved in the validation part of the figure, to the
other part’s disadvantage.
The question that we posed in the beginning of this section was if the
IPZ model is sufficient for control design purposes or do we need to use a
PPZ model instead? To answer that question we start by quoting [Ljung,
1999]: “Feedback control is both forgiving and demanding in the sense
that we can have good control even with a mediocre model, as long as it
is reliable in certain frequency ranges. Loosely speaking, the model has to
be reliable around the cross-over frequency ( | the bandwidth of the
closed loop system), and it may be bad where the closed loop sensitivity
function is small.”
Figure 3.19 shows the open-loop Bode plot for the models in equation
(3.56) and (3.57). Both in the magnitude and the phase plot, there is a
large discrepancy for low frequencies but the disparity between the
models is negligible in the high frequency region. Also around the cross-
over frequency, the two different process models have similar
appearance. This means that the IPZ-model is adequate and well suited
for control design usage. We also note that compared to the PPZ-model,
the IPZ-model in general has larger phase lag. Thus, we can, by using the
IPZ-model for design, be sure to have a stable controller even if the
process is of PPZ-type.
72
Chapter 3. Black-box Models and Controller Structures
3
10
Magnitude (abs)
0
10
-3
10 -5 -4 -3 -2 -1 0 1
10 10 10 10 10 10 10
-20
Phase (deg)
-40
-60
-80
-100 -5 -4 -3 -2 -1 0 1
10 10 10 10 10 10 10
Frequency (rad/s)
Figure 3.19 Bode plot for the IPZ (dashed) and PPZ (dotted) models.
73
Chapter 3. Black-box Models and Controller Structures
3.6 The differential pressure loop
ª 29 s 1 1.7 s º
ª yP º « 0.186 s (1.9s 1) e 0 » ª uP º
«y » « » (3.58)
¬ dP ¼ «0.0395 122 s 1 e 2 s 343s 1 2.6 s » «¬udP »¼
0.0029 e
«¬ s (2.3s 1) s (3.1s 1) »¼
Steam header
3
2 PC
Flash steam 1
Dryer group
PDC
Flash steam
Figure 3.20 A P&ID of the dryer group where the TITO experiment was performed. The
pressure controller PC primarily uses flash steam (from another group operating at higher
pressure), secondly steam from the compressor, and thirdly live steam from the header.
74
Chapter 3. Black-box Models and Controller Structures
300
280
yP (kPa)
260
240
220
2200 2400 2600 2800 3000 3200 3400
16
14
12
u P (%)
10
6
2200 2400 2600 2800 3000 3200 3400
80
70
ydP (kPa)
60
50
40
2200 2400 2600 2800 3000 3200 3400
26
24
22
u dP (%)
20
18
16
14
2200 2400 2600 2800 3000 3200 3400
Time (s)
Figure 3.21 The step response experiment of the pressure and differential pressure loops.
The identified model is given in (3.58). Taken from a fine paper machine.
75
Chapter 3. Black-box Models and Controller Structures
340
yP (kPa)
320
300
280
2000 2200 2400 2600 2800 3000 3200 3400 3600
14
12
10
u P (%)
4
2000 2200 2400 2600 2800 3000 3200 3400 3600
90
80
ydP (kPa)
70
60
50
2000 2200 2400 2600 2800 3000 3200 3400 3600
8
u dP (%)
2
2000 2200 2400 2600 2800 3000 3200 3400 3600
Time (s)
Figure 3.22 A second step response model of a similar process as in Figure 3.20. Taken
from a fine paper machine.
76
Chapter 3. Black-box Models and Controller Structures
where yP and uP are output and input of the steam pressure process,
respectively, and ydP and udP are output and input of the differential
pressure process, respectively. Most conspicuous of this model is that the
differential pressure valve, udP, has no influence on the pressure, yP. This
is also very clear from the figure. This is a nice property since the
pressure controller can then be tuned without regarding the differential
pressure, ydP. The differential pressure, however, is significantly affected
by the pressure valve, most easily seen in Figure 3.21. This
interconnection has to be regarded when tuning the PDC controller. The
purpose of this loop is to maintain good condensate evacuation and even
short term deviations from the set point will influence the drying effect.
However, since the PC loop has a larger effect on the drying, good
regulation of this is more important and the performance of the PDC loop
is of minor importance. Figure 3.22 shows a similar experiment as in
Figure 3.21, but on a different dryer group. This experiment gives
equivalent results.
It is not a big surprise that also the differential pressure can be
modeled as an IPZ process, just like the pressure. Both are acting on the
same system and must, to some extent, be governed by the same physical
equations. However, it should be noted that in some cases it has shown
sufficient to model it as a first order system. Deep analysis of the
differential pressure loop is not covered in this thesis, and from here on
only the moisture cascade loop will be treated
3.7 Summary
In this chapter a black-box model structure has been presented and some
controller structures have been analyzed. First, the closed loop properties
when using a PID controller are examined. It is then followed by two
different controller structures. The main purpose of this is to avoid
resonances in the frequency response, r ĺ y, as feedback with a simple
PID gives, see Figure 3.8. In Section 3.3 this is accomplished by
canceling process dynamics by feedforward in terms of a two-degree-of-
freedom controller. In Section 3.4 dynamics are instead cancelled by state
feedback. The 2DOF controller is best suited for implementation since
the existing DCS today seldom support state feedback controllers.
A different type of model for the steam pressure has also been
investigated where the integrator is replaced by a pole. The conclusion is,
however, that the integrator model is sufficient for controller tuning
purposes and therefore the natural choice.
77
Chapter 3. Black-box Models and Controller Structures
Finally, a multi-variable model for the pressure and differential pressure
loop has been given. The main conclusion of this is that the pressure loop
can be tuned without regarding the differential pressure.
78
Chapter 4. A Physical Model of a Steam Heated Cylinder
4
4.
Let qs (kg/s) be the mass flow rate of steam into the cylinder, qc (kg/s) be
the condensation rate, and qw (kg/s) be the siphon flow rate. Also, let Vs
(m3) and Vw (m3) be the volume of steam and water in the cylinder, and let
Us (kg/m3) and Uw (kg/m3) be the densities of steam and water. The mass
balances for water and steam are then
d
U sVs qs qc , (4.1a)
dt
d
U wVw qc qw , (4.1b)
dt
d
U susVs qs hs qc hs , (4.2a)
dt
d
U wuwVw qc hs qwhw Qm , (4.2b)
dt
d
dt
mC p , m Tm Qm Q p , (4.2c)
where Qm (W) is the power supplied from the water to the metal, Qp (W)
is the power supplied from the metal to the paper, hs (J/kg) is the steam
enthalpy, hw (J/kg) is the water enthalpy, m (kg) the mass of the cylinder
shell, Cp,m (J/(kgK)) the specific heat capacity of the shell, Tm (K) the
mean temperature of the metal, us (J/kg) and uw (J/kg) are the specific
internal energies of steam and water. For the period of one revolution, the
cylinder is in contact with the paper or dryer fabric one part of the time,
and in contact with the surrounding air the other part. Due to higher
resistance for heat transfer, the energy loss to the air is only a fraction of
the energy flow to the paper or fabric. However, [Janson and Nordgren,
1958] found that the temperature variation during one revolution is
negligible (less that 0.3°C) and the metal temperature is therefore assumed
to be independent of the cylinder rotation. The energy flow to the metal is
given by
80
Chapter 4. A Physical Model of a Steam Heated Cylinder
Paper web
Ts
Tm
Qm Tp
Qp
Steam Condensate
Dryer shell
Figure 4.1 A piece of the cross-section of a drying cylinder, visualizing the assumption on
the temperature profile and the energy flows.
Qm D sc Acyl Ts Tm , (4.3)
where Dsc (W/(m2K)) is the heat transfer coefficient from the steam-
condensate interface to the centre of the cylinder shell, Acyl (m2) is the
inner cylinder area, and Ts (K) the steam temperature. The outer surface
area of the cylinder is assumed to be equal to its inside area. The error is
negligible (less than 5%), because the thickness of the cylinder shell is
much smaller than the outer cylinder diameter. The energy flow to the
paper is given by
Qp D cp AcylK Tm T p , (4.4)
where Dcp (W/(m2K)) is the heat transfer coefficient from the center of the
cylinder shell to the centre of the paper sheet, Ș (unitless) is the fraction of
dryer surface covered by the paper web, and Tp (K) is the paper
temperature. Fraction Ș is between 0.5 and 0.7.
For simplicity, all steam within the cylinder cavity is assumed to be
homogeneous with the same pressure and temperature. From (4.3) and
(4.4), we make the assumption of a temperature gradient in the condensate
layer, cylinder shell, and paper web, as illustrated in Figure 4.1.
81
Chapter 4. A Physical Model of a Steam Heated Cylinder
Equations (4.1), (4.2), and (4.3) are a crude nonlinear model for the steam
and condensate system in the cylinder cavity. To obtain a linear second-
order model, we make a few simplifications.
d
U sV qs qc ,
dt
d
dt
mC p , mTm Qm Q p ,
(4.5)
0 qc hs qw hw Qm ,
Qm D sc Acyl Ts Tm ,
Qp D cp AcylK Tm Tp ,
which are a mass balance for the steam, an energy balance for the metal, a
static energy balance for the water, and an algebraic equation for the
energy flow. Eliminating the variables qc and Qm, the model becomes
82
Chapter 4. A Physical Model of a Steam Heated Cylinder
d
hs U sV qs hs qwhw D sc Acyl Ts Tm ,
dt (4.6)
d
dt
mC p , mTm D sc Acyl Ts Tm D cp AcylK Tm Tp .
Assuming that the steam in the cylinder is saturated, enthalpies hs, hw,
density Us and the temperature Ts, are all functions of the pressure p. The
model can thus be written as
dU s dp
hs ( p )V qs hs ( p ) qw ( p )hw ( p ) D sc Acyl Ts ( p ) Tm ,
dp dt (4.7)
dT
mC p ,m m D sc Acyl Ts ( p ) Tm D cp AcylK Tm T p ,
dt
where the states are pressure p and mean metal temperature Tm. The steam
inlet flow, qs, is the input. The equilibrium gives the relations
0
qs0 hs ( p 0 ) qw ( p 0 )hw ( p 0 ) D sc Acyl Ts ( p 0 ) Tm0 ,
(4.8)
0 D sc Acyl Ts ( p 0 ) Tm0 D cp AcylK Tm0 T p0 .
Hence
q s0 hs ( p 0 ) q w0 hw ( p 0 )
T p0 Tm0 . (4.9)
D cp AcylK
83
Chapter 4. A Physical Model of a Steam Heated Cylinder
dU s d'p
hs ( p 0 )V
dp p p0
dt
§ 0 dhs dh dq dT ·
¨¨ qs qw ( p 0 ) w hw ( p 0 ) w D sc Acyl s ¸¸ 'p
© dp dp dp dp ¹ p p0 (4.10)
D sc Acyl 'Tm hs ( p 0 )'qs ,
d'Tm dTs
mC p ,m D sc Acyl 'p D sc Acyl 'Tm D cp AcylK'Tm ,
dt dp p p0
where the states are expressed in terms of deviations and the equilibrium
point in (4.8) have been used to simplify the expression. Assuming that
dhs dh dq dT
qs0 qw ( p 0 ) w hw ( p 0 ) w D sc Acyl s , (4.11a)
dp dp dp dp
and
D cpK D sc , (4.11b)
dU s d'p dTs
hs ( p 0 )V D sc Acyl 'p D sc Acyl 'Tm hs ( p 0 )'qs ,
dp p p0
dt dp p p0
d'Tm dTs
mC p , m D sc Acyl 'p D sc Acyl 'Tm .
dt dp p p 0
(4.12)
84
Chapter 4. A Physical Model of a Steam Heated Cylinder
side of (4.11b) is generally much larger that the left hand side. This is
further discussed in Section 4.4 and Chapter 8. The inequality in (4.11a)
will be commented and examined later in the simulations, and also in
Chapter 8.
Writing the system in standard state-space form, we find that
x Ax B'qs ,
(4.13)
y Cx D'qs ,
ª dTs º
« Acyl dp Acyl » ª 1 º
« » « dU s »
« h V dU s dU
hsV s » «V »
A D sc «
s
dp dp », B « dp »,
« dTs » « »
« Acyl » « »
« dp Acyl » « »
¬ 0 ¼
« mC mC p ,m »¼
¬ p ,m
(4.14)
C >1 0@, D 0.
s a 22 sz b1 § z O z ·
G ( s ) b1 b1 ¨ ¸, (4.15)
ss a11 a 22 ss O O ©s sO¹
where
85
Chapter 4. A Physical Model of a Steam Heated Cylinder
1 D sc Acyl
b1 , z a22 ,
dU s mC p
V
dp
§ dTs · (4.16)
¨ Acyl ¸
dp Acyl ¸
O a11 a22 D sc ¨ .
¨ dU s mC p ,m ¸
¨ hsV dp ¸
© ¹
Note that the model (4.15) has an IPZ structure. The essential parameters
of the model are
x Cylinder volume V
x Cylinder mass m
x Specific heat capacity of metal Cp,m
x Area of the cylinder surface Acyl
x Steam properties hs , dU s dp, dTs dp
x Heat transfer coefficient Dsc
All parameters, except the heat transfer coefficient, Dsc, are known
beforehand, either by machine specifications or from a physics handbook
which includes a steam table and heat capacities. The heat transfer
coefficient depends on both amount of condensate and its degree of
turbulence, and is very difficult to predict. Therefore it is used to fit the
model to the measured data. Note that it is only the last two items in the
parameter list that depend on the operating point.
The following assumptions have been made in the development of the
model
x No blow-through steam
x The steam in the cylinder is saturated
x Paper temperature is constant
x The thermal dynamics of the condensate is fast compared to the
cylinder shell
x The condition (4.11)
The pole O and the zero z are both proportional to the heat transfer
coefficient Dsc. For large s the transfer function (4.15) is approximated by
86
Chapter 4. A Physical Model of a Steam Heated Cylinder
b1
G ( s) | , (4.17)
s
where b1 does not depend on Dsc. For small s the transfer function can be
approximated by
b1 z
G(s) | , (4.18)
Os
where
b1 z hs
, (4.19)
O dTs dU
mC p , m Vhs s
dp dp
does not depend on Dsc. Therefore neither the initial part of a step response
nor the slope of the asymptote depend on the heat transfer coefficient.
Finally, the relations between the black-box and grey-box parameters are
hs mC p , m
kv , T1 ,
dTs dU s D A
mC p , m Vhs sc cyl
dp dp
dU (4.20)
mC p , m hsV s
dp
T2 .
§ dTs dU s ·
D sc Acyl ¨¨ mC p , m Vhs ¸
© dp dp ¸¹
Note that the velocity gain, kv, above is not normalized by the measuring
or any actuator range.
As can be noticed, the grey-box model does not explain the time delay
often seen in the black-box model. It is important to remember that we
often are dealing with sampled measurements and it has been observed, in
practice, that the time delay of the system often is close to the sampling
time. Moreover, we have seen from the derivation of the grey-box model
that there are neglected dynamics in the model, which may possibly give
rise to an estimation of the dead-time that is larger than the true value.
87
Chapter 4. A Physical Model of a Steam Heated Cylinder
Remark 1
A unit analysis of the expressions in (4.20) shows that T1 and T2 are given
in seconds, while kv is given in Pa/kg. Since kv is a velocity gain, one
would expect its unit to include seconds. However, the inflow of steam is
given in kg/s and the time unit is therefore canceled.
Remark 2
The energy flow to the metal (4.3) can be derived as the solution of
Fourier’s law of heat conduction
Q kAT , (4.21)
x1 T1
³ Qdx
x0
³
kA dT
T0
(4.22)
or equivalently
T0 T1 k
Q kA A(T0 T1 ). (4.23)
x1 x0 'x
88
Chapter 4. A Physical Model of a Steam Heated Cylinder
Remark 3
The steady-state solution to (4.1) and (4.2) is
qs qc qw
Qm qc (hs hw ) (4.24)
Qm Qp
As expected, both the mass flows and energy flows are equal. The energy
flow to the metal, Qm, is given by the condensation rate times the
difference in enthalpy due to condensation. The drying is, in other words,
powered by the latent heat of vaporization of the steam.
0
10
-1
10
-3 -2 -1 0
10 10 10 10
-45
Phase (deg)
-65
-90 -3 -2 -1 0
10 10 10 10
Frequency (rad/s)
Figure 4.2 Frequency properties for different values of the heat transfer coefficient Dsc =
500 (dotted), 1000 (dashed), and 2000 (solid).
dimensions and a steam table. It can also be seen that, a higher heat
transfer coefficient yields a lower steam pressure gain, since there is a
larger heat transfer to the cylinder and a higher condensation rate.
The linearized model has also been simulated with a step in the input
signal, shown in Figure 4.3. As the figures illustrate, a higher heat transfer
coefficient gives a lower steam pressure, at a given time instance, due to a
larger heat transfer to the cylinder shell. This has also been pointed out in
[Nelson and Gardner, 1996]. A greater energy flow through the cylinder
shell gives a higher cylinder temperature even though the steam has a
lower temperature, at any time index. The effect of a larger Dsc, in basic
terms, dominates over the effect of a lower steam temperature. It is
essential to remember that there are two effects taking place in this
pressureítemperature process. The heat transfer coefficient varies
between the simulations and the steam temperature is changing with time,
and the cylinder temperature depends on both.
In (4.11a) an inequality that depends on the operating point and
cylinder dimensions was utilized to make a significant simplification.
Using values from this example, we can examine its justification. Apart
90
Chapter 4. A Physical Model of a Steam Heated Cylinder
6000
4000
' p (Pa)
2000
0
0 10 20 30 40 50 60 70 80 90 100
0.6
0.4
' T (K)
0.2
0
0 10 20 30 40 50 60 70 80 90 100
Time (s)
Figure 4.3 Step responses for different values of the heat transfer coefficient Dsc = 500
(dotted), 1000 (dashed), and 2000 (solid). The upper graph shows the steam pressure and
the lower shows the cylinder temperature.
from the steam properties, we also need an expression for the derivative of
siphon flow rate, qw, with respect to the cylinder pressure. Some
experimental values are given in [Stenström and Svanquist, 1991] and
using those, we find that the right hand side of (4.11a) is 10 to 20 times
larger than the left hand side. The next section shows that the model has a
good fit to experimental data.
To evaluate the accuracy of the grey-box model it has been calibrated and
validated against measurements from a steam- and condensate system.
The experiments have been carried out on a paperboard machine and
signals have been measured with a sampling time of 1 s. The cylinder data
is
91
Chapter 4. A Physical Model of a Steam Heated Cylinder
Since the input signal (steam input flow) is not manipulated directly,
neither measured, a model for a steam valve has to be added to (4.13). A
simple approach is to assume a static linear relationship between the
controller signal, uc, and the steam input flow, qs, namely
qs d v uc , (4.25)
1
D sc , (4.27)
1 G cyl
D scc Ocyl
is given, where Gcyl is the distance into the cylinder where the temperature
is equal to the mean cylinder temperature, and Ocyl is the thermal
conductivity of the cylinder shell. In this example, the cylinder thickness
is 25 mm (the mean temperature, Tm, occurs in the middle of the cylinder
shell), and the thermal conductivity is 50 W/(mK). The heat transfer
coefficient through the condensate is then
92
Chapter 4. A Physical Model of a Steam Heated Cylinder
425
p (kPa)
420
415
0 50 100 150 200 250 300 350 400 450 500 550
60
55
u c (%)
50
45
0 50 100 150 200 250 300 350 400 450 500 550
Time (s)
Figure 4.4 Calibrated model (dotted) and measured data (solid).
375
370
p (kPa)
365
360
0 200 400 600 800 1000 1200 1400 1600 1800 2000
55
50
u c (%)
45
40
35
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Time (s)
Figure 4.5 Validation of model (dotted) against measured data (solid).
93
Chapter 4. A Physical Model of a Steam Heated Cylinder
In [Karlsson, 2000], typical values of the heat transfer coefficient are
given. They depend strongly on the condensate thickness and turbulence,
and can vary between 1000 and 4000 W/(m2K), where 2000 is a nominal
value. Nevertheless, the fact that our estimated value from the model is
within that range, gives support for the legitimacy of the model.
The total mass flow rate of steam, during the experiment, to the drying
section is 85 ton/h. The machine has 93 cylinders, so the average steam
flow per cylinder is 0.254 kg/s. The parameter dv, in (4.26), and the
average valve opening gives the steam flow to the particular cylinder in
the model, namely 0.154 kg/s. The steam flow is likely to vary a great
deal between different drying groups but by comparing the two values, we
know that also the second calibration parameter is realistic.
In (4.29) the grey-box model is compared with the corresponding
black-box model, adjusted on the same data set. The velocity gain in the
models is normalized with the measuring range of the pressure gauge. The
black-box model is not shown in the figures but it gives a slightly better fit
to the data, since it has more degrees of freedom (three parameters to
adjust instead of two). There is also a difference in the parameters of the
transfer functions. The time delay is equal to the sampling time and comes
from the identification procedure (pem.m in System Identification
Toolbox for Matlab).
G grey ( s ) 0.00243
50.1s 1 e s
s 20.4s 1
(4.29)
Gblack ( s ) 0.00176
77.8s 1 e s
s 21.3s 1
The grey-box model has also been validated graphically by using the
control signal values to simulate an output. The model output is then
compared with the measured steam pressure. Figure 4.5 shows such an
evaluation. The excitation in the control signal is generated by a series of
steps in the set point (closed loop), which is not shown in the figure to
keep it clear.
94
Chapter 4. A Physical Model of a Steam Heated Cylinder
model this can be accomplished by disregarding the assumption from
(4.11b), saying that
D cpK D sc . (4.11b)
It was previously discussed that for high moisture contents and high fabric
tensions, this assumption might not be valid. The system matrix in (4.14)
now becomes
ª dTs º
« D sc Acyl dp D sc Acyl »
« »
« h V dU s hsV
dU s »
s
Ac « dp dp », (4.30)
« dT »
« D sc Acyl s »
« dp Acyl (D sc D cpK ) »
« mC mC p ,m »
¬ p ,m ¼
and the other matrices are unchanged. By examining time series where the
IPZ-structure is sufficient with cases where it is not, an explanation to the
modeling problem is found. If Dsc >>Dcp then model (4.13) is adequate
and Dsc § Dcp means that (4.30) is a better structure. This has also been
verified by simulation of the primary DAE system (4.1)í(4.3) in
Modelica, see Chapter 8. Closer examination of (4.30) shows that
Dsc >>Dcp gives a system with one fast pole and one close to the origin.
When Dcp is increased, the slow pole moves along the real negative axis
towards the other pole, and it can then no longer be regarded as an
integrator.
It can be shown, knowing that all factors in the elements of (4.30) are
positive, that the eigenvalues of A´ are real and negative. This is nice since
it would be unsatisfying to have a model that is unstable for some
combinations of physical parameters, when it is clear that the process is
stable in reality. The characteristic equation of (4.30) is
where a11, a12, a21, and a22 are the elements of matrix A´. Identification of
the parameters gives
95
Chapter 4. A Physical Model of a Steam Heated Cylinder
dTs
D sc Acyl
dp Acyl (D sc D cpK )
a11 a22 0,
dU s mC p , m
hsV
dp
(4.32)
2 dTs
KD scD cp Acyl
dp
a11a22 a12 a21 ! 0,
dU s
hsVmC p , m
dp
and since the coefficients of the polynomial (4.31) are all positive, the
system is stable. The roots of (4.31) are
and since both a12 and a21 are positive, the solution has no complex parts.
The relation between the position of the poles and the physical parameters
is a bit more complicated than in (4.16).
It can also be shown that the initial dynamics of the modified model in
(4.30) are equal to (4.17). For large s we have
1
G(s) | . (4.34)
dU s
V s
dp
The low frequency properties will be different for the two models
however, since (4.13) contains an integrator and has no steady state gain.
For small s the modified model becomes
hs (D sc D cpK )
lim G ( s ) . (4.35)
s o0 dT
D sc AcylKD cp s
dp
96
Chapter 4. A Physical Model of a Steam Heated Cylinder
4.5 Summary
97
Chapter 5. A Tuning Method for IPZ Models
5.
5
A Tuning Method for IPZ
Models
From practical experience, it has been seen that there can be significant
disturbances in the drying section, which makes disturbance attenuation
an important issue. Since the prevalent machine design standard is to
connect all the dryer groups to the same steam header, a disturbance in
one group easily affects the other groups. In this section, a simple tuning
rule for both PI and PID control is presented. Many tuning methods for PI
and PID control have been proposed previously but most of them
exclusively suppose a first-order system with dead time or an integrating
process with dead time. Since this does not fit very well to the IPZ-
process, a new method is necessary. The design goal is to obtain good
load disturbance response. The tuning rule is based on the four process
parameters of the IPZ transfer function and is therefore denoted as IPZ-
tuning. The idea is to provide operators and control engineers at the paper
mill with a simple tuning tool, which can be used without any complex
optimization calculations. To give the user the option to balance between
robustness and performance, the tuning rule has a design parameter. This
parameter is derived from the Nyquist stability theory and is
dimensionless, which is a nice property since it is then independent of unit
selections. The IPZ-tuning rule is tested and evaluated on an industrial
paper machine. It is also compared with a few other design methods.
98
Chapter 5. A Tuning Method for IPZ Models
5.1 A design method based on optimization
The word optimize comes from the Latin word optimus which means the
best (maximize comes from maximus which means the highest).
Optimization methods are an important tool in nearly all engineering
domains, and the control field is indeed no exception. Many control
methods are based on optimization of some vital criteria, subject to one or
a few constraints. However, it is important to consider how the
optimization problem is set up, since this obviously very much affects the
result and it is not necessarily a good controller that comes out of an
optimal solution. If the optimization problem is not correctly formulated,
the optimal controller might not even be stable. The formulation of the
problem also affects how easily it is solved. There exist very powerful
numerical tools that solve linear or quadratic programming problems,
while non-convex optimization is much more demanding due to less
effective numerics and the existence of several local optimal solutions.
The derivation of the IPZ-tuning rule is based on a design proposal in
[Åström, et al, 1998] and [Åström and Hägglund, 2004]. The idea is to
maximize the integral gain, defined as ki = kc / Ti, subject to a robustness
constraint. Due to the form of the constraint, the optimization problem is
non-convex. Therefore, much effort is put on finding a simple relation
between the process parameters and the optimal solution of controller
parameters, so that the user is relieved from the optimization issue.
Instead, he or she will get an approximate solution from just a few button
pushes on a pocket calculator.
By maximizing the integral gain, the absolute value of the integrated
error (IE) of a step load disturbance is minimized. This can easily be seen
by using the nomenclature in Figure 3.6 and writing the error as
E (s) R( s) Y ( s)
§ P( s) P( s )C ff ( s ) ·
R( s ) ¨¨ D( s) R( s ) ¸¸ (5.1)
© 1 P ( s )C c ( s ) 1 P ( s )C c ( s ) ¹
§ P ( s )C ff ( s ) · P( s)
¨1 ¸
¨ 1 P ( s )C ( s ) ¸ R ( s ) 1 P ( s )C ( s ) D( s ).
© c ¹ c
Assume a step unit load disturbance and that the set point is zero. Evaluate
the integral of the error by the final value theorem
99
Chapter 5. A Tuning Method for IPZ Models
f
P( s) 1
IE ³ e(t )dt
0
lim
s o0 1 P ( s )Cc ( s ) s
P( s) 1
lim 2
(5.2)
s o0 1 Ti s TiTd s s
1 P( s) kc
Ti s
Ti P ( s ) Ti
lim .
s o0 Ti s kc (1 Ti s TiTd s 2 ) P ( s ) kc
Positive area
Control error
Negative area
Time
Figure 5.1 A characteristic response to a load disturbance on the IPZ process, controlled
by a PID controller. It is crucial add a constraint to the optimization problem to avoid an
oscillatory solution, since equally large positive and negative areas cancel each other by
the IE criteria.
100
Chapter 5. A Tuning Method for IPZ Models
Im
R0
Re
–1
Figure 5.2 The robustness constraint of the design method. By varying the radius of the
circle, the degree of robustness is changed.
1
S (s) , (5.3)
1 L( s )
where L(s) is the loop transfer function. The maximum sensitivity, Ms, is
then given by
1
Ms max S (iZ ) max . (5.4)
Z Z 1 L(iZ )
Since |1+L(iZ)| is the distance from a point on the Nyquist curve to the
critical point í1, the shortest distance from the Nyquist curve to the point
í1 is thus 1/Ms. Therefore, we get
1
Ms , (5.5)
R0
and we can then use Ms as our design parameter, when deriving the tuning
rule. The nice thing about the maximum sensitivity function is that it
connects the open loop Nyquist curve with a closed loop property. It is
also dimensionless which is a nice property. The disadvantage is that it
can be difficult to relate to for the unfamiliar user. A common question,
e.g., can be what Ms = 1.2 means in practice. Simulations and
familiarization is one answer to that issue.
101
Chapter 5. A Tuning Method for IPZ Models
Figure 5.3 The complete Nyquist curve for the IPZ-process and a PID controller if the
dead time is assumed to be zero. Otherwise, there will be the typical circular appearance in
the origin.
The shape of the Nyquist curve is changed by varying the controller
parameters, in order to avoid the Ms-circle while achieving a high integral
gain. The derivative time Td, affects the curve close to the origin (high
frequencies), the integral time Ti affects low frequencies, while the gain kc
affects all frequencies equally. However, a few properties of the Nyquist
curve are independent of the controller parameters. Since both the IPZ-
process and a PID controller contain one integrator each, the Nyquist
curve starts at the phase lag íS. At high frequencies the phase is íS/2, if
we assume no dead time. The complete Nyquist plot then appears as in
Figure 5.3.
Note that this design method does not give any suggestion about the
feedforward parameters E or Ȗ of the controller, see (3.19), since these
parameters have no influence on the sensitivity function or the disturbance
rejection.
The design method puts emphasis on disturbance rejection. In
[Slätteke, et al, 2002] it was shown that for a certain value of the
maximum sensitivity (Ms=1.2), maximizing the integral gain or
maximizing the bandwidth of the closed loop system (from set point to
measurement) gave essentially the same controller parameters. That
means that emphasis is put on both the regulator and servo problem. This
does not always hold though, even if the resulting controllers from the two
criteria prove to be close.
102
Chapter 5. A Tuning Method for IPZ Models
Remark
Apart from IE, a few other minimization criteria have been proposed in
different contexts, e.g. the integrated absolute error (IAE), integrated
square error (ISE), and integrated absolute time weighted error (IATE). A
disadvantage with ISE is that it penalizes large errors and gives a small
but long-term error. IAE have the advantage to avoid oscillatory
responses. However, if the system is well damped, IE and IAE give
similar results. Keep in mind that IE = IAE, if the error is positive. This is
also noted in [Åström, et al, 1998]. IE is therefore chosen as the
optimization criterion, since it also is computationally effortless.
A large number of different IPZ models, see (3.1), have been tuned
according to the method described in the previous section, with respect to
different Ms-values. The Ms-values have been chosen to give practically
reasonable and sound controller settings. The models have been chosen to
represent a large variety of cylinder pressure processes encountered at
different paper machines, producing different grades and qualities. The
models cover many different cases, from L << T2 to L > T2, T1 from 50 to
800, and T2 from 2 to 400. However, process parameters T1, T2, and L are
presumed to be positive. This implies that, e.g. non-minimum phase
systems are not included. Furthermore, T1 is always larger than T2 but
cases where T2 is almost as large as T1 are also included in the batch of
models. The purpose of this is to cover the whole range of models found
in the industry. This is also of interest in other process areas since a PPZ
structure appears when two first order models are added (parallel paths in
integrated plants) and if the two time constants differ in magnitude, it can
be modeled as an IPZ process.
The result for Ms = 1.1 is plotted in Figure 5.4, where each circle
represents one set of process parameters and the corresponding optimal
controller parameters. In the plot, the controller gain and integral time
have been drawn in combination with the process parameters to give a
linear dimensionless relationship. The advantage with dimensionless
parameters is that they are independent of the selection of units. Whether
the unit of the velocity gain kv, for example, is given in kPa/(%s) or
simply sí1, does not affect the plot in Figure 5.4.
Even though there is some experimental work in finding combinations
that give a clear relationship between the parameters, there is also a
structure in the selection. Due to their units, kv and kc belongs together, in
103
Chapter 5. A Tuning Method for IPZ Models
0.07
0.06
0.05
0.04
kv kc L
0.03
0.02
0.01
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
(T 2 +L /3)/T 1
1.0
0.9
0.8
T i / (T i+L )
0.7
0.6
0.5
0.4
combination with a parameter that is given in seconds (T1, T2, Ti, or L). As
will be obvious shortly, it must be L because of the influence it has on the
controller gain. Moreover, larger T1 compared to T2 give more lead action
in the process, which must results in a smaller controller gain, given a
104
Chapter 5. A Tuning Method for IPZ Models
certain robustness. Therefore the ratio between T1 and T2 (together with a
correction term L/3, that is found ad hoc) is related to kc. Also, it seems
natural to expect Ti depend on the ‘time constant’ T2, compare e.g. lambda
tuning [Morari and Zafiriou, 1989], and here it is also affected by L in a
way that the expressions on both axes, in Figure 5.4, ranges between zero
and one.
The figure shows an apparent linearity. The least squares method has
been used to get an equation for this relationship, given in the controller
parameters. Obviously, linearity is not a necessity. Curves of higher
degrees are possible to use, as long as the desired parameters can be
analytically solved for.
The chosen Ms-values are 1.1, 1.2, 1.3, and 1.4. The resulting
controller parameters are given by
L
T2
3, T 6T2 L
Ms 1.1 : k c 0.09 i 4L ,
T1 k v L T2 21L
L
T2
3, T 23T2 4 L
Ms 1.2 : k c 0.16 i 5L ,
T1 k v L 9T2 105 L
(5.6)
L
T2
3, T 100T2 17 L
Ms 1.3 : k c 0.23 i L ,
T1 k v L 11T2 94 L
L
T2
3, T 88T2 15 L
Ms 1.4 : k c 0.28 i L .
T1 k v L 12T2 85 L
105
Chapter 5. A Tuning Method for IPZ Models
30
Solid: M = 1.1
s
Dashed: M = 1.2
25 s
Dotted: M = 1.3
s
Dash-dotted: M = 1.4
s
20
i
T
15
10
0
0 10 20 30 40 50
T
2
Figure 5.5 The integral time, Ti, as function of the time constant T2, for the PI tuning rule.
For small T2, the two variables are almost equal, otherwise Ti is smaller that T2. The thin
solid line shows Ti = T2, as a comparison. This is also suggested in a tuning method by
[Nelson and Gardner, 1996].
If the time delay in the IPZ-model is very short, the gain will become
increasingly large and the integral time short, in the formulas above. This
is theoretically correct, since the constraint is still fulfilled for the control
loop. By considering the root locus of the feedback loop and assuming
L = 0, shown in Figure 5.6, we can see that an increasing gain does not
move any of the closed loop poles into the right half plane. Neither by
varying the integral time, the loop will become unstable at any time as
long as it is positive. There is no limitation because of system dynamics. It
has infinite gain margin and any combination of positive kc and Ti can be
achieved. But, high controller gain amplifies measurement noise and gives
large control signals, which will saturate actuators. Limitations are instead
given by [Åström, 2000b]
x Sensor noise
x Unmodeled process dynamics
x Actuator saturation
x Admissible control signal variations
x Sensor and actuator resolution
106
Chapter 5. A Tuning Method for IPZ Models
We can conclude that the dead time of the process is directly related to the
stability of the feedback system. However, it is unusual to obtain a dead
time equal to zero, when identifying the steam pressure process, partly
because of the sample time in today’s discrete DCS-systems. Thus, this
should not impose large restrictions to the usefulness of the method.
In Figure 5.7 and Figure 5.8 simulations of an IPZ-tuned PI controller
for different values of Ms and with ȕ = 1 or ȕ = 0. The process is given by
(200 s 1) s
P ( s ) 0.01 e . (5.7)
s (20 s 1)
With ȕ = 0 the control signal smoother and the overshoot smaller for the
step in set point compared to the case with ȕ = 1. This can also be
accomplished with a more sophisticated set point filter. Remember that
the IPZ tuning is based on optimization of disturbance rejection and can
with advantage be combined with e.g. the 2DOF-controller described in
Section 3.3.
Remark 1
By looking at the root locus in Figure 5.6, an interesting observation can
be made. If the dead time is not dominating the process dynamics, higher
controller gain makes the closed loop system less oscillatory. This
counter-intuitive phenomenon can also be seen when controlling a pure
integrating process with a PI or PID controller.
1 1 1
T2 Ti T1
Figure 5.6 The root locus for the IPZ-process in connection with a PI controller for the
case T1 > Ti > T2. Good tuning normally requires T1 >> Ti § T2.
107
Chapter 5. A Tuning Method for IPZ Models
1.5
1.0
r, y
0.5
0
0 20 40 60 80 100 120 140 160 180 200
3
uc
0
0 20 40 60 80 100 120 140 160 180 200
Time (s)
Figure 5.7 Closed loop response of (5.7) and an IPZ-tuned PI controller for Ms = 1.1
(solid), Ms = 1.2 (dashed), Ms = 1.3 (dotted), and Ms = 1.4 (dash-dotted). ȕ = 1.
1.0
r, y
0.5
0
0 20 40 60 80 100 120 140 160 180 200
1.5
uc
0.5
0
0 20 40 60 80 100 120 140 160 180 200
Time (s)
Figure 5.8 Simulation with the same process and controller as in Figure 5.7 but with ȕ = 0.
108
Chapter 5. A Tuning Method for IPZ Models
Remark 2
The root locus also shows why the open loop zero in the process can not
be canceled by a PI controller, as discussed in Chapter 3. It is only by
having infinite controller gain that one of the feedback poles will reach the
process zero in í1 / T1.
Remark 3
The feature of the IPZ-tuning, to give a non-realizable controller when the
process dead time is zero, is not unique for this method. Many other
tuning methods for minimum phase systems have the same property, see
e.g. [Ziegler and Nichols, 1942], [Chien, et al, 1952], [Cohen and Coon,
1953], [Ho, et al, 1995], [Poulin and Pomerleau, 1999], and [Visioli,
2001].
Remark 4
Small values of T2 will generally give a short integral time, Ti by the IPZ-
tuning. Then it is also important to consider the sample time so that the
process is sampled sufficiently fast [Åström and Hägglund, 2005].
The same technique used in the previous section was also used to obtain
tuning parameters for the PID controller. Using (3.20a) the high frequency
gain for the open loop system is
§ 1 · 1 iZT1 kc kvT1Td
lim kc kv ¨¨1 Td iZ ¸¸ e iZL . (5.8)
Z of
© Ti i Z ¹ iZ (1 iZ T2 ) T2
Physically this means that the open loop system has infinite bandwidth
which is unrealistic. To obtain a finite bandwidth, (3.20b) is used instead,
where N is chosen to 10. The result for Ms = 1.2 is shown in Figure 5.9.
The integral time proves to have a more complex appearance that in the
case of PI control, compare to Figure 5.4. This is also the case for the
derivative time and therefore second-degree curves have been fitted to the
optimal parameters, except for the controller gain. The resulting controller
parameters are given in (5.9). We see that the equations are a bit more
complex than compared to PI control (5.6), but have the same property of
giving infinite controller gain and zero integral time for zero dead time.
In Figure 5.10 and Figure 5.11 simulations of an IPZ-tuned PID
controller for different values of Ms and with ȕ = 1 or ȕ = 0.
109
Chapter 5. A Tuning Method for IPZ Models
0.18
0.16
0.14
0.12
kv kc L 0.1
0.08
0.06
0.04
0.02
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
(T 2+L/3)/T 1
1.0
0.9
0.8
T i / (T i+L)
0.7
0.6
0.5
0.4
0.3
0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
T 2 / (T 2+L)
1.0
0.8
0.6
Td / T2
0.4
0.2
0
0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
T 2 / (T 2+L)
Figure 5.9 The relationship between the PID controller parameters and process parameters
for Ms = 1.2.
110
Chapter 5. A Tuning Method for IPZ Models
3T2 L
kc
22k v T1 L
3T2 L
kc
12k v T1 L
3(3T2 L)
kc
20k v T1 L
111
Chapter 5. A Tuning Method for IPZ Models
1.5
1.0
r, y
0.5
0
0 20 40 60 80 100 120 140 160 180 200
4
uc
1.0
r, y
0.5
0
0 20 40 60 80 100 120 140 160 180 200
2
uc
0
0 20 40 60 80 100 120 140 160 180 200
Time (s)
Figure 5.11 Simulation with the same process and controller as Figure 5.10 but with ȕ = 0.
112
Chapter 5. A Tuning Method for IPZ Models
5.4 Stability regions
ki
C ( s ) kc k d s, (5.11)
s
§ k ·
G (iZ )C (iZ ) r (Z )cos M (Z ) i sin M (Z ) ¨ k c i i ik d Z ¸ 1. (5.12)
© Z ¹
Identifying the real and imaginary parts, we find that the boundary of the
stability region can be represented by
cos M (Z ) Z sin M (Z )
kc , ki Z 2kd . (5.13)
r (Z ) r (Z )
(1 Z 2T12 ) 0.5
r (Z ) kv ,
Z (1 Z 2T22 ) 0.5
(5.14)
S
M (Z ) arctan(ZT1 ) arctan(ZT2 ) ZL.
2
Figure 5.12 and Figure 5.13 show examples of stability regions and IPZ-
tuning. The IPZ-tuning gives parameters rather distant from the boundary
113
Chapter 5. A Tuning Method for IPZ Models
4 4
3 3
2 2
ki
ki
1 1
0 0
0 5 10 15 0 5 10 15
kc kc
Figure 5.12 Stability regions for the IPZ-process controlled by a PI controller. Left:
kv = 0.005, T1 = 50, T2 = 5, L = 3; Right: kv = 0.001, T1 = 200, T2 = 1, L = 2. The points
indicate the position of the IPZ-tuning for the different Ms-values in (5.6). The indications
are moving from the origin to larger kc and ki with larger Ms.
10 10
M s = 1.1 M s = 1.2
8 8
6 6
ki
ki
4 4
2 2
0 0
0 5 10 15 0 5 10 15
kc kc
10 10
M s = 1.3 M s = 1.4
8 8
6 6
ki
ki
4 4
2 2
0 0
0 5 10 15 0 5 10 15
kc kc
Figure 5.13 Stability region for the IPZ-process controlled by a PID controller. The
process is given by kv = 0.005, T1 = 50, T2 = 5, L = 3. The stability region forms a three
dimensional surface in the kcíkiíkd plane. Here, it is visualized in two dimensional plots
by looking at the kcíki plane for the optimal kd (constant in each plot) given by the IPZ-
tuning rule.
114
Chapter 5. A Tuning Method for IPZ Models
of stability. With higher robustness (lower Ms-values), the tuning moves
closer to the origin. With lower robustness, the proposed tuning moves in
direction to the maximum ki-value in the stability region. This implies that
if we assume a stable controller, there is lower boundary for the
achievable IE, for a given process. Also, note that, if the dead time of the
IPZ-process is zero, the stability region fills the whole first quadrant, and
any and ki is feasible.
To verify the tuning rule it has been tested on a real paper machine, for a
PI controller and different Ms-values. It is the first group in a fluting
machine consisting of five cylinders. The process transfer function has
been obtained from a step response and is given by
P ( s ) 0.0196
51.6s 1 e 1.30 s . (5.15)
s 7.79 s 1
The controller settings are calculated for all the values of the maximum
sensitivity function in (5.6) and the obtained values are given in Table 5.1.
The results are presented in Figure 5.14 and Figure 5.15 by a series of
closed loop step responses.
The step responses in the figures and the table show very well the
difference between various values of the tuning parameter Ms. A higher
value provides a more quick response and tighter control, but the valve
must also work harder. This is essential to remember since an aggressive
control signal can introduce disturbances in adjacent dryer groups through
the steam header and must be considered when tuning a whole drying
section of a paper machine.
115
Chapter 5. A Tuning Method for IPZ Models
195
190
185
0 100 200 300 400 500 600 700 800 900 1000
55
Valve position (%)
50
45
40
35
0 100 200 300 400 500 600 700 800 900 1000
Time (s)
200
Pressure (kPa)
195
190
185
800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800
55
Valve position (%)
50
45
40
35
800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800
Time (s)
Figure 5.14 Experiment on a real paper machine showing a closed loop step response,
using IPZ-tuning for a PI controller. Ms = 1.1 (above) and Ms = 1.2 (below). Controller
parameters ȕ = 1.
116
Chapter 5. A Tuning Method for IPZ Models
200
Pressure (kPa)
195
190
185
0 100 200 300 400 500 600 700 800 900 1000
55
Valve position (%)
50
45
40
35
0 100 200 300 400 500 600 700 800 900 1000
Time (s)
200
Pressure (kPa)
195
190
185
800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800
55
Valve position (%)
50
45
40
35
800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800
Time (s)
Figure 5.15 Experiment on a real paper machine showing a closed loop step response,
using IPZ-tuning for a PI controller. Ms = 1.3 (above) and Ms = 1.4 (below). Controller
parameter ȕ = 1.
117
Chapter 5. A Tuning Method for IPZ Models
Another important issue is the steam load to the boiler house. Large
variations in the paper machines steam demand can put out the recovery
boiler and this has to be considered when tuning the drying section too.
For Ms = 1.4, a high frequency component becomes evident, most
easily seen at t = 1000 í 1400. This is probably due to actuator
nonlinearities or disturbances (as a result of the tuning, some disturbance
frequencies are amplified by 40 %).
118
Chapter 5. A Tuning Method for IPZ Models
1.0
y
0.5
0
0 50 100 150 200 250 300 350 400
2
uc
-2
0 50 100 150 200 250 300 350 400
Time (s)
Figure 5.16 Comparison between a PI controller (solid) and a PID controller (dotted)
regulating (5.16). Both controllers are tuned for Ms = 1.2, see Table 5.2. The other
controller parameters are chosen as ȕ = 1, Ȗ = 0, N = 10, see (3.20b). There is a set point
change at t = 10, and a load disturbance starting at t = 250 s.
1 50s 3s
P ( s ) 0.01 e . (5.16)
s (1 15s )
The figure clearly shows that the PID controller puts back the process
value to the set point more quickly than the PI controller during the load
disturbance. The process also reacts faster to a set point change with the
PID controller. However, in both those cases the control signal is slightly
more aggressively used by the PID, which might be a disadvantage due to
cross couplings to adjacent steam groups and increased variations in the
steam demand to the boiler house.
Figure 5.17 shows the gain of the transfer function from r to y and also
the sensitivity function. Figure 5.18 shows the corresponding Nyquist
curves. The PID controller gives a somewhat larger peak in the amplitude
curve for Gry. However, this can be adjusted by the feedforward
parameters ȕ and Ȗ of the PID controller, or by a more general set point
119
Chapter 5. A Tuning Method for IPZ Models
1.5
G ry 1.0
0.5
0 -4 -3 -2 -1 0
10 10 10 10 10
1.5
1.0
G ny
0.5
0 -4 -3 -2 -1 0
10 10 10 10 10
Frequency (Hz)
Figure 5.17 Amplitude diagram for the transfer functions from set point r to process
output y, Gry, and from noise n to process output y, Gny. Gny is equal to the sensitivity
function and Gry is equal to the complementary sensitivity function, if ȕ = 1 and Ȗ = 1. PI
í solid, and PID í dotted. Note that the maximum sensitivity is equal for the to
controllers as the tuning prescribes.
1.0 1.0
0.8 0.8
0.6 0.6
0.4 0.4
0.2 0.2
0 0
-0.2 -0.2
-0.4 -0.4
-0.6 -0.6
-0.8 -0.8
-1.0 -1.0
-2.0 -1.5 -1.0 -0.5 0 0.5 -2.0 -1.5 -1.0 -0.5 0 0.5
Figure 5.18 Nyquist curves for the marked controllers in Table 5.2 and the process
in (5.16), PI í left, PID í right. An often occurring difference between the two optimized
Nyquist curves is that for the PID controller the curve touches the Ms-circle in two points
but only once for the PI controller. This can also be seen in the sensitivity function in
Figure 5.17.
filter, so this is actually not an issue of the feedback loop. The sensitivity
function shows that some disturbances are amplified more by the PID
controller. By looking at the spectrum of pressure disturbances in
Figure 2.18, it can be seen that this extra peak in the sensitivity function
120
Chapter 5. A Tuning Method for IPZ Models
coincides with the variations found in that specific case study. Therefore,
the gain of higher performance by PID control is put in the shade by the
amplification of specific disturbances (the problem should evidently be
resolved by finding the root cause of the disturbances and resolve it, if
possible).
The IPZ-tuning, presented in this chapter, is based on stability of the
closed loop system and disturbance rejection, and does not give any
recommendations on the feedforward settings. However, in Chapter 3 a
two-degree-of-freedom controller structure with set point feedforward is
given that can be combined with the PI or PID controller. Then the
regulation problem and servo problem are separated in a nice way and a
design procedure is given for both.
The conclusion is that PID control of the IPZ-process is more effective
than PI control but it will give a controller with slightly more aggressive
use of the control signal, and also a bit more over-shoot in the set point
response for the standard PID structure. It is therefore difficult to
dogmatically say that the one or the other is better. It simply depends on
how sensitive the steam system is to disturbances and how severe the
cross-couplings between the different steam groups are.
This chapter has presented a tuning method for a PI and PID controller
regulating an IPZ-process. Many different tuning methods have been
presented previously, but they mostly suppose a first or second-order
process, see e.g. [Chien, et al, 1952], [Cohen and Coon, 1953], [Zhuang
and Atherton, 1993], [Ho, et al, 1996], and [Skogestad, 2003], or an
integrating process, see e.g. [Poulin and Pomerleau, 1999], [Wang and
Cai, 2001], [Visioli, 2001], and [Chidambaram and Sree, 2003]. Since the
IPZ-process is not a self-regulating process, it is most natural to use the
latter. However, it is not obvious how to match the IPZ-process to a pure
integrator and these methods do not always give a satisfactory result.
Nevertheless, there a few examples, which are modifications of the
Ziegler-Nichols method that give reasonable tuning, e.g. [Tyreus and
Luyben, 1992], and [Poulin and Pomerleau, 1999], but these methods
have no design parameter which is a deficiency.
In this section, a few design methods found in the literature are
examined to see how well they are compared to IPZ-tuning. The purpose
is not to discredit these methods but to demonstrate the need for a new
121
Chapter 5. A Tuning Method for IPZ Models
tuning method for the steam pressure process. All examples are based on
the process model
1 100s s
P ( s) 0.05 e . (5.17)
s (1 20s )
122
Chapter 5. A Tuning Method for IPZ Models
1
y
0.5
0
0 50 100 150
2.5
2
1.5
uc
1
0.5
0
0 50 100 150
Time (s)
Figure 5.19 Evaluation of Ziegler-Nichols frequency response method for PI control (solid
line) given the process in (5.17). Also shown is the IPZ-tuning for Ms = 1.2 (dotted) and
Ms = 1.3 (dashed).
1
y
0.5
0
0 50 100 150
2.5
2
1.5
uc
1
0.5
0
0 50 100 150
Time (s)
Figure 5.20 Evaluation of Ziegler-Nichols frequency response method for PID control
(solid line) given the process in (5.17). Also shown is the IPZ-tuning for Ms = 1.2 (dotted)
and Ms = 1.3 (dashed).
123
Chapter 5. A Tuning Method for IPZ Models
y 1
0.5
0
0 50 100 150
1.5
1
uc
0.5
0
0 50 100 150
Time (s)
Figure 5.21 Evaluation of Ziegler-Nichols step response method for PI control (solid line)
given the process in (5.17). Also shown is the IPZ-tuning for Ms = 1.2 (dotted) and
Ms = 1.3 (dashed).
1
y
0.5
0
0 50 100 150
2.5
2
1.5
uc
1
0.5
0
0 50 100 150
Time (s)
Figure 5.22 Evaluation of Ziegler-Nichols step response method for PID control (solid
line) given the process in (5.17). Also shown is the IPZ-tuning for Ms = 1.2 (dotted) and
Ms = 1.3 (dashed).
124
Chapter 5. A Tuning Method for IPZ Models
T0
PI: kc 0.45k0 , Ti , (5.19)
1 .2
and
T0 T0
PID: kc 0.6k 0 , Ti , Td . (5.20)
2 8
The ultimate gain and oscillation period can also be calculated directly
from the transfer function. Assume the IPZ process and a P controller. The
condition for the point of instability is
1 sT1 sL
kc e 1. (5.21)
s (1 sT2 )
By letting s = iZ, and equating the complex and real parts, the following
controller settings and Ms-values are obtained
125
Chapter 5. A Tuning Method for IPZ Models
Figure 5.23 The two parameters used in the Ziegler-Nichols step response method.
0.9'uc
PI: kc , Ti 3.33L, (5.24)
aL
and
1 . 2 'u c
PID: kc , Ti 2 L, Td 0 .5 L , (5.25)
aL
where ǻuc is the size of the step. Like in the case of the frequency method,
(5.17) is examined by a simulation. The obtained settings are
and
126
Chapter 5. A Tuning Method for IPZ Models
5.7.2 Tyreus-Luyben’s modified ZN tuning rule
The Ziegler-Nichols tuning rules are derived to give decay ratio of ¼. For
many process control system this is too aggressive, which led [Tyreus and
Luyben, 1992] to derive a more conservative PI tuning rule. The method
is obtained by maximizing the closed loop resonant frequency given a
maximum complementary sensitivity function of 2 dB ( § 1.26 ). The
result is given by
k0
PI: kc , Ti 2.2T0 , (5.28)
3.22
Figure 5.24 shows the simulation results and the method gives a
reasonable good tuning. However, the maximum complementary
sensitivity function is different (in this case 1.16) from the value
prescribed by the method, since it is derived for a different process. Both
the set point response and disturbance rejection is faster than the IPZ-
tuning, but it is also less robust (larger Ms) and the control signal is much
more aggressive. As a comparison, Figure 5.25 shows IPZ-tuning for
Ms = 1.65. The disturbance rejection is then better for IPZ-tuning than
Tyreus-Luyben’s tuning rule. The time to reach a new set point is similar
for the two methods even though the size of the over-shoot is different.
Remember that the over-shoot should be dealt with the feedforward part
of the controller, if necessary.
As for the case of Ziegler-Nichols, the Tyreus-Luyben method lacks a
tuning parameter.
k sL
P( s ) e , (5.30)
s
127
Chapter 5. A Tuning Method for IPZ Models
1.0
y
0.5
0
0 20 40 60 80 100 120 140 160 180 200
2
uc
0
0 20 40 60 80 100 120 140 160 180 200
Time (s)
Figure 5.24 Evaluation of Tyreus-Luyben’s method for PI control (solid line) given the
process in (5.17). Also shown is the IPZ-tuning for Ms = 1.2 (dotted) and Ms = 1.3
(dashed).
1.5
1.0
y
0.5
0
0 20 40 60 80 100 120 140 160 180 200
2
uc
0
0 20 40 60 80 100 120 140 160 180 200
Time (s)
Figure 5.25 Evaluation of Tyreus-Luyben’s method for PI control (solid line) given the
process in (5.17). Also shown is the IPZ-tuning for Ms = 1.65 (dotted) (kc = 1.57,
Ti = 4.09).
128
Chapter 5. A Tuning Method for IPZ Models
1
y
0.5
0
0 20 40 60 80 100 120 140 160 180 200
2
uc
0
0 20 40 60 80 100 120 140 160 180 200
Time (s)
Figure 5.26 Evaluation of the AMIGO tuning rule for PI control (solid line) given the
process in (5.17). Also shown is the IPZ-tuning for Ms = 1.2 (dotted) and Ms = 1.3
(dashed).
1
y
0.5
0
0 50 100 150
2.5
2
1.5
uc
1
0.5
0
0 50 100 150
Time (s)
Figure 5.27 Evaluation of the AMIGO tuning rule for PID control (solid line) given the
process in (5.17). Also shown is the IPZ-tuning for Ms = 1.2 (dotted) and Ms = 1.3
(dashed).
129
Chapter 5. A Tuning Method for IPZ Models
the proposed tuning rule is
0.35
PI: kc , Ti 13.4 L, (5.31)
kL
0.45
PID: kc , Ti 8L, Td 0.5L. (5.32)
kL
The next step is to match the IPZ-process to a pure integrator with time
delay. By letting k be equal to the maximum slope of a unit step response
of the IPZ-process, (3.5) gives
k vT1
k . (5.33)
T2
The time delay, L, is equal to the time delay of the true process. Assuming
(5.17), the controller settings and obtained maximum sensitivity are
and
Figure 5.26 and Figure 5.27 show the simulations. We observe the
AMIGO rule gives a pleasant performance for this process. There is less
over-shoot in the set point response, compared to the IPZ-tuning and the
controllers also bring back the process output to the set point nicely at a
load disturbance.
However, if we instead let the process be given by
1 100s 3s
P ( s ) 0.01 e (5.36)
s (1 10 s )
130
Chapter 5. A Tuning Method for IPZ Models
1.0
y
0.5
0
0 50 100 150 200 250 300
1.5
1.0
uc
0.5
0
0 50 100 150 200 250 300
Time (s)
Figure 5.28 Evaluation of the AMIGO tuning rule for PI control (solid line) given the
process in (5.36). Also shown is the IPZ-tuning for Ms = 1.3 (dashed).
and
Figure 5.28 shows the simulation for (5.37) compared with one case of
IPZ-tuning. The AMIGO tuning does not give a satisfactory result at all
for this process. Note that the two controllers in the figure have the same
robustness measure but still very different performance.
131
Chapter 5. A Tuning Method for IPZ Models
Nelson-Gardner’s pole placement rule for PI control
In [Nelson and Gardner, 1996] a tuning rule for the IPZ process is given,
that is derived from simple pole placement. It can also be found in one of
the exercises in [Sell, 1995], written by Bill Bialkowski, but not as
detailed. In both references the idea is given in words without carrying out
the calculations and the time-delay is disregarded. Let the process be
given by
where the time delay have been approximated by the first-order Taylor
series. The controller is given by
1 sTi
C ( s) k c . (5.40)
sTi
Now, let Ti = T2. This implies that the process pole is canceled by the zero
of the controller. This can be compared with the IPZ-tuning rule, where
the integral time follows the parameter T2 in the sense that larger values of
T2 gives larger values of Ti, and for small T2 the two are almost equal, see
Figure 5.5. The closed loop system is now given by
(k v k c ( L T1 )) 2 4k v k c T2 0 (5.43)
132
Chapter 5. A Tuning Method for IPZ Models
1
y
0.5
0
0 20 40 60 80 100 120 140 160 180 200
2
uc
0
0 20 40 60 80 100 120 140 160 180 200
Time (s)
Figure 5.29 Evaluation of Nelson-Gardner’s pole placement rule (5.44) for PI control
(solid line) given the process in (5.17). Also shown is the IPZ-tuning for Ms = 1.2 (dotted)
and Ms = 1.3 (dashed).
4T2
kc , Ti T2 . (5.44)
k v ( L T1 ) 2
2
(5.45)
T1 L
The simulation result is shown in Figure 5.29 and we can immediately see
that this is not a satisfying controller. We can also see that it is a rather
robust controller, since the Ms-value is fairly small. Further investigations
show that this method gives a quite robust controller with Ms-values
133
Chapter 5. A Tuning Method for IPZ Models
below 1.05 for short time delays (L § 1). Larger time delays or shorter T2
give Ms § 1.1í1.2, and relatively fine simulation results. The disadvantage
with this method is the lack of any design parameter but it gives a
reasonable tuning.
PID control
The poles of a second order process with one zero, can be placed
independently of one another by a PID controller. By assuming that the
time delay can be neglected, the IPZ-process fits into this model structure.
Then the process is given by
(1 sT1 )
P( s) kv , (5.47)
s (1 sT2 )
1 Ti s Ti Td s 2
C ( s) k c . (5.48)
sTi
Gry ( s )
kv kc (1 Ti s TiTd s 2 )(1 T1s )
.
Ti (kv kcT1Td T2 ) s 3 Ti (1 kv kcT1 kv kcTd ) s 2 kv kc (T1 Ti ) s kv kc
(5.49)
( s Z1 )( s 2 2]Z 0 s Z 02 ) (5.50)
134
Chapter 5. A Tuning Method for IPZ Models
kv kc
° Z1Z 02
° k v k c T1Ti Td T2 Ti
°° k v k c (T1 Ti )
® 2Z 0 Z1] Z 02 (5.51)
k k T
° v c 1 i dT T T T
2 i
° Ti (1 k v k c T1 k v k c Td )
° 2Z 0 ] Z1
°¯ k v k c T1Ti Td T2 Ti
Even though any Z0, Z1, and ȗ can be realized there should be some
relation between the process and the desired closed loop poles, to give a
controller with satisfactory performance. E.g. requesting a bandwidth that
is too high for the closed loop system gives very high controller gain. This
injects much noise into the loop, can be damaging to actuators, and might
make the system unstable due to neglected process dynamics.
From (5.52) the controller settings can easily be calculated for a given
closed loop response. As an example, we select to place the poles in
3 3 3
Z0 , Z1 , ] . (5.53)
2T1 4T2 4
135
Chapter 5. A Tuning Method for IPZ Models
21 21
20 20
Steam flow (kg/s)
If this is a good tuning depends on how well we can neglect the time
delay, but also on the other process parameters. Large T1, for example,
tends gives very small kc.
Remark
Trying to place any of the closed loop poles in í1/T1 gives infinite
controller gain. This matches what was observed with PI control.
Controller parameters Ti and Td will have finite values though.
136
Chapter 5. A Tuning Method for IPZ Models
2.5
2.0
1.5
ki
ZN step
1.0
ZN freq
0.5
IPZ Ź TL
ƅ AMIGO
0
0 1 2 3 4 5 6 7
kc
Figure 5.31 Stability region for the process in (5.17) controlled by a PI controller. The
different tuning rules examined in this section are also indicated, ż Ziegler-Nichols
frequency method, + Ziegler-Nichols step test method, Ź Tyreus-Luyben’s modified rule,
ƅ the AMIGO tuning rule, Ƒ Nelson-Gardner, and Ɣ IPZ-tuning (Ms = 1.1í1.4). The
dashed line shows the contour that the IPZ-tuning follows. For large Ms values, the tuning
method is moving towards the peak of the stability region, namely the maximum ki.
To see the difference between the evaluated tuning methods from another
perspective, the proposed controller settings by the different methods for
PI control of (5.17) are indicated in the stability region for the process, see
Figure 5.31. The two Ziegler-Nichols methods distinguish themselves by
grouping together away from the other methods. All four settings from the
IPZ-tuning are shown in the figure. It was previously concluded that the
AMIGO tuning rule and Tyreus-Luyben’s rule give satisfactory controller
settings for this set of process parameters. Therefore, it seems reasonable
that their controller settings are in the vicinity of the IPZ-tuning. The
example of Nelson-Gardner’s pole placement is also in that neighborhood
but closer to the origin, which gives a much more robust control.
Remember that it is only the IPZ tuning among the evaluated methods that
have a design parameter and that give good performance for a wide range
of different parameters of the IPZ model.
Finally, a robustness test is given in Figure 5.32, showing a simulation
of the different methods tuned for (5.17) but where the time delay is
137
Chapter 5. A Tuning Method for IPZ Models
7 12
x 10 x 10
2 2
ZN frequency ZN step
0
y 0 y
-2
-2 -4
0 50 100 150 200 0 50 100 150 200
2 2
TL AMIGO
y 1 y
1
0 0
0 50 100 150 200 0 50 100 150 200
2 2
NG IPZ Ms = 1.3
1
y y
1
0
-1
0
0 50 100 150 200 0 100 15050 200
Time (s) Time (s)
Figure 5.32 Robustness analysis of the tuning methods used in comparison with the IPZ
tuning. The time delay in (5.17) is changed from 1 to 2 seconds.
5.8 Summary
In this chapter, a tuning method for PI and PID control of the IPZ model
has been presented. It is based on optimized load disturbance rejection
subject to a robustness constraint. To calculate the controller parameters,
the method requires the four process parameters of the IPZ-model and a
design parameter that is defined by the user. The process parameters can
easily be obtained from an open loop step response. Since the dead time
has a significant effect on the method, good quality dead time estimation
138
Chapter 5. A Tuning Method for IPZ Models
is essential. The advantage with the IPZ-tuning is that the equations are
simple and it is easy to use. From experiments on several paper machines
at different mills, the method has been validated and proved to perform
well.
A few other tuning methods found in the literature have been
compared to the proposed tuning method. The conclusion is that there is a
need for a new method, since neither of the other give good result when
they are evaluated for different process parameters.
Since the method is based on optimized load disturbance rejection it
can with advantage be combined with a set point filter. Examples of this
are the ȕ and Ȗ factors in (3.19), or the more general 2DOF-controller
described in Section 3.3.
The IPZ structure can be used in many other process areas, apart from
modeling the steam pressure in cylinders [Forsman, 2005]. Therefore the
tuning method has a much wider field of applications than solely drying
section control.
139
Chapter 5. A Tuning Method for IPZ Models
140
Part Modeling and control
of paper moisture in
2 the drying section
6
6.
The drying section is enclosed inside a drying hood. The main purposes of
the hood are to create a controlled environment for the drying process,
improve energy utilization, and also to establish good working conditions
in the machine room. The exhaust air removes the evaporated water from
the paper web while preheated dry air is added to the hood by the supply
air. Traditionally, these two variables are used to control the humidity and
pressure inside the hood. This chapter proposes using the supply air, in
combination with the steam cylinders, to control the moisture in the sheet.
Trials have shown that the supply air has a fairly large and fast impact on
moisture. This suggests that using the supply air would give a closed loop
system with higher performance compared to conventional steam pressure
control. However, the drying capacity of the steam cylinder process is
much larger than the drying capacity of the air process and both processes
are therefore essential to obtain a well functioning control system. The
two manipulated variables are combined in a mid-ranging structure to
control the moisture.
The proposed control structure is evaluated by simulations of a
physical model of the drying section.
143
Chapter 6. Enhanced Moisture Control Using the Air System
6.1 A literature review of drying section models
The model used for simulation and analysis in this chapter and Chapter 7
is developed at the Department of Chemical Engineering, Lund
University, see [Karlsson and Stenström, 2005a] and [Karlsson and
Stenström, 2005b]. In this thesis, the model is not described in detail but a
general overview is given. It is adapted to a board machine at a paper mill
in Sweden. It consists of 93 steam cylinders divided in 12 groups. It is
assumed that the machine speed is 430 m/min and the dry basis weight is
267 g/m2.
144
Chapter 6. Enhanced Moisture Control Using the Air System
One cylinder
Output data
Input data
Time Time
Figure 6.1 Illustration of the simulation method, where one cylinder with a following free
draw is simulated at a time.
The model is built on basic physical relations, in terms of mass and energy
balances, and algebraic equations. The continuity equation forms a basis
for all the balances. Written as a partial differential equation (PDE) it is
given by
w)
J k , (6.1)
wt
145
Chapter 6. Enhanced Moisture Control Using the Air System
Figure 6.2 The model used to simulate the closed loop steam pressure dynamics.
simulation of the next cylinder, by following the flow of the paper. In this
way, the number of nodes in each simulation is greatly reduced, and the
memory requirement is reduced from gigabytes to megabytes. Simulations
shown in this chapter still take 1í2 days on a standard computer. The
simulation technique will have some implications for simulation of the
closed loop system, as it will be shown later.
146
Chapter 6. Enhanced Moisture Control Using the Air System
0.00178(97.37s+1) - 1.0 s Range y: (-100 - 450) 0.00350(44.69s+1) - 1.00 s Range y: (-100 - 450)
Steam group 1 G(s) = e Range u: (0 - 100) Steam group 2 G(s) = e Range u: (0 - 100)
s(7.97s+1) s(2.61s+1)
18 Loss function = 0.31596 96 Loss function = 0.38311
16 94
y 92
14 y
90
12
88
10
86
60 70
u 58 u 68
56 66
500 550 600 650 700 750 800 850 900 0 50 100 150 200 250 300
Time (s) Time (s)
60 166
58 164
y 56 y
162
54
52 160
50 158
64 50
u 62 48
u
60 46
58 44
400 450 500 550 600 650 700 750 0 50 100 150 200 250 300 350
Time (s)
Time (s)
93 266
92
262
91
50
48
u u 45
44
40 40
450 500 550 600 650 700 750 0 100 200 300 400 500
Time (s) Time (s)
0.00270(40.31s+1) - 1.0 s Range y: (-100 - 450) 0.00369(78.64s+1) - 1.96 s Range y: (-100 - 450)
Steam group 7 G(s) = e Range u: (0 - 100) Steam group 8 G(s) = e Range u: (0 - 100)
s(6.86s+1) s(13.23s+1)
Loss function = 0.23028 Loss function = 0.21070
226
198
224
196 222
y y 220
194
218
192 216
214
190
28
64
u u 26
62
60 24
700 800 900 1000 1100 1200 100 150 200 250 300 350 400 450
Time (s) Time (s)
Figure 6.3 Graphical output from the modeling tool [Wallén, 2000] used to identify the
steam pressure dynamics included in the simulation model. The solid lines are process
signals and dotted lines the obtained models (also given as transfer functions). The unit of y
is kPa (gauge) and the unit of uc is %. The figure shows group 1í8.
147
Chapter 6. Enhanced Moisture Control Using the Air System
0.00324(75.96s+1) - 2.42 s Range y: (-100-450) 0.00362(74.83s+1) - 1.89 s Range y: (-100 - 450)
Steam group 9 G(s) = e Range u: (0-100) Steam group 10 G(s) = e Range u: (0 - 100)
s(21.23s+1) s(15.73s+1)
Loss function=0.16867 Loss function = 0.30897
258
250
254
246
250
y y
242
246
238
242
234
238
42 42
u 38 u 38
34
34
50 100 150 200 250 300 350 400 450 500 550 600 650 700 750 800
Time (s) Time (s)
0.00139(110.35s+1) - 1.43 s Range y: (-100 - 450) 0.00138(108.94s+1) - 1.70 s Range y: (-100 - 450)
Steam group 11 G(s) = e Range u: (0 - 100) Steam group 12 G(s)= e Range u: (0 - 100)
s(24.95s+1) Loss function = 0.22603 s(24.26s+1) Loss function = 0.16917
326
324 342
322
340
y 320 y
338
318
316 336
314
334
58 52
u 54 u
48
50 44
0 50 100 150 200 250 300 350 400 450 500 550 600 650 700
Time (s) Time (s)
Figure 6.4 Graphical output from the modeling tool [Wallén, 2000] used to identify the
steam pressure dynamics included in the simulation model. The solid lines are process
signals and dotted lines the obtained models (also given as transfer functions). The unit of y
is kPa (gauge) and the unit of uc is %. The figure shows group 9í12.
Figure 6.5 An example of blow box configuration in both a single-tier (above) and two-
tier (below) machine. The oval areas indicate the region where most evaporation occurs.
The active air volume is essentially bounded by machine equipment and fabrics. The large
arrow indicates the direction of the paper and a few blow boxes are indicated by smaller
arrows. By courtesy of Metso Paper.
148
Chapter 6. Enhanced Moisture Control Using the Air System
Evaporation
Exhaust Air
System
boundary
Cylinder
Wet Air
Paper
Dry Air
Supply
Air Actuator
Figure 6.6 A picture of the air system used in the model.
in the dynamics, see Figure 6.6, and the volume was set to 10 m3. With a
nominal total supply air flow of 50 m3/s divided to 93 cylinders, this
corresponds to a residence time of 18.6 s. Dry air is mixed with
recirculated moist air and in this way, the dew point of the supply air is
manipulated while the flow rate is constant. The idea is to prevent
runnability problems due to web flutter in case of aggressive use of the
supply air.
6.3 A prestudy
To investigate the potential to use the supply air as an actuator, a step test
has been performed on a board machine, see Figure 6.7. The machine
speed is 724 m/min and the basis weight is 134 g/m2. The step test is
performed on the last third of the drying section (the after dryer), while
both the moisture and basis weight controller are put in manual mode (the
steam pressure controllers are still in automatic). The air flows are
manipulated by variable speed fan motors (some machines have air flow
dampers instead). Measurements prior to the experiment showed that
67 % control signal corresponds to 41.4 m3/s supply air and 70.3 m3/s
exhaust air. Both the supply air and exhaust air are changed equally
during the step test to maintain the air balance in the hood. Because of the
instrumentation at the particular mill, it is not possible to automatically
log the signal to the fans, therefore the control signal shown in the figure
has been recreated afterwards from notes taken during the experiment. In
addition, the change in supply air is done 5í10 s after the change in
exhaust air, since it is done manually by one person. The moisture and
149
Chapter 6. Enhanced Moisture Control Using the Air System
65
exhaust air ( C)
o
Dew point
60
55
50
100
supply air (%)
Control signal
50
0
0 200
600 400 800 1000 1200
Time (s)
Figure 6.7 Results of the prestudy where it is shown that the supply air has a large impact
on the sheet moisture, measured at the reel-up. The dew point of the exhaust air is also
shown.
dew point are logged in two different systems and are synchronized at a
later stage by the time tags from the internal system clocks. Therefore, it
is considered somewhat uncertain to make a dynamic model from the
experiment result, but it still clearly gives an idea of the potential to use
the air system as an actuator for the moisture controller.
The experiment gives a promising result. The response is distinct, the
change in sheet moisture is rather significant and the time constant is
somewhere between 15í25 s if it is modeled as a first order process, see
Figure 6.7. This also agrees well with the residence time of the air system
in the model described in Section 6.2.
6.4 Mid-ranging
150
Chapter 6. Enhanced Moisture Control Using the Air System
Pulp tower
Water
QC
Figure 6.8 The consistency in the pulp tower is typically 8í12 %. Dilution to pumpable
suspension takes place at the bottom of the tower by injection of dilution water. An
agitator helps the mixing.
[Shinskey, 1978], [Allison and Isaksson, 1998] and [Allison and Ogawa,
2003]. A common example in the pulp and paper industry is consistency
control of the pulp that is pumped from the storage tower, see Figure 6.8.
The controller has two actuators, a large dilution valve and a smaller valve
for fine adjustments. The main purpose of the water through the large
valve is to make the pulp pumpable. Nevertheless, it is also important to
avoid letting too much water being injected at the bottom of the tower to
prevent the small valve from being closed, and vice versa, since the small
valve often has a smaller operating range and finer resolution. One
solution is to let the small valve return to its midpoint or target value in
steady state. This is known as mid-ranging. In general, the effect of the
inputs differs significantly in both range and speed, and sometimes in
cost. When there is a cost tied to the inputs, the optimal target value for
the expensive input is often not the midpoint.
151
Chapter 6. Enhanced Moisture Control Using the Air System
u2
C(s) P2(s)
r2
r1 u1 y
Ȉ kc Ȉ P1(s) Ȉ
1
Figure 6.9 The structure of the simple mid-ranging controller. The set point r2 is the bias of
the P-controller. In practice, C(s) is often chosen as a PI controller
conclusion is that the MPC solution is superior compared to the other two,
and the valve position controller comes in second hand.
Simple Mid-ranging
This is probably the simplest way to implement mid-ranging and it
consists of only one feedback loop but two controllers in parallel with the
same set point, see Figure 6.9. The controller manipulating u1 is a
P controller, whereas the controller manipulating u2 is a controller with
integral action. When the integral action in C removes the control error,
the output u1 is the offset of the P controller. By adjusting the offset, the
target value for u1 is set. The disadvantage with this solution is that P1 is
limited to be controlled by a P controller. This gives a slow response to
both set point changes and disturbance rejections, since the P controller
gives a steady-state error which is removed by C through a much slower
process P2.
152
Chapter 6. Enhanced Moisture Control Using the Air System
r2 u2
Ȉ C2(s) P2(s)
CD(s)
1
r1 y
Ȉ C1(s) Ȉ
u1 P1(s) Ȉ
1
Figure 6.10 The structure of the valve position controller.
P2 ( s )
CD ( s) , (6.2)
P1 ( s )
Hybrid Mid-ranging
An approach, originally used to control pressure of a steam header, is the
hybrid mid-ranging [Love, 1994]. It consists of one PI controller, two
saturations, a low pass filter, and three gains, see Figure 6.11. The
parameter Rm, should be chosen as
100k1
Rm , (6.3)
k1 k 2
where k1 and k2 are the steady-state gains of P1 and P2, respectively. The
filter in combination with the direct term creates the mid-ranging function.
The transfer function from the output of the PI controller to the
153
Chapter 6. Enhanced Moisture Control Using the Air System
r2
r1 100
Ȉ PI Filter 1 Ȉ
Rm
Ȉ P1(s)
u1
Rm
100
1
100 u2 y
Ȉ
100 Rm
P2(s) Ȉ
1
Figure 6.11 The structure of the hybrid mid-ranging controller.
summation after the filter is in effect a low pass filtered derivative. The
purpose of the saturation block preceding the filter is to allow an offset in
u1 when u2 is reaching its limit. By this, the controller combines mid-
ranging (normal operation) with split-range (saturated u2) in a nice way.
The purpose of the other saturation block is to immediately engage control
signal u2 when u1 reaches its limit.
The PI controller can, without losing the mid-ranging function, be
replaced by any SISO controller.
Mid-range MPC
In [Allison and Isaksson, 1998] a mid-ranging controller is implemented
in the MPC (model predictive control) approach, see Figure 6.12. By
letting the minimized cost function be given by
H p 1 2 H u 1 2
ª r1 (k i ) yˆ (k i | k ) º ª 'u1 (k i ) º
J (k ) ¦
i 0
« r (k i) u (k i) »
¬ 2 1 ¼
¦
i 0
« 'u (k i ) » ,
¬ 2 ¼ R
(6.4)
Q
ªlarge 0 º ªsmall 0 º
Q « 0 , R , (6.5)
¬ small»¼ « 0
¬ large»¼
154
Chapter 6. Enhanced Moisture Control Using the Air System
r1 u1
P1(s)
r2
MPC
u2 y
P2(s) Ȉ
155
Chapter 6. Enhanced Moisture Control Using the Air System
’set point
air flow actuator’ ’set point steam pressure’
r2 u2
Ȉ C2(s) P2(s) ’moisture
disturbance’
CD (s)
’moisture’
1 n
r1 y
Ȉ C1(s) Ȉ
u1 P1(s) Ȉ Ȉ
function of the cascade system in the condensate system and for paper
quality reasons, as described in Section 2.3.
Since the paper drying process includes a significant transport dead-
time, the two controllers C1 and C2 are based on the IMC concept, which
is well known as an effective dead-time compensator for a stable process
with long time-delays [Morari and Zafiriou, 1989]. The transfer function
of an IMC can be written as
H f ( s) P ( s)
C ( s) , (6.6)
1 H f ( s) P ( s) P(s)
0.1365 40 s
P2 ( s ) e . (6.8)
61.5s 1
156
Chapter 6. Enhanced Moisture Control Using the Air System
9.5
8.5
8.0
7.5
60
Air flow actuator (%)
50
40
30
20
-50 0 50 100 150 200 250 300 350 400
Time (s)
Figure 6.14 Step response of the air flow í moisture process (solid) and the obtained
model (dotted), also given in (6.7). The dash-dotted line indicates the final value.
14
12
Moisture (%)
10
460
Pressure set point (kPa)
440
420
400
-50 0 50 100 150 200 250 300 350 400 450 500
Figure 6.15 Step response of the steam pressure í moisture process (solid) and the
obtained model (dotted), also given in (6.8).
157
Chapter 6. Enhanced Moisture Control Using the Air System
The model from air flow actuator to moisture, P1, has two complex
conjugated poles to capture the dynamics. The physical explanation for
the overshoot is a combination of both dynamics in the paper and
cylinder. When the dry air flow is decreased, the amount of moisture in
the air around the sheet is increased. This reduces the evaporation of water
from the paper and the paper moisture increases. This also decreases the
paper temperature which makes the energy flow to the paper to increase,
and this in turn will increase the evaporation. However, the dynamics in
the cylinder is much slower than the air-paper process and the increased
evaporation due to increased energy flow to the paper lags behind the
increase in sheet moisture due to reduced air flow, which gives an over-
shoot in the paper moisture.
Notice the resemblance between the prestudy experiment in Figure 6.7
and the open-loop simulation in Figure 6.14. The moisture in Figure 6.7
reaches a new steady-state promptly in a manner that does not appear as a
process with only one or two real poles. However, it is difficult to
distinguish the over-shoot in the response because of disturbances.
P1 has both faster dynamics and shorter time delay compared to P2.
The design parameter for the IMC is the filter Hf. For C1 it is chosen as
1
H f 1 (s) 2
, (6.9)
§ W1 ·
¨1 s ¸
© Z ¹
where IJ1 = 1 (nominal value) implies that the closed-loop poles are real
and placed at the same distance from the origin as the open-loop poles, see
Figure 6.16. A value of IJ1 greater than one makes the closed-loop system
slower, and vice versa. Apart from the two poles, P1 also has a zero.
However, it has been observed that a simpler model, with only two
complex conjugated poles and no zero, is sometimes sufficient for the
airflow í moisture process. Therefore, a second order filter is chosen in
(6.9), to make the controller C1 realizable in both cases.
Ideally, the only term affecting C2 is CD and the controller is tuned for
that process. This can be seen by inserting C2 into the structure in
Figure 6.13, assuming perfect process models. The purpose of C2 is to
slowly restore the signal u1 to its desired value. In practice, it can not be
assumed to have an ideal decoupling filter, and therefore it is important
that C1 and C2 do not interfere with each other. The closed-loop time
constant of C2 is therefore related to the fast loop and the IMC filter is
chosen as
158
Chapter 6. Enhanced Moisture Control Using the Air System
1
H f 2 ( s) . (6.10)
W
1 2 s
Z
The value of W2 is set to be significantly larger (at least five to ten times)
thanW1, to separate the controllers C1 and C2 in frequency.
The mid-ranging controller is compared to the case when only the
steam pressure in the cylinders are used to control the sheet moisture, here
denoted as steam pressure control. It is assume that the steam pressure
control is based on an IMC tuned for the process in (6.8) and the
corresponding filter is chosen as
1
H f 3 (s) , (6.11)
1 61.5W 3 s
Figure 6.16 The open-loop poles () and the closed-loop double pole (*) of the air flow í
sheet moisture process. The parameter Z is the distance between the open loop poles and
the origin.
159
Chapter 6. Enhanced Moisture Control Using the Air System
1
Compl. sens. function
0.8
0.6
0.4
0.2
0 -4 -3 -2 -1
10 10 10 10
1.5
Sensitivity function
0.5
0 -4 -3 -2 -1
10 10 10 10
Frequency (Hz)
Figure 6.17 Frequency plot for the mid-ranging system (dotted line) and steam pressure
control (solid line). The sensitivity function represents the transfer function from
disturbance n to sheet moisture y, see Figure 6.13. The complementary sensitivity
function is, for this control structure, equal to the transfer function from set point r to
output y.
160
Chapter 6. Enhanced Moisture Control Using the Air System
Remark
The air flow í moisture process can also be modeled with two real poles
and slow zero. The overshoot is then due to the lead action of the zero. In
this work, the process was identified by using System Identification
Toolbox in Matlab which gave a model with two complex conjugated
poles.
161
Chapter 6. Enhanced Moisture Control Using the Air System
10.5
Mid-ranging, W1 = 0.7, W2 = 5.0
Mid-ranging, W1 = 1.2, W2 = 5.0
10.0 Steam pressure control
9.5
Sheet moisture (%)
9.0
8.5
8.0
7.5
0 100 200 300 400 500 600
Time (min)
Figure 6.18 Comparison between the mid-ranging control system (W2 = 5.0) and steam
pressure control. The steam pressure control should be compared to mid-ranging with
W1 = 1.2.
10.5
Mid-ranging, W1 = 1.2, W2 = 5.0
Mid-ranging, W1 = 1.2, W2 = 10.0
10.0 Steam pressure control
Sheet moisture (%)
9.5
9.0
8.5
8.0
7.5
0 100 200 300 400 500 600
Time (s)
Figure 6.19 Moisture in paper for two different W2.
162
Chapter 6. Enhanced Moisture Control Using the Air System
Figure 6.19 shows that the response in sheet moisture is practically
independent ofW2. However, there is a slight difference at the second half
of the simulation. Larger values of W2 reduce the fluctuation since the
steam system is less aggressive in its attempt to ‘mid-range’ the air
system. Figure 6.20 and Figure 6.21 show the corresponding dew point of
the air and steam pressure in the lead group. It is evident that both the
steam and air process becomes less aggressive with larger W2. However,
the larger W2 is, the more the control system becomes ‘single-loop’ and the
advantages of mid-ranging are reduced.
Figure 6.21 also shows the steam pressure in the lead group for steam
pressure control. The steam pressure is not increased as rapidly by the
mid-ranging control as by steam pressure control. In the mid-ranging case,
the fast moisture transients are handled by the air system and the steam
system is only used to restore the air system in steady-state. Less
variations in steam pressure is advantageous since it reduces the injection
of disturbances in the steam and condensate system, which has negative
effect on both steam production and other steam users, see also
Section 2.4.
Remark 1
One drawback of using the dry air flow (and thus the dew point) to control
the sheet moisture, is the risk of reduced efficiency in the heat recovery. A
well optimized drying section ventilation has a dew point close to the
allowed maximum. For the proposed control technique to function well,
the average dew point then needs to be decreased which leads to higher
energy costs. This should be taken into account by weighing the gain of
reduced variability in sheet moisture against increased energy usage
[Lindell and Stenström, 2004], when evaluating the control principle for a
specific drying section.
Remark 2
Physically, the dew point is a driving force for the evaporation of water in
the sheet. A low dew point implies a low vapor partial pressure in the air,
and high difference in vapor pressure between sheet and air. It can be
interpreted as a low dew point pulls out the moisture in the sheet. This is
the opposite of an increase in steam pressure which increases the
evaporation by increasing the vapor pressure in the sheet, and pushes out
the moisture in the sheet. In [Karlsson and Stenström, 2005b] it is shown
that a high vapor pressure inside the sheet can cause delamination
163
Chapter 6. Enhanced Moisture Control Using the Air System
63
62
Dew point (°C)
61
60
59
W1 = 1.2, W2 = 5.0
W1 = 1.2, W2 = 10.0
58
0 100 200 300 400 500 600
Time (s)
Figure 6.20 Change in dew point (mid-ranging).
465
460
Steam Pressure set point (kPa)
455
450
445
164
Chapter 6. Enhanced Moisture Control Using the Air System
(different pulp layers are separated) problems in board machines. Since
the mid-ranging control presented here, gives smaller variations in steam
pressure it should reduce this problem.
Remark 3
An alternative control configuration is to let the dew point be the inner
part of a cascade loop with the moisture control. Controller C1 then gives
a set point to the dew point controller, which is included in P1, compare
with Figure 6.13. However, the dew point is normally measured in the
exhaust air, far from the sheet. Experiments in [Forsman and Birgerson,
1999] show that the time constant in the process from air flow actuator to
dew point is around 1.5 minutes due to the large volume of air in the
hood. This would give a slow moisture control loop, unless the dew point
is measured close to the paper (in the pocket). The experiment in
Figure 6.7 gives a time constant for the actuator í dew point process
around 30 seconds. Also in this case the dew point needs to be measured
closer to the sheet to be useful in a cascade control configuration.
2 2.5 m
Figure 6.22 The air balance in the hood. The large unfilled arrows indicate the exhaust air.
To provide uniform air flows around the cylinders, there is a false ceiling. The smaller
arrows indicate the pressure inside the hood. Below the doors of the hood, there is an under
pressure to prevent moist air to leak into the machine room. The height where the air
pressure inside the hood equals the outside pressure, is called zero-level, [Karlsson, 2000].
165
Chapter 6. Enhanced Moisture Control Using the Air System
C3(s) Ȉ
r3
r2 u2
> Ȉ C2(s) P2(s)
CD (s)
1
1
r1 y
Ȉ C1(s) Ȉ
u1 P2(s) Ȉ
1
Figure 6.23 By allowing the dew point to adjust the air flow actuator set point, it is
possible to set an upper constraint on the dew point. The output from the selector is fed
back to the integral mode of C3 to avoid windup.
air is used to control the zero-level, and the exhaust air to control the dew
point [Forsman and Birgerson, 1999]. The dew point is often measured in
the exhaust air channel.
The proposed control strategy in this chapter uses the supply air to
control the sheet moisture. That means that the exhaust air must be used to
control the zero-level. It is vital to have an upper constraint on the dew
point to prevent condensation that might cause dripping on the sheet and
corrosion on machine units. A possible solution, in the mid-ranging
structure, is shown in Figure 6.23. By adding an extra controller, C3, that
controls the dew point in the exhaust air through a selector that is shared
with set point r2, the air flow actuator set point is increased if the dew
point exceeds set point r3.
6.6 Summary
A new control strategy has been developed for the multi-cylinder drying
section, to reduce sheet moisture variability. The air system in the blow
boxes are used to control sheet moisture, in combination with the steam
system. A previously developed physical model for the drying section is
used in simulations, where the performance of the control system is
compared with a conventional moisture control.
166
Chapter 6. Enhanced Moisture Control Using the Air System
The result is that by combining the conventional steam control loop with
the supply air system, it is possible to significantly improve the
disturbance rejection of the sheet moisture feedback loop, without
changing the level of robustness. This is shown in both in simulations and
by observing that the bandwidth for the disturbance rejection is more than
twice as large for the mid-ranging control. The physical reason for this is
the fast response of the dewatering rate to changes in air dew point.
There are two design parameters in the selected mid-ranging structure,
W1 and W2. The performance of the closed loop system is determined by W1
and W2 determines how fast the steam system should restore the air system
to its original level. Ideally, these two parameters can be chosen
independently but in practice it is recommended to let W2 5W1. Letting
W1 = 1, means that the closed loop poles have the same distance to the
origin as the open loop poles.
167
Chapter 7. Feedforward from a Paper Surface Temperature Measurement
7
7.
The paper machine is often run with constant set points for all quality
variables for a long period of time, typically several hours or even days.
The primary function of the control system is then to maintain the desired
properties in the presence of disturbances. In Chapter 6, it is shown that
the air in the surrounding hood can beneficially be integrated with the
steam pressure control system by mid-ranging. The mid-ranging control
system proved to be far more efficient than using only the steam pressure
as manipulated variable. By using a combination of both the fast but low-
capacity air system and the slow but high-capacity steam system, the mid-
ranging control system was able to significantly reduce the effect of
moisture disturbances.
Due to the long time constants and time delays in the process, there is a
large potential gain in disturbance rejection by introducing feedforward to
the control system. It is essential that the disturbance signal is measured
early in the process to make the feedforward effective. Installation of a
moisture scanner between the press and drying section is possible and the
moisture measurement can be used for feedforward control. However, the
space between the press and drying section is often very small and a
scanner is also a fairly large investment for the mill.
168
Chapter 7. Feedforward from a Paper Surface Temperature Measurement
Therefore, a new signal to be used for feedforward control of the paper
moisture content is presented. The signal is based on the paper surface
temperature in a few positions in the machine direction. The drying
process generates large moisture gradients in the thickness direction of the
paper and the surface becomes dry before the center of the sheet. When
the paper surface becomes dry the surface temperature quickly increases
to a value above 100°C. The general idea is therefore to estimate the
position where the paper surface becomes dry by measuring the surface
temperature. The possibility to measure the signal is verified by
measurements on a real paper machine. The control system is evaluated
by simulations of the drying section model presented in Chapter 6. Similar
work for the wire section, where new feedforward signals have been
proposed, can be found in [Larsson and Gustafsson, 1998] and [Li, et al,
2001].
70
Inlet moisture content: 61.39 %
Inlet moisture content: 62.12 %
60
50
Moisture content (%)
40
Node 2
30
20
Node 1
10
0
0 10 20
50 30 40
60 70 80 90 100
Cylinder no.
Figure 7.1 Moisture content at the first two nodes in the thickness direction. Note that
when the first node reaches a moisture content of 10 %, the second node is still at 60 %.
169
Chapter 7. Feedforward from a Paper Surface Temperature Measurement
free water is equal to the latent heat of vaporization. The bound water,
however, is tightly bound to the fibres by hydrogen bonds, and requires
more energy to evaporate compared to the free water. This means that the
paper is more difficult to dry in presence of only bound water. The
proposed feedforward structure uses this indirectly by measuring the
surface temperature.
As described in Section 6.2, the discretization of the paper sheet model
in the thickness direction is set to 10 nodes. A steady-state simulation of
the moisture content at the two first nodes seen from the lower surface, is
shown in Figure 7.1. The first layer reaches the bound water region at
cylinder number 25. This can be seen since the evaporation rate changes
significantly. The same event happens to the second fibre layer at cylinder
number 43. The paper surface then contains a low moisture content at the
same time as a lot of free water is still present in the middle of the paper
structure.
To validate the model in terms of simulated paper surface temperature,
the temperature has been measured on the paper machine the model is
adapted for. By using an IR-camera which gives the surface temperature
in a two-dimensional window, temperature values along a line in the
130
120
110
Paper temperature (°C)
100
90
80 115
70 110
60 105 Simulated
100 Measured
50
0 20 40 60 80 100
Cylinder no. 95
90
85
80
Simulated
75
Measured
15 20 25 30 35 40 45 50
Cylinder no.
Figure 7.2 Comparison between measured and simulated paper surface temperature. The
most important section for the feedforward structure in enlarged and the 100°C-level is
indicated.
170
Chapter 7. Feedforward from a Paper Surface Temperature Measurement
Steam
Condensate (a) (b)
p1 p1
Cylinder shell
Different surface
temperatures.
Figure 7.3 Schematic picture that shows the changed properties as the paper surface dries
out. A wet surface (a) will act differently compared to a dry surface (b), given the same
external conditions.
machine direction are singled out and the results are shown in Figure 7.2.
The simulations show both the temperature for the free draw and the part
where the paper is in contact with the cylinder. Since the paper in the
contact zone is not visible due to the dryer fabrics, the temperature is only
measured in the free draw. The simulated values show very good
agreement with the measured values. Observe that this is a pure validation
without any adjustment of the model, at an operating point different from
the one where the model was fitted. Paper temperatures exceeding 100°C
are first observed in the free draw section after cylinder number 27 and
this phenomenon appears more frequently after cylinder number 35.
The layer of surface fibres without any surrounding free water show
totally different properties when it comes to both evaporation and heat
transfer mechanisms, compared to fibres surrounded by free water. Since
the evaporation requires a large amount of energy, the absence of free
water will cause the paper temperature to increase. This can also be
considered from another viewpoint. If the surface temperature is above
100°C, there can not exist any free water on the surface since it would
instantly evaporate to the surrounding air. The important event is thus the
formation of a dry fibre layer since the paper surface temperature will
change as the free water is removed. An illustration of the phenomenon is
shown in Figure 7.3.
171
Chapter 7. Feedforward from a Paper Surface Temperature Measurement
(a) (b)
Suggested measuring
position
172
Chapter 7. Feedforward from a Paper Surface Temperature Measurement
150 300
140 Group 4
130
Constant 250
Paper surface temperature (°C)
120 temperature
zone
90
80
Sudden increase 150
Group 2
70 in temperature
Constant
temperature
60 zone
50 100
5 10 15 20 25 30 35 40 45 50
Cylinder no.
Figure 7.5 The temperature curve obtained by measuring the surface temperature as
indicated in Figure 7.4. The steam pressures in the different groups are also shown.
w in = 62.12 %
5
w in = 61.39 %
Free draw temperature difference (°C)
Changed
-1 peak position
5 10 15 20 25 30 35 40 45
Cylinder no.
Figure 7.6 A temperature curve calculated by taking the difference between two adjacent
temperature measurements, as indicated in Figure 7.4.
173
Chapter 7. Feedforward from a Paper Surface Temperature Measurement
two adjacent measurements are shown in Figure 7.6. All the differences in
free draw paper temperature that deviate from zero are either due to a new
steam group or because the paper surface has dried out. The paper surface
temperature differences were fitted to cubic splines and the position the
maximum of these splines (originated from a dried-out surface) will
hereafter be called “peak position”.
Two different cases with different inlet paper moisture ratios are
shown in Figure 7.6. The paper temperature is affected differently at the
two positions where the two first fiber layers dry out but shows the same
behavior for the positions of all other peaks (representing the heat-up
periods). The peak at cylinder number 43 is a reaction of the free water
removal in the second node and does not show large variation in position
(cylinder number).
Note that there might well be a difference in peak position and the
position where the paper surface actually dries out. The maximum of the
temperature difference shows where the dry surface causes maximum
effect and not the exact location of a dry surface. For instance, if position
(a) in Figure 7.4 show 1 °C higher temperature than the previous
measuring position and position (b) show 6 °C higher temperature than
(a), the surface may still have dried out just before (a). The effect of the
dry surface is however radically increased as the paper reaches (b). The
physical meaning of the signals, however, is of minor importance as long
as it indicates changes in paper moisture. Both the peak position and the
sheet moisture at the reel up will be affected by a disturbance in inlet
moisture content and the peak position can therefore be used in
feedforward control.
To measure the peak position, an IR camera should be positioned in
the free draws around the position where the surface becomes dry. It is
crucial to measure at the same position in each free draw to obtain
relevant comparisons of the temperatures. Here, the chosen position is the
middle of the free draw. The peak position may however be difficult to
measure online with the current available measuring technique. The
measuring devices must e.g. withstand the humid environment inside the
hood and problems with fibre coating on the device must be dealt with.
Also, the measuring depth could cause problems if several fibre layers are
visible by the measuring equipment. During the experiments it was also
observed that both the viewing angle and distance to the paper is
important for the result.
174
Chapter 7. Feedforward from a Paper Surface Temperature Measurement
r2 u2
Ȉ C2(z) P2(z) ”peak position”
10
7
0 50 100 150 200 250 300
26.5
Peak position (cyl.)
26.0
25.5
25.0
0 50 100 150 200 250 300
Time (s)
Figure 7.8 The experiment used for identification of H. The solid lines are the result of the
simulation and the dotted line is the obtained simplified model given in (7.1).
175
Chapter 7. Feedforward from a Paper Surface Temperature Measurement
”Temperature
peak position”
Hff,3
Ȉ P1,3
H
Hff,4
u1 y
Ȉ P1,4 Ȉ Ȉ Ȉ
Hff,5
Ȉ P1,5
Figure 7.9 A block diagram of the feedforward structure used in the simulations.
0.185 z 6
H ( z) . (7.1)
1 1.528 z 1 0.599 z 2
176
Chapter 7. Feedforward from a Paper Surface Temperature Measurement
9.0
8.0
Air flow actuator (%)
50
40
30
20
0 50 100 150 200 250 300 350 400
Time (s)
Figure 7.10 Step responses from the five air flow actuators to moisture content at the reel
up. The two dashed lines are the first two air groups which are not used for the
feedforward control.
H H H
H ff ,3 , H ff , 4 , H ff ,5 . (7.3)
3P1,3 3P1, 4 3P1,5
The control signals of the last three air groups now might be different and
therefore C2 was set to control the actuator signal of the first two groups
(u1). Another possible solution is to let C2 control the mean value of all
177
Chapter 7. Feedforward from a Paper Surface Temperature Measurement
five air flow actuators or to introduce multivariable control in the steam
system so that each air group is mid-ranged by a steam group.
Just as in the case of H, the air groups are identified from step responses,
shown in Figure 7.10, and are given by
0.00411 0.00198 z 1 6
P1,3 ( z ) z ,
1 1.609 z 1 0.685 z 2
0.00476 0.00224 z 1 4
P1, 4 ( z ) z , (7.4)
1 1.490 z 1 0.590 z 2
0.00537 0.00214 z 1 1
P1,5 ( z ) z .
1 1.537 z 1 0.639 z 2
7.3 Simulations
11.0
10.5
10.0
No control - open loop
Steam pressure control
Sheet moisture (%)
Mid-ranging only
9.5 Feedforward + mid-ranging
9.0
8.5
8.0
7.5
0 100 200 300 400 500 600
Time (s)
Figure 7.11 Comparison of different control systems.
178
Chapter 7. Feedforward from a Paper Surface Temperature Measurement
26.5
25.5
25.0
8.5
Sheet moisture (%)
8.0
7.5
8.5
8.0
7.5
0 100 200 300 400 500 600
Time (s)
Figure 7.12 Peak position and sheet moisture during the closed loop simulation.
63
62.5
62
61.5
61
Dew point ( C)
o
60.5
60
59.5
59
Mid-ranging only
Feedforward + mid-ranging
58.5
0 100 200 300 400 500 600
Time (s)
Figure 7.13 Change in dew point.
179
Chapter 7. Feedforward from a Paper Surface Temperature Measurement
465
460
Steam pressure (kPa)
455
450
445
440
Steam pressure control
Mid-ranging only
Feedforward + mid-ranging
435
0 100 200 300 400 500 600
Time (s)
Figure 7.14 Comparison of the steam pressure in the lead group by the different control
system.
0.8
Steam pressure control
Mid-ranging only
0.7 Feedforward + mid-ranging
0.6
0.5
Amplitude
0.4
0.3
0.2
0.1
0 -3 -2 -1 0
10 10 10 10
Frequency (rad/s)
Figure 7.15 Frequency plot, showing how the disturbances in inlet moisture propagate, for
the different control systems.
180
Chapter 7. Feedforward from a Paper Surface Temperature Measurement
The most apparent result is that the disturbance is greatly attenuated. After
150 seconds the decreasing rate of the moisture content restoration is
reduced for the feedforward system. The reason for this can be seen in
Figure 7.12, where the peak position is shown for the same simulation.
When the control system compensates for the disturbance by increasing
the control signal, it also slightly affects the peak position. The
feedforward signal is thus affected but this dynamic relation from u1 to the
peak position in Figure 7.7 is neglected. By considering this in (7.2) the
feedforward could probably be slightly improved. For comparative
reasons Figure 7.11 also shows the moisture deviation when all control
systems are switched off (open loop).
The dew points in P1,3 are shown in Figure 7.13 and the steam pressures
are shown in Figure 7.14. Notable for the two feedforward control
systems is that the air flow actuators are not restored to their original
positions by the mid-ranging part. This was also noted in Section 7.2.
Multivariable control of the air and steam groups could be introduced to
give maximum manoeuvring range to the individual air flow actuators, if
desired.
The different control systems are also compared by evaluating their
disturbance rejection in the frequency domain, see Figure 7.15. The figure
is generated by identifying the process from inlet moisture content to
moisture content at the reel-up and the figure simply shows the Bode plot
of the identified model. The feedforward control shows to be superior.
Due to the low pass characteristics of the process, all systems attenuate
high frequencies similarly. For low frequencies the difference between the
systems is considerable. Since there are both slow dynamics and long time
delays in the process, the feedback systems can only attenuate extremely
slow disturbances (note that 10í3 rad/s is a variation with a period time of
almost two hours). The feedforward can beneficially use its prior
information about the disturbances acting on the system.
The feedforward assumes linear systems while the system is nonlinear.
With the presented technique to obtain the control parameters it is
therefore not possible to completely eliminate the effect of disturbances.
Remark
The implementation of the feedforward does not require a mid-ranging
control system. It might as well be combined with any other moisture
control structure. The mid-ranging controller is chosen here to give
natural connection to Chapter 6.
181
Chapter 7. Feedforward from a Paper Surface Temperature Measurement
7.4 Summary
A new measured signal, called the peak position, has been presented that
can be used for feedforward in moisture control. The signal is based on
measurements of paper surface temperature in several free draw zones. By
observing the position at which the paper dries out on the surface it is
possible to compensate for moisture disturbances before they are
measured by the moisture scanner at the reel-up.
With the addition of the peak position in combination with a process
model, it is possible to significantly reduce the transients from inlet
moisture disturbances. From the studied step disturbance it is clear that the
feedforward signal greatly improves the performance of the paper
moisture control system.
The peak position may be difficult to measure with the current
available technique. The potential increase in performance is, however,
large and the signal may be replaced by any other suitable signal that
indicates disturbances in moisture content. Equipment that measure the
surface moisture content are being developed [Boström, 2002]. This
signal would be very suitable for feedforward control. What type of signal
that is finally chosen is of minor importance, as long as it can be used as a
disturbance indication.
182
Chapter 8. Object-Oriented Modeling and Control of the Moisture
88.
Object-Oriented Modeling
and Predictive Control of the
Moisture Content
d
U sV s q s q c q bt ,
dt (8.1)
d
U wV w q c q w .
dt
184
Chapter 8. Object-Oriented Modeling and Control of the Moisture
The energy balances for steam, water, and metal are
d
U s u sV s (q s q bt )hs q c hs ,
dt
d
U w u wV w q c h s q w h w Q m , (8.2)
dt
d
dt
mC p , m Tm Q m Q p ,
where Qm is the power supplied from the water to the metal, Qp is the
power supplied from the metal to the paper, hs is the steam enthalpy, hw is
the water enthalpy, m the mass of the cylinder shell, Cp,m the specific heat
capacity of the shell, Tm the mean temperature of the metal, us and uw are
the specific internal energies of steam and water. The steam and water
volumes add up to the total cylinder volume,
V Vs Vw . (8.3)
Qm D sc Acyl Ts Tm , (8.4)
Qp D cp AcylK Tm T p , (8.5)
where Tp is the paper temperature, Dcp the heat transfer coefficient from
the cylinder shell to the paper, and Ș is the fraction of dryer surface
covered by the paper web. In Chapter 4, Dsc is used to calibrate the model
against plant data. Here, both Dcp and Dsc are possible candidates for that
purpose. Experiments have shown that Dsc depends on both condensate
thickness, machine speed, and the number of spoiler bars, see [Pulkowski
and Wedel, 1988] and [Wilhelmsson, 1995]. However, the condensate has
a turbulent behavior and the heat transfer coefficient has proven to be
difficult to model, see Figure 8.1. Therefore Dsc is used as a free variable
to calibrate the model with. Empirical models for Dcp have been
185
Chapter 8. Object-Oriented Modeling and Control of the Moisture
15 spoiler bars
3000
25 spoiler bars
30 spoiler bars
No spoiler bars
2500
Heat transfer coefficient (W/m K)
2
2000
1500
1000
500
0
0 2 4 6 8 10 12 14 16
Condensate thickness (mm)
Figure 8.1 Steady-state measurements of how the number of spoiler bars affect the heat
transfer coefficient for the condensate [Pulkowski and Wedel, 1988]. It also depends on
machine speed and bar size, and is difficult to model due to the turbulent behavior.
is obtained, where Dcp(0) varies between 200 and 500 W/(m2K). It is well
known that Dcp depends on other things, e.g. the web tension, and surface
smoothness of both paper and cylinder, but this is omitted here.
Since the steam flow to the cylinder cannot be manipulated directly, a
valve model is also needed. From [Thomas, 1999] we have
qs Cv f v ( xv ) ( psh p) U s , (8.7)
186
Chapter 8. Object-Oriented Modeling and Control of the Moisture
1.0
R v = 200
R v = 50
0.8
R = 10
v
fv ( xv ) 0.6
0.4
0.2
0
0 0.2 0.4 xv 0.6 0.8 1.0
Figure 8.2 Equal percentage valve characteristic.
x v 1
f v ( xv ) Rv . (8.8)
Equations (8.1) í (8.9) are a crude nonlinear model for the steam-cylinder
process. By choosing p, Vw, and Tm as state variables and using partial
187
Chapter 8. Object-Oriented Modeling and Control of the Moisture
derivatives, the system can be rewritten into a third order state equation
(most steps are omitted here), where f1, f2, and f3 are defined as the right
hand side of the equations.
dp dV
e11 e12 w Cv f v ( xv ) ( psh p) U s qw qbt f1 ,
dt dt
dp dVw
e21 e22 Cv f v ( xv ) ( psh p ) U s hs (8.10)
dt dt
qbt hs qw hw D sc Acyl (Ts Tm ) f 2 ,
dTm
e33 D sc Acyl (Ts Tm ) (D cp (0) 955u ) AcylK (Tm T p ) f3 ,
dt
where
dU s dU
e11 (V Vw ) Vw w ,
dp dp
e12 Uw Us ,
dU s dh dU dh
e21 hs (V Vw ) U s (V Vw ) s hwVw w U wVw w V ,
dp dp dp dp
e22 U w hw U s hs ,
e33 mC p .
(8.11)
In the rewritings of the energy balances above the specific internal energy
has been eliminated by the definitions us = hs p/ȡs and uw = hw p/ȡw.
Using f1, f2, and f3 the system can be further rewritten into an explicit state
form.
188
Chapter 8. Object-Oriented Modeling and Control of the Moisture
Figure 8.3 A Simulink model of (8.12). The submodel f2 is also opened, which in turn
contains submodels.
model Cylinder
equation
Ms = rhos*Vs;
der(Ms) = qs – qc - qbt; % (8.1)
Mw = rhow*Vw;
der(Mw) = qc - qw;
Es = rhos*us*Vs;
us = hs – p/rhos;
der(Es) = (qs – qbt)*hs – qc*hs;
Ew = rhow*uw*Vw; % (8.2)
uw = hw – p/rhow;
der(Ew) = qc*hs – qw*hw – Qm;
Em = m*Cp*Tm;
der(Em) = Qm - Qp;
V = Vs + Vw; % (8.3)
Qm = alpha_sc*Acyl*(Ts - Tm); % (8.4)
Qp = alpha_cp*Acyl*eta*(Tm - Tp); % (8.5)
end Cylinder;
Figure 8.4 Code segment of the Modelica model of the steam cylinder. In addition,
equations for steam properties are required to give a complete simulation model.
189
Chapter 8. Object-Oriented Modeling and Control of the Moisture
In this form the model can be directly implemented and simulated in e.g.
Simulink, see Figure 8.3. By using Modelica instead the tedious and error
prone procedure of transforming the system to explicit form is avoided
and we let the simulation environment decide the state realization. Since
the transformation of equations is automated, it is also easier to change the
model at a later stage. Equations (8.1) í (8.5) are put into the simulation
environment as they are, see Figure 8.4.
ptot KM w §¨ ptot pv , a ·¸
qevap ln , (8.13)
Rg T p ¨ ptot pv , p ¸
© ¹
x
p v,a p tot . (8.14)
x 0.62
p v, p M pv0 , (8.15)
where pv0 is the partial vapor pressure for free water. This is given by
Antoine’s equation
190
Chapter 8. Object-Oriented Modeling and Control of the Moisture
1500
1.0 o
0.5
30 C
o 500
o
60 C
o
90 C
0 0
0 5 10 15 20 25 0 5 10 15 20 25
Moisture content (%) Moisture content (%)
Figure 8.5 Sorption isotherm, ij, and heat of sorption, ǻHs given by (8.17) and (8.21).
§ ·
¨ 10.127 1690 ¸
¨ T 43.15 ¸
© p ¹
p v 0 10 . (8.16)
As long as capillary transport can bring new water to the paper surface,
the vapor partial pressure at the paper surface is equal to the partial
pressure for free water. When the paper becomes dryer a correction factor
called sorption isotherm, ij, is invoked which has a value between zero
and one, see Figure 8.5. The sorption isotherm of a paper web depends on
its composition and temperature. It is not very well investigated when
compared to other materials [Pettersson and Stenström, 2000], but
[Heikkilä, 1993] gives an empirical expression for paper pulp, namely
where u is the moisture ratio (kg moisture/kg fiber). Also, let vx be the
speed of the paper web (m/s), dy the width of the paper web (m), Axy the
area of the dryer surface covered by paper (m2), and g the dry basis weight
(kg/m2). Then the mass balance of moisture in the paper web can be
written as
d (ugAxy )
d y v x guin Axy qevap d y v x gu. (8.18)
dt
Figure 8.6 shows a schematic picture of the mass flows in the model. To
model the energy balance, introduce
191
Chapter 8. Object-Oriented Modeling and Control of the Moisture
Axy qevap
vx
d y v x guin d y v x gu
Paper web
Figure 8.6 The mass balance for moisture in the paper web. The shaded area is the
cylinder wall
C p , fiber uC p , w
C p, p , (8.19)
1 u
where Cp,p, Cp,fiber, and Cp,w is the specific heat capacity for the paper, fiber
and water, respectively (J/kg·K). As we can see, Cp,p is a weighted sum of
the heat capacities of the parts. From [Wilhelmsson, 1995] we have Cp,fiber
= 1256 J/(kg·K).
Also, let Tp be the paper temperature and ǻH be the amount of energy
needed to evaporate the water. Analogously to the discussion about the
mass balance, if the web is wet enough this energy is equal to the latent
heat of vaporization for free water. When the paper becomes dryer an
extra amount of energy ǻHs (the heat of sorption) is necessary besides the
latent heat of vaporization for free water. The heat of sorption can be
derived from the sorption isotherm by thermodynamic theory and this
relation is known as the law of Clausius-Clapeyron
R g ª d (ln M ) º
'H s « », (8.20)
M w «¬ d (1 / T p ) »¼
M 1
'H s 0.10085u 1.0585T p2 R g . (8.21)
M wM
The amount of energy required to evaporate the water from the surface of
the web is the given by
192
Chapter 8. Object-Oriented Modeling and Control of the Moisture
Axy qevap ( ΔH vap + ΔH s )
d y v x g (1 + u in ) C p , p T p ,in d y v x g (1 + u )C p , pT p
Paper web
Qp
Figure 8.7 The energy balance of the paper web. The shaded area is the cylinder wall.
where ΔHvap is the latent heat of vaporization, equal to 2260 kJ/kg (at
atmospheric pressure).
Reference [Pettersson and Stenström, 2000] investigates some sorption
isotherms found in the literature. Many of those give a heat of sorption
that goes to infinity as u goes to zero. This is physically unrealistic since
the bond energy between the last fraction of water and a cellulose fiber
must be finite. From [Heikkilä, 1993], a finite heat of sorption at the
origin which matches the hydrogen bond energy between water−fiber is
given and is therefore found to be most appropriate. The heat loss in the
paper due to mass evaporation is dominating the heat conduction and
radiation, which therefore can be neglected. In addition, since water is an
incompressible medium there is no pressure volume work on the
surroundings and we write the energy balance as a change in enthalpy.
The energy balance of the paper web is thus modeled as
d ( g (u + 1) Axy C p , p T p )
= Q p + d y v x g (1 + u in )C p , p T p ,in
dt (8.23)
− Axy q evap (ΔH vap + ΔH s ) − d y v x g (1 + u )C p , p T p ,
see Figure 8.7. In addition, we let the heat transfer coefficient from the
cylinder to the paper web depend on the moisture content in the web.
To summarize the complete drying section model, it is given by the
balances in (8.1), (8.2), (8.18), and (8.23), together with the algebraic
relations given in (8.3)−(8.9), (8.13)−(8.17), (8.19), (8.21), and (8.22).
Equation (8.1) and (8.2) define the dynamics of a cylinder, and one set of
these equations are needed for each cylinder in the drying section. The
dynamics of a lumped paper web model is given by (8.18) and (8.23). By
connecting a series of these equations a discretized model is obtained,
193
Chapter 8. Object-Oriented Modeling and Control of the Moisture
60
50
40
Sheet moisture (%)
30
20
10
0
5 10 15 20 25 30 35
Cylinder number
Figure 8.8 Validation by comparing steady-state simulations with moisture
measurements on a fine paper machine.
where the outflow of paper of one component becomes the inflow of the
next component.
The two partial models in Section 8.1, are validated separately in [Slätteke
and Åström, 2005] and [Persson, 1998]. To assure that the combined
model is satisfying it has been validated against steady-state data taken
from a paper machine producing fine paper [Stenström, et al, 1994]. The
machine is running at 708 m/min with a basis weight of 80 g/m2. The
paper is over dried to a final moisture content of only 0.012 kg/kg, since
this is a predryer which is immediately followed by a size press (a unit
where starch is applied to the surface to obtain strength and water
resistance) where a certain amount of rewetting occurs.
The moisture content is measured in 11 different positions, which can
be seen in Figure 8.8 together with the simulation result. The heat transfer
coefficient Įsc is used as fitting parameter with the same value in all
cylinders and 1100 W/(m2·K) is found to give the best fit by visual
inspection. The agreement between model and measurements is good even
194
Chapter 8. Object-Oriented Modeling and Control of the Moisture
Pres s
Pope
Header pressure
HeaderPres...
F
P
k={1000e3}
ID
PID i P
Figure 8.9 The Modelica model of the drying section used in the simulations. The
component of the second dryer group is expanded, showing the cylinders, steam valve, and
PID controller. It is assumed that the temperature and moisture are constant at a single
cylinder due to the high machine speed, hence modeled as one control volume.
195
Chapter 8. Object-Oriented Modeling and Control of the Moisture
8.3 Open loop simulations
550
500
D cp = 500
450
Steam pressure (kPa)
150
140
D cp = 4000
130
75
Valve position (%)
70
65
60
55
100 150 200
Time (s)
Figure 8.10 Pressure response for a step in valve position for two different values of Dcp.
Heat transfer coefficient Dsc = 4000.
196
Chapter 8. Object-Oriented Modeling and Control of the Moisture
0.20
qs
q c+ q bt
0.16
510
Steam pressure (kPa)
500
490
480
470
460
90 100 110 120 130 140 150 160 170 180 190 200
Time (s)
Figure 8.11 Open loop step response in pressure with Dsc = 4000 W/(m2·K). Heat transfer
coefficient Dcp is set to a constant value 500 W/(m2·K).
450
Q p, Q m (kW)
Qm
400 Qp
350
250
Energy derivatives (kJ/s)
d( U su sV s )/dt
200
d( mC p,mT m )/dt
150 d( U wu wV w )/dt
100
50
0
-50
100 150 200 250
Time (s)
Figure 8.12 Open loop step response in pressure with Dsc = 4000 W/(m2·K). Heat transfer
coefficient Dcp is set to a constant value 500 W/(m2·K).
Dcp = 500 W/(m2·K), the process can be modeled with an integrator plus
some dynamics for the given time frame (in reality we know that the
integrator is a slow pole). However, for the response with Dcp = 4000
W/(m2·K) an IPZ model is not sufficient. Note that since the heat transfer
coefficients are different in the two figures, the steady-state solutions are
197
Chapter 8. Object-Oriented Modeling and Control of the Moisture
also different and the two simulations start from different steady-state
steam pressures. The conclusion is unaffected if the two responses start at
the same pressure.
Figure 8.11 gives an explanation for the characteristic appearance of
the open loop step response of the steam pressure. When the steam valve
position is increased, it will increase the steam pressure inside the cylinder
and consequently the saturation temperature of the steam and the
condensation rate. However, the increasing condensation rate lags behind
the increasing steam inlet flow as the condensate layer heats up to the new
steady state temperature. Therefore, there will be a fast initial build-up in
steam pressure, before the steam consumption has reached its new value.
This is also noted in Chapter 3.
Figure 8.12 shows the energy flows to the cylinder and paper, and also
the energy derivatives. In the derivation of the linear model in Chapter 4,
it is simplified by assuming that the thermal dynamics of the water is fast
compared to the dynamics of steam and metal. The energy balance is
therefore replaced by a static model by defining its derivative as zero.
9.5
Sheet moisture (%)
9.0
8.5
325
Steam pressure (kPa)
320
315
310
305
300
80
Valve position (%)
60
40
20
0
0 100 200 300 400 500 600 700 800 900 1000
Time (s)
Figure 8.13 Simulation of set point change in steam pressure and a web break.
198
Chapter 8. Object-Oriented Modeling and Control of the Moisture
Figure 8.12 shows that this assumption is valid. Since the cylinder shell
acts as a thermal filter Qp is increasing more slowly that Qm.
Figure 8.13 shows a simulation where the steam pressure controllers
are in closed loop while the moisture is uncontrolled. The pressure is
controlled by a PI controller tuned according to the method in Chapter 5
with Ms = 1.2. At t = 100 s there is a steam pressure set point change. The
valve position returns to its initial position when the set point is reached
due to the integral action in the process. Note that the characteristic over-
shoot in steam pressure does not result in an over-shoot in the moisture.
At t = 500 s there is a web break and the moisture quickly reaches a level
around 60 % since the sheet has lost the energy source for its evaporation,
and the curve is out of the scale in the figure. Since the main energy sink
is lost for the cylinders, the control valve is almost closed to keep the
desired steam pressure.
Figure 8.14 shows two simulations of moisture response for a set point
change in steam pressure: One with a medium heat transfer coefficient
and one with a high heat transfer coefficient. Modeling the two responses
as a first order process, they are given by
0.109 16.1s
P(s) e D sc 2000 W/m 2 K (8.24)
46.2 s 1
and
0.0917 16.6 s
P( s) e D sc 4000 W/m 2 K (8.25)
21.2s 1
In the last decade, MPC (model predictive control) has received large
attention in the process industry. It has been described as one of the most
199
Chapter 8. Object-Oriented Modeling and Control of the Moisture
8.6
8.2
8.0
Dsc = 2000 W/m2K
7.8
6.2
Sheet moisture (%)
6.0
5.8
385
380
375
370
450 500 550 600 650 700 750 800
Time (s)
Figure 8.14 Response in moisture to step change in steam pressure set point for two
different values of heat transfer coefficient Dcp.
200
Chapter 8. Object-Oriented Modeling and Control of the Moisture
Setpoint
moisture PID- Steam system
uc1 controller First part
r1
MPC Dryer
Moisture
y
where y is the moisture content, uc1 is the steam pressure set point to the
first 50 cylinders, and uc2 the set point to the last 10 cylinders. The PID-
controllers in the inner loops are tuned according to the IPZ-tuning
201
Chapter 8. Object-Oriented Modeling and Control of the Moisture
derived in Chapter 5. The response from uc2 has a significantly shorter
time delay but also smaller process gain. The advantage of the high gain
from uc1 and fast dynamics from uc2 is utilized in a mid-ranging MPC
structure.
The simplified model for the singe-loop case in Figure 8.15, obtained
correspondingly to (8.26), is given by
0.11 13s
Y e Uc , (8.27)
48s 1
2 2
49 3
ª r1 (k i ) yˆ (k i | k )º ª 'uc1 (k i ) º
J (k ) ¦
i 0
« r (k i) u (k i) »
¬ 2 c2 ¼
¦
i 0
« 'u (k i )» , (8.28)
¬ c2 ¼ R
Q
where r1 is the set point for the moisture, and r2 the set point for uc2.
Notation ǔ(k + i | k) denotes the i step-ahead prediction and ǻ is the
difference operator. Since the MPC formulation is inherently discrete, the
cost function is given as a summation. The idea of (8.28) is to let uc2 take
care of the variations in paper moisture and let uc1 be positioned at a level
where uc2 has an adequate control range in both directions in steady-state.
In this way, the first part of the drying section serves as the base level of
the drying while the last part controls the moisture. This is done by
choosing appropriate weighting matrices, Q and R. For the simulations in
this chapter they are chosen as
ª500 0 º ª100 0 º
Q « 0 10 3 », R « 0 10 4 » . (8.29)
¬ ¼ ¬ ¼
202
Chapter 8. Object-Oriented Modeling and Control of the Moisture
7.6
y (%) 7.4
7.2
7.0
450
u c1 (kPa)
440
430
600
u c2 (kPa)
500
400
0 200 400 600 800 1000 1200
Time (s)
Figure 8.17 Simulation of mid-range MPC. There is a set point change at t = 100 s and a
disturbance in inlet moisture content to the drying section from 60.0 % to 60.5 % at
7.6
7.4
y (%)
7.2
7.0
450
uc (kPa)
440
430
0 200 400 600 800 1000 1200
Time (s)
Figure 8.18 Simulation of single loop MPC. A set point change occurs at t = 100 s and a
disturbance in inlet moisture content from 60.0 % to 60.5 % at t = 600 s.
203
Chapter 8. Object-Oriented Modeling and Control of the Moisture
1.5
Magnitude (abs)
0.5
0 -4 -3 -2 -1
10 10 10 10
Frequency (Hz)
Figure 8.19 The sensitivity function for single loop MPC (dotted) and mid-range MPC
(solid), using the linearized model (8.26). The maximum sensitivity is chosen to 1.25. The
difference in limit frequency for the two systems is around a factor two.
full response of a change in ǻuc1 and ǻuc2) and a small control horizon to
limit the computational effort for the MPC. There is also a rate constraint
for the control signals
204
Chapter 8. Object-Oriented Modeling and Control of the Moisture
7.6
Moistuure (%)
7.4
7.2
7.0
6.8
20
Steam consumption (kg/s)
19
18
17
16
15
0 200 400 600 800 1000 1200
Time (s)
Figure 8.20 Comparison between mid-range MPC (solid) and single loop MPC (dotted),
showing the moisture content and steam consumption.
50 3
¦ >r (k i) yˆ (k i | k )@ Q ¦
2 2
J (k ) 'uc (k i ) R , (8.31)
i 1 i 0
and the chosen weights that give the sensitivity in Figure 8.19, and
constraint are
Figure 8.19 also indicates the difference in performance between the two
controllers, at least in the non-constrained case. The mid-range MPC has a
limit frequency almost twice as large as the limit frequency of the single
loop MPC. Since paper machines often are run from several hours to days
with the same set point, disturbance rejection is important and it therefore
makes sense to look at this property.
Figure 8.20 shows the moisture response for the two different
controllers together with the steam flow in the header. The difference in
performance is apparent. The mid-range MPC has both better set point
following and disturbance rejection. However, the transient steam
205
Chapter 8. Object-Oriented Modeling and Control of the Moisture
7.6
Moisture (%)
7.4
7.2
7.0
6.8
25
Steam consumption (kg/s)
20
15
10
0 200 400 600 800 1000 1200
Time (s)
Figure 8.21 Comparison between mid-range MPC (solid) and single loop MPC (dotted),
showing moisture content and steam consumption when tuned to give similar performance
in disturbance rejection.
consumption is twice as large for the mid-range MPC. This is the price
being paid for the extra performance.
In Figure 8.21, the single loop MPC has been tuned to give similar
performance in disturbance rejection as the mid-range MPC, shown in
Figure 8.17 and Figure 8.20. The output of the single loop MPC is in this
case much more aggressive than it was in Figure 8.18 and therefore the
transient steam consumption is much larger. This to be able to match the
performance of the mid-range MPC. The maximum value of the
sensitivity is also larger, hence the control system is less robust, see
Figure 8.22.
The advantage with the proposed control structure is that, in general, it
does not require any rebuild of the drying section. The cylinders are
normally divided into different groups even though all groups follow the
same set point. Therefore, it is simply a matter of changing the controller
software and a majority of the main system vendors have the possibility to
include a MPC package into their DCS system. However, it is an
advantage if the flash steam from the last group is recirculated through a
thermo compressor, see Figure 2.5. If the flash steam is reused by another
group it is important to let the steam pressure of the last group be
constrained so that its pressure never falls below the steam pressure of the
206
Chapter 8. Object-Oriented Modeling and Control of the Moisture
1.5
Magnitude (abs)
1
0.5
0 -4 -3 -2 -1
10 10 10 10
Frequency (Hz)
Figure 8.22 The sensitivity function for single loop MPC (dotted) and mid-range MPC
(solid), when tuned to give similar performance
receiving group. This guarantees that the steam flows in the intended
direction.
This mid-ranging MPC has two essential advantages compared to the
mid-ranging structure used in Chapter 6, see Figure 6.10. Firstly, the MPC
does not require any decoupling filter, CD, since the complete multi-
variable process model is given in the cost function. The MPC therefore
automatically balances the two process inputs to get a desired output and
there is no need for the extra ‘fix’ with a decoupling filter between the two
process inputs as in Figure 6.10. Secondly, the MPC has a systematic
implicit way to treat constraints. For the structure in Chapter 6, this has to
be taken care of explicitly and not as a part of the tuning phase. A
disadvantage can be that it is a bit more difficult to tune a mid-range MPC
compared to e.g. the IMC or Dahlin controller used by many mills today.
8.5 Summary
207
Chapter 8. Object-Oriented Modeling and Control of the Moisture
The model is linked to a MPC toolbox in Matlab to evaluate a new
strategy to control the moisture in the drying section. The strategy utilizes
the possibility to divide the drying section into two parts. By controlling
this multi-variable process with a mid-range structure, the performance of
the closed loop system is greatly improved. The mid-ranging MPC is
compared with single-loop MPC and the evaluation is done by comparing
performance in terms of sensitivity function, disturbance rejection and
steam consumption. An important advantage of the mid-ranging MPC is
that it does not require any rebuild of the physical process.
The advantage of this model compared to the one used in Chapter 6
and Chapter 7 is the simulation time, which is a few minutes instead of a
few days. The model used in the previous two chapters is capable of
simulating microscopic phenomena like water and vapor transport, and
moisture gradients, inside the paper sheet. The model described in this
chapter neglects this and assumes a homogeneous paper in the thickness
direction. For the purpose of the feedback control investigated in this
thesis, this assumption sufficient. However, the feedforward control in
Chapter 7 requires a model that takes differences in the thickness direction
into account to observe the peak position.
208
Chapter 9. Conclusions
9.
9
Conclusions
9.1 Summary
In this thesis, modeling and control of the paper machine drying section is
presented. The treatment is divided into two main parts, control of the
steam and condensate system, and moisture control.
209
Chapter 9. Conclusions
option to balance between robustness and performance the tuning rule has
a design parameter. A nice property of the tuning rule is that the controller
settings are easily calculated from the process parameters of the IPZ
model. It is successfully verified from experiments on real paper machine
at different mills. The IPZ structure can be used in many other process
areas, apart from modeling the steam pressure in cylinders. Therefore the
tuning method has a much wider field of applications than solely drying
section control.
The IPZ model is used to examine some properties of the closed loop
system and to investigate two different control structures. The main
objective is to avoid resonance peaks in the closed-loop frequency
response from set point to process output, and this is accomplished by a
set-point feedforward controller (two-degree-of-freedom controller) and a
state feedback controller. The two types of structures solve the objective
differently but they roughly give the same performance. However, the
two-degree-of-freedom controller is recommended since it is best suited
for implementation in today’s industrial control systems. As this controller
structure separates the servo and regulation problem it also fits well to the
IPZ tuning. The set-point feedforward design gives the desired response
from set-point to output, while a PID controller based on IPZ tuning takes
care of the disturbance rejection in the feedback loop.
By making some assumptions and a linearization, it is shown that the
black-box model can be derived from a few simple mass and energy
balances. This means that the process parameters of the black-box model
are directly related to physical quantities as cylinder dimension, valve
characteristics, and steam properties. The model is validated against
measurements from a real paper machine and the heat transfer coefficient
from steam to cylinder is used as calibration parameter.
Moisture control
Traditionally the moisture in the paper sheet is controlled by manipulating
the pressure in the steam cylinders. To improve the performance of the
closed loop moisture control system both feedforward and feedback is
investigated. A new measurement signal, called the peak position is
introduced. It is based on measurements of the paper surface temperature
and indicates the position where the surface becomes dry. It is shown that
peak position can significantly improve the control performance by
feedforward.
The feedback is carried out by two types of mid-ranging controllers.
The first uses the air supply system in combination with the steam system.
210
Chapter 9. Conclusions
The advantage of the fast impact from the supply air and the high capacity
of the steam system are combined to enhance the performance of the
closed loop system. The second mid-ranging structure utilizes similar
advantages by combining the last drying group with the remaining part of
the dryers to improve performance. Two different of mid-ranging
controllers are used. The first mid-ranging structure combines two single-
loop controllers, in this case IMC, together with a decoupling filter and a
method to tune the system is given. The second structure is implemented
in a MPC. The advantages of the mid-ranging MPC are that the
decoupling filter is not needed and that constraints are explicitly taken
into account.
From mass and energy balances, a dynamic simulation package of a
drying section is implemented in the object-oriented modeling language
Modelica. By drag-and-drop features, a simulation model of practically
any drying section can easily be built. It can be used to investigate effects
of a rebuild, different controller structures, and different ways to operate
the machine. All this can be done for web breaks, grade changes, start-
ups, and normal run.
211
Chapter 9. Conclusions
The differential pressure loop has only been treated briefly. A deeper
analysis and a tuning rule, like the one for the pressure loop, are desired.
This loop will probably benefit from a feedforward structure since it is
greatly affected by the operation of the pressure loop.
From the linear model of the cylinder dynamics there is a potential of
making a recursive identification of the heat transfer coefficient for fault
detection with respect to condensate evacuation. This would be very
beneficial since the heat transfer coefficient have a large impact on the
dynamics of the moisture process and consequently its control
performance.
A systematic analysis of moisture control where different control
paradigms as MPC, Smith-predictor, IMC, PID, deadbeat control,
minimum-variance, or LQG, are compared would be valuable. Every
control system vendor use their control structure but no thorough
comparison have been published in the literature.
212
Appendix A. Glossary
Appendix A
Glossary
213
Appendix A. Glossary
Caliper: The thickness of a sheet of paper or board. Also known as
thickness or bulk. Usually measured in nanometers.
Chemical pulp: Pulp in which wood fibres have been separated by
chemical, rather than mechanical, means.
Chemi-thermomechanical pulp (CTMP): Chemi-mechanical pulp
produced by treating wood chips with chemicals (usually sodium
sulphite) and steam before mechanical defibration.
Coated paper: The uniform application of a coating yields a more
even and more closed surface of printing papers, which is suitable for the
reproduction of fine screen artwork. The coating is applied in separate
coaters or in the paper machine.
Conditioned basis weight: dsd
Corrugated board: Formed fluting that is faced with a liner on both
sides.
Curl: The degree of curvature, measured when paper is positioned
flat.
Dry end: Final part of the paper machine from the drying section
onwards.
Dry weight: Mass of paper sheet excluding the moisture.
Fabric tension: The tension of the fabric creates contact pressure
between the wet web and cylinder. Higher contact pressure gives higher
heat transfer.
Filler: A material or substance that is added to the furnish to smooth
out the spaces between fibers, enhancing the printing properties of the
paper.
Fine paper: A broad term including printing, writing, and cover
papers, as distinguished from wrapping papers and paper not generally
used for printing purposes, which are generally referred to as coarse
papers.
Fines: Small fragments of fibres produced, for example, in the course
of beating or refining.
Fluting: Paper that has been formed into the flutes that make up the
ridged part of the corrugated board between the liners.
Formation: The arrangement of fibers in a sheet of paper; can be seen
by holding it up to a light source.
Freeness: A measure of the drainability of an aqueous suspension of
stock, determined and expressed as specified in a standard method of test.
Furnish: The ingredients or constituents of which the paper is made.
Fibrous furnish concerns only the types and proportions of the fibres
present.
214
Appendix A. Glossary
Grade: (1) A class or level of quality of a paper or pulp which is
ranked, or distinguished from other papers or pulps, on the basis of its
use, appearance, quality, manufacturing history, raw materials, or a
combination of these factors. Some grades have been officially identified
and described; others are commonly recognized but lack official
definition. (2) With reference to one particular quality, one item (q.v.)
differing from another only in size, weight, or grain; e.g., an offset book
paper cut grain long is not the same grade as the same paper cut grain
short.
Fiber orientation: The direction of the fibers in paper.
Kraft paper: A paper made essentially from wood pulp produced by
a modified sulfate pulping process. It is a comparatively coarse paper
particularly noted for its strength, and in unbleached grades is primarily
used as a wrapper or packaging material. It can be watermarked, striped,
or calendered, and it has an acceptable surface for printing. Its natural
unbleached color is brown but by the use of semibleached or fully
bleached sulfate pulps it can be produced in lighter shades of brown,
cream tints, and white. In addition to its use as a wrapping paper, it is
converted into such products as: grocery bags, envelopes, gummed
sealing tape, asphalted papers, multiwall sacks, tire wraps, butcher wraps,
waxed paper, coated paper, as well as specialty bags and sacks.
Liner: A paper that is used as the facing material in the production of
corrugated and solid fibre shipping containers.
Mechanical pulp: Pulp consisting of fibres separated entirely by
mechanical rather than chemical means.
Newsprint: A lightweight paper, made mainly from mechanical wood
pulp, engineered to be bright and opaque for the good print contrast
needed by newspapers. Newsprint also contains special tensile strength
for repeated folding. It does not includes printing papers of types
generally used for purposes other than newspapers such as groundwood
printing papers for catalogs, directories, etc.
Opacity: The ability of a sheet of paper to prevent light transmission
through it. Opacity prevents print that is on one side of a sheet of paper
from showing through to the other side.
Packaging papers: These papers are used to wrap or package
consumer and industrial products such as grocer's bags and sacks,
shopping and merchandise bags, and multiwall shipping sacks used for
shipping such products as cement, flour, sugar, chemicals and animal
food. "Specialty" packaging papers are used for cookies, potato chips, ice
cream, and similar products.
215
Appendix A. Glossary
Paper Stock: A mix of pulp fibers, water, additives, chemicals and
dyes that will be pumped onto the paper machine to form paper.
Picking: Fibers in the paper which tend to pull away from the surface
during the drying process.
Recovery boiler: In wood pulping, a unit for concentrating black
liquor to a stage where the residual carbon is then burned out and the
inorganic sodium salts melted and recovered.
Refiner: A machine, usually equipped with discs or with a cone and
plug, intended for the treatment of fibrous materials in an aqueous
medium to give them some of the properties needed for the manufacture
of pulp or paper with the necessary characteristics.
Retention: Proportion of fibre and filler retained on the paper
machine wire.
Sheet (or Web): A continuous length of paper, such as paper when it
is on the paper machine or on roll-feed equipment.
Sizing: This process can either be applied on the surface of the sheet
or in the sheet: in the first case starch is applied to the surface to increase
its strength and to resist the penetration of oil-based inks (this process is
carried out at the size press, which is about two-thirds of the way down
the dry end); in the second case chemicals are added to the stock at the
pulping stage before the sheet is formed: this is called internal or engine
sizing and its purpose is to stop penetration of water-based inks into the
sheet.
Smoothness: It is a measured property of paper that describes or rates
the flatness and evenness of a sheet's surface.
Stock: The wet pulp before it is fed onto a paper making machine, or
during the paper making processes before it becomes a sheet of paper;
contains around 99% water and 1% fibre.
Stock preparation: Collective term for all treatment necessary for the
preparation of the stock before it reaches the making machine.
Tissue: Collective term for papers of a grammage of less than 30 g/m2
that differ in application and composition but have the common feature of
being thin. Examples of different types of tissue papers include sanitary
grades such as toilet, facial, napkin, towels, wipes, and special sanitary
papers. The extremely thin Japanese tissue papers are sometimes
produced in grammages as small as 6 to 8 g/m2. Other examples of tissue
papers are decorative and laminated tissue papers and crepe papers, often
used in gift wrapping and to decorate. Desirable characteristics here are
appearance, strength, and durability.
Wet end: First part of the paper machine up to the drying section.
White water: The filtrate from the wet end of the paper machine
216
Appendix A. Glossary
Winding: An operation in which the paper or board from the paper
machine is slit and wound into the roll widths ordered by the customer.
Woodfree paper: Paper consisting of chemical pulp fibres. It does
not contain any mechanical pulp beyond a permissible content of 5 % by
mass.
217
Appendix B. Conservation Balance for Compartmental Models
Appendix B
Conservation Balance for
Energy in Compartmental
Models
The models developed in Chapter 4 and Chapter 8 are based on mass and
energy balances, together with a few constitutive relations. The balance
equations come from the fundamental conservation principle in physics,
that mass and energy are neither destroyed nor created but simply change
their form. Putting up a mass balance is rather straightforward but the
energy balance is not always so evident, particularly not for gas systems.
In this appendix the energy balance for the steam and condensate in the
cylinder model is discussed. It will be shown why change in internal
energy depends on in- and outflow of enthalpy, see Figure A.1.
control volume
d
qin hin U i qout hout
dt
Figure A.1 The general energy balance for steam and condensate, used in this thesis. q is
mass flow, h specific enthalpy, and Ui internal energy. The change in internal energy
depends on inflow and outflow of enthalpy.
218
Appendix B. Conservation Balance for Compartmental Models
B.1 The energy balance
x internal energy Ui
x kinetic energy Ek
x potential energy Ep
Observe that we here use internal energy Ui [J] and not specific internal
energy ui [J/kg]. Energy flows into the system in terms of
heat output
work output
energy leaving the control volume with the outgoing fluid
work done by the outgoing fluid
219
Appendix B. Conservation Balance for Compartmental Models
where qin and qout are mass flows [kg/s], and ek and ep are specific energies
[J/kg].
Mechanical work
This is denoted as Wm [J/s] and can be both shaft work and expansion
work. Shaft work is the work done by the fluid or surrounding through a
moving part (e.g. a pump rotor). Expansion work accounts for work done
in expansion or contraction of the control volume. This term is often not
present especially not when the control volume is fixed. Work done on the
system is positive work.
Flow work
This term relates to the work done on the fluid as it moves into and out of
the control volume. This is expressed as
where pin and pout are pressures [Pa] at inlet and outlet, and vin and vout are
specific volumes [m3/kg] at inlet and outlet. To see this, consider a unit
volume of fluid across the boundary at the input, see Figure A.2. Let the
element have a length l [m/kg] and the cross-sectional area of the inlet be
A [m2] Then the energy required to push the element across the boundary
is given by the force times length,
Similarly it can be shown that the flow work at the outlet is given by
pout vout
d
(U i Ek E p ) qin (ui ,in ek ,in e p ,in ) qout (ui ,out ek ,out e p ,out )
dt
qin pin vin qout pout vout Qc ,in Qc ,out Wm ,in Wm ,out
220
Appendix B. Conservation Balance for Compartmental Models
The sum of specific internal energy and the pv term is called specific
enthalpy, and is given by the symbol h (h = u + pv). We can therefore
write the energy balance as
d
(U i Ek E p ) qin (hin ek , in e p ,in ) qout (hout ek , out e p , out )
dt
Qc , in Qc , out Wm , in Wm , out
When the steam cylinder model is developed in Chapter 4, the kinetic and
potential energies are neglected, which is applicable to most systems in
the process industry. Also, there are evidently no mechanical work terms
acting on the system and the internal energy of the control volume is
written as specific internal energy times mass. It is then simplified to
d
( U uiV ) qin hin qout hout Qc , in Qc , out
dt
control volume
inlet
Figure A.2 A small section at inlet to the system. There is an expenditure in energy to
push the fluid across the boundary.
221
Appendix B. Conservation Balance for Compartmental Models
If we are dealing with an incompressible liquid, there is no pressure-
volume work done on the system, and the energy balance can be further
simplified to
d
( U hV ) qin hin qout hout Qc , in Qc , out
dt
This is the case when considering the water phase (moisture) in the paper
web.
222
Appendix C. Solution to the One Dimensional Heat Equation
Appendix C
Solution to the One
Dimensional Heat Equation
wT
Uc ( O T ) k , (C.1)
wt
wT w 2T
a 2 0, (C.2)
wt wx
O
a . (C.3)
Uc
223
Appendix C. Solution to the One Dimensional Heat Equation
T2
T1
Steam side Paper side
T0
T
L
x
Figure C.1 At t = 0, the temperature on the steam side is changed from T1 to T2. The two
steady-state solutions are shown in the figure.
wT w 2T
° a 0, 0 x L, t ! 0
° wt wx 2
°T ( x,0) T x
° 1 (T0 T1 )
L (C.4)
®
°T (0, t ) T , t !0
2
°
° wT ( L, t )
° O Q
¯ wx
224
Appendix C. Solution to the One Dimensional Heat Equation
2
f § 1· S2
Q 4(T1 T2 ) a ¨© k 2 ¸¹ L2
t
§ 1·S
T ( x , t ) T2
O
x ¦
k 1 ( 2k 1)S
e sin ¨ k ¸ x. (C.5)
© 2¹ L
Q DA(T1 T2 ), (C.6)
225
Appendix C. Solution to the One Dimensional Heat Equation
Figure C.2 The temperature as a function of both position x and time t (given in seconds).
T1 = 120°C, T2 = 130°C, L = 0.03 m, c = 500 J/(kg°C), Ȝ = 50 W/(m2°C), U = 7000 kg/m3,
Q = 10 kW.
130
128
126
t = 500
Temperature ( C)
124 t = 100
o
t = 50
122 t=5
t = 20
t=1
120
t=0 t = 10
118
116
114
0 0.5 1.0 1.5 2.0 2.5 3.0
Shell thickness (cm)
Figure C.3 The temperature as a function of position x for different points of time.
226
Appendix C. Solution to the One Dimensional Heat Equation
125
124
123
122
121
Temperature ( C)
o
120
119
118
117
116
115
114
0 20 40 60 80 100 120 140 160 180 200
Time (s)
Figure C.4 The temperature at the cylinder surface (paper side) as a function of time. The
response is similar to a step response of a higher order system.
128
127
126
125
Mean temperature ( C)
124
o
123
122
121
120
119
118
117
0 20 40 60 80 100 120 140 160 180 200
Time (s)
Figure C.5 The mean temperature of the cylinder shell as a function of time. The response
is similar to a step response of a first-order system.
227
Appendix C. Solution to the One Dimensional Heat Equation
Heat conduction in an infinitely long rod
Regard the heat equation for an infinitely long rod
wT w 2T
° a 2 0, x , t ! 0
® wt wx (C.7)
°
¯T ( x,0) g ( x)
This equation can be solved by the Fourier or Laplace transform and the
solution is given by the convolution formula
f
T ( x, t ) G
g ( x, t ) ³ G( x D ) g (D )dD ,
f
(C.8)
where
1 2
G ( x, t ) e x 4at
. (C.9)
4S at
G is called the Green function for the heat equation. If we assume that the
initial condition is given by the Dirac pulse, the solution to (C.7) is
f
1 2
T ( x, t ) G
G ³ G( x D )G (D )dD e x 4at
. (C.10)
f
4S at
Figure C.6 shows this solution as a function of both time and space. The
integral of (C.10) over the whole space is proportional to the thermal
energy in the rod, and evaluates to
f
1 2
³ e x 4at
dx 1, (C.11)
f
4S at
independently of t and a. This because no energy leaves the rod and the
Dirac pulse has unit area.
The implication of (C.10) is that the energy impulse in the origin
affects the temperature in each point x for all t > 0. In other words, all
points in space will sense the impulse instantaneously, and the model
228
Appendix C. Solution to the One Dimensional Heat Equation
(C.1) leads to an infinite propagation velocity for the heat at the starting
point. This is due to the averaging effect of the integral in (C.8). Loosely
speaking, the heat equation does not admit any discontinuous solutions.
This non-physical property of the heat equation has been extensively
discussed in the literature, first explicitly solved by [Cattaneo, 1948]. The
theory has become known as extended thermodynamics, and is a field of
research still today [Herwig and Beckert, 2000], [Pulko, et al, 2002], and
[Belevich, 2004]. [Joseph and Preziosi, 1989] gives a good overview of
the field. Often the drawback of infinite velocity is overcome by adding
an extra term to (C.1) involving a derivative of second order in time. It
then becomes a hyperbolic function which can be regarded as a “wave
equation for heat” (with a certain speed of the heat wave).
Figure C.6 Temperature as a function of both position x and time t (given in seconds), for
the normalized case a = 1.
229
References
References
Allison, B. J., and S. Ogawa (2003): “Design and tuning of valve position
controllers with industrial applications.” Transactions of the Institute
of Measurements and Control, 25(1), pp. 3í16.
230
References
Åström, K. J. (1967): “Computer control of a paper machine í an
application of linear stochastic control theory.” IBM Journal of
Research and Development, 11(4), pp. 389í405.
231
References
Bergström, J., and G. Dumont (1998): ”An object oriented framework for
developing dynamic models of a paper machine.” Dynamic Modeling
Control Applications for Industry Workshop, IEEE Industry
Applications, pp. 63í69.
232
References
Carlberg, M. (ed) (1989): Teknisk uppkäftighet, Skogsindustrierna and
Bokförlaget Bra Böcker.
Chang, D.-M., C.-C. Yu, and I.-L. Chien (2000): “Design sensor
trajectory for control: application to sheet-forming processes.” AIChE
Journal, 46(8), pp. 1581í1592.
Chen, S.-C. (1995): “Modelling of paper machines for control: theory and
practice.” Pulp and Paper Canada, 96(1), pp. 17í21.
233
References
Cohen, G. H., and G. A. Coon (1953): “Theoretical consideration of
retarded control.” Transactions of American Society of Mechanical
Engineers, 75, pp. 827í834.
234
References
FFIF (2004): “Paper and wood, yearbook 2004, statistics 2003.” Finnish
Forest Industries Federation, Download can be found at
www.forestindustries.fi.
235
References
Heaven, E. M., I. M. Jonsson, T. M. Kean, M. A. Manness, and R. N.
Vyse (1994): “Recent advances in cross machine profile control.”
IEEE Control System Magazine, 14(5), pp. 35í46.
236
References
Karlsson, I. (2002): ”Förfarande vid torkning av ett banformigt material.”
Swedish patent, 0004851-2 (SE518049).
237
References
Larsson, J. E., and T. Gustafsson (1998): “Paper machine dry line
position control during grade changes.” Proceedings of the American
Control Conference, pp. 2997í3001, Philadelphia, PA.
Ljung, L. (1999): System identification - Theory for the user, 2nd edition,
Prentice Hall.
238
References
Mattsson, S. E., H. Elmqvist, and M. Otter (1998): “Physical system
modeling with Modelica.” Control Eng. Practice, vol. 6, pp.
501í510.
239
References
Nelson, D., and T. Gardner (1996): “Optimizing paper machine control -
a case study.” Pulp and Paper Canada, 97(11), pp. 43í49.
Nissan, A., and W. Kaye (1955): “An analytical approach to the problem
of drying thin fibrous sheets on multicylinder machines.” Tappi
Journal, vol. 38, no. 7, pp. 385í398.
240
References
Poulin É., and A. Pomerleau (1999): “PI settings for integrating processes
based on ultimate cycle information.” IEEE Transactions on Control
Systems Technology, 7(4), pp. 509í511.
Rao, M., Q. Xia, and Y. Ying (1994): Modeling and advanced control for
process industries: application to paper making processes, Springer-
Verlag, New York.
241
References
Sell, N. J. (ed) (1995): Process control fundamentals - for the pulp and
paper industry, Tappi Press.
Skogestad, S. (2003): “Simple analytic rules for model reduction and PID
controller tuning.” Journal of Process control, 13, pp.291í309.
242
References
Stenström, S., B. Wilhelmsson, L. Nilsson, R. Krook, and R.
Wimmerstedt (1994): “Measurement of reference experimental
drying data for the multi-cylinder paper dryer.” Proceedings of the 9th
International Drying Symposium (IDS’94), Gold Coast, Australia, pp.
1179í1186.
Sun, X., K. Yi, and Y. Sun (2000): “The modeling and control of basis
weight and moisture content.” Proceedings 3rd World Congress on
Intelligent Control and Automation, Heifei, P. R. China, pp.
3724í3728.
243
References
Wallén, A. (2000): Tools for autonomous process control, PhD-thesis,
Department of Automatic Control, Lund Institute of Technology,
Sweden.
Wang, Y.-G., W.-J. Cai (2001): “PID tuning for integrating processes
with sensitivity specification.” Proceedings Conference on Decision
and Control, pp. 4087í4091, Orlando, FL.
Xia, Q, M. Rao, and J. Qian (1993): “Model algorithmic control for paper
machines.” Proceedings of the IEEE Conference on Control
Applications, pp. 203í208, Vancouver, Bristish Columbia.
244
List of Symbols
List of Symbols
245
List of Symbols
n í Noise
p Pa Steam pressure inside cylinder
psh Pa Steam pressure in header
ptot Pa Standard pressure (101.325 kPa)
pv0 Pa Partial vapor pressure for free water
pv,a Pa Partial pressure for water vapor in the air
pv,p Pa Partial pressure for water vapor at paper surface
P(s) í Process transfer function
qbt kg/s Blow through steam
qc kg/s Condensation rate
qevap kg/(m2·s) Evaporation rate
qs kg/s Inflow of steam to cylinder
qw kg/s Outflow of condensate
Q í Weighting matrix for MPC
Qm W Energy flow to cylinder shell
Qp W Energy flow to paper
r í Set point
R í Weighting matrix for MPC
Rg J/(mole·K) Gas constant ( § 8.31)
Rv í Valve rangeability
Td s Derivative time of the PID controller
Ti s Integral time of the PID controller
Tm K Temperature of cylinder shell
Tp K Temperature of paper
u kg/kg Moisture ratio
uc í Controller output
us J/kg Internal energy of steam
uw J/kg Internal energy of condensate
246
List of Symbols
V m3 Volume inside cylinder
Vs m3 Volume of steam
3
Vw m Volume of condensate
vx m/s Speed of paper sheet
w % Moisture content
x kg/kg Water content in air
xv í Valve opening
y í Process output
247
Department of Automatic Control
ISSN 0280-5316
ISRN LUTFD2/TFRT--1075--SE