FEMM 5

Часть I · Интерфейс FEMM 5, геометрия, сетка и математические основы → § 1 · Рабочая среда и документы

Занятие 1.1

Новый интерфейс FEMM 5 и расположение инструментов

Воздушная катушка: N = 100 витков · I = 1 А; основные размеры, модуль магнитной индукции b; Ψ = 95,85 мкВб.
Воздушная катушка1 — Обмотка · 2 — r₂ = 30 ммРасчётный образец. Осесимметричное сечение; показана одна сторона обмотки. Размер r₁ отсчитывается от оси до внутренней поверхности, r₂ — до наружной.Открыть SVG целикомРезультат отдельноИсходная учебная модель

Цель и задачи

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

Задачи

  • Найти режимы построения и выбора объектов.
  • Проверить тип задачи, единицы и частоту.
Что нужно вспомнить перед началом
  • Начальные знания о токе, магнитном поле и единицах СИ.

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

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

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

Инженерный вопрос. Какие исходные данные определяют смысл будущего результата катушки?

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

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

Результат этого этапа — проверенная постановка или воспроизводимая процедура перед серией расчётов.

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

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

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

  3. Ответить на вопрос по своим результатам. Сопоставьте запись с повторно открытой моделью; объясните одно последствие неверного параметра.

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

Результат этого занятия. Краткий описание модели с объяснением, где задан ток и какое тело представляет сечение.

Данные именно этого опыта. Запишите тип задачи, единицы, источник, число витков и объёмный множитель подготовленной катушки.

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

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

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

Независимый контроль. Восстановите исходные данные по повторно открытому файлу и проверьте, достаточно ли сохранённого описания для воспроизведения модели другим исполнителем.

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

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

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

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

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

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

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

Расчёт магнитных систем электромагнитов классическим и численным методом (ОНПУ,2010)

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

Лабораторные работы по методу: Трансформатор от потокосцеплений до схемы замещения и нагрузочной характеристики; Асинхронный двигатель от поля в зазоре до параметров и механической характеристики; Генератор с постоянными магнитами от конструкции ротора до ЭДС и нагрузочного напряжения; Выбор обмотки и расчет момента в кольцевом зазоре; Контрольный опыт: от условия задачи к инженерному выводу; Опыт Эрстеда: момент стрелки, эквивалентный ток и движение; Измерительная рамка: от поля провода к ЭДС; Аналогии полей: проводимость, дуальная сеть и диффузия; От магнитной цепи к карте \Psi (i,\theta ) и моменту SRM

Получаемые величины. \Phi, Bсерд и Bзазора; тяговая сила; секущая индуктивность; отклонение магнитной цепи от поля и численная устойчивость

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

Страницы источника. PDF 1–48 полностью: 1–3 титул/авторы/содержание; 4–13 теория; 14–34 FEMM 4.0.1; 35–43 П-электромагнит и сравнение; 44–45 выводы/литература; 46–48 варианты.

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

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

Знакомство с 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. (1.1) · Источник находится в свойствах, а не на рисунке
  2. (1.2) · Единицы и подпись результата
  3. (1.3) · От параметров интерфейса к размерной модели
  4. (1.4) · Численный пример

Как устроен курс

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

Кому адресован курс

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

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

Как построен каждый урок

Все уроки курса следуют одной логике изложения, которую полезно знать заранее — она подсказывает, где искать ответ на возникший вопрос.

  1. Явление. Урок начинается с физической ситуации: заряд в диэлектрике, ток в проводнике, нагретая стенка. Это отвечает на вопрос «зачем».
  2. Геометрический образ. Понятие рисуется прежде, чем записывается: градиент — как направление крутого склона, поток — как «количество линий», проходящих сквозь поверхность, ротор — как вращение маленького колёсика.
  3. Уравнения. Даётся строгое определение (блок «Определение») и вывод (блок «Утверждение») с явным указанием всех обозначений и единиц.
  4. Допущения и границы применимости. Для каждой формулы оговаривается, в каких условиях она верна: статика или квазистатика, линейная или нелинейная среда, гладкое поле или поле с особой точкой.
  5. Визуализация. Раздел «Лаборатория» описывает, что показано в интерактивной сцене, какие параметры можно менять и на какие вопросы она отвечает.
  6. Пример. Разобранный пример с числами в СИ (блок «Пример»), результат которого следует непосредственно из приведённых формул — его можно и нужно пересчитать самостоятельно.
  7. Интерпретация. Что означает полученное число для инженера: допустимо ли оно, как оно изменится при изменении параметров, какие типичные ошибки (блок «Типичные ошибки») приводят к неверному результату.

