Dynamics and Control of Chemical Processes Prof. Davide Manca Report of Lab 3 "Design of The Control System"
Dynamics and Control of Chemical Processes Prof. Davide Manca Report of Lab 3 "Design of The Control System"
System Representation
Solution
In our system we have one controlled variable which is level of the second tank and one
manipulated variable that is the outlet flow rate of second tank. The aim of the control system
is to keep constant the level in the second tank in the case of perturbation. We will assume
that the system is in steady state when the step disturbance occurs.
First step in this problem is to evaluate system’s characteristic parameters by defining open
loop dynamics of the controlled variable (liquid level in the second tank). Then, we can find the
time delay, gain of the system and characteristic time to reach the new steady state.
Dynamic of the open loop: (Without Controller)
𝑑ℎ1 ℎ1 −ℎ2 ℎ1 (𝑠)−ℎ2 (𝑠)
1) 𝐴1 = 𝐹𝑖 − 1) 0 = 𝐹𝑖 (𝑠) −
𝑑𝑡 𝑟1 𝑟1
𝑑ℎ2 ℎ1 −ℎ2 @Steady state ℎ1 (𝑠)−ℎ2 (𝑠)
2) 𝐴2 = − 𝐹𝑜 2) 0 = − 𝐹𝑜 (𝑠)
𝑑𝑡 𝑟1 𝑟1
We can easily evaluate the values of height in each tank at steady state condition.
• ℎ2 (𝑠) = 6.6 𝑚
• ℎ1 (𝑠) = 𝐹𝑖 ∗ 𝑟1 + ℎ2 (𝑠) = 17.88 𝑚
Flowrate after step disturbance will be:
(𝑠,𝑛𝑒𝑤) (𝑠,𝑜𝑙𝑑)
ℎ2 −ℎ1
Gain of the system calculated by: 𝐾= (𝑛𝑒𝑤) (𝑜𝑙𝑑) = 0.696 𝑠/𝑚2 where time delay
𝐹𝑖 −𝐹𝑖
and characteristic time to reach new steady state are equal to 𝑡𝑑 = 8.8 𝑠 and 𝜏 = 115.6 𝑠.
Now, by using Cohen-Coon controller tuning method it is possible to estimate the system
parameters of proportional control.
1 𝜏 𝑡𝑑
𝐾𝑐 = (1 + ) = 19.35
𝐾 𝑡𝑑 3𝜏
The general form of the model of the proportional control is: 𝑐 = 𝑐𝑠 + 𝐾𝑐 (ℎ − ℎ𝑠𝑝 )
Where is the bias 𝑐𝑠 , i.e. the value of the manipulated variable when the controlled variable is
at the set-point. In this case a valve on the inlet flowrate is manipulated and the net result is
the modification of the second tank outlet flowrate.
Assuming a linear relationship between the valve opening degree and the outlet flowrate, the
controller model can be re-written as:
𝐹 = 𝐹0𝑠 + 𝐾𝑐 (ℎ − ℎ𝑠𝑝 )
Exercise 2.
Again, it is requested to model and compute the dynamics of the same system, but with PI
controller.
Solution
By using Cohen-Coon controller tuning method it is possible to estimate the system parameters
of proportional-integral control.
1 𝜏 𝑡𝑑
𝐾𝑐 = ∗ ∗ ( 0.9 + ) = 17.11
𝐾 𝑡𝑑 12𝜏
3𝑡
𝑡𝑑 ∗ (30 + 𝜏𝑑 )
𝜏𝐼 = = 25.28 𝑠
20𝑡𝑑
(9 + 𝜏 )
Assuming a linear relationship between the valve opening degree and the outlet flowrate, the
controller model can be re-written as:
𝐾𝑐 𝑡
𝐹 = 𝐹0𝑠 + 𝐾𝑐 (ℎ − ℎ𝑠𝑝 ) + ∫ (ℎ(𝑡) − ℎ𝑠𝑝 )𝑑𝑡
𝜏𝐼 0
So, the model of the system becomes:
𝑑ℎ1 ℎ1 − ℎ2
1) 𝐴1
= 𝐹𝑖 −
𝑑𝑡 𝑟1
𝑑ℎ2 ℎ1 − ℎ2 𝐾𝑐 𝑡
2) 𝐴2 = − [𝐹0𝑠 + 𝐾𝑐 (ℎ − ℎ𝑠𝑝 ) + ∫ (ℎ(𝑡) − ℎ𝑠𝑝 )𝑑𝑡]
𝑑𝑡 𝑟1 𝜏𝐼 0
As a result, the following graph obtained for PI controller.
Exercise 3.
It is asked to solve the previous problem with a setpoint change in the height of the second to
8.6 m.
The response of each P and PI control action due to the set point change illustrated above.
Again the plots have the same structure in comparison with plot of old set point, where PI
control action affect level to reach the steady state more quickly.