CC BY
155
2010
НАУЧНЫЙ ВЕСТНИК МГТУ ГА серия Аэромеханика и прочность
№ 151
УДК 553.65.11.32:681.3:629.7.015
МОДЕЛИРОВАНИЕ ОБТЕКАНИЯ ЛОПАСТЕЙ НЕСУЩЕГО ВИНТА С РАЗЛИЧНЫМИ ЗАКОНЦОВКАМИ*
Б. С. КРИЦКИЙ
Описывается методика моделирования процесса обтекания лопастей несущего винта с законцовками произвольной формы, основанная на нелинейной теории тонкой несущей поверхности в нестационарной постановке. Приводятся результаты моделирования.
Ключевые слова: моделирование, обтекание лопастей винта, законцовки, нелинейная нестационарная вихревая теория.
1. Введение
Лопасти несущего винта при горизонтальном полете вертолета движутся по сложным криволинейным траекториям. Сечения лопастей при этом обтекаются с переменными по радиусу и азимутальному положению лопастей углами атаки и числами Маха и Рейнольдса. В наибольшей степени эти параметры изменяются в концевой части лопасти. При больших скоростях полета концевые сечения отступающей лопасти обтекаются со срывом потока, а наступающей - с трансзвуковыми скоростями. В этом случае на лопастях винта резко возрастает переменная часть шарнирного момента и, как следствие, увеличиваются нагрузки в системе управления. В силу этих особенностей на небольшом концевом участке лопасти реализуется значительная
часть тяги и крутящего момента несущего винта. Поэтому компоновка законцовки лопасти за-
метно влияет на аэродинамические характеристики несущего винта.
2. О методике моделирования обтекания лопастей
С целью расчетного исследования аэродинамических характеристик лопастей несущего винта сложной формы, включая законцовки различной геометрии, использовался метод расчета, основанный на вихревой теории тонкой несущей поверхности в нелинейной нестационарной постановке [1, 2]. Лопасти винта (рис. 1) заменяются бесконечно тонкими базовыми поверхностями Б; (1 - количество базовых поверхностей). При обтекании лопастей за ними имеется развитый след с поверхностями тангенциального разрыва скорости О] (] - количество поверхностей тангенциального разрыва скорости).
Среда рассматривается как идеальная несжимаемая. Везде вне лопастей винта Б; и их следов О] течение является безвихревым, т.е. для потенциала возмущенных скоростей Ф(х, у, 2, 1) справедливо уравнение Лапласа:
12Ф 12Ф 12Ф „ . \ гс \
+17=°- (х' у 2) г (8' и ч)'
Если W* - скорость движения точек несущей поверхности, обусловленная как поступательным движением, так и вращательным, а также маховым движением и деформацией, то в этих точках выполняется граничное условие непротекания:
(УФ - W*) •п = 0, (х, у, х) е Б;.
*Работа выполнена при поддержке РФФИ (проект № 08-08-00984-а).
Рис. 1. Схематизация лопастей несущего винта и следа
При переходе через поверхность вихревого следа С] соблюдаются условия непрерывности давления и нормальной составляющей скорости:
Р- = Р+; (УФй)_ = (УФй)+ , (X, у, 2) £ С].
Здесь индексы "_" и "+" относятся к разным сторонам поверхности С].
На тех кромках несущей поверхности, с которых стекают вихревые поверхности С], выполняется гипотеза Чаплыгина-Жуковского о конечности скоростей:
Р- = р+; (УФй)_ = (УФй )+, (х, у, 2) £ Ь].
Здесь Ь] - линия схода поверхностей тангенциального разрыва скорости.
На бесконечном удалении от винта, а также его следа возмущения затухают, поэтому
НшУФ = 0, Я = ^1 х2 + у2 + г2 .
Я®¥
Задача заключается в нахождении потенциала Ф(х, у, 2, 1) возмущенных скоростей на лопастях, а также во всем пространстве.
В рассматриваемой модели линии схода вихревых пелен постулируются всегда на острых задних кромках лопастей винта. При этом комлевые и передние кромки считаются закругленными. Поэтому положение зоны отрыва на лопастях определяется в процессе расчета обтекания винта.
Деформация вихревой пелены, моделирующей след, по времени может быть определена равенством:
(Х,^) = (Х,^1) + ] ^0 (Х,Л,я) dт,
т
где X, Л, С _ координаты точек пелены в момент времени т; (Х1, Ль СО _ в момент времени х1; ’0(Х, Л, С) _ компоненты безразмерной относительной скорости среды.
Для определения нагрузок на несущих поверхностях Б; используется интеграл Коши-Лагранжа. Численный метод решения задачи для несущего винта в нелинейной нестационарной постановке согласно методу дискретных вихрей заключается в дискретизации по пространству и времени. Непрерывные вихревые слои, которыми моделируются базовые поверхности лопастей винта и их вихревые следы, заменяются системами вихревых рамок, а непрерывный по времени процесс изменения граничных условий и параметров течения заменяется ступенчатым. Значения кинематических параметров остаются неизменными в рамках одного временного шага. На каждом временном шаге, начиная с первого, после решения системы линейных алгебраических уравнений, определяются напряженности всех вихревых рамок системы лопастей и следа за ними. По найденным напряженностям вихревых рамок с использованием интеграла Коши-Лагранжа определяются нагрузки на лопастях. Суммированием аэродинамической нагрузки по панелям определяются рас-
пределенные и суммарные характеристики винта. Проверка достоверности методики расчета осуществлялась на ряде этапов. Сопоставлялись результаты расчета с данными трубного и летного экспериментов, установлено удовлетворительное согласование результатов расчета с данными физических экспериментов по определению суммарных и распределенных аэродинамических характеристик несущего винта, а также полей скоростей в окрестности несущего винта [1 - 3].
3. Некоторые результаты моделирования
На первом этапе моделировался процесс обтекания лопастей невозмущенным потоком без их вращения. Исследовались лопасти без геометрической крутки с различными законцовками (рис. 2): прямоугольной лопасти Я, со стреловидной законцовкой 2, с законцовкой двойной стреловидности X и с законцовкой сложной формы Б. На рис. 2 приведено также распределение коэффициента нормальной силы сечения по длине лопасти в концевых сечениях, начиная с относительного радиуса 0,6. Форма законцовок существенно влияет на нагружение лопасти в концевых сечениях. Анализ формы следа за лопастью, поля скоростей в окрестности законцовок лопастей позволяет объяснить механизм такого аэродинамического нагружения лопастей.
X
Рис. 2. Аэродинамическое нагружение концевых сечений лопастей различной формы
На втором этапе исследовалось обтекание потоком этих же лопастей в системе вращающегося винта. Моделировалось обтекание четырехлопастного несущего винта. Относительная скорость движения винта У/юнЯ равнялась 0,2 (V - скорость поступательного движения несущего винта, юнЯ - окружная скорость концов лопастей).
Особенность обтекания заключается в том, что интенсивность вихревых жгутов, сходящих с концов лопастей, зависит от формы законцовок. Поле скоростей в области концевых сечений за лопастью также зависит от формы законцовки. Это, в свою очередь, оказывает влияние на распределение нагрузки по длине лопасти.
На рис. 3 в качестве примера показаны места пересечения концевых вихревых жгутов от впереди идущих лопастей с лопастями № 1 и № 4 для двух типов законцовок - прямоугольной формы Я и сложной формы Б. Анализ распределения коэффициента нормальной силы сечения по длине лопасти показал, что у прямоугольной лопасти вихревой жгут более концентрированный, и он создает пики нагрузки в комлевой части лопасти и на относительном радиусе ~ 0,65 (рис. 4). Лопасть с законцовкой Б нагружается более плавно, а в концевых сечениях меняется по закону, обусловленному формой законцовки.
При определении коэффициента нормальной силы сечения скоростной напор подсчитывался по окружной скорости конца лопасти.
При другом азимутальном положении лопасти (лопасть № 4) вихревые жгуты проходят ближе к концевым сечениям лопасти, а большая ее часть обтекается менее возмущенным потоком. Поэтому коэффициент нормальной силы сечения лопастей Я и Б № 4 протекает более
плавно (рис. 5). Отличия наблюдаются в концевых сечениях, там, где большая интерференция вихрей и лопасти.
а б
Рис. 3. К взаимодействию лопастей с концевыми вихревыми жгутами: а - прямоугольные законцовки Я; б - законцовки сложной формы Б
Рис. 5. Распределение Су по радиусу 4-й лопасти
Таким образом, форма концевой части лопасти существенно влияет на аэродинамические характеристики винта, отработка ее компоновки требует учета многих факторов.
ЛИТЕРАТУРА
1. Белоцерковский С. М., Локтев Б. Е., Ништ М.И. Исследование на ЭВМ аэродинамических и аэроупру-гих характеристик винтов вертолетов. - М.: Машиностроение, 1992.
2. Крицкий Б. С. Математическое моделирование аэродинамики винтокрылого летательного аппарата // Научный Вестник МГТУ ГА, серия Аэромеханика и прочность, № 5 9, 2003.
3. Крицкий Б. С. Моделирование ближнего следа за несущим винтом // Научный Вестник МГТУ ГА, серия Аэромеханика и прочность, № 37, 2001.
THE SIMULATION OF FLOW AROUND A LIFTING ROTOR BLADES WITH DIFFERENT WINGTIPS
Kritsky B. S.
Describe the method of modeling of flow around a rotor blades with the wingtips of arbitrary shape, based on the nonlinear theory of thin lifting surface in non-stationary formulation of problem. The results of the modeling are presented.
Сведения об авторе
Крицкий Борис Сергеевич, 1949 г.р., окончил ВВИА им. Н.Е.Жуковского (1976), доктор технических наук, профессор, главный научный сотрудник ЦАГИ, профессор ВВИА им. Н. Е. Жуковского и Ю.А. Гагарина, автор более 120 научных работ, область научных интересов - численные методы в аэрогидродинамике, аэродинамика винтокрылых летательных аппаратов.
