FEMM 5

Часть I · Интерфейс FEMM 5, геометрия, сетка и математические основы → § 6 · Первый полный расчёт

Занятие 6.1

Полная последовательность расчёта воздушной катушки в графическом интерфейсе

Цель и задачи

Цель. Связать смену тока катушки с изменением рассчитанного поля.

Задачи

  • Построить сетку, решить задачу и открыть результат.
  • Снять точечную индукцию и потокосцепление обмотки.
Что нужно вспомнить перед началом

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

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

Расчётная лабораторная работа · Согласование возбуждения и преобразование тока в магнитный отклик

Инженерный вопрос. Что должно измениться после смены тока и нового решения?

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

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

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

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

  1. Обосновать постановку и способ измерения. Различите ток цепи, число витков, ампер-витки и потокосцепление. Знаки двух сторон одной обмотки должны соответствовать одному физическому току.

  2. Индукция в центре. Выведите, какие величины сохраняются при одинаковых ампер-витках и какие дополнительно зависят от числа витков.

  3. Полный расчёт и смена тока. Выполните действия практикума по его исходной таблице. Запишите ток и поле в одной контрольной точке для исходного и изменённого режимов.

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

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

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

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

Оформление наблюдений. Ток, витки и их знак, поле в одной точке, потокосцепление и указание линейности материала.

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

Проверка из практикума. Проверьте отношение B_z(2 \,\text{А})/B_z(1 \,\text{А})\approx 2 и B_z(-1 \,\text{А})/B_z(1 \,\text{А})\approx -1; при отклонении сначала проверьте актуальность результата и параметры источника.

Независимый контроль. Сравните согласованные варианты с законом Ампера и законом масштабирования, не задавая один источник независимо дважды.

Критерий выполнения. B_{\text{ц}} практически сохраняется при том же NI; \Psi удваивается, потому что удвоено число витков, сцепляющихся с тем же потоком.

Точность и расхождения. Сравнивайте отношения при одинаковой геометрии и материале; отличие нелинейного режима от линейной пропорции сначала объясняют изменением рабочей точки.

Основание сравнения. Опора — закон Ампера и определения цепи и потокосцепления из занятия. Масштабирование по току линейно только для линейного материала при неизменной геометрии.

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

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

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

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

Введение: поле как основа электромеханического преобразования

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

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

Получаемые величины. B и H с единицами; \Phi и \Psi по фактической обмотке; секущая L=\Psi /i; ЭДС при заданном изменении потокосцепления

  1. Задать устройство, источники и нормали поверхностей.
  2. Предсказать направление поля и силы до расчёта.
  3. Получить локальное поле и интегралы в согласованной геометрии.
  4. Проверить единицы, знаки и закон Фарадея; объяснить проектную пользу каждой величины.

Страницы источника. PDF 1–13 полностью: 1–2 преобразование энергии; 3–8 величины и законы; 9–13 задачи по картинам поля.

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

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

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

