* SPICE netlist
G11 0 N004 N003 0 1
C11 N004 0 {L}
E11 N005 0 N004 0 1
V11 N005 N006 0
F11 N003 0 V11 1
Vs Vs 0 PULSE(-1 1 1m .10u .10u 2m 4m)
R12 N006 0 1
D31 N022 0 r
D32 0 N022 r
R11 Vs N003 .5
D21 N011 0 s
G21 0 N009 N008 0 1
C21 N009 0 {L}
E21 N010 0 N009 0 1
V21 N010 N011 0
F21 N008 0 V21 1
R22 N011 0 1
R21 Vs N008 .5
G31 0 N020 N019 0 1
C31 N020 0 {L}
E31 N021 0 N020 0 {.7/Is}
V31 N021 N022 0
F31 N019 0 V31 {Is/.7}
R32 N022 0 1
R31 Vs N019 .5
R1 Vs N002 .5
L1 0 N002 {L}
R3 Vs 3 .5
L3 4 3 {L}
B3 4 0 V=if(abs(I(L3))<Is,0,V(3,4)*(Ls/L-1))
R2 Vs 1 .5
L2 2 1 {Ls}
B2 2 0 V=if(abs(I(L2))>Is,0,V(1,2)*(L/Ls-1))
R9 Vs N014 .5
L9s N015 N014 {Ls}
W1 N015 N016 V1 s
V1 N016 N017 0
W2 N015 N016 V2 s
V2 0 N017 0
L9 N016 N015 {L-Ls}
R6 Vs N007 .5
L6 0 N007 Flux={L-Ls}*limit(x,{-Is},{Is})+{Ls}*x
R7 Vs N012 .5
L7 N013 N012 Flux={(L-Ls)*Is}*pwrs(tanh(pwrs(x/{Is},{n})),{1/n})
R8 Vs N018 .5
L8 0 N018 Flux={L*Is*2/Pi}*atan(x/{Is*2/Pi})
R5 Vs N001 .5
L5 0 N001 Flux={L}*x
L7s N013 0 {Ls}
R4 Vs 5 .5
L4 6 5 {Ls}
B4 6 0 V=V(5,6)*(L/Ls-1)/(1+(I(L4)/Is)**6)
.tran 0 44m 40m uic
+ Is=3p)
+ Vfwd={Is} Vrev={Is})
.param L=1m Ls=10u Is=.7
.param n=2

