Vectorized Mathematical Model of a Slip-Ring Induction Motor †

The paper deals with the modeling of a slip-ring induction motor. Induction motors are very often used in industry and their suitable model is needed to reduce control and operating costs. The identification process of self and mutual inductances of the stator and rotor, and mutual inductances between them in the function of the rotor rotation angle is presented. The dependence of each inductance on the rotor rotation angle is determined experimentally. The inductance matrix is then formulated. Taking the magnetic energy of the inductances and kinetic energy of the rotor into account, the Lagrange function is defined. Next, the motor motion equations are obtained. After making some algebraic transformations and using the dimensionless variables, the motion equations of electric circuits and of the mechanical equation are written separately in the forms facilitating their solution. The solution was obtained using the Simulink model for the stator and rotor currents in the form of vectors. The simulation was controlled by MATLAB script. The results of the simulation are presented in the form of basic variables time courses and compared with some values calculated with the use Steinmetz model of induction motor. The work is followed by two appendices, which contain procedures for determining the inverted inductance matrix.


Introduction
AC induction motors definitely have advantage over DC ones. They are more reliable and their purchase and maintenance costs are significantly lower [1]. This is mainly due to a simple design of the rotor and that the stator is the only element which is connected to a power supply, as well as due to economies of scale. The popularity of induction motors is confirmed by the data published in [2], where it was stated that induction motors constitute 95% of driving devices and consume up to 40-50% of the total produced electricity.
Furthermore, AC induction motors can be used in difficult operating conditions, including areas with high levels of dust, chemically aggressive atmosphere and even in explosion hazard zones as speed control and positioning systems. That is why the induction motors are commonly used in different types of industrial drive systems with velocity or position as controlled variable. It is expected that in the next decade up to 50% of all electric motors will contain induction motors powered from power electronic systems [3].
The driving torque of the induction motor results from the stator and rotor interactions. There are two types of rotors: caged and wound [2][3][4]. The design and construction of a cage rotor is simpler, but its mathematical description is more complex than a wound rotor [5]. A wound rotor has windings connected in a star. The connection point of these windings is usually isolated. The other ends of these windings are led to the rings, on which the graphite brushes slide. External resistors can be connected through these brushes. These resistors can facilitate starting of the motor and shape its operating characteristics [3]. The heat generated in the external resistors does not directly affect the motor interior temperature. As a result, the insulation of the windings and bearings age more slowly.
Slip-ring induction motors are used when high starting torque or low starting current is required. They are particularly suitable for driving systems of the machines with high inertia loads. Currently, ring motors with powers up to 20 MW are produced.
The use of induction motors with slip rings with a wounded rotor circuit connected to a variable external resistance allows speed regulation in a considerable range. However, the thermal losses of resistors associated with low speed motor operation are a serious problem. Slip ring motor drives are often using systems to recover electricity from the rotor circuit, which is rectified and returned to the power supply by means of a variable frequency drive [6].
An alternating current power system, developed by George Westinghouse, was introduced as electric power transmission system in the late 1880s and early 1890s. The first AC commutator-free induction motors were independently invented by Galileo Ferraris and Nikola Tesla, respectively in 1885 and in 1887. The three-phase circuit was first applied by Michal Doliwo-Dobrowolski in 1889. First mathematical model of the induction motor was published in 1897 by Charles Steinmetz. He proposed T-equivalent circuit model. The model is a single-phase circuit of a multiphase induction motor. On the basis of Steinmetz equivalent circuit analysis, it is possible to determine many useful relationships between circuit parameters, current, voltage, speed, power and torque. They describe how electrical input variables are transferred into mechanical output in induction motor. The Steinmetz diagram has been widely used in electrical engineering in the unchanged form for over a hundred years. In [6,7] the author uses the Steinmetz diagram, also for three-phase induction motors. The elements of the single-phase equivalent diagram are recalculated from the given parameters of three-phase circuit of the induction motor.
Some extension of the Steinmetz diagram was accomplished in [8,9] by introducing an ideal rotating transformer (IRTF) between magnetizing inductance and load resistance. The use of this element facilitated the modeling of electric machines. In the chapter on inductive machines, a universal model of winding stream connections was introduced, which enables the transformation to a two-phase circuit in frame d-q-0 that leads to a simplified machine model with IRTF. This model is the basis for a universal model of a magnetic field-oriented machine, which enables analysis and facilitates understanding of the dynamics of inductive machines. This model is the basis for the development of field-oriented control [8].
Futhermore, in [4], chapter 3 considers the induction motor scheme developed by Steinmetz. In this diagram, the resistance of losses in iron is additionally introduced in parallel to magnetization inductance and the resistance of rotor windings is distinguished. The analysis of different types of working areas and working characteristics in these areas was carried out. Chapter 4 presents, between others, a three-phase model of an induction motor and the equations of magnetic fluxes and inductance matrixes. Next, the Park transformation was applied, and four scalar differential equations of ordinary electrical part were obtained for coordinates d-q-n (equivalent d-q-0). Schemes in Simulink for scalar variables solving differential equations and sample diagrams of dynamic processes in the induction motor are also presented.
Summarizing these publications review it can be concluded that the T-equivalent model is the most commonly used one. The element of this model, which is defined as the quotient of resistance and slip, it is a non-linear circuit part. It should be stressed that according to [6] the Steinmetz model is single phased and valid only in steady-state balanced circuit condition. However in drive systems there is often an alternating mechanical load and a three-phase power supply is used, which is not always symmetrical. In addition, the asymmetry of stator and rotor circuits and the non-sinusoidality of currents are associated with the creation of a potential difference between the central points of the winding stars and supply voltages or load resistance. This phenomenon is not taken into account in the publications in question. Especially since the transition from three-phase variables to orthogonal coordinates d-q-0 is associated with a common reference point.
The magnetic streams coupling the stator and rotor windings flow twice through the air gap. Hence it can be assumed that the magnetic fluxes in the induction motor are proportional to the currents, and the stator and rotor currents should be assumed as state variables, respectively, and so applied in [5,10,11]. However, in a three-phase system without a neutral conductor there are only two independent currents. It means that the state equations of the stator and rotor circuits should be of second order. However, the state coordinates of these circuits need not be orthogonal. It can be two of three phase currents.
The purpose of the analysis is to determine the mathematical model of the slip ring induction motor. This model should make it possible to analyze both the influence of motor parameters and power quality disturbances occurring in the motor supply circuit on the output torque of the motor and the impact of dynamic load moments on the motor shaft on its supply system. The models presented above are not sufficient for these purposes The basis of the induction motor model is the inductance matrix. The form of this matrix is defined in [4,5,10,12]. The first chapter proposes a method of measuring the elements of this matrix and their values are determined. Direct measurements of the inductance matrix elements and checking its structure can be done only for the slip ring induction motor. Therefore, the paper includes an analysis of the slip ring induction motor. In practice, the model of this motor is often used as a cage induction motor model [3].
Taking into account the above remarks, the equations of the state of the electric part of the motor circuits have been recorded in the vector-matrix form of the 6th order and then transformed into the 4th order. Non-dimensional variables resulting from the equations and time scaling were used. The results of model simulation with dimensionless variables were converted into physical variables. The simulation was conducted for a balanced motor. Thanks to that, it was possible to use the Steinmetz model to verify the obtained results.
An original procedure for inverting the inductance matrix was developed. It allowed to present the equations in a form facilitating the solution of model equations in MATLAB-Simulink system. The operating diagram of the analyzed model is much smaller and simpler than the one presented in [4,13] and allows to conduct simulation experiments to study power quality disturbances and dynamic mechanical loads of the motor.

