Двухобмоточный магнитопровод: первичное возбуждение1 — N₁ · 2 — N₂Расчётный образец. Линейная сталь μᵣ = 1000; глубина 20 мм. b и h — наружные размеры сердечника между противоположными поверхностями.Открыть SVG целикомРезультат отдельноИсходная учебная модель
Цель и задачи
Цель. Выбрать способ измерения, соответствующий нужной величине.
Если одно из понятий незнакомо, вернитесь к предыдущему занятию и воспроизведите его численный пример. После этого переходите к новой постановке.
Место занятия в инженерной методике
Расчётная лабораторная работа · Выбор наблюдаемых величин и проверка интегрального закона
Инженерный вопрос. Почему точка, профиль, контур и область дают разные результаты?
Эти действия входят в подготовку расчёта любого устройства: они устанавливают связь источника, геометрии и наблюдаемой величины. Их результат позволяет выбрать допустимую постановку и способ извлечения будущей характеристики.
Последовательность решения, отчёт и критерии результата
Результат получают из наблюдаемых величин и проверяют в условиях данного практикума.
Последовательность решения
Обосновать постановку и способ измерения. Определите ориентацию контура или поверхности, охватываемые источники и геометрический множитель. Закон связывает измеряемую величину с независимым источником.
Поток и потокосцепление. Проверьте знаки и единицы интегрального закона на разобранной геометрии.
Четыре способа измерения. Выполните действия практикума по его исходной таблице. Извлеките предусмотренные практикумом точечные, контурные и объёмные величины.
Результат этого занятия. Таблица величин, единиц и областей измерения без смешения плотности и интеграла.
Данные именно этого опыта. Извлеките предусмотренные практикумом точечные, контурные и объёмные величины.
Оформление наблюдений. Контур или область интегрирования, ориентация, точечные и интегральные величины раздельно, источник и невязка закона.
Как обосновать результат
Проверка из практикума. Каждая величина имеет единицу и определение области; для полной энергии не пропущен воздух.
Независимый контроль. Сопоставьте интеграл с источником или суммой потоков, а не с повторным чтением того же интеграла другим именем.
Критерий выполнения. Другой учащийся может воспроизвести измерения по таблице без угадывания выбранных областей.
Точность и расхождения. Нормируйте невязку на физический масштаб источника или потока, объявленный до сравнения. Для баланса около нуля приводите также абсолютную невязку в единицах закона.
Основание сравнения. Используйте интегральный закон и ориентацию контура из теории занятия. Плотность, интеграл по сечению и полный объёмный результат сравнивают после учёта геометрического множителя.
Вывод для устройства
Объясните, даёт ли выбранный интеграл полный параметр устройства или только локальную характеристику его сечения.
Границы вывода. Простой эталон проверяет конкретное действие или физический закон; перенос на устройство требует заново согласовать геометрию, материалы, режим и объёмный множитель.
Дополнительный расчётный метод
Векторная алгебра и векторный магнитный потенциал
Провод с током даёт количественный эталон B(r), A(r), циркуляции и потока. Через него проверяют калибровку потенциала, локальные производные и допустимость скалярного магнитного потенциала.
Получаемые величины.B_r и B\tau прямого провода; разности A и поток на единицу длины; циркуляция H; ошибка численной производной и сеточной аппроксимации
Различить скалярное и векторное поле и вычисляемые операции.
Сопоставить аналитическое поле внутри и снаружи провода с точками и контурами FEMM.
Изменить отсчёт A, сохранив B и разности потенциала.
Проверить скачок скалярного магнитного потенциала вокруг тока и объяснить необходимость разреза области.
Страницы источника. PDF 1–31 полностью: 1–12 операции и потенциалы; 13–23 поле прямого провода; 24–28 поток, Биo–Савар, уравнения потенциала; 29–31 разностная схема.
Переход от распределения поля к инженерным характеристикам требует выбора измеряемой величины и соответствующего интеграла. Значения в точках, интегралы по контуру и объёмные интегралы рассматриваются с учётом их физического смысла и единиц.
Соглашения об обозначениях и единицах
\(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\). При постоянных магнитах и взаимном возбуждении сначала учитывают ненулевое смещение потокосцепления.
Положительные направления, оси и специальные индексы определяются в конкретной задаче. Единицу берите из условия или заголовка таблицы; перед подстановкой в формулу СИ переводите миллиметры в метры, миллигенри в генри и градусы в радианы.
Численные результаты обычно показаны с четырьмя значащими цифрами; для сравнения близких значений в таблицах сходимости сохранено больше цифр. Полные значения остаются в исходных расчётных данных. Малые остатки вычисления координат на осях показаны как ноль; малые физические величины сохраняются в показательной записи.
Решение полевой задачи — это функции E(r), B(r), J(r), T(r): бесконечно много чисел. Инженеру же нужны несколько: ёмкость конденсатора, индуктивность катушки, сопротивление шины, потери, сила на якорь. Все эти величины — интегралы поля по некоторой области: кривой, поверхности или объёму. Переход «поле → интеграл → параметр» и есть мост между теорией поля и теорией цепей, где те же C, L, R выступают уже как готовые элементы схемы.
Каждая интегральная величина требует четырёх уточнений, без которых её число лишено смысла: что интегрируется, по какой области, с какой ориентацией (нормали или обхода) и при каких допущениях результат не зависит от произвола выбора области. Этот урок — систематический каталог таких уточнений; предыдущие модули дали отдельные определения, здесь они собраны и сопоставлены.
Потоки, заряд, ток
Поток электрического смещения и заряд
Ψэ=∫SD⋅ndS[Кл],Q=∫VρdV[Кл].
Для замкнутой поверхности S с внешней нормалью теорема Гаусса даёт Ψэ=Q — свободный заряд внутри. Отсюда способ измерить заряд электрода по полю: охватить электрод любой замкнутой поверхностью, лежащей в диэлектрике, и проинтегрировать нормальную составляющую D. Результат не зависит от выбора поверхности, пока она не захватывает других зарядов.
Магнитный поток и потокосцепление
Φ=∫SB⋅ndS[Вб],Ψ=k=1∑NΦk[Вб].
Поверхность S опирается на контур витка, нормаль согласована с положительным направлением тока в нём по правилу правого винта. Так как ∇⋅B=0, поток не зависит от выбора натянутой поверхности — только от контура. Потокосцепление обмотки — сумма потоков всех витков; равенство Ψ=NΦ справедливо лишь при одинаковых потоках витков (см. поток и потокосцепление). Через векторный потенциал Φ=∮LA⋅dl.
Ток
I=∫SJ⋅ndS[А].
Поверхность S пересекает проводник; её ориентация задаёт знак тока. В стационарном режиме (∇⋅J=0) ток один и тот же через любое сечение одного проводника. При переменных полях это верно лишь в МКС-приближении, где ток смещения пренебрежим; иначе поверхности, охватывающие заряд, дают разные значения I.
Напряжение
Напряжение
UAB=∫ABE⋅dl[В]
— линейный интеграл вдоль ориентированного пути от A к B. В статическом поле ∇×E=0, интеграл не зависит от пути и равен разности потенциалов φA−φB. В переменном поле ∮E⋅dl=−dΦ/dt=0, и напряжение между двумя точками зависит от пути: два пути A→B отличаются на ЭДС контура, образованного ими. Поэтому «напряжение на зажимах» переменного устройства имеет смысл только вместе с указанием, как проложены измерительные провода, а в теории цепей — при допущении, что вне элементов переменный поток отсутствует.
Ёмкость
Ёмкость двухэлектродной системы
Для двух проводников с зарядами +Q и −Q в линейном диэлектрике
C=UQ[Ф],
где U — напряжение между ними. Линейность среды гарантирует, что отношение не зависит от U. Через энергию: C=2Wэ/U2=Q2/(2Wэ).
При n проводниках заряд каждого линейно зависит от всех потенциалов: Qi=∑jcijφj. Матрица коэффициентов cij симметрична (cij=cji — следствие взаимности, доказываемой через теорему Грина), диагональные элементы положительны, внедиагональные отрицательны. Эквивалентно, систему описывают частичными ёмкостями между каждой парой проводников и между каждым проводником и землёй — именно они входят в схему замещения многопроводной линии. Единица везде — фарад.
Индуктивности
Собственная индуктивность
Для контура (обмотки) с током I в линейной среде
L=IΨ[Гн],эквивалентноL=I22Wм,
где Ψ — потокосцепление собственного поля, Wм — энергия магнитного поля во всём пространстве. Два определения совпадают только в линейной системе; в нелинейной вводят секущую и дифференциальную индуктивности (см. следующий урок). Для тонких проводов определение через поток требует оговорки, где проходит контур внутри провода; определение через энергию её не требует и автоматически включает внутреннюю индуктивность.
Взаимная индуктивность
Потокосцепление контура 1, создаваемое током контура 2:
M12=I2Ψ12[Гн].Теорема взаимности:M12=M21. Для тонких контуров в немагнитной среде это видно из формулы Неймана
M12=4πμ0∮L1∮L2∣r1−r2∣dl1⋅dl2,
симметричной относительно перестановки контуров [Калантаров]. Знак M зависит от выбора положительных направлений обхода в обоих контурах.
Для n обмоток Ψi=∑jLijIj; матрица индуктивностей симметрична, Lii=Li, Lij=Mij. Энергия Wм=21∑i,jLijIiIj положительна при любых токах, откуда M122≤L1L2 и коэффициент связи k=M12/L1L2≤1.
Сопротивление, мощность, импеданс
Сопротивление
R=IU=I2P[Ом],G=R1[См],
где U — напряжение между электродами проводника, I — ток через любое его сечение, P=∫VJ⋅EdV — мощность потерь Джоуля. Два определения совпадают при линейном законе Ома J=σE; определение через мощность удобнее в численных моделях, так как объёмный интеграл сходится быстрее линейного.
Мощность потерь
P=∫VJ⋅EdV=∫VσE2dV[Вт]
— мгновенная мощность. Для гармонических полей с амплитудными фазорами средняя мощность
P=21∫VRe{J˙⋅E˙∗}dV=21∫Vσ∣J˙∣2dV
(см. фазоры). Пропуск множителя 21 — ошибка вдвое.
Импеданс
Для гармонического режима отношение комплексных амплитуд напряжения и тока
Z=I˙U˙=R+jX[Ом].
Активная часть связана со средними потерями: R=2P/∣I˙∣2; реактивная — со средней запасённой энергией: для магнитной системы X=ωL, где L=4⟨Wм⟩/∣I˙∣2. Поскольку напряжение в переменном поле зависит от пути, Z привязан к определённому расположению зажимов. При скин-эффекте R и L зависят от частоты (вихревые токи).
Энергия и тепловой поток
Энергия поля (линейная среда)
Wэ=21∫VE⋅DdV,Wм=21∫VB⋅HdV[Дж],
интеграл — по всему пространству, где поле отлично от нуля. Множитель 21 справедлив лишь для линейных сред; общие выражения ∫H⋅dB рассматриваются в уроке об энергии и коэнергии.
Тепловой поток
P=∫Sq⋅ndS[Вт]
— мощность, уходящая через поверхность S в направлении нормали. Для замкнутой поверхности в стационарном режиме равен суммарному тепловыделению внутри: ∮q⋅ndS=∫VqVdV. Это тепловой аналог теоремы Гаусса и основа проверки энергетического баланса.
Силы и моменты — тоже интегральные величины: F=∫VfdV (объёмная плотность силы) или поверхностный интеграл тензора натяжений по поверхности, охватывающей тело. Им посвящён следующий урок.
Сводная таблица
Величина
Интеграл
Область
Ориентация
Единица
Допущение независимости
Ψэ, Q
∫D⋅ndS
замкнутая поверхность
внешняя нормаль
Кл
нет других зарядов внутри
Φ
∫B⋅ndS
поверхность на контуре
правый винт с обходом
Вб
∇⋅B=0 (всегда)
I
∫J⋅ndS
сечение проводника
по выбору
А
∇⋅J=0
U
∫E⋅dl
путь A→B
от A к B
В
∇×E=0 (статика)
C, L, M, R
отношения
—
согласование знаков M
Ф, Гн, Ом
линейность
W
21∫E⋅DdV, 21∫B⋅HdV
всё пространство
—
Дж
линейность
P
∫J⋅EdV
проводник
—
Вт
—
P (тепло)
∫q⋅ndS
поверхность
нормаль
Вт
стационарность для баланса
Плоские и осесимметричные модели
В плоской (двумерной) модели поле не зависит от координаты z, и все объёмные и поверхностные интегралы вырождаются: ёмкость, индуктивность, сопротивление, энергия, потери получаются на единицу длины — C′ [Ф/м], L′ [Гн/м], R′ [Ом/м], Wℓ [Дж/м], P′ [Вт/м] (штрих у энергии не используем: W′ в курсе — коэнергия). Полные значения — умножением на глубину модели d: C=C′d, L=L′d, W=Wℓd, P=P′d. С сопротивлением нужна осторожность: если ток течёт вдоль оси z (плоская магнитная модель шины), модель даёт R′ в Ом/м и R=R′d; если же ток течёт в плоскости модели (плоская задача проводимости), результат имеет вид R′ в Ом·м и R=R′/d — проводник тем шире, чем больше глубина. В осесимметричной модели элемент объёма dV=2πrdrdz, элемент поверхности dS=2πrdl, и множитель 2πr появляется под знаком каждого интеграла. Подробно — в уроке о размерности моделей.
Разобранные примеры
Погонная индуктивность коаксиальной линии двумя способами
Внутренний проводник радиуса a, внешний — радиуса b, ток I течёт по внутреннему и возвращается по внешнему; между ними немагнитная среда, μ=μ0. По закону полного тока при a<r<b: H=I/(2πr), B=μ0I/(2πr).
Через поток. Контур тока — это проводник и обратный проводник; натянем на него радиальную полоску длиной 1 м от r=a до r=b. Поток через неё
Φ′=∫ab2πrμ0Idr=2πμ0Ilnab,L′=IΦ′=2πμ0lnab.Через энергию. На единицу длины Wℓ=∫ab2μ0B22πrdr=4πμ0I2lnab, и L′=2Wℓ/I2=2πμ0lnab — тот же результат. Оба способа игнорируют поле внутри проводников (внутреннюю индуктивность); при постоянном токе она добавляет μ0/(8π) от внутреннего провода, при сильном скин-эффекте исчезает.
Числа: a=1 мм, b=3,5 мм, ln3,5=1,253: L′=2⋅10−7⋅1,253=2,51⋅10−7 Гн/м =0,25 мкГн/м.
Погонная ёмкость и проверка связи с индуктивностью
Та же линия с диэлектриком ε=εrε0. Заряд τ на единицу длины внутреннего проводника по теореме Гаусса даёт D=τ/(2πr), E=τ/(2πεr), напряжение U=∫abEdr=2πετlnab, откуда
C′=Uτ=ln(b/a)2πε.
Перемножим с L′: L′C′=μ0ε — логарифмы сократились. Это общее свойство линий с однородным заполнением: 1/L′C′=1/με — скорость распространения волны. Числа при εr=2,3 (полиэтилен): C′=2π⋅8,854⋅10−12⋅2,3/1,253=1,02⋅10−10 Ф/м =102 пФ/м; L′C′=2,51⋅10−7⋅1,02⋅10−10=2,56⋅10−17 с²/м², и μ0ε0εr=1,2566⋅10−6⋅8,854⋅10−12⋅2,3=2,56⋅10−17 с²/м². Скорость 1/L′C′=1,98⋅108 м/с =c/2,3.
Взаимная индуктивность двух соосных витков
Виток 1 радиуса a с током I1 создаёт на оси на расстоянии z индукцию Bz=2(a2+z2)3/2μ0I1a2. Если виток 2 радиуса b≪a расположен соосно на этом расстоянии, поле по его площади почти однородно и Ψ21≈Bzπb2, откуда
M=I1Ψ21=2(a2+z2)3/2μ0πa2b2.
Вычислять M «наоборот» — поток малого витка через большой — было бы гораздо труднее, но взаимность гарантирует тот же результат. Числа: a=10 см, b=1 см, z=0 (витки в одной плоскости): M=2⋅10−31,2566⋅10−6⋅π⋅10−2⋅10−4=1,97⋅10−9 Гн ≈2 нГн. Точное выражение через эллиптические интегралы приведено в [Калантаров]; приближение малого витка годится при b≲0,2a.
Распространённые заблуждения
Поток через «какую-то» поверхность. Поток Φ определён для контура; поверхность — вспомогательна. Но для Ψэ и теплового потока поверхность существенна, а для заряда она должна быть замкнутой.
Ψ=NΦ всегда. Нет: только при одинаковых потоках всех витков. В катушке с рассеянием крайние витки сцеплены с меньшим потоком.
Напряжение как разность потенциалов в переменном поле.∫E⋅dl зависит от пути, если контур из двух путей охватывает переменный поток. Теория цепей обходит это, помещая весь поток внутрь индуктивностей.
L=Ψ/I в насыщенном сердечнике. Отношение зависит от I и не равно 2W/I2; нужно уточнять, о какой индуктивности речь.
Забытая глубина модели. Плоская модель выдаёт погонные величины; C, L, W умножаются на d, а сопротивление — умножается или делится в зависимости от того, течёт ток вдоль оси z или в плоскости. Пропущенный множитель даёт ошибку на порядки.
Энергия по части пространства.W=21∫B⋅HdV — по всему полю; обрезание области расчёта искажает L, особенно для открытых катушек без сердечника.
Вопросы и упражнения
Коаксиальная линия: a=2 мм, b=6 мм, воздух. Найдите L′ и C′ и проверьте, что L′C′=μ0ε0.
Для той же линии длиной 5 м найдите полные L и C. Как изменится C, если заменить воздух на εr=4? А L?
Докажите, что для линейной системы Ψ/I=2W/I2, исходя из W=∫0ΨidΨ при Ψ=Li. Почему доказательство не проходит при L=L(i)?
В двухпроводной системе измерены коэффициенты c11, c22, c12. Выразите через них ёмкость C между проводниками при условии, что их заряды равны и противоположны (Q1=−Q2).
Два соосных витка радиусов 10 см и 1 см раздвинуты по оси на z=10 см. Во сколько раз уменьшится M по сравнению со случаем z=0?
Плоская магнитная модель медной шины сечением 10×100 мм² (ток вдоль оси z, ρэ=1,72⋅10−8 Ом·м) дала R′=17,2 мкОм/м и L′=0,3 мкГн/м. Найдите R и L шины длиной 4 м. Какое из значений вы ожидаете более надёжным и почему (вспомните о концевых эффектах и о замыкании тока)?
4. Из Q1=c11φ1+c12φ2=−Q2=−(c12φ1+c22φ2) выражаем φ2 через φ1, затем C=Q1/(φ1−φ2)=c11+c22+2c12c11c22−c122.
5. Множитель (a2+z2)3/2 вырастает от a3 до (2a2)3/2=22a3: M уменьшится в 2,83 раза.
Итоги
Каждая интегральная величина определяется четырьмя уточнениями: интегрируемое поле, область, ориентация, допущение независимости от выбора области.
Заряд и поток смещения — по замкнутой поверхности с внешней нормалью; магнитный поток — по поверхности, натянутой на контур, с нормалью по правому винту; ток — по сечению; напряжение — по пути, и в переменном поле оно зависит от пути.
C=Q/U, L=Ψ/I=2W/I2, M12=M21, R=U/I=P/I2, Z=U˙/I˙ — отношения, не зависящие от возбуждения только в линейных системах; для многопроводных систем — симметричные матрицы.
Энергия и мощность — объёмные интегралы по всему полю; с амплитудными фазорами средняя мощность содержит множитель 21.
Плоские модели дают погонные величины (C, L, W умножаются на глубину; R — умножается при токе вдоль оси z и делится при токе в плоскости), осесимметричные — множитель 2πr.
Для катушки нужен ответ о потокосцеплении, для изоляции — о напряжённости, для электромагнита — о силе. Карта поля помогает объяснить результат, но не заменяет значение выбранной физической величины. До запуска решателя следует указать контрольную точку, линию или область и физическую единицу. Тогда сеточная проверка будет оценивать именно полезную величину. Показатель «максимальная B» без координаты, способа сглаживания и материала не является достаточным описанием: максимум может относиться к искусственному острому углу.
Точечное поле и сглаживание
На линейном треугольном элементе потенциал интерполируется линейно, а его производные кусочно постоянны. Показанное B может быть дополнительно сглажено для удобства просмотра. Сглаживание не добавляет новых физических данных и иногда скрывает скачок около границы материала. Для количественного сравнения используется один режим обработки и фиксированная координата. В точке на самой границе следует уточнять, с какой стороны берётся значение; у идеального острого угла локальное поле может иметь сингулярный характер.
Поток через поверхность
В плоскопараллельной модели отрезок сечения, вытянутый вдоль Depth, образует поверхность. Поток через неё равен интегралу нормальной B по линии, умноженному на глубину. Знак задаётся нормалью, связанной с направлением обхода. Если окно обработки результатов уже учитывает Depth, повторно умножать результат не нужно. В осесимметрии поверхность образуется вращением линии и содержит вес 2\pi r. Интеграл через стержень измеряет поток этого стержня, но не автоматически потокосцепление всех витков.
Для плоской задачи d в м, \mathrm dl в м, B в Тл, итог \Phi в Вб; нормаль n единичная.
Потокосцепление и результаты цепи
Для тонких замкнутых витков \Psi=\sum_k\int_{S_k}\mathbf B\cdot\mathbf n_k\,\mathrm dS=\sum_k\oint_{C_k}\mathbf A\cdot\mathrm d\mathbf l. Контуры ориентированы по положительному току, а нормали — правилом правой руки. Это сумма интегралов отдельных витков; один произвольный контур по сердечнику не заменяет геометрию всей обмотки.
Потокосцепление суммирует потоки, сцепляющиеся с отдельными витками. В распределённой обмотке витки могут сцепляться с разными потоками; для неё полезен расчёт через потенциал и определённую схему намотки. В FEMM значение \Psi получают в Circuit Properties. Оно зависит от витков и их знаков, тогда как локальное поле при отключённой измерительной обмотке может быть неизменным. Поэтому ошибку \Psi часто нужно искать в соединении обмоток, а не в сетке железа.
\Psi=\sum_{q=1}^{N}\Phi_q=\sum_{q=1}^{N}\oint_{C_q}\mathbf A\cdot\mathrm d\mathbf l
Для тонких замкнутых витков \Phi_q — поток через поверхность с границей C_q; ориентация согласована с положительным током. Единица потокосцепления — Вб.
Интегралы по области
Для энергии выделяется вся область, где хранится учитываемое поле, включая воздух. Выделение только железа пропускает энергию рассеяния и зазора. Для силы по взвешенному тензору выделяется тело, на которое действует сила, а не весь магнитопровод без разбора. Для потерь выбирают материал с соответствующим механизмом и корректной частотой. Прежде чем сравнивать два интеграла, проверьте, относятся ли они к одному телу и одной глубине. Название интеграла не отменяет условий применимости его формулы.
W_{\text{м}} в Дж; w_{\text{м}} в \text{Дж}/\text{м}^3; формула плотности ограничена линейной изотропной обратимой связью B=\mu H.
Минимальная таблица результата
В отчёт включаются параметры источника, частота, глубина или осесимметрия, число элементов и определения измерений. Для B указывают координату и компоненту, для потока — поверхность и знак, для \Psi — цепь и витки, для энергии — выбранные области. Можно добавить карту и профиль B вдоль линии, но они должны пояснять числовые ответы. Такой формат позволяет повторить измерение и отличить изменение физики от изменения визуальной шкалы. Для каждого вопроса укажите нужный тип измерения: точка, контур или интеграл по области.
Ориентация измерения сохраняется в таблице
Когда контур проводится слева направо, его нормаль выбирается правилом ориентации, и поток получает соответствующий знак. Разворот того же контура меняет знак, но не модуль интеграла. Аналогично, смена начала обмотки меняет знак её потокосцепления. Для сравнения двух методов нужно согласовать именно эти ориентации, иначе числа будут противоположными при физически одинаковом поле. Энергия не зависит от направления обхода контура, а знак компоненты силы определяется положительным направлением выбранной оси. В таблице результатов следует указывать не только значения, но и определения величин. Например: поток через поперечное сечение левого стержня, нормаль вверх; потокосцепление цепи первой обмотки с указанным началом; сила на отдельный якорь по оси y.
От локального поля к наблюдаемой величине
Каждая интегральная величина требует собственной поверхности, объёма или электрической цепи. Цветовая шкала индукции не является измерением потокосцепления.
Для потока выберите ориентированную поверхность. В плоской постановке её элемент равен длине участка контура, умноженной на глубину; знак определяется выбранной нормалью.
Для потокосцепления сложите потоки всех витков. У распределённой обмотки разные витки могут охватывать разные потоки, поэтому произведение одного локального потока на число витков требует обоснования.
Для энергии интегрируйте плотность по всем областям поля, включая воздух. Уточните, является результат мгновенной статической энергией или величиной, определённой для гармонической задачи.
Границы применимости. При сравнении двух моделей должны совпадать область интегрирования, число витков, глубина и амплитудное соглашение. Максимум поля у идеального острого угла может не иметь конечного сеточного предела.
Пояснения к официальному руководству
Подробности, которые следует учитывать при настройке и проверке модели. Все учебные пояснения.
Сглаживание поля: что восстанавливается и что остаётся приближением
В плоской задаче линейный потенциал на треугольнике даёт постоянное по элементу производное поле. На общей вершине несколько соседних элементов поэтому предлагают разные значения B. Карта без сглаживания показывает этот кусочно-постоянный результат. Сглаживание создаёт узловое представление производного поля для интерполяции и отображения; исходные узловые значения потенциала и сама сетка от этого не становятся решениями более высокого порядка.
Общий метод восстановления по группе соседних элементов мог бы аппроксимировать значения поля в центрах треугольников плоскостью по наименьшим квадратам. Однако на нерегулярной сетке такая задача может оказаться плохо обусловленной. Для магнитного поля FEMM использует более простое среднее: вклад каждого соседнего элемента взвешивается обратно пропорционально расстоянию от его центра тяжести до рассматриваемого узла. Поэтому большие и маленькие элементы влияют не одинаково, а восстановление не является простым средним всех ближайших значений.
Каждому соседнему треугольнику приписывается вес, обратный расстоянию от его центра до рассматриваемого узла. Сглаживание делает поле непрерывнее на рисунке, но не создаёт новую сетку.
Обозначения и единицы
\mathbf B_e — индукция в центре элемента, Тл; \mathbf x_e и \mathbf x_{\mathrm{node}} — координаты центра и узла, м; w_e — вес, м⁻¹; восстановленная \mathbf B_{\mathrm{node}} имеет единицу Тл.
Условия применения. На интерфейсах действуют специальные правила; нельзя усреднять несовместимые предельные значения через физический разрыв.
В PDF правило задано словами «inversely proportional to the distance»; нормированная формула раскрывает это описание.
Для электростатики процедура другая: в окрестности узла аппроксимируют узловые напряжения плоскостью и используют её градиент для восстановления электрического поля и смещения с учётом материала. Специальные случаи на границах и поверхностях обрабатываются отдельно. Магнитную формулу усреднения B поэтому нельзя выдавать за точное описание любого поля любого модуля.
На границе материалов физически непрерывны определённые компоненты, а не весь вектор. Без поверхностного заряда непрерывна нормальная D; при разных \varepsilon нормальная E обычно разрывна. Аналогично для границы раздела магнитных сред без поверхностного тока непрерывны нормальная B и тангенциальная H. Сглаживание изображения не должно скрывать эти физически обусловленные скачки. Сравнивайте значения с одной выбранной стороны границы раздела и фиксируйте режим сглаживания.
У острых углов поле может иметь сингулярность. Сглаживание способно изменить показанный максимум, но не создать конечное физическое значение для идеальной острой кромки. Для оценки изоляции указывают реальный радиус скругления; для интегральных сил и энергии отдельно проверяют сходимость нужного интеграла. Сравнение сглаженного и несглаженного профиля полезно как диагностический опыт, но не является самостоятельной строгой оценкой погрешности.
Дано. Однородная B=0{,}5 \,\text{Тл} нормально пересекает стержень шириной 10 \,\text{мм} и глубиной 20 \,\text{мм}. Все 100 витков сцепляются с этим потоком.
Создайте плоскопараллельную модель в миллиметрах с Depth=20 мм и Frequency=0. Постройте два концентрических квадрата: наружный x,y=\pm 30 \,\text{мм}, внутренний x,y=\pm 20 \,\text{мм}. Область между ними — один линейный магнитопровод с \mu_r=1000 и \sigma=0.
Обмотка 1 охватывает левый стержень: прямоугольники x=-40\ldots -30 и x=-20\ldots -10 \,\text{мм}, в обоих y=-10\ldots 10 \,\text{мм}. Назначьте цепь «Обмотка 1» типа series и соответственно +100 и -100 витков.
Обмотка 2 охватывает правый стержень: прямоугольники x=10\ldots 20 и x=30\ldots 40 \,\text{мм}, y=-10\ldots 10 \,\text{мм}. Для положительного взаимного потокосцепления назначьте цепь «Обмотка 2» и соответственно -50 и +50 витков. Подпишите на схеме выбранный положительный обход потока.
Замкните внешнюю квадратную область x,y=-150\ldots 150 \,\text{мм} и назначьте Prescribed A=0. Метки воздуха необходимы во всех независимых замкнутых воздушных областях, включая окно. Материал обмоток: \mu_r=1, \sigma=0, без дополнительного источника J.
Задайте размер элемента в железе и обмотках не более 2 \,\text{мм}. Сначала возбудите только одну цепь, затем другую; каждый раз сохраняйте отдельную копию задачи и оба потокосцепления из View → Circuit Properties.
При I_{1}=1 \,\text{А},I_{2}=0 снимите B в середине левого стержня. Проведите поперёк этого стержня контур и измерьте нормальный поток; подпишите направление контура.
Из результатов цепи снимите \Psi_{1} и \Psi_{2}. Для энергии выделите всё железо, обмотки и воздух, а затем отдельно только железо; сравните значения.
Постройте профиль B вдоль линии через окно и стержни. Сравните сглаженный и несглаженный вид, не меняя исходную модель.
Ожидаемый результат. Поток стержня и \Psi/N близки при небольшом рассеянии, но не обязаны совпадать точно. Полная энергия не меньше энергии только железа.
Как проверить. Каждая величина имеет единицу и определение области; для полной энергии не пропущен воздух.
Обработка и анализ результатов
Почему точка, профиль, контур и область дают разные результаты?
Извлеките предусмотренные практикумом точечные, контурные и объёмные величины.
Приведите результаты к одинаковым единицам и ориентации; вычислите невязку независимых способов измерения.
С чем сравнить. Используйте интегральный закон и ориентацию контура из теории занятия. Плотность, интеграл по сечению и полный объёмный результат сравнивают после учёта геометрического множителя.
Итог работы. Таблица величин, единиц и областей измерения без смешения плотности и интеграла.
Пошаговая настройка расчётного примера
Откройте отдельную страницу с последовательными действиями, точными параметрами, накопительными файлами этапов и программами с русскими комментариями.
Рисунок 5.1. Визуализация численных данных FEMM 5: Двухобмоточный магнитопровод: первичное возбуждение. Это экспорт результатов поля, а не снимок окна программы.
Двумерная модель активной части. Модель стали с постоянной магнитной проницаемостью не описывает насыщение. Идеализированные обмотки с нулевой проводимостью не учитывают все потери. Сходимость силы и момента требует отдельной проверки.
Ориентиры в интерфейсеРисунок 5.2. Интерфейс FEMM 5: осесимметричная катушка, геометрические объекты, метки и воздушная область. Открыть в полном размереРисунок 5.3. Поле воздушной катушки: осесимметричный магнитостатический эталон FEMM 5. Открыть в полном размереРисунок 5.4. Результаты цепи Coil: заданный ток и рассчитанное потокосцепление. Открыть в полном размереРисунок 5.5. Интегральная энергия всей магнитной системы; выделение включает катушку и окружающий воздух. Открыть в полном размереРисунок 5.6. Измерение B и H в точке (−25; 15) мм внутри левого стержня. Открыть в полном размереРисунок 5.7. Ориентированный контур поперёк стержня и интеграл нормальной индукции. Открыть в полном размереРисунок 5.8. Профиль |B| через стержни и окно: сглаживание выключено. Контур от (−45; 15) до (45; 15) мм, 400 точек. Position — расстояние от начала контура в мм; |B| — магнитная индукция в Тл. Единицы восстановлены по модели и выбранной величине: само окно FEMM их не печатает. Открыть в полном размереРисунок 5.9. Профиль |B| через стержни и окно: сглаживание включено. Контур от (−45; 15) до (45; 15) мм, 400 точек. Position — расстояние от начала контура в мм; |B| — магнитная индукция в Тл. Единицы восстановлены по модели и выбранной величине: само окно FEMM их не печатает. Открыть в полном размере
Самостоятельная работа
Подготовьте таблицу B, \Phi, \Psi, W_{\text{м}} для одной модели с точным описанием всех измерений.
Критерий выполнения: Другой учащийся может воспроизвести измерения по таблице без угадывания выбранных областей.