FEMM 5

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

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

Цель. Определить отношение Rac/Rdc, внутреннюю реактивность, распределение тока и потери. По заданным количественным ограничениям выбрать число изолированных жил и частоту вращающегося поля.

Связь с проектированием. Выбрать минимальное число транспонированных жил с Rac/Rdc≤1,5 среди исследованных вариантов; для проводящего цилиндра выбрать частоту с максимальным моментом при потерях не более 50 Вт/м.

Задачи

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

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

Глубокий стержень: 8×40 мм, \sigma =20{,}5 МСм/м, Ipeak=1280 А, f=0, 50, 400 Гц. Идеальные магнитные стенки создают одномерный эталон вытеснения тока.

Пакет меди: 304 мм², габарит 8×42 мм, две колонки и n=1, 2, 4, 8, 16, 32, 64 ряда. Температура 120 °C; \sigma =58/[1+0{,}00393(120-20)] МСм/м. Общий амплитудный ток 6000 А, частота 400 Гц. Геометрия меди и изоляции одинакова в паре сравниваемых соединений.

Цилиндр радиусом 10 мм имеет \sigma =20{,}5 МСм/м и нулевую механическую скорость. Фазированный слой радиусов 19…20 мм создаёт двухполюсное вращающееся поле. Внешняя граница радиусом 160 мм, активная длина 1 м.

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

Вытеснение тока — задача магнитной диффузии

В проводнике переменное поле создаёт вихревую составляющую электрического поля. Она меняет распределение J и потерь даже при сохранённом полном токе. Характерная длина проникновения уменьшается с частотой, проводимостью и проницаемостью.

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

В FEMM токи — комплексные амплитуды. Поэтому средняя мощность AC равна половине R|Ipeak|². При постоянном токе коэффициента 1/2 нет. Сопротивление по мощности нужно согласовать с действительной частью напряжения цепи.

\delta_s=\sqrt{\frac1{\pi f\mu\sigma}},\quad \xi=h/\delta_s
\delta s — глубина проникновения в однородное полупространство; h — глубина проводника.
k_R=\xi\frac{\sinh2\xi+\sin2\xi}{\cosh2\xi-\cos2\xi},\quad k_X=\frac3{2\xi}\frac{\sinh2\xi-\sin2\xi}{\cosh2\xi-\cos2\xi}
kX относится только к внутренней реактивности стержня.
R_{ac}=2P/|I_{peak}|^2,\quad R_{dc}=l/(\sigma S),\quad L_{int,ac}=4W_{av,int}/|I_{peak}|^2
Для DC энергия равна LI ^{2}/2; средняя магнитная энергия гармонического поля — L|Ipeak|^{2}/4.

Литература по методу: Расчёты вихревых полей в электромеханике · PDF 1–8 полностью: 1–3 транспозиция и добавочные потери; 3–6 клетка и вытеснение; 6–8 перемешиватель и источник расчётов.; Транспонированный стержень: выбор числа элементарных проводов · PDF 1: физика; 2–3: цепь и расчётные соотношения; 4–5: измерения и варианты

Дробление и транспозиция задают разные условия токов

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

Идеальная полная транспозиция означает, что каждая элементарная ветвь за полный виток занимает одинаковые положения. В двумерной активной части её предельный эффект задаётся равными токами всех 2n жил. При этом Ie=Iобщ/(2n), а не Iобщ.

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

\widehat I_e=\frac{\widehat I_\Sigma}{2n},\qquad \sum_e\widehat I_e=\widehat I_\Sigma
В режиме равных токов Series-цепь задаёт один Ie в каждой жиле; в Parallel-цепи задаётся общий ток.
Z_{eq}=\frac{V_{parallel}}{I_\Sigma}\quad\text{или}\quad Z_{eq}=\frac{V_{series}}{(2n)I_\Sigma}
Напряжение Series-цепи суммируется по жилам; для общего эквивалента используют среднее напряжение одной физической параллельной ветви.

Литература по методу: Транспонированный стержень: выбор числа элементарных проводов · PDF 1: физика; 2–3: цепь и расчётные соотношения; 4–5: измерения и варианты; Расчёты вихревых полей в электромеханике · PDF 1–8 полностью: 1–3 транспозиция и добавочные потери; 3–6 клетка и вытеснение; 6–8 перемешиватель и источник расчётов.

