FEMM 5

От поля к интегральным величинам

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

Предварительные понятия

Зачем сворачивать поле в число

Решение полевой задачи — это функции E(r)\uE(\ur), B(r)\uB(\ur), J(r)\uJ(\ur), T(r)T(\ur): бесконечно много чисел. Инженеру же нужны несколько: ёмкость конденсатора, индуктивность катушки, сопротивление шины, потери, сила на якорь. Все эти величины — интегралы поля по некоторой области: кривой, поверхности или объёму. Переход «поле → интеграл → параметр» и есть мост между теорией поля и теорией цепей, где те же CC, LL, RR выступают уже как готовые элементы схемы.

Каждая интегральная величина требует четырёх уточнений, без которых её число лишено смысла: что интегрируется, по какой области, с какой ориентацией (нормали или обхода) и при каких допущениях результат не зависит от произвола выбора области. Этот урок — систематический каталог таких уточнений; предыдущие модули дали отдельные определения, здесь они собраны и сопоставлены.

Потоки, заряд, ток

Поток электрического смещения и заряд

Ψэ=∫SD⋅n dS  [Кл],Q=∫Vρ dV  [Кл].\Psi_{\mathrm э}=\int_S\uD\cdot\un\,\dS\ \ [\text{Кл}],\qquad Q=\int_V\rho\,\dV\ \ [\text{Кл}]. Для замкнутой поверхности SS с внешней нормалью теорема Гаусса даёт Ψэ=Q\Psi_{\mathrm э}=Q — свободный заряд внутри. Отсюда способ измерить заряд электрода по полю: охватить электрод любой замкнутой поверхностью, лежащей в диэлектрике, и проинтегрировать нормальную составляющую D\uD. Результат не зависит от выбора поверхности, пока она не захватывает других зарядов.

Магнитный поток и потокосцепление

Φ=∫SB⋅n dS  [Вб],Ψ=∑k=1NΦk  [Вб].\Phi=\int_S\uB\cdot\un\,\dS\ \ [\text{Вб}],\qquad \Psi=\sum_{k=1}^N\Phi_k\ \ [\text{Вб}]. Поверхность SS опирается на контур витка, нормаль согласована с положительным направлением тока в нём по правилу правого винта. Так как ∇⋅B=0\nabla\cdot \uB=0, поток не зависит от выбора натянутой поверхности — только от контура. Потокосцепление обмотки — сумма потоков всех витков; равенство Ψ=NΦ\Psi=N\Phi справедливо лишь при одинаковых потоках витков (см. поток и потокосцепление). Через векторный потенциал Φ=∮LA⋅dl\Phi=\oint_L\uA\cdot\dl.

Ток

I=∫SJ⋅n dS  [А].I=\int_S\uJ\cdot\un\,\dS\ \ [\text{А}]. Поверхность SS пересекает проводник; её ориентация задаёт знак тока. В стационарном режиме (∇⋅J=0\nabla\cdot \uJ=0) ток один и тот же через любое сечение одного проводника. При переменных полях это верно лишь в МКС-приближении, где ток смещения пренебрежим; иначе поверхности, охватывающие заряд, дают разные значения II.

Напряжение

Напряжение

UAB=∫ABE⋅dl  [В]U_{AB}=\int_A^B\uE\cdot\dl\ \ [\text{В}] — линейный интеграл вдоль ориентированного пути от AA к BB. В статическом поле ∇×E=0\nabla\times \uE=0, интеграл не зависит от пути и равен разности потенциалов φA−φB\varphi_A-\varphi_B. В переменном поле ∮E⋅dl=−dΦ/dt≠0\oint\uE\cdot\dl=-\dd\Phi/\dd t\ne0, и напряжение между двумя точками зависит от пути: два пути A→BA\to B отличаются на ЭДС контура, образованного ими. Поэтому «напряжение на зажимах» переменного устройства имеет смысл только вместе с указанием, как проложены измерительные провода, а в теории цепей — при допущении, что вне элементов переменный поток отсутствует.

