FEMM 5

Магнитный поток, потокосцепление и ЭДС индукции

Вычислять поток как поверхностный интеграл нормальной составляющей B. Определять потокосцепление как сумму потоков отдельных витков. Знать, когда допустимо Ψ = NΦ. Связать изменение потокосцепления с ЭДС при согласованной ориентации

Две обмотки на магнитопроводе: N₁ = 100 · N₂ = 50 · I₁ = 1 А; основные размеры, поле при первичном возбуждении; Ψ₂ = 6,83 мВб.
Двухобмоточный магнитопровод: первичное возбуждение1 — N₁ · 2 — N₂Расчётный образец. Линейная сталь μᵣ = 1000; глубина 20 мм. b и h — наружные размеры сердечника между противоположными поверхностями.Открыть SVG целикомРезультат отдельноИсходная учебная модель
Предварительные понятия

От поля к напряжению на зажимах

Распределённое магнитное поле само по себе не измеряется вольтметром. Прибор на зажимах катушки показывает напряжение, а оно определяется скоростью изменения одной интегральной величины — потокосцепления. Между полем B(r,t)\uB(\ur,t) и показанием вольтметра стоят три шага: поток через поверхность, сумма потоков по виткам, производная по времени. На каждом шаге есть соглашение об ориентации, и каждое из них влияет на знак ответа. Этот урок проходит все три шага с той тщательностью, которой они заслуживают: ошибки в знаке ЭДС и в различении Φ\Phi и Ψ\Psi — одни из самых частых в практике.

Магнитный поток

Магнитный поток

Поток индукции через ориентированную поверхность SS с единичной нормалью n\un: Φ=∫SB⋅n dS.\Phi=\int_S\uB\cdot\un\,\dS . Единица — вебер, 1 Вб=1 Тл⋅м2=1 В⋅с1\ \text{Вб}=1\ \text{Тл}\cdot\text{м}^2=1\ \text{В}\cdot\text{с}. Поток — скаляр со знаком: смена направления нормали меняет знак Φ\Phi.

Поверхность должна быть ориентирована: для поверхности, натянутой на замкнутый контур LL, направление n\un и направление обхода LL связывают правилом правого винта (обход против часовой стрелки, если смотреть с конца n\un). Для однородного поля, перпендикулярного плоской площадке площади SS, Φ=BS\Phi=BS; в общем случае интеграл приходится вычислять — в лаборатории ниже это делается численной квадратурой.

Поток через контур не зависит от натянутой поверхности

Пусть S1S_1 и S2S_2 — две различные поверхности с одной и той же границей LL и нормалями, согласованными с обходом LL. Объединение S1S_1 и S2S_2 — замкнутая поверхность, причём на одной из них согласованная нормаль совпадает с внешней, а на другой — противоположна ей. По теореме Остроградского–Гаусса и ∇⋅B=0\nabla\cdot \uB=0: ∮B⋅nвнеш dS=Φ1−Φ2=∫V∇⋅B dV=0⇒Φ1=Φ2.\oint\uB\cdot\un_{\text{внеш}}\,\dS=\Phi_1-\Phi_2=\int_V\nabla\cdot \uB\,\dV=0\quad\Rightarrow\quad\Phi_1=\Phi_2 . Следовательно, можно говорить о потоке, сцепленном с контуром: он один для всех поверхностей, опирающихся на этот контур. Для электрического смещения в области с зарядами такой независимости нет — там ∇⋅D=ρ≠0\nabla\cdot \uD=\rho\ne0.

То же утверждение получается через векторный потенциал: по теореме Стокса Φ=∫S∇×A⋅n dS=∮LA⋅dl,\Phi=\int_S\nabla\times \uA\cdot\un\,\dS=\oint_L\uA\cdot\dl , и правая часть вообще не содержит поверхности. Для плоской задачи с A=Azez\uA=A_z\vec e_z контур — длинный прямоугольник глубиной dd с продольными сторонами в точках P1P_1 (ток «на нас») и P2P_2 (ток «от нас»): торцевые стороны лежат в плоскости xyxy и не дают вклада (A⊥dl\uA\perp\dl), а продольные дают Φ=d [Az(P1)−Az(P2)],\Phi=d\,[A_z(P_1)-A_z(P_2)] , в согласии с результатом предыдущего урока. Для проводников конечного сечения AzA_z усредняют по сечению каждого проводника. Плоская модель при этом игнорирует лобовые части обмотки — это надо помнить, сравнивая расчёт с измерением.

