Nonlinear Digital Simulation Models of Switched Reluctance Motor Drive

: The paper deals with nonlinear simulation models of a drive consisting of the four-phase 8 / 6 doubly salient switched reluctance motor (SRM), the four-phase dissymmetric bridge power converter and the closed-cycle rotor speed control strategy carried out by the pulse width modulation (PWM) with variable angle and combined control scheme with the PI algorithm. All presented considerations are based on a MATLAB-SIMULINK platform. The nonlinear mathematical model of the analyzed SRM drive was obtained as a combination of the two dimensional (2D) ﬁnite element model (FEM) of the motor and the nonlinear model of the electrical network of the power supply circuit. The main model and its seven sub-modules, such as the controller module, one phase simulation module, rotor position angle transformation module, power system module, phase current operation module, “subsystem” module, and electromagnetic torque of one phase operation module, are described. MATLAB functions store the magnetization curves data of the motor obtained by the 2D FEM electromagnetic ﬁeld calculations, as well as the data of magnetic co-energy curves of the motor calculated from the magnetization curves. The 2D specimen insert method is adopted in MATLAB functions for operating the ﬂux linkage and the magnetic co-energy at the given phase current and rotor position. The phase current waveforms obtained during simulations match with the tested experimentally phases current waveforms at the same rotor speed and the same load basically. The simulated rotor speed curves also agree with the experimental rotor speed curves. This means that the method of suggested nonlinear simulation models of the analyzed SRM drive is correct, and the model is accurate.


Introduction
The switched reluctance motor (SRM) drive, described comprehensively in [1], has been developed for mining equipment [2], electric vehicles [3], high speed equipment [4], different generators [5,6], wind generators [7], high speed and high power applications [8,9], flywheel energy storage applications [10], linear transportation [11], automotive applications [12], and so on, due to favorable conditions in the design of the motor and power supply, ease of four quadrants operation, high start torque with low start current, high efficiency within vast rotor speed ranges, high dependability with the independence of magnetic paths for each phase and the independence of circuits in each phase. The performance of switched reluctance motors is better than that of the hysteresis direct torque control of permanent-magnet synchronous motor system [13]. The SRM drive generally requires closed-cycle rotor speed control, where the opening angle and the turn-off angle of the power similar to the one proposed in this paper for the improvement of the low-speed performance of SRMs is presented.
The main goal of the paper was to develop the nonlinear simulation models of the advanced SRM drive using the MATLAB platform combined with the 2D electromagnetic field FE analysis and SIMULINK software. The nonlinear electrical meshwork model of the 4-phase asymmetric bridge SRM power supply circuit and the closed-cycle rotor speed control strategy have been realized by the variable angle pulse width modulation control scheme together with the PI algorithm. Finally, the simulation solution of phase current waveforms and rotor speed curves were compared with the measured waveforms for a real four-phase 8/6 poles SRM motor.
This developed procedure is similar to the algorithm given in [36], where the multi-objective finite element model (FEM) optimization of a PMSM motor was performed with the help of MATLAB and ANSYS-Maxwell scripting capabilities. This result was used after that in paper [37] for the synthesis of the dead-beat current controller for this motor.

System Components and Mathematical Model
The analyzed SRM drive prototype consists of a four-phase 8/6 poles SRM and its four-phase asymmetric bridge power supply. The cross-section of the SRM under consideration has been shown in Figure 1. One phase winding of the motor is composed by two coils placed on the opposite stator poles. The motor windings are supplied by the current from the bridge power supply.
Energies 2020, 13, x FOR PEER REVIEW 3 of 17 procedure similar to the one proposed in this paper for the improvement of the low-speed performance of SRMs is presented. The main goal of the paper was to develop the nonlinear simulation models of the advanced SRM drive using the MATLAB platform combined with the 2D electromagnetic field FE analysis and SIMULINK software. The nonlinear electrical meshwork model of the 4-phase asymmetric bridge SRM power supply circuit and the closed-cycle rotor speed control strategy have been realized by the variable angle pulse width modulation control scheme together with the PI algorithm. Finally, the simulation solution of phase current waveforms and rotor speed curves were compared with the measured waveforms for a real four-phase 8/6 poles SRM motor.
This developed procedure is similar to the algorithm given in [36], where the multi-objective finite element model (FEM) optimization of a PMSM motor was performed with the help of MATLAB and ANSYS-Maxwell scripting capabilities. This result was used after that in paper [37] for the synthesis of the dead-beat current controller for this motor.

