"""Проверить два кольцевых зазора чашечного магнита.

Среда: встроенный MicroPython FEMM 5.
Самостоятельный аналитический контроль; полевой решатель не вызывается.
"""
import math
mu0 = 4 * math.pi * 1e-7
# Обе площади считаем в м²; наружный полюс — кольцо, а не полный круг.
sc = math.pi * .008**2
so = math.pi * (.030**2 - .022**2)
gap, mmf = .001, 300.0
reluctance = gap / mu0 * (1 / sc + 1 / so)
phi = mmf / reluctance
# Силы двух полюсов складываются при открытии диска общей координатой z.
force_z = -phi * phi / (2 * mu0) * (1 / sc + 1 / so)
print("Поток, Вб:", phi)
print("Индукции центрального и наружного полюсов, Тл:", phi / sc, phi / so)
print("Осевая сила, Н:", force_z)