Оглавление

Введение

В современной авиации, возобновляемой энергетике и электротранспорте происходит переход от традиционных трёхфазных электрических машин к многофазным системам с пятью, шестью и более фазами. Это не академическая экзотика — многофазные генераторы уже используются на борту современных широкофюзеляжных самолётов, в мегаваттных ветроустановках морских ветропарков, в электрических силовых установках БПЛА большой грузоподъёмности.

Однако при моделировании таких систем разработчик сталкивается с неожиданной проблемой: классические инструменты MATLAB/Simulink не содержат готовых блоков многофазных машин. Нельзя просто взять стандартную модель трёхфазного двигателя и «растянуть» её на произвольное количество фаз. В этой статье мы разберём математический аппарат моделирования многофазных электрических машин и покажем, почему эта задача оказывается связанной с глубокими результатами теории чисел XVII-XIX веков — теоремой Гаусса-Венцеля о построении правильных многоугольников и простыми числами Ферма.

Многофазные машины в современной технике

Прежде чем погружаться в математику, важно понять, зачем вообще нужны многофазные системы. Трёхфазные машины работают прекрасно уже больше века — стандарты отработаны, компоненты доступны, инженеры имеют огромный опыт проектирования и эксплуатации. Что заставляет индустрию усложнять системы, переходя к большему числу фаз?

Авиационные генераторы: когда надёжность критична

Рассмотрим бортовой генератор современного широкофюзеляжного самолёта. Типичная мощность такого генератора составляет 250 кВА при напряжении 230 В переменного тока. Для трёхфазной системы при таких параметрах фазный ток достигает 360 А. Это предъявляет жёсткие требования к силовым компонентам: необходимы кабели большого сечения, мощные выпрямительные диоды или управляемые выпрямители на базе IGBT-модулей с высоким номинальным током, развитая система охлаждения. В авиации, где каждый килограмм массы критичен, это серьёзная проблема.

Переход к шестифазной конфигурации снижает фазный ток до 210 А — уменьшение на 40%. Это позволяет использовать более лёгкие кабели, стандартные силовые полупроводниковые модули, упростить систему отвода тепла. Однако экономия веса — не главное преимущество.

Ключевое достоинство многофазных генераторов — возможность разделения обмоток на независимые группы. Шестифазная машина может быть выполнена как две независимые трёхфазные обмотки, сдвинутые на 30 электрических градусов (рис. 1). При выходе из строя одной обмотки (короткое замыкание, обрыв, повреждение изоляции) система управления обнаруживает неисправность, отключает повреждённую группу и продолжает работу на оставшейся. Мощность снижается до 60-70% номинальной, но критически важные системы самолёта продолжают получать питание до безопасной посадки.

Рисунок 1 - Обмотки шестифазной электрической машины
Рисунок 1 - Обмотки шестифазной электрической машины

Это принципиальное отличие от трёхфазной системы. В трёхфазном генераторе отказ одной фазы приводит к резкой несимметрии токов, опасным перегрузкам оставшихся фаз, высокому риску каскадного отказа. Аварийная ситуация развивается быстро, времени на реакцию мало. В многофазной системе отказ одной фазы из шести — штатная ситуация, предусмотренная алгоритмами управления. Резервирование встроено в саму архитектуру, без дублирования оборудования.

 Аналогичная логика применима к электрическим силовым установкам беспилотных летательных аппаратов большой грузоподъёмности. БПЛА, предназначенные для доставки грузов или длительного патрулирования, всё чаще оснащаются многофазными электродвигателями именно из соображений отказоустойчивости. Потеря одной фазы не должна приводить к падению аппарата — достаточно обеспечить возможность управляемой посадки или возврата на базу с пониженной мощностью.

Ветрогенераторы: когда мощности растут

В современных морских ветропарках устанавливаются турбины мощностью 10-15 МВт. При таких масштабах трёхфазная конфигурация создаёт серьёзные инженерные проблемы. Фазные токи достигают тысяч ампер, что требует параллельного соединения силовых полупроводниковых модулей с тщательной балансировкой токов между ними, массивных токопроводов, сложных систем охлаждения.

Многофазная архитектура естественным образом распределяет мощность между большим числом фаз. Девятифазный генератор на 15 МВт имеет фазные токи примерно втрое меньше, чем эквивалентный трёхфазный. Это позволяет использовать стандартные промышленные силовые модули без параллельного включения, упрощает конструкцию, повышает надёжность.

В условиях морской установки, где техническое обслуживание затруднено и дорого, а простои ведут к значительным финансовым потерям, повышение надёжности имеет прямое экономическое выражение. Отказоустойчивость многофазной системы означает, что турбина может продолжать работать при повреждении одной или даже нескольких фаз до планового визита сервисной бригады, вместо аварийной остановки с вызовом специального судна.

Снижение пульсаций момента

Трёхфазная электрическая машина создаёт электромагнитный момент с пульсациями шестой гармоники. Для большинства промышленных применений эти пульсации на уровне 2-3% от среднего момента не критичны и практически не влияют на работу привода. Однако существуют применения, где даже такие небольшие пульсации недопустимы.

В прецизионных металлообрабатывающих станках пульсации момента главного привода передаются через механическую систему на инструмент и обрабатываемую деталь, вызывая микровибрации, которые ухудшают качество обработки. В медицинском оборудовании — томографах, центрифугах — вибрации создают помехи для измерительных систем. В научных установках типа синхротронов или прецизионных телескопов требования к плавности вращения ещё жёстче.

Пятифазная машина создаёт пульсации десятой гармоники с амплитудой в 2.5 раза меньшей, чем трёхфазная. Шестифазная — двенадцатой гармоники с амплитудой втрое меньшей. Кроме того, высокочастотные пульсации лучше демпфируются инерцией механической системы. Пульсация на частоте 1200 Гц гасится упругостью валов и муфт гораздо эффективнее, чем на 600 Гц. В результате акустический шум также снижается, что важно для применений, где уровень шума регламентируется санитарными нормами или требованиями комфорта.

