Беспроводная передача энергии через магнитно-связанные индуктивные катушки

в 10:55, , рубрики: беспроводная энергия, воздушный трансформатор, катушка индуктивности, магнитное поле, мощность, схемотехника, физика, Электроника для начинающих

Введение

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

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

Магнитное поле одиночного витка с током

Рассмотрим магнитное поле одиночного витка с током. Найдем магнитное поле витка в любой точке пространства. Почему необходимо подобное рассмотрение? Потому что почти во всех книгах, по крайней мере в тех, которые удалось отыскать автору статьи, решение данной задачи ограничевается лишь нахождением одной компоненты магнитного поля вдоль оси витка $B_z(z)$, в то время как в этой статье мы отыщем закон для магнитного поля во всем пространстве.

Магнитное поле одиночного витка с током
Иллюстрация к закону Био-Савара-Лапласа

Для нахождения магнитного поля, воспользуемся законом Био-Савара-Лапласа (смотри Википедия — Закон Био-Савара-Лапласа). На рисунке видно, что центр системы координат $O$ совпадает с центром витка. Контур окружности витка обозначен как $C$, а радиус окружности — как $a$.По витку течет ток $I$. $vec{r}$ — это переменная-радиус-вектор из начала координат в произвольную точку витка. $vec{r}_0$ — это радиус-вектор в точку наблюдения. Еще нам понадобится полярный угол $varphi$ — угол между радиус-вектором $vec{r}$ и осью $OX$. Расстояние от оси витка до точки наблюдения обозначим за $rho$. И наконец, $mathrm{d}vec{r}$ — элементарное приращение радиус-вектора $vec{r}$.

Согласно закону Био-Савара-Лапласа, элемент контура с током $mathrm{d}vec{r}$ создает элементарный вклад в магнитное поле, который дается формулой

$mathrm{d}vec{B}(vec{r}_0)=frac{mu_0 I}{4pi} cdot frac{[,mathrm{d}vec{r} times (vec{r}_0-vec{r})]}{|vec{r}_0-vec{r}|^3}$

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

$vec{r}_0=(rhocos{varphi}, rhosin{varphi}, z) overset{varphi=0}{rightarrow} (rho, 0, z)$

$vec{r}=(acos{varphi}, asin{varphi}, 0)$

$mathrm{d}vec{r}=(-asin{varphi}, acos{varphi}, 0),mathrm{d}varphi$

$vec{r}_0-vec{r}=(rho -acos{varphi}, -asin{varphi}, z)$

$[mathrm{d}vec{r} times (vec{r}_0-vec{r})]=begin{vmatrix} vec{e}_x& vec{e}_y& vec{e}_z\ -asin{varphi},mathrm{d}varphi& acos{varphi},mathrm{d}varphi& 0\ rho -acos{varphi}& -asin{varphi}& z end{vmatrix}=(azcos{varphi}, azsin{varphi}, a^2 -arhocos{varphi}),mathrm{d}varphi$

$|vec{r}_0-vec{r}|^3=left(rho^2 + a^2 + z^2 -2rho acos{varphi}right)^{frac{3}{2}}$

Для того чтобы найти результирующее магнитное поле, нужно проинтегрировать по всему контуру витка, то есть

$vec{B}(vec{r}_0)=int_C{,mathrm{d}vec{B}(vec{r}_0)}$

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

$B_z(rho, z)=frac{mu_0 I}{4pi} int_0^{2pi}{frac{left( a^2 -rho acos{varphi}right),mathrm{d}varphi}{left(rho^2 + a^2 + z^2 -2rho acos{varphi}right)^{frac{3}{2}}}}$

$B_r(rho, z)=frac{mu_0 I}{4pi} int_0^{2pi}{frac{a,z,,mathrm{d}varphi} {left(rho^2 + a^2 + z^2 -2rho acos{varphi}right)^{frac{3}{2}}}}$

Для нахождения абсолютного значения магнитного поля необходимо просуммировать компоненты по теореме Пифагора $B=sqrt{B_r^2 + B_z^2}$.

Продемонстрируем полученное решение на примере витка радиуса $a=0.1$ (м) и $I=1$ (А).

Амплитуда аксиального магнитного поля
Амплитуда аксиальной компоненты магнитного поля

Амплитуда радиального магнитного поля
Амплитуда радиальной компоненты магнитного поля

Абсолютная амплитуда магнитного поля
Абсолютная амплитуда магнитного поля

