CC BY
66
zing complicated H-plane structures by means of generalized scattering matrices. For triangular boundaries, a specific field representation is proposed. It is based on the domain-product technique and ensures a good convergence of the used series.
Numerical examples demonstrating reliability and efficiency of the approach are presented.
Key words: waveguide, H-plane structure, scattering, field representation, scattering matrix.
УДК 621.372.81 Логачева Л. М.1, Куцак С. В.2, Бондарев В. П.3
1 Старший преподаватель Запорожского национального технического университета 2Аспирант Запорожского национального технического университета 3Канд. физ.-мат. наук, доцент Запорожского национального технического университета
ДИФРАКЦИЯ ВОЛНЫ Н10 В ПРЯМОУГОЛЬНОМ ВОЛНОВОДЕ НА СКАЧКЕ ИМПЕДАНСА УЗКИХ СТЕНОК
На основе проекционного метода дается строгое решение задачи дифракции основной волны Н10 на стыке регулярного и нерегулярного (узкие стенки образованы вторичными волноводами) прямоугольных волноводов. На основе численных результатов приводятся и анализируются зависимости дифракционных характеристик системы от ее параметров.
Ключевые слова: дифракция, импедансные стенки, проекционный метод, прямоугольный волновод.
ВВЕДЕНИЕ
Связь волноводов через отверстия различной формы широко используется при создании сложных волноводных узлов техники СВЧ. Примером такого соединения волноводов различного сечения через отверстия связи может быть поглощающий фильтр гармоник, описанный в [1, 2]. Расчет электрических характеристик такого фильтра весьма затруднен ввиду большого числа вспомогательных волноводов и отверстий связи. Это обстоятельство вызывает необходимость создания расчетной модели, которая позволила бы провести анализ факторов, влияющих на параметры такого рода структур. Известные до настоящего времени методы расчета поглощающих фильтров основаны на применении элементов теории цепей [2], а представленные результаты выполнены в одномодовом приближении.
В представленной работе приводится решение дисперсионного уравнения для нерегулярного волновода с импедансными узкими стенками, которое позволяет в многомодовом приближении решить задачу дифракции основной волны на скачке импеданса узких стенок. В качестве математического аппарата используется импедансный метод, развитый в работах [3, 4] и успешно применяемый для решения широкого класса задач дифракции и распространения радиоволн [5, 6]. Использование эквивалентных граничных условий импедансного типа дает возможность
© Логачева Л. М., Куцак С. В., Бондарев В. П., 2010
заменить последовательность ответвляющихся по узкой стенке вторичных волноводов гладкой поверхностью, на которой задано граничное условие Щукина - Леонтовича, и при этом исключить из рассмотрения поле в боковых волноводах. В общем случае введенный поверхностный импеданс является комплексной величиной, что позволяет учесть потери в поглощающем фильтре.
Решение задачи дифракции проводится методом частичных областей с привлечением проекционного метода, который успешно применяется при решении самых разнообразных волноводных задач [7, 8]. Используемый нами вариант проекционного метода аналогичен развитому в работе [9], поскольку и в данном случае в качестве базисных функций, в отличие от других вариантов метода, используется полная система ортогональных функций, адекватных рассматриваемой задаче.
1. ПОСТАНОВКА ЗАДАЧИ И МЕТОД
РЕШЕНИЯ
Рассмотрим прямоугольный волновод с идеально проводящими широкими стенками и узкими стенками, представляющими собой импеданс конечной длины I. Введем прямоугольную систему координат, начало которой совпадает с одной из вершин прямоугольника рис. 1.
Рис. 1. Волновод с импедансными узкими стенками
Примем следующие обозначения: а - размер широкой стенки основного волновода, Ь - размер узкой стенки основного волновода. Считаем, что на вход волновода падает основной тип волны Н10.
Первоначально задача сводится к записи дисперсионного уравнения, которое позволит определить спектр возбуждаемых волн в прямоугольном волноводе с узкими импедансными стенками. Полученный набор собственных поперечных чисел и функций позволит решить дифракционную задачу и проанализировать значения коэффициентов отражения и прохождения.
Для определения постоянных распространения запишем решение уравнения Гельмгольца для данной системы (зависимость от времени в виде множителя еш опускается):
Ey = Asinyx + B cosyx,
Hz =
(yA cosyx - yBsinyx),
(І)
(2)
где у
поперечное волновое число;
2n
к = — - волновое число свободного пространства; А
в - продольное волновое число; А - длина волны в свободном пространстве; ю - частота волны в свободном пространстве; ц0 - магнитная проницаемость вакуума; А, В - постоянные коэффициенты.
Учитывая граничное условие Щукина - Леонтови-ча на узких импедансных стенках:
[ П, E] = Zs[ П[ П, H]],
(З)
где п - нормаль, направленная вглубь импедансной стенки; 25 - в общем случае комплексный импеданс стенок, получим систему двух однородных уравнений:
B = -jZ°0^-;
ЮД0
jZ “y
A sinya + Bcosya = J—LL (A cosya - Bsinya), ЮД0
(4)
лу0 луа
где - значения поверхностного импеданса при
х =0 и х = а, соответственно.
Исключив из системы (4) постоянные коэффициенты, приходим к дисперсионному уравнению вида
у (2° + 2“,) соБуа + у (1 + у2.?0 ^а) Бшуа = 0, (5)
У
где используется следующая нормировка: 7 = ,,
к
2° = ?-, 2“, = 2; 20 = - сопротивление свободна 20 V е0
ного пространства.
Решение уравнения (5) позволяет определить набор поперечных собственных чисел ут в нерегулярном волноводе и соответствующие им продольные волновые числа
вт = л/к2 - 7л,
где т = 1, 2, 3, ... - номер корня.
Из системы (4) определяются собственные поперечные функции нерегулярного волновода, которые, с точностью до постоянного множителя, можно записать в виде
Фт(X) = Бшух -у^уСОБух. (6)
2. СИСТЕМА ФУНКЦИОНАЛЬНЫХ УРАВНЕНИЙ ВОЛНОВОДА
Предположим, что на границу раздела (2 = 0) в положительном направлении оси 2 со стороны регулярного волновода падает волна Н10 с компонентой
(7)
. ( П x) -j'ГlZ Ey = E0sin ^ —J e ,
где Е0 - амплитуда падающей волны; Гі = ^к2 -
продольное волновое число падающей волны.
Составляющие поля в различных областях волновода (рис. 1) можно записать в виде:
j-'i І) і о * mnx Jrmz , p.-* /Qv
Ey = Ey + у Rmsin----------------e , (z < 0); (S)
a
m = І
Ef] = у (Amejm + Ame^-"^m(x)
m = І
(0 < Z < l);
-j'PmZ , 4- „j'ßm( z - l\
(9)
. mnx -J^m(z-1)
У Tmsin—e m , (z < l); (І0)
a
2 ( m n
a
где Гт = к -1 — І - продольное волновое число
в регулярном волноводе; вт - продольное волновое
m = І
число в нерегулярном волноводе; Ят, Тт, Ат, Ат -постоянные коэффициенты.
Индексы (1), (2) и (3) обозначают соответствующие области волновода на рис. 1.
Для компоненты Нх имеем следующее представление поля:
НХГ) = НХ + — V Г я 8Іп—еГш\
ЮЦо
Н2 = V (Ате
£і Ю^о
(^ < 0); (11)
^2 - ^ - А+те ^ )Ф— (X ),
(0 < 2 < і); (12)
«ДО -
(2 < I).
(13)
Последовательность коэффициентов Ят и Тт подлежит определению.
Используя условие сшивания полей на границах раздела (2 = 0 и 2 = I) регулярного и нерегулярного волноводов, приходим к следующим системам функциональных уравнений относительно искомых коэффициентов Ят, Тт, Ат , Ат:
І (ЯтЭт—х- (А+т + Ате^т‘)Фт(х)| = -Е08Іп—X; (14)
т =1 а а
а
~ ґ тп , -/в і \ пх
І (Я-Г-^ІП—X + (Ат + Ат Є — ')втФт(х)) = ЕоГ^ІД—,
т =1
а
І ([л--е*т + Ат^Фт(х) - Т-зіП"—“хІ = 0;
т = 1
]Т ((а- + А+тЄ^ті) ■ РтФт(х) + Т—Г—ЙП-^) = 0.
т = 1 а
(15)
І (Я—0—„ - (А+т + А- е т )Г„т) = -Е>8„1;
-/в—'
І (Я—Г—О—п + (А- - А- Є — )в-^„т) = Е0ГА1;
(17)
-/в і
£ (Тт8т„ - (Ате т + Ат Ж™) = 0;
т = 1
£ (ТтГтКп + (Ат -А;^1")вт^™) = 0,
т=1
2а п
где ¥пт = -1 Фт (х ёх, п= 1, 2, 3, ....
а0 а
Складывая поочередно второе уравнение с четвертым и первое с третьим в системе (17), получим:
~ --/в I
£ ((Ят + Тт)-Гт -8-п + Ат + Ат )( 1-е "^вт^™) =
т=1
= Е0ГЛ1; (18)
~ л. -/в I
£ ((Ят + Тт)-8т„ - А + Ат )( 1+ е -)■ ¥„т) = -).
т=1
Отнимая поочередно от первого уравнения третье и от второго уравнения четвертое в системе (17), получим:
-/в і
£ ((Я— - Т-) ■ 8—„ + (А- + А-)( 1 - Є - ) ■ ^„-) = -)
(19)
-/в-К
І ((Я- - Т-) ■ Г— ■§-„ - (А- - А-)( 1+ Є - )в-Лт) =
т = 1
= Е0Г10„1.
Введя обозначения
тр і гр ________ _ ґ л + . л - ___ ґ •» - вт1 __ і ,
Ят Т т Хт; I Ат Ат Ст; I 1 Є ^т;
Я— Т— У—, IАт -Ат = О—,
1 о. -/в—1
1 + е = gm,
Заметим, что входящие в системы (14) и (15) ко- получим следующие системы:
—тх
ординатные функции 8Іп—а“ образуют на интервале
[0, а] ортогональную систему, и для нее справедливо следующее соотношение:
г • птх . кпх , а о І 8ІП-------8ІП-------ах = - 0-
■>д а а 2
(„ = 1, 2, 3, ...),
(16)
где 8тп - символ Кронекера.
Условие (16) легко позволяет трансформировать системы функциональных уравнений (14)-(15) в бесконечную систему линейных неоднородных алгебраических уравнений второго рода:
(ХтГт0т„ + Ст^т^т^„т) Е0Г10„1;
(Хт 0т„ CmgmF„m) Е00„1,
(ут 0т„ + ^т^т^пт') Е00„1;
(утГт0т„ DmgmвmF„m) Е0Г10„1,
(20)
т = 1
т = 1
т = 1
т = 1
а
т = 1
т = 1
т = 1
которые, после исключения коэффициентов хт и ут, сводятся к системам уравнений относительно коэф -фициентов Ст и Вт.
£ Ст ^тГп + 1т вт ] ' ?пт = 2Е0Гп5п1,
т = 1
(п = 1, 2, 3.), (22)
£ От [Т^т + вmgm ] ' ?пт = ^ЕдГАь
т =1
(п = 1,2,3, ...). (23)
Коэффициенты Яп и Тп определяются из следу-
ющих выражений:
Яп = £(А+т + Ате^'вт/)Епт - Е08п1,
т=1
(п = 1, 2, 3, ...), (24)
Тп = £ (Ате^ + Ат )Е„т,
т=1
(п = 1, 2, 3, ...). (25)
Очевидно, что матричные элементы и свободные члены систем уравнений (20) и (21) удовлетворяют условию квадратичной сходимости, так что решение последних может быть получено с любой степенью точности методом редукции [6], причем это решение единственно.
m і
Рис. 2. Зависимость модуля коэффициента прохождения |Г| волны Н10 от отношения а/Х
при параметрах l/a =10 и Z° = Z = Zs:
кривая 1 для Zs = 0, 05j; 2 - для Zs = 2j; 3 - для Zs = -2j
3. АНАЛИЗ ЧИСЛЕННЫХ РЕЗУЛЬТАТОВ
Системы уравнений (20) и (21) исследовались при
различных параметрах волновода: а/А, I/а, 2°, 2°. Устойчивость решения и сходимости алгоритма приводятся в табл. 1. Результаты приводятся для симметричной структуры при 2° = /, 2а = /. Для оценки сходимости решения рассматриваются три значения порядка редукции N = 3, 9, 12.
Приведенные данные свидетельствуют о быстрой сходимости алгоритма для симметричной структуры. В случае несимметричной структуры сходимость ухудшается из-за появления четных типов волн. Некоторые результаты численного анализа исследуемой структуры изображены на рис. 2-5.
Ввиду того, что исследуемая система симметрична, в волноводе будут распространяться только нечетные моды.
Анализ численных результатов для модуля коэффициента прохождения | Т основной волны Н10 (рис. 2) показывает, что в случае индуктивного реактивного импеданса (23 = 2/) значение модуля коэффициента прохождения незначительно изменяется с ростом частоты. При емкостном импедансе 2 = -2/) с ростом
m
Рис. 3. Модуль коэффициента прохождения волны Н30 при параметрах l/a =10 и Z° = Z = Zs
Таблица 1. Зависимость модуля коэффициента прохождения от порядка редукции
Параметры системы Модуль коэффициента прохождения
l/a a/A N = 3 N = 9 N = 12
5 0,8 0,9999787 0,9997871 0,9997194
1,2 0,9999977 0,9999955 0,999990б
10 0,8 0,9984375 0,999542 0,999б289
1,2 0,9994745 0,9995б1б 0,9995853
Рис. 4. Модуль коэффициента прохождения волны Н10 при параметрах l/a = 10 и Z Z = Zs
частоты наблюдаются резонансные явления, которые связаны с перераспределением энергии между гармониками внутри нерегулярной части волновода. При импедансе Zs = 0,05j модуль коэффициента прохождения практически равен 1, что приближает данную структуру к регулярному волноводу (идеальные стенки) и подтверждает правильность математической модели.
На рис. 3 приводятся зависимости модуля коэффициента прохождения волны Н30 от а/Х. Анализ зависимостей рис. 3 показывает, что при малых значениях поверхностного импеданса модуль коэффициента прохождения имеет относительно небольшие значения и практически не резонирует с ростом частоты. При Z° = Za = 2j наблюдаются резонансные явления, что при данном значении реактивности импеданса свидетельствует о влиянии высших типов волн. Максимальные значения модуля коэффициента прохождения волны Н30 при Z° = Zа = -2j свидетельствует об интенсивной трансформации волны Н10 в волну Н30.
На рис. 4 представлены модули коэффициента прохождения |Г| основной волны Н10 в случае, когда активная часть поверхностного импеданса не равняется нулю (узкие стенки с потерями). Анализ характеристик показывает, что с ростом частоты коэффициент прохождения также растет. Это связано со слабой трансформацией основной волны в высшие типы волн. Следует также отметить, что в случае индуктивной реактивности комплексного импеданса в диапазоне а/Х = 0,6 + 1,2 модуль коэффициента прохождения имеет постоянное минимальное значение.
Рис. 5. Модуль коэффициента прохождения волны Н30 при параметрах l/a = 10 и Z Z = Zs
На рис. 5 представлена зависимость модуля коэф-фициента прохождения |7| от a/Х для волны Н30. Из рисунка видно, что на величину коэффициента прохождения существенное влияние оказывает тип реактивности. При индуктивном характере Zs модуль коэффициента прохождения |Г| с увеличением a/X превышает модуль коэффициента прохождения |Г| при емкостной реактивности.
ВЫВОДЫ
В настоящей работе методом сшивания получено решение для задачи дифракции основной волны H10 на стыке регулярного и нерегулярного волновода. В рассмотренном варианте проекционного метода в качестве базисных функций использовались собственные функции регулярного и нерегулярного волноводов. Численное решение полученного дисперсионного уравнения позволило определить значения поперечных чисел и постоянных распространения нерегулярного волновода.
Приведенные характеристики коэффициента прохождения свидетельствуют о весомом вкладе высших типов волн в энергию поля системы, что подтверждает адекватность рассмотрения многомодового режима для такого рода структур.
В свою очередь, представленная импедансная модель поглощающего фильтра дает возможность анализировать различные апертуры вспомогательных волноводов и на основе моделирования провести оптимизацию характеристик фильтра.
СПИСОК ЛИТЕРАТУРЫ
1. Вольман В. И. Фильтр поглощающего типа на высокий уровень мощности / В. И. Вольман, Л. М. Логачева // Радиотехника. - 1979. - Т. 34, № 1. - С. 25-28.
2. МатейД. Л. Фильтры СВЧ, согласующие цепи и цепи связи / Д. Л. Матей, Л. Янг, Е. М. Т. Джонс ; [пер. с англ. под ред. Л. В. Алексеева]. - М.: Связь, 1972. -Т. 2. - 495 с.
3. Марков Г. Т. Возбуждение электромагнитных волн / Г. Т. Марков, А. Ф. Чаплин. - М. ; Л. : Энергия, 1967. -367 с.
4. Миллер М. А. Исследование понятия поверхностного импеданса в теории поверхностных волн (обзор) / М. А. Миллер, В. И. Таланов // Изв. вузов. Радиофизика. - 1961. - Т. 4, № 5. - С. 795-830.
5. Арсеньев Е. М. Поверхностные волны над трапециевидно гофрированной поверхностью / Е. М. Арсеньев, В. А. Калошин // Журнал радиоэлектроники. -2008. - № 4. - С. 127-135.
6. Горобец Н. Н. Приближенное решение дисперсионного уравнения для цилиндрического волновода с импе-дансными границами / Н. Н. Горобец, А. П. Удовенко // Радиофизика и радиоастрономия. - 2005. - Т 10, № 2. -С. 172-177.
7. Петрова Ю. Ю. Решение уравнения Гельмгольца в многосвязных волноводных областях: дис. ... канд. физ.-мат. наук : 01.01.03 / Петрова Юлия Юрьевна. -М., 2006. - 114 с.
8. Моденов В. П. Проекционные методы в теории волноводов / В. П. Моденов // Радиотехника и электроника (Россия). - 2005. - Т. 50, № 2. - С. 203-207.
9. Богданов Ф. Г. Дифракция волны Н10 на симметричных диэлектрических стержнях конечной длины / Ф. Г. Богданов // Изв. вузов. Радиофизика. - 1983. -Т. 26, № 2. - С. 246-250.
Надійшла 02.06.2009 Після доробки 03.09.2009
Логачова Л. М, Куцак С. В., Бондарєв В. П.
ДИФРАКЦІЯ ХВИЛІ #10 В ПРЯМОКУТНОМУ ХВИЛЕВОДІ НА СТРИБКУ ІМПЕДАНСА ВУЗЬКИХ СТІНОК
На основі проекційного методу дається строге рішення задачі дифракції основної хвилі Н10 на стику регулярного і нерегулярного (вузькі стінки утворені вторинними хвилеводами) прямокутних хвилеводів. На основі чисельних результатів приводяться та аналізуються залежності дифракційних характеристик системи від її параметрів.
Ключові слова: дифракція, імпедансні стінки, проекційний метод, прямокутний хвилевід.
Logacheva L. M., Kutsak S. V., Bondaryev V. P.
Н10 WAVE DIFFRACTION IN A RECTANGULAR WAVEGUIDE BY NARROW WALLS IMPEDANCE BUMP
Rigorous solution of the principal H10 wave diffraction problem at the joint of regular and non-regular (narrow walls are formed by secondary waveguides) rectangular waveguides has been obtained using the method of projection matching. On the basis of numerical results the dependence of the system diffraction characteristics on its parameters is determined and analyzed.
Key words: diffraction, impedance walls, method of projection matching, rectangular waveguide.
УДК 621.372.81
Самойлик С. С.1, Бондарев В. П.2
1 Аспирант Запорожского национального технического университета 2Канд. физ.-мат. наук, доцент Запорожского национального технического университета
ЭНЕРГЕТИЧЕСКИЕ ХАРАКТЕРИСТИКИ ПРОХОДНОГО РЕЗОНАТОРА НА ПРЯМОУГОЛЬНЫХ ВОЛНОВОДАХ С ЧАСТИЧНЫМ ДИЭЛЕКТРИЧЕСКИМ ЗАПОЛНЕНИЕМ
Расчет параметров проходного прямоугольного резонатора с частичным диэлектрическим заполнением при сильной связи выполнен методом частичных пересекающихся областей с использованием интегральных уравнений макроскопической электродинамики. Этот метод позволяет учитывать как параметры каждого из диэлектриков, такие как диэлектрическая проницаемость и тангенс угла потерь вещества, из которого изготовлена неоднородность, так и его геометрическую форму и место расположения в резонаторе. Проведенный численный анализ позволил вычислить мощности, поглощенные неоднородностями, амплитуды прошедшей через структуру и отраженной волн в зависимости от параметров сложного диэлектрического заполнения.
Ключевые слова: проходной резонатор, коэффициент отражения, коэффициент прохождения.
ВВЕДЕНИЕ ний, а необходимо математическое моделирование
Моделированию процессов, происходящих в СВЧ этих щюцесс°в. °дн°й из теркмых пр°блем при раз-
резонаторах, заполненных диэлектрическим матери- работке таких устройств является трудность формиро-
алом, посвящены исследования [1, 2]. Обзор работ по- вания заданных распределений поля в диэлектри-
казывает, что для детального анализа наблюдаемых ческом материале, помещенном в резонатор. Наиболее
явлений недостаточно экспериментальных исследова- предпочтительным является аналитическое решение
© Самойлик С. С., Бондарев В. П., 2010