Итак, практические преимущества многофазных систем очевидны и убедительны. Но как их моделировать? Здесь начинаются сложности.

Фундаментальная проблема моделирования

При попытке смоделировать многофазную электрическую машину в стандартных инструментах типа MATLAB/Simulink разработчик обнаруживает отсутствие готовых решений. Библиотека Simscape Electrical содержит блоки трёхфазных асинхронных машин, синхронных машин с постоянными магнитами, синхронных машин с электромагнитным возбуждением, но ничего для систем с количеством фаз больше трёх.

Можно попытаться построить модель «с нуля», описывая каждую обмотку отдельно. В такой модели необходимо прописать уравнения для токов и напряжений каждой обмотки, учесть взаимоиндукцию между всеми обмотками статора, взаимодействие с ротором, зависимость магнитных потоков от угла поворота. Для n-фазной машины получается система из нескольких десятков связанных нелинейных дифференциальных уравнений с периодически изменяющимися коэффициентами.

Такая модель симулируется крайне медленно. Численное интегрирование системы с быстро меняющимися коэффициентами требует малого шага интеграции для обеспечения устойчивости и точности. Модель занимает огромную схему, трудна в отладке, практически непригодна для синтеза систем управления. Время расчета одной секунды реального времени может занимать минуты или даже часы процессорного времени.

Однако для трёхфазных машин давно известен элегантный способ обойти эту проблему — координатные преобразования, которые радикально упрощают математическое описание. Эти методы подробно описаны в доступной форме в статье и «Векторное управление для асинхронного электродвигателя на пальцах», где принципы объясняются наглядно без излишнего формализма.

Суть подхода состоит в том, что система с n фазными переменными (токи, напряжения, потокосцепления) преобразуется в эквивалентную двухкоординатную систему. Первое преобразование — трансформация Кларка — переводит n фазных величин в две ортогональные компоненты α и β в неподвижной системе координат. Второе — трансформация Парка — переводит эти компоненты во вращающуюся систему координат d-q, синхронизированную с ротором машины или с вектором потокосцепления статора.

В результате этой двойной трансформации система дифференциальных уравнений с периодически изменяющимися коэффициентами превращается в систему с постоянными коэффициентами. Более того, при установившихся режимах переменные в системе d-q становятся постоянными величинами. Вместо n синусоидальных токов, зависящих друг от друга сложным образом, мы получаем два постоянных тока — один отвечает за магнитный поток (как ток возбуждения в машине постоянного тока), другой за момент (как ток якоря). Управление машиной переменного тока становится столь же простым, как управление машиной постоянного тока — это был настоящий прорыв в теории электропривода.

Для трёхфазных систем эта методология отработана десятилетиями. Существуют готовые формулы, стандартные программные библиотеки, тысячи реализаций в промышленных контроллерах. Но можно ли обобщить этот подход на произвольное число фаз? Оказывается, здесь возникает тонкая математическая проблема, уходящая корнями в классическую геометрию и теорию чисел.

Математическая модель многофазной машины

Построение математической модели начинается с формулировки основных предположений. Эти предположения одновременно упрощают задачу до разрешимого уровня и сохраняют все существенные физические эффекты.

Первое предположение — полная симметрия системы. Все n обмоток статора идентичны по конструкции: одинаковые активные сопротивления R, индуктивности рассеяния Lσ, число витков. Обмотки расположены равномерно по окружности статора, так что электрический угол Θn между магнитными осями соседних обмоток составляет точно 2π/n радиан (рис. 2). Это не произвольное математическое упрощение, а требование к реальной конструкции машины. Любая асимметрия приведёт к паразитным моментам, дополнительным потерям, вибрациям и ухудшению характеристик.

Рисунок 2 - Диаграмма n-фазной электрической машины
Рисунок 2 - Диаграмма n-фазной электрической машины

Второе предположение — соединение обмоток звездой с изолированной нейтральной точкой. Это даёт фундаментальное ограничение: сумма мгновенных значений всех фазных токов равна нулю:

i₁(t) + i₂(t) + ... + i_n(t) = 0

Это не специальное условие, а прямое следствие первого закона Кирхгофа. Если нейтраль изолирована, ток не может втекать или вытекать из узла соединения обмоток — он может только циркулировать между фазами. Данное ограничение уменьшает число независимых переменных с n до (n-1), что критично для возможности перехода к двухкоординатному представлению. Без этого ограничения потребовалось бы три координаты вместо двух.

Третье предположение — представление распределённых обмоток как сосредоточенных. В реальной машине обмотка каждой фазы занимает несколько пазов статора и создаёт пространственно распределённое магнитное поле со сложной конфигурацией. Однако для задач моделирования динамики достаточно заменить распределённую обмотку эквивалентной сосредоточенной с определённой магнитной осью — направлением максимума создаваемого поля. Это стандартное приближение теории электрических машин, точность которого подтверждена многолетней практикой сопоставления расчётов и экспериментов.

Четвёртое предположение — линейность магнитной системы в первом приближении. Реальная машина обладает насыщением магнитопровода, гистерезисом, различными нелинейными эффектами. Однако для базовой модели эти явления не учитываются, что позволяет применять принцип суперпозиции: результирующее магнитное поле есть векторная сумма полей от всех обмоток и от постоянных магнитов (если они есть). Учёт нелинейностей может быть добавлен на следующем этапе уточнения модели, но фундаментальная структура преобразований остаётся той же.

При сделанных предположениях ключевая идея состоит в переходе от n скалярных фазных переменных к векторному представлению. Для мгновенных значений фазных токов i₁(t), i₂(t), ..., in(t) определяется комплексный пространственный вектор тока:

i(t) = Σ i_k(t)·exp(j·2π(k-1)/n)

где j — мнимая единица, exp — комплексная экспонента. Геометрический смысл прозрачен: каждый фазный ток представляется вектором в комплексной плоскости, направленным под углом, соответствующим пространственному положению k-й обмотки. Результирующий вектор — их векторная сумма.

Если токи образуют симметричную систему синусоидальных величин:

i_k(t) = I·cos(ωt - 2π(k-1)/n)

то после подстановки в формулу и применения тригонометрических тождеств получается, что вектор имеет постоянную амплитуду и равномерно вращается с угловой скоростью ω. Это и есть математическое описание вращающегося магнитного поля, которое является основой работы любой машины переменного тока.

Однако для численного моделирования и реализации в микроконтроллерах удобнее работать не с комплексными числами, а с вещественными величинами. Комплексный вектор раскладывается на две ортогональные вещественные компоненты — проекции на оси α и β неподвижной системы координат:

i = i_α + j·i_β

где iα и iβ — вещественные функции времени. Переход от n фазных токов к двум компонентам (iα и iβ) и есть суть трансформации Кларка. Вместо n переменных получаем две, что радикально упрощает математическое описание и делает возможным эффективное численное моделирование. Система из n дифференциальных уравнений сводится к системе из двух уравнений.

Обобщенная трансформация Кларка

Прямая трансформация Кларка для n-фазной системы записывается в матричной форме. Вектор двух компонент [iα, iβ]T получается умножением матрицы преобразования размерности 2×n на вектор-столбец фазных токов [i₁, i₂, ..., in]T.

Матрица прямого преобразования имеет вид:

Структура матрицы отражает геометрию системы. Первая строка содержит косинусы углов, соответствующих пространственному положению каждой обмотки — это проекции единичных векторов фаз на ось α. Вторая строка содержит синусы тех же углов со знаком минус — это проекции на ось β. Коэффициент 2/n обеспечивает нормировку, при которой сохраняется мгновенная мощность — важнейшее физическое требование.

Обратная трансформация восстанавливает n фазных токов из двух компонент (iα, iβ). Матрица обратного преобразования имеет размерность n×2:

Эти матрицы обладают рядом критически важных свойств. Первое — сохранение мгновенной мощности. Сумма произведений фазных токов на фазные напряжения равна сумме произведений компонент α-β на соответствующие напряжения:

P_n = i₁u₁ + i₂u₂ + ... + i_n·u_n = i_α·u_α + i_β·u_β = P_{αβ}

Это следствие унитарности преобразования. Физически это означает, что трансформация координат не вносит и не изымает энергию из системы — она лишь меняет способ математического представления. Любое нарушение этого свойства привело бы к нефизичным результатам моделирования.

Второе свойство — взаимная обратность в определённом смысле. Последовательное применение прямой и обратной трансформации возвращает исходный вектор без искажений (с точностью до компоненты нулевой последовательности, которая отсутствует при изолированной нейтрали). Математически:

T_{n/αβ} · T_{αβ/n} · [i₁, i₂, ..., i_n]^T = [i₁, i₂, ..., i_n]^T

Это критично для систем управления, где цикл «измерение токов → прямая трансформация → вычисления → обратная трансформация → формирование напряжений» выполняется тысячи раз в секунду. Если преобразования не точно обратны, ошибки будут накапливаться, приводя к деградации качества управления.

Третье свойство — инвариантность относительно вращения. Вращение вектора в n-фазной системе соответствует вращению в системе α-β с той же угловой скоростью и фазой. Физически это означает, что вращающееся магнитное поле корректно представляется в новой системе координат без искажения его динамики.

Казалось бы, всё готово для моделирования: формулы записаны, свойства проверены, можно реализовывать в коде и моделировать любое количество фаз. Однако именно здесь проявляется тонкая проблема.

Критический вопрос: точность коэффициентов

Матрицы трансформации содержат тригонометрические коэффициенты вида cos(2πk/n) и sin(2πk/n). Для реализации в программе эти числа нужно представить в памяти компьютера. Возникает вопрос: можно ли записать их в виде точных формул, или придётся использовать приближённые численные значения?

Для n = 3 имеем углы 120° и 240°. Соответствующие функции:

cos(120^o)=-1/2; sin(120^o)=(√3)/2

Это точные значения. Квадратный корень из трёх — иррациональное число, но при реализации достаточно один раз вычислить √(3) = 1.7320508... с максимальной точностью (обычно 15-16 значащих цифр для типа double) и использовать эту константу. Прямая и обратная трансформации будут точно взаимообратными с точностью до машинного эпсилон (около 10-16).

Для n = 4 ещё проще — углы 90°, 180°, 270° дают:

cos(90^o)=0; sin(90^o)=1; cos(180^o)=-1; sin(180^o)=0; cos(270^o)=0; sin(270^o)=-1

Только целые числа 0, ±1. Никаких иррациональностей, максимальная точность представления.

Для n = 5 углы 72° и 144° связаны с золотым сечением φ = (1+√5)/2:

cos(72^o)=(√5-1)/4; sin(72^o)=√(10+2√5)/4cos(144^o)=-(√5+1)/4; sin(144^o)=√(10-2√5)/4

Формулы сложнее, содержат вложенные радикалы, но это всё ещё точные выражения через квадратные корни. Можно предвычислить эти константы один раз при инициализации программы и использовать их с полной точностью.

А что с n = 7? Углы составляют приблизительно 51.43°, 102.86° и так далее. Величина cos(2π/7) не выражается через квадратные корни. Численно:

cos(2π/7)≈0.6234898018587335305...

Это бесконечная непериодическая десятичная дробь. При реализации приходится использовать приближение с конечным числом знаков.

