Другие журналы

научное издание МГТУ им. Н.Э. Баумана

НАУКА и ОБРАЗОВАНИЕ

Издатель ФГБОУ ВПО "МГТУ им. Н.Э. Баумана". Эл № ФС 77 - 48211.  ISSN 1994-0408

Анализ дисперсионных характеристик погонных параметров элементарных отрезков сверхширокополосных балансных печатных щелевых антенн бегущей волны

# 04, апрель 2013
DOI: 10.7463/0413.0548171
Файл статьи: Виленский_P.pdf (1480.17Кб)
авторы: Виленский А. Р., Чернышев С. Л.

УДК.  621.396.677

Россия, МГТУ им. Н.Э. Баумана

chernshv@bmstu.ru

 

Введение

Сверхширокополосные (СШП) печатные щелевые антенны бегущей волны (ПЩАБВ) в настоящее время широко применяются как в качестве одиночных направленных излучателей, так и в составе антенных решёток разнообразных СШП радиоэлектронных систем. Распространённым типом ПЩАБВ являются так называемые балансные ПЩАБВ, например балансные антенны Вивальди [1, 2]. В таких антеннах электромагнитная энергия излучаемой прямой бегущей волны концентрируется в основном в области щелевого зазора, формируемого двумя слоями металлизации несущей диэлектрической подложки. В связи с тем, что обычно ПЩАБВ применяются для излучения сверхкоротких СШП импульсов, актуально рассматривать вопрос об искажении формы излучаемого сигнала внутри структуры антенны. Одним из основных источников таких искажений служит дисперсия погонных характеристик элементарных отрезков ПЩАБВ, рассматриваемой с позиции теории плавных нерегулярных линий передач (НЛП). Электродинамический анализ линий передач (ЛП) с неоднородным диэлектрическим заполнением свидетельствует о том, что в общем случае в таких ЛП могут распространяться только гибридные волны [3]. Гибридные волны обладают дисперсией постоянной распространения и характеристического сопротивления.

В настоящей работе проведён электродинамический анализ дисперсионных характеристик погонных параметров отрезков ЛП балансных ПЩАБВ различных конфигураций. Дисперсионные уравнения получены сведением исходной системы интегральных уравнений, связывающей токи и напряженности электрического поля в поперечном сечении ЛП, к однородной системе линейных алгебраических уравнений (СЛАУ) при помощи спектрального метода Галёркина [4]. Для отыскания компонент двумерной тензорной функции Грина используется метод спектральных иммитансов [5]. Требование равенства нулю определителя системы позволяет отыскать значение постоянной распространения, которая затем используется для вычисления распределения напряженностей электромагнитного поля (ЭМП) в поперечном сечении ЛП и расчёта волнового сопротивления. Следует отметить, что спектральный метод Галёркина ранее применялся для анализа классических щелевых ЛП, результаты могут быть найдены в работах [6, 7].

1. Получение характеристического уравнения. Математическая постановка задачи

А. Описание полей в спектральной области

На рисунке 1,а показано поперечное сечение исследуемой ЛП. Далее будем называть такую ЛП двухполосковой. На рисунке толстыми линиями обозначена металлизация подложки, которую мы в дальнейшем будем считать бесконечно тонким идеально проводящим слоем. Диэлектрическая подложка имеет толщину h и относительную диэлектрическую проницаемость εr. Данная волноведущая структура являются открытой как в направлении оси х, так и в направлении оси у. Две границы раздела сред вакуум-диэлектрик имеют координаты y = ±h/2 соответственно.

а)                                                                    б)

Рисунок 1 – Поперечные сечения двухполосковой ЛП (а) и базовой слоистой структуры (б)

Рассмотрим базовую слоистую структуру – рисунок 1б. На границах раздела сред 1-2 и 2-3 в общем случае существуют произвольные распределения электрических поверхностных токов − Jx,z(x,y), которые обозначены как Jx,z(x,y=-h/2) = Jx,z1(x) и Jx,z(x,y=h/2) = Jx,z2(x).

Рассматриваемое поперечное сечение будем считать сечением регулярной ЛП, в которой вдоль оси z распространяется электромагнитная волна, с постоянной распространения β = β(f), где f − частота электромагнитного поля. Запишем однородные уравнения Гельмгольца для продольных составляющих векторов напряженности электрического и магнитного полей для каждой из подобластей в виде [8]:

 ,  i = 1,2,3,                             (1)

