Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

INDUCTION OF SOLID CONDUCTING CYLINDER WITH AZIMUTH DENSITY OF VORTEX CURRENT IN LOW-FREQUENCY RANGE

Sapogin V.G. 1 Prokopenko N.N. 2 Manzhula V.G. 3 Sapuncov N.E. 1 Nesturina Е.Е. 1
1 Taganrog Institute of Technology Southern Federal University
2 Institute of service and business (branch) Federal State Budget Educational Institution of Higher Professional Education Don State Technical University ISB (branch) FSBEI of HPE «DSTU»
3 Non-state institute of higher education «Institute of communicative technologies»
The analytical method of calculation of induction of solid nonmagnetic conducting cylinder by arbitrary geometry size with azimuth density of vortex current under existence of resistance loss and diamagnetic properties of Foucault currents has been proposed. The radial distributions of azimuth component of vortex electric field’s strength vector, densities of Foucault currents and specific heat power, being emanated in local domain of conducting cylinder, have been obtained. The formula for insertion inductance of solid cylinder has been obtained from energy consideration. The inductance of conducting cylinder comes to zero for two characteristic frequencies of the system and has positive and negative values. The estimates, executed for brass and electronic silicon film, point to that the large values of inductance can be achieved in the region of small resistance loss and minor influence of Foucault currents’ diamagnetism under supplementary execution of condition h/R << 1.
insertion inductance
induction
vortex electric field
alternating current
magnetic field flux
1. Каlantarov P.L., Ceitlin L.А. Raschet inductivnostei: spravochnaia kniga. L.: Energoatomizdat, 1986. 488 р.
2. Manzhula V.G. Isklyuchenie strukturnoj, funktsional’noj i skhemotekhnicheskoj izbytochnosti pri sinteze analogovykh sistem v korpuse. Nauchno-tekhnicheskij vestnik Povolzhya. 2011. no. 2. pр. 123–127
3. Manzhula V. G. Funktsionalno integrirovannaya mikroehlektronnaya sistema zashhity na osnove bystrodejstvuyushhego datchika temperatury. Datchiki i sistemy. 2012. no. 7. pр. 18–22.
4. Savelev I.V. Kurs obschei phiziki. T.2. M.: Nauka, Glavnaia redakcia fiziko-matematicheskoi literatury. 1988. 496 р.
5. Sapogin V.G., Krutchinskiy S.G., Prokopenko N.N., i dr. Integralnie induktivnosti I transformatory analogovykh mikposkhem SVCH – diapazona. Shakhty: GOU VPO UYRGUES, 2010. 273 р.
6. Sapogin V.G., Marchuk V.I., Manjula V.G. Promezhutochnyy otchet po grantu RFFI Teoreticheskie osnovy proektirovaniya integralnykh inductivnostei dla slojnyh funktsionalnykh blokov i IP – module system svyazi I telekommunicatsiy novogo pokoleniya. Proekt no. 12-08-00654а. 2012, 79 р.
7. Sapogin V.G., Prokopenko N.N., Marchuk V.I., Manjula V.G. Pogonnaya induktivnost tsilindricheskikh provodnikov s aksialnoi plotnostyu toka. Inzhenernyy vtstnik Dona. no. 4, 2012. www.ivdon.ru/magazine/ archive / n4t1y2012/1264.
8. Aoki et al. USPAP. Distributed circular geometry power amplifier architecture. US 2005/0030098 A1; feb. 10, 2005.
9. El-Sharawy. USP. Monolithic balanced RF power amplifier. US 6, 424, 227 B1; jul. 23, 2002.
10. Lee et al. USPAP. Systems and method for power amplifiers with voltage bousting multi-primary transformers. US 2008/0164941 A1; Jul. 10, 2008.
11. Komijani et al. USPAP. Reconfigurable distributed active transformers. US 2003/0169113 A1; Sep. 11, 2003.
12. Sanderson. USPAP. High Q monolithic inductors for use in differential circuits. US 2006/0220737 A1; Oct. 5, 2006.

Современные планарные технологии создания индуктивностей или трансформаторов в частотном диапазоне до 10 Ггц используют геометрию квадрата либо октаэдра [8–10]. Но технология многоугольников обладает одним существенным недостатком: при скачкообразном изменении ориентации токопровода в пространстве скин-эффект создаёт благоприятные условия для выбрасывания электрического заряда на подложку. Применение кольцевой геометрии должно ослабить выбрасывание электрического заряда на подложку и тем самым позволит увеличить на указанных частотах как параметр добротности, так и коэффициент передачи трансформаторов [2, 3].

Известные методы расчёта индуктивностей [1] в микрометровом диапазоне геометрических размеров могут приводить к отрицательным значениям индуктивности. Этот недостаток геометрической теории потребовал создания принципиально новых электродинамических физико-математических моделей, которые ориентированы на потребности планарной микросхемотехники.

