Permeability is about 2. This is hopeless for making motors and the like. Permeability of electric steel (found in transformers) is about 4000, for comparison. Permeability of inductor gap-less core material is typically about 50.
Printing the core
This part is reasonably straightforward. I chose a relatively thin and small toroidal core. Small to not need to spend too much wire for the test coil, thin to not get uneven field distribution through its section, and toroidal because it is very easy geometry to calculate.
I think it's important to choose solid concentric infill. This is because gaps across B-field influence results seriously, and need to be minimized to measure the material properties, not properties of a gapped core. Concentric infill results in a serious gap that slices the core in two rings, but gap is along B-field (not across), so it can easily be corrected for by lowering the cross-section area in calculations.
middle-circle radius R = 8 mm = 8e-3 m
core cross-section area S = 23.38 mm^2 = 23.38e-6 m^2
Making a test inductor
Simply wound a testing coil. It is quite important to distribute the turns evenly, to avoid magnetic field leak to outside and skew the results. Other than that, the more copper - the better chances of achieving saturation.
Here's my winding, I think I did a good job in terms of even-ness.
The coil has N = 124 turns. It measures 0.527 Ohms of DC resistance, and has an inductance of 21 uH.
Measuring the inductor
To measure B-H curve, I drive the coil with square wave generator. I can obtain B field by integrating voltage across ideal part of inductor, and H field by measuring current through coil.
On this circuit:
* to the left is high-current square-wave generator
* DUT is the inductor. I approximate it as an ideal winding around the core (made of superconductor) plus DC resistance of the coil. I assume that skin effect and proximity effect are negligible, since I'm driving it at quite low frequencies (around 20KHz).
* 0.1 ohm resistor is a current sense resistor
* scope probes are connected to points labeled Ch.1 and Ch.2. Scope's ground is connected to the point where the blue ground symbol is (so, the generator MUST NOT be grounded!)
So, on Channel1, I get total voltage waveform across the inductor, and on Channel2, I get an inverted value of current.
* capacitor isolates DC component of square wave from biasing the inductor. The capacitor is a large electrolytic cap. The more - the better.
This is how the circuit looks when it's built:
And this is what I get on the scope:
Note that current waveform (cyan) is inverted due to the way probes are connected.
Calculating material properties
H-field can be directly calculated from coil current by using Ampere's circuital law:
i.e. integral of H-field along a contour equals total current passing through that contour.
Choosing the contour to be the center-circle of the core, we get:
H * 2*pi*R = N * i
where R is radius of core's center circle, N is the number of turns in the coil, and i is the current through the coil. So, H field waveform is:
H(t) = N*i(t) / (2*pi*R)
Deriving B-field is quite a bit more involved, and that's what I capture voltage waveform for.
Voltage across the ideal part of the coil is tied to the rate of change of magnetic flux by law of electromagnetic induction:
dФ/dt = U(t) / N
where Ф is B-field flux through the section of the coil toroid, which equals B*S assuming B is constant throughout core section area S, and the coil is tight.
U(t) is the voltage across the ideal part of the coil, which can be calculated by subtracting resistance dropout R_DC*i(t) from the total voltage across the real coil:
U(t) = U_real(t) - i(t)*R_DC
B(t) = integral over time( (U_real(t) - i(t)*R_DC) / N ) / S + offset
'offset' is a constant, that should be chosen so that time-averaged B-field is zero.
On this plot, B and H field wavefoms are shown. H was multiplied by permeability of vacuum mu0, to bring it into a comparable scale...Read more »