FEMM 5

Часть II · Теория четырёх физических модулей FEMM → § 11 · Токопроводимость

Занятие 11.1

Постоянный ток: потенциал, закон Ома и сохранение заряда

Прямой медный проводник. Аналитический эталон. b и h — размеры выделенной геометрии между крайними поверхностями по координатам; пропорции сечения сохранены.
Прямой медный проводник1 — Расчётная областьАналитический эталон. b и h — размеры выделенной геометрии между крайними поверхностями по координатам; пропорции сечения сохранены.Открыть SVG целикомРезультат отдельно

Цель и задачи

Цель. Проверить ток и потери прямого проводника по независимым измерениям.

Задачи

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

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

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

Расчётная лабораторная работа · Определение сопротивления, распределения тока и джоулевых потерь

Инженерный вопрос. Согласуются ли ток вывода, интеграл плотности тока и закон Ома?

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

Лабораторная работа:

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

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

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

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

  2. Постоянный ток в прямом проводнике — аналитический эталон. Проверьте закон Ома и потери в последовательных участках при одних токах и размерах.

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

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

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

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

Данные именно этого опыта. Сохраните напряжение, токи выводов, интеграл тока через сечение и объёмные потери.

Оформление наблюдений. Напряжения, токи выводов, интеграл плотности тока через сечение и потери по областям; затем сопротивления и баланс мощности.

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

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

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

Критерий выполнения. I\approx 580\,\mathrm{А}, R\approx 172{,}413793 мкОм, P\approx 58\,\mathrm{Вт}; токи двух выводов согласованы по знаку и балансу.

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

Основание сравнения. Опора — закон Ома и джоулевы потери из занятия. Удельная проводимость, длина и площадь должны соответствовать тому же режиму; эквивалентный контактный слой не описывает автоматически физику реальной поверхности.

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

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

Границы вывода. Эквивалентные слои контакта и идеальные проводники требуют обоснованных свойств. Двумерная активная длина не включает автоматически реальные присоединения и торцевое распределение тока.

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

Уравнения Максвелла и материальные соотношения

Законы Максвелла определяют выбор измерения в FEMM. Пример с электродами показывает, почему сопротивление зависит от всей краевой задачи, включая форму контактов.

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

Получаемые величины. циркуляция H и охваченный ток; баланс потоков B; ток и сопротивление по электродам; проверка режима квазистатики

  1. Задать ориентацию контура и поверхности.
  2. Сопоставить интеграл с источником и локальные производные с материальным законом.
  3. Для одного проводящего тела изменить расположение электродов при прежнем материале.
  4. Сравнить сопротивление и потери; определить, какие члены полной системы нужны выбранному модулю FEMM.

Страницы источника. PDF 1–20 полностью: 1–2 постановка; 3–15 ротор, циркуляция, дивергенция и поток; 16–20 заряд, закон Ома, материальные уравнения и ток смещения.

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

Аналоговые модели электромагнитного поля

Сходство краевых задач позволяет связать магнитную проводимость с проводимостью токового аналога. Коэффициенты масштаба переводят измеряемую величину аналога в параметр исходного устройства.

Лабораторные работы по методу: Аналогии полей: проводимость, дуальная сеть и диффузия

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

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

Страницы источника. PDF 1–18 полностью: 1–6 прямая аналогия; 7–12 обратная и токовводы; 13–15 сеточные аналоги; 16–18 нелинейные аналоги.

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

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

Распределение постоянного тока определяется законом Ома J=\sigma E=-\sigma \nabla V и условием сохранения заряда \nabla\cdot J=0. При заданных условиях на выводах решение для скалярного потенциала позволяет получить плотность тока в проводящих областях.

Соглашения об обозначениях и единицах
  • \(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. (11.1) · Вывод контрольного соотношения
  2. (11.2) · Объект и физический смысл эталона
  3. (11.3) · Объект и физический смысл эталона
  4. (11.4) · Проверка единиц на эталоне курса

Поле стационарного электрического тока

От электростатики к току

В электростатике проводник был областью, где поле отсутствует: заряды, однажды перераспределившись, останавливаются, и напряжённость внутри обращается в нуль. Стоит, однако, подключить к проводнику источник, поддерживающий разность потенциалов, — и картина меняется. Заряды движутся непрерывно, внутри проводника существует электрическое поле, а само движение зарядов сопровождается выделением тепла. Поле, при котором распределение токов и зарядов не меняется во времени, называют стационарным полем электрического тока. Оно описывается тем же потенциалом φ\varphi, что и электростатическое, но подчиняется другому материальному уравнению и другим граничным условиям.

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

Плотность тока и сила тока

Плотность электрического тока

Вектор J\uJ, направленный вдоль движения положительных зарядов и численно равный заряду, проходящему в единицу времени через единицу площади, перпендикулярной направлению движения. Единица: А/м². Если носители с объёмной плотностью заряда ρ\rho движутся со средней скоростью v\vec v, то J=ρv\uJ=\rho\vec v.

Сила тока II через ориентированную поверхность SS с единичной нормалью n\un — это поток вектора J\uJ: I=∫SJ⋅n dS.I=\int_S\uJ\cdot\un\,\dS . Единица — ампер. Ток — скалярная величина, знак которой определяется выбором нормали: ток «в направлении n\un» положителен. Для провода с равномерным распределением тока по сечению AA формула вырождается в I=JAI=JA, но в общем случае J\uJ меняется по сечению, и делить ток на площадь нельзя — именно это и демонстрирует лаборатория ниже.

Локальный закон Ома

В большинстве проводников при умеренных полях плотность тока пропорциональна напряжённости.

Локальный (дифференциальный) закон Ома

J=σE,E=ρэJ,ρэ=1σ.\uJ=\sigma\uE,\qquad \uE=\rho_{\mathrm э}\uJ,\qquad \rho_{\mathrm э}=\frac1\sigma . Здесь σ\sigma — удельная электрическая проводимость (См/м), ρэ\rho_{\mathrm э} — удельное сопротивление (Ом·м). Для изотропной среды σ\sigma — скаляр, и J∥E\uJ\parallel\uE; в анизотропной среде σ\sigma — тензор.

Проводимость металлов огромна по сравнению с проводимостью электролитов и грунтов: медь σ=5,8⋅107\sigma=5{,}8\cdot10^7 См/м (ρэ=1,72⋅10−8\rho_{\mathrm э}=1{,}72\cdot10^{-8} Ом·м), алюминий 3,5⋅1073{,}5\cdot10^7 См/м (2,86⋅10−82{,}86\cdot10^{-8} Ом·м), влажный грунт — порядка 10−210^{-2} См/м, сухой — 10−310^{-3} См/м и менее. Диапазон охватывает десять порядков, поэтому в задачах с несколькими средами ток практически целиком идёт по металлу, а металлические электроды в слабопроводящей среде можно считать эквипотенциальными.

Удельное сопротивление металлов растёт с температурой: ρэ(T)=ρэ0[1+αρ(T−T0)]\rho_{\mathrm э}(T)=\rho_{\mathrm{э}0}[1+\alpha_\rho(T-T_0)], для меди и алюминия αρ≈4⋅10−3\alpha_\rho\approx 4\cdot10^{-3} К⁻¹. Эта зависимость делает задачу о токе и задачу о нагреве связанными — к ней мы вернёмся в связанных задачах.

Уравнение непрерывности

Заряд сохраняется. Если из замкнутой поверхности SS вытекает ток ∮SJ⋅n dS\oint_S\uJ\cdot\un\,\dS, то заряд QQ внутри убывает с той же скоростью: ∮SJ⋅n dS=−dQ/dt\oint_S\uJ\cdot\un\,\dS=-\dd Q/\dd t. Применяя теорему Остроградского–Гаусса и стягивая объём к точке, получаем локальную форму.

Уравнение непрерывности

∇⋅J=−∂ρ∂t.\nabla\cdot \uJ=-\frac{\partial\rho}{\partial t}. В стационарном режиме распределение заряда не меняется, ∂ρ/∂t=0\partial\rho/\partial t=0, и ∇⋅J=0,∮SJ⋅n dS=0.\nabla\cdot \uJ=0,\qquad \oint_S\uJ\cdot\un\,\dS=0 . Линии тока не имеют ни начала, ни конца: они замыкаются через источник. Интегральная форма — первый закон Кирхгофа: сумма токов, втекающих в любой объём (узел), равна нулю.

Из ∇⋅J=0\nabla\cdot \uJ=0 следует простое и важное правило: вдоль трубки тока — трубки, боковая поверхность которой образована линиями тока, — ток через любое поперечное сечение одинаков. Там, где трубка сужается, плотность тока возрастает.

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

В стационарном режиме ∂B/∂t=0\partial\uB/\partial t=0, поэтому ∇×E=0\nabla\times \uE=0 и, как в электростатике, E=−∇φ\uE=-\nabla \varphi. Подставляя в уравнение непрерывности закон Ома, получаем уравнение для потенциала.

Уравнение стационарного поля тока

∇⋅(σ ∇φ)=0.\nabla\cdot (\sigma\,\nabla \varphi)=0 . В однородной среде (σ=const\sigma=\mathrm{const}) оно переходит в уравнение Лапласа Δφ=0\Delta\varphi=0. Проводимость в однородной области выпадает из уравнения: распределение потенциала определяется только геометрией и граничными условиями, а σ\sigma лишь масштабирует ток.

Граничные условия имеют ясный физический смысл.

  • Электроды. Проводимость электрода (металл) на много порядков выше проводимости среды, поэтому электрод эквипотенциален: φ=φk=const\varphi=\varphi_k=\mathrm{const} на его поверхности (условие Дирихле). Напряжение между электродами U=φ1−φ2U=\varphi_1-\varphi_2.
  • Изолированная поверхность. Ток не может выйти в непроводящую среду: Jn=0J_n=0, то есть ∂φ/∂n=0\partial\varphi/\partial n=0 (условие Неймана). Линии тока касаются такой границы, эквипотенциали подходят к ней под прямым углом.
  • Граница двух проводящих сред. Из ∇⋅J=0\nabla\cdot \uJ=0 следует непрерывность нормальной составляющей J1n=J2nJ_{1n}=J_{2n}, из ∇×E=0\nabla\times \uE=0 — непрерывность касательной E1t=E2tE_{1t}=E_{2t}. Отсюда σ1E1n=σ2E2n\sigma_1E_{1n}=\sigma_2E_{2n}: нормальная составляющая E\uE терпит скачок, на границе существует поверхностный заряд. Закон преломления линий тока — tg⁡α1/tg⁡α2=σ1/σ2\operatorname{tg}\alpha_1/\operatorname{tg}\alpha_2=\sigma_1/\sigma_2 [Бессонов].

Решив краевую задачу, находим E=−∇φ\uE=-\nabla \varphi, J=σE\uJ=\sigma\uE, затем ток любого электрода интегрированием J⋅n\uJ\cdot\un по его поверхности.

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

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

R=UI,G=1R=IU.R=\frac UI,\qquad G=\frac1R=\frac IU . Здесь UU — напряжение между электродами, II — полный ток, стекающий с одного электрода на другой. Единицы: Ом и сименс (См). В линейной среде RR не зависит от UU — это свойство системы (материал и геометрия), а не материала.

Из поля сопротивление извлекают несколькими равноправными способами:

  1. через ток: I=∫S1σE⋅n dSI=\int_{S_1}\sigma\uE\cdot\un\,\dS по поверхности электрода, R=U/IR=U/I;
  2. через мощность (см. ниже): R=P/I2=U2/PR=P/I^2=U^2/P;
  3. в одномерном приближении — суммированием сопротивлений тонких слоёв.

Одномерное приближение для проводника переменного сечения

Пусть ток течёт вдоль оси xx стержня, сечение которого A(x)A(x) меняется медленно. Тогда плотность тока почти однородна по сечению, J(x)=I/A(x)J(x)=I/A(x), E(x)=I/(σA)E(x)=I/(\sigma A), и U=∫0ℓE dx=I∫0ℓdxσ(x)A(x),R=∫0ℓdxσ(x)A(x).U=\int_0^\ell E\,\dd x=I\int_0^\ell\frac{\dd x}{\sigma(x)A(x)},\qquad R=\int_0^\ell\frac{\dd x}{\sigma(x)A(x)} . Это сумма сопротивлений элементарных слоёв dR=dx/(σA)\dd R=\dd x/(\sigma A), соединённых последовательно. Для постоянного сечения R=ℓ/(σA)R=\ell/(\sigma A).

Граница применимости. Формула справедлива лишь при малом наклоне образующей (∣dA/dx∣≪A|\dd A/\dd x|\ll\sqrt{A}, то есть при плавном сужении). При резком изменении сечения, ступеньке или изгибе плотность тока перераспределяется по сечению, линии тока искривляются, и реальное сопротивление больше одномерной оценки: ток не успевает «заполнить» широкую часть и обтекает углы.

Потери Джоуля

Работа поля над зарядами в единице объёма за единицу времени равна J⋅E\uJ\cdot\uE. Эта мощность переходит в тепло.

Объёмная плотность потерь Джоуля

pДж=J⋅E=σE2=J2σ,P=∫VpДж dV.p_{\mathrm{Дж}}=\uJ\cdot\uE=\sigma E^2=\frac{J^2}{\sigma},\qquad P=\int_V p_{\mathrm{Дж}}\,\dV . Единицы: Вт/м³ и Вт. В тепловой задаче pДжp_{\mathrm{Дж}} играет роль объёмного источника qVq_V.

Полная мощность равна UIUI

Подставим E=−∇φ\uE=-\nabla \varphi и воспользуемся тождеством ∇⋅(φJ)=φ ∇⋅J+J⋅∇φ\nabla\cdot (\varphi\uJ)=\varphi\,\nabla\cdot \uJ+\uJ\cdot\nabla \varphi, в котором ∇⋅J=0\nabla\cdot \uJ=0: P=−∫VJ⋅∇φ dV=−∫V∇⋅(φJ) dV=−∮Sφ J⋅n dS.P=-\int_V\uJ\cdot\nabla \varphi\,\dV=-\int_V\nabla\cdot (\varphi\uJ)\,\dV=-\oint_S\varphi\,\uJ\cdot\un\,\dS . На изолированных поверхностях Jn=0J_n=0; остаются электроды с постоянными потенциалами φ1\varphi_1, φ2\varphi_2. Ток, выходящий из объёма через электрод 1 (внешняя нормаль), равен −I-I, через электрод 2 — +I+I. Поэтому P=φ1I−φ2I=UI=I2R=U2/RP=\varphi_1I-\varphi_2I=UI=I^2R=U^2/R. Мощность, распределённая по объёму, в сумме совпадает с мощностью, которую «видит» источник.

Аналогия с электростатикой

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

Электростатика Поле тока
E=−∇φ\uE=-\nabla \varphi E=−∇φ\uE=-\nabla \varphi
D=εE\uD=\varepsilon\uE J=σE\uJ=\sigma\uE
∇⋅D=0\nabla\cdot \uD=0 ∇⋅J=0\nabla\cdot \uJ=0
Q=∮D⋅n dSQ=\oint\uD\cdot\un\,\dS I=∮J⋅n dSI=\oint\uJ\cdot\un\,\dS
C=Q/UC=Q/U G=I/UG=I/U

Соотношение RC=ε/σRC=\varepsilon/\sigma

Пусть одна и та же система электродов помещена сначала в однородный диэлектрик ε\varepsilon, затем в однородную проводящую среду σ\sigma. Потенциал φ\varphi в обоих случаях удовлетворяет уравнению Лапласа с одинаковыми граничными условиями, а по теореме единственности — совпадает. Тогда C=∮εE⋅n dSU,G=∮σE⋅n dSU⇒GC=σε,RC=εσ.C=\frac{\oint\varepsilon\uE\cdot\un\,\dS}{U},\qquad G=\frac{\oint\sigma\uE\cdot\un\,\dS}{U} \quad\Rightarrow\quad \frac GC=\frac\sigma\varepsilon,\qquad RC=\frac\varepsilon\sigma . Зная ёмкость системы, сразу получаем её сопротивление, и наоборот. Величина ε/σ\varepsilon/\sigma имеет размерность времени — это время релаксации заряда, важное для гармонических процессов. Для меди ε0/σ=8,854⋅10−12/5,8⋅107≈1,5⋅10−19\varepsilon_0/\sigma=8{,}854\cdot10^{-12}/5{,}8\cdot10^{7}\approx1{,}5\cdot10^{-19} с.

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

Лаборатория

Лаборатория содержит два режима. В первом — стержень из двух последовательно соединённых материалов σ1\sigma_1, σ2\sigma_2 с линейно изменяющимся сечением; расчёт ведётся в одномерном приближении J=I/A(x)J=I/A(x), а в показаниях выводится полуугол конусности — мера применимости этого приближения. Во втором — радиальный ток в цилиндрическом слое, для которого показано точное решение. Задаются ток II, размеры и проводимости (кнопка «Медь → алюминий» задаёт σ1=58\sigma_1=58, σ2=35\sigma_2=35 МСм/м); выводятся RR, GG, UU, PДжP_{\mathrm{Дж}}, максимумы JJ и pp, продольное сечение, окрашенное по потенциалу, с линиями тока и графики J(x)J(x), E(x)E(x), φ(x)\varphi(x), p(x)p(x).

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

  • В сужении плотность тока и напряжённость растут, а плотность потерь — как J2J^2: тонкие участки греются непропорционально сильно.
  • На границе материалов JJ непрерывна, а EE скачком меняется в σ1/σ2\sigma_1/\sigma_2 раз; график φ(x)\varphi(x) имеет излом.
  • Чем круче конусность стержня, тем дальше от истины одномерная оценка — следите за выводимым полууглом конусности: при углах в десятки градусов приближение ненадёжно.
  • В цилиндрическом слое J∝1/rJ\propto1/r, и большая часть падения напряжения приходится на область вблизи внутреннего электрода.

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

Медный провод

Провод сечением A=2,5A=2{,}5 мм² =2,5⋅10−6=2{,}5\cdot10^{-6} м², длиной ℓ=100\ell=100 м, σ=5,8⋅107\sigma=5{,}8\cdot10^7 См/м, ток I=16I=16 А. R=ℓσA=1005,8⋅107⋅2,5⋅10−6=0,690 Ом,J=IA=6,4⋅106 А/м2=6,4 А/мм2,R=\frac{\ell}{\sigma A}=\frac{100}{5{,}8\cdot10^7\cdot2{,}5\cdot10^{-6}}=0{,}690\ \text{Ом},\qquad J=\frac IA=6{,}4\cdot10^6\ \text{А/м}^2=6{,}4\ \text{А/мм}^2, E=Jσ=0,110 В/м,pДж=JE=7,06⋅105 Вт/м3,P=I2R=177 Вт.E=\frac J\sigma=0{,}110\ \text{В/м},\qquad p_{\mathrm{Дж}}=JE=7{,}06\cdot10^5\ \text{Вт/м}^3,\qquad P=I^2R=177\ \text{Вт}. Проверка: объём провода V=Aℓ=2,5⋅10−4V=A\ell=2{,}5\cdot10^{-4} м³, pДжV=7,06⋅105⋅2,5⋅10−4=177p_{\mathrm{Дж}}V=7{,}06\cdot10^5\cdot2{,}5\cdot10^{-4}=177 Вт. Напряжённость внутри медного провода — всего десятые доли вольта на метр, но она отлична от нуля.

Радиальный ток в цилиндрическом слое (точное решение)

Проводящий слой между коаксиальными электродами радиусов a<ba<b и длины ℓ\ell (краевые эффекты у торцов не учитываем). По симметрии J=J(r)er\uJ=J(r)\vec e_r, и из постоянства тока через любую цилиндрическую поверхность J(r)=I2πrℓ,E(r)=I2πσrℓ,U=∫abE dr=I2πσℓln⁡ba,R=ln⁡(b/a)2πσℓ.J(r)=\frac{I}{2\pi r\ell},\qquad E(r)=\frac{I}{2\pi\sigma r\ell},\qquad U=\int_a^bE\,\dd r=\frac{I}{2\pi\sigma\ell}\ln\frac ba,\qquad R=\frac{\ln(b/a)}{2\pi\sigma\ell}. Тот же результат даёт RC=ε/σRC=\varepsilon/\sigma с ёмкостью коаксиального конденсатора C=2πεℓ/ln⁡(b/a)C=2\pi\varepsilon\ell/\ln(b/a). При a=1a=1 см, b=5b=5 см, ℓ=0,2\ell=0{,}2 м, σ=1\sigma=1 См/м (электролит): R=ln⁡5/(2π⋅1⋅0,2)=1,609/1,257=1,28R=\ln5/(2\pi\cdot1\cdot0{,}2)=1{,}609/1{,}257=1{,}28 Ом. Потенциал φ(r)=φa−I2πσℓln⁡ra\varphi(r)=\varphi_a-\dfrac{I}{2\pi\sigma\ell}\ln\dfrac ra меняется по логарифму — именно его показывает второй режим лаборатории.

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

Электрод в виде полусферы радиуса aa заглублён в однородный грунт σ\sigma так, что его плоская грань лежит на поверхности земли. Ток II растекается радиально в нижнее полупространство: J=I/(2πr2)J=I/(2\pi r^2), E=I/(2πσr2)E=I/(2\pi\sigma r^2). Потенциал относительно бесконечно удалённой точки φ(r)=∫r∞E dr=I2πσr,R=φ(a)I=12πσa.\varphi(r)=\int_r^\infty E\,\dd r=\frac{I}{2\pi\sigma r},\qquad R=\frac{\varphi(a)}I=\frac1{2\pi\sigma a}. При σ=0,01\sigma=0{,}01 См/м, a=0,5a=0{,}5 м: R=1/(2π⋅0,01⋅0,5)=31,8R=1/(2\pi\cdot0{,}01\cdot0{,}5)=31{,}8 Ом. Заметим, что сопротивление определяется в основном грунтом вблизи электрода: половина всего падения потенциала приходится на слой a<r<2aa<r<2a.

Шаговое напряжение. При токе замыкания I=100I=100 А потенциал поверхности земли φ(r)=100/(2π⋅0,01 r)=1592/r\varphi(r)=100/(2\pi\cdot0{,}01\,r)=1592/r В. Человек с длиной шага 0,80{,}8 м, стоящий на расстоянии 1 м от центра электрода, оказывается между точками с потенциалами 15921592 В и 1592/1,8=8841592/1{,}8=884 В: шаговое напряжение ≈708\approx708 В. Отсюда требования к выравниванию потенциалов на подстанциях.

Усечённый конус в одномерном приближении

Алюминиевый стержень (σ=3,5⋅107\sigma=3{,}5\cdot10^7 См/м) длиной ℓ=0,2\ell=0{,}2 м, радиус линейно растёт от r1=5r_1=5 мм до r2=10r_2=10 мм. С r(x)=r1+(r2−r1)x/ℓr(x)=r_1+(r_2-r_1)x/\ell R=∫0ℓdxσπr2(x)=ℓσπ(r2−r1)(1r1−1r2)=ℓσπr1r2=0,23,5⋅107⋅π⋅5⋅10−3⋅10−2=3,64⋅10−5 Ом.R=\int_0^\ell\frac{\dd x}{\sigma\pi r^2(x)}=\frac{\ell}{\sigma\pi(r_2-r_1)}\Bigl(\frac1{r_1}-\frac1{r_2}\Bigr)=\frac{\ell}{\sigma\pi r_1r_2} =\frac{0{,}2}{3{,}5\cdot10^7\cdot\pi\cdot5\cdot10^{-3}\cdot10^{-2}}=3{,}64\cdot10^{-5}\ \text{Ом}. Среднее геометрическое радиусов, а не среднее арифметическое. Наклон образующей (r2−r1)/ℓ=0,025(r_2-r_1)/\ell=0{,}025 мал, одномерное приближение оправдано.

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

  • «Ток течёт по пути наименьшего сопротивления». Ток течёт по всем путям одновременно и делится между ними пропорционально проводимостям: по пути с проводимостью в десять раз меньшей идёт в десять раз меньший, но не нулевой ток. Именно поэтому ток заземлителя растекается по всему полупространству, а не по одной «лучшей» линии.
  • «Внутри проводника с током E=0\uE=0». Это верно только для проводника в электростатическом равновесии. При протекании тока E=J/σ≠0\uE=\uJ/\sigma\ne0; в меди при 6,4 А/мм² напряжённость равна 0,11 В/м — малая, но именно она гонит ток.
  • «J=I/AJ=I/A в любом сечении». Деление тока на площадь законно только при однородном распределении. У изгибов, ступенек, отверстий и вблизи малых электродов плотность тока сильно неоднородна, а сопротивление выше одномерной оценки.
  • Путаница между ρэ\rho_{\mathrm э} и RR. Удельное сопротивление — свойство материала (Ом·м), сопротивление — свойство тела (Ом). «Сопротивление меди» — бессмысленное выражение без указания размеров.
  • «Потенциал электрода неоднороден вдоль него». Падение потенциала вдоль металлического электрода в слабопроводящей среде пренебрежимо мало; электрод — эквипотенциаль. Исключение — длинные тонкие электроды, когда сопротивление самого металла сравнимо с сопротивлением среды.

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

  1. Сферический электрод радиуса aa полностью погружён в безграничную среду σ\sigma. Найдите сопротивление «относительно бесконечности» и сравните с полусферическим заземлителем. Почему результаты отличаются ровно вдвое?
  2. Два коаксиальных слоя: σ1\sigma_1 при a<r<ca<r<c, σ2\sigma_2 при c<r<bc<r<b. Запишите RR и укажите, в каком слое выделяется большая доля мощности при σ1≫σ2\sigma_1\gg\sigma_2.
  3. По плоской пластине толщиной hh ток течёт между двумя противоположными кромками. Покажите, что для квадратной пластины R=1/(σh)R=1/(\sigma h) не зависит от размера стороны. Как изменится результат для прямоугольника со сторонами ℓ\ell (вдоль тока) и ww?
  4. В лаборатории задайте одинаковые материалы и резко различные сечения на концах стержня. Как ведёт себя p(x)p(x)? Оцените, во сколько раз плотность потерь на узком конце превышает плотность на широком при отношении радиусов 1:3.
  5. Используя RC=ε/σRC=\varepsilon/\sigma и ёмкость плоского конденсатора, найдите сопротивление утечки изоляционного слоя толщиной 2 мм площадью 0,5 м² при σ=10−12\sigma=10^{-12} См/м.
  6. Докажите, что в однородном проводнике со стационарным током объёмный заряд равен нулю, а заряды могут располагаться только на поверхностях и границах сред.

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

1. Для сферы J=I/(4πr2)J=I/(4\pi r^2), φ(a)=I/(4πσa)\varphi(a)=I/(4\pi\sigma a), R=1/(4πσa)R=1/(4\pi\sigma a) — вдвое меньше, чем у полусферы, потому что ток растекается в два полупространства: это две одинаковые проводимости, соединённые параллельно.

4. p=J2/σp=J^2/\sigma, а J∝1/r2J\propto1/r^2, поэтому отношение плотностей потерь равно (r2/r1)4=81(r_2/r_1)^4=81.

5. C=ε0S/dC=\varepsilon_0S/d (при εr=1\varepsilon_r=1) и R=ε0/(σC)=d/(σS)=2⋅10−3/(10−12⋅0,5)=4⋅109R=\varepsilon_0/(\sigma C)=d/(\sigma S)=2\cdot10^{-3}/(10^{-12}\cdot0{,}5)=4\cdot10^9 Ом =4=4 ГОм. Для диэлектрика с εr>1\varepsilon_r>1 результат тот же — ε\varepsilon сокращается.

Итоги

  • Стационарное поле тока описывается потенциалом: E=−∇φ\uE=-\nabla \varphi, J=σE\uJ=\sigma\uE, ∇⋅J=0\nabla\cdot \uJ=0, откуда ∇⋅(σ∇φ)=0\nabla\cdot (\sigma\nabla \varphi)=0; в однородной среде — уравнение Лапласа.
  • Граничные условия: электроды эквипотенциальны (φ=const\varphi=\mathrm{const}), изолированные поверхности непроницаемы для тока (Jn=0J_n=0); на границе сред непрерывны JnJ_n и EtE_t.
  • Сопротивление R=U/IR=U/I — характеристика системы; из поля его находят через ток, через мощность P=∫J⋅E dV=I2RP=\int\uJ\cdot\uE\,\dV=I^2R или в одномерном приближении R=∫dl/(σA)R=\int\dd l/(\sigma A) при плавном изменении сечения.
  • Формальная аналогия J↔D\uJ\leftrightarrow\uD, σ↔ε\sigma\leftrightarrow\varepsilon даёт RC=ε/σRC=\varepsilon/\sigma для любой геометрии электродов в однородной среде.
  • Точные решения — цилиндрический слой R=ln⁡(b/a)/(2πσℓ)R=\ln(b/a)/(2\pi\sigma\ell) и полусферический заземлитель R=1/(2πσa)R=1/(2\pi\sigma a) — показывают, что сопротивление сосредоточено там, где плотность тока велика.

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

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

Физический и методический смысл

При постоянном токе J=\sigma E=-\sigma \nabla V и \nabla\cdot J=0. Скалярный потенциал определяет распределение токов в проводящих областях при заданных выводах.

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

Ток I по определению равен \int_S\mathbf J\cdot\mathbf n\,\mathrm dS, а напряжение в стационарном потенциальном поле — \int_a^b\mathbf E\cdot\mathrm d\mathbf l. Переходы I=JS и U=El требуют однородных нормальных и касательных компонент соответственно. Формулы R=l/(\sigma S) и P=I^2R в примере относятся к явно заданной однородной геометрии; при неоднородном поле мощность находят из объёмного интеграла \int_\Omega\mathbf J\cdot\mathbf E\,\mathrm dV.

Вывод контрольного соотношения

Выведите сопротивление прямого бруска, полный ток и джоулеву мощность.

R=\frac{L}{\sigma S},\quad E_x=\frac{U}{L},\quad J_x=\sigma E_x,\quad I=J_xS,\quad P=UI
(11.1)
Постоянная проводимость и равномерное сечение; контакты, нагрев и вытеснение тока в этом эталоне отсутствуют.
  1. Площадь S=0{,}001\cdot 0{,}01=10^{-5}\,\mathrm{м^2}; сопротивление равно 172,4 × 10⁻⁶ Ом.

  2. E=1\,\mathrm{В/м}, J=58\cdot 10^{6}\,\mathrm{А/м^2}; полный ток 580 А.

  3. Мощность 58 Вт совпадает с интегралом J\cdot E по объёму. В центре потенциал 0,05 В.

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

Связь с опытом в FEMM

Решите прямой медный проводник и проверьте I по выводу и интегралу J через поперечное сечение.

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

Проверка и границы вывода

I\approx 580\,\mathrm{А}, R\approx 172{,}413793 мкОм, P\approx 58\,\mathrm{Вт}; токи двух выводов согласованы по знаку и балансу.

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

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

Связь с общими определениями

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

Пояснения к официальному руководству

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

Рисунки из руководства к этому занятию (4)
  1. М.2.53 · Диалог «Материалы»
  2. М.2.54 · Диалог «Свойство узла»
  3. М.2.57 · Диалог «Библиотека материалов для протекания тока»
  4. М.2.58 · Диалог «Свойства проводника»
Формулы и таблицы руководства к этому занятию

Токовый материал: единицы проводимости и смысл loss tangent

В Current Flow частота задаёт статический или установившийся гармонический режим. Материал содержит \sigma _{x},\sigma _{y} в См/м, относительные \varepsilon _{x},\varepsilon _{y} и отдельные Loss Tangent по направлениям. Эти \sigma отличаются единицей от магнитного диалога, где используется МСм/м. При f=0 ёмкостная составляющая исчезает.

Менеджер токовых материалов с записями Air и Ceramic и кнопками Add, Delete, Modify.

Рисунок М.2.53. Диалог «Материалы»

Материал проводящего диэлектрика описывается проводимостью, диэлектрической проницаемостью и тангенсом потерь. Имя Ceramic не определяет эти величины; откройте свойства и проверьте данные для нужной температуры и частоты.

Обозначения на рисунке

Add
Добавить
Delete
Удалить
Modify
Изменить
Air
Воздух
Ceramic
Керамика

Увеличить рисунок · Источник: Д. Микер, руководство FEMM 5, рис. 2.53, с. 80 [32]. Подпись и пояснение переведены для курса.

Current Flow Material Library: слева Copper и Aluminum, справа Air и Ceramic текущей модели.

Рисунок М.2.57. Диалог «Библиотека материалов для протекания тока»

Библиотека токового модуля хранит собственные свойства материалов. Перенос названия Copper из другой физической библиотеки не гарантирует правильные единицы \sigma: перед расчётом проверьте численное значение и его температурную применимость.

Обозначения на рисунке

Library Materials
Материалы библиотеки
Model Materials
Материалы текущей модели
Copper
Медь
Aluminum
Алюминий

Увеличить рисунок · Источник: Д. Микер, руководство FEMM 5, рис. 2.57, с. 83 [32]. Подпись и пояснение переведены для курса.

Для диэлектрика с потерями руководство трактует введённую \varepsilon как модуль комплексной проницаемости. Loss Tangent — \tan \delta, безразмерное число, а не угол \delta в градусах. Угол получают как \arctan (\tan \delta ). Если в документации изготовителя указана действительная часть \varepsilon ', модуль нужно предварительно восстановить: |\varepsilon |=\varepsilon '√(1+\tan ^{2}\delta ). При малых потерях разница мала, но переносить такую приближённую замену на большой \tan \delta нельзя.

Комплексная проницаемость диэлектрика с потерями

\widetilde\varepsilon=|\varepsilon|\left(\cos\varphi-j\sin\varphi\right)
(М.2.19)
При временном множителе \exp (j\omega t) пассивное диэлектрическое запаздывание соответствует отрицательной мнимой части \varepsilon. Поле Rel. Permittivity в токовом модуле трактуется как модуль относительной проницаемости. Если известна действительная часть \varepsilon ', требуемый модуль равен \varepsilon '√(1+\tan ^{2}\varphi ), а не просто \varepsilon ' при больших потерях.

Обозначения и единицы

  • \widetilde\varepsilon=\varepsilon^{\prime}-j\varepsilon^{\prime\prime}; для пассивного материала \varepsilon^{\prime\prime}\ge0.
  • \varphi — угол запаздывания; в тригонометрических функциях радианы.
  • Абсолютная \varepsilon измеряется в Ф/м; относительная \varepsilon _{r} безразмерна. Одна и та же фазовая запись применима к обеим при согласованной нормировке.

Условия применения. Линейный диэлектрик на заданной частоте, конвенция \exp (j\omega t). Индекс направления опущен; при анизотропии каждое направление имеет свой модуль и угол.

Тильда введена как признак комплексной величины. Звёздочка не используется: в формулах мощности она означает комплексное сопряжение.

Тангенс угла диэлектрических потерь

\tan\varphi=\frac{\sin\varphi}{\cos\varphi}
(М.2.20)
В поле Loss Tangent вводят безразмерное отношение, а не угол в градусах. Для \varepsilon =\varepsilon '-j\varepsilon '' оно равно \varepsilon ''/\varepsilon '. Если задано число t=\tan \varphi, нужный угол вычисляют как \varphi =\arctan t, а модуль проницаемости восстанавливают с учётом этого же t.

Обозначения и единицы

  • \tan\varphi=\varepsilon^{\prime\prime}/\varepsilon^{\prime} — безразмерный коэффициент потерь.
  • Обозначение \delta в полях \tan \delta x, \tan \delta y означает тот же угол потерь, который в формуле руководства назван \varphi.

Условия применения. Для обычного пассивного диэлектрика с положительной \varepsilon ' выбирают ветвь 0\le \varphi <\pi /2. Проводимость \sigma учитывается отдельно: повторное включение тех же потерь в \sigma и \tan \varphi приводит к двойному счёту.

Для каждого направления комплексная проводимость складывается из обычной \sigma и диэлектрического вклада. Не добавляйте повторно эквивалентные диэлектрические потери в \sigma, если они уже учтены \tan \delta. Такой токовый модуль не включает индуктивное магнитное вытеснение тока; увеличение частоты не заменяет вихретоковую магнитную постановку.

Комплексная эффективная проводимость по двум направлениям

\begin{aligned}\sigma_{x,\mathrm{eff}}&=\sigma_x+j\omega\varepsilon_0|\varepsilon_{r,x}|e^{-j\varphi_x},\\\sigma_{y,\mathrm{eff}}&=\sigma_y+j\omega\varepsilon_0|\varepsilon_{r,y}|e^{-j\varphi_y}.\end{aligned}
(М.2.21)
Каждый коэффициент суммирует проводящий и диэлектрический токи вдоль своей оси. Раскрытие экспоненты показывает положительную добавку \omega \varepsilon _{0}|\varepsilon _{r}|\sin \varphi к действительной проводимости и реактивную часть \omega \varepsilon _{0}|\varepsilon _{r}|\cos \varphi. При нулевой частоте остаётся \sigma; при нулевом угле диэлектрический вклад чисто реактивен.

Обозначения и единицы

  • \sigma_x,\sigma_y,\sigma_{x,\mathrm{eff}},\sigma_{y,\mathrm{eff}} — См/м.
  • |\varepsilon_{r,x}|,|\varepsilon_{r,y}| — безразмерные вводимые модули; \varepsilon_0 — Ф/м; \omega — рад/с.
  • \varphi_i=\arctan(t_i), где t_i — введённый Loss Tangent по направлению i.

Условия применения. Линейная квазиэлектростатическая гармоника с диагональными материальными коэффициентами в координатах модели. Не включает индуктивный скин-эффект; одни и те же диэлектрические потери учитываются только один раз.

Сохранены обе строки исходного номера (2.21). Модули относительных \varepsilon и отдельные углы \varphi _{x}, \varphi _{y} обозначены явно, поскольку интерфейс задаёт их независимо; это устраняет неоднозначность общего \varphi в оригинале.

Boundary поддерживает Fixed Voltage, Mixed BC, Surface Current Density, Periodic и Anti-Periodic. Conductor Properties задаёт либо полный ток вывода, либо его напряжение, а вторую величину вычисляет. Раздельные поверхности с одним именем образуют общий эквипотенциальный вывод с суммарным током; геометрический разрыв сам по себе не означает электрической независимости.

Токовое Nodal Property с выбором заданного напряжения или Point Current Density; поле источника подписано A/m.

Рисунок М.2.54. Диалог «Свойство узла»

В плоской постановке сосредоточенный в сечении источник имеет линейную нормировку вдоль глубины: на рисунке указаны А/м. Это не плотность объёмного тока А/м² и не автоматически полный ток именованного вывода. Выбор источника и выбор фиксированного потенциала взаимоисключающие.

Обозначения на рисунке

Specified Voltage Property
Заданный потенциал
Point Current Density Property
Точечный токовый источник
Specified Voltage
Значение потенциала
Current Density, A/m
Линейная плотность токового источника, А/м

Увеличить рисунок · Источник: Д. Микер, руководство FEMM 5, рис. 2.54, с. 80 [32]. Подпись и пояснение переведены для курса.

Токовый Conductor Properties с выбранным Prescribed Current, активным током в амперах и выключенным напряжением.

Рисунок М.2.58. Диалог «Свойства проводника»

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

Обозначения на рисунке

Prescribed Current
Заданный полный ток
Voltage (V)
Потенциал, В
Current (A)
Полный ток, А

Увеличить рисунок · Источник: Д. Микер, руководство FEMM 5, рис. 2.58, с. 84 [32]. Подпись и пояснение переведены для курса.

Проверьте на своей модели

  • При заданных \varepsilon ' и \tan \delta вычислите величину для поля permittivity; отдельно покажите, что \tan \delta не является углом.
  • Сверьте полный ток через проводник и через контрольное сечение; для AC сравнивайте комплексные числа.

Источник: руководство FEMM 5 целиком [32]. § 2.8 · § 2.8.1 · § 2.8.2

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

Постоянный ток в прямом проводнике — аналитический эталон

Дано. Длина 100 мм; высота 1 мм; глубина 10 мм; \sigma =58\cdot 10^{6}\,\mathrm{См/м}; слева 0,1 В, справа 0 В; боковые стороны изолированы.

  1. Площадь S=0{,}001\cdot 0{,}01=10^{-5}\,\mathrm{м^2}; сопротивление равно 172,4 × 10⁻⁶ Ом.

  2. E=1\,\mathrm{В/м}, J=58\cdot 10^{6}\,\mathrm{А/м^2}; полный ток 580 А.

  3. Мощность 58 Вт совпадает с интегралом J\cdot E по объёму. В центре потенциал 0,05 В.

Ответ: R=172{,}4137931 мкОм; I=580\,\mathrm{А}; P=58\,\mathrm{Вт}; V_{\mathrm{mid}}=0{,}05\,\mathrm{В}.

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

Углублённый разбор

Прямой проводник: полный ток, плотность тока и потери

Объект и физический смысл эталона

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

Из закона сохранения заряда и закона Ома следует линейность потенциала. Полный ток получается интегрированием плотности по поперечному сечению. В плоской постановке FEMM глубина является физическим размером: без неё нельзя сравнивать полный ток и мощность.

Существующий эталон курса уже решён в FEMM 5. Видео № 3 можно просмотреть после его выполнения, чтобы сопоставить физическую постановку. Отдельное управление MATLAB для этого опыта не нужно.

\frac{\mathrm d}{\mathrm dx}\left(\sigma\frac{\mathrm dV}{\mathrm dx}\right)=0\Rightarrow V(x)=U\left(1-\frac{x}{l}\right)
(11.2)
Постоянная проводимость \sigma; потенциал слева U, справа 0.
J=\sigma\frac{U}{l},\quad I=JS,\quad R=\frac{l}{\sigma S},\quad P=UI=I^2R
(11.3)
J — А/м²; S — м²; I — А; R — Ом; P — Вт. Формула мощности относится к постоянному току.

Проверка единиц на эталоне курса

Дано. Медь \sigma =58\cdot 10^{6} См/м; l=100 мм, ширина 1 мм, глубина 10 мм; U=0{,}1 В.

  1. Площадь поперечного сечения равна 10⁻⁵ м²; напряжённость E=0{,}1/0{,}1=1 В/м.

  2. Умножим плотность тока на площадь и вычислим мощность.

    J=58\cdot10^6\ \text{А}/\text{м}^2,\quad I=580\ \text{А},\quad P=58\ \text{Вт}
    (11.4)

Ответ. J=58 МА/м², I=580 А, R\approx 0{,}1724 мОм, P=58 Вт.

Проверка. Локальное значение co_getpointvalues()[1] — плотность тока, а не ток в амперах. Интеграл активной мощности должен совпасть с UI. Большой ток объясняется идеализированными размерами и источником напряжения, а не допустимым нагревом проводника.

Границы модели.
  • Результат иллюстрирует стационарную электропроводность; температура и допустимый длительный ток здесь не определяются.

Первоисточники: ComProgExpert. FEMM Tutorial #03: DC conduction in a rectangular cube [7].

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

Решите прямой медный проводник и проверьте I по выводу и интегралу J через поперечное сечение.

Таблица 11.1. Исходные данные практикума
Исходная величинаЗначениеЕдиница
L100мм
Высота1мм
Глубина10мм
\sigma58 × 10⁶См/м
U0,1В
  1. Создайте Current Flow Problem, frequency=0, millimeters, Planar, Depth=10; прямоугольник (0,0)–(100,1).
  2. Материал: \sigma_x=\sigma_y=58000000\,\mathrm{См/м}; диэлектрическую проницаемость оставьте конечной, при f=0 она не определяет ток.
  3. Назначьте торцевым проводникам fixed voltage 0.1 и 0 В. Верх и низ не соединяйте с источником.
  4. Решите задачу и прочитайте токи проводников, V в точке (50,0.5), Jx и интегральную мощность. Проверьте знак тока и баланс двух выводов.

Ожидаемый результат. R=172{,}4137931 мкОм; I=580\,\mathrm{А}; P=58\,\mathrm{Вт}; V_{\mathrm{mid}}=0{,}05\,\mathrm{В}. Для FEMM сравнение действительно только при совпадении идеализации и граничных условий.

Как проверить. Большая плотность тока здесь нужна для численной проверки. Допустимый нагрев проводника этим стационарным электрическим примером не установлен. Прямое численное сравнение возможно только при совпадении всех исходных данных и постановки.

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

Согласуются ли ток вывода, интеграл плотности тока и закон Ома?

  1. Сохраните напряжение, токи выводов, интеграл тока через сечение и объёмные потери.
  2. Сопоставьте токи выводов, интеграл плотности тока и потери по областям; восстановите сопротивление и баланс мощности.

С чем сравнить. Опора — закон Ома и джоулевы потери из занятия. Удельная проводимость, длина и площадь должны соответствовать тому же режиму; эквивалентный контактный слой не описывает автоматически физику реальной поверхности.

Итог работы. Оценки сопротивления и мощности с проверкой знаков и зарядового баланса.

Расчётные исследования по теме

Изучите физическую постановку, выполните программу во встроенном Python FEMM 5 и проверьте результаты по сохранённым моделям и независимым критериям.

  • Электрическое сопротивление тонкой пластины

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

    Решений FEMM: 3 · Численных проверок: 8

  • Распределение тока в проводнике переменного сечения

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

    Решений FEMM: 3 · Численных проверок: 9

  • Ток утечки коаксиальной линии

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

    Решений FEMM: 3 · Численных проверок: 8

  • Проводимость анизотропной плёнки

    Квадрат 20 × 20 мм, толщина 0.1 мм; короткие электроды у противоположных углов. Главные оси проводимости совпадают с осями модели.

    Решений FEMM: 6 · Численных проверок: 14

  • Сужение медной дорожки печатной платы

    Медная фольга толщиной 35 мкм: два широких участка 10 × 10 мм соединены шейкой длиной 20 мм и шириной 2 мм; полная длина 40 мм.

    Решений FEMM: 3 · Численных проверок: 9

  • Резистивный датчик положения: холостой отсчёт

    Резистивная полоска 100 × 10 мм, толщина 0.1 мм. На концах 0 и 1 В; потенциал без нагрузки измеряется в точках 20, 50 и 80 мм.

    Решений FEMM: 3 · Численных проверок: 10

  • Распределение тока между параллельными полосами

    Две параллельные полосы длиной 100 мм и шириной по 10 мм, глубина 1 мм. Общие торцевые электроды соединяют медную и бронзовую области.

    Решений FEMM: 3 · Численных проверок: 9

  • Четырёхзондовое измерение поверхностного сопротивления

    Круглая плёнка радиуса 200 мм и толщины 0.1 мм. Четыре соосные позиции с шагом 10 мм: токовые контакты x=±15 мм, измерительные точки x=±5 мм. Радиус токовых контактов 0.2 мм.

    Решений FEMM: 3 · Численных проверок: 8

  • Кондуктометрический датчик уровня жидкости

    Коаксиальные электроды радиусами 2 и 10 мм, высота чувствительной части 100 мм. Жидкость занимает 20, 50 или 80 мм снизу. Расчёт осесимметричный: горизонтальная координата r, вертикальная z.

    Решений FEMM: 9 · Численных проверок: 28

  • Определение координат резистивного сенсорного экрана

    Однородный резистивный прямоугольник 100 × 60 мм, толщина 0.1 мм. Два режима питания последовательно измеряют координаты идеальной точки касания (30,42) мм.

    Решений FEMM: 6 · Численных проверок: 18

  • Полусферический заземлитель и потенциал поверхности грунта

    Полусферический электрод радиуса 0.1 м лежит в однородном полупространстве грунта. Внешняя эквипотенциальная полусфера имеет радиус 10 м. Удельное сопротивление грунта 100 Ом м; электрод 100 В, удалённая граница 0 В. Модель осесимметричная.

    Решений FEMM: 3 · Численных проверок: 16

  • Растекание тока в осесимметричном ступенчатом токоподводе

    Два соосных медных цилиндра длиной по 20 мм соединены ступенью: радиусы 5 и 1 мм. Ток идёт вдоль оси, боковые поверхности изолированы, напряжение между торцами 10 мВ.

    Решений FEMM: 3 · Численных проверок: 9

Все темы учебного практикума

Окна FEMM 5

Снимки интерфейса FEMM показывают геометрию, настройки и работу с результатами в FEMM 5. Откройте изображение, чтобы рассмотреть подписи и значения в полном размере.

Рабочая область модели «Электрическое сопротивление тонкой пластины». Нативный препроцессор FEMM 5.
Рабочая область модели «Электрическое сопротивление тонкой пластины». Нативный препроцессор FEMM 5. Открыть в полном размере
Сохранённое решение «Электрическое сопротивление тонкой пластины». Карта модуля потенциала |V|, В, в проводящей плёнке; постоянный ток.
Сохранённое решение «Электрическое сопротивление тонкой пластины». Карта модуля потенциала |V|, В, в проводящей плёнке; постоянный ток. Открыть в полном размере
Рабочая область модели «Распределение тока в проводнике переменного сечения». Нативный препроцессор FEMM 5.
Рабочая область модели «Распределение тока в проводнике переменного сечения». Нативный препроцессор FEMM 5. Открыть в полном размере
Сохранённое решение «Распределение тока в проводнике переменного сечения». Модуль потенциала |V|, В, и эквипотенциали проводника переменного сечения при постоянном токе.
Сохранённое решение «Распределение тока в проводнике переменного сечения». Модуль потенциала |V|, В, и эквипотенциали проводника переменного сечения при постоянном токе. Открыть в полном размере
Рабочая область модели «Ток утечки коаксиальной линии». Нативный препроцессор FEMM 5.
Рабочая область модели «Ток утечки коаксиальной линии». Нативный препроцессор FEMM 5. Открыть в полном размере
Сохранённое решение «Ток утечки коаксиальной линии». Карта модуля потенциала |V|, В, в слабопроводящей коаксиальной изоляции; стационарная утечка.
Сохранённое решение «Ток утечки коаксиальной линии». Карта модуля потенциала |V|, В, в слабопроводящей коаксиальной изоляции; стационарная утечка. Открыть в полном размере
Рабочая область модели «Проводимость анизотропной плёнки». Нативный препроцессор FEMM 5.
Рабочая область модели «Проводимость анизотропной плёнки». Нативный препроцессор FEMM 5. Открыть в полном размере
Сохранённое решение «Проводимость анизотропной плёнки». Модуль потенциала |V|, В, и эквипотенциали ортотропной проводящей плёнки при постоянном токе.
Сохранённое решение «Проводимость анизотропной плёнки». Модуль потенциала |V|, В, и эквипотенциали ортотропной проводящей плёнки при постоянном токе. Открыть в полном размере
Рабочая область модели «Сужение медной дорожки печатной платы». Нативный препроцессор FEMM 5.
Рабочая область модели «Сужение медной дорожки печатной платы». Нативный препроцессор FEMM 5. Открыть в полном размере
Сохранённое решение «Сужение медной дорожки печатной платы». Карта модуля потенциала |V|, В, и эквипотенциали печатной дорожки с узким участком; постоянный ток.
Сохранённое решение «Сужение медной дорожки печатной платы». Карта модуля потенциала |V|, В, и эквипотенциали печатной дорожки с узким участком; постоянный ток. Открыть в полном размере
Рабочая область модели «Резистивный датчик положения: холостой отсчёт». Нативный препроцессор FEMM 5.
Рабочая область модели «Резистивный датчик положения: холостой отсчёт». Нативный препроцессор FEMM 5. Открыть в полном размере
Сохранённое решение «Резистивный датчик положения: холостой отсчёт». Карта модуля потенциала |V|, В, вдоль резистивного датчика положения; стационарное питание 1 В.
Сохранённое решение «Резистивный датчик положения: холостой отсчёт». Карта модуля потенциала |V|, В, вдоль резистивного датчика положения; стационарное питание 1 В. Открыть в полном размере
Рабочая область модели «Распределение тока между параллельными полосами». Нативный препроцессор FEMM 5.
Рабочая область модели «Распределение тока между параллельными полосами». Нативный препроцессор FEMM 5. Открыть в полном размере
Сохранённое решение «Распределение тока между параллельными полосами». Карта модуля потенциала |V|, В, в параллельных медной и бронзовой ветвях. Одинаковый градиент потенциала сочетается с различными плотностями тока из-за разной проводимости материалов.
Сохранённое решение «Распределение тока между параллельными полосами». Карта модуля потенциала |V|, В, в параллельных медной и бронзовой ветвях. Одинаковый градиент потенциала сочетается с различными плотностями тока из-за разной проводимости материалов. Открыть в полном размере
Увеличенная область токовых контактов четырёхзондовой модели: контакты радиусом 0,2 мм расположены при x = ±15 мм. Два потенциальных зонда задаются точками считывания решения, без отдельных проводящих областей. Внешняя граница плёнки за кадром; подписи скрыты.
Увеличенная область токовых контактов четырёхзондовой модели: контакты радиусом 0,2 мм расположены при x = ±15 мм. Два потенциальных зонда задаются точками считывания решения, без отдельных проводящих областей. Внешняя граница плёнки за кадром; подписи скрыты. Открыть в полном размере
Сохранённое решение «Четырёхзондовое измерение поверхностного сопротивления». Модуль потенциала |V|, В, возле токовых контактов четырёхзондового измерения. Увеличенная область плёнки; внешняя граница за кадром.
Сохранённое решение «Четырёхзондовое измерение поверхностного сопротивления». Модуль потенциала |V|, В, возле токовых контактов четырёхзондового измерения. Увеличенная область плёнки; внешняя граница за кадром. Открыть в полном размере
Рабочая область модели «Кондуктометрический датчик уровня жидкости». Нативный препроцессор FEMM 5.
Рабочая область модели «Кондуктометрический датчик уровня жидкости». Нативный препроцессор FEMM 5. Открыть в полном размере
Сохранённое решение «Кондуктометрический датчик уровня жидкости». Модуль потенциала |V|, В, в осесимметричном резистивном уровнемере; стационарное питание 1 В. Эквипотенциали скрыты для читаемости узкого разреза.
Сохранённое решение «Кондуктометрический датчик уровня жидкости». Модуль потенциала |V|, В, в осесимметричном резистивном уровнемере; стационарное питание 1 В. Эквипотенциали скрыты для читаемости узкого разреза. Открыть в полном размере
Рабочая область модели «Определение координат резистивного сенсорного экрана». Нативный препроцессор FEMM 5.
Рабочая область модели «Определение координат резистивного сенсорного экрана». Нативный препроцессор FEMM 5. Открыть в полном размере
Сохранённое решение «Определение координат резистивного сенсорного экрана». Модуль потенциала |V|, В, в одном координатном режиме резистивной сенсорной панели; питание 1 В.
Сохранённое решение «Определение координат резистивного сенсорного экрана». Модуль потенциала |V|, В, в одном координатном режиме резистивной сенсорной панели; питание 1 В. Открыть в полном размере
Увеличенный меридиональный разрез полусферического заземлителя радиусом 100 мм в однородном грунте. Дуги 200, 500 и 1000 мм задают области сетки; удалённая граница 10 м за кадром.
Увеличенный меридиональный разрез полусферического заземлителя радиусом 100 мм в однородном грунте. Дуги 200, 500 и 1000 мм задают области сетки; удалённая граница 10 м за кадром. Открыть в полном размере
Сохранённое решение «Полусферический заземлитель и потенциал поверхности грунта». Модуль потенциала |V|, В, и эквипотенциали грунта при стационарном растекании тока с полусферического электрода. Увеличенная область; удалённая граница за кадром.
Сохранённое решение «Полусферический заземлитель и потенциал поверхности грунта». Модуль потенциала |V|, В, и эквипотенциали грунта при стационарном растекании тока с полусферического электрода. Увеличенная область; удалённая граница за кадром. Открыть в полном размере
Рабочая область модели «Растекание тока в осесимметричном ступенчатом токоподводе». Нативный препроцессор FEMM 5.
Рабочая область модели «Растекание тока в осесимметричном ступенчатом токоподводе». Нативный препроцессор FEMM 5. Открыть в полном размере
Сохранённое решение «Растекание тока в осесимметричном ступенчатом токоподводе». Модуль плотности тока |J|, А/м², в увеличенном переходе от тонкого стержня к широкому токоподводу. Шкала отображения 0–5·10⁷ А/м²; эквипотенциали скрыты. Торцы модели за кадром.
Сохранённое решение «Растекание тока в осесимметричном ступенчатом токоподводе». Модуль плотности тока |J|, А/м², в увеличенном переходе от тонкого стержня к широкому токоподводу. Шкала отображения 0–5·10⁷ А/м²; эквипотенциали скрыты. Торцы модели за кадром. Открыть в полном размере

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

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

Критерий выполнения: I\approx 580\,\mathrm{А}, R\approx 172{,}413793 мкОм, P\approx 58\,\mathrm{Вт}; токи двух выводов согласованы по знаку и балансу.

Литература

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

Пояснения и типичные ошибки
Что означает отсутствие источника \nabla\cdot J=0?

Заряд не накапливается внутри стационарной области; входящий и выходящий токи балансируют.

Какое контрольное число и единица проверяют этот опыт?

R=172{,}4137931 мкОм; I=580\,\mathrm{А}; P=58\,\mathrm{Вт}; V_{\mathrm{mid}}=0{,}05\,\mathrm{В}.

При каких предположениях получена приведённая формула?

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

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

  • Токовый модуль не вычисляет индукционное вытеснение тока только потому, что проводник сделан из меди.
  • Постоянная проводимость и равномерное сечение; контакты, нагрев и вытеснение тока в этом эталоне отсутствуют.