LONG STROKE MAGNETIC LEVITATION PLANAR STAGES by Irfan-ur-rab Usman B.A., The University of British Columbia, 2008 B.A.Sc., The University of British Columbia, 2008 M.A.Sc., The University of British Columbia, 2010 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY in THE FACULTY OF GRADUATE AND POSTDOCTORAL STUDIES (Mechanical Engineering) THE UNIVERSITY OF BRITISH COLUMBIA (Vancouver) December 2015 Ā© Irfan-ur-rab Usman, 2015 ii Abstract Modern positioning applications often require long strokes in multiple degrees-of-freedom (DOF). One solution to such requirements is a planar stage capable of simultaneous large strokes in X- and Y-translation. An ideal stage concept is based on a planar motor which uses non-contact multi-axis forces to directly actuate a single moving body, without any connecting elements or bearings structures that would induce additional structural modes as well as excess inertia. This thesis presents the design, analysis and experimental results of two planar motors: 1) a permanent magnet synchronous planar motor, and 2) a permanent magnet asynchronous planar motor. The permanent magnet synchronous planar motor is comprised of multiple one-dimensional (1-D) magnet arrays attached to a moving stage and multiple stationary 1-D coils built as a printed circuit board (https://www.youtube.com/watch?v=-r4Tv7GbB8o). This motor addresses many issues with existing synchronous planar levitation motors, including scalability of the XY stroke with minimal increase in controller or drive complexity, and simplified commutation and control due to natural force/torque decoupling and no coil/magnet array edge effects. Modeling, analysis and design of a prototype are shown with motion control results. To use this levitation motor for high accuracy positioning applications, force and torque characteristics of the 6-DOF stage must be highly linear in all axes in order to minimize controller effort and reduce intrinsic motor disturbances. Novel magnet array designs are presented which self-attenuate both force and torque ripple in all 6-DOF, without additional controller or drive complexity. These array designs are tested via simulation and experiment. The permanent magnet asynchronous planar motor is based on the asynchronous induction effect to simultaneously levitate and propel the motion stage. The advantages of this type of planar iii levitation motor are: 1) passive stability in all axes; and 2) simple slab-type homogenous stator comprised of a basic conducting material such as aluminum or copper. This thesis presents novel modeling and analysis and provides a new analytical expression for the generated levitation force and torque with geometric and material property inputs. Experimental force and torque measurements are carried out and compared to the analytical model. iv Preface The research presented in this thesis has been carried out at the University of British Columbia Precision Mechatronics Laboratory (UBC PML), Department of Mechanical Engineering, under the supervision of Dr. Xiaodong Lu. Chapter 1 is a brief literature review of existing multi-degree of freedom direct drive stages using magnetic field interactions. Chapter 2 is based in part on the work published in ļ· Xiaodong Lu, Irfan Usman, "6D Direct-Drive Technology for Planar Motion Stages", Annals of the CIRP-Manufacturing Technology, Vol. 61, No. 1, pp. 359-362, 2012. In this work we present the novel 6-DOF synchronous levitation motor including field, force and torque modeling as well as commutation laws. Experimental motion tracking results from a prototype are also presented. The author modeled, simulated, designed, built and assembled the 6-DOF levitation motor. The 6-DOF position sensor was based on the work of Niankun Rao [1]. The real-time control computer was based on the work of Kristofer Smeds [2]. The 6-DOF motion controller and the multi-channel PWM current driver were implemented by Dr. Xiaodong Lu. Chapter 2 is also based in part on the work published in ļ· Irfan-ur-rab Usman, Xiaodong Lu, "Force Ripple Attenuation of 6-DOF Direct Drive Permanent Magnet Planar Levitating Synchronous Motors", IEEE Transactions on Magnetics. Vol. 51, No. 12, pp. 1-8, Dec. 2015. v In this work a novel magnet array design is presented for force ripple attenuation. The author modelled, simulated, designed, built and measured the load characteristic of the novel split magnet array. Portions of the work have also been presented in conference proceedings as follows: ļ· I. Usman; X. Lu; "Force ripple attenuation in 6DOF planar levitated motor"; In: 2nd International Conference on Virtual Machining Process Technology, 2013 ļ· I. Usman; X. Lu; "Torque and force ripple attenuation of 6DOF levitating motor"; In: 3rd International Conference on Virtual Machining Process Technology, 2014 Chapter 3 includes work on asynchronous levitation motors that is currently part of a manuscript that has been submitted to IEEE Transactions on Magnetics. In this work, a new analytical model for levitation force and drag torque on a rotating magnet disk over a finite thickness homogeneous conducting slab is presented, and experimental load characterization is carried out. This novel analytical model is useful for informing physical understanding and design intuition, especially when compared to existing numerical methods. The rotating magnet disk can be used to make a 6-DOF asynchronous type planar levitation motor [3]. The author carried out the modeling, analysis, simulation and experimental load characterization presented. vi Table of Contents Abstract .......................................................................................................................................... ii Preface ........................................................................................................................................... iv Table of Contents ......................................................................................................................... vi List of Tables ..................................................................................................................................x List of Figures ............................................................................................................................... xi Notation ....................................................................................................................................... xix Acknowledgements ......................................................................................................................xx Dedication ................................................................................................................................... xxi Chapter 1: Introduction ................................................................................................................1 1.1 Thesis Overview ............................................................................................................. 4 1.2 Thesis Contributions ....................................................................................................... 5 1.3 Machine Architectures for Planar Positioning ................................................................ 6 1.3.1 Serial Machine Architecture ....................................................................................... 6 1.3.2 Parallel Kinematic Architecture .................................................................................. 7 1.3.3 Direct Drive Architecture ........................................................................................... 8 1.4 Multi-DOF Direct Drive Machines ................................................................................. 9 1.4.1 Long Stroke Synchronous Magnetic Planar Levitation Motors ............................... 12 1.4.2 Long Stroke Asynchronous Magnetic Planar Levitation .......................................... 22 Chapter 2: Synchronous Planar Levitation Stage ....................................................................29 2.1 Synchronous Levitation Concept and Working Principle ............................................ 30 2.2 Magnetic Field, Force and Torque Generation Modeling for 2-DOF Motor Element . 34 vii 2.2.1 3D Magnetic Field Distribution for Cuboidal Magnet ............................................. 35 2.2.2 Force and Torque Generation on Magnet Array ....................................................... 43 2.2.3 Principle of Equivalent Force Generation ................................................................. 45 2.2.4 Magnetic Field Modeling for Infinitely Extended Halbach Array ........................... 54 2.2.5 Force Due to Single-Phase Repeating Coil Group ................................................... 63 2.2.6 Force Due to 3-Phase Repeating Coil Group with Commutated Current ................. 68 2.3 Force and Torque Ripple Self-attenuation Using Magnet Array Design ..................... 72 2.3.1 Force Vector Rotation from Offset Commutation Law ............................................ 73 2.3.2 Split Magnet Array Design for Self-attenuation of Force Ripple ............................. 77 2.3.3 Experimental Load Verification for Force Ripple Self-attenuation ......................... 82 2.3.4 Quad-Split Array Design for Self-attenuation of Force and Torque Ripple ............. 87 2.4 6-DOF Motion Stage Prototype Design ........................................................................ 99 2.4.1 Prototype Architecture .............................................................................................. 99 2.4.2 Mover Design.......................................................................................................... 103 2.4.2.1 Magnet Array Geometric and Material Parameters ........................................ 103 2.4.2.2 Mover Magnet Array Configuration ............................................................... 105 2.4.2.3 Mover Mechanical Design .............................................................................. 109 2.4.3 Stator Coil Array PCB Board Design and Wiring Configuration........................... 111 2.4.4 Stator Design ........................................................................................................... 123 2.5 6-DOF Prototype Manufacturing and Assembly ........................................................ 127 2.5.1 Mover ...................................................................................................................... 127 2.5.2 Stator ....................................................................................................................... 133 2.5.3 Metrology Reference Frame ................................................................................... 135 viii 2.5.4 Complete System .................................................................................................... 136 2.6 Motion Tracking Results with 6-DOF Prototype........................................................ 138 2.6.1 Plant Modeling and Control Architecture ............................................................... 138 2.6.2 6-DOF Tracking Results ......................................................................................... 139 Chapter 3: Asynchronous Planar Levitation Stage ................................................................143 3.1 Asynchronous Planar Levitation Stage Concept and Working Principle ................... 144 3.2 Asynchronous Levitation Machine Force and Torque Modeling ............................... 148 3.2.1 Infinitely Wide 2D Moving Magnet Array over Homogeneous Conducting Slab . 148 3.2.2 Analytical Shells Model for Force and Torque of Rotating 3D Magnet Disk over Homogeneous Conducting Slab .......................................................................................... 161 3.3 Experimental Force, Torque and Power for Rotating Magnet Disk ........................... 169 Chapter 4: Conclusion ...............................................................................................................180 4.1 Synchronous Planar Levitation Motor ........................................................................ 182 4.1.1 Contribution I: Motor topology for New Type of 6-DOF Synchronous Planar Levitation Motor ................................................................................................................. 182 4.1.2 Contribution II: Analytical Model of Field, Force, Torque and Commutation Laws for New Type of 6-DOF Synchronous Planar Levitation Motor ........................................ 183 4.1.3 Contribution III: Prototype and Experimental Demonstration of New Type of 6-DOF Synchronous Planar Levitation Motor ................................................................................ 183 4.1.4 Contribution IV: Novel Split and Quad-Split Magnet Array Designs for Force and Torque Ripple Self-attenuation ........................................................................................... 184 4.1.5 Contribution V: Experimental Verification of Novel Split Magnet Array Design for Force Ripple Attenuation .................................................................................................... 185 ix 4.2 Asynchronous Planar Levitation Motor ...................................................................... 185 4.2.1 Contribution I: Novel Analytical Model of Field, Force and Torque of Magnet Disk Rotating Above a Finite Thickness Homogenous Conductor............................................. 185 4.2.2 Contribution II: Evaluation of Novel Analytical Model with Finite Element Simulation and Characterization with Experimental Load Measurements ......................... 186 4.3 Future Work ................................................................................................................ 186 4.3.1 Future Work for 6-DOF Synchronous Planar Levitation Motor ............................ 187 4.3.2 Future Work for 6-DOF Asynchronous Planar Levitation Motor Modeling and Analysis............................................................................................................................... 188 Bibliography ...............................................................................................................................189 x List of Tables Table 2.1 Spacing for Split Array Targeting Different Force Harmonics .................................... 79 Table 2.2 Force Harmonic Coefficients in Z ................................................................................ 87 Table 2.3 Analytical Mean and 6 cpĪ» Ripple Forces for Sub-arrays. ........................................... 93 Table 2.4 Analytical Mean and 6 cpĪ» Ripple Torques for Sub-arrays. ........................................ 94 Table 2.5 Force and Torque Amplitudes. .................................................................................... 97 Table 2.6 Alpha Prototype Motor Parameters. .......................................................................... 122 Table 2.7 Predicted Alpha Prototype Motor Performance at flying height 0.5 mm. ................. 137 Table 3.1 2D Asynchronous Levitation Machine Parameters. .................................................. 161 Table 3.2 3D Rotating Disk Asynchronous Levitation Machine Parameters. ........................... 169 Table 3.3 3D Experimental rotating disk asynchronous levitation machine parameters. .......... 173 xi List of Figures Figure 1.1 Novel 6-DOF synchronous planar levitation motor. ..................................................... 3 Figure 1.2 Asynchronous planar levitation motor mover (onboard battery pack, rotary motor controllers, and stator conducting slab not shown). ........................................................................ 3 Figure 1.3 Serial machine architectures. (a) XY H-drive style planar stage. (b) Dual stage. ......... 7 Figure 1.4 Parallel machine architecture for 6-DOF motion. ......................................................... 8 Figure 1.5 Direct drive machine architecture for 6-DOF positioning stage. .................................. 9 Figure 1.6 Planar motor for very small scale movers (such as 2mm by 2mm by 0.4mm) based on diamagnetic levitation and Lorentz force translation, figure adapted from [29]. (a) Top view of mover magnets. (b) Electromagnetic configuration of diamagnetic planar levitation motor. ...... 11 Figure 1.7 Kimās 6-DOF synchronous levitation motor, figure adapted from [36]. .................... 13 Figure 1.8 Scaling up planar stroke of Kim and Trumperās synchronous levitation motor. ........ 14 Figure 1.9 Moving coil synchronous planar levitation motor (umbilical cable to mover not shown), figure adapted from [39]. Red magnet blocks are magnetized in +z, and blue magnet blocks are magnetized in āz. ......................................................................................................... 15 Figure 1.10 Shifting center of pressure of example 3-phase moving coil levitation motor, given a constant levitation force command (umbilical cable not shown). ................................................ 16 Figure 1.11 Moving magnet 2D chessboard type synchronous planar levitation motor with multiple alternating layers of racetrack style coils; figure adapted from [43]. Red magnet blocks are magnetized in +z, and blue magnet blocks are magnetized in āz. .......................................... 18 Figure 1.12 Moving magnet 2D chessboard type synchronous planar levitation motors. (a) Alternating zones of racetrack coils, figure adapted from [45]. (b) Herringbone coil pattern, figure adapted from [46]. .............................................................................................................. 18 xii Figure 1.13 Scaling up XY stroke of moving magnet 2D chessboard type synchronous planar levitation motor (example using motor type similar to [43]). ....................................................... 19 Figure 1.14 Moving magnet 2D chessboard type synchronous planar levitation motor with alternating layers of straight coils; figure adapted from [54]. ...................................................... 20 Figure 1.15 6-DOF moving magnet stage with axially magnetized cylindrical magnets over array of circular coils; figure adapted from [55]. ................................................................................... 21 Figure 1.16 Circular type 2-DOF induction planar motor, figure adapted from [57]. .................. 23 Figure 1.17 Asynchronous planar motor using electromagnets and back-iron, figures adapted from [4] (out-of-plane bearings not shown). (a) Linear 1-DOF induction motor. (b) 3-DOF configuration of asynchronous planar motor. ............................................................................... 23 Figure 1.18 Asynchronous levitation motor using permanent magnet Halbach array, adapted from [58]. ...................................................................................................................................... 24 Figure 1.19 Permanent magnet asynchronous levitation motor elements showing example magnetization patterns. (a) Circumferential magnet disk. (b) Radial magnet wheel. ................... 25 Figure 2.1 Synchronous planar levitation 6-DOF motor electromagnetic topology. ................... 31 Figure 2.2 Y2 magnetic array and coils cross-section .................................................................. 32 Figure 2.3 Cuboidal permanent magnet with vertical (z-directed) uniform magnetization M shown in red. ................................................................................................................................. 37 Figure 2.4 Example 3D finite extension Halbach-patterned magnet array (magnetization directions shown in red). (a) 3D isometric view of finite Halbach array. (b) 2D view of finite Halbach array showing segment dimensions. ............................................................................... 40 Figure 2.5 Analytical field model, flux density vector components of BFEA for example magnet array at zm = -Ī»/5. ....................................................................................................................... 41 xiii Figure 2.6 Error in flux density vector components for zm = -Ī»/5. Error defined as BFEA-BANALYTICAL. ........................................................................................................................ 42 Figure 2.7 General 3D magnet array interacting with differential conductor element. ................ 44 Figure 2.8 Finite x-directed extension Halbach-patterned magnet array. ..................................... 46 Figure 2.9 Infinite x-directed extension Halbach-patterned magnet array ................................... 46 Figure 2.10 Bz from 3D analytical model along test line (zm = -Ī»5, ym = 0). ............................. 47 Figure 2.11 Equivalent force cases. (a) Ī»-width magnet array interacting with Ī»-repeating coil array. (b) repeating Ī»-width magnet arrays with single coil. (c) ā-width magnet array with single coil................................................................................................................................................. 52 Figure 2.12 (a) Wm-width magnet array interacting with Ī·c-repeating coil array. (b) Superposition showing overlap of each identical magnet array. .................................................. 53 Figure 2.13 2D form of the infinite extension magnet array. ....................................................... 56 Figure 2.14 Magnetization distribution within magnet volume for ā-width Halbach array with 4 magnet segments/Ī» separated by gap g. ........................................................................................ 58 Figure 2.15 Single phase coil array below finite width Halbach array. ........................................ 63 Figure 2.16 Partial harmonic model versus 3D analytical model at different zm (a) Bz, partial harmonic model (PHD), 3D analytical model (3DAM) evaluated at ym = 0 with Dm = 20Ī», and error (PHD-3DAM), for ā-width magnet arrays. (b) Fourier coefficients for error at different zm........................................................................................................................................................ 66 Figure 2.17 Modified 3-phase coil array below finite width Halbach array. ................................ 70 Figure 2.18 Multiple stationary coil array stack-up dimensions, N-layers of identically driven coil arrays. ..................................................................................................................................... 72 xiv Figure 2.19 Offset magnet array from commutation center. ........................................................ 74 Figure 2.20 Force vector rotation due to commutation center offset. (a) Position independent force vector F0. (b) Force ripple vector F6. (c) Force ripple vector Fk Ā± 1 resulting from field harmonic Bk. ................................................................................................................................. 76 Figure 2.21 Split array solution for force ripple attenuation of the 6th harmonic of force ripple due to the 5th field harmonic. ........................................................................................................ 77 Figure 2.22 FEA models for āinfiniteā coil arrays (extended in the x-direction over 18Ī» width). (a) Non-split array. (b) Split array. ............................................................................................... 80 Figure 2.23 Finite element analysis of non-split array and split array over āinfiniteā coil arrays( extended over 18Ī» width). Flying height zf = 0.5 mm, Ixr = 0, Izr = 10 [A], Ī» = 30 mm, Br =1.325 T, Dm = 2Ī»,Wc = 0.1583Ī», tc = 0.213 [mm] and 8 layers of active coils spaced 0.643 [mm] apart in the -zm direction. ........................................................................................ 81 Figure 2.24 Experimental arrays for load measurement, showing dimensions and magnetization directions. ...................................................................................................................................... 83 Figure 2.25 Experimental setup for load measurement. ............................................................... 83 Figure 2.26 Experimental test cases with finite coil array of 4Ī» width. (a) Non-split array. (b) Split array. ..................................................................................................................................... 84 Figure 2.27 Comparison experimental non-split array and split array, uncompensated for zf =0.0247Ī» , Ixr = 0, Izr = 9.5 A, Ī» = 30 mm, Dm = 2Ī»,Wc = 0.1583Ī», Hm = Ī»/4, tc = 0.213 mm and N = 8 layers of active coils spaced tg = 0.643 mm apart in the -zm direction. ................... 86 Figure 2.28 Magnet array designs showing force ripple vectors acting at magnetic centers of pressure. (a) Non-split magnet array. (b) Split magnet array. ...................................................... 89 xv Figure 2.29 Quad-split array for force and torque attenuation. (a) Embodiment 1. (b) Embodiment 2. (c) Sub-array example Halbach magnetization pattern. ...................................... 90 Figure 2.30 Quad āsplit force ripple vectors acting at equivalent moment arms of sub-arrays (LI, LII, RI, RII) . ............................................................................................................................. 91 Figure 2.31 Finite element model schematic showing generalized magnet array in relation to stationary coil array with 8 coil layers spaced apart by 0.643 mm in zm. .................................... 95 Figure 2.32 Simulated magnet arrays showing magnetization patterns. (a) Non-split array. (b) Split array with Ļ = Ī»/10. (c) Quad-split array with Ļ = Ī»/10. .................................................. 95 Figure 2.33 FEA forces and torques for non-split array, split array, and quad-split array; zf =Ī»/60 = 0.5mm, Ixr = 0, Izr = 10 A, with 8 coil layers spaced apart by 0.643 mm in zm. ........ 96 Figure 2.34 Spatial harmonic amplitudes for force and torque ripples from FEA; flying height zf = Ī»/60 = 0.5 mm, Ixr = 0, Izr = 10 A, with 8 coil layers spaced apart by 0.643 mm in zm. (a) Translation force, Fx. (b) Levitation force, Fz. (c) Torque in ym-axis, Ty ........................ 98 Figure 2.35 Prototype architecture schematic showing 6-DOF camera metrology, power amplifier and real-time controller. .............................................................................................. 100 Figure 2.36 Schematic cross-section of alpha prototype. ........................................................... 101 Figure 2.37 Alpha prototype solid model showing 6-DOF camera metrology. ......................... 102 Figure 2.38 Alternate mover magnet array layouts (identical Ī» = 30mm). ............................... 106 Figure 2.39 Alpha prototype mover design, exploded view. ...................................................... 110 Figure 2.40 Isometric cut-view of mover . ................................................................................. 110 Figure 2.41 Alpha prototype mover design, detailed planform. ................................................. 111 Figure 2.42 Coil structure isometric view, single Ī»-group Y-elongated coil traces. .................. 113 Figure 2.43 Eight layer coil structure, showing alternating X- and Y-elongated coil traces. ..... 114 xvi Figure 2.44 Single phase multilayer windings, shown for three layers of Y-elongated traces. .. 115 Figure 2.45 Three phase wiring, shown for a single Ī»-coil group. ............................................. 116 Figure 2.46 Serial connections between two single Ī»-coil groups with 3-phase wiring. ............ 116 Figure 2.47 Extendable coil array PCBs. Inset shows connector PCB with solid pin connectors...................................................................................................................................................... 118 Figure 2.48 Coil structure isometric view, X- and Y-elongated coils. (a) Layer X1: x-elongated traces. (b) Layer Y1: y-elongated traces. .................................................................................... 119 Figure 2.49 Coil array PCB showing coil groupings, X- and Y-elongated coils. ....................... 120 Figure 2.50 Coil array PCB stack-up details. ............................................................................. 121 Figure 2.51 Alpha prototype stator exploded view. .................................................................... 124 Figure 2.52 Alpha prototype stator isometric cut-view. ............................................................. 125 Figure 2.53 Isometric cut-view of metrology reference frame. .................................................. 126 Figure 2.54 Magnet sub-array assembly jig. ............................................................................... 128 Figure 2.55 Magnet sub-array assembly jig with magnet array installed. .................................. 129 Figure 2.56 Magnet array assembly jig....................................................................................... 130 Figure 2.57 Magnet assembly. .................................................................................................... 132 Figure 2.58 Stator assembly steps. .............................................................................................. 134 Figure 2.59 Metrology reference frame assembly. ..................................................................... 135 Figure 2.60 Synchronous planar levitation motor system. ......................................................... 136 Figure 2.61 Controller block diagram. ........................................................................................ 139 Figure 2.62 Elliptical path of mover in XY, at zf = 1 mm. ....................................................... 141 Figure 2.63 Tracking error for elliptical path, zf = 1 mm. ........................................................ 142 Figure 2.64 Mover floating with zf = 2.5 mm. .......................................................................... 142 xvii Figure 3.1 Asynchronous levitation planar motor concept using four levitation disks, magnetization shown in red for magnet disk 1 (central actuated pivot mechanism and onboard battery pack/motor controllers not shown). ................................................................................ 146 Figure 3.2 Levitation and propulsion principle (details of actuated pivot mechanism not shown)...................................................................................................................................................... 147 Figure 3.3 Stationary Halbach array of infinite extent over a translating conducting slab. ....... 150 Figure 3.4 Magnetic pressure on magnet array, 2D analytical model versus 2D FEA. .............. 160 Figure 3.5 Rotating permanent magnet disk asynchronous levitation machine. (a) Motor topology (with magnetization pattern for example similar to that shown in Figure 1.19a). (b) Cylindrical shell element of magnet disk. ` ................................................................................................... 165 Figure 3.6 3D FEA model, with Nm = 8. ................................................................................... 167 Figure 3.7 Circumferential versus straight magnetization of magnet segments in XY plane. ... 168 Figure 3.8 Levitation force and drag torque, 3D FEA versus analytical shells method for zf =2 mm. .......................................................................................................................................... 168 Figure 3.9 Experimental rotating disk asynchronous levitation machine. .................................. 171 Figure 3.10 Load test setup. ........................................................................................................ 172 Figure 3.11 6-DOF load characteristics for rotating disk asynchronous machine (copper). ...... 174 Figure 3.12 6-DOF load characteristics for rotating disk asynchronous machine (aluminum). . 175 Figure 3.13 Force, torque and mechanical power iso-lines at different flying heights (copper). 177 Figure 3.14 Force, torque and mechanical power iso-lines at different flying heights (aluminum)...................................................................................................................................................... 178 Figure 3.15 Error between analytical shells model and experimental load characteristic (copper)...................................................................................................................................................... 179 xviii Figure 3.16 Error between analytical shells model and experimental load characteristic (aluminum). ................................................................................................................................. 179 xix Notation Multiple forms of vector notation are used throughout this thesis. A general vector field š¾ is denoted by bold italics, while a scalar is always un-bolded. The vector field notations are defined as š¾ ā” šš„ ?Ģ? + šš¦šĢ + šš§?Ģ? ā” (šš„,šš¦,šš§) ā” [šš„šš¦šš§] where (?Ģ?, šĢ, ?Ģ?) are unit vectors in the principal directions of a Cartesian frame (š„, š¦, š§), and the scalar values šš„,šš¦,šš§ denote the magnitudes of the vector at a location (š„, š¦, š§) along each unit direction. xx Acknowledgements Special thanks are owed to my parents, who have supported and sacrificed for me throughout my life. I wish to thank in particular my research advisor, Dr. Xiaodong Lu, without whom none of this would be possible. His incredible work ethic and creative prowess were and are a constant source of inspiration. I thank my friends of old (in alphabetical order): Yoyo Au, Daniel Fritter, Howie Wong, and Adrian Yu. I thank Erina Okuda Nesbit. I thank my somewhat ludicrous number of cousins who made me feel like I had an infinitely extended and close knit family array. My colleagues and friends at UBC who throughout my many, many years in the laboratory have influenced me as an engineer and as a person: Elizabeth Jean Hu, Kristofer Smeds, Darya Amin-Shahidi, Richard Graetz, Arash Jamalian, Bill Kengli Lin, Fan Chen, Alexander Yuen, Navio Kwok, Niankun Rao, Eric Buckley, Mark Dyck, Rui Chen, Jian Gao and Keir Maguire. The machine shop and instrumentation technicians who let me āborrowā all the tools and knowledge I needed to get this done: Roland Genshorek, Erik Wilson, Bernhard Nimmervol, Markus Fengler, Sean Buxton, and Glen Jolly. And of course thanks to Samuel Earnshaw who set us a most interesting problem to solve, and James Clerk Maxwell for giving us the fundamental tools to solve it. xxi Dedication To my parents. 1 Chapter 1: Introduction Non-trivial modern industrial positioning machines typically require controlled motion in multiple degrees of freedom (DOF). For example, a planar stage (X-Y table) generates simultaneous X and Y translation, producing one of the most fundamental machine elements in manufacturing. High speed, high acceleration, precision planar motion with some out of plane Z-stroke is greatly sought after in many high performance industrial applications, such as lithography wafer steppers [4]. Further, even lower precision applications may benefit from ultra-long XY planar strokes as an enabling technology for reconfigurable manufacturing, such as for agile manufacturing in modular assembly systems [5] [6]. An ideal stage concept that simultaneously satisfies the need for precision motion, large strokes and high scanning rates and accelerations is the application of non-contact multi-axis forces to a single end effector, without any intervening connecting structures or bearings which would induce structural modes and decrease overall stiffness, accuracy, and achievable acceleration. This thesis presents two methods for applying such non-contact forces to an end effector over large strokes using magnetic field interactions. First, a novel 6-DOF long stroke synchronous planar levitation motor (Figure 1.1) is presented which uses driven current through stationary 1D coil arrays to actuate sets of 1D Halbach magnet arrays attached to a mover chassis. The motor topology naturally generates linear and decoupled 6-DOF actuation of the mover, which allows minimal controller and drive complexity. The coil structure is very well suited to be manufactured as a printed circuit board (PCB) using commercially available technology without further effort. The passive mover (unpowered and untethered) combined with low heat generation due to good power efficiency from the high copper fill factor of the coil design mean this motor is well suited to applications that require 2 simultaneously high accelerations, precision and very large planar strokes. The coil design and novel driving principles allow for ease of scaling up the planar stroke, relative to other synchronous levitation machines. Second, a novel analysis of an asynchronous planar levitation motor (Figure 1.2) is presented which utilizes sets of mechanically rotating Halbach permanent magnet disks over a conducting stator to produce an induced field in the stator which generates both lift and propulsion on the mover. A major advantage of this class of levitation motor is the simple stator: at a minimum, only a homogenous slab of conducting material without further structure is required, making this stator more cost effective than the synchronous levitation motor stator for ultra-long strokes. Another advantage is that the mover is passively stable in all 6 axes under gravity preload, meaning position feedback is not a requirement. However, the power efficiency of the asynchronous levitation motor is lower than that of the synchronous levitation motor especially at higher loads because of the Ohmic loss in the conducting stator plus power losses incurred by the mechanical rotation mechanism of each magnet disk. The levitation force must be opposed by gravity for stable operation, so the work volume of the motor is limited in orientation. And, while all axes are passively stable, they are also highly coupled, complicating 6-DOF control. Further, the mover requires some powered mechanisms to mechanically rotate the magnet disks relative to the stator and thus needs either onboard power (such as a battery) or an umbilical cable to the mover. This class of levitation motor is thus suitable for applications with lower precision requirements but very large planar strokes. 3 Figure 1.1 Novel 6-DOF synchronous planar levitation motor. Figure 1.2 Asynchronous planar levitation motor mover (onboard battery pack, rotary motor controllers, and stator conducting slab not shown). 4 1.1 Thesis Overview This thesis is divided into four main parts: Chapter 1: Introduction This chapter consists of the overall structure of the thesis and a literature review of both classes of synchronous and asynchronous planar levitation techniques. Chapter 2: Synchronous planar levitation motor This chapter comprises the bulk of the research work and includes: the novel concept and working principles of the 6-DOF synchronous levitation motor; modeling and analysis of the magnetic field and force generation; presentation, analysis and simulation of magnet array designs and design methodology to attenuate force and torque ripple in synchronous machines; design and manufacture of a prototype 6-DOF planar levitation stage; experimental motion tracking results for the prototype 6-DOF stage; and experimental load testing of the force ripple-attenuating magnet array designs. Chapter 3: Asynchronous planar levitation motor This chapter presents the conceptual design of a 6-DOF asynchronous planar levitation motor based on rotating permanent magnet disks similar to that first presented in [7] and more recently as a commercial venture by [8], [9]. Modeling and analysis of the field distribution is presented which yields a closed form analytical solution for levitation and drag forces for a 2D infinitely long linear motor geometry, which is then extended to a new model for levitation force and drag torque for a 3D magnet array disk and conductor geometry. This novel 3D analytical model has 5 the potential to give detailed insight into the effect of various geometric and material parameters on the performance of the levitation motor without the need for extensive 3D finite element analysis. An experimental asynchronous levitation motor element is built to verify the model and levitation force and drag torque are measured. Chapter 4: Conclusion This chapter summarized the contributions of this thesis to the state of the art and discusses future work. 1.2 Thesis Contributions The original contributions of this work are: Synchronous planar levitation motor ļ· A novel electromagnetic configuration for a new type of 6-DOF synchronous planar levitation motor. ļ· Modeling and analysis for magnetic field, force and torque generation, and commutation laws for this new type of planar levitation stage. The novel motor topology allows the minimal complexity analytical force model, which greatly simplifies the commutation law for the motor. A simpler commutation law provides benefits in terms of computational cost in a practical implementation of the 6-DOF motion stage. ļ· Magnet array design methodology for force and torque ripple self-attenuation in this synchronous planar levitation motor, which can be generalized to other synchronous machines. 6 ļ· Experimental load testing validation of force and torque ripple attenuation using novel magnet array designs. ļ· Prototype and motion tracking results for 6-DOF synchronous planar levitation motor. Asynchronous planar levitation motor ļ· New analytical model for levitation force and drag torque of asynchronous planar levitation motor. ļ· Experimental force, torque and power measurement and characterization of asynchronous planar levitation motor element. 1.3 Machine Architectures for Planar Positioning Three types of machine architecture can be used to solve multi-DOF positioning problems: 1) serial, 2) parallel, and 3) direct drive architectures. 1.3.1 Serial Machine Architecture Serial architectures employ multiple 1-DOF motion elements and stack them in series to achieve either multi-DOF motion or enhanced performance (precision and stroke) on the end effector. Figure 1.3a shows a widely adopted approach to multi-axis machines: a Y-axis linear stage transfers its Y-direction motion to an X-axis linear stage, producing combined XY motion on the end effector (the X-axis stage) [10]. The dual stage concept shown in Figure 1.3b is often used to achieve high precision over large strokes by using a long stroke, low precision coarse stage to carry a high precision, short stroke fine stage, achieving high precision over a long stroke [11]. The serial machine architecture is the most ubiquitous because single-DOF stages are the simplest and 7 most mature positioning technology. However, as each 1-DOF stage has only a single actively controlled axis and five passively constrained axes, the final moving stage will suffer from decreased stiffness, accumulated error motion, and usually very large moving inertia caused by the serial connection of multiple moving masses. Figure 1.3 Serial machine architectures. (a) XY H-drive style planar stage. (b) Dual stage. 1.3.2 Parallel Kinematic Architecture An alternative approach to multi-DOF positioning is a parallel architecture (Figure 1.4), where 1-DOF motion elements are connected in parallel between inertial ground and the end effector. This approach is a popular research area due to its versatile motion potential [12] but suffers from limited motion range, varying force/velocity transmission ratio, and position-dependent stiffness within the work volume [13]. Serial and parallel approaches can also be combined together, such as a planar stage with large yaw motion [14]. 8 Figure 1.4 Parallel machine architecture for 6-DOF motion. 1.3.3 Direct Drive Architecture The ideal stage architecture is the direct drive approach (Figure 1.5), where there is only one moving stage/end effector without any intermediate motion elements and all actuating forces are directly applied to the mover. In linear stages, the advantages of direct drive over a lead-screw feed drive are well understood, including improved speed, accuracy and acceleration [15]. These benefits are still true in multi-DOF direct drive stages, but such drives have not been widely utilized. Conceptually, 6-DOF direct drive systems do not need any mechanical guiding bearings, and consequently have additional benefits including zero friction and natural isolation from ambient vibration. In this thesis, we examine in particular methods using magnetic field interactions to produce non-contact forces on the moving body. 9 Figure 1.5 Direct drive machine architecture for 6-DOF positioning stage. 1.4 Multi-DOF Direct Drive Machines Many efforts have been made towards the ideal multi-DOF direct drive concept. 2-DOF direct drives have been widely used in optical disk players as fine motion stages for simultaneous focus control and radial tracking of the optical pickup unit, such as in [16], but such drives only have sub-mm strokes. Combined rotary and linear 2-DOF direct drives can be advantageously used in high precision spindles with in-feed capability [17], [18]. One well-known example of 3-DOF direct drive stages is a Sawyer stepper motor [19], which can generate XY and small Īø in-plane motion on a single stage and relies on air bearings to constrain out-of-plane motion. The disadvantages of this motor include inherent cyclic errors as well as constraints imposed by the umbilical cord which powers the moving stage. For smoother planar motion, voice coil motors [20] and linear motors [21] have been used in high precision 3-DOF direct drive stages, but their strokes are very limited. There are also 3-DOF rotary direct drives for versatile joints or robotic āwristsā, such as variable reluctance spherical motors [22] and brushless spherical motors [23]. A magnetic field-gradient based 5-DOF actuator for micro-needle manipulation has been developed 10 that controls three translation and two rotational degrees of freedom of a 0.5mm long needle-like soft magnetic mover using externally applied and controlled magnetic fields over millimeter level work volumes [24]. Intended for micro-surgery, the work volume and size of the mover needle is severely limited. Fully 6-DOF direct drive stages using magnetic field interactions present particular challenges. Earnshaw in [25] showed that 6-DOF passive magnetic levitation of a permanent magnet without dynamic stabilisation was impossible with only ferromagnetic materials (i.e. with materials with permeability š ā„ šš). An exception are systems that include diamagnetic materials where š < šš ( [26], [27]). Recently, planar actuators utilizing stationary sheets of diamagnetic material (high purity graphite) interacting with permanent magnet movers to generate levitation, in combination with stationary coil sets to provide actuation in the XY plane, have been presented for very small scale movers (such as 2mm by 2mm by 0.4mm) in [28] [29] [30] (Figure 1.6). This planar motion system has the potential to be run entirely without position feedback as all degrees of freedom are passively stable, and can be actuated over very large planar strokes potentially meters-long. However, the low force density of diamagnetic levitation makes scaling up to large payload masses problematic [27]. 11 Figure 1.6 Planar motor for very small scale movers (such as 2mm by 2mm by 0.4mm) based on diamagnetic levitation and Lorentz force translation, figure adapted from [28]. (a) Top view of mover magnets. (b) Electromagnetic configuration of diamagnetic planar levitation motor. For 6-DOF magnetic actuation using materials with š ā„ šš, dynamic feedback stabilization is required in at least one axis. 6-DOF direct drive stages based on single-axis voice coil type Lorentz force drives have been presented in [31], [32], [33]. Verma et. al has demonstrated a 6-DOF stage using multiple 2-DOF Lorentz force motors [34]. While highly linear in terms of force characteristics and thus suitable for very high precision applications, these stage topologies are limited to small positioning ranges. In this thesis we study two categories of long stroke 6-DOF planar levitation direct drive stages: i) synchronous machines where the magnetic field of the stator travels synchronously with the magnetic field of the mover; and ii) asynchronous machines where a changing magnetic field 12 interacts with a conducting body to produce an induction field in the conductor that resists and repels the applied dynamically changing magnetic field producing force on the conductor and dynamic field source. The use of permanent magnets in both classes of motor allows for greater force densities compared to pure current-based electromagnets. The advantage of permanent magnet 6-DOF synchronous levitation motors are the potential for good power efficiency, making this class of motor well suited to high precision applications where thermal variation is a significant issue, as well as vacuum-environment compatibility. The disadvantages are i) position feedback is required for stability and proper synchronous commutation; and ii) coil topologies can become complex or difficult to manufacture and drive, making scaling to ultra-long planar strokes problematic. The advantages of 6-DOF asynchronous planar levitation motors are i) that the conducting body has a very simple homogenous structure (in the case where the conducting body is the stator, the planar stroke of the stage can be extended with minimal cost and complexity); and ii) all axes are passively stable, and therefore position feedback is not a requirement for stable levitation and actuation. The biggest disadvantage to asynchronous levitation is power efficiency even with the use of permanent magnets, as there are significant Ohmic power losses incurred by the induced eddy currents in the conducting body. This makes this class of levitation motor suitable for applications requiring ultra-long strokes with lower precision positioning requirements. 1.4.1 Long Stroke Synchronous Magnetic Planar Levitation Motors Kim and Trumper demonstrated a large stroke synchronous 6-DOF direct drive planar stage demonstrated 50Ć50Ć0.4 mm XYZ stroke and small rotations with a 300Ć300 mm mover in [35]. This synchronous planar levitation motor used four sets of 2-DOF motor elements comprised of 13 Halbach array linear motors and stationary coils based on the Lorentz force for both levitation and translation (shown in Figure 1.7). This design is sufficient for scanning microscopy applications [36]; however, most industry applications require a motion stage with planar strokes on the order of several hundred millimetres or more. Comparing the motor topology in Figure 1.8a and Figure 1.8b, itās clear that in order to scale the planar stroke (šš„ Ć šš¦), the mover size (šæ Ć šæ) must increase in proportion. Simultaneously, the magnet fill factor as a proportion of the mover area decreases meaning lower force densities relative to mover weight and thus lower maximum accelerations for a given current density limit. Achieving >100 mm stroke requires either a mover on the order of meters [37] or this drive system can only be used as a fine motion stage on top of another long-stroke coarse stage. Figure 1.7 Kimās 6-DOF synchronous levitation motor, figure adapted from [35]. 14 Figure 1.8 Scaling up planar stroke of Kim and Trumperās synchronous levitation motor. A moving-coil long-stroke planar motor for photolithography has been presented in [38] and [39]. The basic design is similar to that shown in Figure 1.9, with a set of moving 2D coils actuated above a stationary 2D magnet array. The planar stroke only depends on the size of the stationary magnet array, and therefore very high mover area to working area ratios can be achieved. It is capable of meters-long strokes in X and Y, and can achieve accelerations of several gās with appropriate forced cooling of the moving coils. Such a design can be successfully industrialized and used in advanced lithography machines [4] [40] [41]. However, the stage requires cooling fluid and electrical power lines attached to the mover which causes significant position dependent disturbance forces. In addition, the stator is a large exposed permanent magnet array, which limits the stage to a few niche applications. In addition, the centre of magnetic pressure on each coil set shifts with the mover XY position, creating a position dependent torque acting on the mover. This 15 is the result of synchronous locking of the magnetic field of the moving coils to the magnetic field of the stationary bed of magnets as illustrated in Figure 1.10 with an example 3-phase moving coil levitation motor. Individual force vectors (green arrows) act at the centre of each coil, producing a net levitation force (black arrow) acting at the centre of magnetic pressure for each planar position. As the coils traverse in the XY plane, the centre of magnetic pressure shifts, imparting a position dependent torque to the mover. Figure 1.9 Moving coil synchronous planar levitation motor (umbilical cable to mover not shown), figure adapted from [38]. Red magnet blocks are magnetized in +z, and blue magnet blocks are magnetized in āz. 16 Figure 1.10 Shifting center of pressure of example 3-phase moving coil levitation motor, given a constant levitation force command (umbilical cable not shown). Cable-free long-stroke 6-DOF direct drive planar motion stages have been presented in [42] [43] and [44]. These existing solutions are essentially based on a 2D chessboard-type magnet array on the mover, actuated by a magnetic field generated by a stationary array of 2D patterned coils. Binnard [42] uses multiple layers of racetrack type coils with alternating orthogonal orientations interacting with a 2D magnet array, similar to that shown in Figure 1.11. Compter et. al [44] uses a single layer of racetrack coils with zones of alternating orientation as seen in Figure 17 1.12a. Jansen et. al [45] uses a single layer of herringbone pattern coil arrays (Figure 1.12b), and explicitly excites coil arrays inside the mover area and away from the edges of the magnet array in order to minimize magnet array edge effects. These moving magnet synchronous levitation motors have several significant advantages over moving coil topologies [38] and Kim and Trumperās synchronous levitation motor [35]: i) the passive mover doesnāt require power or a cooling umbilical cable; ii) forced cooling architecture for the stationary coils is easier to implement as the heat source is stationary; iii) the planar stroke can be extended much larger than the size of the mover itself, meaning a greater utilization of the machine volume. However, these 2D moving magnet/ 2D coil array designs have several disadvantages as well: high system complexity, poor coil fill factor, coil end effects, and actuating force and torque variation with stage position. Extensive research effort has been expended on modelling and compensation schemes for these intrinsic force and torque disturbances, for example as presented in [43] [46] [47] [48] [49] [50] [51]. In addition, the number of coils and drives required scales quadratically with the planar stroke (or linearly with the work area): the motion stage in Figure 1.13b has twice the XY stroke of the motion stage in Figure 1.13a but four times the number of individual coils. This significantly and in some cases prohibitively raises the manufacturing, drive and control complexity of the stage as the planar stroke increases. 18 Figure 1.11 Moving magnet 2D chessboard type synchronous planar levitation motor with multiple alternating layers of racetrack style coils; figure adapted from [42]. Red magnet blocks are magnetized in +z, and blue magnet blocks are magnetized in āz. Figure 1.12 Moving magnet 2D chessboard type synchronous planar levitation motors. (a) Alternating zones of racetrack coils, figure adapted from [44]. (b) Herringbone coil pattern, figure adapted from [45]. 19 Figure 1.13 Scaling up XY stroke of moving magnet 2D chessboard type synchronous planar levitation motor (example using motor type similar to [42]). Ueda and Ohsaki presented the design, analysis and simulation of a 6-DOF levitation stage using a 2D moving magnet array and stationary sets of straight coils arranged in alternating layers oriented in X- and Y-axes (similar to a proposal by Kim in appendix A, page 319 of [52]), and also at a skew angle to both planar axes [53]. The X-oriented straight conductors actuate the mover in the Y- and Z-directions, and the Y-oriented conductors actuate the mover in the X- and Z-directions. The skewed coil layer is used to help mitigate the torque coupling with other axes. This stage has the potential for high coil fill factor due to the straight 1D coil array patterns, and for large strokes and ease of scalability since the number of individual coil drives does not scale 20 quadratically with planar stroke. However, the 2D magnet array itself is quite complex to manufacture, and there is significant torque coupling between axes. In addition, in the presented design coils that do not produce force are necessarily active meaning the power efficiency of the motor is compromised. No experimental prototype demonstrating full 6-DOF actuation has been disclosed yet. Figure 1.14 Moving magnet 2D chessboard type synchronous planar levitation motor with alternating layers of straight coils; figure adapted from [53]. A moving magnet long stroke 6-DOF stage has recently been presented in [54], using an array of stationary circular type coils interacting with cylindrical magnets attached to a mover 21 (Figure 1.15). The prototype stage has achieved 80Ć80Ć25 mm XYZ stroke with Ā± 40 degrees and Ā± 15 degrees rotations about the in-plane XY axes, and infinite rotation about the Z-axis. A significant advantage of this stage compared to other synchronous levitation stages is the very large rotation range combined with large planar strokes achievable. However, the force characteristic over the work volume is highly position dependent and non-linear, making control and commutation complex. Additionally, since the number of coils is linear with the planar area of the stage, scaling up the planar stroke presents similar issues to the moving magnet synchronous levitation stages using arrays of 2D racetrack or circular type coil arrays. Figure 1.15 6-DOF moving magnet stage with axially magnetized cylindrical magnets over array of circular coils; figure adapted from [54]. 22 In Chapter 2, we present a new long-stroke 6-DOF direct drive synchronous planar levitation motor solution comprised of sets of 1D magnet arrays and stationary 1D arrays of straight coils. This novel stage has the following advantages over existing long stroke 6-DOF asynchronous planar levitation motors: (1) the mover can be much smaller than the planar stroke (ratio of working planar area to mover area is very high) ; (2) the moving magnet topology requires no cables to the moving body; (3) the number of coils increases linearly with stroke (versus quadratically in existing moving magnet 6-DOF stage solutions) and thus can scale easily to meters-long XY strokes with cm range Z-stroke; (4) the coil terminals exit at the perimeter of the stator, greatly simplifying the system; (5) coil fill factor can be very high and thus large force density is achievable; and (6) it is simple to control/commutate and has superior force linearity due to freedom from coil end effects and permanent magnet field edge effects. 1.4.2 Long Stroke Asynchronous Magnetic Planar Levitation Linear induction motors without permanent magnets have long been studied in the context of generating single axis forces alone [55] but can also be used to generate multi-DOF planar motion [56]. A 2-DOF circular type planar motor using electromagnetic coils and a back-iron was presented in [57] (Figure 1.16) that was shown to be extendable to a long stroke planar motor using linear induction motor elements. A 3-DOF asynchronous XY planar motor using linear electromagnets and a back-iron has been presented in [56], similar to that shown in Figure 1.17. These planar motors still require an out-of-plane bearing due to the field interactions with the back-23 iron, which generates an attractive force that overcomes the repulsive levitation force developed by the dynamic field in the conductor for reasonable driving current densities. The back-iron is required in these types of planar motors to increase magnetic field flux densities and thus increase thrust force densities in the XY plane. Figure 1.16 Circular type 2-DOF induction planar motor, figure adapted from [57]. Figure 1.17 Asynchronous planar motor using electromagnets and back-iron, figures adapted from [56] (out-of-plane bearings not shown). (a) Linear 1-DOF induction motor. (b) 3-DOF configuration of asynchronous planar motor. 24 To generate levitation forces, permanent magnets can be used as the primary dynamic field in the conductor, eliminating the requirement for a back-iron to increase field densities. The Inductrack concept was proposed in [58] that used a linear permanent magnet Halbach array attached to a mover translating over a set of close-packed shorted coils to generate a repulsive levitation force counterbalanced by gravity (Figure 1.18). By using permanent magnets as the primary dynamic field source instead of coils, they were able to demonstrate very high force densities but still require an additional mechanism for propulsion to lift-off speeds. Figure 1.18 Asynchronous levitation motor using permanent magnet Halbach array, adapted from [58]. Previously, [7] presented the concept of a permanent magnet disk rotated mechanically over a homogenous conducting body (Figure 1.19a), producing both levitation force and drag torque that are functions of the rotation speed of the disk. In this thesis we refer to this type of magnetization pattern as circumferential. The drag torque was advantageously used to simultaneously produce a propulsion force by tilting the magnet disk rotation axis or by partial 25 overlap of the magnet disk and the conducting body. This 3-DOF motor element can be used as the basis of a 6-DOF asynchronous planar levitation motor as demonstrated in [3]. A commercial venture based on this concept was recently started as seen in [8] [9]. Another permanent magnet based induction motor element has been proposed that uses a radial pattern magnet wheel rotating over a conducting body (Figure 1.19b) which is also capable of simultaneous levitation and propulsion, and is self-stable when combined with a split-sheet conducting guideway [59] [60]. However, the radial magnet wheel has a lower force density compared to the circumferentially-patterned magnet disk in [7] because only a portion of the strong side of the permanent magnet field is close to the conductor. In this thesis we therefore confine ourselves to the rotating circumferentially-magnetized magnet disk similar to that found in [7]. Figure 1.19 Permanent magnet asynchronous levitation motor elements showing example magnetization patterns. (a) Circumferential magnet disk. (b) Radial magnet wheel. 26 Much of the electrodynamic force modeling of induction effects for both magnet disks and magnet wheels is carried out numerically using finite element type methods [59] [61] [62] [63] [64] [60]. While such numerical methods have been shown to be accurate to within a few percent, a closed form analytical solution for the forces and torques produced on the magnet disk is a very useful tool for informing the optimal design of the magnet disk and conducting slab, and establishing physical understanding and intuition as to motor working principles. Many analytical models for linear induction motors have been given in the literature. The levitation and drag force expressions for an infinitely thick and infinitely wide conducting slab with a time- and space-varying surface current as the primary field source over an infinitely permeable back iron was given in [65]. The analytical model for a linear permanent magnet Halbach array of finite horizontal extent over a uniform conductor was presented in [66], though we note that an expression for the field penetration depth into the conductor (a critical parameter for this model) is not given. A semi-analytical model was presented for eddy current damping generated by cuboidal magnets over a finite conducting plate [67]; this model assumes that the conducting plate thickness is small relative to the skin depth effect and therefore the field fully penetrates the conductor for all speeds. The analytical expression for the field in the air gap of a linear induction motor with a back iron was given in [68], but only applies to current-based dynamic field sources. Lubin and Rezzoug [69] have presented a 3D analytical field, force and torque model for an axial flux eddy current torque coupling. The coupling is comprised of a circumferential permanent magnet disk separated by an air gap from a conducting plate, with some relative rotational velocity between the two bodies. The permanent magnet disk has an iron yoke backing, as does the conducting plate. The purpose of the iron yokes is to maximize the flux density acting 27 through the conductor. Lubin and Rezzoug extrapolate from [70] where it was shown that the curvature of the magnet disk generates second order effects and can therefore be neglected as a first approximation. They accordingly unwrap the magnet disk into an equivalent linear asynchronous machine, and solve for the field distribution using the 3D magnetic vector potential in Cartesian coordinates. By assuming infinite relative permeability of the yoke material, they are able to set the Dirichlet boundary conditions at the edge of the air gap as perfect magnetic boundaries (i.e. no field extends beyond the planar extents of the iron yokes) and solve the boundary value problem. Similarly, Koo et. al. [71] used the presence of back iron yokes in an axial flux permanent magnet disk synchronous rotary machine to set the magnetic boundary conditions in the radial direction, and used the 3D magnetic vector potential in cylindrical coordinates to solve for the no load magnetic flux distribution. The iron yokes in both the asynchronous torque coupler and the synchronous rotary motor are well designed for the particular application: induced drag torque in the eddy current coupler is maximized at reasonable slip speeds [69], and the synchronously generated torque is similarly amplified in [71]. However, iron yokes always produce an attractive force between magnet disk and yoke acting against the direction of the induced levitation force. Therefore for applications where we desire to levitate the magnet disk using asynchronously generated fields, the presence of any iron would be detrimental to generation of lifting forces at reasonable power levels and rotational velocities. Therefore neither of the solutions from [69] or [71] can be used for a practical magnet disk asynchronous levitation motor. A closed form analytical solution for Fujii et. alās 2-DOF circular type planar motor was presented in [57], but is not applicable to permanent magnet-based primary field sources, has an iron yoke and also requires that the area covered by the dynamic field source (the eddy current region) be much smaller than the entire disk. This is non-ideal as it limits the achievable force 28 density and acceleration since the ratio of force generation area to total mover area is small. Park et al. in [72] presented an analytical model for force and drag torque on a rotating permanent magnet disk over a homogeneous conducting stator without iron yokes. This model applies only to conductors of infinite thickness, and uses experimentally derived correction factors on the motor pole pitch and stator conductivity in order to fit the analytical model to experimental results. As a design tool, this empirical approach is limited as extensive mapping of correction factors is required by either finite element analysis or (in Park et al.ās approach) building and testing different magnet array disks. In Chapter 3, we derive a novel analytical expression for the levitation force and drag torque acting on a rotating circumferentially-patterned permanent magnet disk over a stationary homogeneous conducting slab that takes as input only the geometric and material parameters of the permanent magnet disk and conducting slab. Conceptually this method can be extended to modeling propulsion as a function of tilt angle and rotation speed similar to that presented in [7]. 29 Chapter 2: Synchronous Planar Levitation Stage We present in this chapter a novel long-stroke 6-DOF synchronous planar levitation motor which has the following advantages over existing synchronous levitation motors: (1) the number of coils increases linearly with stroke (versus quadratically in existing solutions) and thus can scale easily to meters-long X-Y strokes with cm range Z stroke; (2) the coil terminals exit at the perimeter of the stator, greatly simplifying the system mechanical and electrical connection design; (3) coil fill factor can be very high and thus large force densities are achievable; (4) it is easy to control/commutate and has superior force linearity due to freedom from coil end effects and field edge effects; (5) a minimal machine footprint relative to its working volume. It utilizes multiple 1D magnet arrays (where the permanent magnet field varies along a single axis) and multiple 1D coil arrays (where the current distribution varies along a single axis) to produce 6-DOF motion, making the design, analysis, modeling, and control much simpler than the synchronous planar levitation motors based on 2D magnet arrays and 2D coil arrays. The prototype presented here demonstrates a range of 260 mm by 60 mm, with a mover of 185 mm per side, as shown in the video [73]. This chapter consists of the following parts: ļ· Synchronous planar levitation motor concept, topology and working principle. ļ· Detailed force and torque modeling for 2-DOF motor element. ļ· Magnet array designs for attenuation of intrinsic force and torque ripples. ļ· Design and fabrication details for alpha prototype of 6-DOF synchronous planar levitation motor. ļ· Experimental motion tracking results from prototype. 30 ļ· Experimental load test results for magnet array design for force ripple attenuation. 2.1 Synchronous Levitation Concept and Working Principle Figure 2.1 shows the concept of a novel 6-DOF direct drive motor, which is composed of a stator coil assembly and a mover. The mover travels on top of the stator coil assembly, with a controllable air gap between the two parts. Instead of a 2D magnet array, four 1D magnet arrays (X1, Y2, X3, and Y4) are installed on the bottom surface of the mover. The force actuation principle of the planar levitation stage is the interaction between the magnetic fields of the four magnet arrays bonded to the mover chassis and zones of excitation current in the stationary 1-D coil array to produce a Lorentz force acting on the mover. Each magnet array is a Halbach pattern with a fundamental spatial wavelength of š. The coil array is made of stacks of alternating layers of X-oriented and Y-oriented conductors, with end turns and coil connections made on the outside perimeter of the working planar range of the motor. Each conductor can be excited individually. The active coils (shown in colours) carry currents, while inactive coils (in grey) have no current. The zones of active coils will change dynamically based on the X and Y location of the mover. Figure 2.1 shows the set of conductors that act on each of the four magnet arrays and the generated forces. As the mover traverses the planar range, each excitation zone (X1, X3, Y2 and Y4 active coils) moves synchronously with the corresponding magnet array. This method of force actuation means that the coils can generate a magnetic field wave that can continuously travel in the XY plane. 31 Figure 2.1 Synchronous planar levitation 6-DOF motor electromagnetic topology. Figure 2.2 shows a cross-section of the Y2 magnet array. The Y2 magnet array has the basic pattern of a traditional 4-segment-per- š Halbach array, composed of both vertically and horizontally magnetized elements. Generally, the advantages of a Halbach array over an array of only vertically magnetized elements include 41% higher magnetic field strength and fewer higher order harmonic fields (no 3rd, 7th, or even harmonics) [74]. The active Y-coils are immersed in the magnetic field below the magnet array, and a Lorentz force is generated on the coils from the interaction between the excitation current and the magnetic field. A force acts on the magnet array at the location of the excitation current, leading to both a force and a torque generated on the Y2 magnet array. The active Y-coils thus generate a force vector acting on the Y2 magnet array in the XZ plane and a torque in the Y-direction. The excitation current through each conductor is 32 commutated to generate two independent forces that comprise a single force vector (š¹š2 ?Ģ? + š¹š2?Ģ?) in the XZ plane, where ?Ģ? and ?Ģ? are unit vectors in the X and Z directions, respectively. We note here that the active coils extend beyond the edges of the finite width magnet array in the X-direction further than shown in the schematic drawing in Figure 2.2. In section 2.2, this drive methodology will be shown to be a key aspect of the motor linearity. Figure 2.2 Y2 magnetic array and coils cross-section This force vector (š¹š2 ?Ģ? + š¹š2?Ģ?) on magnet array Y2 is independent of the moverās position in the Y-direction, because the active Y-coils have a uniform structure in the Y direction. As the magnet array width is designed as an integer number of š, the active coils in the orthogonal direction (X direction) have both zero interacting force on Y2 magnet array and moreover no net torque is developed by the orthogonal coils (X-coils) on the magnet array Y2 because the Halbach array pattern is symmetric with respect to the geometric center of the array as shown in [75]. It is made mirror-symmetric about its vertical middle axis by setting the width of the end magnet elements to ļ¬/8 instead of the ļ¬/4 used in conventional Halbach array designs. In addition, the Y-elongated magnet arrays Y2 and Y4 are spaced apart in the X direction to significantly reduce the 33 coupling between the Y2 magnet array and Y4 active coils. These features enable Y2 to mainly interact with Y2 active coils, minimizing the coupling with active coils Y4, X1, and X2 that are commutated for other magnet arrays. Due to the natural force decoupling between the magnet arrays, motion control design can be greatly simplified. Thus the magnet array Y2 and its corresponding active coil zone combine to form a 2-DOF motor element. Similarly, magnet arrays X1 and X3 work with the X coils to generate forces (š¹š1, š¹š1) and (š¹š3, š¹š3) respectively; and the magnet array Y4 produces forces (š¹š4 , š¹š4). The four Z-forces combined can be used to adjust the air gap between the motor and the stator and control the moverās pitch and tilt. The X and Y forces combined control 2D translation and yaw. Thus the 8 independently controlled forces enable 6-DOF motion control of the mover, allowing over-actuation of the mover in translation (š„, š¦, š§) and Euler angles (š¼, š½, š¾) (rotation around š„, š¦, š§, respectively). Figure 2.2 also shows the cross-section of the stator coil assembly with multiple layers of X and Y coils interlaced. This coil array can be made cost-effectively with well-established printed-circuit board (PCB) manufacturing technology. Insulator layers (FR4 core or prepreg) are inserted between coils. The coils in the same direction on different layers can be connected in parallel or series, depending on the via design and end termination. Multiple PCBs can be laterally connected side by side in both X and Y directions (similar to floor tiles) to form a larger planar coil array. The board-to-board lateral connections (in X and Y) are made at the edges only by connecting pads or through-holes of adjacent boards. As a result, all the end terminals of the entire planar coil array are at the outermost perimeter for ease of wiring to the drive electronics. Therefore, this design allows the planar motor to be easily extended in both X and Y directions for various applications while the number of required coils increases linearly with the in-plane stroke, 34 instead of quadratically as in 2D-chessboard type prior art designs. The coil end effects are eliminated, because the current on each layer only changes orientation to return at the outermost perimeter, outside of the active positioning range of the mover. 2.2 Magnetic Field, Force and Torque Generation Modeling for 2-DOF Motor Element The goal of this section is to obtain an analytical model for force and torque generation on the 2-DOF motor element comprised of the magnet array Y2 and its corresponding active coil zone. This force and torque model is applied to all four 2-DOF motor elements and a general plant model is developed in the prototype design section 2.4. We are not necessarily interested in describing the exact distribution of electromagnetic fields at every point in space but rather the net force and torque generated on each 2-DOF motor element from the interaction of the permanent magnet field and the current in the stator. A harmonic field model was presented in [52]. This harmonic field model applies to periodic magnet arrays of infinite horizontal extent and sufficient depth as to render the problem 2D; that is, without considering edge effects or 3D leakage. We desire to use this harmonic model in our case where we have finite width magnet arrays and cannot neglect edge effects out of hand, without losing accuracy or generality. We prove the equivalence between force generated between an infinitely extended repeating magnet array and a single coil, and the force generated between a finite magnet array and a repeating coil pattern, which we term the principle of equivalent force generation. We further show that the extent of the active coil zone for the Y2 magnet array allows us to approximate the infinite repeating coil array case. We thereafter use a partial harmonic field model to derive the force on a Halbach array from its interaction with the current distribution in the active coil zone. 35 This section is comprised of the following parts: 1) General 3D magnetic field distribution for a cuboidal magnet, showing how the exact spatial field distribution, while accurate, gives us little insight into the operating principles of the motor. 2) Force and torque generated by interaction of a permanent magnet field and current distribution. 3) Principle of Equivalent Force Generation, which allows us to use the simplified harmonic model of field distribution for an infinitely extended repeating magnet array under certain driving conditions. 4) The harmonic field model derivation for a Halbach array of infinite horizontal extent. 5) Force and torque generation expressions which we will use to generate commutation algorithms for 6-DOF motor control. 2.2.1 3D Magnetic Field Distribution for Cuboidal Magnet Several analytical models for permanent magnets have been presented in [43] including the magnetic surface charge model and harmonic model. Here we use the analysis for cuboidal permanent magnets based on the fictitious magnetic charge method shown in [76] and [77]. All electromagnetic field quantities for the synchronous planar levitation motor can be modelled under the magnetoquasistatic (MQS) assumption as we can assume that for synchronous driving conditions all such fields will propagate very quickly throughout the volume of the planar stage relative to the time rate of change of driving current and magnet array position for practically sized stages [78]. Under the MQS assumption, the differential form of Maxwellās equations can be written as: 36 š» Ć š¬ = āšš©šš” (2.1) š» Ć šÆ = š± (2.2) š» ā ššš¬ = š (2.3) š» ā š© = 0 (2.4) and š© = šš(šÆ + š“) (2.5) Here š¬,šÆ,š© denote the electric field [š/š], magnetic field intensity [š“/š] and magnetic field density [š] vectors in Cartesian coordinates respectively; š± is the current density vector; š is the volumetric charge density; šš is the permeability of free space and is equal to 4š Ć 10ā7 [š“/š]; and šš is the permittivity of free space and is equal to 8.854 Ć 10ā12 [š¹/š]. The magnetization š“ [š“/š] represents the macroscopic effects of many microscopic magnetic dipoles. Further, the current induced in a body is described by š± = š(š¬ + š Ć š©) (2.6) where š is the conductivity of the material, š is the relative velocity of the medium with respect to the š© distribution. A cuboidal 3D magnet segment is shown in Figure 2.3 with a coordinate system (š„, š¦, š§) centered inside the body of the permanent magnet which has dimensions (2š, 2š, 2š). For this 37 analysis the relative permeability of the entire volume (both inside and outside the magnet) is equal to one, which is a reasonable estimation of the relative permeability of rare earth magnets with šš ā 1.01 ā 1.05. A uniform z-directed magnetization š“ = šš?Ģ? is given where the magnetization amplitude is šš = šµš/šš and šµš is the magnet remanence. Using the magnetic surface charge method [76], [78], a fictitious magnetic surface charge everywhere can be defined as: š = āā ā ššš“ (2.7) Figure 2.3 Cuboidal permanent magnet with vertical (z-directed) uniform magnetization M shown in red. Since the magnetization is uniform inside the permanent magnet, it follows that it is divergence free within the magnet volume and only the top and bottom surface of the cuboidal magnet will have magnetic charge. The magnetic surface charge at the top surface is šš”šš = šššš and the 38 magnetic surface charge at the bottom surface is šššš”š”šš = āšššš. In the absence of currents, the curl of the field intensity šÆ is zero and can be solved using the scalar magnetic potential: šÆ = āāš (2.8) From [78], with only the magnetic surface charges as sources, the scalar potential can be formulated as š(š) = ā«š(š)šš£4ššš|š|š (2.9) where š is a point in the (š„, š¦, š§) frame centered on the magnet body. Thus from (2.7), (2.8) and (2.9) the field distribution can be analytically solved as: šµš„ =šµš4šāāā(ā1)š+š+ššæš(š ā š)1š=01š=01š=0 (2.10) šµš¦ =šµš4šāāā(ā1)š+š+ššæš(š ā š)1š=01š=01š=0 (2.11) šµš§ =šµš4šāāā(ā1)š+š+ššš”šš2(ššš š)1š=01š=01š=0 (2.12) where the magnetic flux density vector is š©šššššššššš = šµš„ ?Ģ? + šµš¦šĢ + šµš§?Ģ?, and š = āš2 + š2 + š2 (2.13) š = š„ ā (ā1)šš (2.14) š = š¦ ā (ā1)šš (2.15) 39 š = š§ ā (ā1)šš (2.16) and atan2 is a four quadrant arctangent operator [77]. Since the problem is stated as linear with uniform relative permeability everywhere, we can obtain the full field distribution for a cuboidal magnet pattern by simply translating and rotating the field for each cuboidal magnet segment and using superposition to obtain the net magnetic field. Thus we can solve for the exact spatial distribution of any pattern of cuboidal magnet segments. As an example, we take a Halbach-patterned magnet array as shown in Figure 2.4a made up of five cuboidal magnet segments with X- and Z-directed magnetizations only and finite extensions in each direction. The Halbach pattern has a repeating pattern every š, the fundamental spatial wavelength for the magnet array. Figure 2.4b shows the magnetization pattern with segment dimensions. Each cuboidal magnet segment is assumed to be uniformly magnetized in the direction shown within the magnet volume. This is a close approximation of the real case, however, in reality the local magnetization will be influenced by the magnetic fields generated by the neighbouring magnet segments (the so-called ābucking fieldā). Setting magnet remanence šµš = 1.325 [š], array width šš = š and array depth š·š = 2š, the analytical field š©šššššššššš is calculated as per (2.10)-(2.12). Figure 2.5 shows the vector component of flux density š©šššššššššš at š§š = āš/5. As a partial validation of the analytical 3D field model, commercial physics modeling software COMSOL in 3D mode [79] is used to determine the 3D magnetic field for the same magnet array via finite element analysis (FEA) yielding š©šš¬šØ. The difference between the analytical model and FEA is shown in Figure 2.6, with error defined as š©šš¬šØ ā š©šššššššššš. The worst case error as a percentage of the maximum amplitude of the analytical field for this š§š-level is 1.1% for šµš„, 0.48% for šµš¦, and 1.37% for šµš§. The solutions converge with finer meshing in the finite element 40 model; however due to computational hardware limits the FEA mesh was not further refined. The analytical model coded in the MATLAB environment gave a result in 7.8 seconds on an Intel Core i7-3770 CPU at 3.4 GHz, 12 GB of DDR2 RAM, running Windows 7 64bit. The FEA model took approximately 3 hours to yield a solution on the same computer. Figure 2.4 Example 3D finite extension Halbach-patterned magnet array (magnetization directions shown in red). (a) 3D isometric view of finite Halbach array. (b) 2D view of finite Halbach array showing segment dimensions. 41 Figure 2.5 Analytical field model, flux density vector components of š©šš¬šØ for example magnet array at šš =āš/š. 42 Figure 2.6 Error in flux density vector components for šš = āš/š. Error defined as š©šš¬šØ ā š©šššššššššš. 43 We therefore have a 3D analytical field model that captures essentially all characteristics of the spatial distribution of any cuboidal magnet pattern given uniform permeability through the solution space and a known magnetization distribution. Here we note that this analytical model was determined for magnetizations that are parallel with two faces of the cuboidal magnet geometry. However, despite having an exact analytical expression for the field distribution for a finite Halbach array with a very significant improvement in solution time compared to 3D finite element models, we gain little design insight into the operation of the motor. From the summed analytical expression for the total field it is non-trivial to determine a commutation law for multi-DOF force and torque actuation of the Halbach array using only this form of the field distribution [77]. Sensitivity analyses can be carried out on each parameter to determine an optimized motor design [80]. We will show in later sections that the topology and driving method of this synchronous levitation motor lends itself to a simplified analysis for force and torque generation. In the next section we determine the force and torque generated from interaction between a current distribution and the magnetic field from a general magnet array. 2.2.2 Force and Torque Generation on Magnet Array Once again assuming the relative permeability of all materials in the motor is close to unity, the force generated on the magnet array from the interaction of the current distribution in the stator š± and the magnetic field from the Halbach array can be calculated from the Lorentz force expression alone [81]. As illustrated in Figure 2.7, the force acting on the magnet array due to a differential conductor volume with current density š± is 44 š ššššššš = ā(š± Ć š©š)šš£ (2.17) Here š©š is the magnetic field due to the magnet array alone, š± is the applied current density in the stator and šš£ is the volume of the conductor element. The total force on the magnet array is therefore the integration of (2.17) over the entire conductor volume šš: ššššššš = āā (š± Ć š©ššš)šš£ (2.18) The total torque on magnet array is similarly š»šššššš = āā š Ć (š± Ć š©ššš)šš£ (2.19) where š is the moment arm vector between the magnet array center of gravity and the location of the differential conductor volume. Figure 2.7 General 3D magnet array interacting with differential conductor element. 45 2.2.3 Principle of Equivalent Force Generation From section 2.2.1 and the force and torque expressions (2.18) and (2.19) it can be seen that solving for the magnetic field distribution for a 3D Halbach array (or any repeating pattern magnet array) with finite extension in each dimension yields a complex and non-intuitive analytical expression for the magnetic field spatial distribution, which must then be integrated across the volume of the conductor. While this integration is at least numerically solvable, obtaining the exact spatial distribution for such a magnet array is not the primary goal of this analysis nor is it the most useful magnetic field model for design insight into motor operation. Rather, we are interested in the net force and torque generated on the array via the interaction of a repeating pattern current distribution and the permanent magnet field from the Halbach array. In this sub-section we show that the force generated on a Ī»-wide Halbach magnet array by a coil array repeating with the same period of Ī» is equivalent to the force generated on an infinitely wide Halbach magnet array interacting with a single coil driven by the same current. We later extend this analysis to examine the effect of different spatial wavelengths for the repeating coil pattern versus the width of the magnet array. An example Halbach magnet array is shown in Figure 2.8, with equivalent dimensions to those shown in Figure 2.4 and x-directed width šš = š (hereafter referred to as the š-width magnet array). Note the magnetization pattern is flipped compared to the array in Figure 2.4 though the strong field side remains below the magnet array (š§š < 0). Figure 2.9 shows a Halbach magnet array with the same spatial wavelength and repeating magnetization pattern, with infinite x-directed extension šš = ā (hereafter referred to as the ā-width magnet array); other dimensions are equal to the š-width array in Figure 2.8. 46 Figure 2.8 Finite x-directed extension Halbach-patterned magnet array. Figure 2.9 Infinite x-directed extension Halbach-patterned magnet array The z-directed field component šµš§ at the test line shown in Figure 2.8 (š§š = āš/5, š¦š = 0) is calculated for the š-width magnet array using (2.12) for a š„š ā [ā3š, 3š] (with šµš§ falling to 47 0.0152% of the max|šµš§| at these boundaries). In order to calculate šµš§ over the same x-range for the ā-width magnet array, we note that magnetic field values converge far away from the edges of the magnet array as šš increases. That is, šµš§ for šš = 500š is that same as that for šš =1000š over š„š ā [ā3š, 3š] to within 2.1 Ć 10ā5 ppm (relative to the maximum amplitude of |šµš§|). Therefore for the x-range of interest, we consider šš = 1000š equivalent to a true ā-width magnet array. Figure 2.10 shows šµš§ for both the š-width and ā-width magnet arrays. The š-width magnet array clearly shows edge/leakage effects and no strong periodicity. The ā-width magnet array has a repeating pattern with a fundamental wavelength of š. At š§š = āš/5, the field from the ā-width magnet array is almost entirely a single fundamental. This in no way reduces the generality of the force analysis which applies to force generated by coils at any š§š. Figure 2.10 š©š from 3D analytical model along test line (šš = āš/š, šš = š). 48 For the synchronous planar levitation motor 2D forcer, the force generation case is a repeating pattern of y-directed coils immersed in the magnetic field of a finite width magnet array. In Figure 2.11a, a single Ī»-wide Halbach magnet array M0 interacts with a group of y-directed coil traces Ck (k = āā,ā¦ ,ā2,ā1,0,1,2,ā¦ ,ā) which repeat every Ī» and extend from š¦š ā (āā,ā). The coils are simple line elements and have zero size in the XZ plane, reducing the force expression (2.18) to a single spatial integration along the y-axis. The x-locations of the coil traces in the magnet frame are š„š = šæ for coil C0, š„š = šæ + š for coil C1, š„š = šæ ā š for coil Cā1, and so forth. Correspondingly, assuming that the z-directed flux density distribution from the magnet array M0 is šµš§(š„š) in Figure 2.11a, the field at each coil location is šµš§(šæ) for coil C0, šµš§(šæ + š) for coil C1, šµš§(šæ ā š) for coil Cā1, and so forth. As a result, if each coil Ck is driven by the same current š [š“] the total x-directed force as determined from the x-component of (2.18) on the magnet array is š¹š„,0 = āā« ā šµš§(šæ + šš) āš=āāš šš¦āāā (2.20) The order of operations for the discrete summation and continuous integration are interchangeable. Figure 2.11b shows an alternate case, with a single line-element coil C0 interacting with an infinitely extended magnet array which is the combination of the magnet array M0 and its duplicates Mk (k = āā,ā¦ ,ā2,ā1,0,1,2, ā¦ ,ā) which repeat every š. The magnet array M0 and the coil C0 in Figure 2.11b are at identical locations to those in Figure 2.11a. The magnetic field from magnet array M0 at the location of C0 is šµš§(šæ), from magnet array M1 at the location of C0 is šµš§(šæ + š), from magnet array Mā1 at the location of C0 is šµš§(šæ ā š), and so forth. Therefore the total field at the location of C0 is the superposition of the fields from each magnet 49 array, šµš§(šæ ā ā) + āÆ+ šµš§(šæ ā š) + šµš§(šæ) + šµš§(šæ + š) + āÆ+ šµš§(šæ + ā). The corresponding x-directed force on the magnet array is therefore š¹š„,ā = āā« š ā šµš§(šæ + šš) āš=āāšš¦āāā (2.21) Comparing the forces (2.20) and (2.21), the x-directed force between a group of coil traces spaced š apart in š„š and a single š-wide magnet array is equivalent to the x-directed force between a single coil trace and an infinitely wide magnet array with a magnetization pattern that repeats every š. It can be similarly shown that the z-directed force is also equivalent between the two cases. It is clear that the repeating pattern magnet array in Figure 2.11b is exactly equivalent to the ā -width magnet array as shown in Figure 2.11c. Since the field from the šš-wide magnet array goes to zero very rapidly as š„š ā Ā±ā, the contribution from coils far away from the edges of the magnet array in Figure 2.11a donāt contribute much to the total force on the magnet array. The field is significantly attenuated (below 1% of its maximum value as shown in Fig. 11 for the Ī»-wide array) beyond Ī»/2 from the edges of the magnet array. Practically, only coils immediately under the magnet array or extending beyond the edge of magnet array less than Ī»/2 need to be excited in order to consider the magnetic field as purely periodic from the point of view of net force generation. We therefore make the approximation that so long as the extent of the active coil zone for each magnet array is at least Ī»/2 from the edges of the magnet array in š„š, the principle of equivalent force generation applies. The force equivalence analysis thus far has been for the special case of a repeating coil pattern with a spatial periodicity equal to the width of the magnet array, which is in turn an integer value of the Halbach fundamental wavelength š. A more general case is shown in Figure 2.12a, where the magnet array width is šš and the coil array spatial periodicity is šš. For šš ā¤ šš, the 50 superposition method implies that each identical magnet array will have some overlap š where š =šš ā šš as shown schematically in Figure 2.12b. As determined in the motor topology section 2.1, the width šš of the magnet array should be constrained to an integer value of the Halbach spatial period, and the coil pattern should repeat with a spatial wavelength of š for maximum force generation. Therefore, a well-designed coil and magnet array will be constrained to šš = šš and šš = š where š is preferred to be a real positive integer. From these constraints a general force equivalence can be obtained where a Halbach magnet array of width šš = šš interacting with a repeating coil array with spatial wavelength šš = š has the same force as a single coil interacting with š- ā-width magnet arrays superimposed in the same location. Thus the force expression (2.21) becomes š¹š„,ā = āšā« š ā šµš§(šæ + šš) āš=āāšš¦āāā (2.22) and š¹š§,ā = š ā« š ā šµš„(šæ + šš) āš=āāšš¦āāā (2.23) where š = šš/šš. Therefore the force on the 2D motor element can be determined by analyzing the magnetic field of an ā-width magnet array interacting with a single coil and multiplying by the magnet array width factor š = šš/šš. We note here that for non-integer š, while the mean force generated from the interaction with a repeating coil array will still be proportional to š, additional force harmonics will be introduced which may have to be addressed depending on the performance required. For the general case, the practical excitation zone remains one the order of š/2 away from the edges as the decay rate in š„š is dependent on the Halbach period and not the width of the magnet array. 51 We conclude that the net force generated on the 2-DOF motor element can be appropriately analyzed by collapsing a group of coil traces interacting with a finite width magnet array to a single coil interacting with an infinitely extended magnet array. It remains to determine an appropriate and useful analytical model of the field from an ā-width magnet array. 52 Figure 2.11 Equivalent force cases. (a) š-width magnet array interacting with š-repeating coil array. (b) repeating š-width magnet arrays with single coil. (c) ā-width magnet array with single coil. 53 Figure 2.12 (a) š¾š-width magnet array interacting with š¼š-repeating coil array. (b) Superposition showing overlap of each identical magnet array. 54 2.2.4 Magnetic Field Modeling for Infinitely Extended Halbach Array Sub-section 2.2.3 established that modeling the net force on the 2D motor element as a series of single coils interacting with an infinitely extended array is appropriate so long as the excitation zone extends sufficiently beyond the edges of the finite magnet array. In this sub-section we develop a simple analytical model to analyze the strong-side field of an infinitely extended Halbach array. This analysis is based on the original work carried out in [52] where the harmonic model for a general magnet array was presented. Kim formed the problem as 2D in the XZ plane, assuming that the depth of the magnet array was sufficient to neglect leakage effects. Here we begin with the 2D assumption and compare the resulting harmonic model to the exact 3D analytical model integrated over š¦š ā (āā,ā) for an infinitely extended array as presented in section 2.2.3. Though we assume no current in any region, in Chapter 3 we will use this analysis as a starting point to determine the force from the induction levitation effect on the asynchronous planar levitation motor, with induced currents in the stator region. We therefore solve field distributions using the vector magnet potential rather than with the scalar magnet potential because unlike the scalar magnet potential the vector magnet potential applies to each region regardless of the presence of currents, i.e. where the curl of the field intensity šÆ is not zero. The field density distribution must satisfy Gaussā law for magnetic field (2.4) everywhere, i.e. the field density š© is everywhere solenoidal (divergence free). One solution is š© = [šµš„šµš¦šµš§] = ā Ć šØ (2.24) where šØ is the vector magnet potential in Cartesian coordinates. This automatically satisfies Gaussā law (2.4) since the divergence of the curl of a vector field is always zero. Plugging into Ampereās law (2.2), and substituting field intensity using (2.24), we obtain 55 ā Ć ā Ć šØ = šš(š± + ā Ć š“) (2.25) Using the vector identity ā Ć ā Ć šØ = ā(ā ā šØ) ā ā2šØ and setting the Coulomb gauge ā ā šØ = 0 (i.e. the vector potential is divergence free, which is correct for the MQS assumption), we obtain the vector Poissonās equation: ā2šØ = āšš(š± + ā Ć š“) (2.26) This vector Poissonās equation is equivalent to three scalar Poissonās equations. From the Biot-Savart law and the definition of the vector magnet potential (2.24), it is possible to write the magnetic vector potential in integral form in the absence of free surface currents as [52] : šØ =šš4šā«š± + ā Ć š“š ā šā²šš£ +šš4šā«š“ Ć šš ā šā²šš šš (2.27) Here the volume š is enclosed by the surface š, the source location is šā² and the observer location is š. The second term in (2.27) insures continuity of the vector potential across the boundary between permanent magnet and non-permanent magnet. Boundary conditions are given by (2.2), (2.3) and (2.4), and the continuity of the vector magnet potential across each boundary. Figure 2.13 shows the solution volume broken up into three regions: region I-the volume of air on the top (weak side for a Halbach array) of the magnet array, region II-the permanent magnet array with magnetization š“ = šš„ ?Ģ? + šš§?Ģ? which has only š„š-dependence within the magnet volume (this is assumed to be rare earth magnetic material and therefore šš ā 1), and region III- the air gap between magnet array and current distribution. In the 2D limiting case the vector potential šØ has only a y-directed component, šØ = š“š¦š.Ģ For the general magnetization š“ and in the absence of currents, the vector Poisson equation (2.26) reduces to the scalar equation 56 (ā2ā2š„+ā2ā2š§)š“š¦ = šš(āāšš„āš§+āšš§āš„) (2.28) Figure 2.13 2D form of the infinite extension magnet array. The magnetization distribution š“ can be written as a sum of Fourier coefficients as described by [52]. The field due to each magnetization harmonic can be solved separately and summed for the total field solution, and therefore the vector potential can be treated similarly. Therefore in each region we can assume that the form of the vector potential will satisfy the usual Laplace function criterion that the vector potential will be the multiplication of functions that depend only on a single coordinate: š“š¦ = š¹(š„š)šŗ(š¦š)š»(š§š) (2.29) 57 Given a sinusoidal variation in two of the functions, the third must be dependent on hyperbolic trigonometric functions. For a four-segment per spatial wavelength Halbach array of infinite width in the x-direction, the magnetization distribution is shown in Figure 2.14 [74]. The gap š is introduced between magnet elements to maintain generality and allow for later analysis of the effects of practical arrays. The magnetization amplitude is šš = šµš/šš where šµš is the magnet remanence. We define the characteristic harmonic wavelength ššš = š/š2š where š is the harmonic wave number with š > 1, and fundamental characteristic wavelength šš = š/2š. The Fourier coefficients are calculated as šš„š =1šā« šš„šāšš„š/ššššš„šš/2āš/2 (2.30) šš§š =1šā« šš§šāšš„š/ššššš„šš/2āš/2 (2.31) In the remaining analysis we solve the field for the fundamental component of magnetization. All field components can be similarly solved by substituting harmonic magnetization components calculated as per (2.30), (2.31). The total field is the summation of all field harmonics and the fundamental. The magnetization fundamental harmonic in region II is š“š = [šš„10šš§1] = š š [šš1šāšš„š/šš0š1šāšš„š/šš] (2.32) with the fundamental spatial wavelength š and the magnetization amplitude š1 =4šššcos (š4+š2šš) (2.33) 58 Figure 2.14 Magnetization distribution within magnet volume for ā-width Halbach array with 4 magnet segments/š separated by gap š. Since the Poisson equation (2.28) reduces to a pure Laplace equation in region I, the fundamental component of vector potential in region I can have the form: š“š¦1 = š š([šš¼šš§ššš + šš¼šāš§ššš ]šāšš„š/šš) (2.34) Similarly in region III, the vector potential has the form: 59 š“š¦1 = š š([šš¼š¼š¼šš§ššš + šš¼š¼š¼šāš§ššš ]šāšš„š/šš) (2.35) Since the potential must go to zero at š§š = āā,ā, šš¼ = šš¼š¼š¼ = 0. In region II, the Poisson equation (2.28) is driven by the magnetization term. Since the magnetization is given as uniform in š§š within the magnet volume, (2.28) further reduces to (ā2ā2š„+ā2ā2š§)š“š¦1 = ššāšš§1āš„ (2.36) The solution will have a homogenous part which solves (ā2ā2š„+ā2ā2š§)š“š¦1ā = 0 and allow us to fit boundary conditions and a particular solution driven by (ā2ā2š„+ā2ā2š§)š“š¦1š = ššāšš§1āš„. Since šš§1 has the x-dependence šāšš„š/šš , the particular solution must satisfy (ā2ā2š„+ā2ā2š§)š“š¦1š =āššššššš§1 (2.37) From [52] the vector potential along each boundary in region II has the form š“š¦1 = š“š¦1š + (š“š¦1,š ā š“š¦š)š ššā (š§ššš)š ššā (š»ššš) ā (š“š¦1,š ā š“š¦1š)š ššā (š§š ā š»ššš )š ššā(š»ššš) (2.38) where we assume that š“š¦1 has the x-dependence šāšš„š/šš, and š“š¦1,š and š“š¦1,š are the potentials along the š- and š- boundaries respectively. The particular solution š“š¦1š is obtained by solving (2.37) and noting that since š“š¦1š must also have the x-dependence šāšš„š/šš, then ā2š“š¦1šā2š„=ā1ššš2 š“š¦1š. Further, since the magnetization is uniform in z within the magnet volume, the particular solution is 60 š“š¦1š = ššššššš§1 (2.39) From (2.24), the field components as a function of the vector potential are š© = [šµš„šµš¦šµš§] = š š[ āāš“š¦āš§0āš“š¦āš„ ] (2.40) Using the definition of š“š¦1 in each region, and the definition of flux density (2.43), the value of field densities at each boundary are obtained: š©š = š š[ 1šššš¼šāš»š/šššāšš„š/šš0āššššš¼šāš»š/šššāšš„š/šš] (2.41) š©š = š š[ 1šš{ā(š“š¦1,š ā š“š¦1š) ššš”ā (š»ššš) + (š“š¦1,š ā š“š¦1š)šš šā (š»š/šš)}0āšššš“š¦1,š ] (2.42) š©š = š š[ 1šš{ā(š“š¦1,š ā š“š¦1š) šš šā (š»ššš) + (š“š¦1,š ā š“š¦1š)ššš”ā (š»š/šš)}0āšššš“š¦1,š ] (2.43) š©š = š š[ ā1šššš¼š¼š¼šāšš„š/šš0āššššš¼š¼š¼šāšš„š/šš] (2.44) 61 Boundary continuity conditions are derived from the continuity of vector potential, field continuity due to Gaussā law (2.4), and from (2.2) and (2.5) which in the absence of currents yields ā Ć š© = ššā Ć š“. (2.45) Boundary conditions are summarized as šµš„,š ā šµš„,š = āšššš„1,š (2.46) šµš§,š ā šµš§,š = 0 (2.47) šµš„,š ā šµš„,š = šššš„1,š (2.48) šµš§,š ā šµš§,š = 0 (2.49) š“š¦,š = š“š¦,š (2.50) and š“š¦,š = š“š¦,š (2.51) where the subscript denotes vector amplitude direction and the boundary location. In addition, the field š© as well as the vector potential šØ go to zero at š§š = āā,ā. šš„1,š and šš„1,š are the fundamental components of magnetization in the š„š-direction along the š- and š- boundaries respectively. We are particularly interested in region III below the magnet array, as this is where the stator coils will be located. Solving the system of equations (2.41)-(2.44) using the boundary 62 conditions and relations (2.46)-(2.51), the fundamental component of magnetic field in region III is found to be š©š°š°š°,š = š š[ āšššš1 (1 ā šāš»ššš ) šš§š/šššāšš„š/šš0ššš1 (1 ā šāš»ššš ) šš§š/šššāšš„š/šš] (2.52) For the following sections, the primary harmonics of interest are those present in š©š°š°š°, in particular the fundamental harmonic with characteristic wavelength šš =š2š and the fifth harmonic with characteristic wavelength šš5 =15š2š As discussed in [52], for a four-segment-per- š magnet array with an infinite number of periods, the magnetic field can be modelled as the sum of a base component at 1 cycle-per-Ī» (cpĪ») and harmonic components at 5, 9, and 13 cpĪ». Further harmonics can be added depending on the accuracy desired. We note that at š§š = āš/60, the 13th field harmonic amplitude is less than 0.1% of the fundamental component, and that the harmonics decay faster as a function of the harmonic number. From (2.52) the fundamental flux density is š©š°š°š°,š = [āšµ1šš§š/šš1š šš (š„š/šš)0šµ1šš§š/šš1ššš (š„š/šš)] (2.53) where šµ1 =4šµšcos (š4+š2šš)š(1 ā šāš»ššš ). Repeating the above analysis for the fifth harmonic of magnetization, the 5 cpĪ» component of flux density is š©š°š°š°,š = [šµ5šš§š/šš5š šš (š„š/šš5)0āšµ5šš§š/šš5ššš (š„š/šš5)] (2.54) where šµ5 =4šµšcos (š4+š2šš5)5š(1 ā šāš»ššš5). We note that the full field is the summation of all harmonic contributions plus the fundamental 63 š©š°š°š° = š©š°š°š°,š + ā š©š°š°š°,šāš=2 (2.55) We use this result in the following sections to determine both the force coefficient and force ripple due to the fundamental field component and the fifth field harmonic. We first derive the force ripple due to a single field harmonic and later generalize to force ripple due to all other harmonics. 2.2.5 Force Due to Single-Phase Repeating Coil Group Now that we have determined an appropriate analytical model for the field from an ā-width Halbach magnet array, we can determine the force on the šš-magnet array shown in Figure 2.15 by the group of coils identically driven by the current š¼š and spaced apart by šš = š (implying that the current distribution in the coil array has a spatially periodic current distribution). Each coil trace has width šš and thickness š”š, with the top surface located at š§š = āš§š. For this analysis we assume the depth of the magnet array is sufficient to ignore leakage flux and the integral expression ā« š¼ššš¦šāāāā š·šš¼š. Figure 2.15 Single phase coil array below finite width Halbach array. 64 From the general Lorentz force (2.18) and the equivalent force expressions (2.22) and (2.23), combined with the field expression for the ā-width magnet array (2.55), the total force from the total field on the magnet array due to the group of coils is šš =ššššš·š ā« ā« (š½ššĢ Ć š©š°š°š°)šš„ššš§ššš/2āšš/2āš§šā(š§š+š”š) (2.56) where the current density is š½š =š¼šššš”š[š“š2]. Note that for the case shown in Figure 2.15, šššš= 2. As we will show section 2.2.6, the fundamental component of š©š°š°š° contributes to the desired mean force generation on the mover and all other field harmonics interact with stator currents to generate undesired force and torque ripples. In particular, the 5-cpĪ» field harmonic combined with a 3-phase coil pattern produces a 6-cpĪ» force ripple. Theoretically, the higher harmonic field components can be eliminated by using continuously varying Halbach magnet arrays as described in [74]; however, these are not practical to manufacture. To develop a means to minimize the effect from harmonic field components, here we use a partial harmonic model which includes both the fundamental component and a 5th harmonic component. Thus, given a repeating coil pattern, for the purpose of deriving the net force on the magnet array the spatial variation of the magnetic field below the magnet array can be approximated as the partial harmonic model š©š·šÆš“ = š©š°š°š°,š+š©š°š°š°,š for š§š < 0 for the 2D limiting case: š©š·šÆš“ =[ āšµ1šš§ššš š ššš+ šµ5š5š§ššš š šš(5š) 0šµ1šš§ššš ššš š āšµ5š5š§ššš ššš (5š) ) ] . (2.57) Here we introduce the electrical angle š =š„ššš1[ššš] and šµ1, šµ5 are defined as in (2.53), (2.54). Higher order spatial harmonics decay faster than the fundamental harmonic as š§š decreases (the decay rate being proportional to š). For a change in š§š from š§š = āš/60 to š§š = āš/5, the 65 fundamental field harmonic decreases by 68%, whereas the 5th field harmonic decreases by over 99%. Figure 2.16a shows the comparison between this PHD and the exact 3D analytical model (3DAM) solution for an ā-width magnet array with depth š·š = 20š evaluated at š¦š = 0 (middle of the array in š¦š). The error is defined as PHD-3DAM. The field at š§š = āš/60 is highly non-sinusoidal, while at š§š = āš/5 the fundamental harmonic dominates. The maximum modelling error at š§š = āš/60 is 11.2% relative of the fundamental component, and this error is dominated by the 9th, 13th and 17th harmonics as shown in Figure 2.16b, with 6% due to the 9th harmonic alone. At š§š = āš/5, the maximum modelling error falls to 0.0006% relative to the fundamental component, with the dominant error harmonic at 9- cpĪ» making up essentially all of it. This is the result of the higher harmonics decaying faster than the 9th harmonic. If required, at lower flying heights, 9th, 13th and 17th field harmonics can be added into the PHM to reduce modeling error. However, as shown in the following analysis, with a 3-phase coil arrangement the dominant force components are due to the fundamental and 5th field harmonics for practical stator coil geometries with the effects of the 9th field harmonic entirely eliminated, and therefore the analytical field model is usefully simplified to only these harmonics. 66 Figure 2.16 Partial harmonic model versus 3D analytical model at different šš (a) š©š, partial harmonic model (PHD), 3D analytical model (3DAM) evaluated at šš = š with š«š = ššš, and error (PHD-3DAM), for ā-width magnet arrays. (b) Fourier coefficients for error at different šš. 67 The exact force expression (2.56) can be approximated using the PHM field as šš,š·šÆš“ =ššššš·š ā« ā« (š½ššĢ Ć š©š·šÆš“)šš„ššš§ššš/2āšš/2āš§šā(š§š+š”š) (2.58) Evaluating the spatial integrations across the coil geometry, the field on the magnet array is therefore šš,š·šÆš“ = š¼ššššš[āš¾1ššš ššš§ššš + š¾5 ššš (5š) š5š§ššš }0āš¾1 š šš š šš§ššš + š¾5 š šš(5š) š5š§ššš] (2.59) where magnet and coil geometry are combined in the factors š¾1 = šµ1š·š š ššš (šš2šš) (ššš”š) (1 ā šāš”š/šš) (2.60) š¾5 = šµ5šæš¦š ššš (5šš2šš1) (šš15š”š) (1 ā šā5š”š/šš1) (2.61) The torque developed on the finite magnet array is more complex as there is no analogue to the equivalent force generation principle for torque; rather the exact field for the finite width array as calculated by the 3DAM must be used to determine the torque from each coil individually and summed to find total torque: šš¦(š„š, š§š) = ā šš,šāā āā āā Ć š¹š,š āā āā āā ā āš=āā (2.62) where šš,šāā āā āā is the individual moment arm from the center of gravity of the magnet array to the center of the current distribution of coil trace š (assumed to be at the center of the coil trace for all practical driving conditions) as shown in Figure 2.15, and the force š¹š,š āā āā āā ā is the individual force from coil trace š on the magnet array. 68 2.2.6 Force Due to 3-Phase Repeating Coil Group with Commutated Current From (2.59), a constant current š¼š with the magnet array traversing the X-direction produces a rotating force vector in the XZ plane with fundamental and 5th spatial harmonics. Noting that the force expression is sinusoidal with a fundamental wavelength š (a result that was impossible to determine using only the exact field distribution model 3DAM for the finite width array), a commutation law for the coil group can be determined. Since the force on the magnet array is modulated by the sinusoidally varying field, in order to generate a specified force as the mover traverses the XZ plane, the current š¼š is modulated to have the form š¼š = āš¼š„š šāš§ššš ššš š ā š¼š§ššāš§ššš š šš š (2.63) where š¼š„š and š¼š§š represent the force commands on the magnet array for the š„ and š§ directions respectively. This commutation law is function of the coilās location (š„š, š§š) in the magnet frame; accordingly, another active coil trace located an integer number of š away in the x-direction will be excited with the same current. By changing the current š¼š sinusoidally with electrical angle š such that the fundamental harmonic of current is synchronous with the moverās permanent magnet field at the location of the conductor, š¼š will demodulate the permanent magnet field to produce a useful position independent or mean force, as well as 2, 4 and 6 cpĪ» force ripples when we consider only the fundamental and 5th harmonics of field. To obtain the maximum mean force, it is clear that the coil array spatial period of šš must match the magnet array spatial period š. The mean force and the 2 cpĪ» force ripple are from the modulation of the fundamental harmonic of the magnetic field with the fundamental spatial harmonic of excitation current. The 4 and 6 cpĪ» force ripples in force are due to the modulation of the 5th field harmonic with the fundamental harmonic in the excitation current. If the force ripple due to the 9th field harmonic were calculated, it would 69 produce 8 and 10 cpĪ» force ripples. The majority of the force ripples are eliminated using a modified 3-phase coil configuration as shown in Figure 2.17 with currents through additional conductors commutated as š¼š = āš¼š„ššāš§ššš ššš (š āš3) ā š¼š§ššāš§ššš š š š (š āš3) (2.64) š¼š = āš¼š„ššāš§ššš ššš (š ā2š3) ā š¼š§ššāš§ššš š š š (š ā2š3) (2.65) and š¼šā² = āš¼š, š¼šā² = āš¼š , š¼šā² = āš¼š. This multi-phase coil pattern approximates a sinusoidal current distribution with six steps perĪ». In particular, the 2, 4, 8 and 10 cpĪ» force ripples are canceled out in this configuration. In each magnet spatial period š, six individual coil traces contribute to force generation as shown in Figure 2.17. Because the magnetic field components in x- and z-directions are electrically orthogonal (90 degrees apart in electrical angle), this commutation law allows for the fully independent actuation of a translation (x-directed) and levitation (z-directed) force. The mean position independent force generated is therefore šš = š¾0 [š¼š„š0š¼š§š] (2.66) where the motor force constant is š¾0 = 3ššššš¾1. Thus the vector force in (2.66) is the desired force used for 6-DOF motion actuation. However, due to the existence of the 5th harmonic field from the permanent magnet array, an undesirable force ripple of 6 cpĪ» is also produced. By applying the commutation method (2.63)-(2.65), this force ripple can be calculated as šš = š¾6š4š§ššš [āš¼š„š ššš (6š) ā š¼š§š š šš(6š)0āš¼š„š š šš(6š) + š¼š§š ššš (6š)]] (2.67) 70 where š¾6 = 3ššššš¾5. Other force ripple vectors due to field harmonics š©š can be determined similarly. Figure 2.17 Modified 3-phase coil array below finite width Halbach array. The practical coil stator array is a set of coil traces stacked in š§š. Figure 2.18 shows the general case with the trace spacing in the z-direction š”š. As the āgap factorā (the field decay term from š©š), is of the form šš§š/šš1, the total mean force for š coil arrays spaced š”š apart is simply šššµ = š¾01 ā šāš(š”š+š”š)/šš1 ā šā(š”š+š”š)/šš [š¼š„š0š¼š§š] (2.68) Similarly the total 6 cpš force ripple is šššµ = š¾6š4š§ššš1 ā šā4š(š”š+š”š)/šš1 ā šā4(š”š+š”š)/šš [āš¼š„š ššš (6š) ā š¼š§š š šš(6š)0āš¼š„š š šš(6š) + š¼š§š ššš (6š)]] (2.69) We note that given a continuously varying current sheet, rather than the discretely distributed 3-phase coil arrangement analyzed above, the magnetic field produced by the current distribution 71 would be perfectly locked to the magnet array position and therefore the magnet array would see no variation in magnetic pressure with motion in š„š. A continuously varying current distribution is impractical to manufacture, but increasing the number of phases will decrease the position dependent force component at the expense of current driver complexity. The net torque depends on the net force combined with the location of the center of magnetic pressure. This center depends partly on the discrete nature of the coil distribution, which causes a spatial variation of the center of pressure dependent on the coil pitch (for the 3-phase arrangement, š/6). As a first approximation, the center of magnetic pressure can be assumed to be stationary in the magnet frame and close to the center of the array at the location of the middle of the magnet array and in the center of the current distribution in in š§š for thin š”š; i.e. š„š = 0, š§š = ā(š§š +š”š2) for a single layer of coils. Therefore the mean torque is š»š = š¾š(š»š+š”š2+ š§š)š¼š„ššĢ. In general, the force and torque model may differ from the real motor in several ways. First, the geometric variation of the motor from the model will introduce modeling error. Second, and more significantly depending on the scale of the motor, variations in magnetization strength and magnetization direction (up to 10% and +/- 10 degrees for some typical commercial magnets) will contribute to overall modeling error. One advantage of this motor is that since there is no ferritic or highly permeable material in the system, variations between motors can be limited by careful attention to geometric tolerancing and uniform magnet selection alone. 72 Figure 2.18 Multiple stationary coil array stack-up dimensions, N-layers of identically driven coil arrays. 2.3 Force and Torque Ripple Self-attenuation Using Magnet Array Design The force ripple (2.67) at low flying heights can become a performance limiting factor, especially for high precision applications. As shown in (2.53) and (2.54), higher order spatial field harmonics decay with š§š at a higher rate compared to the fundamental field harmonic, so increased flying height reduces force ripple. Flying higher, however, also decreases the achievable acceleration of the mover since the fundamental field harmonic which produces the mean force is reduced, and there is a maximum current achievable by the current driver as well as a maximum current carrying capacity of the coils. Lee et al. [82] used a higher-order Halbach array with 8 magnetization segments per spatial period to reduce higher harmonics of field. Other strategies to minimize linear motor force ripples include non-rectangular-shaped magnet elements for lower force ripples as in Zhang et al. [83] and Okamoto et al. [84]. Here we ignore such software based controller methods for addressing intrinsic motor disturbances since we assume that it is always better to develop a more linear motor than to attempt correction using software/control methods. Such methods can always be applied in addition if computation cost is not a factor. 73 In this section we present two magnet array designs and a design methodology to reduce force and torque ripple in the 2D motor element. We term these designs self-attenuating because they rely only on magnet array configuration and do not change either commutation laws or controller laws. A whole magnet array is split into sections and each section is offset relative to the commutation center of the excitation current zone. This method of splitting magnet arrays to reduce intrinsic motor disturbances can be generalized to other synchronous machines, in particular other planar levitation synchronous motors, ironless linear motors where the primary source of force ripple is due to field and current interactions, and can also be extended to rotary synchronous machines. 2.3.1 Force Vector Rotation from Offset Commutation Law Since different spatial force harmonics are caused by modulation of different field harmonics, it is useful to examine the effect on each spatial force ripple vector of offsetting the magnet center (š„š, š§š) in the š„š direction without compensating for this offset in the electrical angle š used in the commutated currents. Figure 2.19 shows the true magnetic center of the array at (š„š, š§š) and the offset magnet coordinate system (š„šā², š§šā²) used as the virtual magnet center in the commutation law of (2.63)-(2.65), offset by š„š from the true magnet center at (š„š, š§š). As š„šā² = š„š + š„š, the magnetic field PHM from (2.57) can be re-written as š©š·šÆš“ā² =[ āšµ1šš§ššš š šš(šā² āšš) + šµ5š5š§ššš š šš (5(šā² āšš)) 0šµ1šš§ššš ššš (šā² āšš) āšµ5š5š§ššš ššš (5(šā² āšš)) ] . (2.70) Here we define the virtual electrical angle šā² =š„šā²šš1 and offset electrical angle as šš =š„ššš1. 74 Figure 2.19 Offset magnet array from commutation center. The commutation law (2.63)-(2.65) is modified to be synchronous with the virtual electrical angle šā² and thus š¼š = āš¼š„š šāš§ššš ššš šā² ā š¼š§ššāš§ššš š šš šā² (2.71) š¼š = āš¼š„ššāš§ššš ššš (šā² āš3) ā š¼š§ššāš§ššš š š š (šā² āš3) (2.72) š¼š = āš¼š„ššāš§ššš ššš (šā² ā2š3) ā š¼š§ššāš§ššš š š š (šā² ā2š3) (2.73) Carrying out the same force analysis from 2.2.6 using the shifted field (2.70), the position independent force is ššā² = š¾0 [š¼š„š cos šš + š¼š§šsin šš0š¼š„š š šš šš + š¼š§šššš šš] (2.74) Similarly, the 6-cpš force ripple vector becomes in the true magnet frame (š„š, š§š) ššā² = š¾6š4š§ššš [āš¼š„š ššš (6š ā 5šš) ā š¼š§š š šš(6š ā 5šš)0āš¼š„š š šš(6š ā 5šš) + š¼š§š ššš (6š ā 5šš)]] (2.75) 75 It can be seen from (2.74) that the position independent mean force vector is rotated through the angle šš, while (2.75) shows that the 6-cpš force ripple vector is rotated by 5 times that same angle. From a similar analysis for a general force ripple vector due to the šš”ā field harmonic, the rotation of each of the force ripple vectors can be determined whereby the šš”ā field harmonic produces a (š Ā± 1) harmonic force ripple, rotated through the angle ššš, for š > 1 as shown in Figure 2.20. 76 Figure 2.20 Force vector rotation due to commutation center offset. (a) Position independent force vector šš. (b) Force ripple vector šš. (c) Force ripple vector ššĀ±š resulting from field harmonic š©š. 77 2.3.2 Split Magnet Array Design for Self-attenuation of Force Ripple As shown in Figure 2.21, the magnet array of Figure 2.17 is now split into two sub-arrays in the Y-direction: the array half š£ is offset in the +š„š direction by š/20, and the array half š¢ is offset in the āš„š direction by š/20. The total spacing between array halves is thus š = š/10, which is half the spatial wavelength of the 5th field harmonic responsible for the dominant and undesired 6 cpš force ripple. The coil current is commutated based on the relative position of each coil trace with respect to the center of the whole magnet array according to (2.63)-(2.65). The magnet array width šš used to calculate force coefficients is the total width covered by the magnet volume in the š„š-direction (i.e. for the array shown in Figure 2.21 the magnet array width used to calculate š¾0, š¾6 is šš = 2š, not the total physical extent of the array 2š + š). Figure 2.21 Split array solution for force ripple attenuation of the 6th harmonic of force ripple due to the 5th field harmonic. With the split in the magnet array, the 6 cpš force ripple vector in (2.75) becomes two parts, one acting on each array half. The 6 cpš force ripple on array half š¢ is 78 šš,š =š¾62š4š§ššš[ āš¼š„š ššš (6š +š2) ā š¼š§š š šš (6š +š2)0āš¼š„š š šš (6š +š2) + š¼š§š ššš (6š +š2)]] (2.76) and the 6-cpš force ripple array half š£ is šš,š =š¾62š4š§ššš[ āš¼š„š ššš (6š āš2) ā š¼š§š š šš (6š āš2)0āš¼š„š š šš (6š āš2) + š¼š§š ššš (6š āš2)]] (2.77) The 6-cpš force ripples generated on two array halves have a full 180ā angle and equal magnitude. The resulting sum of (2.76) and (2.77) is šš = šš,š + šš,š = 0, independent of lateral position š„š. The position independent force for the split array design is šš = š¾0 cos (š10) [š¼š„š0š¼š§š] = 0.95š¾0 [š¼š„š0š¼š§š] (2.78) Thus the 6 cpš force ripple is eliminated. In comparison with (2.68), the position independent force from the split array design decreases by 5%. The splitting gap š = š/10 between array halves also increases the XY footprint area of the whole mover by 3.3%. This magnet array splitting method can be generalized to target different force ripple harmonics by spacing the array halves apart by half the spatial wavelength of the field harmonic responsible for the targeted force ripple. Table 2.1 shows the designed gap š for a split array targeting different force ripples present in the 4-segment Halbach array with a 3-phase coil arrangement. In reality, for 4-segment Halbach pattern the 5th harmonic of the permanent magnet array field is much larger than the others, and thus is the main design consideration for the split array. 79 SPACING FOR SPLIT ARRAY TARGETING DIFFERENT FORCE HARMONICS Field Harmonic k Field Harmonic Wavelength Force Ripple Separation Gap š 5 š/5 6 ššš š/10 13 š/13 12 ššš š/26 17 š/17 18 ššš š/34 33 š/33 32 ššš š/66 Field harmonics for a 4-segment Halbach array, force harmonics from actuation by a 3-phase coil array. Table 2.1 Spacing for Split Array Targeting Different Force Harmonics To verify the analytical model and split array force analysis, an FEA was carried out in COMSOL Multiphysics using the AC/DC module in 2D planar mode [79]. The 2D planar mode is used to minimize computation time and to more closely follow the assumptions of the analysis. As before, the leakage effects from the finite extension in š¦š will cause a small reduction in the overall force generated. Figure 2.22 shows the finite element models in schematic form for each case. In order to separate the effects of a finite coil array from the effects of the modulation of higher spatial harmonics of field with the commutated current, an āinfiniteā or extended array was used; the coil array was 18š wide compared to 2š magnet arrays. The coil array width was verified to be wide enough to emulate an infinite repeating current pattern by examining the 1 cpš force ripple: the 1 cpš force ripple amplitude should go to zero as the current sheet goes to an infinite width based on the force equivalence principle. Different coil array widths were attempted until the 1 cpš force ripple amplitude was below the numerical precision of the FEA. 80 Figure 2.22 FEA models for āinfiniteā coil arrays (extended in the x-direction over 18š width). (a) Non-split array. (b) Split array. Figure 2.23 shows the results of a pure levitation (Z) command on the non-split and split arrays. The mean force for the non-split array is 14.57 N compared to the mean force for the split array at 13.85 N for a command of š¼š§š = 10 š“ at a flying height of š§š = š/60 above the top conductor layer. This shows a 5% reduction in the force coefficient as predicted by (2.78). The 6 cpš force ripple amplitude in the non-split array is 0.215 N, and in the split array is 3.65 Ć 10ā6 N. The targeted 6 cpš force ripple is reduced to 0.26 parts per million (ppm) (the numerical 81 precision of the FEA) of the mean force in the split array, from 14756 ppm of the mean force in the non-split array: a reduction of more than four orders of magnitude. The remaining force ripple in the split array is at 12 and 18 cpš, due to the 13th and 17th field harmonics respectively. The amplitude of the remaining force ripple is 260 ppm of the mean force and can be targeted by further coil averaging techniques and combinations of different array splitting as first described in [75]. Figure 2.23 Finite element analysis of non-split array and split array over āinfiniteā coil arrays( extended over 18š width). Flying height šš = š. š [š¦š¦], š°šš = š, š°šš = šš [š], š = šš [š¦š¦], š©š = š. ššš [š], š«š =šš,š¾š = š. ššššš, šš = š. ššš [š¦š¦] and 8 layers of active coils spaced š. ššš [š¦š¦] apart in the āš³š¦ direction. 82 2.3.3 Experimental Load Verification for Force Ripple Self-attenuation In order to demonstrate experimentally the effectiveness of the magnet array splitting method presented in 2.3.2, two magnet arrays were built as shown in Figure 2.24: a non-split array of (2š Ćš»š Ć 2š) dimensions and a split array of ((2š +š10) Ć š»š Ć 2š) dimensions designed to attenuate the 6 cpš force ripple. Figure 2.25 shows the experimental setup designed to measure the load characteristic of two 1-D magnet arrays actuated by a 1-D coil array. The magnet array is attached to an ATI Industrial Automation MINI-45 6DOF load cell that measures three forces and three torques via an adapter plate. The load cell is attached to the spindle housing of a CNC. The CNC spindle is used to translate the magnet array relative to the stationary coils, while a set of linear power amplifiers is used to drive the excitation current through the coils. We did not use the switching amplifier used for the alpha prototype in order to get a cleaner load characterization unsullied by electromagnetic switching noise. The drive capacity of the linear amplifiers is limited to 10 A continuous, so all experimental measurements are done near this limit at 9.5 A to maximize the signal-to-noise ratio of the load cell. 83 Figure 2.24 Experimental arrays for load measurement, showing dimensions and magnetization directions. Figure 2.25 Experimental setup for load measurement. Figure 2.26 shows the two test cases measured, with a flying height š§š = 0.0247š =0.74 mm and š¼š§š = 9.5 A (a pure levitation force is commanded). Since the amplifiers were near their maximum continuous current limit, the flying height was not compensated in the commutation command (i.e. š§š in (2.63)-(2.65) was set to zero). The excitation zone was kept to 84 4š wide (four coil groups) due to experimental limitations of the test coil board, and this was deemed sufficient to meet the requirement to consider the magnetic field purely periodic (based on the principle of equivalent force generation in 2.2.3) and is the practical excitation zone used in the alpha prototype. Figure 2.26 Experimental test cases with finite coil array of 4š width. (a) Non-split array. (b) Split array. Figure 2.27 shows the comparison of the experimental non-split array to the split array. Table 2.2 shows the individual force components of interest, comparing the analytical predictions of each force component using (2.84) to experiment. Due to the variation of actual magnet 85 remanence (Ā±5%), the šµš value in the analytical model is calibrated by matching the analytical levitation force š¹0š§ for the non-split array with the experimentally measured levitating force. All other ratios were then compared without further adjustment of any kind. The key result is that the experimental 6 cpš force ripple for the split array is 10 times smaller than that of the non-split array, from an amplitude of 1.1% of the levitation force for the non-split array to 0.12% of the levitation force for the split array. This indicates that the split array design solution is an effective and practical method of force ripple attenuation in synchronous planar levitation motors. The analytical model predicts that the 6 cpš force ripple should be entirely eliminated. In reality, several practical limitations can cause non-perfect cancelation between the array halves: (a) each magnet element is not necessarily magnetized to an identical magnet remanence; (b) the gap spacing š on the experimental split array could not be manufactured to be exactly š/10; and (c) each magnet segment has a geometric tolerance. There exists a small amount of 1 cpš force ripple, due to the use of 4 š wide coils excited in this experimental setup (as opposed to the infinitely wide repeating coil array assumed by the analytical model) which is predicted by numerical simulation (both FEA and numerical integration using the exact 3D field model from 2.2.1). We note also that the orthogonal force š¹š¦ for both magnet array types is essentially zero, at the level of the noise floor of the sensor. This confirms that the motor has minimal coupling to the orthogonal axis. 86 Figure 2.27 Comparison experimental non-split array and split array, uncompensated for š³š = š. ššššš ,šš±š« = š, šš³š« = š. š š, š = šš š¦š¦,šš¦ = šš,šš = š. ššššš,šš¦ =šš, šš = š. ššš š¦š¦ and N = 8 layers of active coils spaced šš = š. ššš š¦š¦ apart in the āš³š¦ direction. 87 FORCE HARMONIC COEFFICIENTS IN Z F0z % F0z F6z % F0z Non-split Array Analytical 12.33 N 100% 0.153 N 1.24% Experimental 12.33 N 100% 0.137 N 1.11% Split Array Analytical 11.72 N 95% 0 N 0 % Experimental 11.79 N 95.6% 0.014 N 0.11% Br = 1.2 T , adjusted to match experimental and analytical non-split array F0z. All harmonics are compared to the analytical non-split array F0z. Table 2.2 Force Harmonic Coefficients in Z 2.3.4 Quad-Split Array Design for Self-attenuation of Force and Torque Ripple The split magnet array method in 2.3.2 demonstrated the elimination of a targeted spatial harmonic of force ripple to within FEA simulation accuracy; further, experimental characterization in 2.3.3 showed an order of magnitude attenuation of force ripple. However, the issue of torque ripple was not addressed. To meet the positioning and control requirements for high precision stages both torque and force ripple must be attenuated as much as possible. For example, force and torque ripple are discussed in the context of magnetically levitated stages in [45]. In this section, the split array method is extended to a magnet array design that simultaneously self-attenuates targeted spatial harmonics of torque and force ripple without inducing any coupling between any axes or additional distortion at the targeted spatial frequency. We term this new magnet array topology a quad-split array. The quad-split array presented is designed to attenuate both 6 cpš force and torque ripple simultaneously. Finite element analysis (FEA) based simulations are carried out to validate the new magnet array design. 88 Figure 2.28a shows the 6 cpš force ripple vectors (šš,š, šš,š) acting on each half of the non-split magnet array at their respective centers of magnetic pressure for a given š„š position. This force ripple cancels out almost cancels out all net 6 cpš torque on the array due to the modulation of the 5th field harmonic and the current distribution for flying heights, notably š§š =š/60. We note, however, that there still always exists a net 6 cpš torque due to the pitch of the coil array (š/6) and the discrete nature of the coil distribution (approximately a stepwise function in š„š). Figure 2.28b shows the split magnet array with 6 cpš force ripple vectors (šš,š, šš,š) acting on each magnet array half at moment arms (š š, š š). The summation of these force vectors cancels out any net 6 cpš force on the array, as described in 2.3.2. However, an additional torque ripple is induced by the opposition of the 6 cpš force ripple vectors (šš,š, šš,š) as they have a relative angle of 180ā at all positions and are on opposing sides of the magnet array CG. We note that only the z-direction components of force ripple contribute to this torque as the x-directed force components act at the same z-height. If the center of magnetic pressure for each array half is modelled at š š =āšāš2?Ģ? ā (š»š+š”š2+ š§š)?Ģ? and š š =š+š2?Ģ? ā (š»š+š”š2+ š§š)?Ģ? respectively, the additional torque ripple resulting from the 6 cpš force ripple is š»š = āš + š2š¾6š4š§ššš [š¼š„š ššš (6š) + š¼š§š š šš(6š)] šĢ (2.79) where šĢ is a unit vector in the š¦šādirection. This torque ripple will reduce the positioning accuracy of the moving stage and must be addressed for high performance positioning stage applications such as lithography wafer steppers. It is desired to attenuate this 6 cpš torque ripple while simultaneously attenuating force ripple without inducing any additional coupling between axes. 89 Figure 2.29 shows two example embodiments of the quad-split array. Each is comprised of eight identical sub-arrays (šæš¼ , šæš¼š¼ , šæš¼š¼š¼ , šæš¼š, š š¼ , š š¼š¼ , š š¼š¼š¼ , š š¼š). Each sub-array is magnetized in a Halbach-pattern in the XZ-plane and is uniform in the Y-axis. Figure 2.29c shows an example magnetization pattern for the sub-array with horizontally and vertically magnetized elements, with a sub-array width in X and depth in Y of š and magnet height in Z of š»š. In embodiment 1 shown in Figure 2.29a, sub-arrays (š š¼š¼ , š š¼š¼š¼) are shifted by š = š/10 in š„š and sub-arrays (šæš¼š¼ , šæš¼š¼š¼) are shifted by ā š = āš/10 in š„š; the remaining sub-arrays (šæš¼ , šæš¼š , š š¼ , š š¼š) are not shifted. In embodiment 2 shown in Figure 2.29b, sub-arrays (š š¼ , š š¼š) are shifted by š = š/10 in š„š and sub-arrays (šæš¼ , šæš¼š) are shifted by ā š = āš/10 in š„š; the remaining sub-arrays (šæš¼š¼ , šæš¼š¼š¼ , š š¼š¼ , š š¼š¼š¼) are not shifted. As the force and torque analysis for both embodiments are similar, the remaining analysis will be for embodiment 1. Figure 2.28 Magnet array designs showing force ripple vectors acting at magnetic centers of pressure. (a) Non-split magnet array. (b) Split magnet array. 90 Figure 2.29 Quad-split array for force and torque attenuation. (a) Embodiment 1. (b) Embodiment 2. (c) Sub-array example Halbach magnetization pattern. Figure 2.30 shows 6 cpš force ripple vectors for sub-arrays (šæš¼ , šæš¼š¼ , š š¼ , š š¼š¼). From (2.75), the 6 cpš force ripple vector on sub-array šæš¼ is šš,š³š° =š¾68š4š§ššš [āš¼š„š ššš (6š) ā š¼š§š š šš(6š)0āš¼š„š š šš(6š) + š¼š§š ššš (6š)}] (2.80) acting on the moment arm ššæš¼āā āā ā, and the 6 cpš force ripple vector on sub-array šæš¼š¼ is šš,š³š°š° =š¾68š4š§ššš [āš¼š„š ššš (6š ā š) ā š¼š§š š šš(6š ā š)0āš¼š„š š šš(6š ā š) + š¼š§š ššš (6š ā š)}] (2.81) acting on the moment arm ššæš¼š¼āā āā āā ā . The total 6 cpš force ripple is cancelled out between these two sub-arrays. Similarly (š š¼ , š š¼š¼) internally cancels out 6 cpš force ripple. In addition, as shown schematically in Figure 2.29a, sections (šæš¼ , šæš¼š¼) and (š š¼ , š š¼š¼) are mirrored across the š¦šš§š-plane. Thus, the 6 cpš torque ripple induced in the š¦š-axis for section šæš¼ is cancelled out by opposing 91 torque ripple induced by section š š¼. Similarly the 6 cpš torque ripple induced by section šæš¼š¼ is cancelled by the torque induced by section š š¼š¼. The entire array is symmetric with respect to both š„šš§š- and š¦šš§š-planes. This ensures that there is no 6 cpš torque ripple generated on any of the š„š-, š¦š- or š§š-axes on the quad-split array. Figure 2.30 Quad āsplit force ripple vectors acting at equivalent moment arms of sub-arrays (š³š°, š³š°š°, š¹š°, š¹š°š°). Table 2.3 summarizes the mean and 6 cpš ripple forces acting on each sub-array as well as the mean and 6 cpš net forces acting on the quad-split array. Similarly, Table 2.4 summarizes the mean and 6 cpš ripple torques acting on each sub-array as well as the mean and 6 cpš net torques acting on the quad-split array. From this analysis, all net 6 cpš force and torque ripples are shown to cancel out in every axis. The net mean force that remains on the quad-split magnet array is 92 šš = 0.9š¾0 [š¼š„š0š¼š§š] (2.82) From Table 2.4, the remaining mean torque is approximately š»š = 0.9š¾š(š»š + š”š2+ š§š)š¼š„ššĢ (2.83) which can be compensated easily as it is a constant torque across all positions and dependent only on the X-direction force command š¼š„š. If other force and torque ripple harmonics are to be targeted for attenuation, the spacing š should be adjusted as per Table 2.1, where š is determined by the half-wavelength of the targeted field harmonic generating the force and torque ripple. To validate the quad-split array design and examine the remaining torque and force ripple, FEA based simulations are carried out using COMSOL Multiphysics [79] in 2D planar mode using motor parameters similar to the final prototype (described in later sections). Figure 2.31 shows a finite element model schematic, with an 18š-wide active coil zone used to minimize edge effects from the finite width of the magnet arrays as described in 0. Figure 2.32 shows the three magnet array types modelled. Figure 2.33 shows the force and torque results of the three FEA simulations, with a flying height š§š = š/60 = 0.5 mm and š¼š§š = 10 A (a pure levitation force command) over an X-direction traverse of š. The flying height component in the commutation law (2.63)-(2.65) is not compensated (š§š = 0) to mimic a maximum current capacity of 10 A on each drive channel. 93 Table 2.3 Analytical Mean and 6 cpš Ripple Forces for Sub-arrays. 94 Table 2.4 Analytical Mean and 6 cpš Ripple Torques for Sub-arrays. 95 Figure 2.31 Finite element model schematic showing generalized magnet array in relation to stationary coil array with 8 coil layers spaced apart by 0.643 mm in šš. Figure 2.32 Simulated magnet arrays showing magnetization patterns. (a) Non-split array. (b) Split array with š = š/šš. (c) Quad-split array with š = š/šš. 96 Figure 2.33 FEA forces and torques for non-split array, split array, and quad-split array; šš = š/šš =š. šš¦š¦, š°šš = š, š°šš = šš š, with 8 coil layers spaced apart by 0.643 mm in šš. Figure 2.34 shows the harmonic amplitudes of the forces and torques from each FEA. The 6 cpš force and torque ripples are summarized in Table 2.5 along with the mean levitation force values. From Figure 2.34a and Figure 2.34b, both the split array and quad-split array show an almost complete attenuation of 6 cpš force ripple compared to the non-split array (approximately four orders of magnitude). This is on the order of the numerical precision of the finite element analysis. Figure 2.34c shows that the split array induces a torque ripple at 6 cpš which is almost 97 six times larger than the non-split array, while the quad split array has attenuated the 6 cpš torque ripple to 5.1% of the non-split array and 0.88% of the split array. From symmetry, 6 cpš torque ripple on the other axes are zero. Thus the quad-split magnet array is able to attenuate both force and torque ripple at a targeted spatial frequency without additional coupling between axes. This is at the cost of a 10% loss in mean force generation and corresponding increase in mover area In addition, this attenuation of torque is partial validation for the approximation of the magnetic center of pressure for each sub-array, as discussed at the end of section 2.2.6. FORCE AND TORQUE AMPLITUDES F0z (N) % F6z (N) % T6y (Nm) % NON-SPLIT 30.1 100 0.453 100 0.00127 100 SPLIT 28.6 95.0 0.000135 0.003 0.007327 577 QUAD-SPLIT 27.2 90.4 0.0000636 0.0014 0.0000652 5.1 š§š =Ī»60= 0.5mm, š¼š„š = 0, š¼š§š = 10 A, with 8 coil layers spaced apart by 0.643 mm in š§š. Table 2.5 Force and Torque Amplitudes. 98 Figure 2.34 Spatial harmonic amplitudes for force and torque ripples from FEA; flying height šš = š/šš =š. š š¦š¦, š°šš = š, š°šš = šš š, with 8 coil layers spaced apart by 0.643 mm in šš. (a) Translation force, šš. (b) Levitation force, šš. (c) Torque in šš-axis, š»š 99 2.4 6-DOF Motion Stage Prototype Design A 6-DOF alpha prototype (i.e. basic proof of concept prototype) is designed to test the synchronous planar levitation principle. The mover under gravity preload is passively stable in three axes (š§, š¼, š½) and inherently unstable in three axes (š„, š¦, š¾). We therefore require position feedback and control in at least three axes. In order to test the prototype, three additional technologies were developed by colleagues in the UBC Precision Mechatronics Laboratory: a 6-axis position sensing stereo-vision camera was developed by Dr. Xiaodong Lu and Niankun Rao [85]; a 48-channel PWM 3-phase power amplifier to drive the stator currents; and a real time computer Tsunami developed by Kristofer Smeds and presented in [2] to accept the position feedback, execute the digital controller, output current commands to the power amplifier and interface with a host user interface on a non-real time PC. 2.4.1 Prototype Architecture The prototype basic architecture is shown in schematic form in Figure 2.35. The system components are a mover with magnet arrays attached, a coil array (built as a printed circuit board) attached to a stationary base frame, a stationary metrology reference frame attached to the same base frame, and a stereo-vision based position sensor. There are four infrared light-emitting diodes (IR-LEDS) attached to the corners of the mover, and another four IR-LEDS attached to the stationary metrology reference frame with known relative positions. The stereo-vision camera outputs the six-axis relative position (š„, š¦, š§, š¼, š½, š¾) (where (š¼, š½, š¾) are rotations about (š„, š¦, š§) respectively) between the metrology frame and the mover based on the 3-DOF (š„, š¦, š§) position of 100 each LED in the position sensor coordinate frame. This six-axis position is fed back to the real time computer and processed by the digital controller which outputs current commands to the power amplifier. The power amplifier drives current through each coil trace in the coil array, actuating the mover in 6-DOF. Not shown is the host computer which runs the user interface to the real time computer. As shown in the schematic cross section of the alpha prototype in Figure 2.36, a stiffening layer is used to support the coil array PCB. This stiffening layer should be electrically insulating to minimize eddy current losses. The mover chassis should be able to support both the actuating forces applied by the stator current to the magnet arrays as well as the forces between the magnet arrays themselves (these forces are analyzed in [86]). Figure 2.35 Prototype architecture schematic showing 6-DOF camera metrology, power amplifier and real-time controller. 101 Figure 2.36 Schematic cross-section of alpha prototype. Figure 2.37 shows the solid model of the alpha prototype. To obtain the necessary field of view using the 6-DOF position sensor, the stereo-vision cameras must be placed approximately a meter away from the plane of the LEDs (i.e. the top of the mover) in the orientation shown. The camera frame is extruded aluminum. This necessarily means that the camera will be subject to vibration as the camera frame is relatively un-stiff. Hence the stationary metrology reference allows a differential measurement in all six axes, meaning the position feedback is decoupled from camera frame vibration. 102 Figure 2.37 Alpha prototype solid model showing 6-DOF camera metrology. 103 2.4.2 Mover Design 2.4.2.1 Magnet Array Geometric and Material Parameters Many different optimization cost functions can be assembled from the force coefficient expression. Here we choose to maximize acceleration based on a prescribed current density maximum, assuming the mass of the mover chassis is negligible relative to the magnet arrays. From the force expression (2.68), the net force coefficient for each magnet array is š¾ =12šµšcos (š4 +š2šš)š(1 ā šāš»ššš )š·ššššš š ššš (šš2šš) (ššš”š) (1 ā šāš”ššš)1 ā šāš(š”š+š”š)šš1 ā šāš”š+š”ššš [šš¼š] (2.84) where š¼š is the amplitude of the force commands š¼š„š , š¼š§š. These force commands can be interpreted as current commands given a constant flying height. The mass for each 2-DOF motor element is š = ššššššš”š»šš·ššš. Therefore acceleration per force command unit is simply š = š¾/š. Halbach Magnet Array Geometry We choose a four segment per wavelength Halbach array to minimize cost and manufacturing complexity. Force is proportional to (1 ā šāš»ššš ) and acceleration is proportional to 1š»š(1 ā šāš»ššš ). To maximize both we define a cost function šš» =1š»š(1 ā šāš»ššš )2. This cost function is maximum at š»š = š/5, implying a rectangular cross-section since each magnet segment width will be š/4. However, non-square cross-section magnet arrays would require three different types of magnets for a Halbach array mirror-symmetric about the magnet array š¦šš§š plane; therefore we design the magnet arrays as š»š = š/4. This configuration requires only two types of arrays for symmetric Halbach-pattern arrays and is relatively close to š/5 š/4 Ć š/4 Ćš·š (width and height) and š/8 Ć š/4 Ć š·š segments magnetized in the š»š direction. The gap 104 between segments š should be kept as close to zero as possible. Adhesive two part epoxy will be used to bond magnet segments together and this necessitates a 50-100 micron gap. We examine coil parameters as they relate to acceleration and power in 2.4.3. The choice of magnet spatial period is driven by cost and maximum load capacity required, as well as the relative achievable geometric tolerances of each magnet segment. The cost of magnet segments scales more than linearly by volume; hence we have a push to minimize the overall size of the magnet array and the spatial wavelength š. We also want some reasonably large force capacity as well as insensitivity to machining and fabrication tolerances on the order of 50-100 micron. We therefore choose š = 30 šš as a compromise between the objectives of cost, manufacturability and actuator performance. Further, a ~2:1 ratio between š·š:šš yields the largest magnet area per mover (highest magnetic fill factor) and therefore the highest acceleration. Permanent Magnet Material Rare earth magnets offer very high remanence; however, the higher the remanence the lower the overall operating temperature of the magnetic material before permanent losses in magnetization occur (at significantly lower than the Curie temperature). In addition, in the Halbach configuration the operating temperature may be less than an isolated magnet segment as the bucking fields from neighboring magnet segments will exert a de-magnetizing field, meaning permanent magnetization losses occur at an even lower temperature. An N44SH grade magnet material is chosen which yields a remanence šµš = 1.325 [š] (from a magnetic energy product of 44 MGOe and a relative permeability šš of 1.027). 105 2.4.2.2 Mover Magnet Array Configuration The analysis from 2.3 allows us to design the magnet array configuration of the mover based on the specific needs of our application. As the alpha prototype is a basic proof of concept, we prioritize simple manufacturability and minimal cost. Figure 2.38 shows several mover topologies (note the view is from the bottom of the mover, or the strong side of the Halbach arrays). All array force coefficients are normalized to the high force basic large mover shown in Figure 2.38a. a. High Force Basic Large Mover The high force basic large mover is a set of four 2š Ć š»š Ć 4š Halbach arrays. Each magnet array requires only 2 types of magnet segments. The 6 cpš force ripple on the mover is at a maximum. This basic mover gives the largest mean force and acceleration of all of the mover topologies considered here. This topology has the highest acceleration and simplest structure and is therefore suitable for applications that are not high precision but benefit from high acceleration. b. High Force, Low Force Ripple Large Mover This mover topology uses four 2š Ć š»š Ć 4š Halbach arrays with the centers of each pair of (X1, X3) and (Y2, Y4) arrays separated by 4š +š12. For operating conditions where each pair of 106 Figure 2.38 Alternate mover magnet array layouts (identical š = šššš). 107 magnet arrays is driven by similar force commands with near zero tip/tilt (š¼, š½), the additional š12 offset insures that the 6 cpš force ripple cancels out on the CG of the mover without compromising mean force generation. However, this will also induce a large 6 cpš torque ripple on the mover. Intrinsic force disturbance is minimized in this topology while maintaining high acceleration and simplicity of manufacture. The induced torque ripple is still an issue. c. Low Force and Torque Ripple Large Mover Using four split arrays as discussed in 2.3.2, this mover uses four (2š +š10) Ć š»š Ć 4š split Halbach arrays which generate internal 6 cpš torque ripples. Note the magnitude of the torque ripple induced from each magnet array pair is smaller than the induced torque ripple in Figure 2.38b by a factor of 4š+š12š+š10= 3.717. To further cancel this torque ripple, the magnet array pairs (X1, X3) and (Y2, Y4) arrays are separated by 4š +š12 to make the torque ripple on each array out of phase by 180ā under similar force commands for each magnet array pair. This mover topology requires twice the number of non-square magnet segments compared to either Figure 2.38a or b and has a larger planform area. The maximum force on the mover is 95% of that generated by Figure 2.38a or b. This mover topology has no 6 cpš intrinsic force or torque ripple at the cost of 5% maximum acceleration and a marginal increase in manufacturing complexity and planform area. 108 d. High Accuracy Low Internal Stress Large Mover Four quad-split arrays of (2š +š5) Ć š»š Ć 4š are arranged as discussed in 2.3.3 to eliminate 6 cpš torque and force ripples at the magnet array level. The advantages of this mover topology are: i) magnet array force commands can be very different and internal cancellation of 6 cpš force and torque ripple will still occur; ii) because the internal force and torque cancellations on each magnet array occur over much smaller moment arms compared to Figure 2.38b or c, cyclic internal stresses acting on the mover chassis will be significantly reduced. This renders this topology very suitable for high precision applications. The cost of this force and torque linearity and low internal stresses leading to vibrations is a 10% reduction in maximum force generated compared to Figure 2.38a and a much more complex assembly requiring four different magnet segments as shown. e. High Force Basic Small Mover This mover topology is shown to illustrate that multiple sizes of mover using the same magnetic wavelength š can be actuated on the same motor stator. With 25% of the mover area, the total maximum force on the mover is Ā¼ of that generated on the high force basic mover in Figure 2.38a, but with Ā¼ of the magnet weight the acceleration is the same. The high force, low force ripple mover topology in Figure 2.38b is chosen for the alpha prototype as a compromise between maximum acceleration, cost and intrinsic force disturbance. The intrinsic torque disturbance was not initially considered to be an issue until later analysis. 109 2.4.2.3 Mover Mechanical Design Figure 2.39 shows an exploded view of the mover design. The mover chassis is a lightweight aluminum honeycomb structure. The magnet arrays are attached to the bottom of the mover chassis using adhesive epoxy. A non-magnetic central spacer is used as part of the assembly process (described in the manufacturing section 2.5.1) and left in place after assembly. The plastic side plates are to protect the honeycomb core which, while very strong in sheer, is yet very thin aluminum foil and easily damaged by side loads. Figure 2.40 shows an isometric cut view of the mover with internal details. The mover chassis is a commercially available aluminum honeycomb sandwich panel with solid aluminum top and bottom skins of 1.59 mm thickness and a honeycomb core of 1/8ā hex cell density and 0.002ā aluminum foil thickness. The core thickness is 50.8 mm. This was chosen as a compromise between stiffness and weight, with the aluminum skins thicker than optimal (as calculated in [52]) to allow tapping with M2 threads. The IR-LEDs are powered by a lithium-ion battery and driven by a constant current circuit as shown and attached directly to the top of the mover chassis using M2 socket head cap screws [1]; a battery powered solution is chosen to negate the need for power cables to the mover. The central spacer block is aluminum. Magnet array and mover topology dimensions are shown in Figure 2.41, with š = 30 šš. Each magnet array is made up of two subarrays with dimensions 2š Ć š»š Ć 2š, with each sub-array made of magnet segments with a depth of š·š2= 60 šš to minimize cost and risk of damaging the magnet segments due to too large an aspect ratio. The four magnet arrays combined weigh 1.68 kg and the overall weight of the mover is š = 2.3 kg. 110 Figure 2.39 Alpha prototype mover design, exploded view. Figure 2.40 Isometric cut-view of mover . 111 Figure 2.41 Alpha prototype mover design, detailed planform. 2.4.3 Stator Coil Array PCB Board Design and Wiring Configuration The coil array PCB must implement the coil structure presented in 2.1. A single set of 3-phase coil traces elongated in the Y-direction is shown in Figure 2.42 with zoomed inset showing coil geometry. Ideally šš = š/6 for the 3-phase coil arrangement, but a practical insulation gap is required and therefore šš =š6ā šš where šš is the insulating gap dimension. Figure 2.43 shows an example of an eight layer coil array with both X- and Y-elongated coil trace layers alternating 112 in Z, with four layers each. Note that the top layer is X1. From (2.68) we define array force coefficients for each magnet array pair as š¾š =12šµšcos (š4 +š2šš)š(1 ā šāš»ššš )š·ššššš š ššš (šš2šš) (ššš”š) (1 ā šāš”ššš)1 ā šāš(š”š+š”š)šš1 ā šāš”š+š”ššš [šš¼š] (2.85) š¾š = šā(š”š+š”ššš )/ššš¾š [šš¼š] (2.86) where š¾š is the force coefficient for the (X1, X3) and š¾š is the force coefficient for (Y2, Y4) and we can define an insulation thickness in the z-direction as š”ššš . From (2.85), minimal separation thickness š”š provides the highest acceleration. The separation thickness between layers of coil traces acting in the same direction is š”š = 2š”ššš + š”š. For the ideal case of zero insulation thickness, the minimum separation thickness is still š”š = š”š. From this structural limitation we see that simply maximizing the copper thickness š”š is non-optimal as the difference between the force coefficients will be too large (with š¾šš¾š= šā(š”š+š”ššš )/šš). The Ohmic power loss in a single š-coil group is š = š¼š2š + š¼š2š + š¼š2š [š] (2.87) with the resistance š = šš ššš”šššæ, šš the resistivity of the coil material, and šæ the total length of the coil trace in each phase for each layer. Therefore the Ohmic loss is directly proportional to the number of coil layers š. From (2.86) and (2.88) we can define a force/power ratio that is related to the motor electrical efficiency, šš = š¾š„/š. We note that the total effective board thickness is š(2š”š + š”ššš ). We can see that there are diminishing returns for additional coil layers; this is 113 intuitively true because each additional coil trace layer is immersed in a field decaying in š§š as šš§š/šš . Figure 2.42 Coil structure isometric view, single š-group Y-elongated coil traces. 114 Figure 2.43 Eight layer coil structure, showing alternating X- and Y-elongated coil traces. Coil trace end turns are implemented outside the active region of the coil array as shown for a single phase in Figure 2.44 (three layers only are shown). The three phase schematic is shown in Figure 2.45 for a single š-coil group. Note that as per standard three phase motors, only two currents are actively controlled, with the third current automatically generated by the wye- or star-center configuration. Because the commutation law and the physical placement of the coil traces mean that š¼š, š¼š and š¼š have a 60ā electrical angle between them, the standard three phase wiring 115 is modified as shown and the star center insures that š¼š ā š¼š + š¼š = 0, i.e. the three phase currents balance. The current controller actively commands only two currents, š¼š and āš¼š with š¼š automatically generated. Figure 2.44 Single phase multilayer windings, shown for three layers of Y-elongated traces. To minimize the number of drivers required for the alpha prototype, series connections between different š-coil groups are desired. Figure 2.46 shows an example of two š-coil groups connected in series with three-phase wiring implemented. 116 Figure 2.45 Three phase wiring, shown for a single š-coil group. Figure 2.46 Serial connections between two single š-coil groups with 3-phase wiring. 117 In addition to proving the basic synchronous levitation motor concept, it was desired that the coil array PCB be able to demonstrate at least one method of scaling up the planar stroke. Figure 2.47 shows a multiple coil array PCB tile arrangement. Y-elongated coil traces on each board (I, II, III) are attached in series using connector PCBs on the bottom side of the coil array, with individually soldered pins. This allows the use of multiple coil array PCBs without interfering with the working area of the mover and minimal increase in the number of current drivers. We mention this additional functional requirement to explain the difference in layout geometry for the connection points for the X-elongated and Y-elongated coil trace layers, as well as the choice of the X-elongated traces as the first layer (making the Y-stroke have a higher maximum acceleration; the Y-direction is the long stroke axis in the multi-tile configuration). Figure 2.48 shows the coil layer structure. As shown, the via connections for the Y-elongated traces are optimized for minimum Y-width (to minimize the non-uniform coil trace region) and minimal interference with the X-elongated traces. The via connections for the X-elongated traces are outside the working region of the multi-tile coil array stator and are therefore laid out more simply, in straight lines. Figure 2.49 shows the overall coil array PCB layout, with 16 Y-elongated š-coil groups and 9 X-elongated š-coil groups. Groups of the same colour are in series (e.g. šŗš1 is in series with šŗš9). Series groups are chosen based on the stroke of the mover and the need to dynamically shift the active coil zones for each magnet array. The active coil zones for each array are four š -coil groups wide. The Y-elongated coil traces provide actuation in the longer stroke for the alpha prototype, the X-direction. A minimum of eight individually addressable Y-elongated coil groups allows the independent actuation of (Y2, Y4); with the series connection shown it is possible to traverse the entire X-stroke of the stator with only the drivers necessary for these eight coil groups. The X-elongated coil traces provide actuation in the short axis (Y-direction) which is less than the 118 length of the mover and therefore coil groups (šŗš1, šŗš2, šŗš3) actuate only magnet array X1 and coil groups (šŗš7, šŗš8, šŗš9) actuate only magnet array X3 across the whole Y-stroke of the mover. These groups are therefore wired in series to save on the number of drives. Figure 2.47 Extendable coil array PCBs. Inset shows connector PCB with solid pin connectors. 119 Figure 2.48 Coil structure isometric view, X- and Y-elongated coils. (a) Layer X1: x-elongated traces. (b) Layer Y1: y-elongated traces. 120 Figure 2.49 Coil array PCB showing coil groupings, X- and Y-elongated coils. A 16 layer PCB is designed with 8 layers of X-elongated coil traces and 8 layers of Y-elongated coil traces. Copper thickness of 6 oz corresponding to 210 šš chosen based on available 121 heavy copper PCB manufacturing capability and the trade-off between different force coefficients between planar axes and maximum acceleration in the strong planar axis. The minimum insulation thickness available for this weight of copper was š”ššš = 127 šš. The minimum insulation gap for this copper thickness based on a standard copper etching process is šš = 254 šš, making the coil width šš = 4.746 šš. Figure 2.50 shows the final board stack-up, with corresponding separation thickness š”š = 464 šš. Figure 2.50 Coil array PCB stack-up details. 122 Motor parameters are summarized in Table 2.6. From (2.85), (2.86) the force coefficients for each magnet array are š¾š = 3.28 [šš¼š], š¾š = 3.06 [šš¼š]. The maximum force achievable for equal current is 7% lower for the weak stroke direction. We note that the stack-up fill factor for the coil trace array is 63%, and the area fill factor for each coil trace layer is 94%, significantly higher than achievable through typical motor coil windings using round wire. In addition, the PCB feature creation technology automatically allows for alignment tolerances on the order of 50 micron without further effort. ALPHA PROTOTYPE MOTOR PARAMETERS Parameter Symbol Value Units Halbach wavelength š 30 mm Magnet remanence šµš 1.325 T Magnet height š»š 7.5 mm Magnet depth š·š 120 mm Magnet array width šš 60 mm Halbach segment gap š 50 šm Coil trace width šš 4.746 mm Coil trace thickness š”š 0.210 mm Coil separation thickness š”š 0.464 mm Coil layer number in each direction š 8 -- Table 2.6 Alpha Prototype Motor Parameters. 123 2.4.4 Stator Design Figure 2.51 shows the exploded view of the stator assembly. The coil array PCB is supported by a ceramic stiffener made from an electrically insulating quartz-based ceramic. The stator base frame is a lightweight and stiff honeycomb of the same stack-up dimensions as the mover chassis in 2.4.2.3. The coil array PCB is bonded to the ceramic stiffener which in turn is bonded to the top surface of the base frame. The metrology reference frame is a removable part that is installed on the base frame using a kinematic mount. The kinematic mount divorces the metrology frame from internal stresses that may be induced on the stator base frame and coil array by (for example) Ohmic heating. The mounting blocks for the kinematic mount as well as the electrical terminal blocks used to strain relieve the coil array wiring (not shown) are attached to the top surface of the base frame using M2 and M3 socket head cap screws. Figure 2.52 shows detailed dimensions of the stator components. The thickness of the ceramic stiffener insures that the field strength at the top surface of the aluminum skinned base frame is 0.6% of the field strength at the top surface of the coil array, providing low eddy current losses without unduly increasing the overall stator height. 124 Figure 2.51 Alpha prototype stator exploded view. 125 Figure 2.52 Alpha prototype stator isometric cut-view. Figure 2.53 shows the details of the metrology reference frame. The monolithic block is solid aluminum. Three 12mm diameter silicon nitride ceramic balls are bonded to cones machined into the bottom surface of the monolithic block. These ceramic balls mate with commercially available v-groove blocks (VB-375-SM from Bal-tecā¢ [87]) mounted to the base frame to form a minimally constrained kinematic mount. The weight of the monolithic block is sufficient to provide a preload for the mount. Kinematic constraints allow for repeatable and stress free installation of the metrology reference frame. The height of the reference frame is designed so that 126 the reference blockās IR-LEDs are in the same focal plane as the IR-LEDs on the mover for a flying height of 1 mm. This minimizes position feedback error due to lens distortion. The IR-LEDs and constant current circuit PCB are identical to those used on the mover, with the exception of the lithium-ion battery which is unnecessary since the metrology frame is stationary and can be powered by a benchtop power supply. Figure 2.53 Isometric cut-view of metrology reference frame. 127 2.5 6-DOF Prototype Manufacturing and Assembly 2.5.1 Mover Sub-arrays are built using the assembly jig shown in Figure 2.54. To partially stabilize the magnet segments in the Halbach array, a non-magnetic spacer/datum plate is placed on top of a magnetic steel base, with datum plane blocks attached to the spacer plate. The attraction force between the magnet segments and the magnetic base help stabilize the magnet segments as they are brought close together. Three planes constrain the sub-array in XYZ, and lateral and vertical clamping screws are used to push the magnet segments against the datum planes. Loctite E-120HP two part epoxy [88] was used to bond each magnet segment together, and plastic spacing sheets (not shown) were used to isolate the epoxy from the assembly jig. Figure 2.55 shows the actual assembly jig with magnet sub-array installed. 128 Figure 2.54 Magnet sub-array assembly jig. 129 Figure 2.55 Magnet sub-array assembly jig with magnet array installed. To assemble the entire mover, a magnet array assembly jig was designed (Figure 2.56). A non-magnetic stainless steel spacer/datum plate is used as the datum surface on which the mover is built and is sandwiched between the sub-arrays and the magnetic steel base plate. The magnetic steel base plate helps stabilize the sub-arrays when they are brought close together in the XY plane, and provides more than gravity preload to push the magnet sub-arrays onto the spacer/datum plate. The magnet arrays are placed strong side facing this plate. Screw clamp devices are used to laterally constrain the magnet sub-arrays while the epoxy cures. The same epoxy [88] is used. The magnetic steel base plate also has M12 tapped holes through which screws can be threaded and used to separate the non-magnetic spacer/datum plate from the magnetic base plate after the epoxy cures. 130 Figure 2.56 Magnet array assembly jig. Figure 2.57 shows the mover assembly process. In Figure 2.57a, the sub-arrays are brought close together in a stable fashion, and the lateral screw clamps begin to push them together. In Figure 2.57b the lateral screw clamps are fully tightened and the magnet sub-arrays locations are 131 checked. The mover chassis is epoxied to the weak side of the magnet arrays in Figure 2.57c. Once the epoxy has cured, the magnet arrays and the mover chassis are an integral unit; however, now the force between the mover and the magnetic steel base plate is very large. Separation screws threaded into the base plate are used to force the spacer/datum plate apart from the base plate against the large normal force between mover and base plate until the mover can be detached from the spacer/datum plate (Figure 2.57d). After the mover is removed from the magnet array assembly jig, the IR-LED and battery circuit PBCs are test fit (Figure 2.57e). Figure 2.57f shows the bottom (strong side) of the mover with flux sensitive paper showing |šµš§|. 132 Figure 2.57 Magnet assembly. 133 2.5.2 Stator The coil array PCB was manufactured by a commercial vendor using a standard heavy copper PCB process. Figure 2.58 shows the coil array PCB as it arrived from the vendor. Figure 2.58b shows the end turn wiring partially installed using 12AWG solid core wire. Using solder paste and conventional solder, care was taken to insure that solder flowed all the way through the via connections to create minimal electrical resistance between the solid core wire connection and the coil traces. Figure 2.58c shows the ceramic stiffener epoxied to the coil array PCB after the end turn wiring is installed using [88]. Figure 2.58d shows the now integral PCB/stiffener unit being epoxied to the honeycomb base frame. The motor is to be installed on the cement floor, so a rubber damping layer is added to the bottom of the base frame to both provide grip and some vibration isolation (Figure 2.58e). The complete stator with the kinematic mount v-groove blocks installed using M2 screws and three-phase wiring complete is shown in Figure 2.58f. The coil array structure is visible from the edge of the PCB (Figure 2.58f). 134 Figure 2.58 Stator assembly steps. 135 2.5.3 Metrology Reference Frame The metrology reference frame monolithic block is machined from a single block of 7071-T6 aluminum. Figure 2.59a shows the top side of the monolithic block with features for attachment of the IR-LED PCBs and the constant current circuit PCB. Figure 2.59b shows the bottom side of the monolithic block with cone features to accept the kinematic mount silicon nitride balls. Figure 2.59c shows the bottom side of the monolithic block after the silicon nitride balls have been bonded to the block using [88]. Figure 2.59d shows the complete metrology reference frame with IR-LEDs and constant current circuit PCB installed. Figure 2.59 Metrology reference frame assembly. 136 2.5.4 Complete System The complete motor system is shown in Figure 2.60, with metrology reference frame installed. The position sensor frame is shown in the background. Not pictured are the 48-channel switching amplifier and real time computer. The maximum force and acceleration parameters for a flying height of 0.5 mm and maximum current density of 20 A/mm2 is shown in Table 2.7. Due to amplifier limitations, these values were not experimentally tested beyond 1g accelerations. Figure 2.60 Synchronous planar levitation motor system. 137 PREDICTED PERFORMANCE OF ALPHA PROTOTYPE MOTOR Parameter Symbol Value Units Mover mass š 2.3 kg Maximum current density š½š 20 A/mm2 Maximum levitation height š§šššš„ 11.6 mm Maximum levitation force š¹š§ššš„ 228 N Maximum x-translation force š¹š„ššš 118 N Maximum y-translation force š¹š¦ššš 110 N Maximum levitation acceleration šš§ššš„ 9.9 g Maximum x-translation acceleration šš„ššš 5.2 g Maximum y-translation acceleration šš¦ššš 4.9 g Table 2.7 Predicted Alpha Prototype Motor Performance at flying height 0.5 mm. 138 2.6 Motion Tracking Results with 6-DOF Prototype 2.6.1 Plant Modeling and Control Architecture The eight actuating forces shown in Figure 2.1 are related to the mover translation and Euler angles via the following dynamic equation: [ š?Ģ?š?Ģ?š?Ģ?š½š„?Ģ?š½š¦?Ģ?š½š§?Ģ? ] =[ 0 0 1 0 0 0 1 01 0 0 0 1 0 0 00 1 0 1 0 1 0 1ā(š§) š 0 š ā(š§) āš 0 āš0 āš āā(š§) š 0 š āā(š§) āšš 0 āš 0 āš 0 š 0 ] [ š¹š¦1š¹š§1š¹š„2š¹š§2š¹š¦3š¹š§3š¹š„4š¹š§4] (2.88) where š is the mover mass, (š½š„, š½š¦, š½š§) are the moments of inertia around three principal axes, (š, š) are the X-distance and Y-distance between the magnet array Y2 centre and the mover center of gravity (CG) as defined in Figure 2.41, and ā(š§) is the Z-distance between the mover CG and the stator coil. The ā(š§) term represents the dynamic coupling between X-Y translation and pitch/roll motion. For example, acceleration in Y will also induce the mover to rotate around the X axis, because the motor forces š¹š1 and š¹š3 act laterally below the mover, instead of ideally through the mover CG. In addition, this coupling coefficient varies with the mover Z position. To eliminate motion coupling and control the mover in 6-DOF, a multi-input-multi-output (MIMO) motion controller is implemented in the digital controller (Figure 2.61). The MIMO controller output (š¼š„š, š¼š§š) for each magnet array is demodulated into 3-phase currents (š¼š, š¼š, š¼š) for each array, and these current commands are dynamically mapped to the stator X and Y coils in the coil switching block based on the mover position. The coordinate transform block converts the position sensor output into the mover CG translation and rotation coordinates. 139 The MIMO controller allows the implementation of 6 independent SISO control loops. Due to sensor feedback limitations which are beyond the scope of this work, a loop transmission cross-over frequency of 50 Hz was implemented for each axis, with approximately 100Hz closed loop -3dB bandwidth. The limitation is purely due to the feedback sensor and not the drive electronics or the levitation motor itself. The approximate minimum dynamic stiffness within the closed loop bandwidth is calculated as ššš¦ššššš,ššš = (2šššš¤)2š, with the closed loop bandwidth ššš¤ = 100 š»š§ and the mass of the mover š = 2.3 šš, giving ššš¦ššššš,ššš =0.9 š/šš. Within the bandwidth of the system the dynamic stiffness will be higher, with infinite DC stiffness due to integrator action within the load capacity of the motor. Figure 2.61 Controller block diagram. 2.6.2 6-DOF Tracking Results The sensor has a position update rate of 8 kHz and a root-mean-square (RMS) error in the three translation axes of 2.0 šš, 2.1 šš, 5.1 šš in š„, š¦, š§ respectively, and regulation RMS error in the three rotational axes of 0.0014ā, 0.0013ā, 0.0004ā in š¼, š½, š¾ as discussed in [1]. For a stationary 140 position command in all 6 axes, the mover has position regulation RMS errors equal to this sensor noise floor. Motion tracking results for a reference elliptical path are shown in Figure 2.62; the reference path is for 260 mm X-stroke and 60 mm Y-stroke at a flying height of 1mm above the coil array PCB surface. Figure 2.63 shows the 6-DOF tracking error for this trajectory. Tracking RMS error in the three translation axes is 7.49 šš, 5.02 šš, 9.20 šš in š„, š¦, š§ respectively, and RMS error in the three rotational axes of 0.0066ā, 0.0082ā, 0.0029ā in š¼, š½, š¾ respectively. The alpha prototype succeeds in its purpose as a proof of concept; therefore we leave the minimization of this tracking error using this position sensor to future work. In theory, the achievable precision of the synchronous levitation motor is largely limited by the sensor noise floor. If nanometer-level precision is desired, laser interferometers with a target mirror block can be used as in a lithography machine. The achievable translation stroke for the alpha prototype is 260 šš Ć 60 šš Ć 2.5 šš (Figure 2.64 shows the mover floating at š§š = 2.5 šš). Small rotational strokes are achievable, on the order of a Ā±5 degrees in š¾, with rotation in tip/tilt only limited by the flying height. 141 Figure 2.62 Elliptical path of mover in XY, at šš = š šš. 142 Figure 2.63 Tracking error for elliptical path, šš = š šš. Figure 2.64 Mover floating with šš = š. š šš. 143 Chapter 3: Asynchronous Planar Levitation Stage We present in this chapter an analysis of levitation and drag torque of an asynchronous planar levitation motor element using rotating permanent magnet arrays. The conceptual design of a levitation and propulsion motor based on rotating permanent magnets was first presented in [7] with a stable 6-DOF planar levitation stage described in [3] and utilized in [9] [8]. The advantage of asynchronous levitation methods over synchronous magnetic levitation are: i) very simple and passive stator structure requiring at a minimum only an electrically conducting medium; ii) passively stable in all axes under gravity preload, therefore requiring no feedback for stable levitation and actuation. This work focusses on modeling and analysis of the field distribution and generation of levitation and drag forces for a linear 2D motor geometry of infinite horizontal extent, which is extended to an analytical expression for force and torque for rotating 3D magnet disk geometries. The novel analytical expression for force and torque can yield design insight into the motor topology without the need for sensitivity studies using 3D FEA or experimentation. This chapter consists of: ļ· 6-DOF asynchronous levitation motor basic concept and working principle. ļ· Modeling and novel analysis for a limiting 2D case. ļ· Extension of the analytical expression for force, drag and power from 2D geometry to 3D geometry using superposition of fields of cylindrical shell elements. ļ· Experimental load test of single magnet disk over a conducting stator and comparison to the analytical shells model. 144 3.1 Asynchronous Planar Levitation Stage Concept and Working Principle The 6-DOF asynchronous planar levitation stage is comprised of a mover with multiple magnet disks with independently controllable rotation speeds and tilt/tip angles attached to a single mover chassis, levitating over a stationary conducting slab. Figure 3.1 shows a particular embodiment of the planar stage with four rotating magnet disks attached to a mover chassis (note that details of the tilt/tip mechanism are not shown, nor are the onboard power source or controllers for the rotary motors). Each magnet disk is a circumferential Halbach array with the strong field side on the āz face of the disk, and is driven by a rotary motor. In typical operation, pairs of magnet disks are driven in counter-rotation. That is, magnet disk 1 is typically driven in an opposing rotation to magnet disk 3; and similarly magnet disk 2 is driven in opposition to the rotation of magnet disk 4. Each rotating magnet disk generates an asynchronous induction field in the conducting body which acts to repel and counteract the motion of the disk; thus when the disk plane is parallel with the conducting slab, each magnet disk is capable of generating a z-directed levitating force and a drag torque in the axis of rotation. When the disk plane is not parallel with the face of the conducting slab, each disk will produce a z-directed levitation force, a drag torque around the axis of rotation, and a propulsion force in the XY plane. Figure 3.2a shows the principle of propulsion and levitation in the zero tilt/tip case. Magnet disk 1 rotates at š1 in the +z-direction and generates a drag torque šš§1 in the āz-direction, and a levitation force š¹š§1. Similarly magnet disk 3 rotates at š3 in the āz-direction and magnet disk 3 generates a drag torque šš§3 in the +z-direction and a levitation force š¹š§3. When š1 = š3, the drag torques on each disk have equal magnitudes and the torque contribution on the mover due to the pair of magnet disks is šš§1 + šš§3 = 0, and the levitation force contribution is š¹š§1 + š¹š§3 and is always in the +z-direction (because the induced field in the conductor always acts to repel the 145 permanent magnet field from the magnet disk). For š1 < š3, |šš§1| < |šš§3| and a +z-directed torque is generated on the mover as well as a +y-directed torque. Similarly, for š1 > š3, |šš§1| >|šš§3| and a net āz-directed torque is produced on the mover and a -y-directed torque. In each driving case, the levitation forces on each magnet disk are in the +z-direction and will always produce a levitation force contribution in the +z-direction. Figure 3.2b shows the principle of propulsion and levitation in the non-zero tilt/tip case. Magnet disk 1 is rotated by š½ in the -y-direction while spinning about its own axis at š1 in the +z-direction, producing drag torque drag torque šš§1, levitation force š¹š§1 and translation force š¹š¦1. The translation force is produced because the left-hand side of the disk is at a lower flying height over the conducing slab, producing a larger induction field acting on that half of the disk; thus a differential in-plane force is produced. Similarly, magnet disk 3 is rotated by š½ in the +y-direction while spinning about its own axis at š3 in the -z-direction, producing drag torque drag torque šš§3, levitation force š¹š§3 and translation force š¹š¦3. For positive values of š½, a +y-directed force is generated on the mover. Thus by controlling the tilt angle š½ in the y-axis and the rotational velocities š1, š3, the pair of magnet disks 1,3 is capable of generating controllable levitation force, translation force in y-direction, and torques in the y-axis and z-axis. Similarly by controlling the tip angle š¼ in the x-axis and the rotational velocities š2, š4, the magnet disks 2,4 are capable of generating a levitation force, a translation force in the x-direction, and torques in the x-axis and z-axis. Combining all four magnet disks allows 6-DOF actuation of the mover. We note that the axes are coupled. Under gravity preload, the mover is passively stable in all axes. 146 Figure 3.1 Asynchronous levitation planar motor concept using four levitation disks, magnetization shown in red for magnet disk 1 (central actuated pivot mechanism and onboard battery pack/motor controllers not shown). 147 Figure 3.2 Levitation and propulsion principle (details of actuated pivot mechanism not shown). 148 3.2 Asynchronous Levitation Machine Force and Torque Modeling We desire to develop a purely physics based analytical model that can be used as a design tool to develop optimized asynchronous levitation planar machines, without having to rely on experimentally derived correction factors. Here we focus on levitation and drag torque of the magnet disk over a widely extended stationary conducting slab, and leave propulsion modeling to future work. A 2D asynchronous levitation machine comprised of an infinitely extended linear Halbach array and an infinitely extended conducting slab with finite thickness is analyzed. Analytical expressions for the levitation force and drag force are derived that are functions only of material parameters and magnet array and stator geometry. The 2D analytical model is then extended to the case of a rotating 3D circumferential Halbach array over a conducting slab of finite thickness by first splitting the problem into cylindrical shell elements with infinitesimally thin radial thickness, solving for the fields generated by each shell element using the 2D model assuming no radial flux leakage, and superimposing the generated fields from each cylindrical shell element; we term this the analytical shells method. This analytical shells method provides a useful design tool that gives levitation force and drag torque for any geometry of magnet array and stator conductor. 3.2.1 Infinitely Wide 2D Moving Magnet Array over Homogeneous Conducting Slab Linear induction motors have been well analyzed in the literature [65], [55]. The asynchronous machine shown in Figure 3.3 has a permanent magnet array with periodic magnetization distribution š“, relative permeability šš = 1 and thickness š»š over a homogeneous conducting slab of conductivity š, relative permeability šš = 1 and thickness š·š. Both magnet array and conducting slab are of infinite horizontal extent, and depth in š¦š is considered sufficient such that 149 the problem is 2D. We note that only the relative motion between the permanent magnet source and the conducting body is relevant to induced forces, therefore we are free to choose to make either one or the other body stationary for ease of analysis (this relative velocity is the slip speed as discussed in the literature). Here we define the relative motion between the magnet array and the conducting slab using the magnet frame of reference š„š, š¦š, š§š with a stationary magnet array and the conductor translating in the š„š-direction at š£š„ [šš ]. We are interested in the field solutions after steady state conditions have been established and therefore MQS assumptions can be applied. Five regions are defined with boundaries as shown in Figure 3.3: region I is the air above the magnet array (š§š ā (š»š, ā)), region II is the permanent magnet region (š§š ā [0, š»š]), region III is the air gap between the magnet array and the conducting slab (š§š ā (āš§š ā š·š , āš§š)), region IV is the conducting slab, and region V is the air below the conducting slab (š§š ā (āā,āš§š āš·š]). As in 2.2.4, we define a magnetic vector potential šØ such that the magnetic flux density everywhere is š© = ā Ć šØ (2.24). The vector potential expression in the absence of surface currents is once again (2.27). The magnetization distribution is defined as in 2.2.4. Since the problem is 2D, the vector potential has only a single component in the š¦š-direction so šØ = š“š¦š,Ģ everywhere. From Gaussā law (2.4), the definition of the flux density (2.24) and Ampereās law (2.2) and the absence of currents in the air regions (I, III, V), we obtain the Laplace equations for regions I, III and V as ā2šØ = 0 (3.1) 150 For region II, assuming the permanent magnetic material is electrically insulating (relative to the conducting slab this is a good assumption), from (2.26) we obtain Figure 3.3 Stationary Halbach array of infinite extent over a translating conducting slab. 151 ā2šØ = āšš(ā Ć š“) (3.2) just as in 2.2.4. The conducting slab region will have induced currents because every part of the conducting slab will experience a changing magnetic field as the slab moves relative to the permanent magnet array. Re-arranging (2.6) and plugging in (2.2), we get the electric field within the conducting slab purely as a function of the total field distribution and the velocity of the conducting slab: š¬ =1ššš(ā Ć š©) ā š Ć š© (3.3) Substituting (3.3) into Faradayās induction law (2.1), we get ā Ć {1ššš(ā Ć š©) ā š Ć š©} = āšš©šš” and rearranging, 1šššā Ć ā Ć š© = āšš©šš”+ š Ć š© Applying the vector identity (ā Ć (ā Ć š©) = ā(ā ā š©) ā ā2š©) and Gaussā law (2.4), we obtain ā1šššā2š© = āšš©šš”+ š Ć š© (3.4) We have assumed steady state conditions and a stationary magnet array, with only the conductor moving with some velocity š = š£š„ ?Ģ?. Therefore šš©šš”= 0, and (3.4) reduces to the vector Poisson equation ā1šššā2š© = š Ć š© (3.5) From (3.5) two scalar Poissonā equations arise: 152 š2šµš„šš„š2 +š2šµš„šš§š2 = āšššš£š„ššµš§šš§š (3.6) and š2šµš§šš„š2 +š2šµš§šš§š2 = šššš£š„ššµš§šš„š (3.7) The magnetization distribution for a general four segment Halbach array is periodic in š„š and uniform within the magnet volume of each segment as depicted in Figure 2.14. We use the fundamental component of a four segment Halbach array as defined in (2.32) for the following analysis. This analysis can be identically repeated for any spatial harmonic of magnetization by substituting the magnetization fundamental component for any other magnetization harmonic component. The full field solution will be the superposition of all solutions for each harmonic. As will be shown in later sections, the fundamental component of magnetization is responsible for the vast majority of the force generation on asynchronous levitation machines, so we focus on the solution using just this Fourier component. Within the conductor volume, we can guess the form of the magnetic field density as š©š¼š = [šµš„š¼ššµš¦š¼ššµš§š¼š] = š š [{šš¼ššš¼šš§š + šš¼ššāš¼šš§š}šāšš„š/šš0{šš¼ššš¼šš§š + šš¼ššāš¼šš§š}šāšš„š/šš] (3.8) The field everywhere must obey (2.4) and therefore we can immediately find that since ā ā š©š¼š =0, āššµš„š¼ššš„š=ššµš§š¼ššš§š and therefore šš¼š =š¼šššššš¼š (3.9) and 153 šš¼š = āš¼šššššš¼š (3.10) So (3.9) becomes š©š¼š = š š[ {š¼šššššš¼ššš¼šš§š āš¼šššššš¼ššāš¼šš§š}šāšš„š/šš0{šš¼ššš¼šš§š + šš¼ššāš¼šš§š}šāšš„š/šš ] (3.11) Evaluating the scalar Poissonās equation (3.6) using (3.12), we obtain the relation (ā1šš2+ š¼š2)šµš„š¼š = āšššš£š„ššµš§š¼ššš§š (3.12) Evaluating (3.13) using (3.12), we find that the š§š-direction field decay term š¼š depends only on the material properties, the magnetization harmonic characteristic wavelength ššand the velocity of the conducting slab: š¼š =1ššā1 ā ššššš£š„šš (3.13) Just as in 2.2.4, boundary conditions are obtained from the continuity of vector potential, field continuity based on Gaussā law (2.4), Ampereās law (2.2) and (2.5). Boundary conditions are summarized below: šµš„,š ā šµš„,š = āšššš„,š (3.14) šµš§,š ā šµš§,š = 0 (3.15) šµš„,š ā šµš„,š = šššš„,š (3.16) 154 šµš§,š ā šµš§,š = 0 (3.17) šµš„,š ā šµš„,š = 0 (3.18) šµš§,š ā šµš§,š = 0 (3.19) šµš„,š ā šµš„,ā = 0 (3.20) šµš§,š ā šµš§,ā = 0 (3.21) š“š¦,š = š“š¦,š (3.22) š“š¦,š = š“š¦,š (3.23) where the subscript denotes both the direction of the field component and the boundary location. As in 2.2.4, the field š© as well as the vector potential šØ go to zero at š§š = āā,ā. We know the vector potential must be continuous everywhere and can confirm this separately, but only use this property to evaluate the boundaries (š, š) and (š, š). We guess that š“š¦ in region I has the form š“š¦ = š š([šš¼šš§ššš + šš¼šāš§ššš ]šāšš„š/šš) (3.24) Similarly in regions III and V, we guess that š“š¦ has the form: 155 š“š¦ = š š([šš¼š¼š¼šš§ššš + šš¼š¼š¼šāš§ššš ]šāšš„š/šš) (3.25) š“š¦ = š š([šššš§ššš + šššāš§ššš ]šāšš„š/šš) (3.26) The vector potential in region II, the permanent magnet region, is defined as per (2.38) for the fundamental spatial component as š“š¦ = š“š¦š + (š“š¦,š ā š“š¦š)š ššā (š§ššš)š ššā (š»ššš) ā (š“š¦,š ā š“š¦š)š ššā (š§š ā š»ššš )š ššā(š»ššš) (3.27) where š“š¦š is the particular solution to the Poisson equation (3.2): š“š¦š = ššššššš§1 (3.28) From (2.24) and the definition of the potential amplitude š“š¦ in each region other than the conductor, the field along each boundary can be calculated as below: š©š = š š[ 1šš{āšš¼šš»ššš + šš¼šāš»ššš }šāšš„š/šš0āššš{šš¼šš»ššš + šš¼šāš»ššš }šāšš„š/šš] (3.29) š©š = š š[ 1šš{ā(š“š¦,š ā š“š¦š) ššš”ā (š»ššš) + (š“š¦,š ā š“š¦š)šš šā (š»š/šš)}0āšššš“š¦,š ] (3.30) 156 š©š = š š[ 1šš{ā(š“š¦,š ā š“š¦š) šš šā (š»ššš) + (š“š¦,š ā š“š¦š)ššš”ā (š»š/šš)}0āšššš“š¦,š ] (3.31) š©š = š š[ 1šš{āšš¼š¼š¼ + šš¼š¼š¼}šāšš„š/šš0āššš{šš¼š¼š¼ + šš¼š¼š¼}šāšš„š/šš] (3.32) š©š = š š[ 1šš{āšš¼š¼š¼šāš§ššš + šš¼š¼š¼šš§ššš}šāšš„š/šš0āššš{šš¼š¼š¼šāš§ššš + šš¼š¼š¼šš§ššš}šāšš„š/šš] (3.33) š©š = š š [āšš¼ššš{šš¼ššāš¼šš§š ā šš¼ššš¼šš§š}šāšš„š/šš0{šš¼ššāš¼šš§š + šš¼ššš¼šš§š}šāšš„š/šš] (3.34) š©š = š š [āšš¼ššš{šš¼ššāš¼š(š§š+š·š) ā šš¼ššš¼š(š§š+š·š)}šāšš„š/šš0{šš¼ššāš¼š(š§š+š·š) + šš¼ššš¼š(š§š+š·š)}šāšš„š/šš] (3.35) š©ā = š š[ 1šš{āšššā(š§š+š·š)šš + ššš(š§š+š·š)šš }šāšš„š/šš0āššš{šš¼šā(š§š+š·š)šš + šš¼š(š§š+š·š)šš }šāšš„š/šš] (3.36) 157 Solving the system of equations from the boundary conditions (3.15)-(3.24) using the field definitions (3.30)-(3.37) and the additional boundary conditions š© = 0, šØ = 0 at š§š = āā,ā, we find that the complex valued coefficient šš¼š is šš¼š = 2šµš (1 ā šāš»ššš )1(š¼ššš + 1) + (āš¼ššš + 1)(š¼ššš ā 1š¼ššš + 1)šā2š¼šš·š (3.37) where šµš =4šµššcos (š4+š2šš). The field at the conducting slab boundary is š©š = š š[ { āšš¼ššš[1 āš¼ššš ā 1š¼ššš + 1šā2š·šš¼š]0[1 +š¼ššš ā 1š¼ššš + 1šā2š·šš¼š]} šš¼ššāš§š/šššāšš„š/šš] (3.38) For a conducting slab of infinite thickness (š·š ā ā) the field along the top of the conducting slab simplifies to š©š = š š [{āšš¼ššš01} 2šµš1 ā šāš»šššš¼ššš + 1šāš§š/šššāšš„š/šš] (3.39) We can partially validate the field solution by setting the conducting slab velocity to zero (š£š„ = 0) and comparing the field distribution to that of a static Halbach array (since the conducting slab has a permeability šš, under static conditions it is considered air). It is easily shown that (3.39) becomes the field solution (2.52) for the stationary infinitely extended Halbach array with š§š =āš§š for š£š„ = 0. To evaluate the force on the conducting slab, we use Maxwellās stress tensor on the boundary shown in Figure 3.3. The total force on any volume is the integration of the Maxwell stress tensor over the surface of the volume enclosed (Chapter 8 of [65]): 158 š = āÆ š ā ššš (3.40) with the stress tensor under MQS conditions defined as š = [šš„š„ šš„š¦ šš„š§šš¦š„ šš¦š¦ šš¦š§šš§š„ šš§š¦ šš§š§] (3.41) Each term within the stress tensor (for a system with uniform permeability šš) is ššš = ššš»šš»š āšš2š»šš»ššæšš (3.42) where š is the direction of the tensor pressure and š is the normal to the surface, and the Kronecker delta šæšš is šæšš = [1 š = š0 š ā š] (3.43) and š»šš»š = š»š„2 + š»š¦2 + š»š§2 (3.44) Examining the boundary drawn around the conducting slab, we note that the bottom surface of the boundary is at š§š ā āā with š© ā 0 and will contribute nothing to the force on the slab. The boundaries at š„š = āā,ā will cancel out due to symmetry in the limit. Conveniently, this allows us to evaluate the total force on the conducting slab by evaluating only the field at the surface of the conductor, š©š. The total force on the conductor is therefore ššššš =12šš[šµš„ššµš§š012(šµš§š2 ā šµš„š2 )] ššš·š (3.45) 159 where šš is the width of the boundary and š·š is the depth of the array and conductor in the š¦š direction. Because the field is periodic in š„š, the mean force can be evaluated using the identity for spatial averaging of two periodic functions [65]: the spatial average of the product of two periodic functions šŗ = š š(?Ģ?šāšš„ššš ), š¾ = š š(?Ģ?šāšš„ššš ) where ?Ģ?, ?Ģ? are the complex valued amplitudes of each function is < šŗš¾ >=12š š[?Ģ??Ģ?ā] (3.46) Here we use the notation <> to mean a spatial averaging operation, and ?Ģ?ā is the complex conjugate of ?Ģ?. Therefore (3.46) becomes ššššš =12ššš š [?Ģ?š„š?Ģ?š§šā012(?Ģ?š§š?Ģ?š§šā ā ?Ģ?š„š?Ģ?š„šā )] ššš·š [š] (3.47) The magnetic pressure developed on the magnet array is therefore š·šššššš = [šš„0šš§] = ā12ššš š [?Ģ?š„š?Ģ?š§šā012(?Ģ?š§š?Ģ?š§šā ā ?Ģ?š„š?Ģ?š„šā )] [Pa] (3.48) The field amplitudes used to evaluate (3.48) are simply the complex valued amplitudes of the vector components in (3.38) with šš¼š defined as in (3.37): ?Ģ?š„š = 2šµš (1 ā šāš»ššš ) šāš§š/ššāšš¼ššš[1 āš¼ššš ā 1š¼ššš + 1šā2š·šš¼š](š¼ššš + 1) + (āš¼ššš + 1)(š¼ššš ā 1š¼ššš + 1)šā2š¼šš·š (3.49) ?Ģ?š§š = 2šµš (1 ā šāš»ššš ) šāš§š/šš1 +š¼ššš ā 1š¼ššš + 1šā2š·šš¼š(š¼ššš + 1) + (āš¼ššš + 1)(š¼ššš ā 1š¼ššš + 1)šā2š¼šš·š (3.50) 160 We can evaluate the closed form analytical solution for levitation and drag forces (3.48) by comparing it to the results from a finite element simulation in 2D. We use [79] in 2D mode, using an analytical array of 15š width, and a copper conducting slab with a thickness of š·š =0.055š = 2 mm (well below the limit where we would start considering the conductor thickness infinite). Model parameters are listed in Table 3.1. We determine the force on the slab via the Lorentz force integrated through the volume of the conducting slab. Figure 3.4 shows the comparison between the analytical magnetic pressure and average magnetic pressure from the 2D FEA. The error decreases as the mesh for the 2D FEA is refined; for the case shown the maximum error is 1.4% between the 2D analytical model and the FEA. The error decreases as velocity š£š„ increases. The effects of field harmonics other than the fundamental should be less than this error. Figure 3.4 Magnetic pressure on magnet array, 2D analytical model versus 2D FEA. 161 2D ASYNCHRONOUS LEVITATION MACHINE PARAMETERS Parameter Symbol Value Units Halbach wavelength š 36.13 mm Magnet remanence šµš 1.325 T Magnet height š»š 7.5 mm Halbach segment gap š 2.682 mm Conductor slab thickness š·š 2 mm Conductor slab conductivity š 5.998Ć 107 S/m Flying height š§š 1.7 mm Table 3.1 2D Asynchronous Levitation Machine Parameters. 3.2.2 Analytical Shells Model for Force and Torque of Rotating 3D Magnet Disk over Homogeneous Conducting Slab In this section we extend the linear 2D asynchronous machine model to a 3D rotating permanent magnet asynchronous motor element. Figure 3.5a shows the motor element topology with a circumferential Halbach patterned magnet disk at some height over a stationary conducting slab, with some relative rotational velocity between the two bodies. The magnet disk has a magnet height of š»š, inner radius š š, and outer radius š š. The conducting slab is below the magnet disk at a flying height of š§š, and the conductor has a thickness in the z-direction of š·š. The conductor dimensions in the XY plane are large compared to the outer diameter of the magnet disk. The magnet array coordinate frame is fixed to the magnet disk and expressed in cylindrical coordinates (šš, šš, š§š). The magnetization distribution š“ is now a periodic function of šš only and uniform in the other axes. The disk rotates at š around the z-axis relative to the stator coordinate frame 162 expressed in Cartesian coordinates (š„, š¦, š§). We define the Halbach pole number šš as the number of full wavelengths of the Halbach array. Figure 3.5b shows a differential cylindrical shell element of the magnet disk at radius šš. If we imaging cutting the circumferential shell element at some point and unwrapping it to form a linear magnet array, we can form a conceptual equivalence between the 2D linear asynchronous levitation machine discussed in 3.2.1 and the cylindrical shell element, with coordinate transforms: š„š = šš3š·šš3š· (3.51) š¦š = šš3š· (3.52) š§š = š§š3š· (3.53) where the coordinates (š„š, š¦š, š§š) are identical to the magnet frame used in 3.2.1 for the 2D infinitely extended case. We note that this shell element thus unwrapped has no end effects in š„š due to circular closure of the magnet array pattern in 3D. We can thus model the field distribution from each shell as the result of an infinitely extended 2D linear asynchronous machine with infinitesimal depth šš using (3.49) and (3.50) determine the field distribution at the top surface of the conductor assuming for the moment no flux leakage in the radial direction. Minimal leakage in the radial direction is a reasonable assumption away from the inner and outer radius of the magnet disk. The total field is obtained by superimposing the field contributions of each shell element. To obtain the total force on the conductor we can use Maxwellās stress tensor, and define a closed surface with the top surface at the conductor face š§š3š· = āš§š from šš3š· ā (š š, š š) and šš3š· ā163 (0, 2š); bottom surface at š§š3š· = āā with the same radial and angular extents; and cylindrical faces at šš3š· = š š , š š from š§š3š· ā (āš§š , āā). As per the analysis in 3.2.1, the surface at š§š3š· =āā will contribute zero force, and assuming no radial flux leakage in the first approximation, the cylindrical faces will also contribute no force. Therefore we can determine the differential force and torque generated on each shell element from the average magnetic pressures calculated in (3.48). By integrating the differential force and torque contributions from all shells along the radius of the magnet array, an approximation of the total force and drag torque on the magnet disk can be calculated, and is equivalent to evaluating the total Maxwell stress tensor. The Halbach spatial wavelength is a function of the radius šš and the Halbach pole number šš as š =2ššš3š·šš (3.54) The shell element has a tangential speed at every point šš3š·š; the 2D analytical model in 3.2.1 is defined based on a stationary magnet array and traveling conductor, therefore š£š„ = āšš3š·š (3.55) The differential levitation force in the z-direction on the shell element is obtained by plugging in the now radially dependent wavelength (3.55) and relative conductor velocity (3.56) into (3.48): šš¹š§ = 2ššš§šš3š·šš (3.56) The total force on the magnet disk is the integration š¹š§ = ā« šš¹š§š šš š= 2š ā« šš§šš3š·ššš šš š (3.57) 164 Similarly the differential torque around the z-axis is ššš§ = 2ššš„šš3š·2 šš (3.58) and the total torque is šš§ = ā« ššš§š šš š= 2š ā« šš„šš3š·2 ššš šš š (3.59) The integrations in (3.57) and (3.59) can be carried out numerically. Since each is only an integration along a single axis, compared to the 3D numerical integration required by FEA, this model still has the advantage in terms of computation time. In addition, even without solving the integrations analytically, the relative shapes of the force and torque curves can be determined by the closed form solutions to the 2D case (3.48) for parameter optimization. 165 Figure 3.5 Rotating permanent magnet disk asynchronous levitation machine. (a) Motor topology (with magnetization pattern for example similar to that shown in Figure 1.19a). (b) Cylindrical shell element of magnet disk. ` We compare this analytical shells model to a 3D finite element study. An example magnet disk is designed with parameters listed in Table 3.2. We note that the total force producing area is much larger than the inner diameter of the magnet disk, so no assumptions about magnet geometry are being used as was the case for the analytical model in [57]. Figure 3.6 shows the FEA model 166 created in [79], with a conductor diameter 10% wider than the outer diameter of the magnet array disk (increasing the conductor diameter made the FEA simulation time excessively long). This array is designed as a practically buildable magnet disk, and therefore each of the magnet segments magnetized in the XY plane (all others are magnetized in the Z-direction) are straight magnetized and not circumferentially magnetized as shown in Figure 3.7. We wish to evaluate the model in the context of a practical design and not as a descriptive model for an impractical case. Figure 3.8 compares the levitation force and torque from the analytical shells model (with integrations in (3.57) and (3.59) carried out numerically with 104 points) to the 3D FEA at different rotational speeds (where rotational speed in revolutions-per-minute (rpm) is 60|š|2š). The maximum error is under 15%, and is due in part to the assumption of no leakage flux at the inner and outer radii of the magnet array, the relatively small XY extent of the conductor slab, and the straight magnetization of the magnet segments. For engineering design applications, this level of error in the magnitude of the force and torque is acceptable, especially as the shape of the load characteristic (e.g. the saturation rate of the levitation force) is accurate to less than 1%. This is predicted from the form of the expressions (3.58), (3.60) because the rotational velocity dependent terms are easily separable from the magnitude terms. We conclude that the analytical shells model can be used to design practical magnet array disks. 167 Figure 3.6 3D FEA model, with šµš = š. 168 Figure 3.7 Circumferential versus straight magnetization of magnet segments in XY plane. Figure 3.8 Levitation force and drag torque, 3D FEA versus analytical shells method for šš = š šš. 169 3D ROTATING DISK ASYNCHRONOUS LEVITATION MACHINE PARAMETERS Parameter Symbol Value Units Outer magnet array radius š š 100 mm Inner magnet array radius š š 20 mm Magnet height š»š 6.35 mm Magnet remanence šµš 1.2347 T Halbach segment gap š 0 mm Magnet array pole number šš 8 -- Conductor slab thickness š·š 25.4 mm Conductor slab conductivity š 5.998Ć 107 S/m Flying height š§š 2 mm Table 3.2 3D Rotating Disk Asynchronous Levitation Machine Parameters. 3.3 Experimental Force, Torque and Power for Rotating Magnet Disk An experimental magnet array disk (Figure 3.9) is designed and built to evaluate the analytical shells model. We note that the magnet segments are not true arc geometries, but rather cube magnet segments which were more easily available. We approximate the gap between magnet segments for each shell element as a linear function of the radius šš and the minimum and maximum gaps (šš, šš): š =šš ā ššš š ā š š(šš ā š š) + šš (3.60) 170 In addition the magnet segments are coated with a thin layer of conductive nickel. This will affect the asynchronous field as nickel has a conductivity of 14.3 MS/m (compared to 60 MS/m for copper) and the nickel coating is approximately 100 micron thick. The experimental load characterization setup is shown in Figure 3.10, with a CNC spindle used to locate the magnet disk in XYZ and rotate the magnet disk. A conducting slab is attached with adhesive to the top of an electrically insulated spacer which is bonded to an adapter plate and an ATI Industrial Automation MINI-45 6-DOF load cell as shown. The thickness of the insulating layer was determined by spinning the magnet array at the maximum test speed condition of 7700 rpm over the adapter plate and load cell alone and increasing the height of the disk until the measured forces and torques were below the noise floor of the sensor, then increasing the height a further 25%. This insures that only the effects of the conducting slab will be measured and there is no extraneous coupling with the load sensor and other apparatus. Two different conducting slab materials with different thicknesses were tested: copper and aluminum. The test setup parameters and magnet disk design parameters are listed in Table 3.3. 171 Figure 3.9 Experimental rotating disk asynchronous levitation machine. 172 Figure 3.10 Load test setup. 173 EXPERIMENTAL ROTATING DISK ASYNCHRONOUS LEVITATION MACHINE PARAMETERS Parameter Symbol Value Units Outer magnet array radius š š 26.30 mm Inner magnet array radius š š 16.70 mm Magnet height š»š 6.35 mm Magnet remanence šµš 1.2347 T Halbach segment gap, maximum šš 3.79 mm Halbach segment gap, minimum šš 1.21 mm Magnet array pole number šš 4 -- Conductor slab thickness š·š (Copper) 25.4 mm š·š (Aluminum) 21.5 mm Conductor slab conductivity š (Copper) 5.998Ć 107 S/m š (Aluminum) 3.5Ć 107 S/m Table 3.3 3D Experimental rotating disk asynchronous levitation machine parameters. Figure 3.11 and Figure 3.12 shows the 6-DOF load characteristic of the magnet disk over copper and aluminum respectively. The disk generates large levitation force š¹š§ and drag torque šš§ as predicted. The torques šš„, šš¦ are non-zero and detectable at higher rotational speeds, and some very small in-plane XY forces on the order of the load sensor noise floor. This is because the magnet array disk rotation axis is out of perpendicularity to the surface of the conducting slab for each case to within 250 micron, and therefore there will be a net torque generated on an axis in the 174 XY-plane as well as some net translation force. The torque is at least an order of magnitude less than the drag torque šš§. This is in line with our conceptual understanding of the asynchronous levitation motor element as discussed in 3.1. The levitation force generated on the copper slab versus the aluminum slab for a given height and rotation speed is larger by approximately the ratio of the conductivity of copper to aluminum; this agrees with the findings from the experimental work in [62]. Figure 3.11 6-DOF load characteristics for rotating disk asynchronous machine (copper). 175 Figure 3.12 6-DOF load characteristics for rotating disk asynchronous machine (aluminum). Figure 3.13 and Figure 3.14 shows a comparison of the experimental load characteristic with derived mechanical power and the analytical shells model at iso-lines of constant flying height (with integrations in (3.57) and (3.59) carried out numerically with 104 points) for the copper and aluminum cases respectively. The levitation force shows good agreement at lower flying heights, while the drag torque shows good agreement at higher flying heights. Figure 3.15 shows the error 176 at all measured points between the analytical model and the experimental measurement. The worst case error in levitation force prediction is 17.82% for the copper slab; in drag torque -29.8%; the power, since it is derived from the spindle speed and measured drag torque, ššš¤šš = |š|šš§, has the same worst case error at -29.8%. Similarly, the worst case error for levitation force prediction is 16.4% for the aluminum slab and 26.3% for the drag torque and power. This is relatively good agreement with the experimental case, given the assumptions about circumferential magnetization, arc geometry, leakage effects and the presence of the nickel coating. In addition, the magnet disk chassis is aluminum with conductivity of 36.9 MS/m, which will affect the induced field. The magnet disk field may also have been affected by the magnetic spindle chuck. We conclude that the analytical shells method is a useful design tool for real magnet array disks, without any need for empirical data fitting. 177 Figure 3.13 Force, torque and mechanical power iso-lines at different flying heights (copper). 178 Figure 3.14 Force, torque and mechanical power iso-lines at different flying heights (aluminum). 179 Figure 3.15 Error between analytical shells model and experimental load characteristic (copper). Figure 3.16 Error between analytical shells model and experimental load characteristic (aluminum). 180 Chapter 4: Conclusion In this concluding chapter, I outline the original contributions of this thesis to the state of the art and suggest future work based on the presented research. Non-contact 6-DOF direct drives have the potential to greatly improve motion stage performance by eliminating friction, bearing/guiding errors and the inertia of intermediate motion elements that are used in traditional stage designs. Magnetic field interactions can be used to apply these non-contact forces to a single mover body. In this thesis I have explored two classes of magnetic non-contact multi-DOF planar motion stages: i) synchronous planar levitation stages and ii) asynchronous planar levitation stages. Permanent magnet-based synchronous levitation stages have the potential to be highly power efficient if the motor topology allows for high copper fill factors and utilization of the permanent magnet field. The novel 6-DOF motor topology presented in this thesis achieves these design goals in addition to minimizing controller and commutation complexity through the novel design which allows for the highly decoupled, linear 6-DOF actuation of a single mover body with very simple control and commutation architecture. Compared to other synchronous planar levitation motors, this novel 6-DOF actuator can be scaled up in planar stroke with minimal increase in the controller complexity or number of current drives required. The high copper fill factor of the coil design, in addition to increasing the overall power efficiency of the motor, also produces very high force densities and thus accelerations. One limitation to scaling up this synchronous motor technology is that this technology relies on high performance 6-DOF metrology solutions (such as planar encoders and multi-axis laser interferometers) for large planar work volumes. The nature of magnetic levitation systems 181 means that the final positioning performance is in large part determined by the performance of the sensor itself. Force generation capabilities are also limited by the cooling solution of the motor, due to excessive heat dissipation at high current working conditions. However, the heat source in the described synchronous moving magnet planar motor is stationary and thus the required active cooling method is easier to incorporate into existing stages than moving coil motors with a necessarily moving heat source. The development of the stationary cooling architecture is beyond the scope of the work in this thesis, but should be relatively simple. While synchronous planar levitation motors have intrinsically better power efficiency compared to asynchronous induction type planar levitation motors, even the minimal complexity of driving alternating 1D coil structures can become prohibitive for ultra-long stroke planar applications. Asynchronous planar levitation motors have several advantages over synchronous levitation techniques: simple homogeneous stator structures (essentially just slabs of conducting material) allow for cost effective ultra-long planar work areas, and the passively stable nature of the 6-DOF actuation means feedback sensors are not required for stable motion. The disadvantages of this type of levitation motor are: i) the power efficiency, especially at higher loads, is limited because of the Ohmic losses induced in the conducting stator plus the power loss in the rotation mechanism for each magnet disk; ii) the coupled nature of all axes means the 6-DOF motion control is not simple; iii) the mover itself must have some active means of mechanically rotating the permanent magnet disks relative to the motor slab, as well as tilting them in tip/tilt axes. These latter requirements mean that either the mover uses some battery to power onboard rotary motors, or an umbilical cable must be attached to the mover to provide power. Nevertheless it may be beneficial to choose asynchronous levitation techniques for certain applications where the complexity of synchronous driving and/or the cost of the stator may be the limiting factor. 182 4.1 Synchronous Planar Levitation Motor Contributions I, II and III have been in part published in [89]. Parts of contributions IV and V have been published as [90]. 4.1.1 Contribution I: Motor topology for New Type of 6-DOF Synchronous Planar Levitation Motor I presented a novel long stroke moving magnet 6-DOF synchronous planar levitation motor that uses symmetric linear 1D magnet arrays interacting with coil distributions in alternating layers of straight 1D coil patterns manufactured as printed circuit boards. This motion stage has the following advantages compared to existing 6-DOF synchronous levitation planar stages: the mover can be much smaller than the planar stroke and therefore the utilization of the work volume can be much higher; the number of coils increases linearly with planar stroke, compared to quadratically for existing 6-DOF synchronous levitation motors; the coil fill factor can be very high making this motor very power efficient with high force density (critical requirements for high precision applications); the driving method (with coil excitation zones deliberately wider than the magnet arrays) minimizes magnet array edge effects; the coil topology (where all end turns and coil return paths are on the perimeter of the working zone) minimizes coil edge effects; and this natural force linearity combined with the natural force decoupling of each 2-DOF motor element makes control and commutation of the stage very simple. 183 4.1.2 Contribution II: Analytical Model of Field, Force, Torque and Commutation Laws for New Type of 6-DOF Synchronous Planar Levitation Motor In this thesis I have presented new analysis and modeling of the relevant magnetic field, force, torque, and commutation laws of this new 6-DOF motion stage. The novel motor topology and driving method allows the application of an analytical force model of minimal complexity, which further allows us to generate quite simple commutation laws and feedback motion controllers based on a rigid body plant model. Simple commutation laws with low computational cost for 6-DOF actuation are highly desired for practical implementations of the motion stage. 4.1.3 Contribution III: Prototype and Experimental Demonstration of New Type of 6-DOF Synchronous Planar Levitation Motor An experimental alpha prototype has been built and successfully demonstrated 6-DOF planar actuation over a 60mm by 260 mm planar stroke with 11.6 mm maximum out of plane stroke and small rotations around all axes. Position feedback control has been implemented with 50 Hz cross-over frequencies for all axes with ~100 Hz closed loop bandwidth and 0.9 N/šš minimum dynamic stiffness. This control bandwidth is limited by the position sensor performance. Regulation positioning results (i.e. with a stationary reference position command in all axes) show RMS error in the three translation axes of 2.0 šš, 2.1 šš, 5.1 šš in š„, š¦, š§ respectively, and RMS regulation error in the three rotational axes of 0.0014ā, 0.0013ā, 0.0004ā in š¼, š½, š¾. These results are consistent with the noise floor of the position sensor. For an XY elliptical reference path of 260 mm X-stroke and 60 mm Y-stroke at a flying height of 1mm above the coil array PCB surface, the tracking RMS error in the three translation 184 axes is 7.49 šš, 5.02 šš, 9.20 šš in š„, š¦, š§ respectively, and RMS tracking error in the three rotational axes of 0.0066ā, 0.0082ā, 0.0029ā in š¼, š½, š¾ respectively. The prototype presented in this thesis was designed for in-plane accelerations of 5.2g in the strong axis and 4.9g in the weak axis, and out of plane (z-directed) accelerations of 9.9g for a maximum current density of 20A/mm2 and a 0.5 mm flying height, with a total mover mass of 2.3 kg. The maximum flying height of 11.6 mm is determined by the ability to exert 1g levitation acceleration at this maximum current density. This current density limit may be exceeded through the use of forced cooling solutions. 4.1.4 Contribution IV: Novel Split and Quad-Split Magnet Array Designs for Force and Torque Ripple Self-attenuation I modeled sources of force and torque ripple arising from the interaction of higher order field harmonics of the permanent magnet field of each array and the commutated current in the discrete 1D coil pattern. I presented novel magnet array designs to minimize targeted force and torque ripples by splitting and offsetting sections of the magnet array relative to a common commutation center. These designs attenuate force and torque ripple without changing either commutation or control algorithms, and without increasing the complexity of the coil design. Finite element simulations of the novel quad-split magnet array showed that both force and torque ripples at the targeted spatial frequency were attenuated to within the numerical accuracy of the simulation for force and by four orders of magnitude for torque; this is at the cost of 10% mean force generation. This method of splitting magnet arrays can be generalized to other synchronous machines, in particular other synchronous planar levitation motors, ironless linear motors where 185 the primary source of force ripple is due to field and current interactions, and can also be extended to rotary synchronous machines. 4.1.5 Contribution V: Experimental Verification of Novel Split Magnet Array Design for Force Ripple Attenuation Experimental load testing of the novel split array was carried out and showed an order of magnitude reduction in the dominant 6-cpš force ripple at the expense of 5% mean force generation. This verifies the practicality of the split array methodology for targeted force ripple attenuation. 4.2 Asynchronous Planar Levitation Motor These contributions have been submitted as a paper to IEEE Transactions on Magnetics and is currently under review. 4.2.1 Contribution I: Novel Analytical Model of Field, Force and Torque of Magnet Disk Rotating Above a Finite Thickness Homogenous Conductor I presented a novel analytical model of a rotating permanent magnet disks over a stationary homogenous conducting slab to produce actuation forces and torques; this magnet disk can be used in a 6-DOF asynchronous planar levitation motor. While there is significant force and torque coupling between axes, the advantages of this asynchronous topology are a combined levitation and propulsion system and a very simple stator structure, as well as passive stability of all axes. The novel analytical expression for levitation force and drag torque is based on an integration of differential shell elements equivalent to multiple 2D infinitely extended asynchronous levitation 186 motors. This force and torque model allows for quick optimization of magnet disk and conducting slab design, and can offer physical and design insight into motor operation based on the form of the expression. This is advantageous compared to existing force and torque models in literature, which are primarily based on numerical finite element methods. Of the two relevant analytical models presented in the literature, one was restricted to non-permanent magnet sources [57] and limited in terms of applicable geometry, and the other used parametric empirical curve fitting to fit the model to experiment with two parameters making it less useful for design [72]. 4.2.2 Contribution II: Evaluation of Novel Analytical Model with Finite Element Simulation and Characterization with Experimental Load Measurements I presented comparisons of this analytical model to a 3D finite element model showing a worst case error of 15% in both levitation force and drag torque. I built and measured a real experimental magnet disk with cuboidal magnet segments and showed a worst case force and torque modeling error of under 18% in levitation force and under 30% in drag torque for a copper conducting slab, and under 17% for levitation force and under 27% for drag torque for an aluminum conducting slab. This analytical model is therefore useful in the design of an asynchronous planar levitation motor using rotating circumferential magnet disks. 4.3 Future Work Common future work for both classes of planar levitation motor should include work on multiple simultaneous movers in the same work volume. 187 4.3.1 Future Work for 6-DOF Synchronous Planar Levitation Motor An important part of the future work on this class of levitation motor is the study of different motor topologies using the same 2-DOF motor elements. Since the 2-DOF motor is a combination of 1D magnet arrays and 1D coils, it is easy to conceptualize the transformation of these 1D patterns on the surface of flat planes onto curved surfaces e.g. spherical surfaces, or curved in plane (non-straight coils). This transformability of the 2-DOF motor has already been implemented for a rotary version of this synchronous levitation motor by Mark Dyck in [91] as the rotary table for a micro-machining application (a topology first disclosed in [75]). Where the synchronous planar levitation motor presented in this thesis has large translation strokes and small rotation range, the rotary version has infinite rotation range in š and small translations in the XY plane. We can conceptualize other such versions of this synchronous levitation motor, possibly for applications where large rotations around multiple axes with small translations are required, similar to applications for the spherical motor presented in [22]. Other possible future work on this synchronous levitation motor include a study of the dynamic switching algorithm of the active coil zones and the effect on mover actuation force; minimal constraint actuation where only three of the quadrants are working (which would have applications with multiple movers partially overlapping in one of the in-plane axes); over-actuation of a large mover with more than four 2-DOF motor elements attached, allowing the control of large structural modes (as outlined in [75]); and topology, control and commutation methods for a synchronous planar motor with simultaneous large rotations and translations. Force and torque ripple self-attenuation of multiple simultaneous harmonics should be analyzed in the context of magnet array designs presented by this author in [75]. An analysis of the effects of skewed and patterned coil shapes on force and torque ripple attenuation has yet to be 188 carried out. Multiple such coil designs have been presented by Dr. Xiaodong Lu and the author in [75]. 4.3.2 Future Work for 6-DOF Asynchronous Planar Levitation Motor Modeling and Analysis Future modeling work for the asynchronous 3-DOF motor element should include expanding the novel analytical shells model to include the effects of tilt/tip in order to model propulsion in the XY plane. In particular, the in-plane translation velocity should be included, similar to the work done numerically on rotating and translating magnet wheels in [61]. Combining a full 3-DOF analytical model with combinations of multiple magnet disks, a 6-DOF dynamic model can be produced. The 6-DOF characterization of this planar motor can be used to linearize and uncouple the axes of the stage for simplified control and actuation. Further, the effects of a conductive magnet disk chassis (a reasonable and common configuration for a practical planar motor) should be modeled, as well as the effect of a conductive coating on the permanent magnet. A more detailed analysis should also include additional field harmonics and their contribution to total force and torque. 189 Bibliography [1] N. Rao, A novel high-speed stereo-vision system for real-time position sensing (Master's Thesis), Vancouver: University of British Columbia, 2011. [2] K. Smeds, High-performance real-time motion control for precision systems (Master's Thesis), Vancouver: University of British Columbia, 2011. [3] N. Fujii, G. Hayashi and Y. Sakamoto, "Characteristics of magnetic lift, propulsion and guidance by using magnet wheels with rotating permanent magnets," in Industry Applications Conference, 2000. Conference Record of the 2000 IEEE, 2000. [4] F. de Jong, B. van der Pasch, T. Castenmiller, B. Vleeming, R. Droste and F. van de Mast, "Enabling the lithography roadmap: an immersion tool based on a novel stage positioning system," Optical Microlithography XXII Proceedings of SPIE, vol. 7274, no. 1, pp. 72741S1-72741S10, 2009. [5] J. Heilala and P. Voho, "Modular reconfigurable flexible final assembly systems," Assembly Automation, vol. 21, no. 1, pp. 20-30, 2001. [6] Z. M. Bi, S. Y. Lang, M. Verner and P. Orban, "Development of reconfigurable machines," International Journal of Production Research, vol. 46, no. 4, pp. 967-992, 2008. [7] N. Fujii, K. Naotsuka, K. Ogawa and T. Matsumoto, "Basic characteristics of magnet wheels with rotating permanent magnets," in Industry Applications Society Annual Meeting, Conference Record of the 1994 IEEE, 1994. 190 [8] D. G. Henderson, "Magnetic levitation of a stationary or moving object". US Patent 20140265690A1, 2014. [9] Hendo, "Hendo Hover," 2015. [Online]. Available: http://hendohover.com/. [Accessed 3 June 2015]. [10] K. Erkorkmaz, J. Gorniak and D. Gordon, "Precision machine tool XāY stage utilizing a planar air bearing arrangement," CIRP Annals, vol. 59, no. 1, pp. 425-428, 2010. [11] H. Shinno, H. Yoshioka and H. Sawano, "A newly developed long range positioning table system with a sub-nanometer resolution," CIRP Annals, vol. 60, no. 1, pp. 403-406, 2011. [12] M. Weck and D. Staimer, "Parallel kinematic machine tools-current state and future potentials," CIRP Annals, vol. 57, no. 1, pp. 671-683, 2002. [13] G. Pritschow, "Parallel kinematic machines (PKM)," CIRP Annals, vol. 49, no. 1, pp. 275-280, 2000. [14] W. In, S. Lee, J. Jeong and J. Kim, "Design of a planar-type high speed parallel mechanism posiitioning platform with the capability of 180 degrees orientation," CIRP Annals, vol. 57, no. 1, pp. 421-424, 2008. [15] Y. Altintas, A. Verl, C. Brecher, L. Uriarte and G. Pritschow, "Machine tool feed drives," CIRP Annals, vol. 60, no. 2, pp. 779-796, 2011. [16] T. Akiba, O. Ueda, M. Koyama and S. Sasanuma, "Optical pickup device having a magnetic circuit whose poles are configured along a common path within one pole within a focus coil". US Patent 5663841, 1997. 191 [17] X. Lu, P. Matthew, U. Irfan-ur-rab, M. Benito and S. Kristofer, "Rotary-axial spidnles for ultra-precision machining," CIRP Annals, vol. 58, no. 1, pp. 323-326, 2009. [18] I.-u.-r. Usman, M. Paone, K. Smeds and X. Lu, "Radially biased axial magnetic bearings/motors for precision rotary-axial spindles," IEEE/ASME Transactions on Mechatronics, vol. 16, no. 3, pp. 411-420, 2011. [19] B. Sawyer.US Patent 3376578, 1968. [20] H. Shinno, H. Yoshioka and K. Taniguchi, "A newly developed linear motor-driven aerostaic x-y planar motion table system for nano-machining," CIRP Annals, vol. 56, no. 1, pp. 369-372, 2007. [21] W. Gao, S. Dejima, H. Yanai, K. Katakura, S. Kiyono and Y. Tomita, "A surface motor-driven planar motion stage integrated with an XY\thetaZ surface encoder for precision positioning," Precision Engineering, vol. 28, no. 3, pp. 329-337, 2004. [22] K. Lee, R. Roth and J. Zhou, "Dynamic modeling and control of a ball-joint-like variable-reluctance spherical motor," Journal of Dynamic Systems, Measurement, and Control, vol. 118, no. 1, pp. 29-40, 1996. [23] E. h. Weck, T. Reinartz, G. Henneberger and R. De Doncker, "Design of a spherical motor with three degrees of freedom," CIRP Annals, vol. 49, no. 1, pp. 289-294, 2000. [24] M. P. Kummer, J. J. Abbot, B. E. Kratochvil, R. Borer, A. Sengul and B. J. Nelson, "Octomag: an electromagnetic system for 5-DOF wireless micromanipulation," IEEE Transactions on Robotics, vol. 26, no. 6, pp. 1006-1017, 2010. 192 [25] S. Earnshaw, "On the Nature of the Molecular Forces which Regulate the Constitution of the Luminiferous Ether," Transactions Of the Cambridge Philosophical Society, vol. 7, no. 1, pp. 97-112, 1842. [26] M. Simon and A. Geim, "Diamagnetic levitation: flying frogs and floating magnets," Journal of Applied Physics, vol. 87, no. 9, pp. 6200-6204, 2000. [27] R. E. Pelrine, "Diamagnetic levitation," American Scientist, vol. 92, no. 5, pp. 428-435, 2004. [28] R. Pelrine, A. Wong-Foy, B. McCoy, D. Holeman, R. Mahoney, G. Myers, J. Herson and T. Low, "Diamagnetically levitated robots: an approach to massively parallel robotic systems with unusual motion properties," in 2012 IEEE International Conference on Robotics and Automation, 2012. [29] SRI International, "Magnetically actutated micro-robots for advanced manipulation applications," 2014. [30] SRI International, "MicroFactory for Smart Manufacturing," 2015. [Online]. Available: http://www.sri.com/work/projects/microfactories-for-smart-manufacturing. [Accessed 6 July 2015]. [31] R. Hollis, S. salcudean and A. Allan, "A six-degree-of-freedom magnetically levitated variable compliance fine-motion wrist: design, mdoeling and control," IEEE Transactions on Robotics and Automation, vol. 7, no. 3, pp. 320-332, 1991. 193 [32] S. Verma, W. Kim and J. Gu, "Six-axis nanopositioning device with precision magnetic levitation technology," IEEE Transactions on Mechatronics, vol. 9, no. 2, pp. 384-391, 2004. [33] Z. Zhang and C.-H. Menq, "Six-axis magnetic levitation and motion control," IEEE Transactions on Robotics, vol. 23, no. 2, pp. 196-205, 2007. [34] S. Verma, H. Shakir and W.-J. Kim, "Novel electromagnetic actuation scheme for multiaxis nanopositioning," IEEE Transactions on Magnetics, vol. 42, no. 8, pp. 2052-2062, 2006. [35] W. Kim and D. Trumper, "High-precision magnetic levitation stage for photolithography," Precision Engineering, vol. 22, no. 2, pp. 66-77, 1998. [36] M. Holmes, R. Hocken and D. Trumper, "The long-range scanning stage: a novel platform for scnaned-probe microscopy," Precision Engineering, vol. 24, no. 3, pp. 191-209, 2000. [37] I. Etxaniz, A. Izpizua, M. SanMartin and J. Arana, "Magnetic levitated 2D fast drive," IEEJ Transactions on Industry Applications, vol. 126, no. 12, pp. 1678-1681, 2006. [38] J. Compter, "Electro-dynamic planar motor," Precision Engineering, vol. 28, no. 2, pp. 171-180, 2003. [39] P. C. M. Frissen, J. C. Compter, M. J. M. Renkens, G. A. J. De Fockert and R. J. M. Coolen, "Displacement device". US Patent 6847134 B2, 2005. [40] ASML, "TWINSCAN NXT:1950i," ASML, [Online]. Available: http://www.asml.com/asml/show.do?lang=EN&ctx=46772&dfp_product_id=822. [Accessed 29 August 2015]. 194 [41] R.-H. M. Schmidt, "Ultra-precision engineering in lithographic exposure equipment for the semiconductor industry," Schmidt, R. H. M. (2012). Ultra-precision engineering in lithographic exposure equipment for the semiconductor industry. Philosophical Transactions of the Royal Society of London A: Mathematical, Physical and Engineering Sciences, vol. 370, no. 1973, pp. 3950-3972, 2012. [42] M. Binnard, "Planar motor with linear coil arrays". US Patent 6445093 B1, 2000. [43] J. Jansen, C. van Lierop, E. Lomonova and A. Vandenput, "Modeling of magnetically leviated planar actuators with moving magnets," IEEE Transactions on Magnetics, vol. 43, no. 1, pp. 15-25, 2007. [44] J. Compter, P. Frissen and J. van Eijk, "Displacement device". US Patent 7948122, 2011. [45] J. Jansen, C. Lierop, E. Lomonova and A. Vandenput, "Magnetically levitated planar actuator with moving magnets," IEEE Transactions on Industry Applications, vol. 44, no. 4, pp. 1108-1115, 2008. [46] J. Rovers, J. Achterberg, M. J. C. Ronde, J. W. Jansen, J. C. Compter, E. A. Lomonova, C. M. M. van Lierop and M. J. G. van den Molengraft, "The deformation of the moving magnet plate of a commutated magnetically levitated planar actuator," Mechatronics, vol. 23, no. 1, pp. 233-239, 2013. [47] J. M. M. Rovers, M. Stock, J. W. Jansen, C. M. M. van Lierop, E. A. Lomonova and Y. Perriard, "Real-time 3D thermal modeling of a magnetically levitated planar actuator," Mechatronics, vol. 23, no. 2, pp. 240-246, 2013. 195 [48] J. Rovers, J. Jansen and . E. Lomonova, "Force and torque errors due to manufacturing tolerances in planar actuators," IEEE Transactions on Magnetics, vol. 48, no. 11, pp. 3116-3119, 2012. [49] J. Rovers, J. Jansen and E. Lomonova, "Multiphysical analysis of moving-magnet planar motor topologies," IEEE Transactions on Magnetics, vol. 49, no. 12, pp. 5730-5741, 2013. [50] J. Rovers, J. Jansen, J. Compter and E. Lomonova, "Analysis method of dynamic force and torque distribution in the magnet array of a commutated magnetically levitated planar actuator," IEEE Transactions on Industrial Electronics, vol. 59, no. 5, pp. 2157-2166, 2012. [51] C. van Lierop, J. Jansen, A. Damen, E. Lomonova, . P. van den Bosch and A. Vandenput, "Model-based commutation of a long-stroke magnetically levitated linear actuator," IEEE Transactions on Industry Applications, vol. 45, no. 6, pp. 1982-1990, 2009. [52] W.-j. Kim, "High-precision planar magnetic levitation," Massachusetts Institute of Technology, Boston, 1997. [53] Y. Ueda and H. Ohsaki, "Six-degree-of-freedom motion analysis of a planar actuator with a magnetically levitated motover by six-phase current controls," IEEE Transactions on Magnetics, vol. 44, no. 11, pp. 4301-4304, 2008. [54] P. Berkelman and M. Dzadovsky, "Magnetic levitation over large translation and rotation ranges in all directions," IEEE/ASME Transactions on Mechatronics, vol. 18, no. 1, pp. 44-52, 2013. 196 [55] M. Poloujadoff, The theory of linear induction machinery, 1 ed., Oxford: Clarendon Press, 1980. [56] M. Kumagai and R. L. Hollis, "Development and control of a three DOF planar induction motor," in 2012 IEEE International Conference on Robotics and Automation , 2012. [57] N. Fujii and M. Fujitake, "Two-dimensional drive characteristics by circular-shaped motor," IEEE Transactions on Industry Applications, vol. 35, no. 4, pp. 803-809, 1999. [58] R. F. Post and D. Ryutov, "The inductrack: a simpler approach to magnetic levitation," IEEE Transactions on Applied Superconductivity, vol. 10, no. 1, pp. 901-904, 2000. [59] J. Bird and T. Lipo, "Characterisitcs of an electrodynamic wheel using 2-D steady-state model," IEEE Transactions on Magnetics, vol. 43, no. 8, pp. 3395-3405, 2007. [60] J. Bird and T. A. Lipo, "Modeling the 3-D rotational and translational motion of a Halbach rotor above a split-sheet guideway," IEEE Transactions on Magnetics, vol. 45, no. 9, pp. 3233-3242, 2009. [61] J. Bird and T. Lipo, "A 3-D steady-state magnetic charge finite element model of an electrodynamic wheel," IEEE Transactions on Magnetics, vol. 44, no. 2, pp. 253-265, 2008. [62] N. Fujii and M. Chida, "Three dimensional force of magnet wheel with revolving permanent magnets," IEEE Transactions of Magnetics, vol. 33, no. 5, pp. 4221-4223, 1997. [63] K. Ogawa, Y. Horiuchi and N. Fujii, "Calculation of electromagnetic forces for magnet wheels," IEEE Transactions on Magnetics, vol. 33, no. 2, pp. 2069-2072, 1997. 197 [64] N. Fujii, Y. Ito and T. Yoshihara, "Characteristics of a moving magnet rotator over a conducting plate," IEEE Transactions on Magnetics, vol. 41, no. 10, pp. 3811-3813, 2005. [65] H. H. Woodson and J. R. Melcher, Electromechanical Dynamics Part II: Fields, Forces, and Motion, 1 ed., New York: John Wiley & Sons, 1968. [66] S.-M. Jang, S.-H. Lee and S.-s. Jeong, "Characteristic analysis of eddy-current brake system using the linear Halbach array," IEEE Transactions on Magnetics, vol. 38, no. 5, pp. 2994-2996, 2002. [67] K. J. W. Pluk, T. van Beek, J. Jansen and E. Lomonova, "Modeling and measurements on a finite rectangular conducting plate in an eddy current damper," IEEE Transactions on Industrial Electronics, vol. 61, no. 8, pp. 4061-4072, 2014. [68] S.-M. Jang, Y.-S. Park, S.-Y. Sung, K.-B. Lee, H.-w. Cho and D.-J. You, "Dynamic characteristics of a linear induction motor for predicting operating performance of magnetic levitation vehicles based on electromagnetic field theory," IEEE Transactions on Magnetics, vol. 47, no. 10, pp. 3673-3676, 2011. [69] T. Lubin and A. Rezzoug, "3-D analytical model for axial-flux eddy-current couplings and brakes under steady-state conditions," IEEE Transactions on Magnetics, vol. 51, no. 10, pp. 1-12, 2015. [70] O. de la Barriere, S. Hlioui, H. Ben Ahmed, M. Gabsi and M. LoBue, "3-D formal resolution of Maxwell equations for the computation of the no-load flux in an axial flux permanent-magnet synchronous machine," IEEE Transactions on Magnetics, vol. 48, no. 1, pp. 128-136, 2012. 198 [71] M.-M. Koo, J.-Y. Choi, Y.-S. Park and S.-M. Jang, "Influence of rotor overhang variation on generating performance of axial flux permanent magnet machine based on 3-D analytical method," IEEE Transactions on Magnetics, vol. 50, no. 11, pp. 1-5, 2014. [72] J. H. Park, J. H. Choi, D. H. Kim and Y. S. Baek, "Parameteric design of the levitation mechanism for maglev planar transportation vehicle," IEEE Transactions on Magnetics, vol. 40, no. 4, pp. 3069-3071, 2004. [73] X. Lu and I.-u.-r. Usman, "6-DOF Magentically-Levitated Precision Planar Motion Stage (Planar Motor) [video file]," 2012. [74] D. Trumper, M. Williams and T. Nguyen, "Magnet arrays for synchronous machines," in IEEE Industry Applications Society Annual Meeting, 1993. [75] X. Lu and I.-u.-r. Usman, "Displacement devices aand methods for fabrication, use and control of same". PCT Patent Patent WO2013059934 A1, 2013. [76] G. Akoun and J.-P. Yonnet, "3D analytical calculation of the forces exerted betwen two cuboidal magnets," IEEE Transactions on Magnetics, Vols. MAG-20, no. 5, pp. 1962-1964, 1984. [77] J. Jansen, E. Lomonova, A. Vandenput and J. Compter, "Design tool for a 6-dof planar motor with moving permanent magnets and standstill coils," in Proceedings of the 4th International Symposium on Linear Drives for Industry Applications, Birmingham, 2003. [78] H. A. Hause and J. R. Melcher, "Introduction to Electroquasistatics and Magnetoquasistatics," in Electromagnetic Fields and Energy, Englewood Cliffs, NJ, Prentice Hall, 1989, pp. 67-70. 199 [79] M. COMSOL, "AC/DC Module; Version 4.3b (x64 Build)," 2013. [80] S. Woo and D.-G. Gweon, "Design and optimization of long stroke planar motion maglev stage using copper strip array," International Journal of Precision Engienering and Manufacturing, vol. 16, no. 3, pp. 479-485, 2015. [81] H. A. Haus and J. R. Melcher, Electromagnetic Fields and Energy, 1 ed., Englewood Cliffs: Prentice-Hall, Inc., 1989. [82] H.-W. Lee, S. Lee, S.-H. Won and J. Lee, "Optimal design of high=precision maglev system using simulation-based DOE and FEM," IEE Proceedings Electric Power Applications, vol. 153, no. 5, pp. 773-779, 2006. [83] Y. Zhang, Z. Yang, M. Yu, K. Lu, Y. Ye and X. Liu, "Analysis and design of fouble-sided air core linear servo motor with trapezoidal permanent magnets," IEEE Transactions on Magnetics, vol. 47, no. 10, pp. 3236-3239, 2008. [84] Y. Okamoto and N. Takahashi, "Minimization of driving force ripple of linear motor for rope-less elevator using topology optimization technique," Journal of Materials Processing Technology, vol. 181, no. 1-3, pp. 131-135, 2007. [85] X. Lu, N. Rao and I.-u.-r. Usman, "Siz-axis position measurement system for levitated motion stages," CIRP Annals, vol. 62, no. 1, pp. 507-510, 2013. [86] H. Allag, J. Yonnet and M. Latreche, "3D analytical calculation of forces between lienar Halbach-type permanent-magnet arrays," in Advanced Electromechanical Motion Systems & Electric Drives Joint Symposium, Lille, 2009. 200 [87] Bal-tec, "Bal-tec," 2015. [Online]. Available: https://www.precisionballs.com/index.php. [Accessed 11 August 2015]. [88] Loctite, "Technical Data Sheet Loctite EA E-120HP," 2014. [89] X. Lu and I. Usman, "6D direct-drive technology for planar motion stages," Annals of the CIRP-Manufacturing Technology, vol. 61, no. 1, pp. 359-362, 2012. [90] I.-u.-r. Usman and X. Lu, "Force ripple attenuation of 6-DOF direct drive permanent magnet planar levitating synchronous motors," IEEE Transactions on Magnetics, vol. 51, no. 12, pp. 1-8, 2015. [91] M. Dyck, Magnetically leviated six degree of freedom micro-machining rotary table (Master's Thesis), Vancouver: University of British Columbia, 2014.
You may notice some images loading slow across the Open Collections website. Thank you for your patience as we rebuild the cache to make images load faster.
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Long stroke magnetic levitation planar stages
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Long stroke magnetic levitation planar stages Usman, Irfan-ur-rab 2015
pdf
Notice for Google Chrome users:
If you are having trouble viewing or searching the PDF with Google Chrome, please download it here instead.
If you are having trouble viewing or searching the PDF with Google Chrome, please download it here instead.
Page Metadata
Item Metadata
Title | Long stroke magnetic levitation planar stages |
Creator |
Usman, Irfan-ur-rab |
Publisher | University of British Columbia |
Date Issued | 2015 |
Description | Modern positioning applications often require long strokes in multiple degrees-of-freedom (DOF). One solution to such requirements is a planar stage capable of simultaneous large strokes in X- and Y-translation. An ideal stage concept is based on a planar motor which uses non-contact multi-axis forces to directly actuate a single moving body, without any connecting elements or bearings structures that would induce additional structural modes as well as excess inertia. This thesis presents the design, analysis and experimental results of two planar motors: 1) a permanent magnet synchronous planar motor, and 2) a permanent magnet asynchronous planar motor. The permanent magnet synchronous planar motor is comprised of multiple one-dimensional (1-D) magnet arrays attached to a moving stage and multiple stationary 1-D coils built as a printed circuit board (https://www.youtube.com/watch?v=-r4Tv7GbB8o). This motor addresses many issues with existing synchronous planar levitation motors, including scalability of the XY stroke with minimal increase in controller or drive complexity, and simplified commutation and control due to natural force/torque decoupling and no coil/magnet array edge effects. Modeling, analysis and design of a prototype are shown with motion control results. To use this levitation motor for high accuracy positioning applications, force and torque characteristics of the 6-DOF stage must be highly linear in all axes in order to minimize controller effort and reduce intrinsic motor disturbances. Novel magnet array designs are presented which self-attenuate both force and torque ripple in all 6-DOF, without additional controller or drive complexity. These array designs are tested via simulation and experiment. The permanent magnet asynchronous planar motor is based on the asynchronous induction effect to simultaneously levitate and propel the motion stage. The advantages of this type of planar levitation motor are: 1) passive stability in all axes; and 2) simple slab-type homogenous stator comprised of a basic conducting material such as aluminum or copper. This thesis presents novel modeling and analysis and provides a new analytical expression for the generated levitation force and torque with geometric and material property inputs. Experimental force and torque measurements are carried out and compared to the analytical model. |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Date Available | 2020-07-31 |
Provider | Vancouver : University of British Columbia Library |
Rights | Attribution-NonCommercial-NoDerivs 2.5 Canada |
DOI | 10.14288/1.0221435 |
URI | http://hdl.handle.net/2429/55916 |
Degree |
Doctor of Philosophy - PhD |
Program |
Mechanical Engineering |
Affiliation |
Applied Science, Faculty of Mechanical Engineering, Department of |
Degree Grantor | University of British Columbia |
GraduationDate | 2016-02 |
Campus |
UBCV |
Scholarly Level | Graduate |
Rights URI | http://creativecommons.org/licenses/by-nc-nd/2.5/ca/ |
AggregatedSourceRepository | DSpace |
Download
- Media
- 24-ubc_2016_february_usman_irfanurrab.pdf [ 15.41MB ]
- Metadata
- JSON: 24-1.0221435.json
- JSON-LD: 24-1.0221435-ld.json
- RDF/XML (Pretty): 24-1.0221435-rdf.xml
- RDF/JSON: 24-1.0221435-rdf.json
- Turtle: 24-1.0221435-turtle.txt
- N-Triples: 24-1.0221435-rdf-ntriples.txt
- Original Record: 24-1.0221435-source.json
- Full Text
- 24-1.0221435-fulltext.txt
- Citation
- 24-1.0221435.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
data-media="{[{embed.selectedMedia}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
https://iiif.library.ubc.ca/presentation/dsp.24.1-0221435/manifest