Advanced
Optimal Design of Permanent Magnetic Actuator for Permanent Magnet Reduction and Dynamic Characteristic Improvement using Response Surface Methodology
Optimal Design of Permanent Magnetic Actuator for Permanent Magnet Reduction and Dynamic Characteristic Improvement using Response Surface Methodology
Journal of Electrical Engineering and Technology. 2015. May, 10(3): 935-943
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 : October 23, 2013
  • Accepted : November 12, 2014
  • Published : May 01, 2015
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Hyun-Mo Ahn
Dept. of Electrical Engineering, Dong-A University, Korea. (hmahn@donga.ac.kr)
Tae-Kyung Chung
School of Electrical and Electronics Engineering, Chung-Ang University, Korea. (tkchung@cau.ac.kr)
Yeon-Ho Oh
Power Apparatus Research Center, HVDC Research Division, Korea Electrotechnology Research Institute, Korea. (yhoh@keri.re.kr)
Ki-Dong Song
Power Apparatus Research Center, HVDC Research Division, Korea Electrotechnology Research Institute, Korea. (yhoh@keri.re.kr)
Young-Il Kim
ENTEC Electric & Electronic Co., Ltd., Korea. (hrkoh@entecene. co.kr)
Heung-Ryeol Kho
ENTEC Electric & Electronic Co., Ltd., Korea. (hrkoh@entecene. co.kr)
Myeong-Seob Choi
ENTEC Electric & Electronic Co., Ltd., Korea. (hrkoh@entecene. co.kr)
Sung-Chin Hahn
Corresponding Author: Dept. of Electrical Engineering, Dong-A University, Korea. (schahn@dau.ac.kr)

Abstract
Permanent magnetic actuators (P.M.A.s) are widely used to drive medium-voltage-class vacuum circuit breakers (V.C.B.s). In this paper, a method for design optimization of a P.M.A. for V.C.B.s is discussed. An optimal design process employing the response surface method (R.S.M.) is proposed. In order to calculate electromagnetic and mechanical dynamic characteristics, an initial P.M.A. model is subjected to numerical analysis using finite element analysis (F.E.A.), which is validated by comparing the calculated dynamic characteristics of the initial P.M.A. model with no-load test results. Using tables of mixed orthogonal arrays and the R.S.M., the initial P.M.A. model is optimized to minimize the weight of the permanent magnet (P.M.) and to improve the dynamic characteristics. Finally, the dynamic characteristics of the optimally designed P.M.A. are compared to those of the initially designed P.M.A.
Keywords
1. Introduction
In recent years, the growing prospect of the electric power industry in a green and eco-friendly manner has gained much attention. Vacuum circuit breakers (V.C.B.s), which use vacuum as an insulation medium, are replacing gas circuit breakers (G.C.B.s) for use in medium-voltage-class power interception [1 - 3] .
Switch actuators for V.C.B.s typically use spring, solenoid, or permanent magnetic actuator (P.M.A.) mechanisms, and medium-voltage-class V.C.B.s have been using spring-driven mechanism for decades. Such mechanisms, however, have drawbacks including complicated structures, reliance on periodic maintenance, and proneness to oil leaks. P.M.A.s, by contrast, have advantages in terms of simpler structure, improved operating time, easy motion control, and high reliability. Accordingly, manufacturers seeking to reduce the V.C.B. size and to improve the operating time have begun to prefer P.M.A. over spring-driven mechanism [4 , 5] .
Fig. 1 shows the structure of the P.M.A. developed for use in a V.C.B. It has four important electromechanical components: the permanent magnets (P.M.s), a movable plunger, trip coils, and static yokes. When an excitation current is supplied through the trip coil, the plunger moves upwards or downwards. The mechanical energy generated due to movement of plunger is transmitted to a vacuum interrupter (V.I.). When the plunger reaches fully open or closed position, it stops to be held still by P.M.s only. The P.M.s, which together with the trip coil compose the P.M.A., are able to generate rapid continuous motion with more efficiency than a spring-type actuator [6 - 8] .
PPT Slide
Lager Image
Structure of P.M.A. model
The response surface method (R.S.M.) is an effective statistical technique that can be used to improve the performance of electrical devices, enabling the easy development of objective functions and the reduction of computing time [9 , 10] .
As most types of electric machines have many design parameters, it is often necessary to search extensively in order to determine which particular variable most affects relevant design results. Because the R.S.M. is generally used to optimize two or three design variables [11] , design of experiments (D.O.E.s) for seven design variables can be developed using tables of mixed orthogonal arrays to assess large numbers of design variables and their levels effectively [12] , with a response value for each parameter combination calculated by means of finite element analysis (F.E.A.). Using the R.S.M., the D.O.E.s according to tables of mixed orthogonal arrays are analyzed. A Pareto chart can be used in order to find the design variables with the largest effects on the response function [13] . The P.M.A. design optimization process proposed here is aimed at minimizing the weight of the P.M. at improving the dynamic characteristics of an initially designed P.M.A.
This paper describes a dynamic characteristic analysis of a P.M.A. using coupled finite element method (F.E.M.) and the design optimization of the P.M.A. using tables of mixed orthogonal arrays and the R.S.M. Using coupled F.E.M., electromagnetic and mechanical characteristics such as the holding force and the action complete time are calculated and used as response values for P.M.A. design optimization. In order to verify the validity of the F.E.A., the calculated dynamic characteristics of an initial P.M.A. model are compared with no-load test results. Finally, the optimal P.M.A. model designed by using the R.S.M. is compared with the initial P.M.A. model in order to verify the validity of the optimally designed model.
2. Characteristic of V.C.B. Mechanism
The force characteristics of the actuators used for driving a vacuum interrupter (V.I.) are shown in Fig. 2 [6 , 14] . It can be seen that the V.I. requires a high level of thrust at the end of stroke. Although the spring actuator produces high levels of thrust in the initial operation, this drops by the end of the stroke; in addition, spring actuators have complicated structures and face limitations to continuous operation. Although the solenoid actuator generates a similar force profile as that required by the V.I., this doesn’t reach the required force of V.I. at the end of the stroke [15] . By contrast, the P.M.A. matches the force characteristics of the V.I. closely during the initial operation of the stroke and also produces high levels of thrust at the end of the stroke, making it a more efficient V.I. driver than either the spring actuator or the solenoid.
PPT Slide
Lager Image
Force characteristic of V.C.B. mechanism
3. Numerical Analysis Model
As shown in Fig. 3 , the P.M.A. model used in this study consists of trip coils, plunger, yokes, and P.M.s, which are placed in the center of the P.M.A. between the yoke and the plunger. The holding force generated by the P.M.s maintains the position of the plunger. At that time, when an excitation current flows in the upper or lower coil, the P.M.A. performs an open or close operation. The nonlinear B-H characteristic of plunger and yoke is shown in Fig. 4 .
PPT Slide
Lager Image
P.M.A. Model (a) side-view (b) front-view
PPT Slide
Lager Image
B-H curve of M-19 steel
4. Numerical Analysis and Experimental Validation
Generally, a P.M.A. can be adequately represented in a mathematical manner through the following three equations: (a) the electric circuit equation with voltage source, (b) the magnetic field equation for magnetic flux and energy, and (c) the equation of motion for velocity and force [5 , 8 , 16] .
- 4.1 Governing equation
To analyze the dynamic characteristics of the P.M.A. model, the electric circuit equation should be coupled to the magnetic field equation. When a source voltage is charged to the capacitor for the driving source, the electric circuit equation can be expressed as
PPT Slide
Lager Image
where, C is the capacitance of capacitor for the driving source [F], i ( t ) is the excitation current [A], R is the resistance [Ω], and L is the inductance [H]. The current density generated by the source capacitor and the P.M. are fundamental to determining the dynamics of the drive plunger, and the magnetic field equation for the P.M.A. is formulated as
PPT Slide
Lager Image
where ν is the magnetic reluctivity,
PPT Slide
Lager Image
is the magnetic vector potential, σ is the conductivity [S/m],
PPT Slide
Lager Image
is the current density [A/m 2 ], V is the electric potential, μ is the permeability, and
PPT Slide
Lager Image
is the magnetization vector [A/m]. In terms of the P.M.A. model, the most important parameter from the magnetic field equation above is the magnetic energy, which operates the plunger. The magnetic force can be calculated using the principle of virtual work and is given by
PPT Slide
Lager Image
where Fmag is the magnetic force [N], z is the motional direction [m], and v is the volume [m 3 ]. This calculated magnetic force is used to build an equation of motion from which plunger acceleration, velocity, and position can be extracted; this is expressed as
PPT Slide
Lager Image
where M is the plunger mass [kg], 2 2 ( d 2 z / dt 2 ) is its acceleration in the direction of motion [m/s 2 ], B is the attenuation constant [N/m/s], ( dz / dt ) is the plunger velocity in the direction of motion [m/s], K is the spring constant [N/m], and f is the frictional force [N].
- 4.2 Comparison of simulation and experiment
In order to determine the holding force, it is necessary to understand the distribution of magnetic flux density in the absence of excitation current. Using F.E.M., the magnetic flux density distribution can be numerically calculated, as shown in Fig. 5 . In this figure, most of the magnetic flux flow occurs on the lower part of the P.M.A. model because the upper part has relatively high magnetic reluctance due to the air gap between the plunger and the upper yoke. With the plunger in the position show, the holding force calculated by F.E.A. is 7777.2 [N], which agrees closely with the value obtained from no-load test.
PPT Slide
Lager Image
Magnetic flux density at no-load condition
In order to accurately analyze the dynamic characteristics of the P.M.A., an understanding of the electromagnetic field-circuit link is necessary. To analyze this, an external electric circuit was coupled to the electromagnetic-circuit equation, as show in Fig. 6 .
PPT Slide
Lager Image
Drive system
As seen in Fig. 7 , the excitation current calculated from simulation is compared with experimental value of no-load test from t=0 to t=70 [ms] and the maximum excitation current calculated by simulation is 107.3 [A], which represents a pattern similar to the current wave obtained from the no-load test. Dynamic characteristic of the plunger was calculated using a coupled electromagnetic-mechanical F.E.A. Fig. 8 shows a comparison of the operating time between the experimental and the F.E.A. values. The plunger action complete time calculated by F.E.A. is 54.8 [ms], which shows good agreement with the experimental value obtained from no-load test.
PPT Slide
Lager Image
Comparison of excitation current of simulation and no-load test
PPT Slide
Lager Image
Comparison of operating time of simulation and noload test
The dynamic characteristics of initial P.M.A. by F.E.A. are compared with the experimental results measured from no-load test and are shown in Table 1 . This table summarizes the close agreement of results discussed above. The P.M.A. design optimization described in this paper is performed based on the results of F.E.A.
Dynamic characteristics of initial P.M.A. model
PPT Slide
Lager Image
Dynamic characteristics of initial P.M.A. model
5. Design Optimization
- 5.1 Response surface method[17]
R.S.M., an effective statistical technique for improving electrical device performance, is used for design optimization of the P.M.A. model. R.S.M. uses a statistical method to develop an approximate response formula expressing the correlation between the design variables and the response variables [10] . Schematically, the response η to k design variables can be expressed as
PPT Slide
Lager Image
If a secondary regression model is used for function approximation, this equation can be re-formulated as follows [10 , 18 , 19] :
PPT Slide
Lager Image
where, β represents the regression coefficients, the indexed x variables are the design variables, ε is the random error, and k is the number of design variables. The least-squares method (L.S.M.) can be used to predict the unknown coefficients. In matrix notation of the fitted coefficients and response model can be expressed as follows [11] :
PPT Slide
Lager Image
where,
PPT Slide
Lager Image
is the vector containing the unknown coefficients. To optimize the design of initial P.M.A. model, we propose the optimal design process shown in Fig. 9 .
PPT Slide
Lager Image
Design optimization process
- 5.2 Design variables and levels
In R.S.M., the selection of design variables is very important [9] ; however, as most electric machines have a large number of design parameter, D.O.E.s of seven design variables made by using tables of mixed orthogonal arrays is used in order to perform an effective analysis in this study. The seven dimensions selected as design variables are shown in Fig. 10 . The design variables have three levels, and these are all shown in Table 2 in which second level of design variables contains values from initial P.M.A. model.
PPT Slide
Lager Image
P.M.A. design variables
Design variables and levels
PPT Slide
Lager Image
Design variables and levels
- 5.3 Optimal design result
The mixed orthogonal arrays of D.O.E.s shown in Table 3 , which represent design variables considered at each level and response values obtained from F.E.A., are used to perform an R.S.M. design optimization of the P.M.A. model. The response functions are the holding force, maximum thrust, action complete time, and P.M. weight, and the response values for each parameter combination are calculated by F.E.A.
Table of mixed orthogonal array L36(21x37)
PPT Slide
Lager Image
Table of mixed orthogonal array L36(21x37)
In order to visualize the impact of each design variable on the response function, the main effects of each function are represented with the analysis of variance (A.N.O.V.A.), as shown in Fig. 11 .
PPT Slide
Lager Image
Main effect analysis for (a) holding force; (b) max thrust force; (c) complete time; (d) P.M. weight
It can be seen that the holding force and the max thrust are maximized at dv 3 (1), dv 5 (1), and dv 7 (1). Conversely, the action complete time is minimized at dv 3 (1), dv 5 (1), and dv 7 (1). The P.M. weight is minimized at dv 1 (-1), dv 2 (-1), and dv 7 (-1).
In order to draw a response surface based on the experimental data from Table 3 , the response function is extracted by the R.S.M. in accordance with the aim of this study to minimize the weight of the P.M. and to improve the dynamic characteristics of an initially designed P.M.A. In other words, the proposed P.M.A. optimal design process attempts to improve the action complete time (y T ) while minimizing PM weight (y W ) maintaining the holding force (y HF ) and maximum thrust (y TF ). The two fitted second-order polynomial response functions of seven design variables that result from this can be expressed as
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
The adjusted coefficients of determination of the holding force, maximum thrust, action complete time, and P.M. weight are 99.9 [%], 99.2 [%], 95.1 [%], and 100 [%], respectively, and the response surfaces for each response function are expressed using (8)-(11), as shown in Fig. 12 , which shows the contour lines and optimal points for each design variable.
PPT Slide
Lager Image
Contour lines and optimal points around each response function
In order to investigate the influence of significant parameters on the design results, the magnitude and importance of the impact of each design variable on the response functions is analyzed by means of the Pareto chart shown in Fig. 13 . It can be seen that the design variable dv 7 is most significant for the holding force, maximum thrust, and action complete time. Because design variables such as dv 1 , dv 2 , and dv 7 represent the dimensions of the P.M., these three variables come to the effective design variables determining P.M. weight. The Pareto chart shows absolute value of the effects of each design variable. The optimal P.M.A. model obtained from R.S.M. is shown in Table 4 .
PPT Slide
Lager Image
Pareto chart of the standardized effects for each response function.
Optimal design results for P.M.A. model
PPT Slide
Lager Image
Optimal design results for P.M.A. model
Dynamic characteristics such as the holding force, maximum thrust, action complete time, and P.M. weight are compared between the initial and optimal models in Table 5 .
Comparison of initial model and optimal model
PPT Slide
Lager Image
Comparison of initial model and optimal model
As shown in Table 5 , the action complete time has been shortened from 54.8 [ms] to 54.2 [ms], and the P.M. weight decreases by 16.7 [%]. Even though the action complete time of the P.M.A. becomes a little faster, we think that the performance of P.M.A. is improved by the proposed method.
6. Conclusion
In this paper, a design optimization process to minimize P.M. weight and to improve the dynamic characteristics of an initially designed P.M.A. model was proposed. In order to verify the validity of numerical analysis based on coupled F.E.M. used in this process, the dynamic characteristics of P.M.A. were compared with experimental results from no-load test and show good agreements.
In order to optimize the initial P.M.A. model, an optimal design process was performed by the proposed R.S.M. using the mixed orthogonal arrays table based on seven design variables and their three levels. In order to verify the validity of the proposed optimizing process, the optimal P.M.A. model was then compared with the initial P.M.A. model. As the optimal designed P.M.A. model shows improved performance relative to the initial model, we conclude that the proposed approach is efficient to improve the performance of the optimal designed P.M.A. and to reduce the number of experiments by the proposed R.S.M. using the mixed orthogonal array table. We therefore believe that this process will be useful in the design optimization of power apparatuses.
Acknowledgements
This work was supported by the Human Resources Development program (No.20134010200550) and the Global Excellent Technology Innovation (No. 10043245) of the Korea Institute of Energy Technology Evaluation and Planning (KETEP) grant funded by the Korea government Ministry of Trade, Industry and Energy.
BIO
Hyun-Mo Ahn He received the B.S. and M.S. degrees in electrical engineering from Dong-A University, Busan, Korea, in 2009 and 2011, respectively. He is currently working toward the Ph.D. degree. His research interests are multi-physics analysis and numerical analysis of power apparatus and electric machines.
Tae-Kyung Chung He received the B.S., M.S., and Ph.D. degrees in electrical engineering from Seoul National University in 1981, 1983, and 1987, respectively. From 1986 to 1987, he worked at Daewoo Heavy Industries as a senior researcher. Since 1988, he has been with the School of Electrical and Electronics Engineering at Chung-Ang University in Korea as a professor. In 1991, he was with Harvey Mudd College in California as a visiting scholar. His main research fields are design of electric motors and actuators.
Yeon-Ho Oh He received the B.S. and M.S. degrees in electrical engineering from Dong-A University, Busan, Korea, in 1991 and 1993, respectively. He is currently a Research Engineer with the Power Apparatus Research Center, HVDC Research Division, Korea Electrotechnology Research Institute (KERI). He is now Team Leader for the estimation of the power apparatus with numerical analysis.
Ki-Dong Song He received the B.S. and M.S. degrees in electrical engineering from InHa University, Incheon, Korea, in 1988 and 1990, respectively. He got his Ph.D. degree in electrical engineering from Busan National University in 2003. He is currently a Research Engineer with the Power Apparatus Research Center, HVDC Research Division, Korea Electrotechnology Research Institute (KERI). His research interests focus on analysis and design for AC Power Devices and HVDC circuit breakers.
Young-Il Kim He received the B.S. degrees in mechanical engineering from Korea Aerospace University, Goyang, Korea, in 1990. He is currently an Executive Director of Power Apparatus R & D Institute in ENTEC Electric & Electronic Co., Ltd.
Heung-Ryeol Koh He received the B.S. degrees in electrical engineering from Hanyang University, Seoul, Korea, in 1990. He is currently a Senior Researcher of Power Apparatus R & D Institute in ENTEC Electric & Electronic Co., Ltd.
Myeong-Seob Choi He received the B.S. degrees in physics and M.S. degrees in electrical engineering from Soonchunhyang University, Asan, Korea, in 2002 and 2004, respectively. He is currently a Senior Researcher of Power Apparatus R & D Institute in ENTEC Electric & Electronic Co., Ltd.
Sung-Chin Hahn He received the B.S., M.S., and Ph.D. degrees in electrical engineering from Seoul National University, Seoul, Korea. He is currently a Professor with the Department of Electrical Engineering, Dong-A University, Busan, Korea, and senior member of the Korean Institute of Electrical Engineers (KIEE). His current research covers multi-physics analysis and design of power apparatus and electric machines.
References
Shipp D. D. , Dionise T. J. , Lorch V. , MacFarlane B. G. 2011 “Transformer Failure Due to Circuit-Breaker-Induced Switching Transients” IEEE Trans. Ind. Appl. 47 (2) 707 - 718    DOI : 10.1109/TIA.2010.2101996
Shipp D. D. , Dionise T. J. , Lorch V. , MacFarlane W. G. 2012 “Vacuum Circuit Breaker Transients During Switching of an LMF Transformer” IEEE Trans. Ind. Appl. 48 (1) 37 - 44    DOI : 10.1109/TIA.2011.2175430
Iturregi A. , Torres E. , Zamora I. , Abarrategui O. “High Voltage Circuit Breaker: SF6 vs. Vacuum” Proc. the 7nd Int. Conf. Renewable Energies and Power Quality 2009 1 - 6
Park K. N. , Son J. W. , Hahn S. C. “Dynamic Characteristic Analysis of Permanent Magnetic Actuators Coupled Electromagnetic-mechanical Finite Element Method” Proc. 2010 International Conference on Electrical Machines and Systems 1706 - 1709
Fang S. , Lin H. , Ho S.L. 2009 “Transient Co-Simulation of Low Voltage Circuit Breaker with Permanent Magnet Actuator” IEEE Trans. Mag. 45 (3) 1242 - 1245    DOI : 10.1109/TMAG.2009.2012579
Hong S. K. , Ro J. S. , Jung H. K. 2014 “Optimal Design of a Novel Permanent Magnetic Actuator using Evolutionary Strategy Algorithm and Kriging Meta-model” J Electr. Eng. Technol. 9 (2) 471 - 477    DOI : 10.5370/JEET.2014.9.2.471
Fang S. , Lin H. , Ho S. L. , Wang X. , Jin P. , Liu H. 2009 “Characteristics Analysis and Simulation of Permanent Magnet Actuator with a New Control Method for Air Circuit Breaker” IEEE Trans. Mag. 45 (10) 4566 - 4569    DOI : 10.1109/TMAG.2009.2024895
Brauer J. R. 2006 Magnetic actuators and sensors Wiley-Interscience
Hong D. K. , Woo B. C. , Koo D. H. , Kang D. H. 2008 “Optimum Design of Transverse Flux Linear Motor for Weight Reduction and Improvement Thrust Force Using Response Surface Methodology” IEEE Trans. Mag. 44 (11) 4317 - 4320    DOI : 10.1109/TMAG.2008.2002474
Khuri A. , Mukhopadhyay S. 2010 “Response Surface Methodology: Advanced Review” WIREs Comp. Stat. 2 (2) 128 - 149    DOI : 10.1002/wics.73
Hong D. K. , Woo B. C. , Ahn C. W. 2010 “Optimum Design for Improvement of PM-type Longitudinal Flux Linear Motor Using the Statistical Methods” International Journal of Modern Physics B 24 (15n16) 2821 - 2826    DOI : 10.1142/S0217979210065696
Omekanda A. M. 2006 “Robust Torque and Torque-perinertia Optimization of a Switched Reluctance Motor using the Taguchi method” IEEE Trans. Ind. Appl. 42 (2) 473 - 478    DOI : 10.1109/TIA.2006.870031
Hong D. K. , Woo B. C. , Koo D. H , Lee K. C. 2010 “Electromagnet Weight Reduction in a Magnetic Levitation System for Contactless Delivery Applications” Sensors 10 (7) 6718 - 6729    DOI : 10.3390/s100706718
Slade P. G. 2008 The vacuum interrupter - theory, design, and application CRC Press
Khan S.H. , Cai M. , Grattan T.V. , Kajan K. , Honeywood M. , Mills S. 2007 “Computation of 3-D Magnetic Field Distribution in Long-Lifetime Electromagnetic Actuators” IEEE Trans. Mag. 43 (4) 1161 - 1164    DOI : 10.1109/TMAG.2007.892277
Ro J. S. , Hong S. K. , Jung H. K. 2013 “Characteristic Analysis and Design of a Novel Permanent Magnetic Actuator for a Vacuum Circuit Breaker” Electric Power Applications, IET 7 (2) 87 - 96    DOI : 10.1049/iet-epa.2012.0176
Box G. E. P. , Hunter J. S. , Hunter J. S. 2005 Statistics for Experimenters: Design Innovation, and Discovery Wiley-Interscience
Hasanien H. M. , Abd-Rabou A. S. , Sakr S. M. 2010 “Design Optimization of Transverse Flux Linear Motor for Weight Reduction and Performance Improvement Using Response Surface Methodology and Genetic Algorithms” IEEE Trans. Energy Convers. 25 (3) 598 - 605    DOI : 10.1109/TEC.2010.2050591
Xie J. , Kang D. H. , Woo B. C. , Lee J. Y. , Sha Z. H. , Zhao S. D. 2012 “Optimum Design of Transverse Flux Machine for High Contribution of Permanent Magnet to Torque Using Response Surface Methodology” J Electr. Eng. Technol. 7 (5) 745 - 752    DOI : 10.5370/JEET.2012.7.5.745