Design of PI Controllers for Hydraulic Control Systems

The paper proposes a procedure for design of PI controllers for hydraulic systems with long transmission lines which are described by models of high order. Design is based on the combination of the IE criterion and engineering specifications (settling time and relative stability) as well as on the application of D-decomposition. In comparison with some known results, the method is of graphical character, and it is very simple (solving nonlinear algebraic equations is eliminated). The paper presents the algorithm of software procedure for design of the controller. The method is compared with other methods at the level of simulation, and its superiority is shown. By applying theNyquist criterion, it is shown that themethod possesses robustness in relation to nonmodelled dynamics.


Introduction
The first heuristic rules for selection of parameters of the PI controller were given by Ziegler and Nichols in 1942 [1].Since then, numerous papers have been devoted to determination of parameters of the controller, according to different types of assumed controlled processes.It is estimated that more than 95% of controllers used in industry are PID controllers, and most of them are PI controllers [2].The performed analysis connected with the efficiency of control contours in industrial production has shown that PI controllers are used to a large extent, but they are commonly tuned in a bad manner [3,4].According to [4], only 20% of controllers in industry operate well, while 30% of them have poor performance due to the incorrectly tuned controller.In [5], based on research in process industry, the managers and engineers reach a common conclusion that tuning of parameters of PI and PID controllers is an important and difficult problem that deserves more attention.According to [6], the PI controller is the most frequently used control algorithm in process industry.The reason for that is its relatively simple structure, which can be easily understood and implemented in practice.Despite its wide application, there is still a need for improved tuning of parameters of the PI controller [7].The heuristic method proposed by Ziegler and Nichols has the great advantage because it requires very little information about the process.However, this method has significant disadvantages, and it gives very bad damping (it is usually  = 0.2) [7].Small damping leads to degradation of the system and very bad performances from the aspect of relative stability and robustness.
Most PI controllers in process industry operate in the regime of set point value.Therefore, it is of principal importance to efficiently solve the problem of load disturbance rejection [8].Reference [9] considers optimal tuning of controllers.It is shown that the IE criterion (integrated error) is in direct relation with integral gain of controllers.Those results are further expanded in [10], where it is proposed to have optimal load disturbance rejection with constraints on the function of system sensitivity and the weighting coefficient of the given value.The problem is reduced to solv nonlinear algebraic equations.The Newton-Raphson method is used for solving these equations.The initial conditions are the critical factor in this method.
Reference [11] presents a simple procedure for tuning of the PI controller.The problem is reduced to minimization of the IE criterion with the constraint that the Nyquist curve of the loop transfer function has a tangent in its left semiplane, parallel to an imaginary axis.The result of the procedure is explicit formulas for proportional and integral gains.Reference [12] proposes a new procedure for tuning of ideal PID controllers in series with the first-order noise filter.
The problem is reduced to solving two nonlinear algebraic equations.
Reference [13] gives the procedure for synthesis of the P controller based on the generalized Nyquist criterion for hydraulic systems with long transmission lines described by models of high order.Such a type of controllers cannot eliminate the static error of control.That weakness is eliminated by the PI controller.The method that enables design of the PI controller for systems of high order is Ddecomposition.The D-decomposition method was developed by Neȋmark [14,15].It was considerably expanded by Mitrović establishing a strong link between the values of tuning parameters of the characteristic polynomial and the appearance of the transition process expressed through the corresponding degree of relative stability of the system [16].The D-decomposition method was fully generalized in the algebraic method developed by Siljak in [17][18][19].For efficient interpretation of results obtained by the D-decomposition method it is also necessary to have appropriate graphical interpretation which, for the systems of high order, was not possible without the corresponding software support.Intensive development of computer technique has renewed the interest in this methodology [20].Let us emphasize that modern control theory allows design of controllers, for processes of high order, whose order will be equal to the order of the process [21][22][23][24].Implementation of such controllers in industry is very complex and expensive.From that aspect, there is a need for design of controllers of low order, which are most frequently present in industrial practice, for control of processes of high order [25][26][27].
This paper proposes the procedure for design of the PI controller for systems of high order.The results from [9], which establishes a direct relation between the IE criterion and the integral gain (the higher the integral gain, the smaller the value of the IE criterion), are the starting point.The result is extended by introducing engineering specifications (settling time and relative stability).It results in a simple and efficient procedure for design of the PI controller for systems of high order.The method is basically graphical, and unlike [10,12], it is not necessary to solve any nonlinear algebraic equations.The complete methodology of design is given in the second chapter of this paper.The third chapter gives the results of simulation for the hydraulic system with long transmission lines.A MATLAB program was created for graphical interpretation of results, and its flow is presented by the algorithm given in this chapter.The fourth chapter presents the results of comparative analysis of the proposed method for design of the PI controller and other methods, showing the superiority of the method proposed in this paper.