Урок завершается вопросами для самопроверки (ответы к части из них скрыты и раскрываются по щелчку) и кратким перечнем «Итоги» со ссылкой на следующий урок. Заголовок, цели обучения и список пререквизитов формируются автоматически и показаны над текстом.

Карта курса

Курс разбит на десять модулей. Уроки, снабжённые лабораторией, отмечены знаком ⚙.

M0. Введение

M1. Математический язык полей

M2. Уравнения поля и краевые задачи

M3. Электростатика

M4. Поля электрического тока

M5. Магнитостатика

M6. Гармонические магнитные поля

M7. Интегральные величины и инженерные параметры

M8. Тепловые поля

M9. Связанные задачи и численное моделирование

Модули M1 и M2 — фундамент: без них невозможно читать остальные. Модули M3–M6 можно изучать в указанном порядке, модуль M8 (тепловые поля) опирается только на градиент и дивергенцию из M1 и может быть пройден сразу после них. Модуль M9 собирает всё вместе и объясняет, как те же уравнения решаются численно.

Соглашения об обозначениях и единицах

Единицы и запись чисел

Все формулы и примеры курса записаны в Международной системе единиц (СИ): метр, килограмм, секунда, ампер, кельвин. Производные единицы — вольт, тесла, вебер, ватт и другие — выражаются через основные без дополнительных множителей, поэтому в формулы не входят коэффициенты вида 4π4\pi или 1/c1/c, характерные для гауссовой системы. Числа записываются в русском формате: десятичный разделитель — запятая (0,50{,}5; 1,2⋅10−31{,}2\cdot10^{-3}), разряды групп не отделяются. Электрическая и магнитная постоянные: ε0=8,854⋅10−12\varepsilon_0=8{,}854\cdot10^{-12} Ф/м, μ0=4π⋅10−7\mu_0=4\pi\cdot10^{-7} Гн/м ≈1,2566⋅10−6\approx1{,}2566\cdot10^{-6} Гн/м [CODATA].

Температура в формулах — абсолютная, в кельвинах; градусы Цельсия допускаются только при описании «бытовых» температур («стенка нагрета до 80 °C»), а в законе Стефана–Больцмана и других формулах излучения — только кельвины. Буква ρ\rho без индекса обозначает объёмную плотность заряда (Кл/м³); удельное электрическое сопротивление записывается как ρэ\rho_{\mathrm э} (Ом·м), плотность массы — ρm\rho_m (кг/м³).

Векторы и операторы

Векторные величины набираются полужирным шрифтом: E\uE, B\uB, J\uJ, q\uq. Модуль вектора — обычным: E=∣E∣E=|\uE|. Единичные базисные векторы — ex,ey,ez\vec e_x,\vec e_y,\vec e_z в декартовой системе и er,eα,ez\vec e_r,\vec e_\alpha,\vec e_z в цилиндрической; радиус-вектор точки — r\ur, единичная нормаль — n\un. Дифференциальные операторы записываются словами и полностью равнозначны записи через оператор «набла» ∇\nabla: ∇φ≡∇φ,∇⋅F≡∇⋅F,∇×F≡∇×F,Δφ≡∇2φ≡∇⋅∇φ.\nabla \varphi\equiv\nabla\varphi,\qquad \nabla\cdot \vec F\equiv\nabla\cdot\vec F,\qquad \nabla\times \vec F\equiv\nabla\times\vec F,\qquad \Delta\varphi\equiv\nabla^2\varphi\equiv\nabla\cdot \nabla \varphi . Оператор ротора в англоязычной литературе обозначается ∇×\nabla\times; в курсе принято традиционное для русской школы ∇×\nabla\times. Дифференциал записывается прямой буквой: dx\dd x, dl\dl, dS\dS, dV\dV.