System Components and Mathematical Model
The analyzed SRM drive prototype consists of a four-phase 8/6 poles SRM and its four-phase asymmetric bridge power supply. The cross-section of the SRM under consideration has been shown in Figure 1. One phase winding of the motor is composed by two coils placed on the opposite stator poles. The motor windings are supplied by the current from the bridge power supply. The topology of power circuit is shown in Figure 2. It consists of two major switches and two diodes in every phase of the power supply system.  The topology of power circuit is shown in Figure 2. It consists of two major switches and two diodes in every phase of the power supply system. procedure similar to the one proposed in this paper for the improvement of the low-speed performance of SRMs is presented. The main goal of the paper was to develop the nonlinear simulation models of the advanced SRM drive using the MATLAB platform combined with the 2D electromagnetic field FE analysis and SIMULINK software. The nonlinear electrical meshwork model of the 4-phase asymmetric bridge SRM power supply circuit and the closed-cycle rotor speed control strategy have been realized by the variable angle pulse width modulation control scheme together with the PI algorithm. Finally, the simulation solution of phase current waveforms and rotor speed curves were compared with the measured waveforms for a real four-phase 8/6 poles SRM motor.
This developed procedure is similar to the algorithm given in [36], where the multi-objective finite element model (FEM) optimization of a PMSM motor was performed with the help of MATLAB and ANSYS-Maxwell scripting capabilities. This result was used after that in paper [37] for the synthesis of the dead-beat current controller for this motor.

System Components and Mathematical Model
The analyzed SRM drive prototype consists of a four-phase 8/6 poles SRM and its four-phase asymmetric bridge power supply. The cross-section of the SRM under consideration has been shown in Figure 1. One phase winding of the motor is composed by two coils placed on the opposite stator poles. The motor windings are supplied by the current from the bridge power supply. The topology of power circuit is shown in Figure 2. It consists of two major switches and two diodes in every phase of the power supply system.   The alterable angle PWM combined control scheme with the PI algorithm is used for the closed-cycle rotor speed drive control. The opening angle of the major switches, θ 1 , and the turn-off angle of the main switches, θ 2 , at the two rotor speed values are given in Table 1, where (θ 1 and θ 2 = 0 are defined as the rotor position while the axis of the rotor slot is in accordance with that of the stator pole of the conducting phase). Table 1. θ 1 and θ 2 at the different rotor speed ranges.
The developed mathematical model of this SRM drive is in fact the nonlinear model of the electrical circuit of the power supply system combined with the two 2D FE motor models. Neglecting mutual inductances, it becomes [134]: where: where: N-each phase turn numbers of winding, L-active length of motor iron core, S-stator pole area at one side, ∆ i and ∆ j -area of the split cell at the right and left side of the stator poles, A i and A j -magnetomotive force vector of the split cell at the right and left side of the stator poles, separately.
The typical equation for the magnetic vector potential (together with the homogeneous boundary conditions) describes the magnetic field distribution within the SRM: where: Energies 2020, 13, 6715 5 of 17 γ-the reciprocal of permeability, J-stator phase winding current density, D 2 -outer diameter of stator, d 2 -diameter of rotor bore, A-magnetic vector potential. where: γ-the reciprocal of permeability, J-stator phase winding current density, D2-outer diameter of stator, d2-diameter of rotor bore, A-magnetic vector potential.
where US is the DC supply voltage of the power supply, and UT is the on-state voltage drop of the major switch.
(2) During the course of commutation, the major switches in one phase circuit are closed 2 , , , , , where, UD is the diode on-state voltage drop.
(3) While the power supply does not work, the current in phase is naturally continuous and one phase major switch is turned off, , , , , , The system equation of the mechanical motion can be written as [1,12,35]: where J is the sum of inertia moments of the motor and loads, Te is the electromagnetic torque of the motor, TL is the torque of the loads, H is the viscous coefficient of friction, t is the time period and θ is the rotor position.
where U S is the DC supply voltage of the power supply, and U T is the on-state voltage drop of the major switch.
(2) During the course of commutation, the major switches in one phase circuit are closed where, U D is the diode on-state voltage drop.
(3) While the power supply does not work, the current in phase is naturally continuous and one phase major switch is turned off, The system equation of the mechanical motion can be written as [1,12,35]: where J is the sum of inertia moments of the motor and loads, T e is the electromagnetic torque of the motor, T L is the torque of the loads, H is the viscous coefficient of friction, t is the time period and θ is the rotor position.