Потокосцепление обмотки

Обмотка из NN витков — это один длинный проводник, образующий NN почти замкнутых петель. ЭДС на её зажимах складывается из ЭДС отдельных витков, а каждая из них определяется потоком через свой виток. Поэтому естественная величина для обмотки — сумма потоков.

Потокосцепление

Ψ=∑k=1NΦk,\Psi=\sum_{k=1}^N\Phi_k , где Φk\Phi_k — поток, сцепленный с kk-м витком (при одном и том же, согласованном с направлением тока, правиле ориентации для всех витков). Единица — вебер, как и у потока. Через векторный потенциал Ψ=∮весь проводA⋅dl\Psi=\oint_{\text{весь провод}}\uA\cdot\dl — интеграл по всему проводнику обмотки.

Приближение равных потоков

Если все витки сцеплены с одним и тем же потоком Φ\Phi, то Ψ=NΦ.\Psi=N\Phi . Это выполняется тогда и только тогда, когда через каждый виток проходят одни и те же линии индукции: витки одинаковы по форме и размеру, плотно уложены, а поле не имеет составляющих, замыкающихся между витками. Условия нарушаются при:

  • рассеянии — часть линий поля замыкается, пройдя лишь через часть витков (длинная катушка в собственном поле, катушка на сердечнике с зазором);
  • разных размерах витков — конические, многослойные и распределённые обмотки, где Φk\Phi_k зависит от радиуса витка;
  • большом шаге намотки, когда виток уже не является плоским контуром в однородном поле.

В этих случаях вводят эффективное число витков Nэф=Ψ/Φ1N_{\text{эф}}=\Psi/\Phi_1 или считают Ψ\Psi прямым суммированием.

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

N=10N=10 витков, радиусы растут от r1=10r_1=10 мм до r10=19r_{10}=19 мм с шагом 1 мм, однородное поле B0=0,1B_0=0{,}1 Тл вдоль оси. Поток витка Φk=B0πrk2\Phi_k=B_0\pi r_k^2, и Ψ=B0π∑krk2=0,1⋅π⋅2185⋅10−6=6,86⋅10−4 Вб,\Psi=B_0\pi\sum_kr_k^2=0{,}1\cdot\pi\cdot2185\cdot10^{-6}=6{,}86\cdot10^{-4}\ \text{Вб}, где ∑rk2=102+112+⋯+192=2185\sum r_k^2=10^2+11^2+\dots+19^2=2185 мм². Для сравнения: NΦ1=10⋅0,1⋅π⋅10−4=3,14⋅10−4N\Phi_1=10\cdot0{,}1\cdot\pi\cdot10^{-4}=3{,}14\cdot10^{-4} Вб — в 2,2 раза меньше; NΦ10=1,13⋅10−3N\Phi_{10}=1{,}13\cdot10^{-3} Вб — в 1,65 раза больше. Оценка по среднему радиусу 14,5 мм даёт 6,61⋅10−46{,}61\cdot10^{-4} Вб — на 3,8 % ниже точного значения, потому что r2‾>rˉ 2\overline{r^2}>\bar r^{\,2}. Для конической обмотки понятие «поток катушки» неоднозначно, а потокосцепление определено однозначно.

Закон электромагнитной индукции

Закон Фарадея для контура

ЭДС, индуцируемая в замкнутом контуре LL, равна взятой с обратным знаком скорости изменения сцепленного с ним потока: e=∮LE⋅dl=−dΦdt,для обмоткиe=−dΨdt.e=\oint_L\uE\cdot\dl=-\frac{\dd\Phi}{\dd t},\qquad\text{для обмотки}\quad e=-\frac{\dd\Psi}{\dd t}. Соглашение об ориентации. Выбирается положительное направление обхода контура (для обмотки — положительное направление тока ii). Нормаль n\un к поверхности, через которую считается Φ\Phi, берётся по правилу правого винта относительно этого обхода. ЭДС ee положительна, если она стремится создать ток в положительном направлении.

