CC BY
43
Индукционные свойства проводящего цилиндра с азимутальной напряжённостью вихревого электрического поля
В.Г. Сапогин, Н.Н. Прокопенко, В.И. Марчук, А.С. Будяков
Введение
Современные планарные технологии создания индуктивностей или трансформаторов в частотном диапазоне до 10 ГГц используют геометрию квадрата, либо октаэдра [1-5]. Но технология многоугольников обладает одним существенным недостатком: при скачкообразном изменении ориентации токопровода в пространстве скин-эффект создаёт благоприятные условия для выбрасывания электрического заряда на подложку. Применение кольцевой геометрии должно ослабить выбрасывание электрического заряда на подложку и тем самым позволит увеличить на указанных частотах, как параметр добротности, так и коэффициент передачи трансформаторов.
Известные методы расчёта индуктивностей [6] в микрометровом диапазоне геометрических размеров могут приводить к отрицательным значениям индуктивности. Этот недостаток геометрической теории потребовал создания принципиально новых электродинамических физико-математических моделей, которые ориентированы на потребности планарной микросхемотехники.
В [7-10] был найден подход к решению задач кольцевой геометрии, в основе которого лежит закон Био-Савара-Лапласа, либо уравнения магнитостатики. Полученные там результаты дают неплохое совпадение с экспериментальными данными и могут быть применимы для частот переменного тока, при которых отсутствует фазовый сдвиг между током и потоком.
Ниже излагается оригинальный подход к решению задач цилиндрической геометрии, в которых между током и потоком может существовать конечный фазовый сдвиг. Подход объединяет в себе
возможности двух- и трёхмерных задач и не ограничивает возможный диапазон изменения частоты переменного тока.
Математическая модель расчёта
Математическая модель расчёта построена на решении первого уравнения Максвелла для цилиндрической симметрии исследуемых электромагнитных полей. В простейшей постановке задачи переменное однородное аксиальное магнитное поле формирует поток через цилиндр, при котором возникающее вихревое электрическое поле имеет одну азимутальную компоненту.
Тогда, внешнее магнитное поле имеет одну единственную компоненту В = (0,0, Вг). Его зависимость во времени имеет вид
Вг = В0со$, Ш, (1)
где В0=еотг не зависит от координат системы и в расчётах ограничено в радиальном направлении внешним радиусом цилиндра Я.
На рис. 1 представлена ориентация проводящего цилиндра по отношению к направлению поля Вг. Радиус отверстия цилиндра обозначен через г1, а произвольная толщина цилиндра обозначена через к.
Магнитное поле, имеющее одну компоненту В2, порождает в цилиндре вихревое электрическое поле, напряженность которого Е = (0, Еф ,0) также
имеет одну компоненту в случае, если вектор индукции В направлен точно по оси г.
Распределение вихревого электрического поля в пространстве кольца описывается первым уравнением Максвелла
оЕ = -—. (2)
дг
Переписанное в проекциях для полей В2 и Еф, оно имеет вид
(гогЕ) г . (3)
дг
Рис. 1. Ориентация проводящего цилиндра по отношению к осевому магнитному полю
Предполагая, что вихревое электрическое поле совпадает по фазе со скоростью изменения магнитного поля
еф = е0(г , (4)
из уравнения (3) получим уравнение, связывающее Е0(г) и В0
1 д
(Е) = юВ0. (5)
г дг
Его решение для Е0(г) имеет вид
Е0(г) = юВ0 г /2 + С1/ г. (6)
В решении (6) избавимся от особенности при г ^ 0, полагая С1=0. Из (6) видно, что азимутальная компонента вектора напряженности вихревого электрического поля - линейная функция радиуса с масштабом
Е* =юВ0 Я /2 = пуВ0 Я. (7)
Распределение вихревого поля в однородном проводящем цилиндре создает вихревые токи Фуко, плотность тока которых имеет такое же направление в пространстве и рассчитывается из локального закона Ома
Уф= Еф / Р, (8)
где р - удельное сопротивление цилиндра.
Из (8) видно, что плотность тока Фуко в однородной проводящей среде р=сот1 также является линейной функцией радиуса и также зависит от времени, как и Еф
у ф= Jo(r Ш , (9)
гдеу0(г) связано с масштабом плотности тока у* = пуВ0 Я / р (10)
соотношением
Л(г) = У* г / Я. (11)
Из соотношения (10) для масштаба плотности тока видно, что при прочих равных условиях он может достигать больших значений на высоких частотах для проводящей среды с малым значением удельного сопротивления р. Это будет приводить к существенному нагреву проводящих цилиндров на высоких частотах.
Удельная тепловая мощность, выделяемая в локальной области проводника при протекании токов Фуко, рассчитывается из закона Джоуля Ленца:
г 2
Ру = РУ'ф = Р* ^у^п2 Ш, (12)
Я
где р* = р* = п2 V2 В02 Я2 / р (13)
- масштаб объемной плотности тепловой мощности. Как видно из (13), масштаб плотности уже пропорционален квадрату частоты и для проводящей среды с произвольным значением удельного сопротивления оказывается наиболее значительным на внешней границе цилиндра. Это будет приводить к выгоранию его внешней боковой поверхности.
Поскольку напряженность вихревого поля Еф зависит линейно от радиуса, то элементарная электродвижущая сила (ЭДС) индукции также будет зависеть от радиуса и момента времени. Направления векторов напряжённости и плотности тока на силовой линии представлены на рис. 2. Из рисунка видно, что элементарная индукция квадратично зависит от радиуса
Рис. 2. Направление векторов напряжённости и плотности тока в цилиндре
Ше = Е • = Еф гйф
Е * г2Шф Я
б1п Ш.
(14)
В (14) учтено, что вектор совпадает с вектором Е.
На одной замкнутой силовой линии вихревого поля электродвижущая
сила индукции зависит от ее радиуса г2
е(г, ?) = е * —-б1п Ш
Я
(15)
(16)
где е* = 2пЯЕ* - масштаб электродвижущей силы индукции.
Усредненное по радиусу значение ЭДС сохраняет ту же зависимость от времени
< е >
1 я — /
— г •>
Я - г1
Я 2
г 1 е *
е* —— бш ш1аг = —
Я
2
3
'1
2
1+^+^
V Я Я2 ,
б1п Ш
(17)
Полный ток, индуцированный в кольце и пересекающий площадь поперечного сечения = ШНШг, может быть получен из интегрирования для фиксированного момента времени
I = | = и
í 2 л 1 - ^ Я ;
Б1П Ш,
где I* = ]*НЯ /2 = пуВ0НЯ2 /(2р)
(18) (19)
масштаб индукционного тока. Соотношения (18) и (19) позволяют определить независящее от времени усредненное значение интегрального омического сопротивления цилиндра индукционному току
< Я >=<^ = Я
= Я (1 + тх/Я + Г12/Я2) (1 - Г2/ Я2) :
(20) (21)
где Я* =в* /(3/*) - масштаб омического сопротивления цилиндра.
На рис. 3 представлен график зависимости приведенного омического сопротивления <Я>/Я* от приведенного радиуса полости г1/Я. Из него видно, что приведенное усредненное значение омического сопротивления при г1 ^ 0 совпадает с единицей, а для тонких цилиндров при г1 / Я ^ 1 стремится к + да.
Это свойство усредненного интегрального омического сопротивления позволит при реализации тонких цилиндров искусственно уменьшать тепловое энерговыделение проводника на высоких частотах изменением приведённого радиуса полости.
Приведённый радиус полости
Рис 3. Зависимость приведённого омического сопротивления от приведённого радиуса полости цилиндра
Вносимая индуктивность цилиндра
Рассчитаем значение вносимой индуктивности из энергетических соображений. Силовые линии магнитного поля Вг пронизывают весь объем цилиндра и вносят в него магнитную энергию W, значение которой зависит от времени
г В2
w (г , t) = | —^dv. (22)
V 0
При усреднении за период изменения поля эта энергия имеет конечное значение:
1 T 02 л, т>2_г,2;
'Т = _
1г 2 , гВо dV В0 пЯ2 к < ж >т = -1ООБ2 юtdt| —-= —-. (23)
Т 0 V 2М-0 0
Поскольку эта же энергия магнитного поля создает вихревое электрическое поле, которое формирует индукционный ток, будем считать, что усредненное значение магнитной энергии (23) совпадает с усредненным за период значением энергии, запасенной индуктивностью
Ь 2
< ж >Т = - < /2 >Т, (24)
•2 ~ где под < / >Т понимается усреднение, выполненное за такой же
период изменения тока
^2
(25)
2
< / > Т = Т 2
• 2 С ,.2\2 I*
1 - ^
Я у
Подставляя в (24) соотношения (23) и (25), получим для вносимой индуктивности выражение
Ь = Ь* (1 - Г!2/Я2 )-2, (26)
где Ь* = 4р 2/(ц 0 V 2V) (27)
- масштаб индуктивности, а
V = пЯ2 к (28)
- объем цилиндра с полостью.
На рис. 4 представлена зависимость приведенной индуктивности Ь/Ь* от приведенного радиуса полости цилиндра г1/Я. Как видно из рисунка, ее значение совпадает со значением 1 при г1 ^ 0. Для радиально тонких
цилиндров при выполнении условия г1 / Я ^ 1 индуктивность резко
увеличивается по сравнению с масштабом.
Это свойство вносимой индуктивности позволяет ее увеличивать за счет изменения радиальной толщины при реализации тонких колец или тонкостенных удлиненных цилиндров. В этом случае они нагреваются меньше.
Рис. 4. Зависимость приведённого значения индуктивности от приведённого радиуса цилиндрической полости
Оценки вносимой индуктивности
Поскольку современные материалы от проводников до собственных полупроводников имеют удельное сопротивление, изменяющееся в диапазоне
15,5 нОмм<р<2,5 кОмм,
то формулы для расчета вносимой индуктивности (26) и (27) дают широкие возможности для ее изменения.
В таблице 1 рассчитаны значения частоты переменного тока V в герцах, которые для объёма бесполого цилиндра с радиусом Я=1 см и высотой к=1 см, дают значения индуктивности в интервале от 100 мГн до 1 нГн для
различных проводящих материалов. В таблице 1 введены обозначения: материал КЭФ - кремний электронный легированный фосфором с удельным сопротивлением р=0,1 Ом-см; германий взят с удельным сопротивлением р=47 Ом-см; кремний взят с удельным сопротивлением р=230 кОм-см).
Из таблицы 1 видно, что большие значения индуктивности достигаются на малых частотах. Для металлов эти частоты заключены в диапазоне от сотых герца (медь) до герц (нихром). Для полупроводников такие значения индуктивности достижимы на частотах больше 1 кГц. Малые значения индуктивности достигаются на высоких частотах. Для металлов это килогерцы и десятки килогерц, а для полупроводников - десятки и тысячи гигагерц.
Таблица №1
Значения частоты переменного тока V в герцах для различных проводящих материалов
Материал/ индуктивность Медь Вольфрам Никель Нихром КЭФ Германий Кремний
100 мГн 4,9310-2 1,56-10-1 1,95-10-1 3,18 3,18-103 1,5-106 7,32-109
10 мГн 1,56-10-1 4,93-10"1 6,18-10"' 10,1 1,01-104 4,73-106 2,32-1010
1 мГн 4,9310-1 1,56 1,95 31,8 3,18-104 1,5-107 7,32-1010
100 мкГн 1,56 4,93 6,18 101 1,01-105 4,73-107 2,32-10п
10 мкГн 4,93 15,6 19,5 318 3,18-105 1,5-108 7,32-Ю11
1 мкГн 15,6 49,3 61,8 1,01-103 1,01-106 4,73-108 2,32-1012
100 нГн 49,3 156 195 3,18-103 3,18-106 1,5-109 7,32-1012
10 нГн 156 493 618 1,01-104 1,01-107 4,73-109 2,32-1013
1 нГн 493 1560 1950 3,18-104 3,18-107 1,5-1010 7,32-1013
Экспериментальные исследования эффекта вносимой индуктивности для катушек сантиметрового геометрического размера и выше возможно введением в них осесимметричных полых цилиндрических сердечников с различной проводимостью [7,9].
Выводы
Предложен аналитический метод расчёта индукционных и индуктивных свойств цилиндра произвольной высоты с азимутальной плотностью вихревого тока.
Метод основан на решении первого уравнения Максвелла для цилиндрической симметрии исследуемых электромагнитных полей, в которых между током и потоком может существовать конечный фазовый сдвиг.
Получены радиальные распределения азимутальной компоненты напряжённости вихревого электрического поля, плотности токов Фуко и удельной тепловой мощности, выделяемой в локальной области проводника.
Определены зависимости среднего значения ЭДС в цилиндре, индукционного тока и интегрального омического сопротивления от приведённого радиуса полости цилиндра.
Получена формула для вычисления вносимой индуктивности проводящего цилиндра.
Её зависимость от приведённого радиуса полости проводящего цилиндра, квадрата удельного сопротивления, квадрата частоты переменного тока и объёма позволяет изменять вносимую индуктивность в широком диапазоне значений.
Оценки, проведённые для существующих проводников и полупроводников, указывают на то, что больших значений индуктивности можно добиться на низких частотах для материалов с малым удельным сопротивлением и малых значений индуктивности - на СВЧ и КВЧ для материалов с большим удельным сопротивлением.
Заключение
В заключение обсудим пределы применимости построенной физико-математической модели. Как известно, это можно сделать только в том случае, если удаётся продвинуться в дальнейшем обобщении и понимании полученных результатов.
Вдумчивый читатель скажет, что при возникновении токов Фуко в цилиндре обязательно появится дополнительное переменное неоднородное магнитное поле (поле отклика вещества В1), которое будет иметь две компоненты B1r и B1z. и влияние которого на полученные вычисления
необходимо учитывать.
Читатель, конечно, прав и такой учёт при вычислении вносимой индуктивности должен быть сделан. В дальнейшем нужно учесть и потери энергии суммарного магнитного поля, возникающие за счёт джоулева нагрева проводящего кольца.
Но, поскольку учёт этих явлений является дальнейшим расширением рамок полученной физико-математической модели, то удовлетворим это замечание введением некоторого математического условия, которое ограничит рамки предложенного исследования.
< B12 >T << < В02 >T. (29)
Что означает следующее: энергия внешнего магнитного поля, вносимая в проводящий цилиндр за период, оказывается значительно большей, чем энергия магнитного поля токов Фуко, вносимая в проводник за то же самое время. Как показано в дальнейшем, даже частичный учёт диамагнетизма токов Фуко, приводит к тому, что условие (29) будет выполняться в достаточно широкой части низкочастотной области.
Работа выполнена при поддержке Российского фонда фундаментальных исследований, грант №12-08-00654/12 (2012-2013 г.г.).
Литература:
1. El-Sharawy. USP. Monolithic balanced RF power amplifier. US 6, 424, 227 B1; jul. 23, 2002.
2. Komijani et al. USPAP. Reconfigurable distributed active transformers. US 2003/0169113 A1; Sep. 11, 2003.
3. Aoki et al. USPAP. Distributed circular geometry power amplifier architecture. US 2005/0030098 A1; feb. 10, 2005.
4. Sanderson. USPAP. High Q monolithic inductors for use in differential circuits. US 2006/0220737 A1; Oct. 5, 2006.
5. Lee et al. USPAP. Systems and method for power amplifiers with voltage bousting multi-primary transformers. US 2008/0164941 A1; Jul. 10, 2008.
6. Калантаров П.Л. , Цейтлин Л.А. Расчёт индуктивностей: справочная книга. - Л.: Энергоатомиздат, 1986. - 488 с.
7.Сапогин В.Г., Крутчинский С.Г., Прокопенко Н.Н., и др. Интегральные индуктивности и трансформаторы аналоговых микросхем СВЧ - диапазона. - Шахты: ГОУ ВПО «ЮРГУЭС», 2010. - С. 273.
8.Сапогин В.Г., Прокопенко Н.Н., Марчук В.И., Манжула В.Г. Погонная индуктивность цилиндрических проводников с аксиальной плотностью тока. Инженерный вестник Дона. №4(часть 1), 2012. http://www.ivdon.ru/magazine/archive/n4t1y2012/1264
9.Сапогин В.Г., Марчук В.И., Манжула В.Г. Промежуточный отчёт по гранту РФФИ «Теоретические основы проектирования интегральных индуктивностей для сложных функциональных блоков и IP - модулей систем связи и телекоммуникаций нового поколения». Проект № 12-08-00654а. 2012, - С. 79.
10. Сапогин В.Г., Прокопенко Н.Н., Будяков А.С. О методах измерения малых индуктивностей для радиотехнических систем на кристалле. Инженерный вестник Дона. №1, 2013.
http://www.ivdon.ru/magazine/archive/n1y2013/1566
