Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Кураев А.А., Попкова Т.Л., Синицын А.К. Электродинамика и распространение радиоволн.pdf
Скачиваний:
262
Добавлен:
02.05.2014
Размер:
4.84 Mб
Скачать

V

= 1

. Коэффициент при первых производных по времени определяет

ф

εa µa

 

 

 

пространственное затухание волнового процесса в среде с джоулевыми потерями.

Таким образом, из УМ непосредственно следует, что электромагнитные процессы имеют волновой характер и, следовательно, любое ЭМП представимо как некоторое наложение электромагнитных волн.

5.2. Электродинамические потенциалы

Установим далее формулы, выражающие напряженность электрическо-

го поля E иr индукцию магнитного поля B через электродинамические по-

тенциалы A (векторный) и Ф (скалярный). Для этого обратимся к системе уравнений Максвелла.

Четвертое уравнение Максвелла утверждает, что поле B чисто вихревое:

r

(5.8)

div B = 0.

Из (5.8) следует, что поле B следует определить как

 

B = rot A ,

(5.9)

так как div rot A = 0 .

 

Воспользуемся теперь вторым уравнением Максвелла

 

r

 

rot E = −

B

.

(5.10)

 

 

t

 

Подставляя в (5.10) выражение B (5.9) и имея в виду, что координаты и время в используемой лабораторной системе отсчета независимы, получим

r

 

 

 

 

 

 

A

= 0 .

(5.11)

rot E +

 

 

 

t

 

 

Из (5.11) с учетом того, что rot grad Ф = 0,следует

 

r

r

 

 

 

A

gradФ.

(5.12)

E = −

 

t

 

 

 

31

Введенные таким образом функции координат и времени Ar(rr,t), Ф(rr,t)

rтрадиционным следованием знаков), как известно, и являются векторным r r

(A)и скалярным (Ф) потенциалами электромагнитного поля E, B . Заметим, что такое определение потенциалов является очевидным

образом неоднозначным. Действительно, поле E, B не изменяется, если r

вместо исходных A , Ф использовать A, Ф, определенные следующим образом:

r

r

=Ф

f (x, y, z,t)

.

(5.13)

A

= A + grad f (x, y, z,t), Ф

 

t

 

 

 

 

 

 

Таким образом, векторный потенциал определен с точностью до градиента произвольной функции f (x, y, z,t), а скалярный – с точностью до

производной той же функции по времени f (x, y, z,t). Это означает, что

t

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

Налагаемые на электродинамические потенциалы условия или связи называются условиями калибровки потенциалов. Из (5.13) очевидно, что

 

всегда можно f (x, y, z,t) выбрать так, что Ф (x, y, z,t)= 0 , (т.е. чисто

переменное электромагнитное поле E, B , когда f

0 , может быть в

t

r

 

принципе описано полностью одним векторным потенциалом A ).

Сформулируем теперь уравнения, определяющие электродинамические

 

 

 

 

r

 

 

 

 

потенциалы A , Ф. Для этого используем первое и третье уравнения

 

Максвелла

 

 

 

 

 

 

1

 

 

 

 

r

r

 

 

r

 

E

 

 

 

rot B = εa

 

+δ ,

(5.14)

 

µa

t

 

 

 

 

 

 

 

 

 

r

 

ρ

 

 

 

 

div E =

 

 

.

 

 

(5.15)

 

εa

 

 

 

 

 

 

 

 

 

 

Здесь εa , µa - соответственно диэлектрическая и магнитная проницаемости пространства, которые при записи (5.14), (5.15) считались постоянными, δ - плотность электрического тока, ρ - плотность электрического заряда.

32

Подставляя в уравнение (5.14) выражения (5.9) и (5.12), а в уравнение

(5.15) – (5.12), получим уравнения, определяющие A , Ф:

 

 

 

 

r

 

 

 

 

r

 

 

2 A

 

 

r

 

rot rot A + εa µa

 

 

 

= −εa µa

 

 

grad Ф + µaδ ,

(5.16)

t2

t

 

r

 

 

 

 

 

 

 

 

2

A

 

ρ

 

 

 

 

 

Ф = −div

 

 

 

.

 

 

 

t

 

εa

 

 

 

Наиболее простой вид принимает уравнения (5.16), (5.17), если использовать условие калибровки Лоренца:

divAr +εaµa Фt = 0 .

При условии (5.18) получим

 

 

 

 

 

 

 

 

 

 

2

µ

2 А

= −µ

 

 

Аε

a t2

 

δ

 

 

a

 

 

 

a

 

 

2

 

 

 

2Ф

 

ρ

 

 

Фεaµa

 

= −

 

.

t2

εa

При записи (5.19) использовано векторное тождество:

rot rot A = grad div A2

A .

(5.17)

(5.18)

(5.19)