В [5–7] был найден подход к решению задач кольцевой геометрии, в основе которого лежит закон Био‒Савара‒Лапласа либо уравнения магнитостатики. Полученные там результаты дают неплохое совпадение с экспериментальными данными и могут быть применимы для частот переменного тока, при которых отсутствует фазовый сдвиг между током и потоком.

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

Физико-математическая модель расчёта

Предположим, что сплошной проводящий немагнитный цилиндр находится во внешнем однородном переменном магнитном поле, направленном по оси z, sapog01.wmf. Его зависимость во времени имеет вид

sapog02.wmf (1)

где B0 = const не зависит от цилиндрических координат системы и ограничено в радиальном направлении внешним радиусом цилиндра R, а в аксиальном направлении произвольной высотой цилиндра, обозначенной через h.

При такой ориентации переменное магнитное поле может порождать в цилиндре вихревое электрическое поле, напряженность которого sapog03.wmf имеет одну компоненту. Впервые экспериментальное доказательство существования такой возможной ориентации вихревого электрического поля было обнаружено при создании бетатронов – индукционных ускорителей электронов. В них как раз ускорение электронов и осуществляется вихревым электрическим полем, силовые линии которого представляют собой концентрические окружности. Они формируются электромагнитом специальной формы (см., например, [4]).

Компоненты полей связаны между собой первым уравнением Максвелла, записанным в проекциях

sapog04.wmf (2)

Полагая в (2), что фаза вихревого электрического поля совпадает по фазе со скоростью изменения магнитного поля

sapog05.wmf (3)

из уравнения (3) получим уравнение, связывающее E0(r) и B0,

sapog06.wmf (4)

Его решение для E0(r) имеет вид

sapog07.wmf (5)

В решении (5) избавимся от особенности при r → 0, полагая C1 = 0. Это условие приводит к тому, что при B0 = 0 E0 = 0. Из (2) видно, что азимутальная компонента вектора напряженности вихревого электрического поля – линейная функция радиуса r с масштабом

sapog08.wmf (6)

В масштабе напряжённости учтено соотношение ω = 2πf. Тогда (5) с учётом (6) имеет вид

sapog09.wmf при 0 ≤ r ≤ R. (7)

Вихревое электрическое поле в проводящем цилиндре создает вихревые токи Фуко, плотность тока которых имеет такое же направление в пространстве sapog10.wmf. Плотность токов Фуко рассчитывается из дифференциального закона Ома

sapog11.wmf (8)

где r – удельное электрическое сопротивление проводящего цилиндра.

Из (8) видно, что плотность тока Фуко в однородной проводящей среде r = const также является линейной функцией радиуса и также зависит от времени, как и Eφ

sapog12.wmf (9)

где j0(r) связано с масштабом плотности тока

sapog13.wmf (10)

соотношением

sapog14.wmf (11)

Из соотношения (10) видно, что при прочих равных условиях масштаб тока может достигать больших значений на высоких частотах для проводящей среды с малым значением удельного сопротивления r. Это будет приводить к нагреву проводящего цилиндра на высокой частоте.

Удельная тепловая мощность, выделяемая в локальной области проводника при протекании токов Фуко, рассчитывается из закона Джоуля Ленца:

sapog15.wmf (12)

где sapog16.wmf (13)

‒ масштаб объемной плотности тепловой мощности.

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

Тепловая мощность, выделяемая во всём цилиндре, находится из (12)

sapog17.wmf (14)

где масштаб тепловой мощности имеет вид

sapog18.wmf (15)

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

sapog19.wmf (16)

где sapog20.wmf (17)

‒ масштаб индукционного тока.

Частичный учёт магнитных свойств токов Фуко

Вихревая плотность токов Фуко (9) создаёт в пространстве цилиндра собственное магнитное поле sapog21.wmf. Это магнитное поле отклика вещества для исходной ориентации вихревых токов в цилиндрической системе координат имеет две проекции B1r и B1z. Направление проекции B1z во всех точках цилиндра противоположно внешнему магнитному полю В0. Предполагая, что компонента B1r не зависит от координаты z, распределение собственного магнитного поля можно рассчитать из второго уравнения Максвелла, записанного в проекциях (частная производная компоненты B1r по переменной z опущена):

sapog22.wmf (18)

Интегрируя (18) для начального условия B1z(0) с учётом (9), получим

sapog23.wmf (19)

где sapog24.wmf (20)

‒ масштаб индукции собственного магнитного поля. Его можно связать с индукцией внешнего магнитного поля безразмерным параметром β

sapog25.wmf (21)

который указывает на влияние диамагнетизма токов Фуко

