Next Article in Journal
Optimizing the Geometric Parameters of a Straight-Through Labyrinth Seal to Minimize the Leakage Flow Rate and the Discharge Coefficient
Next Article in Special Issue
Optimization of Powertrain in EV
Previous Article in Journal
Effects of SiO2 Nanoparticle Dispersion on The Heat Storage Property of the Solar Salt for Solar Power Applications
Previous Article in Special Issue
Experimental Study of Harmonic Influence on Electrical Energy Metering
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Relative Stability of Electrical into Mechanical Conversion with BLDC Motor-Cascade Control

Department of Power Electronics and Energy Control Systems, AGH University of Science and Technology, 30-059 Krakow, Poland
*
Author to whom correspondence should be addressed.
Energies 2021, 14(3), 704; https://doi.org/10.3390/en14030704
Submission received: 4 December 2020 / Revised: 18 January 2021 / Accepted: 24 January 2021 / Published: 29 January 2021

Abstract

:
The conversion of the electrical energy into the mechanical is usually realized by a motor, power electronics and cascade control. The relative stability ( Θ -stability), i.e., the displacement of its eigenvalues of this system is analyzed for a drive with a BLDC motor. The influence of changing the basic parameters of the motor and power supply system on the drive operation is considered. 4th order closed-loop transfer-function of the cascade control is presented, where boundaries of the transfer-function coefficients are used. The cascade system which uncertainty of the resistance, inductance, flux and gain parameters is analyzed. Theoretical calculations for the cascade control, simulations and laboratory tests are included in the article.

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 ( H ) 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 S θ 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 S θ .
For characteristic polynomial M ( s ) , which can be written in a short form:
M ( s ) = k = 0 n a k s k = 0
where coefficients a k have real value, and coefficient a n > 0 .
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):
s = ω 0 e j Ψ = ξ ω 0 + j ω 0 1 ξ 2 ξ = cos ( Ψ ) = c o n s t
where: 0 ξ 1 and ω 0 is natural undamping frequency.
For value ξ =0 control system described by (1) is stable, because S θ area, contained all left-half complex variable s (Figure 1). Make substitution as follows:
s = ω e j θ = ω e j Ψ π 2
to (1) received complex polynomial of the variable ω in the form:
M 1 ( ω ) = k = 0 n a k e j k θ ω k
The roots of (1) are located from left side of straight line l 1 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:
s = ω e j θ = ω e j Ψ π 2
to (1) we can receive the complex polynomial:
M 2 ( ω ) = k = 0 n a k e j k θ ω k
The roots of real polynomial (1) are located on the left side of straight line l 2 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:
M 1 ( ω ) = k = 0 n a k e j k θ ω k = k = 0 n b k + j c k ω k
M 2 ( ω ) = k = 0 n a k e j k θ ω k = k = 0 n b k j c k ω k
where: b k and c k are coefficients in the following form:
b k = a k c o s k Θ c k = a k s i n k Θ
where: k = 0 , 1 , 2 , , n
Polynomial (1), whose roots lie in S θ area, is relatively θ -stable if and only if M 1 ( ω ) or M 2 ( ω ) areHurwitz polynomial.
To check, do a polynomials M 1 ( ω ) or M 2 ( ω ) , 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 M 1 ( ω ) are conjugate with correspodning roots of with M 2 ( ω ) .
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:
D ( ω ) = M 1 ( ω ) · M 2 ( ω ) = i = 0 2 n d i ω i
where: d i are real coefficients:
d 0 = b 0 2 + c 0 2 d 1 = 2 b 0 b 1 + c 0 c 1 d 2 = b 1 2 + c 1 2 + 2 b 0 b 2 + c 0 c 2 d 2 = 2 b 1 b 2 + c 1 c 2 + 2 b 0 b 3 + c 0 c 3 d 2 n 1 = 2 b n 1 b n + c n 1 c n d 2 n = b n 2 + c n 2
Using (9), the following relationships are obtained:
b i 2 + c i 2 = a i 2 , b i b j + c i c j = a i a j cos j i θ
After substituting (12) to (11) the coefficients d i of polynomial D ( ω ) are in the form:
d 0 = a 0 2 d 1 = 2 a 0 a 1 c o s θ d 2 = a 1 2 + 2 a 0 a 2 c o s 2 θ d 3 = 2 a 1 a 2 c o s θ + 2 a 0 a 3 c o s 3 θ d 2 n 1 = 2 a 2 n 1 a n c o s θ d 2 n = a n 2
Thus, the Routh array has 2 n + 1 rows:
ω 2 n | d 2 n d 2 n 2 d 2 d 0 ω 2 n 1 | 2 n d 2 n ( 2 n 2 ) d 2 n 2 2 d 2 · | · · · · | · · · · | · · · ω 0 | d 0
Condition of the θ -stability, using the generalization Routh criterion, is no change of sequence sign:
d 2 n , 2 n d 2 n , , d 0
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]:
u s A u s B u s C = L s L μ 0 0 0 L s L μ 0 0 0 L s L μ d i s A / d t d i s B / d t d i s C / d t + R s 0 0 0 R s 0 0 0 R s i s A i s B i s C + e A e B e C
where L μ 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:
R = 2 R S , L = 2 ( L S L μ )
The transistors of voltage source inverter (Figure 2) are controlled by the pulse width modulation (PWM) method, with modulation factor d = 0 ÷ 1 , so U ( t ) = d U d .
Finally, the differential equations can be written:
U ( t ) = R I ( t ) + ψ e ω ( t ) + L d I ( t ) d t J d ω ( t ) d t = M e ( t ) M m ( t ) M e ( t ) = ψ e I ( t ) E ( t ) = ψ e ω ( t )
where E is back EMF of the two windings, J is moment of inertia, ψ e is permanent magnets flux, M e is electromagnetic torque, M m 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:
R = R T 0 [ 1 + α ( T b T 0 ) ] lub R = R T b N ( 1 + α Δ T )
where:
T 0 = 293 K⇒ 20 C,
T b —operating temperature,
T b N —rated operating motor temperature, depend on the insulation class,
α —temperature coefficient of resistance in Cu α = 4 · 10 3 1 K .
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:
-
T = L R , where T is electromagnetic time-constant,
-
B = J R ψ e 2 , 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 K α ). 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: M e = ψ e I . The system works in a closed-loop, so that the gains of measuring are also included:
  • Y is the gain of the current measuring,
  • K T is the gain of the angular speed measurement ,
  • K α 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 K p .