Ёмкость

Ёмкость двухэлектродной системы

Для двух проводников с зарядами +Q+Q и −Q-Q в линейном диэлектрике C=QU  [Ф],C=\frac{Q}{U}\ \ [\text{Ф}], где UU — напряжение между ними. Линейность среды гарантирует, что отношение не зависит от UU. Через энергию: C=2Wэ/U2=Q2/(2Wэ)C=2W_{\mathrm э}/U^2=Q^2/(2W_{\mathrm э}).

При nn проводниках заряд каждого линейно зависит от всех потенциалов: Qi=∑jcijφjQ_i=\sum_j c_{ij}\varphi_j. Матрица коэффициентов cijc_{ij} симметрична (cij=cjic_{ij}=c_{ji} — следствие взаимности, доказываемой через теорему Грина), диагональные элементы положительны, внедиагональные отрицательны. Эквивалентно, систему описывают частичными ёмкостями между каждой парой проводников и между каждым проводником и землёй — именно они входят в схему замещения многопроводной линии. Единица везде — фарад.

Индуктивности

Собственная индуктивность

Для контура (обмотки) с током II в линейной среде L=ΨI  [Гн],эквивалентноL=2WмI2,L=\frac{\Psi}{I}\ \ [\text{Гн}],\qquad\text{эквивалентно}\qquad L=\frac{2W_{\mathrm м}}{I^2}, где Ψ\Psi — потокосцепление собственного поля, WмW_{\mathrm м} — энергия магнитного поля во всём пространстве. Два определения совпадают только в линейной системе; в нелинейной вводят секущую и дифференциальную индуктивности (см. следующий урок). Для тонких проводов определение через поток требует оговорки, где проходит контур внутри провода; определение через энергию её не требует и автоматически включает внутреннюю индуктивность.

Взаимная индуктивность

Потокосцепление контура 1, создаваемое током контура 2: M12=Ψ12I2  [Гн].M_{12}=\frac{\Psi_{12}}{I_2}\ \ [\text{Гн}]. Теорема взаимности: M12=M21M_{12}=M_{21}. Для тонких контуров в немагнитной среде это видно из формулы Неймана M12=μ04π∮L1∮L2dl1⋅dl2∣r1−r2∣,M_{12}=\frac{\mu_0}{4\pi}\oint_{L_1}\oint_{L_2}\frac{\dl_1\cdot\dl_2}{|\ur_1-\ur_2|}, симметричной относительно перестановки контуров [Калантаров]. Знак MM зависит от выбора положительных направлений обхода в обоих контурах.

Для nn обмоток Ψi=∑jLijIj\Psi_i=\sum_jL_{ij}I_j; матрица индуктивностей симметрична, Lii=LiL_{ii}=L_i, Lij=MijL_{ij}=M_{ij}. Энергия Wм=12∑i,jLijIiIjW_{\mathrm м}=\tfrac12\sum_{i,j}L_{ij}I_iI_j положительна при любых токах, откуда M122≤L1L2M_{12}^2\le L_1L_2 и коэффициент связи k=M12/L1L2≤1k=M_{12}/\sqrt{L_1L_2}\le1.

Сопротивление, мощность, импеданс

Сопротивление

R=UI=PI2  [Ом],G=1R  [См],R=\frac{U}{I}=\frac{P}{I^2}\ \ [\text{Ом}],\qquad G=\frac1R\ \ [\text{См}], где UU — напряжение между электродами проводника, II — ток через любое его сечение, P=∫VJ⋅E dVP=\int_V\uJ\cdot\uE\,\dV — мощность потерь Джоуля. Два определения совпадают при линейном законе Ома J=σE\uJ=\sigma\uE; определение через мощность удобнее в численных моделях, так как объёмный интеграл сходится быстрее линейного.

Мощность потерь