Simulation Models
The nonlinear simulation model of the SRM drive developed by MATLAB-SIMULINK platform is given in Figure 4. In this model, the turn-on angle of the major switches in the power converter, θ 1 , is fixed at "Xon", the turn-off angle of the major switches in the power supply system, θ 2 , is fixed at "Xoff ", which is shown for different rotor speeds ranges in Table 1. "Us" is the DC supply voltage of the power supply, "θ" is the rotor position, "ω" is the angular speed of the motor, "n*" is the set speed, "n" is the practical rotor speed, "D" is the duty ratio of the pulse width modulation signal, "n" is the given rotor speed, "Te" is the total electromagnetic torque of the motor, "Ta", "Tb", "Tc", and "Td" are the electromagnetic torque in A phase, B phase, C phase, D phase, separately, "TL" is the load, "J" is the sum of all system inertia moments, "H" is the viscous coefficient of friction, "ia", "ib", "ic", and "id" are the phase current in A phase, B phase, C phase, D phase, separately. The "PI controller" module is responsible for the closed-cycle rotor speed control. The expanding module of "PI controller" is shown in Figure 5, where, "Kp" is the proportion coefficient, "K I " is the integral coefficient, the input "in1" is the difference among the given rotor speed and the real rotor speed, the output "out1" is the duty ratio of the pulse width modulation semaphore, and the "saturation" is used to limit the duty ratio between 0.000 and 1.000.

