Theory and Practice – Making Use of the Barkhausen Effect
Total Page:16
File Type:pdf, Size:1020Kb
learn invent impact Theory and Practice – Making use of the Barkhausen Effect David C. Jiles Anson Marston Distinguished Professor Palmer Endowed Chair Department of Electrical & Computer Engineering Iowa State University Workshop on Large Fluctuations and Collective Phenomena in Disordered Materials May 18, 2011 Summary • Stochastic/deterministic Barkhausen model – Equations describing the phenomenon – Comparison of theory and experiment • Variation of emissions with other factors such as stress – Variation of Barkhausen signal amplitude – Theoretical predictions – inverse law • Variation of signal with frequency and distance – Possibilities for depth profiling of properties – Challenges 2 www.ece.iastate.edu learn invent impact The variation of magnetization with magnetic field looks deceptively simple Yet behind it lies a large number of very complicated interrelated mechanisms which are hard to model 1 2 3 4 A further problem is that magnetic behavior is inhomogeneous and does not “scale” easily as dimensions change The result has been a number of different models which do not fit together very well www.ece.iastate.edu learn invent impact Langevin-Weiss Model • An array of magnetic moments at temperature T and in a magnetic field H will distribute themselves among the available energy states according to the probability distribution P P = P0 exp (− µ0mH / kBT ) • Integrating over all possible angles and normalizing gives M = Nm [coth (µ0mH / kBT )− (kBT / µ0mH )] • Weiss extended this to include internal coupling proportional to the magnetization He = aM with the result M = Nm [coth (µ0m(H +αM )/ kBT )− (kBT / µ0m(H +αM ))] which decribes the anhysteretic magnetization curve for ferromagnets www.ece.iastate.edu learn invent impact Anhysteretic magnetization curves with different temperatures and different anisotropies M M Low temperature Ms + Ms 1.0 High temperature -15 -10 -5 5 10 15 H H 3D 2D -1.0 1D "Superparamagnetic magnetization equation in two dimensions," D.C. Jiles, S.J. Lee, J. Kenkel and K. Metlov. Applied Physics Letters , 77, 1029, 2000. www.ece.iastate.edu learn invent impact Energy Dissipation and Hysteresis • Energy is only dissipated by irreversible magnetization processes • A model with dissipation proportional to change in M gives good results M M n επ Epin ()M = ∫dM irr =µok∫dM irr 2m 0 0 • Energy output (change in magnetostatic energy) must equal energy input minus energy dissipated due to any losses such as hysteresis µ Mirr dH e = µ Man dH e − µ k dM 0 ∫ 0 ∫0 ∫ irr www.ece.iastate.edu learn invent impact Irreversible Changes in Magnetization • Differentiating the energy equation lead to a differential for the rate of change of magnetization dM irr 1 = ()M an − M irr dH e k • This can be written in the form of a differential with respect to the applied field H dM irr 1 = ()Man −Mirr dH k −α()Man −Mirr www.ece.iastate.edu learn invent impact Isotropic Model of Hysteresis • The irreversible component of magnetization varies according to the differential equation dM irr 1 = ()Man −Mirr dH k −α()Man −Mirr • The reversible component of magnetization varies as dM dM dM rev = c an − irr dH dH dH • Therefore hysteresis can be represented in terms of Ms, a, α, k, and c using the equation dM 1 (M − M ) c dM = an + an dH ()1+ c k −α ()M an − M ()1+ c dH "Determination of theoretical parameters for modelling bulk magnetic hysteresis properties using the theory of ferromagnetic hysteresis," D.C. Jiles, J.B. Thoelke and M.K. Devine. IEEE Trans. Mag. www.ece.iastate.edu learn invent impact Comparison of Model and Experimental Measurement • Comparison of measured and modelled hysteresis curves of cobalt modified gamma iron oxide material After P. Andrei and A. Stancu, “Hysteresis in particulate recording media. Experiment and simulation with Preisach and Jiles-Atherton models,” Journal of Magnetism and Magnetic Materials , 206 , 160-164, 1999. www.ece.iastate.edu learn invent impact Other magnetic properties - Barkhausen effect 1.5 H Virr 1.0 0.5 V 0.0 -0.5 Voltage Max/ Voltage -1.0 Barkhausen -1.5 0Time (s) 0.006 Barkhausen signal shown as voltage against time. The time dependence of magnetic field H, voltage in flux coil V, and Changes in Barkhausen emissions in envelope of irreversible component of voltage in flux coil V irr carbon steel as a result of applied are also shown on the same time scale. tensile stress. www.ece.iastate.edu learn invent impact Stochastic Process Barkhausen model Model assumes Barkhausen activity is proportional to the rate of irreversible change in magnetization • • <M >: Average discontinuous dM BE ' disc ∝ Mirr = χirr H MBE = N < Mdisc > change in magnetization dt due to Barkhausen jump Avalanches: Number of Barkhausen events N(t n) in a given time interval is correlated with number of events in the previous time interval tn-1 N(tn) = N(tn−1)+∆N(tn−1) δrand : Random number ∆N(tn−1) =δrand N(tn−1) (-1.47 < δrand < 1.47) • • ' ' ' MBE (tn) =< Mdisc > χrr HN (tn−1)+δrand N (tn−1) www.ece.iastate.edu learn invent impact Comparison of model with experiment www.ece.iastate.edu learn invent impact Other Factors – Effects of Stress • Applied stress can be treated in most respects like an effective magnetic field which changes the anisotropy of the material 3 σ ∂λ Hσ = 2 µo ∂M T • Later this was extended to cover the case of a uniaxial stress at an arbitrary direction to the applied magnetic field 3 σσσ 2 2 ∂λ∂λ∂λ Hσσσ()θθθ = ()cos θθθ − νννsin θθθ 2µµµo ∂∂∂M T σσσ is the stress, θθθ is the angle between the stress axis and the direction of H σσσ, and ν is Poisson’s ratio. www.ece.iastate.edu learn invent impact Modelled Effects of Stress on Magnetic Moment Orientation Positive magnetostriction stress applied along the y axis www.ece.iastate.edu learn invent impact Detection of Stress Specimens calibration tested under load in Instron Tensile Test machine BN sensor in the middle of gauge length of sample Typical Magnetising and Analysis Conditions : Magnetising frequency : f = 300 Hz Magnetising voltage : V = 12 volt No of bursts : 20 Analysing bandwidth frequency : fa = 20-1250 KHz (about nearest 100 µm) Smoothing parameter: sp = 100 Sampling frequency: fs =2.5 MHz 15 www.ece.iastate.edu learn invent impact Measured Barkhausen Results Core Surface Strip #3#20#36 - -Depth Depth = = 22.9mm 11.2mm0.2mm 1.5 RMSRMS = =0.187 0.0290.010 volt volt 1.0 0.5 0.0 0.0 -0.5 -0.5 BarkhausenBarkhausen effecteffect signalsignal (V) (V) -1.0 Barkhausen effect signal (V) -1.0 -1.5 0.00 0.02 0.04 0.06 0.08 0.10 -1.5 0.00 0.02 0.04Time (sec) 0.06 0.08 0.10 Time (sec) 16 www.ece.iastate.edu learn invent impact Effects of Stress on Barkhausen Emissions Applied stress causes a change in anisotropy energy of magnetic moments 3 E = − σλ (cos 2 θ − v sin 2 θ ) σ 2 This can be expressed as an equivalent field 1 ∂Eσ 3 σ ∂λ Hσ = − = µ0 ∂M 2 µ0 ∂M σ The total field is then the sum of magnetic field, exchange field and “stress-equivalent field” H σ 3λ σ H H H M H s M M e = + σ +α = + 2 + α µoMs The magnetization is then an anisotropic function of this total field H + Hσ +αM a M an ()H,σ = M s coth − a H + Hσ +αM www.ece.iastate.edu learn invent impact Effects of Stress on Barkhausen Emissions dM 1 1 ′ an χ ()σ = ()H,σ ≅ Ms = Ms dH 3a − ()α +α M 3 σ s λsσ 3a − 2 +αMs µoMs This predicts how the rate of change of magnetization with field depends on stress and so can be used to calculate stress. However there is a much easier way... 3λsσ 3a−α + 2 Ms 1 3a−()α +α M µ M = σ s = o s χ′(σ) Ms Ms 1 1 3λsσ which predicts a straight line graph = − 2 χ′(σ) χ′ )0( µoMs of 1/χ’ against σ 18 www.ece.iastate.edu learn invent impact Detection of Tensile Stress - test results Applying tensile stress with monotonically increasing load 0 M P a 1 0 0 0 50 M P a 100 M P a 9 0 0 150 M P a 200 M P a 8 0 0 250 M P a 300 M P a 7 0 0 350 M P a 400 M P a 450 M P a 6 0 0 500 M P a 550 M P a 5 0 0 600 M P a 650 M P a 4 0 0 700 M P a 750 M P a 3 0 0 800 M P a BN peak amplitudepeak (mp) BN 850 M P a 2 0 0 900 M P a 950 M P a 1 0 0 Yield strength = 993MPa 0 -100 -80 -60 -40 -20 0 20 40 60 80 100 Magnetic Field Strength H (in % ) 19 www.ece.iastate.edu learn invent impact Dependence of Barkhausen emissions on stress Fig. 1. Envelope curves of the rectified MBN bursts for carburized SAE 9310 Fig. 2. MBN Peak Amplitude for carburized specimen for different amplitudes of applied SAE 9310 specimen as a function of applied stress using a tensile test machine (residual stress. stress value from XRD measured before stress application was -805MPa). 20 www.ece.iastate.edu learn invent impact Effects of Stress on Differential Susceptibility and Barkhausen Emissions • Barkhausen voltage V MBE is known from previous work* to be proportional to the differential susceptibility • Therefore, a similar linear expression should hold for reciprocal Barkhausen voltage as a function of stress as for the reciprocal differential susceptibility 1 1 3b′σ = − VMBE (σ) VMBE )0( µo • This suggests that the most useful calibration curve for Barkhausen effect as a function of stress is the reciprocal plot www.ece.iastate.edu learn invent impact Comparison of Model equation with Measurements www.ece.iastate.edu learn invent impact Conclusions • Magnetic properties, including permeability and Barkhausen effect depend on other external factors such as stress.