При таком соглашении знак «минус» выражает правило Ленца. Пусть поток растёт, dΦ/dt>0\dd\Phi/\dd t>0. Тогда e<0e<0: ЭДС гонит ток против положительного направления обхода. Такой ток по правилу правого винта создаёт собственное поле, направленное против n\un, то есть против нарастающего внешнего потока. Индуцированный ток препятствует изменению потока, которое его породило — вот содержание знака «минус», а не утверждение, что «ЭДС отрицательна».

Полезно разобраться и с другой записью, встречающейся в теории цепей: напряжение на индуктивности uL=L di/dt=dΨ/dtu_L=L\,\dd i/\dd t=\dd\Psi/\dd t — со знаком плюс. Противоречия нет: uLu_L — это напряжение, которое внешний источник должен приложить к катушке, чтобы преодолеть ЭДС самоиндукции, uL=−eu_L=-e. Знак зависит от того, рассматриваем ли мы катушку как источник (ee) или как приёмник (uLu_L). В поле же ∇×E=−∂B/∂t\nabla\times \uE=-\partial\uB/\partial t — дифференциальная форма того же закона с тем же минусом.

Гармонический режим

ЭДС при синусоидальном потоке

Пусть NN одинаковых витков площади SS пронизываются однородным полем B=B0cos⁡ωtB=B_0\cos\omega t, нормальным к плоскости витков. Тогда Ψ=NSB0cos⁡ωt=Ψmcos⁡ωt\Psi=NSB_0\cos\omega t=\Psi_m\cos\omega t и e=−dΨdt=NSB0 ωsin⁡ωt=ωΨmcos⁡(ωt−90°).e=-\frac{\dd\Psi}{\dd t}=NSB_0\,\omega\sin\omega t=\omega\Psi_m\cos(\omega t-90°). Амплитуда ЭДС Em=ωΨmE_m=\omega\Psi_m; ЭДС отстаёт по фазе от потокосцепления на четверть периода (в комплексной форме E˙=−jωΨ˙\cplx E=-\ii\omega\cplx\Psi). Действующее значение E=Em2=2πf NΦm2=4,44 fNΦm,E=\frac{E_m}{\sqrt2}=\frac{2\pi f\,N\Phi_m}{\sqrt2}=4{,}44\,fN\Phi_m , поскольку 2π/2=4,4432\pi/\sqrt2=4{,}443. Это «формула трансформаторной ЭДС».

Из неё видно, что ЭДС пропорциональна частоте и амплитуде потока, но не самому потоку в данный момент: в моменты, когда Ψ\Psi максимально, e=0e=0, а когда Ψ\Psi проходит через нуль, ∣e∣|e| максимальна. Отсюда же следует связь между напряжением и индукцией в сердечнике трансформатора: при заданном напряжении и числе витков амплитуда потока Φm=E/(4,44fN)\Phi_m=E/(4{,}44fN) фиксирована, и индукция Bm=Φm/SB_m=\Phi_m/S определяется сечением.

Число витков трансформатора

Сердечник с сечением S=25S=25 см² =2,5⋅10−3=2{,}5\cdot10^{-3} м², допустимая амплитуда индукции Bm=1,5B_m=1{,}5 Тл, f=50f=50 Гц, требуемое действующее напряжение обмотки E=230E=230 В. Φm=BmS=3,75⋅10−3 Вб,E1=4,44⋅50⋅3,75⋅10−3=0,833 В/виток,N=2300,833≈276.\Phi_m=B_mS=3{,}75\cdot10^{-3}\ \text{Вб},\qquad E_1=4{,}44\cdot50\cdot3{,}75\cdot10^{-3}=0{,}833\ \text{В/виток},\qquad N=\frac{230}{0{,}833}\approx276 . Величина «вольт на виток» E1=4,44fΦmE_1=4{,}44f\Phi_m — удобный параметр сердечника: все обмотки на нём имеют одинаковую ЭДС на виток (если пренебречь рассеянием).

Виток в поле прямого провода