Simulation Models
The nonlinear simulation model of the SRM drive developed by MATLAB-SIMULINK platform is given in Figure 4. In this model, the turn-on angle of the major switches in the power converter, θ1, is fixed at "Xon", the turn-off angle of the major switches in the power supply system, θ2, is fixed at "Xoff", which is shown for different rotor speeds ranges in Table 1. "Us" is the DC supply voltage of the power supply, "θ" is the rotor position, "ω" is the angular speed of the motor, "n*" is the set speed, "n" is the practical rotor speed, "D" is the duty ratio of the pulse width modulation signal, "n" is the given rotor speed, "Te" is the total electromagnetic torque of the motor, "Ta", "Tb", "Tc", and "Td" are the electromagnetic torque in A phase, B phase, C phase, D phase, separately, "TL" is the load, "J" is the sum of all system inertia moments, "H" is the viscous coefficient of friction, "ia", "ib", "ic", and "id" are the phase current in A phase, B phase, C phase, D phase, separately. The "PI controller" module is responsible for the closed-cycle rotor speed control. The expanding module of "PI controller" is shown in Figure 5, where, "Kp" is the proportion coefficient, "KI" is the integral coefficient, the input "in1" is the difference among the given rotor speed and the real rotor speed, the output "out1" is the duty ratio of the pulse width modulation semaphore, and the "saturation" is used to limit the duty ratio between 0.000 and 1.000.  One phase simulation module of the main model is shown in Figure 6. It contains the power system module "converter", the phase current computation module "calculating i", the one phase electromagnetic torque calculation module "calculating T", the rotor place angle conversion module "θ conversion". The phase difference of the analyzed SRM is 15°. The rotor position angle conversion module is used to convert the real rotor position of every phase as the relative rotor position. The expanding rotor position angle conversion module is indicated in Figure 7. The input "ln1" is the absolute rotor place angle, and the output "out1" is the relative rotor position of the certain phase. "Constant" is one rotor period, θr, such as one rotor period of four-phase 8/6 is 60°. "Constant1" is 0° for A phase, 15° for B phase, 30° for C phase, and 45° for D phase.   One phase simulation module of the main model is shown in Figure 6. It contains the power system module "converter", the phase current computation module "calculating i", the one phase electromagnetic torque calculation module "calculating T", the rotor place angle conversion module "θ conversion". The phase difference of the analyzed SRM is 15 • . The rotor position angle conversion module is used to convert the real rotor position of every phase as the relative rotor position. The expanding rotor position angle conversion module is indicated in Figure 7. The input "ln1" is the absolute rotor place angle, and the output "out1" is the relative rotor position of the certain phase. "Constant" is one rotor period, θ r , such as one rotor period of four-phase 8/6 is 60 • . "Constant1" is 0 • for A phase, 15 • for B phase, 30 • for C phase, and 45 • for D phase. One phase simulation module of the main model is shown in Figure 6. It contains the power system module "converter", the phase current computation module "calculating i", the one phase electromagnetic torque calculation module "calculating T", the rotor place angle conversion module "θ conversion". The phase difference of the analyzed SRM is 15°. The rotor position angle conversion module is used to convert the real rotor position of every phase as the relative rotor position. The expanding rotor position angle conversion module is indicated in Figure 7. The input "ln1" is the absolute rotor place angle, and the output "out1" is the relative rotor position of the certain phase. "Constant" is one rotor period, θr, such as one rotor period of four-phase 8/6 is 60°. "Constant1" is 0° for A phase, 15° for B phase, 30° for C phase, and 45° for D phase.   One phase simulation module of the main model is shown in Figure 6. It contains the power system module "converter", the phase current computation module "calculating i", the one phase electromagnetic torque calculation module "calculating T", the rotor place angle conversion module "θ conversion". The phase difference of the analyzed SRM is 15°. The rotor position angle conversion module is used to convert the real rotor position of every phase as the relative rotor position. The expanding rotor position angle conversion module is indicated in Figure 7. The input "ln1" is the absolute rotor place angle, and the output "out1" is the relative rotor position of the certain phase. "Constant" is one rotor period, θr, such as one rotor period of four-phase 8/6 is 60°. "Constant1" is 0° for A phase, 15° for B phase, 30° for C phase, and 45° for D phase.    angles of the major switches, a pulse width modulation duty ratio, position angle of the machine rotor, phase currents respectively, and an output "out1", which is the phase voltage "U k ". Proposed "MATLAB Function" is used to calculate the phase voltages, and contains switches "turn-on" and "turn-off" rule. This rule is based on the PWM duty ratio. The "MATLAB Function" includes also opening angle and major switches turn-off angle. While θ 1 ≤ θ ≤ θ 2 , if the PWM semaphore is "1", the phase voltage is expressed by Equation (4), if the PWM semaphore is "0", the phase voltage is denoted as Equation (6). While θ > θ 2 , the phase voltage is expressed by Equation (5) if the phase current exists, and the phase voltage is zero if the phase current does not exist. While θ < θ 1 , the voltage is zero if the phase current does not exist, and the voltage is expressed by Equation (5) if the current in phase exists.
Energies 2020, 13, x FOR PEER REVIEW 8 of 17 Figure 8 presents the module of the expanding power converter, where the import of module "ln1"~"ln6", are the DC supply voltages, respectively. This module contains also turn-on and turnoff angles of the major switches, a pulse width modulation duty ratio, position angle of the machine rotor, phase currents respectively, and an output "out1", which is the phase voltage "Uk". Proposed "MATLAB Function" is used to calculate the phase voltages, and contains switches "turn-on" and "turn-off" rule. This rule is based on the PWM duty ratio. The "MATLAB Function" includes also opening angle and major switches turn-off angle. While θ1 ≤ θ ≤ θ2, if the PWM semaphore is "1", the phase voltage is expressed by Equation (4), if the PWM semaphore is "0", the phase voltage is denoted as Equation (6). While θ > θ2, the phase voltage is expressed by Equation (5) if the phase current exists, and the phase voltage is zero if the phase current does not exist. While θ < θ1, the voltage is zero if the phase current does not exist, and the voltage is expressed by Equation (5) if the current in phase exists. From Equation (1), we have: Based on the Equation (9), the simulation model of the phase current calculation module is shown in Figure 9. The import of module, "ln1"~"ln3", are the rotor angular velocity, the phase voltage, the converted relative rotor place, respectively, and the export of module "out1" is the instantaneous phase current. "Subsystem" is the module for calculating ∂Ψ/∂i and ∂Ψ/∂θ, which are gained by the magnetization curves of the motor.  From Equation (1), we have: Based on the Equation (9), the simulation model of the phase current calculation module is shown in Figure 9. The import of module, "ln1"~"ln3", are the rotor angular velocity, the phase voltage, the converted relative rotor place, respectively, and the export of module "out1" is the instantaneous phase current. "Subsystem" is the module for calculating ∂Ψ/∂i and ∂Ψ/∂θ, which are gained by the magnetization curves of the motor.
Energies 2020, 13, x FOR PEER REVIEW 8 of 17 Figure 8 presents the module of the expanding power converter, where the import of module "ln1"~"ln6", are the DC supply voltages, respectively. This module contains also turn-on and turnoff angles of the major switches, a pulse width modulation duty ratio, position angle of the machine rotor, phase currents respectively, and an output "out1", which is the phase voltage "Uk". Proposed "MATLAB Function" is used to calculate the phase voltages, and contains switches "turn-on" and "turn-off" rule. This rule is based on the PWM duty ratio. The "MATLAB Function" includes also opening angle and major switches turn-off angle. While θ1 ≤ θ ≤ θ2, if the PWM semaphore is "1", the phase voltage is expressed by Equation (4), if the PWM semaphore is "0", the phase voltage is denoted as Equation (6). While θ > θ2, the phase voltage is expressed by Equation (5) if the phase current exists, and the phase voltage is zero if the phase current does not exist. While θ < θ1, the voltage is zero if the phase current does not exist, and the voltage is expressed by Equation (5) if the current in phase exists. From Equation (1), we have: Based on the Equation (9), the simulation model of the phase current calculation module is shown in Figure 9. The import of module, "ln1"~"ln3", are the rotor angular velocity, the phase voltage, the converted relative rotor place, respectively, and the export of module "out1" is the instantaneous phase current. "Subsystem" is the module for calculating ∂Ψ/∂i and ∂Ψ/∂θ, which are gained by the magnetization curves of the motor.  Finally, we have: Based on Equation (10), the expanding "Subsystem" module is shown in Figure 10.
In the module, "in1" is the instantaneous current value, "in2" is the instantaneous converted relative rotor position value, "out1" is ∂Ψ/∂i, and "out2" is ∂Ψ/∂θ. In the "Subsystem" the current incremental quantity ∆i is calculated from "Constant1" and "Constant2", but the rotor place incremental quantity ∆θ is calculated from "Constant3" and "Constant4". Figure 11 shows magnetization map of analyzed motor, computed by the electromagnetic 2D FEM. These curves are stored in MATLAB Functions by data, such as "MATLAB Fcn1", "MATLAB Fcn2", "MATLAB Fcn3" presented in Figure 10. In order to compute the flux linkage at the given rotor position and phase current, in view of the magnetization curves, the "MATLAB Function" with the 2-D specimen insert method has been adopted.
Based on Equation (10), the expanding "Subsystem" module is shown in Figure 10.
In the module, "in1" is the instantaneous current value, "in2" is the instantaneous converted relative rotor position value, "out1" is ∂Ψ/∂i, and "out2" is ∂Ψ/∂θ. In the "Subsystem" the current incremental quantity Δi is calculated from "Constant1" and "Constant2", but the rotor place incremental quantity Δθ is calculated from "Constant3" and "Constant4". Figure 11 shows magnetization map of analyzed motor, computed by the electromagnetic 2D FEM. These curves are stored in MATLAB Functions by data, such as "MATLAB Fcn1", "MATLAB Fcn2", "MATLAB Fcn3" presented in Figure 10. In order to compute the flux linkage at the given rotor position and phase current, in view of the magnetization curves, the "MATLAB Function" with the 2-D specimen insert method has been adopted.  According to Equations (9) and (10) and Figures 6-10, the 2D FE model of the electric machine and the nonlinear electrical meshwork model of the power system have been integrated. The magnetic co-energy of the motor can be computed as follows: The magnetic co-energy curves of the motor are calculated by Equation (11) based on the magnetization curves in Figure 11. They are shown in Figure 12. Finally, we have: Based on Equation (10), the expanding "Subsystem" module is shown in Figure 10.
In the module, "in1" is the instantaneous current value, "in2" is the instantaneous converted relative rotor position value, "out1" is ∂Ψ/∂i, and "out2" is ∂Ψ/∂θ. In the "Subsystem" the current incremental quantity Δi is calculated from "Constant1" and "Constant2", but the rotor place incremental quantity Δθ is calculated from "Constant3" and "Constant4". Figure 11 shows magnetization map of analyzed motor, computed by the electromagnetic 2D FEM. These curves are stored in MATLAB Functions by data, such as "MATLAB Fcn1", "MATLAB Fcn2", "MATLAB Fcn3" presented in Figure 10. In order to compute the flux linkage at the given rotor position and phase current, in view of the magnetization curves, the "MATLAB Function" with the 2-D specimen insert method has been adopted.  According to Equations (9) and (10) and Figures 6-10, the 2D FE model of the electric machine and the nonlinear electrical meshwork model of the power system have been integrated. The magnetic co-energy of the motor can be computed as follows: The magnetic co-energy curves of the motor are calculated by Equation (11) based on the magnetization curves in Figure 11. They are shown in Figure 12. According to Equations (9) and (10) and Figures 6-10, the 2D FE model of the electric machine and the nonlinear electrical meshwork model of the power system have been integrated. The magnetic co-energy of the motor can be computed as follows: The magnetic co-energy curves of the motor are calculated by Equation (11) based on the magnetization curves in Figure 11. They are shown in Figure 12. The equation for calculating electromagnetic torque generated by each phase of the double salient reluctance machine can be written as Based on the Equation (12), the simulation model of one phase electromagnetic torque computation module is shown in Figure 13. In the module, "ln1" is the instantaneous current value, "in2" is the instantaneous converted relative rotor place value, and the export of module "out1" is the one phase electromagnetic torque instantaneous value. "Constant1" and "Constant2" are the rotor position increment, Δθ.
The magnetic co-energy curves data of the motor in Figure 12 is stored in MATLAB Functions by data, such as "MATLAB Fcn1", "MATLAB Fcn2" in Figure 13. During calculating for computing magnetic co-energy at the given rotor place and phase current, in view of the magnetic co-energy curves in MATLAB Functions the 2D specimen insert method has been adopted.

