FEMM 5

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

Зубцовые контуры: проводимости, коронки и взаимная связь при смещении

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

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

Задачи

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

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

Матрица двух конечных зубцов рассчитывается при \mu _{r}=10000, зазорах 1, 2, 4 мм и единичных ампер-витках. Противоположные токи сторон одной катушки имеют согласованные знаки.

Серия встречных зубцов использует антипериодическую ячейку с шагом зубцов 30 мм: соседний зубцовый контур имеет противоположную МДС. Магнитный период — 60 мм. Смещение 0…60 мм с шагом 3 мм исследуется двумя независимыми возбуждениями.

Для коронок ширина шлица 20 мм, зазор 5 или 20 мм, глубина 10, 20, 30, 50 мм. Ток 1000 А расположен в нижнем слое толщиной 0,5 мм; опорная линия находится на 1 мм выше дна.

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

От распределения A к матрице индуктивностей

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

В линейной системе потокосцепление является линейной функцией токов. Независимые возбуждения дают столбцы матрицы L. Симметрия энергетического оператора требует L_{12}=L_{21}; знак взаимного элемента зависит от направления сторон обмоток.

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

\Psi_k=N_kl(\overline A_{k+}-\overline A_{k-}),\quad \boldsymbol\Psi=\mathbf L\mathbf i,\quad W=\tfrac12\mathbf i^T\mathbf L\mathbf i
Усреднение выполняется по соответствующим площадям сторон, а не по произвольной линии внутри паза.
L_{12}=L_{21},\qquad |M|\le\sqrt{L_{11}L_{22}}
Взаимность и неотрицательная энергия проверяют согласованность цепей и знаков. Для насыщенной системы применяют локальную дифференциальную матрицу.

Литература по методу: Потокосцепление и ЭДС обмотки: сравнение независимых способов · PDF 1–6: геометрия и прямые методы; 7–12: спектр; 12–13: вращающийся вектор; 14–16: программа и варианты; Зубцовый контур и взаимная проводимость при движении · PDF 1–2: определение; 3–4: программа и движение; 5–6: сравнение токовых моделей

Коронка — отдельная поправка к определённому потоку

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

Разность A этих линий даёт поток выбранной области. Из его проводимости вычитают вклад прямоугольной равномерной части h'/b. Полученная поправка является знаковой: небольшое отрицательное значение означает уменьшение относительно выбранной прямоугольной опоры, а не отрицательную полную магнитную проводимость.

Для глубокого шлица и \beta =1 две приведённые формулы алгебраически эквивалентны. При малой глубине поле опорной линии может оставаться неравномерным; это проверяется серией по глубине и сетке.

\lambda_{\rm кор}=\frac{A(0,-h\prime)-A(0,\delta)}{\mu_0I}-\frac{h\prime}{b}
Активная длина сокращается. Опорная линия должна лежать выше токового слоя в достаточно равномерной части поля.
\lambda_V=\frac{\ln[(\delta/b)^2+1/4]+4(\delta/b)\arctan[b/(2\delta)]}{2\pi}
Формула проводимости коронок напротив гладкого сердечника.
\lambda_{IS}=\frac1{2\rho}-\frac1\pi\left[\ln\frac2{\sqrt{1+\rho^{-2}}}+\frac{\arctan(\rho^{-1})}{\rho}\right],\quad\rho=\frac b{2\delta}
Вариант формулы при полном шаге обмотки \beta =1; сравнивается с тем же определением поправки.

Литература по методу: Проводимость рассеяния по коронкам зубцов · PDF 1: определение; 2–3: формулы и численные примеры; 4: геометрическая чувствительность; 5: программа

От потенциала к потоку и реальной обмотке

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

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

\Phi=l_z\Delta A_z,\qquad \Psi_k=l_z\sum_jN_{kj}\langle A_z\rangle_{S_j}
Знаки разностей и N согласованы с нормалями и направлениями обмотки; усреднение соответствует распределённой стороне катушки.
\boldsymbol\Psi=\mathbf L\mathbf i,\quad L_{kj}\simeq L_{jk},\quad L_{\rm sec}=\Psi/i,\quad L_{\rm diff}=\partial\Psi/\partial i
Матрица постоянна только для линейного материала при фиксированной геометрии.

Литература по методу: Векторная алгебра и векторный магнитный потенциал · PDF24–31; Математическая модель поля и электрических цепей машины · PDF9–11; Магнитные поля в электрических машинах — введение (файл Гончарова) · PDF1–20