P=∫VJ⋅E dV=∫VσE2 dV  [Вт]P=\int_V\uJ\cdot\uE\,\dV=\int_V\sigma E^2\,\dV\ \ [\text{Вт}] — мгновенная мощность. Для гармонических полей с амплитудными фазорами средняя мощность P=12∫VRe⁡{J˙⋅E˙∗}dV=12∫V∣J˙∣2σdVP=\tfrac12\int_V\Re\{\cplx{\uJ}\cdot\cplx{\uE}^*\}\dV=\tfrac12\int_V\frac{|\cplx{\uJ}|^2}{\sigma}\dV (см. фазоры). Пропуск множителя 12\tfrac12 — ошибка вдвое.

Импеданс

Для гармонического режима отношение комплексных амплитуд напряжения и тока Z=U˙I˙=R+jX  [Ом].Z=\frac{\cplx U}{\cplx I}=R+\ii X\ \ [\text{Ом}]. Активная часть связана со средними потерями: R=2P/∣I˙∣2R=2P/|\cplx I|^2; реактивная — со средней запасённой энергией: для магнитной системы X=ωLX=\omega L, где L=4⟨Wм⟩/∣I˙∣2L=4\langle W_{\mathrm м}\rangle/|\cplx I|^2. Поскольку напряжение в переменном поле зависит от пути, ZZ привязан к определённому расположению зажимов. При скин-эффекте RR и LL зависят от частоты (вихревые токи).

Энергия и тепловой поток

Энергия поля (линейная среда)

Wэ=12∫VE⋅D dV,Wм=12∫VB⋅H dV  [Дж],W_{\mathrm э}=\tfrac12\int_V\uE\cdot\uD\,\dV,\qquad W_{\mathrm м}=\tfrac12\int_V\uB\cdot\uH\,\dV\ \ [\text{Дж}], интеграл — по всему пространству, где поле отлично от нуля. Множитель 12\tfrac12 справедлив лишь для линейных сред; общие выражения ∫H⋅dB\int\uH\cdot\dd\uB рассматриваются в уроке об энергии и коэнергии.

Тепловой поток

P=∫Sq⋅n dS  [Вт]P=\int_S\uq\cdot\un\,\dS\ \ [\text{Вт}] — мощность, уходящая через поверхность SS в направлении нормали. Для замкнутой поверхности в стационарном режиме равен суммарному тепловыделению внутри: ∮q⋅n dS=∫VqV dV\oint\uq\cdot\un\,\dS=\int_Vq_V\,\dV. Это тепловой аналог теоремы Гаусса и основа проверки энергетического баланса.

Силы и моменты — тоже интегральные величины: F=∫Vf dV\vec F=\int_V\vec f\,\dV (объёмная плотность силы) или поверхностный интеграл тензора натяжений по поверхности, охватывающей тело. Им посвящён следующий урок.

Сводная таблица

Величина Интеграл Область Ориентация Единица Допущение независимости
Ψэ\Psi_{\mathrm э}, QQ ∫D⋅n dS\int\uD\cdot\un\,\dS замкнутая поверхность внешняя нормаль Кл нет других зарядов внутри
Φ\Phi ∫B⋅n dS\int\uB\cdot\un\,\dS поверхность на контуре правый винт с обходом Вб ∇⋅B=0\nabla\cdot \uB=0 (всегда)
II ∫J⋅n dS\int\uJ\cdot\un\,\dS сечение проводника по выбору А ∇⋅J=0\nabla\cdot \uJ=0
UU ∫E⋅dl\int\uE\cdot\dl путь A→BA\to B от AA к BB В ∇×E=0\nabla\times \uE=0 (статика)
CC, LL, MM, RR отношения — согласование знаков MM Ф, Гн, Ом линейность
WW 12∫E⋅D dV\tfrac12\int\uE\cdot\uD\,\dV, 12∫B⋅H dV\tfrac12\int\uB\cdot\uH\,\dV всё пространство — Дж линейность
PP ∫J⋅E dV\int\uJ\cdot\uE\,\dV проводник — Вт —
PP (тепло) ∫q⋅n dS\int\uq\cdot\un\,\dS поверхность нормаль Вт стационарность для баланса

