CC BY
220
УДК 629.735.33.014.16:533.662.6
РАСЧЕТ ИНДУКТИВНЫХ СКОРОСТЕЙ НЕСУЩЕГО ВИНТА НА ОСНОВЕ МОДЕЛИ С НЕПРЕРЫВНО РАСПРЕДЕЛЕННОЙ ЦИРКУЛЯЦИЕЙ ВИХРЕВОГО СЛОЯ
Б.Л. АРТАМОНОВ
Статья представлена доктором технических наук, профессором Ципенко В.Г.
Установлена связь погонных циркуляций радиального и азимутального вихревых слоев с циркуляцией присоединенного вихря, изменяющейся как по радиусу лопасти, так и по её азимутальному положению. Для типовых законов распределения циркуляции присоединенных вихрей вычислены три компонента индуктивной скорости в характерных точках диска винта. В том числе в точках, лежащих на хорде лопасти при её различных азимутальных положениях и режимах работы винта, что представляет практический интерес для расчета его аэродинамических характеристик.
Ключевые слова: индуктивная скорость, несущий винт, циркуляция, верхний слой.
В работе [1] описана линейная вихревая модель воздушного винта, представляющая собой пространственную геликоидальную пелену заданной геометрии, покрытую непрерывно распределенным двухкомпонентным вихревым слоем. Элементами дискретизации пелены являются треугольные панели, произвольно ориентированные в пространстве. Разработаны метод, алгоритмы и программы расчета трех составляющих вектора индуктивной скорости от системы произвольно ориентированных площадок, покрытых вихревым слоем, интенсивность которого линейно изменяется по поверхности площадки. Решение получено в элементарных функциях. Рассмотрим применение этой модели для несущего винта вертолета.
Зададимся законом изменения циркуляции присоединенных вихрей по диску винта Г(г,щ) и установим её связь с погонной циркуляцией сошедшего вихревого слоя у(г,щ). Это даст возможность решить прямую задачу аэродинамического расчета винта - определение поля индуктивных скоростей по заданному закону циркуляции.
Закон изменения циркуляции присоединенных вихрей получим с учетом её изменения по азимутальному положению лопасти. Воспользуемся уравнением, связывающим циркуляцию в сечении лопасти с его геометрическими и кинематическими параметрами [2]
Г(т ,щ) = 1 су (Г ,Щ)Ь (г )Ж (г ,щ), (1)
где Ь (г ) - закон изменения хорды по радиусу лопасти;
с (г,щ) - закон изменения коэффициента подъемной силы по диску винта;
Ж(Т, щ) - суммарная скорость потока в сечении лопасти.
Для трапециевидных лопастей зависимость Ь (г) имеет вид
Ь = Ь (г) = Ь Л - (л- 1)г ] / [0,7 + 0,3л], (2)
где Ь7 - хорда в характерном сечении лопасти на радиусе г = 0,7Я;
П - сужение лопасти.
При вычислении индуктивной скорости в первом приближении учтем изменение су только по азимуту в следующем виде
Су(г ,Щ) ~ СуО ■ /(Щ) ,
где с - среднее значение коэффициента подъемной силы по диску винта;
/(щ) - закон изменения коэффициента подъемной силы по азимутальному положению лопасти, характерный для шарнирного несущего винта.
Для их вычисления воспользуемся формулами Л.С. Вильдгрубе [3], полученными для шарнирного несущего винта:
3с
c =_3ст_• (3)
У0 п , 1 о 2\' V~V
У ка-, (1 + 1,2j )
f(щ) = (1 + j)/(1 + 2,67^тщ+2jj sin2щ), (4)
где к - коэффициент концевых потерь;
j - коэффициент режима работы винта j = V cosaB .
При вычислении скорости воздушного потока в сечении лопасти W(r, щ) будем считать, что
W(r ,щ) « Üx (r, щ) = r + (V cosaB + vlx)sin^ , (5)
а поле индуктивных скоростей vlx равномерно распределено по диску винта.
Применяя уравнение связи (1) в рамках лопастной теории, следует, по рекомендациям Прандтля, ввести поправку на конечность числа лопастей через функцию В.И. Шайдакова f (r, щ), сводящую циркуляцию на конце лопасти в ноль [4]
кл(1-7 )W
f (r ,щ) = -arccos(e 27ÜÜ ). (6)
ж
После чего закон (1) принимает вид
r(r , Щ) =1 cy (r, Щ)Ъ (r )ÜX (r, y}^ (r, щ). (7)
Для перехода от циркуляции присоединенных вихрей Г(г , щ) к погонной циркуляции свободного вихревого слоя /(r, 5) необходимо вычислить градиенты закона Г(г , щ) по направлениям образующих поверхности S:
/г (г, 5)= d Г(г , 5); / (r, 5)= d Г(г , щ) (8)
dr as
в её дискретных узловых точках 5у(ху, yy, zy).
Первое уравнение характеризует интенсивность радиального вихревого слоя за счет изменения циркуляции присоединенных вихрей по радиусу лопасти, второе - интенсивность тангенциального слоя за счет изменения циркуляции присоединенных вихрей по азимуту.
Перейдем в выражении (8) от производной по длине дуги ds к производной по азимутальному положению лопасти dщ. Поскольку на пелене имеет место соотношение
ds =yj dx 2+dy 2+dz2, то, дифференцируя выражения для координат пелены [1] по углу
x = r sin^-5) -5Vx(r ,5); y = -5Vy(r ,5); z = Fcos^-5) ,
получим
ds r 2+Vx2+Vy2 + 2Vx мп(ще -5)d5 .
Теперь соотношение (8) для /(r, 5) можно представить в виде
ñ,i . _9 1 _ —dr(r5) = Vr 2+Vx2+Vy2 + 2VxSin(щл -5) d5
1 (9) /5 (r ,5).
Ф
¡Г 2+^2+Ку2 + 2Кх81И(^л -5) Связь между углами 5 в (9) и у в законе Г(г, у) устанавливается соотношением
У] = Ул -
Отметим, что для корректного описания типового закона изменения циркуляции (7) точки
по радиусу лопасти г следует распределять неравномерно со сгущением к её концу. Хороший результат дает синусоидальное распределение с координатами
ТП
r = r + (1 -r )sin2(-) , 0 < i < Nr.
1 о 0J v2 fl"
(10)
Приведем результаты расчетов по описанному алгоритму типового закона циркуляции применительно к двухлопастному винту о7 = 0,06, кл = 2, работающему на режиме V = 0,1, ав = 0 и ст = 0,01. На рис. 1 показаны исходный закон циркуляции Г(ги соответствующие ему значения погонных циркуляций радиального /г (r, $) и тангенциального у9 (r, $) вихревых слоев в пределах одного оборота винта с шагом по азимуту А\|/ = = 20°.
Рис. 1. Исходный закон циркуляции Г(г ,щ) и соответствующие ему законы погонных циркуляций радиального /г(г,3) и тангенциального^(г ,3) вихревых слоев
В соответствии с постановкой задачи погонная циркуляция вихревого слоя в узловых точках пелены у(Т,3) должна быть задана в базовой декартовой системе осей координат ОХУ1, связанной с плоскостью и осью вращения винта
В проекциях на оси базовой системы осей координат получим:
(11)
У y (Г,3,)=У
x4+l "X4 -+7ÁKA) "x4
S ij+i -Su »Vi+ij -Su
Уцп +УМАУ
S ij+i -Su 1+1 j -Su
"Z4 "Zy
S ij+i -Su Si+ij -Su
(12)
Найденное поле циркуляций в узловых точках вихревого слоя даёт возможность для вычисления индуктивной скорости Упел(хА, у¡,~z ¡ ) в произвольной точке А пространства воспользоваться формулами (40) и (35) из работы [1]. Практический интерес для расчета аэродинамических характеристик винта представляют скорости в плоскости вращения лопастей при уА = 0.
Будем предполагать, что ось лопасти совпадает с линией начала вихревой пелены (3 = 0). Чтобы обеспечить возможность вычисления индуктивной скорости непосредственно на оси лопасти ( хА = r cosy, yA = 0, zA = r sin y ), скорректируем форму первого участка вихревого следа, обеспечив подход продольных векторов yr (7[, <9П) по перпендикуляру к ней. Для этого координаты конечных точек первого участка пелены должны быть пересчитаны следующим образом:
xi,i =
cos A3
cos(y -A3); yn = VA3;
cos A3
sin(y -A3).
(13)
На рис. 2 показаны в развертке по радиусу законы изменения по диску винта трех компонентов индуктивной скорости ух (г ,щ), V (г, щ), у2 (г, щ) от вихревого следа, сходящего с лопасти
на азимуте щл = 0, при числе витков геликоидальной спирали поб = 10. Следует обратить внимание на значительные «всплески» индукции в концевых сечениях лопасти на азимутах щ « 300°-360°, что свидетельствует о моделировании зарождающегося концевого вихря, сходящего с концевой части лопасти.
Оценим рациональную длину вихревого следа, обеспечивающую необходимую точность расчетов. На рис. 3 показана зависимость по числу учтенных витков вихревого слоя средней по диску индуктивной скорости, вызываемой вихревым следом трехлопастного винта, от каждой лопасти, находящейся на азимуте щЛу = 2п(у-1)/к
V^ = 1 | | ^пелу (г , щ, щЛу )лШщ , (14)
0 0
и от всех лопастей вместе
Vp =
1 кл 2г Г- - -
_ 2 J J УпелУ (r, ¥, ¥"v )rdrdy.
Я v=1 0 0
(15)
х
r
z
r
0.15
0.1
0.05
-0.05
-0.1
-0.15
0 0 2 0 4 0 6 ^ 0 ёг 1
-О
-20 40 60 -80 -100 -120 -140 160 180 200 220 240 260 280 300 -320 340 360
0.12 0.1 0.08 0.06 0.04 0.02 О
-0.02 ( -0.04 -0.06 -0.08 -0.1
Т (
1 1 ^
о гГ
—Г 1
-О
-20 40 60 -80 -100 -120 -140 160 180 200 220 240 260 280 300 -320 340 360
ч»
0.08 0.06 0.04 0.02
О
(
-0.02 -0.04 -0.06 -0.08
у Л— \
0 0 2 0 4 0 6
-О -20 40 60 -80 -100 -120 -140 160 180 200 220 240 260 280 300 -320 -340 360
Рис. 2. Законы изменения компонентов вектора индуктивной скорости в плоскости диска винта о7 = 0,06, кл = 2, работающего на режиме V = 0; ав = 0 и ст = 0,01, вызванной вихревой пеленой, сошедшей с лопасти, находящейся на азимуте = 0 (в развертке по радиусу)
В заключение приведем законы распределения индуктивной скорости по диску винта, полученные на основе описанной выше лопастной модели с непрерывно распределенной погон-
ной циркуляциеи по поверхности вихревого слоя, подсчитанные суммированием индуктивных скоростей от вихревых пелен, сошедших со всех лопастей
(16)
0
-0.005 -0.01 -0.015 -0.02 -0.025
01 > ; ( 5 > ( 8
\
п
об
♦ 1 лопасть ■ 2 лопасть а 3 лопасть --Винт в целом
ПуСР
Рис. 3. Средняя по диску индуктивная скорость, вызываемая вихревым следом от каждой лопасти и винтом в целом для разного числа витков геликоидальной спирали от поб = 1 до поб = 8 для трехлопастного винта о7 = 0,06, кл = 3 при V = 0,1; ав = 0 и ст = 0,01
На рис. 4 приведены зависимости компонентов вектора индуктивной скорости \ (г ,щ),
V (г, щ), \ (г ,щ) в плоскости диска винта о7 = 0,06, кл = 3, работающего на режиме V = 0; ав = 0
и ст = 0,01. Видно, что продольная и боковая составляющие вектора индуктивной скорости на два порядка меньше, чем нормальная к плоскости диска винта. Скорости V (г, щ) на азимутах
между лопастями значительно меньше, чем в районе лопастей. Выводы
1. Установлена связь законов распределения погонной циркуляции двухкомпонентного непрерывного вихревого слоя с геометрией лопасти и режимом работы шарнирного несущего винта.
2. Для типового закона циркуляции Г(г,щ) вычислено поле индуктивных скоростей в плоскости диска винта на режиме косого обтекания. Показаны значительные «всплески» индукции в концевых сечениях лопасти на азимутах щ « 300° - 360°, что свидетельствует о моделировании зарождающегося концевого вихря, сходящего с концевой части лопасти.
3. Установлено, что для корректного моделирования закона изменения индуктивной скорости в плоскости диска винта следует учитывать не менее 4-5 витков его вихревой системы.
4. Показано, что продольная и боковая составляющие вектора индуктивной скорости на два порядка меньше, чем нормальная к плоскости диска винта. Скорости V (г, щ) на азимутах между лопастями значительно меньше, чем в районе лопастей.
У=1
Рис. 4. Законы изменения компонентов вектора индуктивной скорости в плоскости диска
винта а7 = 0,06, кл = 3, работающего на режиме V = 0; ав = 0 и ст = 0,01, вызванной вихревой пеленой, сошедшей со всех лопастей винта (в развертке по радиусу)
ЛИТЕРАТУРА
1. Артамонов Б.Л. Вихревая модель воздушного винта с непрерывно распределенной циркуляцией вихревого слоя // Статья в настоящем Вестнике.
2. Шайдаков В.И. Обобщенная дисковая вихревая теория и методы расчета индуктивных скоростей несущего винта вертолета // Проектирование вертолетов. - М.: МАИ, 1977. - Вып. 406.
3. Баскин В.Э., Вильдгрубе Л.С., Вождаев Е.С., Майкапар Г.И. Теория несущего винта. - М.: Машиностроение, 1973.
4. Шайдаков В.И., Маслов А.Д. Аэродинамическое проектирование лопастей воздушного винта. - М.: МАИ, 1995.
CALCULATION OF INDUCED VELOCITIES OF ROTOR ON THE BASE OF THE MODEL WITH CONTINUOUSLY DISTRIBUTED CIRCULATION OF VORTEX LAYER
Artamonov B.L.
The relationship of linear circulation of radial and azimuthal vortex layers with circulation of bound vortex was created. This circulation of bound vortex varies depending on radius and azimuthal location of blade. Three components of induced velocity in the characteristic points of rotor disk area were calculated for the typical laws of circulation distribution of bound vortexes. Among other characteristic points there are the points lying in the blade chord, when this blade is situated in different azimuthal locations and operating modes of rotor. It is of practical interest in rotor characteristics computation.
Keywords: induced velocity, rotor, circulation, upper layer.
Сведения об авторе
Артамонов Борис Лейзерович, 1947 г.р., окончил МАИ (1972), кандидат технических наук, старший научный сотрудник, заместитель заведующего кафедрой проектирования вертолетов МАИ (Национального исследовательского университета), автор более 180 научных работ, область научных интересов - аэромеханика винтокрылых летательных аппаратов вертикального взлета и посадки.