Simulated and Measured Results
The developed SRM (see Appendix A) drive prototype with the closed-cycle control speed conducted with use of the variable angle of the PWM combined control scheme with the PI algorithm is simulated by the proposed nonlinear simulation models on MATLAB-SIMULINK, and it is also tested experimentally at the same conditions. The photography of the SRM with the torque/speed meter and the load is shown in Figure 14, where the motor is on the left, the torque/speed meter is in the middle, and the load is on the right. The picture of the 4-phase asymmetric bridge power unit with the Intel 8XC196KC MCU controller (high performance member of the MCS 96 microcontroller family) is indicated in Figure 15, where on the right the power converter is shown, the Intel controller The equation for calculating electromagnetic torque generated by each phase of the double salient reluctance machine can be written as Based on the Equation (12), the simulation model of one phase electromagnetic torque computation module is shown in Figure 13. The equation for calculating electromagnetic torque generated by each phase of the double salient reluctance machine can be written as Based on the Equation (12), the simulation model of one phase electromagnetic torque computation module is shown in Figure 13. In the module, "ln1" is the instantaneous current value, "in2" is the instantaneous converted relative rotor place value, and the export of module "out1" is the one phase electromagnetic torque instantaneous value. "Constant1" and "Constant2" are the rotor position increment, Δθ.
The magnetic co-energy curves data of the motor in Figure 12 is stored in MATLAB Functions by data, such as "MATLAB Fcn1", "MATLAB Fcn2" in Figure 13. During calculating for computing magnetic co-energy at the given rotor place and phase current, in view of the magnetic co-energy curves in MATLAB Functions the 2D specimen insert method has been adopted.

