28.11.2014 Views

Online proceedings - EDA Publishing Association

Online proceedings - EDA Publishing Association

Online proceedings - EDA Publishing Association

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Fig. 5. Lower view of the FEM model with permanent magnets in the<br />

opposite configuration and air volume.<br />

11-13 <br />

May 2011, Aix-en-Provence, France<br />

<br />

VI. NUMERICAL CALCULATION OF STATIC LEVITATION<br />

The experimental levitation distance of the proof mass<br />

with 10mm and 1mm was assumed for the<br />

calculation of the magnetic force that was successively<br />

compared to the gravity force in order to verify the<br />

equilibrium. As described in the next section, the<br />

experimental levitation distance of the mid-plane of the<br />

proof mass in the configuration given is L, <br />

1.0687mm; the distribution of the magnetic field from the<br />

FEM model was calculated in correspondence to the vertical<br />

positions 1mm ∆ in the discretized domain, where<br />

∆ 0.5mm is the mesh size. The magnetic forces in<br />

horizontal directions are self-balanced, that reduces the<br />

equilibrium calculation at the vertical direction.<br />

Only one half of the proof mass was considered due to the<br />

symmetry of the levitating system; this part of the graphite<br />

volume was divided in several portions (1mm wide) as<br />

described by Fig. 7. The magnetic force was computed on<br />

the nodes situated along the longitudinal axis of each<br />

portion and then extended to the volume of the entire slice.<br />

Finally, the total magnetic force acting on the proof mass<br />

was calculated by adding all the contributions.<br />

The magnetic forces are calculated in correspondence to<br />

the nodes , situated along the longitudinal axis of each<br />

graphite slice (), as indicated in Fig. 7.<br />

Fig. 7. Division of the proof mass in portions and nodal magnetic force<br />

distribution.<br />

Fig. 6. Distribution of the magnetic field components at 1mm <br />

L, calculated from FEM simulation.<br />

The central finite difference method was used to compute<br />

the derivative of the magnetic flux density reported in Eq.<br />

(10). By introducing<br />

<br />

<br />

(15)<br />

it results, in the discretized domain,<br />

,<br />

<br />

,,<br />

∆<br />

(16)<br />

where the vertical index L, 1mm corresponds to<br />

coordinate of the experimental levitation height of the proof<br />

mass mid-plane and 1 L, ∆, as represented in<br />

Fig. 8.<br />

100

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!