FEMM 5

Часть V · Электрические машины → § 37 · Турбогенераторы: от параметрической геометрии до характеристик, динамики и дефектов

Занятие 37.11

Характеристика холостого хода и линия воздушного зазора

Цель и задачи

Цель. Оценить насыщение по характеристике холостого хода генератора.

Задачи

  • Построить характеристику ХХ по основной гармонике.
  • Определить линию воздушного зазора по объявленному линейному опыту.
Что нужно вспомнить перед началом

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

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

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

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

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

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

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

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

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

  2. Контрольный аналитический пример. Проверьте уравнение выбранного генераторного режима и знаки мощности до построения серии.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Получаемые величины. 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: варианты исходных данных

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

Применение конечно-элементных моделей при учебном проектировании синхронных машин

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

Лабораторные работы по методу: Трансформатор от потокосцеплений до схемы замещения и нагрузочной характеристики; Асинхронный двигатель от поля в зазоре до параметров и механической характеристики; Генератор с постоянными магнитами от конструкции ротора до ЭДС и нагрузочного напряжения; Выбор обмотки и расчет момента в кольцевом зазоре; Асинхронный двигатель: механическая характеристика и вращающееся поле; Синхронный генератор: от полевых параметров к внешней и угловой характеристикам

Получаемые величины. основная ЭДС и её фазовый угол; фазные \Psi; U, I, \cos \varphi при согласованной нагрузке; поправки длины и насыщения

  1. Проверить обмоточную таблицу и период модели.
  2. В одном режиме получить гармонику поля и фактические потокосцепления.
  3. Восстановить ЭДС обоими способами с общей скоростью и нормировкой.
  4. Согласовать внешнюю цепь и сравнить прогноз на независимой точке нагрузки.
  5. Объяснить влияние неплоской конструкции и эффективной длины.

Страницы источника. PDF 1–14 полностью: 1–5 модели и определение ЭДС; 6–9 таблицы и выводы; 10–11 литература и геометрия; 12–14 английское резюме и авторы.

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

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

Характеристика E0(If) строится при неизменной частоте вращения и нулевых токах статора. Сопоставление нелинейной зависимости с линейной оценкой раскрывает насыщение и влияние принятой учебной характеристики B(H).

