Электромагнитные волны в направляющих структурах

advertisement
МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ УКРАИНЫ
Национальный аэрокосмический университет им. Н.Е. Жуковского
«Харьковский авиационный институт»
С.Н. Барсуков
ЭЛЕКТРОМАГНИТНЫЕ ВОЛНЫ В НАПРАВЛЯЮЩИХ СТРУКТУРАХ
Учебное пособие
Харьков «ХАИ» 2007
3
С.Н. Барсуков
ЭЛЕКТРОМАГНИТНЫЕ ВОЛНЫ В НАПРАВЛЯЮЩИХ СТРУКТУРАХ
2007
4
УДК 621.372.8(075.8)
Электромагнитные волны в направляющих структурах/ С.Н. Барсуков.
– Учеб. пособие. – Харьков: Нац. аэрокосм. ун-т «Харьк. авиац. ин-т»,
2007 – 59 с.
Изложена теория электромагнитных процессов в направляющих
структурах СВЧ, в частности, прямоугольных и круглых волноводах,
коаксиальных и полосковых линиях передачи. Рассмотрены особенности и типы направляемых волн, их основные характеристики. Приведена методика энергетического анализа направляющих структур.
Для студентов факультета радиотехнических систем летательных аппаратов при усвоении части 2 лекционного курса «Электродинамика и распространение радиоволн».
Ил. 49. Табл.3. Библиогр.: 9 назв.
Рецензенты: д-р техн. наук, проф. А.В. Полярус,
канд. техн. наук, доц. Е.А. Милькевич
©Национальный аэрокосмический университет им. Н.Е. Жуковского
«Харьковский авиационный институт», 2007 г.
5
1. НАПРАВЛЯЕМЫЕ ВОЛНЫ
Направляемая волна – это волна, которая распространяется
вдоль заданного направления. Приоритетность направления обеспечивается направляющей системой. К направляющим системам относятся: металлические плоскости, металлические или диэлектрические
стержни, металлические трубки и т.д. Конструктивно эти системы
реализуются в виде волноводов, коаксиальных и полосковых линий.
1.1. Особенности направляемой волны. Волновое уравнение
Волна, которая распространяется в направляющей структуре, называется направляемой волной. Направляемая волна – неоднородная: в плоскости фронта волны ее компоненты E и H не постоянны,
поскольку изменяются по закону стоячей волны. В такой волне поверхности равных фаз и равных амплитуд не совпадают. Изменения
амплитуды волны определяются необходимостью выполнения граничных условий на поверхности, которая ограничивает эту структуру.
На рис. 1.1 показана произвольная односвязная направляющая
система, где ΔS – поперечное сечеX
ние направляющей системы, S –
S
плоскость фронта (XOY), Z – направление распространения.
Z
В общем виде комплексная амY
плитуда направляемой волны имеет
два сомножителя:
Y
Рис. 1.1
E& m = E& 0 ( x, y )e − jhz ,
(1.1)
где E& 0 ( x, y ) – амплитудный множитель, то есть функция распределения напряженности поля в плоскости поперечного сечения волны;
e − jhz
– фазовый множитель;
h – продольное волновое число.
6
Волновое уравнение (Гельмгольца) для монохроматической волны в неограниченном пространстве имеет следующий вид:
∇ 2 E& m + k 2 E& m = 0,
k=
где волновое число
2π
ω
= ω εμ =
.
ϑ
λ
В ограниченной направляющей структуре решение уравнения
имеет вид (1.1), что позволяет упростить исходное уравнение.
Разобьем пространственный лапласиан на два слагаемых, тогда
получим
∂2 &
2&
Em + k 2 E& m = 0,
∇ s Em +
2
∂z
где
2
∂2
2 ∂
– двумерный лапласиан в плоскости фронта вол+
∇S =
2
2
∂x
∂y
ны.
Вторая производная по продольной координате с учетом (1.1) заменяется в соответствии со следующим равенством:
∂2 &
Em = (− jh )2 E& m = −h 2 E& m .
2
∂z
Тогда волновое уравнение примет такой вид:
∇ 2s E& m + k 2 − h 2 E& m = 0.
(
)
2
2
Если в уравнение ввести обозначение g = k − h – поперечное волновое число, то получим
∇ 2s E& m + g 2 E& m = 0,
(1.2)
а если сократить его на экспоненту, то –
∇ 2s E& 0 ( x, y ) + g 2 E& 0 ( x, y ) = 0.
Таким образом, волновое уравнение для направляемой волны
зависит только от поперечных пространственных координат.
1.2. Взаимосвязь поперечных и продольных компонент волны
Особенность направляемой волны (1.1) позволяет выразить поперечные компоненты волны через продольные. Поэтому взаимосвязи
E& X , E&Y = f1 E& Z , H& Z ,
H& X , H& Y = f 2 E& Z , H& Z
можно получить исходя из уравнений Максвелла, записанных в комплексной форме для диэлектрической среды:
(
(
)
7
)
rot H& = j ωε E& ;
rot E& = − j ωμ H& .
Запишем компоненты ротора и приравняем составляющие векторов в левой и правой частях уравнений:
→
→ →
X 0 Y0 Z 0
rotE& =
∂ ∂
∂x ∂y
∂
;
∂z
E& x E& y E& z
Заменим производную
второе равенства:
∂E& z ∂E& y
&
= − jωμH& x ,
−
rot x E =
∂z
∂y
⎛ ∂E&
∂E& ⎞
rot y E& = −⎜⎜ z − x ⎟⎟ = − jωμH& y ,
∂z ⎠
⎝ ∂x
∂E& y ∂E& x
&
−
= − jωμH& z .
rot z E =
∂x
∂y
∂
= −( jh )
и выберем первое и
∂z
∂E& Z
&
&
jωμH X + jhEY = −
, μ → −ε
∂y
∂E& Z
&
&
jωμH Y − jhE X =
, E↔H
∂x
∂ H& Z
&
&
− j ωε E X + jh H Y = −
,
∂y
∂H& Z
&
&
.
− jωεEY − jhH X =
∂x
В связи с симметрией уравнений Максвелла вторую пару уравнений получим с использованием принципа эквивалентных замен. Из
полученных уравнений выберем пару уравнений, в которые входят
компоненты H& x , E& y . Решим выбранную систему уравнений относительно этих переменных. Для этого умножим обе части уравнений на
соответствующие коэффициенты, а затем сложим эти уравнения:
Тогда
∂E& Z
&
&
, jh
jωμH X + jhEY = −
∂y
∂H& Z
&
&
, − jωμ
jhH X + jωεEY = −
∂x
jωε ,
− jh .
∂E& Z
∂H& Z
+ jωμ
,
∂y
∂x
∂E&
∂H& Z
− k 2 H& X + h 2 H& X = − jωε Z + jh
.
∂x
∂y
− h 2 E&Y + k 2 E&Y = − jh
8
Учитывая, что
k 2 − h 2 = g 2 , получим
⎛
j ⎜
∂ E& Z
∂ H& Z
E& Y =
+ ωμ
⎜− h
∂y
∂x
g2 ⎜
⎝
⎞
⎟
⎟,
⎟
⎠
(1.3)
⎛
⎞
⎜
&
&
∂ EZ
∂HZ ⎟
j
&
⎜ − ωε
⎟.
+h
HX =−
2⎜
∂
∂
y
x
⎟⎟
g ⎜
⎝
⎠
Используя эквивалентные замены μ → −ε , E → H , имеем еще
пару уравнений:
j ⎛
∂H&
∂E& ⎞
H& Y =
⎜⎜ − h Z − ωε Z ⎟⎟ ,
∂y
∂x ⎠
g2 ⎝
(1.4)
∂H& Z
∂E&
j ⎛
+h Z
E& X = −
⎜⎜ ωμ
∂y
∂x
g2 ⎝
⎞
⎟⎟.
⎠
Эта система из четырех уравнений определяет взаимосвязь поперечных компонент направляемой волны через продольные. В свою
очередь, продольные составляющие H& z , E& z определяются из решения одномерных (скалярных) волновых уравнений
∇ 2S H& Z + g 2 H& Z = 0
∇ 2S E& Z + g 2 E& Z = 0,
с учетом граничных условий Eτ = 0, H n = 0 (краевая задача Дирихле).
1.3. Типы направляемых волн
Проанализируем особенности волнового уравнения в двухсвязных и односвязных направляющих структурах.
1. В первом варианте волновое уравнение
∇ 2s E& + g 2 E& = 0
удовлетворяется в частном случае при двух условиях:
∇ 2s E& = 0,
g = 0.
Равенство нулю плоского лапласиана соответствует статическому
полю. Действительно, из уравнения Максвелла
∂B&
&
,
rotE = −
∂t
9
∂B&
= 0 получим
а при статическом поле
∂t
rotE& = 0.
Тогда, используя тождество, имеем
rotrotE& = graddivE& − ∇ 2 E& = 0,
причем в свободном от зарядов пространстве divE& = 0, следователь2
но, ∇ E& = 0.
Вывод: этот тип волн существует в направляющих системах, по
которым можно передавать постоянный ток (двусвязные системы). В
2
2
2
таких структурах g = 0, поэтому волновое число k = g + h совпадает с волновым числом для неограниченного пространства, т.е.
k = h. В этом случае волна является поперечной ( EZ = 0, H Z = 0 ),
т.е. линии векторов E и H лежат в плоскости поперечного сечения
волны. Такая волна обозначается как ТЕМ-волна (transvers – поперечный).
В односвязных направляющих системах распространяются волны
поперечно-магнитного (TH) или поперечно-электрического (TE) типа.
2. Электрическая (E-волна), или TM (TH), – это поперечномагнитная волна, в которой отсутствует продольная магнитная составляющая H Z = 0.
3. Магнитная (H-волна), или TE, – это поперечно-электрическая
волна, в которой продольная составляющая E Z = 0 .
В таких структурах плоский лапласиан и поперечное волновое
число отличны от нуля:
∇ 2s E& ≠ 0 ,
g ≠ 0.
Постоянная распространения h отличается от волнового числа
неограниченного пространства k , то есть
h = k2 − g2 .
(1.5)
В предельном случае при
k = g = kK ,
h=0
исчезает волновой процесс. Введя понятие критической длины волны,
постоянную распространения можно представить в таком виде:
2
⎛ λ ⎞
2
2
⎟⎟ ,
h = k − k K = k 1 − ⎜⎜
(1.6)
λ
⎝ K⎠
где k =
2π
λ
, kK =
2π
λK
– текущее волновое число и его критическое
значение соответственно.
10
Критическая длина волны λ K – это такое максимальное значение
длины волны, измеренной в свободном
Волновой режим
Отсечка
пространстве, при которой прекращается
распространение волны вдоль направλ ляющей системы.
λK
В зависимости от длины волны возможны следующие режимы направляюРис. 1.2
щей структуры (рис. 1.2):
а) режим распространения: λ < λ K , h – действительное число,
E (t ) = E0 (x, y) cos(ωt − hz);
б) режим отсечки:
λ > λK , h = − jh1, E (t ) = E0 ( x, y)e− h z cosωt.
1
1.4. Характеристики дисперсных волн
Волны E- и H-типов обладают дисперсией, то есть фазовая скорость ϑφ = f (ω ) , причeм в отличие от ТЕМ-волны в этих волнах
h ≠ k , следовательно, они имеют конечное
h
k
g
Рис. 1.3
значение λ K .
К основным характеристикам дисперсных волн относятся следующие параметры:
1. Продольное волновое число
2
⎛ λ ⎞
⎟⎟ ,
h = k 1 − ⎜⎜
⎝ λK ⎠
причем в направляющей системе h < k (рис. 1.3), так как волновое
число неограниченного пространства k раскладывается на продоль2
2
2
ное h и поперечное g волновые числа, т.е. k = g + h .
2. Увеличенная длина волны в направляющей системе
λB =
2π
=
h
λ
⎛ λ
1 − ⎜⎜
⎝ λK
⎞
⎟⎟
⎠
2
,
(1.7)
где λ – длина волны в неограниченном пространстве, то есть ТЕМволны. Эффект увеличения длины волны в направляющей системе
изображен на рис. 1.4, где показана взаимосвязь длины волны ТЕМтипа λ с длиной волны в направляющей структуре λB . На основе
геометрических построений соотношение длин волн описывается
следующей формулой:
11
π
λ = λB cos α = λB cos( − ϕ ) = λB sin ϕ .
2
П ТЕМ
λ
П ТЕМ
ПЕ
ϕ
α
Z
λB
Рис. 1.4
3. Фазовая скорость, которая определяется через продольное
волновое число:
ϑφ =
ω
=
h
ω
⎛ λ
k 1 − ⎜⎜
⎝ λK
⎞
⎟
⎟
⎠
2
,
или ϑ φ =
ϑ0
⎛ λ
1 − ⎜⎜
⎝ λK
⎞
⎟⎟
⎠
2
,
(1.8)
где ϑ0 – скорость в неограниченном пространстве. Следовательно,
фазовая скорость увеличивается, т.е. волна становится «быстрой».
Эта скорость возрастает в сравнении со скоростью в неограниченном
пространстве: ϑф > ϑ0 .
На рис. 1.4 также показано изменение фазовой скорости
λ
ϑφ = в =
T
где
ϑ0 =
λ
T
λ
T sin ϕ
=
ϑ0
,
sin ϕ
.
4. Групповая скорость, которая находится из условия взаимосвязи скоростей
С2
2
ϑ Г ϑφ = С ,
где С = ϑ0 , т.е. ϑ Г =
,
ϑφ
следовательно,
⎛ λ
ϑ Г = ϑ0 1 − ⎜⎜
⎝ λK
2
⎞
⎟⎟ .
⎠
(1.9)
5. Волновое сопротивление Е-волны, которое определяется
отношением поперечных компонент:
12
ZC =
EY
.
HX
Используя уравнения представлений поперечных компонент через продольные (1.3), (1.4) при H z = 0 , имеем
⎛ λ
=
Z CE =
1 − ⎜⎜
ωε ωε
⎝ λK
h
k
2
⎞
⎟⎟ .
⎠
Поскольку волновое число и волновое сопротивление для неограниченного пространства соответственно равны:
k = ω εμ ,
ZC =
μ
,
ε
то
⎛ λ
Z CE = Z C 1 − ⎜⎜
⎝ λK
2
⎞
⎟⎟ .
⎠
(1.10)
Волновое сопротивление для H-волны
Z CH =
ωμ
h
ZC
=
⎛ λ
1 − ⎜⎜
⎝ λK
⎞
⎟⎟
⎠
2
.
(1.11)
Уменьшение поперечной составляющей Е-волны по сравнению
с ТЕМ (рис. 1.5) приводит к уменьшению волнового сопротивления
Z CE в связи с уменьшением поперечной электрической компоненты,
так как
ZC =
E
EX
ЕZ
EX
.
H
П ТЕМ
ПЕ
Z
Замечание: при увеличении длины волны
(λ>λK) Z c
направляемой волны становится
мнимым:
π
E E& X
E
E j2
Z =
= jZ C = Z C e .
H&
Рис. 1.5
При этом электрическая и магнитная компоненты волны сдвинуты на 90 o , что приводит к мнимой величине вектора Пойнтинга П . В
этом случае отсутствует перенос энергии вдоль направляющей структуры.
13
2. ПРЯМОУГОЛЬНЫЙ ВОЛНОВОД
Прямоугольный волновод является направляющей структурой
сантиметрового диапазона волн.
y
Геометрия волновода с выбранной
системой координат приведена на рис. 2.1.
z
Поперечное сечение волновода лежит в
b
плоскости XOY, размеры поперечного се- x
чения – a × b , ось Z соответствует проa
дольной оси волновода.
Рис. 2.1
Условия рассмотрения:
1) волновод неограничен вдоль продольной оси Z;
2) стенки волновода – идеально проводящие;
3) среда внутри волновода – идеальная диэлектрическая.
Определим компоненты волны, распространяющейся внутри волновода. Задачу решим для двух типов волн: Е-волны и Н-волны.
2.1. Волны в прямоугольном волноводе
2.1.1. Волна электрического типа (Е-волна)
Продольная компонента E& z определяется из решения волнового
уравнения
∇ 2s E& z + g 2 E& z = 0
(2.1)
для внутренней области структуры 0 ≤ x ≤ a, 0 ≤ y ≤ b (рис. 2.2) с
учетом четырех граничных условий
y
⎧ x = 0, x = a,
(2.1а)
E& z = 0 при ⎨
b
y
=
0
,
y
=
b
.
⎩
Поскольку E& z является тангенциальной
EZ
x
компонентой относительно ограничивающих
a
плоскостей, то на каждой из них она должна обРис. 2.2
ращаться в нуль.
Для решения уравнения используем метод разделения переменных (Фурье), то есть разобьем уравнение (2.1) на два: по координатам
X и Y . Так как общая форма записи направляемой волны
E& z = E0 ( x, y )e − jhz ,
то решение исходного уравнения можно записать в следующем виде:
E& z = X ( x) ⋅ Y ( y )e − jhz .
Подставим E& z в исходное уравнение (2.1):
14
(2.2)
∂ 2 E& Z
+
∂x 2
∂ 2 E& Z
∂y 2
+ g 2 E& Z = 0 .
Тогда с учетом координатного разбиения амплитуды по x и y
получим
X ′′ ⋅ Y + Y ′′ ⋅ X + g 2 XY = 0 .
Разделим обе части этого уравнения на XY:
X ′′ Y ′′
2
X
+
Y
+ g = 0.
Представим поперечное волновое число как g
волновое уравнение распадается на два:
⎧ X ′′
2
⎪⎪ X + g X = 0;
⎨
⎪Y ′′ + gY2 = 0,
⎪⎩ Y
2
= g 2X + gY2 , тогда
⎧⎪ X ′′ + g 2 ⋅ X = 0;
X
⎨
⎪⎩Y ′′ + gY2 ⋅ Y = 0.
или
Решения однородных уравнений запишем в следующем виде:
X = A cos( g X x) + B sin( g X x);
Y = C cos( gY y ) + D sin( gY y ) .
Тогда (2.2) примет такую форму записи:
E& Z = [ A cos g X x + B sin g X x ]⋅ [C cos gY y + D sin gY y ]⋅ e − jhz .
(2.3)
Определим постоянные величины A, B, C, D с использованием
четырех граничных условий (2.1а). Применим первую пару граничных
условий:
x = 0:
A ⋅ [C cos gY y + D sin gY y ] = 0;
x = a:
откуда
[A cos( g X a) + B sin( g X a)]⋅ [C cos( gY y) + D sin( gY y)] = 0 ,
⎧ A = 0;
⎨
⎩ B sin( g X a ) = 0.
Для следующей пары граничных условий имеем:
C ⋅ [ A cos g X x + B sin g X x ] = 0;
y = 0:
y = b:
откуда
Так как
[A cos( g X x) + B sin( g X x)]⋅ [C cos( gY b) + D sin( gY b)] = 0,
⎧C = 0;
⎨
⎩ D sin( gY b) = 0.
sin( g X a) = 0, sin( gY b) = 0,
⎧ g X a = mπ , m ≠ 0;
⎨
⎩ gY b = nπ , n ≠ 0.
15
то
Следовательно, коэффициенты C = 0 , A = 0 , а поперечные волновые числа
mπ
⎧
g
=
,
⎪⎪ X
a
⎨
⎪ g = nπ .
⎪⎩ Y
b
(2.4)
С учетом найденных коэффициентов ( C = 0 , A = 0 ), введя обозначение BD = E0 , из (2.3) получим окончательное решение волнового уравнения
E& Z = E0 sin ( g X x ) sin ( gY y )e − jhz .
(2.5)
Распределение E& z вдоль координат X и Y при m = 1, n = 1 показано на рис. 1.7, откуда убеждаемся в соблюдении граничных условий Eτ = 0 на стенках волновода.
Из уравнений взаимосвязи поперечных и продольных компонент (1.3) и (1.4) с Y
учетом H& Z = 0 определяются поперечные
компоненты:
EZ
h ∂ E& Z
⋅
E& X = − j
;
2 ∂x
g
h ∂ E& Z
⋅
E& Y = − j
;
2 ∂y
g
ωε ∂ E& Z
&
⋅
HY = − j
;
2 ∂x
g
ωε ∂ E& Z
&
⋅
.
HX = j
2 ∂y
g
X
Рис. 2.3
С учетом известной компоненты E& z (2.5) имеем
g
E& X = − jh X E& 0 cos( g X x) ⋅ sin( gY y ) ⋅ e − jhz ,
g2
g
E& Y = − jh Y E& 0 sin( g X x) ⋅ cos( gY y ) ⋅ e − jhz ,
g2
g
H& X = jωε Y E& 0 sin( g X x) ⋅ cos( gY y ) ⋅ e − jhz ,
g2
g
H& Y = − jωε X E& 0 cos( g X x) ⋅ sin( gY y ) ⋅ e − jhz .
g2
16
(2.6)
2.1.2. Волна магнитного типа (Н-волна)
Волновое уравнение для продольной компоненты волны имеет
вид
Внутренняя
область
∇ 2S H& Z + g 2 H& Z = 0.
структуры
(2.7)
ограничена
0 ≤ x ≤ a, 0 ≤ y ≤ b.
y
координатами
Граничные условия запишем предварительно через электрические компоненты E& X и
b
E& Y E& . В соответствии с рис. 2.4 тангенциальные
Y
компоненты E& X и E&Y должны обратиться в
a x
E& X
нуль в следующих координатах:
Рис. 2.4
⎧ x = 0,
E& Y = 0 при ⎨
⎩ x = a;
⎧ y = 0,
E& X = 0 при ⎨
⎩ y = b.
Решение уравнения имеет вид
H& Z = H Z ( x, y )e − jhz = H Z ( x) ⋅ H Z ( y ) ⋅ e − jhz = X ⋅ Y ⋅ e − jhz .
Аналогично случаю Е-волны, запишем общее решение для
Н-волны:
(2.8)
H& Z = ( A cos g X x + B sin g X x) ⋅ (C cos gY y + D sin gY y ) ⋅ e − jhz .
Перейдем от граничных условий для компонент E& X , E&Y к граничным условиям относительно H& Z . Использовав взаимосвязь продольных и поперечных компонент (1.3), (1.4) с учетом E& Z = 0 , получим
h ∂ H& Z
&
⋅
;
HX =−j
2
∂
x
g
h ∂ H& Z
&
⋅
;
HY = − j
2
∂
y
g
ωμ ∂ H& Z
&
⋅
;
EY = j
2
∂
x
g
ωμ ∂ H& Z
&
EX = − j
⋅
.
2
∂
y
g
Перейдем к граничным условиям, выраженным через H& Z :
E& X = 0,
то есть
∂ H& Z
= 0,
∂y
17
y = 0,
y = b;
E& Y = 0,
то есть
∂ H& Z
= 0,
∂x
x = 0,
x = a.
Вычислим производные от общего решения (2.8) по координатам
X и Y и подставим граничные условия:
x = 0:
∂ H& Z
= (− Ag X sin g X x + Bg X cos g X x) ⋅ (C cos gY y + D sin gY y ) = 0,
∂x
B = 0;
откуда
y = 0:
∂ H& Z
= ( A cos g X x + B sin g X x) ⋅ (−CgY sin gY y + DgY cos gY y ) = 0,
∂y
откуда
D = 0.
Подставляя в (2.8) B = 0 , D = 0 и используя вторую пару граничных условий, имеем:
− Ag X sin g X x = 0, то есть
mπ
gX =
;
a
g X a = mπ ,
или
− CgY sin gY y = 0,
nπ
gY =
.
b
gY b = nπ
или
x = a:
y = b:
то есть
С учетом найденных коэффициентов B = 0 , D = 0 , а также введя
обозначение AC = H 0 , получим окончательное решение волнового уравнения
H& Z = H 0 cos( g X x ) cos( gY y )e − jhz . (2.9)
Y
HZ
X
Очевидно, что на всех стенках структуры тангенциальная компонента
H& Z
Рис. 2.5
принимает максимальное значение, что иллюстрирует рис. 2.5.
Поперечные компоненты волны определим из уравнений взаимосвязи (1.3) и (1.4) при E& Z = 0 :
18
h ∂ H& Z
&
;
HX =−j
⋅
2
∂
x
g
h ∂ H& Z
&
;
HY = − j
⋅
2
∂
y
g
ωμ ∂ H& Z
&
;
EY = j
⋅
2
∂
x
g
ωμ ∂ H& Z
&
.
EX = −j
⋅
2
∂
y
g
Подставляя продольную компоненту (2.9) в эти уравнения, имеем
такие составляющие волн:
g
H& X = jh X H 0 sin( g X x) ⋅ cos( gY y ) ⋅ e − jhz ;
g2
g
H& Y = jh Y H 0 cos( g X x) ⋅ sin( gY y ) ⋅ e − jhz ;
g2
(2.10)
g
E& X = jωμ Y H 0 cos( g X x) ⋅ sin( gY y ) ⋅ e − jhz ;
g2
g
E& Y = − jωμ X H 0 sin( g X x) ⋅ cos( gY y ) ⋅ e − jhz .
g2
В данных выражениях введены поперечные волновые числа
mπ
nπ
, gY =
.
g 2 = g X 2 + gY 2 , g X =
a
Взаимосвязь всех волновых чисел такова:
2
2
2
b
k =h +g ,
где h – продольное волновое число; k – волновое число неограни-
ченного пространства.
2.1.3. Особенности волн
1. Волновод – это многоволновая направляющая система с типами волн Emn, Hmn, например: H01, H10, H11, H12, H21, …;
E11,
E12,
E21, … .
Типы волн E00, H00, E0n, Em0 не существуют, так как при такой комбинации индексов все проекции векторов обращаются в нуль.
Существующие типы волн разделяют на два класса:
• Emn – электрические типы волн (Е-волны), или поперечномагнитные волны (ТМ), у которых отсутствует продольная
магнитная компонента H Z = 0;
19
• Hmn – магнитные типы волн (Н-волны), или поперечноэлектрические волны (ТЕ), у которых E Z = 0.
Следовательно, в электрических волнах Emn обязательно присутствует продольная электрическая компонента E Z ≠ 0 , а в магнитных
волнах Hmn обязательно присутствует продольная магнитная компонента H Z ≠ 0 :
H& Z = 0; H mn → H& Z ≠ 0, E& Z = 0.
Emn → E& Z ≠ 0,
2. В поперечном сечении волновода поле имеет характер стоя-
чей волны. Числовые индексы в обозначении типа волны (m, n ) определяют количество полуволн стоячей волны вдоль соответствующей
стенки волновода: m – вдоль широкой стенки (координата X ); n –
вдоль узкой стенки (координата Y ).
На рис. 2.6 показано распределение продольной составляющей
&
E z для волны типа Е21, т.е. при m = 2 , n = 1 .
Для волн типов E0n, Em0 нулевой индекс
соответствует равномерному распределе- y
нию поля вдоль соответствующей координаты.
E21
Так как для Е-волны E Z ≠ 0, то при
равномерном распределении поля нарушается граничное условие Eτ = 0. Например,
x
Рис. 2.6
для волн типов Е10 и Е01 продольная составляющая поля вдоль соответствующей координаты неизменная (рис. 2.7).
y
y
E01
E10
x
x
Рис. 2.7
Тогда в координатах x = 0, y = b 2 и x = a, y = b 2 для волны
типа Е01 нарушается граничное условие EZ = Eτ ≠ 0, а для волны
типа Е10 – на широких стенках волновода. Это подтверждает невозможность существования волн типов E0n, Em0.
20
3. В продольном сечении волновода волновое число свободного пространства (ТЕМ-волны) раскладывается на продольную и поперечную составляющие (рис. 2.8):
k 2 = h2 + g 2 ,
(2.11)
где h – волновое число бегущей волны;
g – волновое число стоячей волны.
y
g
П тем
g
k
gx
h
Рис. 2.9
Рис. 2.8
gy
x
В поперечном сечении волновода поперечное волновое число
раскладывается на составляющие вдоль соответствующей поперечной координаты (рис. 2.9):
g 2 = g x2 + g 2y .
(2.12)
Составляющие поперечного волнового числа определяют изменение фазы стоячей волны вдоль соответствующей координаты:
ψx = gxx,
ψy = g y y.
Для симметричного выполнения граничных условий полное изменение фазы вдоль соответствующей стенки должно быть кратно π:
mπ = g x a,
nπ = g y b ,
что соответствует значениям волновых чисел
gx =
mπ
,
a
gy =
nπ
.
b
Волновые числа связаны с соответствующими длинами стоячих
волн:
gx =
откуда
т.е.
2π
,
λcx
gy =
2a
λcx = ,
m
a
m=
,
λcx /2
2π
,
λcy
2b
λcy = ,
n
b
n=
.
λcy /2
21
Поэтому второй формальный способ определения поперечных
волновых чисел следует из их определения:
gx =
2π
,
λcx
m=
a
λcx /2
,
4. По формальному правилу
1
m
=
.
λcx 2a
т.е.
λk определяется через критическое
значение волнового числа k k . Критический случай соответствует нулевому значению продольного волнового числа h = 0 , т.е. из выраже2
2
2
ния k = h + g имеем
k k = g = g x2 + g 2y .
Подставляя k k в выражение для
λk =
2π
=
kk
2π
g x2 + g 2y
λk , получаем
2π
=
2
⎛ mπ ⎞ ⎛ nπ ⎞
⎜
⎟ +⎜
⎟
⎝ a ⎠ ⎝ b ⎠
2
,
откуда
λk =
2
2
⎛m⎞
⎛n⎞
⎜ ⎟ +⎜ ⎟
⎝a⎠
⎝b⎠
2
.
(2.13)
Каждый тип волны, определяемый индексами m, n , имеет свое
значение λk .
5. Поскольку волновой режим соответствует действительному
значению h , т.е.
2
⎛ λ ⎞
h = k 1 − ⎜⎜ ⎟⎟ ,
⎝ λk ⎠
то при заданной длине волны λ (измеряемой в неограниченном пространстве) и заданных размерах волновода распространяются только
волны, удовлетворяющие условию λ < λk , причем
mn
λ
λв =
⎛ λ
1 − ⎜⎜
⎝ λk
.
⎞
⎟⎟
⎠
2
По мере уменьшения λ, начиная с низшей волны, могут распространяться волны высших типов.
Основная волна – это низший тип волны с наибольшим значением λk .
22
Определим основную волну для волновода R100 с размером поперечного сечения 10×23 мм:
Н10:
λk10 =
2
2
⎛m⎞ ⎛n⎞
⎜ ⎟ +⎜ ⎟
⎝ a ⎠ ⎝b⎠
Н01: λk 01 = 2b = 20 мм ;
2
= 2a = 46 мм ;
Н20: λk 20 = a = 23 мм .
Диаграмма типов волн приведена на рис. 2.10, откуда следует,
что Н10 – это низший тип. Основная волна Н10 имеет простейшую
структуру поля.
Вырожденные волны E11 и Н11, Е21 и Н21 имеют одинаковые λk .
Размеры волновода определяН01
ются из условий существования
Н20
только основной волны
Н10
λ
a < λ < 2a ,
λ > 2b ,
тогда его размеры будут
2b a 2a
Рис. 2.10
λ
2
< a < λ,
0<b<
λ
2
.
Размер широкой стенки определяет диапазон рабочих длин волн.
Теоретически отношение длин волн на краях диапазона составляет 2
единицы, а реальное отношение – 1,6….1,9 единиц.
6. Параметры волны h , λ В , ϑгр , ϑФ , Z C соответствуют параметрам дисперсных направляемых волн, рассмотренных выше.
7. Компоненты векторов поля основной волны Н10 находятся из
уравнений (2.10) при m = 1, n = 0 :
H& z = H 0 cos( g x x)e − jhz ;
h
H& x = j
H 0 sin( g x x)e − jhz ;
gx
(2.14)
ωμ
E& y = − j
H 0 sin( g x x)e − jhz .
gx
Выразим амплитуды компонент волны через амплитудное значение электрической компоненты. Обозначим амплитудный множитель
компоненты при E& y как
23
ωμ
gx
H 0 = E0 .
Тогда амплитудный множитель компоненты H& x будет
H x0 =
h
h gx
h
H0 =
⋅
E0 =
E0 .
gx
g x ωμ
ωμ
Отношение поперечных компонент определяет волновое сопротивление
E
ZC 0
ωμ
ωμ
Z CH = 0 =
,
=
=
2
2
HX0
h
⎛ λ ⎞
⎛ λ ⎞
⎟⎟
k 1 − ⎜⎜
λ
⎝ K⎠
⎟⎟
1 − ⎜⎜
λ
⎝ K⎠
поэтому
HX0 =
E0
.
Z CH
Амплитудный множитель компоненты H& z будет
g
E
g
E
π
H Z 0 = H0 = H X 0 x = 0 ⋅ x = 0 ⋅ .
h
Z CH h
Z CH ah
Компоненты основной волны, выраженные через E0 , имеют вид
E& y = − jE0 sin( g x x)e − jhz ,
E
H& x = j 0 sin( g x x)e − jhz ,
Z CH
(2.15)
E
π
H& z = 0 ⋅ cos( g x x)e − jhz ,
Z CH ah
а в стандартной форме записи –
⎛ πx ⎞
E& = − jY0 E0 sin ⎜ ⎟e − jhz ,
⎝a⎠
E ⎡
π
⎛ πx ⎞
⎛ πx ⎞⎤
H& = 0 ⎢ jX 0 sin ⎜ ⎟ + Z 0 cos⎜ ⎟⎥ e − jhz .
ah
⎝a⎠
⎝ a ⎠⎦
Z CH ⎣
(2.16)
2.1.4. Структура поля основной волны
Графическое представление изменений компонент волны изображено на рис. 2.11. Вдоль координаты Y компоненты поля не изменяются, так как n = 0 .
24
⎛ πx ⎞
H x`~ sin⎜ ⎟
⎝ a ⎠
⎛ πx ⎞
E y ~ sin⎜ ⎟
⎝ a ⎠
y
x
x
а
⎛ πx ⎞
H z ~ cos⎜ ⎟
⎝ a ⎠
y
y
x
а
а
Рис. 2.11
Z
Изменения компонент поля (2.15) вдоль продольной координаты
E y ~ sin(hz ) ,
H x ~ sin(hz ) ,
H z ~ cos(hz )
приведены на рис. 2.12.
Z
Z
HZ
х
х
y
х
х
х
λВ
2
EY, (-HX)
•
•
x
х
•
•
•
•
•
X
0
E
H
Рис. 2.12
Структуру поля (см. рис. 2.12) строят, используя графические
представления компонент волны вдоль соответствующих координат.
Сплошные линии соответствуют линиям вектора E , а пунктирные
– вектора H .
Структура поля основной волны имеет такие особенности:
1. Линии вектора E лежат в плоскости XOY поперечного сечения, так как эта волна является поперечно-электрической (ТЕ-волна),
а линии вектора H – в плоскости XOZ продольного сечения.
25
2. Электрическая компонента E максимальна в середине широкой стенки x = a 2 . Вдоль координаты Y поле постоянно.
3. Вдоль координаты Z максимумы E смещены на λВ/4 относительно центра семейства магнитных линий.
4. Структура поля и параметры волны не зависят от размера узкой стенки.
5. По внутренней поверхности волновода протекают поверхностные токи плотностью
jS = n0 × H ,
[
]
где n0 – внутренняя единичная нормаль (рис. 2.13).
Линии тока ортогональны линиям H и показаны на рис. 2.14.
y
H
n0 j s
js
H
Z
jCM
Рис. 2.13
jS
Рис. 2.15
Рис. 2.14
Линии тока, протекающего по внутренней поверхности волновода, имеют следующие особенности:
– линии тока пересекают ребра волновода;
– в середине широкой стенки x = a 2 линии тока замыкаются по
нормали к широкой стенке, переходя в ток смещения (рис. 2.15)
jCM =
∂D
;
∂t
– если щели в стенках волновода не
нарушают линий поверхностного тока, то
Z
такие щели – неизлучающие: продольная
щель в середине широкой стенки; поперечная щель в узкой стенке (рис. 2.16).
Поля высших типов волн строят, исходя из простейших типов. Базовыми
26
Рис. 2.16
структурами служат поля E11, H11 (рис. 2.17). Сложные структуры
Emn , H mn получают путем повторения структуры поля вдоль соответствующей координаты.
E11
y
H11
x
y
x
Рис. 2.17
2.2. Энергетические характеристики волновода
2.2.1. Вектор Пойнтинга в энергетических задачах
Задачи анализа, передачи и поглощения мощности решают с использованием соответствующих составляющих векE
тора Пойнтинга.
П
H
Рис. 2.18
Определим различные формы представления
вектора.
1. Мгновенный действительный вектор
П (t ) = E (t ) × H (t )
определяет направление передачи мощности (рис. 2.18). Его модуль
2
численно равен поверхностной плотности мощности, Вт м , или
энергии электромагнитной волны, переносимой через единичную по2
верхность за единицу времени, Дж с ⋅ м .
2. Комплексная форма мгновенного действительного вектора
Пойнтинга вытекает из представления действительной гармонической
функции в виде комбинации двух комплексно-сопряженных функций:
Ε0e j (ωt − kz ) + Ε0e − j (ωt − kz )
Ε (t ) =
= Ε0 cos(ωt − kz ).
(
)
2
Это соотношение можно записать в сокращенной форме, используя мгновенный комплексный вектор E& :
Ε& + Ε& ∗
,
Ε (t ) =
2
аналогично
27
& +Η
&∗
Η
,
Η (t ) =
2
& = Ε& e jωt – мгновенный комплексный вектор;
где Ε
m
Ε& m = Ε& 0 e − jkz – комплексная амплитуда.
Тогда действительный вектор Пойнтинга примет такой вид:
[(
)]
]
)(
1 & &*
E + E × H& + H& * =
4
1
= E& × H& + E& * × H& + E& × H& * + E& * × H& * ,
4
П (t ) = E (t ) × H (t ) =
[
где
[
]
(
1 &* & & & * 1 ⎡ &* &
E × H + E × H = ⎢ E × H + E& * × H&
4
4⎣
)* ⎤⎥⎦ =
= (a + jb ) + (a + jb )* = a + jb + a − jb = 2a =
[
]
[
1
1
*
= Re E& * × H& = Re E& m
× H& m
2
2
]
– слагаемые, не зависящие от времени;
1 & & &* & * 1 &
*
* − 2 jωt
E × H + E × H = Em × H& m e 2 jωt + E& m
× H& m
e
4
[
] 4 [(
(
)
)
]
– слагаемые, зависящие от времени.
С учетом представления этих слагаемых получим вектор П (t ) ,
выраженный через комплексные амплитуды поля:
1
1
*
*
*
* − 2 jωt
П (t ) = Re E& m
× H& m
+ E& m × H& m e 2 jωt + E& m
× H& m
e
. (2.17)
2
[
] 4 [(
(
)
)
]
3. Вектор, не зависящий от времени (среднее значение вектора),
можно найти в результате усреднения (2.17) по периоду его изменения:
1T
Пср = ∫ П (t )dt .
T0
T
2 jωt
Поскольку ∫ e
dt = 0, то в результате усреднения получим
0
1
*
Пср = Re E& m × H& m
,
(2.18)
2
где
[
]
Пср – усредненная за период поверхностная плотность мощ-
ности, которая не зависит от времени.
Для ортогональных векторов E& и H& запишем модуль вектора:
28
[
]
1
*
Пср = Re E& m ⋅ H& m
.
2
Так как E& m = Z c H& m , то
1
2
Пср = Re ⎡ Z c H& m ⎤ .
⎥⎦
2 ⎢⎣
(2.19)
2.2.2. Мощность, передаваемая по волноводу
Мощность, передаваемая по волноводу, определяется потоком
среднего значения вектора Π ср через поперечное сечение волновода (рис. 2.19):
[
]
1
& ∗ Ζ dS ,
Ρ = ∫ Π ср dS = Re ∫ Ε& × Η
0
2 S
S
x
S
y
Рис. 2.19
где S – площадь поперечного сечения волновода;
Ζ 0 – единичная нормаль к
z0
плоскости поперечного сеП СР чения.
Выразим общую форz
мулу мощности через поперечные
составляющие
волны:
Ε& = Ε& x X 0 + Ε& yY0 + Ε& z Z 0 = Ε& ⊥ +Ε& z Z 0 ,
& =Η
& +Η
& z Ζ0 ,
Η
⊥
& − составляющие векторов волны, лежащие в плоскости
где Ε& ⊥ , Η
⊥
поперечного сечения.
Используя выражение для комплексной амплитуды вектора
& = Ε& × Η
& ∗ , выделяя проекцию вектора Пойнтинга на проПойнтинга Π
дольную координату волновода (продольная составляющая вектора)
& z , получаем
&τ = Π
Π
0
& +z Η
&
&∗
& τ = Ε& ⊥ + z0Ε& z × Η
& ∗
Π
⊥
0 z ⋅ z0 = Ε ⊥ × Η ⊥ ⋅ z0 .
Здесь учтено, что смешанное произведение с двумя совпадаю& ⋅ z = 0.
щими векторами равно нулю, т.е. z0 × Η
⊥ 0
(
)(
)
[
]
& = Z [Η
&
Подставим Ε
⊥
c ⊥ × z 0 ], тогда
& × z ]× Η
&∗ ⋅z .
& τ = Zc [ Η
Π
⊥
0
⊥ 0
Применив
формулу
двойного
a × b × c = b (a c ) − c a b , получим
( )
29
векторного
произведения
& ×z ×Η
& * = z (Η
&
&*
&* &
& 2
Η
⊥
0
⊥
0 ⊥ ⋅ Η ⊥ ) − Η ⊥ ( Η ⊥ ⋅ z0 ) = z0 Η ⊥ .
Второе слагаемое равно нулю, как скалярное произведение орто& * (Η
&
гональных векторов, т.е. Η
⊥ ⊥ ⋅ z0 ) = 0 (в соответствии с рис. 2.19).
Тогда мощность, передаваемая по волноводу, определяется как
1
& ⊥ 2 dS .
(2.20)
Ρ = Re ∫ Z c Η
2 S
При действительном Z c (волновой процесс) знак Re опускается:
Ρ=
1 EH
1 1
2
2
ZC ∫ Η ⊥
dS =
Ε⊥
dS ,
∫
EH
2
2Z
S
S
(2.21)
C
где
2
Η⊥
= Η 2x + Η 2y − составляющая вектора в поперечном сечении
волновода;
Ζ cE (или Ζ cH ) – волновое сопротивление в зависимости от типа
волны;
S – площадь поперечного сечения волновода.
Аналогично можно получить выражение для передаваемой мощности P через продольные компоненты волны Ε z , Η z [9, с. 256].
2.2.3. Потери мощности в волноводе
Потери мощности в волноводе в основном определяются конечной проводимостью стенок. При протекании поверхностных токов
часть передаваемой мощности рассеивается в виде тепла (джоулевы
потери). Введем параметры, которыми описывается процесс энергетического обмена между электромагнитным полем и стенками волновода.
1. В идеальном металле протекает поверхностный ток js
(рис. 2.20), определяемый магнитной компонентой поля:
(2.22)
n0 × Η1 = js ,
где j s =
dI s
– линейная плотность тока, а I s – поверхностный ток.
dl
2. На поверхности реального проводника появляется небольшая
тангенциальная составляющая, которая определяется граничным условием Леонтовича
&τ ,
Ε& τ = Z c Η
или
30
[
]
& .
Ε& = Z c n0 × Η
(2.23)
3. Нормальная составляющая вектора Пойнтинга Π ⊥ определяет поток энергии, направленной внутрь стенок волновода
(рис. 2.21):
Π ⊥ = Ετ × Ητ .
4. Характеристики реального металла – диэлектрическая проницаемость ε& , волновое число k& , волновое сопротивление Z c – имеют
комплексный характер. Комплексная диэлектрическая проницаемость
в среде с потерями вводится в соответствии с уравнением Максвелла
& = jωε Ε& + σΕ& = Ε& ( jωε + σ ) = jωε&Ε& ,
rot Η
ε& = ε − j
где
σ
= ε (1 − jtgΔ ) ,
ω
(2.24)
а тангенс угла диэлектрических потерь
tgΔ =
Тогда волновое число
σ
.
ωε
(2.25)
k& = ω ε&μ = β − jα .
n0
IS
n0
js
dl
(2.26)
Hτ
Eτ
Hτ
П⊥
Рис. 2.20
Рис. 2.21
В среде с потерями коэффициенты затухания и фазы таковы:
k0
k0
2
2
α=
где
2
⋅
β=
1 + tg Δ − 1;
2
k 0 = ω εμ .
⋅
1 + tg Δ + 1,
В хорошо проводящей среде при tgΔ >> 1 эти коэффициенты
совпадают:
ω εμ
k
⋅
α = β = 0 ⋅ tgΔ =
2
2
σ
ωσμ
=
.
2
ωε
(2.27)
Волновое сопротивление такой среды с учетом комплексного характера диэлектрической проницаемости
Zc =
μ
μ
μ
=
=
=
− εjtgΔ
ε&
ε (1 − jtgΔ)
31
j
μωε 1 + j μω
=
⋅
,
εσ
σ
2
α
σ
Z c = (1 + j ) .
или
(2.28)
5. Внутри металла существует быстро затухающее поле
− jk&z
− j ( β − jα ) z
−αz − jβ z
&
Ε( z ) = Ετ e
= Ετ e
= Ετ e
e
,
где тангенциальная составляющая на поверхности Ετ = Ε(0) .
n0
Eτ
d
n0
E(z)
j(Z)
jm
z
z
Рис. 2.23
Рис. 2.22
Введем параметр глубины проникновения поля d . Примем, что
на глубине от поверхности z = d поле уменьшается в e раз
(рис. 2.22):
Ετ
−αd
−1
Ε(d ) = Ετ e
т.е.
= Ετ e ;
Ε( d )
αd = 1.
= e,
Глубина d определяется следующей формулой:
d=
1
α
=
2
ωσμ
.
(2.29)
Тогда волновое сопротивление (2.28) примет такой вид:
Z c = (1 + j )
1
.
σd
(2.30)
6. Внутри металла протекает объемный ток j ( z ) (рис. 2.23), который распределен по нормали к поверхности в соответствии с законом
изменения электрической компоненты поля
j ( z ) = σΕ& ( z ) = σΕτ e −α (1+ j ) z ,
(2.31)
причем на поверхности jm = σEτ .
7. Ток, который протекает в поверхностном слое толщиной d и
распределен в нем равномерно (рис. 2.24), называется эквивалентной
плотностью тока. Для получения этого тока проинтегрируем объемный ток (2.31) по всей глубине:
32
∞
&j = &j e −α (1+ i ) z dz = &j
э ∫ m
m
0
n0
1
d
1
;
α (1 + i )
jЭ
Рис. 2.24
писать
(2.32)
Поскольку в приповерхностном слое
Ε& τ = &jэ Z э , то эквивалентное сопротивление
Zэ =
• • •
&
&j = σΕτ = Ε& σd .
э
τ
α (1 + i )
1+ i
1+ i
.
σd
(2.33)
8. Поверхностный ток – это ток, который
протекает в бесконечно тонком слое (по поверхности). Для поверхностного тока можно за-
Ε& τ = Z S &jS ,
(2.34)
где Z S – поверхностное сопротивление,
&j – плотность поверхностного тока.
S
9. Комплексная амплитуда нормальной составляющей вектора
Пойнтинга направлена внутрь стенок волновода (см. рис. 2.21):
& = 1 [ Ε& × Η
& * ] = − 1 Ε& ⋅ Η
& τ* ] ⋅ n0 .
Π
⊥
τ
τ
τ
[
Поскольку
2
2
&τ ,
Ε& τ = Z S Η
то
2
& = −1Z Η
&
Π
⊥
S τ ⋅ n0 .
2
(2.35)
Комплексная мощность потерь определяется как поток нормальной составляющей вектора Пойнтинга через поверхность волновода:
1
(2.36)
Ρ& = Z S ∫ Ητ2 dS .
2 S
Средняя мощность потерь, которая определяет реальные потери,
вычисляется в соответствии со следующим выражением:
1
(2.37)
Ρср = Re Ρ& = RS ∫ Ητ2 dS ,
2 S
где RS – действительная часть поверхностного сопротивления,
RS = 1 σd .
2.2.4. Затухание волн прямоугольного волновода
Мощность потерь определяется действительной частью потока
нормальной составляющей комплексного вектора Пойнтинга через поверхность волновода:
33
[
& dS = Re 1 Ε& × Η
& * ]⋅ dS ,
Ρср = Re ∫ Π
∫ τ
⊥
τ
2
S
S
1
Ρср = RS ∫ Ητ2 dS .
2 S
Определим потери мощности на элементе длины волновода Δz ,
т.е. погонную мощность потерь (рис. 2.25):
Δz
1
1
2
ΔΡср = RS ∫ Ητ dS = RS . ∫ ∫ Ητ2 dldz;
2 S
2
0 L
ΔΡср =
где
∫
1
RS ⋅ Δz ∫ Ητ2 dl ,
2
L
Δz
− интеграл по периметру
L
поперечного сечения волновода.
Потери на единицу длины
волновода
Δl
z
Рис. 2.25
Ρср1 =
ΔΡср
R
= S ∫ Ητ2 dl .
2 L
Δz
Найдем взаимосвязь коэффициента затухания
α = f (Ρср1 ).
(2.38)
α и Ρср1 :
Направляемая затухающая волна имеет такой вид:
Ε& = Ε& 0 ( x, y )e −αz ⋅ e − jβ z ;
& =Η
& ( x, y )e −αz ⋅ e − jβz .
Η
0
Средняя мощность, передаваемая по волноводу,
1
1
Ρср = ∫ Π ср dS = Re ∫ E& × H& ∗ z0 dS = Re⋅ e − 2αz ∫ E& 0 × H& 0∗ z0 dS ,
2 S
2
S
S
где z0 – вектор единичной нормали, направленный внутрь стенок.
[
Обозначим Ρ0 =
]
[
[
1
& * ]⋅ z dS − мощность в сечении
Re ∫ Ε& 0 × Η
0
0
2 S
z = 0 . Вдоль оси волновода мощность уменьшается по закону
Ρср = Ρ0 e − 2αz ,
а изменение средней мощности
dΡср
dz
]
= (−2α ) ⋅ Ρ0 e − 2αz = (−2α ) ⋅ Ρср .
Отсюда коэффициент затухания, неп м ,
34
dΡср
1 dz
1 Ρср1
= ⋅
α =− ⋅
,
2 Ρср
2 Ρср
1 Ρср1
,
2 Ρср
α= ⋅
(2.39)
где Ρср − средняя мощность, передаваемая по волноводу;
Ρср1 = ΡΠ − средняя мощность потерь на единицу длины.
1
Это выражение представляет собой энергетическое определение
коэффициента затухания.
Так как
Ρср ( z ) = Ρ0e − 2αz ,
то после логарифмирования получим
(−2αz ) = ln
Ρср ( z )
Ρ0
,
а окончательное выражение для коэффициента затухания, неп м , в
логарифмическом масштабе имеет следующий вид:
1 1
z 2
α = ⋅ [ ln
Ρ0
].
Ρср ( z )
(2.40)
Единицы измерения коэффициента затухания в масштабе десятичных логарифмов:
Ρ
L ДБ = 10 lg 2 [дБ].
Ρ1
Ρ
LБ = lg 2 [Б] ,
Ρ1
Для перехода к неперам выполним следующее преобразование:
Ρ2
2αz
α = 10 lg
где
Ρ1
= 10 lg e
= 20 ⋅ (αz ) lg e = 20 ⋅ 0,434(αz ) = 8,68(αz ),
α имеет размерность неп м .
Из этого соотношения вытекает взаимосвязь непера и децибела:
1неп = 8,68дБ.
Взаимосвязь линейного и десятичного логарифмического масштаба представлена в табл. 2.1.
Таблица 2.1
L [дБ]
3
6
9
10
20
N [раз]
2
4
8
10
100
Получим выражение для коэффициента затухания, выраженное
через компоненты волны, исходя из энергетического определения затухания
35
1 Ρср1
.
2 Ρср
α= ⋅
Подставим в эту формулу выражения для мощностей
(
R
Ρср = ΡΠ1 = S ∫ Ητ2 dl ,
1
2 L
)
1
Ρср = Re ∫ E& 0 × H& 0* dS ,
2 S
тогда
Rs ∫ Ητ2 dl
1
L
,
& * z dS
2 Re ∫ Ε& 0 × Η
0 0
α= ⋅
[
]
(2.41)
S
что соответствует электродинамическому определению коэффициента затухания.
Для уменьшения затухания применяют материалы с большими
значениями удельной проводимости σ : медь, латунь (90%Cu). Слой
окисла увеличивает затухание, поэтому используют защитное покрытие: лак, Au, Ag.
Шероховатость поверхности увеличивает путь тока и, как следствие, затухание. В сантиметровом диапазоне волн глубина проникновения составляет d = 0,1 мкм , поэтому чистота поверхности соответствует 8 - 10 классам.
2.2.5. Коэффициент затухания основной волны
Исходным выражением для анализа является энергетическое определение коэффициента затухания
1 ΡΠ1
.
2 Ρср
α= ⋅
(2.42)
Сначала вычислим среднюю мощность, передаваемую по волноводу, через электрическую компоненту волны:
1 1
Ρср = ⋅
E⊥2 dS .
∫
2 ZH S
(2.42а)
C
Выпишем все компоненты основной волны Н10 прямоугольного
волновода и укажем их на рис. 2.26:
36
⎡
⎢
& z = Η 0 cos( g x x);
⎢Η
⎢
& x = jhΗ 0 g x sin( g x x);
⎢Η
⎢
g2
⎢
⎢Ε& = − jΗ (ωμ ) g x sin( g x).
0
x
⎢⎣ y
g2
y
x
z
b
E& y
H& z
H& x
a
Рис. 2.26
Электрическая компонента, которая лежит в плоскости поперечного сечения волновода,
Ε ⊥ = Ε y = Ε y 0 sin( g x x),
где введено обозначение для амплитудного значения
g
1
Ε y 0 = Η 0ωμ x = ωμ
⋅ Η0,
2
g
g
x
а поперечное волновое число g
2
n = 0.
Кроме того,
gx =
(2.43)
= g x2 + g 2y = g x2 , так как g y = 0 при
mπ π
= .
a
a
После выполнения соответствующих подстановок
(k =
ω
2π 1
, или ω = kυ =
)
⋅
λ
υ
εμ
запишем формулу для амплитуды
Ε y0 = H 0
ωμ
gx
= H0
2π
λ
⋅
a μ
a
μa
=2
H 0 = 2Z c0 H 0 ,
λ ε
λ
εμ π
1
⋅
(2.44)
где Z c 0 – волновое сопротивление среды в неограниченном пространстве.
Вычислим интеграл для расчета средней мощности (2.42а)
2
Ε
ba
1 − cos 2α
y
2
2
2 πx
dxdy = sin 2 α =
= 0 ab,
∫ E⊥ dS = Ε y 0 ∫ ∫ sin
a
2
2
00
S
тогда средняя мощность, передаваемая по волноводу,
2
2
Ε
Ε
1 1
1
y
y
(2.45)
Ρср = ⋅
⋅ 0 ab = ⋅ ab ⋅ 0 .
H
2 ZH
2
4
ZC
C
Подставляя волновое сопротивление для Н-волны
37
Z CH =
Z C0
⎛ λ
1 − ⎜⎜
⎝ λΚ
⎞
⎟⎟
⎠
2
и амплитудное значение E y 0 (2.44), получаем
4Z C2
3
⎛ λ
21
2a b
0 a 2
( ) = Η0
Ρср = Η 0 ⋅ ab ⋅
⋅ Z C0 ⋅ 1 − ⎜⎜
4
λ2
Z CH λ
⎝ λΚ
2
⎞
⎟⎟ . (2.46)
⎠
Далее рассчитаем потери мощности на единицу длины волновода
R
ΡΠ1 = S ∫ Ητ2 dl.
(2.47)
2 L
При вычислении интеграла по периметру поперечного сечения
волновода L на каждой плоскости нужно выбирать составляющие
вектора H , параллельные (тангенциальные) данной плоскости
(см. рис. 2.26). Для боковых поверхностей (плоскость x0y), которые
соответствуют координатам x = 0 , x = a , тангенциальная компонента
Hτ2 = H z2 , а на основаниях (плоскость x0z) с координатами y = 0 ,
y=b –
Hτ2 = H x2 + H z2 .
Вычислим интеграл по контуру поперечного сечения
a
⎤
⎡b 2
2
2
2
∫ Ητ dl = 2⎢ ∫ Η z dy x = 0 + ∫ (Η x + Η z )dx ⎥ =
⎦⎥
0
L
⎣⎢0
2
a⎡
⎡b
⎤ ⎤
2
2 πx
2 ⎛ πx ⎞ h
2 ⎛ πx ⎞
= 2Η 0 ⎢ ∫ cos
dy
+ ∫ ⎢cos ⎜ ⎟ +
sin ⎜ ⎟⎥dx ⎥ =
2
x
=
0
a
a
⎝
⎠
⎝ a ⎠⎥⎦ ⎥⎦
⎢⎣0
gx
0 ⎢⎣
⎡
a 1 h2 ⎤
2
= 2Η 0 ⎢b + + ⋅
⋅ a ⎥.
2
2
2
g x ⎥⎦
⎢⎣
2 2
Выразим отношение h g x через длину волны, используя выражение для продольного волнового числа
h=
2π
λb
1−
= 2π
38
λ
λ
λk
,
а также значение критической длины волны
h
2π a
=
⋅ =
g x λв π
⎛ λ
2a 1 − ⎜⎜
⎝ λΚ
⎞
⎟⎟
⎠
λk = 2a :
2
2
⎛ λ ⎞
⎛λ
⎛λ ⎞
⎟⎟ = ⎜ Κ
= ⎜ K ⎟ 1 − ⎜⎜
⎝ λ
⎝ λ ⎠
⎝ λΚ ⎠
λ
2
⎞
⎟ − 1.
⎠
Окончательно значение интеграла примет вид
⎡ a a ⎡⎛ 2a ⎞ 2 ⎤ ⎤
⎛
2a 3 ⎞⎟
2
2⎢
2
⎜
∫ Η τ dl = 2Η 0 b + 2 + 2 ⋅ ⎢⎜ λ ⎟ − 1⎥ ⎥ = 2Η 0 ⎜ b + 2 ⎟ ,
⎢
⎥⎥
⎢⎝ ⎠
⎣
L
⎦⎦
⎣
⎝
λ ⎠
а искомая мощность (2.47) будет
2a 3 ⎟⎞
2 ⎜⎛
ΡΠ1 = Η 0 b +
⋅R .
2 ⎟ S
⎜
⎝
λ ⎠
(2.48)
Теперь найдем коэффициент затухания (2.42) с учетом вычисленных интегралов (2.46) и (2.48):
⎞
2a 3 ⎛⎜ bλ2
2a 3
2
+ 1⎟
)
RS Η 0 (b +
2
3
⎜
⎟
2
λ ⎝ 2a
1
1
⎠ .
λ
α= ⋅
= ⋅ RS ⋅
2
2
2
2
3
⎛ λ ⎞
⎛ λ ⎞
a 3b
2a b
⎟⎟
⎟⎟
1 − ⎜⎜
Z
Η0
⋅ Z c0 1 − ⎜⎜
2
2 c0
λ
λ
λ
λ
⎝ Κ⎠
⎝ Κ⎠
После сокращений коэффициент затухания основной волны окончательно принимает вид
2
⎡
⎤
⎛ λ ⎞ ⎛ 2b ⎞
⎢1 + ⎜ ⎟ ⎜ ⎟⎥
RS ⎢⎣ ⎝ 2a ⎠ ⎝ a ⎠⎥⎦
⋅
,
α=
2
Z c0
⎛λ ⎞
b 1− ⎜ ⎟
⎝ 2a ⎠
где Z c 0 =
RS =
(2.49)
μ
– волновое сопротивление неограниченной среды;
ε
1
– поверхностное сопротивление стенок волновода.
σd
Частотная зависимость коэффициента затухания основной волны
для разных соотношений поперечного сечения приведена на
рис. 2.27. На рис. 2.28 показаны изменения коэффициентов затухания
основной волны и волны типа Е11.
39
b
= 0,1
a
α
α
b
= 0,5
a
b
=1
a
f
f кр
E11
H10
f
f кр
Рис. 2.27
Рис. 2.28
На основании приведенных зависимостей можно сделать следующие выводы:
1. При уменьшении частоты (увеличении длины волны) в окрестности критического значения коэффициент затухания резко возрастает и в пределе стремится к бесконечности.
2. При увеличении частоты наблюдается медленный рост коэффициента затухания, что связано с увеличением поверхностного сопротивления
RS =
1
=
σd
1
σ
2πfμ
ωμ
=
.
2σ
2σ
=
2
(2.50)
ωσμ
На высоких частотах предельное выражение коэффициента затухания (2.49) имеет вид
RS
.
bZ c0
α=
3. Для стандартного медного волновода асимптотическое значение
α≅
0,7
λ λ
⋅10 − 3
дБ
,
м
−2
следовательно, при λ = 10 м коэффициент затухания составляет
α = 0,7 дБ м , т.е. на расстоянии l = 10 м затухание α = 7 дБ (волна
затухает в 5 раз).
4. Коэффициент затухания уменьшается с увеличением периметра поперечного сечения:
α~
l
2( a + b )
⎡1 1⎤
=
= 2 ⎢ + ⎥.
ab
ab
⎣a b⎦
5. При постоянном периметре коэффициент затухания минимален при соотношении размеров b a = 1,18 , однако при b a = (0,5...2)
затухание меняется незначительно, поэтому выбирают b = 2a. При
40
b a > 0,5 ближайшим высшим типом волны является волна Н01, что
сокращает рабочий частотный диапазон волновода, поэтому ограничиваются соотношением b a ≤ 0,5 . Следует отметить, что волны низших типов имеют наименьшее затухание (см. рис. 2.28).
6. Волновод в сравнении с коаксиальным кабелем имеет меньшее затухание. Например, на частоте f = 3000 мГц , что соответст-
с 3 ⋅ 108
= 0,1 м , затухание в волноводе совует длине волны λ = =
f 3 ⋅ 109
ставляет α0 = 0,02 дБ/м, а в коаксиальной линии – α0 = 0,8 дБ/м. В
миллиметровом диапазоне волн прямоугольные волноводы практически не используются из-за большого затухания.
7. Полоса частот волновода ΔF = ±(20...25)% Fср , т.е. отношение
частот рабочего диапазона составляет Fmax/Fmin = 1,5.
8. Предельная мощность, передаваемая по волноводу, ограничена возможностью пробоя. Пробой происходит в области наибольшей
напряженности поля. Вид пробоя – искровой разряд при атмосферном давлении или тлеющий – при пониженном давлении. При пробое
происходит полное отражение волны (выход из строя генераторной
лампы).
Максимальная мощность пробоя определяется величиной передаваемой мощности (2.45) при предельном значении напряженности
поля:
⎛ λ
1 − ⎜⎜
ab 2
⎝ λk
Pпроб =
Eпр
4
Z с0
⎞
⎟⎟
⎠
2
.
(2.51)
Для воздуха Ε пр = 30 кВ см , что при λ = 3,2 см и поперечном сечении волновода 10x20 мм2 на основной волне дает значение предельной мощности Рпр = 990 кВт. Обычно допустимая мощность со-
ставляет Ρдоп = (0,2...0,3)Ρпр .
Неоднородности в волноводе уменьшают величину предельной
мощности.
41
3. КРУГЛЫЙ ВОЛНОВОД
Поперечное сечение такого волновода имеет вид окружности радиусом а. Такую структуру удобно представить в цилиндрической системе координат, изображенной на рис. 3.1, где
r – радиальная координата; θ – азимутальный
угол; Z – продольная координата.
0
θ
r
a
Начало координат соответствует центру
окружности – точке 0.
Пусть волновод работает на волне Е-типа.
Рис. 3.1
Определим продольную составляющую электрического поля E& z путем решения соответствующего волнового уравнения
∇ r2θ E& z + g 2 E& z = 0,
(3.1)
где 0 < r ≤ a , 0 < θ ≤ 2π – пространственная область уравнения;
∇ 2rθ – двумерный лапласиан в плоскости поперечного сечения.
Граничное условие следующее:
EZ = 0 при r = a .
Сначала получим выражение для двумерного лапласиана в
цилиндрической системе координат, а затем определим лапласиан в
криволинейной системе координат x1, x2, x3. Так как
∇ 2ϕ = divgradϕ ,
то предварительно запишем вражения для градиента и дивергенции:
gradϕ = ∇ϕ =
1 ∂ϕ
1 ∂ϕ
1 ∂ϕ
e1 +
e2 +
e3 ;
h1 ∂x1
h2 ∂x2
h3 ∂x3
⎤
1 ⎡ ∂
∂
∂
(
)
(
)
(
)
h
h
A
h
h
A
h
h
A
+
+
2 3 1
1 3 2
1 2 3 ⎥.
h1h2 h3 ⎢⎣ ∂x1
∂x2
∂x3
⎦
Подставляя вместо A1 , A2 , A3 проекции градиента
1 ∂ϕ
1 ∂ϕ
1 ∂ϕ
A1 =
A2 =
A3 =
,
,
,
h1 ∂x1
h2 ∂x2
h3 ∂x3
divA = ∇A =
получаем
divgradϕ = ∇ 2ϕ =
=
1 ⎡ ∂ ⎛ h2 h3 ∂ϕ ⎞ ∂ ⎛ h1h3 ∂ϕ ⎞ ∂ ⎛ h1h2 ∂ϕ ⎞⎤
⎜
⎟⎥,
⎜
⎟+
⎜
⎟+
⎢
h1h2 h3 ⎣ ∂x1 ⎜⎝ h1 ∂x1 ⎟⎠ ∂x2 ⎜⎝ h2 ∂x2 ⎟⎠ ∂x3 ⎜⎝ h3 ∂x3 ⎟⎠⎦
42
где h1 , h2 , h3 – коэффициенты Ламэ, которые определяют скорость
перемещения точки вдоль соответствующей координатной линии. Эти
числа являются коэффициентами линейных искажений.
Например,
дифференциальное
изменение
радиальной
координаты dr (рис. 3.2) приводит к соответствующему изменению
координатной линии на величину dl1 = dr . Поскольку линейные
искажения отсутствуют, то h1 =
dl1
= 1. Изменение азимутальной
dr
координаты dθ (рис. 3.3) приводит к соответствующему изменению
координатной линии ( r = const ) на величину dl2 . Тогда dl2 = rdθ , т.е.
h2 =
dl2
= r.
dθ
dr
dθ
r
dl2
r
dl1
Рис. 3.2
Рис. 3.3
Следовательно, для цилиндрической системы координат имеем
X1 X2 X3
r
θ
z
hi 1
r
1
Подставляя найденные коэффициенты hi и конкретные
цилиндрические координаты в лапласиан, выраженный в обобщенной
криволинейной системе, получаем
1 ⎡ ∂ ⎛ ∂ϕ ⎞ ∂ ⎛ 1 ∂ϕ ⎞⎤
∇ 2rθ = ⎢ ⎜ r
⎟+
⎜
⎟.
r ⎣ ∂r ⎝ ∂r ⎠ ∂θ ⎝ r ∂θ ⎠⎥⎦
Тогда волновое
принимает вид
уравнение
в
цилиндрических
координатах
1 ∂ ⎛ ∂E& z ⎞ 1 ∂ 2 E& z
+ g 2 E& 2 = 0,
⎜r
⎟+
r ∂r ⎝ ∂r ⎠ r 2 ∂ϕ 2
или
∂ 2 E& z 1 ∂ 2 E& z
1 ∂E& z 1
+ ⋅r ⋅
+
+ g 2 E& z = 0,
r ∂r r
∂r 2
r 2 ∂ϕ 2
т.е.
43
∂ 2 E& z
1 ∂E& z
1 ∂ 2 E& z
+
+
+ g 2 E& z = 0.
r ∂r
∂r 2
r 2 ∂ϕ 2
(3.2)
В соответствии с методом разделения переменных запишем решение в виде
E& z = R (r ) ⋅ Φ (ϕ ) ⋅ e − jhz ,
которое подставим в волновое уравнение (3.2):
∂R 1
∂ 2Φ
∂2R 1
+ g 2 RΦ = 0.
+ ⋅Φ⋅
+ ⋅R⋅
Φ
2 r
2
2
∂r
∂ϕ
r
1
Отсюда после умножения на
получим
RΦ
1 ∂ 2 R 1 1 ∂R 1 1 ∂ 2 Φ
+ ⋅ ⋅
+
+ g 2 = 0.
R ∂r 2 r R ∂r r 2 Φ ∂ϕ 2
Отделим переменные по r и ϕ :
∂r
1 ∂ 2Φ
r 2 ∂ 2 R r ∂R
2 2
.
⋅
+ ⋅
+g r =− ⋅
R ∂r 2 R ∂r
Φ ∂ϕ 2
Два независимых уравнения равны, если каждое из них равно
2
константе m :
r 2 ∂ 2 R r ∂R
⋅
+ ⋅
+ g 2r 2 = m2 ;
2 R ∂r
R
∂r
1 ∂ 2Φ
−
= m2.
Φ ∂ϕ 2
Так как функции R (r ), Ф(ϕ ) зависят от одной переменной, то перейдем от частных производных к обычным.
В первом уравнении получим единичный коэффициент при стар2
шей производной. Для этого разделим уравнение на r R :
d 2 R 1 dR
R
+
+ g 2r 2 − m2
= 0.
2
2
r dr
dr
(
)r
Окончательно имеем систему уравнений в канонической форме:
d 2 R 1 dR ⎜⎛ 2 m 2 ⎞⎟
⋅ R = 0;
(3.3)
+
+ g −
2
2 ⎟
r dr ⎜
dr
⎝
d 2Φ
dϕ
2
44
r ⎠
+ m 2 Φ = 0.
(3.3а)
Решив эти уравнения, получим следующее:
1. Φ (ϕ ) = A cos(mϕ ) + B sin( mϕ ) = A1 cos m(ϕ − ϕ 0 ),
при этом начало отсчета можно выбрать произвольно, т.е.
ϕ0 = 0 , то-
гда Ф(ϕ ) = A1 cos(mϕ ) .
2. Решение уравнения Бесселя (3.3) является суперпозицией цилиндрических функций 1-го и 2-го рода:
R(r ) = CI m ( gr ) + DN m ( gr ) ,
I m ( x ) – функция Бесселя m -порядка;
N m ( x ) – функция Неймана m -порядка.
где
Общее решение уравнения имеет вид
E& z = [CI m ( gr ) + DN m ( gr )]A1 cos(mϕ ) ⋅ e − jhz .
(3.4)
Графики первых функций Бесселя 1-го и 2-го рода изображены на
рис. 3.4.
Ii
y
1
I0(x)
N0(x)
Ni
I1(x)
N1(x)
0,5
I2(x)
0,6
0,2
x
1
2
3
4
5
6
7
8
x
2
9
4
6
8
-0,5
Рис. 3.4
Из физического смысла задачи в центре волновода r = 0 , x → 0 ,
N ( x ) → −∞ . Эта функция не удовлетворяет физическим условиям,
поэтому полагаем D = 0 . Тогда общее решение имеет вид
E& z = E0 I m ( gr ) ⋅ cos(mϕ ) ⋅ e − jhz .
(3.5)
Поперечное волновое число определяем из граничного условия
r = a,
E& z = 0,
тогда из формулы (3.5) имеем I m ( ga ) = 0, т.е. ga = ν mn , следовательно
ν
g = mn ,
a
45
(3.6)
где
m – порядок функции Бесселя, m = 0,1,2,…;
n – номер корня, n =1,2,3,…;
ν mn – корни функций Бесселя.
Значения корней функций Бесселя приведены в табл. 3.1.
Таблица 3.1
Номер корня n
Порядок функции m
1
3,832
7,016
10,173
0
2,405
5,520
8,654
1
2
3
2
5,135
8,417
11,620
0
Волны круглого волновода Emn нумеруются в соответствии с индексацией корней ν mn , т.е. порядком функции m и номером корня n ,
что определяет типы волн E01, E02, E11 E12 E21, ….
Особенности волн в круглом волноводе:
1. Компоненты волны зависят от волновых чисел, которые взаи2
2
2
мосвязаны соотношением k = g + h .
2. Критическая длина Е-волны определяется из указанного соотh = k 2 − g 2 . В граничном (критиченошения волновых чисел:
ском) случае продольное волновое число h = 0, тогда волновое число
неограниченного пространства совпадает с поперечным, т.е.
ν
2π
2π ν mn
слеk K = g = mn . Поскольку k K =
, то из уравнения
=
a
λk
λK
a
дует, что
λEK =
2πa
ν mn
Критическая длина Н-волны
λH
K =
.
(3.7)
2πa
,
qmn
(3.8)
где qmn – корни производных функций Бесселя.
Значения корней производных функций Бесселя представлены в
табл. 3.2.
Таблица 3.2
Порядок функции m
1
2
1,84
3,05
5,33
6,7
8,54
9,96
Номер корня n
1
2
3
0
3,83
7,02
10,17
46
3. Параметры волн соответствуют направляемым волнам при соответствующем значении λk . Например, продольное волновое число
2
⎛g⎞
h = k 1− ⎜ ⎟ ,
⎝k⎠
или
⎛ λ
h = k 1 − ⎜⎜
⎝ λK
2
⎞
⎟⎟ .
⎠
4. Основная волна определяется из сравнительной диаграммы
длин волн (рис. 3.5).
Основная волна соответствует наибольшему значению критической длины волны, т.е. минимальному корню ν mn ( qmn ).
Из табл. 3.1 и 3.2 следует,
Е11
что таким минимальным корнем
Н01
является q11, что соответствует
Н21
0
волне Н11, т.е. H11 – основная
волна. Критическая длина волны для этого типа максимальна
Е01
Н11
1,64а
λ
2,06а 2,61а 3,41а
Н11:
λK =
Рис. 3.5
Е01:
λK =
2πa
= 3,41a ,
1,84
2πa
= 2,61a и т.д.
2,4
Из диаграммы типов волн следует, что при работе на основной
волне рабочий диапазон длин волн заключен в диапазоне
3,41а…2,61а, что меньше, чем в прямоугольном волноводе.
5. Поскольку E ~ cos(mϕ ) ,
ϕ = 2π , то m соответствует числу це-
лых стоячих волн, которые укладываются вдоль окружности волновода.
Волны при m = 0 не изменяются в зависимости от ϕ , т.е. азимутально однородны. Волны при m = 1 имеют максимум продольной составляющей EZ при ϕ = 0 , ϕ = π .
Так как EZ ~ I m ( gr ) , то в критическом случае ga = ν mn , в свою
очередь,
gk =
2π
λcm
связано с λcт – длиной стоячей волны вдоль ко-
ординаты r , следовательно,
2π
λcm
a = ν mn ,
или
47
a
λcm / 2
ν
= mn .
π
Как следует из табл. 3.1 и 3.2, n определяет число максимумов
продольной составляющей EZ на радиусе а.
6. Структура волны типа Е01 азимутально-однородная ( m = 0 , по
E
× × ×
×
H
E01
. . .
. .
. . .
× × ×
× × ×
. .
Рис. 3.6
координате ϕ – равномерная) (рис. 3.6). Изменение компоненты E z
вдоль радиальной координаты r (3.5) определяет соответствующая
функция Бесселя
r⎞
⎛
E z ~ I m ( gr ) = I m ⎜ν mn ⎟
a⎠
⎝
с учетом представления поперечного волнового числа (3.6). Для волны типа Е01 это выражение конкретизируется ( m = 0 , n = 1 ):
r⎞
⎛
E z ~ I 0 ⎜ν 01 ⎟.
a⎠
⎝
Указанная зависимость E z (r ), изображенная на рис. 3.7, воспроизводит первую полуволну функции I 0 ( x ) (см. рис. 3.4) в интервале
0 ≤ x ≤ ν 01 .
Волна Е11 имеет два максимума E z по координате ϕ ( m = 1, одна стоячая волна) и один максимум вдоль координаты r . Изменение
компоненты E z (r ) показано на рис. 3.8, что соответствует функции
Бесселя первого порядка I1 ( x ) (см. рис. 3.4) в интервале 0 ≤ x ≤ ν 11 :
r⎞
⎛
E z ~ I1⎜ν11 ⎟.
a⎠
⎝
I0
r
r
rm
Рис. 3.7
Рис. 3.8
48
Координату максимума rm находим из производной функции Бесселя I1′ ( xm ) = I1′ (q11 ) = 0 , т.е. q11 = 1,84 , тогда
r
r
xm = ν11 m = 3,83 m = q11 = 1,84 ,
a
a
откуда
rm 1,84
=
= 0,48 .
a 3,83
Волна Н01 является азимутально-равномерной, поперечноэлектрической, т.е. линии Е – концентрические окружности (рис. 3.9).
H01
• •
• •
× ×
× ×
• •
• •
× ×
× ×
Рис. 3.9
H 110
• • •
• •
• • •
× × ×
× ×
×
Рис. 3.10
Рис. 3.11
По структуре простейшие волны – это волны типов Е01, Н01.
Структура основной волны представлена на рис. 3.10. На
рис. 3.11 показано, как основная волна прямоугольного волновода
0
H10 преобразуется в волну типа H11
при его деформации.
49
0
Основная волна H11 применяется редко из-за неустойчивости
плоскости поляризации. Неоднородности волновода приводят к повороту плоскости поляризации.
0
Волна типа E01 в силу круговой симметрии используется во вращающихся сочленениях.
0
Волна типа H 01 имеет наименьшие потери, причем при повышении частоты затухание уменьшается. Это связано с тем, что поверхностный ток имеет только поперечные составляющие, так как возбуждается продольной составляющей вектора H Z , поэтому линии поверхностного тока имеют кольцевой характер. При повышении частоты продольная составляющая вектора уменьшается, т.е. волна приближается к типу ТЕМ, что приводит к уменьшению плотности поверхностного тока. При соотношении a λ > 3...4 затухание этой вол0
ны намного меньше затухания волны типа H11 .
0
Круглые волноводы на волне типа H 01 имеют затухание
(1…2) дБ/км в коротковолновой части сантиметрового и миллиметрового диапазонов. Распространение волн других типов предотвращают, применяя дисковую или спиральную структуру волновода. Подавляемые волны имеют продольную составляющую тока, поэтому при
данной конструкции волновода такая структура тока эффективно подавляется.
4. КОАКСИАЛЬНАЯ ЛИНИЯ ПЕРЕДАЧИ
Эта линия представляет собой двухсвязную систему, состоящую
из двух проводников: центрального и экранирующего. Она изображена на рис. 4.1, где b = d 2 – радиус внутреннего проводника, а
a = D 2 – радиус внешнего проводника. В этой структуре существуют
волны трех типов: E, H и TEM. Основной волной является ТЕМ-волна.
Для такой волны линии векторов Ε и Η лежат в
D
плоскости поперечного сечения.
Компоненты волны в цилиндрической системе координат имеют следующий вид:
Ι&
⎧&
ϕ
, b ≤ r ≤ a;
Η
=
0
⎪⎪
2πr
⎨
&
⎪Ε& = r Z Ι ,
0 c
⎪⎩
2πr
d
(4.1)
Рис. 4.1
50
где Z c =
μ
– волновое сопротивление ТЕМ-волны.
ε
U&
Волновое сопротивление линии Z л = , при этом напряжение
Ι&
между проводниками
a
a
Ι& a dr
Ι&
Ι&
a
&
&
&
U = Zc
ln ,
U = ∫ Εdr = Z C
∫ = Z C ln r b ,
2π b
2π b r
2π
b
тогда
U& Z c a
Zл = =
ln ,
Ι& 2π b
Zл =
μ0
a
μ a 1
ln =
ln .
ε b 2π ε r ε 0 b
1
2π
Учитывая соотношение констант
μ0
= 120π ,
ε0
окончательно получаем
60
Zл =
εr
ln
a
[Ом],
b
(4.2)
где ε r – относительная диэлектрическая проницаемость внутренней
среды.
Коэффициент затухания
RS ∫ Ητ2 dl
1
L
α= ⋅
.
& &*
2
[
]
Re ∫ Ε × Η dS
S
Вычислим интегралы, входящие в это выражение:
Ι2
Ι2
2
2
2
⋅ 2πa +
⋅ 2πb =
а) ∫ Ητ dl = ∫ Η a dl + ∫ Η b dl =
2
2
(2πa)
(2πb)
L
L1
L2
2
=
[
I ⎛1 1⎞
⎜ + ⎟;
2π ⎝ a b ⎠
]
2π a
& ×Η
& * z dS = Z
б) ∫ Ε
0
c ∫ ∫
S
Ι2
0 b (2πr )
2
I2 a
= Zc
ln .
2π b
rdϕdr = Z c
Ι2
a
⋅ 2π ⋅ ln r =
b
(2π ) 2
Используя значения интегралов, получаем коэффициент затухания
51
1 R
α= ⋅ S
2 Zc
⎛1 1⎞
⎜ + ⎟
a b⎠
⋅⎝
.
a
ln
b
(4.3)
Особенности коэффициента затухания следующие:
1. Так как a > b, то большая энергия теряется в центральном
проводнике.
2. Если увеличивать радиус центрального проводника b , то коэффициент затухания уменьшится, так как уменьшится плотность поΙ2 1
верхностного тока и рассеиваемая мощность Ρв ~
⋅ .
2π b
3. Если увеличивать b , то уменьшится сопротивление линии
Z л , тогда при постоянной мощности
Ι2 a
Ρ = Ζc
ln = Ζ л Ι 2
2π
b
нужно увеличивать ток, что вызовет увеличение затухания. Следовательно, есть оптимальное значение радиуса:
dα
a
= 0 → ≅ 3,6.
db
b
Оптимальный выбор размеров a и b дает такие значения
волновых сопротивлений: Ζ л = 77 Ом для воздуха и Ζ л = 50 Ом для
заполнения полиэтиленом с ε r = 2,35 . Стандартные кабели имеют
сопротивления 50 и 75 Ом.
При длине волны λ<10 см затухание в линии превышает затухание в волноводе. Поэтому в сантиметровом диапазоне применяются
только короткие отрезки коаксиальной линии.
4. Коаксиальная линия работает на основной волне типа ТЕМ.
Качественное обоснование критической длины волны приведено на
рис. 4.2.
rСР
2а
2b
πrСР
Рис. 4.2
Если разрезать структуру по диаметру и “разогнуть”, то получим
сечение прямоугольного волновода с размером широкой стенки, равным половине длине окружности:
52
l=
где rc =
2πrc
= πrc ,
2
a+b
– средний радиус.
2
Значение критической длины волны для основного типа
λk ≅
π
2
(D + d ) = π (a + b ).
(4.4)
Условие существования ТЕМ-волны: λ > π (a + b ) . Поэтому при
уменьшении рабочей длины волны нужно уменьшать диаметры проводников.
Затухание электромагнитной волны определяется как джоулевыми потерями в проводниках, так и рассеиванием энергии в диэлектрике. Коэффициент затухания, неп м ,
α = αп + αд
содержит две компоненты затухания:
α п – в проводнике;
α д – в диэлектрике.
При измерении коэффициента затухания в логарифмическом
масштабе, дБ/м,
α = 8,686(α п + α д ).
Коэффициент затухания в диэлектрической среде для ТЕМволны определяется известным выражением
k
2
α д = tgΔ.
С учетом представления волнового числа в диэлектрике
k = 2π λд имеем
αд =
π
π
tgΔ =
ε tgΔ,
λ0 r
λд
(4.5)
ε r выражено через длину волны в вакууме λ0 и относительную диэлектрическую проницаемость диэлектрика ε r .
Параметр α д может быть выражен также через параметры ди-
где
λд = λ0
электрической среды:
1
α д = Z cσ ,
2
где
(4.5а)
Z c = μ ε – волновое сопротивление диэлектрической среды;
σ – удельная проводимость диэлектрика.
Действительно,
53
αд =
σ σ μ
1 σ
1
π
π σ
πσ
tgΔ =
⋅
=
= ⋅
=
.
εμ =
ε 2 ε
λд
λд ωε λ ε 2π 2 εϑф 2
д
T
(4.5б)
Коэффициенты затухания могут быть рассчитаны также на основе теории длинных линий. Эквивалентная Г-образная схема элемента длины линии определяет постоянную распространения
γ = ZY ,
Y = G + jωC – погонные значения последова-
где Z = R + jωL,
тельного сопротивления и параллельной проводимости, причем потери в проводниках и диэлектрике определяются соответственно параметрами R и G .
Преобразованное выражение постоянной распространения
R
G
1+
jωL
jωC
G (ωC ) << 1 дает приблина высоких частотах при R (ωL ) << 1,
γ = jω LC 1 +
женную формулу
⎛
1
R
1
G ⎞
1 R
1
+ ⋅
γ = jω LC ⎜⎜1 + ⋅
+ ⋅ GZ 0 = jβ + α ,
⎟⎟ = jω LC + ⋅
j
L
j
C
Z
2
ω
2
ω
2
⎝
⎠
0 2
где Z 0 = L C – волновое сопротивление линии.
Таким образом, коэффициенты затухания определяются следующими формулами:
- в диэлектрике
1
- в проводнике
α д = GZ 0 ;
2
(4.6)
1 R
,
2 Z0
(4.7)
αп = ⋅
где
Z
⎛b⎞
Z 0 = c ln⎜ ⎟,
2π ⎝ a ⎠
R = Rs
1 ⎛1 1⎞
⎜ + ⎟.
2π ⎝ b a ⎠
(4.8)
5. ПОЛОСКОВЫЕ ЛИНИИ ПЕРЕДАЧИ
Полосковые линии передачи относятся к двухсвязным системам.
Эти линии содержат два металлических проводника в виде тонких
лент, которые электрически изолированы друг от друга. Конструктивно они выполняются как симметричные или несимметричные линии.
54
Симметричная полосковая линия содержит центральный ленточный проводник, расположенный симметрично между плоскими металлическими экранами, которые электрически соединены между собой (рис. 5.1).
3
1
2
Рис. 5.2
Рис. 5.1
Несимметричная полосковая линия (рис. 5.2) представляет собой диэлектрическую пластину (подложку 1), расположенную между
металлическим заземленным экраном 2 и узким металлическим проводником 3 на противоположной плоскости. Несимметричная линия с
диэлектриком при большом значении относительной диэлектрической
проницаемости получила названия микрополосковой линии (МПЛ). В
такой линии толщина положки и ширина полоскового проводника много меньше длины волны в свободном пространстве. Малые размеры
линии позволяют использовать ее в интегральных СВЧ-схемах.
Основная волна в симметричной полосковой линии является
поперечной волной (ТЕМ-волна), которая, в частности, не имеет дисперсии.
В микрополосковой линии, строго говоря, распространяются
гибридные волны, которые содержат продольные электрические и
магнитные компоненты. Обоснование гибридных волн приведено в
[2]. Это связано с более сложной структурой поля, чем в симметричных линиях, так как только часть поля концентрируется в заполненном диэлектриком пространстве. Некоторая часть поля формируется
слоистой структурой диэлектрик–воздух. Однако малость продольных
компонент поля в сравнении с поперечными позволяет рассматривать
в качестве основной волны квази-ТЕМ-волну. Структура поля такой
волны показана на рис. 5.3. На низких частотах
E
эта модель дает приемлемую точность. С увелиH
чением частоты возрастает доля продольных составляющих поля, т.е. волна проявляет дисперсные свойства. При этом волновое сопротивление
и фазовая скорость начинают зависеть от частоты.
Рис. 5.3
Фазовая скорость ТЕМ-волны определяется
электродинамическими параметрами среды, в которой распространяется эта волна:
ϑф =
1
εμ
=
ϑ0
ϑ
≈ 0 ,
ε r μr
εr
55
(5.1)
где ϑ0 – скорость распространения волны в вакууме;
ка.
ε r – относительная диэлектрическая проницаемость диэлектри-
В МПЛ фазовая скорость квази-ТЕМ-волны рассчитывается через параметр эффективной диэлектрической проницаемости ε rэ :
ϑ0
.
ε rэ
ϑф =
(5.1а)
Этот параметр характеризует соотношение между энергиями
поля, концентрирующегося в диэлектрике и воздухе. Как следует из
теории длинных линий, скорость перемещения фронта волны вдоль
линии
1
LC
ϑф =
определяется ее погонными параметрами: емкостью и индуктивностью. Сравнивая ее с фазовой скоростью в линии при отсутствии диэлектрического заполнения
ϑ0 =
1
,
LC B
где C В – погонная емкость линии с воздушным заполнением, получаем
ϑ0
C
=
= ε rэ .
ϑф
CB
(5.2)
Этот параметр входит во все соотношения, используемые при
расчетах МПЛ.
Линия с широким проводником близка к плоскому конденсатору,
где практически вся энергия поля концентрируется в диэлектрике, поэтому ε rэ ≈ ε r . При очень узкой проводящей полоске энергия поля
распределяется почти поровну между воздухом и диэлектриком, поэтому ε rэ близка к полусумме проницаемостей диэлектрика и воздуха:
ε +1
ε rэ ≈ r
.
2
Диапазон изменения этого параметра составляет
ε r +1
2
< ε rэ < ε r .
В литературе приводятся расчетные соотношения при нулевой
толщине полоски [5]
56
ε r + 1 ε r − 1 ⎛ 10h ⎞ −1 2
ε rэ =
+
⎜1 +
⎟
w
2
2
⎝
и при конечной толщине полоски t
⎠
ε r + 1 ε r − 1 ⎛ 10h ⎞ −1 2
+
ε rэ =
− Q,
⎟
⎜1 +
2
2
w
⎝
где
ε −1 t h
Q= r
⋅
,
4 ,6
wh
h – толщина диэлектрика;
w – ширина центральной полоски;
⎠
(5.3)
(5.3а)
ε −1 t h
Q= r
⋅
.
4,6
wh
В работе [7] приведены иные варианты аппроксимирующих выражений для ε rэ . На частотах до 10 ГГц дисперсией (частотной зави-
ε rэ можно пренебречь. В [5] указаны источники, в которых
рассмотрены дисперсионные модели ε rэ . Там же приведено простое
симостью)
выражение для учета частотной зависимости.
Длина волны, распространяющейся в свободном пространстве,
ϑ
λ0 = 0
f
для квази-ТЕМ-волны при наличии диэлектрика, что соответствует
длине волны в МПЛ
λд =
ϑф
f
.
Сравнивая эти длины волн, имеем
λд ϑф
1
=
=
,
λ0 ϑ0
ε rэ
откуда длина волны в МПЛ
λд =
λ0
.
ε rэ
(5.4)
Рабочий диапазон частот МПЛ ограничивается сверху и снизу. Со
стороны верхних частот он определяется критическими частотами
ближайших высших типов волн. Наименьшие критические частоты
имеют волна типа НЕ1 и поверхностная волна, распространяющаяся
вдоль границы раздела диэлектрик–воздух. Для МПЛ, реализуемых в
интегральном исполнении, критическая частота поверхностной волны,
ГГц, [6]
57
fk =
где
h – толщина подложки, мм.
75
h εr −1
,
(5.5)
МПЛ применяют до частот 30 000 МГц. Физические ограничения
вызваны увеличением потерь в линии, а также требованиями высоких
допусков на их изготовление. Кроме того, на волнах короче 1 см преимуществ МПЛ в сравнении с волноводами с точки зрения габаритных
показателей нет. На частотах f > 14 ГГц преимущественно используются МПЛ в интегральном монолитном исполнении. В области низких
частот МПЛ могут применяться для передачи постоянного тока. Однако этот диапазон рекомендуется ограничивать частотой f >100 МГц.
Это связано с повышенными габаритами резонансных элементов на
более низких частотах.
Волновое сопротивление МПЛ определяется в результате решения уравнения Лапласа в частных производных численным методом.
Алгоритмы расчетов приведены в [7]. Приближенные расчетные формулы в замкнутой форме даны в [5]. Для симметричной полосковой
линии в частном случае бесконечно тонкого внутреннего проводника
получено выражение для волнового сопротивления [5]. Этот параметр
вычисляется через эллиптические интегралы. Выражение получено
на основе электростатического анализа методом конформных отображений. В этом же источнике приведено аппроксимированное выражение через элементарные функции.
Волновое сопротивление типовых МПЛ находится в пределах
20…125 Ом. Снизу значение сопротивления ограничено потерями на
излучение и преобразованием в высшие моды.
В несимметричной структуре полосковой линии трудно реализовать линию с большим волновым сопротивлением. Для этого нужно
сильно уменьшать ширину полоски, что резко увеличивает потери в
проводнике. Кроме того, требования на допуски при их изготовлении
становятся нереализуемыми. Уменьшение параметра ε также приводит к увеличению волнового сопротивления, однако при этом возрастают потери на излучение.
Затухание волны обусловлено потерями в проводниках и диэлектрике. Затухание в диэлектрике для полосковой линии определяется
классической формулой для ТЕМ-волны
k
2
α д = tgΔ =
где
π
π
tgΔ =
ε tgΔ,
λд
λ0 r
λ0 – длина волны в вакууме,
λд – длина волны в диэлектрике.
Эти длины связаны очевидным соотношением
58
(5.6)
λд =
λ0
.
εr
В логарифмическом масштабе затухание, дБ/м, определяется по
следующей формуле:
α д = 8,686π
εr
ε
tgΔ = 27,3 r tgΔ.
λ0
λ0
(5.6а)
В микрополосковой линии учитывается эффект перераспределения энергии электромагнитного поля между диэлектриком и окружающей средой, поэтому формула для затухания модифицируется
[5]:
α д = 27,3
εr
εr −1
ε − 1 tgΔ
⋅ rэ
⋅
.
ε rэ
λ0
(5.7)
Потери в диэлектрике обычно намного меньше потерь в проводниках, особенно при использовании высококачественной керамики
(кварц, поликор). При использовании полупроводниковой подложки
затухание существенно возрастает [6].
Потери в проводниках описываются достаточно громоздкими выражениями. Для МПЛ в замкнутой форме они приведены в [5], [6].
Простое приближенное выражение коэффициента затухания для несимметричной линии приведено в работе [6]:
α пр ≈
8,7 Rs
,
Zc w
(5.8)
где Z c – волновое сопротивление линии,
w – ширина центрального проводника.
Ток в проводниках концентрируется в слое, толщина которого определяется глубиной проникновения поля (толщиной скин-слоя).
Толщина проводника должна быть не менее 3 - 5 скин-слоев. Материал проводников должен иметь высокую электропроводность, малую
величину температурного коэффициента сопротивления, а также хорошие технологические свойства (адгезия к подложке).
Материал подложки определяет потери в диэлектрике и длину
волны в линии. Для МПЛ выбирают материал с высоким значением
ε r (~ 10 единиц), малыми потерями, постоянством ε в широком частотном диапазоне. В монолитных схемах используют Si (ε = 11,7),
GaAs (ε = 13,3). Кварцевая подложка характеризуется высокой стабильностью ε, бериллиевая керамика имеет высокую теплопроводность, что обеспечивает хороший теплоотвод.
59
Геометрические размеры МПЛ влияют на ее основные параметры. В то же время заданному значению волнового сопротивления может соответствовать разная геометрия линии.
Ширина основания а выбирается из условия a > 3w , где w – ширина токопроводящей полоски. В этой области пространства сосредоточена большая часть энергии поля.
Уменьшение ширины полоски приводит к уменьшению габаритов
линии, но при этом возрастают потери в проводнике. Волновое сопротивление линии монотонно уменьшается при увеличении параметра
w h , где h – толщина подложки. Уменьшение размера h обеспечивает меньшие потери на излучение, но для обеспечения постоянства
волнового сопротивления нужно уменьшать размер w , что приведет к
увеличению потерь в проводнике. Поэтому минимальное значение h
определяется допустимой величиной потерь в проводнике.
Библиографический список
1. Никольский В.В. Электродинамика и распространение радиоволн. – М.: Наука, 1973. – 607 с.
2. Григорьев А.Д. Электродинамика и техника СВЧ. – М.: Высш.
шк., 1990. – 335 с.
3. Гольдштейн Л.Д., Зернов Н.В. Электромагнитные поля и волны.
– М.: Сов. радио, 1971. – 664 с.
4. Петров Б.М. Электродинамика и распространение радиоволн. –
М.: Горячая линия-Телеком, 2003. – 558 с.
5. Гупта К., Гардж Р., Чадха Р. Машинное проектирование СВЧустройств: Пер. с англ. – М.: Радио и связь, 1987. – 432 с.
6. Малорацкий Л.Г., Явич Л.Р. Проектирование и расчет СВЧэлементов на полосковых линиях. – М.: Сов. радио, 1972. –
232 с.
7. Фуско В. СВЧ-цепи. Анализ и автоматизированное проектирование: Пер. с англ. – М.: Радио и связь, 1990. – 288 с.
8. Барсуков С.М. Теорія електромагнітного поля: Навч. посібник. –
Х.: ХАИ, 1994. – 62 с.
9. Вольман В.И., Пименов Ю.В. Техническая электродинамика. –
М.: Связь, 1971. – 487 с.
60
Оглавление
1.
Направляемые волны...................................................................
1.1. Особенности направляемой волны. Волновое уравнение
1.2. Взаимосвязь поперечных и продольных компонент волны
1.3. Типы направляемых волн………………………………………
1.4. Характеристики дисперсных волн…………………………….
2. Прямоугольный волновод.............................................................
2.1. Волны в прямоугольном волноводе…………………………..
2.1.1. Волна электрического типа (Е-волна)………………..
2.1.2. Волна магнитного типа (Н-волна)……………………...
2.1.3. Особенности волн………………………………………..
2.1.4. Структура поля основной волны………………………
2.2. Энергетические характеристики волновода………………...
2.2.1. Вектор Пойнтинга в энергетических задачах………..
2.2.2. Мощность, передаваемая по волноводу……………..
2.2.3. Потери мощности в волноводе…………………………
2.2.4. Затухание волн прямоугольного волновода…………
2.2.5. Коэффициент затухания основной волны……………
3. Круглый волновод………………………………………………….….
4. Коаксиальная линия передачи………………………………….…..
5. Полосковые линии передачи………………………………………..
Библиографический список………………………………………………..
61
3
3
4
6
8
11
11
11
14
16
21
24
24
26
27
30
33
39
47
51
57
Барсуков Сергей Николаевич
ЭЛЕКТРОМАГНИТНЫЕ ВОЛНЫ В НАПРАВЛЯЮЩИХ СТРУКТУРАХ
Редактор А.Н. Емленинова
Св. план, 2007
Подписано в печать
Формат 60х84 1/16. Бум. офс. № 2. Офс. печ.
Усл. печ. л. 3,2. Уч.-изд. л. 2,43. Т. 100 экз. Заказ
. Цена свободная
Национальный аэрокосмический университет им. Н.Е. Жуковского
«Харьковский авиационный институт»
61070, Харьков-70, ул. Чкалова, 17
http://www.khai.edu
Издательский центр «ХАИ»
61070, Харьков-70, ул. Чкалова, 17
izdat@khai.edu
62
Download