Казалось бы, современные компьютеры имеют высокую точность представления вещественных чисел (15-16 значащих цифр для double) — разве этого недостаточно? Проблема в том, что при моделировании системы управления с обратной связью преобразования выполняются в цикле с высокой частотой. Если частота управления 10 кГц, за одну секунду выполняется 10 000 циклов «прямая трансформация → вычисления → обратная трансформация». За минуту — 600 000 циклов.

Когда коэффициенты приближённые, прямая и обратная трансформации не являются точно взаимообратными. При каждом цикле возникает малая ошибка. Для систем с точными коэффициентами (n = 3, 4, 5, 6) эта ошибка на уровне машинного эпсилон (10-16) и не накапливается. Для систем с приближёнными коэффициентами (n = 7, 9, 11...) ошибка на уровне 10-15 — 10-16 на каждом цикле медленно накапливается, приводя к «дрейфу» системы координат, появлению паразитных гармоник, нарушению энергетического баланса.

Возникает фундаментальный вопрос: для каких n тригонометрические функции cos(2πk/n) и sin(2πk/n) представимы точными формулами через квадратные корни? Это не праздный вопрос — от ответа зависит, какие многофазные системы можно моделировать с высокой точностью, а какие потребуют специальных мер для компенсации накопления ошибок.

Теорема Гаусса-Венцеля: неожиданная связь с геометрией

Ответ на поставленный вопрос даёт замечательная теорема из теории чисел, установленная на рубеже XVIII-XIX веков. История этого открытия начинается в Древней Греции и связана с классической задачей построения правильных многоугольников.

Античные математики умели строить циркулем и линейкой правильные треугольник (n=3), четырёхугольник (n=4), пятиугольник (n=5) и шестиугольник (n=6). Построение пятиугольника, связанное с золотым сечением, считалось особенно красивым и было символом пифагорейской школы. Но правильный семиугольник (n=7) не поддавался. Многие поколения геометров пытались найти построение, но безуспешно. Постепенно сформировалось предположение, что семиугольник вообще невозможно построить циркулем и линейкой, однако строгого доказательства не было на протяжении двух тысячелетий.

Прорыв произошёл 30 марта 1796 года. Девятнадцатилетний студент Карл Фридрих Гаусс, будущий «король математиков», записал в свой дневник открытие, которое потрясло математический мир: правильный 17-угольник можно построить циркулем и линейкой.

Рисунок 3 - Карл Фридрих Гаусс
Рисунок 3 - Карл Фридрих Гаусс

Это было совершенно неожиданно. Семнадцать — простое число, причём достаточно большое, и никто не предполагал, что такое построение возможно. Более того, Гаусс не ограничился единичным результатом. Он сформулировал общий критерий, точно описывающий, какие правильные n-угольники можно построить циркулем и линейкой, а какие нельзя принципиально. Это открытие произвело такое впечатление на юного Гаусса, что он решил посвятить свою жизнь математике, а не филологии, как планировал ранее. По легенде, Гаусс просил выгравировать правильный 17-угольник на своём надгробии, но каменотёс отказался, сочтя задачу слишком сложной.

Теорема Гаусса-Венцеля: правильный n-угольник можно построить циркулем и линейкой тогда и только тогда, когда

n = 2^k · p₁ · p₂ · ... · p_s

где k ≥ 0 — целое число, а p₁, p₂, ..., ps — различные простые числа Ферма.

Числа Ферма определяются формулой:

F_m = 2^{2^m} + 1

где m = 0, 1, 2, 3, ... — неотрицательные целые числа.

Вычислим первые числа Ферма:

F₀ = 2^{2^0} + 1 = 2^1 + 1 = 3 - простоеF₁ = 2^{2^1} + 1 = 2^2 + 1 = 5 - простое F₂ = 2^{2^2} + 1 = 2^4 + 1 = 17 - простоеF₃ = 2^{2^3} + 1 = 2^8 + 1 = 257 - простоеF₄ = 2^{2^4} + 1 = 2^{16} + 1 = 65537 - простоеF₅ = 2^{2^5} + 1 = 2^{32} + 1 = 4294967297 = 641 × 6700417 - составное

На сегодняшний день известно только пять простых чисел Ферма: 3, 5, 17, 257, 65537. Числа F₅, F₆, F₇ и все проверенные дальше (а проверены числа до F₃₂) оказались составными. Остаётся открытым вопрос, существуют ли другие простые числа Ферма — это одна из знаменитых нерешённых проблем теории чисел. Большинство математиков склоняется к гипотезе, что их больше нет, но доказательства не существует.

Почему эта теорема о построении многоугольников связана с нашей задачей моделирования электрических машин? Связь глубокая и неочевидная. Операции, допустимые при построении циркулем и линейкой — проведение прямой через две точки, проведение окружности с центром в точке через другую точку, нахождение точек пересечения — алгебраически соответствуют четырём арифметическим операциям (сложение, вычитание, умножение, деление) и извлечению квадратного корня. Никакие другие операции при построении циркулем и линейкой недопустимы.

Следовательно, координаты точки, построенной циркулем и линейкой начиная с рациональных координат, выражаются формулой, содержащей только рациональные числа, четыре арифметические операции и квадратные корни (возможно, вложенные). Такие числа называются конструктивными.

Вершины правильного n-угольника, вписанного в единичную окружность, имеют координаты (cos(2πk/n), sin(2πk/n)) для k = 0, 1, ..., n-1. Если n-угольник можно построить циркулем и линейкой, то эти координаты — конструктивные числа, представимые через радикалы. Если нельзя — то эти числа не представимы формулами, содержащими только квадратные корни, как бы сложно мы их ни вкладывали друг в друга.

Рисунок 4 - построение многоугольника по теореме Гаусса-Венцеля
Рисунок 4 - построение многоугольника по теореме Гаусса-Венцеля

