A New Approach in Analytical Analysis of Eddy Currents in Laminated Core

Similar documents
Design and Application of Fault Current Limiter in Iran Power System Utility

Lecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b.

CALCULATION OF 2D-THERMOMAGNETIC CURRENT AND ITS FLUCTUATIONS USING THE METHOD OF EFFECTIVE HAMILTONIAN. R. G. Aghayeva

Homework 7 Due 18 November at 6:00 pm

Chapter 6. Electromagnetic Oscillations and Alternating Current

Magnetic field generated by current filaments

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

EE 5337 Computational Electromagnetics

INFLUENCE OF IMPERFECT BONDING ON INTERFACE OF MAGNETO- ELECTRO-ELASTIC HETEROSTRUCTURES: SH WAVES DISPERSION RELATIONS

Vortex Shedding on Combined Bodies at Incidence to a Uniform Air Stream T. Yavuz x, Y. E. Akansu xx, M. Sarıo lu xxx, and M.

An inductance lookup table application for analysis of reluctance stepper motor model

A Design Method for USBL Systems with Skew Three-element Arrays

State-Space Model for a Multi-Machine System

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

A Finite volume simulation of a supersonic laminar Flow: Application to a compression corner

Vertical Sluice Gate Discharge Coefficient

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

IPMSM Inductances Calculation Using FEA

Physics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2

Efficient Macro-Micro Scale Coupled Modeling of Batteries

Table of Common Derivatives By David Abraham

Qubit channels that achieve capacity with two states

Simple Relations between a Uniaxial Medium and an Isotropic Medium

Negative-Index Refraction in a Lamellar Composite with Alternating. Single Negative Layers

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Sensors & Transducers 2015 by IFSA Publishing, S. L.

1 Lecture 20: Implicit differentiation

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

Basic Differentiation Rules and Rates of Change. The Constant Rule

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann

Study on aero-acoustic structural interactions in fan-ducted system

Influence of Radiation on Product Yields in a Film Boiling Reactor

A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation

Chapter 4. Electrostatics of Macroscopic Media

5-4 Electrostatic Boundary Value Problems

Nonlinear Dielectric Response of Periodic Composite Materials

An Approach for Design of Multi-element USBL Systems

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

EFFECTS OF RADIATION AND THERMAL DIFFUSION ON MHD HEAT TRANSFER FLOW OF A DUSTY VISCOELASTIC FLUID BETWEEN TWO MOVING PARALLEL PLATES

EXACT TRAVELING WAVE SOLUTIONS FOR A NEW NON-LINEAR HEAT TRANSFER EQUATION

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1

The Second Order Contribution to Wave Crest Amplitude Random Simulations and NewWave

Open Access An Exponential Reaching Law Sliding Mode Observer for PMSM in Rotating Frame

6. Friction and viscosity in gasses

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges

Lecture Notes: March C.D. Lin Attosecond X-ray pulses issues:

Prep 1. Oregon State University PH 213 Spring Term Suggested finish date: Monday, April 9

1. The electron volt is a measure of (A) charge (B) energy (C) impulse (D) momentum (E) velocity

Applications of First Order Equations

Maxwell s Equations 5/9/2016. EELE 3332 Electromagnetic II Chapter 9. Maxwell s Equations for static fields. Review Electrostatics and Magnetostatics

Capacitance and Dielectrics

Computed Tomography Notes, Part 1. The equation that governs the image intensity in projection imaging is:

Dynamics of the synchronous machine

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann

The effect of nonvertical shear on turbulence in a stably stratified medium

G j dq i + G j. q i. = a jt. and

EE 370L Controls Laboratory. Laboratory Exercise #7 Root Locus. Department of Electrical and Computer Engineering University of Nevada, at Las Vegas

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

Optimal LQR Control of Structures using Linear Modal Model

A General Analytical Approximation to Impulse Response of 3-D Microfluidic Channels in Molecular Communication

Implicit Differentiation

1 M3-4-5A16 Assessed Problems # 1: Do 4 out of 5 problems

CE 579: STRUCTRAL STABILITY AND DESIGN

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.

Situation awareness of power system based on static voltage security region

Approaches for Predicting Collection Efficiency of Fibrous Filters

On Using Unstable Electrohydraulic Valves for Control

SiC-based Power Converters for High Temperature Applications

Simulation of Angle Beam Ultrasonic Testing with a Personal Computer

