Preview

Вестник Концерна ВКО «Алмаз – Антей»

Расширенный поиск

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

https://doi.org/10.38013/2542-0542-2022-1-25-32

Полный текст:

Аннотация

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

Для цитирования:


Жизневский А.К. Сравнительный анализ численных методов пересчета электрического поля на плоскость, не параллельную исходной. Вестник Концерна ВКО «Алмаз – Антей». 2022;(1):25-32. https://doi.org/10.38013/2542-0542-2022-1-25-32

For citation:


Zhiznevsky A.K. Comparative analysis of numerical methods to calculate electric field distribution from the initial plane to a plane non-parallel to the initial one. Journal of «Almaz – Antey» Air and Space Defence Corporation. 2022;(1):25-32. https://doi.org/10.38013/2542-0542-2022-1-25-32

ВВЕДЕНИЕ

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

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

Уменьшение временных затрат на поиск пространственного положения облучателя, близкого к оптимальному, при котором все радиотехнические характеристики ФАР удовлетворяют заданным требованиям, может быть достигнуто путем замены ряда промежуточных экспериментов расчетами, осуществляемыми по специально разработанным алгоритмам.

При разработке этих алгоритмов необходимо решить ряд электродинамических задач прикладной направленности.

Одна из таких задач – пересчет амплитудно-фазового распределения, измеренного на плоской поверхности сканирования P0 (рис. 1), параллельной раскрыву рупорного облучателя, на плоскость входов фазовращателей проходной ФАР P1 при различных углах поворота Θ облучающего рупора.

Рис. 1. Иллюстрация постановки задачи

В более общем виде эту задачу можно сформулировать следующим образом: «пересчет известного амплитудно-фазового распределения поля на исходной плоскости P0 на плоскость P1, не параллельную исходной». Решение этой задачи позволяет избежать проведения дополнительных измерений поля на плоскости сканирования PI0, используя одно измерение, полученное при сканировании на плоскости P0, для всех возможных углов отклонения Θ.

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

Для простоты будем рассматривать только линейную поляризацию вектора напряженности электрического поля и разместим его вдоль оси Оy, что не приведет к потере общности, так как результаты для случая линейной поляризации вектора E вдоль оси Ox и круговой поляризации могут быть получены поворотом системы координат (СК). Размеры плоскости сканирования P0 должны быть таковы, чтобы при рассматриваемом положении облучателя амплитуда поля на ее краях составляла не более –30 дБ от максимального значения [1].

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

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

Решетка диполей расположена в плоскости Oxy. P0 – плоскость, параллельная решетке диполей, находящаяся на расстоянии z0 = 5λ (где λ – длина волны) от начала координат. P1 – плоскость, центр которой находится на расстоянии z1 = 100λ от решетки диполей. Она повернута вокруг оси Оy на угол Θy, равный 10 градусам. Размеры плоскостей P0 и P1 одинаковы и равны 124,5λ на 127λ (250 на 255 точек с шагом 0,5λ).

Электрическое поле от модельной решетки диполей рассчитывается на плоскость P0 (рис. 3). Затем с плоскости P0 пересчитывается на повернутую плоскость P1 каждым из исследуемых методов. Для проведения оценки результатов в качестве эталона используется поле, рассчитанное на повернутой плоскости P1 непосредственно от решетки диполей. Размер модельной дипольной решетки 20 элементов (вдоль оси Ox) на 8 элементов (вдоль оси Oy) излучателей. Длина каждого диполя равняется половине длины волны. Расстояние между диполями – половина длины волны. Диполи ориентированы вдоль оси Oy, что соответствует вектору E вдоль оси Oy. Возбуждение решетки – равноамплитудное с перепадом фазы на 180 градусов при x = 0, что позволяет получить разностную диаграмму направленности (ДН) от поля на плоскости для оценивания точностных характеристик рассматриваемых численных методов.

Рис. 3. Электрическое поле на плоскости P0,
полученное от модельной решетки диполей

1. Интерполяционный метод

Для расчета электрического поля на плоскости P1, не параллельной исходной P0 и повернутой к ней на угол Θy, можно воспользоваться интерполяционным методом (рис. 4).