(5.20)

Однако в случае, когда необходимо выделение содержащей разрыв на границахr источников квазистатической составляющей электрического поля

E , можно воспользоватся кулоновской калибровкой потенциалов:

 

r

 

 

 

 

 

 

(5.21)

div A = 0.

 

 

 

 

 

В этом случае уравнения (5.16), (5.17) принимают вид

 

 

2 r

 

r

 

 

Ф

 

 

 

µa

2 A

 

r

r

(5.22)

A εa

t2

= −µa

δ εa

 

= −µаδ,

 

 

 

 

 

t

 

 

2Ф = −

ρ

.

 

 

 

 

(5.23)

 

 

 

 

 

εа

 

 

 

 

33

Соответственно квазистатическая часть электрического поля Е1

определяется только через Ф, а динамическая Е2

r

- только через А :

r

r

Аr r r

r

 

E1 = − Ф,

Е2 = −

 

, E = E1

+ E 2 .

(5.24)

t

 

 

 

 

 

В дальнейшем мы будем использовать калибровку Лоренца (5.18), поскольку уравнения (5.19), (5.20) в этом случае наиболее просто связывают

электродинамические потенциалы с их источниками - δ и ρ .

5.3.Электричеcкий вектор Герца или поляризационный потенциал

Определим векторный потенциал через электрический вектор Герца Пe следующим образом:

A = εa µa Пtе .

Скалярный потенциал Ф определим через Пe так, чтобы выполнялось условие калибровки Лоренца (5.18) :

εaµa Фt = −div Ar = −εaµa t div Пrе,

тогда Ф = −div Пrе.

Подчеркнем, что указанные определения возможны только при t 0,

т.е. для чисто переменных полей. В общем случае, когда имеются статиче-

ские и стационарныеr

компоненты ЭМП, нужны все четыре составляющие

потенциалов, т.е. и Ar, и Ф.

 

r

определяются только через один вектор

При введении Пe поля E, H

Пre :

 

 

 

 

 

 

r

1

r

 

Пrе

 

 

H =

 

rot A = εarot

 

,

(5.25)

µa

 

 

 

t

 

 

34

r

r

 

 

 

 

 

2

 

e

r

 

A

grad Ф = −ε

 

µ

 

П

(5.26)

E = −

a

a

 

 

+ grad div Пе .

 

t

 

 

t2

 

 

 

Выясним физический смысл Пe . Для этого установим его источник. Обратимся к понятию поляризованности среды и вектору поляризации p .

r

r

= −qсвяз = −∫ ρсвязdV .

(5.27)

pdS

S

 

V

 

Здесь qсвяз - связанный заряд, пересекающий S при поляризации.

Уравнение (5.27) можно переписать, используя теорему Остроградского – Гаусса:

r

div pdV = −ρсвязdV .

V V

Ввиду произвольности V получаем дифференциальный эквивалент

(5.27):

div prcm = −ρсвяз.

Если действуют сторонние источники, создающие рrcm и ними, очевидно, должна существовать идентичная (5.28) связь:

div prcm = −ρcm .

(5.28)

ρст , то между

(5.29)

Дифференцируя по времени (5.29) и, используя уравнение непрерывности, имеем

r

 

 

ρcm = divδrcm .

 

 

div p

= −

(5.30)

 

t

 

 

 

 

t

 

 

 

 

Из уравнения (7.6) с учетом тождества div rot a 0 получаем

 

r

r

 

 

 

 

r

r

r

 

 

p

cm

 

 

 

 

δcm =

 

+ rot a

= δcmP

+δcm3 .

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

rP

 

 

 

r

 

 

 

 

 

 

 

 

pcm

 

 

3

 

Здесь δcm

=

 

 

 

- разомкнутая плотность тока, образующая ρcm , δст

-

 

t

 

 

 

 

 

 

 

 

 

 

замкнутая, вихревая плотность тока, для которой divδcm3

= 0.

 

Рассмотрим далее процессы, связанные только с δстP . Для этого случая

 

УМ принимают вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

35

r rot H

r rot E

 

 

r

 

 

r

 

= εa

E

+

 

Pcm

,

(5.31)

t

 

 

 

 

t

 

 

 

 

r

 

 

 

 

= −µa

H

.

 

(5.32)

t

 

 

 

 

 

 

 

Рассматриваемый случай может соответствовать, например, возбуждению ЭМП с помощью системы диполей, образованных разомкнутыми проводниками с возбужденным в них от источника сторонним током.

Подставим в (5.31) выражения H , E (5.25) и (5.26):

 

 

 

 

 

 

 

Пr

е

 

 

 

 

 

 

 

2 Пr е

r

 

 

Рr

ε

a

rot rot

 

 

 

= ε

a

 

 

ε

a

µ

a

 

 