Главный источник обозначений для нас — учебники [Бессонов] и [Демирчян]; при расхождениях с ними в тексте даётся явное пояснение. Полный перечень символов с единицами и ссылками на уроки собран в глоссарии.

Гармонические процессы и комплексные амплитуды

Значительная часть курса (модули M4, M6, M7) посвящена полям, изменяющимся во времени по синусоидальному закону. Для них применяется метод комплексных амплитуд, и здесь важно зафиксировать соглашения, поскольку в литературе они различаются.

Соглашение о комплексных амплитудах

Мнимая единица обозначается j\ii, j2=−1\ii^2=-1. Гармоническая величина x(t)=Xmcos⁡(ωt+ψ)x(t)=X_m\cos(\omega t+\psi) с амплитудой XmX_m, угловой частотой ω=2πf\omega=2\pi f и начальной фазой ψ\psi представляется комплексной амплитудой (фазором) X˙=Xmejψ,x(t)=Re⁡{X˙ejωt}.\cplx X=X_m e^{\ii\psi},\qquad x(t)=\Re\{\cplx X e^{\ii\omega t}\}. Временная зависимость — ejωte^{\ii\omega t} (не e−jωte^{-\ii\omega t}); комплексные амплитуды — амплитудные, а не действующие значения. Действующее значение X=Xm/2X=X_m/\sqrt2.

Два пояснения. Во-первых, буква j\ii выбрана вместо принятой в математике ii потому, что в электротехнике i(t)i(t) — мгновенный ток; при записи jωLI˙\ii\omega L\cplx I путаница исключена (полужирная J\uJ — плотность тока — с прямой j\ii не смешивается). Во-вторых, выбор ejωte^{\ii\omega t} означает, что дифференцирование по времени сводится к умножению на jω\ii\omega: dx/dt↔jωX˙\dd x/\dd t\leftrightarrow\ii\omega\cplx X. В курсах теоретической физики, например [Ландау] и [Джексон], принята зависимость e−iωte^{-i\omega t}; их формулы переходят в наши заменой i→−ji\to-\ii (комплексным сопряжением). Так, коэффициент затухания поля в проводнике записывается у них через −i=(1−i)/2\sqrt{-i}=(1-i)/\sqrt2, а у нас — через j=(1+j)/2\sqrt{\ii}=(1+\ii)/\sqrt2, тогда как сама глубина проникновения δ\delta вещественна и от соглашения не зависит. При сопоставлении с литературой на это следует обращать внимание.

Из амплитудного соглашения вытекает запись средней мощности: P=12Re⁡{U˙I˙∗}P=\tfrac12\Re\{\cplx U\cplx I^*\}, где звёздочка — комплексное сопряжение, а множитель 12\tfrac12 появляется именно потому, что амплитуды не действующие.

Соглашения об ориентации

Многие величины в теории поля — поток, циркуляция, ЭДС — имеют знак, зависящий от выбора направления. Чтобы знаки в курсе были согласованы с лабораториями, принимаются три правила.

Правила ориентации

  1. Нормаль к замкнутой поверхности — внешняя. Поток ∮SF⋅n dS\oint_S\vec F\cdot\un\,\dS положителен, если поле «выходит» из объёма.
  2. Правило правого винта. Для поверхности с краем направление обхода контура LL и нормаль n\un связаны так: если смотреть с конца вектора n\un, обход совершается против хода часовой стрелки. Иначе: если пальцы правой руки указывают направление обхода, большой палец даёт n\un.
  3. ЭДС индукции e=−dΨ/dte=-\dd\Psi/\dd t, где потокосцепление Ψ\Psi вычислено с нормалью, согласованной по правилу 2 с положительным направлением обхода контура (которое принимается и за положительное направление ЭДС и тока).

