Next Article in Journal
Development of Tea Tree Oil Based Nanoemulgel Loaded with Azithromycin for Enhancing the Antibacterial Activity
Next Article in Special Issue
Thermal-Hydraulic Characteristics of Carbon Dioxide in Printed Circuit Heat Exchangers with Staggered Airfoil Fins
Previous Article in Journal
Effects of Bio-Coal Briquette for Residential Combustion on Brown Carbon Emission Reduction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Control of Technological Processes

by
Anatoliy M. Tsirlin
1,* and
Alexander I. Balunov
1,2
1
System Analysis Research Center, Ailamazyan Program Systems Institute of RAS, 152021 Pereslavl-Zalessky, Russia
2
Department of Cybernetics, Yaroslavl State Technical University, 150023 Yaroslavl, Russia
*
Author to whom correspondence should be addressed.
Processes 2023, 11(6), 1835; https://doi.org/10.3390/pr11061835
Submission received: 16 May 2023 / Revised: 9 June 2023 / Accepted: 12 June 2023 / Published: 16 June 2023
(This article belongs to the Special Issue Advances in Process Optimization and Control)

Abstract

:
The paper formulates conditions under which the roots closest to the imaginary axis (critical roots) of the characteristic equation of a linearized system are real for the maximum possible degree of stability of the closed-loop control system of a technological process with pure delay. For the parameters of the controllers corresponding to the maximum degree of stability, these roots are multiples. Their multiplicity order is one more than the number of coefficients in the transfer function of the controller. It is demonstrated that for a typical technological control object, these conditions are satisfied for all “serial” control laws. This allowed for obtaining analytical expressions for optimal settings and limiting degrees of stability as functions of object parameters for typical dynamic characteristics of technological processes. The paper considers the problem of robust stability for control systems with an object containing pure delay. It has been proven that in the maximum stability problem, the operations of maximizing over controller parameters and minimizing over the set of possible object parameters can be interchanged. Therefore, selecting robust settings amounts to determining the minimum of the maximum stability over the set of possible object parameter values. Controllers with such settings are suitable, without modification, for a whole class of technological processes.

1. Introduction

The stability of dynamic systems and the selection of controller parameters are the most studied topics in control theory. A considerable number of studies were devoted to this theme in the 1950s and 1960s. A.M. Letov [1] paid significant attention to this area, using methods based on the Lyapunov function defined on the state space of the dynamic system. These methods allow for the nonlinearity of equations to be taken into account, but they are only effective for finite-dimensional systems.
Frequency methods are applicable only to linear systems, but they allow for considering systems with delay, which is essential for technological processes as control objects. The dynamics of these objects are mainly determined by heat and mass transfer processes. They are stable, characterized by distributed parameters, and their impulse transfer functions do not change sign and tend to zero or to some constant if the object contains an integrating element. The transfer functions of such objects contain pure delay, and their characteristic equation does not have complex roots. The magnitude and phase of their frequency response decrease monotonically with increasing frequency, with the magnitude tending to zero and the phase to minus infinity (the “monotonicity” effect). For brevity, such objects will be referred to as technological.
Tsypkin introduced the concept of the stability degree of a linear system as the distance from the imaginary axis to the nearest (critical) root of the system’s characteristic equation. If each of the roots is assigned an index ν , then the stability degree is the absolute value of the maximum over ν of the real part of the root (in a stable system, the real parts of all roots are negative). He used frequency methods in his formulated problem of maximizing the stability of linear systems [2,3]. This topic was also the subject of his dissertation.
The problem of maximizing stability was significantly developed by Shubladze and his colleagues [4,5].
Another feature of controllable technological processes is the change in their dynamic properties when the raw material composition and flow intensity are altered. Therefore, researchers have paid significant attention to the problem of object control with varying parameters and the synthesis of robust controllers that, with fixed settings, can control a whole class of objects or one object within a wide range of parameter changes. For finite-dimensional problems, methods for synthesizing such systems with an estimation of the permissible range of object parameters was developed in the work of Polyak and Shcherbakov (see [6]). However, these methods are not applicable to objects with delays, which includes the majority of technological processes.
The selection of controller parameters based on the condition of maximum stability indirectly ensures the robustness of the system. If the dependence of the maximum stability on the parameters of the object is obtained, and the region in which these parameters can vary is known, then the most “unfavorable” combination of parameters from the region of their possible values is selected. At such parameters, the maximum stability is minimal. This is a guaranteed degree of stability of the system throughout the range of changes in its characteristics.
It is significantly easier to solve these problems if an analytical dependence of the maximum stability on the object’s parameters is found. As shown below, such a dependence can be easily obtained when the roots closest to the imaginary axis (critical roots) are real. In [7], the maximum stability in this case is called aperiodic, and in the case when the critical roots are complex—oscillatory. It is not known in advance which of these cases is true. The controller parameters are chosen based on the condition of maximum aperiodic stability, and then it is checked whether the synthesized system has complex roots that are closer to the imaginary axis than real roots. The paper shows under what conditions the answer to this question is negative, and therefore the relation between the object’s parameters and the controller settings, found by the condition of maximum aperiodic stability, can be used to calculate the system.
The problem of calculating the limit aperiodic stability is considered by plotting a hodograph of the extended frequency response of the open-loop system with maximum aperiodic stability.
An important feature of technological processes is that, despite their enormous variety, they are similar to each other in their dynamic characteristics, and in most cases, their dynamics in the Fourier or Laplace transform domain can be approximated by an aperiodic or integrating element with pure delay. This makes it possible to use a range of typical “serial” control laws consisting of proportional (P), integral (I), proportional-integral (PI), and proportional-integral-differential (PID) controllers in control systems.
The paper presents formulas for selecting parameters of typical controllers based on the condition of maximum system stability for objects with delay, single-loop, and two-loop systems. A methodology for selecting robust settings is proposed, and a real-time optimization system is considered to support the optimality conditions of the technological system.

