INTRODUCTION
The contribution of this article is study of two different control strategies for a class of chemical process industries. The model of the distillation column used in this plant, which is characterized by an underactuated structure. According to the diversity and complexity of these systems, it is important to emphasize that none of the technique proposed and developed for fully actuated systems can be applied directly to any underactuated system. Therefore, it is meaningful to develop control methods for this class of systems, and more precisely to develop the most optimized controller design. The control objective for systems characterized by the fact that there are more degrees of freedom than actuators, is to obtain a desirable behavior of several output variables by simultaneously manipulating several inputs channels. Under-actuated systems are less sensitive to modelling errors, so it has to be controlled in its original form to obtain robust stability and performance[1]. Garcia and Morari (1985), and Morari and Zafirious (1989) presented the concept of Internal Model Control (IMC) and have already proved the effectiveness of the framework for robust control of different kinds of Singular-Input Singular Output (SISO) systems. The IMC structure is composed of three principal parts: the process, the internal model in parallel to the process and the controller, which is the dynamic inversion of the model[2]. Recently, many methods have been proposed to control MIMO systems with time delays. For example, Xianku and Hongshuai (2015) proposed a Closed-loop Gain Shaping Algorithm (CGSA) using Padé approximation, which is sufficiently accurate in view of stability analysis[3]. Qibing, Beiyan, Qi and Guanfei (2014) introduced a design method of decoupling internal model control; the basic idea is to realize the decoupling of the controller of non-square processes by inserting some compensation Relative Normalized Gain Array (RNGA). An equivalent transfer function matrix is introduced to approximate the pseudo-inverse of the process transfer function matrix, which makes the design of decoupling internal model control simple and easy to calculate[4]. Guanfei and Qi (2013) integrated IMC with disturbance controllers by choosing different forms of external input/output disturbance. When this disturbance is applied directly as input to the controller, the design of controller needs to compensate the effect of slow dynamic poles by adding some constraints[5]. Pamela, Jebarajan and Baby introduced an approach to regulate the heater power in those systems, which must control the temperature in food processing, pharmaceuticals and in polymerization. The objective is to control the system with both PI controller and IMC structure and to analyze its performance parameters[6]. Qi, Laing, Feng and Si-wen (2012) proposed a novel design IMC controller based on Singular Value Decomposition (SVD), this approach uses SVD in the inverse of the steady-state gain matrix of process[7]. The last decades have shown an increasing interest in the control of under-actuated systems, many IMC methods can be introduced to achieve considerable results on the control for these kinds of systems. For this class of system, the number of inputs is smaller than the number of outputs, which means that the transfer function matrix is not square, and the issue of inversion exists. As a result, the problem of inversion was solved in many literature researches previously mentioned[4-6] by decoupling the internal model control based on the Relative Normalized Gain Array (RNGA), Singular Value Decomposition (SVD) and an Equivalent Transfer Function (ETF) matrix. The solution brought by these methods require complicated calculations and many instructions to implement[7]. Otherwise, our proposed approach; Virtual Inputs is used to augment the system inputs inserting a certain number of virtual columns to the model in order to make it a square matrix. These virtual columns have no influence on the response of the system, and they will be eliminated after. The control of under-actuated systems is challenging, hence the necessity to identify the control technique with less interactions and better performances such as overshoot, setting time etc. Based on internal model control for discrete under-actuated systems a comparison based on the inversion technique of the model will be made between Singular Value Decomposition (SVD) and Virtual Inputs (VI) to realize the internal model controller. In this article, section 2 summarizes the fundamental control problem to be solved. Section 3 describe the general IMC structure. Section 4 details the design of the internal controller based on Singular Value Decomposition. Section 5 details the design of the internal controller based on Virtual Inputs. At the end, we show how our proposed methodology based on Virtual Inputs at distillation column could give good results compared to the other method described previously.
PROBLEM FORMULATION
The process to be controlled is assumed linear and discrete-time governed by the following equation[8]
Wheren x Î ℝn i is the system state, u Î ℝm is the vector of manipulated variables, y Î ℝ n is the vector of system outputs, is the time delays, F, H and C are matrices of corresponding dimensions[8]. For the IMC design technique, it is convenient to express this model in frequency response form. Taking Z transforms of (1), the following input-output model is obtained:
Where is the system transfer
matrix of dimension
, the number
of control inputs is equal to m, the number of outputs is equal to n, making it a
rectangular matrix, which has the form
![indecs-19-561-g10.gif](/xml/21173/indecs-19-561-g10.gif)
![indecs-19-561-g11.gif](/xml/21173/indecs-19-561-g11.gif)
(i) If the model contains time delays, its inverse involves predictive terms, which make the controller unrealizable.
(ii) If the zeros of the transfer matrix of the model outside the complex unit circle yield an unstable perfect controller.
(iii) System equipped with a perfect controller is extremely sensitive to modelling errors and time delays.
(iv) The direct model inversion is also impossible in the case of underactuated systems.
In fact, the model must provide an accurate description of the process dynamics and characteristics. Therefore, the model expression must be very close to that of the plant. For underactuated systems, the number of control inputs is less than the number of outputs and therefore we will have a rectangular matrix that is invertible. This represents the major problem encountered. In this article, to resolve this problem, it is proposed to develop some methods of inversion in the case of under-actuated systems.
INTERNAL MODEL CONTROL DESCRIPTION
The development of the IMC structure has progressed in recent years in order to design an optimal feedback controller. In this section we present the general IMC structure and we describe its basic principles and properties. Due to the fact that the traditional IMC methods cannot solve the control problem of non-square systems, we introduce two design technique to realize the internal controller but, when using a matrix to describe a non-square system, the issue of inversion often emerges.
THE GENERAL IMC DESCRIPTION
The general IMC structure of multivariable systems adopted in this article is shown
in Figure 1, where ,
and
are the output vector of the process, the
output of internal model and the control variable, respectively. r(z) is the
set-point vector, d(z) is the disturbance,
and
represent the transfer function matrix of
process and its model,
is the
transfer matrix of the IMC controller.
We begin by reviewing the properties of the IMC structure. This structure is equivalent to a conventional feedback loop with controller. From Figure 1 the inputs vector u(z) and the system outputs vector y(z) are expressed by
Property 1. Dual Stability; Assuming an ideal model![indecs-19-561-g21.gif](/xml/21173/indecs-19-561-g21.gif)
![indecs-19-561-g7.gif](/xml/21173/indecs-19-561-g7.gif)
![indecs-19-561-g22.gif](/xml/21173/indecs-19-561-g22.gif)
![indecs-19-561-g17.gif](/xml/21173/indecs-19-561-g17.gif)
![indecs-19-561-g7.gif](/xml/21173/indecs-19-561-g7.gif)
![indecs-19-561-g26.gif](/xml/21173/indecs-19-561-g26.gif)
Property 2. Perfect Control; Assume that the controller, yields a closed-loop stable IMC loop, this
controller is equivalent to the inverse model to achieves perfect set-point
satisfaction despite any disturbance. Furthermore,
is often not realizable.
When
is verified, transfer
function (7) becomes
Property 3. Zero Offset; Assume that the steady state gain controller is equal to the inverse model
gain and the closed-loop system in Figure 1 is stable. The key to apply the IMC
structure is the controller who would yield the best output response possible.
However, as in the under-actuated case, this perfect controller cannot be
implemented for previous reasons. We will discuss in the next section the controller
design using two different methods; Singular Value Decomposition (SVD) and Virtual
Input (VI). These methods are based on the transfer function matrix of the model of
the process.
IMC STRUCTURE BASED ON SINGULAR VALUE DECOMPOSITION
The application of the internal model structure to under-actuated systems is considered like our main target. In this section, we will describe the design phase of the internal controller based on SVD and the implementation at the level of the structure IMC. This approach of design can realize the decoupling of under-actuated processes and eliminate the unrealizable factors by inserting compensated terms. Meanwhile, a non-diagonal filter is designed based on SVD matrix theory. We discuss with more details this approach in the following.
The IMC structure of an under-actuated systems based on SVD is shown in Figure 2, where CSVD(z) is the IMC controller.
IMC structure based on SVD method. The SVD of the transfer matrix G(z) of the under-actuated systems is as follows: Where,![indecs-19-561-g33.gif](/xml/21173/indecs-19-561-g33.gif)
![indecs-19-561-g34.gif](/xml/21173/indecs-19-561-g34.gif)
![indecs-19-561-g35.gif](/xml/21173/indecs-19-561-g35.gif)
![indecs-19-561-g36.gif](/xml/21173/indecs-19-561-g36.gif)
![indecs-19-561-g37.gif](/xml/21173/indecs-19-561-g37.gif)
![indecs-19-561-g15.gif](/xml/21173/indecs-19-561-g15.gif)
![indecs-19-561-g38.gif](/xml/21173/indecs-19-561-g38.gif)
![indecs-19-561-g39.gif](/xml/21173/indecs-19-561-g39.gif)
![indecs-19-561-g40.gif](/xml/21173/indecs-19-561-g40.gif)
![indecs-19-561-g41.gif](/xml/21173/indecs-19-561-g41.gif)
![indecs-19-561-g43.gif](/xml/21173/indecs-19-561-g43.gif)
![indecs-19-561-g43.gif](/xml/21173/indecs-19-561-g43.gif)
In the nominal case, equation
(10) becomes
![indecs-19-561-g46.gif](/xml/21173/indecs-19-561-g46.gif)
![indecs-19-561-g47.gif](/xml/21173/indecs-19-561-g47.gif)
![indecs-19-561-g48.gif](/xml/21173/indecs-19-561-g48.gif)
It is obvious that once
![indecs-19-561-g51.gif](/xml/21173/indecs-19-561-g51.gif)
![indecs-19-561-g52.gif](/xml/21173/indecs-19-561-g52.gif)
![indecs-19-561-g48.gif](/xml/21173/indecs-19-561-g48.gif)
![indecs-19-561-g53.gif](/xml/21173/indecs-19-561-g53.gif)
![indecs-19-561-g54.gif](/xml/21173/indecs-19-561-g54.gif)
![indecs-19-561-g55.gif](/xml/21173/indecs-19-561-g55.gif)
![indecs-19-561-g55.gif](/xml/21173/indecs-19-561-g55.gif)
![indecs-19-561-g57.gif](/xml/21173/indecs-19-561-g57.gif)
![indecs-19-561-g62.gif](/xml/21173/indecs-19-561-g62.gif)
![indecs-19-561-g59.gif](/xml/21173/indecs-19-561-g59.gif)
![indecs-19-561-g62.gif](/xml/21173/indecs-19-561-g62.gif)
![indecs-19-561-g76.gif](/xml/21173/indecs-19-561-g76.gif)
![indecs-19-561-g62.gif](/xml/21173/indecs-19-561-g62.gif)
![indecs-19-561-g63.gif](/xml/21173/indecs-19-561-g63.gif)
![indecs-19-561-g66.gif](/xml/21173/indecs-19-561-g66.gif)
![indecs-19-561-g70.gif](/xml/21173/indecs-19-561-g70.gif)
![indecs-19-561-g71.gif](/xml/21173/indecs-19-561-g71.gif)
![indecs-19-561-g72.gif](/xml/21173/indecs-19-561-g72.gif)
![indecs-19-561-g76.gif](/xml/21173/indecs-19-561-g76.gif)
![indecs-19-561-g74.gif](/xml/21173/indecs-19-561-g74.gif)
![indecs-19-561-g61.gif](/xml/21173/indecs-19-561-g61.gif)
![indecs-19-561-g74.gif](/xml/21173/indecs-19-561-g74.gif)
![indecs-19-561-g76.gif](/xml/21173/indecs-19-561-g76.gif)
![indecs-19-561-g77.gif](/xml/21173/indecs-19-561-g77.gif)
![indecs-19-561-g61.gif](/xml/21173/indecs-19-561-g61.gif)
![indecs-19-561-g79.gif](/xml/21173/indecs-19-561-g79.gif)
![indecs-19-561-g55.gif](/xml/21173/indecs-19-561-g55.gif)
![indecs-19-561-g80.gif](/xml/21173/indecs-19-561-g80.gif)
![indecs-19-561-g81.gif](/xml/21173/indecs-19-561-g81.gif)
![indecs-19-561-g80.gif](/xml/21173/indecs-19-561-g80.gif)
![indecs-19-561-g82.gif](/xml/21173/indecs-19-561-g82.gif)
![indecs-19-561-g55.gif](/xml/21173/indecs-19-561-g55.gif)
According to the equation (17), the robustness of the often cannot meet the requirement, for
solving this problem, this approach need to improve the filter F(z). The filtering
structure is designed in general, which makes control system bear the capacity of
high-dimensional decoupling and fast response. First, we use the SVD in the inverse
of the transfer functions matrix of process
. Then, we can use the term after
decomposition to improve the filter
and to obtain the controller based on
SVD.
Step 1: Use the SDV in the inverse of the steady-state gain matrix of
Step 2: Let satisfy the following
formula
Step 3: The improved internal model controller is as follows
The robustness of the system can be greatly enhanced by adding a
filter in the feedback loop,
and the filter time constant can be set to half of the maximum delay time in this
loop[15].
IMC STRUCTURE BASED ON VIRTUAL INPUTS (VI)
Focusing now on the inversion method of the under-actuated systems based on Virtual Inputs. To successfully apply our approach, firstly we need to modify the IMC basic structure mentioned in Figure 1, so that it becomes applicable to underactuated systems with more outputs than control inputs. Secondly, we design an approximate inverse of the model plant which is inspired by the studies of[16,17] in the case of MIMO systems and over-actuated systems[18].
The modified IMC structure presented in Figure 3 is characterized by two more blocks with respect to the basic IMC structure presented in Figure 1[19]. The first new block is the Virtual Inputs Augmentation VIA (z) which is used to augment the system inputs inserting virtual![indecs-19-561-g90.gif](/xml/21173/indecs-19-561-g90.gif)
![indecs-19-561-g91.gif](/xml/21173/indecs-19-561-g91.gif)
The inserting columns mentioned previously can be chosen as first-order transfer functions, which verify the stability criterion, and in order to simplify the study and avoid inversion problems[19]. The studied system can be shown through the following equation[18, 19]:
When the model and the system match perfectly![indecs-19-561-g93.gif](/xml/21173/indecs-19-561-g93.gif)
![indecs-19-561-g94.gif](/xml/21173/indecs-19-561-g94.gif)
![indecs-19-561-g98.gif](/xml/21173/indecs-19-561-g98.gif)
![indecs-19-561-g99.gif](/xml/21173/indecs-19-561-g99.gif)
![indecs-19-561-g100.gif](/xml/21173/indecs-19-561-g100.gif)
![indecs-19-561-g101.gif](/xml/21173/indecs-19-561-g101.gif)
![indecs-19-561-g101.gif](/xml/21173/indecs-19-561-g101.gif)
![indecs-19-561-g103.gif](/xml/21173/indecs-19-561-g103.gif)
![indecs-19-561-g104.gif](/xml/21173/indecs-19-561-g104.gif)
![indecs-19-561-g107.gif](/xml/21173/indecs-19-561-g107.gif)
![indecs-19-561-g108.gif](/xml/21173/indecs-19-561-g108.gif)
- Monovariable system where
- Overactuated system where
- Underactuated system where
SIMULATION RESULTS AND COMPARISON
To analyze the comparative study on the above control technique, we considered a chemical process industry. In order to test the control effect of discrete underactuated internal model controller based on SVD and VI, a side-stream distillation control problem suggested by Doukas and Lyben will be used[20].
Figure 5 shows the side-stream distillation scheme that serves to separate benzene, toluene and xylene. The problem posed in this case was to control the concentrations of four impurities in three product streams with only three manipulated variables: reboiler duty, reflux ratio and side stream flow rate.
Doukas and Lyben[20] treated the distillation column as a![indecs-19-561-g113.gif](/xml/21173/indecs-19-561-g113.gif)
![indecs-19-561-g113.gif](/xml/21173/indecs-19-561-g113.gif)
In this section, we evaluate the controller performance using the above-mentioned approach.
SIMULATION RESULTS USING THE SVD CONTROLLER
Starting with the internal controller based on Singular Value Decomposition, and
using the whole approach seen in section 4. Using equation (23), can be obtained by
Using equation (24), is given as follows
Using equation (25),
![indecs-19-561-g128.gif](/xml/21173/indecs-19-561-g128.gif)
![indecs-19-561-g131.gif](/xml/21173/indecs-19-561-g131.gif)
The expression of the inverse of
![indecs-19-561-g136.gif](/xml/21173/indecs-19-561-g136.gif)
![indecs-19-561-g137.gif](/xml/21173/indecs-19-561-g137.gif)
The final expression of the internal controller based on SVD described by the equation (26) is the following
SIMULATION RESULTS USING THE VIRTUAL INPUT CONTROLLER
Dealing now with our proposed approach, the Virtual Input methods applied on the same
system studied previously. Considering the same underactuated system with three
control and four outputs. The system transfer matrix G(z) is given by (33). Let us
consider the case of perfect modeling , The augmented model transfer function is
of dimension
. The system and
the model outputs are expressed respectively by (58) and (59).
The model transfer functions
![indecs-19-561-g154.gif](/xml/21173/indecs-19-561-g154.gif)
![indecs-19-561-g155.gif](/xml/21173/indecs-19-561-g155.gif)
![indecs-19-561-g156.gif](/xml/21173/indecs-19-561-g156.gif)
![indecs-19-561-g157.gif](/xml/21173/indecs-19-561-g157.gif)
![indecs-19-561-g158.gif](/xml/21173/indecs-19-561-g158.gif)
![indecs-19-561-g159.gif](/xml/21173/indecs-19-561-g159.gif)
![indecs-19-561-g160.gif](/xml/21173/indecs-19-561-g160.gif)
![indecs-19-561-g161.gif](/xml/21173/indecs-19-561-g161.gif)
![indecs-19-561-g162.gif](/xml/21173/indecs-19-561-g162.gif)
![indecs-19-561-g163.gif](/xml/21173/indecs-19-561-g163.gif)
![indecs-19-561-g165.gif](/xml/21173/indecs-19-561-g165.gif)
![indecs-19-561-g166.gif](/xml/21173/indecs-19-561-g166.gif)
![indecs-19-561-g168.gif](/xml/21173/indecs-19-561-g168.gif)
![indecs-19-561-g169.gif](/xml/21173/indecs-19-561-g169.gif)
![indecs-19-561-g170.gif](/xml/21173/indecs-19-561-g170.gif)
![indecs-19-561-g171.gif](/xml/21173/indecs-19-561-g171.gif)
The gain matrix relative to
is given by:
COMPARISON ANALYSIS
After using both Singular Value Decomposition and Virtual Inputs approaches, it seems to be quite interesting making a comparison between them showing their effectiveness in terms of stability, robustness, precision and tracking signal. From the data presented in table 1, it is relevant to note that the desired specifications of the closed-loop responses are not met with the SDV controller and this is expected because of the problem of the interaction which are too strange as illustrated by the decoupling of the controller and the simulation results obtained with the SVD controller in Figure 6. SVD method has a problem if we do not make the right choice of the parameter of the filter, in this case, the approximation cannot be made with effectiveness and the system can be diverging. Added to that this approach towards the under-actuated process is to square the system by make the decoupling of the model affecting the characteristic of the system. Unfortunately, this operation can decrease the performance of the system and making it so poor by neglecting some information’s.
The internal model controller using a Virtual Inputs as inversion approach of the
model , better performances
in terms of rise time, settling time and in terms of error dynamics precision are
obtained, as shown by table1, these results are satisfactory because the steady
state is zero for all responses as shown in Figure 7. This mean that the set-point
tracking is ensured. The internal model control using the Virtual Input approach is
having the benefits of small overshoot, faster tracking features, and less
interaction compared to the Singular Value Decomposition method as indicated by the
Figures 6 and 7 and by the table below. The Virtual Input method may be fails when
the system has a disturbance acting on the outputs of the system. For this reason,
with load disturbance the response needs to be controlled to attain robustness and
performance.
In order to verify the robustness of the proposed internal model controller, we added a white noise with a variation of 0,2 acting on the outputs of the system in the interval of time [100s, 150s, 200, 250s]. The influence of this noise is observed in Figure 8.
It is shown that the closed loop system is not sensitive to noise. We conclude that the method it offers strong robustness in case of the given perturbation and good tracking features.CONCLUSION
To control an under-actuated industrial process demands not only to know all the characteristic of his model, but also to achieve the desired performance when there exist load disturbances. Both Virtual Input and SVD approaches are applied to distillation column process, to prove the controller achievability for this industrial process and that it is the most efficient approach in terms of rapid response, set-point tracking and disturbance rejection. Through a comparative analysis, the Virtual Inputs approach avoids the complex calculation, such as calculate the inverse of the matrix, the controller structure is simple, has a few tuning parameters, and is easy to be accepted by operators. The simulation results show that our proposed method (Virtual Inputs) ensures suitably the set-point tracking and the disturbance rejection for the external disturbance as illustrated in distillation control problem suggested by Doukas and Lyben and gives better results compared to SVD approach. The modified internal model control scheme gives more degree of freedom in the controller design technique in order to improve the performance of the controlled system. As future works, it will be interesting to handle with unstable multivariable non-square systems and non-linear systems as proposed in order to design an internal controller with more flexibility and to apply the proposed Virtual Input controller to real processes.