sapog26.wmf (22)

где sapog27.wmf (23)

– масштаб частоты системы. Из (22) видно, что диамагнетизм слабо проявляет себя при b → 0 и его влияние велико на высоких частотах, когда b >> 1. При значении b = 1 масштабы индукций магнитных полей одинаковы.

Поскольку современные материалы от проводников до собственных полупроводников имеют удельное сопротивление, изменяющееся в диапазоне 15,5 нОм·м < ρ < 2,3 кОм·м, то формула для расчета масштаба частоты (23) даёт широкие возможности для ее изменения.

В табл. 1 приведены значения масштаба частоты для восьми проводящих материалов.

Таблица 1

Масштаб частоты для материалов с различным удельным сопротивлением для радиуса цилиндра R = 2 мм.

№ п/п

Материал

Удельное сопротивление (Ом⋅м)

Масштаб частоты f*, Гц

1

Медь

15,5⋅10–9

1,96⋅103

2

Вольфрам

48,9⋅10–9

16,2⋅103

3

Никель

61,4⋅10–9

7,78⋅103

4

Латунь (марг.)

2,12⋅10–7

26,9⋅103

5

Нихром

1,0⋅10–6

1,27⋅105

6

KEF (кремний электронный)

1,0⋅10–3

1,27⋅108

7

Германий

0,47

5,95⋅1010

8

Кремний

2,3⋅103

2,91⋅1014

Как видно из табл. 1, масштаб частот современных материалов лежит в широком диапазоне частот от 2⋅103 до 3⋅1014 Гц.

Энергетический баланс системы

Средняя за период Т изменения поля магнитная энергия, поступающая в цилиндр, находится из интегрирования

sapog28.wmf (24)

где результирующая индукция магнитного поля системы

sapog29.wmf (25)

Тогда < W > состоит из двух слагаемых: средней за период энергии внешнего поля

sapog30.wmf (26)

и средней энергии магнитного поля отклика вещества

sapog31.wmf (27)

Интегрирование за период скалярного произведения sapog32.wmf обращается в нуль из-за наличия фазового сдвига π/2 между индукциями B1 и B0. В связи с этим между переменными полями B1 и B0 возникает своеобразная интерференция, которая и суммирует обе энергии

sapog33.wmf (28)

Вносимая за период магнитная энергия (28) расходуется на нагрев проводящего цилиндра Wq и ту часть энергии, которая запасается в индуктивность системы

sapog34.wmf (29)

Заметим, что в соотношении (29) не учтена магнитная энергия, которая поставляется в систему радиальной компонентой магнитного поля отклика вещества B1r. Её влияние на процессы дополнительной генерации энергии магнитного поля в пространстве проводящего цилиндра будут учтены в отдельной работе.

Тепловая энергия, расходуемая на нагрев цилиндра за период, вычисляется интегрированием по времени (14)

sapog35.wmf (30)

Среднее за период значение квадрата тока найдём из (16)

sapog36.wmf (31)

Индуктивные свойства сплошного немагнитного цилиндра

Подставляя в (28) и (29) вычисленные энергии, получим для приведённой индуктивности цилиндра соотношение

sapog37.wmf (32)

где sapog38.wmf (33)

‒ масштаб индуктивности, который может принимать большие значения для аксиально тонких плёнок при выполнении условия h/R << 1.

Как видно из (32), индуктивность проводящего цилиндра обращается в нуль на двух характеристических частотах F1 и F2 (F1 < F2):

sapog39.wmf (34)

sapog40.wmf (35)

Нули индуктивности делят весь частотный интервал на три части: низкочастотный (выполнено условие 0 < β < β1), среднечастотный, в котором индуктивность принимает отрицательные значения (выполнено условие β1 < β < β2) и высокочастотный, в котором она снова принимает положительное значение (выполнено условие β2 < β < ∞).

Как видно из (32), в низкочастотном диапазоне индуктивность при β → 0 (диамагнетизм токов и омические потери небольшие) может принимать положительные, большие по сравнению с масштабом значения, зависящие от частоты по закону

sapog41.wmf (36)

Этот результат представляет интерес для создания больших индуктивностей на низких частотах в малом объёме. При b << 1 формула (36) переходит в формулу

sapog42.wmf (37)

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

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

Исследуем возможности увеличения индуктивности в низкочастотной области за счёт уменьшения частоты переменного поля f.

Сплошной цилиндр из марганцовистой латуни будет иметь положительную индуктивность (r = 2,12⋅10–7 Ом×м, (см. табл. 2)) при геометрических размерах R = 2 мм, h = 1 см (выполнено условие h/R > 1). Масштаб частоты f* = 2,69⋅104 Гц. Масштаб индуктивности L* = 1,58 нГн. Первый нуль индуктивности находится на первой характеристической частоте F1 = 0,1605f* = 4,3 кГц. Приведём значения индуктивности цилиндра, которые могут быть реализованы на частотах, меньших первой характеристической частоты в целое число раз.

Таблица 2

Зависимость индуктивности от частоты для условия h/R > 1.

F(частота)

L/L*

L (микрогенри)

F1/2

77,3

0,122

F1/3

232

0,366

F1/4

464

0,734

F1/5

775

1,22

F1/6

1,16⋅103

1,84

F1/7

1,63⋅103

2,57

F1/8

2,17⋅103

3,43

F1/9

2,80⋅103

4,41

F1/10

3,49⋅103

5,51

Как видно из табл. 2, увеличение индуктивности за счёт уменьшения омических потерь и влияния диамагнитных свойств токов Фуко может быть в сорок пять раз больше в диапазоне изменения частоты переменного магнитного поля от 0,43 до 2,15 кГц.

Сделаем оценки индуктивности, представляющие интерес для планарных технологий, при выполнении условия h/R << 1. Для тонкой плёнки h = 1 мкм и радиусом R = 1 см, выполненной из электронного кремния (см. табл. 3), имеем: масштаб частоты f* = 5,07⋅1010 Гц. Масштаб индуктивности L* = 39,5 нГн. Первый нуль индуктивности находится на первой характеристической частоте F1 = 0,1605f* = 8,14 ГГц. Приведём значения индуктивности цилиндра, которые могут быть реализованы на частотах, меньших первой характеристической частоты в целое число раз в диапазоне от 0,81 Ггц до 4,07 Ггц.

Применение плёнки из того же материала с той же толщиной, но радиусом в 1 см даёт масштаб частоты f* = 5,07⋅106 Гц. Масштаб индуктивности L* = 0,39 мГн. Первый нуль индуктивности находится на первой характеристической частоте F1 = 0,1605f* = 0,81 МГц. Значения индуктивности изменяются в диапазоне от 30 мГн до 1,4 Гн в интервале изменения частоты от 81 кГц до 0,4 МГц.

Таблица 3

Зависимость индуктивности от частоты для условия h/R << 1

F(частота)

L/L*

L (микрогенри)

F1/2

77,3

3,05

F1/3

232

9,17

F1/4

464

18,4

F1/5

775

30,6

F1/6

1,16⋅103

45,9

F1/7

1,63⋅103

64,3

F1/8

2,17⋅103

85,7

F1/9

2,80⋅103

110

F1/10

3,49⋅103

138

Выводы

  • Предложен аналитический метод расчёта индуктивных свойств сплошного проводящего немагнитного цилиндра произвольной высоты и радиуса с азимутальной плотностью вихревого тока при наличии омических потерь и диамагнитных свойств токов Фуко.
  • Получены радиальные распределения азимутальной компоненты напряжённости вихревого электрического поля, плотности токов Фуко и удельной тепловой мощности, выделяемой в локальной области проводящего цилиндра.
  • Из энергетических соображений получена формула для вычисления вносимой индуктивности сплошного проводящего цилиндра.
  • Показано, что индуктивность проводящего цилиндра обращается в нуль для двух характеристических частот F1 и F2.
  • Нули индуктивности делят интервал частот на три части: низкочастотный, среднечастотный и высокочастотный.
  • На низких и высоких частотах индуктивность проводящего цилиндра принимает положительные значения, а на средних частотах индуктивность отрицательна.
  • Показано, что вносимая индуктивность зависит от частоты и имеет большие абсолютные значения в диапазоне низких частот.
  • Оценки, проведённые для латунного цилиндра в диапазоне низких частот, при выполнении условия h/R > 1 указывают на то, что больших положительных значений индуктивности латуни на уровне 1 мкГн, отличающихся в 45 раз, можно добиться в интервале частот от 0,43 кГц до 2,15 кГц.
  • Оценки, проведённые для плёнки из электронного кремния в диапазоне низких частот, при выполнении условия h/R << 1 указывают на то, что больших значений индуктивности на уровне 50 мкГн, отличающихся в 45 раз, можно добиться в интервале частот от 0,81 ГГц до 4,07 ГГц.

Работа выполнена при поддержке Российского фонда фундаментальных исследований, грант № 12-08-00654/12 (2012-2013 гг.).

Рецензенты:

Бубнов В.А., д.т.н., профессор, заведующий кафедрой естественнонаучных дисциплин, ГБОУ ВПО «Московский городской педагогический университет», г. Москва;

Крутчинский С.Г., д.т.н., профессор кафедры систем автоматического управления, ФГАОУ ВПО «Южный федеральный университет», г. Ростов-на-Дону.

Работа поступила в редакцию 16.12.2013.