Optimal Design of Multi-DOF Deflection Type PM Motor by Response Surface Methodology

Journal of Electrical Engineering and Technology.
2015.
May,
10(3):
965-970

This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/)which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

- Received : May 03, 2014
- Published : May 01, 2015

Download

PDF

e-PUB

PubReader

PPT

Export by style

Article

Metrics

Cited by

TagCloud

This paper uses response surface methodology as the optimization method of torque of multi-DOF deflection type PM motor. Firstly, the application of Taguchi algorithm selects structural parameters affecting the motor torque largely which simplifies the optimization process greatly. Then, based on the central composite design (CCD), response surface equation numerical model is constructed by the finite element method. With the aid of experiment design and analysis software, the effects of the interaction among factors on the index are analyzed. The results show that the analytical method is efficient and reliable and the experimental results can be predicted by response surface functions.
Three-dimensional structure of the motor
By controlling the energization of the stator coils, the motor can generate a continuous rotation torque in order to achieve continuous motion, which is the core of the electrifying strategy.
R
r), height of the rotor (
H
), air gap length between stator and rotor (
σ
) and the distance between the double stator (
L
). Parameter values of the variables can be shown from
Table 1
.
Motor optimization parameters and parameter values of the variables
According to Taguchi algorithm, this paper selects orthogonal table
L
_{9}
(3
^{4}
) and establishes the experimental matrix, where
L
represents an orthogonal matrix, “9” represents the number of tests, “3” represents parameter variables number, “4” represents the number of parameters
[6]
. According to the experimental matrix, the maximum torque is got by finite element method, as shown in
Table 2
.
Orthogonal table and the results of finite element analysis
To get the influence rate of each parameter, Taguchi method calculates the average value of each parameter, and then analyzes the variance
[9]
. The influence rate of each optimization parameters calculated by the following formula:
where
xi
represent
R
r,
H
,
σ
,
L
respectively,
T
represents the maximum torque of the new multi-DOF deflection type PM motor,
m_{xi}
(
T
_{i}
) represents the average value of
T
.
Optimization of all the variables is impractical and may not be necessary for real applications
[10]
.
Table 3
shows that the influence rates of four parameters are 17.5%, 49.0%, 27.8%, 5.7% respectively, which can provide the basis for an economic solution.
Influence rate of optimization parameters on motor torque
R
r), height of the rotor (
H
) and air gap length between stator and rotor (
σ
). This paper uses multiple linear regression method to establish the response surface model. A second-order function is used to approximate Response Surface model, which can be described as:
where 𝛽
_{0}
, 𝛽
_{i}
, 𝛽
_{ii}
are estimated by a regression from experiments,
ε
is the error estimate
[11]
.
Response surface model in terms of the observation can be expressed by the following matrix:
Based on least squares estimation, the vector of estimator
b
can be found as follow:
The fitted regression model is:
It is convenient to transform the natural design variables to coded variables, which are usually defined to be dimensionless quantities with mean zero and the same spread or standard deviation
[12]
. The coding rule is described as follow:
where
x
_{i}
represent the permanent magnet rotor diameter (
R
r), height of the rotor (
H
) and air gap length between stator and rotor (
σ
).
The appropriate selection of the points for experiments can reduce the variance of the coefficients, which makes the response surface more reliable. The central composite design (CCD) is applied as the experimental design method to estimate the regression. Model of CCD consists of eight runs at the corners of a cubic (±1, ±1, ±1), plus six runs at the star points, (0,±
α
,0) , (0,0,±
α
) , and one run at the cubic centre (0, 0, 0). CCD model can be described as follow:
In order to verify the accuracy of the obtained fitting equation, it is necessary to analyze the variance of response surface model.
Table 4
is called an analysis-of-variance (ANOVA) table, where
N
presents total number of experiments;
p
presents the number of unknown parameters.
The CCD model
Analysis of variance
The total variation in a set of data is called the total sum of squares (
SST
). In the following formula (7)-(11),
Y
_{u}
represents response value,
represents the average value of test dates, and
represents Predicted value of response surface model.
The sum of squares due to regression (SSR) reflects the random error and system error. The formula is
The sum of squares (SSE) unaccounted for by the fitted model is
The ratio between mean square regression and mean square residual is
F
- statistic, which is used to test the computed statistics. It is define as
The coefficient of determination
R
^{2}
is used to determine the precision of the regression equation from CCD. It is shown as
α
is 1.682 to guarantee the rotation and robustness of the model. The design parameters are provided in
Table 5
.
Application of CCD designing parameter variables
By the finite element software, the maximum torque for different test sets are obtained, as shown in
Table 6
.
Experimental arrangement and its dates
Based on least squares estimation, equation of mathematical models of response surface for the motor torque can be shown in (12)
where
X
1,
X
2,
X
3 represent the permanent magnet rotor diameter (
Rr
), height of the rotor (
H
) and air gap between stator and rotor (
σ)
from CCD respectively.
P
value is less than 0.01 indicating a good fit of the model, and the lack of fit value
P
is equal to 0.0849 (greater than 0.05). Furthermore, the coefficient of determination
R
^{2}
is equal to 0.9939 (greater than 0.8000). The analysis of variance indicates that linear relationship between response value and variables is significant and the data regularity can be reflected by the response model accurately.
X
1,
X
2,
X
3,
X
1
X
2,
X
2
X
3,
X
1
^{2}
,
X
2
^{2}
,
X
3
^{2}
influence the motor torque greatly and the remaining items are not significant.
Analysis of variance
Response surface for effects of the rotor diameter and the height of rotor on response value
Response surface for effects of the rotor diameter and the air gap length on response value
Response surface for effects of the height of rotor and the air gap length on response value
Fig. 3
shows that when air gap length is 20 mm, the interaction between the rotor diameter and the height of rotor is remarkable. Compared with the rotor diameter, magnitude of response surface changes more greatly with increasing height of the rotor, which indicates that the height of rotor has bigger influence on response surface than rotor diameter. The bigger the rotor height is, the larger the toque of motor is. However, with the rotor diameter increasing, motor torque increases first and then decreases. The reason is that when the volume of motor is constant, with the increase of rotor diameter, the stator coil and the stator core of motor corresponding decrease, resulting in reducing flux. So, when the rotor diameter is at 15-19 mm, the response value of motor torque is bigger.
Fig. 4
indicates that when the rotor height is 10 mm, the interaction between the rotor diameter and the air gap length is not significant. It can be seen from the graph that the rotor diameter has bigger influence on response surface than the air gap length. With the increase of air gap length, the torque of multi-DOF deflection type PM motor is smaller. However, with the rotor diameter increasing, motor torque increases first and then decreases. When the rotor diameter is at 15-17.5 mm and air-gap length is at 0.70-1.35 mm, motor torque is bigger.
Fig. 5
shows that when the rotor diameter is at 10 mm, the interaction between the height of rotor and the air gap length is obvious. It can be seen from the graph that the rotor height has bigger influence on response surface than the air gap length. The bigger the rotor height is, the larger the toque of motor is. However, with the air gap length increasing, motor torque decreases accordingly.
By Design-Expert 8.0.5.0 software, when
Rr
=16.45 mm,
H
=14 mm,
σ
=1.18 mm, torque of the new multi-DOF deflection type PM motor reaches maximum value. While the motor air gap is larger, the magnetic flux leakage is greater, so the efficiency of the motor will be reduced. Taken together, the optimal motor structure is adjusted as:
Rr
=16.45 mm,
H
= 14 mm,
σ
=10.7 mm.
Experiment date for test validation
Table 8
shows that Relative errors between predicted values and measured values are less than 5%.The torque of improved motor is increased by 29.4%, improving the performance of the motor.
Rr
), height of the rotor (
H
) and air gap length between stator and rotor (
σ
). Based on the central composite design (CCD), response surface equation numerical model is constructed by the finite element method. The analysis of variance indicates that linear relationship between response value and variables is significant and the data regularity can be reflected by the response model accurately. With the aid of software Design-Expert 8.0.5.0, when
Rr
=16.45 mm,
H
=14 mm,
σ
=0.7 mm, torque of this motor reaches maximum value and is increased by 29.4% compared with the initial model. The results prove the effectiveness and correctness of the presented optimal design scheme, which provides the reference for further related optimal design of multi-DOF motors or actuators.
Zheng Li received his B.Sc. and Ph.D. degrees in electrical engineering and power electronics and electric drive from Hefei University of Technology, Hefei, China, in 2002 and 2007, respectively. Since 2007, he has been a Lecturer, Associate Professor and Professor with the School of Electrical Engineering, Hebei University of Science and Technology. He is the author of more than 70 published papers. His current research interests include design, analysis, and control of novel motors and actuators, intelligent control, and power electronics.
Lu Zhang received her B.Sc. degree in Control Technology and Instruments from the Agricultural University of Hebei, China. Her research interests are novel multi-DOF motor and its control technology.
Qingqing Lun received her B.Sc. degree in Electrical Engineering and Automation from Hebei University of Science and Technology, China. Her research interests are special PM motors’ design and optimization.
Hongbo Jin received her B.Sc. degree in Mechatronic Engineering from Central South University, China. Currently, he is the Dean and Senior Engineer of Department of Novelty Consulting Center, Hebei Science and Technology Information Institute. His research interests are special motor drive and control.