Statics, Quasistatics, and Transmission Lines

Outline. Calculus for the Life Sciences II. Introduction. Tides Introduction. Lecture Notes Differentiation of Trigonometric Functions

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

Physics 115C Homework 4

1 Introuction To ene everthing precisel we consier a specic moel. Let us concentrate our attention on Bose gas with elta interaction (quantum Nonlinea

MEASUREMENT OF LOSSES IN DIRECT TORQUE CONTROLLED PERMANENT MAGNET SYNCHRONOUS MOTOR BY ANALYSIS AND EXPERIMENTAL SETUP

water adding dye partial mixing homogenization time

10. Magnetism. ) it is. S G appropriate to call the magnetic pole

Application of Filippov Theory to the IEEE Standard Anti-windup PI Controller

Summary of the Class before Exam1

Math 342 Partial Differential Equations «Viktor Grigoryan

CHAPTER 32. Answer to Checkpoint Questions

Investigation of local load effect on damping characteristics of synchronous generator using transfer-function block-diagram model

and from it produce the action integral whose variation we set to zero:

Characterization of lead zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy

Impurities in inelastic Maxwell models

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Lecture 10 Notes, Electromagnetic Theory II Dr. Christopher S. Baird, faculty.uml.edu/cbaird University of Massachusetts Lowell

Lagrangian and Hamiltonian Mechanics

Sph4c Chapter 2 Simple Machines LoRusso

Static Equilibrium. Theory: The conditions for the mechanical equilibrium of a rigid body are (a) (b)

Hyperbolic Functions

1 Heisenberg Representation

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom

Interaction force in a vertical dust chain inside a glass box

To understand how scrubbers work, we must first define some terms.

Electromagnet Gripping in Iron Foundry Automation Part II: Simulation

A USER-FRIENDLY TOOL FOR SIMULATING THE TIME-DEPENDENT FIELD DISTRIBUTION IN PE INSULATION ON THE BASIS OF A PHYSICAL APPROACH

Linear First-Order Equations

Characteristics of Dark Current and Photocurrent in Superlattice Infrared Photodetectors

Transcription:

J. Basic. Appl. Sci. Res., (7)741-745, 1 1, TextRoa Publication ISSN 9-434 Journal of Basic an Applie Scientific Research www.textroa.com A New Approach in Analtical Analsis of E Currents in Laminate Core Jaber Merrikhi Ahangarkolaei *, Yaser Merrikhi Ahangarkolaei an Mir Maji Temoori Department of Electrical Engineering, Firoozkooh Branch, Islamic Aza Universit, Firoozkooh, IRAN ABSTRACT Conventional methos of e-currents moeling are base on solving iffusion equation, which its variable is the magnetic fiel intensit; an then calculating e currents istribution b using the Curl of magnetic fiel intensit. Using the Curl operator necessitates calculating the erivative of fiel intensit, which is obtaine b some numerical fiel analsis methos such as Finite Elements Metho (FEM). Taking erivative causes man isasters, such as iscontinuit in the current ensit istribution. In this paper, for the first time, iffusion equation is written for the electric fiel strength as the main variable. Base on the analtical formula erive from the propose metho, the same results are achieve. Simplicit of the propose metho enables us to emplo it in numerical analsis such as finite element metho. KEY WORDS: e current, laminate core, power loss, skin effect. INTRODUCTION Ferromagnetic materials, such as iron, nickel an their allos are wiel use for builing cores of magnetic evices [1]. In the esign of ac machines an rives, the preiction of power losses in magnetic cores is important []. The etermination of e currents in a given magnetic configuration is alwas a complicate task. As Depicte in Fig.1, the cores are usuall laminate to reuce e-current loss, however this loss can still be significant, especiall at the higher frequencies prouce b pulse-with moulation waveforms in tpical ac rives [1]. Fig.1: Laminate magnetic core Design engineers usuall preict lamination loss using two methos []. One wa is to use manufacturer s atasheets. Usuall in the atasheet of the core, a cluster of curves can be foun, which show the core loss per unit volume uner sine-wave excitation versus the maximum flux ensit for ifferent values of frequenc, but such curves are often not available or in practical applications are often ifferent from those uner which the atasheets ha been obtaine. The other metho is to use a simple formula for e-current loss as a function of lamination thickness, but it is onl vali for frequencies of about 6 Hz or less. In the case of most ac rives, significant frequenc content in the kilohertz an megahertz regions is require. At such high frequencies, e-current loss is much greater than hsteresis loss an, thus, core loss is almost entirel ue to e currents []. The subject of core losses in AC machines an in particular laminations has been a subject of research since the as of Steinmetz [3]. Steinmetz propose approximate formulas for the harmonic case which are quite similar to those currentl use toa. His work was empiricall erive an supporte from a theoretical viewpoint b Bertotti [4]. It is assume in [3] that the magnetic flux istribution in the sheet is uniform, making the metho applicable onl in the frequenc range where the skin effect is negligible or to the magnetic cores with low conuctivities or thin sheets [5]. This latter formulation i not prove accurac enough which motivate Bertotti to propose an aitional term, the excess loss term which accounte for the ifference between measurements an theor [6]. * Corresponing Author: Jaber Merrikhi Ahangarkolaei, Department of Electrical Engineering, Firoozkooh Branch, Islamic Aza Universit, Firoozkooh, IRAN. E-mail: [jmerrikhi@iaufb.ac.ir] 741

