FEMM 5

Лабораторная работа

Пазовый зазор: коэффициент Картера и комплексная передача гармоник

Цель. Определить добавочное магнитное сопротивление, обусловленное раскрытием паза, и коэффициент Картера. При одинаковом возбуждении оценить изменение амплитуд и фаз пространственных гармоник.

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

Задачи

  1. Сверить гладкий зазор с решением уравнения Лапласа.
  2. Получить k\delta для четырёх раскрытий и отделить сеточную погрешность.
  3. Сопоставить комплексные коэффициенты одинаковых гармоник пазовой и гладкой моделей.

Исходные данные и условия опыта

Плоская задача, активная длина 1 м, зубцовое деление 40 мм, зазор 2 мм, глубина паза 20 мм. Раскрытие 0, 4, 8 или 12 мм. Воздух имеет \mu _{r}=1; поверхности идеальных сердечников задают H\tau =0.

В серии Картера задан средний поток Bср=0,8 Тл через постоянные значения A на боковых сечениях. В гармонической серии одинаковый токовый слой толщиной 0,3 мм расположен у верхнего сердечника; J_{1}=10 А/мм², J_{3}/J_{1}=0{,}3∠0,6 рад. Проводимость равна нулю, поэтому частота 50 Гц задаёт фазоры без вытеснения тока.

Теория расчётного метода

Почему раскрытие увеличивает требуемую МДС

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

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

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

W=\frac{B_{\rm ср}^{2}\tau\delta_{\rm эфф}l}{2\mu_0},\qquad k_\delta=\frac{2\mu_0W}{B_{\rm ср}^{2}\tau l\delta}
W — статическая энергия всего воздушного зазора и паза; размеры в метрах.
k_{\delta,\rm conf}=\frac{\tau}{\tau-\frac{4\delta}{\pi}[\rho\arctan\rho-\frac12\ln(1+\rho^2)]},\quad \rho=\frac{b_0}{2\delta}
Конформная оценка для прямоугольного достаточно глубокого раскрытия напротив гладкого сердечника; влияние соседних раскрытий и конечной глубины проверяют полем.

Литература по методу: Гладкий зазор: аналитическая проверка поля и независимых источников ошибки · PDF 1–4: границы; 5–9: слои и периодичность; 9–12: количественная программа; Чётное и нечётное поля паза: Картер и разностная проводимость · PDF 1: смысл разложения; 2–3: варианты, интегралы, k\delta и проверка; Чётное поле паза и контроль комплексной потенциальной функции · PDF 1–3: отображения; 4–6: силовая линия; 6–8: потенциалы и программа

Гармоника сохраняет амплитуду и фазу

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

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

Дробление источника и дробление конечных элементов — разные приближения. Поэтому сначала сравнивают 24 и 48 участков токового слоя, затем отдельно уточняют сетку.

\widehat B_{y,\nu}(y)=\frac{\mu_0\widehat J_\nu\sinh(k_\nu t)\cosh(k_\nu y)}{k_\nu\sinh(k_\nu\delta)},\quad k_\nu=2\pi\nu/\tau
Модуль гармоники в области ниже верхнего слоя, 0\le y\le \delta -t. Пространственная производная задаёт синусную фазу B относительно косинусного J.
\widehat C_\nu=\widehat B_{\nu,\rm паз}/\widehat B_{\nu,\rm гладк}
Сравниваются одинаковый порядок, координата выборки, шаг и фазовая система отсчёта.

Литература по методу: Влияние раскрытия пазов на гармоники воздушного зазора · PDF 1: допущения; 2–4: расчётные выражения; 5–6: парные модели; 7: первичная литература; Гладкий зазор: аналитическая проверка поля и независимых источников ошибки · PDF 1–4: границы; 5–9: слои и периодичность; 9–12: количественная программа; Нечётное поле паза как аналитический эталон рассеяния · PDF 1–3: постановка и формула; 4–7: модели и сохранение величин; 7–8: программа

Уравнение поля, материалы и границы

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

На границе материалов сохраняется нормальная составляющая B; тангенциальная H непрерывна при отсутствии поверхностного свободного тока. Внешняя граница и симметрия являются физическими условиями задачи. Их изменение проверяется отдельно от уменьшения конечных элементов: сходимость внутри неверной области не устраняет ошибку постановки.

\nabla\times\mathbf H=\mathbf J,\quad \nabla\cdot\mathbf B=0,\quad \mathbf B=\nabla\times(A_z\mathbf e_z)
Статическая/квазистатическая система без существенного тока смещения.
\nabla\cdot(\nu\nabla A_z)=-J_z,\qquad B_x=\partial_y A_z,\quad B_y=-\partial_x A_z
Здесь \nu =1/\mu; производные берутся по координатам в метрах, если A_z задан в Вб/м.

Литература по методу: Уравнения Максвелла и материальные соотношения · PDF3–20; Векторная алгебра и векторный магнитный потенциал · PDF1–12,24–31; Математическая модель поля и электрических цепей машины · PDF1–11

Разделение переменных и эталон гладкого зазора

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

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