Знак «минус» в третьем правиле — закон Ленца, но только при согласованной ориентации: если нормаль выбрать противоположной, знак формально сменится, а физика, разумеется, нет. В лаборатории потокосцепления стрелка положительного направления обхода показана на сцене явно.

Как пользоваться лабораториями

Интерактивные лаборатории открываются прямо в тексте урока в месте, где они нужны. Трёхмерную сцену можно вращать (перетаскивание мышью или одним пальцем), сдвигать (перетаскивание правой кнопкой или двумя пальцами) и масштабировать (колесо прокрутки или щипок). На панели параметров справа (на узких экранах — снизу) каждый параметр представлен ползунком, числовым полем и единицей измерения; кнопка «Сброс» возвращает исходные значения, а заранее подобранные варианты — например, набор полей с положительной, отрицательной и нулевой дивергенцией — выбираются из списка. Пробник — перемещаемая точка, в которой выводятся локальные значения: компоненты поля, модуль градиента, производная по направлению. У каждой сцены есть текстовая альтернатива, пригодная для чтения с экрана и для случаев, когда трёхмерная графика недоступна; все элементы управления доступны с клавиатуры.

Происхождение данных

Каждая лаборатория снабжена пометкой о том, откуда взяты показываемые числа:

  • аналитическое решение — вычислено по замкнутым формулам, приведённым в уроке;
  • численное решение — получено дискретизацией (например, методом конечных элементов) с указанием метода и размера сетки; такие данные имеют погрешность, которая обсуждается в уроке о сходимости (модуль M9);
  • качественная иллюстрация — картина передаёт характер явления, а не количественный результат;
  • предвычисленные данные — рассчитаны заранее для набора параметров и воспроизводятся без пересчёта (например, временная шкала нестационарного прогрева). В графиках происхождение каждой кривой указано в легенде (например, «Кранк–Николсон» и «аналитический ряд»), а качественные элементы — линии поля, оттенки фона — изображаются полупрозрачно.

Линии поля — не траектории

Линии поля (силовые линии, линии тока, линии магнитной индукции) — это геометрическое представление направления вектора: касательная к линии в каждой точке совпадает с направлением поля. Они не являются траекториями каких-либо частиц: заряд, обладающий инерцией, в общем случае не движется вдоль линии E\uE; линии B\uB никуда не «текут», хотя в анимациях их иногда рисуют движущимися. Густота линий отражает модуль поля лишь при специальном построении (в плоских задачах — при равном шаге потока между соседними линиями), поэтому выводы о величине поля в курсе всегда делаются по числам, а не по рисунку.

Что входит в курс и что остаётся за его рамками

Курс охватывает электростатику, поля постоянного и гармонического тока, магнитостатику, квазистационарные магнитные поля (вихревые токи, скин-эффект), тепловые поля в твёрдых телах и методы перехода от распределённого поля к сосредоточенным параметрам — ёмкости, индуктивности, сопротивлению, тепловому сопротивлению. Отдельно рассматривается, как эти задачи решаются численно и как проверять результат.

За рамками курса остаются электромагнитные волны, излучение антенн и волноводы: соответствующие явления упоминаются только там, где нужно указать границы применимости статических и квазистатических приближений (условие малости размеров системы по сравнению с длиной волны). Конвекция и тепловое излучение рассматриваются не как самостоятельные механизмы переноса, а только как граничные модели — коэффициент теплоотдачи α\alpha и закон Стефана–Больцмана на поверхности тела; движение жидкости и перенос излучения в объёме не изучаются. Магнитные материалы описываются однозначной кривой B(H)B(H) без учёта гистерезиса; потери на гистерезис обсуждаются качественно.

Глоссарий и литература

