Regular paper

PERMANENT MAGNET HOMOGENEOUSLY MAGNETIZED ALONG ITS AXES Ana N. MLADENOVIĆ1 Abstract: The paper presents the magnetic field calculation of permanent magnet, homogeneously magnetized along its axes. Method used in the paper is based on a system of equivalent magnetic dipoles. The results that are obtained using this analytical method are compared with results obtained using program FEMLAB. Magnetic field and magnetic flux density distributions of permanent magnet are shown in the paper. Keywords: Magnetic field, Permanent magnet, Magnetic dipole.

INTRODUCTION To determine the magnetic field components in vicinity of permanent magnets, starts from presumption that magnetization, M , of permanent magnet is known. The following methods are useful in practical calculation: a) Method based on determining distribution of microscopic Ampere's current; b) Method based on Poisson's and Laplace's equations, determining magnetic scalar potential; and c) Method based on a system of equivalent magnetic dipoles. Magnetic field inside and outside the permanent magnet, if magnetization of permanent magnet is known, can be calculated using equivalent system of volume and surface microscopic Ampere's currents, which are determined as J a (r ') = rot M (r ') , and

(1)

J sa (r ') = M (r ') × nˆ ,

(2)

B (r ) = rot A(r ) .

(4)

Inside a permanent magnet, magnetic field can be determined using relation H = B µ0 − M . (5) Outside a permanent magnet, magnetic field can be determined using relation H = B µ0 . (6) The second method is based on determining magnetic scalar potential ϕ m . Inside a permanent magnet magnetic scalar potential satisfies Poisson's equation ∆ϕ m = div M .

(7)

Magnetic field vector can be presented as H = − grad ϕ m .

(8)

Because outside a permanent magnet M = 0 , magnetic scalar potential, ϕ m 0 satisfies Laplace's equation, ∆ϕ m 0 = 0 , (9) where H 0 = − grad ϕ m 0 .

(10)

The third method that is mentioned in the paper for magnetic field calculation is based on superposition of results obtained for elementary magnetic dipoles.

where nˆ is unit vector of outgoing normal (Fig.1).

Fig.2 - Elementary magnetic dipole

Elementary magnetic dipole (Fig.2) has magnetic moment (11) d m = M dV ′ .

Fig.1 - Permanent magnet

These currents produce magnetic vector potential: µ A(r ) = 0 4π

∫ J (r ′) a

dV ′ µ0 + R 4π

V

where R = r − r ′ .

∫J S

sa ( r ′)

d S′ , R

(3)

This magnetic moment produces, at field point P , elementary magnetic scalar potential 1 Rdm 1 RM (12) d ϕm = = dV ′ , 3 4π R 4π R 3

Magnetic flux density is 1

University of Nis, Faculty of Electronic Engineering, Aleksandra Medvedeva 14, 18000 Nis, Serbia&Montenegro, E-mail: [email protected]

where R = r − r ′ is distance from the point where the magnetic field is being calculated to elementary source, and R = r − r ′ . After integration magnetic scalar potential is obtained as 1 R d m 1 RM ϕm = = d V ′ . (13) 4π 4π R 3 R3



V



V

ϕ m1 = −

M 4π

h 2 R2 l

∫ ∫∫ ((x − x')

h − R1 0 2

(x − x ′) d x ′ d y ′ d z ′ 32 2 + ( y − y ')2 + (z − z ')2 )

(19

) The second part of the magnetic circuit is homogeneously magnetized in positive direction of x-axis, ) M = Mx . (20)

PROBLEM DEFINITION

Scalar product is

Permanent magnet that is observed in the paper is homogeneously magnetized along its axes. This is magnetic circuit that is made of ferromagnetic material and it consists of four parts. Each of these parts is magnetized in direction shown in the Fig.3. Dimensions of the permanent magnet are also presented in the Fig.3.

Substituting expressions (15) and (21) in (13), magnetic scalar potential ϕ m 2 can be calculated

ϕm 2 =

M 4π

RM = M ( x − x ′) .

h 2 − R1 l