Прямоугольный контур длиной ℓ=1\ell=1 м лежит в одной плоскости с длинным проводом; ближняя сторона на расстоянии a=1a=1 см, дальняя — b=11b=11 см. Поле B=μ0i/(2πr)B=\mu_0i/(2\pi r), поток Φ=∫abμ0i2πrℓ dr=μ0iℓ2πln⁡ba=2⋅10−7⋅ln⁡11⋅i ℓ=4,80⋅10−7 i  Вб при ℓ=1 м.\Phi=\int_a^b\frac{\mu_0i}{2\pi r}\ell\,\dd r=\frac{\mu_0i\ell}{2\pi}\ln\frac ba=2\cdot10^{-7}\cdot\ln11\cdot i\,\ell=4{,}80\cdot10^{-7}\,i\ \ \text{Вб при }\ell=1\text{ м}. При i=Imcos⁡ωti=I_m\cos\omega t, Im=100I_m=100 А, f=50f=50 Гц: Φm=4,80⋅10−5\Phi_m=4{,}80\cdot10^{-5} Вб, Em=ωΦm=314,2⋅4,80⋅10−5=15,1E_m=\omega\Phi_m=314{,}2\cdot4{,}80\cdot10^{-5}=15{,}1 мВ. Тот же поток получается как ℓ [Az(a)−Az(b)]\ell\,[A_z(a)-A_z(b)] с Az=−(μ0i/2π)ln⁡rA_z=-(\mu_0i/2\pi)\ln r — проверка формулы для плоских задач. Это простейшая модель наводки от силовой шины на соседнюю петлю.

Лаборатория

Лаборатория позволяет пройти всю цепочку от поля к ЭДС. Источник поля — либо однородное B0\uB_0 с задаваемыми модулем и направлением, либо поле кругового витка с током, вычисляемое аналитически через эллиптические интегралы. Поверхность интегрирования — диск, у которого меняются радиус, положение на оси и ориентация нормали; поток через него находится численной квадратурой. Обмотка задаётся числом витков NN, радиусом первого витка r1r_1, шагом, конусностью κ\kappa (радиус растёт от витка к витку: rk=r1[1+κ(k−1)/(N−1)]r_k=r_1[1+\kappa(k-1)/(N-1)]), положением и ориентацией оси. Выводятся BnB_n в центре диска, поток Φ\Phi через диск, потоки отдельных витков Φk\Phi_k, потокосцепление Ψ=∑Φk\Psi=\sum\Phi_k, для сравнения NΦ1N\Phi_1 и относительное отклонение. В гармоническом режиме всё поле умножается на cos⁡ωt\cos\omega t (задаются частота ff и фаза ωt\omega t, есть анимация); выводятся e(t)=−dΨ/dte(t)=-\dd\Psi/\dd t, амплитуда Em=ωΨmE_m=\omega\Psi_m и действующее значение Em/2E_m/\sqrt2, строятся графики Ψ(t)\Psi(t) и e(t)e(t) за период; стрелка на обмотке показывает положительное направление обхода.

Что проверить:

Поток, потокосцепление и индуктивность