2. Calculation of Optimal Controller Settings

Let us consider a linear single-loop automatic control system (Figure 1).
Here, W 0 ( p ) и W r ( S , p ) are the transfer functions of the object and the controller, respectively, where S is the vector of controller parameters. The transfer function of the closed-loop system is
W s ( p ) = W 0 ( p ) 1 + W 0 ( p ) W r ( S , p ) ;
assume that roots of the characteristic equation do not coincide with the zeros of the object’s transfer function.
The controller usually has a standard PID structure
W r ( S , p ) = S 2 p + S 1 + S 0 p .
If S 2 = 0 , then the controller is proportional-integral (PI), if S 2 and S 1 are zero, it is integral (I), and if S 2 and S 0 are zero, it is proportional (P).
An integral component ensures an astatic character of the transition processes. In this case, the characteristic equation has no zero roots. If there is no integral component ( S 0 = 0 ), then the characteristic equation has a root at the origin, and the system, when subjected to step-like input, does not return to the equilibrium state and exhibits a “residual static error”. In the latter case, when calculating the maximum stability degree, all roots of the characteristic equation but zero are taken into account.
Based on the made assumptions, the characteristic equation of System (1) can be written equivalently as
p W 0 ( p ) + S 2 p 2 + S 1 p + S 0 = 0 , p ν = ρ ν ± i ω ν , ν = 1 , 2 , ,
Choosing the value of vector S aims to ensure that the roots of the characteristic equation are located to the left of the imaginary axis while the distance from the nearest roots η * = | ρ * | is maximized. This requirement is formalized as a minimax.
Problem A:
ρ * ( S ) = max ν ρ ν ( S ) min S .
The dependency ρ ν ( S ) is determined by solving Equation (3), where the real parts of all roots are known to be negative. If the maximization problem with respect to ν has no solution, then instead of a maximum, the exact upper bound ρ ν shall be found.
As only the coefficients at low powers of p in the characteristic equation depend on the controller settings, any changes in these settings do not affect the sum of the real parts of its roots
ρ S = ν ρ ν ,
which, for the nth power equations, according to Vieta’s formulas, is equal to the coefficient at p n 1 with the opposite sign. Thus, decreasing the real parts of the roots closest to the imaginary axis leads to an increase in ρ ν for the remaining roots. At the same time, the vector S primarily affects the roots located closer to the imaginary axis in the p [1] plane. Therefore, the maximum stability
ρ * = min S ρ * ( S )
corresponds to the case where, for several real roots or several pairs of complex roots, the values of the real parts are the same.
The number of such “critical” roots, if task (4) has a solution, is one more than the number of adjustable controller coefficients. This type of solution structure is typical for minimax problems and, in the case of the problem of uniform approximation, is known as the “Chebyshev alternance” principle.
For the case where all critical roots are real, problem (4) is relatively easy to solve, and in some cases, it can be solved analytically. The limit of aperiodic stability can be found from the condition of multiplicity of the ( m + 1 ) th critical root
d m ( d ρ ) m ρ W 0 ( ρ ) ρ * = 0 ,
where m is the number of required controller settings (in Equation (3), m = 3 ).
For S 0 = 0 , the left-hand side of Equation (3) can be reduced by p. At the same time, a root with ρ 0 = 0 , and changing settings S 2 , S 1 does not change ρ 0 , which leads to a steady-state error. This leads to abrupt changes of ρ * ( S ) at S 0 = 0 .
If Equation (6) can be solved when S 0 0 , then the corresponding optimal settings are found as
S 2 * = 0.5 d 2 d ρ 2 ρ W 0 ( ρ ) ρ * ,
S 1 * = 2 S 2 * ρ * + d d ρ ρ W 0 ( ρ ) ρ * ,
S 0 * = ρ * S 1 * + S 2 * ρ * + 1 W 0 ( ρ * ) .
For S 0 = 0 , settings S 1 and S 2 , the PD-controller is chosen in such a way that the roots of the equation
1 W 0 ( p ) + S 2 p + S 1 = 0
satisfy requirement (4). The limit of aperiodic stability (taking into account all non-zero roots) and the corresponding settings are found from the conditions:
d 2 d ρ 2 1 W 0 ( ρ ) ρ * = 0 ,
S 2 * = d d ρ 1 W 0 ( ρ ) ρ * ,
S 1 * = ρ * S 2 * + 1 W 0 ( ρ * ) .
The settings found in this way, corresponding to the limit of aperiodic stability of the system, are a solution to problem (4) only if all complex roots of the characteristic equation of the system for the settings selected in this way are to the left of the line ρ a * ± i ω . This leads to
Problem B: Under which conditions does the limit of aperiodic stability ρ a * coincide with the maximum stability of the system (i.e., critical roots are real)?

