This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Hybrid FEM-Analytical Force and Torque Models
of a Reaction Sphere ActuatorLeopoldo Rossini, Emmanuel Onillon, Olivier Chetelat, and Yves Perriard
Abstract—This paper presents a hybrid FEM-analyticalmodel for the magnetic flux density, the force and torque ofa Reaction Sphere (RS) actuator for satellite attitude control.The RS is a permanent magnet synchronous spherical actuatorwhose rotor is magnetically levitated and can be acceleratedabout any desired axis. The spherical actuator is composed ofan 8-pole permanent magnet spherical rotor and of a 20-coilstator. Due to the highly complex geometry of the sphericalrotor, consisting of 8 bulk permanent magnet poles withtruncated spherical shape adjusted on the back-iron structurewith truncated octahedral shape, a pure analytical approachis not possible. Therefore, in this article we adopt a hybridapproach in which FEM or measured derived values are com-bined with other boundary conditions on a known analyticalstructure to derive expressions for the magnetic flux density,the force, and the torque. The Laplace equation is solved byexploiting powerful properties of spherical harmonic functionsunder rotation to derive closed-form linear expressions forall possible orientations of the rotor. The proposed modelsare experimentally validated using a developed laboratoryprototype and with finite element simulations.
Index Terms—Spherical actuator, attitude control, electro-magnetic modeling, electromagnetic forces, magnetic levitation.
I. INTRODUCTION
A Reaction Sphere (RS) has been proposed as a po-
tential alternative to Attitude and Orbit Control System
(AOCS) actuators based on Reaction Wheels (RWs) or
Control Moment Gyroscopes (CMGs) [1]. The RS can be
categorized into the family of electromagnetic spherical
actuators, mainly proposed in the robotic literature [2]–[4].
The RS consists of an 8-pole permanent magnet spherical
rotor that is magnetically levitated and can be accelerated
about any axis by a 20-pole stator with electromagnets. A
schematic illustration of the RS is depicted in Fig. 1. A
RS laboratory prototype has been manufactured to validate
force and torque analytical models [5]. In this prototype,
the 8 permanent magnet poles of the spherical rotor have
been discretized using a mosaic of 728 cylindrical magnets
to approximate the desired fundamental spherical harmonic
of the magnetic flux density. More recently, a spherical
rotor optimized to improve its manufacturability has been
designed and manufactured [6]. As illustrated in Fig. 2, this
rotor has 8 bulk permanent magnet poles with truncated
spherical shape that are parallel magnetized and adjusted
on the back-iron structure, which has truncated octahedral
This work was supported in part by the European Space Agency(ESA), under the GSTP contract 20227/06/NL/SFe and the NPI contract40000101313/10/NL/PA, and by maxon motor ag.
L. Rossini, E. Onillon, and O. Chetelat are with the SystemsDivision of the Swiss Center for Electronics and Microtechnology(CSEM), Rue Jacquet-Droz 1, CH-2000 Neuchatel, Switzerland (email:[email protected]).
Y. Perriard is Professor at the Integrated Actuators Laboratory (LAI),Ecole Polytechnique Federale de Lausanne (EPFL), Lausanne, Switzerland(email: [email protected]).
shape. The optimization of design parameters was carried out
to minimize the magnetic flux density distortion with respect
to the fundamental spherical harmonic of degree 3 and order
2 (octupole), which has interesting properties relevant to the
actuator control [5].
In this article, the magnetic flux density, the force and the
torque models of the optimized spherical rotor are derived
and proposed in Section II and Section III, respectively.
As firstly proposed in [5], the procedure to derive these
models exploits powerful properties of spherical harmonic
functions under rotation so that force and torque models for
all possible orientations of the rotor are expressed in closed
form as linear combination of mutually orthogonal force and
torque characteristic matrices. However, contrary to [5], a
pure analytical approach is no longer possible due to the
highly complex geometry of the optimized rotor. Therefore,
in this article we adopt a hybrid approach in which Finite
Element Method (FEM) or measured derived values are com-
bined with other boundary conditions on a known analytical
structure to derive expressions for the magnetic flux den-
sity, the force, and the torque. Subsequently, procedures to
compute the spherical harmonic decomposition coefficients
offline and online are proposed in Section IV and Section
V, respectively. As it will be illustrated, these coefficients
fully determine the magnetic flux density, the force, and the
toque in stator coordinates for any possible orientation of
the rotor. Finally, experimental and numerical verifications
of the proposed models using the developed prototype and
finite element simulations are proposed in Section VI.
II. MAGNETIC FLUX DENSITY MODEL
For the derivation of the magnetic flux density model,
we use the rotor and stator schematics in Fig. 3 with
the respective relevant dimensions. We are interested in
computing the magnetic flux density in the airgap (region
1), which is enclosed in the stator (region 2), which in turn
is surrounded by air (region 3).
Fig. 1: Schematic illustration of the 3-D motor of the RS,
which is composed of a 8-pole rotor (cube) and a 20-pole
The matrix A (Γ) has dimension Nm × 7 and contains the
spherical harmonic basis of degree 3 as
A (Γ)T=
12R
03 (γ1)
12R
03 (γ2) . . . 1
2R03 (γNm
)R1
3 (γ1) R13 (γ2) . . . R1
3 (γNm)
R23 (γ1) R2
3 (γ2) . . . R23 (γNm
)R3
3 (γ1) R33 (γ2) . . . R3
3 (γNm)
I13 (γ1) I13 (γ2) . . . I13 (γNm)
I23 (γ1) I23 (γ2) . . . I23 (γNm)
I33 (γ1) I33 (γ2) . . . I33 (γNm)
,
(26)
where Γ is the parameter vector of sensor spherical angular
coordinates
Γ = [γ1, γ2, . . . , γNm] , (27)
with γk = (θk, φk), k = 1, 2, . . . , Nm. Finally, using (23),
the force characteristic matrix KF in (21) can be readily
calculated as
K3F =
x1
2·K0
F,R + x2 ·K1F,R + x3 ·K
2F,R + x4 ·K
3F,R
+ x5 ·K1F,I + x6 ·K
2F,I + x7 ·K
3F,I (28)
where xk is the kth entry of the solution vector x, Km
F,R =
2Re{
Km
F,3
}
, and Km
F,I = −2Im{
Km
F,3
}
, with m = 1, 2, 3.
Although not reported, the torque characteristic matrix can
be derived as in (28). Finally, notice that the matrices Km
F,R,
Km
F,I, Km
T,R, and Km
T,I are constant and are computed offline.
VI. VERIFICATION WITH FEM SIMULATIONS AND
LABORATORY PROTOTYPE
A. Magnetic Flux Density Model
In this section, the magnetic flux density model is verified
using FEM simulations and experimental measurements for
various values of the maximum number of harmonics Nh
taken into account. The setup for the rotor magnetic flux
density measurement is reported in Fig. 5. The magnetic flux
density is measured with a Gaussmeter (MAGNETPHYSIK
FH55) at J = 8 inclination angles θ = [15, 25, ..., 85]degrees. The rotor is supported by a rigid guiding axis and
connected to an AC electric motor so that, for each angle θ,
a total of K = 600 points equally distributed on the interval
φ ∈ [0, 360] degrees are recorded. The magnetic flux density
is measured at I = 4 radial distances r = [92, 94, 96, 98]mm. Given values of the radial component of the measured
magnetic flux density B1r,meas (r, θ, φ) and values of the
developed model BNh
1r (r, θ, φ) computed as in (20), we
define the Mean Normalized Relative Error (MNRE) as
MNRE(%) =
I∑
i=1
J∑
j=1
K∑
k=1
ǫ (ri, θj , φk)
IJK· 100 %, (29)
where
ǫ (ri, θj , φk) =
∣
∣
∣BNh
1r (ri, θj , φk)−B1r,meas (ri, θj , φk)∣
∣
∣
maxj,k
|B1r,meas (ri, θj , φk)|.
(30)
Fig. 5: Setup for rotor magnetic flux density measurements.
The MNRE can be computed similarly for the remaining tan-
gential components of the flux density, which are evaluated
using FEM simulations.
In Fig. 6, we report the measured magnetic flux density
compared to FEM simulations and values derived using the
developed model at inclination angles θ equal to 55, 65, and
75 degrees. Values computed using the developed model are
reported for maximum spherical harmonics degrees Nh = 3and Nh = 19. As can be observed, measurements are in good
agreement with the expected simulated values. Furthermore,
the magnetic flux density profiles are sinusoidal although
higher order harmonics are observed for θ equal to 65 and 75
degrees. These fluctuations are due to the gaps at the vertexes
of the truncated octahedron structure (see Fig. 2). Finally,
notice that the profiles computed using the hybrid FEM-
analytical model with Nh = 3 and Nh = 19 are in good
agreement with measured values for θ equal to 55 and 65
degrees. For θ equal to 75 degrees, where fluctuations due to
the vertexes occur more singificantly, with Nh = 3 the model
can only capture the fundamental component of the profile,
which is completely approximated using Nh = 19. In Fig.
(7) we report the MNRE values computed using (29) as a
function of the maximum spherical harmonic reconstruction
degree Nh. As expected, the higher the number of harmonics
taken into account in the model is, the smaller the MNRE
becomes.
B. Force and Torque Models
In this section, force and torque direct models are verified
with FEM simulations for 5 randomly-generated orientations
of the rotor. Specifically, for each randomly-generated rotor
orientation, parametrized using ZYZ Euler angles α, β, and
γ, 3 different current vectors i are computed using the
inverse model of (16) as i =[
K3F
]+F +
[
K3T
]+T [5].
Matrices K3F
and K3T
are computed as in (21) using Nh = 3while F and T are randomly generated forces and torques
with norm equal to 98.1 N and 0.2 Nm, respectively. Finally,
+ is the matrix pseudo-inverse operator. Subsequently, these
15 current vectors i are applied to the finite element model,
to compute FFEM and TFEM, and to the developed hybrid
FEM-analytical force and torque model for various values of
Nh, to compute FNh and T
Nh . Notice that FNh and TNh
are computed with (16) using KF and KT defined as in
(21) with Nh ∈ I. To compare simulated forces FFEM and
698
0 45 90 135 180−500
−250
0
250
500
φ (deg)
Br(m
T)
θ = 55 deg
Model, Nh = 3Model, Nh = 19MeasurementsFEM
0 45 90 135 180−500
−250
0
250
500
φ (deg)
Br(m
T)
θ = 65 deg
Model, Nh = 3Model, Nh = 19MeasurementsFEM
0 45 90 135 180−500
−250
0
250
500
φ (deg)
Br(m
T)
θ = 75 deg
Model, Nh = 3Model, Nh = 19MeasurementsFEM
Fig. 6: Measured magnetic flux density compared to FEM
simulations and values derived using the developed model
at θ = 55 (top), 65 (middle), and 75 (bottom) degrees.
torques TFEM to those computed using the analytical model
FNh and T
Nh , we use the norm relative error defined, for
the force, as
∣
∣
∣
∣
‖FFEM‖−‖FNh‖‖FFEM‖
∣
∣
∣
∣
·100%. The expression for the
torque norm relative error is computed similarly. Moreover,
the angle between FFEM and FNh and between TFEM and
TNh is also used for comparison.
The mean and standard deviation of the norm relative
errors as a function of the maximum spherical harmonic
degree Nh taken into account in force and torque models are
reported in Fig. 8 (top). The mean and standard deviation of
the angle errors are reported in the same figure (bottom). As
it can be noticed, the norm relative errors and angle errors,
3 5 7 9 11 13 15 17 190
1
2
3
4
5
6
7
MNRE(%
)
Maximum spherical harmonic degree Nh
MNRE Br
MNRE Bθ
MNRE Bφ
Fig. 7: MNRE values as a function of the maximum spherical
harmonic reconstruction degree Nh
together with their variances, can be made smaller by in-
creasing the number of harmonics Nh taken in consideration
in the models.
3 5 7 9 11 13 15 17 190
2
4
6
8
10
12
14
16
Maximum spherical harmonic degree Nh
Meannorm
error(%
)
ForceTorqueσ forceσ torque
3 5 7 9 11 13 15 17 190
1
2
3
4
5
6
7
8
9
Maximum spherical harmonic degree Nh
Meanan
gleerror(deg)
ForceTorqueσ forceσ torque
Fig. 8: Force and torque mean norm (top) and mean angle
(bottom) errors as a function of the maximum spherical
harmonic degree Nh. Standard deviations are also depicted
taking into consideration the 15 configurations.
VII. CONCLUSIONS
In this article we presented the development of magnetic
flux density, force, and torque hybrid FEM-analytical models
for the optimized RS rotor. The development of the magnetic
699
flux density model uses the general analytical solution to
the Laplace equation within the RS airgap with analytical
and simulated/measured boundary conditions to determine
the particular solution. Then, based on the derived magnetic
flux density model, analytical expressions for the force and
torque were derived. The calculation of force and torque
models relies on the Lorentz force law and exploits the
superposition principle.
The spherical harmonics analysis performed on the sim-
ulated magnetic flux density in Section IV showed that,
as desired, the magnetic flux density is mainly based on
its principal spherical harmonic of degree 3 and order 2.
However, higher order components of degree 7, 9, 11, 13,
15, 17, and 19 also appear in the spectrum and are taken
into consideration in the developed hybrid FEM-analytical
models.
Experimental verifications on the magnetic flux density
have been performed and showed a good correspondence
between measured data and values computed using the
developed magnetic flux density model. It was illustrated
that the error between measured and analytical values can
be reduced by increasing the maximum spherical harmonic
degree Nh taken into account in the model.
Finally, FEM simulations have been performed to validate
the developed force and torque models. As expected, the av-
erage error between the analytical and simulated forces and
torques decreases by increasing the number of harmonics
taken into consideration in the model.
Experimental measurements to validate the developed
force and torque models together with the online technique
to estimate the spherical harmonics coefficients (as described
in Section V) are currently ongoing.
REFERENCES
[1] O. Chetelat, “Torquer apparatus,” U.S. Patent 2 010 007 303, Jan. 14,2010.
[2] L. Yan, I.-M. Chen, C. K. Lim, G. Yang, and K.-M. Lee, Design,
Modeling and Experiments of 3-DOF Electromagnetic Spherical Ac-
tuators. Reading, MA: Springer, 2011.[3] A. Gofuku, R. Sasaki, T. Yano, and M. Shibata, “Basic experimental
results of a 14-12 spherical motor,” in Power Electronics, Electrical
Drives, Automation and Motion (SPEEDAM), 2012 International
Symposium on, june 2012, pp. 832 –836.[4] W. Chen, L. Zhang, L. Yan, and J. Liu, “Design and control of
a three degree-of-freedom permanent magnet spherical actuator,”Sensors and Actuators A: Physical, vol. 180, no. 0, pp. 75 –86, 2012. [Online]. Available: http://www.sciencedirect.com/science/article/pii/S092442471200235X
[5] L. Rossini, O. Chetelat, E. Onillon, and Y. Perriard, “Force andtorque analytical models of a reaction sphere actuator based onspherical harmonic rotation and decomposition,” IEEE/ASME Trans.
Mechatronics, to appear.[6] L. Rossini, S. Mingard, A. B. ans E. Forzani, E. Onillon, and
Y. Perriard, “Rotor design optimization for a reaction sphere actuator,”in International Conference on Electrical Machines and Systems
(ICEMS), Sapporo, Japan, 2012.[7] E. P. Furlani, Permanent Magnet and Electromechanical Devices:
Materials, Analysis and Applications. San Diego: Academic Press,2001.
[8] L. Hu, K.-M. Lee, and X. Fu, “A method based on measured boundaryconditions for reconstructing the magnetic field distribution of an elec-tromagnetic mechatronic system,” IEEE/ASME Trans. Mechatronics,vol. 15, no. 4, pp. 595 – 602, 2010.
[9] D. W. Ritchie and G. J. L. Kemp, “Fast computation, rotation, andcomparison of low resolution spherical harmonic molecular surfaces,”Journal of Computational Chemistry, vol. 20, no. 4, pp. 383 – 395,1998.
[10] D. Healy, Jr., D. Rockmore, P. J. Kostelec, and S. S. B. Moore,“Ffts for the 2-sphere - improvements and variations,” The Journal
of Fourier Analysis and Applications, vol. 9, pp. 341 – 385, 1996.