Design of the PI Controller for Systems of High Order
In [9], the problem of disturbance rejection is reduced to minimization of the following IE criterion:  As the result, the following relation is obtained: The result is extended in [10] by introducing constraints on the sensitivity function and the weighting factor of the set point value in the optimization problem.The problem belongs to the category of nonconvex optimization and is reduced to solving nonlinear algebraic equations by iterative methods.
This paper also considers the IE criterion, but instead of constraints from [10], engineering constraints are introduced on (i) relative stability (damping of the closed loop) and (ii) settling time.
The solution of the problem lies in the intersection between the previous two conditions.In accordance with the minimum of the IE criterion, the additional condition is that the integral gain should have the maximum value.This is the basis for development of a simple graphical method based on D-decomposition.
The transfer function of the PI controller is The transfer function of the process is represented in the form: The characteristic equation of the automatic control system from Figure 1 is determined by the equation: By connecting (7) and ( 8), the final expression for the characteristic equation of the automatic control system in the complex domain is obtained as follows: Taking into account (9), it is necessary to express the complex number  in a suitable form and use it for establishing the relation between the damping degree  and the variable parameters of the controller,   and   , contained in the characteristic equation ( 9) for the automatic control system.This is how the area from the "" plane below the straight line  = const.(Figure 2) is mapped in the area of the corresponding damping coefficient represented by the curve  = const., in the parameter plane of tuning parameters of the controller (  ,   ) as follows: where By connecting ( 9) with ( 12), the characteristic equation of the automatic control system obtains the form: where (,   ) and (,   ) represent the real and imaginary parts of the polynomial  1 (,   ).The real and imaginary parts of the polynomial  1 (,   ) are obtained by connecting ( 8) with ( 12) as follows: where   and   are Chebyshev functions of the first and second kinds for which the following recurrent equations hold: where (,   ) and (,   ) represent the real and imaginary parts of the polynomial (,   ), and they are determined based on the following equations: By connecting ( 13) with ( 14) and ( 18), it is obtained that By separating the real and imaginary parts from (20) and by solving it, the following system of equations is obtained: Mathematical Problems in Engineering By solving the system of equations at   ̸ = 0, 0 ≤  < 1 the expressions for the parameters   and   of the PI controller are obtained as follows: For   = 0, the system of ( 21) gives the solution: Equation ( 25) represents a singular straight line which, together with the curve described by (22), represents the closed contour of the region of possible solutions, in the plane of variable parameters   ,   shown in Figure 3.
For  = 1, the system of ( 21) is reduced to one equation described by Equations ( 22) through (26) completely define the parameter (  ,   ) plane for the selected value of the damping coefficient  = const., 0 ≤  ≤ 1, at the change of undamped frequency   = (0,  max ).The parameter plane with its stability limits is presented in Figure 3.
From Figure 3, three characteristic areas particularly stand out for the following: (i)  = 0, the system is at the stability limit and has oscillatory character, the area below the curve  = 0 is the stability area, and the area above the curve is the nonstability area; (ii)  = 1, the system has critical damping, and the area below the curve  = 1 is the area in which the system has monotonous (aperiodic) character.The critically damped response is desired in a lot of applications because it represents the fastest aperiodic response possible; (iii)  = const., the curve for the required degree of system damping from which the variable parameters of the controller (  ,   ) can be chosen.gain   .Selection of the maximum integral gain will not always yield good results regarding robustness defined by phase margin (  ), gain margin (  ), and robustness measure   [27].This problem is particularly noticeable in systems of high order.In the character of response of the closed loop, it is reflected in the increased overshoot () and the increased settling time   .
In order to solve this problem, the area of selection of parameters of the controller (  ,   ) has been narrowed by taking out the area in which the system will possess good settling time (  ) on the corresponding curve  = const.
In order to have a system with good settling time, it is necessary that all real parts of the poles of the transfer function of the closed loop should have a location to the left of the straight line   = const.(from Figure 2).The area from the "" plane to the left of the straight line   = const.(Figure 2) is thus mapped in the area of the corresponding settling time  min = const., in the parameter plane of tuning parameters of the controller (  ,   ) as follows: Based on ( 12) and ( 27), it is clear that Here it is necessary to find the minimum value of undamped frequency that will satisfy the condition: In order to have all real parts of the complex number  = − to the left of the straight line   , it is necessary to find the variable parameters (  ,   ) that will satisfy those conditions.This is realized by replacing the complex number  = − min in (9), so that the characteristic equation now obtains the form: The value of undamped frequencies   min and   max , whose location is identified in Figure 2, is read from the graphs   (,   ) = 0 and   (,   ) = 0, as it is shown in Figure 4.The graphical interpretation of (30),  min = const., in the parameter plane (  ,   ) of the required damping coefficient  = const., is shown in Figure 5.
The straight line  min = const.intersects the curve of the required damping coefficient  = const.in the parameter plane at the point  with the coordinates (  ,   min ).At the point , the integral gain has a minimum value.The point  with the desired coordinates (  ,   ) lies at the part of the curve  = const.between the points  and , at the interval of undamped frequencies   = (  min ,   max ).The lower the order of the process, the closer the point  to the point (  ,   max ), and the higher the order of the process, the closer the point  to the point (  ,   min ).In this narrow range, by selecting the appropriate step, it is very easy to find the point  with the parameters of the controller that will guarantee quality operation of the automatic control system.

Control of a Process with a Long Transmission Line
In order to show the efficiency of the proposed method of design of the PI controller, simulation in the program package MATLAB for the transfer function of the process  1 () has been carried out as follows: (31) The transfer function described by (31) represents a mathematical model of a real process used in industry.It is a pump controlled hydromotor, where the variable flow pump and the hydromotor of constant flow are connected by means of a long transmission line.The complete identification of this transfer function was done in paper [13].The motive for taking exactly this transfer function of the process for testing the proposed method is its high order-it is of the tenth order.When the system is of high order, as in the case of the transfer function  1 (), then the poles which are closest to the imaginary axis in the complex plane have the decisive influence on the character of the transition process.The proposed methodology allows location of such poles of the transfer function at a safe distance from the imaginary axis; that is, it allows the closed loop of the system to have all its poles to the left of the straight line   , which is presented in Figure 2.
Based on the programme created in MATLAB and according to the procedure proposed in the paper, the parameters of the PI controller can be determined for any transfer function of the process described by (4), so that the closed loop of the system could possess the required damping coefficient and good settling time.
The flow of the programme written in MATLAB is represented by the algorithm shown in Figure 6.
The result of operation of the programme is the graphical interpretation of parameters of the controller shown in Figure 7, for different values of the damping coefficient 0 ≤  ≤ 1 of the closed loop of the automatic control system.The stability limit is marked by the curve  = 0, from which it is possible to read the value of the critical proportional gain  16) and ( 17)) Data input (4) (  ),   (  ) (( 22)-( 26))   = 30.39.The values of the maximum integral gain   max and the corresponding proportional gain   for different values of the damping coefficient  are read from Figure 7.The values of parameters of the controller obtained from Figure 7 are shown in Table 1.
The responses of the automatic control system presented in Figure 8 are recorded for the designed values of parameters of the PI controller from Table 1 as well as for the critical value of the proportional gain   = 30.39 at  = 0. From Figure 8 it can be seen that high overshoots and high values of settling time are obtained in the step response, even for higher values of the damping coefficient ( = 0.6, 0.8, and 1.0) of the closed loop.The values of overshoot and settling time are shown in Table 1.
Figure 9 shows how the designed controller, in accordance with the minimum of the IE criterion, rejects the load caused by the action of disturbance  = 1 at the value of reference   = 0. From Figure 9 it is possible to conclude that good results from the aspect of load disturbance rejection can be expected for the values of the damping coefficient  ≥ 0.4.Figure 9 also shows that in an uncontrolled process the load of the system which is caused by disturbance will not Step response 2.5  be rejected.However, according to this criterion, satisfactory values of overshoot, settling time, and stability margin will not be obtained, which is best seen from Table 1.
In order to eliminate the undesired values of increased overshoot and settling time, for the chosen damping coefficient  = 0.6, the settling time has been estimated according to  min =   min .The values of the minimum undamped frequency   min and   max are read from Figure 10.
The minimum value of the integral gain   min = 29.75 and the corresponding proportional gain   = 10.73 is thus obtained in the parameter plane, as it can be seen from Figure 11.Going along the part of the curve shown in Figure 11, from the point of minimum integral gain toward the point of maximum integral gain, with the appropriate step (automatically generated by the program) and recording the response, it is very easy to reach the point on the curve with the parameters of the controller that will give best performances for the given process.
Figure 12 shows the step response of the closed loop with the parameters of the controller read from Figure 11.It can be seen from Figure 12 that for the designed values of parameters of the controller   = 10.78 and   = 57.15, the closed loop of the automatic control system possesses good performances (the overshoot of 9.98% and the settling time of 72.8 ms).

Comparison with Other Methods
The PI controller designed by means of the methodology presented in this paper for the process described by (31) has been compared with the Ziegler-Nichols (ZN) [1] and Tyreus-Luyben (TL) [28]   of parameters of controllers in industrial practice.The results of comparison in the step response for the described process are shown in Figure 13.The complete performances from the aspect of response, relative stability and robustness are presented in Table 2.
Based on the results presented in Table 2, it can be seen that the Ziegler-Nichols method gives a very oscillatory response characterized by the overshoot of 62.3% with the settling time of 148 ms and very bad robustness (phase margin 28 ∘ , gain margin 1.89, and stability margin 0.48).The Tyreus-Luyben method gives a well-damped response characterized by the overshoot of 13.6%, but with a very high value of settling time of 195 ms.From the aspect of robustness, the Tyreus-Luyben method gives good results.The method proposed in this paper gives a well-damped response characterized by the overshoot of 9.98% with considerably better settling time of 72 ms in comparison with the ZN and TL methods.From the aspect of robustness, the   proposed method also gives a better result (phase margin 56.5 ∘ , gain margin 2.74, and stability margin 0.64) too.
The results of comparison from the aspect of load disturbance rejection are presented in Figure 14.From Figure 14 it is possible to see that the proposed method for the values of the damping coefficient ( = 0.5 and  = 0.6) quite satisfactorily rejects the load caused by the action of disturbance  = 1 in comparison with the Ziegler-Nichols method and does it considerably better in relation to the Tyreus-Luyben method.
Remark.It is well known that the Ziegler-Nichols frequency method cannot be applied to the processes of the first and second orders because such processes cannot be brought to the oscillatory stability limit by changing the proportional gain.Also, the Ziegler-Nichols method and some other similar methods cannot be applied to unstable processes.The proposed methodology does not have any restrictions referring to the order of the system.It can be very efficiently applied to low-order processes.Another advantage of this method is that it can also be applied to unstable processes of low order.

Conclusion
This paper has developed an efficient and simple graphical methods for design of the PI controller, which achieves high performances for a broad range of linear processes.For design, it is necessary to know the transfer functions of the process.The processes of high order have been considered in the paper.In comparison with the procedures for tuning of the PI controller proposed in literature, the method described in this paper is characterized by great simplicity and clear engineering specifications.The results of simulations show good robustness in relation to unmodelled dynamics as well as superiority over some other methods of tuning of controllers.The proposed method is suitable for on line real-time implementation and for auto tuning of the PI controller.The methodology of design of the PI controller presented in this paper is simple and can easily be adopted by industry.The extension of this method to design of the PID controller is underway.

𝜉 = c o n s t a nFigure 2 :
Figure 2: Area with the required settling time and relative stability.

Figure 3 Figure 3 :
Figure 3 allows reading the maximum value of integral gain at the point ,   =   max , and the corresponding proportional

Figure 4 :Figure 5 :
Figure 4: Change of the proportional and integral gains of the controller as a function of undamped frequency for the required damping coefficient.

Figure 6 :
Figure 6: Algorithm for determination of parameters of the PI controller.

Figure 7 :
Figure 7: Parameter plane for different values of the damping coefficient with the stability limits.

Figure 8 :Figure 9 :Figure 10 :
Figure 8: Response of the system for the reference value   = 1 at   max , for changes of the damping coefficient 0 ≤  ≤ 1.

Figure 11 :
Figure 11: Parameter plane for the required damping degree and required settling time.

Figure 12 :
Figure 12: System response at the defined damping coefficient  = 0.6 and the required settling time.

Figure 13 :
Figure13: Comparative presentation of the system responses, to the action of reference   = 1, for the proposed method (PM), the Ziegler-Nichols (ZN) method, and the Tyreus-Luyben (TL) method.

Figure 14 :
Figure14: Comparative presentation of the system response to the action of load disturbance  = 1, for the proposed method (PM), the Ziegler-Nichols method (ZN) and the Tyreus-Luyben method (TL).

Table 1 :
Results of design of the PI controller with performances.
methods for tuning of parameters of the PI controller.Those are two very frequent methods of tuning

Table 2 :
Comparative presentation of the results of design of the PI controller with performances for three methods.