Радиотехника и электроника, 2020, T. 65, № 2, стр. 145-157

Сверхкороткоимпульсные излучатели на основе длинной линии

К. В. Горбачев a, Ю. И. Исаенков a, А. В. Ключник a, В. И. Мижирицкий a, В. М. Михайлов a, Е. В. Нестеров a*, В. А. Строганов a

a Объединенный институт высоких температур Российской академии наук
125412 Москва, ул. Ижорская, 13, стр. 2, Российская Федерация

* E-mail: nst@ihed.ras.ru

Поступила в редакцию 26.11.2018
После доработки 04.03.2019
Принята к публикации 10.03.2019

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

Аннотация

Представлены результаты расчетно-теоретических и экспериментальных исследований сверхкороткоимпульсных излучателей с излучающими элементами на основе вибратора и дискоконусной антенны, источником возбуждения которых являются колебания токов в длинной линии, коммутируемой ключом с малым временем коммутации τ. Определены временные характеристики поля излучения и его спектр. Показано, что в зависимости от параметра l/(cτ) (l – длина плеча вибратора) реализуются различные режимы излучения, отличающиеся частотными характеристиками и формой диаграммы направленности. Приведено качественное объяснение основных закономерностей излучения на основе расчета спектра собственных колебаний (в спектральной области) и на основе суммирования вкладов от элементарных токов, текущих по поверхности излучателя (во временной области).

ВВЕДЕНИЕ

В настоящее время широкое распространение получили излучающие системы, формирующие сверхкороткие и сверхширокополосные импульсы излучения (далее – СКИ-излучатели и СШП-излучатели). Они используются для локации и для исследований электромагнитной совместимости радиоэлектронной аппаратуры (РЭА). Подробный обзор СШП-излучателей и их применений представлен в [1]. Под сверхкороткими в данной работе будем понимать импульсы наносекундной и субнаносекундной длительности, что соответствует нескольким полупериодам колебаний электрического поля.

Классификация типов излучателей основана на их спектре, который в свою очередь определяется формой импульсов, формируемых источником напряжения. В “классических” СШП-излучателях источник питания формирует короткие однополярные видеоимпульсы, спектр излучения которых начинается от нуля, а граничная частота спектра определяется временем переключения используемых ключевых устройств (разрядников, полупроводниковых ключей, обострителей). Для формирования биполярных импульсов используются отрезки длинных линий и несколько ключей [1, 2]. В схемах с радиальной формирующей линией возможно формирование как однополярных, так и биполярных импульсов напряжения на нагрузке [3].

В данной работе рассмотрен другой тип импульсных источников, основанных на коммутации длинной линии [4, 5], спектр излучения которых имеет несколько выраженных резонансов, гармоник основной частоты. Такие излучатели имеют ряд характерных особенностей, которые ранее детально не рассматривались.

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

1. ПРИНЦИП РАБОТЫ И СХЕМА СКИ-ИЗЛУЧАТЕЛЯ

Излучатели на основе коммутируемой формирующей линии впервые подробно были рассмотрены в работах [6, 7]. В качестве источника энергии в таких излучателях используется заряженная до напряжения U0 емкость отрезка длинной линии, или емкость системы колебательных контуров. Нагрузкой является излучающий элемент, например вибратор. Таким образом, СКИ-излучатель – это устройство, в состав которого входит как источник сверхкоротких высоковольтных импульсов, так и излучающий элемент (антенна). Особенностью рассматриваемых излучателей является то, что их нельзя разделить на независимые части – генератор и антенну – и конструкция антенны существенным образом влияет на параметры формируемых импульсов.

Принципиальная схема СКИ-излучателя с питанием от длинной линии представлена на рис. 1. Отрезок длинной линии (формирующая линия) длиной L имеет волновое сопротивление W, излучатель имеет входной импеданс Z, а ключ имеет небольшое активное сопротивление r. Зарядное напряжение подается через нагрузку, импеданс которой значительно больше волнового сопротивления линии: Zн$ \gg $ W.

Рис. 1.

Принципиальная схема СКИ-излучателя.

Работа СКИ-излучателя аналогична работе формирующей линии в рассогласованном режиме. После срабатывания ключа по линии бежит волна, которая преломляется и отражается в сечении АА'. Преломленная волна уходит на вибратор, формируя на нем распределение токов и поле излучения в пространстве. Отраженная в сечении АА' волна напряжения, распространяясь обратно по линии, отражается со сменой полярности от ключа с небольшим сопротивлением, и далее процесс повторяется. Таким образом, в длинной линии и на вибраторе возникают затухающие колебания, при этом часть мощности рассеивается на активных сопротивлениях линии и ключа, а часть расходуется на излучение (рассеянием мощности на импедансе ZL будем пренебрегать). В пространстве формируется “цуг” излучения, имеющий вид затухающих по амплитуде пиков (полуволн) различной полярности. Форма цуга определяется процессами отражений волны токов и напряжений в линии и в общем случае имеет достаточно сложный вид. При определенных условиях цуг может иметь регулярную форму, близкую к затухающей синусоиде.

Определяющую роль в генерации излучения играют параметры ключа, а именно его коммутационная характеристика, задающая длительность фронта тока первого импульса и производную от тока dI/dt, и его резистивная характеристика в замкнутом состоянии, влияющая на эффективность излучателя и длительность “цуга”.

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

2. СПЕКТР СОБСТВЕННЫХ КОЛЕБАНИЙ

Исследование физических процессов в СКИ-излучателе начнем с анализа коэффициента передачи и спектра собственных колебаний.

Можно показать, что возникающие при коммутации ключа процессы определяются переходной характеристикой цепи (реакции цепи на ступеньку напряжения на входе системы). Мы определим коэффициент передачи K(ω) из которого могут быть получены переходная и импульсная характеристики. Частотная зависимость коэффициента передачи определяется в первую очередь спектром собственных колебаний, соответствующих полюсам K(ω).

Для исследования спектра собственных колебаний ограничимся случаем, когда нагрузкой длинной линии является симметричный вибратор (СВ). Будем считать, что после срабатывания ключ имеет небольшое, чисто активное сопротивление r. Входной импеданс вибратора с хорошей точностью может быть представлен импедансом эквивалентной длинной линии: Z(ω) = RΣ(ω) + iX(ω) [5]. Полуволновому резонансу соответствует частота ωv = 2πfv . Вблизи полуволнового резонанса мнимая часть импеданса обращается в ноль. Из граничных условий для токов и напряжений в длинной линии несложно найти систему уравнений для определения падающей и отраженной волн, а условие равенства нулю детерминанта этой системы дает частоты собственных колебаний системы:

(1)
$\begin{gathered} \det = \left( {1 + \frac{r}{W}} \right)\left( {1 + \frac{Z}{W}} \right)\exp \left( {ikL} \right) - \\ - \,\,\left( {1 - \frac{r}{W}} \right)\left( {1 - \frac{Z}{W}} \right)\exp \left( { - ikL} \right) = 0. \\ \end{gathered} $

Если на одном конце линии имеем режим, близкий к короткому замыканию, а на другом – к холостому ходу, уравнение (1) можно решить разложением по малым параметрам r0/W $ \ll $ 1 и W/|Z| $ \ll $ 1. Решение уравнения (1) ищем в виде

(2)
$\begin{gathered} f = {{f}_{n}} + \delta {{f}_{n}},\,\,\,\,{{f}_{n}} = {{{{f}_{0}}\left( {2n + 1} \right)} \mathord{\left/ {\vphantom {{{{f}_{0}}\left( {2n + 1} \right)} 4}} \right. \kern-0em} 4}, \\ {{f}_{0}} = {c \mathord{\left/ {\vphantom {c L}} \right. \kern-0em} L},\,\,\,\,n = 0,1,... \\ \end{gathered} $

Подставляя (2) в (1) и разлагая по малым параметрам до членов второго порядка по δfn, получаем из (1) квадратное уравнение вида

(3)
$A{{x}^{2}} + x + iB = 0,$

где x = δω/f0, A = Z  '(fn)f0/Z(fn), B = W/Z(fn) + r/W, Z '′(ω) = ∂Z(ω)/∂ω . Два корня этого уравнения дают поправки δfn к собственным частотам.

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

(4)
$\begin{gathered} \delta {{f}_{{n1}}} \approx - iB = \left( {\frac{2}{\pi }} \right) \times \\ \times \,\,{{f}_{0}}\left\{ {i\left( {\frac{r}{W} + \frac{{W{{R}_{\sum }}\left( {{{\omega }_{n}}} \right)}}{{{{{\left| {Z\left( {{{\omega }_{n}}} \right)} \right|}}^{2}}}}} \right) + \frac{{WX\left( {{{\omega }_{n}}} \right)}}{{{{{\left| {Z\left( {{{\omega }_{n}}} \right)} \right|}}^{2}}}}} \right\}. \\ \end{gathered} $

Квадратичный член в (2) имеет смысл удерживать лишь при условии, что коэффициент A $ \gg $ 1. Это может иметь место вблизи полуволнового резонанса, и именно в этой области частот второй корень уравнения (3) также оказывается вблизи действительной оси, что соответствует появлению второго выраженного резонанса. При A $ \gg $ 1 два корня уравнения (3) имеют вид

(5)
$\delta {{f}_{{n1}}} \approx - iB,\,\,\,\,\delta {{f}_{{n2}}} \approx - {{{{f}_{0}}} \mathord{\left/ {\vphantom {{{{f}_{0}}} {\left( {2\pi A} \right)}}} \right. \kern-0em} {\left( {2\pi A} \right)}}.$

Полученные поправки к частоте дают затухание собственных колебаний (мнимая часть δfn) и сдвиг резонансной частоты относительно своего невозмущенного значения (действительная часть δfn).

Проанализируем сначала собственные волны, соответствующие первой моде δfn1. Как следует из (4), затухание этой моды колебаний складывается из двух частей – диссипативного затухания на сопротивлении ключа и потерь на излучение. Потери на излучение сравниваются с потерями на ключе при

(6)
$W = Wc = {{\left( {{r \mathord{\left/ {\vphantom {r {{{R}_{\sum }}}}} \right. \kern-0em} {{{R}_{\sum }}}}} \right)}^{{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}}}}\left| Z \right|.$

Оценим порядок величины Wc вблизи полуволнового резонанса при сопротивлении ключа 0.3 Ом: Wc ~ (0.3 × 70)1/2 = 4.5 Ом. При W > Wc превалируют потери на излучение, а при W < Wc – потери в ключе.

На низких частотах (по сравнению с полуволновым резонансом), при ω → 0, имеем RΣ → 0, |Z| → ∞ (импеданс вибратора близок к импедансу емкости). В результате Wc → ∞ и потери на излучение оказываются меньше потерь на сопротивлении ключа. Практически вся энергия в этом случае рассеивается на ключе.

Энергетическая эффективность рассматриваемой системы определяется отношением потерь на излучение к омическим потерям. При волновом сопротивлении линии, значительно меньшем сопротивления вибратора, эффективность растет пропорционально W 2. При увеличении волнового сопротивления линии КПД преобразования запасенной в системе энергии в энергию излучения стремится к единице, а абсолютная величина излучаемой мощности имеет максимум вблизи Wc.

Смещение резонанса относительно fn минимально на частотах, где реактивная составляющая импеданса вибратора равна нулю. При этом даже на частотах, где Xn) = 0, смещение резонанса отлично от нуля и определяется поправками более высокого порядка.