Рис. 4. Интерполяционный метод

Для реализации данного способа необходимо найти амплитудно-фазовое распределение (АФР) на каждой из промежуточных плоскостей P1Пр, P2Пр, ..., PnПр, которое создается электрическим полем исходной плоскости P0. Число плоскостей n зависит от расстояния между ними. Величина dдис выбирается из условия, что расстояние между соседними точками пересечения искомой плоскости P1 и параллельных плоскостей P1Пр, P2Пр, ..., PnПр (узловыми точками) не должно быть больше чем λ/2.

Следовательно, dдис ≤ λ sinΘy / 2. В общем случае, когда поворот плоскости P1 происходит вокруг оси Oz или вокруг двух и более осей, необходимо произвести интерполяцию из узловых точек на регулярную сетку [2]. В рассматриваемом случае точки пересечения совпадают с регулярной сеткой плоскости P1 и интерполяция не требуется.

2. Метод эквивалентных поверхностных токов (метод Гюйгенса – Кирхгофа)

Реализация данного метода основана на применении в расчетах электрического поля эквивалентных электрических и магнитных диполей. Электрические и магнитные поля заменяются соответственно на магнитные и электрические моменты.

Известны формулы [3] для электрического поля, создаваемого электрическим и магнитным диполем:

(1)

(2)

где эквивалентные электрический и магнитный моменты пропорциональны магнитному и электрическому полям соответственно [3]:

(3)

(4)

где HSm(x, y, z0)вычисляется с использованием аппарата теории плоских волн [1], [4], [5].

В формулах (1), (2) диполи ориентированы вдоль оси Oz.

Так как рассматривается линейная поляризация вектора напряженности электрического поля Ey, то эквивалентный магнитный диполь также направлен вдоль оси Oy, а эквивалентный электрический диполь (соответствующий Hx) – вдоль оси Ox. Для учета этого условия в расчетах применяются повороты СК. Формулы (1)–(4) позволяют рассчитать поле на плоскости P1, зная поле на плоскости P0.

3. Метод поворотной трансформации спектра плоских волн

В данном методе задаются две СК [6]. Первая СК выступает в роли исходной СК (x, y, z), в которой находится излучающая плоскость P0 с координатами (x, y, z0). Во второй СК (, ŷ, ) находится искомая плоскость P1 с координатами (, ŷ, 1).

Единичные векторы r = (x, y, z) и = (, ŷ, ) могут быть найдены с помощью матрицы поворота T следующим образом:

(5)

(6)

Матрица поворота T (7) вокруг трех осей в трехмерном пространстве представляет собой произведение трех матриц поворота вокруг каждой из осей, а в случае если поворот производится вокруг одной оси, то используется одна соответствующая повороту матрица: (8), (9) или (10).

(7)

где Wxx), Wyy), Wzz) – матрицы поворотов вокруг осей Ox, Oy и Oz соответственно, θx, θy, θz – углы поворота вокруг осей Ox, Oy и Oz против часовой стрелки для наблюдателя, смотрящего с конца оси.

Матрицы вращения вокруг каждой из осей могут быть описаны следующим образом:

(8)

(9)

(10)

При этом матрицы вращения обладают следующими характеристиками:

(11)

где Wζ –1(χ) и WζT(χ)– обратная и транспонированная матрицы вращения соответственно.

В общем виде эти матрицы могут быть представлены следующим образом:

(12)

(13)

Если электрическое поле E(x, y, z) с волновым числом k задано на плоскости P0, то электрическое поле на произвольных плоскостях, параллельных Oxy, задается в следующем виде [1], [6]:

(14)

где kx, ky, kz – проекции волнового вектора, kz = . а F(kx, ky) – спектр плоских волн электрического поля E(x, y, z), вычисляемый как:

(15)

Зададим новую СК так, чтобы ее оси были параллельны одноименным осям системы Oxy, расстояние от ее начала до плоскости P1 равнялось нулю, а центр плоскости P0 имел координаты (0, 0, –(z1 – z0)). Это позволит обнулить z в формуле (14), при этом в формуле (15) будет произведена подстановка d = –(z1 – z0).

Электрическое поле (x, y, z) на повернутой плоскости P1 находится путем подстановки (6) и (12) в (14):