2.1. Conditions of Optimality of the Limit of Aperiodic Stability

Suppose (8) and (9) give ρ a * , S 0 * , S 1 * , and S 2 * . They correspond to the frequency response of the open-loop system W c ( i w ) = W 0 ( i ω ) W r ( S * , i ω ) and its extended frequency response
W c ( i ω , ρ a * ) = W 0 ( i ω ρ a * ) W r ( S * , i ω ρ a * ) .
Figure 2 shows the location of the roots of the characteristic equation of the system in cases where ρ a * corresponds to the maximum stability (a) and when it does not correspond (critical roots are complex) (b).
To determine which case is true, let us construct a Nyquist plot of an extended frequency response of an open-loop system W c ( i ω , ρ a * ) . It is clear that due to condition (3) at ω = 0 , this plot will pass through the point ( 1 , i 0 ) (Figure 3). If all the characteristic equation roots of the system lie to the left of the line ρ a * ± i ω , then when ω changes from zero to infinity, the Nyquist plot W c ( i ω , ρ a * ) will not encircle the point ( 1 , i 0 ) (Figure 3a).
Let us denote the frequencies with the phase of frequency W c ( i ω , ρ a * ) equal to π as critical frequencies ω 0 ν . In this case, ω 00 = 0 ; ω 01 corresponds to a rotation of the plot by 2 π , ω 02 by 4 π , and so on.
At the frequency ω 00 | W c ( i 0 , ρ a * ) | = 1 .
The properties of conformal mapping determine the following.
Statement:
The limit of aperiodic stability is the maximum possible stability of the system only if ν > 0
| W c ( i ω 0 ν , ρ a * ) | 1 .
Particularly, for optimality ρ a * , it is sufficient for the modulus of | W c ( i ω , ρ a * ) | to decrease monotonically while ω increases. Conversely, the limit of aperiodic stability is not optimal if | W c ( i ω , ρ a * ) | increases monotonically while the frequency increases. It is typical for technological processes to exhibit a monotonic decrease of the modulus of the extended frequency response of an open-loop system with frequency.

2.2. Typical Technological Process Control System

A typical industrial subject of control is defined as [8,9]. The dynamics of most processes can be approximately described by a transfer function
W 0 ( p ) = k e p τ T p + 1 ,
as it is close to the dynamics of the subject in question.
To calculate its three coefficients from the curve obtained after applying a step of Δ at the input, a tangent is drawn to this curve at the inflection point. The distance from the point of intersection of this tangent with the x-axis to the origin is τ . The ratio of the steady-state deviation to Δ is k. The tangent of the slope of the tangent is 1 / T .
Let us note that from the transfer function (16) after substituting k = T Θ and subsequent limit transition at T we obtain an integrating object with delay
W 0 ( p ) = e p τ Θ p .
By directing T to zero in (16), we obtain a pure delay object
W 0 ( p ) = k e p τ .
For each of the noted typical objects, analytical expressions for the achievable aperiodic stability ρ a * and corresponding controller settings, expressed through the object parameters, were obtained using (7)–(9) and (11)–(13) (see also [2,5]).
Let us determine if the limit of aperiodic stability reaches the maximum possible stability for systems with typical objects (16)–(18) and I, PI, and PID-controllers included in the negative feedback.

2.2.1. The Pure Delay Object

I-controller:
ρ a * = 1 τ .
The only value of the setting parameter
S 0 = 1 k τ e 1
corresponds to ρ a * . The modulus of the extended frequency response of the open-loop system
| W c ( i ω , ρ a * ) | = 1 ω 2 τ 2 + 1
decreases monotonically from 1 at ω = 0 to zero at ω →∞. In this case, the limit of aperiodic stability is the maximum possible.
PI-controller:
ρ a * = 2 τ .
The corresponding setting parameters
S 0 = 4 k τ e 2 , S 1 = 1 k e 2 .
The modulus of the extended frequency response of the open-loop system is 1. Thus, the critical roots are the real roots and any number of complex roots at critical non-zero frequencies. The limit of stability is ρ a * .
PID-controller:
ρ a * = 3 τ .
The corresponding setting parameters
S 0 = 27 2 k τ e 3 , S 1 = 5 k e 3 , S 2 = τ 2 k e 3 .
The modulus of the extended frequency response of the open-loop system
| W c ( i ω , ρ a * ) | = τ 6 τ 2 ω 2 2 τ 2 + 4 ω 2 ω 2 τ 2 + 9 .
It is 1 at ω = 0 and increases monotonically with the increase of ω . Thus, for such a system, the limit of aperiodic stability is not the maximum possible.

