Download Article (PDF)
Total Page:16
File Type:pdf, Size:1020Kb
Proceedings of 2012 International Conference on Mechanical Engineering and Material Science (MEMS 2012) Analysis of 3D Magnetic Field of Magnetic Suspension Spherical Reluctance Motor Hu Daotian Zeng Li College of Mechanical Engineering Yangzhou University College of Mechanical Engineering Yangzhou University Yangzhou University,YZU Yangzhou University,YZU Yangzhou, China,18252718210 Yangzhou, China,13665203656 [email protected] [email protected] Abstract—It is a very complicated problem that analyzing and equator line of spherical rotor. stator 1 and stator 3 keep coaxial calculating the magnetic field of the magnetic suspension with the x-axis; the function is to drive the rotor to rotate around spherical reluctance motor, because the structure of the stator x-axis and produce magnetic suspension force that makes the and rotor, the distribution of magnetic field present rotor stabilize in the x-axis direction, Stator 2 and stator 4 keep three-dimension. By calculating 3-D magnetic field and analyzing coaxial with the y-axis, which is aimed at driving the rotor to the magnetic field of the magnetic suspension spherical rotate around the y-axis and producing magnetic suspension reluctance motor, Then we draw the magnetic flux density force that makes the rotor stabilize in the y-axis direction; distribution curve of the salient pole stator and rotor in different Stator 5 is located in the top of spherical rotor and keeps coaxial position, when the motor stator is in different current. with z-axis, its function is to drive the spherical rotor to rotate around z-axis and produce stable magnetic suspension force in Keywords-magnetic suspension; finite element method; 3D the z-axis direction. magnetic field; spherical reluctance motor The surface of spherical rotor of magnetic suspension I INTRODUCTION spherical motor has mutual orthogonal grooves. The stator includes six salient poles and one ladder circle. Based on Multiple degrees of freedom motor refers to a kind of reluctance motor principle, there are multiphase windings motor which has two or three rotating degrees, the motor can which drive the stator to rotate in salient poles of stator. After rotate round fixed-point axis of the space. The spherical motor each phase winding is charged with electricity , it not only belongs to multiple degrees of freedom (MDOF) motor. produces electromagnetic torque which drives the rotor to rotate, MDOF motor has a wide application prospect in robots, but also provides radial magnetic suspension force for rotor.[4] multiple coordinates mechanical processing center, aerospacecraft, electric gyroscope, all-round tracking antenna, turret turntable , prosthetic human body , medical equipment , camera work station, panorama photography work station, blender, globe valve, and other multiple degrees of freedom equipments.[1] It is too difficult to simplify magnetic field of spherical motor by using 2D magnetic field because it has complex boundary conditions, anisotropy, nonlinearity and other characteristics , and that it is a typical 3D magnetic field, we must further discuss its calculation problem. Currently, the subject about the design and analysis of spherical motor is still being studied at home and abroad.[2-3] Fig.1 The structure of the magnetic suspension spherical reluctance motor This study plans to use 3D finite element analysis(FEA)method to calculate and analyze the magnetic III. ANALYSIS OF MAGNETIC FIELD OF MAGNETIC field of MDOF magnetic suspension spherical reluctance SUSPENSION SPHERICAL RELUCTANCE MOTOR motor in order to offer some basis for the design and A. The Basic Principle of 3D Magnetic Field Calculation of calculation of MDOF magnetic suspension spherical reluctance Magnetic Suspension Spherical Reluctance Motor motor. In the finite element of 3D magnetic field, three II. THE BASIC STRUCTURE OF MAGNETIC SUSPENSION components of the magnetic vector potential A are Ax, Ay, Az. [5-6] SPHERICAL RELUCTANCE MOTOR In the magnetic field of 3D linear media, the magnetic A 0 The stator and rotor 3D structure of magnetic suspension vector potential A satisfy following demands under : spherical reluctance motor are shown in Fig.1,the rotor is a unit 2 JA ,In the magnetic field V; salient pole sphere that has grooves, and it locates in the middle of the motor; Four stators are symmetrically distributed on the The China Natural Science Foundation for this project. Item Number: 50975249.and Natural Science Foundation of Universities in Jiangsu Pro for this project. Item Number: 08KJB460008. © 2012. The authors - Published by Atlantis Press 5 AC/ k k k k s ,on the boundary S; P P P P ()()]}i j A j i A J Nk dxdydz Corresponding conditions variational problem is: x z xj y z yj y i 1 2 Note that its value is not zero only when i is the node of FA()[ | A | J A ] dxdydz min,AC | s this unit; Therefore, unit coefficient matrix corresponding to 2 one unit is SSS , among the equation is third If the whole magnetic field is discrete into k0 units and p0 k ij k ij k nodes. In the rectangular coordinate system,the magnetic matrix. S is composed of S square matrix that the k ij k vector potential Ai of any node i has three components, which are called Axi ,Ayi, Azi; Therefore,each component of the number is p0 p 0 .The same magnetic vector potential of p0 nodes can be lined into the argument, k k k k k k T , among the NNNNNNNk ... following column vector: x1y 1z 1 xp0yp 0zp 0 equation, [][AAAAAA x1 y 1 z 1 x 2... xp ] k k , k k , Nxi Jx Pj dxdydz Nyi Jy Pj dxdydz If each unit has p0 nodes,the shape function corresponding vk vk to node i is Pk ,and the magnetic vector potential A of any i Nk J Pk dxdydz, i 1,2,..., p . p zi z j 0 0 k v point of each unit can be expressed as APA i i . Three k i1 We should solve the coefficient matrix of above unit, then components of the magnetic vector potential A are respectively realize further expansion , getting the finite element equation expressed as that is SAN . Solution of above equation is [A]. p p p APAAPAAPA 0 k , 00k , k x i1 i xi y i11i yi z i i zi B. Establishment of Model and the Decision of Boundary Energy functional F(A) of the whole magnetic field can be Condition expressed as the sum of each unit energy functional: Numerical analysis for electromagnetic field of motor k FAFAFA 0 mainly use the finite element method that is the most effective i1 k and the most widely used at present, boundary element method, The unit energy functional is and finite difference method. ANSYS is a kind of finite 1 1 element analysis software ,which is widely used at FA()[ |A |2 J A ] dxdydz { present.[7]The basic principle of analysis for electromagnetic k 2 2 kkfield using ANSYS is that processing object first should be divided into finite unit(including some nodes),then magnetic Az Ay 2 AxA z 2 Ay Ax 2 [( )()()] potential or electric potential of the each node under certain y z z x x y boundary conditions and initial conditions can be calculated (J A J A J A)} dxdydz according to magnetic vector potential and electric scalar x x y y z z potential. According to extreme value theory,the variational problem is equivalent to the following equation: TAB.1 THE STRUCTURAL PARAMETERS OF MAGNETIC SUSPENSION SPHERICAL RELUCTANCE MOTOR F k0 F F k F the radius of spherical top of stator salient pole/mm 51 k 0 , 0 k 0 , k 1 k 1 Axi Axi Ayi Ayi unilateral air gap of between stator surface and rotor surface/mm 1 F k F 0 k 0 k 1 Azi Azi the radius of outside spherical surface of rotor/mm 50 In the above equation, any unit k after discretization has p0 the radius of inner spherical surface of stator/mm 20 nodes. Then the ith node of this unit can be showed as k k k k the radius of outside spherical surface of stator/mm 38 F p 1 P P P P k 0 { [( i j i j )A j1 xj Axi k y y z z the length of air gap/mm 1 Pk Pk Pk Pk the number of stator groove 6 ()()]}j i A j i A J Nk dxdydz x y yi x z zj x i the number of groove in rotor latitude direction 7 k k k k F p 1 P P P P k 0 { [( i j i j )A the number of groove in rotor longitude direction 12 j1 yj Ayi k x x z z k k k k the number of stator coil 50 P Pj Pj P ()()]}i A i A J Nk dxdydz polar altitude of stator/mm 25 x y xj y z zj y i coil current/A I k k k k F p 1 P P P P k 0 { [( i j i j )A A j1 x x z z zj Based on magnetic potential and electric potential, we can zi k get other physical quantities of electric field such as the 6 distribution of line of magnetic force, magnetic induction, the motor becomes larger, and the magnetic flux density of electric displacement flux, energy of electromagnetic field, section where it is in the air gap becomes smaller and smaller. magnetic force, torque, inductance, capacitance and so on But when it enters into the next magnetic pole, the magnetic through dealing with the solution. Due to the symmetry of flux density of section starts to become larger again.