Соглашения об обозначениях и единицах
  • \(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. (6.1) · Что показывают линии потенциала
  2. (6.2) · Независимый эталон конечного соленоида
  3. (6.3) · Первый закон линейности
  4. (6.4) · Как получается контрольное поле катушки
  5. (6.5) · Как получается контрольное поле катушки
  6. (6.6) · Численный пример
  7. (6.7) · Численный пример

От геометрии к физическому решению

В этом занятии используется та же осесимметричная катушка, чтобы новый результат не сопровождался ненужной сменой геометрии. Сечение обмотки занимает радиусы 20\text{–}30 \,\text{мм} и высоту 200 \,\text{мм}, а сто витков равномерно распределены по этому сечению. Ток положителен в азимутальном направлении. В модели отсутствуют ферромагнитные детали, поэтому связь между источником и полем линейна. Внешняя полуокружность радиусом 300 \,\text{мм} выбрана для первого расчёта; влияние размера внешней области ещё предстоит проверить.

Зачем сетку сначала посмотреть

Автоматическая сетка пригодна для начального решения, однако она может быть слишком редкой у кромок обмотки и слишком плотной далеко от неё. Важно увидеть треугольники по всему воздуху и внутри катушки, убедиться в отсутствии неразбитых щелей и тонких случайных областей. Для исходного опыта установите в обмотке ограничение размера элемента 2 \,\text{мм}, а в окружающем воздухе 10 \,\text{мм}. Ошибка линейного решателя 10^{-8} характеризует решение алгебраической системы, а не точность пространственной модели.

Что показывают линии потенциала

В осесимметричном магнитном поле силовые линии в меридиональном сечении связаны с функцией потока 2\pi rA_\theta. Карта линий помогает увидеть замыкание поля, но число изображённых линий задаётся настройкой визуализации. Оно не равно числу витков и не является непосредственным измерением потока. Цветовая карта индукции также зависит от выбранной шкалы. Для сравнения двух токов сохраните одинаковую шкалу и одновременно снимите числовое значение в одной и той же точке.

\Phi(r,z)=2\pi rA_\theta(r,z)
(6.1)
Функция магнитного потока \Phi измеряется в Вб; A_\theta — азимутальная компонента потенциала в \text{Вб}/\text{м}; r — радиус в м.

Точка, контур, область и цепь отвечают на разные вопросы

Точечное значение характеризует локальное поле. Линейный интеграл измеряет поток через выбранную поверхность либо циркуляцию напряжённости. Интеграл по области даёт энергию, объём или силу при соответствующей формуле. Результаты цепи описывают обмотку целиком: ток, напряжение и потокосцепление. Значение B в центре соленоида нельзя автоматически подставить вместо среднего B по каждому витку. Поэтому первое проверяемое число — индукция в центре, а потокосцепление записывается отдельно для последующих занятий.

Независимый эталон конечного соленоида

Тонкостенный соленоид с равномерной намоткой имеет известное поле в центре. Формула учитывает конечную длину и потому лучше бесконечного соленоида. В данном случае реальная учебная обмотка имеет радиальную толщину 10 \,\text{мм}; использование среднего радиуса 25 \,\text{мм} даёт близкую, но приближённую оценку. При сравнении следует назвать эту идеализацию. Дополнительное расхождение может быть обусловлено внешней границей и сеткой. Несколько процентов допустимы только как начальная учебная оценка, а не как доказательство завершённой верификации.

B_\text{ц}\simeq\frac{\mu_0NI}{2\sqrt{R_\text{ср}^2+(l/2)^2}}
(6.2)
B_{\text{ц}} — индукция в центре в Тл; R_{\text{ср}} и l — средний радиус и длина в м; \mu_0=4\pi\cdot 10^{-7} \,\text{Гн}/\text{м}.

Первый закон линейности

Повторите решение при 2 \,\text{А} и при -1 \,\text{А}. В линейной задаче положительное удвоение тока удваивает индукцию и потокосцепление, а смена знака меняет направление поля. Энергия, напротив, зависит от квадрата тока и остаётся положительной. Такой опыт проверяет источник, знак и выбранную точку лучше визуального впечатления от карты поля. При каждом изменении сохраните отдельную исходную модель и заново запустите решатель; просмотр прежнего результата не заменяет пересчёта.

\mathbf B(2I)=2\mathbf B(I),\qquad W_\text{м}(2I)=4W_\text{м}(I)
(6.3)
Законы справедливы для линейной магнитостатической модели без постоянных магнитов.

Почему именно конечный соленоид удобен для первого опыта

Его поле можно оценить без численного решателя. Представьте намотку как набор круглых витков, равномерно распределённых вдоль оси. Поле каждого витка в центральной точке выражается законом Био — Савара; интегрирование по продольной координате от -l/2 до l/2 даёт формулу конечного соленоида. Средний радиус заменяет конечную радиальную толщину намотки одним эквивалентным слоем. Поэтому близость результата метода конечных элементов к формуле подтверждает согласованность источника и геометрии, но небольшая разность естественна. В сложной модели двигателя пришлось бы одновременно проверять схему обмотки, свойства магнитопровода и условия симметрии. Сначала проверяется простое поле, затем на нём изучаются остальные инструменты.

Как получается контрольное поле катушки

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

  1. Поле отдельного кругового витка на оси следует из закона Био — Савара. Расстояние до его плоскости обозначено z.

    B_z(z)=\frac{\mu_0IR^2}{2(R^2+z^2)^{3/2}}
    (6.4)
  2. Просуммируйте витки вдоль катушки: число витков на единицу длины постоянно. В центре интеграл даёт приведённый в уроке эталон.

    B_{\text{ц}}=\frac{N}{l}\int_{-l/2}^{l/2}\frac{\mu_0IR^2}{2(R^2+z^2)^{3/2}}\,\mathrm dz=\frac{\mu_0NI}{2\sqrt{R^2+(l/2)^2}}
    (6.5)
  3. Сравните осевую компоненту в одной и той же точке. Затем удвойте ток: линейность поля и квадратный рост энергии проверяют разные части расчётной цепочки.

Таблица 6.1. Обозначения и единицы
ВеличинаФизический смыслЕдиница СИ
BМагнитная индукцияТл
NЧисло витков с учётом принятого направления1
IТок обмотки; в статике постоянныйА
\mu_0Магнитная постояннаяГн/м
RСредний радиус тонкой катушким
lДлина катушким

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

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

Индукция в центре

Дано. N=100; I=1 \,\text{А}; l=0{,}200 \,\text{м}; R_{\text{ср}}=0{,}025 \,\text{м}.

  1. Вычислим геометрический знаменатель.

    2\sqrt{0{,}025^2+0{,}100^2}=0{,}206155\,\text{м}
    (6.6)
  2. Учитываем конечную длину.

    B_\text{ц}\simeq\frac{4\pi\cdot10^{-7}\cdot100}{0{,}206155}=6{,}0956\cdot10^{-4}\,\text{Тл}
    (6.7)

Ответ: B_{\text{ц}}\approx 0{,}610 \,\text{мТл} при 1 \,\text{А}; при 2 \,\text{А} ожидается около 1{,}219 \,\text{мТл}.

Проверка: Для бесконечного соленоида получилось бы 0{,}628 \,\text{мТл}. Более короткая катушка имеет меньшую центральную индукцию.

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

Полный расчёт и смена тока

Таблица 6.2. Исходные данные практикума
Исходная величинаЗначениеЕдиница
Тип задачиОсесимметричная (Axisymmetric)
Длина катушки200мм
Внутренний / наружный радиус20 / 30мм
Число витков100
Ток цепи1А
Частота0Гц
Радиус внешней полуокружности300мм
Относительная проницаемость воздуха и обмотки1
Удельная электрическая проводимость учебной обмотки0См/м
  1. Создайте магнитную задачу: «Файл → Новая» (File → New), затем «Магнитная задача» (Magnetics Problem). В «Постановке задачи» (Problem Definition) выберите миллиметры (millimeters), осесимметричность (Axisymmetric), частоту 0 \,\text{Гц} и точность линейного решателя 10^{-8}.
  2. Постройте прямоугольник обмотки по точкам (20; -100), (30; -100), (30; 100), (20; 100). Горизонтальная координата — радиус r, вертикальная — координата z. Геометрия занимает только r \ge 0.
  3. Замкните внешнюю область правой полуокружностью радиуса 300 \,\text{мм} с центром (0; 0) и отрезком оси r=0. Обмотка должна полностью находиться внутри области. Суммарный угол двух четвертей внешней дуги — 180^\circ.
  4. В «Свойствах материала» (Material Properties) создайте «Воздух» с относительной проницаемостью 1 и проводимостью 0 и «Учебную обмотку» с теми же параметрами. Разместите метку области (Block Label) в обмотке и отдельную метку воздуха вне её.
  5. В «Свойствах цепи» (Circuit Properties) создайте последовательную цепь (series) «Катушка», задайте ток 1 \,\text{А}. В свойствах метки обмотки назначьте эту цепь и число витков (Turns)+100. Не добавляйте к ней независимо объёмную плотность тока.
  6. Создайте границу «A=0» типа Prescribed A, оставьте коэффициенты A_{0}, A_{1}, A_{2} равными нулю и назначьте её внешней дуге и оси. На оси это соответствует регулярности азимутального потенциала; на удалённой дуге — усечению бесконечной области.
  7. Сохраните исходную модель; создайте и осмотрите сетку, выполните расчёт и откройте результат. В режиме точки снимите B в (0; 0), в результатах цепи (View → Circuit Properties) запишите ток, напряжение и потокосцепление. Сохраните таблицу с единицами отдельно от цветовой карты.
  8. Выполните ещё два расчёта при 2 \,\text{А} и -1 \,\text{А}. Снимайте B_z в одной точке (0;0); запишите также \Psi и W_{\text{м}} с указанием выбранных областей интегрирования.

Ожидаемый результат. Аналитическая оценка B_z(0;0)\approx 0{,}610 \,\text{мТл} для тонкого конечного соленоида. Числа FEMM должны быть измерены учеником; заранее вычисленным полевым результатом эта оценка не является.

Как проверить. Проверьте отношение B_z(2 \,\text{А})/B_z(1 \,\text{А})\approx 2 и B_z(-1 \,\text{А})/B_z(1 \,\text{А})\approx -1; при отклонении сначала проверьте актуальность результата и параметры источника.

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

Что должно измениться после смены тока и нового решения?

  1. Запишите ток и поле в одной контрольной точке для исходного и изменённого режимов.
  2. Сравните результаты при согласованном возбуждении; проверьте знаки и разделите ток витка, ампер-витки и потокосцепление.

С чем сравнить. Опора — закон Ампера и определения цепи и потокосцепления из занятия. Масштабирование по току линейно только для линейного материала при неизменной геометрии.

Итог работы. Две строки результатов и вывод о соответствии изменения поля возбуждению.

Пошаговая настройка расчётного примера

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

Все расчётные примеры курса

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

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

Таблица 6.3. Результаты расчётного образца
ВеличинаРезультатЕдиница
Частота0Гц
Ток цепи Катушка1А
Потокосцепление цепи Катушка95,85 × 10⁻⁶Вб
Энергия поля47,92 × 10⁻⁶Дж
Модуль B в контрольной точке608,5 × 10⁻⁶Тл
Координаты контрольной точки (r, z)0,001; 0мм
Визуализация численных данных FEMM 5: Воздушная катушка. Это экспорт результатов поля, а не снимок окна программы.
Рисунок 6.1. Визуализация численных данных FEMM 5: Воздушная катушка. Это экспорт результатов поля, а не снимок окна программы.

Осесимметричная учебная модель. Независимо проверены изменения сетки и внешнего радиуса; аналитическая модель тонкостенного соленоида служит отдельным приближённым ориентиром.

Ориентиры в интерфейсе
Версия FEMM 5 и дата сборки: окно «О программе».
Рисунок 6.2. Версия FEMM 5 и дата сборки: окно «О программе». Открыть в полном размере
Окно New Document: для магнитной задачи выберите Magnetics Problem.
Рисунок 6.3. Окно New Document: для магнитной задачи выберите Magnetics Problem. Открыть в полном размере
Интерфейс FEMM 5: осесимметричная катушка, геометрические объекты, метки и воздушная область.
Рисунок 6.4. Интерфейс FEMM 5: осесимметричная катушка, геометрические объекты, метки и воздушная область. Открыть в полном размере
Постановка осесимметричной магнитостатической задачи: миллиметры, частота 0 Гц и точность.
Рисунок 6.5. Постановка осесимметричной магнитостатической задачи: миллиметры, частота 0 Гц и точность. Открыть в полном размере
Учебная обмоточная область: относительная проницаемость 1 и нулевая проводимость; источник задаётся цепью.
Рисунок 6.6. Учебная обмоточная область: относительная проницаемость 1 и нулевая проводимость; источник задаётся цепью. Открыть в полном размере
Последовательная цепь Coil: ток 1 А; число витков задаётся отдельно в метке области.
Рисунок 6.7. Последовательная цепь Coil: ток 1 А; число витков задаётся отдельно в метке области. Открыть в полном размере
Свойства метки катушки: материал Winding, цепь Coil и 100 витков.
Рисунок 6.8. Свойства метки катушки: материал Winding, цепь Coil и 100 витков. Открыть в полном размере
Граничное условие A0: нулевой заданный магнитный потенциал на оси и удалённой дуге.
Рисунок 6.9. Граничное условие A0: нулевой заданный магнитный потенциал на оси и удалённой дуге. Открыть в полном размере
Треугольная конечно-элементная сетка катушки в воздухе. В областях проводников элементы мельче, чем у внешней границы расчётной области.
Рисунок 6.10. Треугольная конечно-элементная сетка катушки в воздухе. В областях проводников элементы мельче, чем у внешней границы расчётной области. Открыть в полном размере
Поле воздушной катушки: осесимметричный магнитостатический эталон FEMM 5.
Рисунок 6.11. Поле воздушной катушки: осесимметричный магнитостатический эталон FEMM 5. Открыть в полном размере
Результаты цепи Coil: заданный ток и рассчитанное потокосцепление.
Рисунок 6.12. Результаты цепи Coil: заданный ток и рассчитанное потокосцепление. Открыть в полном размере

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

Повторите опыт при N=200 и I=0{,}5 \,\text{А}. Сравните B_{\text{ц}} и \Psi с исходными N=100, I=1 \,\text{А}.

Критерий выполнения: B_{\text{ц}} практически сохраняется при том же NI; \Psi удваивается, потому что удвоено число витков, сцепляющихся с тем же потоком.
Первоисточники и дополнительное чтение по занятию
  1. [33] — 24–26, 37–52
  2. [20] — 10–22 скачанного ENSEM_TP_FEMM.pdf
Пояснения и типичные ошибки
Почему бесконечный соленоид завышает B_{\text{ц}}?

Он не учитывает выход потока у торцов.

Можно ли получить \Psi из одной центральной точки B?

Только при отдельно обоснованной однородности и одинаковом сцеплении всех витков.

Что меняется при смене знака I?

Поле и потокосцепление меняют знак; энергия линейной катушки сохраняется.

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

  • Сравнивать цвет без фиксированной шкалы.
  • Принимать тонкостенную оценку за фактически выполненный FEMM-расчёт.
  • Снять результат старого тока после редактирования модели.