2.2.2. An Integrating Object with Delay

PI-controller:
ρ a * = 2 2 τ .
Optimal controller parameters
S 0 = Θ τ 2 2 2 1 3 e 2 2 , S 1 = Θ τ 2 2 1 e 2 2 .
Modulus W c ( i ω , ρ a * )
| W c ( i ω , ρ a * ) | 0.1177 + 0.6864 ω 2 τ 2 0.1177 + 0.6864 ω 2 τ 2 + ω 4 τ 4 .
It monotonically decreases with increasing frequency, therefore, the limit of aperiodic stability is the maximum possible. It is important to note that this modulus does not depend on the object’s parameter Θ .
PID-controller:
ρ a * = 3 3 τ .
Corresponding controller parameters
S 0 = Θ τ 2 3 3 3 1 3 3 2 e 3 3 , S 1 = Θ τ 3 3 2 3 e 3 3 ,
S 2 = Θ 2 3 3 1 3 3 2 2 e 3 3 .
The modulus of the extended frequency response of the open-loop system
| W c ( i ω , ρ a * ) | 2.5844 + 0.8615 ω 2 τ 2 2.5844 + 3.2151 ω 2 τ 2 + ω 4 τ 4 .
It is easy to see that the limit of aperiodic stability is the maximum possible, and | W c ( i ω , ρ a * ) | is not dependent on Θ .

2.2.3. Aperiodic Object with Delay

Parameters of the aperiodic object with delay are determined as a result of the experiment.
I-controller:
ρ a * = 1 + 1 2 T 1 + 1 4 T 2 .
Corresponding controller parameters
S 0 = 1 2 T ρ a * e ρ a * .
Modulus W c ( i ω , ρ a * ) becomes
| W c ( i ω , ρ a * ) | = c T 2 ω 4 + ω 2 ( b 2 2 a T ) + a 2 ,
where
a = ρ a * ( T ρ a * 1 ) , b = 1 2 T ρ a * , c = | S 0 e ρ a * | .
Thus, | W c ( i 0 , ρ a * ) | = 1 . Both the numerator and denominator at ω = 0 yield the same expression
2 T 2 T 1 + 1 4 T 2 .
With the increase in frequency | W c ( i ω , ρ a * ) | tends towards zero because
b 2 > 2 a T .
Indeed,
b 2 = 4 T ρ a * ( T ρ a * 1 ) + 1 ,
2 a T = 2 T ρ a * ( T ρ a * 1 ) .
As the modulus of the extended frequency response of the open-loop system monotonically decreases with the increase of frequency, the limit of stability is considered aperiodic.
PI-controller:
η a * = 1 2 T + 2 τ 1 4 T 2 + 2 τ 2 .
Optimal setting parameters:
S 0 = η k e η τ ( 1 T η ) ,
S 1 = 0.14 e τ T k 4 T τ + 1 .
The module of the extended frequency response of the open-loop system monotonically decreases with the increase of ω , therefore, the limit of aperiodic stability is the maximum possible.
PID-controller:
η * = 1 2 T + 3 τ 1 4 T 2 + 3 τ 2 .
The corresponding setting parameters
S 0 = η 3 τ 2 k e η τ T τ T 2 η + 1 2 ,
S 1 = e η τ k ( 3 T τ η 2 + τ 2 η 2 τ η T τ 2 η 3 1 ) ,
S 2 = τ e η τ k τ η 2 + 2 T η T τ T τ η 2 2 1 .
The modulus of frequency response of the open-loop system decreases monotonically from 1 to 0 while ω increases, thus, η * is the maximum possible.

3. The Robust Stability and Settings of Technological Process Control Systems

