SciELO - Scientific Electronic Library Online

 
 issue103Molecular dynamics simulations of the mechanical behavior of nanostructured and amorphous Al80Ti15Ni5 alloyHeavy metal removal by biopolymers-based formulations with native potato starch/nopal mucilage author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Revista Facultad de Ingeniería Universidad de Antioquia

Print version ISSN 0120-6230On-line version ISSN 2422-2844

Rev.fac.ing.univ. Antioquia  no.103 Medellín Apr./June 2022  Epub Feb 17, 2022

https://doi.org/10.17533/udea.redin.20201010 

Artículo original

Control scheme selection and optimal tuning in industrial process control using factorial experiment design

Selección de estrategia y sintonización óptima de control industrial usando un diseño de experimento factorial

David Roberto Acosta-Villamil1 

José Fernando Noguera-Polania1 

Arnaldo Verdeza-Villalobos2 

Blanca Luz Foliaco-Romero1 

Adriana Fernanda Rincón-Montenegro1 

Marco Enrique Sanjuan-Mejía1 

1Departamento de Ingeniería Mecánica, Universidad del Norte. Km. 5 vía Puerto Colombia, Bloque k, Piso 8. C.P. 0000. Barranquilla, Colombia.

2Departamento de Ingeniería Industrial, Universidad Simón Bolívar. Carrera 59 # 59-65. C. P. 1086. Barranquilla, Colombia.


ABSTRACT

In this study, a novel experimental approach for the optimal selection of an actuator-based control strategy is presented. The proposed approach is a two-stage method: first, a two-level factorial experiment design with n factors (2n) was applied to compare different control schemes. Schemes comparison was carried out in terms of energy consumption and closed-loop performance. For the best relative scheme, a Central Composite Face-centered (CCF) design was completed obtaining the controller parameters that optimize the performance in terms of the Integral Absolute Error (IAE) while operating in a region of low energy consumption. The proposed approach was experimentally tested using real data obtained from a laboratory prototype plant. Some experimental tests illustrating the suitability of our method are shown at the end of this article.

Keywords: Optimal control tuning; process control; design of experiments; factorial design

RESUMEN

En este estudio se presenta un nuevo enfoque experimental para la selección óptima de una estrategia de control basada en el actuador. El enfoque propuesto es un método de dos etapas: primero se aplica un diseño de experimento factorial de dos niveles con n factores (2n) para comparar diferentes esquemas de control. La comparación de esquemas se lleva a cabo en términos de consumo de energía y rendimiento de circuito cerrado. Para el mejor esquema relativo, se completa un Diseño Central Compuesto Centrado en las Caras (CCF, por sus siglas en inglés) obteniendo parámetros de controlador que optimizan el rendimiento, en términos del Error absoluto integral (IAE, por sus siglas en inglés), mientras operan en una región de bajo consumo de energía. El enfoque propuesto se probó experimentalmente utilizando datos reales obtenidos de una planta prototipo de laboratorio. Algunas pruebas experimentales que ilustran la idoneidad de nuestro método se muestran al final de este artículo.

Palabras claves: Sintonía óptima de control; control de procesos; diseño de experimentos; diseño factorial

1. Introduction

Controllers are a vital part of industrial processes, and plants around the world use predominantly proportional-integral-derivative (PID) controllers to keep important process variables at desired setpoints. It is estimated that over 90% of control loops employ PID control, quite often with the derivative gain set to zero, i.e., PI control[1]. Although academic research in control theory has developed tools for tuning controllers and measuring performance, only about 32% of control loops have acceptable performance, and 36% of the processes are operating in open-loop because of problems with the controller[2].

Tuning equations are a tool to adjust controller parameters, which, due to their versatility and generalization capabilities, provide a way to optimize the response of processes finding specific parameters according to performance requirements. For example, to preserve the integrity of actuators, some parameters can be found to minimize the variation of the controller output[3]. Tuning equations for a PID controller based on Internal Model Control (IMC), where the IAE is minimized for a second-order inverse response system with a time delay have been developed[4]. Also, tuning equations for PID controllers to be implemented in photobioreactors, by optimizing the IAE and variance for both controller output and process variables have been investigated[5]. For instance, tuning equations to optimize the process performance using PID controllers for heat exchangers and thermal processes have been used[6].

