Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Конспект лекций Электродинамика и РРВ

.pdf
Скачиваний:
333
Добавлен:
10.05.2015
Размер:
3.16 Mб
Скачать

можно определить и постоянную распределения Г; для нахождения коэффициента затухания надо лишь отделить ее мнимую часть. В частности, для Е- и Н- волн волноовода с идеально проводящей оболочкой χ2 > 0. Пусть внутренняя среда является поглощающей, тогда

Г = k2 (1itg) − χ2 = 2 ik 2tg.

(16)

Разделение вещественной и мнимой частей приводит к следующим формулам:

 

 

1

 

4

 

4

 

2

 

 

2

 

2

 

Г = k

2

(

 

/ k

 

+ tg

∆ +

 

/ k

 

) ,

′′

 

 

1

 

 

 

 

 

 

 

 

 

 

 

(17)

 

 

 

4

 

4

 

2

 

2

2

 

Г = k

2

(

 

/ k

 

+ tg

 

∆ −

/ k

 

) .

При k отсюда получаются выражения

 

 

 

 

 

 

1

 

 

2

 

 

 

 

 

 

 

 

Г

= k 2 ( 1+ tg ∆ +1,

 

 

 

 

 

 

 

 

 

(18)

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

′′

 

 

 

 

2

 

 

 

 

 

 

 

 

Г

= k 2 ( 1+ tg ∆ −1,

 

 

 

 

 

 

 

 

 

 

справедливые для Т – волны.

Периодические структуры

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

Примеры одномерно–периодических или продольно–периодических

структур даны на рис.2. Структуры могут быть открытыми (а, б, в, г) и экранированными ( д, е).

Рис.2. Примеры одномерно-периодических структур

Свободные электромагнитные поля в одномерно–периодических структурах подчиняются так называемой теореме Флоке, выражающей следующее свойство комплексных амплитуд векторов Е и Н:

Еm (x, y, z + Λ) = Еm (x, y, z)eiϕ,

(19)

Hm (x, y, z + Λ) = Hm (x, y, z)eiϕ

,

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

В силу теоремы Флоке (19) нетрудно построить следующие периодические по z функции:

E (x, y, z) = Em (x, y, z)eiγz ,

(20)

H (x, y, z) = Hm (x, y, z)eiγz

,

γ = ϕ/Λ, где γ - специально веденный параметр. Периодичность записанных векторных функций следует из того, что множитель exp(iγz) компенсирует «естественный» фазовый сдвиг, выражающий согласно (19) на отрезке длинной Λ.

Функции E и H можно разложить в ряды Фурье, выразив коэффициенты разложения через соответствующие интегралы, например,

 

 

 

 

 

 

 

 

 

E (x, y, z) = En (x, y)e−ϕ(2

πn / Λ) z ,

 

(21)

где

 

 

n=−∞

 

 

 

 

 

 

zE (x, y, z)ei(2πn / Λ) z dz =

 

 

zE

 

 

 

 

E (x, y) =

1

1

 

 

(x, y, z)ei(γ+2

π/ Λ) zdz.

(22)

Λ

Λ

 

n

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

z

 

 

 

 

Точно также можно представить H . Введя множитель exp(-iγz), снова

перейдем от E и H к Εm и

Hm соответственно:

 

 

 

 

 

 

Em (x, y, z) = E n (x, y)ei(γ+2

πn / Λ) z ,

n=−∞

 

(23)

 

 

 

Hm (x, y, z) = Hm (x, y)ei(χ+2

πn / Λ) z .

n=−∞

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

тудами E0 exp(n z), H0 exp(n z) и постоянными распространения

 

Гn = γ + n

2π

,

n = 0,±1,±2,,±∞.

(24)

Λ

 

 

 

 

Эти волны, называемые пространственными гармониками, имеют следующие фазовые скорости

vф(n) = ω Гn = ω (γ + 2πn Λ)

(25)

и одну общую групповую скорость

vгр = dω n = dω dγ.

(26)

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

Периодические структуры, для которых Λ << λ, называют частыми: на расстоянии длины в однородной среде укладывается большое число пространственных периодов.

Пусть, например, рассматривается ребристая структура, однородная в направлении y (рис.3).

Рис.3. Ребристая структура

Если она частая, то можно ожидать, что по отношению к верхнему полупространству (x < 0) граница структуры x = 0 проявляет некоторые усредненные свойства, а волновой процесс, распространяющийся в направлении z, представляется главным образом нулевой гармоникой разложений (23).

В пазах структуры может существовать поле, в общих чертах, показанные на рис.3. Волновой процесс в целом можно считать Е-волной. В первом приближении поле внутри паза не изменяется в направлении z, имея характер стоячей T-волны по х, так что

E

m

= z

i2Asin k(x d),

H

m

= y

 

2A

cos k(x d)

(27)

 

 

 

0

 

 

 

0

Z

 

(0 < x <d); выполняется условие Em = 0 при x = d. На основании (27) находим следующее соотношение между комплексными амплитудами Em и Hm на плоскости x = 0:

Em (0) = ZS (Hm (0)×x0 ), ZS = i tg kd .

(28)

Поэтому, пренебрегая площадью ребер при x = 0, можно сказать, что эта граница структуры проявляет себя как поверхность, характеризуемая импедан-

сом ZS.

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

Лекция № 18. Прямоугольный металлический волновод

Среди полых волноводов наиболее распространен прямоугольный волновод – металлическая труба прямоугольного поперечного сечения (рис.1).

Рис. 1. Прямоугольный металлический волновод

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

В прямоугольном волноводе с идеально проводящей оболочкой могут существовать только волны классов E и H. Рассматривая E–волны, мы должны решить первую краевую задачу для прямоугольного контура

2 E +χ2E = 0,

E = 0 на L .

(1)

z

z

z

 

 

В декартовой системе координат уравнение Гельмгольца имеет вид:

2E

+

2E

2

E = 0.

(2)

z

 

z + χ

x2 y2

z

 

 

 

Применение метода разделения переменных начинается с предположе-

ния, что неизвестное решение Ez можно представить в виде произведения

функций разных координат: Ez (x, y)= X (x)Y (y). Подстановка этого представ-

ления в (2) дает:

 

 

 

 

 

d 2 X Y + d 2Y

X + χ2 XY = 0,

 

dx2

dy2

 

 

 

 

а после деления на XY получаем уравнение:

 

 

1

d 2 X +

1

d 2Y = −χ2

,

(3)

 

X

 

 

dx2

Y dy2

 

 

где слагаемые слева – функции разных аргументов. Они, таким образом, независимы, а, следовательно, каждое слагаемое равно константе; обозначив эти

константы −χ2x и −χ2y , получаем вместо (3) следующие два обыкновенных дифференциальных уравнения:

d 2 X + χ2 X = 0,

d 2Y 2Y = 0,

(4)

dx2

x

dy2

y

 

 

 

 

которые эквивалентны уравнению (3) при:

 

χ2x + χ2y 2 .

(5)

Общие решения уравнений (4) можно выразить, как известно, в тригонометрической форме следующим образом:

X = Acosχx x + Bsin χx x, Y = C cosχy y + Dsin χy y.

(6)

Здесь введен ряд неопределенных констант, неопределенными являются

также χx и χy.

Граничное условие Ez = 0 на L означает следующие требования:

Ez = 0

(7)

при x = 0, y = 0; x = a, y = b. Мы должны иметь, в частности, Ez (0, y) = 0 . Под-

ставив в (6) x = 0, видим, что это возможно только при А = 0. Потребовав далее выполнения равенства Ez (x,0)= 0 , точно также убеждаемся, что С = 0, то есть:

Ez = N sin χx x sin χy y,

(8)

где N = BD – неопределенная константа. Остается наложить на (8) условия Ez (a, y) =0 и Ez (x,b)= 0 . Они выполняются при χxa = mπ (m = 0,1,2…) и χyb = nπ (n = 0,1,2…), то есть:

χx =

mπ

,

χy =

nπ

,

(9)

a

b

 

 

 

 

 

где m = 1,2,3… и n = 1,2,3…(нулевые значения m и n исключаются, так как при этом Ez = 0 ).

Итак, решение задачи (1) для прямоугольного волновода дает согласно

(8), (9).

 

Emn = Emn sin mπsin nπ,

 

 

 

z

 

0

a

b

 

 

 

 

 

 

 

 

 

2

mπ 2

nπ 2

 

 

χmn

=

 

+

, m =1,2,3...; n =1,2,3...,

 

 

 

 

 

a

 

b

 

(10)

 

 

 

 

 

 

 

 

 

 

где Emn

- неопределенные коэффициенты. Зная эти собственные функции Emn

и

0

 

 

 

 

 

z

 

соответствующие им собственные значения χ2mn , путем подстановки (10) в соотношения:

Emn = z

Emn i(Г

mn

χ

mn

)2

 

Emn emn z

,

m

0

z

 

 

z

(11)

Hmmn = −(iωε0ε χmn2 )

(z0 × Ezmn )emn z ,

 

Выразим электромагнитные поля Е–волн:

E mm n = E0m n {z 0 sin (χmx n x )sin (χmy n y )

 

 

 

 

i

Г

2

 

[x 0 χmx n cos (χmx n x )sin (χ y y )+

 

 

 

 

 

 

 

 

 

 

 

 

m n

 

 

 

 

 

 

 

 

 

 

 

 

χ

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+y 0 χmy n sin (χmx n x )cos (χmy n y )]}e m n z ,

 

 

 

 

 

 

(12)

 

 

 

 

 

 

 

m n

 

 

 

 

E 0m n Гm n

 

 

 

 

m n

 

 

 

m n

 

 

 

 

 

m n

 

 

 

 

 

 

 

 

 

 

H m

= i

 

 

 

 

 

 

 

 

 

[x 0 χ y

sin (χ x

 

x )cos (χ y

 

y )

 

 

 

 

 

 

 

 

 

E

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z m n

 

 

χm n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y 0 χmx n cos (χmx n x )sin (χmy n y )]}e m n z ,

 

 

 

 

 

 

 

 

где

 

 

Гmn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ZmnE =

 

= Z 1(fкрmn

f )2 = Z 1(λ λкрmn )2 ,

 

Z0 =120π

 

µ ε,

(13)

 

 

 

ωε

ε

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Гmn = k

 

 

1(fкрmn

 

f )2

= k

1(λ λкрmn )2 ,

 

 

 

 

 

 

(14)

причем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mn

 

 

 

 

c

 

 

 

mπ

2

nπ

2

 

 

mn

 

 

 

 

 

mπ 2

nπ 2

 

 

 

fкр

=

 

 

 

 

 

 

 

 

 

 

 

+

 

 

,

 

λкр

= 2

 

 

 

 

 

+

 

 

 

.

(15)

 

 

2π εµ

 

 

 

 

 

 

 

 

 

 

 

 

a

 

b

 

 

 

 

 

 

 

 

 

a

 

b

 

 

Анализ H–волн требует решения краевой задачи:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 H +χ2H = 0,

 

 

∂H ∂ν = 0 на L .

 

 

 

 

 

 

(16)

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

z

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Представление искомого решения в виде Hz (x, y)= X (x)Y (y) приводит

к соотношению:

 

Acos(χ

 

x)+ Bsin (χ

 

x)

 

 

(

 

 

 

 

)

 

 

(

 

 

 

)

 

 

H

z

(x, y)=

x

x

C cos

χ

y

y

+ Dsin

χ

y

y

.

(17)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На решение (17) следует наложить условие ∂Hz

∂ν = 0

на L , что озна-

чает:

 

 

 

 

 

 

 

 

 

 

 

 

 

∂Hz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 0

 

при x =0, x = a,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(18)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂Hz

= 0

 

при

 

y = 0, y = b.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дифференцируя Hz (x, y) по х, а затем по у и требуя обращения в нyль

соответствующих производных при х = 0 и у = 0 , находим , что В = 0 и D = 0, то есть

Hz (x, y)= N cos(χx x) cos(χy y),

(19)

где N – неопределенная константа.

Налагая теперь условие обращения в нуль тех же производных при х = а и у = b соответственно, приходим к прежним выражениям (9).

Поэтому:

H zmn = H 0mn cos (χmnx x) cos (χmny y ),

2 mπ 2

nπ 2

mn

)