Для контура в собственном поле Ψ\Psi пропорционально току (в линейной среде): Ψ=Li\Psi=Li, где LL — собственная индуктивность; для двух контуров Ψ12=Mi2\Psi_{12}=Mi_2 — взаимная индуктивность. Через векторный потенциал L=Ψ/i=1i∮A⋅dlL=\Psi/i=\dfrac1i\oint\uA\cdot\dl; для плоской обмотки с NN витками и глубиной dd при усреднении AzA_z по сечениям прямого и обратного проводников L=Ndi[⟨Az⟩+−⟨Az⟩−].L=\frac{Nd}{i}\bigl[\langle A_z\rangle_+-\langle A_z\rangle_-\bigr]. Систематическое изложение этих величин — в уроке об интегральных величинах, расчётные формулы для многих конфигураций — в справочнике [Калантаров].

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

  • «Поток через катушку». Поток определён для поверхности, а не для тела. Выражение имеет смысл лишь после указания поверхности (или контура, на который она натянута). Для обмотки с рассеянием разные витки сцеплены с разными потоками, и единственная корректная величина — Ψ\Psi.
  • «ЭДС пропорциональна потоку». ЭДС пропорциональна скорости изменения потокосцепления. Сильное постоянное поле не даёт ЭДС вовсе; слабое, но быстро меняющееся — даёт большую. Поэтому Em=ωΨmE_m=\omega\Psi_m растёт с частотой.
  • «Ψ=NΦ\Psi=N\Phi всегда». Это приближение, верное для плотной обмотки одинаковых витков без рассеяния. Для конической катушки из примера ошибка превышает 100 %.
  • Знак ЭДС. Минус в e=−dΨ/dte=-\dd\Psi/\dd t — утверждение о направлении относительно выбранного обхода, а не о том, что ЭДС «отрицательна». Без явного указания положительного направления обхода и нормали знак ee не определён. Напряжение на индуктивности uL=+dΨ/dtu_L=+\dd\Psi/\dd t — та же физика в другом соглашении.
  • «Поток зависит от формы поверхности». Для поля B\uB — нет, в силу ∇⋅B=0\nabla\cdot \uB=0. Это позволяет выбирать поверхность, удобную для расчёта (например, плоскую), но нужно следить, чтобы её граница была в точности нужным контуром.

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

  1. Квадратная рамка со стороной 20 см вращается с частотой 50 об/с в однородном поле B=0,2B=0{,}2 Тл вокруг оси, лежащей в её плоскости перпендикулярно полю. Найдите Ψ(t)\Psi(t) для N=50N=50 витков и амплитуду ЭДС.
  2. Цилиндрическая катушка радиуса rr из NN витков находится в однородном поле B0B_0, направленном под углом θ\theta к её оси. Запишите Ψ\Psi. Почему шаг намотки в этом случае не влияет на результат?
  3. Две концентрические обмотки на одном стержне без рассеяния: N1=100N_1=100, N2=400N_2=400. Поток Φm=2\Phi_m=2 мВб при f=50f=50 Гц. Найдите действующие ЭДС обеих обмоток и их отношение.
  4. Обмотка намотана в два слоя: 20 витков радиуса 10 мм и 20 витков радиуса 12 мм в однородном осевом поле. Найдите Nэф=Ψ/Φ1N_{\text{эф}}=\Psi/\Phi_1, где Φ1\Phi_1 — поток внутреннего витка.
  5. Докажите, что при изменении направления обхода контура на противоположное ЭДС ee меняет знак, а правило Ленца остаётся в силе.
  6. В плоской модели длиной dd две шины сечением w×hw\times h с токами ±i\pm i. Поясните, почему для вычисления Ψ\Psi нужно усреднять AzA_z по сечению шины, а не брать значение в её центре, и в каком случае разница существенна.

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

1. Ψ=NBScos⁡ωt\Psi=NBS\cos\omega t с S=0,04S=0{,}04 м², ω=2π⋅50=314,2\omega=2\pi\cdot50=314{,}2 с⁻¹: Ψm=50⋅0,2⋅0,04=0,4\Psi_m=50\cdot0{,}2\cdot0{,}04=0{,}4 Вб, Em=ωΨm=125,7E_m=\omega\Psi_m=125{,}7 В, e=125,7sin⁡ωte=125{,}7\sin\omega t В.

3. E1=4,44⋅50⋅100⋅2⋅10−3=44,4E_1=4{,}44\cdot50\cdot100\cdot2\cdot10^{-3}=44{,}4 В; E2=4⋅44,4=177,6E_2=4\cdot44{,}4=177{,}6 В; отношение E2/E1=N2/N1=4E_2/E_1=N_2/N_1=4 — коэффициент трансформации.

4. Ψ=B0π(20⋅102+20⋅122)⋅10−6=B0π⋅4880⋅10−6\Psi=B_0\pi(20\cdot10^2+20\cdot12^2)\cdot10^{-6}=B_0\pi\cdot4880\cdot10^{-6}, Φ1=B0π⋅100⋅10−6\Phi_1=B_0\pi\cdot100\cdot10^{-6}, Nэф=48,8N_{\text{эф}}=48{,}8 вместо N=40N=40.

Итоги

Далее — магнитные цепи, насыщение и постоянные магниты: как свести расчёт потока в сердечнике с зазором к закону Ома для магнитной цепи и когда это допустимо.

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