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

Кураев АА, Попкова ТЛ, Синицын АК 2004

.pdf
Скачиваний:
382
Добавлен:
15.06.2014
Размер:
5 Мб
Скачать

A&mne (0) = (e&mnee&mne+ ) Kzmn / Ktmn (0),

&e

 

ee+

(8.97)

Bmn (0) = (emn + emn ) W / Ktmn (0),

 

&

&

 

A&kмl(0) = (e&kмl+ + e&kмl) W / Ktkl (0),

B&kмl(0) = (e&kмl+ e&kмl) Kzkl / Ktkl (0),

Безразмерная мощность, переносимая волновым полем через поперечное сечение волновода в выбранных переменных имеет вид

 

 

2

&

e

&

e*

2

&

м

&

м*

 

(8.98)

P(z)= −α ∑ ∑χmnreal (Asmn Bsmn )+ χklreal (AsklBskl ) .

s

mn

 

 

 

 

 

kl

 

 

 

 

 

 

Исходя из представления (8.96) на регулярных участках, а также в точках волновода, где az=bz=qz=hz=0, мощности прямой и обратной волн в выбранных безразмерных переменных выражаются следующим образом:

 

 

 

 

 

 

 

 

 

 

j

 

&e

 

 

 

 

 

j

&e

*

P

 

 

 

 

 

 

 

 

 

dA

 

 

 

 

dB

 

+

±

= −α ∑ ∑χmn real Asmn ±

 

 

 

 

 

Bsmn ±

 

 

 

 

 

 

 

 

 

2

 

 

 

&e

 

 

 

 

smn

 

&e

 

 

 

 

smn

 

 

 

 

s

mn

 

 

 

 

Kzmn

 

 

dz

 

 

 

 

Kzmn

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

&м

 

 

 

 

j

&м

*

 

 

 

 

 

 

2

 

&м

dA

&м

 

dB

 

 

 

 

 

 

 

 

 

 

 

 

 

skl

 

 

 

 

 

skl

 

 

 

 

 

 

 

+χkl real Askl ±

 

 

 

 

Bskl ±

 

 

 

 

 

 

 

 

 

 

(8.99)

Kzkl

 

dz

 

Kzkl

dz

 

 

 

 

kl

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Эффективность взаимодействия электронного потока с возбуждаемым электромагнитным полем оценивается величиной электронного КПД

 

1

Ne

γi (z)γi (0)

 

 

ηe =

;

(8.100)

 

1 γi (0)

 

Ne i=1

 

 

или волнового КПД

ηv =

P(z)P(0)

;

ηv± =

P± (z)P± (0)

;

G0 (γ0 1)

G0 (γ0 1)

 

 

 

 

Расчет эффективности по разным формулам позволяет контролировать погрешность вычислений, т.к. при отсутствии потерь в стенках волновода из-

121

за сохранения баланса энергии должно выполняться условие ηe =ηv , а также на регулярных участках ηv =ηv+ +ηv.

Ниже приведены примеры моделирования приборов и устройств СВЧ на основе развитой теории.

ЛБВ-О на нерегулярном гофрированном волноводе прямоугольного се-

чения. Схема ЛБВ приведена на Рис.8.3. Формируемый электронной пушкой

Рис. 8.3.

ленточный электронный поток, проходя в зазоре волновода с гофрированны-

 

y

 

 

Dyp

 

 

 

x

Dxm

Dym

Dxp

Рис.8.4

ми стенками (Рис. 8.4), взаимодействует с усиливаемой Е11-волной и другими связанными с ней волнами, после чего осаждается на коллекторе. Ввод сигнала и вывод СВЧ энергии осуществляется через согласованныеr переходы. В

области взаимодействия создано однородное магнитное полеF = F0 zr0 , обес-

печивающее минимизацию поперечных составляющих скорости пучка. Профиль гофра задавался в виде:

D

= −D ;

D

ym

= −D

yo

h

ym

(z)sin2 (n π z / L) ;

(8.101)

xm

xo

 

 

 

 

 

 

 

ym

 

D

= D ;

D

yp

= D

yo

+ h

yp

(z)sin2 (n π z / L) .

 

xp

xo

 

 

 

 

 

 

yp

 

В математической модели ЛБВ-О приняты трехмерныеr уравнения движения крупных частиц (8.89) с начальными условиями βi0 = β0 zr0 , rr0i = 0 . Начальные условия для амплитуд (8.87) задавались исходя из того, что на вход