1. Introduction

With the development of technology, engineering calculation is more complex which takes more time. Many engineering problems are nonlinear, so the optimization of engineering calculation forced scholars to find new and reliable method of mathematical programming. Response surface methodology is a progressive approximation optimization method which can efficiently solve this nonlinear optimization problem.
The response surface method combines statistics and mathematics
[1]
. In recent years, due to the high speed development of statistics, the response surface method is no longer confined to the biological and chemical industry, in the ecology, food science and other new areas also show its strong advantages
[2]
.
Traditionally, multi-DOF motion is often completed by the combination of many single-DOF motion motors and mechanical transmission mechanisms. This not only results in low precision devices, but also decreases dynamic and static performance. Meanwhile, many of the applications of reduction gears increase nonlinear friction and energy losses
[3]
. Motors which have simple structure and can achieve multi-DOF motion become the best alternative means, therefore, research and application of multi-DOF motion of motor has attracted wide attention of scholars at home and abroad
[4
,
5]
.
This paper uses response surface methodology as the optimization method of torque of multi-DOF deflection type PM motor. Based on Taguchi algorithm, appropriate optimization parameters are screened out. Based on the central composite design (CCD), the response surface equation numerical model is constructed by the finite element method. By Design-Expert 8.0.5.0 software, the optimal motor structure is obtained. Results of the improved motor torque are analyzed and compared with the initial model.
2. Application of Taguchi Method