A_k(x,y)=(c_ke^{ky}+d_ke^{-ky})e^{jkx}
Плоская гармоника при k\ne 0; реальные поля получают действительной частью.
A_\nu(r,\theta)=(a_\nu r^\nu+b_\nu r^{-\nu})e^{j\nu\theta},\qquad A_0=a_0+b_0\ln r
Кольцевой линейный зазор; коэффициенты следуют из заданных условий на двух радиусах.

Литература по методу: Аналитические методы: назначение эталонов и допущения · PDF1–19; Разделение переменных: поле плоского и кольцевого зазора · PDF1–24

Конформное преобразование и показатель поля у угла

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

В преобразовании Шварца–Кристоффеля показатели множителей определяются внутренними углами многоугольника. Возле угла \alpha локальная координата связана с исходной степенью \alpha /\pi. Отсюда следует показатель поля \pi /\alpha -1. Для входящего угла 270^\circ поле растёт как r^(-1/3); измерение показателя по нескольким расстояниям позволяет объяснить рост сеточного максимума в L-образном проводнике и необходимость конечного радиуса скругления.

\frac{dF}{dz}=\frac{dF/dw}{dz/dw},\qquad \frac{dz}{dw}=C\prod_k(w-a_k)^{\alpha_k/\pi-1}
Первая формула задаёт масштаб комплексного поля, вторая — производную отображения Шварца–Кристоффеля.
z-z_0\sim(w-w_0)^{\alpha/\pi},\qquad |\nabla U|\sim r^{\pi/\alpha-1}
Локальный асимптотический закон; коэффициент и конечный диапазон расстояний задаются полной краевой задачей.

Литература по методу: Конформные преобразования и комплексная потенциальная функция · PDF1–17,26–30; Преобразование Шварца–Кристоффеля · PDF1–13

Чётное и нечётное поле: от зубцовой зоны к Картеру

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

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

\mathbf B=\mathbf B_{\rm even}+\mathbf B_{\rm odd},\qquad k_C=\delta_{\rm eff}/\delta
Разложение использует одну линейную геометрию и согласованные возбуждения.
\delta_{\rm eff}=\frac{2\mu_0 W}{B_{\rm av}^2\tau l_z}
Энергетическое определение для выбранного зубцового деления при фиксированном среднем потоке в воздушной области.

Литература по методу: Зубчатость зазора, Картер и зубцовые контуры · PDF1–20,34–54; Hrabovcová, Rafajdus, Makyš. Analysis of Electrical Machines (2020) · PDF21–31

Модели для работы

Пазовый зазор: коэффициент Картера и комплексная передача гармоник

Плоская задача, активная длина 1 м, зубцовое деление 40 мм, зазор 2 мм, глубина паза 20 мм. Раскрытие 0, 4, 8 или 12 мм. Воздух имеет \mu _{r}=1; поверхности идеальных сердечников задают H\tau =0. В серии Картера задан средний поток Bср=0,8 Тл через постоянные значения A на боковых сечениях. В гармонической серии одинаковый токовый слой толщиной 0,3 мм расположен у верхнего сердечника; J_{1}=10 А/мм², J_{3}/J_{1}=0{,}3∠0,6 рад. Проводимость равна нулю, поэтому частота 50 Гц задаёт фазоры без вытеснения тока. Идеальные линейные сердечники исключают насыщение и потери в стали. Плоская периодическая область не содержит кривизну, скос и торцы. Комплексное отношение при совместном источнике не равно универсальной передаче отдельной гармоники.

Модели и сохранённые решения серии
Результаты расчёта
ВеличинаРезультатЕдиница
Коэффициент Картера, b=8 мм1,0991
Эффективный зазор2,198мм

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

1 · Постановка и нормировка потока

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

  1. Запишите геометрический зазор, раскрытие, деление и активную длину. Покажите поверхности идеального железа и два боковых сечения с заданным A.
  2. Для опыта Картера вычислите разность A из Bсрτ; для гармонического опыта объявите нормировку J, пространственные фазы и координату выборки.
  3. До расчёта назначьте допустимое увеличение МДС и список гармоник. Зафиксируйте, что в каждой из двух серий остаётся постоянным.

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

Подробное выполнение: 8.5 · Источники тока, витки, знаки и связи между токами и напряжениями на выводах; 3.2 · Условия Дирихле, Неймана, Робина и открытая граница

2 · Эталон и коэффициент Картера

Гладкая задача проверяет единицы и способ получения энергии.

  1. Рассчитайте гладкий зазор на сетках 0,60; 0,30; 0,15 мм и сопоставьте k\delta с единицей.
  2. При раскрытии 8 мм повторите три сетки; для 4 и 12 мм получите сравнительные точки.
  3. Вычислите эффективный зазор из энергии; сопоставьте с конформной формулой.

Записать в протокол. Таблица b_{0}, h, W, δэфф, k\delta и отклонение от аналитической оценки.

Подробное выполнение: 3.2 · Условия Дирихле, Неймана, Робина и открытая граница; 3.3 · Независимое уточнение сетки, области и положения; 37.3 · Зубцово-пазовая структура статора и контроль тиражирования

3 · Передача гармоник