Соглашения об обозначениях и единицах
  • \(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. (37.97) · Линейная опора и коэффициент насыщения
  2. (37.98) · Вывод характеристики ХХ, линейной опоры и интерполяции
  3. (37.99) · Вывод характеристики ХХ, линейной опоры и интерполяции
  4. (37.100) · Вывод характеристики ХХ, линейной опоры и интерполяции
  5. (37.101) · Вывод характеристики ХХ, линейной опоры и интерполяции
  6. (37.102) · Вывод характеристики ХХ, линейной опоры и интерполяции

Условия опыта ХХ

Обмотка статора разомкнута: iA=iB=iC=0. Изменяется только ток возбуждения, частота и геометрия неизменны. Для каждой точки нужна основная гармоника потокосцепления или эквивалентная проверенная процедура её извлечения; одно случайное мгновенное значение \Psi не является амплитудой.

Одноценная B(H) без гистерезиса не описывает остаточный магнетизм, восходящую и нисходящую ветви реального опыта. Нулевая ЭДС при If=0 — свойство этой модели, а не опровержение остаточного напряжения машины.

Линейная опора и коэффициент насыщения

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

Коэффициент насыщения ksat=If,actual/If,linear при заданном E показывает дополнительную потребность в возбуждении. Отношение вертикальных ординат при одном If имеет иной смысл. При интерполяции избегайте экстраполяции за рассчитанный диапазон.

E_0(I_f)=\frac{\omega_e\Psi_1(I_f)}{\sqrt2},\qquad k_{sat}(E_*)=\frac{I_{f,\mathrm{нелин}}(E_*)}{I_{f,\mathrm{линия}}(E_*)}
(37.97)
If — А, E_{0} — действующая фазная ЭДС В; скорость одинакова для обеих кривых.
  1. При E*=230 В нелинейная модель требует 6 А, линия зазора — 4 А.

  2. Получаем ksat=6/4=1{,}5; это горизонтальное сравнение при общей ЭДС.

Согласованность серии

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

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

Вывод характеристики ХХ, линейной опоры и интерполяции

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

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

Связь с программой. В шаге 10 operating() выполняется при current=0 для ряда If. Первая гармоника определяется усреднённым пространственным вектором; выбранный малотоковый наклон можно проверить отдельной таблицей или линейным вариантом B(H).

  1. В линейном уравнении поля K a=If f1 матрица K постоянна; следовательно a и линейный функционал \Psi пропорциональны If.

    \mathbf a=I_f\mathbf K^{-1}\mathbf f_1\quad\Longrightarrow\quad\Psi_1=c_\Psi I_f
    (37.98)
  2. Переход к действующей ЭДС даёт линейную опору с наклоном kE, пропорциональным скорости.

    E_{gap}=k_E I_f,\qquad k_E=\frac{\omega_e c_\Psi}{\sqrt2}
    (37.99)
  3. Нелинейная характеристика определяется полем каждой точки. При обратной интерполяции целевая ЭДС должна находиться между двумя рассчитанными ординатами.

    I_f(E_*)\approx I_{f,1}+\frac{E_*-E_1}{E_2-E_1}(I_{f,2}-I_{f,1}),\quad E_1\le E_*\le E_2
    (37.100)
  4. Для линейной опоры If,gap=E*/kE; деление нелинейного тока на этот ток даёт коэффициент насыщения при общей ЭДС.

    k_{sat}(E_*)=\frac{k_E I_{f,nonlinear}(E_*)}{E_*}
    (37.101)
  5. При неизменном магнитостатическом поле \Psi не зависит от назначенной кинематической скорости; ЭДС масштабируется частотой.

    \frac{E_0(f_2,I_f)}{E_0(f_1,I_f)}=\frac{f_2}{f_1},\qquad\Psi_1(f_2,I_f)=\Psi_1(f_1,I_f)
    (37.102)
Таблица 37.32. Обозначения и единицы
ВеличинаФизический смыслЕдиница СИ
c_\PsiЛинейный коэффициент потокосцепления по IfГн
k_EНаклон линейной ЭДСВ/А
I_{f,nonlinear}Возбуждение по насыщенной ХХА
k_{sat}Объявленный коэффициент насыщения1
E_0,E_*ЭДС ХХ и общий выбранный уровень сравненияВ действующее
\widehat\Psi_1Амплитудное потокосцепление основной гармоникиВб
I_fПостоянный ток возбужденияА
fЧастота движения поляГц

Границы применимости. Пропорциональность частоте относится к статическому полю без скоростных вихретоковых эффектов.

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

Контрольный аналитический пример

Дано. При E*=230 В нелинейная характеристика требует If=6 А, линия воздушного зазора — If=4 А.

  1. Сравниваются токи при одной ЭДС и скорости.

  2. ksat=6/4=1{,}5.

Ответ: Коэффициент насыщения ksat=1,5 в данном аналитическом примере.

Проверка: Это самостоятельный аналитический эталон. Сравнение с FEMM допустимо при совпадении всех входов, геометрии и определения величины.

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

Шаг 10: Характеристика холостого хода и линия воздушного зазора

Таблица 37.33. Исходные данные практикума
Исходная величинаЗначениеЕдиница
Учебная модельдвухполюсный цилиндрический ротор; 24 открытых паза статора, 12 открытых пазов ротора; отдельные немагнитные клинья
СредаFEMM 5 · встроенный MicroPython · магнитостатика позиций
Частота изменения положения50Гц
  1. Выполните run-10.py; проверьте нулевые токи статора во всех точках.
  2. Постройте E0(If) по первой гармонике, отдельно подпишите фазное и линейное напряжение.
  3. Укажите способ получения линии зазора; при необходимости выполните отдельную линейную серию.
  4. Выберите E внутри общего диапазона двух кривых и определите два If интерполяцией.
  5. Повторите одну насыщенную точку с более плотным угловым набором и сеткой.

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

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

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

Где полевая характеристика отклоняется от линии воздушного зазора?

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

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

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

Встроенный MicroPython: Характеристика холостого хода и линия воздушного зазора · Порядок запуска

Выполнить независимый шаг 10

Сохраните reproduce.py и выбранный run-NN.py в одну рабочую папку и прочитайте README.md. Укажите её абсолютный путь в TG_ROOT; пример на странице использует C:/FEMM5/turbogenerator. Запускайте через редактор Python FEMM 5. Каждый шаг восстанавливает необходимые входы без запуска предыдущих; внешний pyFEMM не требуется.

Python · встроенный MicroPython FEMM 5
"""Характеристика холостого хода и линия воздушного зазора.

Встроенный MicroPython FEMM 5; отдельный воспроизводимый шаг учебного турбогенератора.
Результат относится к адаптированной геометрии, а не к промышленной машине 340 МВт.
"""
# Укажите рабочий каталог; сохраните рядом reproduce.py и прочитайте README.md.
TG_ROOT = r"C:/FEMM5/turbogenerator"
STAGE = "10"
CASE = "standard"
# Программа заново строит нужные предпосылки и сохраняет модели, таблицы и квитанцию шага.
exec(open(TG_ROOT + "/reproduce.py").read(), globals())

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

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

Таблица 37.34. Результаты расчётного образца
ВеличинаРезультатЕдиница
Рассчитанные точки характеристики61
Максимальное напряжение ХХ в серии208,7В
График сохранённых численных данных учебного образца; условия и формулы приведены в тексте. Это самостоятельная научная визуализация, а не снимок интерфейса FEMM.
Рисунок 37.7. График сохранённых численных данных учебного образца; условия и формулы приведены в тексте. Это самостоятельная научная визуализация, а не снимок интерфейса FEMM.

Нелинейная 2D магнитостатика активной части и отдельная редуцированная динамика. Проверка сетки относится к напряжению ХХ. Уточнение углового шага и баланс мощности приведены отдельно на этапе 15.

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

Получите характеристики при двух скоростях и проверьте пропорциональность ЭДС частоте при одинаковых If и магнитостатических входах; насыщение сравнивайте по \Psi 1.

Критерий выполнения: Токи статора равны нулю, скорость фиксирована, ЭДС действующая фазная; линия зазора имеет явное определение, интерполяция не выходит за диапазон.
Первоисточники и дополнительное чтение по занятию
  1. [65] — С. 1–12; ХХ, КЗ, регулировочные и U-образные характеристики; в статье рассмотрен турбогенератор 325 МВт.
  2. [105] — С. 26–32; Угловые функции, гармоники и фазовые соотношения; нормировка курса определена отдельно.
  3. [33] — Магнитная плоская задача, цепи mi_/mo_, встроенный Python и постпроцессор.
Пояснения и типичные ошибки
Почему в модели нет восходящей и нисходящей ветвей ХХ?

Одноценная безгистерезисная B(H) не хранит предысторию намагничивания.

Можно ли вычислять коэффициент насыщения как отношение любых двух ординат?

Нет. Приведённое определение сравнивает токи при одной и той же ЭДС; другое отношение нужно назвать и определить отдельно.

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

  • Изменение скорости не изменяет магнитостатическое поле при тех же токах в данной постановке.
  • Безгистерезисная ХХ не является экспериментальной петлёй реальной машины.