Введение
Думаю, что многие из читателей видели хотя бы один ролик на популярных видеосервисах, где электричество передается через пустое пространство при помощи индуктивных катушек.
В этой статье мы хотим обратиться к первоосновам процесса беспроводной передачи энергии с помощью магнитного поля. Начав с рассмотрения простейшей индуктивной катушки, и вычисления ее индуктивности, мы постепенно перейдем к теории электрических цепей, в рамках которой, будет показан и обоснован способ передачи максимальной мощности при прочих равных условиях. Итак, начнем.
Магнитное поле одиночного витка с током
Рассмотрим магнитное поле одиночного витка с током. Найдем магнитное поле витка в любой точке пространства. Почему необходимо подобное рассмотрение? Потому что почти во всех книгах, по крайней мере в тех, которые удалось отыскать автору статьи, решение данной задачи ограничивается нахождением лишь одной компоненты магнитного поля и лишь вдоль оси витка — , в то время как мы отыщем закон для магнитного поля во всем пространстве.
Иллюстрация к закону Био-Савара-Лапласа
Для нахождения магнитного поля, воспользуемся законом Био-Савара-Лапласа (смотри Википедия — Закон Био-Савара-Лапласа). На рисунке видно, что центр системы координат совпадает с центром витка. Контур окружности витка обозначен как , а радиус окружности — как .По витку течет ток . — это переменная-радиус-вектор из начала координат в произвольную точку витка. — это радиус-вектор в точку наблюдения. Еще нам понадобится полярный угол — угол между радиус-вектором и осью . Расстояние от оси витка до точки наблюдения обозначим за . И наконец, — элементарное приращение радиус-вектора .
Согласно закону Био-Савара-Лапласа, элемент контура с током создает элементарный вклад в магнитное поле, который дается формулой
Теперь остановимся подробнее на переменных и выражениях, входящих в формулу. С учетом аксиальной симметрии задачи можем записать
Для того чтобы найти результирующее магнитное поле, нужно проинтегрировать по всему контуру витка, то есть
После подстановки всех выражений и некоторых тождественных преобразований получаем выражения для аксиальной и радиальной компоненты магнитного поля соответственно
Для нахождения абсолютного значения магнитного поля необходимо просуммировать компоненты по теореме Пифагора .
Продемонстрируем полученное решение на примере витка радиуса (м) и (А).
Амплитуда аксиальной компоненты магнитного поля
Амплитуда радиальной компоненты магнитного поля
Абсолютная амплитуда магнитного поля
Заметим, что для витка произвольной формы, на больших расстояниях , т.е. много больше характерного размера витка, поведение магнитного поля будет стремиться к найденному решению.
Катушка индуктивности. Магнитно-связанные катушки
Теперь, когда мы знаем решение для магнитного поля одного витка, можем найти индуктивность катушки, состоящей из витков. По определению индуктивность — это коэффициент пропорциональности между током в витке и магнитным потоком через площадь сечения витка. Мы пользуемся здесь идеальной моделью катушки, которая безразмерна по направлению своей оси симметрии. Конечно же, на практике такого не бывает. Однако, как приближенные, полученные формулы будут достаточно хороши. Хотя катушки и считаются безразмерными вдоль , необходимо задаться ненулевым радиусом сечения провода. Обозначим его , и пример равным (мм). Иначе при интегрировании магнитного потока подынтегральное выражение обратится в бесконечность.
Индуктивно связанные катушки
На рисунке изображены две магнитно связанные катушки. Пусть первая катушка имеет радиус и содержит витков, а вторая — и соответственно. Тогда для нахождения собственных индуктивностей необходимо вычислить магнитный поток каждой катушки через свое собственное сечение.
Поскольку в катушке много витков, найдем величину, называемую потокосцепление, дважды умножив на количество витков
По определению, индуктивность это коэффициент пропорциональности в формуле . Таким образом, получим собственные индуктивности катушек
Пусть центры катушек разделены расстоянием , лежат на одной оси, и их плоскости витков сориентированы параллельно. Для нахождения взаимной индуктивности, нужно вычислить потокосцепление, образуемое одной катушкой через сечение другой, то есть
Тогда взаимная индуктивность катушек дается выражением
Насколько известно автору, такие интегралы можно взять только численно.
Заметим, что как правило и . Коэффициентом связи катушек называется величина
Исследуем зависимость коэффициента связи катушек от расстояния. Для этого рассмотрим две одинаковые катушки с радиусом витков (м) и количеством витков . При этом собственная индуктивность каждой из катушек составит (мГн).
Коэффициент связи катушек от расстояния между ними
График не изменится, если одинаково изменить число витков в обеих катушках, либо одинаково изменить радиус обеих катушек. Коэффициент связи удобно выражать в процентах. Из графика видно, что даже при расстоянии между катушками в 1 (мм) коэффицент связи меньше 100%. Коэффициент падает до 10% на расстоянии порядка 60 (мм), и до 1% на 250 (мм).
Беспроводная передача энергии
Итак, нам известны индуктивности и коэффициент связи. Теперь воспользуемся теорией электрических цепей переменного тока для поиска оптимальных параметров, при которых передаваемая мощность оказалась бы максимальной. Для понимания этого параграфа читатель должен быть знаком с понятием электрического импеданса, а также с законами Кирхгофа и законом Ома. Как известно из теории цепей, две индуктивно-связанные катушки образуют воздушный трансформатор. Для анализа трансформаторов удобна Т-образная схема замещения.
Воздушный трансформатор и его эквивалентная схема
Передающую катушку слева будем условно называть «трасмиттер», а принимающую катушку справа — «ресивер». Между катушками коэффициент связи . На стороне ресивера находится потребитель, представленный нагрузкой . Нагрузка в общем случае может быть комплексной. Входное напряжение на стороне трансмиттера , а входной ток — . Напряжение, передаваемое на ресивер — , и передаваемый ток . Полный импеданс на стороне трансмиттера обозначим как , а полный импеданс на стороне ресивера .
Предполагается, что на вход схемы подается синусоидальное напряжение .
Обозначим — сопротивления и индуктивности катушек (две собственные и одна взаимная) соответственно. Тогда, согласно теории трансформатора
С другой стороны, согласно нашим обозначениям
где — полные активные сопротивления на стороне трансмиттера и ресивера соответственно, и — полные реактивные сопротивления.
Импеданс связи равен .
Найдем входной ток цепи
где знак обозначает параллельное соединение сопротивлений. Тогда напряжение, переданное на ресивер
И наведенный ток
Можем найти комплексную мощность, переданную в ресивер
Таким образом имеем выражение для комплексной мощности
Выражение для активной компоненты мощности
Выражение для реактивной компоненты мощности
В большинстве практических задач требуется передать максимальную активную мощность, поэтому
Либо, что то же самое
Для удобства введем функцию
и исследуем ее на наличие экстремумов
Откуда получаем систему из двух уравнений
Эта система имеет пять решений, два из которых нефизичны, так как приводят к мнимым значениям величин, которым полагается быть действительными. Три других физических решения приведены ниже вместе с соответствующими формулами для мощности
Решение 1
Мощность
Решение 2 и 3
Мощность для решений 2 и 3
Решение 2 и 3 нужно использовать, когда реактивное сопротивление связи достаточно велико
Когда же это не так, нужно использовать решение 1. Чаще всего в реальных ситуациях окажется мало, поэтому рассмотрим решение 1 несколько подробнее.
Решение 1: . И соответствующая ему активная мощность дается формулой
Из формулы мощности видно, что мощность зависит от реактивного сопротивления связи , а значит и от частоты передачи , и от геометрии взаимного расположения катушек, которая учитывается коэффициентом связи .
Как заметили внимательные читатели, зависимость — нелинейная. Функция достигает максимума при .
Исследование формулы мощности на экстремумы
Максимальная активная мощность при равна
Таким образом, вышеозначенная формула представляет абсолютный теоретический предел переданной активной мощности при любых условиях. При этом для реактивной мощности, переданной в ресивер, имеем
Численное моделирование
Продемонстрировать работу всей вышеизложенной теории можно, выполнив симуляцию SPICE модели нашего устройства из двух связанных катушек.
SPICE модель двух индуктивно-связанных катушек
Симуляция выполнена для коэффициента связи %, что соответствует 25 см удаления между катушками. Параметры катушек те же, что и в предыдущем параграфе, принятые для построения графика .
Получается, что реактивные сопротивления каждой из катушек необходимо скомпенсировать конденсаторами и . То есть настроить каждый из контуров (передающий и принимающий) в резонанс на заданной частоте. Если предположить, что величина нагрузки действительная, то величины емкостей могут быть найдены из формул
Ниже приведены два графика для переданного напряжения и переданной мощности во времени на частоте (кГц).
Переданное напряжение
Переданная мощность
Из рисунков видно, что на расстоянии 25 (см) переданное напряжение оказалось приблизительно в 2.5 меньше входного, а переданная пиковая мощность — приблизительно в 4 раза меньше мощности, потребляемой от входа, что согласуется с полученными формулами.
В заключении опишем, какие меры можно предпринять для увеличения передаваемой мощности:
- увеличить количество витков в катушках
- увеличить радиус витков
- увеличить частоту передачи
- уменьшить расстояние между катушками
- ввести магнитный сердечник, принадлежащий обеим катушкам (замкнутый либо открытый)
- ввести незамкнутый магнитный сердечник, принадлежащий лишь катушке-ресиверу
Пожалуй, написание этой статьи накладывает на автора обязательство изготовить и протестировать такую систему из двух катушек в лабораторных условиях, но это уже совсем другая история. Благодарю за внимание.
Литература
- Сивухин, Д. В. «Общий курс физики. Т. 3: Электричество и магнетизм.» (1990).
- Бессонов, Лев Алексеевич. Теоретические основы электротехники. Электромагнитное поле. Общество с ограниченной ответственностью Издательство ЮРАЙТ, 2012.
- Лаврентьев, М. А., and Б. В. Шабат. «Теория функций комплексной переменной.» (1972).
Комментарии (42)
DEM_dwg
30.10.2017 15:06Привет!
Вопрос, такой а можно ли будет использовать в качестве провода соединяющего катушки, медный лист. Т.е. сделать питание машинки с помощью одной шины которая будет соприкасаться с медным листом.
DEM_dwg
30.10.2017 15:31Т.е. от связи катушек можно будет полностью избавиться???
Но зачем же тогда они у вас на схеме соединены?938MeV Автор
30.10.2017 16:02Эта схема — эквивалентная. Она используется только для расчета. Реальная схема — та, что выше. Точно так же и в симуляциях. Симулятор выдает ошибку, если нет общей земли. Но в реальности общая земля не нужна.
kirillaristov
30.10.2017 15:55Что это за необычная запись уравнений такая?
938MeV Автор
30.10.2017 16:03Каких именно уравнений?
kirillaristov
31.10.2017 11:21При первой загрузке вместо красивых формул подгрузилась их текстовая реализация. Смотрел через opera mini.
Victor_koly
30.10.2017 16:23А можете теперь для указанных соотношений индуктивностей объяснить, при каких условиях выполняется это:
x_1 = x_2 =-x_3
Или может там для возможности такого импеданса нужны какие-то емкости ещё в цепи?938MeV Автор
30.10.2017 16:25Было немножко неясно в статье написано. Уже исправлено. Нужно добавить в контур каждой катушки по конденсатору, настроенному в резонанс с собственной индуктивностью катушки. Более точно: нужно избавиться от реактивного сопротивления, какое бы оно ни было, и в ресивере, и в трансмиттере.
FarsenaL
30.10.2017 16:38+1Серьёзный аналитический подход. Не пробовали с FEM-расчетом сравнить и проверить на сколько сложнее/проще, точнее/ресурсозатратнее получается?
938MeV Автор
30.10.2017 16:49Не пробовал, но думаю, что FEM был бы полезен для конкретных геометрий катушек и их взаимных положений. С помощью FEM можно аккуратно посчитать индуктивности и коэффициент связи. А дальше все то же самое, насколько я понимаю.
Alyoshka1976
30.10.2017 16:46+2Все очень интересно, но исправьте, пожалуйста, фамилию немецкого электротехника в тексте статьи на правильную.
vanxant
30.10.2017 18:14С пунктом 3 рекомендаций (про увеличение частоты) точно не ошиблись? В формуле для p2 в числителе вторая степень x3, и соответственно, частоты, а в знаменателе уже четвёртая.
Ну и в целом эта рекомендация кажется контр-интуитивной. С ростом частоты растёт сопротивление катушек, а значит и потери на них. Разве нет?938MeV Автор
31.10.2017 02:20Это очень хороший вопрос! Если повысить частоту, то реактивные сопротивления возрастут. Возрастет и сопротивление связи x3. Если удастся сделать его достаточно большим (x32 > r1r2), то можно будет использовать решение #2 или #3, и уже передавать активную мощность (p2 = |u1|2/4r1), что вообще не зависит от x3!
HOMPAIN
31.10.2017 02:33На индуктивностях потерь не бывает(не считая излучения магнитного поля, но оно мизерное на таких частотах). Они схожи по смыслу с конденсаторами, только накапливают ток, а потом выдают его обратно.
За один цикл колебания поля передаётся одинаково количество энергии, не зависимо от частоты. Соответственно чем выше частота, тем больше циклов на единицу времени, больше переданной энергии. А потери на сопротивлении катушки всегда постоянны.
То что с ростом частоты растёт импеданс катушки по факту не на что не влияет. Никто не запрещает просто увеличить напряжение, что бы через цепь тёк тот же ток.GermanRLI
31.10.2017 16:19В идеальных индуктивностях, которые в приблизительных расчетных моделях — не бывает.
В реальных индуктивностях вылезают межвитковые ёмкости (в результате у реальной катушки есть собственная резонансная частота, после которой импеданс, наоборот падает), активное сопротивление, на высоких частотах — ещё и скин-эффект, который активное сопротивление увеличивает, на больших токах и, следовательно сильных полях — насыщение в сердечнике, ну и нагрев провода, который ещё больше увеличит активное сопротивление.
Так что с ростом частоты — потери в катушке растут.
armature_current
31.10.2017 02:25Ваша конечная формула для мощности на приемнике в пределе при выполнении советов 1 и 3 дает обратную зависимость от взаимоиндуктивности в квадрате (если произведение активных сопротивлений r1 и r2 будет гораздо меньше индуктивного сопротивления х3), т.е. выводы противоречат полученному выражению максимума мощности. И с каждым последующим увеличением частоты или количества витков, что в обоих случаях приводит к увеличению индуктивного сопротивления еще больше будет приводить к противоречию.
938MeV Автор
31.10.2017 02:40Все верно. Дело в том, что p2(x3) — это функция, у которой есть экстремум. Я предположил, что в большинстве реальных ситуаций, мы подбираемся к этому экстремуму слева, то есть должны увеличивать x3. Наверно, такое предположение слегка опрометчиво. Добавлю в статью исследование функции p2(x3).
armature_current
31.10.2017 09:03Спасибо за оперативные пояснения. Вопрос в продолжение дискуссии. Верно ли утверждение, что при условии принятых допущений (в силу указанных элементов схемы замещения) максимальное значение передаваемой мощности через трансформатор не превышает 2,5%? Я вроде бы смотрю, но в схеме замещения нигде не вижу оговорок про магнитные свойства среды, в которой осуществляется передача энергии и естественно прихожу к выводу, что полученное выражение распространяется на все возможные варианты, в том числе и на обычный трансформатор с магнитным сердечником. Уточните пожалуйста конечный результат, при каких условиях он получен
938MeV Автор
31.10.2017 09:48По формулам получается p2 = |u1|2/4r1. Если я нигде не ошибся — выходит что так. Это максимально возможная мощность. Магнитные свойства среды не влияют, если пренебречь потерями на перемагничивание и вихревыми токами в среде. Добавлением магнитного сердечника (в первом приближении) вы просто увеличиваете индуктивную связь k.
Иными словами, формулы должны работать для любого трансформатора как с, так и без магнитного сердечника.
Sdima1357
А резонанс не поможет?
Sdima1357
Из Википедии:
"
Использование резонанса несколько увеличивает дальность передачи. При резонансной индукции передатчик и приёмник настроены на одну частоту. Производительность может быть улучшена ещё больше путём изменения формы волны управляющего тока от синусоидальных до несинусоидальных переходных формы волны. "
mayorovp
Резонас тут уже используется. Неявно и с самого начала. Потому что при разных частотах активно использующееся тут понятие "импеданс" не имеет смысла.
AntonSt
По-моему, "импеданс" будет иметь смысл при любой фиксированной частоте, не обязательно резонансной.
mayorovp
Читайте комментарий выше, там все написано.
Sdima1357
то mayorovp
«Резонас тут уже используется»
Спасибо, я понял, что он возникает неявно при оптимизации параметров, просто в статье этого не было указано явно. (сейчас указано)
938MeV Автор
У меня эта цитата из Википедии вызывает одни сомнения. Уж не знаю, как несинусоидальные формы сигналов могут улучшить «производительность».
Sdima1357
И я не знаю :(
Sdima1357
Видимо вот отсюда:
РАЗРАБОТКА МОДУЛЯ СИСТЕМЫ БЕСКОНТАКТНОЙ ПЕРЕДАЧИ ЭНЕРГИИ ПИТАНИЯ ИМПЛАНТАТОВ
ftf.kubsu.ru/htmlfiles/raspis/tom2.pdf
Nataly75
Как раз на столе аж 4 пояса от пульсометров. Частота 5 кГц. Стреляют метра на полтора. 3в батарейки хватает года на 2 (в зависимости от использования).
artskep
Возможно это отзвук борьбы правок фриков.
Буквально недавно на местном форуме на такого натолкнулся — там и «вихри» и «катушки Мишина», которые лечат рак… И там было что-то про про М-фазу, особый смысл меандра…
Вобщем полный набор фрика в одном флаконе вкупе с полной безграмотностью в измерениях (подключить параллельно генератор, осциллограф и «катушку», и смотреть на графики — я, конечно, уже давно забыл курс метрологии, но уверен, что при такой схеме меня бы сразу на пересдачу направили бы).
938MeV Автор
Да. Катушки с конденсаторами как раз и настраиваются в резонанс. С одной стороны, кажется очевидным. Но вот выясняется, что есть и другие решения.