Скин-эффект и эффект близости: как литцендрат решает проблему потерь

Скин-эффект — это вытеснение переменного тока к поверхности проводника, следующее из уравнений Максвелла: плотность тока экспоненциально спадает вглубь с характерной длиной δ = √(2/ωμσ), называемой глубиной скин-слоя. Эффект близости — это дополнительные вихревые токи, наведённые в проводнике соседними проводниками с переменным полем; в многослойной обмотке его вклад растёт как квадрат числа слоёв. Литцендрат подавляет оба эффекта за счёт разделения медного сечения на множество изолированных транспонированных жил.

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

Откуда берётся скин-эффект: пять строк от Максвелла

Внутри проводника с переменным током работают три уравнения из четырёх максвелловских — закон Фарадея, закон Ампера без тока смещения (на частотах до сотен МГц в меди им можно пренебречь) и закон Ома в локальной форме:

Подставим одно в другое, применим тождество ротора и получим для синусоидального режима с частотой ω уравнение Гельмгольца для плотности тока:

Решение для полубесконечного проводника убывает в глубину по экспоненте:

Сама глубина скин-слоя при этом:

Ток экспоненциально падает и одновременно отстаёт по фазе: на каждой глубине δ амплитуда ослабляется в e ≈ 2,72 раза, а фаза отстаёт на 1 радиан. На глубине 5δ от поверхности плотность тока составляет уже менее 1% от поверхностной.

Это фундаментально, и отсюда следует неочевидное: не существует никакой «границы» между проводящей областью и непроводящей. Ток есть везде, просто в глубине он ничтожно мал. Инженерное приближение «ток течёт в слое толщиной δ» — удобное, но неточное.

Точная мощность потерь учитывает экспоненциальный хвост и даёт на единицу площади поверхности:

Эта формула — основа всех последующих расчётов.

Круглый провод: точное бесселево решение

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

Отсюда получается точное отношение AC-сопротивления к DC:

На практике пользуются аппроксимацией через безразмерный параметр Q = d/δ — отношение диаметра жилы к скин-слою:

Именно из этой таблицы следует рабочий критерий выбора жилы в литцендрате: Q ≤ 1, то есть d ≤ δ. Не «до 2δ», как часто пишут в справочниках — это верхняя аварийная граница, а практический оптимум глубже.

Возьмём медный провод диаметром 1 мм на частоте 1 МГц. Скин-слой δ ≈ 66 мкм, Q ≈ 15. Плотность тока на оси провода составляет примерно 0,01% от поверхностной — средняя часть сечения электрически мертва. Провод всё ещё «проводит», но лишь тонкой трубкой у поверхности. Заменим его литцендратом из 100 жил по 0,1 мм (Q = 1,5, F_R ≈ 1,03) — и ток снова течёт по всей площади меди.

Эффект близости: уравнение Доуэлла и почему слои убивают

Скин-эффект — это реакция провода на собственное поле. Эффект близости — реакция на чужое поле, созданное соседями. И вот тут физика становится недоброй.

В 1966 году P. L. Dowell опубликовал работу «Effects of eddy currents in transformer windings» (Proc. IEE, vol. 113), где вывел аналитическое решение для многослойной обмотки из плоской фольги. Идея простая: заменил круглые провода в слое эквивалентной фольгой той же медной площади, решил одномерную задачу, получил формулу для AC-сопротивления m-слойной обмотки:

В модели есть один важный коэффициент, который Доуэлл ввёл для перехода от фольги к круглому проводу — коэффициент заполнения слоя обмотки η (в англоязычной литературе — porosity factor):

Это не то же самое, что привычный из расчёта трансформаторов коэффициент заполнения окна k_зап = S_медь/S_окно — тот двумерный (доля меди по площади всего окна), а η — одномерный, описывает плотность укладки витков по одной координате (вдоль слоя). И не имеет ничего общего с геологическим «коэффициентом пористости» из СНиПов и ГОСТ 25100, несмотря на похожее английское название — там это про объём пор в грунте.