Таким образом, теорема Гаусса-Венцеля точно указывает, для каких n коэффициенты трансформации Кларка могут быть записаны в виде точных формул через квадратные корни. Для этих n реализация преобразований в программе будет иметь максимальную точность, ограниченную только разрядностью вещественной арифметики. Для остальных n придётся использовать приближённые численные значения со всеми вытекающими проблемами накопления ошибок.

Практически применимые конфигурации

Применяя теорему Гаусса-Венцеля к задаче моделирования многофазных машин, составим таблицу для практически интересных значений n ≤ 12:

n

Разложение

Точные формулы

Применимость

3

20·3

✓

Стандарт

4

22

✓

Отличная

5

20·5

✓

Отличная

6

21·3

✓

Отличная

7

7 (не Ферма)

✗

Проблематична

8

23

✓

Отличная

9

32

✗

Проблематична

10

21·5

✓

Хорошая

11

11 (не Ферма)

✗

Проблематична

12

22·3

✓

Хорошая

Из таблицы видно, что наиболее перспективными для практического применения являются конфигурации с n = 3, 4, 5, 6, 8, 10, 12 фаз. Именно эти системы можно моделировать с максимальной точностью без специальных мер по компенсации накопления ошибок.

Системы с n = 7, 9, 11, 13, 14, 15 фаз проблематичны с точки зрения моделирования. Это не означает, что их нельзя моделировать вообще — численные методы достаточно развиты для работы с приближёнными коэффициентами. Однако требуются дополнительные меры: периодическая ренормализация для предотвращения дрейфа, более аккуратный выбор численных методов интегрирования, контроль энергетического баланса. Всё это усложняет разработку и снижает вычислительную эффективность модели.

Интересно, что эти математические ограничения хорошо согласуются с практикой. В промышленности действительно наиболее распространены системы с 3, 5, 6 и 12 фазами. Семифазные системы практически не встречаются, хотя с точки зрения распределения мощности они были бы вполне логичны. Возможно, трудности моделирования и управления сыграли свою роль в выборе конфигураций.

Реализация в системе моделирования Engee

Теперь рассмотрим, как всё это реализуется в практической задаче моделирования многофазной машины в современной российской системе моделирования Engee.

Структура модели многофазного синхронного двигателя

Разработанная модель представляет собой полнофункциональную систему n-фазного синхронного двигателя с постоянными магнитами (СДПМ) и векторным управлением (рис. 5). Ключевая особенность модели — параметрическая настройка числа фаз через скрипт инициализации, что позволяет исследовать различные конфигурации без перестроения схемы.

Рисунок 5 - Модель системы векторного управления n-фазного СДПМ в Engee
Рисунок 5 - Модель системы векторного управления n-фазного СДПМ в Engee

Модель включает в себя следующие основные компоненты:

Блок n-фазного синхронного двигателя с постоянными магнитами (рис. 6), математическое описание которого реализовано в системе координат d-q. Это ядро модели не зависит от числа фаз и описывается универсальными дифференциальными уравнениями электромагнитных и механических процессов.

Рисунок 6 - Структура блока СДПМ
Рисунок 6 - Структура блока СДПМ

Блоки прямого и обратного преобразования Кларка, которые осуществляют переход между n-фазной системой координат и двухкоординатной системой α-β. Эти блоки используют предвычисленные матрицы трансформации, которые рассчитываются в скрипте инициализации в зависимости от заданного числа фаз.

Блоки прямого и обратного преобразования Парка, обеспечивающие переход между неподвижной системой координат α-β и вращающейся системой d-q, синхронизированной с ротором машины. Эти преобразования зависят от текущего электрического угла поворота ротора.

Систему векторного управления, построенную по каскадному принципу. Внешний контур содержит ПИ-регулятор скорости, который формирует задание на момент (ток по оси q). Внутренние контуры включают два независимых ПИ-регулятора токов для осей d и q, обеспечивающих точную отработку заданных токов.

n-фазный силовой инвертор с широтно-импульсной модуляцией, который преобразует n фазных напряжений, сформированных системой управления, в последовательности импульсов для управления силовыми транзисторами.

Датчики фазных токов и положения ротора, обеспечивающие обратные связи для системы управления.

Ядро модели: уравнения в осях d-q

Центральный блок модели — математическое описание синхронной машины в системе координат d-q. Это ядро не зависит от числа фаз и описывается одинаковыми уравнениями для любой конфигурации.

Уравнения напряжения обмоток связывают напряжения и токи по продольной (d) и поперечной (q) осям:

u_d = R_s·i_d + L_d·di_d/dt - ω_e·L_q·i_qu_q = R_s·i_q + L_q·di_q/dt + ω_e·(L_d·i_d + ψ_{PM})

Электромагнитный момент определяется выражением:

T_{em} = (3/2)·p·[ψ_{PM}·i_q + (L_d - L_q)·i_d·i_q]

Механическое уравнение движения описывает динамику ротора:

J·dω_m/dt = T_{em} - T_{load} - B·ω_m

где Rs — активное сопротивление обмотки фазы, Ld и Lq — индуктивности по продольной и поперечной осям, ψPM — потокосцепление от постоянных магнитов, p — число пар полюсов, ωe = p·ωm — электрическая угловая скорость, J — момент инерции ротора, B — коэффициент вязкого трения.

Эти уравнения универсальны! Переход от трёхфазной машины к пятифазной или двенадцатифазной не требует изменения этого блока — меняются только трансформации на входе и выходе.

Предварительное вычисление матриц трансформации в скрипте

Ключевая особенность реализации в Engee — матрицы трансформации Кларка предварительно вычисляются в скрипте инициализации один раз перед запуском модели. Это обеспечивает максимальную производительность симуляции.

Скрипт инициализации (код приведен ниже) выполняет следующие функции:

Задание числа фаз nphases как основного параметра конфигурации. Изменение этого единственного параметра автоматически перенастраивает всю модель на другое число фаз без необходимости модификации схемы.