New tuning approaches for PID controllers include the use of genetic algorithms to obtain controller parameters[7], and metaheuristic algorithms like swarm optimization, which use as a performance index the IAE[8]. In other swarm intelligence algorithms, the estimated mathematical models of systems or processes can be used to modify the parameters of PID controllers. These adaptive control strategies can be implemented in commercial programmable logic controllers (PLC)[9, 10].

No doubt, control loops are the most critical components in automation systems. Energy and raw material consumption, operation safety, product quality, and thus, the financial profitability are directly or indirectly linked to the performance of process control systems[11]. Hence, a vital interest for any company must be to reach, restore, and maintain top performance of control loops. However, most of the process control systems currently used are complex, usually comprising many hierarchical levels; thus, it is almost impossible for plant personnel to maintain them on top performance. Additionally, determining a feasible operating region within acceptable variable ranges and optimal performance is a difficult task in process engineering; therefore, several works investigate the correlation between performance criteria and controller parameters or structure. For instance, using the knowledge from the feasible operating region to design the controller structure for a flash distillation process, where the improved design outperforms the classical configuration[12].

The use of the design of experiments (DoE) methods, such as the response surface methodology, as a design tool, is widely accepted to be a powerful approach to simultaneously optimize complex processes in terms of various response variables, for instance, system performance and energy consumption[13, 14]. The response surface methodology can also be used together with simulation studies to investigate the effects that factors have upon response variables[15].

In this work, a novel experimental-based approach to reach an optimal performance is proposed. The core of our method is DoE. The flow control process from the TecQuipment's Process Trainer is used in a two-stage method: first, a two-level factorial experiment design with n factors (2n) was applied to compare different control schemes. The scheme comparison is carried out in terms of energy consumption and closed-loop performance. For the best relative scheme, a central composite face-centered (CCF) design was completed obtaining controller parameters which optimize performance, in terms of the integral absolute error (IAE) index, while operating in a region of lowest energy consumption. Notice that it is not necessary to develop a complex first principle model of the process to be controlled, enabling the use of our approach in industrial process plants. This work’s main motivation was to propose an experimental approach to find a controller parameter within an acceptable search space, which yielded an optimal system response regarding variables related to the system outputs. In the case presented, using the response surface methodology and defining the limits of the search space with a traditional tuning method, optimal controller parameters were found with respect to the response variables IAE and energy consumption.

2. Methodology

In this section, a case study is presented, where a control scheme is selected, and the optimal tuning of a flow control loop is performed, for which experimental data collected from plant operating conditions.

2.1 System description

The prototype plant CE117 Process Trainer from TecQuipment (see Figure 1] was utilized to validate the proposed approach.

Figure 1 TecQuipment's process trainer 

The CE117 Process Trainer is a fully integrated laboratory-scale process control plant equipped with valves, pumps, power supplies, and ancillaries to allow the implementation of flow, level, temperature, and pressure control strategies. More precisely, in the flow control loop, the prototype plant consists of a reservoir tank from which water is lifted by a hydraulic pump (P1) into a process tank; in this line, a servo-controlled proportional valve (V1) and a flow transmitter (FT) are situated. Additionally, after the pump, there is a process loop bypass valve (V2), which will be used within the framework of this study to simulate disturbances due to leaks. The piping and instrumentation diagram (P&ID) describing the process under study is presented in Figure 2. The drivers of the components of the plant use analog signals with a range from 0.0 to 10.0 Volts. Further technical data are presented in Table 1.

Figure 2 Piping and instrumentation diagram of the process 

2.2 Control Schemes

Two actuator-based control strategies were considered for the flow control loop of the process: in the first one, the actuator of the final control element is the pump (P1). In contrast, in the second one, both the pump (P1) and the proportional valve (V1) are considered as actuators. In both cases a proportional-integral (PI) controller was considered.

Control schema 1: in this scheme, the hydraulic pump is driven by the PI controller while the proportional valve is held open with its driver at a constant value of 8.0 V; in Figure 3, a piping and instrumentation diagram of this scheme is presented.

Figure 3 P&ID of the pump-based control strategy 

Control schema 2: in this case, both the pump and the proportional valve drivers are manipulated by PI controllers, with flow being the controlled variable in both cases. The piping and instrumentation diagram of this strategy is shown in Figure 4.

Figure 4 P&ID of the pump and valve-based control strategy 

2.3 Benchmark Process

To evaluate the control strategies regarding both energy consumption and closed-loop performance, a setpoint signal consisting of a sequence of three steps with amplitudes of 5.0 L/min, 2.0 L/min, and 3.0 L/min, respectively, with a total duration of 150 seconds was established as the benchmark. Figure 5 shows the setpoint signal. To assess energy consumption, a Belkin WeMo Insight Switch is used to obtain two instantaneous power measurements during each cycle.

Figure 5 Setpoint sequence for flow control loops 

The power measurements are averaged and multiplied by the process duration to calculate energy consumption (Ec) in [W/h]; this is one of the process variables analyzed in the DoE. Considering that the energy consumption of the pump is significantly higher than that of the valve, the controller for the valve is tuned and kept invariant throughout this work. In contrast, the controller for the pump is investigated using the proposed methodology.

2.4 Tuning for the Valve PI-Controller

A process step test was carried out as follows: 1) with the system in open loop, a step change in the manipulated variable was applied. It is important to remark that the magnitude of the change should be large enough to generate a measurable change in the transmitter signal, in this case, a change in the valve driver is made from 2.0 V to 9.0 V, while keeping the pump driver at 8.0 V. 2). The response of the transmitter output signal is recorded, making sure that the resulting plot of the controlled variable versus time covers the entire test period from the introduction of the step test until the system reaches a new steady state. Using the output response of the plant, i.e., liquid flow as a time series, the dynamic behavior of the system can be identified. In this case, the maximum rate of change can be observed as soon as the input step was applied, with a small time-delay. This response can be modeled as a first-order system plus dead time; for this reason, a linear model was identified with a First-Order-Plus-Dead-Time (FOPDT) structure using the MATLAB System Identification Toolbox. The model is described by the transfer function in (1), and the model and real plant response are plotted in Figure 6.

Figure 6 FOPDT model and actual responses for the proportional valve 

The model was then used to find the parameters for the PI Controller of the valve. With the help of the interactive PID-Tuner interface in Matlab, the gains in the controller are tuned to obtain closed-loop stability, adequate performance, and robustness[16]. In a parallel configuration, the gains for the PID are set to kc = 3.0 and ki = 1.9, which results in a not too aggressive control of the flow via the valve. For a unit step, the closed-loop model had a rise time under 3 seconds, a settling time under 10 seconds with an overshoot of less than 3%.

A FOPDT model is obtained for the pump (2), with a procedure like the one used to obtain the valve model, varying the pump driver from 2.0 V to 7.0 V.

For the pump, tuning is accomplished using the rules of Internal Model Control[17]. Before presenting this procedure, the DoE methodology must be introduced, considering that the value of the tuning parameter is selected according to the experimental conditions described in the DoE methodology.

2.5 DoE Methodology

A full factorial (24, n = 1) experiment with two levels for each factor is conducted initially to screen for significant factors. The experimental conditions are shown in Table 1. Based on the results, a Central Composite Face-centered (CCF) is carried out and examined in section 3, to determine optimum areas of operation from which a new operating point is proposed and tested in the real plant.

Table 1 Conversion details for components of the system 

The coded (1) and (-1) PI-controller gains k c and k i , in a parallel configuration, are obtained by varying the parameter λ from the tuning rules by internal model control[17-19]. Both values for the parameter λ (0.1 and 2.0) yield acceptable process closed-loop performance. The equations for a PI-Controller are given by Equation (4), and ki = kcτ in the PI parallel configuration. The values of kp = 0.760, τ = 1.560 and t0 = 0.200 refer to the FOPDT model in (3) and (4).

Thus, the natural values for Factors B and C, are found and shown in Table 2 for the complete set of Factors.

