Advanced
Optimal Design of Multi-DOF Deflection Type PM Motor by Response Surface Methodology
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
Copyright © 2015, The Korean Institute of Electrical Engineers
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
Author
Metrics
Cited by
TagCloud
About the Authors
Zheng Li
Corresponding Author: School of Electrical Engineering, Hebei University of Science and Technology, China. (Lzhfgd@163.com)
Lu Zhang
School of Electrical Engineering, Hebei University of Science and Technology, China. (jiuselu.fly@163.com)
Qingqing Lun
School of Electrical Engineering, Hebei University of Science and Technology, China. (jiuselu.fly@163.com)
Hongbo Jin
Hebei Science and Technology Information Institute, China. (jinhbo@163.com)

Abstract
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.
Keywords
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
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.
- 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 ( 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
PPT Slide
Lager Image
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
PPT Slide
Lager Image
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:
PPT Slide
Lager Image
where xi represent R r, H , σ , L respectively, T represents the maximum torque of the new multi-DOF deflection type PM motor, mxi ( 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
PPT Slide
Lager Image
Influence rate of optimization parameters on motor torque
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 ( 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:
PPT Slide
Lager Image
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:
PPT Slide
Lager Image
Based on least squares estimation, the vector of estimator b can be found as follow:
PPT Slide
Lager Image
The fitted regression model is:
PPT Slide
Lager Image
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:
PPT Slide
Lager Image
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.
PPT Slide
Lager Image
The CCD model
Analysis of variance
PPT Slide
Lager Image
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,
PPT Slide
Lager Image
represents the average value of test dates, and
PPT Slide
Lager Image
represents Predicted value of response surface model.
PPT Slide
Lager Image
The sum of squares due to regression (SSR) reflects the random error and system error. The formula is
PPT Slide
Lager Image
The sum of squares (SSE) unaccounted for by the fitted model is
PPT Slide
Lager Image
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
PPT Slide
Lager Image
The coefficient of determination R 2 is used to determine the precision of the regression equation from CCD. It is shown as
PPT Slide
Lager Image
- 3.2 Analysis of the results
- 3.2.1 Response surface model
The value of α 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
PPT Slide
Lager Image
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
PPT Slide
Lager Image
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)
PPT Slide
Lager Image
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.
- 3.2.2 The analysis of variance
From Table 7 , the model 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
PPT Slide
Lager Image
Analysis of variance
- 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
Response surface for effects of the rotor diameter and the height of rotor on response value
PPT Slide
Lager Image
Response surface for effects of the rotor diameter and the air gap length on response value
PPT Slide
Lager Image
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.
- 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
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.
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 ( 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.
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
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.
References
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