Version 4
SHEET 1 1300 680
WIRE 160 176 144 176
WIRE 640 176 464 176
WIRE 688 176 640 176
WIRE 464 192 464 176
WIRE 144 208 144 176
WIRE 144 208 64 208
WIRE 416 208 144 208
WIRE 640 208 640 176
WIRE 64 288 64 208
WIRE 144 288 144 208
WIRE 464 288 464 272
WIRE 640 288 640 272
WIRE -192 304 -192 272
WIRE 64 384 64 368
WIRE 144 384 144 368
WIRE -192 416 -192 384
FLAG 144 384 0
FLAG 64 384 0
FLAG 416 256 0
FLAG 464 288 0
FLAG 640 288 0
FLAG -192 416 0
FLAG -192 272 0
FLAG 688 176 VB
IOPIN 688 176 Out
FLAG 160 176 VL
IOPIN 160 176 Out
SYMBOL ind 128 272 R0
WINDOW 3 39 66 Left 2
WINDOW 39 39 94 Left 2
SYMATTR Value Hc={Hc} Bs={Bs} Br={Br} A={A} Lm={Lm} Lg={Lg} N={N}
SYMATTR SpiceLine Rser=0
SYMATTR InstName L1
SYMBOL current 64 368 R180
WINDOW 123 0 0 Left 2
WINDOW 39 0 0 Left 2
WINDOW 0 -28 87 Left 2
WINDOW 3 39 -50 Right 2
SYMATTR InstName I1
SYMATTR Value PWL(.0 0 {1*ip} {ip} {2*ip} 0 {3*ip} {-ip} {4*ip} 0 {5*ip} {ip})
SYMBOL g 464 176 R0
WINDOW 3 20 97 Left 2
SYMATTR Value {1/(A*N)}
SYMATTR InstName G1
SYMBOL cap 624 208 R0
WINDOW 3 42 41 Left 2
SYMATTR Value 1
SYMATTR InstName C1
SYMBOL bi -192 304 R0
SYMATTR InstName BH
SYMATTR Value I=I(L1)*{N/Le}
TEXT 24 456 Left 2 !.tran 0 {ip*tfac}
TEXT -96 -136 Left 2 ;Nonlinear Inductor Test, Display of L(I) and B(I) or B(H)\nInductance: VL vs. I(L1), since I1 has always dI/dT=1A/sec\nDisplay VB vs. I(L1) or vs. I(BH) to show mag. flux density B(I) or B(H)\nL(I) = (dB(I)/dI) * A * N ,  B = mu0 * muR * (N / Le) * I
TEXT -80 56 Left 2 !.param ip =  .01
TEXT -80 24 Left 2 !.param Hc=16 Bs=.44 Br=.1 A=25u Lm=20m Lg=0 N=1000
TEXT -80 88 Left 2 !.param tfac = 5
TEXT -88 -8 Left 2 ;Parameters for user-edit:
TEXT 128 56 Left 2 ;ip : peak current & PWL time-interval
TEXT 128 88 Left 2 ;tfak: 5 ->full hysteresis loop, 1 -> initial to B(ip)
TEXT 616 24 Left 2 ;see help for L
TEXT -96 136 Left 2 !.param Le = {(Lm+Lg)**2/(Lm+Lg)} ;effective magn. length
TEXT 400 -168 Bottom 2 ;Posted to http://groups.yahoo.com/group/LTspice by <gmilder@...>
RECTANGLE Normal 784 112 -96 -32