In recent years, researchers have paid significant attention to the problem of controlling objects with varying dynamic characteristics and synthesizing robust controllers that can control a whole class of objects or a single object over a wide range of its parameters, loads, and others without reconfiguration. Methods for synthesizing such systems with estimation of permissible range of possible object parameters were developed in studies [6,7,8,9,10,11,12,13,14] and others. Most of these studies are devoted to linearized systems; their characteristic equation’s left-hand side is a polynomial of the form
P n ( a , p ) = a 0 + a 1 p + + a n p n ,
where the polynomial coefficients can take values belonging to some set V a . A polynomial (29) is said to be robustly stable if it is stable for any a V a . Its coefficients a V a should be positive—this is necessary but not sufficient.
In [6,10], the problem of robust stability of the polynomial is solved specifically when the set V a is a parallelepiped delimited by interval constraints
a i ̲ a i a i ¯ , i = 0 , , n .
In [10], four polynomials are constructed, and their values a are selected in such a way that their stability guarantees the stability of P n ( a , p ) . In [11], using the Mikhailov stability criterion based on constraints (30), a system is obtained with characteristics determined by the vector a 0 V a and deviations | a i a 0 | γ α i . It is shown there that it is necessary to construct only one Mikhailov plot of this system to determine whether the original polynomial is robustly stable for a given γ , and at what maximum γ this stability is maintained.
For technological processes with pure delay in their transfer functions, the obtained results are not applicable. These technological processes include those in the chemical, metallurgical, food industries, energy, and others.
The denominator of the transfer function of the closed-loop control system of such objects is not a polynomial, and all results of the automatic control theory (ACT), which are based on the properties of polynomials, are not applicable to such systems. Stability criteria such as Routh–Hurwitz, Mikhailov, logarithmic frequency characteristic methods, state-space methods, and others are not applicable as well. These control system features of technological processes have been repeatedly emphasized by Rotach [15].
The monotonicity of a modulus and phase of technological linearized objects is generally valid, even for open-loop control systems, which, as shown below, simplifies the solution of the problem of robust stability and the selection of robust controller settings.
The characteristic equation of the closed-loop control system is
1 + W ( S , p ) = 0 .
In this equation, W ( S , p ) is the transfer function of the open-loop system, which is equal to the product of the transfer function of the control object W 0 ( p ) and the transfer function of the controller W R ( S , p ) , S V S , where the feedback is negative.
A sufficient condition for the stability to be equal to x 0 is to fulfill the constraint imposed on the modulus of the extended frequency response of the open-loop system:
| W ( S ( x 0 ) , i ω x 0 ) | < 1 ,
where, as shown in the first section, S ( x 0 ) is the solution of the system of equations containing the derivatives with respect to x 0 of the function W ( S , x 0 ) (conditions for the multiplicity of the real root closest to the imaginary axis).
Below are the conditions for the existence of robust controller settings in technological process control systems, equations defining the boundary of the robust D-decomposition, and a methodology for selecting controller parameters for a limited and closed set V a of possible values of the transfer function parameters a of the object. The conditions are specified for systems with typical technological objects and controllers.
Let us consider single-loop systems, where the transfer function of the object depends on the coefficients a V a . We will assume that the open-loop systems are stable or neutral, and that the modulus M and phase φ satisfy the monotonicity conditions
M ( a , S , ω k ) M ( a , S , ω k + 1 ) ,
where ω k is the solution of the equation
φ ( a , S , ω k ) = π ( 1 + 2 k ) , k = 0 , 1 , , φ ω < 0 a V a .
For a technological object and any controller, for which frequency response modulus and phase do not increase with frequency, these conditions are met by default. However, if the controller contains a differentiating element, the monotonicity conditions can be verified by inequality (31).
According to the Nyquist criterion, a dynamic system with feedback is stable if the open-loop system is stable, and the Nyquist plot of the open-loop system W ( a , S , i ω ) = W R ( S , i ω ) W 0 ( a , i ω ) does not encircle the point ( 1 , i 0 ) while varying ω from zero to infinity.
For systems with frequency response modulus that monotonically decreases with increasing ω , this means that the conditions
φ ( a , S , ω 0 ) = π , M ( a , S , ω 0 ) < 1
are satisfied at the first intersection w ith real axis at frequency ω 0 > 0 . For brevity, let us denote the frequency response modulus of the open-loop at its intersection with the negative real axis as M π ( a , S , ω ) .
The expression r ( a , S ) = 1 M π ( a , S , ω ) is referred to as stability margin. Condition (3) can be expressed in the form of inequality M π ( a , S , ω 0 ) < 1 .
Systems satisfying conditions (32) and (33) have a Nyquist plot of the open-loop system, as shown in Figure 4.
Let us provide several definitions related to robust systems.
Definition 1.
1. A closed-loop system is said to be robustly stable if there exists a permissible vector S such that r ( a , S ) > 0 a V a .
For systems that satisfy conditions (32) and (33), this implies the existence of a vector of controller parameters S * V S such that the minimum over S of the maximum M π ( a , C , ω ) over a is less than 1:
min S V S max a V a M π ( a , S , ω ) < 1 .
2. A robust D-partition in the parameter space of the controller is defined as the set of all parameter values for which the closed-loop system is robustly stable. Therefore, the system is robustly stable if the set defined by the D-partition is not empty.
The boundary of the region of robust stability in the parameter space of the controller S is defined by the condition
min S V S max a V a M π ( a , S , ω ) = 1 .
The set V s may include non-negative values of S.
3. Let the system be robustly stable at a certain value of S. Let us define the robust degree of stability η as a non-negative number such that the system that has an extended frequency response of the open-loop system W ( a , S , i ω η ) is stable for all a V a except for the set of values a = a * , for which it is on the stability boundary.
For this system,
max a V a , ω M ( a , S , η , ω ) = 1 / φ ( a , S , η , ω ) = π .
Here, M ( a , S , η , ω ) and φ ( a , S , ω ) are the modulus and phase of the extended frequency response of the open-loop system.
The region bounded by the extended frequency response is a mapping of all the points in the root locus plane of the closed-loop system, lying to the right of a line parallel to the imaginary axis with an x-axis of η . Due to the properties of conformal mapping, this region expands with the growth of η . Therefore, the module of the extended frequency response M ( a , S , η , ω 0 ( η ) ) increases with the growth of η for each fixed value of the phase of the extended frequency response
M π ( a , S , η , ω ) η > 0 .
From inequality (38), it follows that for the system to be robustly stable, it is necessary and sufficient for there to exist such η > 0 for which conditions in (37) are satisfied.
The problem of selecting the controller parameters based on the conditions of maximum robust stability can be formulated as follows, taking into account the introduced definitions:
η * = max S min a V a , ω η ( S , a ) / M π ( a , S , η , ω ) = 1 .
The problem is that the function η ( S , a ) cannot be expressed in analytical form.
On the set determined by the condition M π ( a , S , η , ω ) = 1 , the derivatives
η / a = M π / a M π / η ; η / S = M π / S M π / η
due to condition (38) are opposite in sign to the derivatives of M π over these variables, which leads to an equivalent form of the problem of maximizing robust stability:
min S max a V a , ω M π ( a , S , η * , ω ) = 1 .
Here, η * is known not to exceed (see [16])
η 0 = min a V a , ω max S M π ( a , S , η 0 , ω ) = 1 .
From the inequality η * η 0 , it follows that for the system to be robustly stable, it is necessary to have a positive value η 0 of problem (41).
For processes with the response of the open-loop system monotonically dependent on the frequency, the stability and its margin monotonically depend on each other. This leads to the equation η * и η 0 . This lets us use expressions for maximum stability obtained in the previous section when calculating robust settings.
Problems (39) and (41) are equivalent, and η * = η 0 when the solution of the internal problem in (39) a * does not depend on S. That is, the minimum stability of the system over a is achieved at the same value of the object parameter vector for any controller settings. It happens when the function M π ( a , S , η , ω ) has a multiplicative or additive form.
In the first case,
M π ( a , S , η , ω ) = F a ( a , η , ω ) F S ( S , η , ω ) .
Let us denote the maximum of F a over a V a and the minimum F S over S as F a * ( ω , η ) and F S * ( ω , η ) . The condition of the maximum over ω of the product of these functions being equal to 1 determines the maximum possible robust stability. The order of the maximum and minimum operations does not affect the form of the functions F a * и F S * , which means that problems (39) and (41) are equivalent, and the robust stability can be found by initially selecting the controller settings based on the minimum condition M π (or maximum η ) for any admissible object parameters, and then finding the maximum M π (or minimum η * ( a ) ) over a V a . Similar reasoning is applicable to the additive form of the function M π .
The modulus of the frequency response of the open-loop control system is a product of the moduli of the frequency characteristics of the object, which depend only on the parameters a, and the controller, which depends on the vector S. Therefore, this expression is multiplicative, and the achievable stability margin does not depend on the order of operations of finding the maximum of this modulus over a and the minimum over S. Because the stability and its margin are monotonically dependent on each other, the same statement holds for the stability itself. This implies an algorithm for calculating robust controller settings.
It should be noted that choosing the controller parameters based on the maximum robust stability condition is more natural than choosing based on the maximum stability margin condition as η * is directly related to the duration of the transition in the system [17].