+ grad div П

е

+

ст

.

 

 

 

 

2

 

 

 

t

 

 

 

 

 

 

 

t

 

 

 

t

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

Используя в левой части записанного уравнения тождество rot rot Пе = grad div Пе 2 Пе , получаем

 

 

 

 

 

2

r

е

 

r

 

 

2

r

е

 

П

 

pcm

 

 

 

П

 

+εа µа

 

 

 

 

 

= 0 .

 

 

 

t2

 

εa

t

 

 

 

 

 

 

 

 

В случае, который мы рассматриваем t 0 , и поэтому окончательно имеем

 

 

r е

 

r

 

 

r

 

 

 

 

2

 

2 П

е

 

pcm

 

 

 

 

 

П

εа µа

 

 

= −

 

.

 

(5.33)

 

t

 

εa

 

 

 

 

 

 

 

 

 

 

Полученное волновое уравнение для

Пe

(5.33) по структуре не отлича-

 

 

 

 

r

 

r

 

 

 

 

ется от уравнения для E,

H , A, Ф. Источником же Пe , как следует из (5.33),

является сторонняя поляризация pcm , откуда и название этого вектора – поляризационный потенциал. Простая связь Пe и pcm весьма полезна при решении конкретных задач. Например, если prcm = Z0 pcm (т.е. prcm имеет только Z – компоненту) из (7.9) очевидно, что для описания возбуждаемого ЭМП достаточно использовать только одну компоненту Пе Z0 Пze . Т.е. в этом случае векторная задача сводится к скалярной! Заметим, что в этом случае

ЭМП имеет только пять компонент: в H отсутствует в соответствии с (5.25) продольная компонента Hz. Такое поле можно назвать поперечно – магнитным.

36

5.4. Фиктивные магнитные точки и заряды. Перестановочная двойственность УМ.

Магнитный вектор Герца

Обратимся теперь к замкнутым токам. Как известно из решения задач магнитостатики, магнитное поле рамки с током i соответствует полю маг-

нитного диполя, т.е. магнитный момент M замкнутого тока определяется как

Mr = iµa Snr,

где Snr - вектор-площадка, охватываемая током i. С другой стороны, можно представить, что

Mr = iµa Snr = rlqm ,

где qm – магнитные заряды диполя длиной l.

В свою очередь, im = dqdtm , где im – разомкнутый магнитный ток.

Таким образом, замкнутый ток i создает поле, подобное тому, которое возбуждает разомкнутый магнитный ток im. Такие эквивалентные qm и im называются фиктивными магнитными зарядами и токами. В ряде случаев такая замена эффективна, поскольку позволяет не решать непосредственно поставленную задачу, а воспользоваться уже решенной с соответствующими переобозначениями.

При введении фиктивных магнитных токов и зарядов УМ переписываются следующим образом:

r

 

 

r

 

 

 

 

E

 

 

 

 

rot H =εa

,

 

(5.34)

t

 

 

 

 

 

 

r

 

 

Hr

r

 

rotE

= −µa

 

 

 

 

δm ,

(5.35)

 

t

 

 

 

 

 

 

 

r

 

 

 

 

 

 

(5.36)

div D = 0 ,

 

 

 

 

 

r

= ρm .

 

 

 

div B

 

 

(5.37)

Очевидно от системы (5.34)...(5.37) можно перейти к исходным УМ путем следующих переобозначений:

r r

r r

r

ρm → −ρ .

E H ,

H E, εa → −µa ,

µa → −εa , δm → −δ ,

37

Это свойство называется перестановочной двойственностью УМ. Это свойство позволяет, как указывалось выше, не решать ряд задач непосредственно, а воспользоваться переобозначениями в решении, относящемся к “двойственной” задаче. В частности, очевидно, что задачу возбуждения ЭМП замкнутыми токами решать нет смысла как самостоятельную, т.к. она формулируется как задача возбуждения ЭМП разомкнутыми магнитными токами, т.е. в виде двойственной уже решенной задаче возбуждения ЭМП разомкнутыми электрическими токами.

Для разомкнутых магнитных токов, как и в электрическом случае, имеет место связь

rP

Jrcm

 

 

 

 

 

 

δm =

 

,

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

где Jrcm - намагниченность среды, создаваемая δmP .

 

Второе УМ при этом записывается в виде

 

r

 

 

r

 

r

 

 

 

H

 

Jcm

 

 

rot E

= −µa

.

(5.38)

t

 

 

 

 

 

t

 

ца Пr тПо: аналогии с “электрической” задачей введем магнитный вектор Гер-

r

 

 

 

 

Пr т

 

 

E

= −µ

a

rot

 

 

,

(5.39)

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

2 r т

 

r

 

H = −εa µa

П

+ grad div Пт .

(5.40)

 