* LTspice schematic (.asc)
Version 4
SHEET 1 1616 1152
WIRE 504 -272 504 -304
WIRE 504 -160 504 -192
WIRE 456 -256 456 -304
WIRE 456 -304 376 -304
WIRE 312 -304 312 -272
WIRE 312 -192 312 -160
WIRE 456 -208 456 -160
WIRE 376 -160 376 -208
WIRE 376 -272 376 -304
WIRE 376 -304 312 -304
WIRE 264 -208 264 -160
WIRE 264 -304 184 -304
WIRE 264 -304 264 -256
WIRE 184 -160 184 -192
WIRE 184 -272 184 -304
WIRE -408 -448 -408 -464
WIRE 664 -304 664 -272
WIRE 664 -192 664 -160
WIRE 664 -304 632 -304
WIRE 504 -304 552 -304
WIRE 808 432 808 480
WIRE 808 544 808 576
WIRE 744 576 744 544
WIRE 744 480 744 432
WIRE 744 432 808 432
WIRE 152 -304 184 -304
WIRE 40 -304 72 -304
WIRE -408 -544 -408 -576
WIRE -408 -576 -360 -576
WIRE 664 432 744 432
WIRE 664 -64 744 -64
WIRE 744 80 744 48
WIRE 744 -32 744 -64
WIRE 504 464 504 432
WIRE 504 576 504 544
WIRE 456 480 456 432
WIRE 456 432 376 432
WIRE 312 432 312 464
WIRE 312 544 312 576
WIRE 456 528 456 576
WIRE 376 576 376 528
WIRE 376 464 376 432
WIRE 376 432 312 432
WIRE 264 528 264 576
WIRE 264 432 184 432
WIRE 264 432 264 480
WIRE 184 576 184 544
WIRE 184 464 184 432
WIRE 664 432 664 464
WIRE 664 544 664 576
WIRE 664 432 632 432
WIRE 504 432 552 432
WIRE 152 432 184 432
WIRE 40 432 72 432
WIRE 504 -32 504 -64
WIRE 504 80 504 48
WIRE 456 -16 456 -64
WIRE 456 -64 376 -64
WIRE 312 -64 312 -32
WIRE 312 48 312 80
WIRE 456 32 456 80
WIRE 376 80 376 32
WIRE 376 -32 376 -64
WIRE 376 -64 312 -64
WIRE 264 32 264 80
WIRE 264 -64 184 -64
WIRE 264 -64 264 -16
WIRE 184 80 184 48
WIRE 184 -32 184 -64
WIRE 664 -64 664 -32
WIRE 664 48 664 80
WIRE 664 -64 632 -64
WIRE 504 -64 552 -64
WIRE 152 -64 184 -64
WIRE 40 -64 72 -64
WIRE -376 -304 -344 -304
WIRE -264 -304 -232 -304
WIRE -152 -304 -120 -304
WIRE -376 176 -344 176
WIRE -264 176 -248 176
WIRE -248 176 -232 176
WIRE -120 320 -120 288
WIRE -152 176 -136 176
WIRE -120 176 -120 208
WIRE -136 176 -120 176
WIRE -376 -64 -344 -64
WIRE -264 -64 -248 -64
WIRE -248 -64 -232 -64
WIRE -120 80 -120 48
WIRE -152 -64 -136 -64
WIRE -120 -64 -120 -32
WIRE -136 -64 -120 -64
WIRE 40 176 72 176
WIRE 152 176 184 176
WIRE 280 256 280 336
WIRE 264 176 280 176
WIRE 392 256 376 256
WIRE 296 256 280 256
WIRE 392 336 376 336
WIRE 296 336 280 336
WIRE 392 336 392 256
WIRE 616 176 600 176
WIRE 280 256 280 176
WIRE 392 176 392 256
WIRE 376 176 392 176
WIRE 504 176 520 176
WIRE 392 176 424 176
WIRE 280 176 296 176
WIRE -808 -64 -776 -64
WIRE -696 -64 -664 -64
WIRE -584 -64 -552 -64
WIRE -808 176 -776 176
WIRE -696 176 -664 176
WIRE -584 176 -552 176
WIRE -808 432 -776 432
WIRE -696 432 -664 432
WIRE -584 432 -552 432
WIRE -808 -304 -776 -304
WIRE -696 -304 -664 -304
WIRE -584 -304 -552 -304
WIRE -552 256 -552 272
WIRE -376 432 -344 432
WIRE -264 432 -248 432
WIRE -248 432 -232 432
WIRE -120 576 -120 544
WIRE -152 432 -136 432
WIRE -120 432 -120 464
WIRE -136 432 -120 432
FLAG -408 -448 0
FLAG 312 -160 0
FLAG 376 -160 0
FLAG 456 -160 0
FLAG 504 -160 0
FLAG 664 -160 0
FLAG 184 -160 0
FLAG 264 -160 0
FLAG 808 576 0
FLAG 744 576 0
FLAG 40 -304 Vs
IOPIN -576 -112 In
FLAG -360 -576 Vs
IOPIN -976 -384 Out
FLAG 744 80 0
FLAG 312 576 0
FLAG 376 576 0
FLAG 456 576 0
FLAG 504 576 0
FLAG 664 576 0
FLAG 184 576 0
FLAG 264 576 0
FLAG 40 432 Vs
IOPIN -576 624 In
FLAG 312 80 0
FLAG 376 80 0
FLAG 456 80 0
FLAG 504 80 0
FLAG 664 80 0
FLAG 184 80 0
FLAG 264 80 0
FLAG 40 -64 Vs
IOPIN -576 128 In
FLAG -376 -304 Vs
IOPIN -992 -112 In
FLAG -120 -304 0
FLAG -376 176 Vs
IOPIN -992 368 In
FLAG -120 320 0
FLAG -248 176 3
FLAG -136 176 4
FLAG -376 -64 Vs
IOPIN -992 128 In
FLAG -120 80 0
FLAG -248 -64 1
FLAG -136 -64 2
FLAG 40 176 Vs
IOPIN -576 368 In
FLAG 616 176 0
FLAG -808 -64 Vs
IOPIN -1424 128 In
FLAG -552 -64 0
FLAG -808 176 Vs
IOPIN -1424 368 In
FLAG -552 272 0
FLAG -808 432 Vs
IOPIN -1424 624 In
FLAG -552 432 0
FLAG -808 -304 Vs
IOPIN -1424 -112 In
FLAG -552 -304 0
FLAG -376 432 Vs
IOPIN -992 624 In
FLAG -120 576 0
FLAG -248 432 5
FLAG -136 432 6
SYMBOL g 312 -288 R0
SYMATTR InstName G11
SYMATTR Value 1
SYMBOL cap 360 -272 R0
SYMATTR InstName C11
SYMATTR Value {L}
SYMBOL e 504 -288 R0
SYMATTR InstName E11
SYMATTR Value 1
SYMBOL voltage 536 -304 R270
SYMATTR InstName V11
SYMATTR Value 0
SYMBOL f 184 -272 R0
SYMATTR InstName F11
SYMATTR Value V11
SYMATTR Value2 1
SYMBOL voltage -408 -560 R0
SYMATTR Value PULSE(-1 1 1m .10u .10u 2m 4m)
SYMATTR InstName Vs
SYMBOL res 648 -288 R0
SYMATTR InstName R12
SYMATTR Value 1
SYMBOL diode 728 480 R0
SYMATTR InstName D31
SYMATTR Value r
SYMBOL diode 792 544 M180
SYMATTR InstName D32
SYMATTR Value r
SYMBOL res 56 -288 R270
SYMATTR InstName R11
SYMATTR Value .5
SYMBOL Misc\\xvaristor 728 -48 R0
SYMATTR InstName D21
SYMATTR Value s
SYMATTR Prefix D
SYMBOL g 312 -48 R0
SYMATTR InstName G21
SYMATTR Value 1
SYMBOL cap 360 -32 R0
SYMATTR InstName C21
SYMATTR Value {L}
SYMBOL e 504 -48 R0
SYMATTR InstName E21
SYMATTR Value 1
SYMBOL voltage 536 -64 R270
SYMATTR InstName V21
SYMATTR Value 0
SYMBOL f 184 -32 R0
SYMATTR InstName F21
SYMATTR Value V21
SYMATTR Value2 1
SYMBOL res 648 -48 R0
SYMATTR InstName R22
SYMATTR Value 1
SYMBOL res 56 -48 R270
SYMATTR InstName R21
SYMATTR Value .5
SYMBOL g 312 448 R0
SYMATTR InstName G31
SYMATTR Value 1
SYMBOL cap 360 464 R0
SYMATTR InstName C31
SYMATTR Value {L}
SYMBOL e 504 448 R0
SYMATTR InstName E31
SYMATTR Value {.7/Is}
SYMBOL voltage 536 432 R270
SYMATTR InstName V31
SYMATTR Value 0
SYMBOL f 184 464 R0
SYMATTR InstName F31
SYMATTR Value V31
SYMATTR Value2 {Is/.7}
SYMBOL res 648 448 R0
SYMATTR InstName R32
SYMATTR Value 1
SYMBOL res 56 448 R270
SYMATTR InstName R31
SYMATTR Value .5
SYMBOL res -360 -288 R270
SYMATTR InstName R1
SYMATTR Value .5
SYMBOL ind2 -136 -288 M270
SYMATTR InstName L1
SYMATTR Value {L}
SYMATTR Type ind
SYMBOL res -360 192 R270
SYMATTR InstName R3
SYMATTR Value .5
SYMBOL ind2 -136 192 M270
SYMATTR InstName L3
SYMATTR Value {L}
SYMATTR Type ind
SYMBOL bv -120 192 R0
SYMATTR InstName B3
SYMATTR Value V=if(abs(I(L3))<Is,0,V(3,4)*(Ls/L-1))
SYMBOL res -360 -48 R270
SYMATTR InstName R2
SYMATTR Value .5
SYMBOL ind2 -136 -48 M270
SYMATTR InstName L2
SYMATTR Value {Ls}
SYMATTR Type ind
SYMBOL bv -120 -48 R0
SYMATTR InstName B2
SYMATTR Value V=if(abs(I(L2))>Is,0,V(1,2)*(L/Ls-1))
SYMBOL res 56 192 R270
SYMATTR InstName R9
SYMATTR Value .5
SYMBOL ind2 280 192 M270
SYMATTR InstName L9s
SYMATTR Value {Ls}
SYMATTR Type ind
SYMBOL csw 296 256 R270
SYMATTR InstName W1
SYMATTR Value V1 s
SYMBOL voltage 408 176 M90
SYMATTR InstName V1
SYMATTR Value 0
SYMBOL csw 296 336 R270
SYMATTR InstName W2
SYMATTR Value V2 s
SYMBOL voltage 616 176 R90
SYMATTR InstName V2
SYMATTR Value 0
SYMBOL ind2 392 192 M270
SYMATTR InstName L9
SYMATTR Value {L-Ls}
SYMATTR Type ind
SYMBOL res -792 -48 R270
SYMATTR InstName R6
SYMATTR Value .5
SYMBOL ind2 -568 -48 M270
SYMATTR InstName L6
SYMATTR Value Flux={L-Ls}*limit(x,{-Is},{Is})+{Ls}*x
SYMATTR Type ind
SYMBOL res -792 192 R270
SYMATTR InstName R7
SYMATTR Value .5
SYMBOL ind2 -568 192 M270
SYMATTR InstName L7
SYMATTR Value Flux={(L-Ls)*Is}*pwrs(tanh(pwrs(x/{Is},{n})),{1/n})
SYMATTR Type ind
SYMBOL res -792 448 R270
SYMATTR InstName R8
SYMATTR Value .5
SYMBOL ind2 -568 448 M270
SYMATTR InstName L8
SYMATTR Value Flux={L*Is*2/Pi}*atan(x/{Is*2/Pi})
SYMATTR Type ind
SYMBOL res -792 -288 R270
SYMATTR InstName R5
SYMATTR Value .5
SYMBOL ind2 -568 -288 M270
SYMATTR InstName L5
SYMATTR Value Flux={L}*x
SYMATTR Type ind
SYMBOL ind -568 160 R0
SYMATTR InstName L7s
SYMATTR Value {Ls}
SYMBOL res -360 448 R270
SYMATTR InstName R4
SYMATTR Value .5
SYMBOL ind2 -136 448 M270
SYMATTR InstName L4
SYMATTR Value {Ls}
SYMATTR Type ind
SYMBOL bv -120 448 R0
SYMATTR InstName B4
SYMATTR Value V=V(5,6)*(L/Ls-1)/(1+(I(L4)/Is)**6)