Плоские и осесимметричные модели

В плоской (двумерной) модели поле не зависит от координаты zz, и все объёмные и поверхностные интегралы вырождаются: ёмкость, индуктивность, сопротивление, энергия, потери получаются на единицу длины — C′C' [Ф/м], L′L' [Гн/м], R′R' [Ом/м], WℓW_\ell [Дж/м], P′P' [Вт/м] (штрих у энергии не используем: W′W' в курсе — коэнергия). Полные значения — умножением на глубину модели dd: C=C′dC=C'd, L=L′dL=L'd, W=WℓdW=W_\ell d, P=P′dP=P'd. С сопротивлением нужна осторожность: если ток течёт вдоль оси zz (плоская магнитная модель шины), модель даёт R′R' в Ом/м и R=R′dR=R'd; если же ток течёт в плоскости модели (плоская задача проводимости), результат имеет вид R′R' в Ом·м и R=R′/dR=R'/d — проводник тем шире, чем больше глубина. В осесимметричной модели элемент объёма dV=2πr dr dz\dV=2\pi r\,\dd r\,\dd z, элемент поверхности dS=2πr dl\dS=2\pi r\,\dd l, и множитель 2πr2\pi r появляется под знаком каждого интеграла. Подробно — в уроке о размерности моделей.

Разобранные примеры

Погонная индуктивность коаксиальной линии двумя способами

Внутренний проводник радиуса aa, внешний — радиуса bb, ток II течёт по внутреннему и возвращается по внешнему; между ними немагнитная среда, μ=μ0\mu=\mu_0. По закону полного тока при a<r<ba<r<b: H=I/(2πr)H=I/(2\pi r), B=μ0I/(2πr)B=\mu_0I/(2\pi r).

Через поток. Контур тока — это проводник и обратный проводник; натянем на него радиальную полоску длиной 1 м от r=ar=a до r=br=b. Поток через неё Φ′=∫abμ0I2πrdr=μ0I2πln⁡ba,L′=Φ′I=μ02πln⁡ba.\Phi'=\int_a^b\frac{\mu_0I}{2\pi r}\dd r=\frac{\mu_0I}{2\pi}\ln\frac ba,\qquad L'=\frac{\Phi'}{I}=\frac{\mu_0}{2\pi}\ln\frac ba . Через энергию. На единицу длины Wℓ=∫abB22μ02πr dr=μ0I24πln⁡baW_\ell=\int_a^b\frac{B^2}{2\mu_0}2\pi r\,\dd r=\frac{\mu_0I^2}{4\pi}\ln\frac ba, и L′=2Wℓ/I2=μ02πln⁡baL'=2W_\ell/I^2=\frac{\mu_0}{2\pi}\ln\frac ba — тот же результат. Оба способа игнорируют поле внутри проводников (внутреннюю индуктивность); при постоянном токе она добавляет μ0/(8π)\mu_0/(8\pi) от внутреннего провода, при сильном скин-эффекте исчезает.

Числа: a=1a=1 мм, b=3,5b=3{,}5 мм, ln⁡3,5=1,253\ln3{,}5=1{,}253: L′=2⋅10−7⋅1,253=2,51⋅10−7L'=2\cdot10^{-7}\cdot1{,}253=2{,}51\cdot10^{-7} Гн/м =0,25=0{,}25 мкГн/м.

Погонная ёмкость и проверка связи с индуктивностью

Та же линия с диэлектриком ε=εrε0\varepsilon=\varepsilon_r\varepsilon_0. Заряд τ\tau на единицу длины внутреннего проводника по теореме Гаусса даёт D=τ/(2πr)D=\tau/(2\pi r), E=τ/(2πεr)E=\tau/(2\pi\varepsilon r), напряжение U=∫abE dr=τ2πεln⁡baU=\int_a^bE\,\dd r=\frac{\tau}{2\pi\varepsilon}\ln\frac ba, откуда C′=τU=2πεln⁡(b/a).C'=\frac{\tau}{U}=\frac{2\pi\varepsilon}{\ln(b/a)} . Перемножим с L′L': L′C′=μ0εL'C'=\mu_0\varepsilon — логарифмы сократились. Это общее свойство линий с однородным заполнением: 1/L′C′=1/με1/\sqrt{L'C'}=1/\sqrt{\mu\varepsilon} — скорость распространения волны. Числа при εr=2,3\varepsilon_r=2{,}3 (полиэтилен): C′=2π⋅8,854⋅10−12⋅2,3/1,253=1,02⋅10−10C'=2\pi\cdot8{,}854\cdot10^{-12}\cdot2{,}3/1{,}253=1{,}02\cdot10^{-10} Ф/м =102=102 пФ/м; L′C′=2,51⋅10−7⋅1,02⋅10−10=2,56⋅10−17L'C'=2{,}51\cdot10^{-7}\cdot1{,}02\cdot10^{-10}=2{,}56\cdot10^{-17} с²/м², и μ0ε0εr=1,2566⋅10−6⋅8,854⋅10−12⋅2,3=2,56⋅10−17\mu_0\varepsilon_0\varepsilon_r=1{,}2566\cdot10^{-6}\cdot8{,}854\cdot10^{-12}\cdot2{,}3=2{,}56\cdot10^{-17} с²/м². Скорость 1/L′C′=1,98⋅1081/\sqrt{L'C'}=1{,}98\cdot10^8 м/с =c/2,3=c/\sqrt{2{,}3}.

Взаимная индуктивность двух соосных витков

Виток 1 радиуса aa с током I1I_1 создаёт на оси на расстоянии zz индукцию Bz=μ0I1a22(a2+z2)3/2B_z=\dfrac{\mu_0I_1a^2}{2(a^2+z^2)^{3/2}}. Если виток 2 радиуса b≪ab\ll a расположен соосно на этом расстоянии, поле по его площади почти однородно и Ψ21≈Bzπb2\Psi_{21}\approx B_z\pi b^2, откуда M=Ψ21I1=μ0πa2b22(a2+z2)3/2.M=\frac{\Psi_{21}}{I_1}=\frac{\mu_0\pi a^2b^2}{2(a^2+z^2)^{3/2}} . Вычислять MM «наоборот» — поток малого витка через большой — было бы гораздо труднее, но взаимность гарантирует тот же результат. Числа: a=10a=10 см, b=1b=1 см, z=0z=0 (витки в одной плоскости): M=1,2566⋅10−6⋅π⋅10−2⋅10−42⋅10−3=1,97⋅10−9M=\dfrac{1{,}2566\cdot10^{-6}\cdot\pi\cdot10^{-2}\cdot10^{-4}}{2\cdot10^{-3}}=1{,}97\cdot10^{-9} Гн ≈2\approx2 нГн. Точное выражение через эллиптические интегралы приведено в [Калантаров]; приближение малого витка годится при b≲0,2ab\lesssim0{,}2a.

Распространённые заблуждения

  • Поток через «какую-то» поверхность. Поток Φ\Phi определён для контура; поверхность — вспомогательна. Но для Ψэ\Psi_{\mathrm э} и теплового потока поверхность существенна, а для заряда она должна быть замкнутой.
  • Ψ=NΦ\Psi=N\Phi всегда. Нет: только при одинаковых потоках всех витков. В катушке с рассеянием крайние витки сцеплены с меньшим потоком.
  • Напряжение как разность потенциалов в переменном поле. ∫E⋅dl\int\uE\cdot\dl зависит от пути, если контур из двух путей охватывает переменный поток. Теория цепей обходит это, помещая весь поток внутрь индуктивностей.
  • L=Ψ/IL=\Psi/I в насыщенном сердечнике. Отношение зависит от II и не равно 2W/I22W/I^2; нужно уточнять, о какой индуктивности речь.
  • Забытая глубина модели. Плоская модель выдаёт погонные величины; CC, LL, WW умножаются на dd, а сопротивление — умножается или делится в зависимости от того, течёт ток вдоль оси zz или в плоскости. Пропущенный множитель даёт ошибку на порядки.
  • Энергия по части пространства. W=12∫B⋅H dVW=\tfrac12\int\uB\cdot\uH\,\dV — по всему полю; обрезание области расчёта искажает LL, особенно для открытых катушек без сердечника.

Вопросы и упражнения

  1. Коаксиальная линия: a=2a=2 мм, b=6b=6 мм, воздух. Найдите L′L' и C′C' и проверьте, что L′C′=μ0ε0L'C'=\mu_0\varepsilon_0.
  2. Для той же линии длиной 5 м найдите полные LL и CC. Как изменится CC, если заменить воздух на εr=4\varepsilon_r=4? А LL?
  3. Докажите, что для линейной системы Ψ/I=2W/I2\Psi/I=2W/I^2, исходя из W=∫0Ψi dΨW=\int_0^\Psi i\,\dd\Psi при Ψ=Li\Psi=Li. Почему доказательство не проходит при L=L(i)L=L(i)?
  4. В двухпроводной системе измерены коэффициенты c11c_{11}, c22c_{22}, c12c_{12}. Выразите через них ёмкость CC между проводниками при условии, что их заряды равны и противоположны (Q1=−Q2Q_1=-Q_2).
  5. Два соосных витка радиусов 10 см и 1 см раздвинуты по оси на z=10z=10 см. Во сколько раз уменьшится MM по сравнению со случаем z=0z=0?
  6. Плоская магнитная модель медной шины сечением 10×10010\times100 мм² (ток вдоль оси zz, ρэ=1,72⋅10−8\rho_{\mathrm э}=1{,}72\cdot10^{-8} Ом·м) дала R′=17,2R'=17{,}2 мкОм/м и L′=0,3L'=0{,}3 мкГн/м. Найдите RR и LL шины длиной 4 м. Какое из значений вы ожидаете более надёжным и почему (вспомните о концевых эффектах и о замыкании тока)?

Ответы к пунктам 1, 4 и 5

1. ln⁡3=1,099\ln3=1{,}099; L′=2⋅10−7⋅1,099=2,20⋅10−7L'=2\cdot10^{-7}\cdot1{,}099=2{,}20\cdot10^{-7} Гн/м; C′=2π⋅8,854⋅10−12/1,099=5,06⋅10−11C'=2\pi\cdot8{,}854\cdot10^{-12}/1{,}099=5{,}06\cdot10^{-11} Ф/м; L′C′=1,11⋅10−17L'C'=1{,}11\cdot10^{-17} с²/м² =μ0ε0=\mu_0\varepsilon_0 (1,2566⋅10−6⋅8,854⋅10−12=1,113⋅10−171{,}2566\cdot10^{-6}\cdot8{,}854\cdot10^{-12}=1{,}113\cdot10^{-17}).

4. Из Q1=c11φ1+c12φ2=−Q2=−(c12φ1+c22φ2)Q_1=c_{11}\varphi_1+c_{12}\varphi_2=-Q_2=-(c_{12}\varphi_1+c_{22}\varphi_2) выражаем φ2\varphi_2 через φ1\varphi_1, затем C=Q1/(φ1−φ2)=c11c22−c122c11+c22+2c12C=Q_1/(\varphi_1-\varphi_2)=\dfrac{c_{11}c_{22}-c_{12}^2}{c_{11}+c_{22}+2c_{12}}.

5. Множитель (a2+z2)3/2(a^2+z^2)^{3/2} вырастает от a3a^3 до (2a2)3/2=22 a3(2a^2)^{3/2}=2\sqrt2\,a^3: MM уменьшится в 2,832{,}83 раза.

Итоги

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

Применение в FEMM 5