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,
и
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
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
If
, then the controller is proportional-integral (PI), if
and
are zero, it is integral (I), and if
and
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 (), 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
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:
The dependency
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
which, for the
nth power equations, according to Vieta’s formulas, is equal to the coefficient at
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
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
th critical root
where
m is the number of required controller settings (in Equation (
3),
).
For
, the left-hand side of Equation (
3) can be reduced by
p. At the same time, a root with
, and changing settings
does not change
, which leads to a steady-state error. This leads to abrupt changes of
at
.
If Equation (
6) can be solved when
, then the corresponding optimal settings are found as
For
, settings
and
, the PD-controller is chosen in such a way that the roots of the equation
satisfy requirement (
4). The limit of aperiodic stability (taking into account all non-zero roots) and the corresponding settings are found from the conditions:
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
. This leads to
Problem B: Under which conditions does the limit of aperiodic stability 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
, and
. They correspond to the frequency response of the open-loop system
and its extended frequency response
Figure 2 shows the location of the roots of the characteristic equation of the system in cases where
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
. It is clear that due to condition (
3) at
, this plot will pass through the point (
) (
Figure 3).
If all the characteristic equation roots of the system lie to the left of the line , then when ω changes from zero to infinity, the Nyquist plot will not encircle the point (
Figure 3a).
Let us denote the frequencies with the phase of frequency equal to as critical frequencies . In this case, corresponds to a rotation of the plot by by , and so on.
At the frequency .
The properties of conformal mapping determine the following.
Statement:
The limit of aperiodic stability is the maximum possible stability of the system only if Particularly, for optimality , it is sufficient for the modulus of to decrease monotonically while increases. Conversely, the limit of aperiodic stability is not optimal if 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
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 .
Let us note that from the transfer function (
16) after substituting
and subsequent limit transition at
we obtain an integrating object with delay
By directing
T to zero in (
16), we obtain a pure delay object
For each of the noted typical objects, analytical expressions for the achievable aperiodic stability
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:
The only value of the setting parameter
corresponds to
. The modulus of the extended frequency response of the open-loop system
decreases monotonically from 1 at
= 0 to zero at
→∞. In this case, the limit of aperiodic stability is the maximum possible.
PI-controller:
The corresponding setting parameters
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 .
PID-controller:
The corresponding setting parameters
The modulus of the extended frequency response of the open-loop system
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:
Optimal controller parameters
Modulus
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:
Corresponding controller parameters
The modulus of the extended frequency response of the open-loop system
It is easy to see that the limit of aperiodic stability is the maximum possible, and
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:
Corresponding controller parameters
Modulus
becomes
where
Thus,
. Both the numerator and denominator at
yield the same expression
With the increase in frequency
tends towards zero because
Indeed,
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:
Optimal setting parameters:
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:
The corresponding setting parameters
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
where the polynomial coefficients can take values belonging to some set
. A polynomial (
29) is said to be robustly stable if it is stable for any
. Its coefficients
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
is a parallelepiped delimited by interval constraints
In [
10], four polynomials are constructed, and their values
a are selected in such a way that their stability guarantees the stability of
. In [
11], using the Mikhailov stability criterion based on constraints (
30), a system is obtained with characteristics determined by the vector
and deviations
. 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
In this equation,
is the transfer function of the open-loop system, which is equal to the product of the transfer function of the control object
and the transfer function of the controller
, where the feedback is negative.
A sufficient condition for the stability to be equal to
is to fulfill the constraint imposed on the modulus of the extended frequency response of the open-loop system:
where, as shown in the first section,
is the solution of the system of equations containing the derivatives with respect to
of the function
(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 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
. We will assume that the open-loop systems are stable or neutral, and that the modulus
M and phase
satisfy the monotonicity conditions
where
is the solution of the equation
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 does not encircle the point while varying from zero to infinity.
For systems with frequency response modulus that monotonically decreases with increasing
, this means that the conditions
are satisfied at the first intersection w ith real axis at frequency
. For brevity, let us denote the frequency response modulus of the open-loop at its intersection with the negative real axis as
.
The expression is referred to as stability margin. Condition (3) can be expressed in the form of inequality .
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 .
For systems that satisfy conditions (32) and (33), this implies the existence of a vector of controller parameters such that the minimum over S of the maximum over a is less than 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 conditionThe set 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 is stable for all except for the set of values , for which it is on the stability boundary.
For this system,
Here,
and
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
increases with the growth of
for each fixed value of the phase of the extended frequency response
From inequality (
38), it follows that for the system to be robustly stable, it is necessary and sufficient for there to exist such
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:
The problem is that the function
cannot be expressed in analytical form.
On the set determined by the condition
, the derivatives
due to condition (
38) are opposite in sign to the derivatives of
over these variables, which leads to an equivalent form of the problem of maximizing robust stability:
Here,
is known not to exceed (see [
16])
From the inequality
, it follows that
for the system to be robustly stable, it is necessary to have a positive value 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 и . This lets us use expressions for maximum stability obtained in the previous section when calculating robust settings.
Problems (
39) and (
41) are equivalent, and
when the solution of the internal problem in (
39)
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
has a multiplicative or additive form.
Let us denote the maximum of
over
and the minimum
over
S as
and
. 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
и
, 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
(or maximum
) for any admissible object parameters, and then finding the maximum
(or minimum
) over
. Similar reasoning is applicable to the additive form of the function
.
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
and corresponding maximum stability
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 and its corresponding (critical) parameter values are found. If the obtained minimum value is positive, the system is robustly stable, and the corresponding settings are the desired ones.
It should be noted that the set can be any closed and bounded set, and 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
Here,
S is the vector of controller parameters.
The limit of stability is determined by expression (
25):
The constraints on
T and
highlight the set
of their possible values.
Figure 5 shows how critical values
and
are determined in this case. For these values,
is minimal (usually corresponding to the maximum of the ratio
).
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
The calculation of the limiting degree of stability and the corresponding parameters of the regulator using the formulas given in the Table gives
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:
Robust settings correspond to the values
The transient process in a closed system with and robust settings is shown in the same figure, curve (b).