Вычисление матрицы прямого преобразования Кларка Tαβ/n размерностью 2×n. Эта матрица содержит косинусы и синусы углов, соответствующих пространственному положению каждой обмотки, с нормировочным коэффициентом 2/n для сохранения мощности.

Вычисление матрицы обратного преобразования Кларка Tn/αβ размерностью n×2. Эта матрица восстанавливает n фазных величин из двух ортогональных компонент α и β.

Проверку взаимной обратности матриц. Произведение Tn/αβ · Tαβ/n должно давать единичную матрицу с точностью до машинного эпсилон, что подтверждает корректность вычислений.

Задание параметров машины: электрических (сопротивления, индуктивности, потокосцепление магнитов), механических (момент инерции, трение, число пар полюсов) и номинальных режимов работы.

Настройку параметров силового инвертора: напряжение звена постоянного тока, частота ШИМ-модуляции.

Настройку коэффициентов регуляторов: пропорциональные и интегральные коэффициенты для регулятора скорости и регуляторов токов, а также ограничения по токам и напряжениям.

using LinearAlgebra
function clarke_transform_nphase(n_phases::Int; verbose=true)
    # Инициализация матриц
    T_Clarke = zeros(2, n_phases)
    T_Inverse = zeros(n_phases, 2)  
    # Заполнение матриц
    for k in 1:n_phases
        angle = 2π*(k-1)/n_phases
        T_Clarke[1, k] = cos(angle)
        T_Clarke[2, k] = sin(angle)
        T_Inverse[k, 1] = cos(angle)
        T_Inverse[k, 2] = sin(angle)
    end
    # Нормализация
    T_Clarke = (2/n_phases) * T_Clarke
    if verbose
        is_constructible = check_gauss_wantzel(n_phases)
        println("\n========================================")
        println("Конфигурация: $n_phases-фазная система")
        println("========================================")
        println(is_constructible ? "✓ Удовлетворяет теореме Гаусса-Ванцеля" : 
                                   "✗ НЕ удовлетворяет теореме Гаусса-Ванцеля")
        println("Угол между фазами: $(round(360/n_phases, digits=2))°")
        println("\nНорма ошибки ортогональности: $(norm(T_Clarke * T_Inverse - I(2)))")
    end 
    return T_Clarke, T_Inverse
end
function check_gauss_wantzel(n::Int)
    fermat_primes = [3, 5, 17, 257, 65537]
    temp_n = n
    while mod(temp_n, 2) == 0
        temp_n = div(temp_n, 2)
    end
    temp_n == 1 && return true
    for fp in fermat_primes
        if mod(temp_n, fp) == 0
            temp_n = div(temp_n, fp)
            mod(temp_n, fp) == 0 && return false
        end
    end 
    return temp_n == 1
end

# Параметры синхронного генератора с постоянными магнитами
n = 6; # Число фаз
Pnom = 800e3;       # Номинальная мощность генератора, Вт (800 кВт)
wnom = 39.27;       # Номинальная угловая скорость вращения, рад/с
Mnom = Pnom/wnom;   # Номинальный момент, Н·м
Inom = 2040;        # Номинальный ток, А
Lsd = 0.0023;       # Индуктивность по продольной оси (d-axis), Гн
Lsq = 0.002;       # Индуктивность по поперечной оси (q-axis), Гн
Rs = 0.077;         # Активное сопротивление статора, Ом
zp = 8;             # Количество пар полюсов
Psi_pm = 2/3*(Mnom/(zp*Inom)); # Потокосцепление постоянных магнитов, Вб
J = 113;            # Момент инерции ротора, кг·м²

# Параметры преобразователя
Unom = 380;         # Номинальное выходное напряжение, В
Fpwm = 10000;       # Частота ШИМ (широтно-импульсной модул��ции), Гц
Tpwm = 1/Fpwm;      # Период ШИМ, с

# Примеры использования
T_Clarke, T_Inverse = clarke_transform_nphase(n)

Для конфигураций с числом фаз n = 3, 4, 5, 6, 8, 10, 12 тригонометрические коэффициенты в матрицах представимы через квадратные корни из рациональных чисел. Это обеспечивает максимальную точность вычислений — прямая и обратная трансформации являются точно взаимообратными с точностью до машинного эпсилон (около 10-16 для типа double). Для этих конфигураций долговременное моделирование абсолютно стабильно без накопления ошибок округления.

Блоки преобразований в модели

В модели Engee трансформации реализованы в виде отдельных функциональных блоков, использующих предвычисленные матрицы.

Блок прямого преобразования Кларка принимает вектор n фазных токов [i₁, i₂, ..., in]ᵀ и выполняет матричное умножение на предвычисленную матрицу Tαβ/n, получая на выходе два ортогональных компонента [iα, iβ]ᵀ. Это преобразование сокращает число переменных с n до двух, радикально упрощая дальнейшие вычисления.

Блок прямого преобразования Парка преобразует компоненты iα и iβ из неподвижной системы координат во вращающуюся систему d-q, используя текущий электрический угол поворота ротора θelec. Формулы преобразования включают умножение на матрицу вращения с элементами cos(θelec) и sin(θelec). В результате получаются компоненты id и iq, которые в установившемся режиме являются постоянными величинами.

Блоки регуляторов токов работают с постоянными величинами id и iq, что кардинально упрощает синтез управления. Два независимых ПИ-регулятора формируют напряжения ud и uq для отработки заданных токов id_ref и iq_ref.

Блок обратного преобразования Парка преобразует вычисленные напряжения ud и uq обратно в неподвижную систему координат, получая uα и uβ. Это зеркальное преобразование относительно прямого, использующее ту же матрицу вращения, но транспонированную.