Traditionally, if four different parameters have to change and each parameter has three variations, 81 experiments needed to be performed
[6]
. The traditional methods are often time-consuming and laborious .The application of Taguchi algorithm selects structural parameters affecting the motor torque largely which simplifies the optimization process greatly.
- 2.1 Initial structure data of the motor

While achieving a continuous movement, the motor can complete the tilting movement within a certain angle
[7]
. The novel multi-DOF deflection type PM motor has a 6-pole permanent magnet rotor of the “drum” shape. The height of the rotor is 10mm and the outer diameter of the rotor is 31 mm, the inner diameter of rotor is 10 mm. Magnetic N pole and S pole are arranged alternately. The motor contains two layers of stator coils; each layer has 8 stator coils which can be controlled independently. The height of the stator is 17.5 mm and the outer diameter of the stator is 90 mm, the inner diameter of the stator is 78mm, the distance between the two stators is 5.5 mm. By giving different current to stator coils, the rotor and stator realize the multi-degree of freedom movement around a spherical bearing maintaining a constant gap (5.5 mm).The cross-section of this motor is shown in
Fig. 1
. For the PM motor, the electromagnetic torque generated by the interaction between the current-carrying coils and PM poles making the motor achieve the z-rotation, x-rotation and y-rotation
[8]
.
PPT Slide

Lager Image

- 2.2 Determining the optimal structure parameters of the motor