3.1. Algorithm for Selecting Robust Controller Settings for Technological Processes

1. The dependences of the controller’s optimal settings S * ( a ) and corresponding maximum stability η * ( a ) on the parameters of the transfer function of the control object are found based on the conditions of the proximity of real roots of the closed-loop control system to the imaginary axis (Appendix A Table A1).
2. The minimum value of the function min a V a η * ( a ) and its corresponding (critical) parameter values a * are found. If the obtained minimum value η * ( a * ) is positive, the system is robustly stable, and the corresponding settings S * ( a * ) are the desired ones.
It should be noted that the set V a can be any closed and bounded set, and η * ( a ) is a continuous and bounded from the function below, which guarantees the existence of a minimum.

3.2. Robust Control System for an Aperiodic Object with Pure Delay

As an example, let us consider a system consisting of an aperiodic object with delay and a PI-controller.
The transfer function of the open-loop system is
W ( a , p , S ) = K e p τ T p + 1 W R ( S , p ) , a = ( K , T , τ ) .
Here, S is the vector of controller parameters.
The limit of stability is determined by expression (25):
η * = 1 2 T + 2 τ 1 4 T 2 + 2 τ 2 .
The constraints on T and τ highlight the set V 0 of their possible values. Figure 5 shows how critical values T * and τ * are determined in this case. For these values, η * is minimal (usually corresponding to the maximum of the ratio τ / T ).
After substituting the critical object parameters into expressions (26) and (27), the robust controller settings are determined.

