Advanced
Development of a Simulator of a Magnetic Suspension and Balance System
Development of a Simulator of a Magnetic Suspension and Balance System
International Journal of Aeronautical and Space Sciences. 2010. Sep, 11(3): 175-183
Copyright ©2010, The Korean Society for Aeronautical Space Science
  • Published : September 15, 2010
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Dong-Kyu Lee
Jun-Seong Lee
Jae-Hung Han
jaehunghan@kaist.edu
Yoshiyuki Kawamura
Sang-Joon Chung

Abstract
The increased demand for a higher performing magnetic suspension and balance system (MSBS) resulted in an increase in costs for the efforts necessary for achieving an improved MSBS. Therefore, MSBS performance should be predicted during the design in order to reduce risk. This paper presents the modeling and simulation of an MSBS that controls 6-degree of freedom (DOF)of an aerodynamic body within the MSBS. Permanent magnets and electromagnets were modeled as coils, and this assumption was verified by experimental results. Finally, an MSBS simulator was developed, predicting that the MSBS is able to contain the model within a bounded region as well as measure external forces acting on the body during wind tunnel tests.
Keywords
1. Introduction
Wind tunnel tests are often used observe the air flow around an aerodynamic body as well as measure the aerodynamic forces acting on an aerodynamic body.Oftentimes, mechanical supports are employed in order to fix the aerodynamic body at the center of the test section.However, these mechanical supports cause disturbances to the air flow around the aerodynamic body, decreasing the accuracy of the experimental results. Magnetic suspension and balance systems (MSBSs) provide ideal solutions for these disturbance problems. An MSBS has the ability to levitate the magnetized model at the center of the test section of a wind tunnel without any mechanical support, such as sting or strut, by using magnetic forces and moments( Figs. 1 and 2 ). Thus, an aerodynamic model, which is partially or fully composed of ferromagnetic material, can be suspended in the test section by an MSBS. One point to note is that the remaining portions of the aerodynamic model not composed of ferromagnetic materials must be made of a rigid, electrically non-conducting material in order to avoid side effects caused by electromagnetic fields. Meanwhile,an MSBS can easily change the position and attitude of the model as well as measure the external forces and moments simultaneously acting on the model while refraining from physical contact. Thus, an MSBS eliminates support interferences during wind tunnel tests (Covert, 1988).
The basic idea of MSBSs was originally proposed by Holmes
Lager Image
The arrangement of the electromagnets of the 60 cm magnetic suspension and balance system in JAXA (Higuchi et al. 2008).
(1937). He suggested that magnetic forces can be used to reduce frictional forces acting on bearings. Subsequently,an MSBS was successfully developed for wind tunnel tests by Laurenceau et al. (1957). As a result, MSBSs were mostly utilized to eliminate support interferences during wind tunnel tests. The advantage of eliminating support interferences is most helpful during wind tunnel tests under super sonic conditions. Using an MSBS, Owen et al. (2007)measured the lift and drag coefficients of a cone under a supersonic regime, thus eliminating shock interferences caused by stings or supports in other systems. Sawada et al.(2004) used the MSBS installed on the JAXA subsonic wind tunnel in order to measure the drag coefficients of a sphere.Higuchi et al. (2006, 2008) observed the air flow around several cylinders with particle image velocimetry during wind tunnel tests with the same MSBS. Additionally, Higuchi et al. (2006, 2008) also measured drag coefficients and the base pressure of these cylinders.
Lager Image
A magnetically supported cylinder model in the 60 cm magnetic suspension and balance system in JAXA (Higuchi et al.2008).
Regardless of the advancements made in regards to MSBSs,further development of an improved MSBS is still needed,requiring greater effort and increasing costs. In particular, the size of an electric power supply system and cooling system for the electromagnet is substantially increasing as demands of specific MSBS functions increase, thus increasing the costs of MSBS research and development efforts (Covert, 1988).Consequently, developing an MSBS without performance prediction carries more risks. Performance should then be predicted during the design process. In this study, magnetic fields generated by electromagnets and permanent magnets of an MSBS were modeled under certain assumptions. The calculated results were compared to experimental results as a means of verifying the assumptions. Relationships between state/control variables and magnetic forces/moments were also obtained in state space form so that the 6-degree of freedom (DOF) control simulation was performed with the system matrix of the MSBS. 6 independent proportionaldifferential controllers were used for each DOF in the simulation. Finally, 6-DOF control simulator was developed and the performance of the MSBS was successfully predicted with step-response simulation and an external force measurement simulation.
2. Mathematical Formulation
- 2.1 6-DOF equation of motion
The equation describing the motion of a rigid body in the body fixed system can be expressed as follows:
Lager Image
where
Lager Image
represent the mass of the body,the linear and angular velocities, the linear momentum and the angular momentum, respectively.
Lager Image
are the external forces and external moments acting on the body. For an aerodynamic model located in the MSBS with no air flow,the gravitational force and the magnetic forces and moments are the only external forces and moments that exist. These forces can be rewritten as
Lager Image
where F MSBSx , F MSBSy , F MSBSz denote the X, Y and Z components of the magnetic force, and similarly F MSBSx , F MSBSy, F MSBSz are the X, Y and Z components of the magnetic moment,respectively. [φ][θ][ψ] is the transformation matrix that converts the force and moment expressions from the MSBS fixed frame to the body fixed frame. Therefore, if the magnetic forces and moments generated by the MSBS were appropriately obtained, the equation of motion of the aerodynamic model can be completely defined and can be expressed in state space form by applying the small disturbance theory.
- 2.2 Magnetic force and moment
The magnetic force exerted on the point charge q , that is moving with the velocity of
Lager Image
is
Lager Image
if a magnetic field
Lager Image
exists. The magnetic field can be generated by another moving charge q 1 as follows:
Lager Image
where
Lager Image
is the position vector from q 1 to q,
Lager Image
is the velocity vector of q 1
Lager Image
is known as
Lager Image
With Eq. (3), the magnetic force acting on a conducting wire
Lager Image
generated by magnetic field
Lager Image
can be obtained. If a conducting wire
Lager Image
with cross section A is placed inside the magnetic field
Lager Image
and N point charges of quantity q move in the conducting wire with the velocity
Lager Image
then the force acting on the conducting wire is
Lager Image
Since
Lager Image
Lager Image
are parallel to each other, Eq. (5) can be reduced to
Lager Image
Lager Image
in Eq. (6) means the electric current I, and
Lager Image
can be simply expressed as
Lager Image
In the same context, when the electric current I 1 flows through a conducting wire
Lager Image
located at
Lager Image
from the origin and the electric current I 2 flows through a conducting wire dl 2 located at r 2 from the origin, the magnetic force acting on the conducting wire
Lager Image
from the origin, the magnetic force acting on the conducting wire
Lager Image
generated by magnetic field from a conducting wire
Lager Image
can be written as
Lager Image
Lager Image
Magnetic force interaction of two current circuits.
When these conducting wires
dl and dl 2 form closed curves, as shown in Fig. 3 , the magnetic force acting on closed curve No. 2 can be obtained by integrating Eq. (8) as follows:
Lager Image
form closed curves, as shown in Fig. 3 ,the magnetic force acting on closed curve No. 2 can be obtained by integrating Eq. (8) as follows:
Lager Image
Lager Image
Magnetic field generated by an arbitrary current circuit
Meanwhile, the magnetic field
Lager Image
generated by an arbitrary current circuit through which electric current I 1 flows ( Fig. 4 )can be written as follows.
Lager Image
Thus the magnetic field generated by an electromagnet can be calculated by Eq. (11) if it is modeled as a circular coil(Christy et al., 1993).
Lager Image
where N is the number of turns of the coil.
3. Modeling and Simulation
- 3.1 System description
The MSBS studied in this paper comprised a 33 cm × 40 cm test section. The system was composed of two large permanent magnets and eight electromagnets, as shown in Fig. 5 . The outer and inner diameters of the electromagnets were 180 mm and 140 mm, respectively. The thickness of the magnet was 20 mm, and the number of turns was 255. The large permanent magnets were made of samarium-cobalt,of which its remanence was 1.066 T. The large permanent magnets possessed dimensions of 150 mm as the outer diameter, 50 mm as the inner diameter, and 25 mm in thickness.
This study employed a cylinder type object, as shown in Fig. 6 , as the aerodynamic model. The mass of the model was measured by an electronic scale and the moment of inertia of the model was obtained from a 3D CAD program.The values of the mass and the moment of inertia of the aerodynamic model are shown in Table 1 . The aerodynamic model contained several small permanent magnets in
Lager Image
Schematic diagram of magnetic suspension and balance system(MSBS).
Lager Image
Moment of inertia of cylinder model obtained by using CADprogram.
order to react to the magnetic field generated by the MSBS.The small permanent magnets were made of neodymium with a remanence of 1.39 T. The permanent magnets were cylindrical in shape, possessing dimensions of 35mm in diameter and 40 mm in thickness.
Mass and moment of inertia of the cylinder model obtainedfrom CAD program
Lager Image
Mass and moment of inertia of the cylinder model obtainedfrom CAD program
The MSBS had the ability to countervail the weight of the aerodynamic model with the magnetic forces generated by the two large permanent magnets. It controlled the position and attitude of the model through active control using eight electromagnets. The changes in position and attitude of the model were continuously detected by several laser sensors, and a microprocessor provided the control signals in real-time to generate appropriate input currents for each electromagnet.
3.2 Magnetic field
Eq. (11) was used to calculate the magnetic fields generated by the electromagnet. The electromagnet was modeled as a circular coil with 255 turns and an average diameter of 160 mm. The calculated magnetic fields were compared to the experimental results for various input currents. As shown in Fig. 7 , the magnetic field in the z-direction parallel to the axial direction of the electromagnet measured at the 20 mm above the center of the electromagnet (X = 0 mm, Y = 0 mm, Z= 20 mm) and it is consistent with the simulated result.
Lager Image
Results taken from an electromagnet magnetic field simulationin comparison to experimental results X = 0 mm Y = 0 mm Z= 20 mm.
Lager Image
Results taken from an electromagnet magnetic field simulationin comparison to experimental results X = -95 mm Y = -95mm Z = 93.5 mm.
Figure 8 shows the magnetic field in a different position(X = -95 mm, Y = -95 mm, Z = 93.5 mm). Slight differences exist between the simul ated and measured magnetic fields illustrated in Fig.8 . This discrepancy was caused by the limit of the sensor resolution. The experimental results, measured by a Gauss meter with a resolution of 0.1 mT, were bounded within 0.1 mT. Through the comparisons, the conclusion follows that electromagnets can be modeled as circular
The major parameters for modeling the permanent magnets as circular coils
Lager Image
The major parameters for modeling the permanent magnets as circular coils
coils.
The large permanent magnets for the MSBS and the small permanent magnets used in the aerodynamic model were modeled as coils with some assumptions on dimensions arranged in Table 2 . The large permanent magnet was modeled as a coil and the diameter of the modeled coil was the mean value of the inner and outer diameters of the large permanent magnet. The small permanent magnet was modeled as a coil which possessed same diameter with the small permanent magnet. The number of turns was assumed as one, and each magnitude of the input current for each modeled coil was calculated using Eq. (12) with the kno wn thickness, the number of turns and the remanence of each permanent magnet (Christy et al., 1993).
Lager Image
The magnetic fields generated by the two large permanent magnets in the MSBS were modeled as a system with two coils at the top and bottom. The calculation was performed along the axis of the permanent magnets. The results from the calculation were compared to the experimental results.Figure 9 shows the magnetic field distribution for the entire domain, and Fig. 10 represents an enlargement of the same results for the operating domain. The calculated and measured results are the same on the operating range determined to be near the middle of the two permanent magnets, indicating that the large permanent magnets can be modeled as coils.It was observed that the difference between the measured and simulated results increases as the z-directional position approaches the permanent magnet. This difference originates from the assumption that permanent magnets may be modeled as coils. The magnetic fields generated by the small permanent magnet were calculated in the same manner as the large permanent magnets, represented in Figs. 11 and 12.The calculated magnetic fields correspond to the measured magnetic fields in the operating domain, indicating that the small permanent magnet can be also modeled as a coil.
- 3.3 Magnetic force
Based on the results obtained from the aforementioned
Lager Image
Results from the magmetic field simulation compared to the experimental results for two large permanent magnets entire domain.
Lager Image
Results from the magnetic field simulation compared to theexperimental results for two large permanent magnets operatingdomain.
experiments, it is believed that the magnetic fields generated by the electromagnets and the permanent magnets can be calculated by modeling these magnets as coils, allowing researchers to obtain the magnetic forces and moments existing between these magnets using Eq. (9). The two large permanent magnets of the MSBS countervail the weight of the aerodynamic model in the MSBS. Calculations show that the weight equilibrium position of the model is directly below 40 mm in z-direction from the center of the MSBS. After calculating the magnetic forces and moments for various ranges of state variables and input variables, it was observed that the relationships between the magnetic forces/moments and the state/input variables appear linear around the weight equilibrium position. Thus the force/moment derivatives with respect to the state/input variables could be extracted as shown in Table 3 , while relatively small amounts of change were neglected. These derivatives were
Lager Image
Results from the magnetic field simulation campared to the experimental results for small permanent magnet entire do-main.
Lager Image
Results from the magnetic field simulation compared to theexperimental results for a small permanent magnet operatingdomain.
used to obtain the equation of motion of the aerodynamic model in state space form so that the simulation of motion of the model can be performed efficiently.
- 3.4 6-DOF control
The equation of motion of the aerodynamic model in the MSBS, which was obtained in state space form, was imported into MATLAB Simulink to simulate 6-DOF control, as shown in Fig. 13 . Limiting factors, such as the operating range and the actuator delay of electromagnets, were modeled. The sensor noise was also modeled in the simulation block.Six independent proportional-differential controllers were employed to control each DOF, and the control gains were obtained in accordance to classical control theory.In the end, the 6-DOF control simulator of the MSBS was developed to allow performance predictability of the MSBS. As shown in Fig. 14 , a step response simulation was
Lager Image
Block diagram for the simulation of 6-degree of freedom (DOF)control of the magnetic suspension and balance system.
performed and the root mean square errors (RMSEs) for each DOF were calculated in order to quantitatively evaluate the performance. The calculated RMSEs for X, Y, Z, θX, θY and θZ were 0.0851 mm, 0.0842 mm, 0.1146 mm, 0.0911°,0.1411° and 0.0975° respectively. These results ensure that the MSBS possesses sufficient performance needed to fix the aerodynamic model around the weight equilibrium position inside of the MSBS.
- 3.5 External force measurement
An MSBS not only acts as contactless actuator but also contactless balance, facilitating the measurement of external forces acting on an aerodynamic body within the system while forcing the position and the attitude of the body to follow the demanded values. When wind tunnel tests using an MSBS are performed, the equation of motion of the aerodynamic body can be expressed as follows.
Lager Image
where m is the mass, or moment of inertia, of the body, x is the displacement of the body and i denotes the direction
The changes of the magnetic forces and moments with respect to the changes of the state variables and the input variables
Lager Image
The changes of the magnetic forces and moments with respect to the changes of the state variables and the input variables
of the equation, such as X, Y, Z, θ x , θ Y , θ z ., F ; i,t is the external force in the i-direction and t is time. During the experiment, x can be calculated from x i,t' which is measured by the position sensor of the MSBS. F iMSBS, t is proportional to the input current of the electromagnet as well as to the control commands obtained from the microprocessor of the MSBS.Thus F iAerodynamic ,t can be calculated from Eq. (13). This process was simulated with the 6-DOF control simulator of the MSBS while the static drag force was assumed to be the only aerodynamic force acting on the cylinder model. The demanded position and attitude remained constant during the wind tunnel tests. For this case, the model remained at a fixed position with a certain attitude allowing the left side of Eq. (13) to be canceled out so that the experimental results may be represented by the mean values. Therefore, the mean values of the aerodynamic forces and the magnetic forces and moments possessed the same magnitudes. Fig. 15 illustrates the simulated results, showing the relationships between the mean values of the drag force versus the mean values of the control commands for 6-DOF control in the body frame (see Fig. 6 ). It is observed that one of the mean control commands increased proportional to the drag force while the other values were near the zero value, indicating that 6 component forces acting on the body, such as aerodynamic forces, can be measured from the input commands.
4. Conclusions
This study developed a 6-DOF control simulator for an MSBS. The MSBS contained eight electromagnets and two permanent magnets as well as an aerodynamic model. The aerodynamic model remained at the weight equilibrium position within the MSBS. In addition, the model comprised
Lager Image
Step response of 6-DOF control simulation of the magnetic suspension and balance system.
Lager Image
Control commands for 6-DOF control with respect to various drag forces acting on the aerodynamic body
of several permanent magnets to handle magnetic forces and moments generated by the MSBS. The two permanent magnets countervail the weight of the model and the eight electromagnets control the 6-DOF of the model.These magnets were modeled as coils in order to calculate the magnetic forces and moments. Subsequently, the eight electromagnets control the 6-DOF of the model.These magnets were modeled as coils in order to calculate the magnetic forces and moments. Subsequently, the relationships between the magnetic forces/moments and the state/input variables were obtained in order to extract the force and moment derivatives. The equation of motion of the model was obtained in state space form and the 6-DOF control simulation was performed using this equation. The simulator developed in this study successfully predicted the performance of the MSBS. The step response simulation result indicated that the MSBS can hold the aerodynamic model around a specific position. Furthermore, wind tunnel test simulations confirmed that the MSBS can measure the external forces, such as drag forces, acting on the aerodynamic body during the tests.
Acknowledgements
The authors gratefully acknowledge the financial support provided by the Agency for Defense Development under the contract UD090082JD and KAIST Institute for Design of Complex Systems. The first and second authors would like to thank the Brain Korea 21 Project in 2010.
References
Covert E. E 1988 Magnetic suspension and balance systems IEEE Aerospace and Electronic Systems Magazine 3 14 - 22
Christy R. W , Milford F. J , Reitz J. R 1993 Foundations of Electromagnetic Theory 4th ed.Reading Addison-Wesley Pub. Co. MA
Higuchi H , van Langen P , Sawada H , Tinney C.E 2006 Axial flow over a blunt circular cylinder with andwithout shear layer reattachment Journal of Fluids andStructures 22 949 - 959
Higuchi H , Sawada H , Kato H 2008 Sting-free measurements on a magnetically supported right circular cylinder aligned with the free stream Journal of Fluid Mechanics 596 49 - 72
Holmes F. T 1937 Axial magnetic suspensions Reviewof Scientific Instruments 8 444 - 447
Owen A. K , Owen F. K 2007 Hypersonic free flight measurement techniques 22nd International Congress on Instrumentation in Aerospace Simulation Facilities Pacific Grove CA 1 - 11
Sawada H , Kunimasu T , Suda S 2004 Spheredrag measurements with the NAL 60cm MSBS Journal ofWind Engineering 98 129 - 136    DOI : 10.1016/j.jfluidstructs.2006.04.020