Если одно из понятий незнакомо, вернитесь к предыдущему занятию и воспроизведите его численный пример. После этого переходите к новой постановке.
Место занятия в инженерной методике
Расчётная лабораторная работа · Построение карты отклика и локальных параметров
Инженерный вопрос. Зависит ли отклик одной оси от тока другой оси?
От геометрии и обмоток переходят к потокосцеплениям и ЭДС, затем к характеристикам ХХ, КЗ, нагрузки и регулирования. Эти характеристики определяют требуемое возбуждение и допустимый электрический режим.
Карты dq с перекрёстным насыщением. Выполните действия практикума по его исходной таблице. Сохраните токи d и q, оба потокосцепления и локальные производные контрольных точек.
Ответить на вопрос по своим результатам. Проверьте карту в независимых точках; сопоставьте интерполяцию или производную с прямым расчётом и повторите при меньшем шаге.
Что включить в отчёт
Результат этого занятия. Карты с демонстрацией перекрёстного насыщения и диапазона применимости.
Данные именно этого опыта. Сохраните токи d и q, оба потокосцепления и локальные производные контрольных точек.
Оформление наблюдений. Координаты каждой точки, отклик, шаг производной и способ интерполяции; контрольные точки не должны совпадать с узлами определения параметра.
Как обосновать результат
Проверка из практикума. Прямое и обратное преобразования восстанавливают фазные данные; низкотоковый момент имеет независимый ориентир.
Независимый контроль. Сопоставьте прогноз с прямым расчётом в независимой точке и измените шаг производной или плотность карты отдельно от сетки.
Критерий выполнения. Оси, нормировка и ток возбуждения указаны; получены две модели и два уровня возбуждения; шаг производной проверен.
Точность и расхождения. Разделите невязку полевого решения, ошибку интерполяции и чувствительность производной к приращению. Вблизи нулевого отклика сохраняйте абсолютную оценку с единицей.
Основание сравнения. Используйте определения карты, локальной производной и коэнергии из занятия. Вывод относится к проверенному диапазону; экстраполяция и переход через границу карты требуют отдельного обоснования.
Вывод для устройства
Объясните, в каком диапазоне карта может заменить прямые полевые расчёты в модели цепи или движения и почему за её пределами нужен новый расчёт.
Границы вывода. Каждая характеристика требует своих независимых условий и скорости. Двумерный срез, заданные токи и статическое поле не подтверждают автоматически режим на сеть, устойчивость и полную мощность реального генератора.
Дополнительный расчётный метод
Коэффициенты поля возбуждения явнополюсного генератора
Коэффициенты поля связывают геометрию полюса с основной ЭДС, полным потоком и потребной МДС возбуждения.
\tau — длина полюсного деления на окружности расточки; для неперекрывающихся полюсов необходимо b_p\le \tau. Это необходимое, но не достаточное условие геометрической пригодности.Численный пример применения метода
Для явнополюсной машины Di=6,50 м, p=24, ширина полюсного наконечника 306 мм. Определим полюсное деление и перекрытие.
Число полюсов равно 2p=48.
Полюсное деление \tau=πDi/(2p)=425{,}424 мм.
Относительное перекрытие полюсным наконечником \alpha =306/425{,}424=0{,}71928. Межполюсное расстояние по дуге равно 119,424 мм.
Ответ.\tau \approx 425{,}4 мм; \alpha \approx 0{,}7193. Эти величины задают геометрию рабочего зазора перед расчётом формы поля.
Страницы источника. PDF 1–2: определения; 3–4: модели и насыщение; 5: варианты исходных данных
Продольная и поперечная проводимости явнополюсной машины
Коэффициенты k_d и k_q объясняют различие продольного и поперечного магнитного отклика и помогают перейти к индуктивным параметрам явнополюсной машины.
Таблицы потокосцеплений описывают влияние возбуждения, токов статора и положения ротора на насыщенную синхронную машину. Продольные и поперечные составляющие позволяют исследовать режимы турбогенератора и гидрогенератора в общей системе координат.
Соглашения об обозначениях и единицах
\(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\). При постоянных магнитах и взаимном возбуждении сначала учитывают ненулевое смещение потокосцепления.
Положительные направления, оси и специальные индексы определяются в конкретной задаче. Единицу берите из условия или заголовка таблицы; перед подстановкой в формулу СИ переводите миллиметры в метры, миллигенри в генри и градусы в радианы.
Численные результаты обычно показаны с четырьмя значащими цифрами; для сравнения близких значений в таблицах сходимости сохранено больше цифр. Полные значения остаются в исходных расчётных данных. Малые остатки вычисления координат на осях показаны как ноль; малые физические величины сохраняются в показательной записи.
Продольную ось d выбирают вдоль оси поля ротора; q сдвигают на электрические 90^\circ. Для симметричного фазного набора преобразование координат выделяет составляющие, удобные для описания синхронной машины. Оно не меняет физические токи, а выражает их в другой системе. В линейном ориентире \Psi_d = L_d i_d + \Psi_{\text{возб}}, \Psi_q = L_q i_q. При насыщении \Psi_d и \Psi_q зависят от обоих токов и тока возбуждения, поэтому нужна карта и локальная матрица Якоби.
У неявнополюсного турбогенератора разница между d и q обычно меньше, чем у явнополюсного гидрогенератора; реальное отношение определяется геометрией и данными, а не названием машины. Карты строят при фиксированной температуре и выбранном токе возбуждения либо расширяют третьей осью. Из напряжений v_d = R_si_d - \omega_{\text{эл}} \Psi_q и v_q = R_si_q + \omega_{\text{эл}} \Psi_d видно, как скоростная ЭДС появляется через производную потокосцепления при вращении координат. Статическое FEMM-напряжение не заменяет эти члены.
Амплитудно-инвариантная система dq с множителем 2/3; это момент синусоидальной модели основной гармоники без явной зависимости энергии от угла при фиксированных осевых токах. Токи положительны внутрь статора. В насыщенном случае используются обе согласованные карты потокосцеплений; зубцовый момент требует угловой проверки.
Порядок проверки результата
Численный пример ниже является аналитическим учебным эталоном. Его значения вычислены по указанным формулам; они не обозначают результат решения новой модели методом конечных элементов. Практикум связывает эталон с полевой задачей: сначала воспроизведите уже проверенный результат базового раздела, затем измените только один исследуемый параметр. Именно такое сравнение позволяет отделить ошибку программы от влияния геометрии, сетки и физического приближения. Сохраните исходную модель, рабочую копию и таблицу с заголовками величин. В подписи каждой величины укажите единицу, положительное направление и способ получения; для гармонических величин дополнительно укажите амплитудную нормировку. Если условие задачи не выполнено, остановите серию и сохраните причину отказа вместе с последним корректным состоянием.
Пояснение программы
Продольное потокосцепление состоит из индуктивного слагаемого и потокосцепления возбуждения; поперечное определяется поперечной индуктивностью. Формула момента использует амплитудные токи и выбранную нормировку преобразования трёх фаз. Члены напряжения, обусловленные вращением, показаны отдельно в комментариях: неподвижное решение магнитного поля не содержит их автоматически.
Таблицы генератора и перекрёстные производные
Возбуждение и токи статора изменяются совместно; каждая таблица соответствует заданной температуре и геометрии. Пространственное приближение здесь — основная синусоидальная гармоника без явной зависимости энергии от угла при фиксированных осевых токах. Полный момент пазовой машины дополнительно проверяется угловой картой.
Зафиксируйте преобразование фазных токов в осевые. Для каждой пары осевых токов восстановите фазные токи и выполните полевое решение.
Получите оба потокосцепления и вычислите момент по общей билинейной формуле. Сравните его с полевым интегралом в выбранных контрольных состояниях.
Частные производные таблиц по двум токам формируют локальную матрицу. Ненулевые смешанные элементы показывают перекрёстное насыщение и требуют совместного решения цепи.
Потокосцепление обмотки; тот же физический смысл, что у \Psi
Вб
L
Индуктивность при явно указанной линеаризации
Гн
T
Электромагнитный момент
Н·м
p
Число пар полюсов
1
i_d,i_q
Осевые токи в выбранной амплитудной нормировке
А
Границы применимости. Статические параметры, зависящие от токов, не следует дифференцировать как константы. При использовании локальной модели с динамическими индуктивностями нужно сохранить всю локальную матрицу. Формула момента через две осевые компоненты описывает синусоидальную модель основной гармоники без явной зависимости энергии от угла при фиксированных осевых токах. Зубцовый момент и пространственные пульсации проверяют по полной угловой полевой карте.
Та же формула раскрывается как 3[\Psi_{\text{возб}} i_q + (L_d - L_q)i_d i_q], показывая возбуждённую и реактивную части линейного ориентира.
Ответ: Момент линейного аналитического эталона равен +14{,}1 \,\text{Н}\cdot \text{м}. При положительной скорости и токах внутрь статора это двигательный знак момента. Пример проверяет алгебру общей карты; для генераторного режима в той же пассивной ориентации электромагнитный момент должен тормозить ротор.
Проверка: Знак момента меняется при изменении знака i_q; число пар полюсов относится к электрическому углу.
Выберите ранее рассмотренные модели турбогенератора и явнополюсного генератора; укажите оси d/q и число пар полюсов.
Выведите обратное преобразование токов и задавайте физические фазные токи для каждой точки карты.
Во встроенном Python считывайте три фазных \Psi и преобразуйте их в выбранную d/q-нормировку.
В каждой модели постройте низкотоковый ориентир и насыщенную карту при одном токе возбуждения.
Получите локальные производные, проверяя шаг тока. Сравните зависимость L_d и L_q без принудительного приравнивания.
Добавьте второй ток возбуждения и отделите его влияние от изменения геометрии.
Ожидаемый результат. Карты и оси заданы однозначно; скоростные члены восстановлены из физических \Psi, а не из статического напряжения.
Как проверить. Прямое и обратное преобразования восстанавливают фазные данные; низкотоковый момент имеет независимый ориентир.
Обработка и анализ результатов
Зависит ли отклик одной оси от тока другой оси?
Сохраните токи d и q, оба потокосцепления и локальные производные контрольных точек.
Проверьте карту в независимых точках; сопоставьте интерполяцию или производную с прямым расчётом и повторите при меньшем шаге.
С чем сравнить. Используйте определения карты, локальной производной и коэнергии из занятия. Вывод относится к проверенному диапазону; экстраполяция и переход через границу карты требуют отдельного обоснования.
Итог работы. Карты с демонстрацией перекрёстного насыщения и диапазона применимости.
Код выполняется во встроенном Python FEMM 5. Он проверяет приведённый аналитический вывод; получение полевых данных выполняется отдельным практикумом с указанной моделью.
Python · встроенный MicroPython FEMM 5
"""Насыщение и таблицы потокосцеплений по продольной и поперечной осям генератора.Среда: встроенный MicroPython FEMM 5.Аналитический учебный пример; полевой решатель не вызывается.Параметры, единицы величин и порядок вывода пояснены по стадиям."""# Линейные продольная и поперечная индуктивности, Гн; потокосцепление возбуждения, Вб.# Токи по осям d и q заданы амплитудными значениями, А.pole_pairs,ld_h,lq_h,excitation_wb=2,0.08,0.05,1.0id_a,iq_a=-2.0,5.0# Ось d направлена вдоль возбуждения; ось q — поперёк. Потокосцепления выражены в Вб.lambda_d=ld_h*id_a+excitation_wblambda_q=lq_h*iq_a# Коэффициент 3/2 соответствует выбранному преобразованию трёх фаз и амплитудным токам.torque_nm=1.5*pole_pairs*(lambda_d*iq_a-lambda_q*id_a)# Порядок вывода: продольное потокосцепление, Вб; поперечное потокосцепление, Вб; момент, Н·м.print(lambda_d,lambda_q,torque_nm)# Скоростное напряжение добавляется отдельно:# vd = resistance*id_a - omega_electrical*lambda_q# vq = resistance*iq_a + omega_electrical*lambda_d
Пошаговая настройка расчётного примера
Откройте отдельную страницу с последовательными действиями, точными параметрами, накопительными файлами этапов и программами с русскими комментариями.
Рисунок 30.1. Визуализация численных данных FEMM 5: Синхронная машина с явнополюсным ротором. Это экспорт результатов поля, а не снимок окна программы.
Двумерная модель активной части. Модель стали с постоянной магнитной проницаемостью не описывает насыщение. Идеализированные обмотки с нулевой проводимостью не учитывают все потери. Сеточная проверка момента этой рабочей точки выполнена отдельно; её критерии, методы и границы приведены в протоколе. Токи статора заданы как мгновенное состояние. Это поле при токовом возбуждении, а не автоматически согласованная нагрузка генератора. Скоростную ЭДС восстанавливают отдельно; напряжение неподвижного решения её не содержит. У острых кромок модели стали с постоянной проницаемостью возможны большие локальные значения B. Они не описывают реальный материал и сами по себе не доказывают насыщение или повреждение. Проверяйте нелинейную характеристику материала и сходимость поля в физически заданной области.
Ориентиры в интерфейсеРисунок 30.2. Геометрия: цилиндрический ротор с токовым возбуждением; неподвижный снимок, скоростная ЭДС отдельно Открыть в полном размереРисунок 30.3. Распределение магнитной индукции: цилиндрический ротор с токовым возбуждением; неподвижный снимок, скоростная ЭДС отдельно Открыть в полном размереРисунок 30.4. Геометрия: явнополюсный ротор с токовым возбуждением; отдельная геометрия Открыть в полном размереРисунок 30.5. Распределение магнитной индукции: явнополюсный ротор с токовым возбуждением; отдельная геометрия Открыть в полном размереРисунок 30.6. Синхронная машина с явнополюсным ротором: выбранная сетка воздушного зазора. Открыть в полном размереРисунок 30.7. Синхронная машина с явнополюсным ротором: момент ротора на выбранной сетке. Открыть в полном размереРисунок 30.8. Встроенная консоль FEMM: Насыщение и таблицы потокосцеплений по продольной и поперечной осям генератора. Аналитическая проверка программы занятия. Показана консоль исходной редакции программы, зафиксированной при съёмке. Текущая программа с русскими пояснениями проверена отдельно во встроенном Python FEMM; этот снимок не является новым запуском текущей редакции. Открыть в полном размере
Самостоятельная работа
Сопоставьте локальные d/q-матрицы турбо- и гидрогенератора в низкотоковой и насыщенной точках; объясните различие.
Критерий выполнения: Оси, нормировка и ток возбуждения указаны; получены две модели и два уровня возбуждения; шаг производной проверен.