t 2

 

 

 

 

 

 

 

 

Подставляя (5.39), (5.40) в (5.38), получаем при t 0 следующее вол-

новое уравнение для Пr т :

 

2

r

т

 

2 Пr т

 

Jcm

 

 

 

 

П

 

εа µа

 

= −

 

.

(5.41)

 

 

t 2

µa

 

 

 

 

 

 

 

 

Структура уравнения (5.41) полностью идентична (5.33), что непосредственно может быть использовано при решении двойственных задач путем переобозначений.

38

Опять можно отметить, что если специфика задачи такова, что при Jrcm = Zr0 Jcm , достаточно использовать только одну компоненту Пr т =Z0 Пzm . Компонентой Пzm определяются в соответствии с (5.39) ЭМП, в которых отсутствует компонента Ez. Такие поля называются поперечно-электрическими.

5.5. Граничные условия для Пzm и Пze на идеально

проводящих продольных и поперечных поверхностях

Для большой группы задач, относящихся к регулярным волноводам и резонаторам, ЭМП может быть представлено в виде суперпозиции попереч- но-электрических (ТЕ) и поперечно-магнитных (ТМ) полей (Самарский А.А., Тихонов А.Н. «О представлении полей в волноводе в виде суммы полей ТЕ и ТМ, ЖТФ 1948 – т.18, №7 – с.959-970»). Поскольку эти поля выражаются со-

ответственно через Пzm и Пze, важнейшее значение приобретает формули-

ровка граничных значений (ГУ) для этих функций на идеально проводящих продольных и поперечных стенках указанных устройств.

Введем обобщенно-цилиндрическую систему координат q1,q2 ,Z , где q1,q2 - обобщенные криволинейные в общем случае координаты в попереч-

ной к Z плоскости Q3.

Координатные поверхности в этой системе определяются следующим образом :

Q1 : q1 = Const; Q2 : q2 = Const; Q3 : Z = Const .

Соответственно направляющие вектора системы определяются как

er1 dqr1 Q1, er2 dqr2 Q 2 , Z0 .

ГУ для Пze .

Будем исходить из того, что тангенциальная составляющая E на идеально проводящей поверхности S равна нулю. Таким образом,

Ez

Q , Q

2

= 0, Eq , q

2

Q

= 0.

 

1

1

 

3

Запишем E , выраженное через Z0 Пze :

r

 

µ

 

2 Пe

r

 

 

Пe

E = −ε

a

a

 

z

Z

0

+ grad

z .

 

 

t

2

 

 

 

 

 

 

 

 

 

 

 

z

39

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полагая процессы гармоническими

 

 

 

 

 

= jω , имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

r

 

 

 

 

 

Пe

k 2

=ω2ε

 

µ

 

 

,

 

 

 

 

 

 

E = k 2Z

0

Пe + grad

 

 

z ,

a

a

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пe

r 1

 

 

 

Пe

 

r

1

 

 

 

 

Пe

r

 

2 Пe

 

grad

 

 

z

= e

 

 

 

 

 

z

+ e

 

 

 

 

 

 

 

 

 

z

+ Z

 

 

z

,

 

 

 

 

q

 

2 h

 

q

 

 

 

 

 

 

z

1 h

 

 

 

z

 

 

 

 

2

 

 

z

 

0

z

2

 

 

 

 

 

1

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

h1, h2 – соответственно метрические коэффициенты Ламэ обобщенных

координат q1, q2.

Таким образом, имеем

 

 

 

 

 

 

2

 

Пe

 

 

 

 

 

 

 

 

 

E

z

Q1

, Q2

= k 2

+

 

 

 

Q1 , Q2

= 0 .

 

2

 

 

 

z

 

z

 

 

 

 

 

 

 

 

 

 

 

 

Это условие должно выполняться при любых Z, q1 или Z, q2, следовательно, необходимо потребовать

Пze

 

Q , Q

= 0 .

(5.42)

 

 

 

1

2

 

Иначе говоря, на идеально проводящих боковых стенках Q1, Q2 Пze об-

ращается в ноль.

Условие на поперечной стенке Q3 имеет вид

 

 

 

 

 

 

 

 

1

 

 

 

 

 

Пe

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Eq1 , q2

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

Q3 = 0.

 

 

Q3

= h

 

 

 

q

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1, 2

 

 

 

 

1, 2

 

 

 

 

 

Таким образом, если

 

 

0 , получаем

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1, 2

 

 

 

 

 

Пe

 

 

 

 

=

Пe

 

 

 

= 0 .

 

 

 

 

 

 

(5.43)

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

Q

 

rz

Q

 

 

 

 

 

 

z

 

3

 

n

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ГУ для Пzm .

 

 

 

 

 

 

 

 

 

 

 

 

 

v

Будем исходить из ГУ для нормальной составляющей H на идеально

