A Semi-Analytical Method For The Design of Coil-Systems For Homogeneous Magnetostatic Field Generati
A Semi-Analytical Method For The Design of Coil-Systems For Homogeneous Magnetostatic Field Generati
A Semi-Analytical Method For The Design of Coil-Systems For Homogeneous Magnetostatic Field Generati
net/publication/216676054
CITATIONS READS
31 2,070
1 author:
Marco A. Azpurua
Universitat Politècnica de Catalunya
76 PUBLICATIONS 678 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
EMC-STD ( 21NRM06) Metrology for emerging electromagnetic compatibility standards View project
All content following this page was uploaded by Marco A. Azpurua on 27 May 2014.
M. A. Azpúrua*
Instituto de Ingenierı́a, Laboratorio de Electromagnetismo Aplicado,
Centro de Ingenierı́a Eléctrica y Sistemas, Caracas, Venezuela
1. INTRODUCTION
¡ ¢ Bmax − Bmin
H Bmax , Bmin , B̄ = 1 − , (4)
B̄
where Bmax and Bmin are the maximum and minimum values taken by
the magnitude of the magnetic flux density within the working volume,
respectively. Thus, H is measured as the maximum deviation of the
magnitude of magnetic flux density in relation to the average value of
~ within V .
the magnitude of B,
∗ z − hi −(y − b)
Bix2 = 2 2
q , (10)
(x − ai ) + (z − hi ) (x − ai )2 + (y − bi )2 + (z − hi )2
∗ z − hi y−b
Bix3 = 2 2
q , (11)
(x + ai ) + (z − hi ) (x + ai )2 + (y − bi )2 + (z − hi )2
Progress In Electromagnetics Research B, Vol. 37, 2012 177
∗ z − hi −(y + b)
Bix4 = 2 2
q , (12)
(x + ai ) + (z − hi ) 2 2
(x + ai ) + (y + bi ) + (z − hi )2
and,
µ0 Ii ¡ ∗ ∗ ∗ ∗
¢
Biy (x, y, z) = Biy1 + Biy2 + Biy3 + Biy4 , (13)
4π
∗ , B∗ , B∗
where, Biy1 ∗ are given by (14), (15), (16) and (17)
and Biy4
iy2 iy3
respectively.
∗ z − hi x+a
Biy1 = 2 2
q , (14)
(y − bi ) + (z − hi ) (x + ai )2 + (y − bi )2 + (z − hi )2
∗ z − hi x−a
Biy2 = 2 2
q , (15)
(y − bi ) + (z − hi ) (x − ai )2 + (y − bi )2 + (z − hi )2
∗ z − hi x−a
Biy3 = 2 2
q , (16)
(y + bi ) + (z − hi ) (x − ai )2 + (y + bi )2 + (z − hi )2
∗ z − hi x+a
Biy4 = 2 2
q , (17)
(y + bi ) + (z − hi ) (x + ai )2 + (y + bi )2 + (z − hi )2
and finally,
µ0 Ii ∗ ∗ ∗ ∗
Biz (x, y, z) = (Biz1 + Biz2 + Biz3 + Biz4 ), (18)
4π
∗ , B∗ , B∗
where, Biz1 ∗ are given by (19), (20), (21) and (22)
and Biz4
iz2 iz3
respectively.
∗ (x − ai )(y − bi )
Biz1 = q (β1 + β2 ) , (19)
(x − ai )2 + (y − bi )2 + (z − hi )2
∗ (x − ai )(y + bi )
Biz2 = q (β1 ∗ + β2 ) , (20)
2 2 2
(x − ai ) + (y + bi ) + (z − hi )
∗ (x + ai )(y − bi )
Biz3 = q (β1 + β2 ∗ ) , (21)
2 2 2
(x + ai ) + (y − bi ) + (z − hi )
∗ (x + ai )(y + bi )
Biz4 = q (β1 ∗ + β2 ∗ ) , (22)
2 2 2
(x + ai ) + (y + bi ) + (z − hi )
where β1 , β2 , β1 ∗ and β2 ∗ are given by,
1
β1 = (23)
(y − bi ) + (z − hi )2
2
178 Azpúrua
1
β2 = (24)
(x − ai ) + (z − hi )2
2
1
β1 ∗ = (25)
(y + bi ) + (z − hi )2
2
1
β2 ∗ = (26)
(x + ai ) + (z − hi )2
2
12
Magnitude of the Magnetic Flux Density, B (mT)
Imax 0.5 Imax 0.25 Imax 0.1 Imax Imax 0.5 Imax 0.25 Imax 0.1 Imax
12
B: 9.94 B: 10
10 z: 0.5 z: 0
B: 9.946 B: 10.01
z: 0.5 z: 0 10
8
8
6 B: 5.004
z: 0 6 B: 5.002
z: 0
4
B: 2.502 4
B: 2.501
z: 0
z: 0
2 B: 1.001 B: 1
z: 0
2
z: 0
0 0
-1 - 0. 5 0 0.5 1 -1 - 0. 5 0 0.5 1
Position on the longitudinal axis, z (m) Position on the longitudinal axis, z (m)
(a) (b)
Figure 5. Magnitude of the magnetic flux density of the designed coil
system. (a) (x, y) = (0, 0). (b) (x, y) = (±0.25, ±0.25) m.
184 Azpúrua
8. VALIDATION OF RESULTS
11 12
10 11
9 10
8 9
7 8
6 7
5 6
4 5
FEM Simulated B FEM Simulated B
Synthesis Results Synthesis Results
3 4
-1 - 0. 5 0 0.5 1 -1 - 0. 5 0 0.5 1
Posición sobre el eje longitudinal z [m] Position on the longitudinal axis, z (m)
(a) (b)
9. CONCLUSION
An alternative method for the design of coils systems for the generation
of magnetostatic fields has been proposed in this work. The
spatial distribution of the magnetostatic field is estimated superposing
the magnetic induction numerically computed from the analytical
expression of the magnetic field generated by each coil, obtained using
the Biot-Savart’s law and the current filament method. The estimated
magnetic field using this method, coincides with the calculated using
FEM, within the working volume.
Compared with other approaches used in the design of coils
systems, the proposed method combines the speed of computation
associated with the evaluation of the analytical expressions previously
known (for the most common geometries used in the construction of
windings) and the possibility of synthesizing the desired field within
Progress In Electromagnetics Research B, Vol. 37, 2012 187
REFERENCES