Доклады БГУИР_DoкLADY BGUIR
2019, № 1 (119) 2019, N0.1 (119)
УДК 621.396.96
СПОСОБ ИЗМЕРЕНИЯ ЧАСТОТ СЛЕДОВАНИЯ ЛОПАСТЕЙ ВИНТОВ ЛЕТАТЕЛЬНОГО АППАРАТА НА ОСНОВЕ СВЕРТКИ СПЕКТРА «ВТОРИЧНОЙ» МОДУЛЯЦИИ
С Р. ГЕЙСТЕР1, Т Т. НГУЕН2
1 Научно-исследовательский институт Вооруженных Сил Республики Беларусь, Республика Беларусь 2Белорусский государственный университет информатики и радиоэлектроники, Республика Беларусь
Поступила в редакцию 7 декабря 2018
Аннотация. Предлагается новый способ измерения частот следования лопастей винтов летательного аппарата (в частности, несущего (НВ) и рулевого (РВ) винтов вертолета) для цепей адаптации радиолокационного датчика (РЛД) с обращенным синтезом апертуры антенны. В основе способа лежит свертка спектра «вторичной» модуляции сигнала с одновременной компенсацией доплеровской частоты отраженного сигнала. Это позволяет в типовых для РЛД условиях быстрого изменения доплеровской частоты обеспечить длительное когерентное накопление и, соответственно, высокую точность измерения частот вращения винтов.
Ключевые слова: частота вращения, частота следования лопастей, несущий винт, рулевой винт, вертолет.
Abstract. A novel measurement method for repetition frequencies of the blades of an aerial vehicle (in part, helicopter's main and tail rotor) for adaptation circuits of the radar sensor with an inverse synthesis of antenna aperture is proposed. The method based on the convolution of the spectrum of the «secondary» modulation of the signal with simultaneous compensation of the Doppler frequency of the reflected signal. This allows in typical conditions for the radar sensor rapid change of Doppler frequency to provide a long coherent accumulation and, accordingly, high accuracy of measurement of the rotational frequencies of the rotors.
Keywords: rotational frequency, repetition frequency of the blades, main rotor, tail rotor, helicopter.
Doklady BGUIR. 2019, Vol. 119, No. 1, pp. 68-74
Measurement method for repetition frequencies
of the blades of an aerial vehicle on the base of the convolution
of the «secondary» modulation spectrum
S.R. Heister, T.T. Nguyen
Постановка задачи
Ошибка = / NL измерения частоты вращения ¥у винта с NL лопастями
определяет качество формирования радиолокационного изображения (РЛИ) винта и напрямую связана с ошибкой 5^. измерения частоты следования импульсов «вторичной» модуляции
¥Ьу = NLFV (или частоты следования лопастей) для приближающихся (удаляющихся) лопастей в отраженном сигнале (ОС) винта. Для определения частоты ^ необходимо вначале точно измерить частоту следования приближающихся (удаляющихся) лопастей, а затем
разделить ее на предполагаемое число лопастей NL.
Обычно при распознавании самолетов частоту следования импульсов «вторичной» модуляции FLV, формируемых лопастями винта (лопатками компрессора или турбины), определяют по спектру принятого сигнала [1]. Это возможно при двух условиях: первое -большое количество лопастей (лопаток) и большая частота вращения вала компрессора (турбины) или тянущего винта; второе - малое изменение радиальной скорости летательного аппарата (ЛА) за время измерения, что напрямую связано с изменением ракурса наблюдения.
Оба условия нереализуемы при наблюдении ЛА в РЛД и, в частности, вертолетов на предельно малых высотах: частоты вращения НВ вертолета находятся в пределах FV = 3...5 Гц при числе лопастей NL = 2...8; дальность действия РЛД RRLD небольшая (500...1000 м). Поэтому ракурс наблюдения в ходе полета ЛА относительно РЛД, ведущего мониторинг в секторе АРГг^ > (20°...40°) под углом установки около 45°, будет быстро изменяться [2].
Изменение радиальной скорости ЛА относительно РЛД в ходе когерентного накопления ОС в узкополосных фильтрах быстрого преобразования Фурье (БПФ) вызывает «размывание» спектра «вторичной» модуляции и приводит или к неприемлемым ошибкам, или к принципиальной невозможности измерения частоты FLV . Именно это является главной причиной невозможности применения традиционного способа измерения частоты вращения винта. Далее рассматривается новый способ измерения частоты следования импульсов «вторичной» модуляции FLV (частоты следования приближающихся (удаляющихся) лопастей винта), устраняющий влияние быстрого изменения радиальной скорости, обеспечивающий свертку спектра «вторичной» модуляции (ССВМ) ОС винта и когерентное накопление свертки.
Физические основы когерентного накопления свертки спектральных составляющих «вторичной» модуляции отраженного сигнала винта
При монохроматическом зондирующем сигнале (МХЗС) отраженный от винта сигнал ирь V (^ включает [3, 4] последовательности амплитудно-фазомодулированных импульсов
«вторичной» модуляции от приближающихся (ПЛ) ир рь V (^ и удаляющихся (УЛ) ир V ^) лопастей, следующих с одинаковым периодом Тр = 1 / FLV . Спектр комплексной огибающей ОС винта формируется относительно доплеровской частоты FdopK корпуса ЛА и при МХЗС
(см. рис. 1) имеет огибающую, определяемую спектром одиночного импульса «вторичной» модуляции, формируемого ПЛ, и спектром одиночного импульса «вторичной» модуляции, формируемого УЛ. Под огибающей на частотах Frq(п) = FdopK + nNLFV, п = ±1,±2,...,±Nrq
располагаются дискретные составляющие, обусловленные периодичностью ОС ПЛ и УЛ винта, где - число анализируемых спектральных составляющих. Составляющие
«вторичной» модуляции в левой части амплитудно-фазочастотных спектров (АФЧС) ОС НВ (или РВ) обозначим как п, а в правой части - как п.
На основе моделирования с использованием моделей ОС, описанных в [3], проведены исследования разностей фаз между соседними спектральными составляющими АФЧС ОС НВ и ОС РВ зависшего вертолета на примере Ми-2 (Vt = 0 м/с). Формировались и исследовались свертки (п+1) и ^п+1^п. Исследования позволили выявить следующие важные свойства АФЧС
«вторичной» модуляции ОС НВ и ОС РВ.
Первое - разности фаз между соседними спектральными составляющими для каждой (правой или левой) части АФЧС ОС НВ (или РВ) практически одинаковы.
Второе - фазы сверток в левой (п+1) и правой ^п+1^,*п частях спектра ОС НВ Ми-2,
у которого количество лопастей нечетное, отличаются на 180°, а для спектра ОС РВ Ми-2, у которого количество лопастей четное, практически одинаковы.
u
--1
-(n+ 2) -
0+1)
Планерная составляющая
- '«+1 -
I ^я+2
0
F
dopK
f - fo
Рис. 1. Спектр сигнала «вторичной» модуляции винта
Свойства АФЧС комплексной огибающей ОС винта являются физической основой для реализации свертки спектра «вторичной» модуляции ОС винта:
1. Когерентные свойства сверток соседних спектральных составляющих «вторичной» модуляции позволяют реализовать свертку всех спектральных составляющих отдельно для правой и отдельно для левой частей спектра «вторичной» модуляции. Свертка спектра «вторичной» модуляции в спектральной области формируется следующим образом:
для ОС ПЛ
и для ОС УЛ
_PL = Z - n+1-n ;
n=1
N„-1
_ UL = Z- -nS-(n+1) ■ n=1
(1а) (1б)
2. Когерентность сверток правой и левой частей спектра «вторичной» модуляции винта позволяет сформировать общую свертку спектра «вторичной» модуляции (ССВМ) винта:
5
пРи нечетном количестве лопастей _V = _PL - _UL = Z (-n+1-n ---n--(n+1)) ;
n=1
N„-1
при четном количестве лопастей
_V =_PL +-ZJUL = Z (-n+1-n + —-n—-(n+1) ) .
(2а) (2б)
n
2
2
N, -1
N, -1
n=1
Способ измерения частот следования лопастей несущего и рулевого винтов на основе свертки спектра «вторичной» модуляции сигнала
Способ измерения частоты следования лопастей винта может быть реализован в многоканальном устройстве (МКУ) свертки спектра «вторичной» модуляции (ССВМ), в произвольном канале которого вычисляется ССВМ ОС винта для измеряемой частоты Flv izm, устанавливаемой в диапазоне от FLV mm до FLV max с шагом AF. Свертка может быть
получена как во временной, так и в спектральной областях. Измеренное значение частоты FhV отождествляется с частотой настройки канала FLV imm , на выходе которого получен наибольший сигнал. При таком подходе потенциальная ошибка измерения частоты следования лопастей винта лежит в диапазоне ±0,5AF . Частота FLV imm (n) для канала с номером n описывается выражением:
FLv_„m (n) = Flv _mm + nAF, n = , Nf = ceil [(F^ - Fl^ )/AF] .
Структурные схемы каналов МКУ ССВМ с номерами n, в которых реализуется ССВМ ОС УЛ и ПЛ винта для частоты FLV imm (n) во временной области, представлены на рис. 2.
По аналогии с выражениями (1) эти схемы обеспечивают формирование во временной области сверток спектра «вторичной» модуляции ОС приближающихся лопастей
Ta
ZII_PL_V (FLV_zm (n)) = j UP_PL_V (t)U'P_PL_V (t)eXP [i2^FLV_zm (n)t]dt (3а)
0
и удаляющихся лопастей
ZEи_Г _ (п)) =\ иР_ЦЬ^ ^_V (0еХР _ от (п)']^ .
(3б)
ЦР РЬ V (0
^Л РЬ V ^ *т (п)) Цр
иь V (I)
Х^-ехр [-2р _ ^ (и)г]
иь V ( Fьv izm (и))
Х^- еХР _iZm (и)?]
а б
Рис. 2. Структурная схема п -го канала МКУ ССВМ для ОС ПЛ (а) и ОС УЛ (б) винта
В этих схемах одновременно выполняется автоматическая компенсация доплеровской частоты сигнала, отраженного от корпуса (планера) ЛА.
Структурные схемы устройств, реализующих способ измерения частоты следования лопастей НВ р ш и РВ р ¡^ одновинтового вертолета на основе ССВМ ОС винта во
временной области, представлены на рис. 3 и 4 соответственно.
иР
Цр
)
)
Рис. 3. Структурная схема устройства измерения частоты следования лопастей НВ
1
Линия задержки
и„
)
Устройство определения задержки 1-го импульса ОС ПЛ
Устройство определения задержки 1-го импульса ОС УЛ
^ _ РЬШ Устройство удаления ОСНВ МКУ
ССВМ
Т —I р _ ^—; п
Устройство МКУ
К _шш ОС НВ ССВМ
п = 0
| Линия задержки
Р _ ¡V
п = 0
Устройство
поиска максимумов и
определения частоты
II _ РЬ _ ¡V (РЬ ¡V)
Устройство оценки разности фаз и решения о
«четном» или «нечетном» числе лопастей
Оценка
Р
Признак р
иР
Л()
п = Ыг
^и_иь_¡^^ Ь_¡V
Р ¡V)
Рис. 4. Структурная схема устройства измерения частоты следования лопастей РВ
Принцип работы устройства измерения частоты следования лопастей НВ, представленного на рис. 3, следующий. ОС, из которого удалены отражения от корпуса и мешающие отражения (МО), поступает на вход устройства БПФ, на выходе которого формируется спектр. Затем относительно доплеровской частоты планера Р(0орК выделяются спектры ОС ПЛ и УЛ винтов. Эти спектры преобразуются во временные реализации ОС ПЛ иР РЬ v (0 и ОС УЛ ир иЬ v (^) путем обратного БПФ (ОБПФ). Далее отраженные сигналы ПЛ и УЛ винтов запоминаются и поступают на входы одинаковых МКУ ССВМ. На выходах МКУ
ССВМ формируются свертки ^ рЬ ш (р ш гт (п)) и Zп_иЬ_Ж (РЬ_Ж_гт (п)) . Свфтки
поступают в устройство поиска максимумов и определения частоты р ш, где формируется
Т
оценка частоты следования лопастей несущего винта Р МУ . В данном случае оценка частоты Р МУ соответствует частоте Рь ш ш (п) канала, имеющего наибольший квадрат модуля
I I2 I I2 "
свеРтки из рЬ ш (РЬ_Ш_гтг (п))| и ^II_иЬ_(РЬ _МУ _ ггт (п))| . ПРи3нак РЬМУ «четного»
(«нечетного») числа лопастей НВ формируется на выходе устройства оценки разности фаз и решения о «четном» или «нечетном» числе лопастей, в котором вычисляется разность фаз
ЛфZ ш сверток Хи РЬ МУ (Р МУ) и Хи иь МУ (РЬ МУ). Признак «нечетное» число лопастей формируется при модуле разности фаз |ЛфZ МУ|, близком к 180°, а признак «четное» число
лопастей - при модуле разности фаз |ЛфZ МУ |, близком к 0°.
Затем в устройстве, представленном на рис. 4, выполняется процесс измерения частоты следования лопастей рулевого винта РЬ КУ и принятия решения о «четном» («нечетном») числе
лопастей РВ. Отраженные сигналы ПЛ и УЛ винтов, сформированные путем ОБПФ в устройстве на рис. 3, поступают в устройства определения задержек первых импульсов «вторичной модуляции» в ОС ПЛ НВ РЬМУ и в ОС УЛ НВ иьМУ относительно начала
анализа. Затем на основе полученных задержек и периода повторения ТР МУ = 1 / р МУ путем временной селекции из ОС удаляются импульсы ПЛ и УЛ НВ. Полученные сигналы ир РЬ КУ (1) и ир иь КУ (г) для ПЛ и УЛ РВ поступают в два одинаковых МКУ ССВМ,
на выходах которых формируются свертки 2П РЬ RУ (РЬ RУ г?т (п)) и иь RУ(РЬ КУ кт (п)). Частота следования лопастей РВ Рь КУ определяется по наибольшему значению квадрата
I Р I |2 '
модуля свертки иь_КУ (РЬ_КУ _гт (п))| и ^II _РЬ_КУ (РЬ_КУ _ит (п))| . Признак РЩУ «четное»
(«нечетное») число лопастей РВ формируется по аналогии с формированием признака для НВ.
На основе измеренного значения частоты Рь МУ (или Рь КУ) можно получить оценку
частоты вращения НВ (или РВ) делением Рь МУ (или Рь КУ) на предполагаемое число лопастей
НВ (или РВ) и поиском реального сочетания частоты вращения и числа лопастей винта.
Заметим, что разработанный способ может также использоваться для измерения частоты следования лопастей винта двухвинтового вертолета и частоты следования лопастей винта (лопаток компрессора или турбины) самолета.
Результаты экспериментальных исследований способа
Условия исследований: вертолет Ми-2 удалялся от РЛД, зондирующий сигнал представлял собой МХЗС с круговой поляризацией и длиной волны Х = 0,0125м. Квадратурные составляющие принятого сигнала на видеочастоте преобразовывались в цифровой вид с частотой дискретизации 96 кГц. Обработке подвергался преобразованный ОС, из которого были удалены ОС корпуса и мешающие отражения. Время когерентного накопления выбрано равным Та = 1,365с. Границы анализируемого диапазона частот следования лопастей для НВ приняты равными Рь ш тп = 9 Гц и Рь ш тах = 20 Гц, а для РВ -
РЬ_КУ _тгп = 40 Гц и РЬ_КУ _тах = 80 Шаг изменения частоты ЛР = ° 01Гц .
На рис. 5 (для НВ) и рис. 6 (для РВ) представлены зависимости квадратов модулей выходных сверток МКУ ССВМ от анализируемой частоты. Для уменьшения уровней боковых лепестков каналов МКУ ССВМ было использовано весовое окно Хеммига. Анализ зависимостей показывает, что свертка спектра ОС УЛ (или ПЛ) НВ (или РВ) имеет максимум
в районе частоты следования лопастей НВ (или РВ) в диапазоне анализируемых частот. Кроме того, на рис. 6 имеются небольшие дополнительные пики, обусловленные влиянием остатков ОС НВ и расположенные на частотах, кратных частоте следования лопастей НВ.
а б
Рис. 5. Зависимости квадратов модулей выходных сверток МКУ ССВМ от измеряемой частоты,
а - для ОС ПЛ НВ, б - для ОС УЛ НВ
а б
Рис. 6. Зависимости квадратов модулей выходных сверток МКУ ССВМ от измеряемой частоты,
а - для ОС ПЛ РВ, б - для ОС УЛ РВ
Сравнение сверток, представленных на рис. 5, позволяет получить оценку частоты следования лопастей НВ - Рь му = 12,52Гц. Модуль разности фаз между свертками
^и РЬ (Р му) и Хи ш му (р му) составляет |Дф2 ^ = 188,48° и, следовательно, НВ имеет нечетное число лопастей. Исходя из рис. 6 оценка частоты следования лопастей РВ составляет р Яу = 48,97Гц, а число лопастей РВ является четным, так как модуль разности фаз между
свертками Хи _ рь фь _ Яу )и Хи _ Яу фь _ Яу ) составляет |Дф2 = 0,114° .
Полученные результаты измерений соответствуют параметрам вертолета Ми-2 и подтверждают практическую работоспособность разработанного способа.
Заключение
Разработанный способ измерения частоты следования лопастей винта на основе свертки спектра «вторичной» модуляции отраженного сигнала наряду с выполнением основной функции устраняет негативное влияние на измерение быстрого изменения радиальной скорости летательного аппарата. Это свойство обеспечивает принципиальную возможность решения задач измерения параметров винтов ЛА в РЛД, для которого типовым условием наблюдения является быстрое изменение ракурса объекта наблюдения. Разработанный способ обеспечивает также существенное повышение отношения «сигнал винта/помеха» за счет когерентного накопления свертки, что приводит к соответствующему
снижению ошибок измерения искомой частоты. Важным свойством способа является возможность определения признака четности числа лопастей винта. Измеренное значение частоты следования лопастей и признак четности используются для определения частоты вращения винта, знание которой необходимо для построения радиолокационного изображения винта.
Список литературы
1. Гейстер С.Р., Курлович В.И., Шаляпин С.В. Экспериментальные исследования спектральных портретов винтовых и турбореактивных самолетов в радиолокаторе обзора с непрерывным зондирующим сигналом // Радиолокация и радиометрия: сб. ст. из науч.-техн. серии. 2000. Вып. 3. С. 90-96.
2. Гейстер С.Р., Нгуен Т.Т. Облик радиолокационного датчика обнаружения и распознавания низколетящих целей, использующего обращенный синтез апертуры антенны // Наука и воен. безопасность. 2017. № 1 (51). С. 22-27.
3. Гейстер С.Р., Нгуен Т.Т. Математические модели радиолокационного сигнала, отраженного от винта вертолета, в приложении к обращенному синтезу апертуры // Наука и воен. безопасность. 2018. № 3 (57). С. 34-40.
4. Гейстер С.Р., Нгуен Т.Т. Маховые движения и изгибы лопастей несущего винта вертолета в модели отраженного сигнала // Наука и воен. безопасность. 2018. № 4 (58). С. 21-29.
References
1. Gejster S.R., Kurlovich V.I., Shaljapin S.V. Jeksperimental'nye issledovanija spektral'nyh portretov vintovyh i turboreaktivnyh samoletov v radiolokatore obzora s nepreryvnym zondirujushhim signalom // Radiolokacija i radiometrija: sb. st. iz nauch.-tehn. serii. 2000. Vyp. 3. S. 90-96. (in Russ.)
2. Gejster S.R., Nguen T.T. Oblik radiolokacionnogo datchika obnaruzhenija i raspoznavanija nizkoletjashhih celej, ispol'zujushhego obrashhennyj sintez apertury antenny // Nauka i voen. bezopasnost'. 2017. № 1 (51). S. 22-27. (in Russ.)
3. Gejster S.R., Nguen T.T. Matematicheskie modeli radiolokacionnogo signala, otrazhennogo ot vinta vertoleta, v prilozhenii k obrashhennomu sintezu apertury // Nauka i voen. bezopasnost'. 2018. № 3 (57). S. 34-40. (in Russ.)
4. Gejster S.R., Nguen T.T. Mahovye dvizhenija i izgiby lopastej nesushhego vinta vertoleta v modeli otrazhennogo signala // Nauka i voen. bezopasnost'. 2018. № 4 (58). S. 21-29. (in Russ.)
Сведения об авторах
Гейстер С.Р., д.т.н., профессор, главный сотрудник Научно-исследовательского института Вооруженных сил Республики Беларусь.
Нгуен Т.Т., аспирант кафедры информационных радиотехнологий Белорусского государственного университета информатики и радиоэлектроники.
Адрес для корреспонденции
20013, Республика Беларусь,
г. Минск, ул. Бровки, 6,
Белорусский государственный университет
информатики и радиоэлектроники
e-mail: [email protected]
Гейстер Сергей Романович
Information about the authors
Heister S.R., D.Sci., professor, main researcher of Research institute of the Armed forces of Republic of Belarus.
Nguyen T.T., PG student of information radiotechnologies department of Belarusian state university of informatics and radioelectronics.
Address for correspondence
220013, Republic of Belarus, Minsk, P. Brovka st., 6, Belarusian state university of informatics and radioelectronics e-mail: [email protected] Heister Sergey Romanovich