Физический смысл η простой: η → 1 для плотной упаковки витков встык, η → 0 для редкой намотки с большими промежутками. В формулу Доуэлла η входит как множитель √η внутри Δ, и при рыхлой обмотке потери снижаются — магнитное поле «размазывается» по большему объёму и сцепляется с меньшим количеством меди. Это объясняет старый инженерный трюк с прокладочным шнуром в катушках высокодобротных контуров: искусственно разреженная укладка уменьшает эффект близости.

Не пугайтесь гипертригонометрии в самой формуле. Важно вот что: вторая скобка содержит множитель (m² − 1). Потери от эффекта близости растут как квадрат числа слоёв обмотки. Один слой — эффекта близости нет вообще, работает только скин. Два слоя — уже в три раза больше. Пять слоёв — в 24 раза больше, чем у одного.

Чтобы увидеть цифры, не формулы — возьмём типовой пример из справочника Ridley Engineering (2019): обмотка трансформатора из 5 слоёв медного провода 0,3 мм на частоте 100 кГц. Расчёт по Доуэллу даёт средний рост сопротивления 11,6 раз. У самого внутреннего слоя — 27 раз. Увеличим диаметр до 0,6 мм в попытке «добавить меди и снизить потери» — среднее подскакивает до 51 раза, а внутренний слой — до 123. Больше меди, в восемь раз больше потерь. Это контринтуитивный, но строгий результат одномерной электродинамики.

Литцендрат ломает эту катастрофу в корне: заменяя одну толстую жилу на сотню тонких с транспозицией, он переводит обмотку в режим, где d ≪ δ, и первый член в формуле Доуэлла (собственно скин) становится пренебрежимым, а второй (близость, ∝ m²) остаётся, но с маленьким Δ, потому что жила тонкая.

Внутренняя и внешняя близость: почему их надо разделять

В работах Ferreira (1994) и Sullivan—Zhang (COMPEL 2014) эффект близости в литцендрате разделяется на две составляющие, и это разделение принципиально для понимания физики.

• Внешняя близость (external proximity) — поле, созданное соседними витками обмотки, действует на жилу «снаружи пучка». Именно она даёт множитель (m² − 1) в формуле Доуэлла. Борется транспозицией + разреженной намоткой + чередованием (interleaving) первичной и вторичной обмоток в трансформаторе.
• Внутренняя близость (internal proximity) — поле, созданное другими жилами в том же пучке, действует на данную жилу. Она приводит к тому контринтуитивному эффекту: наивное расщепление толстого провода на множество тонких жил может увеличить потери, если d не ушло заметно ниже δ. Борется выбором достаточно тонкой жилы и правильной упаковкой.

Количественно: интегральная P_прокс в литцендратной обмотке — это сумма двух вкладов, внешний растёт как m², внутренний — как n (число жил в пучке). При грамотном проектировании оба члена примерно равны — это и есть оптимум. Сильный перекос в любую сторону означает, что либо слишком мало жил (перегружен внешний вклад), либо слишком много (перегружен внутренний, лак съел сечение меди).

Как это видит добротность: Q(f) и её максимум

Всё сказанное выше — про потери в меди. Но для инженера, проектирующего высокодобротный контур или резонансный преобразователь, важнее интегральный параметр: добротность катушки Q = ωL/R. Это та же самая буква Q, что в безразмерном параметре Q = d/δ, но означает совсем другое — отношение запасённой энергии к потерянной за период. Будем звать её Q_кат, чтобы не путать.

На низких частотах:

Q_кат растёт линейно с частотой — индуктивное сопротивление растёт, а активное держится на уровне DC. Но по мере приближения к области, где Q = d/δ становится порядка единицы, знаменатель начинает раздуваться из-за F_R(f):

Так как F_R растёт по Доуэллу быстрее, чем линейно, после некоторой частоты Q_кат достигает максимума и начинает снижаться. Этот максимум — оптимальная рабочая точка катушки, её и ищут инженеры при проектировании. Для литцендратных катушек индукционных нагревателей максимум Q_кат обычно приходится на 30–100 кГц, для катушек беспроводной зарядки Qi — на 100–200 кГц, для MRI-катушек — на десятки МГц (с использованием специальных конструкций).

Elektrisola публикует прямые измерения Q(f) для своих стандартных литцендратов. Типовая катушка из 120 жил по 0,1 мм на воздушном сердечнике показывает максимум Q_кат ≈ 500 на частоте 200 кГц — в 2–3 раза выше, чем у аналогичной катушки из одножильного медного провода той же DC-эквивалентной площади. Это и есть тот инженерный выигрыш, ради которого литцендрат и делают.

Три частотных диапазона — и где литцендрат реально работает

Современные модели (Wojda & Kazimierczuk, 2012; Geng et al., IET Power Electronics, 2021) делят работу литцендрата на три режима:

Граница между средним и высоким диапазоном — отдельная интересная физика.

Почему литцендрат умирает выше 2 МГц: межжильная ёмкость

Каждая пара соседних жил в литцендрате — это микроконденсатор. Лак толщиной 5–15 мкм с диэлектрической проницаемостью ε_r ≈ 3,5–4 (полиуретан, полиэфиримид) образует плоский зазор. Удельная межжильная ёмкость для типичного литцендрата 0,1 мм получается порядка 1–10 пФ/м на пару соседних жил, а с учётом многосотенного числа пар — десятки нанофарад на метр суммарной ёмкости.

На частоте f реактивное сопротивление этой ёмкости:

На 100 кГц для 10 нФ получаем 160 Ом — в тысячу раз больше омического сопротивления жил, ёмкость не мешает. На 5 МГц для той же ёмкости — всего 3,2 Ома, уже сопоставимо с AC-сопротивлением. Ток, вместо того чтобы течь через лак (который формально изолятор), начинает «переливаться» между жилами через ёмкостную связь. Транспозиция теряет смысл, литцендрат электрически превращается в одну толстую жилу.

К этому добавляются две беды:

• Самоиндуцируемая резонансная частота (SRF) катушки. Паразитная ёмкость обмотки резонирует с индуктивностью на SRF, выше которой компонент перестаёт быть индуктивностью и становится ёмкостью. Для литцендратных катушек SRF обычно 2–10 МГц.
• Ёмкостные коммутационные потери:
  P_C = C · V² · f

  P_C

  ёмкостные потери, Вт

  V

  амплитуда напряжения на ёмкости, В

  C

  паразитная ёмкость обмотки, Ф

  f

  частота переключений, Гц

  На частотах выше 1 МГц и напряжениях 400–800 В эти потери могут превысить выигрыш от снижения омического сопротивления. Парадоксально, но в некоторых флайбек-преобразователях переход с одножильного провода на литцендрат ухудшает общий КПД именно из-за этого.

Практический вывод: рабочий диапазон литцендрата — 10 кГц – 1 МГц. С оговорками — до 2–3 МГц. Выше 5 МГц — уже специальные конструкции с увеличенной толщиной лака или плоские плёночные проводники.

Чем сейчас считают: от Доуэлла к 3D FEA

Уравнение Доуэлла — одномерное приближение. Для большинства трансформаторов его точности достаточно, но расхождение с реальностью может достигать 60% при неудачной геометрии обмотки (Nan & Sullivan, IEEE Trans. Power Electronics, 2009). Это заставило индустрию перейти на численные методы.

Сегодняшний инструментарий выглядит так:

• Метод Ferreira (1994): точное бесселево решение для каждой жилы, дополненное суперналожением магнитных полей соседей. Точнее Доуэлла для редких обмоток, менее точен для плотных.
• Sullivan—Zhang (COMPEL 2014): аналитическая модель с учётом шага скрутки, направлений и многоступенчатости. Реализована в MATLAB, выложена в открытый доступ на сайте Thayer School (Dartmouth).
• 2D FEA с дискретной транспозицией (Roy, Sakly, IEEE 2014): ключевая идея — решают 2D FEA для параллельных жил, затем математически «перемешивают» их в матрице импеданса по заданному протоколу скрутки. Точность в пределах 10% по отношению R_AC/R_DC.
• 3D FEA с литц-моделью в Ansys Maxwell: литцендрат моделируется одним цилиндром с анизотропной проводимостью. Это позволяет считать катушки с тысячами витков без прорисовки каждой жилы.

В 2022 году Lin et al. предложили метод эффективной гомогенизации для transient FEA, который сжимает вычисление 10 000-жильной обмотки до нескольких минут CPU-времени. Для сравнения: прямое моделирование каждой жилы в той же задаче заняло бы несколько дней.

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