Simulated and Measured Results
The developed SRM (see Appendix A) drive prototype with the closed-cycle control speed conducted with use of the variable angle of the PWM combined control scheme with the PI algorithm is simulated by the proposed nonlinear simulation models on MATLAB-SIMULINK, and it is also tested experimentally at the same conditions. The photography of the SRM with the torque/speed meter and the load is shown in Figure 14, where the motor is on the left, the torque/speed meter is in the middle, and the load is on the right. The picture of the 4-phase asymmetric bridge power unit with the Intel 8XC196KC MCU controller (high performance member of the MCS 96 microcontroller family) is indicated in Figure 15, where on the right the power converter is shown, the Intel controller In the module, "ln1" is the instantaneous current value, "in2" is the instantaneous converted relative rotor place value, and the export of module "out1" is the one phase electromagnetic torque instantaneous value. "Constant1" and "Constant2" are the rotor position increment, ∆θ.
The magnetic co-energy curves data of the motor in Figure 12 is stored in MATLAB Functions by data, such as "MATLAB Fcn1", "MATLAB Fcn2" in Figure 13. During calculating for computing magnetic co-energy at the given rotor place and phase current, in view of the magnetic co-energy curves in MATLAB Functions the 2D specimen insert method has been adopted.

Simulated and Measured Results
The developed SRM (see Appendix A) drive prototype with the closed-cycle control speed conducted with use of the variable angle of the PWM combined control scheme with the PI algorithm is simulated by the proposed nonlinear simulation models on MATLAB-SIMULINK, and it is also tested experimentally at the same conditions. The photography of the SRM with the torque/speed meter and the load is shown in Figure 14, where the motor is on the left, the torque/speed meter is in the middle, and the load is on the right. The picture of the 4-phase asymmetric bridge power unit with the Intel 8XC196KC MCU controller (high performance member of the MCS 96 microcontroller family) is indicated in Figure 15, where on the right the power converter is shown, the Intel controller is on the left, and the isolation/enlargement and current/voltage protected circuit board for major switches in the power system is shown in the middle.
In the prototype, the DC voltage of the power supply is 132 V, the frequency of the PWM semaphore is 5.0 kHz, the major switches opening and turn-off angles in the power electronics is fixed based on the rotor speed ranges included in Table 1. The rating output torque is 0.5 Nm and the rating rotor speed is 900 rpm.
Energies 2020, 13, x FOR PEER REVIEW 11 of 17 is on the left, and the isolation/enlargement and current/voltage protected circuit board for major switches in the power system is shown in the middle. In the prototype, the DC voltage of the power supply is 132 V, the frequency of the PWM semaphore is 5.0 kHz, the major switches opening and turn-off angles in the power electronics is fixed based on the rotor speed ranges included in Table 1. The rating output torque is 0.5 Nm and the rating rotor speed is 900 rpm.    is on the left, and the isolation/enlargement and current/voltage protected circuit board for major switches in the power system is shown in the middle. In the prototype, the DC voltage of the power supply is 132 V, the frequency of the PWM semaphore is 5.0 kHz, the major switches opening and turn-off angles in the power electronics is fixed based on the rotor speed ranges included in Table 1. The rating output torque is 0.5 Nm and the rating rotor speed is 900 rpm.    The comparison between simulated and experimental peak value of currents is given in Table 2. As it results from this comparison, the maximum error between the calculated and measured values is not more than 8% for low speeds and not more than 6% for high speeds.   Based on Table 1, the opening angle is fixed at 2.5 • in Figures 16a and 17a, the opening angle is fixed at 0 • in Figures 16b-d and 17b-d. The turn-on angle in Figures 16a and 17a are later than those in Figures 16b-d and 17b-d, so that the ascending rate of the phase current in the former case is lower than the ascending rate of the phase current in the latter. At the same rotor speed, the load in Figures 16c and 17c is larger than the load in Figures 16b and 17b, so that the duty ratio of the PWM semaphore in Figures 16c and 17c is larger than that in Figures 16b and 17b, the phase current in Figures 16c and 17c is bigger than that in Figures 16b and 17b. Figures 16 and 17 prove that at the same rotor speed and load, simulated phase current waveforms agree basically with the experimental ones.
The curves of simulated rotor speeds are given in Figure 18 and the curves of experimentally tested rotor speeds are given in Figure 19. The comparison of simulated and tested rotor speeds is given in Table 3.  The comparison between simulated and experimental peak value of currents is given in Table 2. As it results from this comparison, the maximum error between the calculated and measured values is not more than 8% for low speeds and not more than 6% for high speeds.   Table 1, the opening angle is fixed at 2.5° in Figures 16a and 17a, the opening angle is fixed at 0° in Figures 16b-d and 17b-d. The turn-on angle in Figures 16a and 17a are later than those in Figures 16b-d and 17b-d, so that the ascending rate of the phase current in the former case is lower than the ascending rate of the phase current in the latter. At the same rotor speed, the load in Figures  16c and 17c is larger than the load in Figures 16b and 17b, so that the duty ratio of the PWM semaphore in Figures 16c and 17c is larger than that in Figures 16b and 17b, the phase current in Figures 16c and 17c is bigger than that in Figures 16b and 17b. Figures 16 and 17 prove that at the same rotor speed and load, simulated phase current waveforms agree basically with the experimental ones.
The curves of simulated rotor speeds are given in Figure 18 and the curves of experimentally tested rotor speeds are given in Figure 19. The comparison of simulated and tested rotor speeds is given in Table 3.  (c) Figure 19. Curves of the tested rotor speed: (a) the given rotor speed is n = 400 rpm, and the load is T = 0.05 Nm (ordinate: 230 rpm/div., abscissa: 2.5 s/div.); (b) n = 800r/min and T = 0.05 Nm (ordinate: 230 rpm/div., abscissa: 2.5 s/div.); (c) n = 800r/min and T = 0.50 Nm (ordinate: 570 rpm/div., abscissa: 0.5 s/div.). As can be seen, the curves of the simulated rotor speeds agree with the tested experimentally ones at the same certain rotor speed and load, and the closed-cycle rotor speed control can be implemented by the PI controller which works with high accuracy. All this proves the correctness of the adopted assumptions and the good functioning of the proposed method of SRM control.

