Вихревая модель крыла с непрерывно распределенной циркуляцией вихревого слоя Текст научной статьи по специальности «Механика и машиностроение»

УДК 629.735.33.014.16:533.662.6

ВИХРЕВАЯ МОДЕЛЬ КРЫЛА С НЕПРЕРЫВНО РАСПРЕДЕЛЕННОЙ ЦИРКУЛЯЦИЕЙ ВИХРЕВОГО СЛОЯ

Б.Л. АРТАМОНОВ

Статья представлена доктором технических наук, профессором Ципенко В.Г.

Рассматривается линейная вихревая модель крыла конечного размаха, представляющая собой плоскую прямоугольную пространственную пелену, покрытую непрерывно распределенным вихревым слоем. Элементами дискретизации пелены являются четырехугольные панели, лежащие на ее поверхности. Созданы метод, алгоритмы и программа расчета трех составляющих вектора индуктивной скорости от произвольно ориентированной прямоугольной площадки, покрытой вихревым слоем, интенсивность которого линейно изменяется по поверхности площадки. Решение получено в элементарных функциях. Численным способом решена задача определения закона циркуляции присоединенных вихрей по размаху крыла и расчета его аэродинамических характеристик, основываясь на принятой вихревой модели и гипотезе плоских сечений.

Ключевые слова: крыло, циркуляция, вихревой слой, линейная вихревая модель.

Вихревая модель трапециевидного крыла конечного размаха разрабатывается на основе теории несущей линии и гипотезы плоских сечений с целью создания численного метода расчета его аэродинамических характеристик на основе взятых из эксперимента аэродинамических характеристик профилей. Предлагаемая вихревая модель позволяет избежать особенностей при вычислении индуктивных скоростей как на самой вихревой пелене, так и в контрольных точках поверхности, создающей подъемную силу.

В вихревой теории крыла конечного размаха ключевым является вопрос о выделении особенности при расчете индуктивных скоростей, вызываемых системой вихрей на самой вихревой поверхности. Подробно эта проблема изложена в работе академика Б.Н. Юрьева [1], где анализируются различные численные приемы решения этой задачи как в рамках дискретной, так и в рамках непрерывной модели вихревой пелены. Однако точного решения этой задачи им получено не было.

Практический интерес представляет расчет аэродинамических характеристик крыла в широком диапазоне углов атаки, включающем в себя как до-, так и засрывные режимы обтекания. Решению этой задачи в рамках линейной вихревой модели крыла были посвящены работы Г.Ф. Бураго [2], Г.И. Майкапара [3], Г.В. Каменкова [4], Т.К. Сиразетдинова [5].

В работах [2; 3] задача определения закона распределения циркуляции по размаху крыла решалась методом последовательных приближений в предположении дискретного распределения циркуляции. При такой постановке задачи основное затруднение заключалось в неоднозначности и неустойчивости решения при резком изменении коэффициента подъемной силы профиля на закритических углах атаки. В работах [4; 5] был применен прием разложения закона циркуляции в ряд и сведения задачи к системе линейных уравнений. В этом случае при плавном характере зависимостей коэффициента подъемной силы профиля по углу атаки и циркуляции по размаху крыла уже в первом приближении получаются решения, хорошо совпадающие с экспериментальными данными. Для неплавных законов высокая точность решения не гарантируется.

С.Г. Бураго разработал численный метод расчета распределения подъемной силы по размаху нестреловидного крыла на закритических углах атаки. В соответствии с теорией Прандтля предполагалось, что на концах крыла циркуляция и коэффициент подъемной силы равны нулю. Метод расчета был опробован как на гладких зависимостях сур(а), так и на "обвальных" характеристиках в зоне срыва. На некоторых углах атаки наблюдался расходящийся итерационный процесс, несмотря на применение коэффициента релаксации.

Общая постановка задачи

Рассмотрим плоское прямое трапециевидное крыло размахом L, установленное под углом ф к вектору скорости набегающего невозмущенного потока V. Расположим крыло вдоль оси OZ правой декартовой системы координат так, чтобы его середина совпадала с началом координат (рис. 1). Закон изменения хорды такого крыла Ь по размаху

ъ в относительных величинах Ь = 2Ь / Ь, ъ — 2ъ / Ь будет иметь вид

Ь(2) = —Ц |-(Л"1)2, (1)

Ц1 + Л) "

где г| = Ь(0)/Ь(1) - сужение крыла; X = Ь/Б -удлинение крыла; S - площадь крыла.

