Electronics > Projects, Designs, and Technical Stuff
K_fe: Core Loss Coefficient of an Inductor
TimNJ:
Hi all,
I'm working through the design of a resonant AC inductor for an LLC power supply. This is basically my first time doing this.
I'm trying out the design procedure outlined in Ch 15. of 'Fundamentals of Power Electronics' by Erickson and Maksimovic. There is a parameter K_fe, "core loss coefficient" which is necessary to determine the optimal AC flux density for lowest total losses. The book, and countless white-papers I've found online, say "K_fe can be determined from the manufacturer's data" but none of them explain how to do this. (Someone on EDABoard asked this same exact question without getting a good answer: https://www.edaboard.com/showthread.php?379779-How-to-calculate-core-loss-coefficient).
I can determine Beta, the "Core Loss Exponent" by finding the slope of the log-log core loss plot. Does anyone know how to determine K_fe, which has the following units?
I think K_fe depends on B, flux density, so I'm a little confused on how this equation (with output AC peak flux density) also would have a flux density input.
Thank you.
TimNJ:
I decided to re-familiarize myself with some math concepts that I haven't touched in 10 years or so, since high school. Here's what I came up with.
First off, I believe the "core loss coefficient" to be k in the Steinmetz equation.
In a log-log plot, like those typically used to describe a magnetic material's core loss, you can determine the exponent of the power function by determining the slope of the log-log plot. You can do this by evaluating b = log(F1/F0) / log(X1/X0), where (X0, F0) and (X1, F1) are a set of points that exist on the plot.
For the 100KHz, 100C plot, I found b to be 4.16.
Here's where my confusion remains...
I'm taking a flux density of 200mT as a nominal operating value, 100KHz switching frequency. If we want to determine the coefficient k, we can use the following (In this case m is equal to b above):
I believe this allows you to calculate the coefficient. However, I've tried finding the coefficient at different points on the plot, and the numbers seem to vary.
For example, on the 100Khz/100C plot...
At (300mT, 8.1W/cm3), k = 1212W/cm3*Tb
At (200mT, 0.90W/cm3), k = 728W/cm3*Tb
At (100mT, 0.04W/cm3), k = 578W/cm3*Tb
Math and I aren't best friends, but I don't understand why the coefficient would change...surely that can't be correct?
Anyone know what I'm missing here?
Thanks!
floobydust:
<crickets>
I'm not fluent in going down the rabbit hole for core loss calulations, they are super extremely complicated. I find almost always the copper losses are higher.
I found one reference in "Modern Ferrite Technology" to a book by Colonel W.T. McLyman (pic related) first published 1982.
Colonel W.T. McLyman, "Magnetic Core Selection for Transformers and Inductors"
which I would hunt down and see his empirical equation discussed.
He also has: "Transformer and Inductor Design Handbook"
This is about all the help I can give.
T3sl4co1l:
Right, you can derive the exponent from the log-log plot, although when the exponent itself varies say with frequency, that's a bit awkward. You may have to iterate an approximation in that case; but that's alright, it will converge quickly, and getting within 10% of the true answer is better than the real material will likely do anyway. :)
Anyway, the offset of the log-log line is in units of loss (W/cm^3 or whatever) per flux density (tesla). The slope is beta. :)
Um, give or take a +1 to beta I think?
It may help to understand core loss in terms of the Q factor of the core itself. Constant-Q curves on the loss diagram are diagonal lines of slope = 1. A lot of materials have an exponent near there, so that Q essentially depends on frequency alone. I've some powder core datasheets that I've doodled such notes on; typical figures for a Sendust mu=60 core are:
FreqQ1k1005k8325k75100k52250k37500k22
This is a little trickier for gapped ferrite, because you can control the gap. The gap itself is lossless* (air), so as you make more of the core's effective length into airgap, the Q rises proportionally. (A terrible Q < 1 core, with air gap, can still be used quite effectively, if you need to.)
Mind, you aren't varying Bpk by varying gap -- what you are varying is the amount of magnetization required to reach Bpk. Which means for the same inductance, you need more turns, and that's where Bpk and therefore losses go down with increasing gap.
*Not counting possible eddy currents in wires near the gap.
Tim
TimNJ:
floobydust and Tim, thank you for your answers. Being a complicated subject, I'm going to sleep on your responses and see if I can come up with meaningful answers!
For now, I'd like to share my progress in determining K_fe, some problems I ran into, and maybe tips for others attempting to do the same. (Props to my dad, the math teacher, for helping.)
-----
1.) First mistake I made was reading the values off the core loss plot incorrectly. Notice on the x-axis (flux density), each tick is displayed at 100, 200, 300, etc. On the y-axis (core loss), the ticks are 101, 102, 103, etc. In between 102 and 103 are minor tick marks corresponding to 2*102, 3*102, 4*102, etc.
I however, initially read these as 202, 302. 402, etc. which is totally not the same thing! Don't do that!
-----
2.) I didn't really understand what the "core loss exponent" and "core loss coefficient" were, but now it is quite obvious. Core loss fits the power function model, that is y = a*xb. In order to determine the optimum flux density that yields the lowest total power loss, that is Pcu + Pfe (copper + core loss), we need to know exactly what the core loss function looks like. In "Fundamentals of Power Electronics", Kfe corresponds to a, and beta corresponds to b. The optimum flux density is found where Ptot = Pcu + Pfe is minimized, that is the derivative = 0, i.e. at the bottom of a parabolic function.
-----
3.) While b, the exponent, can be found by determining the slope between two points on the core loss curve, I found this method to be rather inaccurate. Since core loss is on a log-log plot, small errors tend to be amplified, especially where the tick markets get compressed. Calculating by hand, using two points, I got lots of different exponent values, between 2 and 4. This error/variance means the calculated coefficient will also be off. tl;dr Doing this entirely by eye/hand does not seem very accurate.
A more accurate method is to select a number of points on the curve (4 to 6, maybe), and fit an equation to it, using Excel. In this case, I chose 4 points, and determined the exponent to be 2.83, and the coefficient to be 27.42W/cm3*Tb. I'm much more comfortable with these values because they are quite similar to the "typical" values noted in the book. (Note, before I was getting coefficient values on the order of 500W/cm3*Tb or more.)
Navigation
[0] Message Index
[#] Next page
Go to full version