Энергия, коэнергия и выбор независимой переменной

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

Равенство W=W'=i\Psi /2 требует линейности и отсутствия независимого возбуждения в рассматриваемой однопортовой задаче. Для нескольких обмоток энергетическое тождество содержит сумму токов на потокосцепления всех независимых портов. При постоянном возбуждении генератора динамический баланс использует отдельно энергию переменных токов по формуле блока d/q и мощность преобразования. Для насыщенной стали интегрируют фактическую характеристику. Сила определяется частной производной при фиксированных токах, а не производной вдоль нагруженной траектории с меняющимся током.

W=\int_0^\Psi i(\psi)\,d\psi,\quad W^\prime=\int_0^i\Psi(j)\,dj,\quad W+W^\prime=i\Psi
Однозначная квазистатическая характеристика без отдельного гистерезисного цикла.
F_x=\left.\frac{\partial W^\prime}{\partial x}\right|_{\mathbf i},\qquad T=\left.\frac{\partial W^\prime}{\partial\theta_m}\right|_{\mathbf i}
Угол \theta m выражается в радианах; положительное направление координаты задано в модели.

Литература по методу: Энергия магнитного поля через векторный потенциал · PDF1–9,19–30; Hrabovcová, Rafajdus, Makyš. Analysis of Electrical Machines (2020) · PDF130–146,184–189

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

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

Эффективный зазор определяется по интегральному потоку или энергии при явно фиксированной МДС либо фиксированном потоке. Коэффициент Картера связывает реальную зубцовую область с гладкой моделью. Он не служит множителем для любого локального пика 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

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

Зубцовые контуры: проводимости, коронки и взаимная связь при смещении

Матрица двух конечных зубцов рассчитывается при \mu _{r}=10000, зазорах 1, 2, 4 мм и единичных ампер-витках. Противоположные токи сторон одной катушки имеют согласованные знаки. Серия встречных зубцов использует антипериодическую ячейку с шагом зубцов 30 мм: соседний зубцовый контур имеет противоположную МДС. Магнитный период — 60 мм. Смещение 0…60 мм с шагом 3 мм исследуется двумя независимыми возбуждениями. Для коронок ширина шлица 20 мм, зазор 5 или 20 мм, глубина 10, 20, 30, 50 мм. Ток 1000 А расположен в нижнем слое толщиной 0,5 мм; опорная линия находится на 1 мм выше дна. Матрица относится к линейной стали \mu _{r}=10000 и объявленным направлениям контуров. Период 60 мм описывает чередование МДС, а не равномерное униполярное возбуждение. Выделенная по двум линиям A проводимость рассеяния отличается от полной проводимости катушки по энергии.

Модели и сохранённые решения серии
Результаты расчёта
ВеличинаРезультатЕдиница
Взаимная индуктивность при совмещении7,675мкГн
Поправка коронок b=20, \delta =5 мм-0,00891

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

1 · Контуры, знаки и выделяемая проводимость

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

  1. Сопоставьте стороны каждой катушки со знаками тока и положительным направлением потокосцепления.
  2. Объявите границы участка рассеяния: опорную линию шлица и касающуюся линию у гладкого сердечника. Не подменяйте этот поток полной энергией.
  3. Для движущихся зубцов запишите шаг 30 мм, антипериодическую МДС соседнего контура и полный магнитный период 60 мм.

Записать в протокол. Схема контуров, направления токов и определение каждой проводимости.

Подробное выполнение: 8.6 · Периодичность, антипериодичность и скользящий воздушный зазор; 8.7 · Одна катушка и потокосцепление без смешения с потоком

2 · Два независимых возбуждения

Один расчёт не определяет матрицу двух катушек.

  1. Возбудите только катушку A, затем только B; извлеките обе пары потокосцеплений.
  2. Вычислите среднее A по каждой стороне и независимо восстановите потокосцепление.
  3. Сопоставьте взаимные элементы и энергию при зазорах 1, 2, 4 мм.

Записать в протокол. L_{11}, L_{22}, M, энергия и невязки.

Подробное выполнение: 8.7 · Одна катушка и потокосцепление без смешения с потоком; 8.8 · Две обмотки и матрица собственной и взаимной индуктивности; 8.9 · Знаки M, взаимность и энергетическая проверка

3 · Зубцы в относительном движении

Положение определяет связь контуров.

  1. Выполните два возбуждения в каждом положении 0…60 мм с шагом 3 мм.
  2. Постройте M(x), сравните M_{12} и M_{21} в каждой точке.
  3. Проверьте смену знака после 30 мм и возврат после 60 мм; укажите физический смысл антипериодической МДС.