Будем предполагать, что сечения крыла образованы профилем, аэродинамические характеристики которого Сур = Сур(а), схр = схр(а) известны и заданы на основе круговой продувки в аэродинамической трубе в диапазоне углов атаки -180° < а < 180°. Принимая гипотезу плоских сечений в скоростной системе осей координат, можно записать выражения для аэродинамических сил, действующих на элемент крыла шириной сЬ:

(2)

где '(г) = -^/у2 +уу2^) - результирующий вектор скорости в сечении крыла; уу - вертикальный компонент индуктивной скорости в сечении крыла; Р(г) = Аг^[уу(2)/У] - угол притекания потока в сечении крыла; а - угол атаки результирующего вектора скорости в сечении крыла

а(2) = Ф(2)-Р(2). (3)

В отличие от обычно применяемых формул [1], здесь в обоих выражениях (2) присутствуют и коэффициент подъемной силы сур и коэффициент сопротивления схр профиля. Это позволяет рассчитывать аэродинамические характеристики крыла в широком диапазоне углов атаки, в том числе и за срывом, если известна круговая продувка профиля.

Интегрируя элементарные аэродинамические силы (2) по размаху крыла и переходя к коэффициентам этих сил, получим:

<4)

2

Сук = — |[сур(а)с08/?(г) + схр(а)81п /3(г)]]¥ (г)Ь(г)с1г;

Схк — '

£

'(7) = + РЙ = Лгс1в |у (7)

где У^^^^ + му2^), Р(2) = АпЛё|у(г)_. (5)

Запишем уравнение связи циркуляции вектора скорости в сечении крыла с его геометрическими, кинематическими и аэродинамическими характеристиками

Г(2) = ^сУР(а)Ь(2)\У(2).

(6)

Будем предполагать, что с крыла сходит плоская вихревая пелена, направленная по вектору скорости невозмущенного потока, который набегает на крыло параллельно оси ОХ.

Погонная циркуляция вихревой пелены y(z) изменяется непрерывно вдоль размаха крыла и численно равна градиенту изменения циркуляции по размаху крыла

y(z) = Af(z). (7)

dz

Стремясь получить аналитические выражения для расчета индуктивной скорости в произвольной точке пространства от плоской вихревой пелены, разобьем ее сечениями zi (1 < i < N) на отдельные полосы и будем предполагать, что погонная циркуляция непрерывна по поверхности пелены и изменяется в пределах полосы (zi, zi+1) по линейному закону:

y(z) = Az + B; Zi<z<zM; l<i<N-l, (8)

где А, В - коэффициенты, зависящие от значений у на краях полосы:

Zi+l — Zi Zi+l — Zi

Отметим, что линейный закон распределения погонной циркуляции вихревого слоя y(z) равносилен квадратичному закону распределения циркуляции r(z) по размаху крыла. Такой прием аппроксимации был впервые предложен профессором А.В. Чесаловым [6].

Индуктивная скорость, создаваемая вихревой полосой с линейно распределенной погонной циркуляцией

Рассмотрим участок вихревой пелены, расположенный в плоскости XOZ правой системы осей координат (рис. 2). Положим, что погонная циркуляция вихревого слоя на границах полосы zi+1) известна и равна (7^ ум), а между этими точками изменяется линейно. Вычислим индуктивную скорость V, создаваемую вихревой полосой в точке Га = {ха,уа^а}.

Воспользуемся уравнением Био-Савара, записанным в векторной форме [7]

= (9)

4тх|г

где dГ - вектор-элемент вихря, создающий индуктивную скорость dv; г - вектор, соединяющий элемент вихря с точкой, в которой вычисляется индуктивная скорость.

Зададим точку вихревой пелены, из которой исходит вектор-элемент dГ радиус-вектором р = {х,0,г}. Тогда, как следует из рис. 2:

' = = ~ (Ю)

Рис. 2. К расчету индуктивной скорости от вихревой полосы

Г=ГА- Р= XÄ-X,yA,ZÄ-Z ; |Г| = V(XA " Х)2 + УА + (ZA " Z)2 .