Inductance Matrix of a Slip-Ring Induction Motor
The three-phase circuits of both the stator and rotor should be considered as part of the motor dynamics analysis. In each phase of each circuit there are windings, which are mutually coupled. The type of couplings depends on the rotor rotation angle in relation to the stator. The couplings between the windings also depend on the properties of the magnetic circuit. In the induction motor, the air gap between the rotor and the stator is important. The magnetic flux generated by the currents of the stator flows through the gap twice. This allows to assume that the motor magnetic circuit is linear and may be described using inductances. That is why the base of the induction motor mathematical model is a matrix of inductances. The matrix may be described experimentally. The separated physical parts of the motor are powered and relevant elements of the inductance and resistance matrix are calculated. One drawback of this approach is that the measurements can be performed only on a real, existing machine. The results of these measurements allow us to determine which parameters are relevant for a given model. In the process of identifying the parameters of the model, the measurements are performed as the function of the rotor rotation angle. The inductance matrix L of an induction motor depends on rotor rotation angle and consists of four submatrices 3 × 3: Ls-the stator matrix of self-and mutual inductances, Msr-the matrix of mutual inductances of the rotor in relation to the stator, Lr-the rotor matrix of self-and mutual inductances and Mrs-the matrix of mutual inductances of the stator in relation to the rotor.
Taking into account the matrix L symmetry [3], it can have 21 different elements. Simultaneous identification of all these parameters causes the parameter measuring system to be rather complex. Therefore, the measurements were carried out in several stages. One of the stator or rotor windings was connected to a 50 Hz AC source and the source current and the voltages of all motor windings were measured. Measurements of both the powered winding and the remaining ones were performed using an eight-canal simultaneous measurement system. The diagram of connections for the case when one winding of the stator is powered is presented in Figure 1. The supplied winding was marked with the thick line.
Therefore, the measurements were carried out in several stages. One of the stator or rotor windings was connected to a 50 Hz AC source and the source current and the voltages of all motor windings were measured. Measurements of both the powered winding and the remaining ones were performed using an eight-canal simultaneous measurement system. The diagram of connections for the case when one winding of the stator is powered is presented in Figure 1. The supplied winding was marked with the thick line.
In the first stage of parameters determination, the supply was connected to the first winding of the stator and both the supply current and voltages: Is1, Us1, Us2 and Us3 were measured. The measured voltages are described by equations: After using the Golay-Sawitzki filter to the currents and voltages signals, the inductances and resistances were determined from the above-mentioned equations with the use of the least squares' method. To determine mutual inductances, it is necessary to use the measurements taken at the supply of chosen one winding- Figure 1. The Ls1 inductance of the winding 1 of the stator is sum of leakage inductance Ls_1 and magnetizing inductance Lsm1: The magnetizing inductance Lsm1 of the winding 1 equals the negative sum of mutual inductances of the remaining windings Lsm12 and Lsm13: In the first stage of parameters determination, the supply was connected to the first winding of the stator and both the supply current and voltages: I s1 , U s1 , U s2 and U s3 were measured. The measured voltages are described by equations: After using the Golay-Sawitzki filter to the currents and voltages signals, the inductances and resistances were determined from the above-mentioned equations with the use of the least squares' method. To determine mutual inductances, it is necessary to use the measurements taken at the supply of chosen one winding- Figure 1.
The L s1 inductance of the winding 1 of the stator is sum of leakage inductance L s_1 and magnetizing inductance L sm1 : Energies 2020, 13, 4015 5 of 20 The magnetizing inductance L sm1 of the winding 1 equals the negative sum of mutual inductances of the remaining windings L sm12 and L sm13 : The measurements were performed winch induction motor of SZUDe36a 2P44 type. It is the three-phase slip-ring and has three pairs of poles. The nominal stator current for the star connection of windings was equal to 4.1 A, the nominal rotor current for the star connection of windings was 21 A and the nominal motor speed was 920 revolutions per minute. The rated power is 1.5 kW. The measurements were executed at the supply of the windings with the current equal to approx. 0.25 of the rated current of given winding.
All stator windings were measured. The averaged measurement results showed that the stator magnetizing inductance was equal to L sm = 0.187 H and the stator leakage inductance to L s_ = 0.0293 H. The average phase resistance of windings of the stator was also measured. It amounted to R s = 10.5 Ω and contained the eddy currents resistance and the windings copper resistance equaling R sDC = 3.7 Ω, (measured with the multimeter).
The measured values of self and mutual inductance coefficients of stator phase windings were very close. The parameters differences for the circuit of different phases in the function of the rotation angle did not exceed 2% of their mean value. Therefore, it was acknowledged that they were equal for all phases of the motor and independent on the rotation angle [4,5]. This applies also to both the leakage and the magnetizing inductances. It means that the matrix Ls is symmetric and the elements of the diagonal are equal and amount to L s . The elements outside the diagonal equal −0.5L sm : where: where: Similar measurements were conducted for the rotor. It was obtained the averaged magnetizing inductance L rm = 3.9 mH and the leakage inductance L r_ = 0.55 mH from (5) and (6). The average phase resistance of the stator windings amounted to R r = 0.523 Ω and contained the resistance introduced by eddy currents and the windings (copper) resistance equaling R rDC = 0.2 Ω.
The matrix Lr has the form similar to that of Ls matrix. The elements of the diagonal are equal and amount to L r . and the elements outside the diagonal are equal to-0.5L rm . where: The mutual inductances between the stator and rotor as well as between the rotor and stator in the function of the rotor rotation angle were determined for the supply given to one winding of the Energies 2020, 13, 4015 6 of 20 stator or rotor, respectively. In case of the supply given to the rotor it is assumed that the voltages on stator windings are described by the equations: .
As in the previous case the signals of currents and voltages were filtered. The elements of the matrix of mutual inductances in the function of the rotor rotation angle were identified using the least squares method. The measurements were carried out for the supply given to the individual windings of the rotor.
The waveforms of coefficients of the mutual inductances of stator windings in relation to the rotor first phase winding versus the rotor rotation angle are presented in Figure 2.
Energies 2020, 13, x FOR PEER REVIEW 6 of 20 As in the previous case the signals of currents and voltages were filtered. The elements of the matrix of mutual inductances in the function of the rotor rotation angle were identified using the least squares method. The measurements were carried out for the supply given to the individual windings of the rotor.
The waveforms of coefficients of the mutual inductances of stator windings in relation to the rotor first phase winding versus the rotor rotation angle are presented in Figure 2. The coefficients of mutual inductances between stator and rotor are periodic functions of the rotor rotation angle. Their period amounts to 120 degrees and results from the number of poles pairs. Their amplitudes differ at most by approx. 2%. The inductances generate the cyclic matrix-circulant [14]: where: 3  cos  2  3  cos  3  cos   3  cos  3  cos  2  3  cos   2  3  cos  3  cos  3 In the next stage the stator winding was supplied and measurements of voltages on rotor windings in the function of the rotation angle were conducted. The voltages are described by an equation analogous to (10). The measured coefficients are also periodic functions of the rotor rotation angle they generate the circulant as well. Their amplitudes differ by 2%.
The amplitude of the mutual inductances among the stator and rotor equals Msr = 0.0275 H. The amplitude of the mutual inductances between the rotor and stator is approx. 4% smaller and amounts to Mrm = 0.0264 H.
In the following part of the chapter it is assumed that the Mrs and Msr inductances are equal to: with the accuracy of about 2% of Mrs and Msr values. The coefficients of mutual inductances between stator and rotor are periodic functions of the rotor rotation angle. Their period amounts to 120 degrees and results from the number of poles pairs. Their amplitudes differ at most by approx. 2%. The inductances generate the cyclic matrix-circulant [14]: where: In the next stage the stator winding was supplied and measurements of voltages on rotor windings in the function of the rotation angle were conducted. The voltages are described by an equation analogous to (10). The measured coefficients are also periodic functions of the rotor rotation angle they generate the circulant as well. Their amplitudes differ by 2%.
The amplitude of the mutual inductances among the stator and rotor equals M sr = 0.0275 H. The amplitude of the mutual inductances between the rotor and stator is approx. 4% smaller and amounts to M rm = 0.0264 H. In the following part of the chapter it is assumed that the M rs and M sr inductances are equal to: with the accuracy of about 2% of M rs and M sr values.
Basing upon the measurements we have concluded that phase resistances of the stator circuits have very close values. This applies also to the rotor resistances per phase. The phase resistances of the stator and rotor circuits are placed on diagonals of the matrices respectively Rs and Rr.
Finally, the inductance matrix of the motor is assumed in the form: The resistances matrix of the motor is also a concatenation of the stator and rotor diagonal resistances matrices: The 0 nxn denotes the zero matrix of n x n dimensions. The phase resistances of the stator and rotor circuits are placed on diagonals of the matrices: In general, in the model the stator phase winding resistances can have different values. The same can be for the rotor winding. However, basing upon the measurements, it was assumed in simulation that phase resistances of the stator circuits are equal. This applies also to the rotor resistances per phase.
The inductance matrix of a slip-ring motor (1) depends on rotor rotation angle and together with motor currents (time derivative of electric charges) determine the magnetic energy, which aggregated with the rotor kinetic energy form the Lagrange function. where: L(ϕ)-the matrix of inductances, dependent on the angle of rotor rotation in relation to a stator, Q-the motor currents column vector, ω = . ϕ-the angular velocity of the rotor, J-the moment of rotor inertia.
The Lagrange function is the difference of kinetic and potential energies and it does not take into account the friction and external forces. Consequently, the Lagrange function does not describe the energy flow [15].
Basing on the Lagrange function the d'Alembert-Lagrange equation may be formulated [15]. Using the equation with virtual velocities as variations of virtual coordinates it may be easy checked that induction motors are holonomic systems. Thanks to that, after using the Euler-Lagrange equations and taking into account forces of friction and external excitation, the motion equations may be obtained in the form: d dt d dt where Uz-denotes the supply voltage vector, Uo-the vector of voltages between s neutral points, R-the windings resistance matrix, T L -the mechanical load torque, T F -the mechanical friction torque After substitution of Lagrange function into the Equations (19) and (20) and some transformations we get The vectors of currents and voltages are of the 6th order and it refer to concatenated three-phase stator and rotor vectors. So: where: Is = [I s1 I s2 I s3 ] T , Ir = [I r1 I r2 I r3 ] T , 1 3x1 = [1, 1, 1] T , U os , denote voltage between center points of star power supply and star winding for the stator, U or -voltage between center points of star winding of the rotor and the star load or power supply of rotor, Us, Ur-are vectors of the phase voltage that supply power to the stator and rotor windings. The motor inductance matrix L(ϕ) and the resistance matrix R are described as (15) and (16) respectively. The voltages supplying the stator are assumed in form: where: Es-the diagonal matrix of magnitudes and time functions of power supply phase voltages.
where: q = 2 · π/3. The mean value of the phase voltages amplitudes is described as: The (21) components are column vectors of voltages of the 6th order. Their first three rows represent the stator equation, the other three are the rotor one.
In order to simplify the writing to the dimensionless form, the voltage and currents of the rotor are transferred to the stator level. For this purpose, the stator equations are divided by the mean amplitude of the phase supply voltages Es, the rotor equations are divided by Es·(M/L sm ). Next the time scaling τ = ω s ·t is used and then the substitution of the value of stator and rotor currents divided by Es/(ω s L sm ) and Es/(ω s M) respectively is performed. After the dimensionless variables' definition: i s = I s /(Es/ω s L sm ), i r = I r /(Es/(ω s M)) r s = Rs/(ω s L sm ), r r = Rr/(ω s L rm ) u os = U os /Es, u or = U or /(Es · (M/L sm )) e s = Es/Es, u r = Ur/(Es · (M/L sm )) , (21) may be written in form. where: A circle above the state variables denotes the time derivative calculated in relation to τ, which describes the time after time scaling.
The three-phase equation of a stator and a rotor may be described by two instantaneous values of currents. It means that the stator and rotor circuit may be also described by two equation system of the second order. Therefore, the stator (or/and rotor) equations of the 1st and 2nd phase will be used only for the analysis. The 3rd phase current is replaced with a negative sum of the current of the first and second phases. The replacing process for the stator and rotor currents is shown as T 32 transformation [14]. where: In order to eliminate voltages between center points of stars of power supply and stator windings, the third phase equation should be subtracted from the first-phase equation and from the second-phase equation of stator. The same should be done with the equations of the rotor. The above subtractions correspond to a premultiplication by a matrix [16] separately for stator and rotor equations: After transformations the equations of the electric part of the motor may be expressed as follows: where: The Λ and Θ matrices are sized 6 × 6 dimensions. After applying the T 23 and T 32 transformations, the Λ 2 and Θ 2 matrices are 4 × 4 dimensions and their elements-submatrices are 2 × 2 dimensions. The values and designations of these sub-arrays for Λ 2 are as follows: where: Similarly, elements of the Θ 2 matrix are in the form: mr s2 = T 23 · r s · T 32 = r s1 + r s3 r s3 r s3 r s2 + r s3 , (41) mr r2 = T 23 · r r · T 32 = r r1 + r r3 r r3 r r3 r r2 + r r3 , Computing elements of matrices (34) and (35) yields possibility calculation of currents derivatives using equation: The Λ 2 −1 -inverted matrix may be determined using the procedures presented in Appendices A and B. The mechanical Equation (22) of the motor may also be transformed. After substitution of wet friction: and time scaling, it has the following form: where: The mathematical model describes the inductive slip-ring motor with symmetrical inductances. The remaining elements of the model may be asymmetric, i.e., they may have unequal phase components.
The tests conducted on the model may be divided into the analysis of the model of the symmetrical arrangement and sensitivity studies of the remaining elements affecting its characteristics. The model uses dimensionless variables and parameters and that is why it is simple and maybe useful in simulation of the slip-ring induction motor.

The Model of the Motor in Simulink
The motor electric and mechanical Equations (44) and (46) respectively, provide the basis for the creation of the vectorized simulation model in Simulink. The vectorization simplifies the diagram of the model and facilitates its use. The model diagram is presented in Figure 3. The physical parameters of the modeled motor are presented in Table 1. Physical variables are converted into dimensionless Simulink input variables and simulation output variables into physical variables of the object in the MATLAB script, which controls the simulation experiment.
The  The physical parameters of the modeled motor are presented in Table 1. Physical variables are converted into dimensionless Simulink input variables and simulation output variables into physical variables of the object in the MATLAB script, which controls the simulation experiment.
The The simplicity of the final Simulink model ought to be emphasized. It results from the application of the transformation of the circuit equation to (44) and the use of vector-matrix notation. The simplicity of the model may be estimated by comparing it with the models presented in [4,13,17]. The obtained model is also more convenient for the arrangement of simulation experiments.

T-Equivalent Induction Motor Model
For a rough check of the modelling results, it was decided to use the possibly simple model of an induction motor. It was decided to use the Steinmetz model. His replacement diagram of a multi-phase induction motor in the form of T-equivalent circuit is presented in Figure 4.

T-Equivalent Induction Motor Model
For a rough check of the modelling results, it was decided to use the possibly simple model of an induction motor. It was decided to use the Steinmetz model. His replacement diagram of a multiphase induction motor in the form of T-equivalent circuit is presented in Figure 4.
Synchronous angular velocities of the stator's magnetic field rotation and the rotor's angular speed were determined as s  , r  respectively. The element on which the output load is generated in this scheme is the resistance with the value Rr'/s. Steinmetz presented the power emitted on this resistance as a sum of electromechanical output power and thermal power. To determine these powers, he introduced the following equation: Multiplying this equation by 2 3 r I   the following equation is obtained in which the Pgap is air gap power, Pem-electromechanical output power and Pr-the heat power generated on the resistance. The above active power of rotor is respectively equal to: To simplify further analysis, the IEEE recommends using Thevenin's claim. The result is the diagram shown in Figure 5. The following components are used in the diagram: • R s , X s -stator resistance and leakage reactance of stator, • X m -motor magnetizing reactance, • R r , X r -rotor resistance and leakage reactance of rotor transformed to stator side, This is a single-phase model of a multiphase induction motor. The model of the motor is valid in steady-state balanced circuit condition [6]. The slip of the motor is defined as Synchronous angular velocities of the stator's magnetic field rotation and the rotor's angular speed were determined as ω s , ω r respectively. The element on which the output load is generated in this scheme is the resistance with the value R r '/s. Steinmetz presented the power emitted on this resistance as a sum of electromechanical output power and thermal power. To determine these powers, he introduced the following equation: Multiplying this equation by 3 · I r 2 the following equation is obtained in which the P gap is air gap power, P em -electromechanical output power and P r -the heat power generated on the resistance. The above active power of rotor is respectively equal to: To simplify further analysis, the IEEE recommends using Thevenin's claim. The result is the diagram shown in Figure 5. In this diagram there is a voltage source with VTE voltage and ZTE impedance in series with rotor reactance and load resistance. These variables are equal to: Using the above result, it is possible to determine the stator current from the circuit in Figure 4: The electromechanical output power can be determined from (52) and on this basis the electromechanical moment may be calculated.
It should be stressed that the above dependencies result from the scheme adopted by Steinmetz and the assumption that the individual phase circuits of a three-phase induction motor can be treated as working independently on torque.

