NONLINEAR FINITE ELEMENT METHOD IN MAGNETISM
|
|
- Gordon Carson
- 5 years ago
- Views:
Transcription
1 POLLACK PERIODICA An International Journal for Engineering and Information Sciences DOI: /Pollack Vol. 4, No. 2, pp (2009) NONLINEAR FINITE ELEMENT METHOD IN MAGNETISM Miklos KUCZMANN Laboratory of Electromagnetic Fields, Department of Telecommunication Széchenyi István University, Egyetem tér 1, H-9026 Győr, Hungary Received 9 February 2009; accepted 26 May 2009 Abstract: The paper presents and compares three potential formulations to solve nonlinear static magnetic field problems by applying the fixed-point technique and the Newton-Raphson scheme. Nonlinear characteristics have been handled by the polarization method in the two algorithms. The proposed combination of Newton-Raphson scheme and the polarization formulation result in a very effective nonlinear solver, because only the derivate of the characteristics, i.e. only the permeability or the reluctivity has to be used. That is why, this method can be prosperous to solve nonlinear problems with hysteresis, and it is faster than the classical fixed-point method. Keywords: Finite element method, Fixed-point method, Newton-Raphson method, Nonlinear magnetics 1. Introduction Nonlinear electromagnetic field problems can only be solved by iterative algorithms. There are two widely used techniques, the fixed-point method, and the Newton- Raphson scheme. The first one is known as a stable, but very slow algorithm, while the second one has quadratic convergence speed close to the solution of nonlinear equations, however Newton-Raphson scheme is sometimes divergent [1]-[7]. The nonlinear characteristics of ferromagnetic materials can be handled by the polarization formulation [1]-[5]. This formulation is used to prescribe the nonlinear constitutive relations between the magnetic field intensity H and the magnetic flux density B as HU ISSN Akadémiai Kiadó, Budapest
2 14 M. KUCZMANN and B = µ H + R, (1) H = ν B + I, (2) where R and I are the nonlinear residual terms determined iteratively, moreover µ and ν are permeability-like and reluctivity-like quantities. In other words, the nonlinear output of the hysteresis models can be split into a linear part and a nonlinear part. In the case of fixed-point method, µ and ν are usually constant [1]-[5], and µ max + µ µ = min, (3) 2 ν max + ν ν = min. (4) 2 Here µ max, µ min, ν max, and ν min are the maximum and the minimum slope of the nonlinear characteristics. This selection results in global convergence of the fixed-point technique. The other possibility for representation of nonlinear behavior is the application of nonlinear permeability or nonlinear reluctivity [6]-[8], i.e. B = [ µ ( H )]H or H = [ ν ( B) ]B. The main disadvantage of this method is that the derivatives of the elements of tensors [ µ () ] or [ ν () ] have to be determined, and these functions must be monotonous. In general, [ µ () ] and [ ν ( ) ] is not monotonous. This is the reason why, in the case of hysteresis, the classical Newton-Raphson technique can be difficult to use. This paper presents a combination of the polarization formulation (1) or (2) and the Newton-Raphson method, which results in an algorithm as simple as the fixed-point technique, and only the characteristics have to be monotonous, neither the permeability nor the reluctivity. 2. Governing nonlinear equations Let suppose the domain of interest = 0 m, ( 0 and m are filled with air and with magnetic material, respectively) is surrounded by the boundary Γ, which can be decomposed into two parts (refer to Fig. 1), i.e. Γ = Γ H ΓB, where boundary conditions are prescribed [9], [10]. Nonlinear static magnetic field problems can be described by the following Maxwell s equations [9], [10]
3 NONLINEAR FINITE ELEMENT METHOD IN MAGNETISM 15 H = J, B = 0, in, (5) where H, B and J are the magnetic field intensity, the magnetic flux density and the source current density. Source current density is placed only in the air region. The following boundary conditions are prescribed on the boundary [9], [10] H n = K, B n = b, on on ΓH, ΓB, (6) where K and b are the surface current density, and the charge density of fictitious magnetic surface charges, moreover n is the outer normal unit vector of the domain. The constitutive relation between the magnetic field intensity and the magnetic flux density is linear in air, B = µ 0H or H = ν 0B, where µ 0 and ν 0 = 1 µ 0 are the permeability and the reluctivity of vacuum. In the nonlinear region filled with ferromagnetic material, the constitutive laws in the form of (1) or (2) can be used. Fig. 1. Structure of a static magnetic field problem 2.1. The magnetic vector potential by nodal FEM The magnetic vector potential A can be introduced from the second equation in (5), i.e. [9], [10] B = A. (7) Substituting it to (2) and (5), and enforcing Coulomb gauge implicitly, result in the following partial differential equations [9], [10] ( ν 0 A ) ( ν 0 A) = J, in 0, (8)
4 16 M. KUCZMANN ( ν A ) ( ν A) = I, in m. (9) The following boundary conditions are appended to these equations [9], [10], ( ν A + I ) ( ν A) n = K, 0 n = K, on Γ H, (10) A n = 0, on Γ H, (11) n A = α, on Γ B, (12) ν A = 0, ν0 A = 0, on Γ B. (13) Here the term α is a constant representing b in (6) [9]. In this case, the magnetic vector potential is approximated by nodal shape functions The magnetic vector potential by edge FEM The other possibility is using the ungauged formulation of magnetic vector potential, which is approximated by the edge shape functions [9], [10], [11], i.e. the unknowns are associated to the edges of the mesh. The source current density must be approximated by the impressed current vector potential T as T = J, in, (14) since J = 0, i.e. the following partial differential equations are valid [9], [10], [11], ( ν 0 A ) = T, in 0, (15) ( ν A ) = T I, in m, (16) with the boundary conditions (10) and (12), and K = 0 according to the first equation in (23) The reduced magnetic scalar potential The source current density is approximated by using the impressed current vector potential T as in (14), and from the first Maxwell s equation in (5) [9], [10], [11] H = T Φ (17)
5 NONLINEAR FINITE ELEMENT METHOD IN MAGNETISM 17 can be obtained, where Φ is the so called reduced magnetic scalar potential, and it is approximated by nodal FEM. Substituting it back to (1) and to (5) results in the partial differential equations, ( µ 0 Φ ) = ( µ 0T ), in 0, (18) ( µ Φ ) = ( µ T ) R, in m, (19) and in the following boundary conditions are valid Φ = 0, on Γ H, (20) [ µ ( T Φ ) + R] ( T Φ ) n = b, µ 0 n = b, on ΓB. (21) 2.4. The impressed current vector potential The impressed current vector potential can be obtained by the solution of the partial differential equation [10] T = J, in, (22) with the following boundary conditions T n = K, T n = 0, on on ΓH, ΓB. (23) This ungauged formulation can be solved by the usual finite element procedures (see in [10]). 3. The Newton-Raphson method in the polarization formulation The magnetic vector potential as well as the reduced magnetic scalar potential has to be updated by A and by Φ in every nonlinear iteration steps [5], [6], [7]. The weak form of the potential formulations can be built up by using the weighted residual method and the Galerkin scheme [9], [10]. The nonlinear terms I and R are linearized as and I ( A + A) I( A) + ([ ] [ 1] ν ) ( A) ν d, (24) ( Φ Φ ) R( Φ ) ([ µ ] [ 1] µ ) ( Φ ) R + d, (25)
6 18 M. KUCZMANN where [ ν d ] and [ µ d ] are the differential reluctivity tensor and the differential permeability tensor of the nonlinear characteristics, which are updated in every nonlinear steps. Here µ and ν are selected as in (1) and in (2). Applying (24) and (25) in the weighted residual method, the weak form of the potential formulations can be obtained The magnetic vector potential by nodal FEM After applying the weighted residual method to the partial differential equations and boundary conditions presented in 2.1, the following weak formulation can be obtained, ν 0 W ( A + A) + ν W A 0 + W [ ν d ] m ( A) + ν W ( A + A) + W I m W J W K dγ = 0. Γ 0 H + ν 0 W ( A + A) m m 0 (26) Here W is the vector weighting function [9], [10]. The magnetic vector potential is updated by A, which is the solution of (26) in every iteration steps. After applying the weighted residual method to the partial differential equations and boundary conditions presented in 2.2, a weak formulation similar to (23) can be obtained The magnetic vector potential by edge FEM After applying the weighted residual method to the partial differential equations and boundary conditions presented in 2.2, the following weak formulation can be obtained ν 0 W ( A + A) + ν W A 0 m + W [ ν d ] m ( A) + W I W T = 0, m 0 m (27) where W is the vector weighting function. The magnetic vector potential is updated by A, i.e. the solution of (27) The reduced magnetic scalar potential After applying the weighted residual method to the partial differential equations and boundary conditions presented in 2.3, the following weak formulation can be obtained,
7 NONLINEAR FINITE ELEMENT METHOD IN MAGNETISM 19 µ 0 N ( Φ + Φ ) + µ N Φ + N [ µ d ] 0 m m N R µ N T Nb dγ = 0. m 0 m ΓB ( Φ ) (28) Here N is the scalar weighting function [9], [10]. The magnetic scalar potential is updated by Φ, which is the solution of (28) in every iteration steps. 4. Simulation results Three problems (a modified version of TEAM 10, TEAM 13, and a modified version of TEAM 24 defined by the COMPUMAG Society) have been solved by the above-introduced Newton-Raphson algorithm and by the classical fixed-point technique. The same nonlinear characteristics presented in Fig. 2 have been used, and the curve has been approached by a simple piecewise linear approximation. Originally, this nonlinearity is measured using the material from TEAM Modified version of TEAM problem No. 10 The problem can be seen in Fig. 3 [12]. Steel plates have been placed around a racetrack shaped coil. Only the eighth, of the problem can be analyzed, because of symmetry. Fig. 2. The nonlinear characteristics of the simulated problems Fig. 3. The arrangement of TEAM problem No. 10 Two types of mesh have been analyzed. The first one consists of 8413 tetrahedra, the second one has been built up by finite elements. The first mesh results in 12850, 57278, and unknowns for Φ, vector A, and nodal A formulation, respectively. The second discretization results in 63807, , and unknowns for the same formulation.
8 20 M. KUCZMANN Fig. 4 shows the finer mesh of the problem. The magnetic flux density is driven by the steel plates around the coil as it can be seen in Fig. 5. Fig. 4. The FEM mesh of TEAM 10 Fig. 5. Magnetic flux is driven by the plates Fig. 6 shows the magnetic flux density along a line [12] placed inside the plates (the source current density is given in the figures) and along a line just below the horizontal plate. The three formulations give practically the same results. It is noted that, the nodal A formulation is more sensitive to the density of the finite element mesh. Fig. 6. Distribution of magnetic flux density inside the plates and in air under the horizontal plate The Newton-Raphson method is much faster than the classical fixed-point method as it can be seen in Fig. 7. The measured magnetic flux density in the center of the central plate is 1.67 T [12]. The simulated data are T, T, T by the reduced magnetic scalar potential formulation, by the edge element represented magnetic vector potential and by
9 NONLINEAR FINITE ELEMENT METHOD IN MAGNETISM 21 the magnetic vector potential approximated by nodal finite elements. It can be concluded that the edge element based A formulation gives the best result, moreover B simulated by the reduced magnetic scalar potential is a little larger than the results obtained from the magnetic vector potential, moreover the nodal vector potential formulation is more sensitive to the density of the mesh 4.2. TEAM problem No. 13 Fig. 7. Comparison of the number of iterations by Newton-Raphson method and the fixed-point method This problem is a modified version of TEAM 10 [13], [14]. The U-shaped yokes have been translated as it can be seen in Fig. 8. Only the half of the geometry has been analyzed. Fig. 8. The arrangement of TEAM problem No. 13 The coarse mesh consists of tetrahedra, which results in 40180, , and unknowns for Φ, edge element based A, and nodal A, respectively. The
10 22 M. KUCZMANN dense mesh contains finite elements. It means , , and unknowns for the same potentials. Fig. 9 shows some comparisons between the results simulated by the three, presented potential formulations. The magnetic flux density simulated by the reduced magnetic scalar potential is a little larger than the results obtained from the magnetic vector potential formulations, moreover the nodal vector potential formulation is more sensitive to the density of the mesh, as it was experienced in the last example, too. Comparison between measured and simulated data shows that the vector A formulation is the closest to the measured data. Fig. 9. Distribution of magnetic flux density inside the plates 4.3. Modified version of TEAM problem No. 24 This is a modified version of the problem TEAM 24 (Fig. 10) [15]. The problem contains a motor which rotor has been locked. The source current of the coils is constant, i.e. a static magnetic field problem has been analyzed, and the characteristics of the stator and the rotor are given in Fig. 2. Fig. 10. Model of TEAM problem 24
11 NONLINEAR FINITE ELEMENT METHOD IN MAGNETISM 23 The number of finite elements is 39060, which results in and unknowns for the Φ -formulation and for the edge element based A -formulation, respectively. The magnetic flux is mainly driven by the stator and the rotor steels as it can be seen in Fig. 11, and a comparison between the simulated results is shown in Fig. 12. Fig. 11. The magnetic flux density inside the motor Fig. 12. Magnetic flux density along the path shown in Fig Conclusion It can be concluded that the mentioned Newton-Raphson method combined with the polarization formulation is much faster than the fixed-point technique, however the weak formulations and the applied algorithms are almost the same. Only the permeability or the reluctivity of the material has to be simulated.
12 24 M. KUCZMANN The aim of further research is to find the way to solve nonlinear eddy current field problems, and taking the nonlinearity into account by the vector Preisach model of hysteresis. Acknowledgements This paper was supported by the János Bolyai Research Scholarship of the Hungarian Academy of Sciences (BO/00064/06), by Széchenyi István University ( ), by the Hungarian Scientific Research Fund (OTKA PD 73242), and by Hungarian Science and Technology Foundation (OMFB-00725/2008). References [1] Dlala E., Belachen A., Arkkio A. A fast fixed-point method for solving magnetic field problems in media with hysteresis, IEEE Trans. on Magn, Vol. 44, 2008, pp [2] Dlala E., Arkkio A. Analysis of the convergence of the fixed-point method used for solving nonlinear rotational magnetic field problems, IEEE Trans. on Magn, Vol. 44, 2008, pp [3] Auserhofer S., Bíró O., Preis K. A strategy to improve the convergence of the fixed-point method for nonlinear eddy current problems, IEEE Trans. on Magn, Vol. 44, 2008, pp [4] Saitz J. Newton-Raphson method and fixed-point technique in finite element computation of magnetic field problems in media with hysteresis, IEEE Trans. on Magn, Vol. 35, 1999, pp [5] Chiampi M., Repetto M., Chiarabaglio D. An improved technique for nonlinear magnetic problems, IEEE Trans. on Magn, Vol. 30, 1994, pp [6] O Dwyer J., O Donnell T. Choosing the relaxation factor for the solution of nonlinear magnetic field problems by the Newton-Raphson method, IEEE Trans. Magn, Vol. 31, 1995, pp [7] Koh C. S., Ryu J. S., Fujiwara K. Convergence acceleration of the Newton-Raphson method using successive quadratic function approximation of residual, IEEE Trans. Magn, Vol. 42, 2006, pp [8] Fonteyn K., Belahcen A., Arkkio A. Properties of electrical steel sheets under strong mechanical stress, Pollack Periodica, Vol. 1, 2006, pp [9] Bíró O. CAD in electromagnetism, Advances in electronics and electron physics, Vol. 82, 1991, pp [10] Kuczmann M., Iványi A. The finite element method in magnetism. Akadémiai Kiadó, Budapest, [11] Bíró O. Edge element formulations of eddy current problems, Computer Methods in Applied Mechanics and Engineering, Vol. 160, 1999 pp [12] Preis K., Bárdi I., Bíró O., Magele C., Renhart W., Richter K.R., Vrisk G. Numerical analysis of 3D magnetostatic fields, IEEE Trans. Magn, Vol. 27, 1991, pp [13] Nakata T., Fujiwara K. Summary of results for benchmark problem 13 (3-D nonlinear magnetostatic model), COMPEL, Vol , pp [14] Preis K., Bárdi I., Bíró O., Magele C., Vrisk G., Richter K. R. Different finite element formulations of 3D magnetostatic fields, IEEE Trans. Magn, Vol. 28, 1992, pp [15] Allen N., Rodger D. Description of TEAM workshop problem 24, nonlinear time-transient rotational test rig,
Handling Nonlinearity by the Polarization Method and the Newton-Raphson Technique
Handling Nonlinearity by the Polarization Method and the Newton-Raphson Technique M. Kuczmann Laboratory of Electromagnetic Fields, Department of Telecommunications, Széchenyi István University, Egyetem
More informationNewton-Raphson Method in Nonlinear Magnetics
Newton-Raphson Method in Nonlinear Magnetics M. Kuczmann Széchenyi István University, Department of Telecommunications, Laboratory of Electromagnetic Fields H-9026, Egyetem tér 1, Gyır, Hungary Abstract:
More informationFinite Element Analysis of Single-Phase Induction Motors
Finite Element Analysis of Single-Phase Induction Motors Dániel Marcsa, Miklós Kuczmann Széchenyi István University, Győr, Egyetem tér 1, Hungary, H-9026, email: marcsadaniel@yahoo.co.uk Abstract: The
More informationLABORATORY OF ELECTROMAGNTIC FIELDS SZÉCHENYI ISTVÁN UNIVERSITY FINITE ELEMENT ANALYSIS
FINITE ELEMENT ANALYSIS OF SINGLE-PHASE INDUCTION MOTOR D. MARCSA, M. KUCZMANN Laboratory of Electromagnetic Fields, Department of Telecommunications Széchenyi István University, H-9026, Egyetem tér 1,
More informationThe Newton-Raphson method accelerated by using a line search - comparison between energy functional and residual minimization
Physics Electricity & Magnetism fields Okayama University Year 2004 The Newton-Raphson method accelerated by using a line search - comparison between energy functional and residual minimization Koji Fujiwara
More informationEddy Current Losses in the Tank Wall of Power Transformers
Eddy Current Losses in the Tank Wall of Power Transformers Erich Schmidt Institute of Electrical Drives and Machines, Vienna University of Technology A 14 Vienna, Austria, Gusshausstrasse 25 29 Phone:
More informationSIMULATION OF A PMS MOTOR BY THE HELP OF TWO DIFFERENT DESIGN SOFTWARE TOOLS
Bulletin of the Transilvania Universit of Braşov Series I: Engineering Sciences Vol. 6 (55) No. 1-213 SIMULATION OF A PMS MOTOR BY THE HELP OF TWO DIFFERENT DESIGN SOFTWARE TOOLS G. KOVÁCS 1 M. KUCZMANN
More informationNumerical problems in 3D magnetostatic FEM analysis
Numerical problems in 3D magnetostatic FEM analysis CAZACU DUMITRU Department of electronics, computers and electrical engineering University of Pitesti Str.Tirgu din Vale nr.1 Pitesti ROMANIA cazacu_dumitru@yahoo.com
More informationInduction Motors Simulation by Finite Element Method and Different Potential Formulations with Motion Voltage Term
Induction Motors Simulation by Finite Element Method and Different Potential Formulations with Motion Voltage Term by Dániel Marcsa B.Sc. student in Electrical Engineering supervisor Dr. Miklós Kuczmann,
More informationMethod for determining relaxation factor for modified Newton-Raphson method
Physics Electricity & Magnetism fields Okayama University Year 1992 Method for determining relaxation factor for modified Newton-Raphson method Koji Fujiwara Takayoshi Nakata N. Okamoto Kazuhiro Muramatsu
More informationEddy-Current Effects in Circuit Breakers During Arc Displacement Phase
Eddy-Current Effects in Circuit Breakers During Arc Displacement Phase Olivier Chadebec, Gerard Meunier, V. Mazauric, Yann Le Floch, Patrice Labie To cite this version: Olivier Chadebec, Gerard Meunier,
More informationCharacteristics Analysis of the Square Laminated Core under dc-biased Magnetization by the Fixed-point Harmonicbalanced
Characteristics Analysis of the Square Laminated Core under dc-biased Magnetization by the Fixed-point Harmonicbalanced FEM Author Zhao, Xiaojun, Li, Lin, Lu, Junwei, Cheng, Zhiguang, Lu, Tiebing Published
More informationThe Linear Induction Motor, a Useful Model for examining Finite Element Methods on General Induction Machines
The Linear Induction Motor, a Useful Model for examining Finite Element Methods on General Induction Machines Herbert De Gersem, Bruno Renier, Kay Hameyer and Ronnie Belmans Katholieke Universiteit Leuven
More information338 Applied Electromagnetic Engineering for Magnetic, Superconducting, Multifunctional and Nano Materials
Materials Science Forum Online: 2014-08-11 ISSN: 1662-9752, Vol. 792, pp 337-342 doi:10.4028/www.scientific.net/msf.792.337 2014 Trans Tech Publications, Switzerland Torque Characteristic Analysis of an
More informationKeywords: Electric Machines, Rotating Machinery, Stator faults, Fault tolerant control, Field Weakening, Anisotropy, Dual rotor, 3D modeling
Analysis of Electromagnetic Behavior of Permanent Magnetized Electrical Machines in Fault Modes M. U. Hassan 1, R. Nilssen 1, A. Røkke 2 1. Department of Electrical Power Engineering, Norwegian University
More informationNonlinear analysis of eddy current and hysteresis losses of 3-D stray field loss model (Problem 21)
Physics Electricity & Magnetism fields Okayama University Year 2000 Nonlinear analysis of eddy current and hysteresis losses of 3-D stray field loss model (Problem 21) Norio Takahashi Toshiomi Sakura Zhiguang
More informationFinite Element Modeling of Electromagnetic Systems
Finite Element Modeling of Electromagnetic Systems Mathematical and numerical tools Unit of Applied and Computational Electromagnetics (ACE) Dept. of Electrical Engineering - University of Liège - Belgium
More informationDedicating Finite Volume Method to Electromagnetic Plasma Modeling: Circuit Breaker Application
Dedicating Finite Volume Method to Electromagnetic Plasma Modeling: Circuit Breaker Application Loïc Rondot, Vincent Mazauric, Yves Delannoy, Gérard Meunier To cite this version: Loïc Rondot, Vincent Mazauric,
More informationUpdate On The Electromagnetism Module In LS-DYNA
12 th International LS-DYNA Users Conference Electromagnetic(1) Update On The Electromagnetism Module In LS-DYNA Pierre L'Eplattenier Iñaki Çaldichoury Livermore Software Technology Corporation 7374 Las
More informationAssessing COMSOL Performances on Typical Electromagnetic Problems Faced by Turbo-Generator Manufacturers
1 Assessing COMSOL Performances on Typical Electromagnetic Problems Faced by Turbo-Generator Manufacturers E. A. Badea *1, J. Oesterheld 1, S. Friedel 2 and M. Olsson 2 1 Alstom (Switzerland) LTD, 2 COMSOL
More informationTeam Problem D Static Force Problem
Team Problem 20 3-D Static Force Problem 1. General description The model is shown in Fig.1. The center pole and yoke are made of steel. The coil is excited by a dc current. The ampere-turns are 1000,
More informationDynamic simulation of a coaxial magnetic gear using global ODE's and DAE s and the rotating machinery, magnetic interface
Dynamic simulation of a coaxial magnetic gear using global ODE's and DAE s and the rotating machinery, magnetic interface M. Ostroushko 1, W. Zhang 1, W. M. Rucker 1 1. Institute of Theory of Electrical
More informationworks must be obtained from the IEE
Title Eddy-current analysis using and stop hysterons vector Author(s) Matsuo, T; Osaka, Y; Shimasaki, M Citation IEEE TRANSACTIONS ON MAGNETICS (200 1172-1177 Issue Date 2000-07 URL http://hdl.handle.net/2433/39989
More informationThe magnetic field diffusion equation including dynamic, hysteresis: A linear formulation of the problem
The magnetic field diffusion equation including dynamic, hysteresis: A linear formulation of the problem Marie-Ange Raulet, Benjamin Ducharne, Jean-Pierre Masson, G. Bayada To cite this version: Marie-Ange
More informationCONSIDER a simply connected magnetic body of permeability
IEEE TRANSACTIONS ON MAGNETICS, VOL. 50, NO. 1, JANUARY 2014 7000306 Scalar Potential Formulations for Magnetic Fields Produced by Arbitrary Electric Current Distributions in the Presence of Ferromagnetic
More informationTeam Problem D Non-Linear Magnetostatic Model
Team Problem 13 3-D Non-Linear Magnetostatic Model 1. General Description The model is shown in Fig.1. An exciting coil is set between two steel channels, and a steel plate is inserted between the channels.
More informationADAM PIŁAT Department of Automatics, AGH University of Science and Technology Al. Mickiewicza 30, Cracow, Poland
Int. J. Appl. Math. Comput. Sci., 2004, Vol. 14, No. 4, 497 501 FEMLAB SOFTWARE APPLIED TO ACTIVE MAGNETIC BEARING ANALYSIS ADAM PIŁAT Department of Automatics, AGH University of Science and Technology
More informationTheoretical Formulation of a Time-Domain Finite Element Method for Nonlinear Magnetic Problems in Three Dimensions
Progress In Electromagnetics Research, ol. 153, 33 55, 2015 Theoretical Formulation of a Time-Domain Finite Element Method for Nonlinear Magnetic Problems in Three Dimensions Su Yan 1, 2, * and Jian-Ming
More informationKent R. Davey American Electromechanics 2275 Turnbull Bay Rd. New Smyrna Beach, FL
1 Analytic Analysis of Single and Three Phase Induction Motors Kent R. Davey American Electromechanics 2275 Turnbull Bay Rd. New Smyrna Beach, FL 32168-5941 Abstract- The analysis of single and multiphase
More informationDevelopment and analysis of radial force waves in electrical rotating machines
DOI: 10.24352/UB.OVGU-2017-098 TECHNISCHE MECHANIK, 37, 2-5, (2017), 218 225 submitted: June 20, 2017 Development and analysis of radial force waves in electrical rotating machines S. Haas, K. Ellermann
More informationStudy on problems in detecting plural cracks by alternating flux leakage testing using 3D nonlinear eddy current analysis
Physics Electricity & Magnetism fields Okayama University Year 2003 Study on problems in detecting plural cracks by alternating flux leakage testing using 3D nonlinear eddy current analysis Yuji Gotoh
More informationExamination and Development of RF inductors by using the Finite Element Method
Examination and Development of RF inductors by using the Finite Element Method by Zoltán Pólik B.Sc. student in Electrical Engineering supervisor Dr. Miklós Kuczmann, Ph.D. Associate Professor A thesis
More informationGlobal formulation of 3D magnetostatics using flux and gauged potentials
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2000; 00:1 6 [Version: 2000/01/19 v2.0] Global formulation of 3D magnetostatics using flux and gauged potentials M.
More informationPulse eddy currents using an integral-fem formulation for cracks detection
International Journal of Applied Electromagnetics and Mechanics 33 (2010) 1225 1229 1225 DOI 10.3233/JAE-2010-1242 IOS Press Pulse eddy currents using an integral-fem formulation for cracks detection Gabriel
More informationInvestigation of evolution strategy and optimization of induction heating model
Physics Electricity & Magnetism fields Okayama University Year 1999 Investigation of evolution strategy and optimization of induction heating model Makoto Horii Norio Takahashi Takashi Narita Okayama University
More informationANSYS Based FEM Analysis for Three and Four Coil Active Magnetic Bearing-a Comparative Study
International Journal of Applied Science and Engineering 2013. 11, 3: 277-292 ANSYS Based FEM Analysis for Three and Four Coil Active Magnetic Bearing-a Comparative Study Pabitra Kumar Biswas a,* and Subrata
More informationFINITE ELEMENTS, ELECTROMAGNETICS AND DESIGN
FINITE ELEMENTS, ELECTROMAGNETICS AND DESIGN Edited by S. RATNAJEEVAN H. HOOLE Professor of Engineering Harvey Mudd College Claremont, CA, USA and Senior Fellow, 1993/94 Department of Electrical Engineering
More informationInfluence of magnetic anisotropy on flux density changes in dynamo steel sheets
ARCHIVES OF ELECTRICAL ENGINEERING VOL. 64(1), pp. 81-88 (2015) DOI 10.1515/aee-2015-0008 Influence of magnetic anisotropy on flux density changes in dynamo steel sheets WITOLD MAZGAJ, ZBIGNIEW SZULAR,
More informationHomogenization of the Eddy Current Problem in 2D
Homogenization of the Eddy Current Problem in 2D K. Hollaus, and J. Schöberl Institute for Analysis and Scientific Computing, Wiedner Hauptstr. 8, A-14 Vienna, Austria E-mail: karl.hollaus@tuwien.ac.at
More informationNeural Network Based Vector Hysteresis Model and the Nondestructive Testing Method
Neural Network Based Vector Hysteresis Model and the Nondestructive Testing Method by Miklós Kuczmann MSc in Electrical Engineering supervisor prof. Amália Iványi A thesis submitted to the Budapest University
More informationA 2-Dimensional Finite-Element Method for Transient Magnetic Field Computation Taking Into Account Parasitic Capacitive Effects W. N. Fu and S. L.
This article has been accepted for inclusion in a future issue of this journal Content is final as presented, with the exception of pagination IEEE TRANSACTIONS ON APPLIED SUPERCONDUCTIVITY 1 A 2-Dimensional
More informationGlobal formulation for 3D magneto-static using flux and gauged potential approaches
Global formulation for 3D magneto-static using flux and gauged potential approaches Maurizio Repetto, Francesco Trevisan 1 Dep. of Ingegneria Elettrica, Industriale, Politecnico di Torino C. Duca degli
More informationInverse Problems and Optimal Design in Electricity and Magnetism
Inverse Problems and Optimal Design in Electricity and Magnetism P. Neittaanmäki Department of Mathematics, University of Jyväskylä M. Rudnicki Institute of Electrical Engineering, Warsaw and A. Savini
More information2019 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media,
9 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising
More informationfield using second order edge elements in 3D
The current issue and full text archive of this journal is available at http://www.emerald-library.com using second order edge elements in 3D Z. Ren Laboratoire de GeÂnie Electrique de Paris, UniversiteÂs
More informationElectromagnetic Analysis Applied to the Prediction of Stray Losses in Power Transformer
Electromagnetic Analysis Applied to the Prediction of Stray Losses in Power Transformer L. Susnjic 1), Z. Haznadar 2) and Z. Valkovic 3) 1) Faculty of Engineering Vukovarska 58, 5 Rijeka, Croatia, e-mail:
More information3-D eddy current analysis in moving conductor of permanent magnet type of retarder using moving coordinate system
Physics Electricity & Magnetism fields Okayama University Year 1997 3-D eddy current analysis in moving conductor of permanent magnet type of retarder using moving coordinate system Kazuhiro Muramatsu
More informationPresented at the COMSOL Conference 2010 Paris
Presented at the COMSOL Conference 2010 Paris Fakultät Elektrotechnik und Informationstechnik Institut für Feinwerktechnik und Elektronik-Design Prof. Dr.-Ing. habil. Jens Lienig Robust and Reliability-Based
More informationMedical Physics & Science Applications
Power Conversion & Electromechanical Devices Medical Physics & Science Applications Transportation Power Systems 1-5: Introduction to the Finite Element Method Introduction Finite Element Method is used
More informationDETERMINING STATIC TORQUE CHARACTERISTICS OF PERMANENT MAGNET STEP MOTOR
DETERMINING STATIC TORQUE CHARACTERISTICS OF PERMANENT MAGNET STEP MOTOR Miroslav BJEKIĆ Abstract: The paper analyzes three methods of calculating static torque of permanent magnet step motor. Analytical
More informationEffect of Reaction Plate on Performance of Single-Side Linear Induction Motor in Different Speeds and Frequencies with Finite Element Method
Effect of Reaction Plate on Performance of Single-Side Linear Induction Motor in Different Speeds and Frequencies with Finite Element Method O. Maktabdar MSC Student, University of Birjand, Birjand, Iran,
More information1205. Optimal design of quadratic electromagnetic exciter
1205. Optimal design of quadratic electromagnetic exciter Ji Hyun Choi 1, Chang Hyun Park 2, Jin Ho Kim 3 1, 3 Electric Motor and Actuator Laboratory, Department of Mechanical Engineering, Yeungnam University
More informationAnalytical Model for Magnetic Anisotropy Dedicated to Non-Oriented Steel Sheets
Analytical Model for Magnetic Anisotropy Dedicated to Non-Oriented Steel Sheets Floran Martin, Deepak Singh, Anouar Belahcen, Paavo Rasilo, Ari Haavisto, and Antero Arkkio Aalto University - Departement
More informationExperimental and Theoretical Study of Interlaminar Eddy Current Loss in Laminated Cores
Tampere University of Technology Experimental and Theoretical Study of Interlaminar Eddy Current Loss in Laminated Cores Citation Bikram Shah, B., Rasilo, P., Belahcen, A., & Arkkio, A. (207). Experimental
More informationModeling of tape wound cores using reluctance networks for distribution transformers
No. E-4-AAA- Modeling of tape wound cores using reluctance networks for distribution transformers Amir Baktash, Abolfazl Vahedi Centre of Excellence for Power System Automation and Operation Electrical
More informationEnergy-Based Variational Model for Vector Magnetic Hysteresis
Energy-Based Variational Model for Vector Magnetic Hysteresis L. Prigozhin 1 V. Sokolovsky 1 J. W. Barrett 2 S. E. Zirka 3 1 Ben-Gurion University of the Negev, Israel 2 Imperial College London, UK 3 Dnepropetrovsk
More informationAGENERAL approach for the calculation of iron loss in
IEEE TRANSACTIONS ON MAGNETICS, VOL. 51, NO. 1, JANUARY 2015 6300110 Comparison of Iron Loss Models for Electrical Machines With Different Frequency Domain and Time Domain Methods for Excess Loss Prediction
More informationSOLIDIFICATION SURFACE SPEED CONTROL OF FERROMAGNETIC PIECES USING EDDY CURRENT HEATING
SOLIDIFICATION SURFACE SPEED CONTROL OF FERROMAGNETIC PIECES USING EDDY CURRENT HEATING MIHAI MARICARU, MARILENA STĂNCULESCU, 1 VALERIU ŞTEFAN MINCULETE, 1 FLOREA IOAN HĂNŢILĂ 11 Key words: Coupled eddy
More informationFinite Element Method based investigation of IPMSM losses
Finite Element Method based investigation of IPMSM losses Martin Schmidtner 1, Prof. Dr. -Ing. Carsten Markgraf 1, Prof. Dr. -Ing. Alexander Frey 1 1. Augsburg University of Applied Sciences, Augsburg,
More informationME 608 Numerical Methods in Heat, Mass, and Momentum Transfer
ME 608 Numerical Methods in Heat, Mass, and Momentum Transfer Solution to Final Examination Date: May 4, 2011 1:00 3:00 PM Instructor: J. Murthy Open Book, Open Notes Total: 100 points Use the finite volume
More informationAnalysis and Case Study of Permanent Magnet Arrays for Eddy Current Brake Systems with a New Performance Index
Journal of Magnetics 18(3), 276-282 (2013) ISSN (Print) 1226-1750 ISSN (Online) 2233-6656 http://dx.doi.org/10.4283/jmag.2013.18.3.276 Analysis and Case Study of Permanent Magnet Arrays for Eddy Current
More informationTitle. Author(s)Igarashi, Hajime; Watanabe, Kota. CitationIEEE Transactions on Magnetics, 46(8): Issue Date Doc URL. Rights.
Title Complex Adjoint Variable Method for Finite-Element A Author(s)Igarashi, Hajime; Watanabe, Kota CitationIEEE Transactions on Magnetics, 46(8): 2739-2742 Issue Date 2010-08 Doc URL http://hdlhlenet/2115/45008
More informationScientific Computing
Lecture on Scientific Computing Dr. Kersten Schmidt Lecture 4 Technische Universität Berlin Institut für Mathematik Wintersemester 2014/2015 Syllabus Linear Regression Fast Fourier transform Modelling
More informationINFLUENCE OF MAGNETIC ANISOTROPY ON THE ROTATIONAL MAGNETIZATION IN TYPICAL DYNAMO STEEL SHEETS
Zeszyty Problemowe Maszyny Elektryczne Nr 3/014 (103) 173 Witold Mazgaj, Zbigniew Szular, Adam Warzecha Cracow University of Technology INFLUENCE OF MAGNETIC ANISOTROPY ON THE ROTATIONAL MAGNETIZATION
More informationComputation and Analysis of the DC-Biasing Magnetic Field by the Harmonic-Balanced Finite-Element Method
International Journal of Energy and Power Engineering 2016; 5(1-1): 31-36 Published online October 10, 2015 (http://www.sciencepublishinggroup.com/j/ijepe) doi: 10.11648/j.ijepe.s.2016050101.14 ISS: 2326-957X
More information3D-FEM-Simulation of Magnetic Shape Memory Actuators
3D-FEM-Simulation of Magnetic Shape Memory Actuators M.Schautzgy* 1/2, U.Kosiedowski² and T.Schiepp 1 1 ETO MAGNETIC GmbH, Hardtring 8, 78333 Stockach, Germany 2 University of Applied Science Konstanz,
More informationModelling of Specimen Interaction with Ferrite Cored Coils by Coupling Semi-Analytical and Numerical Techniques
Modelling of Specimen Interaction with Ferrite Cored Coils by Coupling Semi-Analytical and Numerical Techniques A. Skarlatos, E. Demaldent, A. Vigneron and C. Reboud 25 th 28 th of June, Bratislava, Slovak
More informationComparison between Analytic Calculation and Finite Element Modeling in the Study of Winding Geometry Effect on Copper Losses
Comparison between Analytic Calculation Finite Element Modeling in the Study of Winding Geometry Effect on Copper Losses M. Al Eit, F. Bouillault, C. March, G. Krebs GeePs Group of electrical engineering
More informationMulti-Scale FEM and Magnetic Vector Potential A for 3D Eddy Currents in Laminated Media
Multi-Scale FEM and Magnetic Vector Potential A for 3D Eddy Currents in Laminated Media K. Hollaus, and J. Schöberl Institute for Analysis and Scientific Computing, Wiedner Hauptstr. 8, A-14 Vienna, Austria
More informationDESIGN FEATURES AND GOVERNING PARAMETERS OF LINEAR INDUCTION MOTOR
CHAPTER 5 DESIGN FEATURES AND GOVERNING PARAMETERS OF LINEAR INDUCTION MOTOR 5.1 Introduction To evaluate the performance of electrical machines, it is essential to study their electromagnetic characteristics.
More informationIEEE TRANSACTIONS ON POWER DELIVERY, VOL. 22, NO. 1, JANUARY /$ IEEE
IEEE TRANSACTIONS ON POWER DELIVERY, VOL. 22, NO. 1, JANUARY 2007 195 Analysis of Half-Turn Effect in Power Transformers Using Nonlinear-Transient FE Formulation G. B. Kumbhar, S. V. Kulkarni, Member,
More informationModeling and Analysis of Leakage Flux and Iron Loss Inside Silicon Steel Laminations
International Journal of Energy and Power Engineering 6; 5(-): 8-5 Published online November 5, 5 (http://www.sciencepublishinggroup.com/j/ijepe) doi:.68/j.ijepe.s.65.7 ISSN: 36-957X (Print); ISSN: 36-96X
More informationThe initial magnetization curve shows the magnetic flux density that would result when an increasing magnetic field is applied to an initially
MAGNETIC CIRCUITS The study of magnetic circuits is important in the study of energy systems since the operation of key components such as transformers and rotating machines (DC machines, induction machines,
More informationSHAPE DESIGN OPTIMIZATION OF INTERIOR PERMANENT MAGNET MOTOR FOR VIBRATION MITIGATION USING LEVEL SET METHOD
International Journal of Automotive Technology, Vol. 17, No. 5, pp. 917 922 (2016) DOI 10.1007/s12239 016 0089 7 Copyright 2016 KSAE/ 092 17 pissn 1229 9138/ eissn 1976 3832 SHAPE DESIGN OPTIMIZATION OF
More informationMagnetostatic Analysis of Solenoid
Magnetostatic Analysis of Solenoid 1. Introduction 2. Model View 3. Materials 4. Load & Restraint Information 5. Coils Information 6. Force and Torque Information 7. Study Properties 8. Results Table 9.
More informationTime-Harmonic Modeling of Squirrel-Cage Induction Motors: A Circuit-Field Coupled Approach
Time-Harmonic Modeling of Squirrel-Cage Induction Motors: A Circuit-Field Coupled Approach R. Escarela-Perez 1,3 E. Melgoza 2 E. Campero-Littlewood 1 1 División de Ciencias Básicas e Ingeniería, Universidad
More informationDefinition Application of electrical machines Electromagnetism: review Analogies between electric and magnetic circuits Faraday s Law Electromagnetic
Definition Application of electrical machines Electromagnetism: review Analogies between electric and magnetic circuits Faraday s Law Electromagnetic Force Motor action Generator action Types and parts
More informationCOIL GUN is an importance member in electromagnetic
94 IEEE TRANSACTIONS ON PLASMA SCIENCE, VOL. 39, NO. 1, JANUARY 2011 Application of FE-BECM in Field Analysis of Induction Coil Gun Shoubao Liu, Jiangjun Ruan, Yadong Zhang, Yu Zhang, and Yujiao Zhang
More informationNEURAL NETWORK MODEL FOR SCALAR AND VECTOR HYSTERESIS
Journal of ELECTRICAL ENGINEERING, VOL. 54, NO. 1-2, 2003, 13 21 NEURAL NETWORK MODEL FOR SCALAR AND VECTOR HYSTERESIS Miklós Kuczmann Amália Iványi The Preisach model allows to simulate the behaviour
More informationANALYSIS OF FLUX DENSITY AND FORCE ON DOUBLE SIDED LINEAR INDUCTION MOTOR WITH DIFFERENT POLE SHAPES
ANALYSIS OF FLUX DENSITY AND FORCE ON DOUBLE SIDED LINEAR INDUCTION MOTOR WITH DIFFERENT POLE SHAPES Chetan Vasudeva Sanjay Marwaha Gurmeet Singh D.A.V University, Sarmastpur, S.L.I.E.T, Longowal, Sangrur,
More informationA Symmetric and Low-Frequency Stable Potential Formulation for the Finite-Element Simulation of Electromagnetic Fields
A Symmetric and Low-Frequency Stable Potential Formulation for the Finite-Element Simulation of Electromagnetic Fields Martin Jochum, Ortwin Farle, and Romanus Dyczij-Edlinger Abstract A low-frequency
More informationProjectile s Velocity Effect for Voltage Induced at Sensing Coil for Applying to Air Bursting Munition
Journal of Magnetics 18(2), 90-94 (2013) ISSN (Print) 1226-1750 ISSN (Online) 2233-6656 http://dx.doi.org/10.4283/jmag.2013.18.2.090 Projectile s Velocity Effect for Voltage Induced at Sensing Coil for
More informationRobust optimal design of a magnetizer to reduce the harmonic components of cogging torque in a HDD spindle motor
Microsyst Technol (2014) 20:1497 1504 DOI 10.1007/s00542-014-2153-4 Technical Paper Robust optimal design of a magnetizer to reduce the harmonic components of cogging torque in a HDD spindle motor Changjin
More informationAn inductance lookup table application for analysis of reluctance stepper motor model
ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR
More informationMathematical Modelling and Simulation of Magnetostrictive Materials by Comsol Multiphysics
Excerpt from the Proceedings of the COMSOL Conference 008 Hannover Mathematical Modelling and Simulation of Magnetostrictive Materials by Comsol Multiphysics Author: M. Bailoni 1, Y.Wei, L. Norum 3, 1
More informationCoupled electromagnetic, thermal and stress analysis of large power electrical transformers
Coupled electromagnetic, thermal and stress analysis of large power electrical transformers DANIELA CÂRSTEA High-School Group of Railways, Craiova ROMANIA ALEXANDRU ADRIAN CÂRSTEA University of Craiova
More informationMOTIONAL MAGNETIC FINITE ELEMENT METHOD APPLIED TO HIGH SPEED ROTATING DEVICES
MOTIONAL MAGNETIC FINITE ELEMENT METHOD APPLIED TO HIGH SPEED ROTATING DEVICES Herbert De Gersem, Hans Vande Sande and Kay Hameyer Katholieke Universiteit Leuven, Dep. EE (ESAT), Div. ELEN, Kardinaal Mercierlaan
More information2577. The analytical solution of 2D electromagnetic wave equation for eddy currents in the cylindrical solid rotor structures
2577. The analytical solution of 2D electromagnetic wave equation for eddy currents in the cylindrical solid rotor structures Lale T. Ergene 1, Yasemin D. Ertuğrul 2 Istanbul Technical University, Istanbul,
More informationMSFEM with Network Coupling for the Eddy Current Problem of a Toroidal Transformer
MSFEM with Network Coupling for the Eddy Current Problem of a Toroidal Transformer Markus Schöbinger, Simon Steentjes, Joachim Schöberl, Kay Hameyer, Karl Hollaus July 18, 2017 ISTET 2017 Illmenau July
More informationLevitation force analysis of medium and low speed maglev vehicles
Journal of Modern Transportation Volume 2, Number 2, June 212, Page 93-97 Journal homepage: jmt.swjtu.edu.cn DOI: 1.17/BF3325784 1 Levitation force analysis of medium and low speed maglev vehicles Guoqing
More informationOperation of an Electromagnetic Trigger with a Short-circuit Ring
Operation of an Electromagnetic Trigger with a Short-circuit Ring Dejan Križaj 1*, Zumret Topčagić 1, and Borut Drnovšek 1,2 1 Faculty of Electrical Engineering, University of Ljubljana, Ljubljana, Slovenia,
More informationCOMPARISON OF DEFECT DETECTION IN ALUMINUM AND STEEL PLATES USING AN ELECTROMAGNETIC ACOUSTIC TRANSDUCER
The 12 th International Conference of the Slovenian Society for Non-Destructive Testing Application of Contemporary Non-Destructive Testing in Engineering September 4-6, 2013, Portorož, Slovenia COMPARISON
More informationExperimental Tests and Efficiency Improvement of Surface Permanent Magnet Magnetic Gear
IEEJ Journal of Industry Applications Vol.3 No.1 pp.62 67 DOI: 10.1541/ieejjia.3.62 Paper Experimental Tests and Efficiency Improvement of Surface Permanent Magnet Magnetic Gear Michinari Fukuoka a) Student
More informationChap. 1 Fundamental Concepts
NE 2 Chap. 1 Fundamental Concepts Important Laws in Electromagnetics Coulomb s Law (1785) Gauss s Law (1839) Ampere s Law (1827) Ohm s Law (1827) Kirchhoff s Law (1845) Biot-Savart Law (1820) Faradays
More informationFINITE ELEMENT MODELING OF THE BULK MAGNITIZATION OF RAILROAD WHEELS TO IMPROVE TEST CONDITIONS FOR MAGNETOACOUSTIC
FINITE ELEMENT MODELING OF THE BULK MAGNITIZATION OF RAILROAD WHEELS TO IMPROVE TEST CONDITIONS FOR MAGNETOACOUSTIC RESIDUAL STRESS MEASUREMENTS J. P. Fulton and B. Wincheski Analytical Services and Materials,
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationPower Loss Analysis of AC Contactor at Steady Closed State with Electromagnetic-Thermal Coupling Method
Journal of Information Hiding and Multimedia Signal Processing c 2017 ISSN 2073-4212 Ubiquitous International Volume 8, Number 2, March 2017 Power Loss Analysis of AC Contactor at Steady Closed State with
More informationDesign and development of coil casing MRF brake system
Design and development of coil casing MRF brake system R. Siti Lydia 1,*, R. Mokhtar 1, A.B. Muhamad Husaini 1 and B.J. Norhaniza 1 1 System Engineering and Energy Laboratory, Section of Manufacturing,
More informationComparison of Magnetic Field and Iron Loss Models for 2-Dimensional Magnetic Anisotropic Properties
Kiyoshi WAJIMA, Teruyuki TAMAKI Nippon Steel Corporation Comparison of Magnetic Field and Iron Loss Models for 2-Dimensional Magnetic Anisotropic Properties Abstract. In numerical analysis of electromagnetic
More information