Интенсивность элемента вихря dГ вычисляется через погонную циркуляцию вихревого слоя на отрезке ей = е{хД0

dГ = y(z)dS. (11)

Записывая интеграл от (9) по всей вихревой полосе, имеем

1 z'+i dSx r

v = — Jy(z)dzJ—-

z S Г

Раскрывая входящее в (12) векторное произведение и учитывая (10), получим выражение для вектора индуктивной скорости у{ух,уу,у2} через двойной интеграл по поверхности вихревой пелены. В проекциях на оси координат соотношение (12) распадается на три

V. = 0

^ г1+1 00 у Л

дх

0 (Ха-Х) + УА + (?Л-*У дх

3/2

(13)

|*а-*)2+ У А+ ( )2

3/2

После интегрирования будем иметь в относительном виде для компонента уу

_1_

4л-

уу =—\А

ха 1п

-г, + <2(гг)

+ Ая1^(Ха)Уа

4(

атеЩ

Ха

2а -гм +()(гг+\) (г а —гг+\)

+ Уа

атеЩ-

1а - гм

Уа

-атеЩ-

2 А-

+ (Zi+1 -гг)

УаЯ^м)

- aтеtg

Ха

(га - гг)

Уа0(*)

+ (Агл + Л)1п

Уа у

(<2(гг+1 )-ха)(<2&)-ХА)