Example of the Model Waveforms
The results of the simulation were compared with the values determined for the same motor parameters using the Steinmetz model, which allows us to determine the motor currents and powers. However, the full state vector in this model cannot be reproduced as the rotor speed is not available. The simulation model is used to determine the steady state speed as a function of the load torque. The course of this speed during motor start-up for the load torque TL = 1 Nm and TL = 15 Nm is shown in Figure 6. These graphs show that after just one period of run there is an increase in motor speed equal to about 20% of steady state speed. The time of reaching the steady state is equal to about 5 to 10 periods of supply voltage. The slip for the load torque TL = 1 Nm is about 0.025 and for TL = 15 Nm it is 0.435. In this diagram there is a voltage source with V TE voltage and Z TE impedance in series with rotor reactance and load resistance. These variables are equal to: where Z s = R s + jX s . The current I r ' flowing in the load may be calculated as follows: Using the above result, it is possible to determine the stator current from the circuit in Figure 4: The electromechanical output power can be determined from (52) and on this basis the electromechanical moment may be calculated.
It should be stressed that the above dependencies result from the scheme adopted by Steinmetz and the assumption that the individual phase circuits of a three-phase induction motor can be treated as working independently on torque.

Example of the Model Waveforms
The results of the simulation were compared with the values determined for the same motor parameters using the Steinmetz model, which allows us to determine the motor currents and powers. However, the full state vector in this model cannot be reproduced as the rotor speed is not available. The simulation model is used to determine the steady state speed as a function of the load torque. The course of this speed during motor start-up for the load torque T L = 1 Nm and T L = 15 Nm is shown in Figure 6. These graphs show that after just one period of run there is an increase in motor speed equal to about 20% of steady state speed. The time of reaching the steady state is equal to about 5 to 10 periods of supply voltage. The slip for the load torque T L = 1 Nm is about 0.025 and for T L = 15 Nm it is 0.435.  Figure 7 shows the rotor currents for load torque TL = 1 Nm. In steady state they had an amplitude equal to 2 A and a period of 0.82 s which corresponds to a slip equal to 0.025. The current amplitude at the beginning of the start-up was 48 A. In Figure 8, for torque TL = 15 Nm in steady state it was about 5.5 A and the period was about 0.046 s. At the beginning of the start-up process, the current amplitude reached 50 A. The determined periods of current oscillation of the rotor are in accordance with the rotor speed courses in Figure 6. The frequency and the magnitude of the currents waveforms of the rotor decreases with the reduction in the load torque.
Starting currents can be estimated using Steinmetz model dependencies. The observed model data differs from those calculated for s = 1 using (55) by about 5 percent of the later ones. This can be explained by a fairly high acceleration during the first period and a reduction in slip s, which also reduces the rotor current.  Figure 7 shows the rotor currents for load torque T L = 1 Nm. In steady state they had an amplitude equal to 2 A and a period of 0.82 s which corresponds to a slip equal to 0.025. The current amplitude at the beginning of the start-up was 48 A.  Figure 7 shows the rotor currents for load torque TL = 1 Nm. In steady state they had an amplitude equal to 2 A and a period of 0.82 s which corresponds to a slip equal to 0.025. The current amplitude at the beginning of the start-up was 48 A. In Figure 8, for torque TL = 15 Nm in steady state it was about 5.5 A and the period was about 0.046 s. At the beginning of the start-up process, the current amplitude reached 50 A. The determined periods of current oscillation of the rotor are in accordance with the rotor speed courses in Figure 6. The frequency and the magnitude of the currents waveforms of the rotor decreases with the reduction in the load torque.
Starting currents can be estimated using Steinmetz model dependencies. The observed model data differs from those calculated for s = 1 using (55) by about 5 percent of the later ones. This can be explained by a fairly high acceleration during the first period and a reduction in slip s, which also reduces the rotor current. In Figure 8, for torque T L = 15 Nm in steady state it was about 5.5 A and the period was about 0.046 s. At the beginning of the start-up process, the current amplitude reached 50 A. The determined periods of current oscillation of the rotor are in accordance with the rotor speed courses in Figure 6. The frequency and the magnitude of the currents waveforms of the rotor decreases with the reduction in the load torque.
Starting currents can be estimated using Steinmetz model dependencies. The observed model data differs from those calculated for s = 1 using (55) by about 5 percent of the later ones. This can be explained by a fairly high acceleration during the first period and a reduction in slip s, which also reduces the rotor current.  The initial value (for s = 1) determined from relation (56) is about 10.4 A as in these figures. Differences occur for steady state stator current amplitudes. For TL = 1 Nm the calculation shows amplitude of 4.7 A and the simulation in Figure 9 shows a value of about 3.3 A. Similarly, for TL = 15 Nm the calculation gives 6.2 A and the graph in Figure 10 shows 5.5 A.
The stator current waveforms presented in Figures 9 and 10, have two stages: the start-up and steady-state. The shape of the waveforms during the stages is clear. In the case of a smaller load torque the current stabilizes after only four periods of the power supply voltage. The magnitudes of currents in steady-state rise with the increase in load torque value. The stator currents waveforms frequency is independent of the load torque and constant, and equal to the frequency of the supply network. From the rotor current, it is possible to determine the value of the stator current amplitude at the beginning of the motor start-up process for both load torque values. These currents are shown in Figures 9 and 10. The initial value (for s = 1) determined from relation (56) is about 10.4 A as in these figures. Differences occur for steady state stator current amplitudes. For TL = 1 Nm the calculation shows amplitude of 4.7 A and the simulation in Figure 9 shows a value of about 3.3 A. Similarly, for TL = 15 Nm the calculation gives 6.2 A and the graph in Figure 10 shows 5.5 A.
The stator current waveforms presented in Figures 9 and 10, have two stages: the start-up and steady-state. The shape of the waveforms during the stages is clear. In the case of a smaller load torque the current stabilizes after only four periods of the power supply voltage. The magnitudes of currents in steady-state rise with the increase in load torque value. The stator currents waveforms frequency is independent of the load torque and constant, and equal to the frequency of the supply network. The initial value (for s = 1) determined from relation (56) is about 10.4 A as in these figures. Differences occur for steady state stator current amplitudes. For T L = 1 Nm the calculation shows amplitude of 4.7 A and the simulation in Figure 9 shows a value of about 3.3 A. Similarly, for T L = 15 Nm the calculation gives 6.2 A and the graph in Figure 10 shows 5.5 A.
The stator current waveforms presented in Figures 9 and 10, have two stages: the start-up and steady-state. The shape of the waveforms during the stages is clear. In the case of a smaller load torque the current stabilizes after only four periods of the power supply voltage. The magnitudes of currents in steady-state rise with the increase in load torque value. The stator currents waveforms frequency is independent of the load torque and constant, and equal to the frequency of the supply network. The electric torque during the start-up motor for the different load torque values is presented in Figure 11. In the steady state, the electromagnetic torque waveforms are close to the load torque. This is due to the inclusion of the mechanical part in the simulation model. The Steinmetz model for load torque TL = 1 Nm, gives torque value approximately 30% lower than the simulation steady state. The compliance was obtained for the higher value of the load moment. It seems that it may result from the simplification of the T-equivalent model. The T-model is sensitive for the slip value close to zero. In both cases the constant component of the start-up torque appears. It determines the dynamics of the start-up process of the motor. After reaching a specified angular speed, the electric torque converges to the sum of torques of loads and friction.

