#### Abstract

This paper presents a new nonlinear model predictive control (MPC) algorithm for Hammerstein systems subject to constraints on the state, input, and intermediate variable. Taking no account of constraints, a desired linear controller of the intermediate variable is obtained by applying pole placement to the linear subsystem. Then, actual control actions are determined in consideration of constraints by online solving a finite horizon optimal control problem, where only the first control is calculated and others are approximated to reduce the computational demand. Moreover, the asymptotic stability can be guaranteed in certain condition. Finally, the simulation example of the grade transition control of industrial polypropylene plants is used to demonstrate the effectiveness of the results proposed here.

#### 1. Introduction

Hammerstein systems consist of the cascade connection of a static (memoryless) nonlinear function followed from a linear dynamic system. Under certain assumptions (such as fading memory), the Hammerstein model could approximately represent nonlinear dynamics of real systems and has been successfully applied to many kinds of industrial processes such as pH neutralization [1], distillation [2], and polymerization transition [3]. In recent years, constrained control of Hammerstein systems has become one of the most needed and yet very difficult tasks for the process industry.

Model predictive control (MPC) is an effective control algorithm for handling constrained control problems, and various MPC algorithms have been proposed for control of Hammerstein systems with constraints [4β11]. Making use of the geometry structure of Hammerstein model, [4β9] developed two-step MPC schemes for Hammerstein systems with input constraints, where the intermediate variable was firstly obtained by linear MPC and then the actual control was calculated by solving nonlinear algebraic equations, desaturation, and so forth. Although two-step MPC algorithms possess more light computational burden than that where the nonlinearities were incorporated into optimal control problems directly [10, 11], solving nonlinear algebraic equations will inevitably have error and the restricted control is generally different from the desired one for constrained Hammerstein systems [6]. Therefore, the performance and stability properties may be deteriorated in the presence of constraints. Moreover, to the best of authorsβ knowledge, the conversional two-step MPC schemes usually have limited ability to dealing with input constraints. So, [12] presented a new two-step MPC scheme for Hammerstein systems with overall constraints, where the actual controller was determined by solving a finite horizon optimization problem that was described as tracking quadratic, optimal trajectories of the intermediate variable. However, the online computational demand of optimization problem of the scheme is still an open question.

In this paper, we consider the problem of the grade transition control of industrial polypropylene plants and propose an efficient NMPC scheme for Hammerstein systems with overall constraints based on the pole placement method [13]. We firstly obtained desired trajectories of Hammerstein systems by using pole placement to the linear subsystem and then calculated actual control actions by solution of finite horizon optimal control problems. In order to reduce the online computational burden of the optimization problem, only the first control action is calculated and the remaining is approximated by solving nonlinear algebraic equations over predictive horizon. Thus, the performance and stability properties can be guaranteed in the face of overall constraints, with moderate increment of the computational demand. Finally, an example of the grade transition control of propylene plants is used to demonstrate the effectiveness of the proposed algorithm.

#### 2. Pole Placement-Based NMPC Algorithm

Consider the discrete-time Hammerstein systems described by where , , and are the state, intermediate variable and control input, respectively. Vector field represents the nonlinear relationship between the input and intermediate variable, satisfying . It is assumed that is stabilizable and the state is available for feedback control. Without loss of generality, we also assume that the origin is the equilibrium of system (1).

Given constraints of system (1) as follows: where sets , , and are constraints on the state, intermediate variable, and input, respectively. The goal of this paper is to design an efficient NMPC algorithm for the Hammerstein system (1) with the constraint (2).

##### 2.1. Design of NMPC Algorithm

Consider desired poles, denoted by , of the linear subsystem of (1), which represent some desired performance of system (1). Associated with this set of desired poles is a linear controller of the intermediate variable which is generated via Lyapunovβs method without any constraints [13]. Since the controller (3) is unable to be implemented by real control action, here it is labeled as with associated closed-loop states .

Define a finite horizon objective function as follows: where the vector is -step ahead prediction from time instant , integer is the predictive horizon, and matrices and are the weighting matrices of the state and input, respectively. Then a new two-step MPC algorithm of Hammerstein systems with overall constraints is presented below.