3. Experiments

The factorial experimental design and the results for both process variables are presented in Table 3, in the randomized order of the runs. Data analysis was performed using the software Statgraphics Centurion XVII.

Table 2 Experimental conditions 

Table 3 Screening experiment 

3.1 Screening energy consumption

Before discussing the results of the analysis of variance (ANOVA), the assumptions of normality, homogeneity of variance, and independence of the residuals were tested for both process variables (Appendix).

Consequently, after examining the normal probability plot of the ANOVA, interactions with a low level of significance were removed from the analysis of variance sequentially. Factor B, its interactions, and the interaction between Factor C and D are excluded from the ANOVA, thus remaining with the factors shown in the ANOVA results in Table 4.

Table 4 Analysis of variance for E c  

From the ANOVA results, the factors with higher effects are A and D, which can also be seen in Figure 7.

Figure 7 Main effects plot for Ec 

It could be inferred apriori that Factor D plays a significant role in energy consumption. When the bypass valve is opened, to maintain the same flow setpoint, a higher energy input is necessary. For Factor A, the performed analysis suggests that the pump and valve control scheme has lower energy consumption. It furthermore shows that the integral gain on the controller affects both controller schemes differently in terms of energy consumption.

3.2 Screening closed-loop performance

A similar procedure is conducted with the results of the ANOVA using the normal probability plot for the IAE, and sequentially excluding interactions and factors with low levels of significance. It is concluded that based on the conducted experiment, only Factor B has a significant effect on the closed-loop performance in terms of the IAE, as can be seen in the Pareto in Figure 8.

Figure 8 Pareto chart of the effects for the IAE 

3.3 Central Composite Face-centered (CCF) DoE

According to Montgomery[14], the Response surface methodology (RSM) is a set of powerful statistical techniques that can be used to model problems where one or many independent variables are influenced by several other variables. The analysis of problems using this method can optimize these response variables, where the analytical nature between the response variables and the independent variables is unknown. Normally, a low order polynomial in the region of interest defined by the independent variables is used to investigate the relationship. The CCF is a class of DoE and is the most popular design structure to identify second-order models, where the curvature of the independent variable is of interest.

Based on the conclusion drawn from the previous screening experiment, the relation between factors B and C for the pump and valve control scheme is further explored, i.e., Factor A equals 1 and fixing the bypass valve to closed; thus, Factor D equals -1. A CCF design is completed with axial points located at the centers of the faces of the squared region, defined by the experimental conditions used in the screening experiment, with two central points. The new experimental conditions associated with Factors C and D equal to zero are found by interpolating the parameter λ from the tuning rules by internal model control (IMC). For the region between the coded values -1 and 1 and natural values 0.1 and 2.0, respectively, the center yields for the coded value 0, the natural value λ = 0.95. Following the IMC tuning rules with the new parameter λ, the new levels for the factors are found to be kc = 1.210 and ki = 0.774. The CCF design with the values for the output variables is presented in Table 5.

Table 5 Central composite face-centered 

Using the ANOVA results to analyze the response variable energy consumption, first interactions among factors are sequentially eliminated, and the ultimately factor B is removed from the analysis of variance. It is clear that the information extracted from the ANOVA, presented in Table 6, indicates that the Integral Constant of the PI-Controller significantly affects the consumption of energy. That can be explained with the impact of changing the slope with which the response moves towards the reference, when this tracking is slow for small values of the parameter K i , in this case towards the point Factor C equal to 1, consumption increases as seen in Figure 9.

Table 6 CCF-analysis of variance for E c  

Adjusted R-Squared = 68.33 %

Figure 9 CCF-Plot of main effects for Ec 

The results of the ANOVA for the IAE response, after removing the interactions with a low level of significance, are shown in Table 7. The response surface plot resulting from the regression model of the analysis is shown in Figure 10.

Figure 10 CCF-Response Surface for the estimated IAE 

Table 7 CCF-analysis of variance for IAE 

Adjusted R-Squared = 77.10 %

3.4 Analysis of the results

