CC BY
33
2008
Доклады БГУИР
№ 3 (33)
УДК 534.23
ФИЗИЧЕСКИЕ ОСНОВЫ ФОРМИРОВАНИЯ АКУСТИЧЕСКОЙ ВОЛНЫ ВИНТОМ САМОЛЕТА И МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ВРЕМЕННОЙ СТРУКТУРЫ АКУСТИЧЕСКОГО СИГНАЛА ИДЕАЛЬНОЙ ЛОПАСТИ
А.М. ДЖЕКИ, А С. ГЕЙСТЕР
Белорусский государственный университет информатики и радиоэлектроники П. Бровки, 6, Минск, 220013, Беларусь
Поступила в редакцию 26 мая 2008
Описываются физические основы формирования акустической волны воздушным дозвуковым винтом легкомоторного самолета, а также математическая модель распределения давления вдоль вращающейся лопасти винта для зон повышенного и пониженного давления. Представляется оригинальная методика определения и математическая модель временной структуры акустического сигнала, формируемого идеальной лопастью.
Ключевые слова: воздушный винт, акустический сигнал.
Введение
Обнаружение легкомоторных летательных аппаратов (вертолетов и легкомоторных самолетов) на малых и предельно малых высотах является важной задачей. В приграничной зоне мониторинг малых высот осуществляется с целью пресечения ввоза наркотиков и контрабандных товаров, вывоза ценностей, а также исключения незаконного пересечения границы воздушным путем. Решение задачи мониторинга на малых и предельно малых высотах с использованием радиолокационных средств затруднено объективными факторами [1, 2], связанными с негативным влиянием земной поверхности на распространение радиоволн, низкой радиолокационной заметностью и малыми скоростями движения легкомоторных летательных аппаратов (ЛА). Это обусловливает необходимость разработки и создания новых эффективных средств обнаружения, в которых используются естественные физические явления, возникающие при полете легкомоторных ЛА на малых и предельно малых высотах. Одним из таких явлений является звук, порождаемый воздушными винтами.
Физические основы формирования акустического сигнала винтом самолета
На летательном аппарате тяга двигательной установки создается либо непосредственно за счет реакции струи, вытекающей из сопла двигателя, либо с помощью специального аэродинамического преобразователя — воздушного винта. Воздушный винт при вращении отбрасывает своими лопастями массу воздуха, проходящую через плоскость вращения винта, назад с увеличенной скоростью. За счет реактивного действия отбрасываемой массы воздуха на лопасти создается тяга винта [3].
В легкомоторных самолетах и вертолетах используются дозвуковые винты, которые при заданном диаметре имеют такое число оборотов, при котором результирующая скорость, равная сумме окружной скорости вращения винта и скорости полета, во всех сечениях лопастей
винта остается дозвуковой. Так как максимальная окружная скорость винта соответствует концу лопасти, то в данном случае рассматривается концевая скорость
Ук = 2пЯь / Ту, (1)
где Яь — радиус винта; Ту — период вращения винта.
На рис. 1 изображена модель двухлопастного винта, в котором лопасти, представленные прямоугольными пластинами, характеризуются относительно наблюдателя как приближающаяся и удаляющаяся. Воздействие на воздух происходит: во-первых, при "набеге" лопасти и проявляется как в сжатии воздуха на передней стороне лопасти, так и в его перемещении; во-вторых, в разрежении воздуха на обратной стороне лопасти.
Рис. 1. Модель двухлопастного винта
Для определения картины возмущения воздуха, возникающего при создании аэродинамической силы на лопасти, необходимо в первую очередь знать распределение давления по поверхности лопасти. Изменение давления в произвольной точке на поверхности лопасти по сравнению с давлением невозмущенного потока оценивают [3] коэффициентом давления
Р = (р - Р»)Рр (2)
где (р-рК) — избыточное (возмущенное) давление в данной точке на поверхности лопасти; рУ2/2 — скоростной напор (давление) невозмущенного потока плотностью р при скорости движения анализируемой точки поверхности у.
Определение распределения давления по поверхности реальной лопасти является сложным процессом, сочетающим как теоретические, так и экспериментальные исследования. Однако основные особенности формирования акустической волны воздушным винтом можно рассмотреть на примере идеальной лопасти, для которой предполагается отсутствие:
вихревых потоков, возникающих при затекании воздуха из зоны повышенного давления в зону пониженного давления на обратную сторону лопасти;
затекания воздуха из зоны повышенного давления в зону пониженного давления на краю лопасти.
Математическая модель временной структуры акустического сигнала идеальной лопасти
Приближающаяся идеальная лопасть. Оценим распределение давления в продольном сечении приближающейся идеальной лопасти. Для этого избыточное давление под вращающейся лопастью в произвольной точке, расположенной на расстоянии г,_ L от оси винта и движущейся со скоростью Vs_L, со стороны набегающего невозмущенного потока можно представить в виде
АРп (Г _L ) = (Pn (г,_Ь ) _ Р» ) = к
п PVS2_ L
р_а
2
(3)
где кпр_а — коэффициент, учитывающий угол атаки набегающей лопасти в точке.
С учетом того, что скорость движения точки лопасти Vs_L = 2пг,_L / ^ , на рис. 2 для приближающейся идеальной лопасти представлены закон изменения скорости точек V (гL ) на лопасти и закон распределения избыточного давления Арп (г,_ L ) вдоль лопасти.
АРп ( Г _ ь)
V,.
У
^Рп (Г _ Ь)
V,.
Г _ ь
' , _ N
<-►
Ось вращения винта
Лопасть винта
X
X
0
0
Рис. 2. Изменение скорости точек и избыточного давления вдоль лопасти
Звуковая волна в фазе повышенного давления, формируемая приближающейся лопастью, распространяется в соответствии с фронтом повышенного давления. Фронт и уровень избыточного давления определяются квадратичной функцией (3), что свидетельствует об отсутствии у лопасти, являющейся акустической излучающей антенной, единого фазового центра. Прием акустической волны осуществляется в дальней зоне, применительно к которой действует [1] фильтрующее свойство пространства. Пренебрегая временем распространения звуковой волны, это свойство заключается в том, что в данный момент в точку приема приходит волна, соответствующая фрагменту фронта давления, перпендикуляр к которому в данный момент совпадает с направлением на точку приема.
На рис. 3 (вид сверху) представлены пояснения к определению точки излучения на линии фронта избыточного давления для приближающейся идеальной лопасти.
Рис. 3. Определение точки излучения на линии фронта избыточного давления: X0Y — прямоугольная система координат лопасти; XV0YV — прямоугольная система координат, связанная с направлением на точку приема; — угол между касательной к фронту давления для точки лопасти, расположенной на расстоянии rs-L от оси вращения, и осью лопасти; <pL_n — угловое положение оси лопасти относительно оси XV, перпендикулярной направлению на точку приема; (n-<L-tn)) — угловое рассогласование для
точки анализа, расположенной на расстоянии rs-L от оси вращения
Точка лопасти, излучение которой для текущего положения лопасти придет в точку наблюдения, характеризуется касательной к фронту избыточного давления, перпендикуляр к которой соответствует направлению на точку приема. В математическом виде условие приема сигнала повышенного давления, формируемого в точке анализа на расстоянии rs—L , представляется выполнением равенства Лфtn s_L = 0, где ф f s—L = (п — <L-tn ).
Найдем аналитическое описание неизвестного угла фfs—L для каждой точки приближающейся лопасти, что позволит найти форму акустического сигнала в точке приема при вращении лопасти. Для этого воспользуемся коэффициентом давления
L) = ЛР"(r._L ) , (4)
г n—max
характеризующим распределение давления вдоль лопасти (см. выражение (3)) с нормировкой относительно максимального значения Лрп—max , соответствующего точке на конце лопасти:
ЛРи—max = Лр (R ) = kU ^ = kp—а , (5)
V
где VK — концевая скорость лопасти, определяемая выражением (1).
С учетом (3) избыточное давление Арп (гб,_ь) в точке лопасти, находящейся на расстоянии от оси вращения, определяется выражением
АРп (Г ь ) = кП. (6)
ТIV
Коэффициент давления рп ), характеризующий относительное распределение давления вдоль лопасти, является отношением выражений (6) и (5):
Рп Ь ) =
АРп (г,_ь ) = г2
Ар n_max
в?
(7)
Для нахождения угла фу,_ь необходимо взять производную от рп (г,_ь) по [_ь , которая определяется выражением
Рп (г,_ь ) =
ЛРп Ь ) Г
Лг,
= 2-
ВЬ
(8)
С учетом этого искомый угол фуь для любой точки идеальной лопасти, находящейся на расстоянии г,_ь от оси вращения, определяется как
фь = агс^! Рп (Г_ь) I = ^
С г >
2 ._ь
V
В
(9)
ь
С учетом того что угол фу,_ь изменяется монотонно, получим начальное фуN и конечное фук значения этого угла для идеальной приближающейся лопасти, при которых соответственно начинается и заканчивается формирование избыточного давления в направлении на точку приема. С учетом выбранной системы координат (рис. 3)
фуN = агс^ I Р(Вь ) I = агс^
С 2 ^
V Вь J
фу.,_к = агс^ I Р([_N) I = агс^
С [ >
2 . _ д
V
В
(10)
ь J
На основании выражения (9) для текущего углового положения лопасти фь_гп можно
найти точку на лопасти, излучение которой соответствует направлению приема. Расстояние до этой точки определяется выражением
>2
В
г
ь_гп
= Чггё (п_фь_п ),
2
а уровень избыточного давления равен
АРп (Гь_п ) = К _
р2п 2 г^,
гп
Р _а
(11)
(12)
Угловое положение лопасти фь_гп при вращении изменяется равномерно в соответствии со скоростью вращения Qv = 2п / ^ и в выбранной системе координат (см. рис. 3) в момент времени г равно фь_гп (г) = tQv . Абсолютное значение гг времени задержки акустического сигнала при распространении его от лопасти до точки приема в случае пассивной локации не играет большого значения, поэтому в целях упрощения опустим его из рассмотрения. С учетом этого математическая модель временной структуры акустического сигнала Ррг (г),
формируемого приближающейся идеальной лопастью, в точке приема описывается выражением
Рг ($) =
0 при 0 < г < гу.,_д
рп2 в4
КЛ
гг Р_а
т
2 (п_ гQv) при гу,_д < г < t
у ,_к >
(13)
0 при г
у_ к
< г
2
где кг — коэффициент, учитывающий уменьшение уровня давления при распространении на трассе "летательный аппарат - точка приема"; г^ 5 - N = (п-ф ^ 8-К )- — время начала фор-
0.у
мирования акустического давления при угловом положении лопасти ф ь - ш = (п - ф f 8 - N ), которое отсчитывается с момента г = 0, при котором фь-гп = 0;
г^ в-К = (п - фf ^-к )- — время окончания формирования давления при угловом положении
лопасти фь-ш = (п - ф f 3-к ) .
Удаляющаяся идеальная лопасть. При удалении лопасти на ее обратной стороне формируется зона разреженного давления, фронт которой вдоль идеальной лопасти также можно описать квадратичной функцией. По аналогии с приближающейся лопастью можно получить математическую модель временной структуры акустического сигнала рий (г), формируемого удаляющейся идеальной лопастью. Выражение для этого сигнала в точке приема имеет вид
Ри* 0) =
0 при 0<г <г
f.s-NV '
рп2 Я4
-кЛ-а V^2 (^ ) при < г < гг,-ку ; (14)
V
0 при гг,-KV < г
гг р-а 2
tf - К¥
ф f.я-NV
где кр_а — коэффициент, учитывающий угол атаки удаляющейся лопасти; г f ^ - ш
Qv
— время начала формирования акустического сигнала пониженного давления при угловом положении лопасти фь-гп = ф f 3-^¡у , причем время отсчитывается с момента г = 0, при котором
фь-гп = 0; К = фку — время окончания формирования акустического сигнала пони-Qv
женного давления при фь-гп = ф^ s-кv .
С учетом выбранной системы координат (рис. 3)
С г \ Г-Л С 2 ^
ф f.s - Nv = агй§ ^ Р(Г - N) J = агй§ ^ 2Я~ J, ф- ку = агй§ ^ Р( Яь) ^ = агй§
V Яь У
(15)
Формулы (13) и (14) позволяют математически описать временную структуру акустических сигналов, формируемых идеальной приближающейся и удаляющейся лопастью соответственно, в точке приема.
Заключение
Одним из физических явлений, которое может быть использовано для обнаружения легкомоторных летательных аппаратов на малых и предельно малых высотах, является возникновение звуковых колебаний. При полете винтового самолета акустические колебания или акустические волны создаются лопастями воздушного винта. Возле перемещающейся лопасти возникают поля повышенного и пониженного давления, что при круговом движении лопасти приводит к формированию двухполярного акустического сигнала. Описание этих полей является сложной задачей, при решении которой используются как теоретические, так и экспериментальные исследования. В данном случае разработана математическая модель временной структуры акустического сигнала, создаваемого идеальной лопастью. В дальнейшем эта модель
может быть использована для описания реальных лопастей, а также для построения алгоритмов обнаружения летательных аппаратов по их акустическим сигналам.
BASIC PHYSICS OF ACOUSTIC WAVE FORMING PERFORMED BY THE AIRCRAFT PROPELLER AND MATHEMATICAL MODEL OF ACOUSTIC SIGNAL TIME STRUCTURE OF THE IDEAL BLADE
ABDUALLA MOHAMED JAKI, A.S. HEISTER Abstract
Basic physics of acoustic wave forming performed by a small engine plane are described. Mathematical model of pressure distribution along the rotating propeller blade in pressure and depression zones is considered. Original method of determination and mathematical model of acoustic signal time structure formed by the ideal blade are represented.
Литература
1. Охрименко А.Е. Основы радиолокации и радиоэлектронная борьба. М., 1983, Ч. 1.
2. Гейстер С.Р. Системное проектирование и расчет радиолокаторов противовоздушной обороны. Минск, 1999. Ч. 1.
3. Вотяков В.Д. Аэродинамика летательных аппаратов и гидравлика их систем: Учеб.: В 2 ч. 1972. Ч. 1.