Глоссарий содержит все обозначения курса с единицами измерения, кратким определением и ссылками на уроки, где символ вводится. Ссылки на него встроены в текст: например, дивергенция. Библиография включает только реально существующие издания; в тексте они обозначаются фамилией первого автора в квадратных скобках — [Бессонов], [Тамм], [Кочин]. Для векторного анализа рекомендуются [Кочин] и [Борисенко], для строгого изложения интегральных теорем — [Фихтенгольц], для уравнений поля — [Бессонов], [Демирчян] и [Тамм], для теплопроводности — [Лыков] и [Исаченко].

Вопросы для самопроверки

  1. Чем комплексная амплитуда I˙=5ejπ/6\cplx I=5e^{\ii\pi/6} А отличается от действующего значения тока? Чему равно действующее значение?
  2. Запишите d2x/dt2\dd^2x/\dd t^2 для гармонической величины x(t)x(t) через её комплексную амплитуду. Как изменился бы ответ при соглашении e−jωte^{-\ii\omega t}?
  3. Через замкнутую поверхность поток ∮J⋅n dS\oint\uJ\cdot\un\,\dS оказался отрицательным. Заряд внутри объёма растёт или убывает?
  4. Почему пометка «численное решение» в лаборатории сопровождается указанием размера сетки, а пометка «аналитическое решение» — нет?
  5. Является ли линия магнитной индукции траекторией движения какого-либо объекта? Приведите довод.

Ответы к пунктам 1–3

1. Амплитуда тока 5 А, начальная фаза 30°; мгновенный ток i(t)=5cos⁡(ωt+π/6)i(t)=5\cos(\omega t+\pi/6) А. Действующее значение I=5/2≈3,54I=5/\sqrt2\approx3{,}54 А.

2. d2x/dt2↔(jω)2X˙=−ω2X˙\dd^2x/\dd t^2\leftrightarrow(\ii\omega)^2\cplx X=-\omega^2\cplx X. При соглашении e−jωte^{-\ii\omega t} результат тот же ((−jω)2=−ω2(-\ii\omega)^2=-\omega^2), но первая производная сменила бы знак: −jωX˙-\ii\omega\cplx X.

3. Поток плотности тока через внешнюю нормаль — это ток, вытекающий из объёма. Отрицательный поток означает, что ток втекает, и заряд внутри растёт: dQ/dt=−∮J⋅n dS>0\dd Q/\dd t=-\oint\uJ\cdot\un\,\dS>0.

Итоги

  • Каждый урок проходит путь «явление → геометрия → уравнения → допущения → визуализация → пример → интерпретация».
  • Единицы — СИ, числа — в русском формате, векторы — полужирные, операторы ∇\nabla, ∇⋅\nabla\cdot, ∇×\nabla\times равнозначны ∇\nabla, ∇⋅\nabla\cdot, ∇×\nabla\times.
  • Гармонические процессы описываются амплитудными комплексными амплитудами при зависимости ejωte^{\ii\omega t}; мнимая единица — j\ii.
  • Нормаль к замкнутой поверхности внешняя, обход и нормаль связаны правилом правого винта, e=−dΨ/dte=-\dd\Psi/\dd t.
  • Лаборатории помечают происхождение данных; линии поля — средство изображения направления, а не траектории.

Далее — скалярные и векторные поля и системы координат: язык, на котором будут записаны все последующие уравнения.

Уточнения постановки и расчётной методики FEMM 5

Модель состоит из объектов с разными функциями

Расчётная модель создаётся в препроцессоре, а полученное решение исследуется в окне обработки результатов. Каждый геометрический объект выполняет определённую функцию: узел (Node) задаёт точку, отрезок (Segment) соединяет два узла прямой, дуга (Arc Segment) описывает криволинейную границу. Метка области (Block Label) связывает замкнутую часть плоскости со свойствами материала. Поэтому построение контура и назначение материала представляют собой разные этапы: свойства линии не определяют материал её внутренней области, а метка вне замкнутого контура не формирует расчётный объём. Перед решением следует последовательно проверить объекты каждого типа и убедиться, что выделены именно те элементы, свойства которых требуется задать.

Постановка задачи предшествует размерам