The operating conditions that minimize the estimated IAE are shown in Table 8, with their coded values, and in Table 9, their mapped natural values are listed.

Table 8 IAE response optimization for the estimated IAE 

Optimal Value =15.33

Table 9 Natural values of the optimal operating condition 

The found optimal operation conditions are then tested in the real process, obtaining the results for Energy consumption and IAE shown in Table 10.

Table 10 Actual performance of the optimal point 

The closed-loop response of the real system using the optimal operating conditions for the controller gains is shown in Figure 11, and the authors qualify it as an excellent performance. The control effort for the pump and the proportional valve using the controller gains for the optimal operating condition are shown in Figure 12.

Figure 11 The closed-loop response of the system at the optimal operating point 

The closed-loop performance of the experimental condition corresponding to the CCF run number 8 is better in terms of the IAE; However, its energy consumption is higher since consumption is highly correlated to the control effort of the pump. In Figure 13, the control effort for both experimental conditions, CCF-Run 8, and the optimal operating point are presented. Notice that the pump control effort for the experimental condition deemed optimal by the regression model and the response surface methodology is considerably lower than the control effort for the CCF-Run 8.

Figure 12 Control effort for the optimal operating point 

Figure 13 Pump control effort for CCF-Run8 (1, -1, 0, -1) and Optimal point 

4. Conclusions

The response surface methodology (RSM) is a mighty collection of techniques that, if used appropriately, become tools to optimize complex processes where several response variables are affected by numerous identified factors. The statistical information abstracted from the data collection representing a stochastic and even dynamic process allows the methodology’s user to explore the true response function. Regression models manage to extract the natural error of the system and perform system optimization without analytically determining the complicated relationship between response variables and factors of influence.

An objective parameter to compare the energy consumption of two different control architectures was proposed. Following a DoE approach for the assessment, two-level factorial experiments were used for screening Factors, and a comparison of controller schemes was completed in a real process.

Within this work, the RSM was used successfully to minimize two variables of interest, Energy consumption, and closed-loop performance, of a specific process, by exploring the gains in a controller as Factors. The found optimal operating condition is tested in the real system. The regression model predicted at the optimal operating conditions an =15.33, and the real system responded with an IAE=15.92. The difference between them is only 3.8%, with energy consumption considerably lower than most of the operating conditions explored in the data collection process. This is a clear indication of the benefits of the proposed approach.

The regression model from the proposed methodology found controller parameters, which were used as a new experimental condition. Subsequently, using these parameters for the controller, the response showed a lower control effort than all the other experimental conditions. This was because the control effort for the pump in this new experimental condition required a lower voltage causing lower energy expenditure. The proportional gain was set close to one of the boundaries causing a fast response to the step changes followed by a smooth error elimination phase whose rate depends on the integral time.

6. Acknowledgements

This work was supported by Universidad del Norte, Universidad Simón Bolívar and Colciencias through Convocatoria No. 727 de 2015.

Referencias

[1] C. Knospe, “PID control,” IEEE Control Syst. Mag., vol. 26, no. 1, February 2006. [Online]. Available: https://doi.org/10.1109/MCS.2006.1580151Links ]

[2] M. L. Brisk, “Process control: Potential benefits and wasted opportunities,” Aust. J. Electr. Electron. Eng., vol. 2, no. 1, 2005. [Online]. Available: https://doi.org/10.1080/1448837X.2005.11464113Links ]

[3] P. Klán and R. Gorez, “PI controller design for actuator preservation,” IFAC Proc. Vol., vol. 41, no. 2, 2008. [Online]. Available: https://doi.org/10.3182/20080706-5-KR-1001.00981Links ]

[4] F. Castrillón and D. Castellanos, “New tuning rules for PID control lers based on IMC with minimum IAE for inverse response processes,” DYNA, vol. 82, no. 194, November 2015. [Online]. Available: http://dx.doi.org/10.15446/dyna.v82n194.46744Links ]

[5] A. Verdeza and L. Di Mare and M. Sanjuán and A. Bula, “Diseño de ecuaciones de sintonía para controladores PID (Proporcional-Integral-Derivativo) implementados en fotobiorreactores,” Inf. tecnológica, vol. 27, no. 4, 2016. [Online]. Available: http://dx.doi.org/10.4067/S0718-07642016000400013Links ]