Записать в протокол. Матрицы и график взаимной проводимости на магнитный период.

Подробное выполнение: 8.6 · Периодичность, антипериодичность и скользящий воздушный зазор; 17.2 · Таблицы потокосцеплений и локальная матрица индуктивностей; 37.3 · Зубцово-пазовая структура статора и контроль тиражирования

4 · Рассеяние коронок

Границы потока определяют смысл вычисленной проводимости.

  1. Найдите A на опорной линии шлица и на гладкой поверхности напротив его оси.
  2. Получите λкор для двух зазоров и четырёх глубин; сопоставьте две аналитические формулы.
  3. Для глубины 20 мм и зазора 5 мм повторите сетки 0,60; 0,30; 0,15 мм и оцените абсолютную погрешность малой поправки.

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

Подробное выполнение: 4.1 · Наблюдаемые величины и обработка результатов; 3.3 · Независимое уточнение сетки, области и положения; 23.3 · Роторные токи, вытеснение тока и параметры схемы

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

Раскрытие 20 мм; зазоры 5 и 20 мм. Противоположные зубцовые контуры: шаг 30 мм, магнитный период 60 мм, \mu _{r}=10000.

  1. При b=20 мм, \delta =5 мм и h'=19 мм получена λкор=-0,00889997; формула даёт -0,00891301.
  2. При \delta =20 мм, h'=29 мм получена λкор=0,330032; формула даёт 0,330682.
  3. Взаимная индуктивность противоположных контуров M(0)=7{,}67455 мкГн; M(30 мм) меняет знак, M(60 мм) восстанавливается.

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

Проводимость коронок при разных глубинах

\delta, ммh', ммСетка, ммλкор FEMMλкор формула
590,380,18 × 10⁻⁶-0,008913
5190,6-0,009133-0,008913
5190,3-0,0089-0,008913
5190,15-0,008658-0,008913
5290,3-0,009283-0,008913
5490,3-0,009317-0,008913
2090,30,33460,3307
20190,30,33020,3307
20290,30,330,3307
20490,30,330,3307

Взаимная проводимость и смещение

x, ммL_{11}, мкГнM_{12}, мкГнM_{21}, мкГн
014,077,6757,675
313,366,8436,843
612,025,1775,177
910,583,223,22
129,4871,3891,389
159,171122,6 × 10⁻⁶122,6 × 10⁻⁶
189,486-1,388-1,388
2110,58-3,22-3,22
2412,02-5,177-5,177
2713,37-6,844-6,844
3014,07-7,675-7,675
3313,36-6,843-6,843
3612,02-5,177-5,177
3910,58-3,22-3,22
429,487-1,389-1,389
459,171-122,6 × 10⁻⁶-122,6 × 10⁻⁶
489,4861,3881,388
5110,583,223,22
5412,025,1775,177
5713,376,8446,844
6014,077,6757,675

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

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

ВеличинаЕдиницаЧто определяет
L_{11}, L_{22}, M(x)ГнЭнергия и взаимная связь контуров
λкор1Поправка к пазовому рассеянию
  1. Среднее A и потокосцепление цепи согласуются менее чем на 0,2%; энергия согласуется с LI ^{2}/2.
  2. Взаимность и повторение через 60 мм проверяются с допуском 0,2%.
  3. Абсолютная сеточная разность λкор оценивается отдельно от относительной: поправка может быть близка к нулю.

Итоговый вывод. Выберите модель проводимости коронок по отношению раскрытия к зазору и глубине равномерного участка. Объясните знак поправки и влияние положения зубцов на M.

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

Литература

  1. Потокосцепление и ЭДС обмотки: сравнение независимых способов. — 16 с.
  2. Зубцовый контур и взаимная проводимость при движении. — 6 с.
  3. Проводимость рассеяния по коронкам зубцов. — 5 с.
  4. Основные понятия векторной алгебры, используемые при описании электромагнитного поля. Векторный магнитный потенциал. Раздел 3. — 31 с.
  5. Математическое описание электромагнитного поля в электрической машине. Раздел 5. — 15 с.
  6. Магнитные поля в электрических машинах — введение. — 75 с.
  7. Выражение энергии магнитного поля через векторный магнитный потенциал. Раздел 4. — 30 с.
  8. 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].
  9. Магнитное поле в зазоре электрической машины с учётом зубчатости сердечников. Раздел 11. — 54 с.