проводящей поверхности:

 

 

 

 

 

 

 

Hq

 

Q

= 0, Hq

2

 

Q

= 0, H z

 

Q

= 0 .

 

 

 

1

 

 

1

 

 

 

 

 

 

2

 

 

 

 

 

 

 

3

40

Для составляющей Hz имеем

 

 

 

 

 

2

 

Пm

 

 

 

 

 

 

H

z

Q3

= k 2

+

 

 

 

= 0 ,

 

2

 

 

 

z

 

z

Q3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

Откуда ввиду того, что

k 2

+

 

 

 

 

 

 

 

 

0 , получаем

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пm

 

 

Q3

= 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для составляющих Hq

, Hq

2

 

можно записать

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

Пm

 

 

 

 

 

 

 

 

 

 

 

 

1

 

Пm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

= 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q1

 

 

 

Q1

 

h

q

 

 

z

 

 

Q1

 

 

z h

 

q

 

 

 

Q1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

Пm

 

 

 

 

 

 

 

 

 

 

 

 

1 Пm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

z

 

 

 

= 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q2

 

 

 

Q2

 

 

 

h

 

 

q

2

 

z

 

 

 

Q2

 

 

 

z h

 

q

2

 

 

 

 

 

Q2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

Из этих условий при

 

 

 

 

 

 

0 следует

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 Пzm

 

 

 

 

 

 

 

 

 

= 0

.или

 

Пzm

 

 

 

 

Sбок = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q , Q

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h1, 2

 

 

q1, 2

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.44)

(5.45)

ГУ (5.42)..(5.45) вместе с уравнениями (5.33) и (5.41) формулируют краевую задачу для определения ЭМП в регулярных волноводах и резонаторах.

41

ЧАСТЬ 2

НАПРАВЛЯЮЩИЕ СИСТЕМЫ ИЛИ ВОЛНОВОДЫ

ГЛАВА VI

ОБЩАЯ ТЕОРИЯ РЕГУЛЯРНЫХ ВОЛНОВОДОВ

6.1. Типы направляющих систем

Для минимизации потерь в устройствах СВЧ используются направляющие системы или волноводы (ВВ). Волны в ВВ занимают конечное поперечное сечение, они неоднородны в этом сечении и не рассеиваются в окружающее пространство. Назначение ВВ: линии связи; связь элементов СВЧ устройств; электродинамические системы в электронных приборах СВЧ; ускорители; системы СВЧ обработки в машиностроении, деревообработке, пищевой промышленности, сельском хозяйстве.

Рассмотрим некоторые типы ВВ. Вообще говоря, ВВ можно разделить на «открытые», в которых замкнутый внешний экран отсутствует, а поле концентрируется вблизи направляющих структур, и «закрытые», в которых в той или иной форме присутствует экран, ограничивающий поперечное распределение поля. Конечно, это разграничение достаточно условно: в «открытых» ВВ может быть поставлен экран, в «закрытых» - экран в определенных диапазонах может быть «прозрачным».

Начнем с систем, условно относимых к «открытым». На рис.6.1. изображена двухпроводная линия.

Рис. 6.1. Рис. 6.2.

При h << λ ввиду противофазности токов в проводах линии излучение в окружающее пространство весьма мало. На рис. 6.2 показана симметричная полосковая линия: 1 – токонесущая полоска, 2 – экран.

Несимметричная полосковая линия или микрополосковая изображена на рис.6.3. Здесь 1 – токонесущая полоска, 2 – диэлектрическая подложка, 3 – экран. Все поперечные размеры весьма малы по сравнению с длиной волны, излучение во внешнее пространство практически отсутствует. Указанные структуры нашли очень широкое распространение в микромощной технике.

42

1

2

3

Рис. 6.3.

Нетрудно увидеть развитие структуры микрополосок: экранированная микрополосковая линия (рис.6.4), щелевая линия (рис.6.5), копланарная линия (рис.6.6), полосково – щелевая экранированная линия (рис.6.7),

2

Рис. 6.4. Рис. 6.5.

Рис. 6.6. Рис. 6.7.

z

Рис. 6.8. Рис. 6.9.

двущелевая линия (рис.6.8). В коротковолновом диапазоне (вплоть до светового) используются диэлектрические волноводы (рис.6.9), на рабочих часто-

Рис. 6.10.

43

Z

Рис. 6.11.

тах которого на границе диэлектрика имеет место полное внутреннее отражение волн.

К открытым относятся и лучевые ВВ: зеркальный ВВ (рис. 6.10), линзовая линия (рис.6.11).

На рис.6.12 показана экранированная двухпроводная линия – «закрытая» система. Естественное ее развитие – коаксиальная линия (рис.6.13), весьма широко распространенная. Дальнейшее развитие ВВ большой мощности – полые ВВ.

