Радиотехника и электроника, 2022, T. 67, № 9, стр. 835-846

Однополяризационная двухсекционная сверхширокополосная антенная решетка с экранированной первой секцией

С. Е. Банков *

Институт радиотехники и электроники им. В.А. Котельникова РАН
125009 Москва, ул. Моховая, 11, стр. 7, Российская Федерация

* E-mail: sbankov@yandex.ru

Поступила в редакцию 16.10.2020
После доработки 02.08.2021
Принята к публикации 25.12.2021

Полный текст (PDF)

Аннотация

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

1. ПОСТАНОВКА ЗАДАЧИ

Актуальной научно-технической задачей является создание сверхширокополосных (СШП) и сверхдиапазонных (СД) антенных решеток (АР) [1, 2]. Достижение рабочих диапазонов частот больших (8–10) : 1 при помощи традиционных технических решения в этой области сталкивается с рядом трудноразрешимых проблем, отмеченных в [3, 4]. В работе [5] в качестве альтернативной схемы построения СШП АР предложена двухсекционная (ДС) решетка.

Основная идея функционирования ДС АР основана на анализе собственных волн, существующих в бесконечных СШП-решетках, результаты которого изложены в [68]. Суммируя их можно сделать следующий вывод. Проблема расширения рабочего диапазона АР тесно связана с проблемой перехода с несимметричной линии передачи (ЛП) на симметричную. Под несимметричными ЛП мы понимаем такие структуры как коаксиальная, микрополосковая, полосковая и т.д. ЛП. Общей их чертой является то, что токи, текущие по проводникам в противоположных направлениях имеют разную структуру, которая обусловлена разной геометрией проводников, по которым они протекают. Например, в коаксиальной линии токи текут по внутреннему и внешнему проводникам, имеющим разные радиусы. Несимметричные ЛП удобны для построения широкого класса радиоэлектронных устройств. Поэтому элемент СШП АР должен иметь вход в виде такой ЛП.

Под симметричными ЛП обычно понимают такие структуры, как двухпроводную, щелевую и т.д. линии. В них существуют два одинаковых проводника расположенных симметрично относительно некоторой плоскости. Для построения радиоэлектронных устройств такие ЛП неудобны, так как они имеют существенно более высокий уровень радиационных потерь на некоторых нерегулярных элементах. Однако, по той же причине они удобны для построения излучателей, в том числе СШП излучателей. Наиболее распространенным примером может служить антенна Вивальди [2], которая представляет собой плавно расширяющуюся щелевую ЛП.

Переход с несиметричной ЛП на симметричную видимо принципиально включает в себя частотно-зависимые [2] элементы. Можно предложить также конструкцию перехода с диссипативными элементами, но такие технические решения встречаются редко и мы их обсуждать не будем. Наличие частотно-зависимых элементов ограничивает расширение частотного диапазона АР. Например, до определенного предела диапазон антенны Вивальди может расширяться за счет увеличения плавности изменения ширины щелевой линии и соответствующего ее удлинения. Однако, как только диапазон плавного трансформатора сравняется с диапазоном перехода с микрополосковой ЛП на щелевую, входящего в состав антенны, то общий диапазон устройства перестанет увеличиваться.

Таким образом, построение СШП переходов с несимметричной ЛП на симметричную становится решающим фактором для улучшения показателей качества СШП АР. В работе [5] предложен новая схема построения СШП АР, так называемая, ДС АР. Традиционная схема построения СШП АР содержит один трансформатор сопротивления, который обеспечивает изменение сопротивления от стандартного значения в 50 Ом до сопротивления близкого к 376 Ом, соответствующих волновому сопротивлению свободного пространства. Переход, связывающий несимметричную и симметричную ЛП, располагается в точке, где сопротивление близко к 50 Ом.

В ДС АР добавлен второй понижающий трансформатор сопротивления от стандартных 50 Ом до некоторого меньшего значения, которое, как правило, находится в диапазоне 15…25 Ом. В работе [9] показано, что размещение перехода в точке с пониженным сопротивлением существенно расширяет полосу его рабочих частот и делает возможным построение АР с диапазоном (20…30) : 1.

Отметим, что симметричные ЛП относятся к числу высокоомных линий. Так достижение характеристического сопротивления порядка 15…25 Ом в щелевой ЛП требует реализации весьма узких щелей, часто технологически недопустимых. По этой причине в работах [5, 9] в качестве симметричной ЛП взят плоский волновод (ПВ), в котором нужное значение сопротивления достигается относительно просто. В качестве несимметричной ЛП в указанных работах используется микрополосковая линия (МПЛ). При этом МПЛ в каналах АР расположены рядом друг с другом и между ними существует электромагнитная связь. Поэтому, правильно говорить не об отдельных ЛП, а об их системе.

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

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

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

2. ИССЛЕДУЕМАЯ СТРУКТУРА

На рис. 1а, 1б показаны два варианта конструктивного исполнения ДС АР. Первая секция выполнена в виде схемы на основе симметричной полосковой ЛП, у которой экранные проводники электрически связаны при помощи, например металлизированных отверстий. В этом случае осуществляется полная экранировка первой секции.

Рис. 1.

Конструкции двухсекционных решеток с экранированной первой секцией с горизонтально (а) и вертикально ориентированными (б) зондами.

Полосковую ЛП удобно выполнить на подложке с повышенной диэлектрической проницаемостью. При этом ЛП с низким характеристическим сопротивлением будут иметь относительно небольшую ширину.

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

Связь двух секций в варианте на рис. 1а осуществляется при помощи перехода. Он состоит из горизонтально ориентированного зонда, который является продолжением полоскового проводника, присоединенного к проводнику рупора. В состав перехода входит также полость, расположенная сзади рупора и ограниченная двумя металлическими вертикальными стенками. Она выполняет функцию объемного резонатора, который в идеальном случае должен иметь в точке расположения зондов бесконечно большое входное сопротивление. На практике оно зависит от частоты и является фактором, ограничивающим диапазон АР.