Частота вращающегося поля и электромагнитное перемешивание

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

Для двухполюсного поля и неподвижного тела переданная электромагнитная мощность превращается в джоулевы потери, P=ωсинх T. Сопоставление интеграла Лоренца, взвешенного тензора напряжений и мощности даёт независимый контроль.

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

\widehat J_s(\theta)=J_m e^{-j\theta},\quad \omega_{syn}=2\pi f,\quad P_{rotor}=\omega_{syn}|T|\ \ (n=0,p=1)
Фазированный слой из 48 участков приближает бегущую пространственную волну первого порядка.

Литература по методу: Расчёты вихревых полей в электромеханике · PDF 1–8 полностью: 1–3 транспозиция и добавочные потери; 3–6 клетка и вытеснение; 6–8 перемешиватель и источник расчётов.

Диффузия поля и частота в системе проводника

Из законов Фарадея, полного тока и Ома для неподвижной линейной проводящей среды следует уравнение магнитной диффузии. Его гармоническое решение задаёт глубину проникновения. Частота должна соответствовать системе координат рассматриваемого проводника: для ротора основной асинхронной гармоники это частота скольжения.

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

\frac{\partial\mathbf B}{\partial t}=\frac1{\mu\sigma}\nabla^2\mathbf B,\qquad \delta_s=\sqrt{\frac2{\omega\mu\sigma}},\qquad k=\frac{1+j}{\delta_s}
Изотропные постоянные \mu и \sigma, ток смещения несущественен.
f_2=sf_1,\qquad P_{\rm Cu,2}=sP_\delta,\qquad T_e=P_\delta/\Omega_s
Последние равенства относятся к согласованной основной гармонике и принятой схеме асинхронной машины.

Литература по методу: Электродинамические поля и проникновение в проводящий массив · PDF1–23,34–42; Интегральные величины при резко выраженном поверхностном эффекте · PDF1–15; Магнитные поля в электрических машинах — введение (файл Гончарова) · PDF60–74

Скин-эффект, поле рассеяния и дробление проводника

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

Число элементарных проводов выбирают по потерям при заданных температуре, медном сечении и полном токе. Идеальная транспозиция задаёт равные токи; параллельное соединение без транспозиции задаёт равенство напряжений и допускает неравные токи. Эти две электрические модели сравниваются отдельно. Одна только геометрическая нарезка массива на полосы не определяет схему соединения.

\varphi(\xi)=\xi\frac{\sinh2\xi+\sin2\xi}{\cosh2\xi-\cos2\xi},\qquad \psi(\xi)=2\xi\frac{\sinh\xi-\sin\xi}{\cosh\xi+\cos\xi}
Нормированная высота \xi согласуется с геометрией и глубиной проникновения выбранного пазового эталона.
\varphi(\xi)=1+\frac4{45}\xi^4+O(\xi^8),\qquad \psi(\xi)=\frac13\xi^4+O(\xi^8)
Независимые низкочастотные проверки: собственный фактор стремится к 1, добавочный — к 0.
R_{ac}=P_{\rm Cu}/I_{\rm rms}^2,\qquad \widehat I_{\rm strand}=\widehat I_{\rm total}/N
Вторая формула применена только при принятой идеальной транспозиции и одинаковых элементарных проводах.

Литература по методу: Переменное поле в прямоугольном проводнике паза · PDF1–21; Электродинамические поля и проникновение в проводящий массив · PDF24–42; Расчёты вихревых полей в электромеханике · PDF1–6

Поверхностный импеданс как проверка сильного вытеснения

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

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

Z_s=\frac{1+j}{\sigma\delta_s},\qquad P\simeq\frac12\operatorname{Re}Z_s\int_S|\widehat H_t|^2dS
Фазоры пиковые; для RMS множитель 1/2 не используется. Тонкие тела и зоны краёв требуют отдельной проверки.

Литература по методу: Интегральные величины при резко выраженном поверхностном эффекте · PDF1–15; Электродинамические поля и проникновение в проводящий массив · PDF9–23

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

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