122

подается сигнал в виде падающей Е11 - волны с амплитудой e11e+re = e0 , ампли-

туды всех остальных распространяющихся волн задавались равными нулю. Тем самым обеспечивается полное согласование на входном конце лампы. Исходя из условия селекции рабочей E11-волны на базовой частоте, были выбраны безразмерные (отнесенные к λ0 2π ) значения Dxo=10, Dyo=2 . Опти-

мизировались: глубина гофра hym(z)=hyp(z)=hy(z), длина гребенки L, количество периодов nym=nyp=ny, , ток пучка I0 амплитуда сигнала e0. Анализ показывает, что в волноводе при таких размерах распространяется 12 типов волн, однако в результате расчетов выяснилось, что основной "паразитной" волной в рассматриваемой ЛБВ-О является волна НЕ11, которая связана в нерегулярном волноводе с волной Е11 и оказывает существенное влияние на процесс усиления. Все другие распространяющиеся типы волн не возбуждаются.

Был произведен поиск оптимизированных по величине волнового КПД вариантов ЛБВ-О в диапазоне напряжений электронного пучка (0.5≤β00.85). Анализ полученных результатов показал необходимость учета НE11 - волны. Оказалось, что предварительные расчеты по упрощенной модели с учетом только Е11- волны дают значительную погрешность. При β0 >0.7 были най-

дены варианты с КПД 36% на регулярной гребенке и с КПД 57% на нерегулярной. Ниже приведены два варианта.

Вариант1(регулярная гребенка):

β0 =0.75,

I0 =100A,

L =71, ny = 40, hy =1.34,

e0 = 0.0025,

η = 32%,

Ku =13 Дб.

Вариант2: (нерегулярная гребенка)

β0 =0.85,

I0 = 386 A,

L = 88, ny = 40, η = 57%, Кu = 23 Дб.

В этом варианте глубина гофра изменяется в пределах 1.005 < hy <1.345.

За счет оптимального профилирования глубины гофра эффективность лампы может быть увеличена примерно в 1.5-2 раза.

Таким образом, полученная математическая модель позволяет производить расчеты и оптимизацию ЛБВ-О с нерегулярным профилем гофра замедляющей системы. Проведенные исследования показали высокую эффективность ЛБВ-О с замедляющей системой в виде гребенки с оптимально профилированным гофром. Причем, с увеличением напряжения электронного пучка КПД такой лампы возрастает.

Заграждающий фильтр на основе периодического гофра. Рассмотрим расчет заграждающего фильтра для волны Н10 - типа. Геометрия отрезка волновода, выполняющего функции фильтра представлена на Рис. 8.4. Профиль гофра задавался формулами (8.101).

123

При расчете фильтра граничные условия выбираются следующим образом. На входе задаются падающая Н10 волна с амплитудой e&01м+ =1 + j0 и от-

раженная e&01м= e0 e jϕ0 . Амплитуды всех других волн задавались равными

нулю. На выходном конце при z=L задается условие отсутствия встречной волны P(L)=0 и минимизируется величина проходящей мощности P+(L). В

уравнениях (19)-(24) полагается I&C = I&Be = I&=0 .

Задавались значения Dx0 , Dy0 и методом пристрелки подбирались зна-

чения e0, ϕ0, L, hy, ny из условия P(L)=0 и min P+(L).

Расчеты, в частности, показали, что при Dx0 < π, Dy0 < π, (т.е. пока не распространяются волны Н12, Е12), рассматриваемый периодический фильтр (hy=const) при ny=10 обеспечивает 25ДБ затухания. При этом волны Н01, Н11,

Е11 не возбуждаются.

Один из вариантов такого фильтра приведен ниже:

Dx0 = 1.75, Dy0 = 2.75, L =70.5 , nv = 10, hy=1.38, P+(L)/P+(0) =

0.0034.

На Рис. 8.5 приведено распределение мощностей прямой и встречной волн (8.99) P+(z), P(z). Полоса заграждения такого фильтра уменьшается с увеличением ny.

P+, P-

2

 

 

 

 

 

1,5

 

 

 

 

 

 

 

P+

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

0,5

 

 

 

 

 

 

0

 

 

 

 

 

 

-0,5

 

P-

 

 

 

 

