Mathematical modeling of first order process with dead time using various tuning methods for industrial applications

Sathish Kumar Shanmugam1 , Yuvaraj Duraisamy2 , Meenakumari Ramachandran3 , Senthilkumar Arumugam4

1Jansons Institute of Technology, Coimbatore, Tamilnadu, India

2Tamilnadu College of Engineering, Coimbatore, Tamilnadu, India

3Kongu Engineering College, Erode, Tamilnadu, India

4Chettinad College of Engineering and Technology, Karur, Tamilnadu, India

1Corresponding author

Mathematical Models in Engineering, Vol. 5, Issue 1, 2019, p. 1-10. https://doi.org/10.21595/mme.2019.20447
Received 12 December 2018; received in revised form 2 January 2019; accepted 27 January 2019; published 31 March 2019

Copyright © 2019 Sathish Kumar Shanmugam, et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Abstract.

The proposed research involves Mathematical modeling for first order process with dead time using various tuning methods for industrial applications. Different tuning methods are proposed. Proposed method selection depends on plant operating conditions and also depending upon the process dynamics. The PID controller is most widely used for industrial process control. Modeling were developed for modified internal control model [4-6] in this proposed research. The proposed work is the modeling and simulation of three different first order processes with dead time. The standard controller tuning method is used to obtain the steady state response of first order with dead time.

Mathematical modeling of first order process with dead time using various tuning methods for industrial applications

Highlights
  • The controller tuning is performed to stabilize the process output according to the set point changes.
  • The controller tuning is based on the error signal and it is performed to nullify the error.
  • The feed-forward controller takes control action according to the load changes.
  • Dead time is the time taken by the system to produce the response after the input is applied

Keywords: IMC, tuning, transfer function, simulation, analysis.

1. Introduction

The preferred output of a control system is obtained by tuning the controller parameters according to the set point. The first order system is widely used in industrial control system to represent the behavior of process dynamics [1, 2]. The delay takes place in process measurement, controller action, actuator operation and computation is called delay time. Dead time is the time taken by the system to produce the response after the input is applied [3].

2. Block diagram of feedback control system

Fig. 1 shows the block diagram proposed IMC Feedback Control System. The algorithm is necessity oe for tuning. The control action is done by the actuator. The input signal which is attenuated will control the parameter. Sensor is used for controlling the signal measure the process variable from the process output. The control algorithm is used to nullify the error signal to zero [4-6].

Fig. 1. Block diagram of proposed IMC feedback control system

 Block diagram of proposed IMC feedback control system

3. Literature review

• Astrom and Haggland proposed an alternative method of the system. The technique itself create its own test signal based on the estimation [7].

• The robustness of the automatic tuning of controller by Zeigler’s technique is improved by Astrom and Haggland method. In this robust loop shaping rule called MIGO and AMIGO is introduced. MIGO concentrates on maximizes the integral gain to improve the robustness of the system.

• Wallen et al., enhanced the loop shaping based tuning by introducing an additional constraint the ratio of integral time to the derivative time for tuning the PID controller. And also, it makes some conditions the curvature should be negative and it should have a decreasing phase of a curve.

• Zeigler and Nichols proposed the first auto tuning procedure for controllers. The sensitivity, Reset rate and pre-act time are the characteristics of the reaction curve for designing an automatic controller [8, 9].

4. Materials and methods

4.1. Controller tuning

When PID controller has been selected, Different methods are required for tuning. This is called as ‘controller tuning’. Different methods adopted for the proposed tuning are analysed and its results were carried out with the modified structure. The method proposed is Internal Model Control (IMC).

Fig. 2. Proposed IMC schematic circuit

 Proposed IMC schematic circuit

Fig. 3. Proposed IMC schematic circuit with surface change

 Proposed IMC schematic circuit with surface change

Fig. 4. Modified proposed IMC schematic circuit with surface change

 Modified proposed IMC schematic circuit with surface change

Fig. 5. Modified proposed IMC schematic circuit with surface change inner loop

 Modified proposed IMC schematic circuit with surface change inner loop

The standard feedback of proposed IMC circuit is given by:

(1)
g c s = q S 1 - g p s   q S .

The equivalent circuit equation of proposed structure IMC is given by:

(2)
g c s = q S 1 - g ~ p s   q S .

The PID controller tuning equation is given by:

(3)
g c s = k c τ 1 τ D s 2 + τ 1 s + 1 τ 1 s 1 τ F s + 1 .

The proposed method first order process is given by:

(4)
g ~ p s = k p τ p s + 1 .

After filtering in the process is given by:

(5)
q s = q ~ s f s = g p - 1 s f s = τ p s + 1 k p 1 λ s + 1 ,       q s = 1 k p τ p s + 1 λ s + 1 ,

After transforming the equation can be written as:

(6)
g c s = q S 1 - g ~ p s   q S = τ p s + 1 k p ( λ s + 1 ) 1 - k p τ p s + 1 τ p s + 1 k p ( λ s + 1 ) = τ p s + 1 k p λ s .

Multiplying with tuning parameters:

(7)
g c s = τ p k p λ τ p s + 1 τ p s .

Similarly, for second order process:

(8)
g c s = τ 1 + τ 2 k p λ τ 1 τ 2 s 2 + ( τ 1 + τ 2 ) s + 1 ( τ 1 + τ 2 ) s .

4.2. Methods of tuning for servo response

The controller tuning is performed to stabilize the process output according to the set point changes. The feedback controller is used for setting point tracking. The error signal is generated which is the differences between setpoint and process variable. The controller tuning is based on the error signal and it is performed to nullify the error. The controller tuning is performed to stabilize the process output according to the load disturbance. The feedback controller is not suitable for load disturbance tracking. The advanced control mechanism such as feedforward controller is used to track the load disturbance. The feedforward controller takes control action according to the load changes.

4.2.1. Cohen-Coon tuning

This method is also based on a delayed first order rise, and the method tunes the PID gains to achieve quarter-amplitude damping, i.e. each peak in the transient is one quarter of the immediate preceding peak. If the unit step response of the open loop plant is available, then process gain (K), dead time (td) and time constant (T) can be determined. PID controllers are designed to test the system and selected based on the performance and selection criteria. Table 1 shows the Cohen Coon method.

Table 1. Cohen Coon

Controller type
(KC)
T I
T D
PID
  1.49 k m τ m
2τm
1.27 τ m

4.2.2. Ziegler Nichols tuning

Ziegler Nichols method is also called as Ultimate Cycle Method. It is based on adjusting a closed loop until steady state oscillations occur. This method is used to minimizing the total error of the response. To return to desired level operation as soon as possible and to keep the maximum deviation as small as possible this method is used. The integral absolute error is the integral error which defines the stability of the system. Larger the integral error, smaller the stability is. Thus, by minimizing the integral absolute error we get the good robustness and performance. Thus, IAE is also taken measuring parameter. Control objectives focus on the behavior of the system such as to produce a zero steady state error, a fast transient response to a step command, a short settling time and low overshoot. It is also desirable to make the system less sensitive to disturbance. The period of these oscillations is defined as the ultimate period Pu. The controller gains at which these oscillations occur is referred to as the ultimate gain Ku. Ziegler Nichols method are shown in Table 2.

Table 2. Ziegler Nichols method

Controller type
K C
T I
T D
PID Zeigler Nichol’s
K u 1.7
P u 2
P u 8
Ultimate cycle PID developed by mc Millan
1.009 T m k m τ m 1 1 + T m τ m         . 65 2
2 τ m 1 + T m τ m         . 65
0.26TI

4.2.3. FOPTD transfer functions

The PID tuning formulas are applied for the three FOPTD transfer functions to obtain the analysis. The FOPTD transfer function is considered based on ratio of dead-time to the time constant less than one, equal to 1 and greater than 1. Z-N method is time consuming and forces the system to margin if instability. Many other algorithms have been proposed to solve these problems by obtaining critical data (ultimate gain and frequency) under more acceptable conditions. One of these methods is damped oscillation method. Z-N method are shown in Table 3.

Table 3. Z-N method

Controller type
K C
T I
T D
PID
K u / 1.3
2.2 P u
P u / 5.2

The implicit manner based PID control algorithm is simple and robust to handle the model inaccuracies and hence using IMC-PID tuning method a clear trade-off between closed loop performance with damped oscillation and robustness to model inaccuracy is achieved with a single tuning parameter are show in Table 4 for implicit manner.

Table 4. For implicit manner

Controller type
K C
T I
T D
PID
2 / K ( λ + θ / 2 )
2λ+θ
λ θ + ( θ 2 ) / 4 2 T + θ

4.2.4. Integrating process transfer functions

The PID tuning formulas are applied for the three different order transfer functions to obtain the analysis in order to study the integrating process performance. The different order integrating process transfer functions considered for the analysis are.

The PID tuning formulas are applied for the three different order transfer functions to obtain the analysis in order to study the integrating process performance. The different order integrating process transfer functions considered for the analysis are,

(i) Transfer function 1.

The pure integrator process transfer function with low delay consider as example of benchmark process default value of second order system for controlling:

(9)
G p 1 s = 1 s e - 2 s .

(ii) Transfer function 2.

The pure integrator process transfer function with low delay consider as example of benchmark process default value of second order system for controlling:

(10)
G p 2 s = 1 s e - 5 s .

(iii) Transfer function 3.

The higher order integrator process transfer function without delay consider as example of benchmark process default value of second order system for controlling:

(11)
G p 3 s = ( 0.17 s + 1 ) 2 s s + 1 2 ( 0.028 s + 1 ) .    

