1. Introduction
Typical uncertainty control problems are presented for one feedback loop, but this paper presents uncertainty analysis for two feedback loops in the cascade control system. This article describes the change in the eigenvalue placement when the parameters of the mathematical model are changed (are uncertain).
Mathematical model parameters in every plant are varying when it works, so they are uncertain. Additionally, the parameters identification method are not perfect, so model parameters are uncertain also. This, uncertainty analysis is very important in every automatic control system. Uncertainty in theory leads to the following methods:
Classical Nyquist criterion which leads to gain and phase margin [
1,
2],
Hardy space (
) analysis and synthesis [
2,
3],
Kharitonov theorem [
4,
5],
relative stability [
6,
7,
8]
In this paper Kharitonov and relative stability are used to cascade control of Brushless Direct Current (BLDC) motor. The motor model is similar to synchronus (with and without permanent magnets), separately-excited DC and induction motor [
9]. The closed-loop system eigenvalues for nominal transfer-function parameters of a plant are in certain sector in Gauss’s plane, but if parameters of a plant are uncertain, then sector will be larger—it is
relative stability.
The paper presents model of the motor and the cascade system with uncertainty of the resistance, inductance, flux and gain parameters. Moreover, relative stability is analized, thus Kharitonov and -stability theory is used.
Article presents:
relative stability theory which is more stringent than asymptotic stability [
4,
6,
7,
8,
10],
model of BLDC motor with power converter,
uncertainty of model parameters,
cascade control system, where transfer-function of the closed-loop system is analyzed, this 4th order function is the basis of -stability analysis in the paper,
the boundaries of the transfer-function coefficients,
relative stability analysis of the uncertainty closed-loop cascade control,
some experimental results of the BLDC motor speed control.
The asymptotic stability testing, but not relative stability, for electric drives is used among other papers [
11,
12,
13,
14,
15,
16].
2. Idea of Relative Stability and Generalization of Routh Theorem
The theory and method for testing relative stability (
-stability) is described in this section.
-stability [
6] of linear feedback system is described by angle
between imaginary axis Im and sector’s band (
Figure 1a). Sector
determination for control system poles leads to minimal value of relative damping factor
, therefore overshoot of step response.
A development method for study stability from [
10] is article [
6], where has been presented another approach. The linear time-invariant dynamic control system is relatively
-stability, if all roots of characteristic equation are located inside specified area (
angle or damping coefficient
). In literature [
4], this area called
.
For characteristic polynomial
, which can be written in a short form:
where coefficients
have real value, and coefficient
.
If a control system is relatively stable then relative-damping coefficient
is greater than the smallest, real relative-damping coefficient results from polynomial coefficient (
1):
where:
and
is natural undamping frequency.
For value
=0 control system described by (
1) is stable, because
area, contained all left-half complex variable
s (
Figure 1). Make substitution as follows:
to (
1) received complex polynomial of the variable
in the form:
The roots of (
1) are located from left side of straight line
in
s-plane (
Figure 1a), if all roots of complex polynomial (
4) are inside the left-half
variable (
Figure 1b).
On the other hand, applying substitution:
to (
1) we can receive the complex polynomial:
The roots of real polynomial (
1) are located on the left side of straight line
in
s-plane, if all roots of after substitution (
6) are inside the left-half
variable (
Figure 1b).
A complex polynomial (
4) and (
6) are obtained:
where:
and
are coefficients in the following form:
where:
Polynomial (1), whose roots lie in area, is relatively -stable if and only if or areHurwitz polynomial.
To check, do a polynomials
or
, which have complex coefficients, have all roots with negative real parts, can apply generalization of the Hurwitz criterion, which was discussed in [
10], so the Routh criterion generalization can be used, too. The roots of both polynomials (
4) and (
6) have complex roots, which do not conjugate complex pairs. While roots of
are conjugate with correspodning roots of with
.
In order to study only -stability generalization of the Routh criterion, it is necessary to build a modified Routh array, which has real coefficient.
Given polynomials (
7) and (
8) done multiply both polynomials in order to find real coefficients. Thus, the following polynomial is obtained:
where:
are real coefficients:
Using (
9), the following relationships are obtained:
After substituting (
12) to (
11) the coefficients
of polynomial
are in the form:
Thus, the Routh array has
rows:
Condition of the
-stability, using the
generalization Routh criterion, is no change of sequence sign:
in the first column of the Routh array (
14).
3. Mathemathical Model of BLDC Motor
The
Figure 2 shows diagram of a BLDC motor supplies by a voltage source inverter.
The mathematical model of BLDC motor is in the following form [
17]:
where
is mutual inductance between stator coils.
Considering the principle of the motor control (at any moment, the current flows through two stator windings) which result from the design, it is possible to show (present courses of) currents and electromotive forces as in
Figure 3, where
is the position of the rotor.
At any time of the control, the current flows through the two stator windings, thus it is assumed:
The transistors of voltage source inverter (
Figure 2) are controlled by the pulse width modulation (PWM) method, with modulation factor
, so
.
Finally, the differential equations can be written:
where
E is back EMF of the two windings,
J is moment of inertia,
is permanent magnets flux,
is electromagnetic torque,
is a load torque.
4. Uncertainty of Model Parameters
The BLDC motor considered in the paper is characterized by a model error, which is due to the uncertainty of the parameters. These errors result from: parameters identification, the conditions and nature of the motor operation.
One of the parameters is the stator windings resistance. This value, apart from the identification errors, changes as the motor temperature increases during operation, according to the relation:
where:
= 293 K⇒ 20 C,
—operating temperature,
—rated operating motor temperature, depend on the insulation class,
—temperature coefficient of resistance in Cu .
Another parameter that affects the model’s uncertainty is the moment of inertia J, which, in extreme cases can change during operation. For example, a sheet metal reeler. The article omits changes to this parameter.
The inductance of motor windings, which is defined as the ratio of the flux to the current I also contributes to the mathematical model error. In the initial start-up flux and current values change so that their ratio is constant. However, in the further stage of start-up, the flux saturates, while the current arises, which causes the L inductance to decrease.
The above parameters have a particular impact on the values of time constants:
- -
, where T is electromagnetic time-constant,
- -
, where B is electromechanical time-constant,
that directly affect the nature of the drive operation.
5. BLDC Motor—Speed Control System
Angular speed of the BLDC motor is controlled by using a cascade control structure, where the control of the inverter transistors depends on the current position of the rotor (feedback ). In this type of cascade control, the motor current is usually measured in the intermediate circuit of the inverter, the position of the rotor is measured by the Hall sensors placed on the motor stator, and a incremental encoders are used to measure accurate angular speed.
The control structure shown in
Figure 4 consists of a master controller (angular speed control) and the slave controller that is responsible for controlling the electromagnetic torque or quantity proportional to it. In direct current drives, the quantity proportional to the electric torque is the current which is defined as:
. The system works in a closed-loop, so that the gains of measuring are also included:
Y is the gain of the current measuring,
is the gain of the angular speed measurement ,
is the feedback path, measuring the current position of the rotor.
The block diagram of the BLDC drive (
Figure 5) was determined from the mathematical model (
18), from the diagram in
Figure 4 and the inverter gain
.
The PI current controller assumes the form:
and the PI speed controller:
Thus, closed-loop transfer-function is equal:
which leads to:
where:
The article considers the uncertainty of the following parameters:
. By applying the parameter limits (lower or upper) for the denominator coefficients of the closed-system transfer-function, an interval polynomial is obtained in the form:
where:
are the lower and upper ranges of coefficients (
24) and
. Their values are determined according to the following approach:
boundary coefficient
at
:
boundary coefficient
at
:
boundary coefficient
at
:
boundary coefficient
at
s:
boundary coefficient
:
6. Uncertainly Parameters of BLDC Drive
The BLDC motor drive with the following rated parameters is considered:
= 5.25 []—single-phase resistance,
= 0.46 [mH]—single-phase inductance,
E = 2.62 []—electromotive force,
= 24.9 []—torque constant,
= 24.9 []—emf constant,
J = 9.9 · [kgm]—moment of inertia.
The ranges of the parameters uncertainty for the tests according following forms:
Thus, the interval polynomial coefficients are:
If the current controller (
20) parameters are:
and for velocity controller (
21):
This rated closed-loop transfer-function (
23) is:
The obtained transfer-function
, more precisely the characteristic polynomial, is subjected to algebraic studies of the relative
-stability and the value of the relative damping factor is
, which corresponds to the angle
. Upper values determine the boundary sector of relative stability (
Figure 6), “to which the results of changes of individual parameter values will refer”. to which there will be referred the results of the particular parameter value changes.
6.1. Interval of the Stator Resistance
The first uncertain parameter is the resistance of the motor stator, which is in the range:
Applying the minimum resistance value
to the transfer-function coefficients (
23) equation can be written:
Analyzing the denominator of transfer-function (
29), we get the relative damping factor
, corresponding to the angle
. On the other hand, if we use the maximum value of resistance
to the transfer-function coefficients (
23), we obtain the transfer-function of a closed circuit in the form:
where dumping factor and angle .
The
Figure 7 shows sectors of relative stability determined by characteristic polynomial roots, for different resistance values.
After sector analysis with disturbed parameters and referring them to the nominal sector the system is relative -stable for a resistance lower than the nominal one (), because this sector is located inside the nominal sector . When the resistance reaches the maximum value of , study object will lose relative stability because sector is located outside the nominal sector (all roots aren’t inside the sector). Thus, for the system has greater overshoot and oscillations during start-up than for rated parameters. On the other hand, for the the system has the smallest overshoot, which is confirmed by the relative damping factor .
The simulation results (in Matlab environment) obtained for angular velocity (
Figure 8) and current of the inverter (
Figure 9) confirm
-stability analysis.
These results confirm validity of the statements regarding relative damping factor , whose value represents the overshoot value in the real system. Each waveform, for a resistance value of , or , has the lowest overshoot, and for a resistance value of , or , has the highest overshoot.
Thus, the statements of the relative damping factor
, whose the value corresponds to the overshoot in the real system, are confirmed (
Figure 8 and
Figure 9).
6.2. Interval of the Inductance
The range of inductance variation is:
thus, transfer-function for
and
are in the following forms:
After analysing above transfer-functions, we get:
for (
31):
,
for (
32):
.
So -stability is always fulfilled and the simulation studies confirmed above results and conclusions.
6.3. Interval of Inverter Gain
The range of variability inverter parameter
is determined by the voltage on capacitor in DC circuit:
where transfer-function (
23) for
assumes form:
Analysis of the denominator (
33) gives the relative damping factor
corresponding to the angle
.
Applying
substitution, transfer-funtcion (
23) is in the following form:
where for
damping factor eqauls
and angle
.
Figure 10 shows relative stability sectors determined by the roots of characteristic polynomials for different values of the inverter gain
.
Analysis of the received sectors with disturbed parameters and reference to the nominal sector, shows that the system is relative
-unstable for gain value less than nominal (
), because this sector is located outside of the nominal sector. In second case the gain value is
and the system is relative stable because the sector
is within the nominal sector. The simulations (
Figure 11) confirm the studies results. For
, the control system has the lowest overshoot compared to the speed waveforms at rated and minimum
gain value, which is presented in
Figure 11.
6.4. Interval of the Flux Linkage
The last parameter tested is flux linkage
, which takes limit values in the range:
and gets a transfer-function for
:
According to the analysis of the characteristic polynomial (
35), the result is
and
. This polynomial compared to the characteristic transfer-function polynomial (
28) is relative
-unstable (
Figure 12).
The transitions of angular speed, shown in the
Figure 13, based on the cascade control structure (
Section 5) confirms that the system is
-unstable for transfer-function (
35).
7. Kharitonov Theorem
One of the most popular methods of testing uncertain dynamic systems is the Kharitonov theorem. Unfortunately, successive maximum and minimum values of coefficients in the interval characteristic equations are often mutually exclusive (in one equation). Therefore, the theorem very quickly allows to determine asymptotic or relative stability, but the results are approximate.
The interval polynomial is in the form [
18]:
where
are a range of uncertain polynomial coefficients and
.
Theorem 1 (Kharitonov Theorem).
Every polynomial from the polynomial family (36) is stable if and only if the following four Kharitonov polynomials are stable: The stability conditions of the low degrees interval polynomials are given in [
5] and for an uncertain fourth-degree polynomial
The interval polynomial is stable if the following inequalities occur:
and polynomials
are stable, so from Liénarda-Chiparta theorem [
18] stability conditions are in the form
Thus, stability condition for considered interval polynomial are (
42), (
45) and (
46).
After inserting the coefficients of interval polynomials
for Kharitonov’s polynomial (
41) there are obtained four polynomials (
37), (
38), (
39) and (
40) representing the worst conditions of the BLDC drive:
for (
37):
,
for (
38):
,
for (
39):
,
for (
40):
,
Comparing received sectors with a nominal sector
(
28), the drive is
-unstable in each case. The results obtained even show the possibility of loss of asymptotic stability by the control system, while the real results are in the open left Gauss’s half plane.
8. Experimental Results
Low power BLDC motor (rated voltage is 15 V, current is 1.1 A) was used in laboratory stand (Motion Control Kit F240 produced by Technosoft in 1998 year [
19]), thus the experiments were safe. MCK243 is a complete motion kit, including a voltage source inverter and a three-phase permanent magnet synchronous motor with incremental encoder and Hall sensors. The main element of this kit is the TMS320F243 Digital Signal Processor (DSP) for digital motion control applications.
Unfortunately, not all calculations could be proved practically. The following section presents the change of and R. For a better representation of -stability laboratory experiments were realized for an increased range of parameter changes.
Testing the changes of resistance is difficult, so for
additional resistances were used and a start-up was performed. However, for
, the controllers were optimized for a higher stator resistance, and the experiment was carried out for the rated one—this case is an artificial proof of previous calculations. The experimental results are shown in
Figure 14.
Obtained results are similar to presented in
Figure 8.
The next stage of laboratory tests was the change of voltage on the capacitor in the DC circuit, which was considered in the analysis as a change
. Experiments were carried out for 12, 15 and 25 V, again the results for a larger range of parameter changes were observed. The aim was to observe the effects of voltage changes in the DC circuit (
Figure 15)—the step responses for 12 V and 15 V are practically the same, so for 12 V is omitted.
Obtained results are similar to presented in
Figure 11, so laboratory results confirm previous mathematical analysis.
Experimental verification of the remaining theoretical calculations and simulations is difficult on a real motor, but the presented experiments of uncertainty the most important parameters R, showed the correctness of theoretical calculations.
The next experiments are not as simple because changing the inductance of the motor windings is difficult to implement. For this reason, optimization of the controllers for the incorrectly identified inductance
L was done. It was assumed, in the optimization of the controllers where
mH (the inverse of the uncertainty assumptions of
L) and the experiment was performed for the nominal value of
L. The results are shown in
Figure 16.
The obtained results slightly differ from the theory, possibly due to the nonlinearity of the low-power motor tested. In addition, the obtained result can be interpreted as: what happens if L is incorrectly identified and used to adjust controllers.
The next step was to introduce multiple perturbations of the motor model parameters to selection of the controllers settings and additionally lower supply voltage was used. Two experiments were realized/prepared for:
- (1)
V,
- (2)
V.
Experiment 1 shows less damping ratio than the nominal case, so the system is not
-stable. This is similar to the case shown in
Figure 18. On the other hand, experiment 2, where the supply voltage is reduced, shows that the obtained signal is deformed (almost flat part for
t = 0.06–0.14 s) and this is the effect of voltage limitation which is not considered in the calculation.
9. Conclusions
The analytical approach was used to determine the fourth-order transfer function for a cascade control system with a BLDC motor, for which the limits of changes of the numerator and denominator coefficients were determined.The denominator polynomial was used to calculate the relative stability of the closed-loop control system for different cases (numerical calculations).
Simulation researches confirmed the correctness of analytical-numerical calculations from the first part of the article. Later, experimental results were presented.
The article presents -stability test for the cascade control structure of the BLDC motor drive, where the influence of changes in individual parameters (uncertainty) of the motor model and voltage inverter amplification were examined. The biggest influence on the relative stability is the change of resistance, followed by the amplification of the converter and the combined rotor flux. The inductance of stator windings has practically no effect on the -stability of the cascade control structure and this is the most important result of this work.
The results presented in the paper can be extended to electric drives with PMSM or DC motors. Furthermore, the considered closed-loop transfer-function can be used for uncertainty studies by frequency domain methods, although the open-loop transfer-function (Nyquist criterion) is more appropriate in this case.