The PI current controller assumes the form:
G R I ( s ) = K R T R s + 1 T R s = K R T R s + K R T R s = m s + 1 V s
and the PI speed controller:
G R ω ( s ) = K ω T ω s + 1 T ω s
Thus, closed-loop transfer-function is equal:
G z ( s ) = K ω T ω s + 1 T ω s B K p m R s + B K p R B T V s 2 + B V + B K p m Y R s + B K p Y R + V k T 1 J s 1 + K T K ω T ω s + 1 T ω s B K p m R s + B K p R B T V s 2 + B V + B K p m Y R s + B K p Y R + V k T 1 J s
which leads to:
G z ( s ) = b 2 s 2 + b 1 s + b 0 a 4 s 4 + a 3 s 3 + a 2 s 2 + a 1 s + a 0
where:
b 2 = K ω T ω k T K p m 1 R b 1 = K ω T ω k T K p 1 R + K ω k T K p m 1 R b 0 = K ω k T K p 1 R a 4 = T ω J T V a 3 = T ω J V + T ω J K p m Y 1 R a 2 = T ω J K p Y 1 R + K ω T ω k T K t K p m 1 R + T ω V k E 2 R a 1 = K ω T ω k T K t K p 1 R + K ω k T K t K p m 1 R a 0 = K ω k T K t K p 1 R
The article considers the uncertainty of the following parameters: K p , L , R , ψ e . 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:
M ( s ) = a ̲ 4 , a ¯ 4 s 4 + a ̲ 3 , a ¯ 3 s 3 + a ̲ 2 , a ¯ 2 s 2 + a ̲ 1 , a ¯ 1 s + a ̲ 0 , a ¯ 0
where: a ̲ n , a ¯ n are the lower and upper ranges of coefficients (24) and a n > 0 . Their values are determined according to the following approach:
  • boundary coefficient a 4 at s 4 :
    a ̲ 4 = T ω J V L m i n R m a x
    a ¯ 4 = T ω J V L m a x R m i n
  • boundary coefficient a 3 at s 3 :
    a ̲ 3 = T ω J V + T ω J K p m i n m Y 1 R m a x
    a ¯ 3 = T ω J V + T ω J K p m a x m Y 1 R m i n
  • boundary coefficient a 2 at s 2 :
    a ̲ 2 = T ω J K p m i n Y 1 R m a x + K ω T ω k T m i n K t K p m i n m 1 R m a x + T ω V k E m i n 2 R m a x
    a ¯ 2 = T ω J K p m a x Y 1 R m i n + K ω T ω k T m a x K t K p m a x m 1 R m i n + T ω V k E m a x 2 R m i n
  • boundary coefficient a 1 at s:
    a ̲ 1 = K ω T ω k T m i n K t K p m i n 1 R m a x + K ω k T m i n K t K p m i n m 1 R m a x
    a ¯ 1 = K ω T ω k T m a x K t K p m a x 1 R m i n + K ω k T m a x K t K p m a x m 1 R m i n
  • boundary coefficient a 0 :
    a ̲ 0 = K ω k T m i n K t K p m i n 1 R m a x
    a ¯ 0 = K ω k T m a x K t K p m a x 1 R m i n

6. Uncertainly Parameters of BLDC Drive

The BLDC motor drive with the following rated parameters is considered:
  • R s = 5.25 [ Ω ]—single-phase resistance,
  • L s = 0.46 [mH]—single-phase inductance,
  • E = 2.62 [ V 1000 obr min ]—electromotive force,
  • k T = 24.9 [ mNm A ]—torque constant,
  • k E = 24.9 [ mNm A ]—emf constant,
  • J = 9.9 · 10 7 [kgm 2 ]—moment of inertia.
The ranges of the parameters uncertainty for the tests according following forms:
  • R [ 0.5 R N , 1.1 R N ] R [ 5.25 Ω , 11.55 Ω ]
  • L [ 0.7 L N , 1.1 L N ] L [ 0.644 mH , 1.012 mH ]
  • K p [ 0.95 K p N , 1.3 K p N ] K p [ 1.996 V V , 2.7313 V V ]
  • Ψ e [ 0.85 Ψ e N , Ψ e N ] Ψ e 0.0102 Vs rad , 0.012 Vs rad
Thus, the interval polynomial coefficients are:
  • a 4 2.7174 · 10 14 , 9.3943 · 10 14
  • a 3 4.9362 · 10 10 , 5.0621 · 10 10
  • a 2 8.1549 · 10 8 , 2.4617 · 10 7
  • a 1 7.1677 · 10 6 , 2.5386 · 10 5
  • a 0 3.5682 · 10 4 , 1.3 · 10 3
If the current controller (20) parameters are:
m = 0.0879 · 10 3 , V = 0.0123
and for velocity controller (21):
K ω = 8.9474 , T ω = 0.02
This rated closed-loop transfer-function (23) is:
G N ( s ) = 7.87 · 10 8 s 2 + 8.998 · 10 4 s + 0.04479 4.27 · 10 14 s 4 + 4.946 · 10 10 s 3 + 9.807 · 10 8 s 2 + 9.764 · 10 6 s + 4.861 · 10 4
The obtained transfer-function G N ( s ) , more precisely the characteristic polynomial, is subjected to algebraic studies of the relative θ -stability and the value of the relative damping factor is ξ = 0.5 , which corresponds to the angle θ = 30 . 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:
R [ 0.5 R N , 1.1 R N ]
Applying the minimum resistance value ( R m i n ) to the transfer-function coefficients (23) equation can be written:
G R m i n ( s ) = 1.574 · 10 7 s 2 + 1.8 · 10 3 s + 0.08958 8.54 · 10 14 s 4 + 5.019 · 10 10 s 3 + 1.96 · 10 7 s 2 + 1.95 · 10 5 s + 9.72 · 10 4
Analyzing the denominator of transfer-function (29), we get the relative damping factor ξ = 0.693 , corresponding to the angle θ = 43.87 . On the other hand, if we use the maximum value of resistance ( R m a x ) to the transfer-function coefficients (23), we obtain the transfer-function of a closed circuit in the form:
G R m a x ( s ) = 7.155 · 10 8 s 2 + 8.18 · 10 4 s + 0.04072 3.882 · 10 14 s 4 + 4.94 · 10 10 s 3 + 8.92 · 10 8 s 2 + 8.88 · 10 6 s + 4.42 · 10 4
where dumping factor ξ = 0.4524 and angle θ = 26.90 .
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 ( R m i n < R N ), because this sector is located inside the nominal sector θ N . When the resistance reaches the maximum value of ( R m a x ) , study object will lose relative stability because θ R m a x sector is located outside the nominal sector (all roots aren’t inside the θ R N sector). Thus, for R = R m a x the system has greater overshoot and oscillations during start-up than for rated parameters. On the other hand, for the R = R m i n 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 R m i n = 5.25 Ω , or ξ = 0.693 , has the lowest overshoot, and for a resistance value of R m a x = 11.55 Ω , or ξ = 0.4524 , 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:
L [ 0.7 L N , 1.1 L N ]
thus, transfer-function for L m i n and L m a x are in the following forms:
G L m i n ( s ) = 7.87 · 10 8 s 2 + 8.998 · 10 4 s + 0.04479 2.989 · 10 14 s 4 + 4.95 · 10 10 s 3 + 9.81 · 10 8 s 2 + 9.764 · 10 6 s + 4.86 · 10 4
G L m a x ( s ) = 7.87 · 10 8 s 2 + 8.998 · 10 4 s + 0.04479 4.697 · 10 14 s 4 + 4.95 · 10 10 s 3 + 9.81 · 10 8 s 2 + 9.764 · 10 6 s + 4.86 · 10 4
After analysing above transfer-functions, we get:
  • for (28): θ = 30 ,
  • for (31): θ = 29 . 95 ,
  • for (32): θ = 30 . 013 .
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 K p is determined by the voltage on capacitor in DC circuit:
K p [ 0.95 K p N , 1.3 K p N ]
where transfer-function (23) for ( K p m i n ) assumes form:
G K p m i n ( s ) = 7.477 · 10 8 s 2 + 8.548 · 10 4 s + 0.04255 4.27 · 10 14 s 4 + 4.94 · 10 10 s 3 + 9.39 · 10 8 s 2 + 9.276 · 10 6 s + 4.62 · 10 4
Analysis of the denominator (33) gives the relative damping factor ξ = 0.4785 corresponding to the angle θ = 28 . 59 .
Applying ( K p m a x ) substitution, transfer-funtcion (23) is in the following form:
G K p m a x ( s ) = 1.023 · 10 7 s 2 + 1.17 · 10 3 s + 0.05823 4.27 · 10 14 s 4 + 4.97 · 10 10 s 3 + 1.23 · 10 7 s 2 + 1.269 · 10 5 s + 6.32 · 10 4
where for ( K p m a x ) damping factor eqauls ξ = 0.6130 and angle θ = 37 . 81 .
Figure 10 shows relative stability sectors determined by the roots of characteristic polynomials for different values of the inverter gain K p .
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 ( K p m i n < K p N ), because this sector is located outside of the nominal sector. In second case the gain value is K p m a x and the system is relative stable because the sector θ K p m a x is within the nominal sector. The simulations (Figure 11) confirm the studies results. For K p m a x , the control system has the lowest overshoot compared to the speed waveforms at rated and minimum K p gain value, which is presented in Figure 11.

6.4. Interval of the Flux Linkage

The last parameter tested is flux linkage ψ e , which takes limit values in the range:
ψ e [ 0.85 Ψ e N , Ψ e N ]
and gets a transfer-function for ψ m i n :
G ψ m i n ( s ) = 6.69 · 10 8 s 2 + 7.648 · 10 4 s + 0.03807 4.27 · 10 14 s 4 + 4.95 · 10 10 s 3 + 9.387 · 10 8 s 2 + 8.3 · 10 6 s + 4.132 · 10 4
According to the analysis of the characteristic polynomial (35), the result is ξ = 0.4744 and θ = 28 . 32 . 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]:
P ( s ) = [ a ̲ n , a ¯ n ] s n + [ a ̲ n 1 , a ¯ n 1 ] s n 1 + + [ a ̲ 1 , a ¯ 1 ] s + [ a ̲ 0 , a ¯ 0 ]
where [ a ̲ i , a ¯ i ] are a range of uncertain polynomial coefficients and a ̲ n > 0 .
Theorem 1
(Kharitonov Theorem). Every polynomial P ( s ) from the polynomial family (36) is stable if and only if the following four Kharitonov polynomials are stable:
p 1 ( s ) = a ̲ n s n + a ̲ n 1 s n 1 + a ¯ n 2 s n 2 + a ¯ n 3 s n 3 + a ̲ n 4 s n 4 + a ̲ n 5 s n 5 +
p 2 ( s ) = a ¯ n s n + a ¯ n 1 s n 1 + a ̲ n 2 s n 2 + a ̲ n 3 s n 3 + a ¯ n 4 s n 4 + a ¯ n 5 s n 5 +
p 3 ( s ) = a ̲ n s n + a ¯ n 1 s n 1 + a ¯ n 2 s n 2 + a ̲ n 3 s n 3 + a ̲ n 4 s n 4 + a ¯ n 5 s n 5 +
p 4 ( s ) = a ¯ n s n + a ̲ n 1 s n 1 + a ̲ n 2 s n 2 + a ¯ n 3 s n 3 + a ¯ n 4 s n 4 + a ̲ n 5 s n 5 +
The stability conditions of the low degrees interval polynomials are given in [5] and for an uncertain fourth-degree polynomial
P 4 ( s ) = [ a ̲ 4 , a ¯ 4 ] s 4 + [ a ̲ 3 , a ¯ 3 ] s 3 + [ a ̲ 2 , a ¯ 2 ] s 2 + [ a ̲ 1 , a ¯ 1 ] s 1 + [ a ̲ 0 , a ¯ 0 ]
The interval polynomial is stable if the following inequalities occur:
a ̲ 4 > 0 , a ̲ 3 > 0 , a ̲ 2 > 0 , a ̲ 1 > 0 , a ̲ 0 > 0
and polynomials
p 2 ( s ) = a ¯ 4 s 4 + a ¯ 3 s 3 + a ̲ 2 s 2 + a ̲ 1 s 1 + a ¯ 0
p 4 ( s ) = a ¯ 4 s 4 + a ̲ 3 s 3 + a ̲ 2 s 2 + a ¯ 1 s 1 + a ¯ 0
are stable, so from Liénarda-Chiparta theorem [18] stability conditions are in the form
32 = a ¯ 3 a ̲ 2 a ̲ 1 a ¯ 3 2 a ¯ 0 a ¯ 4 a ̲ 1 2 > 0
34 = a ̲ 3 a ̲ 2 a ¯ 1 a ̲ 3 2 a ¯ 0 a ¯ 4 a ¯ 1 2 > 0
Thus, stability condition for considered interval polynomial are (42), (45) and (46).
After inserting the coefficients of interval polynomials a 0 ÷ a 4 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): θ = 0 , ξ = 0
  • for (38): θ = 1 . 83 , ξ = 0.0319
  • for (39): θ = 20 . 55 , ξ = 0.3510
  • for (40): θ = 13 . 64 , ξ = 0.2358
Comparing received sectors with a nominal sector θ G N (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 K p 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 2 R N additional resistances were used and a start-up was performed. However, for 0.5 R N , 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 K p . 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, K p 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 L = 1 / 0.7 L ˙ N = 1.31 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)
R = 0 , 5 R N , L = 0 , 5 L N , ψ = 0 , 8 ψ N , U = 15 V,
(2)
R = 2 R N , L = 1 , 2 L N , ψ = 1 , 2 ψ N , U = 10 V.
The step responses are shown in Figure 17.
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.

Author Contributions

Conceptualization, S.S., G.S., J.G.; methodology, S.S., G.S.; software, S.S., J.G.; validation, S.S., G.S.; formal analysis, S.S., G.S., J.G.; writing—review and editing, S.S., G.S., J.G.; visualization, S.S.; supervision, S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BLDCBrushless Direct Current

References

  1. Dorf, R.C.; Bishop, R.H. Modern Control Systems; Prentice-Hall: Upper Saddle River, NJ, USA, 2008. [Google Scholar]
  2. Skogestad, S.; Postlethwaite, I. Multivariable Feedback Control: Analysis and Design; Wiley: Chichester, UK, 2005. [Google Scholar]
  3. Doyle, J.C.; Francis, B.; Tannenbaum, A. Feedback Control Theory; Macmillan Publishing: New York, NY, USA, 1990. [Google Scholar]
  4. Bhattacharyya, S.P.; Chapellat, H.; Keel, L.H. Robust Control: The Parametric Approach; Prentice-Hall: Upper Saddle River, NJ, USA, 1995. [Google Scholar]
  5. Bialas, S. Resistance Stability of Polynomials and Matrices; Uczelniane Wydawnictwo Naukowo-Techniczne: Cracow, Poland, 2002. (In Polish) [Google Scholar]
  6. Stojic, M.R. Analytic tests for relative stability and strict aperiodicity. IEEE Trans. Educ. 1998, 41, 333–336. [Google Scholar] [CrossRef]
  7. Stojic, M.R.; Vukosavic, S.N. A generaliztion of Kharitonov’s four polynomial concept for robust relative-stability problems. Electron. Energ. 1993, 6, 1–12. [Google Scholar]
  8. Turowicz, A. Geometry of Polynomial Zeros; Wydawnictwo PWN: Warsaw, Poland, 1967. [Google Scholar]
  9. Sieklucki, G. Analysis of the transfer-function models of electric drives with voltage controlled source. Prz. Elektrotech. 2012, 88, 250–255. [Google Scholar]
  10. Stojic, M.R.; Siljak, D.D. Generalization of Hurwitz, Nyquist and Mikhailov stability criteria. IEEE Trans. Autom. Control. 1964, 10, 250–254. [Google Scholar] [CrossRef]
  11. Milivojevic, N.; Krishnamurthy, M.; Gurkaynak, Y.; Sathyan, A.; Lee, Y.; Emadi, A. Stability Analysis of FPGA-Based Control of Brushless DC Motors and Generators Using Digital PWM Technique. IEEE Trans. Ind. Electron. 2012, 59, 343–351. [Google Scholar]
  12. Chen, H.; Bai, H.; Zhang, Y.W.; Shi, Y.H. Stability and damping performance analysis of induction motor. In Proceedings of the 8th International Conference on Advances in Power System Control, Operation and Management (APSCOM 2009), Hong Kong, China, 8–11 November 2009; IET: London, UK, 2009; pp. 1–6. [Google Scholar]
  13. Rigatos, G.; Siano, P.; Marignetti, F.; Cecati, C. Nonlinear Optimal Control of Induction Motors for Maglev Trains Propulsion. In Proceedings of the 2019 IEEE International Conference on Mechatronics (ICM), Ilmenau, Germany, 18–20 March 2019; IEEE: Ilmenau, Germany, 2019. [Google Scholar]
  14. Reddy, N.T.; Basam, V.R.; Pasumarthi, M.R. A New Method for Multi-loop Control Tuning of High Power PM-BLDC Motor Drive. In Proceedings of the 2019 4th International Conference on Electrical, Electronics, Communication, Computer Technologies and Optimization Techniques (ICEECCOT), Mysuru, India, 13–14 December 2019; IEEE: Mysuru, India, 2019. [Google Scholar]
  15. Tashakori, A.; Ektesabi, M. Stability analysis of sensorless BLDC motor drive using digital PWM technique for electric vehicles. In Proceedings of the IECON 2012—38th Annual Conference on IEEE Industrial Electronics Society, Montreal, QC, Canada, 25–28 October 2012; IEEE: Montreal, QC, Canada, 2012. [Google Scholar]
  16. Hou, H.; Yao, W.; Zhang, W. Vector control of single phase brushless DC motor. In Proceedings of the 2016 19th International Conference on Electrical Machines and Systems (ICEMS), Chiba, Japan, 13–16 November 2016; IEEE: Chiba, Japan, 2016. [Google Scholar]
  17. Krishnan, R. Electric Motor Drives. Modelling, Analysis and Control Theory; Prentice-Hall: Upper Saddle River, NJ, USA, 2001. [Google Scholar]
  18. Ackermann, J.; Blue, P. Robust Control: The Parameter Space Approach; Springer: London, UK, 2002. [Google Scholar]
  19. Demirtas, M. Off-line tuning of a PI speed controller for a permanent magnet brushless DC motor using DSP. Energy Convers. Manag. 2011, 52, 264–273. [Google Scholar] [CrossRef]
Figure 1. θ -stability in term generalization of Routh theorem: (a) s-plane (b) ω -plane.
Figure 1. θ -stability in term generalization of Routh theorem: (a) s-plane (b) ω -plane.
Energies 14 00704 g001
Figure 2. Diagram of a BLDC motor supplies by a voltage source inverter.
Figure 2. Diagram of a BLDC motor supplies by a voltage source inverter.
Energies 14 00704 g002
Figure 3. The current (- - -) and electromotive forces (—) of the motor.
Figure 3. The current (- - -) and electromotive forces (—) of the motor.
Energies 14 00704 g003
Figure 4. Cascade control of BLDC drive.
Figure 4. Cascade control of BLDC drive.
Energies 14 00704 g004
Figure 5. Block diagram of the control system for a BLDC motor drive.
Figure 5. Block diagram of the control system for a BLDC motor drive.
Energies 14 00704 g005
Figure 6. Sector of relative θ -stability at rated drive parameters.
Figure 6. Sector of relative θ -stability at rated drive parameters.
Energies 14 00704 g006
Figure 7. Sectors of relative θ -stability for different stator resistance values.
Figure 7. Sectors of relative θ -stability for different stator resistance values.
Energies 14 00704 g007
Figure 8. Speed waveforms during start-up with disturbed stator resistance parameters.
Figure 8. Speed waveforms during start-up with disturbed stator resistance parameters.
Energies 14 00704 g008
Figure 9. Current waveforms during start-up with disturbed stator resistance parameters.
Figure 9. Current waveforms during start-up with disturbed stator resistance parameters.
Energies 14 00704 g009
Figure 10. Sectors of relative θ -stability for different gain inverter values K p .
Figure 10. Sectors of relative θ -stability for different gain inverter values K p .
Energies 14 00704 g010
Figure 11. Speed waveforms during start-up with disturbed gain inverter parameters.
Figure 11. Speed waveforms during start-up with disturbed gain inverter parameters.
Energies 14 00704 g011
Figure 12. Sectors of relative θ -stability for different flux values.
Figure 12. Sectors of relative θ -stability for different flux values.
Energies 14 00704 g012
Figure 13. Speed waveforms during start-up with disturbed flux parameters.
Figure 13. Speed waveforms during start-up with disturbed flux parameters.
Energies 14 00704 g013
Figure 14. Step responses for different stator resistances.
Figure 14. Step responses for different stator resistances.
Energies 14 00704 g014
Figure 15. Step responses for different voltage in the DC circuit— K p interval.
Figure 15. Step responses for different voltage in the DC circuit— K p interval.
Energies 14 00704 g015
Figure 16. Step responses for too large inductance in the optimization of the controllers.
Figure 16. Step responses for too large inductance in the optimization of the controllers.
Energies 14 00704 g016
Figure 17. Step responses for many mistakes in identification.
Figure 17. Step responses for many mistakes in identification.
Energies 14 00704 g017
Figure 18. Sectors of relative θ -stability for Kharitonov’s theorem.
Figure 18. Sectors of relative θ -stability for Kharitonov’s theorem.
Energies 14 00704 g018
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sobieraj, S.; Sieklucki, G.; Gromba, J. Relative Stability of Electrical into Mechanical Conversion with BLDC Motor-Cascade Control. Energies 2021, 14, 704. https://doi.org/10.3390/en14030704

AMA Style

Sobieraj S, Sieklucki G, Gromba J. Relative Stability of Electrical into Mechanical Conversion with BLDC Motor-Cascade Control. Energies. 2021; 14(3):704. https://doi.org/10.3390/en14030704

Chicago/Turabian Style

Sobieraj, Sylwester, Grzegorz Sieklucki, and Józef Gromba. 2021. "Relative Stability of Electrical into Mechanical Conversion with BLDC Motor-Cascade Control" Energies 14, no. 3: 704. https://doi.org/10.3390/en14030704

APA Style

Sobieraj, S., Sieklucki, G., & Gromba, J. (2021). Relative Stability of Electrical into Mechanical Conversion with BLDC Motor-Cascade Control. Energies, 14(3), 704. https://doi.org/10.3390/en14030704

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