-1

 

 

 

 

 

 

 

 

 

 

 

-1,5

 

 

 

 

 

 

-2

 

 

 

 

 

 

0

0,2

0,4

0,6

0,8

1

 

 

 

 

 

 

z/L

 

 

 

Рис 8.5.

 

 

 

Расчет дисперсионных характеристик периодического волновода. При расчете дисперсионных характеристик волн в периодических волноводах искомое решение для поля в волноводе подчиняется условию Флоке:

r

z + kd ),

r

r

r

 

(E&(x, y,

B&(x, y,

z + kd ))= (E&(x, y, z),

B&(x, y, z)) e jkϕ0 .

(8.102)

124

здесь d - период, 0 ϕ0 π - набег фазы нулевой пространственной гармоники поля. Постоянные распространения пространственных гармоник определяются через ϕ0 и d по формуле:

 

 

βφn =

 

d

 

 

 

; n =0,

±1, ±2, ...,

βφn = vφn / c

(8.103)

 

 

ϕ0 + 2nπ

vφn - фазовая скорость гармоники.

 

 

 

 

Hkl)

При расчете зависимости ϕ0(W) для волны определенного типа (Emn или

граничные

 

условия

при

z=0

задаются

следующим

образом:

&+

=1,

&

= eor . Методом пристрелки подбираются eor, ϕ0, для которых вы-

eij

eij

 

полняются условия при z=d , вытекающие из (8.102):

 

&

 

(d )

&

 

 

 

 

jϕ

0

;

&

 

 

&

jϕ0

 

 

 

Aij

= Aij (0) e

 

 

Bij (d )=

Bij (0) e .

 

 

 

 

 

Анализ решений этой задачи, в частности, показал, что для отрезка пе-

риодического

волновода нельзя

одновременно

удовлетворить

условию

(8.102)

и

 

 

условиям

полного

согласования

на входе и

выходе

&

 

 

&

 

(0)

&

 

 

&

(L)=0 .

Следовательно,

решение, удовлетворяю-

A

(0)= B

 

 

= A

(L)= B

 

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

В качестве примера расчета на Рис. 8.6. приведена зависимость ϕ0(W) для периодического гофрированного волновода, используемого в приведенном выше варианте 1 ЛБВ-0. Верхняя кривая соответствует расчету для Е11+НЕ11 волны, нижняя - для Е11 без учета НЕ11. Анализ показывает, что при

W& =1 нулевая пространственная гармоника Е11+НЕ11 имеет замедление βφ0 = 0.71, эта же гармоника Е11-волны имеет замедление βφ0 =0.755 . Таким

образом, усиление в этом варианте реализуется при синхронизме именно с пространственной гармоникой комбинированной Е11+НЕ11 волны.

Развитая строгая теория произвольно-нерегулярного волновода с прямоугольным сечением позволяет корректно решать широкий круг задач анализа и оптимизации электронных приборов и устройств СВЧ. Приведенные в статье примеры демонстрируют недопустимость использования в этих задачах упрощенных теорий, в которых не учитывается связь различных типов волн в нерегулярных периодических и непериодических волноводах.

125

3

 

 

 

 

 

 

ϕ0

 

 

 

 

 

 

2,5

 

E11+HE11

 

 

 

2

 

 

 

E11

 

 

 

 

 

 

 

 

1,5

 

 

 

 

 

 

1

0,85

0,9

0,95

1

1,05

1,1

0,8

 

 

 

 

 

 

W

 

 

 

Рис.8.6.

 

 

 

126

8.9 Т-функции для решения двухточечных задач в теории нерегулярных волноводов

В п.8.4 развита общая теория произвольно-нерегулярных волноводов, основанная на отображении произвольно-нерегулярной внутренней поверхности волновода на регулярный цилиндр. В преобразованной (косоугольной) системе координат решение представляется в виде связанных нормальных волн с использованием проекционной процедуры. При этом амплитуды связанных волн определяются системой ОДУ с переменными коэффициентами, вид которых определяется профилем нерегулярного волновода. Граничные условия к этой системе ставятся в начальном и конечном сечении отрезка нерегулярного волновода (двухточечная задача). Решение этой задачи традиционными методами не встречает затруднений, если рассматриваются только распространяющиеся волны. Как показано ниже, для точного расчета волновода необходим учет и связанных с распространяющимися закритических волн, существенно меняющих характеристики волновода. Однако для закритических волн численное решение граничной (двухточечной) задачи с использованием пошаговых методов типа Рунге-Кутта или Хемминга невозможно из-за их быстрой расходимости (из-за малых ошибок появляются резко возрастающие решения). В этом случае необходимо строить аналитические решения на системе заданных узловых точек, удовлетворяющие граничным условиям краевой задачи и представляющие собой разложение искомых функций в базисе специальных функций, обеспечивающих разрешимость получающейся системы алгебраических уравнений относительно коэффициентов разложения. Такие специальные дискретно определенные Т-функции введены и описаны ниже.

8.9.1 Т-функции

Будем предполагать, что эти функции определены на системе N равноотстоящих точек tm: t0, t1tN-1 интервала LN. Интервал LN определим следующим образом:

LN =

2π

sin(π / N )

= 2N sin(π / N ).

(8.104)

 

 

 

 

 

(π / N )

 

 

В пределе при N → ∞ :

 

 

lim L

N

= 2π lim

sin(π / N )

 

= 2π.

 

 

 

n→∞

 

 

 

(π / N )

 

 

127

Кроме того, tm = mh, h =

LN

= 2sin(π / N ).

(8.105)

N

 

 

 

Определим системы базисных функций k-го порядка первого рода tcNk(ktm) и второго рода tsNk(ktm) следующим образом:

 

 

 

 

 

2πk

 

 

 

πk tm

 

 

 

tc

 

(kt

 

 

t

 

 

 

 

 

 

= cos

 

 

 

,

 

Nk

 

m

 

LN

 

m

 

N sin(π /

 

 

 

 

 

 

 

 

 

 

 

N)

 

k = -2, -1, 0, 1, 2… k – целое число.

 

 

 

 

 

2πk

 

 

 

πk t

m

 

 

ts

Nk

(kt

m

) = sin

 

t

 

= sin

 

 

.

 

 

 

 

 

 

 

LN

 

m

 

N sin(π

/

 

 

 

 

 

 

 

 

 

 

N )

Введем также комплексную Nk-функцию

teNk (ktm ) = tcNk (ktm ) + jtsNk (ktm ).

В пределе N → ∞, tN 0, tm t , имеем

lim tcNk (ktm ) = lim cos(

(π / N )

 

ktm ) = coskt,

sin(π / N )

N →∞

N →∞

 

 

lim tsNk (ktm ) = lim sin(

 

(π / N )

ktm ) = sin kt,

sin(π / N )

N →∞

N →∞

 

 

lim teNk (ktm ) = cos(kt) + j sin kt = e jkt .

N →∞

(8.106)

(8.107)

(8.108)

Таким образом, в пределе, для непрерывного аргумента t, TNk-функции переходят в обычные тригонометрические функции.

В дискретном же варианте, как следует из определений (8.106), (8.107), имеют место иные соотношения:

 

 

tcNk (ktm ) = cos(kxm ) ,

(8.109)

 

 

tsNk (ktm ) = sin(kxm ) ,

(8.110)

где xm =

2π m,

tm = 2sin(π / N ) m .

 

 

N

N

 

Таким образом, аргументы в Т–функциях и тригонометрических функциях k-го порядка различны. Это различие и реализует основное свойство Т

128

функций на дискретном множестве точек tm, приводящее к аналитическим решениям: производные от этих функций на множестве tm с точностью до коэффициента равны самим функциям. Действительно:

tcNk (ktm ) = −Rk tsNk (ktm ) ,

(8.111)

tsNk (ktm ) = Rk tcNk (ktm ) ,

(8.112)

teNk (ktm ) = jRk t&eNk (ktm ) ,

(8.113)

tc′′Nk (ktт) = −Rk2tcNk (ktm ) ,

(8.114)

ts′′Nk (ktт) = −Rk2tsNk (ktm ) ,

(8.115)

te′′Nk (ktт) = −Rk2teNk (ktm ) ,

(8.116)

Rk = sin(πk / N ) / sin(π / N ) .

 

8.9.2 Взаимодействие Hoi волн в гофрированном волноводе с круговым сечением

Воспользуемся общей теорией нерегулярных волноводов, развитой в [26]. В рассматриваемом случае система дифференциальных уравнений для

амплитуд связанных волн H01(A&1 ), …, H0i (A&i ) имеет вид (источники внутри отрезка волновода отсутствуют):

d

2 &

 

 

 

 

 

 

 

 

 

 

 

2

1

 

1

 

 

1 dg

 

2

 

 

 

 

 

 

 

 

 

 

 

 

Ai

 

 

 

 

 

 

µi

 

 

 

2

 

&

 

 

 

 

 

 

 

 

 

 

 

dT 2

 

+

 

1

 

 

 

g 2

3

µ

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

µ

 

 

 

 

 

i

g dT

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

1 dg

 

 

 

 

&

 

 

 

1 dg

2

 

(5)

 

 

1 d

2

g

 

(6)

 

&

 

 

 

 

 

(4) dAj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ

ij

 

 

 

 

 

+

 

 

 

 

 

γ

ij

 

 

 

 

2

γ

ij

A

j

= 0,

(8.117)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j=1 g dT

 

 

 

dT

 

 

 

 

 

 

 

 

 

 

 

g dT

 

 

 

 

 

 

 

 

 

 

 

 

 

g dT

 

 

 

 

 

 

 

 

 

 

 

 

 

ji

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i =1,2, ...

Здесь Т = 2λπ z, z - расстояние вдоль оси волновода, 0 T L0 , λ -

длина волны в свободном пространстве; g(T) = b(z)/b1кр, b(z) – радиус внут-

129

ренней поверхности нерегулярного волновода, b1кр – критический радиус для

волны Н01, b1кр = λ µ1 / 2π, µi (i =1, 2, 3...)

 

- корни производной функции Бес-

селя 1-го рода нулевого порядка (J0(µi )= 0).

 

 

 

 

 

 

 

γij(4)

=

 

 

 

µi µ j

 

2J

0

(µ

i

)

,

 

 

 

 

 

 

 

 

 

µi2 µ2j

 

J0 (µ j

)

 

 

 

 

 

 

 

 

ij

 

 

 

 

 

 

 

 

 

 

 

 

 

γij(5)

 

 

µi µ j (µi 2 + 3µ j 2 )

 

 

2J

0

(µ

)

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

,

 

(µi2

µ2j )2

 

 

 

 

 

J0 (µ j

)

ij

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γij(6)

= −

µi µ j

 

 

 

2J

0

(µ

i

)

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

 

 

 

ij

 

 

 

 

µi2 µ2j J0

(µ j

 

 

 

 

 

 

 

 

 

В случае, когда а) вход и выход отрезка нерегулярного волновода со-

гласованы; б) на входе и выходе выполняются условия: dTdg (0) = dTdg (L0 ) = 0;

в) сигнал подается только с левого конца и только на волне Н01, граничные условия в системе (8.117) имеют вид [26]:

 

 

 

 