Глубокий стержень: 8×40 мм, \sigma =20{,}5 МСм/м, Ipeak=1280 А, f=0, 50, 400 Гц. Идеальные магнитные стенки создают одномерный эталон вытеснения тока. Пакет меди: 304 мм², габарит 8×42 мм, две колонки и n=1, 2, 4, 8, 16, 32, 64 ряда. Температура 120 °C; \sigma =58/[1+0{,}00393(120-20)] МСм/м. Общий амплитудный ток 6000 А, частота 400 Гц. Геометрия меди и изоляции одинакова в паре сравниваемых соединений. Цилиндр радиусом 10 мм имеет \sigma =20{,}5 МСм/м и нулевую механическую скорость. Фазированный слой радиусов 19…20 мм создаёт двухполюсное вращающееся поле. Внешняя граница радиусом 160 мм, активная длина 1 м. Равные токи описывают идеальную полную транспозицию; реальные лобовые соединения, длина и контактные сопротивления отдельно не рассчитаны. Температура меди задана, а не получена из теплового равновесия. Для цилиндра задана нулевая скорость. Результат не определяет скорость расплава и не включает гидродинамику.

Модели и сохранённые решения серии
Результаты расчёта
ВеличинаРезультатЕдиница
kR глубокого стержня при 50 Гц2,5271
Выбранное число жил1281
Rac/Rdc выбранного пакета1,2171
Частота цилиндра при P\le 50 Вт/м200Гц

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

1 · Электрические условия и бюджет потерь

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

  1. Выпишите температуру, проводимость, площадь меди, амплитудный и действующий общий ток. Проверьте формулу Ie=Iобщ/(2n).
  2. Составьте две отдельные схемы: общее напряжение параллельных ветвей и равные токи идеальной транспозиции. Укажите, как напряжение FEMM преобразуется в эквивалентное Z.
  3. Назначьте ограничения Rac/Rdc≤1,5 для пакета и P\le 50 Вт/м для цилиндра. Скорость цилиндра зафиксируйте равной нулю.

Записать в протокол. Две схемы соединений, нормировка фазоров и критерии двух инженерных решений.

Подробное выполнение: 35.1 · Потери в проводниках обмоток: скин-эффект, близость и гармоническая сетка; 34.1 · Магнитные потери, потери в проводниках и граница баланса мощности

2 · Одномерный эталон стержня

Аналитическая проверка выявляет ошибки фазоров и мощности.

  1. Рассчитайте DC и AC стержень; получите профили комплексной J по глубине.
  2. По P найдите kR; по магнитной энергии только проводника найдите kX.
  3. Сопоставьте с формулами Эмде и повторите сетки 0,70; 0,35; 0,175 мм при 50 Гц.

Записать в протокол. \delta s, kR, kX, Z, P и профили J.

Подробное выполнение: 23.3 · Роторные токи, вытеснение тока и параметры схемы; 29.1 · Серия расчётов по частоте и определение параметров асинхронного двигателя; 35.1 · Потери в проводниках обмоток: скин-эффект, близость и гармоническая сетка

3 · Пакет изолированных жил

Одинаковая площадь позволяет сравнивать электрические условия.

  1. Для каждого n рассчитайте Parallel и равные токи Series=6000/(2n) А.
  2. Сохраните комплексный ток каждой жилы, суммарный ток, мощность и полное эквивалентное Z.
  3. Постройте Rac/Rdc от числа жил и выберите первый исследованный вариант с коэффициентом не более 1,5.

Записать в протокол. Таблица 2n, Ie, Rac/Rdc, P и распределение токов ветвей.

Подробное выполнение: 35.1 · Потери в проводниках обмоток: скин-эффект, близость и гармоническая сетка; 34.1 · Магнитные потери, потери в проводниках и граница баланса мощности; 29.2 · Проверка параметров на независимых частотах

4 · Выбор частоты силового воздействия

Силовой максимум нужно сопоставить с тепловым ограничением.

  1. Выполните частотную серию 10…800 Гц для неподвижного цилиндра.
  2. Сравните моменты Лоренца и WST, а также P и 2πfT.
  3. Найдите максимум момента среди исследованных частот; затем выберите частоту при ограничении P\le 50 Вт/м.
  4. Отдельно уточните сетку, число участков источника и внешнюю область при 400 Гц.

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