(гА-г^+У;

-2 -2

где 0(г) = ^¡(гл - г) + ха + уА .

Особый интерес представляет вычисление индуктивной скорости в плоскости XOZ, где находится вихревая поверхность. Для этого положим в полученном выражении уА = 0. Эта операция связана с раскрытием неопределенности типа Ох °о, что приводит к выражению

Уу(хл,7л) = — <! А 4п

хл 1п

7л -71 + 0(7,)

+ (2^1 )

+(А7А +В)1п

(0(7н.1)-хл)(0(70-хл)

-2

(7^71) +ул

(14)

7л -71+1+0(71+1)

Скорости в плоскости начала вихревой пелены находим, полагая хА = 0

0(7^1)00(7,)

(7а - г;)2

Здесь С>0 - функция С>(г) при ха = уА = 0: £>0(гм) = \гА -гг+1|; <20(г,) = \гА -г\.

Подставляя эти значения в (14), а также коэффициенты А и В по формулам (8), после упрощений получим выражение

1

Уу0(7л) = — л(7,+1 -7,)+(л7л +В)1п 4т1

уу0 (7л) =

4%

71+1 71

— —

7л - 71+1

7л "71

(15)

Как видно из (15), на краях полосы при 7А = 7;и 7А = 7;+1 имеет место особенность логарифмического типа. Это не дает возможности вычислить индуктивную скорость непосредственно на краю полосы при отличной от нуля циркуляции в этой точке. Если же циркуляция на краю полосы обращается в ноль, то скорость оказывается величиной конечной.

При расчете характеристик крыла задача упрощается, поскольку соседние полосы имеют общую границу, на которой особенность отсутствует и скорости конечны. Если рассмотреть две соседствующие друг с другом вихревые полоски ( ) и ( ), то индуктивная скорость в

начальной точке линии их соприкосновения 7А = 7; будет равна

V

vy0(zi) =

An

(Yi+l -YiHYih

zi - zi+l

zi — zi-l

(16)

При равной ширине участков z; — zi+1 = z; -z^.

скорость в точке zi не зависит от величины связанной с ней циркуляции и определяется только циркуляцией соседних точек.

Формулы (15) и (16) полностью решают задачу определения индуктивных скоростей во всех сечениях крыла ( -1 < гА <1), в том числе и в тех, где заданы значения погонной циркуляции вихревого слоя (zA = Zi), за исключением точек, расположенных на концах крыла (гА = г1 =-1, гА =1). Если точка гАне совпадает ни с одной из границ вихревых полос Zi, то расчет сводится к суммированию индукции от всех полосок вихревой пелены крыла

_ _ j N-1

Vyo(zA) = —X 4л U

(У,

- ч , Yi+l(zA -Zi)-Yi(ZA — Zi+l) ^ 1+1 Ii)+ - - 111

— —

za - Zi+l

za - Zi

ZA ^ Zi .

(17)

Zi+l - Zi

Если точка zA совпадает с какой-либо из внутренних границ вихревых полос zk (1 < k < N), то расчет сводится к суммированию по формуле (15) индукции от всех полосок вихревой пелены, за исключением двух, примыкающих к zk, и с учетом индуктивного воздействия двух последних по формуле (16)

vy0 (zk ^ = ~~~

4п

k-2

22

i=l

(Y1+1-Y1H^+i(z^zi^y,(z^z^l)ln

zi+l - zi

zk -zi+l

zk zi

n-I ik1

zk zk+l

n +

zk zk—l

(2 < k < N - 1). (18)

i 1 i zi 1 zi

zk -zi+l

zk zi

Определение закона распределения циркуляции по размаху крыла

Расчет поля индуктивных скоростей сводится к выбору рационального закона распределения дискретных сечений Zi по размаху крыла и вычислению в этих точках значений погонной

циркуляции вихревого слоя у;(г;)на основе известных циркуляций Г(г;). В большинстве случаев законы распределения циркуляции по размаху крыла близки к эллиптическим, для корректного описания которых распределение расчетных точек должно сгущаться к концам крыла. Оптимальное их количество должно составлять 8 - 10 на полуразмах [1; 2].

Вычисление значений у; (г;) по известным циркуляциям Г; (г;) согласно (7) связано с численным дифференцированием зависимости Г(г) в точках г = . Применим способ, предложенный А.В. Чесаловым [6], и опишем закон Г^) на участке [ Zi, Zi+1 ] полиномом

Г'(г) = 0,5Аг2 +Вг + С, ^ (19)

Производная Г'(г) будет линейной функцией от г , что соответствует допущению (8), принятому при выводе формул индуктивной скорости от вихревой полоски

Г'(г) = у (г) = Аг + В, Zi < г < Zi+1. (20)

Коэффициенты А, В, С найдем из условий, что парабола проходит через две заданные точки Г; (г;), Г;+1 (г;+1), а в первой из них имеет заданную производную Г[ ):

Г(^) = Г1; Г(^+1) = Г1+1; Г(^) = Г|. (21)

Определение коэффициентов А, В, С по условиям (21) сводится к системе линейных уравнений, решение которой имеет вид:

.(Гы-Г^-Г^-^).

А = 2-

В =-7=-—~Г2-> (22)

c = f | г;Ztzi+1 (zi+1 - i,) + z\ (Г 1+1 - Г,)

1 _ (z+l-z)2 '

Далее, по (20) можно найти значение r'(zi+1) на конце отрезка [ z¡, zi+1 ]

— Г - Г —

r'(zi+1) = y(zi+1) = Azi+1 + В = 2 _1+1 Г;, (23)

Z+i"Z

которое может быть принято в качестве начального условия для построения аппроксимирующего полинома на следующем отрезке [ z¿, zi+J. Тем самым обеспечивается гладкость сплайна

r(z) во всем диапазоне задания точек Г (z) •

Идея описанного приема аппроксимации закона T(z) заключается в том, что при обтекании крыла невозмущенным потоком без скольжения этот закон всегда является симметричной функцией относительно корневого сечения крыла. Поэтому производная P(z) в этом сечении равна нулю, что позволяет построить аппроксимирующие сплайны Г (z), двигаясь влево и вправо от корневого сечения крыла.

Таким образом, алгоритм определения закона распределения циркуляции по размаху крыла сводится к последовательным приближениям. Переход к следующему приближению выполняется одновременно для всех сечений крыла, используя общий коэффициент релаксации ц

f (v+D (- } = f (V) (- } + Kf } _ f (V) (z¡)), 2 < i < N - 1. (24)

Для устойчивой апериодической сходимости на всех углах атаки крыла, включая закрити-ческие, коэффициент релаксации должен лежать в диапазоне 0,05 < ц < 0,2. Описанный алгоритм последовательных приближений устойчиво сходится как в линейной области зависимости cyp= cyp(a), так и на закритических углах атаки.

В качестве примера приведем результаты расчетов характеристик прямого крыла X = 5, П = 1 на базе профиля NACA 0012 в диапазоне углов установки 5° < ф <30°. Как видно из рис. 3, зависимости T(z) и vy(z) являются достаточно гладкими функциями, несмотря на то что по размаху крыла взято всего 17 сечений.

В соответствии с формулами (4) коэффициенты подъемной силы и сопротивления крыла

вычисляются интегрированием по его размаху законов с (z), с (z), b(z), vy (z), W(z) . Поскольку эти законы заданы значениями функций в дискретных точках z¿, вычисление обоих интегралов может быть выполнено методом Симпсона, полагая, что значения подынтегральных функций на концах крыла равны нулю.

На рис. 4 показаны зависимости коэффициента подъемной силы крыла сук(ф, X) и коэффициента сопротивления крыла Схк(ф, X) в диапазоне изменения угла установки -30° < ф < 30° для различных удлинений крыла X = 3, 5, 7, 10. Видно, что качественный характер этих зависимостей полностью соответствует общепринятым представлениям теории крыла конечного размаха [1 - 4; 6; 8].

ГОО

Уу(г)

,0.25 ■

I I I I I 0

045-

■01

0 05

0

Гй

-ф=5°

Ф= 10° Ф= 15° ф=20° ф=25° ф=30°

-1. 0 -0. 8 -0.6 -0. 4 -0. 2 0. 0 0. 2 0. 4 0. 6 0. 8 1. 0

6 -0 4 -0 з-её" 00. 4 и -0

|/

\

\

/ к

\

у V

1

-0 5 - *

У1уй

-ф=5° ф=10° ф=15° ф=20° ф=25° ф=30°

Рис. 3. Законы изменения циркуляции и индуктивной скорости по размаху крыла

-1

0

Сук(ф,

Схк(ф, X)

Рис. 4. Аэродинамические характеристики крыльев различных удлинений

Выводы

1. Математическая модель крыла, построенная на основе вихревой пелены с непрерывно распределенной по ее участкам погонной циркуляцией, дает возможность вычислять индуктивные скорости в произвольной точке пространства, непосредственно на вихревой пелене и на поверхности, с которой сходит вихревой слой.

2. Алгоритм вычисления компонентов индуктивной скорости, сведенный к элементарным функциям, позволяет строить на его основе итерационные методы расчета воздушной нагрузки, распределенной по размаху несущей поверхности крыла, и его интегральные аэродинамические характеристики на режимах как досрывного, так и засрывного обтекания.

3. Работоспособность алгоритма продемонстрирована в широком диапазоне углов атаки профиля (-5° < а < 25°) и чисел Маха (0,3 < М < 0,9), а достоверность подтверждена сопоставлением с результатами ЦАГИ.

ЛИТЕРАТУРА

1. Юрьев Б.Н. Крыло с изменяющейся по размаху циркуляцией: избранные труды. - М.: Изд-во АН СССР, 1961. - Т. 2.

2. Бураго Г.Ф. Опыт расчета распределения подъемной силы по размаху крыла на закритических углах атаки // Труды ВВИА им. Н.Е. Жуковского. - 1944.

3. Майкапар Г.И. К аэродинамическому расчету крыла конечного удлинения // Труды ЦАГИ. - 1944. - № 547.

4. Каменков Г.В. Теория крыла в закритической области // Труды КАИ. - 1947. - Вып. 18.

5. Сиразетдинов Т.К. Крыло конечного размаха на больших углах атаки // Труды КАИ. - 1956. - Вып. 31.

6. Чесалов А.В. Построение поляры Лилиенталя монопланного крыла произвольной формы // Труды ЦАГИ. - 1929. - № 342.

7. Лойцзянский Л.Г. Механика жидкости и газа. - М.: Наука, 1973.

8. Карафоли Е. Аэродинамика крыла самолета // Несжимаемая жидкость. - М.: Изд-во АН СССР, 1956.

9. Миль М.Л. и др. Вертолеты // Расчет и проектирование // Аэродинамика. - М.: Машиностроение, 1966. - Т. 1.

VORTICAL MODEL OF THE WING COVERED WITH CONTINUOUSLY DISTRIBUTED CIRCULATION OF THE VORTICAL LAYER

Artamonov B.L.

The linear vortical model ot the final scope of a wing is exsamined. It representis the flat rectangular spatial veil covered with continuously distributed vortical layer. Elements of digitization of a veil are the quadrangular panels laying on its surface. Method, algorithms and the program of calculation of three making vectors of inductive speed from any guided rectangular platform covered with a vortical layer are created. Its intensity linearly changes on the surface of a platform. The decision is received in elementary functions. The numerical way solves the task of a definition of the law of circulation of the attached whirlwinds in scope of a wing and calculation of its aerodynamic characteristics, being based on the accepted vortical model and a hypothesis of flat sections.

Key words: wing, circulation, vortical layer, linear vortical model.

Сведения об авторе

Артамонов Борис Лейзерович, 1947 г.р., окончил МАИ (1972), старший научный сотрудник, кандидат технических наук, заместитель заведующего кафедрой проектирования вертолетов МАИ (НИУ), автор более 180 научных работ, область научных интересов - аэромеханика винтокрылых летательных аппаратов вертикального взлета и посадки.