2

mn

)

2

,

(20)

χmn =

a

 

+

b

 

= (χx

 

+ (χy

 

 

 

 

 

 

 

 

 

 

 

 

 

где m = (0), 1, 2, …, n = (0), 1, 2, …, Н0mn – неопределенные коэффициенты. Нулевые значения m и n можно брать лишь при сочетании с ненулевыми (поэтому

нули взяты в скобки). Дело в том, что Hz00 есть константа, причем χ002 = 0; это решение не относится к классу Н-волн. Подставляя собственные функции Hzmn

и собственные значения χ2mn в соотношения:

Emn = −(iωµµ

0

χ2

) ( H mn ×z

)emn z ,

m

mn

z

0

 

Hmnm = z0 Hzmn i(Гmn χ2mn ) Hzmn emn z ,

получаем Н–волны:

 

E mm n = iZ mHn H 0m n

Г

m n

[x 0 χmy n cos

(χmx n x )sin (χmy n y )

 

χ

2

 

 

m n

 

 

 

 

y 0 χmx n sin (χmx n x )cos (χmy n y )]e m n z ,

 

 

 

 

 

Г

mn

 

H mmn

= H 0mn z0 cos (χmnx x )cos (χmny y )+ i

 

 

2

 

 

 

 

 

χmn

x0χmnx sin (χmnx x )cos (χy y )+ y 0χmny cos (χmnx x )sin (χmny y ) emn z ,

(21)

(22)

где

Z H

=

ωµµ0 =

 

Z

 

 

=

Z

 

 

,

(23)

 

 

 

 

 

 

 

 

mn

 

Гmn

1

mn

f )

2

 

mn

)

