Основы теории синхронных машин с несколькими обмотками на статоре Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

УДК 622

ОСНОВЫ ТЕОРИИ СИНХРОННЫХ МАШИН С НЕСКОЛЬКИМИ ОБМОТКАМИ НА СТАТОРЕ

Стрижков Игорь Г ригорьевич д.т.н., профессор

Кубанский государственный аграрный университет, Краснодар, Россия

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

Ключевые слова: ТЕОРИЯ СИНХРОННЫХ МАШИН, МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ, МЕТОД ПРЕОБРАЗОВАНИЯ КООРДИНАТ

Общие положения теории синхронных машин

К вопросам теоретического исследования и математического анализа синхронных машин (СМ) различных конструкций в последние десятилетия были привлечены значительные научные силы. Современная теория СМ располагает самыми передовыми достижениями математики, физики и электромеханики. Значительный вклад в теорию СМ внесли российские ученые. Еще в 30-40 годы прошлого века широко известные работы Горева А.А. [1], Постникова И.М., Костенко М.П. и др. создали отечественную школу электромеханики, сохраняющую высокий авторитет в мире и сегодня.

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

Современной электромеханике известен целый ряд синхронных машин, имеющих несколько трехфазных обмоток, действующих в едином

UDC 622

FUNDAMENTALS OF THE THEORY OF SYNCHRONOUS MACHINES WITH MULTIPLE WINDINGS OF THE STATOR

Strizhkov Igor Grigorievich Dr.Sci.Tech., professor

Kuban State Agrarian University, Krasnodar, Russia

The article is an attempt to generalize the theory of the mathematical description of synchronous machines with multiple three-phase windings on the stator and arbitrary location in space of the machine. The features of the transformations of coordinates in the mathematical description of the machine are shown

Keywords: THEORY OF SYNCHRONOUS MACHINES, MATEMATICAL MODELING, TWO REACTION THEORY

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

В круг вопросов, отнесенных к общим вопроса построения математического описания синхронной машины с несколькими электрически и магнитно связанными обмотками, входят описание потокосцеплений машины, уравнения равновесия напряжений в электрических цепях статорных обмоток, уравнения электромагнитной мощности и момента, активной и реактивной мощности, уравнение движения машины. Вопрос о целесообразности обобщения теории СМ с несколькими обмотками на статоре поставлен в научной литературе [2] и определенная работа проделана и по реализации этой проблемы [2, 3, 4]. Однако завершенной формы теория СМ с несколькими обмотками на статоре до настоящего времени не получила. В настоящей работе сделана попытка такого обобщения.

Объект математического описания и принятые допущения

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

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

Ротор вращается в общем случае с произвольной угловой частотой О; положение его относительно статора характеризуется углом у между магнитной осью фазы а1 и продольной осью ротора ё.

Рисунок 1 - Расположение осей обмоток многообмоточной синхронной машины

Как и классической синхронной машине, рассматриваемой машине присущи идеализирующие допущения: отсутствие насыщения магнитной цепи; синусоидальное распределение магнитного поля и МДС вдоль воздушного зазора; симметрия трехфазных обмоток статора; симметрия магнитопровода относительно осей ё и q; отсутствие гистерезиса, потерь в стали и вытеснения тока в проводниках; независимость сопротивлений рассеяния от пространственного положения ротора; постоянство активных сопротивлений обмоток; токи во всех обмотках переменного тока (якорных и демпферных) принимаются синусоидальными.

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

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

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

Потокосцепления обмоток в фазной системе координат

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

В теоретической электромеханике, как правило, применяют симметричные многофазные координатные системы, когда фазные координаты представляют собой симметричную многолучевую звезду. Для трехфазных машин фазные координатные оси обозначают буквами а, Ь, с. Подробно свойства изображающих векторов анализируются в [1, 2, 5]. Здесь рассматриваются свойства этих векторов применительно к случаю расположения на статоре нескольких трехфазных обмоток с несовпадающими магнитными осями.

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

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

Іa + 4 + 4 = 0 (1)

Строго говоря, условие (1) диктуется только необходимостью построения симметричных координатных осей. При токах ia, гь, 4, имеющих равную амплитуду, но временной сдвиг, отличающийся от 2п/3, их также можно представить в виде проекций на временные оси изображающего вектора, имеющего постоянную длину и вращающегося с постоянной угловой скоростью. При этом угол поворота осей должен быть равен сдвигу начальных фаз токов. Однако в рассматриваемой задаче принято условие симметрии токов и рассматривается случай симметричного расположения фазных осей.

Ряд авторов для придания наглядности процедуре представления скалярных величин в виде изображающего вектора [1, 5] совмещают фазные координатные оси с магнитными осями фазных обмоток. Такое совмещение не имеет физического смысла, поскольку для скалярных параметров машины как ток, напряжение, потокосцепление и др., оси координат имеют характер временных и к конструкции машины никакого отношения не имеют. При расположении на статоре нескольких трехфазных обмоток, их параметры могут определяться как проекции

соответствующих изображающих векторов на общие оси а, Ь, с независимо от углового смещения обмоток в пространстве машины.

.................а

Рисунок 2 - Изображающий вектор тока трехфазной обмотки

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

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

трехфазных систем координат для нескольких трехфазных обмоток.

Рассмотрим подобный случай на конкретном примере. Пусть на статоре СМ расположены две трехфазные обмотки, смещенные относительно друг друга на центральный угол в (рис. 3). Обозначим их условно номерами 1 и 2. Начальные фазы токов в фазных обмотках а1 и а2 соответственно ф1 и ф2. Токи каждой из обмоток симметричны.

Расположим оси а1, Ь1, с1 и а2, Ь2, с2 совпадающими с осями соответствующих фазных обмоток. Векторы ^ и ^2 результирующих МДС соответствующих обмоток займут положение под углом у1 и у2 осям а1 и а2. Причем

^ = 1, 2) = ф, - п/2 .

Рисунок 3 - Построение изображающих векторов МДС двух трехфазных обмоток с несовпадающими магнитными осями

Пространственное смещение изображающих векторов

а = ух - (у2 - в) = Ш + ф1 - Ш - ф2 + в = Ф1 - Ф2 + в.

(2)

Уравнение (2) устанавливает связь между угловым смещением изображающих векторов МДС обмоток, смещением начальных фаз токов и угловым смещением обмоток. В рассматриваемом случае результирующая МДС обеих обмоток определяется геометрическим суммированием изображающих векторов ¥1 и ¥2.

Из уравнения (2) следует, что максимальная результирующая МДС двух обмоток при фиксированных амплитудах токов имеет место при совпадении векторов ¥1 и ¥2 по направлению, что имеет место при смещении начальных фаз токов на угол -в. Если эти обмотки включены параллельно, то общий или суммарный ток будет меньше алгебраической суммы токов. И наоборот, совпадающие по начальной фазе токи обмоток 1

и 2 создадут результирующую МДС меньшую, чем алгебраическая сумма МДС обмоток.

Однако, при углах а не превосходящих +15 град. эл., расхождение между алгебраической и геометрической суммами МДС невелико. Так, при равенстве модулей Е1 и Е2 и а = 15 град. эл., алгебраическая сумма превосходит геометрическую на 1,7 %.

Таким образом, изображающие векторы скалярных параметров следует представлять, когда это возможно, в единой трехфазной системе координат, в то время как векторы МДС многофазных обмоток

- в собственных координатах, оси которых совмещены с осями соответствующих фазных обмоток.

Чтобы выявить аналитическое выражение для суммарной МДС нескольких многофазных обмоток с учетом неодинаковости начальных фаз токов и смещения обмоток в пространстве машины, воспользуемся формой представления МДС фазной обмотки при произвольных начальных фазах токов и угловых положениях обмоток, предложенных [5]

Е(в, 0 = ^'ф1т собэ(ю? - фп - в + вп) + F"фlm соб(- Ш + фп - в + вп),

где Еф1т, Е"ф1т - амплитуда соответственно прямой и обратной волны 1-й

гармоники фазной МДС;

фп - начальная фаза тока обмотки;

вп - угол между действительной осью координат комплексной плоскости и осью рассматриваемой обмотки (рис. 4.); в - угловая координата рассматриваемой точки.

В комплексной форме функция МДС обмотки имеет вид:

Е(в, 0 = £'ф1т + Е 'ф1т, где Еф1т = Еф1т ехр[/(ю? - фп)] ехр(/вп) = Е'ф1т ехр(/Ш) ехр[/(вп - фп)]

изображает прямую вращающуюся волну МДС фазы, а комплексная функция

Е 'ф1т= Е'' ф1т exp(-/Ш) exp[/'(вn + фп)]

- обратную волну.

а

Рисунок 4 - Изображающий вектор тока и его проекции на оси

Для представления результирующей МДС трехфазной обмотки, просуммируем МДС всех фазных обмоток

£1т £ 1т £ а1т + £ Ь1т + £ с1т £ ф1т ехр[/ фп)] ехр(/ Рп) +

+ £ 'ф1т ЄХр[/ (Ш - фп - 2п/3)] ехр[/' (Рп + 2п/3)] + £ 'ф1т ехр[/ (Ш - фп - 2п/3)] X Хехр[(Рп - 2п/3)] = 3 £ 'ф1т = £ 'ф1т ЄХр(/' Ш) ехр[(Рп - фп)]. (3)

Уравнение (3) является аналитическим выражением МДС одной трехфазной обмотки. Для определения суммарной МДС N статорных обмоток следует суммировать МДС этих обмоток в комплексной форме.

N N

Ерез = Е Ет = ^/Ю *)Е Е (Дп - Фп )]} . (4)

1 1

Уравнение (4) показывает, что при изображении векторов МДС на векторной диаграмме ЭМДС с общей системой координат, вектор МДС обмотки должен опережать вектор тока этой же обмотки на угол вп. В этом случае угловое смещение обмоток проявляет себя как эквивалентное смещение начальной фазы тока.

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

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

Преобразование координат может быть реализовано различным образом. Для теории СМ фундаментальное значение имеет преобразование переменных, эквивалентное замене координат осей а, Ь, с на

ортогональные, жестко связанные с ротором, т.е. вращающиеся с угловой частотой ротора. При этом одна из осей принимается совпадающей с продольной осью ротора (ё), другая ^) располагается перпендикулярно, т.е. по поперечной оси ротора. В отечественной литературе ось q чаще принимают опережающей ось ё, однако это требование не является обязательным.

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

Токи фаз обмотки Жп представим как проекции изображающего вектора Іп на оси а,Ь,с, предполагая, что условие (1) выполняется:

/сп = 1п соб(5п + р). 0 Преобразование координат осуществляется использованием следующих уравнений:

где [4^о]п coloп(idn, iqп, /оп);

[/а,Ь,с]п colon(ian, ІЬп,> /сп);

[^п] - матрица преобразования координат.

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

Іап Іп с°з5п,

ІЬп !п с°8(§п р),

(5)

(6)

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

Единственным обязательным требованием ко всякой матрице [Ап] является неравенство ее нулю в любой момент времени *. Таким образом,

Наибольшее употребление в отечественной литературе по теории СМ [31 ] нашли преобразования, выполняемые с использованием матрицы

[Ап] Ф 0.

(7)

Рисунок 5 - Преобразование тока обмотки к осям по матрице

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

І = А[ІТ+ІЇ •

Как уже отмечалось, аналогичные выражения можно записать для напряжений, ЭДС и потокосцеплений.

Подставив выражение (5) в уравнение (6) и положив у = О^ - у0, получим

їа = I соб(7-£);

Ід =-1 эшО-^;

І0 = 0.

Эти же соотношения получим из векторной диаграммы по рис. 4. Это означает, что амплитуда тока в относительных единицах при переходе от а, Ь, с координат к координатам ё, q, ° не изменяется.

Обратные преобразования осуществляются по формуле

где

соб у - Бт у 1 [Ак ]= с°8(у-р) - 81п(у-г) 1 соб(у+г) - вт(у+р) 1

(8)

Установленная аналогия аналитического преобразования переменных величин их геометрической интерпретации позволяет заключить, что

где /dn , /дп - векторные величины, причем |^п| = 4п; |/чп| = /чп.

В классической теории СМ широкое распространение получила идея представления переменных параметров в координатах d, q, ° в комплексной форме [2]. Для этого ось d совмещают с вещественной (+1), и ось q с мнимой (+/) осями. При этом вместо вещественных параметров с индексами d и q используются комплексные h(h = /, и, е, у) = hd + ^, являющиеся изображающим вектором переменной величины h.

Преобразование напряжений по аналогии с токами осуществляется по уравнениям:

[ud,q,o]n [А1] [ua,b,c]n,

где [ud,q,o]n colon(udn, иоп-> и°п);

[Ua,b,c]n = C0l0n(u

an, Ubn, исп).

Обратные преобразования

[ua,b,c]n [АК] [^^°]п;

(9)

и = 1и2 + и2 •

где и н \ ёп цп

и и.Ап + идп5 иёп |иД ицп I идп|-

Применение преобразований с матрицей (7) вносит еще одну существенную особенность в уравнения математического описания СМ по типу уравнений Парка-Горева. В уравнениях потокосцеплений обнаруживается относительно малое заполнение матрицы индуктивных сопротивлений, что приводит к облегчению решения системы уравнений Парка-Горева аналитическими методами. Для СМ с одной трехфазной обмоткой уравнения потокосцепления имеют вид [2]

У а " ха 0 Хауё 0 Ха/ "

У, 0 Хч 0 X ауц 0

= Хауё 0 X ауц 0 Х/уа X ^ у а

0 X ауц 0 X уч 0 1 уч

У/ _ к 1 0 х/уа 0 х/ _ _ 1Г _

Как будет показано ниже, применение матриц преобразования переменных параметров (7) и (8) в случае наличия на статоре нескольких обмоток, сдвинутых в пространстве на некоторый угол, отличный от нуля и не кратный п, приводит к появлению взаимной индукции разных обмоток не только в одноименных осях ортогональной системы координат, как в случае соосных обмоток, но и в разноименных осях, что ведет к более полному заполнению матрицы индуктивных сопротивлений. Это привносит определенное усложнение и новизну в теорию преобразования координат по типу Парка-Горева. Рассмотрим этот вопрос более детально.

Желание максимальным образом уменьшить заполнение матрицы индуктивных сопротивлений (индуктивностей), привело ряд авторов к выводу о целесообразности использования разных матриц преобразования

переменных параметров для разных статорных обмоток. Так, в [3] предлагается использовать матрицу вида

С08(/-Ь) СО$(у-0-р) со$,(у-0 + р)

- ^т(у-р) - $\п(у-р-р) - $лп(у-0 +р)

7^5 705 л/05

(10)

где различие в коэффициентах 2/3 и (2/3)12; 1/2 и 2-12 объясняется различным подходом к формированию базисных величин при переходе к относительным единицам.

Аналогичные преобразования реализованы в [4] применительно к СМ с обмотками, смещенными на фиксированный угол п/6. Матрица (10) учитывает пространственный угол смещения обмоток и является матрицей, максимально разрешающей матрицу индуктивных

сопротивлений машины в координатах ё, q. Кроме того, она не равна нулю, вследствие чего математически корректна.

С точки зрения геометрической интерпретации, преобразования по матрице (10) предусматривают поворот вектора тока обмотки Ж в сторону отставания на угол в с последующим проецированием его на оси ё и q. Это можно трактовать как перемещение обмотки Ж до совпадения магнитных осей обеих статорных обмоток (частный случай двойной статорной обмотки впервые рассмотрен автором в статье [6]).

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

поворотом каждого вектора на угол ви. В случае последовательного включения двух или нескольких обмоток статора, расчетный ток в обмотках будет иметь разные значения составляющих по осям d и q, т.е. в последовательных участках цепи будут протекать различные расчетные токи (естественно, что такое несоответствие имеет место только в промежуточных вычислениях в координатах обобщенного вектора; при переходе к естественным координатам а, Ь, с токи примут естественный характер распределения). Вследствие наличия такого недостатка, на наш взгляд весьма существенного, этот способ преобразования координат при решении ряда задач, как, например, задачи анализа статических режимов СМ, не может быть рекомендован как оптимальный. Однако, в задачах сопряженных с использованием численных методов, для которых характерно пренебрежение наглядностью промежуточных результатов вычисления, подобный недостаток вполне допустим, и метод может быть рекомендован для применения.

Ниже представлен анализ преобразований, выполненных как по единой для всех статорных обмоток матрице преобразований, так и с использованием индивидуальных матриц преобразования для каждой обмотки, имеющих вид

Уравнения равновесия напряжения статорных обмоток в фазных координатах можно представить в следующем виде:

Цп = гпЦп + №^/Ж.

(] = а, Ь, с; п = 1, 2, ..., Ы)

С°8(у - Рп) С0»(У - Рп - Р) С0»(У - Рп + Р)

-э!п(у рп) -^п(У-Рп-Р) -^п(У-рп +Р) . (11)

0,5

0,5

0,5

Уравнения равновесия напряжения в осях d, q, о после реализации преобразования по уравнению (9) с матрицей преобразования (7) или (11) (любой из этих матриц), дает уравнения:

^п rnidn + dn/dt ^qndY/dt; )

^п = Гп^п + d¥qn/dt + ^п^уМ .

Потокосцепления в машине с несколькими обмотками на статоре

Анализ индуктивных связей обмоток базируется на известном методе представления постоянных и переменных параметров машины в естественных а, Ь, с и ортогональных вращающихся d, q координатных системах.

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

а) Потокосцепления обмоток в а, Ь, с координатах

В наиболее общем виде потокосцепления как функции токов в обмотках в естественных координатах можно представить в матричной форме следующим уравнением

т = №м,

где А(А=У, 0 = С°1оп[(Йа1, АЫ, Ас1), (^2, Аь2, М, .^(Ак, АьК, Ак), (^, hyq, А)]; [Ь] - по матрице (12).

[ь] =1

2

Ьа1 Ьа1Ъ1 Ьа1с1

ЬЪ1а1 1 Ь 1 с1

Ьс1а1 Ьс1 1 Ьс1

ЬаИа1 ЬМ 1 ЬаИс1

ЬЪЫа1 ЬЪМ1 ЬЪЫс1

ЬсЫа1 ЬсМ 1 Ьс^с1

Ьапа1 ЬапЪ 1 Ьапс1

~о Ь п 1 Ь пс1

Ьспа1 Ьсп 1 Ьспс1

■К 1 ■К с1 ■К

Ьуда1 1 ■К Ьудс1

ы 1

Ьа1аМ Ьа1ЪN Ьа^ . .. Ьа1ап Ьа1 п Ьа1 сп •К а1 •К а1 Ьа^

ЬЪ 1аМ ЬЪ1Ъ N ЬЪ 1cN . .. ЬЪ 1ап Ь 1 п Ь 1сп ЬЪ ^ ЬЪ1yq 1

Ьс1аМ Ьс1ЪN Ьc1cN . .. Ьс1ап Ьс1 п Ьс1 сп Ьc1yd Ьc1yq Ьс^

ЬаЫ ЬaN Ъ N ЬaNcN . .. Ь<$^ап Ь^Ъп ЬаЫсп ЬаN

ЬЪЫ^ % ы ЬЪNcN . .. ЬЪNan Ьът п ЬЪШп ЬЪNyd ЬЪ^ 1

ЬcNaN ЬcNЪN Ь N .. Ь^^ап ЬсШп ^с№сп ЬcNyd ЬcNyq ЬсN

ЬапаЫ ЬапЪN .. Ьап Ьап п Ьапсп Ьanyd Ьanyq Ьап

ЬЪп^ ЬпЪ1Я ЬЪпс^Ы . .. ЬЪ пап п ~о Ь псп ЬЪnyd ~о > ~о

ЬcnaN Ьсп ЪN .. Ьспап Ьсп п п с с Ьcnyq Ь сп

ЬydaN ЬydЪN Ьу^с!^ . п а ■К п ~о ■К п ■К ■к 0 Ь у¥

ЬyqaN Ьyq ЪN ЬyqcN . .. Ьyqan п ~о ■К Ьyqcn 0 ■К 0

X ЬМ . п Ьрп п Ьfyd 0

(12)

Согласно принципу взаимности, взаимные индуктивности с взаи-мообратными индексами равны между собой, т.е. Ьа1Ь1 = ЬЬ1а1; ЬЬ1с2 = Ьс2Ь1;

••• -^опуё Ьydcn, и т.д., но ЬЬ1с2 ф Ьс1Ь2 и т.д.

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

Пространственное положение ротора определяется углом у между продольной осью d и осью фазной обмотки а1 (уа1 = у).

б) Собственные индуктивности фазных обмоток статора

Известно [5], что эти индуктивности являются периодическими функциями угла ух с периодом п. Причем эта функция имеет вид:

Ьх = 4 + ^СОБУх. (13)

Постоянная составляющая 1О представляет собой сумму постоянной составляющей индуктивности 4, соответствующей потоку в воздушном

зазоре, и рассеяния 4, практически не зависящей от положения ротора. Относительно 1О можно сказать, что она определяется параметрами самой фазной обмотки, а также магнитной проводимостью воздушного зазора.

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

7а1 = 7; 7Ь1 = 7 + 2р = 7 - р; 7е1 = 7 + р ;

7ап = 7 - вп; 7Ьп = 7 - вп - р; 7сп = 7 - вп + р. (14)

Поскольку рассматриваемые обмотки имеют одинаковое число витков, уравнения (13) и (14) можно записать в ином виде:

Ьап 1оп + 12п соб2(7 вп);

ЬЬп 1оп + 12п соб2(7 вп р);

Ьсп 1оп + 12п соб2(7 вп+ р);

п = 1, 2, ..., N в1 = 0.

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

Ьап 1оеп + 4п + 12псоб2(7 вп);

ЬЬп 1оеп + 4п + 12псоб2(7 в р); (15)

Ьсп 1оеп + 4п + 12псоб2(7 вп + р).

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

1оп = 1о1 (№п ^обп)2/(^1 ^об1)2;

12п = 121 (№п ^обп)2/(^1 ^об1)2;

4, = /,1 Жп2/^12.

Для обмоток с равными числами витков и обмоточными коэффициентами составляющие коэффициента самоиндукции разных обмоток отличаются только аргументом косинуса.

Уравнения (15) удобно представить в общем виде:

и(к = ап, Ъ п, Сп) = [/о1 + /21соб2(уй- вп)] (Жп./Жъ)2,

где №п_ - эффективное число витков обмотки п (№п- = №п £обп).

В более краткой форме

/п = /1 (Жп-/^1-)2.

В теории СМ принято выражать постоянную и переменную составляющие индуктивности через ее экстремальные значения или, иначе говоря, через индуктивности по продольной оси /,^ и поперечной /д. Поскольку среднее значение синусоидальной величины равно среднему арифметическому его максимального и минимального значений

/оп (4п + /дп)/2.

Амплитуда выражается другой известной формулой

/2п (^п /дп)/2.

Понятно, что для СМ неявнополюсной конструкции = /чк и переменная составляющая взаимной индуктивности обращается в нуль.

в) Взаимные индуктивности фазных обмоток

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

виде составляющих по продольной и поперечной осям. Для этого МДС обмотки, создающую синусоидальную пространственную волну, амплитуда которой совпадает с осью обмотки а1, разложим на две волны, действующие по осям d и д. Амплитуды этих составляющих равны:

Fad = Е собу; Е Бту.

Аналогично могут быть разложены ток и потокосцепление обмотки а1, т.е.

^1 = /а1 ^1 Сову; ^Ч1 = - /а1 /д1 вШу. (16)

Потокосцепление этих волн с любой обмоткой статора зависит от угла между осью этой обмотки и осью d. Например, для обмотки Ъ1, составляющей с осью d угол (у - р), потокосцепление ¥а1Ь1 запишется в виде

^а1Ь1 = ¥Ыа1 = ^1 соб(у - 2п/3) - ^ бш(у - 2п/3),

что после известных преобразований, подробно представленных в [1], дает уравнение взаимной индуктивности

-^а1Ь1 -^Ь1а1 (/d1 + /д1)/4 (/d1 /д1)/2 соБ(2У + п/3);

-^Ь1о1 -^о1Ь1 (/d1 + /д1)/4 (/d1 /д1)/2 соБ(2У п);

^а1с1 с1а1 (/d1 + /д1)/4 (/d1 /д1)/2 соБ(2У + 5п /3).

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

Взаимные индуктивности фазных обмоток № имеет те же особенности, что и обмотки Ж1, но, вследствие смещения этих обмоток

относительно оси начала отсчета угла у, изменяется аргумент косинуса. При этом имеем:

^апЬп -^Ьпап 1оп/2

+ 12п соб2(у - вп + р);

-^Ьпсп L'спЬп 1оп/2

+ 12п СОБ2(у — вп);

^апсп ^спап 1оп/2

+ 12п ООБ2(у — вп — р).

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

Для взаимных индуктивностей фаз обмотки имеем:

ЬапЬп = Ььпап = [- /01/2 — /21 сов2(у — вп + р)] ^/Ж^2;

^Ьпсп = LcnЬn = [— 1о1/2 — 121 сОБ2(У — вп)] (Жп/Ж1)2;

Ьапсп = Ьспап = [— /о1/2 — /21 СОв2(у — вп + р)] ( Ж/Ж^)2.

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

¥а1ап = ^аСОБ(у — вп) — ^т(у — вп),

что с учетом уравнения (16) дает

^а1ап = ^[/«пСОву СОв(у — вп) + /ц^Шу БШ(у — вп)].

После тригонометрических преобразований получим:

La1an 0,5[(/ё1 + /д1)соввп + (/ё1 ^1)сов(2У вп)].

Индуктивная связь обмоток изменяется с изменением угла вы; в то же время для вы = 0 имеем Ьа1аЫ = Ьа1.

Вследствие симметричного расположения пар фазных обмоток a1-aN, b1-bN, c1-cN, взаимные индуктивности последних двух пар обмоток имеют вид:

£b1bn = 0,5[(ld1 + lq1) cosPn + (Zd1 - V) cos(2y - Pn - p)]; )

Lc1cn 0,5[(/d1 + lq1) C0SPn + (/d1 lq1) cos(2Y Pn + p)],

или в более краткой форме записи:

La1an Lana1 lo1 cosPn + l21 C°s(2Y Pn);

Lb1bn Lbnb1 lo1 cosPn + l21 cos(2Y Pn p); (17)

Lc1cn Lcnc1 lo1 cosPn + l21 cos(2y Pn + p).

При необходимости учета взаимной индуктивности по потоку рассеяния соответствующую величину xs1n следует суммировать с членами правой части уравнений (17).

Имея ввиду, что коэффициенты взаимной индуктивности по полю в воздушном зазоре пропорциональны произведению чисел витков обмоток и распространяя результаты по уравнениям (17) на индуктивности произвольных обмоток n и т, получим:

Lanam = [/o1 cos(Pn - Pm) + ll1 cos(2Y - Pn + Pm)] Wn WJWx\

Lbnbm = [lo1 cos(Pn - Pm) + /21 cos(2Y - Pn + Pm - p)] Wn WJWx\

Lcncm = [/o1 cos(bn - Pm) + ll1 cos(2Y - Pn + Pm + p)] Wn WJW\ .

г) Взаимная индуктивность разноименных фазных обмоток

Взаимная индуктивность фаз а1 и bn определяется их взаимным смещением на угол (Pn + p). Для обмоток с одинаковым числом витков можно записать

¥a1bn = Z'a1[/d1 cosy cos(y - Pn - p) + /q1 sinY sin(Y - Pn - p)], что после преобразований дает:

La1bn = 0,25(/d1 + /q1)(V3 sinPn - cosPn) + 0,5(/d1 - /q1) cos(2Y - Pn - p) (18)

или La1bN = 0,5/o1 cos(Pn + p) + /21 cos(2y - Pn - p).

Поскольку пары фазных обмоток b1-cN и c1-aN расположены симметрично относительно пары a1-bn, их можно записать в виде:

Lb1cn Lcnb1 /o1n +

0,5(/d1 - /q1)cos (2y - Pn) ;

Lc1an Lanc1 /o1n + 0,5(/d1 /q1)cos(2Y Pn + p), (19)

где /o1n = /o1 cos(Pn + p) = - 0,25(/d1 + /q1)^V3 sinPn + cosPn).

Уравнения (19) показывают, в частности, что при bn = 0 сопротивления взаимной индукции совпадают с соответствующими значениями сопротивлений обмотки W1.

Взаимная индуктивность обмоток а1 и сп определяется их взаимным расположением под углом 2p + pn, следовательно для потокосцепления ^a1cn можно записать

¥a1cn= ia1 [/d1 cosy cos(y - Pn + p) + /q1 sinY sin(Y - Pn + p)], что после преобразований дает:

La1cn = 0,25(/d1 + /q1)^V3sinPn - cosPn) + 0,5(/d1 + /q 1) cos(2y - Pn + p) =

= /o1 cos(Pn - p) + /21 cos(2Y - Pn + p).

Сопротивления взаимной индукции фазных обмоток b1 и an наиболее просто может быть определено на основании следующих рассуждений. Из рис. 1 видно, что результирующие оси пар обмоток a1-bn и b1-an совпадают, следовательно, совпадают и аргументы косинуса вторых слагаемых по уравнениям (18). Аналогичным образом они совпадают для пар b1-cn и c1-bn, a1-cn и c1-an. Однако углы смещения катушек в каждой паре отличаются. Если в рассмотренном выше случае угол между фазами a1 и bn составляет (Pn + p) и совпадает с аргументом косинуса в последнем

из уравнений (19), то угол между осями фазных обмоток ап и Ь1 составляет (р - вп), вследствие чего постоянная составляющая /оп1 принимает вид:

/о1п = /о1 СОБ(Рп - р) = 0,25(41 + /ц1)( 73 Б1пРп - СОБРп)

Таким образом, можно записать

-^апЬ1 -^Ь1ап 1оп1

+ /21 ООБ(2у - вп - р);

-^Ьпо1 -^о1Ьп /оп1 + /21 СОБ(2у Рп);

-^опЬ1 -^Ь1оп /оп1

+ /21 еоБ(2у - вп + р).

Сопротивления взаимоиндукции между двумя произвольными обмотками Жп и Жт определим с использованием их аналогии с обмотками Ж1 и Жп.

Если результирующая ось обмоток а1 и ап смещена на угол вп/2 по отношению к оси а1, то ось пары ап и ап смещена на угол (вп - вт)/2, а угол смещения самих обмоток равен (вп - вт). В соответствии с этим

-^апат /о1 СоБ(вп - вт) + /21 ОоБ(2у - вп + вт);