Блок обратного преобразования Кларка преобразует два напряжения uα и uβ в n фазных напряжений [u₁, u₂, ..., un]ᵀ путём умножения на предвычисленную матрицу Tn/αβ. Эти напряжения подаются на ШИМ-модулятор для формирования управляющих импульсов инвертора.

Результаты моделирования в Engee

Модель была протестирована для шестифазной системы. Параметры машины соответствовали синхронному двигателю мощностью 800 кВт с номинальной угловой скоростью 39.27 рад/с (приблизительно 375 об/мин). Основные параметры представлены в коде скрипта инициализации.

Электрические параметры двигателя включают индуктивности по продольной и поперечной осям Ld = 0.0023 Гн и Lq = 0.002 Гн соответственно, активное сопротивление обмотки фазы Rs = 0.077 Ом, потокосцепление постоянных магнитов ψPM = 6.67 Вб. Номинальный фазный ток составляет 2040 А, что при шести фазах обеспечивает номинальную мощность 800 кВт. Момент инерции ротора J = 113 кг·м² отражает массивную конструкцию, характерную для мощных тихоходных машин.

Параметры силового преобразователя: номинальное выходное напряжение 380 В, частота ШИМ-модуляции 10 кГц. Высокая частота коммутации обеспечивает низкий уровень высших гармоник в токах и снижает акустический шум.

Сценарий испытаний включал несколько характерных режимов работы, позволяющих оценить динамические и установившиеся характеристики системы. Сначала выполнялся разгон от нуля до номинальной скорости 39.27 рад/с, затем система работала на номинальной скорости в режиме холостого хода без нагрузки на валу. После установления переходных процессов производился наброс номинальной нагрузки с моментом сопротивления 20.4 кН·м, соответствующим номинальной мощности. Далее система продолжала работать под нагрузкой до конца симуляции, что позволило оценить долговременную стабильность.

Результаты моделирования (рис. 7) убедительно подтвердили корректность реализации обобщённой трансформации Кларка для шестифазной конфигурации. Наблюдается быстрая отработка задания по скорости с плавным выходом на номинальный режим без перерегулирования. Переходный процесс при разгоне характеризуется монотонным нарастанием скорости, что свидетельствует о правильной настройке регулятора скорости и достаточной динамике внутренних контуров тока. При набросе нагрузки система демонстрирует хорошее качество регулирования — провал скорости минимален и составляет приблизительно 4 % от заданного значения, что не является критичным.

Рисунок 7 - Переходный процесс по скорости
Рисунок 7 - Переходный процесс по скорости

Фазные токи двигателя имеют симметричную синусоидальную форму со сдвигом 60 градусов (2π/6 радиан) между соседними фазами, как и должно быть для идеально симметричной шестифазной системы. Коэффициент гармонических искажений THD не превышает 3 процентов, что является отличным показателем и обусловлено высокой частотой ШИМ 10 кГц в сочетании с качественной фильтрацией индуктивностями обмоток. Форма токов близка к идеальной синусоиде, высшие гармоники подавлены до незначительного уровня. Это важно для минимизации дополнительных потерь в стали и меди, а также для снижения пульсаций момента.

Распределение мощности между фазами происходит правильно и равномерно. Все шесть фазных токов (рис. 8) имеют одинаковую амплитуду, что подтверждает отсутствие паразитной несимметрии, вносимой системой управления или численными ошибками. При номинальной мощности 800 кВт и номинальном токе 2040 А амплитуда фазного тока составляет около 1020 А (действующее значение 721 А на фазу). Для сравнения, эквивалентная трёхфазная машина той же мощности требовала бы фазных токов около 2050 А действующего значения — на 73% больше. Это наглядно демонстрирует преимущество многофазной конфигурации для снижения требований к силовым компонентам инвертора.

Рисунок 8 - Графики токов всех 6 фаз под нагрузкой
Рисунок 8 - Графики токов всех 6 фаз под нагрузкой
Рисунок 9 - Графики переходных процессов токов во всех 6-ти фазах
Рисунок 9 - Графики переходных процессов токов во всех 6-ти фазах

Система управления работает стабильно без колебаний и автоколебаний во всех режимах — как при холостом ходе, так и под номинальной нагрузкой. Графики компонент тока id и iq показывают, что компонента id поддерживается близкой к нулю (стратегия MTPA — максимальный момент на ампер), а компонента iq пропорциональна моменту нагрузки. Компоненты тока id и iq в установившемся режиме являются практически постоянными величинами с пульсациями менее 1 процента от среднего значения (рис.10, 11). Это ключевой показатель корректности всей цепочки преобразований: измерение шести фазных токов → прямое преобразование Кларка (6 → α,β) → прямое преобразование Парка (α,β → d,q) → регуляторы → обратное преобразование Парка (d,q → α,β) → обратное преобразование Кларка (α,β → 6) → формирование шести фазных напряжений. Любая ошибка в матрицах трансформации или в реализации преобразований проявилась бы как пульсации или дрейф компонент d-q.

Рисунок 10 - Переходные процессы тока id
Рисунок 10 - Переходные процессы тока id
Рисунок 11 - Переходные процессы тока iq
Рисунок 11 - Переходные процессы тока iq

Особенно важным результатом является долговременная стабильность моделирования. Симуляция на интервале 100 секунд модельного времени (что эквивалентно сотням тысяч периодов коммутации) не выявила накопления ошибок или дрейфа системы координат. Энергетический баланс сохраняется с высокой точностью: разность между электрической мощностью, потребляемой из сети постоянного тока инвертором, и механической мощностью на валу двигателя в точности равна потерям в активных сопротивлениях обмоток и потерям на коммутацию в инверторе с точностью лучше 0.01 процента. Это подтверждает, что предвычисленные матрицы трансформации Кларка для шестифазной системы математически корректны, и преобразования действительно являются взаимообратными с точностью до машинного эпсилон.