где поперечный оператор, ,,.

Поперечные компоненты поля в известном виде выражаются через продольные из уравнений Максвелла в дифференциальной форме [8].

Для рассматриваемой двумерной задачи поверхностные токи связаны с электрическим полем структуры через интегральный оператор, ядром которого является двумерная функция Грина многослойной среды [8]:

 

 

где Z(q1q2)(ij)(x) −  компоненты двумерной функции Грина, связывающие q1-компоненту вектора поля на границе раздела i с q2-компонентой вектора поверхностного тока на границе j. Здесь для определённости границе раздела y= h/2 присвоен номер 1, границе = -h/2 – номер 2.

Выражения (2) фактически определяют решение задачи отыскания полей через свёртку распределения токов с функцией Грина. В связи с этим удобно перейти в спектральную область, воспользовавшись преобразованием Фурье по координате х в виде:

,                                                                                         (3)

где черезобозначено преобразование Фурье оригинала f(x).

Тогда, используя свойства преобразование Фурье от свёртки двух функций, выражения (2) могут быть переписаны в спектральной области в виде:

 

где символ "^" обозначает преобразование Фурье исходной функции.

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

 ,     ,                            (5)

где скалярные функции  и определяются следующим образом:

 

где .

Неизвестные коэффициенты ...определяются из граничных условий на границах раздела 1 и 2:

Граница y = -h/2:

Граница y = h/2:

 

Таким образом, при известном значении β возможно отыскать ЭМП структуры. Для получения дисперсионного уравнения нам понадобятся спектральные выражения для компонент функции Грина.

Б. Получение компонент функции Грина

В общем случае компоненты функции Грина могут быть найдены при использовании граничных условий (7). Однако такой подход оказывается вычислительно ёмким. Удобно воспользоваться методом спектральных иммитансов для отыскания искомых выражений [5].

 

 

где                                 

Оставшиеся функции в (8) могут быть найдены в работе [5].

В. Применение спектрального метода Галёркина

Принимая во внимание центральную симметрию задачи, справедливы следующие выражения:

 

 

Таким образом, исходная система уравнений может быть сокращена вдвое. Далее, в соответствии с методом Галёркина, зададимся аппроксимационными выражениями для компонент поверхностных токов, например, на границе раздела 1:

 

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

Переходя в спектральную область в (12), получаем:

 

 

Требуя в соответствии с методом Галёркина равенства нулю скалярного произведения невязки системы с каждой из базисных функций, получаем СЛАУ относительно коэффициентов разложения:

                                                                 (12)

где однородность системы (12) следует из теоремы Парсеваля, а элементы блочной матрицы и компоненты векторов определяются следующим образом

 

К дисперсионному уравнению относительно β0 переходим, требуя равенства нулю определителя системы (12).

                                 (14)

После отыскания β0 можно положить  = 1 А/м. Затем в матрице (12) выделяется базисный минор порядка Nx+Nz-1 и рассчитываются оставшиеся коэффициенты. Далее по найденному распределению поверхностных токов рассчитываются спектральные распределения компонент ЭМП структуры. Наконец, отыскав распределения ЭМП, можно переходить к вычислению характеристического сопротивления двухполосковой ЛП:

  где                             (15)

 ‑ активная мощность, переносимая волной через поперечное сечение ЛП.

Г. Базисные функции

Рассмотрим вопрос выбора типа и количества базисных функций. Для решения рассматриваемой задачи целесообразно пользоваться базисами полной области, функции которых компактно распределены в спектральной области [9]. Двухполосковая ЛП не обладает симметрией в распределении поверхностных токов по полоскам. В связи с этим используемый базис, вообще говоря, должен быть составлен из набора чётных и нечётных относительно центра полоска функций.

Из общефизических соображений в качестве априорных можно отметить только условия на рёбрах полосков:

              (16)

где x0 = (S + W)/2.

Критерием отбора наиболее подходящей системы функций служило условие максимальной точности расчёта при минимальном размере системы. В результате вычислительным путём было найдено, что заданному критерию удовлетворяет система функций вида

 

 