-^ЬпЬт /о1 СоБ(вп - вт) + /21 СоБ(2у - вп + вт - р);

^спст /о1 соБ(вп - вт) + /21 соб(2у - вп + вт + г); (20)

-^Ьпат /о1 СоБ(вп - вт + р) + /21 СоБ(2у - вп + вт - р);

-^спЬт 1о1 СоБ(вп - вт + р) + /21 Со8(2у - вп + вт);

-^апСт /о1 СоБ(вп - вт + р) + /21 Со8(2у - вп + вт + р).

Распространяя полученные результаты на случай обмоток с неодинаковым числом витков, получим

^Ьпат = [/о1 СоБ(вп - вт + р) + /21 СоБ(2у - вп + вт - р)] ЖЖ/ЖД

^пЬт = [/о1 Со8(вп - вт + р) + /ц Сов(2У - вп + вт)] ЖШт/Ж?\ (21)

ДтСт = [/о1 СоБ(вп - вт + р) + /21 СоБ(2у - вп + вт + р)] ЖпЖт/ЖД

В наиболее общем виде взаимные индуктивности двух любых фазных обмоток можно записать

2

Ьт]т0' ^ Ь, С; ] ^ Ь, с) [/о1соБвщт + /21 соб(2У + в^т)] ЖпЖт/Ж1 . (22)

При необходимости учета взаимной индуктивности по потокам рассеяния, она должна быть добавлена в правую часть уравнений (20), (21) и (22) как величина /зпт, не зависящая от положения ротора.

д) Взаимные индуктивности обмоток статора и ротора Известно [5], что частота изменения переменной составляющей взаимной индуктивности неподвижных обмоток, расположенных на статоре, и вращающихся вместе с ротором обмоток возбуждения и демпферной в два раза меньше частоты переменной составляющей взаимных индуктивностей статорных обмоток. Это объясняется тем, что для взаимных индуктивностей статорных обмоток не имеет значения полярность магнитного полюса ротора, в то время как для взаимной индуктивности обмоток, расположенных соответственно на статоре и роторе, максимум наступает при совпадении магнитных осей обмоток. Кроме того, взаимным индуктивностям обмоток статора и ротора свойственно отсутствие постоянной составляющей взаимной индуктивности.

Предполагая, как и прежде, что учитываются лишь основные гармонические пространственного распределения магнитного потока, можно представить взаимные индуктивности в общем виде уравнением

Ь = Ь1 собо,

где Ь1 - взаимная индуктивность обмоток при совпадении их магнитных осей;

о - угол, линейно связанный с углом у, при котором совпадают оси обмоток.

Для взаимной индуктивности фазных обмоток обмотки Wn с обмоткой возбуждения уравнение принимает вид:

Lanf la1f Cos(Y — Pn);

Lbnf = lalf cos(Y - Pn - p); (23)

Lcnf = lalf COS(Y - Pn + p).

Взаимоиндуктивная связь статорных обмоток с продольной составляющей демпферной обмотки (yd) по форме аналогична уравнениям (23), поскольку эта демпферная обмотка расположена на роторе аналогично обмотке возбуждения, следовательно

Lanyd lyd Cos(Y — Pn);

Lbnyd = lyd

cos(Y - Pn - p); (24)

Lcnyd lyd Cos(Y - Pn + p).

Поскольку ось поперечной демпферной обмотки (yq) опережает продольную ось ротора d на угол п/2, взаимные индуктивности обмоток статора с обмоткой yq имеет вид:

Lanyq lyq Cos(Y Pn + п/2) lyq SlH(Y Pn);

Lbnyq = lyq cos(Y - Pn - p + п/2) = -lyq sln(Y - Pn - p); (25)

Lcnyq = lyq COS(Y - Pn + p + п2) = -lyq Sln(Y - Pn + p).

Взаимные индуктивности между обмоткой возбуждения и демпферной yd не зависит от пространственного положения ротора и потому является величиной постоянной: lfyd = ConSt. В то же время взаимная индуктивность обеих обмоток с поперечной демпферной обмоткой yq отсутствует, поскольку оси этих обмоток перпендикулярны: lfyq = lyqf = 0.

Для обмоток с различным числом витков уравнения по типу (23,...,25) принимают вид:

Lanf = lalf COS(Y - Pn) WJW^

Ьь^ = /аИ Соб(у - вп - р) Жп/Жь

Ь^ = /аИ Соб(у - вп + р) Жп/Жь

Ьапуё /аИ соб(У - вп) Жп/Ж1;

ЬЬпуё /а1£ соб(у - вп - р) Жп/Жь

Ьспуа = /аИ соб(у - вп + р) Жп/Жь

Ьапуя /а1уя яп(у - вп) Жп/Ж1;

ЬЬпуя /а1уя

бш(у - вп - р) Жп/Жь

Ьспуя /а1уя бш(у - вп + р) Жп/Жь

Как видим, взаимные индуктивности, определяющие потокосцепление обмоток статора с обмотками ротора, являются периодическими функциями углового положения ротора. Для вращающейся машины угол у в функции времени определяется выражением [2]

г

у = | № + ro,

0

где у0 - значение у при ? = 0.

Для равномерно вращающегося ротора

у = Ы + У0 .

Потокосцепления обмоток в ф д координатах

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

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

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

а) Преобразование координат при равенстве чисел витков статорных обмоток по условию минимального заполнения матрицы индуктивностей в координатах ф д

Потокосцепление обмотки а1:

^а1 ^‘а1Іа1 + -^а1Ь1^Ь1 + -^а1с^'с1 +...+ -^а1аК^аК + -^а1ЬК^ЬК + La1ydiyd +

(26)

Уравнение обратных преобразований токов

[Іа,Ь,о]п [-^] [id,q,o]nN,

(27)

где [Іа,Ь,с]п colon(ian, iЬn, Ісп);

[і^д,о]п colon(idn, iqn, Іоп);

Г с°8(>-Д,) 8Іп(г-Д,) 1

[£] = с^(^- Д -Р) «п(Г-Д -Р) 1 ^^Г-Д +р) ^п(Г- Д + Р) 1

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

^п ^2 (29)

где /б2 - взаимная индуктивность по потоку рассеяния фазных обмоток I и 7 при совпадающих магнитных осях; ву - угол смещения обмоток * и7.