*Algorithm 1 (modified two-step MPC, MTMPC). *This algorithm comprises the following steps:*Step 1. *Given desired poles and off-line compute the gain matrix in (3) such that the eigenvalues of are those specified in vector .*Step 2. *Set , , in the objective function in (4).*Step 3. *With and the current state , solve on-line the nonlinear algebra equations (5) without constraints to obtain desired control actions *Step 4. *With the current state , determine on-line the actual control action as follows:
where is defined by (4).*Step 5. *Implement , set , and go back Step 3.

*Remark 2. *For the algorithm here, only the first control action (i.e., ) is optimized subject to constraints at each time. So, the online optimization problem has only decision variables versus decision variables for general nonlinear MPC, which results in a significant reduction of the online computational demand since the computational demand, in general, grows exponentially with the number of decision variables [14]. Thus, the predictive horizon should be chosen reasonably large to ensure adequate performance. On the other hand, since only the first control is implemented in MPC scheme, it is expected that only enforcing constraints for the first prediction should give excellent results. In addition, the closed-loop performance of MPC systems clearly depends on how desired poles are arranged. The further work will focus on how to determine optimally.

It should be pointed out that the feasibility of optimization problem (6) may not be guaranteed theoretically with respect to arbitrary constraints in (2). However, it is usually ensured in practice via adjusting the intermediate variable constraint. Next, we show that the closed-loop MPC system (1) is asymptotically stable, provided that the feasibility of (6) is hold.

##### 2.2. Analysis of Stability

In terms of Algorithm MTMPC, the MPC law of constrained Hammerstein system (1), (2) is defined as with the closed-loop system

*Definition 3. *Set is called an attractive region of system (8) if for all , the system trajectories when and satisfies
where subset .

Let a set of points be the desired poles of the linear subsystem . Then from inverse Lyapunovβs theorem [13], there is systematic matrices and such that In following, we present the result on our proposed algorithm.

Theorem 4. *Consider the system (1), (2), there exists a nonempty region such that the closed-loop system (8) is asymptotically stable if the inequality (11) is satisfied
**
where function and with . Moreover, the set is an attractive region of system (8).*

*Proof. *Consider the system (1), (2), define a level set of function as
such that , for all , . Note that the set is always nonempty since the constraint in (2) includes the origin as the internal point.

Let function be a Lyapunov function candidate of system (8). Then, for any , we have
where . Substituting (10) into (13) yields
From the inequality (11), it is straight forward to obtain that
Hence, function is a Lyapunov function of system (8), and the system is asymptotically stable in egion . Furthermore, from the definition of attractive regions and set in (12), we derive that is an attractive region of system (8).

This establishes the theorem.

*Remark 5. *From the proof of the above theorem, we know that the stability property of closed-loop system (8) is guaranteed by the linear controller of the intermediate variable, provided that the feasibility of optimization problem (6) and condition (11) holds. This suggests that the stability and optimality can be separated to some extent. Therefore, the design parameter of the objective function in (4) can be tuned freely to achieve more satisfactory performance, regardless the stability of the closed-loop system. In addition, it is pointed out that the attractive region obtained here is merely a theoretical concept and the computation of is not an easy task.

*Remark 6. *Instead of conventional two-step MPC algorithms, where the actual control was determined by solving nonlinear algebraic equations based on desired controllers of the intermediate variable resulted from linear MPC, actual control actions here are derived by solving a finite horizon optimal control problem (i.e., nonlinear MPC), which is defined as tracking the desired state trajectories driven by pole placement state feedback controller. Thus, overall constraints of Hammerstein systems are taken into account when designing the predictive controller, and then performance and stability properties are guaranteed. Meanwhile, the ability to handling constraints is enhanced in the algorithm proposed here. Note that the linear controller of intermediate variables in (3) can be any one that stabilizes the linear subsystems of Hammerstein models with certain desired performance.

#### 3. Simulation Example

Taking an example of polypropylene (PP) grade transition control, we illustrate the effectiveness of the proposed algorithm.