Заметим, что для витка произвольной формы, на больших расстояниях $zgg a$, т.е. много больше характерного размера витка, поведение магнитного поля будет стремиться к найденному решению.

Подсказка...

Для подобных вычислений и построения графиков удобно использовать MathCad 15

Катушка индуктивности. Магнитно-связанные катушки

Теперь, когда мы знаем решение для магнитного поля одного витка, можем найти индуктивность катушки, состоящей из $n$ витков. По определению индуктивность — это коэффициент пропорциональности между током в витке и магнитным потоком через площадь сечения витка. Мы пользуемся здесь идеальной моделью катушки, которая безразмерна по направлению своей оси симметрии. Конечно же, на практике такого не бывает. Однако, как приближенные, полученные формулы будут достаточно хороши. Хотя катушки и считаются безразмерными вдоль $OZ$, необходимо задаться ненулевым радиусом сечения провода. Обозначим его $delta$, и пример равным $delta=0,1$ (мм). Иначе при интегрировании магнитного потока подынтегральное выражение обратится в бесконечность.

индуктивно связанные катушки
Индуктивно связанные катушки

На рисунке изображены две магнитно связанные катушки. Пусть первая катушка имеет радиус $a_1$ и содержит $n_1$ витков, а вторая — $a_2$ и $n_2$ соответственно. Тогда для нахождения собственных индуктивностей необходимо вычислить магнитный поток каждой катушки через свое собственное сечение.

$Phi=iint_S{vec{B}cdotvec{,mathrm{d}S}}=int_0^{2pi}{int_0^{a-delta}{B_z(rho,z)rho,mathrm{d}rho,mathrm{d}varphi}}=2piint_0^{a-delta}{B_z(rho,z)rho,mathrm{d}rho}$

Поскольку в катушке много витков, найдем велечину, называемую потокосцепление, дважды умножив на количество витков

$Psi=frac{1}{2}n^2mu_0 I int_0^{a-delta}{int_0^{2pi}{frac{left( a^2 - rho acos{varphi}right),mathrm{d}varphi}{left(rho^2 + a^2 + z^2 -2rho acos{varphi}right)^{frac{3}{2}}}}rho,mathrm{d}rho}$

По определению, индуктивность это коэффициент пропорциональности $L$ в формуле $Psi=LI$. Таким образом, получим собственные индуктивности катушек

$L_1=frac{1}{2}n_1^2mu_0 int_0^{a_1-delta}{int_0^{2pi}{frac{left( a_1^2 - rho a_1cos{varphi}right),mathrm{d}varphi}{left(rho^2 + a_1^2 -2rho a_1cos{varphi}right)^{frac{3}{2}}}}rho,mathrm{d}rho}$

$L_2=frac{1}{2}n_2^2mu_0 int_0^{a_2-delta}{int_0^{2pi}{frac{left( a_2^2 - rho a_2cos{varphi}right),mathrm{d}varphi}{left(rho^2 + a_2^2 -2rho a_2cos{varphi}right)^{frac{3}{2}}}}rho,mathrm{d}rho}$

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

$Psi_{12}=frac{1}{2}n_1 n_2mu_0 I int_0^{a_2-delta}{int_0^{2pi}{frac{left( a_1^2 - rho a_1cos{varphi}right),mathrm{d}varphi}{left(rho^2 + a_1^2 + z^2 -2rho a_1cos{varphi}right)^{frac{3}{2}}}}rho,mathrm{d}rho}$

Тогда взаимная индуктивность катушек дается выражением

$M_{12}=frac{1}{2}n_1 n_2mu_0 int_0^{a_2-delta}{int_0^{2pi}{frac{left( a_1^2 - rho a_1cos{varphi}right),mathrm{d}varphi}{left(rho^2 + a_1^2 + d^2 -2rho a_1cos{varphi}right)^{frac{3}{2}}}}rho,mathrm{d}rho}$

Насколько известно автору, такие интегралы можно взять только численно.
Заметим, что как правило $Psi_{12}=Psi_{21}$ и $M_{12}=M_{21}=M$. Коэффициентом связи катушек называется величина

$k=frac{M}{sqrt{L_1L_2}}$

Исследуем зависимость коэффициента связи катушек от расстояния. Для этого рассмотрим две одинаковые катушки с радиусом витков $a_1=a_2=0.1$ (м) и количеством витков $n_1=n_2=100$. При этом собственная индуктивность каждой из катушек составит $L_1=L_2=8.775$ (мГн).

График зависимости коэффициента связи двух одинаковых катушек от расстояния между ними
Коэффициент связи катушек от расстояния между ними

График не изменится, если одинаково изменить число витков в обеих катушках, либо одинаково изменить радиус обеих катушек. Коэффициент связи удобно выражать в процентах. Из графика видно, что даже при расстоянии между катушками в 1 (мм) коэффицент связи меньше 100%. Коэффициент падает до 10% на расстоянии порядка 60 (мм), и до 1% на 250 (мм).

Беспроводная передача энергии

Итак, нам известны индуктивности и коэффициент связи. Теперь воспользуемся теорией электрических цепей переменного тока для поиска оптимальных параметров, при которых передаваемая мощность оказалась бы максимальной. Для понимания этого параграфа читатель должен быть знаком с понятием электрического импеданса, а также с законами Кирхгоффа и законом Ома. Как известно из теории цепей, две индуктивно-связанные катушки образуют воздушный трансформатор. Для анализа трансформаторов удобна Т-образная схема замещения.

Трансформатор и его схема замещения
Воздушный трансформатор и его схема замещения

Передающую катушку слева будем условно называть «трасмиттер», а принимающую катушку справа — «ресивер». Между катушками коэффициент связи $k$. На стороне ресивера находится потребитель, представленный нагрузкой $z_L$. Нагрузка в общем случае может быть комплексной. Входное напряжение на стороне трансмиттера $u_1$, а входной ток — $i_1$. Напряжение, передаваемое на ресивер — $u_2$, и передаваемый ток $i_2$. Полный импеданс на стороне трансмиттера обозначим как $z_1$, а полный импеданс на стороне ресивера $z_1$.

Предполагается, что на вход схемы подается синусоидальное напряжение $u_1=u_{1m}sin{omega t}$.

Обозначим $R_{coil,1}, R_{coil,2}, L_{coil,1}, L_{coil,2}, M$ — сопротивления и индуктивности катушек (две собственные и одна взаимная) соответственно. Тогда, согласно теории

$z_1=R_{coil,1} + jomega (L_{coil,1} - M)$

$z_2=R_{coil,2} + jomega(L_{coil,2} - M) + R_{load} + j X_{load}$

С другой стороны, согласно нашим обозначениям

$z_1=r_1 + j x_1$

$z_2=r_2 + j x_2$

где $r_1, r_2$ — полные активные сопротивления на стороне трансмиттера и ресивера соответственно, и $x_1, x_2$ — полные реактивные сопротивления.

Импеданс связи равен $z_3=j omega M=j x_3$.

Найдем входной ток цепи

$i_1=frac{u_1}{z_1 + z_2 || z_3}$

где знак $||$ обозначает параллельное соединение сопротивлений. Тогда напряжение, переданное на ресивер

$u_2=u_1 - i_1 z_1=u_1left(1 - frac{z_1}{z_1 + z_2 || z_3}right)$

И наведенный ток

$i_2=frac{u_2}{z_2}=frac{u_1}{z_2}left(1 - frac{z_1}{z_1 + z_2 || z_3}right)$

Можем найти комплексную мощность, переданную в ресивер

$s_2=u_2 i_2^*=p_2 + jq_2$

Таким образом имеем выражение для комплексной мощности

$s_2=|u_1|^2 z_2 left|frac{z_3}{z_1z_2+z_1z_3+z_2z_3}right|^2$

Выражение для активной компоненты мощности

$p_2=|u_1|^2 r_2 left|frac{z_3}{z_1z_2+z_1z_3+z_2z_3}right|^2$

Выражение для реактивной компоненты мощности

$q_2=|u_1|^2 x_2 left|frac{z_3}{z_1z_2+z_1z_3+z_2z_3}right|^2$

В большинстве практических задач требуется передать максимальную активную мощность, поэтому

$p_2 rightarrow mathrm{max} Rightarrow left|frac{z_3}{z_1z_2+z_1z_3+z_2z_3}right|^2 rightarrow mathrm{max}$

Либо, что то же самое

$left|z_1 + z_2 + frac{z_1z_2}{z_3}right|^2 rightarrow mathrm{min}$

$left|r_1 + jx_1 + r_2 + jx_2 +frac{(r_1 + jx_1)( r_2 + jx_2)}{jx_3}right|^2 rightarrow mathrm{min}$

$frac{1}{x_3^2}|(r_1x_3 + r_2x_3 + r_1x_2 + r_2x_1) + j(x_1x_3 + x_2x_3 + x_1x_2 - r_1r_2)|^2 rightarrow mathrm{min}$

Для удобства введем функцию

$f(x_1,x_2)=(r_1x_3 + r_2x_3 + r_1x_2 + r_2x_1) + j(x_1x_3 + x_2x_3 + x_1x_2 - r_1r_2)$

и исследуем ее на наличие экстремумов

$|f(x_1,x_2)|^2 rightarrow mathrm{min}$

Откуда получаем систему из двух уравнений

$frac{partial|f|^2}{partial x_1}=2mathbb{Re}(f)r_2 + 2mathbb{Im}(f)(x_2 + x_3)=0$

$frac{partial|f|^2}{partial x_2}=2mathbb{Re}(f)r_1 + 2mathbb{Im}(f)(x_1 + x_3)=0$

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

$x_1=-x_3,quad x_2=-x_3, quad p_2=frac{|u_1|^2,x_3^2,r_2}{left(r_1r_2 + x_3^2right)^2}, quad q_2=-frac{|u_1|^2,x_3^3}{left(r_1r_2 + x_3^2right)^2}$

$x_1=frac{1}{r_2}left(sqrt{r_1r_2left(x_3^2-r_1r_2right)}-r_2x_3right), quad x_2=frac{1}{r_1}left(sqrt{r_1r_2left(x_3^2-r_1r_2right)}-r_1x_3right), quad p_2=frac{|u_1|^2}{4,r_1}, quad q_2=frac{|u_1|^2,x_2}{4,r_1,r_2}$

$x_1=-frac{1}{r_2}left(sqrt{r_1r_2left(x_3^2-r_1r_2right)}+r_2x_3right), quad x_2=-frac{1}{r_1}left(sqrt{r_1r_2left(x_3^2-r_1r_2right)}+r_1x_3right), quad p_2=frac{|u_1|^2}{4,r_1}, quad q_2=frac{|u_1|^2,x_2}{4,r_1,r_2}$

Наиболее практичным представляется первое из трех решений $x_1=-x_3,quad x_2=-x_3$. И соответствующая ему активная мощность дается формулой

$p_2=frac{|u_1|^2,x_3^2,r_2}{left(r_1r_2 + x_3^2right)^2}$

Из формулы мощности видно, что мощность зависит от реактивного сопротивления связи $x_3=2pi,f,k,sqrt{L_{coil,1}L_{coil,2}}$, а значит и от частоты передачи $f$, и от геометрии взаимного расположения катушек, которая учитывается коэффициентом связи $k$.

Численное моделирование

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

SPICE модель двух катушек
SPICE модель двух индуктивно-связанных катушек

Симуляция выполнена для коэффициента связи $k=1$%, что соответствует 25 см удаления между катушками. Параметры катушек те же, что и в предыдущем параграфе, принятые для построения графика $k$.

Получается, что реактивные сопротивления каждой из катушек необходимо скомпенсировать конденсаторами $С_1$ и $С_2$. То есть настроить каждый из контуров (передающий и принимающий) в резонанс на заданной частоте. Если предположить, что величина нагрузки действительная, то величины емкостей могут быть найдены из формул $omega L_1=frac{1}{omega C_1}, omega L_2=frac{1}{omega C_2}$.

Ниже приведены два графика для переданного напряжения и переданной мощности во времени на частоте $f=10$ (кГц).

Переданное напряжение
Переданное напряжение

Переданная мощность
Переданная мощность

Из рисунков видно, что переданное напряжение оказалось приблизительно в 2.5 меньше входного, а переданная пиковая мощность — приблизительно в 4 раза меньше мощности, потребляемой от входа.

В заключении опишем, какие меры можно предпринять для увеличения передаваемой мощности:

  1. увеличить количество витков в катушках $n_1, n_2$
  2. увеличить радиус витков $a_1, a_2$
  3. увеличить частоту передачи $f$
  4. уменьшить расстояние между катушками $d$
  5. ввести магнитный сердечник, принадлежащий обеим катушкам (замкнутый либо открытый)
  6. ввести незамкнутый магнитный сердечник, принадлежащий лишь катушке-ресиверу

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

Литература

  1. Сивухин, Д. В. «Общий курс физики. Т. 3: Электричество и магнетизм.» (1990).
  2. Бессонов, Лев Алексеевич. Теоретические основы электротехники. Электромагнитное поле. Общество с ограниченной ответственностью Издательство ЮРАЙТ, 2012.
  3. Лаврентьев, М. А., and Б. В. Шабат. «Теория функций комплексной переменной.» (1972).

Автор: Антон Ткачев

Источник

Поделиться