Рис. 6.12. Рис. 6.13.

На рис.6.14 изображен ВВ круглого сечения, на рис.6.15 – прямоуголь-

ного.

Рис. 6.14. Рис. 6.15.

Находят применение и волноводы специального сечения: эллиптические (рис.6.16), ребристые: П – (рис. 6.17) и Н – образного сечения (рис.6.18).

Рис. 6.16. Рис. 6.17.

44

a

z

 

h

Рис. 6.18. Рис. 6.19.

В СВЧ электронных приборах с черенковским излучением (лампы бегущей и обратной волн) используются замедляющие системы.

На рис.6.19 изображена широко применяемая в приборах средней мощности спиральная замедляющая система. Ее замедление можно оценить

приближенно чисто геометрически:

Vф

= h

πa

, где V - фазовая скорость

 

c

 

ф

распространения основной пространственной гармоники вдоль оси системы Z, c – скорость электромагнитных волн в пустоте. В мощных приборах используются системы типа гофрированных или диафрагмированных волноводов.

Такие же системы используются и в ускорителях заряженных частиц.

6.2. Постановка и схема решения волноводных задач (регулярные ВВ)

ВВ будем называть регулярным в случае, когда конфигурация поперечного сечения ВВ S и свойства заполняющей ВВ среды εa, µa не за-

висят от координаты Z, отождествляемой с осью волновода, т.е.

s , ε&a, µ&a f (Z ) , но может быть εa ,µa = f (q1,q2 ), где q1, q2 – обобщен-

ные координаты в S .

Как указывалось выше, все поля в регулярных волноводах могут быть представлены как суперпозиция волн ТМ (Е) и ТЕ (Н). В свою очередь, ТМ

или Е – волны определяются функцией Пzl , ТЕ или Н – функцией Пzm . Таким образом, совокупность Пzl , Пzm - полностью определяет волны в ВВ.

Введем следующие условия и упрощения.

1. Рассматриваем собственные волны ВВ, которые существуют вне области, занятой источниками, т.е. в волновых уравнениях для Пzl, m Pcm = 0 и

Jcm = 0 .

2.Стенки волновода можно считать идеально проводящими, т.е. на них справедливы ГУ (5.42) и (5.45).

3.Электромагнитные процессы гармонические, т.е. можно использовать метод комплексных амплитуд.

45

При перечисленных условиях волноводные задачи сводятся к следующим краевым задачам для Пzl, m .

2

& l, m

+ K

2 & l, m

= 0,

K

2

= ω

2

εa µa ,

(6.1)

 

Пz

 

Пz

 

 

 

 

 

 

 

 

 

 

 

 

&

&

 

П&zl Sбок = 0 ,

 

 

 

 

 

 

 

(6.2)

& m

 

 

 

 

 

 

 

 

 

 

 

Пz

Sбок = 0.

 

 

 

 

 

 

 

(6.3)

 

r

 

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

 

 

 

Воспользуемся условием регулярности волновода: +Z и -Z направления равноправны, а постоянная распространения волн вдоль Z (Г) постоянна. То-

гда представим искомые решения для П&zl,m в форме

& l,m

&l,m

ψ

l,m

(q1,

q2 )

e

m jΓ Z

.

(6.4)

Пz

= A

 

 

Здесь A&l,m - комплексная амплитуда волны, ψ l,m (q1,q2 ) - функция се-

чения или мембранная функция, Г – постоянная распространения; верхний знак соответствует попутной (относительно z) волне, нижний – встречной. Решение записано в обобщенно – цилиндрической системе координат q1, q2, Z.

Представим

 

2

 

&

l,m

 

2

&

l

,m

 

2

&

l,m

 

 

 

 

 

 

 

 

 

 

 

 

Пz

 

= Пz

 

 

+

Z 2

Пz

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

h П

 

 

h

 

П

 

 

2 П&zl,m =

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

& l,m

 

 

 

1

 

& l,m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

+

 

 

 

 

 

,

 

h h

 

q

 

h

 

 

 

q

q

h

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

1

 

1

 

 

 

 

1

2

2

2

 

 

h1(q1, q2),

h2(q1, q2) – метрические коэффициенты Ламэ для q1, q2.

 

 

 

 

 

2

 

& l,m

 

 

 

 

 

2 & l,m

 

 

 

 

 

 

 

 

 

 

 

 

Учтем, что

 

 

 

Пz

 

 

 

= −Г

Пz

 

. Тогда (6.1) приобретает вид

Z 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

 

 

 

2 Г& 2 )П&Zl,m = 0.

 

 

 

 

 

 

 

 

(6.5)