In propylene polymerization industry, melt index (, g/10βmin) is usually used to identify grade indices of PP homopolymer and both of and concentration of ethylene (, %) in PP denote grade indices of total phase PP [12, 15β18]. This implies that the nature of grade transition process is to operate the response of the variables and in polymer. Then according to the propylene polymerization mechanism, we have the following mathematical model which formulates the dynamics of the grade transition in polypropylene process [12]: where states and are the cumulative (, g/10βmin) and (, %) of PP, respectively; intermediate variables and denote the instantaneous (, g/10βmin) and (, %) of PP, respectively; input, , , denote the reaction temperature , concentration ratio of hydrogen to propylene (%), and concentration ratio of ethylene to propylene (%), respectively, and parameter is the average polymer residence time ( = 2βh). The parameters of model () can be identified on-line by industrial data [15].

The grade transition sequence considered here is ABC, where grades A and B are PP homopolymer and C is total phase PP. It should be pointed out that the feed flow rate of ethylene is always zero (i.e., ) in production of PP homopolymer grades A and B. The grade transition process is defined by the set-point of states and input/output constraints as shown in Table 1, where superscripts , , and denote accumulation, instantaneous value, and increment constraints, respectively. Each grade transition is optimized over 10βh horizon with a sampling interval of 0.5βh.

In the simulation running, the closed-loop desired poles of system (16) are placed at , and the corresponding state-feedback matrix in (3) is computed by . Moreover, the parameters in the objective function (4) are chosen as penalty on the states, , with , penalty on the inputs, and a horizon length of . Assume that grade transition AB takes place at the tenth hour and BC at the thirtieth hour. Then some curves of the grade transition process being controlled by NMPC law are shown in Figures 1 and 2, where the solid line in Figure 1 is the set-point, the dotted line the instantaneous, and the dashed line the cumulative.

**(a) Melt index**

**(b) Concentration of ethylene in polypropylene**

**(a) Reaction temperature**

**(b) Concentration ratio of hydrogen to propylene**

**(c) Concentration ratio of ethylene to propylene**

In order to calculate the quantity of off-specification polymer in grade transition process, we employ the criterion where on-specification polymer is defined as the polymer with and within a Β±5% deviation of desired specification in Table 1. All other polymers beyond the criterion are categorized as off-specification production. Similarly, the transition time is defined as the time interval when the cumulative properties of polymer are out of the specification range Β±5% of the previous grade and lasting until the cumulative properties enter and stay within the specification range Β±5% of the new target grade. Therefore, the time of grade AB transition is 8.5 hours and the time of grade BC transition is 6.5 hours. Compare to the time of 10 hours in actual grade transition operation, this decreases the quantity of off-specification polymer in grade transition process and then increases the productive time of on-specification polymer. Also, Figure 2(a) suggests that there is no change in reaction temperature during the grade transition processes, which is desired in actual operation. Finally, it can be seen that the control laws, intermediate variables, and the states in the grade transition operation do not violate the constraints in Table 1.

#### 4. Conclusion

Together with the pole placement method, the paper presented an efficient MPC algorithm for control of Hammerstein systems subject to constraints on states, inputs, and intermediate variables. Applying pole placement to the linear subsystem of Hammerstein models yielded linear controllers of the intermediate variable and then the actual control action was determined by online solving an optimization problem. The size of the online optimization problem depended only on the number of inputs, not on the predictive horizon, which greatly reduced the computational demand of the algorithm. With condition of feasibility of the optimization problem, the asymptotical stability of the closed-loop system with constraints was guaranteed and the results on simulation example of grade transition control of polypropylene plants showed the effectiveness of the results obtained here.

#### Acknowledgments

This work is supported by the National Natural Science Foundation of China under Grant no. 60904040, Specialized Research Fund for Doctoral Program of Higher Education of China under Grant no. 20093317120002 and Natural Science Foundation of Zhejiang Province under Grant no. Y1100911. Finally, the authors are grateful to the reviewsβ constructive comments. This paper includes revised and extended text of the first authorβs presentation at the 29th CCC, July 27, 2010, Beijiang, China.