"""Холостой ход, установившееся короткое замыкание, нагрузка и мощность.

Среда: встроенный MicroPython FEMM 5.
Аналитический учебный пример; полевой решатель не вызывается.
Параметры, единицы величин и порядок вывода пояснены по стадиям.
"""

import math

# Действующие фазные напряжение и ЭДС, В; продольное и поперечное реактивные сопротивления, Ом.
voltage_rms, emf_rms, xd, xq = 230.0, 260.0, 8.0, 5.0
angle = math.radians(20)
# Первое слагаемое мощности связано с возбуждением, Вт.
excited_w = 3 * voltage_rms * emf_rms / xd * math.sin(angle)
# Второе слагаемое связано с магнитной несимметрией по осям d и q, Вт.
# При равенстве xd и xq оно обращается в ноль.
reluctance_w = 1.5 * voltage_rms**2 * (1 / xq - 1 / xd) * math.sin(2 * angle)
# Порядок вывода: слагаемое возбуждения; реактивное слагаемое; суммарная активная мощность, Вт.
print(excited_w, reluctance_w, excited_w + reluctance_w)