(16)

Векторы f(kx, ky, kz) и (kˆx, kˆy, kˆz) в пространстве Фурье преобразуются один в другой также с помощью матриц T и T–1. Поэтому, используя соотношение = Tf, аргумент экспоненты в выражении (16) можно упростить до 2πi(x+ y ŷ).

При этом электрическое поле на повернутой плоскости P1 будет определяться выражением:

(17)

где проекции волнового вектора kx, ky в уравнении (16) выражаются с использованием соотношения f = T–1 следующим образом:

(18)

(19)

Замена в интеграле (16) переменных kx, ky на x, y осуществляется с помощью применения Якобиана:

(20)

где

(21)

С использованием (20) может быть реализовано преобразование электрического поля к виду:

(22)

где

(23)

Здесь

(24)

Графические иллюстрации по результатам моделирования процесса применения метода трансформации спектра приведены на рисунке 5.

Рис. 5. Спектр плоских волн до трансформации и после

4. Сравнение численных методов

Результаты анализа эффективности использования вышеописанных численных методов приведены на рисунке 6, где представлены азимутальные сечения ДН, рассчитанные с шагом 0,01 градуса.

Основное внимание уделено области одного из двух главных лепестков ДН и области равносигнального направления. Приведенные графики позволяют судить о степени приближения местоположения облучателя в пространстве к его оптимальному значению.

В таблице 1 приведены численные параметры полученных ДН.

Таблица 1. Сравнительный анализ методов

Временные затраты на вычисления по каждому из методов были измерены при различном числе отсчетов на плоскостях P0 и P1. Соответствующий график представлен на рисунке 7. Из этого рисунка также видно, что метод интерполяции, как и ожидалось, предъявляет повышенные требования к объему оперативного запоминающего устройства (ОЗУ) ЭВМ. При достижении количеством элементов порога в 70 000 единиц происходит прекращение дальнейшего расчета из-за переполнения ОЗУ при вызове сильно расходующей память функции интерполяции. Объем доступной ОЗУ ЭВМ при выполнении численного моделирования составлял 11 ГБ.

ВЫВОДЫ

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

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

1. Шубников В. В., Калашников В. С. и др. Синтез характеристик антенн по измерениям в ближней зоне. СПб.: ГУАП, 2016. 309 с.

2. https://octave.sourceforge.io/octave/function/interp2.html

3. Никольский В. В., Никольская Т. И. Электродинамика и распространение радиоволн. Изд. 3-е, перераб. и дополн. М.: Наука, 1989. 544 с.

4. Evans G. E. Antenna measurement techniques // Norwood. 1990.

5. Slater D. Near-field antenna measurements. Artech House. 1991.

6. Matsushima K. Formulation of the rotational transformation of the fields and their application to digital holography // Applied optics. 2008. T. 47. №. 19. C. D110–D116.


Об авторе

А. К. Жизневский
Научно‑образовательный центр акционерного общества «Северо-Западный региональный центр Концерна ВКО «Алмаз – Антей»
Россия

Жизневский Александр Константинович – аспирант Научно-образовательного центра акционерного общества «Северо-Западный региональный центр Концерна ВКО «Алмаз – Антей» – Обуховский завод»; заместитель начальника цеха, АО «ЗРТО». Область научных интересов: измерения в ближней зоне, настройка антенных решеток.

Санкт-Петербург



Рецензия

Для цитирования:


Жизневский А.К. Сравнительный анализ численных методов пересчета электрического поля на плоскость, не параллельную исходной. Вестник Концерна ВКО «Алмаз – Антей». 2022;(1):25-32. https://doi.org/10.38013/2542-0542-2022-1-25-32

For citation:


Zhiznevsky A.K. Comparative analysis of numerical methods to calculate electric field distribution from the initial plane to a plane non-parallel to the initial one. Journal of «Almaz – Antey» Air and Space Defence Corporation. 2022;(1):25-32. https://doi.org/10.38013/2542-0542-2022-1-25-32

Просмотров: 437


Creative Commons License
Контент доступен под лицензией Creative Commons Attribution 4.0 License.


ISSN 2542-0542 (Print)