Таким образом, результаты моделирования в системе Engee полностью подтвердили работоспособность и корректность обобщённой трансформации Кларка для шестифазной конфигурации синхронного двигателя мегаваттного класса. Модель демонстрирует высокую точность, вычислительную эффективность и долговременную стабильность, что делает её пригодной как для научных исследований, так и для практического проектирования систем управления многофазными электрическими машинами.

Моделирование многофазных электрических машин — не тривиальное обобщение трёхфазного случая. Возможность построения эффективной модели связана с глубокими результатами классической математики — теоремой Гаусса-Венцеля о построении правильных многоугольников и теорией простых чисел Ферма.

Ключевые выводы работы

  1. Обобщённая трансформация Кларка позволяет преобразовать систему с n фазными переменными в эквивалентную двухкоординатную систему α-β. Это радикально упрощает моделирование, сводя систему из n дифференциальных уравнений с периодически меняющимися коэффициентами к системе из двух уравнений с постоянными коэффициентами.

  2. Математическая корректность и вычислительная стабильность трансформации зависят от возможности точного представления тригонометрических коэффициентов через квадратные корни. Эта возможность существует тогда и только тогда, когда число фаз представимо в виде n = 2k · p₁ · p₂ · ... · ps, где pi — различные простые числа Ферма (3, 5, 17, 257, 65537).

  3. Практически применимыми для моделирования и управления являются системы с n = 3, 4, 5, 6, 8, 10, 12 фазами. Эти конфигурации обеспечивают долговременную стабильность численного моделирования без специальных мер по компенсации накопления ошибок.

  4. Ядро математической модели — дифференциальные уравнения в системе координат d-q — не зависит от числа фаз. Переход от трёхфазной машины к многофазной требует замены только матриц трансформации на входе и выходе модели, внутренняя динамика описывается одинаково.

  5. Системы управления, разработанные для трёхфазных машин (векторное управление, прямое управление моментом и т.д.), без изменений переносятся на многофазные машины при использовании соответствующих трансформаций координат.

Направления дальнейших исследований

Описанный в статье подход открывает несколько направлений для дальнейшей работы.

Оптимизация для специфических применений. Различные конфигурации многофазных систем имеют свои преимущества. Четырёхфазные системы исключительно просты в реализации благодаря целочисленным коэффициентам. Пятифазные обеспечивают хороший компромисс между снижением фазных токов и сложностью. Шестифазные и дуальные трёхфазные позволяют использовать стандартные трёхфазные компоненты. Необходим детальный сравнительный анализ этих конфигураций для различных классов применений.

Учёт высших гармоник. В реальных машинах присутствуют высшие пространственные гармоники магнитного поля, которые в простейшей модели не учитываются. Многофазные системы позволяют использовать эти гармоники для создания дополнительного момента или, наоборот, компенсировать их влияние. Требуется обобщение трансформации Кларка для учёта нескольких пространственных гармоник.

Несимметричные режимы. Одно из ключевых преимуществ многофазных систем — способность работать при отказе одной или нескольких фаз. Моделирование таких несимметричных режимов требует модификации трансформаций координат. Это активная область исследований с важными практическими приложениями в авиации и других критичных системах.

Многомерные системы управления. Традиционное векторное управление использует только две степени свободы (оси d и q). В системах с числом фаз больше трёх появляю��ся дополнительные степени свободы, связанные с высшими пространственными гармониками. Разработка алгоритмов управления, использующих эти дополнительные каналы, может дать новые возможности для оптимизации характеристик привода.

Силовая электроника. Переход к многофазным системам требует развития специализированных топологий инверторов и методов модуляции. Классические методы пространственно-векторной ШИМ нуждаются в обобщении на произвольное число фаз. Это взаимосвязанная задача: выбор топологии инвертора влияет на возможности реализации алгоритмов управления.

Проблемные конфигурации. Системы с n = 7, 9, 11 фаз в настоящее время избегаются из-за проблем с точностью коэффициентов трансформации. Однако с точки зрения распределения мощности семифазная система выглядела бы весьма привлекательно. Возможно, существуют альтернативные подходы к моделированию и управлению, которые позволили бы эффективно работать с этими конфигурациями.

Заключительные замечания

Многофазные электрические машины переходят из категории исследовательской экзотики в категорию промышленных решений. Этот переход стал возможным благодаря развитию силовой электроники, микропроцессорной техники и, что не менее важно, математических методов моделирования и управления.

Обобщённая трансформация Кларка, основанная на классических результатах Гаусса и Венцеля, обеспечивает математический фундамент для этого перехода. Она позволяет использовать хорошо развитую теорию векторного управления трёхфазными машинами, распространяя её на системы с произвольным (в определённых рамках) числом фаз.

Для практикующих инженеров ключевое сообщение простое: при проектировании многофазных систем следует отдавать предпочтение конфигурациям с n = 3, 4, 5, 6, 8, 10, 12 фазами. Эти числа не случайны — они определяются фундаментальными математическими ограничениями, известными уже более двух веков. Игнорирование этих ограничений приведёт к усложнению разработки, проблемам с моделированием, необходимости применения специальных численных методов.

Для исследователей открывается обширное поле деятельности. Многофазные системы предоставляют новые степени свободы для оптимизации характеристик электропривода. Дополнительные фазы можно использовать для снижения пульсаций момента, компенсации высших гармоник, повышения отказоустойчивости, улучшения энергетических показателей. Полный потенциал этих возможностей ещё предстоит раскрыть.

В более широком контексте эта работа иллюстрирует непреходящую ценность фундаментальной математики. Результаты, полученные из чистого интеллектуального любопытства столетия назад, находят применение в современных технологиях. Это служит напоминанием о том, что инвестиции в фундаментальную науку — это инвестиции в будущие технологии, даже если связь между ними не всегда очевидна сразу.

Полезные ссылки

Подробное видео с объяснением: