FEMM 5

Часть V · Электрические машины → § 24 · Синхронные генераторы

Занятие 24.2

Характеристика ХХ и возбуждение турбогенератора

Цель и задачи

Цель. Получить характеристику холостого хода из серии возбуждений.

Задачи

  • Получить основную гармонику потокосцепления по углам.
  • Построить E_0(I_{\text{в}}) при неизменной скорости.
Что нужно вспомнить перед началом

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

Место занятия в инженерной методике

Расчётная лабораторная работа · Восстановление электрических характеристик синхронного генератора

Инженерный вопрос. Где характеристика отклоняется от линейного роста?

От геометрии и обмоток переходят к потокосцеплениям и ЭДС, затем к характеристикам ХХ, КЗ, нагрузки и регулирования. Эти характеристики определяют требуемое возбуждение и допустимый электрический режим.

Последовательность решения, отчёт и критерии результата

Результат получают из наблюдаемых величин и проверяют в условиях данного практикума.

Последовательность решения

  1. Обосновать постановку и способ измерения. Для характеристики запишите, что фиксировано: скорость, возбуждение, напряжение, нагрузка или мощность. Скоростная ЭДС восстанавливается из пространственного потокосцепления.

  2. Насыщающаяся характеристика ХХ. Проверьте уравнение выбранного генераторного режима и знаки мощности до построения серии.

  3. Ручная характеристика ХХ. Выполните действия практикума по его исходной таблице. Сохраните возбуждение и вычисленную ЭДС каждой точки при одной скорости.

  4. Ответить на вопрос по своим результатам. Сведите потокосцепления, фазоры и мощности к одному режиму; проверьте независимые условия и выбранные направления.

Что включить в отчёт

Результат этого занятия. Характеристика ХХ с объяснением линейной области и насыщения.

Данные именно этого опыта. Сохраните возбуждение и вычисленную ЭДС каждой точки при одной скорости.

Оформление наблюдений. Возбуждение, скорость, фазные потокосцепления, ЭДС и согласованные токи/напряжения в объёме конкретного опыта; независимые условия каждой строки.

Как обосновать результат

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

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

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

Точность и расхождения. Уточняйте шаг возбуждения или угла отдельно от сетки. При КЗ оценивайте абсолютное остаточное напряжение и его отношение к ненулевому масштабу режима.

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

Вывод для устройства

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

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

Дополнительный расчётный метод

Магнитные поля в электрических машинах — введение (файл Гончарова)

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

Лабораторные работы по методу: Трансформатор от потокосцеплений до схемы замещения и нагрузочной характеристики; Асинхронный двигатель от поля в зазоре до параметров и механической характеристики; Генератор с постоянными магнитами от конструкции ротора до ЭДС и нагрузочного напряжения; Выбор обмотки и расчет момента в кольцевом зазоре; Определение коэффициентов поля явнополюсной машины; Асинхронный двигатель: механическая характеристика и вращающееся поле; Синхронный генератор: от полевых параметров к внешней и угловой характеристикам; Потери массивного ротора: частота, гармоники и магнитное смещение; Зубцовые контуры: проводимости, коронки и взаимная связь при смещении; Проводник на переменном токе: сопротивление, транспозиция и частота силового воздействия; Замкнутый стальной виток: проникновение поля, сопротивление и потери; Контрольный опыт: от условия задачи к инженерному выводу; Измерительная рамка: от поля провода к ЭДС; От магнитной цепи к карте \Psi (i,\theta ) и моменту SRM; Зубцовый момент: пусковая помеха и спектр PM-генератора; Нагрузка PM-генератора и рабочая область магнитов

Получаемые величины. \Phi =\Delta Az·depth, B=\Delta Az/\Delta n; L, M; Bp и ошибка обратного восстановления; kf, k_{\Phi}, kd, kq; E0(If), Xd, Xq, U, I,\varphi,\delta; Jротора(f2), T(s), отличие дефектного ротора

Страницы источника. PDF1–75; A_z/поле 1–16; L_{17}–20; FFT20–29; синхронные 30–60; AD60–74

Литература по методу

Коэффициенты поля возбуждения явнополюсного генератора

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

Лабораторные работы по методу: Определение коэффициентов поля явнополюсной машины

Получаемые величины. k_f = B_{1}/B полюс; k_\Phi = Φполный/\Phi 1; \alpha \delta = B ср/B полюс; коэффициент насыщения k_z

  1. Рассчитываем ненасыщенную геометрию при заданном токе возбуждения.
  2. Интегрируем поток на полюсном делении и разложить B по гармоникам.
  3. Вводим зубцы, ярмо, проверенную B(H) и заполнение сталью.
  4. Сравниваем коэффициенты и требуемую МДС при одном целевом потоке.
\tau=\frac{\pi D_i}{2p},\qquad \alpha_{\mathrm{геом}}=\frac{b_p}{\tau}
\tau — длина полюсного деления на окружности расточки; для неперекрывающихся полюсов необходимо b_p\le \tau. Это необходимое, но не достаточное условие геометрической пригодности.
Численный пример применения метода

Для явнополюсной машины Di=6,50 м, p=24, ширина полюсного наконечника 306 мм. Определим полюсное деление и перекрытие.

  1. Число полюсов равно 2p=48.
  2. Полюсное деление \tau=πDi/(2p)=425{,}424 мм.
  3. Относительное перекрытие полюсным наконечником \alpha =306/425{,}424=0{,}71928. Межполюсное расстояние по дуге равно 119,424 мм.

Ответ. \tau \approx 425{,}4 мм; \alpha \approx 0{,}7193. Эти величины задают геометрию рабочего зазора перед расчётом формы поля.

Страницы источника. PDF 1–2: определения; 3–4: модели и насыщение; 5: варианты исходных данных

Литература по методу

Теория и вывод формул

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

Соглашения об обозначениях и единицах
  • \(i(t)\), \(u(t)\) — мгновенные значения; \(\widehat I\), \(\widehat U\) — комплексные амплитуды; \(\underline I_{\text{д}}\), \(\underline U_{\text{д}}\) — комплексные действующие значения. Для синусоиды \(\widehat I=\sqrt{2}\,\underline I_{\text{д}}\); действующее значение тока равно \(|\underline I_{\text{д}}|\). Нормировка обязательно учитывается в формулах мощности.
  • \(\Phi=\int_S\mathbf B\cdot\mathrm d\mathbf S\) — поток через ориентированную поверхность; \(\Psi=\sum_k\oint_{C_k}\mathbf A\cdot\mathrm d\mathbf l\) — потокосцепление последовательных витков с согласованными замкнутыми контурами. Только при одинаковом сцепленном потоке каждого витка \(\Psi=N\Phi\). Обе величины измеряются в веберах: число витков безразмерно. Запись «Вб·виток» подчёркивает суммирование.
  • \(\theta_{\text{м}}\), \(\Omega\) (также \(\omega_{\text{мех}}\)) — механические угол и угловая скорость. Для основной гармоники \(\theta_{\text{э}}=p\theta_{\text{м}}\), \(\omega_{\text{э}}=p\Omega\), где \(p\) — число пар полюсов. Геометрический период зубцов задаётся отдельно. В производных по углу используйте радианы.
  • Статическая индуктивность \(L_{\text{ст}}=\Psi/i\) — наклон секущей из начала координат; это определение используется при \(\Psi(0)=0\) и \(i\ne0\). Динамическая индуктивность \(L_{\text{дин}}=\partial\Psi/\partial i\) — наклон касательной (дифференциальная, инкрементальная); положение и остальные токи фиксированы. Для линейной характеристики без смещения обе равны одной индуктивности \(L\). При постоянных магнитах и взаимном возбуждении сначала учитывают ненулевое смещение потокосцепления.
  • Положительные направления, оси и специальные индексы определяются в конкретной задаче. Единицу берите из условия или заголовка таблицы; перед подстановкой в формулу СИ переводите миллиметры в метры, миллигенри в генри и градусы в радианы.
  • Численные результаты обычно показаны с четырьмя значащими цифрами; для сравнения близких значений в таблицах сходимости сохранено больше цифр. Полные значения остаются в исходных расчётных данных. Малые остатки вычисления координат на осях показаны как ноль; малые физические величины сохраняются в показательной записи.
Указатель формул этого занятия
  1. (24.4) · Последовательность двух серий
  2. (24.5) · Основная гармоника и производная

Последовательность двух серий

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

В начальном линейном диапазоне \Psi_{\text{амп}} пропорционально току возбуждения. При росте насыщения зубцов, ярма и ротора прирост потокосцепления на один ампер уменьшается. Характеристика ХХ описывает магнитное состояние и связь с обмоткой, но не задаёт полноценные переходные параметры генератора. Остаточная магнитная индукция возникает только при подходящей модели материала и истории, а не из произвольного нулевого тока.

E_0(I_{\text{в}})=\frac{2\pi f}{\sqrt2}\Psi_{1,\text{амп}}(I_{\text{в}})
(24.4)
Нижний индекс 1 обозначает основную пространственную гармонику фазного потокосцепления; f соответствует заданным p и n.

Основная гармоника и производная

При редкой ручной серии полезно сначала оценить амплитуду близкой синусоиды по максимумам и минимумам. Более строгая проверка использует коэффициенты синуса и косинуса по полному электрическому периоду. Производная усиливает высокие гармоники: компонент \Psi_h \cos(h \Theta) даёт ЭДС с множителем h \omega_{\text{э}}. Поэтому грубая угловая таблица может хорошо передавать основной поток, но плохо передавать форму ЭДС.

В практикуме 12\text{–}24 положения достаточно для освоения основной характеристики, но это не доказательство точности всех пазовых гармоник. При переходе к нелинейному материалу ток возбуждения сохраняют одинаковым во всех положениях, а полевой расчёт повторяют заново. Сетку особенно проверяют в зубцах и насыщенных перемычках.

Фактически приложенный цилиндрический образец решён магнитостатически: ток цепи «Обмотка возбуждения» = 2 \,\text{А}, мгновенные фазные токи A = 1 \,\text{А}, B = 1 \,\text{А}, C = -2 \,\text{А}; глубина 80 \,\text{мм}. Это один расчёт при фиксированных токах с линейной сталью, а не завершённая характеристика ХХ, КЗ или согласованная нагрузка генератора. Скорость, угол нагрузки и скоростная ЭДС в этом снимке не вычислены. Серии и цепное согласование практикума выполняются отдельно.

\Psi_a=\sum_{h\ge1}\Psi_{h,\text{амп}}\cos(h\Theta+\varphi_h)\ \Rightarrow\ e_a=\sum_{h\ge1}h\omega_{\text{э}}\Psi_{h,\text{амп}}\sin(h\Theta+\varphi_h)
(24.5)
Вывод непосредственно использует \mathrm d\Theta/\mathrm dt = \omega_{\text{э}}; знак зависит от начального фазового направления, которое фиксируется в отчёте.

От насыщения поля возбуждения к характеристике холостого хода

Скорость и схема обмотки сохраняются; изменяется только ток возбуждения. Напряжение сравнивается на одной стороне и для одной выбранной гармоники.

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

  2. Примените закон Фарадея при неизменной скорости. Нелинейность зависимости ЭДС от тока возбуждения отражает насыщение магнитной цепи.

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

Таблица 24.4. Обозначения и единицы
ВеличинаФизический смыслЕдиница СИ
IТок обмотки; в статике постоянныйА
\PsiПотокосцепление обмотки; тот же физический смысл, что у \PsiВб
fЧастота гармоникиГц
U_{\text{д}}Действующее значение напряженияВ
BМагнитная индукцияТл

Границы применимости. Увеличение основной гармоники не описывает рост всех гармоник одинаково. Экстраполяция материала за последние табличные точки не подтверждается одной гладкой характеристикой ЭДС.

Разобранный численный пример

Насыщающаяся характеристика ХХ