2 П&Zl,m + (K

 

 

 

 

 

 

 

 

Учитывая структуру решения (6.4), сокращаем в (6.5) общие множите-

ли A&l,mem jΓZ . В результате исходная трехмерная краевая задача (6.1)..(6.3) редуцируется к двухмерной:

2 l,m

2 l,m

2

& 2

&

2

,

(6.6)

ψ

+ ψ

= 0, = K

Г

 

46

ψ l(l)= 0,

(6.7)

ψrm

(l)= 0 .

(6.8)

n

 

 

Здесь l - контур поперечного сечения волновода.

Поставленная краевая задача (6.6)..(6.8) известна как задача Штурма-

Лиувилля. Ее содержание можно сформулировать так: найти значения 2 , при которых существует нетривиальные (т.е. ненулевые в данном случае) решения (6.6) при ГУ первого рода (6.7) или ГУ второго рода (6.8) на контуре сечения волновода l. Значение v , при которых имеются нетривиальные

решения задачи Штурма-Лиувилля называются собственными значениями, а отвечающие им функции ψv (q1,q2 ) называются собственными функциями.

Каждому собственному значению v соответствует одна собственная функция ψv . Если это не так, то такой случай называется вырожденным. Собственные значения v можно расположить в виде бесконечного упорядоченного ряда:

0 < 1 < 2 <... < v <...

Собственные значения и собственные функции задачи ШтурмаЛиувилля обладают рядом специальных свойств. Остановимся на наиболее важных из них.

1. v – действительные числа. Для доказательства этого воспользуемся первой теоремой Грина:

Ф ψdS + ψ 2ФdS = ψ Фnr dl.

(6.9)

S

S

l

 

Положим Ф =ψ =ψvl,m , причем в соответствии с (6.6) имеет место равенство

2ψvl,m + v2ψvl,m = 0 .

Умножая это равенство скалярно на ψvl,m , получаем

2

 

ψ

v

 

2 = −ψ

2ψ

v

.

(6.10)

 

 

v

 

 

 

v

 

 

 

Учтем также, что в связи с ГУ (6.7) или (6.8) имеем

47

ψv

ψrv

dl

= 0 .

 

(6.11)

l

 

 

n

 

 

 

 

 

 

Учитывая (6.10) и (6.11) из (6.9) (при Ф =ψ =ψvl,m ), получаем

 

2

 

 

 

ψv

 

2 dS

 

 

 

 

 

 

=

S

 

 

 

 

.

(6.12)

 

 

 

 

 

 

 

 

 

v

 

 

 

ψv

 

2 dS

 

 

 

 

 

 

S

Из полученного результата (6.12) непосредственно следует исходное утверждение: справа стоит вещественное, строго положительное число, следовательно, и v - действительное число.

2. Собственные функции ψv (q1,q2 ) взаимно ортогональны на S (l). Докажем это. Пусть известны две собственные функции ψi и ψj, отвечающие собственным значениям i2 и 2j , причем i2 2j . Запишем уравнение (6.6)

сначала для ψi, затем для ψj. Первое уравнение умножим скалярно на ψj, второе – на ψi и из первого результата вычтем второй. В результате получим

ψ

j

2ψ

i

ψ

2ψ

j

= ( 2 2 )ψ ψ

.

 

 

 

 

 

 

 

 

 

 

i

 

 

 

j

 

i

i j

 

 

 

 

 

 

Проинтегрируем получившееся по S (l):

 

 

 

 

( j

2ψ

i

 

i

2ψ

j )

dS

 

(

j

i )

i j

dS

 

.п

(6.13)

 

 

ψ

 

ψ

 

 

= 2 2

ψ ψ

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

Применяя к левой части уравнения (6.13) дважды первую теорему Гри-

на (6.9) (сначала Ф =ψi , ψ =ψ j , затем Ф =ψ j , ψ =ψi ), получаем

 

ψ j

ψri

ψi

ψrj

dl = ( 2j

i2 )

ψiψ j dS .

 

 

(6.14)

l

 

n

 

 

 

n

 

 

 

 

 

S

 

 

 

 

 

В силу ГУ (6.7) или (6.8) левая часть в (6.14) равна нулю (заметим, что

этот результат получится и при ненулевых ГУ).

 

 

 

 

Поэтому при 2j i2

из (6.15) получаем искомый результат

 

 

 

ψi (q1,

q2 ) ψ j (q1,

q2 )dS = 0

при

ij.

 

 

 

 

S (l)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

48

e,m v

3. Собственные функции ψve,m (q1, q2 ) составляют полную на S (l) систему функций, т.е. нельзя найти такую новую функцию ϕ (q1, q2 ), которая была бы ортогональна на S (l) всем функциям ψ . Это означает, что любая функция ϕ (q1, q2 ) может быть представлена на S (l) в виде разложения

вбазисе {ψv (q1, q2 )}.

6.3.Общие свойства электрических (Е) волн

врегулярных волноводах

Весь набор Е-волн определяется системой решений {П& zve } краевой задачи (6.1), (6.2). Т.е. для любой конкретной задачи

& e

& e

&e e

&V Z

 

&

&2

2

e

,

Пz

= Пzv , где Пzv = Avψv e

 

Гv = k

v .

v=1

Краевая задача для ψv2 , v2 сформулирована в предыдущем разделе. Рассмотрим компоненты полей Е-волн:

r

 

r

 

& e

 

& e

 

 

 

 

 

 

&

&2

 

& e

Пzv

 

Пzv

 

& e

&2

2

 

 

Ev = k

Z0 Пzv + grad

 

,

 

= − v

Пzv ,

k

=ω εaµa

 

 

 

 

z

 

z

 

 

 

 

&

&

 

 

 

 

 

 

 

 

 

 

 

Hr&v = jωε&arot (Zr0 П&zve )

Продольная составляющая есть только у электрического поля Е-волны:

&

2

& e 2

&e

m &V Z e

(q1, q2 ).

Evz = v

Пzv = v

Av e

ψv

 

 

& e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отметим, что Пzv с точностью до постоянного (но размерного!) мно-

жителя v2

совпадает с E&vz .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Проанализируем далее поперечные составляющие полей.

 

r

&

&

 

& &

 

&

 

er

ψ e

 

er ψ e

 

 

&

e

e

V

Z

 

 

1

 

v

 

2

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Evt = − v grad t Пzv = − v Av e

 

 

 

 

 

q1

+

 

 

 

 

.

 

(6.15)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h1

 

h2 q2

 

 

r

r

&

r

& e

 

&

&e

 

&

er

 

ψ e

 

 

er

ψ e

&

&

 

 

V Z

 

1

 

 

v

 

2

v

 

Hv = Hvt

= jωεarot (Z0

Пzv )= jωεa Av e

 

 

 

h

 

q

h

q

.

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

1

1

 

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

49

1. Составим скалярное произведение (Er&vt , Hr&vt ):

r&

 

r&

 

 

&e

2

&

 

&

 

 

j2Г&V Z

×

 

 

 

(Evt , H vt )

= (Av

)

Гvωεa e

 

 

 

 

 

 

 

 

 

e

 

 

 

 

 

e

 

 

 

 

e

1

 

 

e

 

1

 

ψv

 

 

1 ψv

 

1 ψv

 

ψv

0 .

×

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q h q

 

 

h

 

q h

 

q

 

h

 

2

 

 

 

 

2

 

 

 

1

1

 

2

 

 

 

 

 

1

 

 

1

2

 

 

 

 

Таким образом, Er&vt Hr&vt . Но как следует из (6.15), Er&vt определяется градиентом собственной функции ψve , иначе говоря, графически поле Er&vt

совпадает с градиентом собственной функции ψve .

r

r

 

 

Поскольку же H vt Evt , линии напряженности магнитного поля явля-

ются линиями уровня ψ e (q , q

2

). Таким образом, мы получили эффективное

v

1

 

правило графического построения полей волны типа Е волновода: 1) ψ e (q , q ) является мембранной функцией E ;