Совместное решение уравнений (26,...,29) и уравнений индуктивности (12) дает:

¥а1 = (/о1 + /в1 + /21 СоБ2у)(/о1 + *11 СоБу + ;Ч1 БШу) +

+ [- /о1/2 + /б1 соБр + /21 соб(2у - р)][/о1 + *'«11 соб(у - р) + ^1 бш(у - р)] +

+ [- /о1/2 + /б2 собр + /21 соб(2у + р)][*о1 + *11 соб(у + р) + /Ч1Б1п(у + р)] + ...

...+ [/о1 соБвк + /б2 соБвк + /21 соб(2у - впк][г'ок + соб(у - вк) +

+ Б1п(у - вк)] + [/о1 соБ(вк + р) + /б2 соБ(вк + р) + /21 соб(2у - вк - р)]Х х[*ок + соб(у - вк - р) + Б1п(у - вк - р)] + [/о1 СоБ(вк - р) +

+ /б2 соБ(вк - р) + /21 соб(2у - вк + р)][*ок + 1«ш соб(у - вк + р) +

Б1П(у - вк + р)] + Ьа1уё *уё + Ьа1уд *уд + Ьа11 ^

Используя в дальнейших преобразованиях известные [5] соотношения:

соБа соБв + соБ(а - р) соБ(в - р) + соБ(а + р) соБ(в + р) = 3/2 соБ(а - в);

Б1па соБв + Б1п(а - р) соБ(в - р) + Б1п(а + р) соБ(в + р) = 3/2 Б1п(а - в);

еоБа + еоБ(а - р) + еоБ(а + р) = 0; Б1иа + Бт(а - р) + Б1и(а + р) = 0, в итоге получим:

+

^ai = id1(3/o1/2 + 3/s1/2 + 3/2i/2) cosy + iq1(3/o1/2 + 3/s1/2 + 3/2i/2) siny +...+ idN(3/oi/2 + 3/si/2 + 3I21/2) cosy + iqN(3/oi/2 + 3/si/2 + 3/21/2) siny +

+ La1ydiyd + ^aiyq^yq + Laifif.

Поскольку Ya1 = Yo + Yd1 cosy + Yq1 siny, имеем:

3 3 N 3 N

Y,i = - + - (/01 + /2, )S '* + - /,2 S U,

2 2 ,=1 2 ,=2

3 3 N 3 N

Y = - /^ + - (/01 - /2. )Siq + ^ /,2 S I',, .

2 2 j=1 2 J=2

Аналогичные преобразования уравнений потокосцепления фазной обмотки аИ дают следующие результаты:

3 3 N 3 N

Y, = -(/„ -U'di + т(/о, + /21 )Si, + -/,2 S;

2 2 J=1 2 J=2

3 3 N 3 N

Y, = -(/„ -/,2)/,i + т(/01 + /21 )Siq, + T/,2 S'* .

2 2 J=1 2 J=2

Матрица индуктивностей (индуктивных сопротивлений) приобретает

вид:

где Ldi = /oi + /si + /21;

Lqi /oi + /si

Ld1 0 (N d Ld 0 Ld 2 0

0 L„ 0 L 2 q 2 0 L 2 q 2

Ldi 0 Ld 2 0 Ld 2 0

0 L,i 0 L 2 q2 0 (N q Lq

Ldi 0 <N d Ld 0 (N d Ld 0

0 Li 0 L 2 q2 0 L 2 q 2

Ld2 /oi + /s2 : + /21;

L q2 /oi + /s2 /21.

б) Преобразования потокосцеплений для машины с разным числом витков статорных обмоток по условию минимального заполнения матрицы индуктивности

При выполнении преобразования учитывается, что индуктивность самоиндукции обмотки пропорциональна квадрату числа витков, а взаимоиндукции - пропорциональна произведению чисел витков взаи-моиндуцирующих обмоток. При этом уравнение (26) можно переписать в виде:

2

^а1 (^ а1^'а1 + ^ а1Ь1^Ь1 + ^ а1о1^о1)Ж1 + ...+ ((^ а^Ы^аЫ + ^ а^Ы^ЬЫ + а^Ысы) К1 + -^а1у6^у6 + La1yq/yq + La1fif,

где Ьга1}(] = а1, Ь1, с1, а2, ...,аЫ, ЬЫ, сЫ) = ^/(К^Ю

Аналогичные вышеприведенным преобразования дают следующие результаты

(/О1 + /';\ + /21 + (4 + /';1 + /21 )К1 I /фК] ;

1=2 _

(4 +/,; -121+(/01 + гл - /21)К I /К ;

1=2 _

* = 3 *д1 2

* =-тдЫ о

(/01+С КI /к + г, к I /ж +(/;1 - 4

1=2 1=2

Ы Ы

(/0. - с к I /к, + /;к I !,к, +(/; - 4 №.

2/

Ы дЫ

1=2

1=2

Матрица индуктивностей

1 L Iw 0 Ld2W1 0 . .. L'd Wi 0 " " Wi

0 L'qWi 0 Lqw . .. 0 L'q 2W1 Wi

Ld 2W2 0 LdW 0 . .. LdW 0 W2

[L] _ 3 _ 2 0 L'q W 0 Lqw . .. 0 L'q 2W2 X W2

Ld 2Wn 0 Ld 2WN 0 . .. LdW 0 Wn

0 L'q W 0 Lq2WN . .. 0 L'qW _ Wn

где L’d1 = /r0i + /'s1 + /Г21; L'd2 = / ol + / s2 + / 21;

Lq1 = / o + /rSl - /Г21; L'q2 = / ol + / s2 - ‘ / 21.

в) Преобразования потокосцеплений и индуктивностей при естественном взаимном расположении изображающих токов

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

Для потокосцепления фазной обмотки а1 запишем:

¥а! = ^1 + ^2 +...+ ^N,

где ^1 - потокосцепление от собственного тока обмотки W\,

¥2, ^N - потокосцепление от тока в обмотках W2 и WN;

Составляющая ^1 преобразуется по тем же уравнениям, что и в п. а. Новый характер принимают преобразования для составляющей ¥n. Рассмотрим их подробнее:

¥n = [(/o1 + ls2) COSPn + /21COS(2y - Pn)](/on + i'dnCOSy + ZqnSiny) + + [(/o1+ /s2) COS(Pn + p) + /21COS(2y - Pn - p)] [i'oN + i'dnCOS(y - p) + ^т(у - p)] + + [(lo1+ /s2) COS(Pn - p) + /210OS(2y - Pn+ p)] [i'oN + i'dnCOS(y + p) + ZqnSin(y + p)]. (31)

После преобразования уравнения (31), получим:

^n = (3/2) COSY [idn(/o1 + /s2 + /21) COSPn + iqn(/o1 + /s2 + /21) SinPn] +

+ Бту [/Д/о! + 42 - /21) СОбРп - г'яп(/о1+ 1а2 - Ы ЭШ^]

или

^п - (3/2) /Д/о1 + /я2 + /21) соэвп + (3/2) /цп(/о1 + /э2 + /21) ЭШ^

^дп (3/2) /яп(/о1 + 1э2 /21) С°8Рп (3/2) /ёп(1о1 + 1э2 /21) БШ-Рп

Матрица индуктивности в этом случае принимает вид:

Ьа 1 0 .. . Ьа 2С0sДv Ь, 2 БШ ^ .. . Ьа 2С0SД Ьа 2 Э1п А

0 ^ .. . - Ь 2«1п 3 Ь 2С0Э РМ .. . - Ь 2Э1п Рп Ь 2С0SД

Ьа 2соэЬ Ьа 2SlnД .. . Ал 0 .. . Ьа2С°5(рп ^) Ьа 2 Э1п(Рп -А )

- Ь 2«1п & Ь 2С0SД .. .0 Ь?1 .. . - Ь 2 Э1п(Рп -Аг ) Ь2 С0э(Рп -Л )

ЬЛ 2С^ Рп Ь2 Э1пРп .. . Ьа2 С0э(Рп -Л ) Ь 2 Э1п(Рп -Л ) .. . Ьл 0

- Ь 2 Э1п /Зп Ь 2С0SД .. . - Ь 2 Э1п(Рп -Р ) Ь2С°5(Рп -Р) .. .0 4,1

(32)

Применительно к обмоткам с разным числом витков, матрица индуктивности имеет вид (33) и приобретает особенности, свойственные матрице (30).

' « 0 . . Ьа2Wlc0s ДN Ьа 2 РN .. Ьа 2 С° Р„ Ьа 2 «1п РN ' ^ ■

0 Ь,1^, . . - ^2^,^ РN Ь,2 ^1С0 РN .. - Ь,2 Д Ь„2№1 С° Д

Ьа2^ С° 3N Ьа 2^ э1п Д . . ЬЛ 0 .. Ьа 2^ c0s(Д„ РN ) Ьа 2 ^ РN ) ^^ N

Ьq'2W N ^ 3N N С° Д. . .. 0 Ьд1 .. - Ь,2 ^ -ДN ) Ь,2 ^ С0S(Д„ - ДN )

Ьа С° Д Ьа 2^„ ^ Д . . Ьа 2 ^ С0^ -ДN ) Ь 2^ РN ) .. 4,1^ 0

- Э1п Р, Ьд2^п С° Д . . - Ь,2РN) Ь,2^ С0Э(Д -Р) .. .0 ь,1»; , . №п -

Ь]=

(33)

Мощность и электромагнитный момент

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

3

2

Р = Па11а1 + иЬ11Ь1 + ис1*с1 + иа 2?'а 2 + иЬ2*Ь 2 + Пс 2*с 2 + ••• +

N

+ UaNІaN + UbNІbN + UcNІcN = Е и п* п •

п=1

7=а,Ь,с

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

3 N

Р = — > и„,1А + и г + 2 ип й (34)

г* / л ап ап дпдп 0 п 0 п V /

2 п-1

Подставив в уравнение (34) значение напряжения иёп, идп, получим:

Р = | X [(л, - ^„Р7 + гпг'ап )г'ап + (Р^„ - ^а„Р7 + К„ ] (35)

где р = сИЖ - оператор дифференцирования.

Преобразование уравнения (35) перегруппировкой слагаемых дает:

N . . N

Р=Ё (р^ + р^чп *чп)+Ё ч- ^ )рг+Ё г« (/і + С) (з6)

, (**2 1 ’2

п

п=1 п=1

Уравнение (36) показывает, что потребляемая статорными обмотками активная мощность Р расходуется, как и у классической СМ, на изменение запаса электромагнитной энергии в машине, на тепловые потери и оставшаяся, основная часть, в виде электромагнитной мощности передается ротору.

Разделив электромагнитную мощность на угловую частоту вращения ротора ю = ру, получим выражение в относительных единицах для электромагнитного момента, создаваемого всеми обмотками статора:

N

т

п=1

= Ё (^Лп-^Лп) (3.101)

Уравнение движения ротора электрической машины вращательного движения имеет вид:

т - тс = Н (С2у/М2),

3.102)

где Н - механическая постоянная вращающихся масс; тс - момент сопротивления рабочей машины.

Применительно к СМ с несколькими обмотками на статоре оно приобретает вид:

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

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

3. Рациональное математическое описание синхронных машин с

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

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

Выводы

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

Литература

1. Горев А. А. Переходные процессы синхронной машины / А. А. Горев. -Л.: Госэнерго-издат, 1950. 551с.

2. Копылов И.П. Математическое моделирование электрических машин. 3-е изд., пере-раб и доп / И.П. Копылов. -М.: Высш. шк., 2001. 327с.

3. Страхов С.В. Уравнения синхронной машины с двумя трехфазными обмотками на статоре/ С.В. Страхов, В.И. Сегал. //Изв вузов. Электромеханика. 1982, № 7. С. 623-628.

4. Крылов О. А. Математическая модель синхронной машины с двумя трехфазными обмотками, сдвинутыми в пространстве / О.А. Крылов // Электричество. 1989, №3. С. 6368.

5. Кононенко Е.В. Электрические машины (Спец. курс) / Е.В. Кононенко, Г.А. Сипай-лов, К. А. Хорьков. -М.: Высш. шк., 1975. 279 с.

6. Стрижков И. Г. Математическое описание синхронного двигателя с двойной якорной обмоткой / И.Г. Стрижков // Повышение эксплуатационной надежности электрооборудования в с.х.: Труды Челяб. гос. агроинж. ун-та. 1991. С. 10-18.