При приближении резонанса линии к полуволновому резонансу вибратора возникает вторая мода колебаний и спектральное распределение принимает характерный “двухпичковый” вид. Как следует из (5), добротность второй моды пропорциональна (RΣX ' – $XR_{\Sigma }^{'}$)/|Z|2. Характерной особенностью данной моды является возможность резкого увеличения добротности при изменении частоты. Сдвиг частоты пропорционален ${{ - ({{R}_{\Sigma }}R_{\Sigma }^{'} + XX{\kern 1pt} ')} \mathord{\left/ {\vphantom {{ - ({{R}_{\Sigma }}R_{\Sigma }^{'} + XX{\kern 1pt} ')} {{{{\left| Z \right|}}^{2}}}}} \right. \kern-0em} {{{{\left| Z \right|}}^{2}}}}.$ В окрестности полуволнового резонанса производные импеданса вибратора положительны, поэтому частота будет сдвигаться в область более низких частот по сравнению с резонансом длинной линии.

Коэффициент передачи для рассматриваемой системы находили из S-матрицы четырехполюсника, состоящего из трех последовательно соединенных четырехполюсников – ключа, линии и вибратора. Численные расчеты частотной зависимости коэффициента передачи полностью соответствуют полученным выше выражениям в соответствующих предельных случаях. Как будет показано ниже, эффективность СКИ-излучателей зависит от взаимного положения резонансов длинной линии и резонансов излучателя. Мы не будем приводить детальный анализ K(ω), так как ниже будут использоваться точные численные решения для токов и напряжений. В качестве примера на рис. 2 приведен результат численного расчета спектра тока в центре длинной линии с волновым сопротивлением W = 30 Ом длиной 15 см, подключенной к симметричному цилиндрическому вибратору с длиной плеча 16 см и диаметром 3 см. В расчете линия заряжалась до единичного напряжения 1 В и коммутировалась через резистор r = 0.3 Ом. Продолжительность (время) коммутации ключа τ (изменение импеданса от бесконечности до r) выбирали равной 0.15 нс. Далее модель вибратора и длинной линии с такими параметрами назовем СВ1.

Рис. 2.

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

В спектре импульса тока видны резонансы длинной линии на частотах 480, 1370 и 2330 МГц. При этом расчетное значение f0 ≈ 500 МГц. Полуволновой резонанс этого вибратора близок к 500 МГц, поэтому низкочастотный резонанс имеет два пика, причем частота низкочастотного пика сдвинута относительно f0 примерно на 180 МГц в область низких частот. Низкочастотная часть спектра тока определяется длительностью цуга излучения T. Заметная доля энергии содержится на частотах вплоть до f ~ 1/T.

3. ХАРАКТЕРИСТИКИ СКИ-ИЗЛУЧАТЕЛЕЙ

Прежде чем переходить к анализу полей, определим основные, наиболее часто используемые характеристики при описании работы СШП- и СКИ-излучателей [1, 4].

Анализ полей будем выполнять в дальней зоне. Для гармонических сигналов граница дальней зоны определяется условиями R > 2D2/λ при D > λ или R > (2…3)λ при D < λ. Для СКИ-сигналов, которые содержат в спектре как низкие, так и высокие частоты, определение границы дальней зоны требует уточнения. Граница дальней зоны СКИ-излучателя определяется по критерию неизменности формы импульса излучения в фиксированном направлении или по критерию синфазности Е- и Н-компонент поля. Оба эти критерия эквивалентны. Для инженерной оценки границ дальней зоны можно использовать условия R ≥ 2D2/λ и R > 4D, где D – размер излучателя, а в качестве λ используется верхняя частота спектра, при этом погрешности, связанные с пренебрежением высшими члена разложения на границе дальней зоны, составят не более 10% (подробнее см. в [1]).

В дальней зоне векторный потенциал $\vec {A}$(R, t) равен интегралу от плотности тока на вибраторе $\vec {J}\left( {r,t} \right)$:

$\vec {A} = {{{{\mu }_{{\text{0}}}}} \mathord{\left/ {\vphantom {{{{\mu }_{{\text{0}}}}} {\left( {4\pi R} \right)}}} \right. \kern-0em} {\left( {4\pi R} \right)}}\int {\vec {J}\left( {r,t - {R \mathord{\left/ {\vphantom {R {c - {{\left( {\vec {n}\vec {r}} \right)} \mathord{\left/ {\vphantom {{\left( {\vec {n}\vec {r}} \right)} c}} \right. \kern-0em} c}}}} \right. \kern-0em} {c - {{\left( {\vec {n}\vec {r}} \right)} \mathord{\left/ {\vphantom {{\left( {\vec {n}\vec {r}} \right)} c}} \right. \kern-0em} c}}}} \right)} dr,$

а напряженность магнитного и электрического полей

$\begin{gathered} \vec {H}\left( {\vec {R},t} \right) = {{{\text{rot}}\vec {A}} \mathord{\left/ {\vphantom {{{\text{rot}}\vec {A}} {{{\mu }_{0}}}}} \right. \kern-0em} {{{\mu }_{0}}}} = {{\left[ {{{{\vec {A}}}_{t}}\vec {n}} \right]} \mathord{\left/ {\vphantom {{\left[ {{{{\vec {A}}}_{t}}\vec {n}} \right]} {c{{\mu }_{0}},}}} \right. \kern-0em} {c{{\mu }_{0}},}} \\ \vec {E}\left( {\vec {R},t} \right) = \left[ {\vec {n}\left[ {{{{\vec {A}}}_{t}}\vec {n}} \right]} \right]. \\ \end{gathered} $

Здесь $\vec {n} = {{\vec {R}} \mathord{\left/ {\vphantom {{\vec {R}} R}} \right. \kern-0em} R},\vec {R}$ – радиус-вектор направленный в точку наблюдения, ${{\vec {A}}_{t}} = {{\partial A} \mathord{\left/ {\vphantom {{\partial A} {\partial t}}} \right. \kern-0em} {\partial t}}$. Форма импульса поля в дальней зоне зависит от направления $\vec {n}$ и не зависит от расстояния R, которое входит в виде комбинации t* = t – R/c и дает только сдвиг начала отсчета времени. Электрическое и магнитное поля в дальней зоне поперечны и синфазны.

Построим в дальней зоне сферу радиусом R с центром в СКИ-излучателе, $\vec {n}$ – единичный вектор направленный из центра сферы в точку наблюдения. Электрическое поле в каждой точке на сфере в направлении $\vec {n}$ можно представить в виде

(7)
$\vec {E}\left( {\vec {n},t} \right) = {{\vec {e}}_{\theta }}{{E}_{\theta }}\left( {\vec {n}} \right){{u}_{\theta }}\left( {\vec {n},t{\text{*}}} \right) + {{\vec {e}}_{\varphi }}{{E}_{\varphi }}\left( {\vec {n}} \right){{u}_{\varphi }}\left( {\vec {n},t{\text{*}}} \right).$

Здесь ui($\vec {n}$, t) – скалярные функции единичной амплитуды, описывающие временную зависимость поля, а ${{E}_{j}}\left( {\vec {n}} \right)$ – пиковые значения компонент поля в направлении $\vec {n}$, θ и φ углы сферической системы координат. Как следует из (7), пространственное распределение излучения (в том числе поляризация) изменяется во времени. В каждый фиксированный момент времени, нормируя каждую компоненту в (7) на ее максимальное значение, можно построить пространственную диаграмму направленности (ДН). Наибольший интерес представляет распределение поля в моменты времени t1t2…, соответствующие максимумам излучения. В общем случае эти моменты времени в разных направлениях разные. Однако для рассматриваемых СКИ-излучателей разброс времен достижения максимумов в разных направлениях значительно меньше длительности полуволны (см. ниже), поэтому для анализа направленных свойств им можно пренебречь. Для анализа пространственного распределения излучения во временной области будем использовать нормированные ДН в моменты времени, соответствующие максимумам излучения:

${{F}_{\theta }}\left( {\vec {n},{{t}_{1}}} \right) = {{{{E}_{\theta }}\left( {\vec {n},{{t}_{1}}} \right)} \mathord{\left/ {\vphantom {{{{E}_{\theta }}\left( {\vec {n},{{t}_{1}}} \right)} {{{E}_{\theta }}\left( {{{{\vec {n}}}_{1}},{{t}_{1}}} \right)}}} \right. \kern-0em} {{{E}_{\theta }}\left( {{{{\vec {n}}}_{1}},{{t}_{1}}} \right)}},$
${{F}_{\theta }}\left( {\vec {n},{{t}_{2}}} \right) = {{{{E}_{\theta }}\left( {\vec {n},{{t}_{2}}} \right)} \mathord{\left/ {\vphantom {{{{E}_{\theta }}\left( {\vec {n},{{t}_{2}}} \right)} {{{E}_{\theta }}\left( {{{{\vec {n}}}_{2}},{{t}_{2}}} \right)}}} \right. \kern-0em} {{{E}_{\theta }}\left( {{{{\vec {n}}}_{2}},{{t}_{2}}} \right)}}, \ldots ,$

где ${{\vec {n}}_{1}}$,$\vec {n}{}_{2}$ – направления максимального излучения. Аналогичным образом могут быть получены ДН для φ-компоненты поля и для полного поля.

Введенные таким образом ДН можно связать с традиционными ДН для гармонических сигналов $\vec {F}\left( {\omega ,\vec {n}} \right)$. Фурье-компонента электрического поля на частоте ω имеет вид

$\vec {E}\left( {\omega ,\vec {n}} \right) = c\left( \omega \right){{\exp \left( { - ikR} \right)} \mathord{\left/ {\vphantom {{\exp \left( { - ikR} \right)} {R \cdot \vec {F}\left( {\omega ,\vec {n}} \right)}}} \right. \kern-0em} {R \cdot \vec {F}\left( {\omega ,\vec {n}} \right)}},$

где c(ω) = Emax(ω)R, Emax(ω) – величина поля в направлении максимума диаграммы на частоте ω . Для θ-компоненты поля получаем

(8)
$\begin{gathered} {{F}_{\theta }}\left( {\vec {n},{{t}_{1}}} \right) = \\ = \int {\left[ {{{{{E}_{{\theta \max }}}\left( \omega \right)} \mathord{\left/ {\vphantom {{{{E}_{{\theta \max }}}\left( \omega \right)} {{{E}_{\theta }}\left( {{{{\vec {n}}}_{1}},{{t}_{1}}} \right)}}} \right. \kern-0em} {{{E}_{\theta }}\left( {{{{\vec {n}}}_{1}},{{t}_{1}}} \right)}}} \right]} {{F}_{\theta }}\left( {\omega ,\vec {n}} \right)\exp \left( {i\omega {{t}_{1}}} \right)d\omega \,. \\ \end{gathered} $

Таким образом, ДН во временной области ${{F}_{\theta }}\left( {\vec {n},{{t}_{1}}} \right)$ есть сумма диаграмм направленности на разных частотах, взятых с соответствующим весом. Из (8) следует, что если излучатель имеет частотно независимую диаграмму направленности во всем значимом частотном диапазоне СКИ-излучения, то все его диаграммы направленности во временной области по модулю совпадают с этой диаграммой направленности и отличаются только фазовыми множителями.

Излучаемая в направлении $\vec {n}$ энергия импульса имеет вид

$\begin{gathered} {{W}_{{{\text{изл}}}}}\left( {\vec {n}} \right) = {{\left( {{{\tau }_{\theta }}{\kern 1pt} {{E}_{\theta }}{{{\left( {\vec {n}} \right)}}^{2}} + {{\tau }_{\varphi }}{\kern 1pt} {{E}_{\varphi }}{{{\left( {\vec {n}} \right)}}^{2}}} \right)} \mathord{\left/ {\vphantom {{\left( {{{\tau }_{\theta }}{\kern 1pt} {{E}_{\theta }}{{{\left( {\vec {n}} \right)}}^{2}} + {{\tau }_{\varphi }}{\kern 1pt} {{E}_{\varphi }}{{{\left( {\vec {n}} \right)}}^{2}}} \right)} {{{Z}_{0}},}}} \right. \kern-0em} {{{Z}_{0}},}} \\ {{\tau }_{\theta }} = {{\tau }_{\theta }}\left( {\vec {n}} \right) = \int {{{u}_{\theta }}{{{\left( {\vec {n},t} \right)}}^{2}}dt} , \\ \end{gathered} $

где Z0 = 120π – волновое сопротивление свободного пространства. Нормируя эту величину на максимальное значение, получаем вполне однозначную диаграмму направленности по энергии ${{F}^{2}}\left( {\vec {n}} \right)$.

Кроме ДН для анализа СКИ-излучателей используют эффективный потенциал (ЭП) и энергетический КПД (эффективность СКИ-излучателя). Эффективный потенциал – произведение максимального пикового значения напряженности поля в дальней зоне на расстояние R:

(9)
${\text{ЭП}} = {{E}_{{\max }}}R.$

В дальней зоне эта величина не зависит от расстояния. Иногда используют также значения эффективного потенциала в разных направлениях. Так как эффективный потенциал прямо пропорционален зарядному напряжению U0, то удобно использовать безразмерный нормированный эффективный потенциал:

(10)
${{{{E}_{{\max }}}R} \mathord{\left/ {\vphantom {{{{E}_{{\max }}}R} {{{U}_{0}}}}} \right. \kern-0em} {{{U}_{0}}}}.$

Энергетический КПД равен отношению излученной энергии к энергии, запасенной перед началом коммутации.

СКИ-излучатели оптимизируются по указанным выше характеристикам за счет выбора параметров длинной линии и параметров излучателя.

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

4. АНАЛИЗ ПОЛЕЙ ИЗЛУЧЕНИЯ

Для расчета полей излучения основной характеристикой является распределение тока по излучающему элементу. “Источником” в данном случае является волна токов, возникающая при коммутации заряженной длинной линии.

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

Для вибраторов, которые практически могут быть использованы в СКИ-излучателях, распределение тока может быть получено только численными методами. В представленных ниже результатах численных расчетов длинная линия заряжалась до напряжения 1 В, а затем срабатывал ключ, формировался переходной колебательный процесс и система разряжалась.

На рис. 3а представлены осциллограммы токов в длинной линии и на поверхности конуса в сечении aa' для излучающего элемента в виде дискоконусной антенны с параметрами: длина b = 30 см, диаметр a = 23 см, длина линии L = 20 см, W = 30 Ом, r = 0.3 Ом, время коммутации τ = 0.15 нс, зарядное напряжение – U0 = 1 В (далее – модель ДК1). Коммутация ключа осуществляется в момент времени t = 0. Сечение дискоконусной антенны вертикальной плоскостью и обозначения размеров показаны на рисунке 3б. Внутренний проводник коаксиальной линии соединен в верхней части с диском. Коммутация внутреннего проводника заряженной коаксиальной линии на внешний проводник (конус) осуществляется ключом К, расположенным в нижней части линии (ключ показан условно).

Рис. 3.

Токи (а) в линии (сплошная кривая) и на поверхности конуса (штриховая) дискоконусной антенны в сечении aa'; дискоконусная антенна, К – ключ (б).

На рис. 4а показаны сечения симметричного цилиндрического вибратора, в которых рассчитывались токи, представленные на рис. 4б. На рис. 4в в растянутом масштабе представлен первый полупериод токов. Сечения 1 и 11, а также 2 и 22 расположены симметрично относительно центра вибратора. Для расчетов использовалась модель СВ1 (см. разд. 2).

Рис. 4.

Сечения вибратора, выбранные для расчета токов (а), токи в сечениях 1, 2 и 3 (б), первая полуволна токов в сечениях 1, 2, 3, 11, 22 (в).

В осциллограммах тока можно выделить два временных интервала – область “вынужденных” колебаний (две первые полуволны) и область свободных колебаний (начиная с третьей полуволны). В области вынужденных колебаний токи в разных сечениях вибратора расфазированы, а первая полуволна тока имеет короткий фронт, примерно равный времени коммутации ключа τ. Сдвиги токов первой полуволны в разных сечениях определяются временем распространения импульса тока по вибратору от одного сечения до другого со скоростью примерно равной скорости света. Начиная с третьей полуволны расфазировка токов по вибратору незначительна.

Наличие первого максимума с крутым фронтом является характерной особенностью рассматриваемых СКИ-излучателей. Так как поле излучения пропорционально производной от тока, то наиболее крутой части фронта импульса тока будет соответствовать первый максимум поля излучения. Второй максимум излучения (он может даже превосходить первый) обусловлен значительным перепадом между первой и второй полуволнами тока.

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

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

Указанные выше особенности распределения тока по излучателю приводят к зависимости формы импульса от направления излучения. Зависимость формы импульса излучения от направления была исследована для СШП-импульсов (см., например, [1]). В спектральной области этот эффект можно объяснить зависимостью формы ДН вибратора от частоты [8]. В частности, при l/λ < 0.5 ДН вибратора дипольного типа, а на более высоких частотах возникают боковые лепестки и максимумы “прижимаются” к оси вибратора.

В данной работе рассмотрен эффект зависимости формы импульса от направления распространения применительно к СКИ-излучателям и приведена его простая физическая интерпретация во временной области.

Для анализа поля излучения во временной области разобьем вибратор на элементы, в пределах которых ток (“элементарный” ток) можно считать не зависящим от координаты (рис. 5). В качестве шага разбиения возьмем длину Δz < cτ.

Рис. 5.

Суммирование вкладов элементарных токов в полное излучение под разными углами к оси вибратора.

Векторный потенциал равен сумме вкладов от элементарных токов:

(11)
$\begin{gathered} \vec {A} = {{{{\mu }_{0}}} \mathord{\left/ {\vphantom {{{{\mu }_{0}}} {\left( {4\pi R} \right)}}} \right. \kern-0em} {\left( {4\pi R} \right)}}\sum {{{{\vec {J}}}_{m}}} \left( {t - {R \mathord{\left/ {\vphantom {R c}} \right. \kern-0em} c} - {{\left( {\vec {n}{{{\vec {r}}}_{m}}} \right)} \mathord{\left/ {\vphantom {{\left( {\vec {n}{{{\vec {r}}}_{m}}} \right)} c}} \right. \kern-0em} c}} \right), \\ m = 1,2,...,{l \mathord{\left/ {\vphantom {l {\Delta z}}} \right. \kern-0em} {\Delta z}}. \\ \end{gathered} $

Как следует из численных расчетов, в пределах первой полуволны форма тока в разных сечениях для основной части вибратора меняется слабо, убывая по амплитуде к его краям. При этом временная задержка между соседними элементарными токами примерно равна времени распространения сигнала Δz/c. В этом приближении для оценки суммы можно произвести замену

${{J}_{m}}\left( {t{\text{*}} - \left( {\vec {n}{{{\vec {r}}}_{m}}} \right)} \right) \to {{C}_{m}}\vec {J}\left( {t{\text{*}} + {{m\Delta z} \mathord{\left/ {\vphantom {{m\Delta z} c}} \right. \kern-0em} c} - m\Delta z{{\cos \left( \theta \right)} \mathord{\left/ {\vphantom {{\cos \left( \theta \right)} c}} \right. \kern-0em} c}} \right),$

а в качестве $\vec {J}\left( t \right)$ использовать ток в центральном сечении плеча вибратора. Вклад в поле излучения от элементарного тока вычисляется по формуле

(12)
$\begin{gathered} \delta {{{\vec {E}}}_{m}} = {{{\vec {e}}}_{\theta }}{{\sin \left( \theta \right)} \mathord{\left/ {\vphantom {{\sin \left( \theta \right)} {\left( {4\pi R} \right)}}} \right. \kern-0em} {\left( {4\pi R} \right)}} \times \\ \times \,\,{{C}_{m}}{{J}_{t}}\left( {t{\text{*}} + {{m\Delta z} \mathord{\left/ {\vphantom {{m\Delta z} c}} \right. \kern-0em} c}\left( {1 - \cos \left( \theta \right)} \right)} \right). \\ \end{gathered} $

Производная от тока Jt(t) локализована на фронте и спаде первой полуволны. На фронте она имеет вид пика амплитудой ~U0/(Wτ) и длительностью ~τ (см. рис. 5).

Полное поле излучения является сумой вкладов элементарных токов, причем результирующее поле определяется временными задержками на распространение вида mΔzcos(θ)/c и временными задержками в распределении тока по вибратору ~mΔz/c.

Если выполнить суммирование полей излучения элементарных токов $\delta {{\vec {E}}_{m}}$, то в зависимости от параметра l/(сτ) получим совершенно разные пространственные распределения поля.

При l/(сτ) $ \ll $ 1 максимум излучения будет находиться в плоскости (xy) центрального сечения вибратора, где пространственные задержки обращаются в ноль.

При l/(сτ) $ \gg $ 1, максимум излучения смещается к оси на угол θ ~ 30°…50°. Чем меньше τ, тем ближе первый максимум прижимается к оси вибратора. Длительность импульса от первой полуволны в плоскости центрального сечения оказывается порядка l/c, а длительность импульса под небольшим углом к оси оказывается меньше, и для очень малых времен переключения она будет стремиться к τ. Такое распределение поля объясняется зависимостью задержек от направления распространения. Действительно, основной вклад в поле в плоскости центрального сечения (плоскость xy, перпендикулярная оси вибратора и проходящая через его центр) от первой полуволны тока вносит центральная часть вибратора от +сτ до –сτ (область “синхронности” элементарных токов). Поля от других частей вибратора приходят с временной задержкой, большей длительности импульса элементарного тока и не увеличивают амплитуду, а формируют “плоскую” часть импульса излучения. При l/(сτ) $ \gg $ 1 размеры области “синхронности” элементарных токов оказываются меньше длины плеча вибратора. Поэтому при увеличении длины плеча вибратора амплитуда максимума поля изменяться не будет.

При уменьшении угла θ временные задержки на распространение импульса от соседних элементарных токов компенсируют в определенной степени задержки токов на поверхности вибратора. При этом размеры области синхронности элементарных токов возрастают и могут приблизиться к длине плеча вибратора. Так как элементарные токи вдоль оси не излучают ($\delta {{\vec {E}}_{m}}\sim \sin \left( \theta \right)$), то максимум поля будет расположен под некоторым углом к оси вибратора. Положение максимума будет смещаться ближе к оси вибратора при уменьшении τ, а величина этого максимума будет расти пропорционально длине плеча вибратора. Так как распределения токов на разных плечах вибратора почти совпадают (см. рис. 4в), то будут формироваться симметричные максимумы в противоположных направлениях.

Амплитуда поля в центральном сечении (xy), которая пропорциональна размеру области синхронности элементарных токов (~τ) и производной от элементарного тока (~1/τ), практически не будет зависеть от длительности фронта. Амплитуда поля под углом к центральному сечению будет расти с уменьшением длительности фронта, так как амплитуда элементарного тока пропорциональна 1/τ, а размеры области синхронности от длительности фронта почти не зависят. При некотором значении времени коммутации τ = τс амплитуды максимумов в разных направлениях сравниваются. Величина τс зависит от конструкции антенны.

Указанные качественные соображения не учитывают ограниченную полосу рабочих частот цилиндрического вибратора, которая начинает сказываться при небольших временах коммутации τ. Для более широкополосных излучателей увеличение амплитуды с уменьшением τ будет более значительным (см. ниже).

На рис. 6 представлены результаты численных расчетов зависимости эффективного нормированного потенциала ER/U0 от продолжительности коммутации τ в плоскости центрального сечения (2) и под углом 45° к оси симметричного вибратора.

Рис. 6.

Зависимость нормированного эффективного потенциала (первая полуволна) от времени коммутации τ под углами θ = 45° (1) и 90° (2).

Результаты расчетов подтверждают представленные выше качественные соображения: зависимость амплитуды максимума при θ = 90° оказывается значительно более плавной, чем при θ = 45°.

В “поздневременной” области цуга колебания токов на всей поверхности вибратора практически синфазны. Временные задержки обращаются в ноль в центральном сечении вибратора. Поэтому ДН для рассмотренных на рис. 3 и 4 вибраторов для этой части цуга излучения будет дипольного типа с максимумом в плоскости центрального сечения.

Для иллюстрации представленных выше качественных соображений на рис. 7 и 8 представлены результаты численных расчетов полей излучения в центральном сечении вибратора и под углом 45° к его оси, а также угловые распределения абсолютной величины электрического поля E(θ) (θ – угол от оси вибратора) на расстоянии 1 м в моменты времени, соответствующие первым двум максимумам излучения. Моделирование проводилось для модели СВ1.

Рис. 7.

Поперечная компонента электрического поля излучения (а) и модуль спектральной плотности электрического поля (б) под углами θ = 90° (сплошная) и θ = 45°(штриховая) к оси вибратора.

Рис. 8.

Распределение абсолютной величины электрического поля E(θ) в плоскости (zy): а – первый максимум излучения (t1 = 0.9 нс, Emax = 0.39 В/м, θmax = 56°), б – второй максимум излучения (t2 = 2.0 нс, Emax = 0.46 В/м, θmax = 54°).

Из рисунков видно, что длительность первой полуволны под углом 45° к оси меньше, чем в плоскости центрального сечения. Спектр излучения под углом 45° смещен в сторону высоких частот. Абсолютная величина первого максимума поля под углом 45° на 10% больше соответствующего значения в плоскости центрального сечения. Длительность первой полуволны (по 3 дБ) в плоскости центрального сечения составила 0.5 нс, а под углом 45° – 0.17 нс, при этом длительность импульса в плоскости центрального сечения примерно соответствует времени распространения волны по плечу вибратора l/c. Расчеты показали также, что пространственное распределение поля в поздневременной области – области свободных колебаний – близко к распределению дипольного типа.

Распределение поля по углу θ демонстрирует формирование максимумов излучения под углом к оси вибратора.

Для дискоконусной антенны эффект изменения формы импульса от направления может быть выражен значительно более сильно. Как уже отмечалось ранее, для симметричного цилиндрического вибратора смещение максимума к оси “подавлено” множителем sin(θ) в элементарном токе. Для конического излучателя фактор sin(θ), подавляющий излучение вдоль оси, в данном случае не работает, так как элементарные токи направлены под разными углами к оси. В результате амплитуда первого максимума под небольшим углом к оси будет значительно превышать амплитуду излучения под углом 90°.

В качестве примера на рис. 9 представлены осциллограммы поперечной компоненты электрического поля на расстоянии 1 м и зависимость абсолютной величины поля E(θ) в момент времени, соответствующий первому максимуму для модели ДК1 (см. рис. 3). На рис. 10 представлены соответствующие этим осциллограммам спектры.

Рис. 9.

Поперечная компонента электрического поля под углами θ = 90° (сплошная кривая), θ = –45° (штриховая) и θ = +45° (пунктирная) к оси конуса (а) и пространственное распределение поля E(θ) в момент первого максимума t1 = 1.5 нс, Emax = 0.76 В/м, θmax = 131° для дискоконусной антенны (б).

Рис. 10.

Зависимость модулей спектральной плотности электрического поля от частоты под углами 45° (штриховая кривая) и 90° (сплошная) к оси дискоконусной антенны.

Видно, что высокие частоты высвечиваются в основном под небольшим углом к оси, а низкие частоты сосредоточены в плоскости (xy). Низкочастотный резонанс соответствует собственным колебаниям системы.

Интересной особенностью СКИ-излучателя является перемещение максимума излучения по угловым координатам в процессе излучения. Этот эффект виден на рис. 5 и 8. В частности, для дискоконусной антенны (см. рис. 8) максимум излучения сначала находится под углом θ = 90°, а затем перемещается на угол θ = –45°. В моменты времени, соответствующие второму и последующим максимумам, максимум расположен под углом θ = 90°.

Перемещение максимума можно объяснить изменениями временных задержек элементарных токов по вибратору. Эти задержки на фронте импульса (первого полупериода) равны Δz/c, а в зоне свободных колебаний (после третьего полупериода) они практически равны нулю.

СКИ-излучатель в общем случае не имеет фазового центра. В частности, существуют интервалы времени, когда поля излучения на сферической поверхности в разных направлениях противофазны (см. рис. 9а). Отсюда, в частности, следует, что использование СКИ-излучателя в качестве облучателя требует специального расчета зеркала.

5. ОПТИМИЗАЦИЯ СКИ-ИЗЛУЧАТЕЛЕЙ

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

В локационных приложениях необходимо получить максимальную амплитуду первого пика при максимально крутых фронтах и минимальной длительности. Уменьшение длительности фронта и первой полуволны, направленной под углом 30°…45° к оси, может быть достигнуто снижением времени срабатывания ключа.

В ряде приложений требуется получить максимальный эффективный потенциал. Увеличить эффективный потенциал можно как за счет уменьшения ширины диаграммы направленности, так и за счет перераспределения основной доли (80% и более) энергии цуга в две первые полуволны. Такое перераспределение обеспечивается выбором соответствующих параметров излучателя.

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

Увеличение направленности одного СКИ-излучателя может быть получено за счет увеличения его размеров, как в продольном, так и в поперечном направлении. Отметим, что ширина ДН первого максимума оказывается заметно меньше, чем ширина ДН в плоскости поперечного сечения. Например, на рис. 8 ширина ДН по уровню 3 дБ составляет 44° для первого максимума под углом 45° к оси, в то время как в плоскости поперечного в максимуме второго полупериода ширина ДН составляет 90°. Это объясняется тем, что более высокие частоты высвечиваются под небольшими углами к оси.

6. РЕЗУЛЬТАТЫ ЭКСПЕРИМЕНТАЛЬНЫХ ИССЛЕДОВАНИЙ

Экспериментальные исследования проводились с двумя типами СКИ-излучателей:

– макет 1 (симметричный вибратор с длиной плеча 15 см, диаметром 6 см, длинная линия с волновым сопротивлением W = 60 Ом, длиной 13 см),

– макет 2 – вибратор над проводящим диском (диаметр вибратора 6 см, длина плеча 17.5 см, длина линии 15 см, волновое сопротивление W = = 60 Ом).

Зарядное напряжение для исследуемых излучателей составляло от 100 до 300 кВ.

Оценки параметров ключа, используемого в экспериментах, показывают, что он может обеспечить длительность фронта первого импульса ~200 пс и имеет сопротивление в состоянии проводимости ~100…300 мОм.

Электрическое поле измерялось в дальней зоне (на расстоянии 2.5 м и более) с помощью преобразователей напряженности импульсного электрического поля – ИППЛ-Л и ИППЛ-Д (производство ВНИИОФИ) с временем нарастания переходной характеристики 50 и 200 пс и длительностью переходной характеристики 5.6 и 25 нс соответственно. Измерения проводились под разными углами к оси излучателя. На рис. 11–12 для макета 1 представлены осциллограммы поля и их спектры, измеренные в плоскости центрального сечения, под углом θ = 90°.

Рис. 11.

Осциллограмма цуга (кривая 1) и его спектр (кривая 2) при θ = 90°.

Рис. 12.

Осциллограмма первых трех полуволн цуга при θ = 90°.

Низкочастотный резонанс в спектре излучения в плоскости поперечного сечения имеет характерную двухпичковую форму, так как длина линии примерно равна длине плеча вибратора. Первый пик имеет максимальную амплитуду, а его ширина меньше ширины второй полуволны. На рис. 12 представлены три первые полуволны цуга. Из рисунка видно, что длительность фронта примерно равна 0.2 нс, а длительность первой полуволны примерно равна l/c и в два раза меньше длительности второй полуволны.

На рис. 13 представлена осциллограмма поля и его спектр, измеренные под углом θ = 50° к оси вибратора. Спектр излучения под углом 50° смещен в более высокочастотную область, длительность первой полуволны заметно меньше, а нормированный эффективный потенциал больше чем в плоскости центрального сечения.

Рис. 13.

Осциллограмма цуга (1) и его спектр (2) при θ = 50°.

На рис. 14 представлены измеренные зависимости нормированного эффективного потенциала ER/U0 от угла θ для макетов 1 и 2. Эти диаграммы построены по первому (основному) пику излучения. Максимум излучения для рассматриваемых макетов направлен под углами 30°…50° к оси.

Рис. 14.

Зависимость ER/U0 от угла θ для макетов 1 и 2.

ЗАКЛЮЧЕНИЕ

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

В результате численного моделирования исследуемых систем получены характерные зависимости поля излучения от времени, которые имеют форму цуга, содержащего короткий первый импульс, определяемый продолжительностью коммутации ключа, и последующую поздневременную часть, соответствующую собственным колебаниям системы. Спектр излучения имеет максимумы на частотах, соответствующих резонансам длинной линии f = c/4L, 3c/4L,… (L – длина линии). При этом низкочастотная часть спектра определяется длительностью цуга излучения.

Форма диаграммы направленности симметричного цилиндрического вибратора существенно зависит от параметра l/(cτ). При l/(cτ) $ \gg $ 1 основной пик излучения имеет фронт длительностью τ и направлен под углом около θ ~ 30°…45° к оси излучателя, с уменьшением τ его амплитуда растет, а угол θ уменьшается, максимум прижимается ближе к оси. При l/(cτ) $ \ll $ 1 максимум излучения лежит в плоскости центрального сечения вибратора (θ = 90°), длительность первой полуволны порядка l/c. Более того, в разные моменты времени направление максимума излучения изменяется, происходит смещение максимума по угловым координатам. Схожий характер поведения демонстрирует форма ДН дискоконусной антенны. Качественный анализ результатов расчета спектра собственных колебаний (в спектральной области) и результатов суммирования вкладов от элементарных токов, текущих по поверхности излучателя (во временной области) позволяет объяснить основные закономерности излучения.

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

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

  1. Беличенко В.П., Буянов Ю.И., Кошелев В.И. Сверхширокополосные импульсные радиосистемы. Новосибирск: Наука, 2015.

  2. Губанов В.П., Ефремов А.Н., Кошелев В.И. и др. // ПТЭ. 2005. Т. 48. № 3. С. 46.

  3. Балдыгин В.А., Григорьев И.Н., Крученов М.Б. и др. // Изв. вузов. Физика. 2018. Т. 61. № 9/2. С. 86.

  4. Фортов В.Е., Исаенков Ю.И., Михайлов В.М. и др. // РЭ. 2013. Т. 58. № 11. С. 1102.

  5. Безруков М.Ю., Горбачев К.В., Исаенков Ю.И. и др. // Труды 17-й Межд. Крымской конф. “СВЧ-техника и телекоммуникационные технологии”, Севастополь, 10–14 сент., 2007. Севастополь: Вебер, 2007. Т. 2. С. 630.

  6. Baum C.E. // Circuit and Electromagnetic System Design Notes. Note 45. September 10, 2000.

  7. Agee F.J., Baum C.E., Prather W.D. et al. // IEEE Trans. 1998. V. PS-26. № 3. P. 860.

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

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