CC BY
121
УДК 623.466.55
ОЦЕНКА РАДИОЛОКАЦИОННЫХ ХАРАКТЕРИСТИК МИШЕНИ -«БОЛЬШОЙ КОРАБЕЛЬНЫЙ ЩИТ»
М.В. Зонов
EVALUATION OF RADAR CHARACTERISTICS OF TARGET -“GREAT SHIP SHIELD”
M.V. Zonov
Предложена модель перемещения кажущегося центра отражения многоточечного надводного объекта на примере мишени - «большой корабельный щит пр. 436-бис». Результатом работы модели является процесс изменения мощности отраженного от мишени сигнала во времени и процесс флюктуации кажущегося центра отражения в пределах продольной оси мишени. Справедливость модели подтверждена путем сравнения результатов работы модели и результатов натурной работы с ракетой по этой мишени (участок прямого подлета ракеты к мишени с 7 км). Для обеспечения корректного сопоставления результатов модели и натурной работы модель отражения дополнена упрощенной моделью угловой автоматики радиолокационной системы ракеты.
Ключевые слова: радиолокация, характеристики мишени.
The paper presents a model of the apparent movement of the center of reflection of the multi-point surface of the on-water object as an example of the target - the «great ship shield». The result of the modeling is the process of changing of the power reflected from the target signal in time and process fluctuation of the apparent center of reflection within the longitudinal axis of the target. The validity of the model is confirmed by comparing the results of the model and results of field work with a missile on the target.
Keywords: radar, characteristics of target.
При анализе результатов испытаний ракет возникла необходимость в оценке радиолокационных характеристик мишени - «большого корабельного щита пр. 436-бис» (БКЩ) на конечном участке полета ракеты, а именно измерение эффективной отражающей поверхности (ЭПР) и получение процесса флюктуации кажущегося центра отражения (КЦО) в пределах продольной оси мишени. Аналогичные измерения параметров мишени были проведены более 30 лет назад в воинской части, эксплуатирующей данную мишень, но в настоящее время воинская часть этими результатами не располагает.
Проведение дополнительных облетов мишени самолетом-лабораторией является весьма дорогостоящим мероприятием. Поэтому для оценки радиолокационных характеристик мишени целесообразнее провести моделирование отражения радиоволн от мишени в сантиметровом диапазоне длин волн.
Оценка радиолокационных характеристик мишени приведена в научно-техническом отчете [1].
Зонов Михаил Владимирович - инженер-конструктор ОАО «ОКБ «Новатор», г. Екатеринбург; [email protected]
Но в данном отчете результаты моделирования приведены только для участков полета ракеты в несколько десятков длин волн, чего явно недостаточно для решения поставленной задачи. Получение результатов работы этой модели для всего конечного участка полета ракеты затруднительно, так как отчет не содержит подробного описания модели, а также не приведена и сама модель.
В связи с этим предложена простая и удобная (в плане вычислительных затрат) модель перемещения КЦО мишени, которая в отличие от указанного моделирования [1] не учитывает сферичности фронта отраженной волны в ближней зоне наблюдения мишени.
Корректность предложенной модели подтверждена путем сравнения результатов расчета и результатов натурной работы с реальной ракетой по БКЩ. Сравнение результатов расчета и натурной работы выполнено по имеющимся в телеметрической информации параметрам юа и ра - угловой скорости и угловому положению по азимуту строительной оси антенны радиолокационной сис-
Zonov Mihail Vladimirovich - design engineer of OSC “EDB “Novator”, Yekaterinburg; [email protected]
темы наведения (РСН) ракеты. Для преобразования расчетного положения КЦО мишени в юа и ва использована упрощенная модель системы угловой автоматики (СУА) РСН.
Согласно данным в [2] мишень БКЩ представляет собой надводную плавающую платформу с расположенными на ней в особом порядке 11 отражающими элементами - уголковыми отражателями (УО) с разными высотами расположения и разными значениями ЭПР.
Как сказано в [3], такой многоточечный объект формирует суммарный кажущийся центр отражения. По мере подлета к мишени изменяется характер ее отражающей способности - ЭПР и смещается положение КЦО. Временные флюктуации угла прихода отраженного от БКЩ сигнала называют угловым шумом или угловым мерцанием объекта. Это явление приводит к угловой ошибке при сопровождении радиолокационной системой наведения ракеты.
Исходными данными для расчета являются телеметрические параметры, полученные с натурной работы ракетой по мишени БКЩ: координаты ракеты - Хр, Zp, Yp, проекции скорости ракеты на оси х и г - Vx, Vz и координаты цели - Хц, Zц. РСН ракеты работает в многочастотном режиме автосопровождения (цикличный перебор из к частот (массив длин волн ^[к]) с фиксированным шагом) при обработке 5 подпачек из п импульсов. Темп считывания исходных данных и шаг расчета соответствуют периоду повторения импульсов в РСН.
Также в качестве исходных данных использованы координаты расположения УО на мишени (массивы трехмерных координат Х[г], Y[i], Z[i]) и ЭПР каждого из УО (массив а[/']), где i = 0.. .11.
На рис. 1 показана структурная схема программы расчета углового положения, скорости перемещения и ЭПР КЦО мишени.
Принимается, что БКЩ расположен перпендикулярно к направлению подлета ракеты, вследствие чего ось ХБКщ сонаправлена с осью Zр (рис. 2).
р.і/і
Модуль расчета углового положения УО —► —► —► —► / \ Модуль / \ ДНА РСН
1' ’' ’' ’ 1'' ’' 11 ’' ’'
Исх. данные:
х[ц, ш т, Модель
ОЙ БКЩ
Н
Исх. данные: Хр, ^р, Гр, Ух, V
Хц, ^ц, Ь|к|
Ша/Т
ва/і
Рис. 1. Структурная схема программы расчета характеристик БКЩ
-БКЩ
ПКР ;
Рис. 2. Ситуационная модель для расчета углового положения КЦО мишени
Координаты центра БКЩ совмещены с координатами цели Хц, 2ц.
Исходные данные Х[і], У[і], 2[і] поступают в модуль, выполняющий расчет углового положения уголкового отражателя БКЩ. Угловое положение і-го УО может быть определено как ' X [і]'
¥ уО И = -a tan
D
(1)
где Х[і] - координата і-го уголкового отражателя по продольной оси БКЩ; Б - текущая дальность от ракеты до центра БКЩ (Хц, 1Ц).
Модуль ДНА РСН содержит аппроксимирующую функцию вида
9а (¥) =
sin(uy)
u¥
(2)
где и = 20, V = 2 - параметры аппроксимации ДНА; у (в рад, см. рис. 2) является аргументом функции диаграммы направленности антенны, соответствующим угловому положению /'-го уголкового отражателя БКЩ из (1) и текущему угловому положению строительной оси антенны ра[/] для у-й подпачки (У = 1.5), и может быть определен как
¥ = ¥УО[/] -Ра[У].
В соответствии с (2) ширина ДНА по уровню половинной мощности составляет 8°.
С помощью функции, аппроксимирующей диаграмму направленности антенны РСН, в модуле ДНА РСН производится корректировка значения ЭПР каждого уголкового отражателя БКЩ.
Корректировка ЭПР /-го уголкового отражателя (о[/]) осуществляется следующим образом:
о[/]* = о[/] 0а(у). (3)
Расчет [3-5] углового положения кажущегося центра отражения БКЩ (ХКцО, ГКцО, ZКцО) осуществляется для текущего положения ракеты при излучении РСН одного импульса с длиной волны Х[к] (к - номер рабочей частоты РСН) путем нахождения углового положения КЦО, образованного первыми двумя отражателями (см. рис. 2, УО1 и УО2), а затем последовательного нахождения углового положения КЦО, образованного следующим уголковым отражателем (УО/) и кажущимся центром отражения, определенным ранее. Процесс продолжается до достижения последнего УОц.
Соотношения для расчета КЦО для БКЩ в соответствии с ситуационной моделью (см. рис. 2) приведены ниже:
1) дальность от ракеты до первого уголкового отражателя (УО1)
ё = ^(( - 7[1] - Хр )2 +((] - 7р )2 +( + X [1] - 7р )2;
2) дальность от ракеты до второго уголкового отражателя (УО2)
¿2 ^(Хц -Щ]-Хр)2 +(2]-Гр)2 +(ц + Х[2]-7р)2;
3) разность дальностей от ПКР до первых двух отражателей
Аё = ¿2 — ¿1;
4) набег фазы сигнала РСН (в рад) на интервале Аё
4пАё
ф =-------;
т Х[к]
5) коэффициент т отношения откорректированных по (3) ЭПР УО1 и УО2
о[1]
о[2]
6) суммарная ЭПР УО1 и УО2 (ЭПР кажущегося центра отражения от УО1 и УО2)
JКЦО
= о[1] +о[2] + 2-
V°[1]* °[2]*
cos ф;
7) расстояние между УО1 и УО2
DI2 = >/Х[1] -X[2]) +(1] - Y[2])2 + ([1] - Z[2])2
8) угол 9 (в рад)
' Ad
9 = arcsin
D
12
9) величина смещения КЦО («центр тяжести»)
AL = --
D12 a
2(cos2 c + b sin2 c) ’
где коэффициенты a, b, c определяются следующим образом:
п2
1 -m q u
a =-------cos 9, b =
1 + m
1 - m 1 + m
c = -
2nAd
10) координаты КЦО (ХКцО J2), образованного УО1 и УО2,
X
КЦО 12
X + X 2 +ALX1 - x2
D
12
После обработки последнего излученного РСН импульса у-й подпачки производится вычисление усредненного XКцО за подпачку по формуле
X
КЦО
X X
n=1
КЦО
64
где ХКцО - положение кажущегося центра отражения БКЩ для п-го излученного импульса РСН.
Модель системы угловой автоматики РСН вычисляет угловое положение антенны РСН в азимутальной плоскости, соответствующее отработанному усредненному за подпачку положению КЦО (ХКцО ).
На рис. 3 приведена упрощенная структурная схема модели СУА РСН по азимуту.
Моделирование звеньев произведено путем численного интегрирования с шагом (А/), равным периоду повторения импульсов РСН, умноженному на число импульсов в подпачке. Коэффициент усиления к = 6, постоянная времени Т = 0,5 с.
Угол линии визирования (фЛВ) вычисляется по формуле
ф ЛВ =-a sin
X
КЦО
D
Рис. 4. Изменение углового положения и угловой скорости перемещения КЦО в зависимости от дальности ракета-БКЩ
Переменной Б в данном случае обозначена дальность от ракеты до КЦО БКЩ:
Б-і/(Хц -^кцэ Хр) + (^Кцэ ) +(^ + хкцо ^р
Угловая скорость антенны РСН юа для у-й подпачки (юа[/]) вычислена как
' кФлв - кРа [У - 1]
«а [У] -Ша [І - 1 + А/
Т
'-«а [ І - 1]
ЭПР БКЩ, кв.м
І І І I
ііі
17500
Iі I
г
1' '11 Т [Г
.її і і Лііі
|І I»'
ПГ1 Гг 1 р I
Т' I
15000
12500
10000
6750 6250 5750 5250 4750 4250 3750 3250 2750 2250 1750 1250 750 250
Дальность, м
Рис. 5. Зависимость расчетной ЭПР кажущегося центра отражения БКЩ от дальности ракета-БКЩ
22500
20000
7500
5000
2500
0
50
м
40
30
- -30
1 -40
Рис. 6. Зависимость смещения КЦО вдоль продольной оси БКЩ ( Л Кц() ) от дальности ракета-БКЩ
Угловое положение антенны РСН по азимуту Ра (Ра[/]) определяется следующим образом:
Ра [ І ] -Ра [ І - 1] + А/«а [ І ].
При расчете Ра[/] и Ша[/] для первой подпачки (/ = 1) принимается ша/-1] = Ра[/-1] = 0, т. е. в начале расчета строительная ось антенны РСН направлена на центр БКЩ (нулевые начальные условия).
Результаты работы программы расчета углового положения, скорости перемещения КЦО и ЭПР БКЩ приведены на рис. 4, 5. На рис. 4 показано изменение углового положения и изменение угловой скорости перемещения КЦО в азимутальной плоскости в зависимости от дальности ракета-БКЩ (расчетные Ша[/] и Ра[/] обозначены как ша_Расч и Ра_Расч). На этом же рисунке для сравнения и подтверждения результатов расчета приведены эти же параметры ша и Ра, взятые из телеметриче-
ской информации (ша_НР и Ра_НР) с натурной работы ракеты (прямолинейный подлет ракеты перпендикулярно продольной оси БКЩ на конечном участке траектории).
На рис. 5 приведена зависимость расчетной ЭПР кажущегося центра отражения БКЩ в зависимости от дальности от ракеты до БКЩ.
На рис. 6 приведена зависимость смещения КЦО вдоль продольной оси БКЩ (ХКцО ) от дальности от ракеты до БКЩ.
Выводы
Проведена оценка характеристик для мишени типа БКЩ пр. 436-бис при прямом подлете ракеты с дальности 7 км, перпендикулярно продольной оси мишени. Амплитуда флюктуаций КЦО вдоль продольной оси БКЩ в среднем составила ±15 м относительно центра БКЩ. Среднее значение ЭПР
БКЩ составило порядка 15 000 м2. Сравнение графиков на рис. 4 указывает на схожесть полученных результатов и результатов натурной работы. Например, значение дисперсии для юа_Расч составило 0,14 (°/с)2, для юа_НР - 0,16 (°/с)2. Математическое ожидание для юа_Расч составило -0,05 °/с, для юа_НР составило -0,08 °/с.
Разница в значениях дисперсии и математических ожиданиях связана с неточностью принятых моделей, возможным несовпадением фаз сигналов в реальной РСН с расчетным набегом фазы в модели БКЩ и отсутствием учета отражений от морской поверхности. Также дополнительный вклад в разницу математических ожиданий вносит имеющаяся статическая ошибка СУА реальной РСН. В целом принятая модель позволяет достаточно достоверно оценить радиолокационные характеристики БКЩ.
Намечен путь по улучшению разработанной программы расчета ЭПР, углового положения и скорости перемещения КЦО мишени, заключающийся в разработке модели системы управления ракеты по курсу (на рис. 1 - блок, обозначенный пунктирной линией). Это позволит выбирать раз-
личные пути подлета ракеты к БКЩ, а также корректировать курс полета ракеты в зависимости от изменения юа, выдаваемого РСН в бортовую аппаратуру ракеты.
Литература
1. Оценка радиолокационных характеристик БКЩ: науч.-техн. отчет. 1-4 этапы НИР / УГТУ-УПИ. Екатеринбург, 2006-2007. Техн. архив ОКБ «Новатор», ТДК № 42947, 43123, 43300, 43391.
2. Инструкция по оборудованию и использованию мишеней при стрельбе ракетами с радиолокационными или тепловыми устройствами самонаведения. - М.: Воениздат МО СССР, 1965.
3. Тверской, Г.Н. Имитаторы эхо-сигналов судовых радиолокационных станций / Г.Н. Тверской, Г.К. Терентьев, И.П. Харченко. - Л.: Судостроение, 1973. - 224 с.
4. Корн, Г. Справочник по математике / Г. Корн, Т. Корн. - М.: Наука, 1974. - 832 с.
5. Моделирование в радиолокации / А.И. Леонов, В.Н. Васенев, Ю.И. Гайдуков и др.; под ред. А.И. Леонова. - М.: Сов. радио, 1979. - 264 с.
Поступила в редакцию 12 октября 2010 г.