Conclusions
The model of slip-ring induction motor presented in this paper was obtained on the basis of inductance matrix in the form specified in the literature. The elements of this matrix were determined experimentally for the selected motor using methods developed by the authors. After the Lagrange function was formulated, motion equations of the motor were determined, dimensionless variables The electric torque during the start-up motor for the different load torque values is presented in Figure 11. The electric torque during the start-up motor for the different load torque values is presented in Figure 11. In the steady state, the electromagnetic torque waveforms are close to the load torque. This is due to the inclusion of the mechanical part in the simulation model. The Steinmetz model for load torque TL = 1 Nm, gives torque value approximately 30% lower than the simulation steady state. The compliance was obtained for the higher value of the load moment. It seems that it may result from the simplification of the T-equivalent model. The T-model is sensitive for the slip value close to zero. In both cases the constant component of the start-up torque appears. It determines the dynamics of the start-up process of the motor. After reaching a specified angular speed, the electric torque converges to the sum of torques of loads and friction.

Conclusions
The model of slip-ring induction motor presented in this paper was obtained on the basis of inductance matrix in the form specified in the literature. The elements of this matrix were determined experimentally for the selected motor using methods developed by the authors. After the Lagrange function was formulated, motion equations of the motor were determined, dimensionless variables In the steady state, the electromagnetic torque waveforms are close to the load torque. This is due to the inclusion of the mechanical part in the simulation model. The Steinmetz model for load torque T L = 1 Nm, gives torque value approximately 30% lower than the simulation steady state. The compliance was obtained for the higher value of the load moment. It seems that it may result from the simplification of the T-equivalent model. The T-model is sensitive for the slip value close to zero. In both cases the constant component of the start-up torque appears. It determines the dynamics of the start-up process of the motor. After reaching a specified angular speed, the electric torque converges to the sum of torques of loads and friction.

