CC BY
131
Системы управления,связи и безопасности №2. 2019
Systems of Control, Communication and Security ISSN 2410-9916
УДК 621.396.96
Метод снижения радиолокационной заметности воздушных судов путем параметрического синтеза оптимальной формы канала воздухозаборников
Кузнецов В. А.
Постановка задачи: задача снижения радиолокационной заметности воздушных судов является актуальной во всем мире, а ее решение состоит из комплекса превентивных мер по снижению эффективной площади рассеяния и, кроме того, сокращению информативности демаскирующих признаков. Существенный вклад в суммарную эффективную поверхность рассеяния в передней полусфере вносит воздухозаборник. Вместе с тем, наличие эффекта вторичной доплеровской модуляции при отражении электромагнитных волн от воздухозаборника обуславливает необходимость снижения его радиолокационной заметности. Целью работы является снижение радиолокационной заметности воздушного судна путем параметрического синтеза Б-образной формы канала его воздухозаборников, оптимальной по критерию минимума эффективной площади рассеяния в заданном диапазоне конструктивно возможных изменений. Используемые методы: при решении задачи снижения радиолокационной заметности используется параметрический синтез формы канала, для оценки эффективной площади рассеяния для учета поляризации падающей электромагнитной волны модифицирован метод геометрической оптики, при этом широко используется математическое моделирование. Новизна: элементами новизны является непосредственно сам метод параметрического синтеза оптимальной по критерию минимальной при всех возможных условиях наблюдения эффективной площади рассеяния формы канала воздухозаборника, а также модификация метода геометрической оптики для учета поляризационных свойств каждого фацета поверхности канала. Результат: применение предложенного метода позволяет снизить эффективную площадь рассеяния канала воздухозаборника воздушного судна во всем диапазоне углов облучения более чем на 90% по сравнению с прямой формой канала, а также синтезировать форму Б-образного канала воздухозаборника с требуемой точностью при параметрически заданных форме входной апертуры и диапазоне конструктивно возможных искривлений геометрии канала. Практическая значимость: представленное решение задачи снижения радиолокационной заметности может быть использовано при оценке эффективной площади рассеяния с учетом реальной нагрузки канала в виде лопаток входного направляющего аппарата и компрессоров авиационного двигателя, а также при разработке макетов воздухозаборников перспективных воздушных судов.
Ключевые слова: радиолокационная заметность, эффективная площадь рассеяния, геометрическая оптика, фацетная модель, воздухозаборник.
Введение
Для снижения боевых потерь авиационной техники во всем мире разрабатывается и применяется целый комплекс превентивных мер по снижению ее радиолокационной (РЛ) заметности (рис. 1), единственным показателем которой является эффективная площадь рассеяния (ЭПР).
Библиографическая ссылка на статью:
Кузнецов В. А. Метод снижения радиолокационной заметности воздушных судов путем параметрического синтеза оптимальной формы канала воздухозаборников // Системы управления, связи и безопасности. 2019. № 2. С. 180-202. DOI: 10.24411/2410-9916-2019-10209. Reference for citation:
Kuznetsov V. A. The radar aircraft visibility reducing method by using air intake duct optimal shape parametric synthesis. Systems of Control, Communication and Security, 2019, no. 2, pp. 180-202. DOI: 10.24411/2410-9916-2019-10209 (in Russian).
DOI: 10.24411/2410-9916-2019-10209
Systems of Control, Communication and Security
ISSN 2410-9916
Рис. 1. Меры по снижению радиолокационной заметности воздушного судна
В общем случае, конструктивные меры можно разделить на четыре направления: специальная форма фюзеляжа (shaping), радиопоглощающие материалы (РПМ или RAM - radar absorbing materials), активные и пассивные способы противодействия (active and passive cancellation) [1], причем к каждой составной части самолета применяется, по возможности, весь комплексный подход.
В открытых зарубежных и отечественных источниках, например, в [1-3], представлено огромное количество диаграмм обратного вторичного излучения (ДОВИ) различных воздушных судов (рис. 2).
а) б) в) г) д)
Рис. 2. Примеры воздухозаборников и диаграмм обратного вторичного излучения самолетов F-15 (а), Су-27 (б), F-22 (в), F-35 (г), Су-57 (д)
В них отмечается существенный вклад воздухозаборников сверхзвуковых самолетов 4 и старше поколений в ЭПР в передней полусфере, и только кромок крыла - у самолетов, выполненных по технологии Stealth. Результаты их анализа показывают, что используемые подходы к моделированию и расчету ДОВИ в большинстве своем не учитывают элементы авиационной силовой установки [4]. Кроме того, расчет суммарных значений ЭПР производится только по ракурсу, без учета разрешающей способности РЛС, что приводит к невозмож-
DOI: 10.24411/2410-9916-2019-10209
Системы управления,связи и безопасности №2. 2019
Systems of Control, Communication and Security ISSN 2410-9916
ности формирования с их помощью дальностно-доплеровских радиолокационных портретов [5].
Особый интерес представляют элементы авиационной силовой установки, поскольку целью снижения их радиолокационной заметности является сокращение набора признаков, используемых в бортовых радиолокационных системах (БРЛС) при решении задач распознавания воздушных целей. В данном контексте речь идет об эффекте вторичной доплеровской модуляции, возникающем при отражении электромагнитных волн (ЭМВ) от движущихся составных частей объектов локации, а именно от вращающихся лопаток компрессора низкого давления, вследствие чего по значениям частот, составляющих спектр такой модуляции, появляется возможность распознать тип авиационной силовой установки, а, значит, и непосредственно саму воздушную цель. Для этого применяются экранирующие щитки, выполняющие функцию клина, Б-образная форма канала воздухозаборников и входной направляющий аппарат (ВНА) двигателя, играющий роль так называемого радар-блокера, а также радиопо-глощающие покрытия (РПП) на всех этих элементах (рис. 3).
а) б) в)
Рис. 3. Пример конструктивных решений по снижению радиолокационной заметности элементов авиационной силовой установки: входной апертуры воздухозаборника (а), Б-образного канала воздухозаборника (б), входного направляющего аппарата (в)
Что касается Б-образной формы канала воздухозаборника сверхзвукового самолета, то разработчикам приходится параллельно решать две задачи: с одной стороны, необходимо обеспечить высокую газодинамическую устойчивость канала для исключения возможности помпажа авиационной силовой установки, и, с другой стороны, удовлетворить требования заказчика по заданному уровню радиолокационной заметности в передней полусфере. Поэтому разработка оптимальной геометрии канала воздухозаборника, удовлетворяющей одновременно двум критериям, а именно минимальной радиолокационной заметности и высокой газодинамической устойчивости, является сложной и актуальной на сегодняшний день научно-технической задачей на стыке нескольких научных направлений. Несмотря на важность удовлетворения обоим критериям, в данной статье внимание уделяется только радиолокационной замет-ности.
DOI: 10.24411/2410-9916-2019-10209
Системы управления,связи и безопасности №2. 2019
Systems of Control, Communication and Security ISSN 2410-9916
Целью данной работы является снижение радиолокационной заметности канала воздухозаборника путем параметрического синтеза S-образной формы канала, оптимальной по критерию минимума ЭПР в заданном диапазоне конструктивно возможных изменений.
Постановка задачи
Как известно, единственным средством обнаружения воздушных целей на больших дальностях является БРЛС. Для сокращения дальности обнаружения самолетов радаром противника применяется технология снижения их радиолокационной заметности. Так, например, на рис. 4 приведены зависимости дальности обнаружения от ЭПР, рассчитанные при заданных импульсной мощности излучения БРЛС Римп=1 и 20 кВт, диаметре антенны 0,9 м, коэффициенте полезного действия антенны 0,7, эффективной площади антенны 0,45 м2, длительности импульса 10 мкс, количестве интегрируемых импульсов 100, частотах сигнала 1, 3 и 10 ГГц, общих потерях 3 дБ. Очевидно, что для увеличения дальности обнаружения малозаметных воздушных целей с ЭПР менее 1 м2 требуется существенно увеличивать энергетику излучения БРЛС, а это, в большинстве случаев, либо невозможно, либо ведет к снижению скрытности и демаскировке самолета-носителя БРЛС.
300
5 250
^ 200 J3
о
о 150 л С га
4100 о га
^ 50
0'.......'........'..........'-1,, ,-,- —г-т =
101 10° 10"1 10"2 10"3 10"4
ЭПР, м.кв.
Рис. 4. Зависимости дальности обнаружения от ЭПР воздушной цели
С другой стороны, уменьшение средних по ракурсам значений ЭПР воздушного судна обеспечивает значительное снижение боевых потерь (рис. 5) [6].
Для оценки максимальной дальности обнаружения воздушной цели необходимо знать ее среднюю или суммарную ЭПР в направлении облучения. В случае одиночной или сосредоточенной цели чаще всего прибегают к моделированию ДОВИ статистическими методами на основе экспериментальных данных [7]. Однако в этом случае отсутствует возможность учесть или оценить вклад воздухозаборников в суммарную ЭПР цели. Для учета отдельных, в том числе движущихся, конструктивных элементов и их электромагнитных свойств наиболее широкое распространение получили фацетные модели представления поверхности объекта сложной формы и лучевая концепция моделирования процессов рассеяния и многократного переотражения ЭМВ [8].
DOI: 10.24411/2410-9916-2019-10209
Systems of Control, Communication and Security
ISSN 2410-9916
Относительные потери Относительные потери ю|- _
010510 20 30 40 5€ ЭПР=15м; ЭПР=15м: заметкостью
^ПР 'г н с рэп
ЭПР = 1 м
Рис. 5. Зависимости относительных боевых потерь от ЭПР самолета
В данной работе фацетная модель канала воздухозаборника представляет собой вытянутую произвольную полость, поверхность которой описывается идеально проводящими плоскими пластинами треугольной формы - фацетами, заданными тремя вершинами с координатами в трехмерном пространстве. Каждый фацет при этом характеризуется вектором нормали и электромагнитными параметрами. В качестве априорной информации может быть задан радиус выходной апертуры R (соответствующий радиусу двигателя), длина канала L, диапазоны конструктивно возможного изменения геометрии канала, электромагнитные характеристики материалов.
Задача синтеза S-образной формы канала воздухозаборника состоит в следующем: необходимо аналитически задавать форму входной апертуры, закон расширения или сужения поперечных сечений вдоль длины канала и характер изогнутости канала в двух плоскостях (вертикальной и горизонтальной). Такое параметрическое задание геометрии формы канала и ее изменение в некотором диапазоне значений с учетом их пространственного ограничения позволит производить расчет значений ЭПР, и, в случае их минимальных значений, запоминать оптимальные параметры геометрии канала с последующим синтезом его поверхности.
Синтез входной апертуры воздухозаборника
На сегодняшний день известно множество используемых форм входной апертуры воздухозаборников: от скошенных квадратных и прямоугольных до треугольных и трапециевидных (рис. 6).
Одним из немаловажных критериев выбора формы входной апертуры канала является ограничение прохождения вертикально или горизонтально поляризованной ЭМВ в канал. Тенденцией последних разработок является размещение апертуры канала в верхней части фюзеляжа для исключения её видимости наземными радиолокационными системами, однако такое конструктивное решение ведет к существенному снижению маневренности воздушного судна. Кроме того, на самолетах F-35 Lightning II и J-20 Chengdu перед входной апертурой канала используют «наплывы» (рисунок 5, в и ж) для исключения попадания ЭМВ, переотраженной от фюзеляжа, путем ее рассеяния.
DOI: 10.24411/2410-9916-2019-10209
Systems of Control, Communication and Security
ISSN 2410-9916
ж) з) и)
Рис. 6. Примеры форм входной апертуры воздухозаборника Су-57 (а), F-22 (б), F-35 (в), ТР-23 (г), X-47B (д), X-47A (е), J-20 (ж), Qaher F-313 (з), проекта ПАК ДА (и)
Моделирование выходной апертуры канала воздухозаборника в трехмерном пространстве в форме окружности не представляет особой сложности. Пусть вытянутая сторона канала направлена по оси У, длина канала Ь=тХ, где т - коэффициент пропорциональности, X - длина волны, а число углов, задающих угловое положение точек апертуры, задано Ыа, тогда система уравнений, с помощью которых рассчитываются координаты точек апертуры выглядит следующим образом:
^ = Я ^(а)
У = ь ,
радиус двигателя, п -а. = а1_1 + За - значение текущего угла, За
На рис. 7 приведены примеры вычисленных согласно (1) координат точек входной апертуры воздухозаборника с различными значениями параметра Жа.
где R=nX -
(1)
коэффициент пропорциональности,
= 360\..e,/ = l.JVa.
а) б) в)
Рис. 7. Пример сформированных массивов точек выходной апертуры канала воздухозаборника радиусом Я=1 м при N„=10 (а), N„=50 (б), N,=100 (в)
Для задания входной апертуры канала в форме выпуклых многоугольников интересное решение нашлось в работе [9] как для прямоугольной, так и для полярной системы координат. Совместное использование параметрического
DOI: 10.24411/2410-9916-2019-10209
Systems of Control, Communication and Security
ISSN 2410-9916
представления и универсального уравнения позволило получить описание различных двумерных фигур в параметрическом виде. Добавление третьей координаты к выражениям из [9] позволило однозначно задать любой прямоугольник в трехмерном пространстве системой параметрических уравнений:
Ч- = ((А + Аз;) А + АА
У = Ь , (2)
^ = ((А1- + А4;) А + АзА
- ключи, определяющие последовательность вычислений,
где А,-A5i
4, =
1
T
_i_
d
T _i_
d
Au =
T
_l_
d
T d
2
A =
T
_i_
d
d
A4l =
T
_i_
d
T d
4
A. = T - d •
T _i_
d
[*] - оператор округления к меньшему,
d=kX - длина диагонали прямоугольника, к - коэффициент пропорциональности, С - угол наклона диагонали относительно горизонта, Т - параметр, изменяемый в пределах от d до 5d с шагом / (Ыа -1).
В случае, если С =п/4, то система уравнений (2) описывает квадраты разного размера в зависимости от величины заданной диагонали ^
Поскольку прямоугольник представляет собой частный случай параллелограмма, то, соответственно, он является одновременно и параллелограммом, который в общем случае может быть получен путем введения в (2) величины деформирования Адеф:
x =
А
деф
d
Ai, + Ai cos(^)
A5l + A2l {Адеф + At COS(^))
+
+A
L3l
d cos(^) +
Адеф
A5l
d
(3)
/
У- = 0
3 =( (А- + Аз-) А + А2гА) ял(С) Очевидно, что если Адеф =0, то система уравнений (3) будет определять
параллелограмм в виде прямоугольника. Если прямоугольник деформировать на величину Адеф одновременно с двух сторон, то будет получена равнобедренная трапеция, которая однозначно задается системой уравнений, полученной из системы уравнений (2) путем введения двусторонней деформации:
3
<
DOI: 10.24411/2410-9916-2019-10209
Systems of Control, Communication and Security
ISSN 2410-9916
X,.
+A
А^ A + a/+ ^ A + A„ (Ад. + A costf)) +
d
4i
r
3i
d А
\
dcos(^) + 2АдеЛ --^A
деф
d
• (4)
У = 0
3 = ( (А + А ) А + )
На рис. 8 приведены некоторые результаты моделирования формы входной апертуры канала воздухозаборника с количеством точек N«=100 и соответствующими параметрами для систем уравнений (1)-(4).
А =0 С=л/4 d=2,02 м
а)
А =0,5 z=n/3,1 d=2,88 м б)
А =0,9 Z=n/8 d=2,04 м в)
А =0,5 z=n/3,1 d=2,56 м
г)
А =1,9 Z=n/2 d=1,42 м
д)
Рис. 8. Синтезированные массивы точек входной апертуры воздухозаборника различной формы при N«=100
Как видно из представленных рисунков, параметрическое задание многоугольников позволяет синтезировать идентичные реальным формы входной апертуры канала воздухозаборников современных воздушных судов со сниженной радиолокационной заметностью, представленные на рис. 6. Следует отметить, что в дальнейшем, при необходимости, «срезанность» или наклон плоскости входной апертуры могут быть реализованы несколькими способами: путем нахождения точек пересечения поверхности канала с некоторой плоскостью среза или изогнутой поверхностью фюзеляжа, либо простейшим наклоном входной апертуры путем домножения координат точек апертуры на заданный косинус или синус требуемого угла наклона или среза.
В соответствии с целью работы и сформулированными для ее достижения задачами, на следующем этапе синтеза необходимо формализовать изогнутость по двум осям и расширение (сужение) промежуточных сечений по всей длине воздухозаборника и синтезировать соответствующую заданным параметрам геометрию поверхности канала в виде множества вершин фацетов.
Синтез прямого канала воздухозаборника
Существует множество способов и алгоритмов задания поверхности вытянутого по оси вращения объекта, в том числе линейными и нелинейными функциями (например, [10]). Однако простейшим в смысле вычислительной
<
DOI: 10.24411/2410-9916-2019-10209
Системы управления,связи и безопасности №2. 2019
Systems of Control, Communication and Security ISSN 2410-9916
сложности является соединение точек входной и выходной апертур отрезками прямой с последующей интерполяцией промежуточных точек канала путем деления отрезков на равные части, число которых соответствует заданному числу сечений Ыь, как показано на рис. 9.
а) б) в)
Рис. 9. Синтезированные прямые каналы воздухозаборников при #¿=15
Основным требованием для адекватного формирования поверхности канала предлагаемым способом является соответствие числа точек, образующих входную и выходную апертуры, заданному значению Ыа. По промежуточным точкам, образованным таким образом, можно рассчитать площадь сечений для оценки газодинамической устойчивости. Очевидно, что чем больше заданное Ыь, тем точнее будет воспроизведена поверхность канала и, вместе с тем, тем дольше будет рассчитываться значение ЭПР канала заданной геометрии. Дальнейшее представление результатов для возможности сравнения и наглядности будет производиться на примере простейшего из рассматриваемых каналов - с круглыми апертурами.
Синтез 8-образного канала воздухозаборника
Одним из наиболее известных способов формирования вытянутых изогнутых объемных тел является задание закона изменения кривизны их строительной оси. В данной работе предлагается следующая формализация изменения координат (х', у\, г) точек центральной оси Б-образного канала воздухозаборника:
DOI: 10.24411/2410-9916-2019-10209
Systems of Control, Communication and Security
ISSN 2410-9916
= ■
A, 2
'cos(^-) - Г
L
y = У + (i -1)
_L_
N
(5)
z. =
A
_z
2
cos(^- ) L
где А х, А 2 - величины смещения центра входной или выходной апертуры от нулевого значения по осям X и 2 соответственно, у'0 - начальная точка цен-
тральной оси изогнутого тела, i
1...N
Очевидно, что в случае воздухозаборника круглой формы, то есть при R = Const, точки промежуточных сечений можно сформировать путем вычисления их координат согласно (1) при замене yi = yl. Однако, в случае задания входной апертуры иной формы, реализация адекватной поверхности объемного тела заданной кривизны невозможна. Один из предлагаемых вариантов правильного линейного формирования точек поверхности по промежуточным сечениям представлен выше при рассмотрении синтеза прямого канала воздухозаборника. На рис. 10 представлены некоторые результаты вычисления координат Na точек промежуточных NL сечений относительно строительной оси S-образной кривизны, вычисленной согласно (5) с соответствующими параметрами. Как видно из рисунка, формализация искривления строительной оси канала (5) неизбежно приводит к смещению входной или выходной апертуры. Поэтому на этапе задания начальных параметров синтеза целесообразно указывать диапазон максимальных смещений, допустимых конструктивно при компоновке воздушного судна. При этом интерес представляет указание смещений в форме соотношений с длиной волны X для последующего анализа полученных результатов.
Следующим этапом синтеза геометрии канала является параметрическое сужение или расширение NL промежуточных сечений вдоль заданной длины канала L. Простейшим решением этой задачи является изменение промежуточных сечений по заранее заданному закону. В данной работе в качестве примера рассматривается функция sin(y), где у - аргумент функции, изменяющийся в требуемом диапазоне, определяющем форму кривой. На самом деле функция может быть любой, обеспечивающей требуемое параметрическое искривление. Вычисление координат (хг", уг", z7") точек промежуточных сечений вдоль всей
длины канала с параметром искривления
A
R
производится согласно выражени-
ям:
х" = (, + х" ) + Ад (, + х" ) sin(^. )
Уi j Уi
(6)
z ..
( + Z" ) + Ar ( zj + Z" ) sin(ai )
<
-X
DOI: 10.24411/2410-9916-2019-10209
Системы управления,связи и безопасности №2. 2019
Systems of Control, Communication and Security ISSN 2410-9916
где I- - порядковый номер угла, задающего угловое положение точки сечения, 7 - - порядковый номер сечения вдоль длины канала Ь.
Рис. 10. Синтезированные точки поверхности Б-образного канала воздухозаборника относительно прямого, радиусом Я=1 м с параметрами искривления:
Лх =-1,0 м (а), Лх =+1,0 м (б), Лг =-1,0 м (в), Лг =+1,0 м (г), Лх=Л2 =-1,0 м (д), Лх=Л2 =+1,0 м (е)
На рис. 11 для примера представлены некоторые результаты вычисления координат точек искривленных сечений канала, сформированных при соответствующих параметрах и диапазоне аргумента у. Как видно на рис. 11, е, согласно предлагаемого метода параметрического синтеза канала можно сформировать произвольную полость с контролируемой с помощью параметров Ыа и ЫЬ точностью задания поверхности канала.
Дальнейшее преобразование полученного массива точек в фацетную модель Б-образного канала воздухозаборника производится с помощью одного из многочисленных методов триангуляции [11].
Пример результатов такого преобразования представлен на рис. 12, откуда видно, что при переходе к фацетной модели поверхности канала появляется возможность вычисления нормали к плоскости каждого фацета, используемой при расчетах ЭПР.
Следует отметить, что для оценки ЭПР нагруженного канала потребуется сформировать фацетную модель заглушки выходной апертуры канала либо «нагрузить» канал реальными фацетными моделями лопаток ВНА и компрес-
DOI: 10.24411/2410-9916-2019-10209
Системы управления,связи и безопасности №2. 2019
Systems of Control, Communication and Security ISSN 2410-9916
сора низкого давления двигателя. В рамках данной статьи рассматривается только «прямая» заглушка.
Рис. 11. Синтезированные точки поверхности Б-образного канала воздухозаборников относительно прямого, радиусом Я=\ м с параметрами искривления: АЛ =+0,2 м, 0...П (а), АЛ =-0,2 м (б), АЛ =+0,2 м, 0.. ,4п (в), А Л =+0,2 м, п/2.п (г), А х =+1,0 м, А г =-2,0 м, А Л =+0,5 м, 0.п (д), Ах =+1,5 м, Аг =+0,5 м, АЛ =+0,3 м, 0.4п, N¿=55 (е)
Рис. 12. Результат синтеза Б-образного канала воздухозаборника на этапах формирования: точек (а), поверхности (б) и треугольных фацетов (в) с параметрами Я=1 м, N¿=15, N«=101, Ах =-1,5 м, Аг =-1,5 м, Ак =+0,7 м, ^=0...л
Следует отметить, что для оценки ЭПР нагруженного канала потребуется сформировать фацетную модель заглушки выходной апертуры канала либо
DOI: 10.24411/2410-9916-2019-10209
Systems of Control, Communication and Security
ISSN 2410-9916
«нагрузить» канал реальными фацетными моделями лопаток ВНА и компрессора низкого давления двигателя. В рамках данной статьи рассматривается только «прямая» заглушка.
Таким образом, с учетом рассмотренных выше параметров синтеза поверхности S-образного канала воздухозаборника можно сформулировать критерий оптимальности параметрического синтеза:
S * = argmin <(max(<{e,p})' d Ддеф, Дх, Az, Дй)
S*eS (R,L,Na, NL )
< < <JTo
,(7)
A <Д <Д w
х min х х max
Д . <Д <Д
z min z z max
Д R min < Д R < Д R max
то есть из всех вариантов S возможных форм канала воздухозаборника с заранее заданными радиусом двигателя R, длиной канала L и детальностью (Na, Nl) фацетной модели требуется подобрать оптимальную геометрию с такими формой апертуры и параметрами искривления канала (d Ддеф, Дх, Д, Дд), удовлетворяющими диапазону конструктивно возможных смещений
(Дхmin'Дхmax'Дzmin'Дzmax' ДRmin. ДRmax) , при которых на всех ракурсах облучения
-90° < в < 90о и 0° < р < 360° обеспечивается минимальная ЭПР о, не превышающая регламентированного техническим заданием значения отз.
Методика оценки радиолокационной заметности канала воздухозаборника
Результаты анализа методов оценки ЭПР произвольных полостей, достаточно подробно представленные в [8], показали, что наиболее известным и широко используемым на сегодняшний день является Shooting and Bouncing Ray (SBR) [12] на основе лучевого подхода к моделированию процессов рассеяния ЭМВ и ее многократного переотражения. Определение координат точек пересечения лучей и треугольников (фацетов) реализовано с помощью алгоритма
[13], а вычисление значений ЭПР канала производится методом физической оптики с учетом электромагнитных свойств каждого фацета. В качестве примера на рис. 13 представлены результаты расчета ЭПР полостей, опубликованные в
[14]. Несмотря на заявляемую авторами работы высокую адекватность результатов, получаемых с помощью такого подхода, его недостатком являются огромные вычислительные затраты при плотности потока мощности излучения У20.
Кроме того, «сбор» многократно переотраженных лучей производится на плоскости входной апертуры или плоскости фронта ЭМВ без учета направления обратного вторичного излучения. В данной работе предложено использовать упрощенное решение [8], основанное на модифицированном методе геометрической оптики с учетом поляризационных свойств каждого фацета [5], верифицированное по экспериментальным данным, представленным в работе
DOI: 10.24411/2410-9916-2019-10209
Systems of Control, Communication and Security
ISSN 2410-9916
[15]. Расчет ЭПР предложено выполнять не по значениям напряженности электрического поля, а по значениям мощности отраженной и падающей ЭМВ с учетом направления обратного вторичного излучения, совпадающего с ракурсом облучения.
а)
б)
Рис. 13. Сравнительная оценка ЭПР круглых прямого и Б-образного каналов при ^=5Х/2, £=10Х, Дх=2,5Х на вертикальной (а) и горизонтальной (б) поляризациях
Следует обратить внимание, что оптимизацию геометрии канала можно производить по минимуму средней ЭПР по всем ракурсам, либо по медианным ЭПР в диапазоне углов (например, +5, +10 градусов), либо по удовлетворению всех измеренных значений ЭПР во всем диапазоне ракурсов требуемому значению ЭПР. Исходя из того, что обычно для каждого конкретного самолета в таблицах справочника указывается всего одно значение, можно сделать вывод о целесообразности поиска такой геометрии канала, которая во всех ракурсах обеспечивает ЭПР меньше заданной в техническом задании. При этом необходимо учесть, что полученное при оценке значение ЭПР канала с нагрузкой -лишь одна из составляющих суммарной ЭПР воздушного судна.
Результаты
Для примера в расчетах использована простейшая фацетная модель канала воздухозаборника, представленная на рис. 9, а. Размеры модели соответствуют соотношениям: £/2^=3,5, ^=33,3Х~1 м, £=233,33Х~7 м, где Л,=0,03 м. Диапазон конструктивно возможных смещений АхтЬ= 0, Ахтах = 1,8 м,
А
-0,3 м, А
0,3 м, А
Я тт
-0,3 м, А
Я тах
произвольно, в следующих пределах: 0 <Ах< 601
-101 < А < 101.
2
-101 <АС < 101
0,3 м в данном случае задан
(8)
На рис. 14 представлены диаграммы ЭПР канала с Б-образным искривлением геометрии канала по оси X. Параметр Ах изменялся в диапазоне от нуля до 100Х с шагом 5Х (0,15 м), при постоянных А2 и АЯ .
DOI: 10.24411/2410-9916-2019-10209
Системы управления,связи и безопасности №2. 2019
Systems of Control, Communication and Security ISSN 2410-9916
1.2
а) б)
Рис. 14. Диаграммы ЭПР S-образного канала воздухозаборника (а) и их максимальные значения ЭПР (б), полученные при в = 0°, 0° < ( < 180°,
0<А < 100А, А = 0, А, = 0
X 5 Z ' R
Как видно из приведенных на рис. 14, а диаграмм ЭПР, полученных при контролируемом с помощью параметра А искривлении канала только в одной плоскости, центральный и единственный пик (А X = 0 и ( = 90°) сначала раздваивается (АX = 20А« R /2,( = 75° и (2 = 90°) с одновременным снижением ЭПР, а при дальнейшем искривлении до А = 100А » 3R /2, центральный максимум значительно уменьшается до значений, близких к нулю, при этом первый максимум вместе со снижением ЭПР смещается в сторону по азимуту на 35° до (( = 55°. Анализ зависимости максимальных значений диаграмм ЭПР от А (рис. 14, б) позволяет сделать вывод о возможности снижения ЭПР канала воздухозаборника до 61,8 % (с 1,1 м2 до 0,42 м2 при Ax=60i) за счет S-образной формы канала при смещениях даже в одной плоскости.
Исходя из предлагаемого метода параметрического синтеза канала, по результатам промежуточной оценки ЭПР при всех азимутальных углах наблюдения 0°<(< 180° (рис. 14, а) необходимо выбрать только максимальные значения, а затем во всем диапазоне изменяемых параметров определить минимальные значения ЭПР из всех максимальных (рис. 14, б). Таким образом, можно обеспечить оптимальное на данном этапе значение параметра искривления: А X = 60 А = 1,8 м при начальном ограничении А X max = 60 А = 1,8 м.
Следует отметить, что результаты расчетов ЭПР S-образного канала при отрицательных АX, как и при отрицательных или положительных смещениях
А по оси Z в том же диапазоне, будут идентичными только в том случае, когда не учитывается поляризация ЭМВ. При необходимости исследования поляризационных характеристик канала требуются дополнительные вычисления для каждого типа смещения.
DOI: 10.24411/2410-9916-2019-10209
Системы управления,связи и безопасности №2. 2019
Systems of Control, Communication and Security ISSN 2410-9916
Вместе с тем, интерес представляет исследование зависимости значений ЭПР канала от степени его кривизны одновременно по двум осям X и 2. Как показано на рис. 15, а, искривление геометрии канала в двух плоскостях приводит к сопоставимым рис. 14, а результатам.
1.2
а) б)
Рис. 15. Диаграммы ЭПР Б-образного канала воздухозаборника (а) и их максимальные значения ЭПР (б), полученные при в = 0°, 0° < р < 180°,
-1001 <Дх< 0, -101<Д2< 101, Д = 0
1.2
а) б)
Рис. 16. Диаграммы ЭПР Б-образного канала воздухозаборника (а) и их максимальные значения ЭПР (б), полученные при в = 0°, 0° < р < 180°, Дх = 0,
Д = 0, -Ш<ДЛ< 101
По аналогии с результатом анализа максимальных значений ЭПР, представленных на рис. 14, б, по графику на рис. 15, б можно определить оптимальное значение параметров искривления строительной оси канала: Дх = -601 = -1,8 м и Д = 21 = 0,06 м, обеспечивающие минимальную радиолокационную заметность Б-образного канала во всем диапазоне ракурсов облучения.
На рис. 16, а представлены диаграммы ЭПР исходного канала с расширением/сужением поперечных сечений с помощью параметра Дк. Из рисунка
DOI: 10.24411/2410-9916-2019-10209
Systems of Control, Communication and Security
ISSN 2410-9916
видно, что только сужение прямого канала до ДЛ < -4Я позволяет уменьшить ЭПР до 0,76 м2. Из графика с максимальными значениями ЭПР канала на рис. 16, б можно сделать вывод, что оптимальное значение Дк составляет -10Х, при котором обеспечивается минимальная радиолокационная заметность.
Следует отметить, что, несмотря на сопоставимость результатов в первом и втором случаях (Д и Д), оптимальные значения параметров искривления
отличаются, что обуславливает необходимость учета всех возможных искажений канала. Такая необходимость исследования одновременного изменения всех параметров Дх ,Дz ,ДR во всех конструктивно допустимых диапазонах приводит к тому, что могут потребоваться значительные вычислительные ресурсы и существенные временные затраты. Так, оценка ЭПР в диапазоне углов в = 0° и 0°<^< 180° при 0<Д< 100А, -Ш<Д2<Ш и -Ш<ДЛ<Ш занимает
около 6 суток при проведении вычислений на ЭВМ с процессором Intel Core i5-4200M CPU 2,5 ГГц, ОЗУ 8 Гб под управлением 64-разрядной операционной системой Windows 7 без визуализации и оптимизации кода в среде MATLAB. Процесс оценки ЭПР в заданных диапазонах повторяется более 1,667 млн. раз. Для каждой реализации геометрии канала с конкретными значениями Д, Д, Дд определяется максимум ЭПР во всем диапазоне условий наблюдения (то есть при в = 0°, 0°<^< 180°). Среди всей выборки максимумов диаграмм ЭПР, приведенной на рис. 17, находится минимальное значение, которое и будет определять оптимальную геометрию канала воздухозаборника согласно предложенного критерия (7) при ограничительных условиях (8). Как видно из рисунка, множество максимальных значений диаграмм ЭПР канала при конкретных значениях Дх, Д обусловлено учетом третьего параметра искажения
канала ДR .
Рис. 17. Максимальные значения ЭПР S-образного канала воздухозаборника, полученные при в = 0°, 0°<^< 180°, 0 < Дх < 601, -10Л<Д2 < 101 и
-Ш<ДЛ < Ш
DOI: 10.24411/2410-9916-2019-10209
Системы управления,связи и безопасности №2. 2019
Systems of Control, Communication and Security ISSN 2410-9916
а) б) в)
Рис. 18. Проекции максимальных значений ЭПР Б-образного канала воздухозаборника на плоскости параметров <г0Дх (а), <г0Дг (б), ст0Дк (в)
Полученная выборка для наглядности может быть также представлена проекциями множества значений на оси этих параметров. Так, на рис. 18 представлены зависимости максимальных значений ЭПР от каждого параметра искажения.
Полученный массив значений максимумов ЭПР имеет размерность 4, поэтому для его представления в трехмерном виде произведена выборка минимумов максимальных значений ЭПР во всем диапазоне изменения ДЛ. На рис. 19 представлены минимумы ЭПР в координатах Дх, Д2, Дк в виде точек с числовыми значениями.
Рис. 19. Минимальные значения ЭПР Б-образного канала воздухозаборника, полученные при в = 0°, 0°<р< 180°, 0 <Д < 601, -101 <Д < 101 и
-101<ДЛ< 101
Анализ приведенных на рисунке значений позволяет определить координаты точки с минимальным значением ЭПР, то есть определить параметры геометрии канала воздухозаборника, оптимальной по критерию минимальной радиолокационной заметности. Результаты поиска экстремума по данным, по-
DOI: 10.24411/2410-9916-2019-10209
Systems of Control, Communication and Security
ISSN 2410-9916
казанным на рис. 19, показали, что минимум ЭПР канала составляет 0,11257 м2, что соответствует параметрам искривления А х = 601 = 1,8 м, А г = -51 = -0,15 м и АЛ= -91 = -0,27 м.
При данных параметрах оптимальная геометрия канала имеет вид, представленный на рис. 20.
а) б)
Рис. 20. Оптимальная по минимуму радиолокационной заметности геометрия канала воздухозаборника на этапе синтеза точек поверхности (а) и после триангуляции с вычисленными нормалями каждого фацета (б)
Азимут, град.
а) Аф = 1°, 10Х
Азимут,
б) аф = 0,1°, m
Азимут
в) Аф = 1°, 11
Азимут,
г) Аф = 1°, М0
Рис. 21. Диаграммы ЭПР канала, рассчитанные с шагом Аф = 1° (а), Аф = 0,1° (б), с плотностью потока мощности излучения 1Х (в) и М0 (г)
Оценка ЭПР канала воздухозаборника со сниженной радиолокационной заметностью производилась вышеуказанным модифицированным методом геометрической оптики [5] с различным шагом дискретизации углов наблюдения и плотностью потока мощности излучения (рис. 21). Как видно из приведенных графиков, уменьшение шага углов наблюдения Аф при остальных равных значениях параметров приводит к незначительному уточнению значений углов ра-
DOI: 10.24411/2410-9916-2019-10209
Системы управления,связи и безопасности №2. 2019
Systems of Control, Communication and Security ISSN 2410-9916
курса облучения и достаточно существенному изменению значений ЭПР канала (рис. 21, а, б). Уплотнение потока мощности излучения с 10Х до Ш0 (рис. 21, а, в, г) ведет к значительному изменению формы диаграммы ЭПР, а также некоторому смещению значений на 3° в области р = 20°. Поэтому для дальнейшего практического применения предлагаемого метода параметрического синтеза оптимальной геометрии требуется верификация полученных результатов экспериментальными данными измерений с последующей калибровкой значений плотности потока излучения. Полученные зависимости и сформулированные выводы сопоставимы с результатами и выводами, опубликованными в [1-3].
Важно также отметить, что перед реализацией и макетированием реальных образцов каналов с рассчитанной оптимальной геометрией в обязательном порядке требуются расчеты по газодинамической устойчивости. Для этого целесообразно использовать двухкритериальный метод параметрического синтеза оптимальной геометрии канала воздухозаборника.
Выводы
Таким образом, в данной работе представлен метод параметрического синтеза формы канала воздухозаборника, оптимальной по критерию минимальной ЭПР во всех ракурсах наблюдения, обеспечивающей снижение его радиолокационной заметности до 90%. Предложенная формализация итеративного изменения параметров искривления геометрии канала позволяет с необходимой точностью воспроизвести форму апертур и поверхность канала в заданных при проектировании конструктивных ограничениях.
Следует отметить, что используемый при расчетах ЭПР модифицированный метод геометрической оптики с учетом поляризации - далеко не самый точный, но самый быстрый, поэтому для получения более достоверных результатов целесообразно применять один из приближенных (асимптотических) методов физической оптики. В случае наличия в канале реальной нагрузки в виде лопаток входного направляющего аппарата и компрессора низкого давления требуется дополнительное использование методов геометрической или физической теории дифракции. Их программные реализации в рамках предложенного подхода и рассмотренного метода не составляют особой сложности.
Как следует из полученных в работе результатов, для расчета ЭПР в X-диапазоне при синтезе оптимальной геометрии канала достаточной плотностью потока мощности излучения является значение, близкое к длине волны. Дальнейшее уплотнение потока излучения до Ш0 и Х/20 требуется лишь для уточнения результатов.
В данной статье, вследствие ограничения по объему материала, не представлены результаты исследований ЭПР канала воздухозаборника по всему диапазону углов места, с различными формами входных апертур, реальными нагрузками канала в виде лопаток входного направляющего аппарата и компрессора низкого давления, с применением РПМ, в различных частотных диапазонах и типах поляризации ЭМВ, а также по газодинамической устойчивости.
DOI: 10.24411/2410-9916-2019-10209
Системы управления,связи и безопасности №2. 2019
Systems of Control, Communication and Security ISSN 2410-9916
Литература
1. Knott E. F., Shaeffer J. F., Tuley M. T. Radar cross section. Scitech Publishing, 2004. 638 p.
2. Сухаревский О. И., Василец В. А., Кукобко С. В. Рассеяние электромагнитных волн воздушными и наземными радиолокационными объектами: монография. Харьков: Харьковский национальный университет Воздушных Сил имени Ивана Кожедуба, 2009. 468 с.
3. Вождаев В. В., Теперин Л. Л. Характеристики радиолокационной заметности летательных аппаратов. М.: Физматлит, 2018. 376 с.
4. Кузнецов В. А. Геометрические эффекты поляриметрического рассеяния электромагнитных волн элементами турбореактивного двигателя воздушного судна со сниженной радиолокационной заметностью // Академические Жуковские чтения: сборник статей по материалам Международной научно-практической конференции. Воронеж: Военно-воздушная академия, 2018. С. 138-142.
5. Кузнецов В. А., Амбросов Д. В. Динамическая модель истребителя F-22 Raptor // Радиолокация, навигация, связь: сборник трудов Международной научно-технической конференции. Воронеж: Воронежский государственный университет, 2019. Т. 3. С. 126-137.
6. Лагарьков А. Н., Погосян М. А. Фундаментальные и прикладные проблемы стелс-технологий // Вестник РАН. 2003. Т. 73. № 9. С. 779-787.
7. Кузнецов В. А. Структура и свойства наземных объектов на изображениях РСА в задачах распознавания // Телекоммуникации. 2012. № 10. С. 31-38.
8. Кузнецов В. А., Гончаров С. А. Структурно-параметрический синтез малогабаритной радиолокационной станции с синтезированной апертурой беспилотного летательного аппарата ближнего действия // Системы управления, связи и безопасности. 2017. № 9. С. 28-72.
9. Белов А. М. Параметрические уравнения четырехугольников // Информационный сайт А.М. Белова [Электронный ресурс]. 01.08.2009. - URL: www.stob2.narod.ru/30s.htm (дата обращения 21.04.2019).
10. Agoston M. K. Computer graphics and geometric modeling: implementation and algorithms. Springer, 2005. 920 p.
11. Скворцов А. В., Мирза Н. С. Алгоритмы построения и анализа триангуляции. Томск: Томский государственный университет, 2006. 168 с.
12. Ling H., Chou R., Lee S. W. Shooting and bouncing rays: calculating the RCS of an arbitrary shaped cavity // IEEE Transactions of Antennas and propagation. 1989. Vol. 37. No. 2. P. 194-205.
13. Moller T., Trumbore B. Fast, minimum storage ray/triangle intersection // Journal of Graphics Tools. 1997. Vol. 2(1). P. 21-28.
14. Chou R., Ling H., Lee S. W. Reduction of the radar cross section of arbitrarily shaped cavity structures // Electromagnetics Laboratory report, 1987. No. 87-6. 123 p.
15. Справочник по радиолокации / под ред. М. И. Скольника, пер. с англ. под общей ред. В. С. Вербы. В 2 книгах. М.: Техносфера, 2004. 672 с.
DOI: 10.24411/2410-9916-2019-10209
Системы управления,связи и безопасности №2. 2019
Systems of Control, Communication and Security ISSN 2410-9916
References
1. Knott E. F., Shaeffer J. F., Tuley M. T. Radar cross section. Scitech Publishing, 2004. 638 p.
2. Sukharevsky O. I., Vasilets V. A., Kukobko S. V., Nechitaylo S. V., Sazonov A. Z. Rasseianye Elektromahnytnykh Voln Vozdushnymy y Nazemnymy Radyolokatsyonnymy Obektamy [The Electromagnetic Wave Scattering By Aerial And Ground Radar Objects]. Kharkiv, Khar'kovskii natsional'nyi universitet Vozdushnykh Sil imeni Ivana Kozheduba, 2009. 468 p. (in Russian).
3. Vozhdaev V. V., Teperin L. L. Kharakteristiki Radiolokatsionnoi Zametnos-ti Letatel'nykh Apparatov [Characteristics of Radar Aircraft Visibility]. Moscow, Fizmatlit, 2018. 376 p. (in Russian).
4. Kuznetsov V. A. Geometricheskie effekty poliarimetricheskogo rasseianiia elektromagnitnykh voln elementami turboreaktivnogo dvigatelia vozdushnogo sudna so snizhennoi radiolokatsionnoi zametnost'iu. Akademicheskie Zhukovskie chteniia: sbornik statei po materialam Mezhdunarodnoi nauchno-prakticheskoi konferentsii [Proceedings of the International Scientific and Practical Conference]. Voronezh, Air Force Academy, 2018, pp. 138-142 (in Russian).
5. Kuznetsov V. A., Ambrosov D. V. Dinamicheskaia Model' Istrebitelia F-22 Raptor. Radiolokatsiia, navigatsiia, sviaz': sbornik trudov Mezhdunarodnoi nauchno-tekhnicheskoi konferentsii [Radiolocation, navigation, communication. Proceedings of the International Scientific and Technical Conference]. Voronezh, Voronezhskii gosudarstvennyi universitet, 2019, no. 3, pp. 126-137 (in Russian).
6. Lagar'kov A. N., Pogosian M. A. Fundamental'nye i Prikladnye Problemy Stels-tekhnologii [Fundamental and Applied Problems of Stealth Technology]. VestnikRossiiskoi akademii nauk, 2003, Vol. 73, no 9, pp. 779-787 (in Russian).
7. Kuznetsov V. A. Struktura i Svoistva Nazemnykh Ob"ektov na Izobrazheniiakh RSA v Vadachakh Raspoznavaniia [Structure and Properties of Ground Objects on Radar Images in Recognition Tasks]. Telekommunikatsii, 2012, no. 10, pp. 31-38 (in Russian).
8. Kuznetsov V. A., Goncharov S. A. Structural-Parametric Synthesis of Short-Range Unmanned Aerial Vehicle Small-Sized SAR. Systems of Control, Communication and Security, 2017, no. 9, pp. 28-72 (in Russian).
9. Belov A. M. Parametricheskie uravneniia chetyrekhugol'nikov. Informatsionnyi Sait A.M. Belova, 2009. Avaiable at: http://www.stob2.narod.ru/30s.htm (accessed 21 April 2019) (in Russian).
10. Agoston M. K. Computer Graphics and Geometric Modeling: Implementation and Algorithms. Springer, 2005. 920 p.
11. Skvortsov A. V., Mirza N. S. Algoritmy Postroeniia i Analiza Trianguliatsii [Algorithms for Constructing and Analyzing Triangulation]. Tomsk, Tomskii gosudarstvennyi universitet, 2006. 168 p. (in Russian).
12. Ling H., Chou R., Lee S. W. Shooting and Bouncing Rays: Calculating the RCS of an Arbitrary Shaped Cavity. IEEE Transactions of Antennas and propagation, 1989, vol. 37, no. 2, pp. 194-205.
13. Moller T., Trumbore B. Fast, Minimum Storage Ray/Triangle Intersection. Journal of Graphics Tools, 1997, vol. 2(1), pp. 21-28.
DOI: 10.24411/2410-9916-2019-10209
Системы управления,связи и безопасности №2. 2019
Systems of Control, Communication and Security ISSN 2410-9916
14. Chou R., Ling H., Lee S. W. Reduction of the radar cross section of arbitrarily shaped cavity structures. Electromagnetics Laboratory report. 1987, no. 87-6, 123 p.
15. Skolnik M.I. Radar Handbook. McGraw-Hill Education, 2008. 1328 p.
Статья поступила 31 мая 2019 г.
Информация об авторe
Кузнецов Виктор Андреевич - кандидат технических наук. Старший преподаватель кафедры эксплуатации бортового авиационного радиоэлектронного оборудования. Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия имени профессора Н.Е. Жуковского и Ю.А. Гагарина» (г. Воронеж). Область научных интересов: системный анализ, радиолокация, распознавание образов, фракталы. E-mail: kuzzviktor@mail.ru
Адрес: 394064, г. Воронеж, ул. Старых Большевиков, д. 54А.
The radar aircraft visibility reducing method by using air intake duct optimal shape parametric synthesis
V. A. Kuznetsov
Purpose: the problem of reducing the radar visibility of aircraft is relevant throughout the world, and its solution consists of a set of preventive actions to reduce the radar cross section and, moreover, to reduce the informativity of the demasking signs. The air intake duct makes a significant contribution to the total effective radar cross section area in the forward hemisphere. However, the secondary Doppler modulation effect presence during the air intake electromagnetic waves reflection process necessitates its radar visibility reduction. The aim of this paper is to reduce the radar aircraft visibility by parametric synthesis of the S-shaped air intake duct. The shape of the air intake duct is optimal by the criterion of the minimum of the effective radar cross section area in a given range of structurally possible changes. Methods: when solving the problem of reducing radar visibility, the parametric synthesis of the channel shape is used. The geometrical optics method is modified to estimate the effective radar cross section area to account for the polarization of the incident electromagnetic wave. The mathematical modeling is widely used. Novelty: the element of novelty is the parametric synthesis method of the effective radar cross section area of the air intake duct shape. The shape in that case is optimal by the criterion of the minimum under all possible conditions for observing. Also the element of novelty is the geometrical optics method modification which allows to take into account the polarization properties of each facet channel surface. Results: The proposed method allows to reduce the effective radar cross section area of the aircraft air intake duct more than 90% in comparison with the direct channel shape, over the entire range of irradiation angles. It also allows to synthesize the shape of the S-shaped air intake duct with required accuracy when an input aperture shape and a structurally possible channel geometry curvature range are parametrically defined. Practical relevance: the presented solution of the radar visibility reducing problem can be used when the estimating process of the effective radar cross section area taking into account the actual air intake duct load. That load has the form of blades of the inlet guide vanes and compressors of an aircraft engine. The method also can help to develop models of air intakes for prospective aircraft.
Key words: radar visibility, radar cross section, geometrical optics, faceted model, air intake duct.
Information about Author
Victor Andreevich Kuznetsov - Ph.D. of Engineering Sciences. Senior lecturer. Chair of aircraft radioelectronic equipment exploitation. Military Educational and Scientific Centre of the Air Force N.E. Zhukovsky and Y.A. Gagarin Air Force Academy» (Voronezh). Field of research: system analysis, radiolocation, pattern recognition, fractals. E-mail: kuzzviktor@mail.ru
Address: Russia, 394064, Voronezh, Starykh Bolshevikov, 54A.
DOI: 10.24411/2410-9916-2019-10209