В конструкции на рис. 1а сокращение габаритных размеров происходит только за счет уменьшения длины первой секции в силу применения диэлектрической подложки.

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

3. ФОРМУЛИРОВКА ГРАНИЧНОЙ ЗАДАЧИ ДЛЯ ПЕРЕХОДА

На рис. 2 показана структура, для которой мы сформулируем и решим граничную задачу. Она представляет собой часть периода двумерной решетки, т.е. часть так называемого канала Флоке [10], к анализу которого сводится задача о бесконечной периодической решетке. Периоды решетки вдоль разных координатных осей равны ${{P}_{{y,z}}}$.

Рис. 2.

Модель перехода с полосковой линии на плоский волновод.

Рассмотрим граничные условия, которым удовлетворяет поле в структуре. На металлических поверхностях при $z = 0,h,H$, $x = 0, - L$ выполняются условия равенства нулю тангенциальных компонент электрического поля. На зонде, выполняется граничное условие следующего вида:

(1)
$U(y) = J(z){{R}_{g}} + E,$
где $E,{{R}_{g}}$ – ЭДС и внутреннее сопротивление сосредоточенного источника, включенного между зондом и нижним экраном плоского волновода, $U$ – напряжение между проводниками расположенными при $z = 0,h$:
(2)
$U = \int\limits_0^h {{{E}_{z}}(y,z)dz} ,$
$J$ – полный ток, текущий по проводнику:

(3)
$J(z) = \int\limits_{ - w/2}^{w/2} {{{I}_{z}}(y,z)dy} .$

Здесь ${{I}_{z}}(y,z)$ – плотность электрического тока на проводнике.

Выражение (1) является аналогом закона Ома и записано при ряде допущений. Мы предполагаем, что параметры $w,h$ удовлетворяют неравенствам

(4)
$kh,kw \ll 1,$
где k – волновое число свободного пространства. При выполнении соотношений (4) можем допустить, что плотность тока и электрическое поле постоянны по оси 0z. Тогда соотношение (1) преобразуется к следующему виду:

(5)
$h{{E}_{z}}(y) = J{{R}_{g}} + E.$

В формуле (5) сохранилась зависимость от координаты $y$. Чтобы избавиться от нее, проведем усреднение (5) на интервале $\left| y \right| < {w \mathord{\left/ {\vphantom {w 2}} \right. \kern-0em} 2}$:

(6)
$\frac{h}{w}\int\limits_{{{ - w} \mathord{\left/ {\vphantom {{ - w} 2}} \right. \kern-0em} 2}}^{{w \mathord{\left/ {\vphantom {w 2}} \right. \kern-0em} 2}} {{{E}_{z}}(y)dy} = J{{R}_{g}} + E.$

Сосредоточенный источник, параметры которого входят в формулу (6), позволяет описывать возбуждение перехода со стороны первой секции АР. Он выполняет функцию порта, при помощи которого переход соединяется с внешним устройством 1. Функцию порта 2 выполняет волновод, который присоединяется к структуре справа (см. рис. 2). Он имеет высоту $h$ и ширину ${{P}_{y}}$.

На поверхностях при $y = {{ \pm {{P}_{y}}} \mathord{\left/ {\vphantom {{ \pm {{P}_{y}}} 2}} \right. \kern-0em} 2}$ выполняются граничные условия периодичности:

(7)
$E({{ - {{P}_{y}}} \mathord{\left/ {\vphantom {{ - {{P}_{y}}} 2}} \right. \kern-0em} 2})\exp ( - i{{\kappa }_{y}}{{P}_{y}}) = E({{{{P}_{y}}} \mathord{\left/ {\vphantom {{{{P}_{y}}} 2}} \right. \kern-0em} 2}),$
где под $E({{ \pm {{P}_{y}}} \mathord{\left/ {\vphantom {{ \pm {{P}_{y}}} 2}} \right. \kern-0em} 2})$ понимаем любую компоненту электромагнитного поля. Здесь ${{\kappa }_{y}}$ – параметр периодического режима, определяющий сдвиг фаз между каналами решетки по оси $0y$. Поскольку решетка двумерная, то она имеет сдвиг фаз и по оси $0z$. Он описывается параметром ${{\kappa }_{z}}$. Указанные параметры связаны с углами сканирования решетки при помощи соотношений
(8)
${{\kappa }_{y}} = k\sin \theta \cos \varphi ,\,\,\,\,{{\kappa }_{z}} = k\sin \theta \sin \varphi ,$
где угол $\theta $ отсчитываем от оси $0x$, а $\varphi $ от оси $0y$.

Отметим, что допущение о постоянстве поля по оси $0z$, сделанное на основании неравенств (4), верно не только на поверхности зонда, но и во всей области $x > 0$.

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

4. РЕШЕНИЕ ГРАНИЧНОЙ ЗАДАЧИ

Для решения поставленной задачи разделим нашу сложную структуру на более простые области и запишем в них поле в форме, удовлетворяющей дифференциальным уравнениям, а также части граничных условий. В качестве частичных областей выберем область 1, расположенную при $ - L < x < 0$, область 2, которой соответствует интервал $0 < x < d$ и область 3, характеризуемую неравенством $x > d$.

Область 1 – это резонатор, область 2 – отрезок прямоугольного волновода длиной $d$ с поперечными размерами $h,{{P}_{y}}$ и область 3 – полубесконечный волновод с такими же поперечными размерами. Особенностью волновода, входящего в рассматриваемую структуру, являются граничные условия (7), установленные на его вертикальных стенках.

Используя неравенство (4), мы допустили, что поле в областях 2 и 3 не зависит от координаты $z$. Предположим также, что компонента электрического поля ${{E}_{z}}$ много больше компонент ${{E}_{{x,y}}}$. Такое допущение типично для представления поля между двумя близко расположенными металлическими проводниками.