Подробное выполнение: 50.2 · Вихревые токи в электромагнитном аппарате: экранирование и предел гармонической модели; 34.1 · Магнитные потери, потери в проводниках и граница баланса мощности; 40.2 · Сила по тензору напряжений Максвелла, взвешенному интегралу и коэнергии

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

6000 А peak, 400 Гц, 120 °C, 304 мм² меди. Условия: Rac/Rdc≤1,5; для цилиндра отдельно P\le 50 Вт/м.

  1. Глубокий стержень при 50 Гц: kR=2,52682; Эмде 2,52659.
  2. Первый исследованный пакет, удовлетворяющий ограничению, содержит 128 жил: Ie=46,875 А peak, Rac/Rdc=1,21709, P=1730{,}79 Вт/м.
  3. Среди исследованных частот максимум момента цилиндра получен при 400 Гц. При ограничении P\le 50 Вт/м выбраны 200 Гц: T=0{,}0389618 Н·м/м, P=48{,}9739 Вт/м.

Выбор относится к исследованной сетке числа жил и частот. Теплоотвод и скорость жидкости отдельно не рассчитаны.

Эталон глубокого стержня

f, ГцСетка, ммkR FEMMkR ЭмдеkX FEMMkX Эмде
00,351111
500,72,5272,5270,60310,5989
500,352,5272,5270,6030,5989
500,1752,5272,5270,60270,5989
4000,357,2127,1970,21640,2084

Одинаковая площадь меди: разные ограничения ветвей

Жил 2nСоединениеIцепи peak, АRac/RdcP, Вт/м
2Общее напряжение60009,40313370
2Равные токи30009,40313370
4Общее напряжение60009,413370
4Равные токи150014,2720300
8Общее напряжение60009,3813340
8Равные токи75025,2935970
16Общее напряжение60009,45913450
16Равные токи37513,6519400
32Общее напряжение60009,57713620
32Равные токи187,54,4046262
64Общее напряжение60009,61313670
64Равные токи93,751,8582642
128Общее напряжение60009,62213683
128Равные токи46,881,2171731

Цилиндр в двухполюсном вращающемся поле

f, ГцT Лоренц, Н·м/мT WST, Н·м/мP, Вт/мНевязка P-\omega T, %
100,002518440,002519070,158280,0260451
500,01237180,01237493,887720,026089
1000,02346390,023469914,74670,026226
2000,03896180,038972248,97390,0267642
2250,04137540,041386658,50880,0269511
2500,04327990,043291868,00240,0271578
4000,04723780,0472521118,7560,0287716
8000,03998780,0400032201,0710,0350827

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

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

ВеличинаЕдиницаЧто определяет
Rac/Rdc, kX, Z1 Ом/мПараметры проводника для схемы замещения
P, IеВт/м, АНагрев и распределение токов
T(f)Н·м/мВыбор частоты электромагнитного воздействия
  1. Стержень: изменение kR по сетке менее 1%; отклонение от Эмде менее 3%, внутренней реактивности менее 2%.
  2. Суммарный ток и \operatorname{Re} Z согласуются с заданным током и мощностью менее чем на 0,2%.
  3. Момент Лоренца, WST и P/ωсинх согласуются менее чем на 2,5%; уточнение сетки и источника менее 1,5%, внешней области менее 4%.

Итоговый вывод. Выберите число жил по Rac/Rdc и частоту по моменту и потерям. Объясните, почему дробление без транспозиции даёт другой результат и почему максимум частоты не означает максимум полезной силы.

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

Литература

  1. Транспонированный стержень: выбор числа элементарных проводов. — 5 с.
  2. Фисенко В. Г. Расчёты вихревых полей в электромеханике. — 8 с.
  3. Электродинамические поля в электрических машинах. Раздел 13. — 42 с.
  4. Расчёт интегральных величин (тока, потока, электрического и магнитного сопротивлений) для тел произвольной формы при резко выраженном поверхностном эффекте. Раздел 15. — 15 с.
  5. Магнитные поля в электрических машинах — введение. — 75 с.
  6. Переменное электромагнитное поле в прямоугольном проводнике в пазу электрической машины. Раздел 14. — 21 с.