УДК 621.396.674.3
Леонюк А.С., Васильев Н.А.
Военная академия связи им. С.М. Буденного
Исследование широкополосных свойств трехлучевой стелющейся
антенны, размещенной в слоистой среде Арктической зоны
Аннотация. Разработка аварийных быстроразвертываемых стелющихся антенн для
коротковолновой радиосвязи в условиях Арктики является сложной задачей синтеза. Существуют
различные подходы к синтезу антенн. Среди них можно выделить параметрический, структурнопараметрический и конструкционный синтез.
В статье приводится один из способов
параметрического синтеза трехлучевой стелющейся антенны, размещенной в слоистой среде и его
результаты, представленные в виде графика зависимости коэффициента стоячей волны от
изменяемых геометрических параметров трехлучевой стелющейся антенны.
Ключевые слова: конструкционный синтез, параметрический синтез, арктическая зона, трехлучевая
стелющаяся антенна, коэффициент стоячей волны, входное сопротивление антенны.
Введение. Для ближней (примерно до
600 – 1000 км)
коротковолновой
радиосвязи пространственной волной
широко применяются горизонтальные, а
также
наклонные
и
согнутые
симметричные вибраторы. Использование
традиционных
антенно-мачтовых
устройств в зоне Островной Арктики
проблематично. Это связано с наличием в
данной физико-географической области
таких факторов природного характера, как
сильные (ураганные) ветры, скорость
которых может достигать 40 м/с, и
промерзший на большую глубину грунт,
структура которого вызывает трудности в
развертывании
антенно-мачтовых
устройств военных радиостанций. В таких
условиях
важное
место
отводится
разработке альтернативных мачтовым
антеннам новых антенных устройств,
обеспечивающих
устойчивое
функционирование линий радиосвязи.
Одним из вариантов антенных устройств
являются безмачтовые (стелющиеся)
антенны,
расположенные
непосредственно на поверхности земли.
Среда распространения радиоволн
состоит
из
верхнего
и
нижнего
полупространств.
Нижнее
полупространство
зоны
Островной
Арктики на протяжении 9…10 месяцев в
году представляет собой сложную
систему, состоящую из снежного покрова,
льда,
вечномерзлого
грунта
и
вызывающую затухание, ослабление и
отражение распространяющихся в ней
радиоволн от границ раздела сред.
Поэтому под поверхностью земли будем
понимать слоистый подстилающий объем,
состоящий из снежного покрова и вечной
мерзлоты.
Разработка
быстроразвертываемых
стелющихся антенн для КВ-радиосвязи в
условиях Островной Арктики является
сложной задачей синтеза. Существуют
различные подходы к синтезу антенн.
Среди
них
можно
выделить
параметрический,
структурнопараметрический и конструкционный
синтез [1].
Параметрический
синтез
(ПС)
позволяет при заданной структуре
проектируемой
антенной
системы
подобрать значения параметров входящих
в
нее
элементов,
при
которых
характеристики
системы
будут
находиться в заданных разработчиком
пределах
При
параметрическом
синтезе
стелющейся антенны, рассматриваемой в
данной статье, необходимо учитывать
влияние
плоскослоистых
структур
(снежный покров, мерзлый грунт), в
которые помещена антенна, на ее
электродинамические характеристики. В
качестве прототипа антенны выбрана
штатная антенна радиостанции средней
мощности Р-166 2ВНС 13/9, с некоторым
отличием, выраженным в отсутствии
мачты, к которой крепятся плечи антенны,
поперечной
перемычки
(шунта),
соединяющей
два
проводника
и
добавлении
еще
одного
плеча.
Концептуальная модель трехлучевой
стелющейся антенны (ТСА) представлена
на рис. 1.
Возбуждающий зазор
l
Рис. 1. Концептуальная модель трехлучевой
стелющейся антенны
Рассмотрим
основные
этапы
параметрического синтеза трехлучевой
стелющейся антенны.
Выбор рабочей полосы частот. Будем
руководствоваться диапазоном частот
выбранной в качестве прототипа антенны,
который составляет 5...15 МГц.
Определение
допустимых
габаритных размеров антенны. Длины
плеч ТСА соответствуют длинам плеч
2ВНC 13/9 .
Выбор численного метода. Расчет
электрических
характеристик
ТСА
целесообразно производить численным
методом. Для этих целей воспользуемся
методом
конечных
разностей
во
временной области (КРВО).
Определение
критериев
ПС.
Критерием ПС выберем КСВ антенны.
Ввиду особенностей размещения ТСА
питание целесообразно осуществлять
коаксиальным кабелем, изготавливаемым
промышленностью
с
волновым
сопротивлением ф 150 Ом . Потребуем,
чтобы КСВ 2,5 при ф 150 Ом .
Ввод переменных и пределов их
изменения. В качестве переменной ПС
выберем угол между плечами антенны
27 , 45 , 63 . Величины
углов
продиктованы
особенностями
метода КРВО, заключающимися в
разбиении
расчетной
области
на
кубические
ячейки.
Расчет
КСВ
целесообразно
производить,
предварительно
определив
входное
сопротивление антенны.
Методика
расчета
входного
сопротивления
трехлучевой
стелющейся антенны, размещенной в
слоистой среде. Рассмотрим этапы
предлагаемой методики.
1)Введение исходных данных. На
данном этапе вводятся: частота входного
сигнала f, его амплитуда A, число
частотных отсчетов N f , шаг по частоте
f , длина плеч вибратора l, размеры
трехмерной
расчетной
области,
геометрические размеры слоев, в которых
расположена антенна и физические
константы (скорость света, абсолютная
диэлектрическая
проницаемость,
абсолютная магнитная проницаемость).
2) Расчет
параметров
модели
возбуждения и определение приращений
по
времени.
Исходными
данными
являются амплитуда гармонического
сигнала A и количество временных шагов
в периоде колебания Nш. Подлежат
вычислению количество синусоид за
время
Nsin f / f ,
воздействия
максимальное число временных шагов
Tstep Nsin Nш , частота дискретизации
сигнала
временного
f д f Nш
шага
и
t 1
величина
fд .
Число
синусоид определяется в зависимости от
времени, необходимого для достижения
установившегося режима, означающего
синусоидальную повторяемость всех
полей
дискретизируемой
области.
Точность и устойчивость алгоритма
пропорциональна количеству временных
шагов в периоде колебания. Амплитуда
гармонического сигнала A может быть
произвольной.
3) Вычисление
параметров
гармонического сигнала. Данный этап
необходим для установления соответствия
между
комплексными
амплитудами,
полученными
после
процедуры
преобразования Фурье (ПФ) и заданной
частотой
входного
воздействия.
Определение порядка распределения
частот в векторе дискретного Фурьеобраза осуществляется при использовании
частоты Найквиста, т. е. максимальной
частоты в спектре сигнала, которая может
быть представлена дискретным ПФ
[2].
Определение
n 1
2 t
комплексной
амплитуды
искомой
величины для заданной частоты сигнала
происходит при учете, что первое
значение в вычисляемом векторе Фурьеобраза соответствует частоте f = 0 Гц.
Таким образом, порядковый номер
комплексной амплитуды в векторе
дискретных значений будет равен:
n f
1.
f
4) Определение
приращений
по
пространству. В трехмерной решетке для
обеспечения пошаговой устойчивости
алгоритма во времени, шаг временной
дискретизации t выбирают согласно
условию
Куранта – Фридрихса – Леви:
x t 3 c .
5) Построение временной и сигнальной
шкал:
t (i) t (i 1) t, i 1, 2,...,T ,
step
где t – вектор значений временных шагов.
Источник возбуждения симметричного
вибратора моделируется следующим
образом:
s A sin(2ft ),
где A – амплитуда входного воздействия;
f – частота входного воздействия.
6) Автоматический ввод в алгоритм
модели антенны в слоистой среде.
Согласно введенным геометрическим
размерам сред, антенны, выбранным
приращениям по пространству x , y и
z аппроксимируется антенна в слоистой
среде. В качестве антенны выбрана ТСА.
7) Составление
матрицы
электрических параметров и присвоение
их величин моделируемым средам. На
данном этапе осуществляется ввод
параметров
всех
задействованных
материалов и присвоение физических
величин, характеризующих их свойства,
объектам модели:
Material(m, p) = 1,
где m – порядковый номер материала; p –
физическая величина, характеризующая
свойства среды (материала); после знака
равенства показано численное значение
данной физической величины.
После того как будут заданы все
материалы, производится присвоение
характеризующих их физических величин
объектам модели:
Index(I, J, K) = m,
где Index(I, J, K) – массив размером
I J K .
8) Ввод идеально согласованных слоев.
Для
имитации
бесконечного
полупространства,
примыкающего
к
границам расчетной области, в методике
использована техника формирования
идеально согласованных слоев (ИСС),
согласно которой каждая составляющая
напряженности электромагнитного поля
(ЭМП) разбивается на две компоненты
[3].
9) Модель
возбуждения
антенны.
Напряжение, прикладываемое к зажимам
11) Вычисление I А t . Ток на зажимах
вертикальная
симметричного вибратора определяется
по закону Ампера в конечно-разностной
форме [4], т. е. численным решением.
Таким образом, получаем множество
мгновенных значений тока на каждом
составляющая вектора напряженности
электромагнитного поля в точке питания
антенны, узле с координатами I, J, K; s –
источник синусоидального возбуждения;
T – номер итерации (рис. 2);
12) ПФ
полученных
мгновенных
значений
тока
и
подводимого
напряжения. В рассматриваемом случае
непрерывные функции I А (t ) и U А (t )
антенны
U АT ,
выражается
через
напряженность электрического поля:
EzT ( I , J , K ) s(T 1) ,
где
–
EzT ( I , J , K )
U АT s(T 1) z,
где z – шаг
пространству.
дискретизации
по
J-1 J
z
K
y
K−1
x
y
Зазор
EzT ( I , J , K )
Рис. 2. Модель возбуждения антенны
10) Расчет методом КРВО E и H
составляющих напряженности ЭМП.
Методом КРВО вычисляются значения
магнитных составляющих напряженности
H , B , H , B
x
x
y
y
H z , Bz на текущем
временном отсчете T во всей расчетной
области.
Далее
методом
КРВО
вычисляются значения электрических
составляющих
напряженности
ЭМП
E ,G , E ,G
x
x
y
y
, Tstep .
и
Результатом
процедуры численного ПФ является
вектор значений напряжения и тока в
комплексном
виде,
т. е.
K+1
Плечи антенны
, Tstep
U Am U АT , T 1, 2,
J+1
, Tstep .
определены на интервале [0, Tstep],
поэтому дискретным аналогом данных
функций
будем
называть
векторы
I Am I АT , T 1, 2,
z
ЭМП
T
временном шаге T I А , T 1, 2,
Ez , Gz во всей расчетной
области на временном шаге T + 1 с
использованием обновленных значений
H x , Bx , H y , By H z , Bz . Процесс вычислений
повторяется и заканчивается при условии,
если T = Tstep. Вычисление методом
конечных разностей во временной
области всех составляющих в данном
случае – решение задачи распространения
электромагнитной волны.
U Am F(U AT , T 1, 2,
I Am F( I AT , T 1, 2,
, Tstep ) ,
, Tstep ) . Комплексные
амплитуды
тока
и
подводимого
напряжения, соответствующие заданной
частоте, определяются согласно пункту 3.
13) Вычисление RА и X А . Получив
значения комплексных амплитуд тока и
напряжения на зажимах вибратора,
рассчитывают входное сопротивление
U Аm
ZА
.
I Аm
Заключение.
По
предложенной
методике произведен расчет входного
сопротивления ТСА с разными углами
между
плечами,
размещенной
в
трехслойной среде (воздух, снежный
покров и мерзлый грунт) с параметрами:
ɛr2 = 1,2; σ2 = 1·10–6 См/м и ɛr3 = 4,2;
σ3 = 1·10–4
См/м,
соответствующими
зимнему периоду в зоне Островной
Арктики. Антенна располагалась на
границе раздела воздуха и снежного
покрова. КСВ рассчитывается согласно
1
выражению: KСВ
1
RА iX А ф
RА iX А ф
RА iX А ф
.
RА iX А ф
Результаты расчета КСВ представлены
на рис. 3.
КСВ
4,5
α=27°
4
α=45°
3,5
α=63°
3
2,5
2
1,5
Полученные
результаты
свидетельствуют о том, выбранная в
качестве критерия характеристика КСВ
находится в заданных пределах для
исследуемого диапазона частот при угле
между проводниками ТСА 45 . В
качестве
критерия
могут
быть
использованы и другие электрические
характеристики антенны, такие как,
например,
коэффициент
усиления,
коэффициент полезного действия. В
дальнейшем предполагается изменять не
только угол между проводниками, но и их
длины.
1
0,5
0
5
6
7
8
9
10
11
12
13
14
f, МГц
Рис. 3. Зависимость КСВ ТСА от частоты
СПИСОК ЛИТЕРАТУРЫ
[1] THE PRINCIPLES OF ANTENNAS CONSTRUCTIVE SYNTHESIS IN DISSIPATIVE MEDIA Borodulin
[2]
[3]
[4]
R.U., Sosunov B.V., Makarov S.B. Lecture Notes in Computer Science. 2017. Т. 10531 LNCS. С. 455-465.
Teukolsky, S. A. Press Numerical recipes in C. The Art of Scientific Computing/ S.A. Teukolsky, T. William//
Cambridge: Cambridge University Press – s.l.: 2-nd edition, 2002.
Berenger, J. P. Three-Dimensional Perfectly Matched Layer for the Absorption of Electromagnetic Waves/ J. P.
Berenger// J. Comp. Phys. – 1996. – Vol. 127. – Pp. 363–379.
Сальников, Д. В. Методика расчета входного сопротивления несимметричного вибратора методом
КРВО/ Д. В. Сальников// Научно-технический сборник: сб. ст./ Труды академии/СПБ., 2005.
Отзывы:
Авторизуйтесь, чтобы оставить отзыв