Выразим далее поле в областях 1–3 через компоненты ${{E}_{z}},{{H}_{z}}$, т.е. в виде суперпозиции электрических и магнитных волн [11]. В рамках сделанных допущений нетрудно показать, что поле в областях 2 и 3 полностью описывается полем электрических волн, у которых ${{H}_{z}} = 0$. В области 1 мы не можем допустить постоянство поля по координате $z$, и поэтому должны учитывать вклад волн магнитного типа. В соответствии со сделанными выше замечаниями записываем следующие соотношения:

(9)
$\begin{gathered} {{E}_{{z1}}} = \sum\limits_{n,m} {A_{{n,m}}^{e}\cos {{\alpha }_{m}}z{\text{sh}}{{\gamma }_{{n,m}}}(x + L)\exp ( - i{{\kappa }_{n}}y)} , \\ {{H}_{{z1}}} = \sum\limits_{n,m} {A_{{n,m}}^{h}\sin {{\alpha }_{m}}z{\text{ch}}{{\gamma }_{{n,m}}}(x + L)\exp ( - i{{\kappa }_{n}}y)} , \\ {{E}_{{z2}}} = \sum\limits_n {\left( {{{B}_{n}}{\text{ch}}{{\beta }_{n}}x + {{C}_{n}}{\text{sh}}{{\beta }_{n}}x} \right)\exp ( - i{{\kappa }_{n}}y)} , \\ {{E}_{{z3}}} = \sum\limits_n {{{D}_{n}}\exp ( - i{{\kappa }_{n}}y - {{\beta }_{n}}(x - d))} , \\ {{\kappa }_{n}} = {{\kappa }_{y}} + {{2\pi n} \mathord{\left/ {\vphantom {{2\pi n} {{{P}_{y}}}}} \right. \kern-0em} {{{P}_{y}}}},\,\,\,\,{{\alpha }_{m}} = {{\pi m} \mathord{\left/ {\vphantom {{\pi m} H}} \right. \kern-0em} H}, \\ {{\gamma }_{{n,m}}} = \sqrt {\kappa _{n}^{2} + \alpha _{m}^{2} - {{k}^{2}}} ,\,\,\,\,{{\beta }_{n}} = {{\gamma }_{{n,0}}}, \\ n = \ldots - 1,0,1, \ldots ,\,\,\,\,m = 0,1,..., \\ \end{gathered} $
где $k$ – волновое число свободного пространства.

Соотношения (9) записаны так, что поля удовлетворяют уравнениям Максвелла, а также граничным условиям на металлических поверхностях и условиям периодичности. Для полного выполнения всех граничных условий нам осталось выполнить их в плоскостях $x = 0$ и $x = d$. Начнем с плоскости $x = 0$.

Электрическое поле при $x = 0$, $0 < z < h$ обозначим ${{E}_{z}}$. При $z > h$ функция ${{E}_{z}}$ равна нулю в силу граничного условия на металле. Другая тангенциальная к плоскости $x = 0$ компонента ${{E}_{y}}$ равна нулю при $x = 0$, $0 < z < H$ в силу упомянутых граничных условий. Ее равенство нулю при $0 < z < h$ следует из принято выше допущения относительно характера поля в областях 2 и 3.

Выразим компоненту электрического поля ${{E}_{y}}$ в области 1 через компоненты ${{E}_{z}},{{H}_{z}}$ при помощи известных соотношений [11] и приравняем ее нулю при $x = 0$, $0 < z < H$. Воспользуемся ортогональностью функций $\sin {{\alpha }_{m}}z$ на интервале $0 < z < H$ и ортогональностью функций $\exp ( - i{{\kappa }_{n}}y)$ при $\left| y \right| < {{{{P}_{y}}} \mathord{\left/ {\vphantom {{{{P}_{y}}} 2}} \right. \kern-0em} 2}$. В результате получаем следующее соотношение между коэффициентами $A_{{n,m}}^{{e,h}}$:

(10)
$A_{{n,m}}^{h} = - \frac{{{{\alpha }_{m}}{{\kappa }_{n}}}}{{k{{W}_{0}}{{\gamma }_{{n,m}}}}}A_{{n,m}}^{e},$
где ${{W}_{0}}$ – волновое сопротивление свободного пространства.

Выразим коэффициенты $A_{{n,m}}^{e}$ через электрическое поле ${{E}_{z}}$. Для этого приравняем представления (9) для компонент ${{E}_{{z1}}}$ и ${{E}_{{z2}}}$ при $x = 0$ функции ${{E}_{z}}$ и еще раз воспользуемся ортогональностью тригонометрических и экспоненциальных функций:

(11)
$\begin{gathered} {{A}_{{n,m}}} = \frac{{{{\varepsilon }_{m}}}}{{{{P}_{y}}H{\text{sh}}{{\gamma }_{{n,m}}}L}}\frac{{\sin {{\alpha }_{m}}h}}{{{{\alpha }_{m}}}}{{e}_{n}}, \\ {{e}_{n}} = \int\limits_{ - {{P}_{y}}/2}^{{{P}_{y}}/2} {{{E}_{z}}(y)\exp (i{{\kappa }_{n}}y)dy} , \\ {{\varepsilon }_{m}} = \left\{ \begin{gathered} 1,{\text{ }}m = 0, \hfill \\ 2,{\text{ }}m > 0. \hfill \\ \end{gathered} \right. \\ \end{gathered} $

В силу граничных условий при $x = 0$, $0 < z < h$ электрическое поле в области 1 должно равняться электрическому полю в области 2. Исходя из этого найдем коэффициенты ${{B}_{n}}$, выражая их через поле ${{E}_{z}}$:

(12)
${{B}_{n}} = {{{{e}_{n}}} \mathord{\left/ {\vphantom {{{{e}_{n}}} {{{P}_{y}}}}} \right. \kern-0em} {{{P}_{y}}}}.$

Таким образом, с помощью соотношений (11) и (12) мы выполнили все граничные условия для электрического поля при$x = 0$. Граничные условия для магнитного поля состоят в непрерывности компонент ${{H}_{{z,y}}}$ при $x = 0$, $0 < z < h$. Полностью выполнить их невозможно из-за того, что мы положили ${{E}_{y}} = 0$, поэтому удовлетворим их только для компоненты ${{H}_{y}}$.

Представим ее в областях 1 и 2 с помощью соотношений из [11]:

(13)
$\begin{gathered} {{H}_{{y1}}} = \sum\limits_{n,m} { - \frac{{i\cos {{\alpha }_{m}}z{\text{ch}}{{\gamma }_{{n,m}}}(x + L)\exp ( - i{{\kappa }_{n}}y)}}{{k{{W}_{0}}{{\gamma }_{{n,m}}}}}} \times \\ \times \,\,\left( {\kappa _{n}^{2} - {{k}^{2}}} \right){{A}_{{n,m}}}, \\ {{H}_{{y2}}} = \sum\limits_n { - \frac{{i{{\beta }_{n}}\exp ( - i{{\kappa }_{n}}y)}}{{k{{W}_{0}}}}\left( {{{B}_{n}}{\text{sh}}{{\beta }_{n}}x + {{C}_{n}}{\text{ch}}{{\beta }_{n}}x} \right)} . \\ \end{gathered} $

При выводе первого равенства $(13)$ мы воспользовались соотношением (10) и опустили верхний индекс в обозначении коэффициентов ${{A}_{{n,m}}}$. Введем обозначение:

(14)
${{Q}_{n}} = \frac{{\left( {\kappa _{n}^{2} - {{k}^{2}}} \right)}}{{{{k}^{2}}{{P}_{y}}Hh}}\sum\limits_m {\frac{{{{\varepsilon }_{m}}}}{{{{\gamma }_{{n,m}}}}}{\text{cth}}{{\gamma }_{{n,m}}}L{{{\left( {\frac{{\sin {{\alpha }_{m}}h}}{{{{\alpha }_{m}}}}} \right)}}^{2}}} .$

Для выполнения граничных условий приравняем соотношения (13) при $x = 0$, $0 < z < h$. Проинтегрируем полученное равенство по $z$ и воспользуемся ортогональностью функций $\exp ( - i{{\kappa }_{n}}y)$. В результате найдем коэффициенты ${{C}_{n}}$:

(15)
${{C}_{n}} = \frac{{{{k}^{2}}{{Q}_{n}}}}{{{{\beta }_{n}}}}{{e}_{n}}.$

Следующий этап решения – выполнение условий при $x = d$. Как следует из разд. 3, при $x = d$, $0 < z < h$, $\left| y \right| < {w \mathord{\left/ {\vphantom {w 2}} \right. \kern-0em} 2}$ протекает электрический ток, имеющий плотность ${{I}_{z}}$. Касательное электрическое поле непрерывно на листке с током, а магнитное терпит разрыв [11]. Отсюда получаем следующие соотношения:

(16)
$\begin{gathered} {{B}_{n}}{\text{ch}}{{\beta }_{n}}d + {{C}_{n}}{\text{sh}}{{\beta }_{n}}d = {{D}_{n}}, \\ {{D}_{n}} + {{B}_{n}}{\text{sh}}{{\beta }_{n}}d + {{C}_{n}}{\text{ch}}{{\beta }_{n}}d = {{{{j}_{n}}k{{W}_{0}}} \mathord{\left/ {\vphantom {{{{j}_{n}}k{{W}_{0}}} {i{{\beta }_{n}}{{P}_{y}}}}} \right. \kern-0em} {i{{\beta }_{n}}{{P}_{y}}}}, \\ {{j}_{n}} = \int\limits_{ - w/2}^{w/2} {{{I}_{z}}(y)\exp (i{{\kappa }_{n}}y)dy} . \\ \end{gathered} $

Подставим в формулы (16) выражения (12), (15):

(17)
$\begin{gathered} \frac{{{{e}_{n}}}}{{{{P}_{y}}}}{\text{ch}}{{\beta }_{n}}d + \frac{{{{k}^{2}}{{Q}_{n}}}}{{{{\beta }_{n}}}}{{e}_{n}}{\text{sh}}{{\beta }_{n}}d = {{D}_{n}}, \\ {{D}_{n}} + \frac{{{{e}_{n}}}}{{{{P}_{y}}}}{\text{sh}}{{\beta }_{n}}d + \frac{{{{k}^{2}}{{Q}_{n}}}}{{{{\beta }_{n}}}}{{e}_{n}}{\text{ch}}{{\beta }_{n}}d = \frac{{{{j}_{n}}k{{W}_{0}}}}{{i{{\beta }_{n}}{{P}_{y}}}}. \\ \end{gathered} $

Выразим с использованием системы (17) коэффициенты ${{e}_{n}}$ через${{j}_{n}}$:

(18)
${{e}_{n}} = k{{W}_{0}}\exp ( - {{\beta }_{n}}d){{{{j}_{n}}} \mathord{\left/ {\vphantom {{{{j}_{n}}} {i{{\beta }_{n}}}}} \right. \kern-0em} {i{{\beta }_{n}}}}\left( {1 + {{Q}_{n}}{{k}^{2}}{{{{P}_{y}}} \mathord{\left/ {\vphantom {{{{P}_{y}}} {{{\beta }_{n}}}}} \right. \kern-0em} {{{\beta }_{n}}}}} \right).$

Теперь мы можем найти электрическое поле ${{E}_{z}}$ при $x = d$ и приступить к выполнению граничного условия (6):

(19)
${{E}_{z}} = \sum\limits_n {\frac{{k{{W}_{0}}\left( {{\text{ch}}{{\beta }_{n}}d + {\text{sh}}({{\beta }_{n}}d){{Q}_{n}}{{k}^{2}}{{{{P}_{y}}} \mathord{\left/ {\vphantom {{{{P}_{y}}} {{{\beta }_{n}}}}} \right. \kern-0em} {{{\beta }_{n}}}}} \right)\exp ( - {{\beta }_{n}}d - i{{\kappa }_{n}}y){{j}_{n}}}}{{i{{\beta }_{n}}{{P}_{y}}\left( {1 + {{Q}_{n}}{{k}^{2}}{{{{P}_{y}}} \mathord{\left/ {\vphantom {{{{P}_{y}}} {{{\beta }_{n}}}}} \right. \kern-0em} {{{\beta }_{n}}}}} \right)}}} .$

Аппроксимируем ток на узкой полоске константой:

(20)
${{I}_{z}}(y) = C.$

Тогда результатом выполнения условия (6) будет определение постоянной C. Подставим выражение (19) в формулу (6) и найдем ее:

(21)
$\begin{gathered} C = \frac{E}{{\frac{{kh{{W}_{0}}}}{{2iw{{P}_{y}}}}\sum\limits_n {\left( {1 + \frac{{{{\beta }_{n}} - {{Q}_{n}}{{k}^{2}}{{P}_{y}}}}{{{{\beta }_{n}} + {{Q}_{n}}{{k}^{2}}{{P}_{y}}}}\exp ( - 2{{\beta }_{n}}d)} \right)\frac{{f_{n}^{2}}}{{{{\beta }_{n}}}} - w{{R}_{g}}} }},\,\,\,{{f}_{n}} = 2\frac{{\sin {{\kappa }_{n}}{w \mathord{\left/ {\vphantom {w 2}} \right. \kern-0em} 2}}}{{{{\kappa }_{n}}}}. \hfill \\ \hfill \\ \end{gathered} $

Теперь не составляет труда найти полный ток $J$ и входное сопротивление ${{Z}_{{i1}}}$ перехода по входу 1:

(22)
${{Z}_{{i1}}} = \frac{{ikh{{W}_{0}}}}{{2{{w}^{2}}{{P}_{y}}}}\sum\limits_n {\left( {1 + \frac{{{{\beta }_{n}} - {{Q}_{n}}{{k}^{2}}{{P}_{y}}}}{{{{\beta }_{n}} + {{Q}_{n}}{{k}^{2}}{{P}_{y}}}}\exp ( - 2{{\beta }_{n}}d)} \right)\frac{{f_{n}^{2}}}{{{{\beta }_{n}}}}} .$

По входному сопротивлению нетрудно определить коэффициент отражения от первого входа ${{S}_{{11}}}$:

(23)
${{S}_{{11}}} = \frac{{{{Z}_{{i1}}} - {{R}_{g}}}}{{{{Z}_{{i1}}} + {{R}_{g}}}}.$

Для полного описания перехода нам необходимо определить все элементы его матрицы рассеяния. Начнем с коэффициента передачи S21 со входа 1 на вход 2. Для этого нам нужно найти амплитуду волны в волноводе при $x > d$, которая возникает в нем при возбуждении входа 1. Это нетрудно сделать c помощью соотношений (9), (17), (18). Используя их, получаем выражение для амплитуды электрического поля прошедшей волны ${{E}_{{zt}}}$:

(24)
${{E}_{{zt}}} = - \frac{{iCk{{W}_{0}}{{f}_{0}}}}{{2{{P}_{y}}{{\beta }_{0}}}}\left( {1 + \frac{{{{\beta }_{0}} - {{Q}_{0}}{{k}^{2}}{{P}_{y}}}}{{{{\beta }_{0}} + {{Q}_{0}}{{k}^{2}}{{P}_{y}}}}\exp ( - 2{{\beta }_{0}}d)} \right).$

Здесь постоянная $C$ определяется формулой (21). Мощность, переносимая прошедшей на вход 2 волной ${{P}_{t}}$ , связана с полем простым соотношением –

(25)
${{P}_{t}} = - \frac{{i{{\beta }_{0}}{{{\left| {{{E}_{{zt}}}} \right|}}^{2}}h{{P}_{y}}}}{{2k{{W}_{0}}}},$

мощность волны, падающей на вход 1 –

(26)
${{P}_{i}} = {{{{{\left| E \right|}}^{2}}} \mathord{\left/ {\vphantom {{{{{\left| E \right|}}^{2}}} {8{{R}_{g}}}}} \right. \kern-0em} {8{{R}_{g}}}}.$

Тогда в соответствии с определением матрицы рассеяния [12] получаем

(27)
${{S}_{{21}}} = \sqrt {\frac{{{{P}_{t}}}}{{{{P}_{i}}}}} \exp (i\arg {{E}_{{zt}}}).$

Оставшиеся два элемента матрицы рассеяния ${{S}_{{12,22}}}$ найдем из свойств взаимного и недиссипативного четырехполюсника [12]:

(28)
$\begin{gathered} {{S}_{{12}}} = {{S}_{{21}}},\,\,\,\,\left| {{{S}_{{22}}}} \right| = \left| {{{S}_{{11}}}} \right|, \\ \arg {{S}_{{22}}} = - \arg {{S}_{{11}}} + 2\arg {{S}_{{12}}} + \pi . \\ \end{gathered} $

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

5. ЭКВИВАЛЕНТНАЯ СХЕМА ПЕРЕХОДА

Проведем анализ соотношения (22) и найдем с его помощью выражения, описывающие элементы эквивалентной схемы. Структуру эквивалентной схемы зададим исходя из физических соображений (рис. 3).

Рис. 3.

Эквивалентная схема перехода, 1 и 2 – точки включения полосковой линии и плоского волновода соответственно.

Входящие в эквивалентную схему элементы имеют простой физический смысл. Сопротивление ${{R}_{w}}$ – это характеристическое сопротивление волновода, образующего вход 2, а сопротивление ${{R}_{g}}$ – это сопротивление входа 1, которое также может быть характеристическим сопротивлением ЛП, формирующей этот вход. Реактивное сопротивление ${{X}_{L}}$, включенное последовательно с резистором ${{R}_{g}}$ описывает реактивное сопротивление зонда, введенного в волновод. Обычно оно имеет индуктивный характер, о чем свидетельствует индекс $L$. Реактивное сопротивление ${{X}_{r}}$ соответствует входному сопротивлению резонатора, приведенному к сечению $x = d$.

Рассматриваем структуру в одноволновом режиме, в котором в волноводе существует одна распространяющаяся волна. Такой режим имеет место при выполнении следующих неравенств:

(29)
$h < {\lambda \mathord{\left/ {\vphantom {\lambda 2}} \right. \kern-0em} 2},\,\,\,\,{{P}_{y}} < {\lambda \mathord{\left/ {\vphantom {\lambda {(1 + \sin \theta \cos \varphi )}}} \right. \kern-0em} {(1 + \sin \theta \cos \varphi )}},$
где $\lambda $ – длина волны в свободном пространстве.

Если соотношения (29) справедливы, то все постоянные ${{\beta }_{n}}$, за исключением ${{\beta }_{0}}$, являются действительными числами. Для параметра ${{\beta }_{0}}$ получаем

(30)
${{\beta }_{0}} = i{{b}_{0}},$
где ${{b}_{0}}$ – действительная величина.

С учетом соотношения (30) запишем выражение для входного сопротивления в новой форме:

(31)
$\begin{gathered} {{Z}_{{i1}}} = i{{X}_{L}} + a(1 + \cos \psi ) + ia\sin \psi , \\ a = \frac{{kh{{W}_{0}}}}{{2{{w}^{2}}{{P}_{y}}}}\frac{{f_{0}^{2}}}{{{{b}_{0}}}},\,\,\,\,{{X}_{L}} = \frac{{kh{{W}_{0}}}}{{2{{w}^{2}}{{P}_{y}}}} \times \\ \times \,\,\sum\limits_n {^{{(0)}}\left( {1 + \frac{{{{\beta }_{n}} - {{Q}_{n}}{{k}^{2}}{{P}_{y}}}}{{{{\beta }_{n}} + {{Q}_{n}}{{k}^{2}}{{P}_{y}}}}\exp ( - 2{{\beta }_{n}}d)} \right)\frac{{f_{n}^{2}}}{{{{\beta }_{n}}}}} , \\ \exp (i\psi ) = \frac{{i{{b}_{0}} - {{Q}_{0}}{{k}^{2}}{{P}_{y}}}}{{i{{b}_{0}} + {{Q}_{0}}{{k}^{2}}{{P}_{y}}}}\exp ( - 2i{{b}_{0}}d). \\ \end{gathered} $

Здесь верхний индекс (0) означает, что суммирование ведется по всем значениям n кроме нулевого. Между тем можно записать выражение для входного сопротивления ${{\tilde {Z}}_{{i1}}}$ по входу 1, используя параметры эквивалентной схемы:

(32)
${{\tilde {Z}}_{{i1}}} = i{{X}_{L}} + \frac{{i{{X}_{r}}{{R}_{w}}}}{{{{R}_{w}} + i{{X}_{r}}}}.$

Сравнивая выражения (31) и (32), получаем

(33)
${{R}_{w}} = 2a,\,\,\,\,{{X}_{r}} = 2a\frac{{1 + \cos \psi }}{{\sin \psi }}.$

Соотношения (31)–(33) полностью определяют все параметры эквивалентной схемы. Важно отметить, что ${{f}_{0}} = w$, поэтому

(34)
${{R}_{w}} = {{kh{{W}_{0}}} \mathord{\left/ {\vphantom {{kh{{W}_{0}}} {{{P}_{y}}{{b}_{0}}}}} \right. \kern-0em} {{{P}_{y}}{{b}_{0}}}}.$

Интересно, что в соотношение для сопротивления волновода (34) входят только его геометрические размеры и углы сканирования. Причем при сканировании в плоскости $\varphi = 90^\circ $ сопротивление не зависит от угла $\theta $, так как от этого угла не зависит постоянная ${{b}_{0}}$.

Знание параметров эквивалентной схемы оказывается весьма полезным при настройке перехода и физической интерпретации численных результатов.

6. ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ПЕРЕХОДА

Обсудим результаты численного исследования перехода. На рис. 4–7 представлены частотные зависимости коэффициента отражения перехода по входу 1. Коэффициент отражения выражен в децибелах. Результаты расчетов получены для нулевых углов сканирования $\theta ,\varphi $.

Рис. 4.

Частотная характеристика перехода при w = = 2.5 (1), 3.5 (2), 4.5 (3), 5.5 (4).

Рис. 5.

Частотная характеристика перехода при d = = 0 (1), 0.5 (2),1 (3), 1.5 (4).

Рис. 6.

Частотная характеристика перехода при L = = 16 (1), 14 (2), 12 (3), 10 (4), 8 (5).

Рис. 7.

Частотная характеристика перехода при H = = 10 (1), 12 (2) и 14 (3).

На рис. 4 зависимости получены для параметров L = 13.5, H = 14, d = 0.5, Py = 15, h = 1, w = 2.5, 3.5, 4.5, 5.5. Здесь и далее все размеры приведены в миллиметрах. Видно, что ширина зонда $w$ достаточно сильно влияет на уровень коэффициента отражения в центральной части рабочего диапазона перехода. Из результатов разд. 5 можно сделать вывод, что этот параметр наиболее сильно сказывается на индуктивности зонда, размещенного внутри волновода.

Нижняя граница рабочего диапазона почти не зависит от ширины $w$, поскольку на низких частотах сопротивление зонда стремится к нулю при любом значении индуктивности. Верхняя граница диапазона меняется при изменении $w$. Причем она располагается тем выше, чем шире зонд.

Для однозначного определения рабочего диапазона необходимо задать предельный уровень коэффициента отражения. Пусть он будет равен –10 дБ. В этом случае, как следует из рис. 4, рабочий диапазон $D \approx 25$.

На рис. 5 показано влияние расстояния d на частотную характеристику коэффициента отражения перехода при L = 13.5, H = 14, d = 0, 0.5, 1, 1.5, Py = 15, h =1, w = 5.5. Поведение кривых 1–4 похоже на поведение кривых на рис. 4. Видно, что положение точки размещения зонда слабо влияет на нижнюю границу диапазона и сравнительно сильно на верхнюю. Также от него зависит уровень коэффициента отражения в центральной части рабочей полосы частот.

На рис. 6 показана частотная зависимость коэффициента отражения перехода от длины резонатора L = 16, 14, 12, 10, 8 при H = 14, d = 0.5, Py = 15, h = 1, w = 4.5. Видно, что от данного параметра зависит как нижняя, так и верхняя границы рабочего диапазона частот. При этом увеличение длины $L$ уменьшает обе границы.

Поведение частотной характеристики коэффициента отражения при изменении размера резонатора показывает основные факторы, определяющие диапазон частот перехода. Рассогласование перехода на низких частотах связано с тем, что реактивное сопротивление резонатора при $f \to 0$ приобретает индуктивный характер и стремится к нулю. По этой причине к нулю стремится активная часть входного сопротивления перехода, а коэффициент отражения растет.

Следует отметить, что величина индуктивности определяется размерами резонатора: длиной $L$ и высотой $H$. При их увеличении индуктивность растет, а нижняя граница полосы рабочих частот смещается вниз.

Вместе с тем увеличение длины $L$ сказывается на поведении перехода на высоких частотах. Под высокими частотами понимаем частоты, на которых размеры резонатора сравнимы с половиной длины волны в свободном пространстве. В этой части диапазона в резонаторе возбуждается его основное колебание. На частоте резонанса входное сопротивление резонатора стремится к нулю и шунтирует активную часть входного сопротивления перехода. При этом наблюдается всплеск коэффициента отражения, определяющий верхнюю границу рабочего диапазона перехода. Появление резонанса на высоких частотах хорошо видно из кривой 1 на рис. 6.

На рис. 7 представлена зависимость характеристик перехода от высоты резонатора H = 10, 12, 14 при L = 14, d = 0.5, Py = 15, h = 1, w = 4.5. Влияние высоты аналогично влиянию его длины: увеличение данного параметра одновременно снижает нижнюю и верхнюю границы рабочего диапазона.

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

Рис. 8.

Модель перехода для электродинамического моделирования.

На рис. 9 представлена частотная характеристика перехода, полученная при L = 14, H = 14, d = 0.5, Py =15, h = 1, w = 4.5 с помощью приближенной и точной моделей. Видно достаточно хорошее совпадение кривых 1 и 2, отличие наблюдается в центральной части рабочего диапазона, в которой коэффициент отражения весьма мал. Границы диапазона приближенная модель описывает с весьма высокой точностью.

Рис. 9.

Частотные характеристики перехода, полученные с помощью приближенной (1) и точной (2) моделей.

7. МОДЕЛИ ПЕРВОЙ И ВТОРОЙ СЕКЦИЙ РЕШЕТКИ

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

При решении задачи синтеза будем опираться на результаты работы [13], в которой рассмотрена задача синтеза ступенчатого и плавного трансформатора сопротивления в диапазоне частот большем 10 : 1. В данной же работе ограничимся ступенчатыми трансформаторами. В рамках методики, развитой в [13], можно получить значения характеристических сопротивлений отрезков ЛП, формирующих трансформатор. Следующий этап синтеза – определение по заданным значениям сопротивлений геометрических размеров ЛП. Задача синтеза решается для некоторых фиксированных значений углов сканирования. В данной работе они равны нулю.

Для первой секции, выполненной на основе стандартной ЛП, например полосковой, определение геометрических параметров не составляет труда. Известно большое число программ [14], а также литературных источников [15], которые позволяют найти по заданному значению сопротивления геометрические размеры ЛП.

Относительно нестандартной задачей является определение геометрических параметров второй секции. Она представляет собой последовательность отрезков плоских волноводов, вместе образующих плоский рупор. В рамках модели АР в виде канала Флоке плоские волноводы, бесконечные по оси $0y$, ограничиваются стенками с условиями (7). Сопротивление такого волновода было найдено выше (34). Это соотношение решает задачу определения высоты плоского волновода по найденному в ходе синтеза характеристическому сопротивлению.

Отдельной задачей является определение характеристического сопротивления последней секции трансформатора – пустого канала Флоке. В работе [16] приведено соотношение для этого параметра.

Анализ двух секций АР проводим стандартными методами теории ЛП СВЧ. Поскольку трансформаторы сопротивления представляют собой цепочки последовательно соединенных четырехполюсников, то для определения их характеристик удобно воспользоваться аппаратом матриц передачи [12].

8. ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ РЕШЕТКИ В РЕЖИМЕ СКАНИРОВАНИЯ

Рассмотрим результаты численного исследования ДС АР, полученные с помощью описанной выше приближенной модели. Синтез первой и второй секций проводили при следующих условиях: нижняя и верхняя частоты рабочего диапазона 0.45 и 10 ГГц соответственно, максимальный уровень коэффициента отражения каждой секции в рабочем диапазоне частот 0.13, число скачков сопротивления в первом и втором трансформаторах 12 и 24 соответственно, длина второй секции 330. Переход имел следующие параметры: L = 14, H = 14, d = 0.5, Py,z = 15, h = 1, w = 4.5.

На рис. 10–13 представлены результаты численного исследования в виде частотных зависимостей коэффициента отражения АР – ${{S}_{{11}}}$. Кривая на рис. 10 получена при $\theta = 0$. Рисунки 11–13 получены для $\varphi = 0^\circ ,\,\,45^\circ ,\,\,90^\circ $ соответственно, и каждый для трех значений $\theta = 17^\circ ,\,\,33^\circ ,\,\,49^\circ $.

Рис. 10.

Частотная характеристика решетки при $\theta = 0$.

Рис. 11.

Частотная характеристика решетки при $\varphi = 0$ и θ = 17° (а), 33° (б) и 49° (в).

Рис. 12.

Частотная характеристика решетки при$\varphi = 45^\circ $ и θ = 17° (а), 33° (б) и 49о (в).

Рис. 13.

Частотная характеристика решетки при $\varphi = 90^\circ $ и θ = 17° (а), 33° (б) и 49° (в).

Представленные на рис. 10–13 данные позволяют дать оценку характеристикам сканирования обсуждаемой СШП АР. В режиме излучения по нормали к плоскости АР (см. рис. 10) антенна имеет диапазон частот 23.75, что даже несколько больше диапазона частот, использованного при синтезе трансформаторов, который равен 22.2. Расширение диапазона происходит несмотря на снижение верхней частоты до 9.5 ГГц за счет понижения нижней границы до 0.4 ГГц.

При увеличении угла места $\theta $ решетка ведет себя по-разному в разных плоскостях сканирования. Отметим, что плоскость $\varphi = 0$ часто называют плоскостью вектора магнитного поля, или $H$-плоскостью, плоскость $\varphi = 90^\circ $ – плоскостью вектора электрического поля, или $E$-плоскостью. Наконец, плоскость $\varphi = 45^\circ $ – диагональной плоскостью.

Наилучшие характеристики решетка демонстрирует в диагональной плоскости (см. рис. 12а–12в), в которой увеличение угла места слабо влияет на нижнюю границу рабочего диапазона, а верхнюю границу сдвигает в область высоких частот. Следует отметить, что при определении рабочего диапазона, в случае если верхняя его граница, определенная по критерию ${{S}_{{11}}} < - 10$ дБ, оказывается больше 10 ГГц, мы полагали ее равной 10 ГГц. Такое определение верхней границы связано с тем, что периоды решетки равны 15 и на частотах, больших 10 ГГц, в ней могут появляться при некоторых углах сканирования побочные дифракционные максимумы. Таким образом, максимальная рабочая частота анализируемой решетки не может быть больше 10 ГГц не в силу ее рассогласования, а из-за нарушения нормальных условий излучения.

Худшие характеристики сканирования решетка имеет в Е-плоскости, в которой увеличение угла сканирования приводит одновременно к увеличению минимальной частоты и уменьшению максимальной.

Зависимость рабочего диапазона D от угла $\theta $ для $\varphi = 0^\circ ,\,\,45^\circ ,\,\,90^\circ $, полученная путем обработки кривых на рис. 10–13 , представлена на рис. 14. Видно, что максимальный сектор сканирования в диагональной плоскости достигает величины $ \pm 50^\circ $ при $D > 20$. В плоскости вектора Н он равен $ \pm 41^\circ $ при том же ограничении диапазона, а в плоскости вектора Е сектор сканирования находится в пределах $ \pm 27^\circ $.

Рис. 14.

Зависимость рабочего диапазона решетки от угла сканирования для φ = 0° (1), 45° (2) и 90° (3).

ЗАКЛЮЧЕНИЕ

Представленные выше результаты позволяют сделать вывод о том, что ДС АР рассмотренного типа является перспективным для создания СШП решеток техническим решением. Полученные характеристики сканирования представляют практический интерес, хотя достигнутые значения диапазона и сектора сканирования не являются предельными, поскольку исследованная структура не была детально оптимизирована по совокупности показателей качества. Их следует рассматривать только как первое приближение в решении задачи построения широкоугольных СШП АР с диапазоном частот большим 10 : 1. Оптимизация такого сложного объекта, как СШП АР, и достижение предельных характеристик выходит далеко за рамки научного исследования.

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

  1. Schaubert D., Kasturi S., Elsallal M.W., Van Cappellen W. // Proc. 2006 1 st Europ. Conf. Antennas and Propagation (EuCAP). EuCAP. Nice. 06–10 Nov. N.Y.: IEEE, 2006. P. 1.

  2. Воскресенский Д.И., Котов Ю.В., Овчинникова Е.В. // Антенны. 2005. № 11. С. 7.

  3. Банков С.Е., Калошин В.А., Нгуен К.З. // Труды 4-й Всерос. микроволновой конф. Москва. Ноябрь 2016. С. 410.

  4. Acedo E., Garcia E., González-Posadas V. // IEEE Trans. V. AP-58. № 1. P. 68.

  5. Бaнкoв C.E. // PЭ. 2018. T. 63. № 6. C. 524.

  6. Банков С.Е., Скородумова Е.А. // РЭ. 2015. Т. 60. № 5. С. 470.

  7. Банков С.Е. // Журн. радиоэлектроники. 2017. № 11. http://jre.cplire.ru/jre/nov17/12/text.pdf.

  8. Банков С.Е., Дупленкова М.Д. // РЭ. 2018. Т. 63. № 1. С. 25.

  9. Банков С.Е., Дупленкова М.Д. // РЭ. 2020. Т. 65. № 3. С. 211.

  10. Банков С.Е., Курушин А.А., Гутцайт Э.М. Решение оптических и СВЧ задач с помощью HFSS. М.: Оркада, 2012.

  11. Марков Г.Т., Чаплин А.Ф. Возбуждение электромагнитных волн. М.: Радио и связь, 1983.

  12. Сазонов Д.М. Антенны и устройства СВЧ. М.: Высш. школа, 1988.

  13. Бaнкoв C.E. // PЭ. 2019. T. 64. № 8. C. 741.

  14. Разевиг В.Д., Потапов Ю.В., Курушин А.А. Проектирование СВЧ устройств с помощью MicowaveOffice. М.: Солон-Пресс, 2003.

  15. Ганстон М.А.Р. Справочник по волновым сопротивлениям фидерных линий. М.: Связь, 1976.

  16. Банков С.Е., Калошин В.А., Ле Н.Т. // РЭ. 2018. Т. 63. № 12. С. 1263.

Дополнительные материалы отсутствуют.