где индексы «е» и «о» обозначают чётные (even) и нечётные (odd) базисные функции соответственно;  – количества чётных и нечётных базисных функций для x- и z-компонент тока соответственно. Преобразования Фурье от выражений (17) отыскиваются численно. На рисунках 2а, 2б приведены графики первых четырёх базисных функций для x- и z-компонент тока соответственно.

 

а)                                                                    б)

Рисунок 2 - Графическое представление первых четырёх базисных функций x‑компоненты тока (а), z- компоненты (б)

 

2. Результаты расчёта

В качестве примера рассмотрим ЛП, для которых S/W>0.02. В балансных ПЩАБВ таким конфигурациям соответствуют области, из которых начинает формироваться излучение. При относительно высоких S/W распределение токов по полоскам имеет слабую асимметрию, поэтому вычислительным путём было найдено, что при , mz = 8 ошибка вычисления β не превышает 2,5 %, а ошибка вычисления Z0 меньше 5 %. Разработанный вычислительный алгоритм был протестирован на подложках с низкими (случай антенн) и высокими (трактовые устройства) εr. Результаты расчёта коэффициента укорочения длины волны приведены на рисунках3 ‑ 5. Параметрами кривых служат значение диэлектрической проницаемости подложки εr, нормированная толщина подложки h/λ и нормированная ширина полоска W/λ. Стартовая точка на оси абсцисс S/λ в каждом случае соответствует значению 2h/λ. Аналогичные зависимости для характеристического сопротивления двухполосковой ЛП продемонстрированы на рисунках 6 ‑ 8.

 

 

Рисунок 3 – Зависимость коэффициента укорочения от нормированного ширины зазора S/λ; ε= 2,22, h/λ = 0,04

 

 

Рисунок 4 – Зависимость коэффициента укорочения от нормированного ширины зазора S/λ; ε= 3,38, h/λ = 0,005

 

 

Рисунок 5 – Зависимость коэффициента укорочения от нормированного ширины зазора S/λ; εr = 9,8, h/λ = 0,02

 

 

Рисунок 6 – Зависимость волнового сопротивления от нормированного ширины зазора S/λ; εr = 2,22, h/λ = 0,04

 

 

Рисунок 7 – Зависимость волнового сопротивления от нормированного ширины зазора S/λ; εr = 3,38, h/λ = 0,005

 

 

Рисунок 8 – Зависимость волнового сопротивления от нормированного ширины зазора S/λ; εr = 9,80, h/λ = 0,02

 

Проведённый анализ позволяет заключить, что наиболее сильно дисперсия постоянной распространения проявляется для геометрий ЛП с малыми значениями S/λ. В то время как двухполосковые ЛП с большими значениями S/λ обладают наибольшей дисперсией характеристического сопротивления. Как и в случае щелевых ЛП дисперсия усиливается при увеличении диэлектрической проницаемости подложки.

Полученные результаты также позволяют сделать другой интересный вывод. Анализ дисперсионных характеристик Z0 для случая электрически толстых подложек свидетельствует о стабилизации величины Z0, начиная с некоторого значения S/λ для широких полосков). Более того, при дальнейшем разнесении полосков наблюдается уменьшение величины Z0. Указанный эффект наблюдался также и при исследовании описанных конфигураций ЛП в САПР конечно-разностным методом. Физически эту закономерность можно объяснить тем, что при сильном разнесении широких полосков распределение поверхностных токов начинает симметрироваться и всё большая часть поперечного поля оказывается перехваченной подложкой.

В заключение сопоставим результаты, полученные с помощью квазистатического подхода [10, 11], с помощью САПР, реализующей конечно-элементный метод, и с помощью предложенного подхода. Объектами анализа выберем двухполосковые ЛП с высокими значениями S/W. Результаты продемонстрированы на рисунках 9 – 12.

 

 

Рисунок 9 – Сопоставление частотных зависимостей коэффициента укорочения; W = 5 мм, S = 5 мм, h = 0,5 мм, εr = 3,38

 

Рисунок 10 – Сопоставление частотных зависимостей характеристического сопротивления; W = 5 мм, S = 5 мм, h = 0,5 мм, εr = 3,38

 

Рисунок 11 – Сопоставление частотных зависимостей коэффициента замедления; W = 3 мм, S = 20 мм, h = 0,5 мм, εr = 9,8

 

Рисунок 12 – Сопоставление частотных зависимостей характеристического сопротивления; W = 3 мм, S = 20 мм, h = 0,5 мм, εr = 9,8

 