Ahangarkolaei et al., 1 Moreover, in the last ecaes, a great number of papers have been evote to obtain approximate e-current losses. Conventional methos of moeling e-currents are base on solving iffusion equation for magnetic fiel intensit. The e currents istribution is then obtaine b taking the Curl of magnetic fiel intensit. Using the Curl operator necessitates calculating the erivative of fiel intensit. Taking erivative causes man problems, such as iscontinuit in the current ensit istribution an amplification of the errors [7]. The objective of this paper is to present a new analtical metho of e current ensit an e current core power loss from Integral form of Faraa's law (Eq.1), that can be use in numerical analsis. E. l (1) t Where H, J, B, E are the magnetic fiel strength, current ensit, magnetic flux ensit, an electric fiel strength respectivel. 1.Calculation of e current in a single sheet For a lamination with thickness an height l as shown in Fig., an assuming constant magnetic permeabilit an conuctivit, electrical fiel onl has the z irection component as a function of x. Fig.: Single lamination an the coorinate sstem. Neglecting the isplacement current, Maxwell s equations an Ohm s law becomes: ( x, t ) ( x, t ) () t x Letting E(t) be sinusoial an assuming the complex exponential form: j t ( x, t) Re( ( x) e ) (3) Where ω is the angular frequenc. Combining the (Eq.3), (Eq.4) equations, we have : (, ) E (, ) z x t j x t (4) x Stu of Fig. will show that, bounar conitions exist for thin sheets with an axiall applie fiel: E E (5) x (6) z z (7) ( x ) E (8) ( x ) E (9) Eq. 4 is a stanar secon orer ifferential equation, which has a general solution an it is convenient to write the solution of ifferential equation in term of hperbolic functions. A sinh( Tx ) B cosh( Tx ) (1) Where T is constant an is equal to: T (1 j) f j (11) Which A an B being arbitrar constant values. Base on the bounar conition, E iels sinh( Tx) E (1) sinh( T ) 74

J. Basic. Appl. Sci. Res., (7)741-745, 1 There remains onl E to be evaluate. Using Faraa s line integral an integrating aroun the path ABCD in Fig., iel: E. l (13) t Where ϕ is the total flux within the line ABCD. If the height is much greater than the thickness of sheet, the electrical fiel E in complex notation form is: j E (14) l In the case shown in Fig., The magnetic fiel intensit is escribe b an orinar secon-orer ifferential equation, which is the Helmholtz equation, H ( x, t) H ( x, t) (15) t x For a sinusoial applie magnetic strength, Eq.5 becomes. H ( x, t ) j H ( x, t ) (16) x Solution of Eq.17 can be foun in previous articles an becomes [8]: cosh( Tx) H ( x) H (17) cosh( T ) Where, H is the bounar conition. H( ) H( ) H (18) The magnetic fiel intensit as a function of x is shown in Fig.3 for H =1(A/m), σ= 1.4e6 an µ r =3 at f=1, an 5 khz. It can be seen that the magnetic fiel intensit ecreases with increasing frequenc f in the center of the lamination. Fig.3: The magnetic fiel intensit as a function of x at f=1, an 5 khz The magnetic flux flowing through the cross section of a single lamination perpenicular to the xz plane is: / / l lh H ( x) s ( ) ( ) cosh( ) H x x l H x x Tx x s / cosh( T ) / 743