Conclusions
The presented nonlinear simulation model of the 4-phase 8/6 poles SRM drive has been developed on a MATLAB-SIMULINK platform, which is combined with the magnetization curves data computed by the 2D electromagnetic field FE analysis (Figures 11 and 12). The nonlinear electrical meshwork model of the 4-phase asymmetric bridge power electronics main circuit and the closed-cycle rotor speed control strategy have been realized by the variable angle pulse width Figure 19. Curves of the tested rotor speed: (a) the given rotor speed is n = 400 rpm, and the load is T = 0.05 Nm (ordinate: 230 rpm/div., abscissa: 2.5 s/div.); (b) n = 800 r/min and T = 0.05 Nm (ordinate: 230 rpm/div., abscissa: 2.5 s/div.); (c) n = 800r/min and T = 0.50 Nm (ordinate: 570 rpm/div., abscissa: 0.5 s/div.). As can be seen, the curves of the simulated rotor speeds agree with the tested experimentally ones at the same certain rotor speed and load, and the closed-cycle rotor speed control can be implemented by the PI controller which works with high accuracy. All this proves the correctness of the adopted assumptions and the good functioning of the proposed method of SRM control.

Conclusions
The presented nonlinear simulation model of the 4-phase 8/6 poles SRM drive has been developed on a MATLAB-SIMULINK platform, which is combined with the magnetization curves data computed by the 2D electromagnetic field FE analysis (Figures 11 and 12). The nonlinear electrical meshwork model of the 4-phase asymmetric bridge power electronics main circuit and the closed-cycle rotor speed control strategy have been realized by the variable angle pulse width modulation combined control scheme with the PI algorithm. It includes one main module and seven sub-modules, such as: controller, one phase simulation module, rotor place angle conversion module, power electronics module, phase current computation module, "subsystem" module and one phase electromagnetic torque computation module (Figure 4). The magnetization curves data of the motor are stored in the Energies 2020, 13, 6715 15 of 17 MATLAB functions, the magnetic co-energy curves data of the motor calculated from the magnetization curves data are also stored in other MATLAB functions. The two dimensions specimen insert method was adopted in MATLAB functions for computing the flux linkage and the magnetic co-energy curves at the given rotor position and the given phase current. The waveforms of the simulated current agree very well with the waveforms of the experimentally tested ones at the same rotor speed and load (Figures 16 and 17). The simulated rotor speed curves also agree basically with the curves of the experimentally tested rotor speeds (Figures 18 and 19). It was proved that the method of the proposed nonlinear simulation models of the analyzed SRM drive is correct, and the final model is accurate. The main achievement of the work is the creation of a full stable algorithm for controlling SRM machines. The limitations of the proposed algorithm are caused by the accuracy of the SR machine FEM model and the quality of the determined torque and magnetic flux values. On the hardware side, the limitations are the speed and accuracy of the hardware used. The results obtained in this paper enable better control of SR machines in different working states. The proposed method can be applied to develop other simulation models of SRM drives with other double salient reluctance motors and different topologies of the control system and power electronics. Simulation models proposed in the paper can contribute to develop semi-physical simulation platforms which can be used for actual time simulation and-as a consequence-facilitate development of the SRM drive. In the future, research will consist in improving the proposed method and automating the connection of individual computational programs: FEM, MATLAB and SIMULINK.
Funding: This research was funded by Xuzhou Science and Technology Plan Project, grant number KH17004.

Conflicts of Interest:
The authors declare no conflict of interest.