Дано. Учебный аналитический закон \Psi_{1,\text{амп}} = 0{,}30 \tanh(I_{\text{в}}/4) Вб; f = 50 \,\text{Гц}.

  1. При I_{\text{в}} = 2 \,\text{А} получаем \Psi_{\text{амп}} = 0{,}138635 \,\text{Вб} и E_0 = 30{,}79699 \,\text{В}.

  2. При 4 \,\text{А}: 0{,}228478 \,\text{Вб} и 50{,}75511 \,\text{В}. При 6 \,\text{А}: 0{,}271544 \,\text{Вб} и 60{,}32202 \,\text{В}.

  3. Двукратный рост тока от 2 до 4 \,\text{А} повышает ЭДС только в 1,648 раза; дальнейший рост от 4 до 6 \,\text{А} даёт ещё меньшее относительное увеличение.

Ответ: Предел ЭДС данной аналитической модели при большом токе равен 66{,}64324 \,\text{В}.

Проверка: При малом I_{\text{в}} наклон E_0^{\prime} = (2\pi\cdot50/\sqrt{2})\cdot 0{,}30/4 = 16{,}66081 \,\text{В}/\text{А}. Это математический эталон, не материал или результат реальной машины.

Практикум в FEMM 5

Ручная характеристика ХХ

Таблица 24.5. Исходные данные практикума
Исходная величинаЗначениеЕдиница
Радиус расточки / наружной поверхности статора60 / 100мм
Пазы: число121
Пазы: диапазон радиусов60–85мм
Пазы: угловая ширина14град
Активная длина80мм
Число витков каждой стороны со знаком40витков
Пары полюсов11
Ротор / вал58 / 18мм
Катушки возбуждения: витки сторон со знаком\pm 200витков
Ток возбуждения2А
  1. Базовый статор: расточка радиуса 60 \,\text{мм}, наружный радиус 100 \,\text{мм}, глубина 80 \,\text{мм}. Создайте 12 пазов между радиусами 60 и 85 \,\text{мм} угловой шириной 14^\circ с шагом 30^\circ. Наружный воздушный квадрат x,y = \pm 170 \,\text{мм} замкните условием A = 0.
  2. Материалы первого эталона: воздух и идеализированная обмоточная область \mu_r = 1, \sigma = 0; линейная сталь \mu_r = 1000, \sigma = 0. Каждый паз имеет 40 витков с учётом знака. Сопротивление реальной меди здесь не рассчитано; оно задаётся отдельным цепным параметром.
  3. Двухполюсная обмотка: последовательность пазов от 0^\circ — A+, A+, C-, C-, B+, B+, A-, A-, C+, C+, B-, B-. Все стороны имеют по 40 витков с учётом знака. Здесь каждая фазная полоса занимает два соседних паза.
  4. Ротор — кольцо железа от внутреннего вала радиуса 18 \,\text{мм} до наружного радиуса 58 \,\text{мм}. Воздушный зазор до расточки равен 2 \,\text{мм}. Катушечные стороны возбуждения — кольцевые сектора радиусов 30\text{–}55 \,\text{мм}, по 30^\circ ширины, с центрами 90^\circ и 270^\circ; назначьте цепь «Обмотка возбуждения», +200 и -200 витков.
  5. Задайте ток возбуждения 2 \,\text{А}. Ротору и его обмоточным областям назначьте общую группу. Для ручного поворота вращайте всю группу; поле возбуждения не заменяется магнитом. Исходная угловая точка определяется как совмещение оси поля с выбранной осью статора.
  6. Выберите токи возбуждения 1; 3; 5 \,\text{А} и нулевые фазные токи. Для каждого выполните 12 положений на одном электрическом периоде, не повторяя конечную точку.
  7. Для каждой серии запишите фазные потокосцепления. Получите основную амплитуду, отдельно сохраните максимумы и минимумы.
  8. При n = 3000 \,\text{об}/\text{мин} вычислите E_0. Сравните исходный линейный вариант и отдельную копию с учебной B\text{–}H-таблицей из 19.5.
  9. Для одной точки удвойте число положений и проверьте основную гармонику. Перепроверьте локальную индукцию на более мелкой сетке.