Ahangarkolaei et al., 1 T T lh sinh( ) sinh( ) lh T tanh( ) (19) T T cosh( ) T Substitution of (Eq.) an (Eq.15) into (Eq.13) an J E prouces: j H sinh( ) sinh( ) T Tx Tx J z ( x) tanh( ) TH () T sinh( T ) cosh( T ) The magnitue of the e current ensit as a function of x is shown in Fig.4, for H =1(A/m), σ= 1.4e6 an µ r =3 at f=1, an 5 khz. Fig. 4: The magnitue of the e current ensit as a function of x at f=1, an 5 khz.e current power loss The e current power loss ensit istribution in a conucting core average with respect of time is: 3 P ( x ) E. J J ( W / m ) (1) Substitution of (Eq.) into (Eq.1) prouces: 1 sinh( Tx) P( x) TH () cosh( T ) The average power loss ensit over the lamination with is the power loss per unit volume, calle the specific core loss. Its average value in the sheet is: / / 1 1 sinh( Tx) TH sinh( T ) T Pv P( x) x TH x (3) / T cosh( T ) cosh ( ) Next, to emonstrate the valiit of Eq.3, the result of Eq.3 is compare with the e current loss from the Curl of magnetic fiel intensit. RESULTS Recall that Eq. 3 is obtaine from integral form of Maxwell's equations. In the latest research Instea of Eq.3, in [8] Bermúez proposes Eq.4 that obtaine from Curl of magnetic fiel intensit. 744

J. Basic. Appl. Sci. Res., (7)741-745, 1 H sinh( ) sin( ) Pv W m cosh( ) cos( ) 3 ( / ) Where f. Bermúez in [8] showe that Eq.4 is exact for an frequenc. In Fig. 5, we show the losses given b the above analtical formulas versus frequenc in a last given lamination. Can be seen, the results of Eq.3 is completel consistent with the results of Eq. 4. (4) Conclusion Fig. 5: Average e-current loss versus f(hz) In this paper, a new analtical consieration of e current loss is evaluate. For this purpose, Integral form of Faraa s law is emploe to calculate e current ensit an e current core power loss. Utilization Steinmetz an Bertotti theoretical methos, a comparison between this metho an their methos are accomplishe. Base on the analtical formula erive from the propose metho, the same results are achieve. Simplicit of the propose metho enables us to emplo it in numerical analsis such as finite element metho. Acknowlegment This research work has been accomplishe at Department of electrical engineering, Firoozkooh branch, Islamic Aza Universit, Firoozkooh, IRAN. Authors woul like to thanks from its financial support. REFERENCES 1. J. Xu, A. Lakhsasi, Z. Yao an V. Rajagopalan,1996. A Practical Moeling Metho for E-Current Losses Computation in Laminate Magnetic Cores. IEEE Inustr Applications Conference :153-1537. John R. Brauer, Zoltan J. Cenes, Bruce C. Beihoff, an Kenneth P. Phillips,. Laminate Steel E-Current Loss versus Frequenc Compute Using Finite Elements. IEEE Transactions On Inustr Applications, Vol. 36, No. 4 : 113-1137 3. C. P. Steinmetz,1984. On the law of hsteresis.proc. IEEE, vol. 7, no. : 196-1 4. G. Bertotti, 1988. General properties of power losses in soft ferromagnetic materials. IEEE Trans. Magn. Vol. 4 : 61-63 5. Ema Dlala, 8. A Simplifie Iron Loss Moel for Laminate Magnetic Cores. IEEE Transactions On Magnetics, Vol. 44, No. 11: 3169-317 6. Yu Zhang, Rui Guan, Pragasen Pilla an Ming-C. Cheng, 9. General Core Loss Moels on a Magnetic Lamination. IEMDC '9. IEEE International Electric Machines an Drives Conference, 9. : 159-1534 7. J.Merrikhi, J.S.Moghani, E.Fallah, 8. A New Approach in Finite Element Analsis of E Current with Flux Skin Effect. Amirkabir International Journal of Science an Technolog : 1-7 8. Alfreo Bermúez, Dolores Gómez, an Pilar Salgao, 8. E-Current Losses in Laminate Cores an the Computation of an Equivalent Conuctivit. IEEE Transactions On Magnetics, Vol. 44, No. 1: 473-4738 745