3. Выводы

В данной работе представлены результаты исследования дисперсионной модели двухполосковой ЛП. Переход в спектральную область позволил отыскать распределения компонент ЭМП в исследуемой структуре, а также построить эффективные вычислительные алгоритмы определения погонных характеристик ЛП при помощи спектрального метода Галёркина. Полученные данные свидетельствуют о том, что для рассмотренных конфигураций линий максимальное расхождение вычисленных значений волнового сопротивления в диапазоне (3–11) ГГц не превышает 5 %. Ошибка вычисления коэффициента укорочения во всех случаях не превысила 2,5 %.

Результаты анализа широких двухполосковых ЛП с высокими значениями межполоскового зазора (S) позволяют сделать вывод о сильной дисперсии волнового сопротивления таких линий. Следовательно, для построения аккуратной модели ПЩАБВ необходимо использовать дисперсионную модель.

Список литературы

1. Zhaofa Z., Guangyou F., Kamarudin M.F.B., Sato M. Analysis of balanced antipodal Vivaldi antenna and its application for mine-like targets detection // IEIC Technical Reports. 2002. Vol. 102, no. 48. P. 13-18.

2. Bourqui J., Okoniewski M., Fear E.C. Balanced antipodal Vivaldi antenna for breast cancer detection // Proc. 2nd Eur. Conf. Antennas Propag. Edinburgh, U.K., Nov. 11-16, 2007. 5 p.

3. Itoh T., Mittra R. Analysis of microstrip transmission lines. Illinois University: Antenna Lab., 1972. 47 p.

4. Itoh T, Mittra R. A technique for computing dispersion characteristics of shielded microstrip  lines (Short papers) // IEEE Trans. Microwave Theory Tech. 1974. Vol. MTT-22, no. 10. P. 896-898. DOI: 10.1109/TMTT.1974.1128375

5. Itoh T.T. Spectral domain immitance approach for dispersion characteristics of generalized printed transmission lines // IEEE Trans. Microwave Theory Tech. 1980. Vol. MTT-28, no. 7. P. 733-736. DOI: 10.1109/TMTT.1980.1130158

6. Janaswamy R., Schaubert D.H. Dispersion characteristics for wide slotlines on low-permittivity substrates (Short paper) // IEEE Trans. Microwave Theory Tech. 1985. Vol. MTT-33, no. 8. P. 723-726. DOI: 10.1109/TMTT.1985.1133064

7. Janaswamy R., Schaubert D.H. Characteristic impedance of wide slot line on low permittivity substrates // IEEE Trans. Microwave Theory Tech. 1986. Vol. MTT-34, no. 8. P. 900-902. DOI: 10.1109/TMTT.1986.1133465

8. МарковГ.Т., ЧаплинА.ФВозбуждениеэлектромагнитныхволнМ.-Л.: Энергия, 1967. 376 с.

9. Вычислительные методы в электродинамике : пер. с англ. / под ред. Р. Митры. М.: Мир, 1977. 485 с.

10. Моделирование малогабаритных сверхширокополосных антенн: коллективная монография / В.Б. Авдеев, А.В. Ашихмин, А.В. Бердышев, С.В. Корочин, В.М. Некрылов, А.В. Останков, Ю.Г. Пастернак, И.В. Попов, А.П. Преображенский; под ред. В.Б. Авдеева, А.В. Ашихмина. Воронеж: Воронежский государственный университет, 2005. 233 с.

11. Фуско В. СВЧ цепи. Анализ и автоматизированное проектирование: пер. с англ. М.: Радиоисвязь, 1990. 288 с.


Тематические рубрики:
Поделиться:
 
ПОИСК
 
elibrary crossref ulrichsweb neicon rusycon
 
ЮБИЛЕИ
ФОТОРЕПОРТАЖИ
 
СОБЫТИЯ
 
НОВОСТНАЯ ЛЕНТА



Авторы
Пресс-релизы
Библиотека
Конференции
Выставки
О проекте
Rambler's Top100
Телефон: +7 (915) 336-07-65 (строго: среда; пятница c 11-00 до 17-00)
  RSS
© 2003-2022 «Наука и образование»
Перепечатка материалов журнала без согласования с редакцией запрещена
 Тел.: +7 (915) 336-07-65 (строго: среда; пятница c 11-00 до 17-00)