Determination of Combination of Input Values of Simulated Motor Using Taguchi Method
Determination of Combination of Input Values of Simulated Motor Using Taguchi Method
ISSN No:-2456-2165
Abstract:- The most popular type of motors in use today For dynamic simulation purposes, it is important to
is squirrel-cage induction motors. All types of motor know first, that the output parameters of a simulated motor
faults can be studied if a motor could be simulated based and a actual motor are the same. Rotational speed (rpm),
on the motor without faults. The goal of this study was to torque, and stator current are selected as output parameters of
determine the input optimum values for making it almost a simulated motor because these parameters are given on the
identical values of outputs f the simulated motor to the nameplate data on the motor. The parameters on the
target values on the nameplate data of the actual motor. nameplate of a motor become the parameters to check output
In this study, the Taguchi method, calculated using parameters of a simulated motor, correct or not. Output
MATLAB, was used to optimize the values input of the parameters from a simulated model should be very close to
simulated motors. For dynamic simulation purposes, the nameplate data on the actual motor if the model is robust.
inputs and outputs of a simulated motor need to be In this study, the parameters on the nameplate of the motor
selected. In this study, stator resistance, stator are called the targets and outputs of the simulated model
inductance, rotor bar resistance, rotor bar inductance, developed in Simulink MATLAB are called the output
rotor end ring resistance, and rotor end ring were selected parameters.
as the input parameters. Rotations per minute (rpm),
torque, and stator current were selected as output An actual motor with 11 kW and 415V is studied in this
parameters; they are also provided on the nameplate data research work. All the parameters of a normally-functioning
on the actual motor. These results allow us to develop actual motor are given in Appendix 1. Inputs of the dynamic
dynamic simulations to study a variety of motor faults simulation model are classified into three categories: 1)
and to compare the results of a simulated motor with and parameters on the nameplate from the manufacturer; 2)
without faults, to better understand and analyze the calculated parameters from the material and the actual
frequency domain, determine the root cause of the dimensions of the motor; and 3) measured parameters in the
vibrations. All code work is shown in the paper and one laboratory. Some parameters such as the number of stator
new formula was developed based on this study. coils/pole/phase, resistances and inductances of stator and
rotor bars and so on are calculated and/or measured
Keywords:- Model, Simulation, Optimization, Modelling, experimentally if not supplied by the manufacturers.
Dynamic Simulation, MATLAB. 1. Data on the nameplate from the manufacturer: nominal
voltage [V], nominal torque [Hm], nominal current [A],
I. INTRODUCTION number of poles, rotational speed [rpm], frequency [Hz],
/𝑌 connection, efficiency, power factor, number of stator
Squirrel induction motors are used mainly in industry. coils/poles/phases.
Experiments on motors, conducted in a laboratory setting, are 2. Calculated data: number of turns in the stator coils, air gap
expensive and limited, because the motor should be broken to average length [m], synchronous speed [rpm], slip [%],
create the fault to be studied. Another motor is again needed rotor angular velocity [rps], rotor mechanical speed [hz],
to be broken to study a different fault. Large capacity motors rotor frequency [hz], stator slot pitch [rad], pole pitch
are generally not allowed to break because of cost and safety [rad], number of stator coils/pole/phase, number of turns
factors. However, conducting research on a simulated motor, in the stator.
with and without faults, saves an enormous amount of time 3. Data from the laboratory: stator layout, number of rotor
and money. Motors of any capacity, with any type of fault bars, number of stator slots, stator slot pitch [slots], rotor
can be studied if a motor could be simulated in Simulink stack length [m], average radius [m], stator resistance
MATLAB or other software. In the present study, an actual [ohm], stator leakage inductance [H], rotor bar resistance
motor without a fault needed to be simulated. After [ohm], rotor bar inductance [H], rotor end ring resistance
developing a dynamic simulation, several types of motor [ohm], rotor end ring inductance [H].
faults were studied, based on the simulated motor without
faults. It is necessary to adequately know the electric, Six parameters such as stator resistance; stator
magnetic, and mechanical behavior of motors without faults inductance; rotor bar resistance; rotor bar inductance; rotor
for a reliable diagnostic method of fault identification. end ring resistance; and rotor end ring inductance from the
input parameters of a simulated motor, are selected as input
parameters, because these parameters are the most influential
TABLE I
INPUT PARAMETERS
Experimental levels
Input parameter Notation
1 (-10%) 2 3 (10%)
Resistance of stator Rs Rs-0. 1*Rs Rs Rs+0. 1*Rs
Resistance of rotor bar Rb Rb-0.1*Rb Rb Rb+0. 1*Rb
Resistance of rotor end-ring Re Re-0.1Re Re Re+0.1Re
Inductance of stator Ls Lls-0. 1*Lls Lls Lls+0.1*Lls
Inductance of rotor bar Lb Llb-0. 1*Llb Llb Llb+0.1*Llb
Inductance of end ring Le Lle-0. 1*Le Lle Lle+0.1*Lle
The values of input parameters of resistance of stator, resistance of rotor bar, resistance of rotor end-ring, inductance of
stator, inductance of rotor bar, inductance of end ring are listed in Table II.
TABLE II
VALUES OF INPUT PARAMETERS
Rotational speed (RPM) Torque (Te) Stator current (Is)
Notation 1 (- 3
1 (-10%) 2 3 (10%) 1 (-10%) 2 3 (10%) 2
10%) (10%)
Rs 2.35 2.608454802 2.87 3.58 3.98 4.37 4.05 4.50 4.95
1.59703E- 1.95193E- 3.31.E-05 6.67.E-05
Rb 1.77448E-05
05 05 3.64E-05 3.6E-05 7.3E-05 7.3E-05
1.37214E- 1.67706E- 1.37.E-06 3.10.E-06
Re 1.5246E-06
06 06 1.51E-06 1.5E-06 3.4E-06 3.4E-06
Ls 0.01458 0.0162 0.01782 0.02178 1.98.E-02 0.02178 0.0033 3.00.E-03 0.0033
1.73696E- 2.12295E- 1.75.E-06 4.44.E-07
Lb 1.92995E-06
06 06 1.93E-06 1.9E-06 4.9E-07 4.9E-07
1.3704E- 1.67493E- 1.04.E-07 2.22.E-07
Le 1.52266E-07
07 07 1.14E-07 1.1E-07 2.4E-07 2.4E-07
An efficient way to optimize the input parameters simultaneously is to plan experiments using L 18 orthogonal arrays [1].
Code to distribute the six input parameters within the L18 orthogonal array was developed specifically for this project by the
author [2] and the results are shown in Table III.
The most influential six input parameters of RPM are determined by values in each row. 18 rows mean 18 experiments.
These 18 experiments are repeated until the minimum differences are achieved between the output parameters of the model, and
the target (parameters on the nameplate of the motors). This is the same processes of the torque (Te), the stator current (Is). The
final best results of targets such as torque, rotational speed, and stator are shown Table IV.
TABLE IV
OUTPUT PARAMETERS
RPM Te Is RPM Te Is
№ №
1450 73.69 20.8 1450 73.69 20.8
1 1452.133 74.615 20.228 10 1449.512 74.302 20.768
2 1449.811 74.337 20.644 11 1451.055 74.483 20.529
3 1447.175 74.008 21.123 12 1448.834 74.223 20.639
4 1449.839 74.312 20.991 13 1448.580 74.166 21.060
5 1447.344 74.010 21.311 14 1448.623 74.204 20.879
6 1447.737 74.106 20.761 15 1447.788 74.067 21.101
7 1448.524 74.159 21.318 16 1447.117 73.987 21.522
8 1446.365 73.916 21.246 17 1446.335 73.858 21.573
9 1445.179 73.726 21.696 18 1446.672 73.951 21.172
The 18 output values obtained experimentally were used to generate ‘the nominal the best’ loss function using [3] and [4].
The analysis of variance (ANOVA) method was used to determine the percentage contribution of each of the parameters
studied; we used this method to quantitatively estimate the input that the parameters made to the overall outputs, expressed as a
percentage, as described previously [1] - [5]. The fundamental techniques of ANOVA involve the total sum of squares, the
variance ratio, and degrees of freedom.
To simplify the analysis, code was developed that generated ‘the nominal the best,’ (S/N) ratios, and sum of squares [6].
The S/N ratio and values for the sums of the squares of, outputs (are torque, rotations per minute, stator current), are
tabulated in Table V.
The evaluated S/N ratios (S/Nij ) are determined by the S/N ratios that were calculated using values produced from the
experimental distributions (see Table III). The evaluated S/N ratio is calculated by (1), which was developed based on this
research study.
∑ni=1 S/Nirj
S/Nij = (1)
nk
Here: S/Nirj − sum of the S/N ratio in the particular column, nk − number of non-zero values of the S/N ratio.
Step 1. Distribution of 1s, 2s, and 3s into orthogonal arrays: The Rs, Rb, Re, Ls, Lb, Le columns of of orthogonal arrays kept in
orth18.txt have been divided into three columns each, matching the three experimental levels. For example: the Column Rs is
divided into columns labeled Rs1, Rs2 and Rs3. All 1s are in Rs1, all 2s are in Rs2, & all 3s are in Rs3. All the rest of the parameters
follow this same process.
Step 2. Replacement by S/N ratio: All 1s in Rs1, Rb1, Re1,Ls1,Lb1,Le1, and all 2s in Rs2, Rb2, Re2, Ls2, Lb2, Le2 and all threes in Rs3,
Rb3, Re3, Ls3, Lb3, Le3 were replaced by the S/N ratios of rotational speed, torque and stator current corresponding to the
experimental distribution of the row inTable V. The evaluated S/N ratio is calculated by (1).
Step 3. Reconstruct a table: The evaluated S/N ratios of RPM are reconstructed into the table with six rows and three columns.
Step 4. Plot a Taguchi graph: Taguchi graphs are plotted by the evaluated S/N ratios of each parameter (each row).
The evaluated S/N ratios of torque, stator current are calculated the same way as the evaluated S/N ratios of rotational speed.
Codes for Steps 1 - 4 were developed by the author [7] and the results are shown in Tables VII-IX.
TABLE VII
DISTRIBUTION OF ORTHOGONAL ARRAYS
Input parameters Distribution of input parameters
R L Rs Rb Re Ls Lb Le
s b e s b e 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3 1 2 3
1 1 1 1 1 1 1 1 1 1 1 1
1 2 2 2 2 2 1 2 2 2 2 2
1 3 3 3 3 3 1 3 3 3 3 3
2 1 1 2 2 3 2 1 1 2 2 3
2 2 2 3 3 1 2 2 2 3 3 1
2 3 3 1 1 2 2 3 3 1 1 2
3 1 2 1 3 3 3 1 2 1 3 3
3 2 3 2 1 1 3 2 3 2 1 1
3 3 1 3 2 2 3 3 1 3 2 2
1 1 3 3 2 1 1 1 3 3 2 1
1 2 1 1 3 2 1 2 1 1 3 2
1 3 2 2 1 3 1 3 2 2 1 3
2 1 2 3 1 2 2 1 2 3 1 2
2 2 3 1 2 3 2 2 3 1 2 3
2 3 1 2 3 1 2 3 1 2 3 1
3 1 3 2 3 2 3 1 3 2 3 2
3 2 1 3 1 3 3 2 1 3 1 3
3 3 2 1 2 1 3 3 2 1 2 1
Based on Table VII the results of the evaluated S/N ratios of rotational speed (RPM), torque (Te) and stator current (Is) were
calculated and are shown in Table VIII.
TABLE VIII
DISTRIBUTION OF S/N RATIO AND S/N RATIOS
S/N Rs Rs Rs3 Rb Rb Rb Re1 Re Re3 Ls Ls2 Ls3 Lb Lb Lb Le1 Le2 Le3
1 2 1 2 3 2 1 1 2 3
Rotational Speed (RPM)
-6.6 - - -6.6 - -6.6 -6.6
6.6 6.6 6.6
14.5 14. 14. 14. 14. 14. 14.
5 5 5 5 5 5
-9.2 - -9.2 -9.2 -9.2 - -9.2
9.2 9.2
15.8 15. 15. 15. 15. 15. 15.
8 8 8 8 8 8
-8.5 - -8.5 - -8.5 - -8.5
8.5 8.5 8.5
-7.9 - -7.9 -7.9 - -7.9 -7.9
7.9 7.9
The evaluated S/N ratios of rotational speed, stator current and torque among all levels can be determined by using the
values from the final row in Table VIII. The evaluated S/N ratios and calculated values are shown in Table IX.
TABLE IX
THE EVALUATED S/N RATIOS AND CALCULATED VALUES
Parameter Evaluated S/N ratios at Degrees Sum of Mean Variance % of Error of
r experimental levels (L) of Squares square ratio variance mean
1 2 3 Freedom (SS_S/Nij) (Msq) Fr ratio Fr square
(DOFr) (eMsq)
Rotation per minute (rpm)
Rss 0.546 -2.076 -9.862 2 351.680 175.840 6.460 29.379
Rb -0.024 -3.294 -8.075 2 196.756 98.378 3.614 16.437
Re -3.839 -2.040 -5.513 2 36.192 18.096 0.665 3.023
27.218
Lls -5.123 0.276 -6.545 2 155.396 77.698 2.855 12.981
Llb -6.757 1.605 -6.240 2 263.476 131.738 4.840 22.010
Lle -6.232 -3.168 -1.993 2 57.469 28.734 1.056 4.801
Torque (Te)
Rss 4.359 7.058 14.350 2 320.576 160.288 18.062 46.361
Rb 5.440 8.314 12.014 2 130.331 65.165 7.343 18.848
Re 9.949 7.304 8.514 2 21.037 10.518 1.185 3.042
8.874
Lls 5.721 7.547 12.499 2 147.594 73.797 8.316 21.345
Llb 8.104 9.755 7.908 2 12.347 6.174 0.696 1.786
Lle 7.978 9.889 7.901 2 15.216 7.608 0.857 2.201
Stator current (Is)
Rss 14.639 15.423 4.558 2 440.587 220.294 4.226 37.132
Rb 11.541 10.778 12.300 2 6.949 3.475 0.067 0.586
Re 7.361 10.646 16.613 2 264.002 132.001 2.532 22.250
52.123
Lls 13.452 11.111 10.056 2 36.263 18.131 0.348 3.056
Llb 11.634 15.324 7.661 2 176.220 88.110 1.690 14.851
Lle 11.088 11.846 11.685 2 1.914 0.957 0.018 0.161
Lb2
Rs1 Ls2
Rb1
Rs3
Ls3
Rb3
Rb2 Re3
Lb1 Lb3 Le1 Le3
Re2 Ls2
Rs2
Rb1 Ls1
Rs1
Re3
Rs2 Lb2
Rs1
Ls1
Rb3
Rb1 Lb1 Le2 Le3
Rb2 Ls2 Le1
Re2
Ls3
Re1 Lb3
Rs3
Here, the optimum conditions were: Rs2 Rb3 Re3 Ls3 Lb3 Le3 in the RPM.
The non-optimum input parameters are shown in Table X and the optimum input parameters are shown in Table XI
TABLE X
NON-OPTIMUM INPUT PARAMETERS AND OUTPUTS CORRESPONDING THE NON-OPTIMUM INPUT
PARAMETERS
Percentage to close to the
Notatio Optimum input Output parameters target
Input parameter
n values RPM Te Is
RPM Te Is
1450 73.69 20.8
Rotational speed (RPM)
Resistance of stator Rs 3.58
Resistance of rotor bar Rb 2.43E-05
Resistance of rotor end-
Re 1.68E-06 97.03
ring 1449.76 74.62 20.18 99.98% 98.76%
%
Inductance of stator Ls 1.80E-02
Inductance of rotor bar Lb 1.93E-06
Inductance of end ring Le 1.14E-07
Torque (te)
Resistance of stator Rs 4.37
Resistance of rotor bar Rb 3.64E-05
Resistance of rotor end-
Re 99.79
ring 1.66E-06 1442.27 73.85 20.76 99.47% 99.79%
%
Inductance of stator Ls 1.78E-02
Inductance of rotor bar Lb 1.74E-06
Inductance of end ring Le 9.36E-08
Stator Current (Is)
Resistance of stator Rs 3.98
Resistance of rotor bar Rb 3.31E-05
Resistance of rotor end-
Re 99.12
ring 1.37E-06 1445.16 74.15 20.62 99.67% 99.39%
%
Inductance of stator Ls 1.98E-02
Inductance of rotor bar Lb 1.75E-06
Inductance of end ring Le 1.04E-07
TABLE XI
OPTIMUM INPUT PARAMETERS AND OUTPUTS CORRESPONDING THE OPTIMUM INPUT PARAMETERS
Percentage to close to the
Notatio Optimum input Output parameters target
Input parameter
n values RPM Te Is
RPM Te Is
1450 73.69 20.8
Resistance of stator Rs 4.50
Resistance of rotor bar Rb 3.31E-03
Resistance of rotor end- 1.00E-04
Re 74.18
ring 973.78 33.87 15.43 67.16% 45.96%
%
Inductance of stator Ls 1.98E-02
Inductance of rotor bar Lb 1.75E-06
Inductance of end ring Le 1.04E-07
“Values at the optimum levels” means optimum inputs to make the minimum difference between outputs and targets. As
shown in Table XI, each of the values at the optimum levels has produced different outputs. As can be seen in the torque and
stator current columns in Table XI, stator current and torque are increasing, but rpm is decreasing compared to the value in the
column for the optimum level of rotational speed. It means that the optimum inputs of rotational speed to make the minimum
difference between output and targets of rotational speed could not make the minimum difference between output and targets of
stator current and torque. From the results of this study, we have found that the Taguchi method is useful when optimizing single
parameter performance characteristics. On the other hand, this method is too complicated a procedure to be used to determine
input parameters for optimizing more than one output, simultaneously.
RPM
500
0
0 0.5 1 1.5 0
0 0.2 0.4 0.6 0.8 1 1.2 1.4
-1000 time,c -500 time,c
150 100 Torque
Torque
100
50
Te
Te
50
0
0 0 0.2 0.4 0.6 0.8 1 1.2 1.4
0 0.5 1 1.5
-50 time,c -50 time,c
100 Stator Current 100 Stator Current
50 50
Is
0
Is
0
0 0.5 1 1.5 0 0.2 0.4 0.6 0.8 1 1.2 1.4
-50
-50
-100
-100 time,c time,c
a) Non optimum input parameter b) optimum input parameter for Te
Fig.2. Outputs of the simulated motor
The relation between each of the six input parameters and three output parameters are studied based on the simulated motor
(Fig.3). The effects on the three output parameters were obtained by changing the six parameters are changed one by one and
keeping rest constant. Normalized values of the six inputs and three outputs are between 0 and 1.
REFERENCES