В окне Problem Definition выбирают систему единиц геометрии, Planar либо Axisymmetric и частоту. Для начала используется магнитостатика: частота равна нулю. Число 20 в координате означает 20 \,\text{мм} только после выбора миллиметров. В плоскопараллельной модели поле рассчитывается по сечению, а глубина определяет объёмные интегралы. В осесимметричной модели сечение вращается вокруг вертикальной оси; горизонтальная координата является радиусом. Это разные физические предположения, а не два варианта отображения одной карты поля.

Источник находится в свойствах, а не на рисунке

Последовательная цепь задаёт ток одного витка; число витков со знаком в метке обмотки определяет полный ток через сечение. Воздушная катушка и окружающий воздух имеют одинаковую магнитную проницаемость, но только катушка связана с источником. Название «Медь» не создаёт возбуждения. Для первого занятия удельная электрическая проводимость учебной обмотки принимается равной нулю: так выделяется поле заданного тока без дополнительных вихретоковых эффектов. Это идеализация многовитковой обмотки, а не утверждение о реальной проводимости меди.

\mathcal F_C=\oint_C\mathbf H\cdot\mathrm d\mathbf l=NI
(1.1)
МДС определяется циркуляцией H и выражается в А. Равенство NI следует из закона Ампера в магнитостатике для контура, сцепленного с N одинаково ориентированными последовательными витками; I — ток одного витка.

Замыкание геометрии и внешняя область

Линии магнитной индукции замкнуты и обычно проходят через окружающий воздух. Поэтому внешнюю область необходимо включить в задачу. Все границы должны образовывать однозначные замкнутые области без щелей и случайных пересечений. Каждая независимая область получает свою метку. Удалённая граница A=0 в первом опыте является приближением бесконечного пространства. Её расстояние проверяется позже отдельно от сетки. Отсутствующее в геометрии воздушное окно нельзя исправить увеличением числа конечных элементов.

Сетка, решение и просмотр — три разных действия

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

Единицы и подпись результата

В интерфейсе числовой ввод обычно использует точку как десятичный разделитель; в русском отчёте применяется запятая. Магнитная индукция измеряется в теслах, поток — в веберах, потокосцепление — в веберах, индуктивность — в генри. Формально число витков безразмерно; уточнение «виток» помогает смысловому различению, но не является новой единицей СИ. Для геометрии удобно использовать миллиметры, однако в аналитических формулах нужно перейти к метрам. Проверьте, как смена масштаба или глубины меняет геометрические множители.

20\,\text{мм}=0{,}020\,\text{м},\qquad S=bd
(1.2)
Площадь поперечного сечения S в \text{м}^2 получается при ширине b и глубине d в метрах.

Как читать модель в понятном порядке

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

От параметров интерфейса к размерной модели

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

  1. Переведите геометрические размеры в метры для ручной оценки. Смена единицы в формуле требует изменения численного значения длины, площади и объёма.

  2. МДС выбранного замкнутого контура определяется циркуляцией \mathcal F_C=\oint_C\mathbf H\cdot\mathrm d\mathbf l. В магнитостатике контур, сцепленный с N одинаково ориентированными последовательными витками, даёт по закону Ампера \mathcal F_C=NI; ток каждого провода при этом остаётся I.

    \mathcal F_C=\oint_C\mathbf H\cdot\mathrm d\mathbf l=NI,\qquad [\mathcal F_C]=\mathrm A
    (1.3)
  3. Для плоской модели вычислите площадь сечения потока как произведение ширины на глубину. Удвоение глубины при неизменном поле удваивает поток, энергию и силу.

Таблица 1.1. Обозначения и единицы
ВеличинаФизический смыслЕдиница СИ
NЧисло витков с учётом принятого направления1
IТок обмотки; в статике постоянныйА
dАктивная глубина плоской моделим
\PhiМагнитный поток через выбранную поверхностьВб

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

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

Проверка размеров и магнитодвижущей силы