Conclusions
The model of slip-ring induction motor presented in this paper was obtained on the basis of inductance matrix in the form specified in the literature. The elements of this matrix were determined experimentally for the selected motor using methods developed by the authors. After the Lagrange function was formulated, motion equations of the motor were determined, dimensionless variables were introduced and the minimal form of motion equations of the electric part of the induction motor model were established. These equations are recorded in the vector-matrix form of the fourth order. These are non-linear equations, with non-linearity resulting from the dependence of the inductance matrix on the angle of rotor rotation.
These equations are original. The form of these equations allowed to formulate a simple operational diagram of the slip-ring induction motor model in Simulink. In order to check the correctness of the model, the T-equivalent diagram of the induction motor was presented. The correspondence between the initial period of the motor start processes and some variables of the steady state vector was obtained. Such possibilities were made possible using this additional model.
Three-phase circuits of the motor stator and rotor are described using differential equations of two-phase currents of the stator and of two-phase currents of the rotor circuits. As a result, the motor circuits are more easily observed and the forms of the equations are simpler.
The model may be easily adapted for the research on the influence of asymmetry of power supply, dynamical load of the motor and influence of unbalance of other components of motor circuits. The preliminary results of tests are promising. The tests conducted on the model may be divided into the analysis of the model of the balanced (symmetrical) arrangement and sensitivity studies of the remaining elements affecting its characteristics as in [14].
The model makes simulation of a slip-ring induction motor simpler than the ones in [9,13,17]. The authors consider a concept of the presented mathematical transformations and the model as a novelty.
From (A4) we get: and after substituting into (A3): Hence The substitution of the above into (A7) yields: Similarly, the remaining submatrices can be obtained as follows: