L=mu*mu0*S*W^2/l, B=mu*mu0*H, где мю - относительная проницаемость, мю0 - таковая от вакуума (4e-7*Pi), H- напряженность поля. Далее, H=I*W/l, где I - собственно ток через обмотку, W - витки, l - средняя длина магнитной линии. Примерно, при малом зазоре, мю=l/g, где g - зазор. Сердечник в рабочей точке должен быть примерно на половине Bs - индукции насыщения.
Как-то вот так.