Advanced
3D Magnetic Analysis of Permanent Magnets in Spherical Configuration
3D Magnetic Analysis of Permanent Magnets in Spherical Configuration
Journal of Electrical Engineering and Technology. 2016. Jan, 11(1): 93-99
Copyright © 2016, 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 : August 07, 2014
  • Accepted : July 19, 2015
  • Published : January 01, 2016
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Yusuf Oner
Corresponding Author: Dept. of Electrical and Electronic Engineering, Pamukkale University, Turkey. (yoner@pau.edu.tr)
Selami Kesler
Dept. of Electrical and Electronic Engineering, Pamukkale University, Turkey. (skesler@pau.edu.tr)

Abstract
The present study aims to increase the amount of surface flux by changing the magnetic directions of a spherical magnet (NdFeB) consisting of four poles. For this purpose, the magnetic directions of quartile spherical slices constituting the spherical magnet are manipulated and their three-dimensional analyses are carried out by using finite-element method via Maxwell environment. The analysis of the magnetic quartile spheres with four different magnetic directions are compared to the each other, and then the quartile sphere with the best surface flux distribution is suggested for rotor structure. It is clear emphasized that the induced torque of the spherical motor, in which such a rotor is used, will be improved as well.
Keywords
1. Introduction
Permanent magnets have unique abilities over conventional electromagnets: they can continuously generate magnetic flux without consumption of input energy and no power for cooling. Especially, the neodymium-iron-boron (NdFeB) magnet developed in 1980s has enlarged the application of permanent magnets in the electrical industry [1 - 2] . The Halbach magnet array was proposed by Halbach in 1980 for undulatory or electron wiggler [3] . If there is a magnet for which the magnetization direction rotates continuously, magnetic flux is constructed via superposition on one side of the magnet while the flux is destructed on the other side. Since continuously rotating magnetization is difficult to manufacture, the ideal array is approximated with a segmented magnet [4] . Bojan Stumberger et al. point out the comparison of torque capability of the three-phase permanent magnet synchronous motors with different permanent magnet arrangements [5] . Moon G. Lee and Dae-Gab Gweon emphasized the multi-segmented trapezoidal magnet array for a high-precision positioning system [6] . More and more analytical approaches were proposed by many authors for calculating the magnetic field created by arc-shaped permanent magnets whose polarizations can be radial [7 - 9]
The design of the spherical magnet can also be the permanent magnet array. Actually, the spherical permanent magnet array (SPMA) is generally used in spherical motors. Spherical motors can be described as motors that have movement capability that can be possessed by one or more than one motor. Unlike other electric motors, their movement capability is not mono-axial; they can move along different axes. Due to this characteristic, they are chosen in industrial for different purposed applications. However, they have not achieved commercial success yet because they do not have a sound structure, and they have a complex electric and magnetic design.
A lot of research about spherical motors have been going on for several decades. There are many different structured spherical motors resembling induction motors [10 , 11] , direct-current motors [12 - 14] , stepper motors [15] , variable reluctance motors [16 , 17] and ultrasonic motors [18 , 19] . Three dimensional nature of the electromagnetic field distribution in the foregoing spherical motors make their electromagnetic and dynamic behavior difficult to analyze, and this has been a significant obstacle to their design optimization and servo application. As a result, the potential benefits of employing spherical motors have not been realized. Therefore, most of these motor forms remain un-commercialized. In recent past, new formulas of the spherical actuators moving in two-free axes are analyzed and approved with the help of experiences [20 , 21] .
And the features of rotors used in these designs are very different from each other. The directions of the magnets in the form of pieces used in a spherical motor to be designed should be selected in such a way that the surfaces at flux values can evenly be distributed [22 - 24] . Hongfeng Li and et al carried out the spherical harmonic analysis of a halbach array permanent magnet spherical motor [25] . Design and analysis of a variable arc permanent magnet array for the spherical motor were reported by Changliang Xia and et al. [26] . Hyun-Jong Park and et al introduced a several studies about performance increasing of permanent magnet spherical motor. [27 - 29] . Bin Li and et al implemented a spherical motor design by using spherical harmonics [30] . Zhe Qian and et al studied on analysis of permanent magnet spherical motor with cylindrical poles and reported effects of that pole structure [31] . Hyung-Woo Lee et al focused on calculation of the produced torque in permanent magnet spherical motors [32] . Guo Jinjun et al investigated the effect of pole layouts of permanent magnets and electromagnets on inclination torque of a multi degrees of freedom spherical motor [33] .
As known, the directions of permanent NdFeB magnets are determined during production, so they possess predetermined magnetic characteristics. So, in this study, the permanent magnet spherical rotor structures are investigated. For this purpose, magnetic analyses of this rotor type are made via simulations by changing magnetic directions of the quartile spherical slices forming the rotor to find out the increase in the surface flux values.
For each rotor structure, analysis is carried out for the permanent magnet material with the same features. In this way the possibility of making comparison is created.
2. Rotor Structure and Magnetic Directions
The spherical magnet introduced in this study consists of combining four-quartile spheres. Every quartile includes the same NdFeB 35 material as a single piece. The permanence value of the NdFeB 35 permanent magnet is Br=1,234 T, and its relative permeability is μr=1,15 [34] . Fig. 1. shows the spherical rotor structure designed with four pieces of permanent magnets, and the B-H curve for the spherical magnet can be seen in Fig. 2 .
PPT Slide
Lager Image
Permanent magnet rotor in spherical structure and Ansoft Maxwell software has been used for the proposed design and simulations [35].
PPT Slide
Lager Image
B-H curve and material specifications of the permanent magnet supported by the manufacturer.
The magnetic directions of the four pieces of quartile spheres forming the spherical magnet are changed in four different ways, and their magnetic analyses are conducted for each changing.
Every variation of the magnetic-field directions results in a change of flux directions inside and on the surface of the spherical magnet, and in this way, the magnetic direction yielding the highest flux value on the surface of the sphere is achieved for the rotor using in the application.
Br (1,234 T) is the value computed in the conditions of a closed loop (magnetic core) magnetic circuit [34 , 36] . This value was calculated to be around 0.6 T on the surface of the magnet (without magnetic core).
In numerical calculations made for the geometry of the spherical motor designed by Howe et al., the flux values on the surface of the spherical rotor were calculated with legendry polynomial. Since the same structure is used, the same calculations are directly used here
3. 3D Magnetic Analyses
3D magnetic analyses are carried out for 4-different modeled spherical permanent magnets. The magnetic dimension of Model-1 is adjusted in such a way as to pass through the middle of the quartile magnetic slice and at α=45°. Two opposite poles have N magnetic direction, and the other two opposite poles have S magnetic direction.spherical coordinate system
Model-2 is adjusted in such a way as to have a magnetic direction of α =90°; the magnetic dimension of Model-3 is adjusted to be α =0°; and the magnetic direction of Model-4 is adjusted to be α =180°. In Fig. 3. the vector images of the flux densities on the surface of the sphere can be seen. While one pole is directed outside, the other one is directed from outside to inside.
PPT Slide
Lager Image
The vector images of the flux intensity of the designed 4-different model spheres are: (a) Model-1 magnetic direction of the magnet α =45°; (b) Model-2 magnetic direction of the magnet α =90°; (c) Model-3 magnetic direction of the magnet α =0°; (d) Model-4 magnetic direction of the magnet α =180°.
At the end of the analyses carried out, when a section is taken from the center of the spherical magnet and then its magnetic flux lines are examined, many different magnetic distributions will be seen, as in Fig. 4 . The purpose of the present study is to create a high magnetic flux density and even flux distribution on the surface of the sphere. In Fig. 4. the magnetic flux distribution occurring on the central area is shown in the form of vector and power lines. By comparing these different analyses, information of the strength of the flux on the surface of the magnet can be obtained.
PPT Slide
Lager Image
(a) Model-1; (b) Model-2; (c) Model-3; (d) Model-4, magnetic flux density on the central section, vector image of the flux and distribution of the flux lines inside the sphere and on the surface of it, respectively.
In line with the architecture of the spherical motor, the external surface of the spherical magnet that can be used should be one pole, and the internal side should be one pole. Therefore, the best flux distribution can be seen in Model-1. Among the other models, Model-1 can be selected for spherical rotor application as it yields the highest level of flux and the smoothest flux. In Model-4, there is no good flux distribution on the surface in spite of high flux density. Therefore, structure of Model-1 which has low harmonic order and good flux distribution can be employed as a rotor for the spherical motors.
Magnetic poles occurring on the surfaces of the magnetic spheres depending on their directions can clearly be seen in Fig. 5 . Here, the sections marked with red represent the N pole, and the sections marked with blue represent the S pole. When we look at the Fig. 5. attached to Model-1, we can see the inner edge become the N pole and the external surface becomes its opposite the S pole. In this way, the external surface of the quartile sphere slice creates a pole extending from the center towards the outside. This is very suitable for the spherical motor architecture and completes movement of the spherical rotor.
PPT Slide
Lager Image
Polarization of the quartile spherical slice in (a) Model-1; (b) Model-2; (c) Model-3; (d) Model-4
4. Calculation of Magnetic Flux Distributions
The magnetic field region generated by the permanent magnet spherical rotor lays on the basis of distribution information, design optimization and the creation of the proper model for the spherical motor [20] . Given the fact that magnetic flux lines extend on the soft iron stator without leakage, the magnetic field region can be divided into two sub-areas. Air gap having μ0 magnetic permeability and coil region which is the permanent magnet surface having μ0 μr magnetic permeability. Therefore,
PPT Slide
Lager Image
Here μr is the magnet relative permeability and M is the magnetization vector based on residue magnetization. The μr value of a permanent magnet having linear demagnetization characteristics is constant and there is a M = Br / μ0 relationship between the magnetization vector, M and the permanent magnetism, Br .
The radial constituent of the magnetization vector can be expressed as follows by extending in the form of Mr spherical harmonics:
PPT Slide
Lager Image
where,
PPT Slide
Lager Image
m is associated with ordinary l . degree of Legendre polynomial and Mlm can be expressed as follows [20] .
PPT Slide
Lager Image
PPT Slide
Lager Image
By considering the radial constituent boundary conditions of the distribution, the flux density on the spherical permanent magnet surface can be calculated by Eq. (5).
PPT Slide
Lager Image
where,
PPT Slide
Lager Image
And, Rm is the radius of the spherical rotor.
5. Results and Discussion
The magnetic flux densities of the rotor are calculated by Eq. (5) for the position of Rm =30 mm, r = 0.031 m and θ =0° for the gap between α =- π and + π . For the radial constituents of the magnetic flux densities, a comparison of the results between calculated by analytical method and obtained by finite elements method is given in Fig. 6 , for designed model-1.
PPT Slide
Lager Image
Comparison of the results for the radial constituents of magnetic flux densities for the θ=0° position of the rotor for the gap between α=-π and +π calculated by analytical method and finite elements method.
In order to measure the fluxes on the spherical magnet surfaces, the mechanism illustrated in Fig. 7 is established, and surface fluxes of the spherical magnet can be depicted graphically. The spherical magnet is mounted to a step motor shaft and rotated step by step by the computer-aided program. At each step, it waits for 2 seconds, and then it rotates one more step. During the waiting time, the magnetic flux value on the surface of the spherical magnet is measured by the sensor of gauss meter and then all of the measured data are transferred into the computer. The obtained data are tabulated by an Excel page. By using the stored data, the flux densities on the surface of the spherical magnet are plotted as shown in Fig. 8 . The Gauss Meter being used for the measurements has a Tesla Unit option.
PPT Slide
Lager Image
Measurement system for flux density of the spherical magnet
PPT Slide
Lager Image
At θ=0° position of the rotor, the comparison of the radial constituents of the magnetic flux densities for the gap between α=-π and +π Model 1-2-3-4.
Simulation results obtained by changing of magnetic field vectors of the magnets for Model-1 among the studied models complies with the results obtained by calculation with a numerical method as shown in Fig. 6 , where surface flux density of the magnet is calculated by a numerical method, and FEM, and measured by an experiment.
From Fig. 8 , it can be seen that the highest flux value on the surface is obtained from Model-1. On the other models, the surface flux values are lower than Model-1 and the fluxes on the surfaces are not homogeneous. Comparisons of calculations and measurements on the spherical magnet surfaces in terms of magnetic flux density are given in Fig. 6 and Fig. 8 . Examining Model-2 and Model-3, the flux values are high, but have too harmonic components. For Model-4, the magnetic flux does not make a sine wave at all. But, the results obtained from Model 1 meet the expectations for a good structure of spherical rotor by permanent magnets.
6. Conclusion
In this study, four-different manipulated spherical magnet are designed and their three-dimensional analyses are realized by using finite-element method via Maxwell environment. The analysis of the magnetic quartile spheres with four different magnetic directions are compared to the each other, and then the quartile sphere with the best surface flux distribution is suggested for rotor structure.
Due to its high and smooth surface flux, which is very important to the architecture of the spherical motor, Model-1 can be used for the spherical motor to get desired torque characteristic. By manipulating different magnet geometries or making different spherical magnet designs, the surface flux values can be improved. In future studies, by using magnets with four different magnetic directions, more available spherical structure can be formed, and more improved flux density on the surface can be achieved, as well.
BIO
Yusuf Oner He received the B.Sc. degree from the Gazi University, Turkey in 1993, the M.Sc. degree from the University of Gazi, Turkey in 1998, all in electrical education. He received the Ph.D. degree in electrical education from the Gazi University, Turkey in 2004. In 1999, he joined the Gazi University, and he was appointed as a Research and Teaching Assistant from 1999 to 2004. He joined the Department of Electrical and Electronic Engineering, Pamukkale University, Denizli, Turkey, in 2007 as a Associate prof. Dr. His research interests include 3D magnetic analysis, magnetic bearing, electrical machine design and analysis, spherical motor.
Selami Kesler He received the BSc, MSc and PhD Degrees from Karadeniz Technical University, Turkey in 1991, in 1998 and in 2006 respectively all in Electrical Engineering. He joined Karadeniz Technical University in 1992 and he worked as a lecturer until 2004. He joined the Department of Electrical Education, Pamukkale University, Denizli, Turkey, in 2004 as a lecturer. He joined the Department of Electrical and Electronic Engineering, Pamukkale University, Denizli, Turkey, in 2007 as a assistant professor. His fields of interest are power electronics, AC drives, dynamic analysis of the electrical machines, intelligent control techniques and advanced DSP’s.
References
Coey J.M.D. 2002 “Permanent Magnet Applications” Journal of Magnetism and Magnetic Materials 248 441 - 456    DOI : 10.1016/S0304-8853(02)00335-9
Choi Jae Seok , Yoo Jeonghoon 2010 “Optimal Design Method for Magnetization Direction of a Permanent Magnet Array” Journal of Magnetism and Magnetic Materials 322 2145 - 2151    DOI : 10.1016/j.jmmm.2010.01.047
Lee. Moon G. , Lee. Sung Q. , Gweon. Dae-Gab 2004 “Analysis of Halbach Magnet Array and its Application to Linear Motor” Mechatronics 14 115 - 128    DOI : 10.1016/S0957-4158(03)00015-1
Ofori-Tenkorang J. , Lang J. H. “A Comparative Analysis of Torque Production in Halbach and Conventional Surface-Mounted Permanent Magnet Synchronous motors” Conference Proceedings IEEE IAS 30th Annual meeting 8-12 October 1995 657 - 663
Stumberger Bojan , Stumberger Gorazd , Hadziselimovic Miralem , Hamler Anton , Gorican Viktor , Jesenik Marko , Trlep. Mladen 2007 “Comparison of Torque Capability of Three-Phase Permanent Magnet Synchrounos Motors with Different Permanent Magnet Arrangement” Journal of Magnetism and Magnetic Materials 316 261 - 264    DOI : 10.1016/j.jmmm.2007.02.113
Lee Moon G. , Gweon. Dae-Gab 2004 “Optimal Design of a Double-Sided Linear Motor with a Multi-Segmented Trapezoidal Magnet Array for a High Precision Positioning System” Journal of Magnetism and Magnetic Materials 281 336 - 346    DOI : 10.1016/j.jmmm.2004.03.052
Babic S. I. , Akyel C. 2008 “Improvement in the analytical calculation of the magnetic field produced by permanent magnet rings,” Progress In Electromagnetics Research C 5 71 - 82
Ravaud R. , Lemarquand G. , Lemarquand V. , Depollier C. 2008 “The three exact components of the magnetic field created by a radially magnetized tile permanent magnet,” Progress In Electromagnetics Research, PIER 88 307 - 319    DOI : 10.2528/PIER08112708
Ravaud R. , Lemarquand G. , Lemarquand V. , Depollier C. 2009 “Discussion about the analytical calculation of the magnetic field created by permanent magnets,” Progress In Electromagnetics Research B 11 281 - 297    DOI : 10.2528/PIERB08112102
Davey. K. , Vachtsevanos. G.J. , ve Power. R. 1987 “An Analysis of Fields and Torques in Spherical Induction Motors” IEEE Transactions on Magnetics 23 273 - 281    DOI : 10.1109/TMAG.1987.1064749
Foggia. A. , Oliver. E. , Chappuis. F. 1988 “New Three Degree of Freedom Electromagnetic Actuator” Conference Record IAS Anual Meeting New York
Hollis. R.L. , Ralph L. 1989 “Magnetically Levitated Fine Motion Robot Wrist with Programmable Compliance” EP Pat. EP0294760 A1
Kaneko. K. , Yamada. I , Itao. K. 1989 “A Spherical DC Servo Motor with Three Degrees of Freedom” Transactions of the ASME Journal of Dynamic System Measurement and Control 111 398 - 402    DOI : 10.1115/1.3153067
Lim C.K. , Yan L. , Chen I.M. 2004 Mechanical design & numerical electromagnetic analysis of a DC spherical actuator Conference on Robotics, Automation and Mechatronics, Proceedings of the IEEE 536 - 541
Lee. K.M. , Kwan C. 1991 “Design Concept Development of A Spherical Stepper for Robotic Applications” IEEE Transactiona on Robotics and Automation 7 175 - 181    DOI : 10.1109/70.68082
Roth. R. , Lee. K.M. 1995 “Design Optimisation of A Three Degrees of Freeedom Variable Reluctance Spherical Wrist Motor” Transactions of the ASME Journal of Engineering for Industry 117 378 - 388    DOI : 10.1115/1.2804344
Zhou. Z. , Lee. K.M. 1996 “Real Time Motion Control of Multi Degree of Freedom Variable Reluctance Spherical Motor” Proceedings of the 1996 IEEE International Conference on Robotics and Automation Minneopolis-Minnesota 2859 - 2864
Zhao C. , Li Z. , Huang W. 2005 Optimal design of the stator of a three-DOF ultrasonic motor Sensors Actuators A 121 494 - 499    DOI : 10.1016/j.sna.2005.03.010
Shigeki T. , Guoqiang Z. , Osamu M. 1996 Development of new generation spherical ultrasonic motor International Conference on Robotics and Automation, Proceedings of the IEEE 2871 - 2876
Wang. J. , Jewell. G. , Howe. D. 1998 “Analysis, Design and Control of Novel Spherical Permanent Magnet Actuator” IEE Proceedings on Electrical Power Applications 145 61 - 71    DOI : 10.1049/ip-epa:19981635
Wang. J. , Jewell. G. , Howe. D. 1997 “Modeling of a Novel Spherical Permanent Magnet Actuator” Proceeding of IEEE International Conference on Robotics and Automation New Mexico 1120 - 1125
Lee S. J. , Jiles D.C. 2000 “Modeling and simulation of a permanent magnet array in elliptical configurations” MSM 2000 Technical Proceedings of the 2000 International Conference on Modeling and Simulation of Microsystems 668 - 671
Oner Y. 2007 “A Permanent Magnet Spherical Rotor Design and Three Dimensional Static Magnetic Analysis” Sensors and Actuators A 137 200 - 208    DOI : 10.1016/j.sna.2007.02.018
Xu Yanliang , Feng Kaijie “Analysis on Toothless Permanent Magnet Machine with Halbach Array” IPEMC 2006
Li Hongfeng , Xia Changliang , Shi Tingna 2007 “Spherical Harmonic Analysis of a Novel Halbach Array PM Spherical Motor” Proceeding of IEEE International Conference on Robotics and Biomimetics Sanya, China 2085 - 2089
Xia Changliang , Xin Jianguo , Li Hongfeng , Shi Tingna 2013 “Design and Analysis of a Variable Arc Permanent Magnet Array for Spherical Motor” IEEE Transactions on Magnetics 49 (4) 1470 - 1478    DOI : 10.1109/TMAG.2012.2231092
Park Hyun-Jong , Lee Ho-Joon , Cho Su-Yeon , Ahn Han-Woong , Lee Ki-Doek , Park Chan-Yeop , Won Sung-Hong , Lee Ju 2013 “A Performance Study on a Permanent Magnet Spherical Motor” IEEE Transactions on Magnetics 49 (5) 2307 - 2310    DOI : 10.1109/TMAG.2013.2238616
Park Hyun-Jong , Cho Su-Yeon , Ahn Han-Woong , Lee Ho-Joon , Won Sung-Hong , Lee Ju 2012 “A Study of Advanced Spherical Motor for Improvement of Multi-DOF Motion” Journal of Electrical Engineering & Technology 7 (6) 926 - 931    DOI : 10.5370/JEET.2012.7.6.926
Lee Ho-Joon , Park Hyun-Jong , Won Sung-Hong , Ryu Gwang-Hyun , Lee Ju 2013 “Improvements of Performance of Multi-DOF Spherical Motor by Doubled Air-gap Feature” Journal of Electrical Engineering & Technology 8 (1) 90 - 96    DOI : 10.5370/JEET.2013.8.1.090
Li Bin , Yu Rujian , Li Hua , Li Guidan 2014 “Design Considerations of a Permanent Magnetic Spherical Motor Using Spherical Harmonics” IEEE Transactions on Magnetics 50 (8)
Qian Zhe , Wang Qunjing , Li Guoli , Guo Xiwen , Hu Cungang , Yan Hui 2013 “Design and Analysis of Permanent Magnetic Spherical Motor with Cylindrical Poles” International Conference on Electrical Machines and Systems Busan, Korea 644 - 649
Lee Hyung-Woo , Kang Dong-Woo , Lee Ju 2010 “Torque Calculation Method of a Permanent Magnet Spherical Motor” Journal of Electrical Engineering & Technology 5 (3) 431 - 434    DOI : 10.5370/JEET.2010.5.3.431
Guo Jinjun , Kim Dong-Hoon , Son Hungsun 2013 “Effects of Magnetic Pole Design on Orientation Torque for a Spherical Motor” IEEE/ASME Transactions on Mechatronics 18 (4) 1420 - 1425    DOI : 10.1109/TMECH.2013.2242900
Oner Y. 2004 “Sabit Mıknatıslı Küresel Motorun Bilgisayar Destekli Üç Boyutlu Manyetik Analizi, Tasarımı ve Uygulaması”, Ph.D Thesis Gazi University, Institıte of Technology
2000 Ansoft: Maxwell 3D, Getting started: a 3D Magnetic Force Problem
Xia Changliang , Li Hongfeng , Shi Tingna 2008 “3D Magnetic Field and Torque Analysis of a Novel Halbach Array Permanent-Magnet Spherical Motor” IEEE Transactions on Magnetics 44 (8)