The transfer function is approximated to integrator process with delay by using half rule approximation in order to work with, and the transfer function is example of benchmark process default value of second order system for controlling:

(12)
G p 3 s = 1 s e - 1.69 s .

5. Results and discussions

5.1. MATLAB software for simulation of FOPDT

Simulink is widely used for analysis of first order process with dead time. The Simulink is used in this proposed work to obtain the response of first order with dead time. Simulink block diagram of FOPTD set point and transfer function for load disturbance are shown in Figs. 6 and 7.

Fig. 6. Simulink block diagram for FOPTD for set point

 Simulink block diagram for FOPTD for set point

Fig. 7. Simulink block diagram of FOPTD transfer function for load disturbance

 Simulink block diagram of FOPTD transfer function for load disturbance

5.1.1. Transfer function 1

The ratio of time delay θ to the time constant τ is less than one:

(13)
G p 1 s = 1 5 s + 1 e - s .

Fig. 8 shows the output of transfer function 1. Fig. 11 shows the output of transfer function 1 setpoint.

Fig. 8. Output of transfer function 1

Output of transfer function 1

5.1.2. Transfer function 2

The ratio of time delay θ to the time constant τ is unity:

(14)
G p 2 s = 4 10 s + 1 e - 10 s .

Fig. 9 shows the output of transfer function 2. Fig. 12 shows the output of transfer function 2 setpoint.

Fig. 9. Output of transfer function 2

 Output of transfer function 2

5.1.3. Transfer function 3

The ratio of time delay θ to the time constant τ is greater than one:

(15)
G p 3 s = 1 7 s + 1 e - 14 s .

Fig. 10 shows the output of transfer function 3. Fig. 13 shows the output of transfer function 3 setpoint.

Fig. 10. Output of transfer function 3

 Output of transfer function 3

Fig. 11. Output of transfer function 1 (set point)

 Output of transfer function 1 (set point)

Fig. 12. Output of transfer function 2 (set point)

 Output of transfer function 2 (set point)

Fig. 13. Output of transfer function 3 (set point)

 Output of transfer function 3 (set point)

5.2. Simulation results of closed loop methods

The various output loop response is shown in Figs. 14 to 18.

Fig. 14. Ziegler Nichols closed loop response

 Ziegler Nichols closed loop response

Fig. 15. IMC closed loop response

 IMC closed loop response

Fig. 16. Integral of square error closed loop response

 Integral of square error closed loop response

Fig. 17. Integral of absolute error closed loop response

 Integral of absolute error closed loop response

Fig. 18. Integral of square time error closed loop response

 Integral of square time error closed loop response

6. Conclusions

The first order with delay process is widely used to represent the process dynamics. The selection of controller and tuning methods are important consideration in the control system design. This proposed work is the tuning and analysis of three different first order with dead time. This analysis shows that the performance is best for the ratio less than 1 for Zeigler Nichols, Cohen coon, ITAE load disturbances and IMC, ITAE, ISE, ISTE, GP set point method. This model produces a better transient response. This proposed work can be extended by using Fuzzy Logic Controller for first order with dead time and performance characteristics can be compared.

References

  1. Ozipneci Burak, Tolbert Leon M. Simulink implementation of induction machine model – a modular approach. IEEE International Electric Machines and Drives Conference, 2003. [Publisher]
  2. Dorf C. R., Bishop R. F. Modern Control System. 10th Edition, Prentice Hall, 2005. [CrossRef]
  3. Dumitrescu A., Fodor D., Jokinen T., Rosu M., Bucurencio S. Modeling and simulation of electric drive system using Matlab/Simulink environment. International Conference on Electric Machines and Drives (IEMD), 1999, p. 451-453. [Publisher]
  4. Le Huy H. Modeling and simulation of electrical drives using Matlab/Simulink and power system block set. 27th Annual Conference of the IEEE Industrial Electronics Society, 2001. [CrossRef]
  5. Nizam Kamarudin M., Rozali Sahazati Md. Implementation of cascade control DC motor. International Conference on Power and Energy, 2008. [CrossRef]
  6. Sondhi S., Hote Y. V. Fractional order PID controller for load frequency control. Energy Conversion and Management, Vol. 85, 2014, p. 343-353. [Publisher]
  7. Debbarma S., Dutta A. Utilizing electric vehicles for LFC in restructured power systems using fractional order controller. IEEE Transactions on Smart Grid, Vol. 8, Issue 6, 2017, p. 2554-2564. [Publisher]
  8. De Keyser R., Muresan C. I., Ionescu C. M. A novel auto-tuning method for fractional order PI/PD controllers. ISA Transactions, Vol. 62, 2016, p. 268-275. [Publisher]
  9. Moosavi S. H. S., Bardsiri V. K. Satin bowerbird optimizer: a new optimization algorithm to optimize ANFIS for software development effort estimation. Engineering Applications of Artificial Intelligence, Vol. 60, 2017, p. 1-15. [Publisher]