Ожидаемый результат. Линейная характеристика имеет постоянный наклон; нелинейная отклоняется от него по мере насыщения. Угловая зависимость связывает магнитное поле с ЭДС.

Как проверить. Проверьте полноту периода, механический/электрический угол и одинаковую скорость во всех точках. Гармоническая характеристика не называется временной петлёй или ударным процессом.

Обработка и анализ результатов

Где характеристика отклоняется от линейного роста?

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

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

Итог работы. Характеристика ХХ с объяснением линейной области и насыщения.

Пошаговая настройка расчётного примера

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

Все расчётные примеры курса

Расчётный образец FEMM 5

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

Таблица 24.6. Результаты расчётного образца
ВеличинаРезультатЕдиница
Частота0Гц
Активная длина80мм
Ток цепи Фаза А1А
Потокосцепление цепи Фаза А0,03417Вб
Ток цепи Фаза В1А
Потокосцепление цепи Фаза В0,05922Вб
Ток цепи Фаза С-2А
Потокосцепление цепи Фаза С-0,1016Вб
Ток цепи Обмотка возбуждения2А
Потокосцепление цепи Обмотка возбуждения1,186Вб
Момент на исходной учебной сетке-0,04961Н·м
Энергия поля1,334Дж
Модуль B в контрольной точке0,1528Тл
Координаты контрольной точки (x, y)56,42; 17,25мм
Момент на проверенной уточнённой сетке-0,04951Н·м
Визуализация численных данных FEMM 5: Синхронная машина с неявнополюсным ротором. Это экспорт результатов поля, а не снимок окна программы.
Рисунок 24.6. Визуализация численных данных FEMM 5: Синхронная машина с неявнополюсным ротором. Это экспорт результатов поля, а не снимок окна программы.

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

Ориентиры в интерфейсе
Геометрия: цилиндрический ротор с токовым возбуждением; неподвижный снимок, скоростная ЭДС отдельно
Рисунок 24.7. Геометрия: цилиндрический ротор с токовым возбуждением; неподвижный снимок, скоростная ЭДС отдельно Открыть в полном размере
Распределение магнитной индукции: цилиндрический ротор с токовым возбуждением; неподвижный снимок, скоростная ЭДС отдельно
Рисунок 24.8. Распределение магнитной индукции: цилиндрический ротор с токовым возбуждением; неподвижный снимок, скоростная ЭДС отдельно Открыть в полном размере
Синхронная машина с неявнополюсным ротором: выбранная сетка воздушного зазора.
Рисунок 24.9. Синхронная машина с неявнополюсным ротором: выбранная сетка воздушного зазора. Открыть в полном размере
Синхронная машина с неявнополюсным ротором: момент ротора на выбранной сетке.
Рисунок 24.10. Синхронная машина с неявнополюсным ротором: момент ротора на выбранной сетке. Открыть в полном размере

Самостоятельная работа

Постройте две характеристики ХХ — линейную и нелинейную — и объясните магнитные области, определяющие изгиб.

Критерий выполнения: ЭДС восстановлена из угловых потокосцеплений, указаны ток возбуждения и скорость каждой серии, основная гармоника проверена сгущением угла.
Первоисточники и дополнительное чтение по занятию
  1. [33] — 24–39; 45–52
  2. [2] — Программа: трансформаторы, индукционные, синхронные машины и машины постоянного тока
  3. Н. Бьянки. PM motor design with FEMM — скачанный отдельный раздел, физика и постобработка — 68–71: механический и электрический углы; старый интерфейс и Lua не назначаются обязательным инструментом
Пояснения и типичные ошибки
Почему E_0(I_{\text{в}}) нельзя получать на разных скоростях?

Скорость прямо множит ЭДС. Изменение наклона тогда смешает механику и насыщение.

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

Она умножает амплитуду h-й гармоники на h; малые высокие компоненты потокосцепления могут заметно влиять на ЭДС.

На что обратить внимание

  • Изменение скорости одновременно с током возбуждения.
  • Амплитуда потокосцепления делится на \sqrt{2} дважды.
  • Редкая серия названа точным спектром пазовой ЭДС.