2.1. The basics of Simulink model
Dynamics of multiphase electric drive is described by nonlinear differential equations with time varying coefficients. Mathematical description of three phase motor in a stationary three phase reference frame due to varying mutual inductances because rotation is seldom used for modelling and analysis. An equivalent circuit of symmetrical six-phase machine with two winding sets together with mathematical model of machine is given in [
24]. Stator of equivalent circuit has two direct current supply voltages
and
as well as
and
in reference frame with
and
axis. Rotor supply voltages
and
are accepted equal to zero for cage type rotor. The d-q mathematical model in various forms is used for analysis of induction motor dynamics. Rotor reference frame was employed in [
5], [
21,
22]. The stationary reference frame, assuming rotor variables transformed to fictitious stator windings also is widely used in analysis of electrical machines. Stationary reference frame is applied in [
16,
18],[
23]. The rotating with synchronous speed reference frame is applied in[
15], [
24,
25,
26] to model six-phase controlled induction drive. The Simulink models elaborated in any reference frame allow analyse control methods and transients in the systems of multiphase drives [
15,
17,
24].
The Simulink libraries do not include six-phase voltage converters. Elaborated six-phase voltage system in model, shown in
Figure 1, is made of two typical three phase generators. Two generators allow getting six-phase voltage system from two three phase voltage sets:
and
. The voltage set
is shifted by 60 electrical degrees with respect to voltage set
. Two sets of three phase voltages are transformed to two phase stationary reference frames using Clarke’s transform for six-phase voltage applied for each three-phase voltage set separately.
The mathematical model of the motor in synchronously rotating reference frame is expressed as voltage drops
and
as well as
and
across two sets of stator windings and single common rotor winding voltages
and
, referred to a stator as this:
where flux linkages aligned with the direct and quadrature are calculated from equivalent circuit [
27] which assumes presence of two similar three phase windings sets in the stator with their leakage inductance
and mutual leakage inductance between two stator winding sets
. Stator of equivalent circuit has two direct current supply voltages
and
in reference frame with
and
axis. Rotor is supplied by voltages
and
which are assumed equal to zero for cage type rotor. The other notations as
denotes stator resistor,
is rotor resistance referred to stator,
is magnetizing inductance and
is mutual leakage inductance between two stator windings,
is stator leakage inductance,
is rotor leakage inductance,
and
are stator
axis current components,
and
are stator
axis current components,
and
are rotor current components in
and
axis accordingly are calculated:
where
,
are stator
axis flux linkage,
,
, are stator
axis flux linkage,
and
are rotor
and
axis flux linkages accordingly.
The electromagnetic torque, delivered by motor is calculated as:
where P is number of poles.
The movement of electric drive is described by equation:
where
is motor rotational velocity,
is rotor inertia,
is electromagnetic torque and
is load torque.
Analysis of multiphase motor dynamics applies all general assumption for electrical machines. Three phase machine transformation from the stationary three phase reference frame
to two-phase stationary reference frame
and synchronously rotating reference frame
is described in [
29], where the stator voltage
axis is aligned with
axis. The transformations are made in this way. If stationary
axes of three phase machine are oriented at angle
, the voltages
and
can be represented in matrix form as:
where
is zero component. In balanced machine
.
Usually, is assumed equal to zero.
Inverse transform of (5) can be made in this way:
The synchronously rotating
axes rotate with synchronous speed
with respect to
axes and form the angle
. The voltages on
axes are transformed into
axes by this way:
According to reference frame transform theory, set of three phase voltages
,
, and
shifted by 120
0 apart, starting with phase voltage
lagging the voltage
by 60
0, could be transformed to stationary reference frame in this way:
Equation (8) can be written as:
Then the inverse transform becomes:
According to Eq.9, sinusoidal variables in reference frame appear as sinusoidal quantities in rectangular reference frame and .
Voltages
and
can be transformed to synchro-nously rotating reference frame like this:
Thus, we considered two transforms for two three phase voltage sets
and
. The voltages are shifted by 120 electrical degrees in each set. These sets are transformed to stationary rectangular reference frame with equations (5) and (8). The variables from rectangular stationary reference frame are transformed to synchronously rotating reference frame using equations (7) and (11). Simulink model of six phase motor elaborated on the (1), (2), (3) and (4) and is given [
28].
The block diagram in
Figure 1 shows elaborated Simulink model of six-phase drive with indirect field-oriented control and PID and PID fuzzy controllers.
Reference speed is compared with motor speed (
Figure 1) in PID or PID fuzzy controller. The controller elaborates reference torque signal for two field-oriented blocks. Block “Field oriented control 1” produces gating pulses for three phase converter tuned to generate three phase voltage set V
A, V
B and V
C. The next block transforms three phase voltages to two phase voltages in stationary reference frame and the last conversion transforms two phase stationary to rotational reference frame, and forms two voltages v
ds1 and v
qs1. Block “Field oriented control 2” elaborates gating pulses to three phase inverter, tuned to produce voltage set V
D, V
E and V
F. Voltage V
D lags voltage V
A by 60 electrical degrees. This voltage set is transformed to two phase stationary reference frame and afterwords to two phase rotating reference frame and gives two voltage v
ds2 and v
qs2 entering as input signals to six-phase motor model. Motor currents i
ds1, i
qs1 and i
ds2, i
qs2 in two phase rotational reference frame are output of motor model. Each pair of currents are converted to two phase stationary reference frame i
ds1, i
qs1 with Eq. 6, where voltages should be replaced by currents and for i
ds2, i
qs2: Eq. 9; afterwords they are transformed to two sets of three phase currents I
A, I
B, I
C and I
D, I
E,I
F, used as feedback for two field oriented control blocks, elaborating gating pulses according to produced torque reference signal and calculated motor currents.
Gating pulses change voltages and finally motor speed in the control system follows the speed reference.
Fuzzy controller is one of intelligent controllers, suitable for systems even having no mathematical model. Fuzzy control is nonlinear control. The six-phase motor d – q model in the rotating reference frame based on the equivalent motor circuit with two stator windings and a single rotor winding is a non-linear complex model, which can have variable parameters and load disturbance. The fuzzy controller can control similar systems. Operation of fuzzy logic controller is based on sets. Fuzzy logic controller consists of fuzzification block, the knowledge base, and the inference engine. Fuzzification block of fuzzy controller transforms crisp inputs into fuzzy sets [
19]. Seven linguistic variables related to the membership functions, are used to describe the fuzzy sets: PB (positive big), PM (positive medium), PS (positive small), ZE (zero), NS (negative small), NM (negative medium), NB (negative big). The designed system controller is based on the error, and the integral error and error derivative as inputs. Each set represents several linguistic variables defining the possible state of the output. The knowledge base consists of the membership functions and the fuzzy rules, which are obtained by knowledge of the system operation according to the environment for example, provided speed reference or application of load. Knowledge base defines the rules represented by IF-THEN rule base, relating the input and output variables. Designed controller comprises 49 rules.
The inference engine, which is based upon the input fuzzy sets, uses the IF-THEN rules in the knowledge base to make the decisions.
Knowledge base relates 49 IF-THEN rules and membership functions to produce the fuzzy output. The defuzzification process converts the fuzzy output variables to crisp output. Designed fuzzy controller was based on per unit membership functions.
2.2. Model of indirect rotor flux orientation
According to [
30] in indirect field-oriented control, rotor flux vector is aligned with the axis
d of the revolving
reference frame. Two reference frames stationary
and revolving at synchronous speed
are shifted by the angle
between two reference frames.
Angle between two reference frames
varies with rotor rotation. Rotor flux
is aligned with the axis d and its projection to the d axis
and
. Then torque, delivered by an induction motor is calculated as:
where p is number of poles and
is torque producing current,
is magnetizing inductance and
is rotor inductance, referred to stator.
Current
is calculated in rotating reference frame from chosen speed reference
and current
from given flux reference
, rotating at synchronous speed
.
The angular position of
is reference frame orientation angle and is calculated as:
where
is measured mechanical speed of the motor,
is rotor time constant,
reference current, calculated from reference flux
,
is slip speed.
Vector control Simulink model of six-phase motor enters the Simulink model of total six-phase drive shown in
Figure 1.
The small power motor parameters were measured experimentally performing locked rotor and no-load tests. They are given in
Table 1.