Совпадающие источники отделяют влияние границы от влияния обмотки.

  1. Рассчитайте гладкую и пазовую модели при 24 и 48 участках одинакового токового слоя.
  2. Возьмите 256 точек на y=1 мм, без повторения конечной точки периода, и выполните комплексное преобразование Фурье.
  3. Для порядков 1 и 3 вычислите \operatorname{Re} C\nu, \operatorname{Im} C\nu, |C\nu | и arg C\nu; объясните смешение гармоник.

Записать в протокол. Профили комплексного By(x), спектр и отношения одинаковых гармоник.

Подробное выполнение: 37.8 · Угловое потокосцепление, ДПФ и ЭДС холостого хода; 4.1 · Наблюдаемые величины и обработка результатов; 8.6 · Периодичность, антипериодичность и скользящий воздушный зазор

4 · Погрешность и выбор раскрытия

Решение принимается по приросту МДС и спектру при объявленном источнике.

  1. Сведите сеточную разность k\delta, ошибку относительно конформной оценки и изменение B_{1} при дроблении слоя в отдельную таблицу.
  2. Выберите раскрытие по заданному пределу МДС и проверьте комплексные C_{1}, C_{3} для той же геометрии.
  3. Укажите, какие дополнительные данные о стали, обмотке и кривизне нужны для переноса результата на конкретную машину.

Записать в протокол. Выбранное раскрытие, запас относительно численной ошибки и границы применения.

Подробное выполнение: 3.3 · Независимое уточнение сетки, области и положения; 37.8 · Угловое потокосцепление, ДПФ и ЭДС холостого хода

Контрольный расчёт и численные результаты

\tau =40 мм, \delta =2 мм, b_{0}=8 мм, глубина паза 20 мм, Bср=0,8 Тл, l=1 м.

  1. По энергии получено k\delta =1{,}09915 и δэфф=2,19829 мм.
  2. При сохранённом потоке потребная МДС возрастает на 9,91466%.
  3. Для одинакового совместного токового слоя комплексные отношения C_{1} и C_{3} приведены в таблице. Они сохраняют амплитуду и фазу.

Результаты собственных полевых моделей FEMM 5 с линейными идеальными сердечниками.

Раскрытие паза и эффективный зазор

b_{0}, ммСетка, ммk\delta FEMMk\delta конформныйδэфф, мм
00,6112
00,3112
00,15112
40,31,0291,0292,059
80,61,11,0992,201
80,31,11,0992,199
80,151,0991,0992,198
120,31,1991,1982,398

Комплексное отношение при общем источнике

Порядок \nu\operatorname{Re} C\nu\operatorname{Im} C\nu|C\nu |arg C\nu, рад
10,9869-0,0010780,9869-0,001092
30,72740,14740,74220,1999

Количественные результаты и их проверка

Заполните таблицу собственными значениями; каждому числу сопоставьте режим и файл решения.

ВеличинаЕдиницаЧто определяет
k\delta и δэфф1 ммДобавочная МДС раскрытия
C\nu1 радИзменение пространственного спектра
  1. Гладкий зазор: k\delta отличается от единицы менее 0,3%.
  2. Для b_{0}=8 мм изменение k\delta при сетке 0,30→0,15 мм менее 1,5%; конформная оценка отличается менее 1%.
  3. Гармоника гладкого зазора отличается от решения Лапласа менее 2,5%; дробление слоя 24→48 меняет B_{1} менее 2,5%.

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

Границы применения результата

Литература

  1. Гладкий зазор: аналитическая проверка поля и независимых источников ошибки. — 12 с.
  2. Чётное и нечётное поля паза: Картер и разностная проводимость. — 3 с.
  3. Нечётное поле паза как аналитический эталон рассеяния. — 8 с.
  4. Чётное поле паза и контроль комплексной потенциальной функции. — 8 с.
  5. Влияние раскрытия пазов на гармоники воздушного зазора. — 7 с.
  6. Уравнения Максвелла. Раздел 2. — 20 с.
  7. Основные понятия векторной алгебры, используемые при описании электромагнитного поля. Векторный магнитный потенциал. Раздел 3. — 31 с.
  8. Математическое описание электромагнитного поля в электрической машине. Раздел 5. — 15 с.
  9. Аналитические методы исследования поля в электрической машине. Общие положения. Раздел 7. — 19 с.
  10. Аналитические методы исследования поля в электрической машине. Метод разделения переменных. Поле в плоском и кольцевом зазоре. Поверхностный ток как источник поверхностного потенциала. Раздел 8. — 33 с.
  11. Конформные преобразования. Комплексная потенциальная функция. Раздел 9. — 30 с.
  12. Конформные преобразования (продолжение). Преобразование Шварца — Кристоффеля. Раздел 10. — 13 с.
  13. Магнитное поле в зазоре электрической машины с учётом зубчатости сердечников. Раздел 11. — 54 с.
  14. Hrabovcová V., Rafajdus P., Makyš P. Analysis of Electrical Machines. — London : IntechOpen, 2020. — ISBN 978-1-83880-207-3. — DOI: 10.5772/intechopen.90135. Главы 1–3 распространяются по лицензии CC BY 3.0 [13].