&

 

 

 

 

 

µ

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dA

 

 

 

i

 

 

 

 

2

&

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

(0)

= j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dT

1

µ1

 

 

 

/ g0

Ai (0),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

(8.118)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&

 

 

 

 

 

 

µi

 

 

 

 

2

&

 

 

 

 

 

 

 

 

 

 

 

dAi

(L

 

 

 

/ g

(L

)

 

 

 

 

 

 

 

 

 

 

dT

) = − j 1

µ

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

1

 

 

 

 

 

L

i

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

µ

i

2

 

2

 

 

 

µ

i

2

 

 

2

 

 

 

µ

1

 

2

2

 

причем, при

 

 

 

 

/ g

>1

 

 

 

 

 

/ g

= − j

 

 

 

/ g

1 . Здесь g0 =

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

µi

 

 

 

 

µ1

 

 

 

 

 

 

µ1

 

 

 

 

 

 

 

 

 

 

 

 

g(0), gL = g(L0).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в виде, удобном для дальнейшего

Рассмотрим

исходную систему (8.117)

численного счета:

d 2 A&i

dT 2

где

+ Q

&

I

 

 

(T )

(T )A

G

 

i

i

j=1

 

i j

 

 

 

ji

 

 

 

 

 

 

 

 

Gij

 

 

&

 

 

 

 

&

 

 

 

dAj

 

 

 

 

 

 

 

 

 

+ Qi j (T )Aj

= 0, i =1, 2,...,

(8.119)

 

dT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 dg

(4)

, i j ,

 

=

 

 

 

 

γ ij

 

g dT

 

 

 

 

 

 

 

130