2

 

 

 

 

 

 

(fкр

 

 

1(λ λкр

 

 

 

а постоянные распространения Гmn и критические частоты

fкрmn по-прежнему

выражаются формулами (14), (15).

Полученные решения показывают, что прямоугольному волноводу свойственно бесконечное множество свободных электромагнитных полей волн классов Е и Н, которые определяются выбором чисел m и n в (10) или, соответственно (20). Говорят, что выбранное решение дает тип поля, или тип волны

Еmn (в классе Е) и Нmn (в классе Н).

При фиксированной частоте f только для некоторых достаточно малых m и n будет выполнено условие f > fкp (15). Поэтому лишь для конечного числа типов поля постоянные распространения Гmn (14) окажутся вещественными. Эти типы поля имеют характер распространяющихся волн, которые переносят энергию. Все остальные типы поля, составляющие бесконечное множество, энергии не переносят и экспоненциально затухают. Передача энергии невозможна, если f < fкр для всех m и n.

Особую роль играет тип волны с наименьшей критической частотой min fкрmn . Если a > b (см. рис.1), то

min fкрmn = fкр10 =

c

 

1

.

(24)

 

 

 

εµ 2a

 

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

волновода. Обычно соблюдается условие f > fкр10 при f < fкрmn для остальных

типов волн. В этом случае в волноводе перенос энергии осуществляется только одной основной волной.

Рассмотрим строение электромагнитных полей. В классе Е простейшим является тип поля Е11 (рис. 2).

Характер этой картины можно определить из следующих соотношений. Поскольку речь идет о волне класса Е, вектор Н вообще не имеет продольной компоненты: все магнитные силовые линии лежат в поперечных плоскостях. Известно, что они всегда должны быть замкнутыми кривыми (среда однородная, так что линии В и Н одинаковы). Из соображений симметрии ясно, что центром семейства замкнутых магнитных силовых линий должна быть средняя точка сечения. Заметим, что в этой же точке лежит максимум продольного тока смещения. В направлении z максимумы Еz и дD/дt сдвинуты на Λ/4.

Рис. 2. Распределение поля волны Е11

На рис.2 показаны два продольных сечения «мгновенного снимка» волны Е11: x = a/2 и y = b/2. Отмечен отрезок структуры длиной Λ/4. Видно, что длина волны в волноводе Λ есть пространственный период поля.

При анализе любых структур в классе Е рассмотренный тип поля Е11 служит «элементарной ячейкой».

Отметим следующее: во–первых, все плоские границы между ячейками могут быть заменены идеально проводящими плоскостями, это не нарушает структуру поля; во–вторых, направления силовых линий в соседних ячейках согласованы таким образом, что на их границах тангенциальные компоненты Е и Н непрерывны; в–третьих, во всех случаях, когда m и n больше единицы, появляются замкнутые электрические силовые линии; семейства замкнутых электроических и магнитных силовых линий как бы сцеплены подобно звеньям це-

пи. В качестве примера на рис.3 представлена структура поля Е32 (m = 3, n = 2). Можно сказать, что поперечное сечение разбито на 3х2 клеток, в каждой из которых воспроизводится поле Е11 .

Перейдем к обсуждению полей класса Н, начав с поля Н11. Оно не является простейшим, но более простые структуры, когда m = 0 или n = 0, будут рассмотрены отдельно. На рис.4 показаны два поперечных и два продольных сечения структуры поля Н11 для момента времени t = 0.

Структура поля Н11 играет роль элементарной ячейки при анализе всех более сложных структур, когда m 1, n 1. В качестве примера на рис.5 пока-

зана структура поля Н32 .

Рис. 3. Структура поля волны Е32

Рис. 4. Структура поля волны Н11

Рис. 5. Структура поля волны Н11

Все рассматривавшиеся выше Е– и Н–волны , как говорят, являются попарно вырожденными: разным собственным функциям Ezmn и Hzmn при одних и тех же m и n соответствуют разные собственные значения χ2mn . Поэтому раз-

личные по структуре полей волны Етп и Нтп имеют одинаковые постоянные распространения Гтп, а следовательно, равные фазовые скорости. Вырождение снимается при переходе к реальному металлическому волноводу. Более того, собственные волны оказываются уже гибридными: типы поля Етп и Нтп связываются в некоторые комбинации с преобладанием одного из них.

С повышением частоты или увеличения поперечных размеров волновода растет отношение Г/χ. В пределе отношение продольных компонент к поперечным стремится к нулю. Волны классов Е и Н переходят в Т–волны.