Дано. N=100; I=1 \,\text{А}; ширина пакета 20 \,\text{мм}; глубина 20 \,\text{мм}.

  1. В магнитостатике контур, охватывающий сто последовательных витков с током 1 \,\text{А} и согласованными знаками, имеет МДС 100 \,\text{А} по закону Ампера.

  2. Перевод размеров в метры необходим до вычисления площади.

    S=0{,}020\cdot0{,}020=4\cdot10^{-4}\,\text{м}^2
    (1.4)

Ответ: \mathcal F=100 \,\text{А}; S=400 \,\text{мм}^2=0{,}0004 \,\text{м}^2.

Проверка: Если ошибочно подставить 20 вместо 0{,}020 \,\text{м}, площадь получится завышенной в миллион раз.

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

Описание подготовленной воздушной катушки

Таблица 1.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. Вместо запуска расчёта последовательно покажите на экране: окно постановки; метку обмотки; материал; цепь; внешний воздух. Запишите, чем отличаются координата узла, число витков и ток цепи.

Ожидаемый результат. Получена замкнутая осесимметричная геометрия с одним заданным источником NI=100 \,\text{А}; расчётный результат на этом занятии не требуется.

Как проверить. Повторно откройте сохранённую модель и восстановите описание её исходных данных без подсказки преподавателя.

Оценка результата

Какие исходные данные определяют смысл будущего результата катушки?

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

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

Итог работы. Краткий описание модели с объяснением, где задан ток и какое тело представляет сечение.

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

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

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

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

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

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

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

Ориентиры в интерфейсе
Версия FEMM 5 и дата сборки: окно «О программе».
Рисунок 1.2. Версия FEMM 5 и дата сборки: окно «О программе». Открыть в полном размере
Окно New Document: для магнитной задачи выберите Magnetics Problem.
Рисунок 1.3. Окно New Document: для магнитной задачи выберите Magnetics Problem. Открыть в полном размере
Интерфейс FEMM 5: осесимметричная катушка, геометрические объекты, метки и воздушная область.
Рисунок 1.4. Интерфейс FEMM 5: осесимметричная катушка, геометрические объекты, метки и воздушная область. Открыть в полном размере
Постановка осесимметричной магнитостатической задачи: миллиметры, частота 0 Гц и точность.
Рисунок 1.5. Постановка осесимметричной магнитостатической задачи: миллиметры, частота 0 Гц и точность. Открыть в полном размере
Учебная обмоточная область: относительная проницаемость 1 и нулевая проводимость; источник задаётся цепью.
Рисунок 1.6. Учебная обмоточная область: относительная проницаемость 1 и нулевая проводимость; источник задаётся цепью. Открыть в полном размере
Последовательная цепь Coil: ток 1 А; число витков задаётся отдельно в метке области.
Рисунок 1.7. Последовательная цепь Coil: ток 1 А; число витков задаётся отдельно в метке области. Открыть в полном размере
Свойства метки катушки: материал Winding, цепь Coil и 100 витков.
Рисунок 1.8. Свойства метки катушки: материал Winding, цепь Coil и 100 витков. Открыть в полном размере
Поле воздушной катушки: осесимметричный магнитостатический эталон FEMM 5.
Рисунок 1.9. Поле воздушной катушки: осесимметричный магнитостатический эталон FEMM 5. Открыть в полном размере

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

Составьте одностраничное описание модели и укажите, где находятся все пять исходных параметров: тип, единицы, частота, ток, витки.

Критерий выполнения: Все значения читаются из повторно открытой модели; единицы и направление тока записаны явно.

Литература

Полные библиографические описания

Пояснения и типичные ошибки
Где задаётся материал?

В свойствах метки замкнутой области, а не у линии контура.

Достаточно ли корректной геометрии для правильного задания источника поля?

Нет. Нужно проверить связь обмотки с цепью и число витков со знаком.

Зачем до построения выбирать единицы?

Одни и те же координаты дают разные реальные размеры при разных единицах.

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

  • Считать область железом только потому, что её контур похож на магнитопровод.
  • Забыть метку воздуха в замкнутом окне.
  • Смешать миллиметры в геометрии с метрами в формулах.