[6] J. Duarte and W. Orozco, “Optimización de sintonización de controladores PID bajo el criterio IAE aplicados a procesos térmicos,” Revista Inge@UAN, vol. 5, no. 10, pp. 35-45, 2015. [ Links ]

[7] J. Shen, “New tuning method for PID controller,” ISA Trans., vol. 41, no. 4, October 2002. [Online]. Available: https://doi.org/10.1016/S0019-0578(07)60103-7Links ]

[8] M. T. Özdemir and D. Öztürk and Ī. Eke and V. Çelik and K. Y. Lee, “Tuning of optimal classical and fractional order PID parameters for automatic generation control based on the bacterial swarm optimization,” IFAC-PapersOnLine, vol. 48, no. 30, 2015. [Online]. Available: https://doi.org/10.1016/j.ifacol.2015.12.429Links ]

[9] O. García, D. Acosta, and C. Diaz, “GPU-Implementation of a sequential monte carlo technique for the localization of an ackerman robot,” in International Conference on Applied Informatics ICAI 2018: Applied Informatics, Bogotá, Colombia, 2018, pp. 309-320. [ Links ]

[10] J. A. López, A. Duque, and A. F. Navas, “Sintonización de un controlador PID en un PLC haciendo uso de inteligencia de enjambres/auto-tuning of a PID control ler implemented in a PLC using swarm intelligence,” Prospectiva, vol. 15, no. 1, January 2017. [Online]. Available: https://doi.org/10.15665/rp.v15i1.679Links ]

[11] M. Bauer, A. Horch, L. Xie, M. Jelali, and N. Thornhill, “The current state of control loop performance monitoring - A survey of application in industry,” J. Process Control, vol. 38, February 2016. [Online]. Available: https://doi.org/10.1016/j.jprocont.2015.11.002Links ]

[12] K. L. Morales and H. D. Álvarez, “Determination and use of feasible operation region in flash distillation control,” Revista Facultad de Ingeniería Universidad de Antioquia, vol. 95, April 2020. [Online]. Available: http://dx.doi.org/10.17533/udea.redin.20190738Links ]

[13] A. S. Comas, A. T. Palacio, S. T. Mendoza, and D. N. Rodado, “Aplicación del diseño de experimentos taguchi para la identificación de factores de influencias en tiempos de impresión 3d con modelado por deposición fundida,” Int. J. Manag. Sci. Oper. Res., vol. 1, no. 1, pp. 43-48, Jan. 2016. [ Links ]

[14] D. C. Montgomery, Design and analysis of experiments, 9th ed. John wiley & sons, 2017. [ Links ]

[15] J. C. Salazar and A. B. Zapata, “Análisis y diseño de experimentos aplicados a estudios de simulación,” DYNA, vol. 76, no. 159, pp. 249- 257, Sep. 2009. [ Links ]

[16] K. J. Åström and T. Hägglund, Advanced PID control . ISA-The Instrumentation, Systems, and Automation Society, 2006. [ Links ]

[17] C. A. Smith and A. B. Corripio, Principles and Practice of Automatic Process Control, 3rd ed. John wiley & sons, 2006. [ Links ]

[18] J. Duarte and et al, “Auto-ignition control in spark-ignition engines using internal model control structure,” J. Energy Resour. Technol., vol. 139, no. 2, March 2017. [Online]. Available: https://doi.org/10. 1115/1.4034026Links ]

[19] H. Zeng and K. Xiao, “A new design of multivariable decoupling internal model controller,” J. Theor. Appl. Inf. Technol., vol. 48, no. 1, pp. 417-422, Feb. 2013. [ Links ]

Received: February 25, 2020; Accepted: October 29, 2020

* Corresponding author: David Roberto Acosta Villamil, e-mail: villamilr@uninorte.edu.co

5. Declaration of competing interest

We declare that we have no significant competing interests, including financial or non-financial, professional, or personal interests interfering with the full and objective presentation of the work described in this manuscript.

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License