3.3. Example

Consider a closed system with a PI controller and an object with a transfer function
W ( p ) = e p 1.5 p + 1 .
The calculation of the limiting degree of stability and the corresponding parameters of the regulator using the formulas given in the Table gives
η * = 0.88 , S 0 = 0.37 , S 1 = 0.56 .
The transient process corresponding to a single perturbation and the maximum degree of stability is shown in Figure 6, curve (a).
Let the parameters of the object be changed within: 0.5 τ 1.5 , 1 T 2 . Robust settings correspond to the values
T = 1 , τ = 1 , 5 , η r = 0.77 , S 0 = 0.25 , S 1 = 0.25 .
The transient process in a closed system with and robust settings is shown in the same figure, curve (b).

4. Main Results

Conditions were obtained under which the limit of aperiodic stability was achieved. The problem of robust stability for single-loop linear feedback systems was considered. In these systems, the modulus and phase of the frequency response of the open-loop system monotonically decrease with frequency, and the object contains delay, so that the denominator of the transfer function of the closed-loop system is not a polynomial. An algorithm for robust selection of controller settings was proposed for control systems and technological objects containing delay.
It was proven that it is possible to choose robust settings for typical industrial controllers by minimizing the set of possible object parameters (not necessarily characterized by interval constraints) for the maximum achievable stability, which depends on these parameters.
It was shown that for positive bounded values of τ * , T * of the transfer function of typical technological objects, the stability η * is positive. This means that there are parameters of typical controllers that ensure robust stability. The obtained conditions allow for solving the inverse problem: to find the range of variation of object parameters for which the robust stability of the control system will be equal or higher than the specified one. This makes it possible to design controllers with fixed settings for an entire class of objects.

Author Contributions

Conceptualization, A.M.T.; methodology, A.M.T. and A.I.B.; validation, A.I.B.; writing—original draft preparation, A.I.B.; writing—review and editing, A.M.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. The settings that correspond to the maximum of aperiodic stability η for a typical technological process.
Table A1. The settings that correspond to the maximum of aperiodic stability η for a typical technological process.
Controller η
S
W ( p ) = k e p τ T p + 1 W ( p ) = e p τ Θ p W ( p ) = k e τ
12345
P η 1 τ + 1 T 1 τ
S 1 0.37 T k τ e τ T 0.37 Θ τ
PI η 1 2 T + 2 τ 1 4 T 2 + 2 τ 2 0.59 τ 2 τ
S 1 e η τ k ( 2 T η + τ η T τ η 2 1 )
0.23 k 2 T τ + 1.6 e 0.9 T τ 1
0.54 Θ τ 0.14 k
S 0 τ η 2 e η τ k T τ T η + 1
0.19 k τ 0.41 T τ + 1.83 e 1.2 T τ + 1
0.08 Θ τ 2 0.54 τ k
I η 1 2 T + 1 τ 1 4 T 2 + 1 τ 2 1 τ
S 0 η k e η τ ( 1 T η ) 0.37 τ k
PD η 1 T + 2 τ 2 τ
S 1 0.14 e τ T k 4 T τ + 1 0.54 Θ τ
S 2 0.14 T k e τ T 0.14 Θ
PID η 1 2 T + 3 τ 1 4 T 2 + 3 τ 2 1.27 τ 3 τ
S 1 e η τ k ( 3 T τ η 2 + τ 2 η 2 τ η T τ 2 η 3 1 )
0.19 k 4.13 T τ + 2.3 e 1.4 T τ 1
0.79 Θ τ 0.25 k
S 0 η 3 τ 2 k e η τ T τ T 2 η + 1 2
0.29 k τ 0.73 T τ + 1.33 e 1.5 T τ + 1
0.19 Θ τ 2 0.67 τ k
S 2 τ e η τ k τ η 2 + 2 T η T τ T τ η 2 2 1
0.1 τ k 2 T τ + 1.25 e 1.1 T τ 1
0.2 Θ 0.025 τ k