In fact, the purpose of optimizing the performance of the motor is to allow the motor to be further improved, for example, motor efficiency, power factor, torque performance, a minimum amount of material, etc. Multi-DOF deflection Type PM Motor as a new special motor, the purpose of the study is to allow the motor to achieve better multi-DOF motion, so the performance index of the motor efficiency and power factor is not the primary consideration. Therefore, the optimization object of this motor is torque characteristics. The ultimate aim is to pick out the optimal motor structure size for maximum torque in a certain space occupied. According to the structure characteristics of the motor, preliminary selection of optimization parameters variable number is 4. The optimization parameters are the permanent magnet rotor diameter (
Motor optimization parameters and parameter values of the variables

PPT Slide

Lager Image

Orthogonal table and the results of finite element analysis

PPT Slide

Lager Image

PPT Slide

Lager Image

Influence rate of optimization parameters on motor torque

PPT Slide

Lager Image

3. Experimental Design and Analysis

- 3.1 Response surface method

Since the motor is used for robot joints, eyes and other small space components, optimization on the motor under the condition of constant volume, has a vital significance. When the motor volume is constant, the torque of the motor is mainly related with the permanent magnet rotor diameter (
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

Analysis of variance

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- 3.2 Analysis of the results

- 3.2.1 Response surface model

The value of
Application of CCD designing parameter variables

PPT Slide

Lager Image

Experimental arrangement and its dates

PPT Slide

Lager Image

PPT Slide

Lager Image

- 3.2.2 The analysis of variance

From
Table 7
, the model
Analysis of variance

PPT Slide

Lager Image

- 3.2.3 Analysis of each factor interaction

Fig. 3
to
Fig. 5
describe the interaction between the various factors in the response surface model. From
Fig.3
and
Fig. 5
, large changes in the magnitude of the response surface can be visually seen which indicates that the interaction between the various factors influences the motor response significantly. In
Fig. 4
, response surface is smooth revealing that the influence between the rotor diameter and rotor diameter on response values is smaller.
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- 3.2.4 Validation of the numerical model

Two groups of representative parameters are selected to verify the reliability of the response surface model. Comparative results are shown in
Table 8
.
Experiment date for test validation

PPT Slide

Lager Image

4. Conclusion

This paper uses response surface methodology as the optimization method of torque of multi-DOF deflection type PM motor. Optimization parameters selected by Taguchi algorithm is the permanent magnet rotor diameter (
Acknowledgements

This work is supported by the National Natural Science Foundation of China (No. 51107031), the Natural Science Foundation of Hebei Province of China (No. E2014208134), the Foundation of China Scholarship Council (No. 2011813032).

BIO

Li J T
,
Liu Z J
,
Jabbar M A
2004
“Design optimization for cogging torque minimization using response surface methodology”
IEEE Transactions on Magnetics
40
(2)
1176 -
1179
** DOI : 10.1109/TMAG.2004.824809**

Vadde K K
,
Syrotiuk V R
,
Montgomery D C
2006
“Optimizing protocol interaction using response surface methodology”
IEEE Transactions on Mobile Computing
5
(6)
627 -
639
** DOI : 10.1109/TMC.2006.82**

Li Z.
2009
“Robust control of PM spherical stepper motor based on neural networks”
IEEE Transactions on Industrial Electronics
56
(8)
2945 -
2954

Lee J
,
Kim D
,
Baek S
2008
“Newly Structured Double Excited Two-Degree-of-Freedom Motor for Security Camera”
IEEE Transaction on Magnetics
44
(11)
4041 -
4044
** DOI : 10.1109/TMAG.2008.2002802**

Shohei I
,
Katsuhro H
,
She M
2010
“Proposal of Electromagnetic Spherical Actuator with 3-DOF”
Computation and Mathematics in Electrical and Electronic Engineering
29
(4)
994 -
1003
** DOI : 10.1108/03321641011044415**

Hwang K Y
,
Rhee S B
,
Yang B Y
2007
“Rotor pole design in spoke-type brushless DC motor by response surface method”
IEEE Transactions on Magnetics
43
(4)
1833 -
1836
** DOI : 10.1109/TMAG.2007.892616**

Li Z
2014
“Magnetic-field computation of a novel 3-DOF deflection-type PM motor with analytical and finite-element methods”
ELEKTROTEHNIŠKI VESTNIK
81
(1-2)
51 -
56

Li Z
,
Zhang Y
,
Wang Q
,
Wang X
2013
“Analytical modeling and analysis of magnetic field for a novel 3-DOF deflection type PM motor”
Proceedings of the CSEE
33
(1)
219 -
225

Kim S I
,
Hong J P
,
Kim Y K
2006
“Optimal design of slotless-type PMLSM considering multiple responses by response surface methodology”
IEEE Transactions on Magnetics
42
(4)
1219 -
1222
** DOI : 10.1109/TMAG.2006.871950**

Lee J Y
,
Kim S I
,
Hong J P
2006
“Optimal design of superconducting motor to improve power density using 3D EMCN and response surface methodology”
IEEE Transactions on Applied Superconductivity
16
(2)
1819 -
1822
** DOI : 10.1109/TASC.2005.869633**

Kim Y K
,
Hong J P
,
Hur J
2003
“Torque characteristic analysis considering the manufacturing tolerance for electric machine by stochastic response surface method”
IEEE Transactions on Industry Applications
39
(3)
713 -
719
** DOI : 10.1109/TIA.2003.811777**

Cheng Y
,
Xu G
,
Zhu D
2006
“Thermal analysis for indirect liquid cooled multichip module using computational fluid dynamic simulation and response surface methodology”
IEEE Transactions on Components and Packaging Technologies
29
(1)
39 -
46
** DOI : 10.1109/TCAPT.2005.848589**

Citing 'Optimal Design of Multi-DOF Deflection Type PM Motor by Response Surface Methodology
'

@article{ E1EEFQ_2015_v10n3_965}
,title={Optimal Design of Multi-DOF Deflection Type PM Motor by Response Surface Methodology}
,volume={3}
, url={http://dx.doi.org/10.5370/JEET.2015.10.3.965}, DOI={10.5370/JEET.2015.10.3.965}
, number= {3}
, journal={Journal of Electrical Engineering and Technology}
, publisher={The Korean Institute of Electrical Engineers}
, author={Li, Zheng
and
Zhang, Lu
and
Lun, Qingqing
and
Jin, Hongbo}
, year={2015}
, month={May}