2)силовые линии Errtv являются линиями градиентаz ψve (q1, q2 );

3)силовые линии Htv являются линиями уровня ψve (q1, q2 ).

2.Введем и рассчитаем волновое сопротивление Wve для волн типа Е:v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

e

 

 

 

1

 

 

 

 

 

e

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&e

V Z

 

 

 

ψv

 

+

 

 

 

 

ψv

 

 

 

 

 

 

 

 

&

 

 

 

 

 

v Av e

 

 

h

 

 

q

 

h

 

q

2

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

Evt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wv

 

=

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

H& vt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

e 2

 

 

 

1

 

 

 

 

 

e 2 2

 

 

 

 

 

 

 

 

 

 

 

 

&e

e

V Z

&

 

 

 

 

ψv

 

 

+

 

 

 

 

ψv

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jA

 

 

 

ωε

 

 

 

 

 

 

 

 

q

 

 

 

h

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

a

h

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Г&v

 

 

 

k&

 

 

 

 

v 2

 

 

0

 

 

 

 

v

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

=

 

 

 

 

1

&

 

=W

 

 

1

&

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ωεa

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ωεa

 

 

k

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

 

 

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где W 0 =

 

µ&a

- волновое сопротивление свободного пространства с

 

 

 

 

 

 

 

 

ε&a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

характеристиками εa, µa. Наиболее просто записывается Wve , когда k& - дейст-

вительное число (тогда k = ω εa µa

=

ω

=

2π ).

 

 

υ

 

λ

В этом случае, учитывая, что v - действительная величина, можно ввести

одно из важнейших понятий в теории ВВ – критическую длину волны v-го типа волн:

50