∫ ∫ ∫ ((x − x')

(21)

(x − x ′) d x ′ d y ′ d z ′ 32 + ( y − y ' )2 + ( z − z ' )2 )

(22)

2

h − − R2 0 2

The third part of the magnetic circuit is homogeneously magnetized in angular direction M = M θˆ .

(23)

Relations between coordinates x , y , z and cylindrical coordinates r , θ and z are x = r cosθ , y = r sinθ and z = z .

(24)

Using these relations, distance from the point where the magnetic field is being calculated to elementary source, can be presented as R = r 2 + r ' 2 −2rr ' cos(θ − θ' ) + ( z − z ' ) 2 .

Fig.3 – Permanent magnet

Outside the permanent magnet magnetic scalar potential, at field point P( x, y, z ) , can be determined using superposition of results obtained for each magnetized part, (14) ϕ m = ϕ m1 + ϕ m 2 + ϕ m 3 + ϕ m 4 , where ϕ m1 and ϕ m 2 are magnetic scalar potentials of

As magnetization has only angular component θ , scalar product RM is obtained as RM = [(x − x')xˆ + ( y − y ') yˆ + (z − z ')zˆ ]M θˆ .

scalar potentials of fourth part. These magnetic scalar potentials can be determined using the expression (13), where R = r − r′ =

( x − x ') 2 + ( y − y ' )2 + ( z − z ' )2

,

R = r − r ′ = (x − x')xˆ + ( y − y ') yˆ + (z − z ')zˆ .

(26)

Relation between unit vectors xˆ , yˆ , zˆ and θˆ is

two ends that are magnetized in different direction, ϕ m 3 is magnetic scalar potentials produced by the part which is magnetized in angular direction and ϕ m 4 is magnetic

(25)

∂y 1 ⎛ ∂x ∂z xˆ + yˆ + θˆ = ⎜ r ⎝∂θ ∂θ ∂θ

⎞ zˆ ⎟ , ⎠

(27)

and the following expressions are also satisfied

(15)

1 ∂x ˆ 1 ∂z 1 ∂y ˆ θˆ xˆ = , θ yˆ = , θ zˆ = . (28) r ∂θ r ∂θ r ∂θ Using relations (24) and (28), the following relations are obtained (29) θˆ xˆ = − sinθ , θˆ yˆ = cosθ , θˆ zˆ = 0 ,

(16)

and scalar product can be presented as

and

The first part of the magnetic circuit is homogeneously magnetized in negative direction of x-axis, ) (17) M = − Mx . RM = − M ( x − x ′) . Scalar product is (18) Substituting expressions (15) and (18) in (13), magnetic scalar potential ϕ m1 is obtained

RM = Mr ' sin( θ − θ ' ) .

(30)

Substituting expressions (25) and (30) in (13), magnetic scalar potential is calculated as

ϕm3 =

M 4π

h R2 2 π

∫ ∫ ∫ [r

R1 h 0 2

r ' sin(θ − θ' ) d r ' d z ' d θ' 2

+ r ' 2 −2rr ' cos (θ − θ') + (z − z ')2

]

32

(31) The fourth part of the permanent magnet is homogeneously magnetized in positive direction of y-axis,

Scalar product is

M = Myˆ .

(32)

RM = M ( y − y ′) .

(33)

Substituting expressions (15) and (33) in (13), magnetic scalar potential ϕ m 4 is obtained ϕm 4 =

M 4π

h 2 R1 l +l0 + a

( y − y ′) d x ′ d y ′ d z ′ 32 2 + ( y − y ') 2 + ( z − z ' ) 2 )

∫ ∫ ∫ ((x − x')



h − R1 l +l0 2

(34) The solutions of integrals presented in the expressions (19), (22), (31) and (34) are very complex. Because of that the expression for magnetic scalar potential (14) is very large and it won’t be shown in the paper, but it is used for determining the components of magnetic field. Magnetic field vector can be expressed as H = − grad ϕ m ,

(35)

therefore its components are Hx = −

∂ϕ ∂ϕ m ∂ϕ , H y = − m , H z = − m . (36) ∂x ∂z ∂y

NUMERICAL RESULTS

Distribution of magnetic flux density outside the permanent magnet is presented in the Fig.4.

Magnetic flux density is obtained using the analytical method for the following dimension of permanent magnet R1 h = 2, R 2 h = 4 , a h = 2, l h = 6 and l 0 h = 0.2 . The Fig.5 presents distribution of magnetic flux density (arrow) and magnetic potential (contour) obtained using FEMLAB software. Comparing these figures, the excellent agreement between analytical method results and FEMLAB results is evident. Table I Normalized magnetic field values along the direction

y = 0, z = 0

x/l 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0

H /M 1.223270 0.515419 0.267753 0.150165 0.087647 0.051969 0.030798 0.019224 0.016823 0.023926

Table II Normalized magnetic field values along the direction

x l = 1.0167, z = 0 y/l 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3

H /M

0.028094 0.054448 0.257539 0.354298 0.195708 0.184541 0.182569 0.042942 0.015348 0.007758 0.006290 0.006057 0.005785

Fig.4 – Distribution of magnetic flux density

In the Table I and Table II, magnetic field values along characteristic directions, for mentioned dimensions of permanent magnet, are presented.

Fig.5 – Distribution of magnetic flux density (FEMLAB)

Fig.6 – Distribution of magnetic flux density

When the length of the part that is magnetized along y-axis, is equal to R1 + R 2 , distribution of magnetic flux density outside the permanent magnet is presented in the Fig.6. It is obtained using the analytical method, for the following dimension of permanent magnet R1 h = 4, R 2 h = 6 , a h = 2, l h = 6 and l 0 h = 0.2 . The Fig.7 presents distribution of magnetic flux density (arrow) and magnetic potential (contour) obtained using FEMLAB software.

y/l 0.9 1.0 1.1 1.2 1.3

H /M

0.015348 0.007758 0.006290 0.006057 0.005785

CONCLUSION

Permanent magnet, homogeneously magnetized along its axes is observed in the paper. Method that is used for magnetic field determination is based on superposition of results that are obtained for elementary magnetic dipoles. The tables with magnetic field values, in different points, in vicinity of permanent magnet, are shown. Magnetic flux density distribution of permanent magnet is also presented in the paper. Magnetic field lines have the same form and the same direction as magnetic flux density lines, outside the magnet. Results obtained by analytical method are satisfactory confirmed using FEMLAB software. REFERENCES

Fig.7 – Distribution of magnetic flux density (FEMLAB)

Comparing Fig.6 and Fig.7 the results obtained by analytical method are satisfactory confirmed using FEMLAB software. In the Table III and Table IV, magnetic field values along characteristic direction, for mentioned dimensions of permanent magnet, are presented. Table III Normalized magnetic field values along the direction

y = 0, z = 0

x/l 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0

H /M 0.438976 0.210030 0.128901 0.085889 0.059395 0.041883 0.029927 0.021712 0.016179 0.012612

Table IV Normalized magnetic field values along the direction

x l = 1.0167, z = 0 y/l 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8

H /M 0.028094 0.054448 0.257539 0.354298 0.195708 0.184541 0.182569 0.042942 (continue)

[1] Textbook of Electromagnetics, Edited by Hermann Uhlmann (Veličković D. M., Chapter 3.2: Magnetostatic field), Technische Universitaet Ilmenau. [2] A. N. Mladenović, S. R. Aleksić: “Rod shaped permanent magnet”, 48. Conference ETRAN 2004, Čačak, 6th-11th Jun, Proceedings of full papers, pp. 233-236. (In Serbian) [3] A. N. Mladenović, Slavoljub R. Aleksić: Toroidal shaped permanent magnet with air gap, International Conference on Applied Electromagnetics PES2005, Niš 23rd-25th May 2005, Proceedings of Extended Abstracts, pp. 23-24. [4] Ana N. Mladenović, Slavoljub R. Aleksić:” Methods for magnetic field calculation”, 11th International Conference on Electrical Machines, Drives and Power Systems ELMA 2005, Sofia, Bulgaria, 15-16 September 2005, Vol.2, pp. 350-354. [5] Ana N. Mladenović: “Toroidal shaped permanent magnet with two air gaps”, International PhD-Seminar “Numerical Field Computation and Optimization in Electrical Engineering”, Proceedings of Full Papers, Ohrid, Macedonia, 20-25 September 2005, pp.153-158. Ana N. Mladenović was born in Niš, Serbia on July 29, 1977. In 1996 she enrolled in the Faculty of Electronic Engineering, University of Niš, Serbia and Montenegro, with the Computer Science as major. She graduated in 2003 with the final thesis in Electromagnetic. She enrolled in the postgraduate studies at the same faculty at the Department of Theoretical Electrical Engineering, in 2004.

Since 2003 she has been working as a teaching assistant at the Department of Theoretical Electrical Engineering, at the Faculty of Electronic Engineering of the University in Niš. Her main research field are permanent magnets and she is author or co-author of 10 papers presented at national and international conferences.

instructions for preparation of extended abstract - PhD Seminar 2006

May 25, 2005 - [5] Ana N. Mladenović: “Toroidal shaped permanent magnet with two air gaps”, International PhD-Seminar “Numerical. Field Computation and Optimization in Electrical Engineer- ing”, Proceedings of Full Papers, Ohrid, Macedonia, 20-25. September 2005, pp.153-158. Ana N. Mladenović was born in Niš,.

363KB Sizes 0 Downloads 344 Views

Recommend Documents

Preparation of Extended Abstract for NS2008 Conference Proceeding
[3] I. Streeter, G.G. Wildgoose, L. Shao, R.G. Compton,. Sens. Actuators B 133 (2008) 462. [4] W. Zhao, J.J. Xu, Q.Q. Qiu, H.Y. Chen, Biosens. Bioelectron. 22 (2006) 649. [5] N. Gibson, O. Shenderova, T.J.M. Luo, S. Moseenkov, V. Bondar, A. Puzyr, K.

instructions for preparation of manuscripts -
The authors wish to express their acknowledgement to the Basque Government and the Provincial Council of. Gipuzkoa. This work has been funded by them as ...

instructions for preparation of papers
sources (RESs) could be part of the solution [1,2,3]. A HPS is ... PV-HPSs are one of the solutions to the ..... PV Technology to Energy Solutions Conference and.

Extended abstract
'DEC Systems Research Center, 130 Lytton Av- enue, Palo-Alto ... assigned to any server, (we call such tasks un- .... (the optimum off-line algorithm) runs task.

Extended Abstract -
the early 1990s, Sony entered the market and secured a leading position due to ...... of Humanities and Social Sciences, Rose-Hulman Institute of Technology.

Extended Abstract
Keywords: Limit angular speed, Von Mises criterion, Annular Disk. ..... It is also quite obvious that for disks without attached masses, failure always occurs at the ...

ASEG Extended Abstract
out a perfectly focused image at the zero correlation lag. However, there are other classes of penalty functions that can be used in the ASM-IDT inversion procedure; e.g., ones that compensate for illumination irregularities (Yang et al., 2012) or mo

instructions for preparation of full papers
simplify the data structure and make it very direct to be accessed so that the ... enterprise and road companies should be used as much as possible to .... transportation data management and analysis and fully integrates GIS with .... Nicholas Koncz,

instructions for preparation of a camera ready manuscript
acoustic field performed in the frequency or the time domain. .... measurement ( s(t) ) in a noise free environment at a range of 5 km and depth of 100 m is.

Instructions for the Preparation of a
Last great milestone in this environment was made by introducing Internet and .... {Class/Individual2} are two dynamic text placeholders which are in process of ...

instructions to authors for the preparation of manuscripts
e-learning system developers is to build such systems that will create individualized instruction .... to answer the research question, the researcher can carry out ...

instructions to authors for the preparation of papers -
(4) Department of Computer Science, University of Venice, Castello 2737/b ... This paper provides an overview of the ARCADE-R2 experiment, which is a technology .... the German Aerospace Center (DLR) and the Swedish National Space ...

instructions for preparation of full papers
combination of bus, subway, and train routes is not an easy job even for local ..... transportation data management and analysis and fully integrates GIS with .... TransCAD Software (2003) TransCAD, Caliper Corporation, Newton MA, USA.

instructions to authors for the preparation of manuscripts
All these capabilities make airships an attractive solutions for the civilian and military ..... transmit video, remote control, and other data exchange tasks. Camera based ... Sensors, Sensor. Networks and Information Processing Conference, 2004.

Instructions for the Preparation of a
Faculty of Natural Sciences and Mathematics and Education, University of Split, ... to be direct application of the information and communication technology.

Work instructions for preparation of an orphan maintenance ...
May 30, 2018 - If the opinion was reached without need for LOQ, go to step 3. 2. After the AR has been sent to the EC delete LOQ and summary of sponsor's ...

Extended Abstract for ICFSMA 2011, Dresden
2 Photonic jet interpretation in wave optics scope. In the following discussion, we consider an unpolarised harmonic plane wave, with a wavelength in vacuum λ, which propagates along the z-axis in a free ..... [6] Lee, S.-H., et al., Characterizing

instructions for preparing usnccm vi final abstract
In this paper, we introduce a framework for probing the effect of resolving particle morphology using discrete element modeling. Specifically, the formulation uses ...

abstract instructions for authors - numiform 2007
Today, many scientists, engineers, companies, governamental and non-governamental agencies agree that hydrogen will be an important fuel in the future. A relevant application of hydrogen energy is related to the problem of air pollution caused by roa

Extended Abstract - Open Repository of National Natural Science ...
ations, whereas the benchmark in the competitive analysis is an optimal mechanism that knows the value distribution. We propose a mechanism Γ1 which obtains constant competitive ratios under various valuation models. We also conduct numerical simula

Instructions to Prepare the Abstract for the International ...
physics simulations”, Proc. Int. Workshop on Coupled Methods in Numerical. Dynamics, Faculty of Mechanical Engineering and Naval Architecture, University.

Instructions for COLING/ACL 2006 Proceedings
knowledge sources used in inferring the temporal ordering of events ... data sparseness, we used temporal reasoning as ..... The mining uses a set of lexical and.

Instructions for Preparing LREC 2006 Proceedings
bombardement ('bombing') instance. 2.1.2. Tokens and parts of speech. One of the most distinctive features of the French. Treebank with relation to the Penn Treebank is its mixed tokenization: words are tokenized as simple and compound words. This si