References

  1. Letov, A.M. Ustoichivost nelineinykh reguliruemykh sistem. In Stability of Nonlinear Regulated Systems; Fizmatgiz: Moscow, Russia, 1962; 484p. (In Russian) [Google Scholar]
  2. Tsypkin, Y.Z. On the upper boundary of stability of single-loop automatic control systems. Autom. Remote Control 1952, XIII, 425–428. [Google Scholar]
  3. Tsypkin, Y.Z.; Bromberg, P.V. On the degree of stability of linear systems. Izv. AN SSSR. Ser. OTN 1945, 6, 1163–1165. (In Russian) [Google Scholar]
  4. Shubladze, A.M. Methods for synthesizing control systems of maximum stability. Autom. Remote Control 1980, XXXXI, 28–37. [Google Scholar]
  5. Zagarii, G.I.; Shubladze, A.M. Synthesis of Control Systems Based on the Criterion of Maximum Stability; Energoatomizdat: Moscow, Russia, 1988; 170p. (In Russian) [Google Scholar]
  6. Polyak, B.T.; Shcherbakov, P.S. Robust Stability and Control; Nauka: Moscow, Russia, 2002; 303p. (In Russian) [Google Scholar]
  7. Tatarinov, A.V.; Tsirlin, A.M. Mathematical Programming Problems with Complex Variables and the Stability Limit of Dynamical Systems; Izvestiya RAN, Teoriya i Sistemy Upravleniya: Moscow, Russia, 1995; pp. 28–33. (In Russian) [Google Scholar]
  8. Dudnikov, E.G. Fundamentals of Automatic Control of Thermal Processes; Gosudarstvennoe Energeticheskoe Izdatel’stvo: Moscow, Russia, 1956; 264p. (In Russian) [Google Scholar]
  9. Gorecki, H.; Metzger, M. Analysis and synthesis of time-lag control systems. IEEE Trans. Syst. Man Cybern. 1974, SMC-4, 594–595. [Google Scholar] [CrossRef]
  10. Kharitonov, V.L. Asymptotic stability of a family of linear differential equations. Differ. Equ. 1978, 1, 2086–2088. [Google Scholar]
  11. Polyak, B.T.; Tsypkin, Y.Z. Robust Nyquist criterion. Autom. Remote Control 1992, XXXXIII, 25–31. [Google Scholar]
  12. Ulyanov, S.; Litvintseva, L.; Dobrynin, V.; Mishin, A. Moskau, Intelligent Robust Control; ProtentLabs: Moscow, Russia, 2011. (In Russian) [Google Scholar]
  13. Tatarinov, A.V.; Tsirlin, A.M. The limiting degree of aperiodic stability of linear systems and the choice of parameters of industrial regulators. Model. Anal. Informatsionnykh Sist. 2012, 19, 87–96. (In Russian) [Google Scholar] [CrossRef]
  14. Popov, A.Y.; Pukhov, S.S.; Tsirlin, A.M. Minimization of the Least Upper Bound of the Real Parts of Quasipolynomial Roots and the Limit Stability Degree of Linear Dynamic Feedback System. Autom. Remote Control 2015, 76, 1137–1155. [Google Scholar] [CrossRef]
  15. Rotach, V.Y. Theory of Automatic Control; ICC MPEI: Moscow, Russia, 2008; 396p. (In Russian) [Google Scholar]
  16. Davydov, E.G. Issledovanie operatsii. In Operations Research; Graduate School: Moskva, Russia, 1990; 384p. (In Russian) [Google Scholar]
  17. Tsypkin, Y.Z. Fundamentals of Automatic Systems Theory; Nauka: Moscow, Russia, 1977; 560p. (In Russian) [Google Scholar]
Figure 1. Control system structure.
Figure 1. Control system structure.
Processes 11 01835 g001
Figure 2. The limit of aperiodic stability may (a) or may not (b) coincide with the maximum stability.
Figure 2. The limit of aperiodic stability may (a) or may not (b) coincide with the maximum stability.
Processes 11 01835 g002
Figure 3. The Nyquist plot of the extended response of an open-loop system for cases where the maximum degree of stability is aperiodic (a) and oscillatory (b).
Figure 3. The Nyquist plot of the extended response of an open-loop system for cases where the maximum degree of stability is aperiodic (a) and oscillatory (b).
Processes 11 01835 g003
Figure 4. Nyquist plot of the open-loop control system.
Figure 4. Nyquist plot of the open-loop control system.
Processes 11 01835 g004
Figure 5. Set V a of possible values of T and τ and the choice of their critical values.
Figure 5. Set V a of possible values of T and τ and the choice of their critical values.
Processes 11 01835 g005
Figure 6. Process examples.
Figure 6. Process examples.
Processes 11 01835 g006
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tsirlin, A.M.; Balunov, A.I. Optimal Control of Technological Processes. Processes 2023, 11, 1835. https://doi.org/10.3390/pr11061835

AMA Style

Tsirlin AM, Balunov AI. Optimal Control of Technological Processes. Processes. 2023; 11(6):1835. https://doi.org/10.3390/pr11061835

Chicago/Turabian Style

Tsirlin, Anatoliy M., and Alexander I. Balunov. 2023. "Optimal Control of Technological Processes" Processes 11, no. 6: 1835. https://doi.org/10.3390/pr11061835

APA Style

Tsirlin, A. M., & Balunov, A. I. (2023). Optimal Control of Technological Processes. Processes, 11(6), 1835. https://doi.org/10.3390/pr11061835

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop