Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Цейтлин Г.М. Аэродинамика и динамика полета самолета с ТРД учебник

.pdf
Скачиваний:
103
Добавлен:
27.10.2023
Размер:
21.97 Mб
Скачать

Возвратившись к выводу уравнения (2.2-1), нетрудно устано­ вить, что правая часть записанной выше пропорции есть отношение

.

A-

D

V

произведении

ат—^-,

которые в соответствии с известным из ме­

ханики принципом Даламбера можно назвать «силами инерции», препятствующими изменению скорости элементарных воздушных частиц в двух динамически подобных явлениях. Переходя к про­ изводной пропорции

 

.„

 

2

dVt

dV2

—jdbi

u . a - 3

d b 3

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

 

 

 

 

 

 

 

IS

 

S

из наличия геометрического и кинематического подобия (j^1

j ^ - ;

V]

V2

^

и заменяя

элементарные

1

к

ли­

•jy-—

"^т/")>

отрезки нормалей

ниям

тока

dti\

и dn2

пропорциональными

им характерными

линей­

ными размерами 1\ и h,

получаем

 

 

 

 

 

 

 

 

PiVi^i

__

 

 

(26)

 

 

 

 

 

(A,

(J.2

 

 

>

Безразмерное отношение «сил инерции» к силам вязкого

трения

 

 

 

 

 

R e ~ i £ U - £

 

(2.7)

называют

ч и с л о м

Р е й н о л ь д е а.

 

 

 

Согласно уравнению

(2.6) в

явлениях,

подобных по силам

тре­

ния, числа Рейнольдса одинаковы. Повторяя рассуждения в обрат­ ной последовательности, можно убедиться, что равенство чисел Рейнольдса (Rei = Re2) является не только необходимым, но и до­ статочным условием подобия аэродинамических явлений по силам трения. Другими словами, число Рейнольдса является критерием подобия аэродинамических явлений по силам трения.

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

ших числах

Re

(соответствующих

реальным скоростям

полета)

вязкость воздуха

проявляется только

в непосредственной

близости

от

поверхности

тела. Коэффициент

сопротивления трения тела

сх?р

— ~~^

с увеличением числа

Re

уменьшается.

 

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

50

давления dP = pdS, то в соответствии с уравнением (2.2-1) можно записать

Разделяя параметры, относящиеся к натуре и модели явления,

и учитывая, что в силу геометрического подобия

= -ф-,

полу-

чаем

 

 

 

Pi

Pi

 

 

Отношение — можно определить из формулы

скорости

звука

a2 = kRT= k ~ : -£p = JL^ Тогда полученное выше равенство, вы­ ражающее условие подобия явлений по силам давления, прини­

мает вид

4 ( т. е. Mj = Ма .

«1

1

2

Таким образом, в явлениях, подобных по силам давления, чи­ сла М должны быть одинаковыми. Иначе говоря, число М являет­ ся критерием подобия аэродинамических явлений по силам дав­ ления.

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

Зависимости аэродинамических коэффициентов от критериев геометрического, кинематического и динамического подобия назы­ вают а э р о д и н а м и ч е с к и м и х а р а к т е р и с т и к а м и т е л . Исследование аэродинамических характеристик крыла и других частей самолета, а также самолета в целом — основной предмет прикладной авиационной аэродинамики.

§ 2.4. Пограничный слой при малых числах М

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

51

Пусть около такой поверхности движется воздушный поток с некоторой скоростью V (рис. 2.4). Так как кроме движения всей воздушной массы относительно тела внутри нее имеется еще хаоти­ ческое (тепловое) движение молекул, то скорость любой молекулы можно представить векторной суммой

VH = V+V.

где Ут — скорость теплового движения данной молекулы в данное мгновение.

За счет скорости VT молекулы отклоняются от общего направле­ ния движения воздушной массы и сталкиваются с неровностями стенки. Так как элементарные участки поверхности, о которые уда­

Рис. 2.4. К объяснению образования погра­ ничного слоя

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

Поскольку на поверхности тела скорость воздуха равна нулю, а на некотором расстоянии от этой поверхности поток движется со скоростью V, то в результате проявления вязкости воздуха здесь образуется переходный слой (см. § 1.5).

Слой воздуха, непосредственно прилегающий к поверхности те­ ла, в котором действуют силы вязкого трения и скорость постепен­

но нарастает

от нуля до

скорости внешнего потока, называют по­

г р а н и ч н ы м

с л о е м .

 

Общие закономерности, действующие в пограничном слое, удобнее всего рассмотреть на примере тонкой плоской пластинки, параллельной линиям тока невозмущенного потока (рис. 2.5). Око­ ло такой пластинки внешний поток во всех точках имеет одинако­ вые скорости и давления, чем исключается воздействие на погра­ ничный слой различных дополнительных факторов.

Как уже говорилось, проявление вязкости воздуха определяется

52

числом Re. Для исследования состояния пограничного слоя в раз­ личных его точках обычно пользуются местным числом Рейнольдса:

при подсчете которого в качестве характерного линейного размера используется расстояние х данной точки от передней кромки пла­ стинки.

Если число М невелико, можно считать, что плотность и темпе­ ратура воздуха в пограничном слое постоянны и, следовательно, не меняется коэффициент кинематической вязкости v. Так как ско-

X

X

Рис. 2.5. Пластинка,

параллельная ли­

Рис. 2.6. Пограничный слой пластинки

ниям

тока

 

рость на внешней границе пограничного слоя пластинки тоже по­

стоянна и равна скорости невозмущенного

потока,

то

при У» —

= const число Re* пропорционально

координате

х,

т. е. увеличи­

вается по линейному

закону

с удалением от

передней кромки.

Вблизи передней

кромки

число Rex

мало. Это значит,

что силы

трения велики по сравнению

с инерционными

силами.

Поэтому

скорость быстро нарастает вдоль нормали к поверхности и уже на

небольшом расстоянии от нее достигает

скорости внешнего пото­

ка — толщина пограничного слоя 8 здесь

мала.

Структура пограничного слоя на передней части пластинки ла­ минарная (слоистая) — перемешивания воздуха между смежными элементарными струйками практически нет.

С удалением от носка пластинки число Re* постепенно увеличи­ вается и силы трения, вызывающие изменения скорости по толщи­ не пограничного слоя, играют все меньшую роль по сравнению с инерционными силами, препятствующими этим изменениям. Следо­ вательно, нарастание скорости с удалением от поверхности замед­ ляется, а толщина пограничного слоя увеличивается.

Выделив в пограничном слое произвольную воздушную частицу (например, в виде кубика, рис. 2.6) и представив себе профиль скоростей, легко обнаружить, что сторона частицы, расположенная дальше от стенки, движется быстрее, чем сторона, расположенная

53

Рис. 2.7. Эпюры скоростей в ламинарном и турбулентном пограничных слоях

ближе к стенке. Это приводит к деформации

и вращению

частиц,

т. е. к образованию вихрей.

 

 

 

С увеличением местного

числа Рейнольдса

(с удалением

от но­

ска пластинки) количество, интенсивность и размеры вихрей воз­

растают и они начинают взаимодействовать между собой. В ко­

нечном счете это приводит к

разрушению ламинарной структуры

и переходу ее в турбулентную

(завихренную), характеризующуюся

интенсивным перемешиванием воздуха по всей толщине погранич­ ного слоя.

Протяженность

Ахт

переходного

(от ламинарного

к турбулент­

ному

состоянию)

участка пограничного

слоя

сравнительно

неве­

лика,

и при решении

практических

задач

этот участок

обычно

 

 

 

сводят

к так называемой

точке

пере­

 

 

 

хода Т.

 

 

 

 

 

 

 

 

 

Характерно, что разрушение

лами­

 

 

 

нарной

структуры

пограничного слоя

 

 

 

происходит

при определенном

значе­

 

 

 

нии местного числа Рейнольдса, кото­

 

 

 

рое называют

критическим:

 

 

 

 

 

Re,кр •

 

 

 

(2.9)

Величина критического числа Рей­ нольдса тем больше, чем лучше обра­ ботана поверхность пластинки и чем меньше турбулентность (завихрен­

ность) внешнего потока. В среднем можно считать ReK P =105 . Если известны параметры невозмущенного потока, то, зная ве­

личину ReK P , легко определить протяженность ламинарного участка пограничного слоя, равную координате точки перехода:

Минимально возможная скорость полета современного самоле­

та составляет примерно

200 км/ч = 55

м/с. Такая

скорость дости­

гается на

посадке,

т.

е.

непосредственно

у

земли,

где

v =

 

 

 

 

 

 

105.] 4 5 . Ю - 6

 

 

= 1,45-10~5

м 2 / с

Для

этого

случая хт

'-^

 

= 0,027. Из

рассмотренного

примера

видно, что даже при минимальной

ско­

рости значение координаты х т намного меньше

 

протяженности

любой из основных частей самолета

в направлении

оси

Ох. Это

значит, что в полете на большей части

внешних поверхностей

само­

лета пограничный слой имеет турбулентную структуру.

 

 

В отличие от ламинарного пограничного слоя, в котором

пере­

нос количества движения по нормали к линиям

тока

осуществляет­

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

54

непосредственно от стенки (в нижней части пограничного слоя) скорость нарастает значительно быстрее. По той же причине затор­ моженные воздушные частицы выносятся значительно дальше от стенки в воздушный поток, вследствие чего толщина турбулентного пограничного слоя больше. Так как интенсивность вихревых дви­ жений мало зависит от местного числа Рейнольдса, то увеличение темпа нарастания скорости в нижней части турбулентного погра­ ничного слоя с увеличением числа Rex происходит значительно медленнее, чем в ламинарном пограничном слое. Сравнение про­ филей скоростей в турбулентном и ламинарном пограничных слоях схематично показано на рис. 2.7.

§ 2.5. Нагрев воздуха в пограничном слое. Кинетический

нагрев конструкции летательного аппарата

Как уже говорилось, в пограничном слое воздух тормозится до полной остановки непосредственно около поверхности тела. По­ скольку кинетическая энергия воздуха при этом переходит в тепло, температура в нижней части пограничного слоя повышается.

Если считать процесс торможения воздуха в пограничном слое адиабатным, то его температуру на поверхности тела можно опре­ делить по формуле (1.23-2):

Т 0 п = Тт + 6,2Ттт.

Прирост температуры пропорционален числу М^:

&Tu = Tn-Tm=0,2rmm.

В реальных условиях главным образом за счет конвективного отвода тепла от нижней области пограничного слоя в менее на­

гретую верхнюю область и внешний поток часть тепла

рассеивает­

ся и действительное

повышение

температуры АГ

оказывается

меньшим,

чем при отсутствии

теплообмена.

Отношение

г—

 

 

 

 

 

 

 

 

 

 

 

ал

называют

к о э ф ф и ц и е н т о м

 

в о с с т а н о в л е н и я

т е м п е ­

р а т у р ы .

Как показывают

теоретические

и экспериментальные

исследования, для ламинарного пограничного слоя г=0,85

и для

турбулентного г=0,9.

Поэтому

повышение

температуры

в

нижней

части пограничного слоя можно

определять

 

по формуле

 

 

 

ЬТ = 0,2гТ

М2 в

0,1757

М 2 .

 

 

(2.10)

 

 

СС

00

'

ОС

 

ОС

 

\

/

График

этой зависимости

показан

на рис. 2.8.

 

 

 

От воздуха тепло передается конструкции летательного аппа­ рата (кинетический нагрев конструкции) и через некоторое время она прогревается до температуры нижней части пограничного слоя. В связи с этим уже в современном самолетостроении возникает ряд новых сложных технических проблем. При температурах 60—80° С, что соответствует числам М.= 1,21,3, теряет свою проч­ ность и становится непригодным для фонарей кабин органическое

55

стекло. В длительном полете с числом М = 2-г-2,5 обшивка самолета прогревается до 200—300° С и алюминиевые сплавы становятся не­ пригодными для ее изготовления. При еще больших числах М вы­ ходят на предел прочности титановые сплавы, а затем и жаропроч­ ные стали, возникает необходимость в отыскании новых, более жа­ ропрочных конструкционных материалов.

Не менее сложной оказывается задача охлаждения кабин эки­ пажа, отсеков с топливными баками и различным оборудованием, чувствительным к высоким температурам. Дело в том, что при

высоких

скоростях

полета системы охлаждения, основанные на

 

 

 

 

непосредственном

использовании

холодного

 

 

 

 

забортного воздуха, непригодны, ибо (в ре­

 

 

 

 

зультате

торможения)

в

воздухозаборни­

 

 

 

 

ках и радиаторах

воздух

нагревается

 

так

 

 

 

 

же, как и в пограничном

слое. Приходится

 

 

 

 

применять

холодильные

установки,

осно­

 

 

 

 

ванные

на

механическом

срабатывании

 

 

 

 

теплоперепада

 

(турбохолодильники)

и

на

 

 

 

 

испарении

специальных

жидкостей.

 

Та­

 

 

 

 

кие

установки

 

имеют

 

определенные

га­

 

 

 

 

бариты

и

вес,

потребляют

значительные

 

 

 

 

мощности.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Одним

из

наиболее

кардинальных

ре-

и

j

g j

Мм шений

проблем,

связанных с

кинетическим

 

 

 

 

нагревом

конструкции

летательного

аппа-

Рис. 2.8. Повышение тем-

рата,

является

увеличение высоты полета,

пературы

 

пограничного

В в е р

х н и

х

С Л О Я х атмосферы,

на

высотах

 

 

 

 

60—80 км и более, воздух сильно

разрежен,

 

 

 

 

в связи с чем коэффициент

теплопередачи

от пограничного слоя к стенке

столь

мал,

что

даже

длитель­

ный полет с числами Л4^8 - М0 не вызовет опасного нагрева

 

кон­

струкции.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Предполагается, что, применяя специальные теплоизолирующие покрытия и системы охлаждения кратковременного действия, в ко­ торых пары охлаждающей жидкости выводились бы непосредствен­ но в поток без конденсации для повторного использования, можно будет обеспечить полет с числами М ~ 4 ч - 6 и на значительно мень­ ших высотах (10—20 км). Разумеется, безопасная продолжитель­ ность такого полета строго ограничена.

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

Поскольку давление по толщине пограничного

слоя не меняется

и остается равным давлению во внешнем потоке

(для тонкой пла­

стинки— давление рк), то с повышением температуры

происхо­

дит изобарное расширение воздуха. Средние значения

расстояния

56

между молекулами и длины их свободного пробега увеличивают­ ся, то же изменение скорости происходит на большем отрезке нор­

мали, и толщина

пограничного слоя

увеличивается.

 

§ 2.6. Сопротивление трения

Разделим пограничный слой пластинки на множество предель­

но тонких слоев

(рис. 2.9). Так как

скорости в смежных слоях

различны, между ними действуют силы вязкого трения, равные

количеству

движения, переносимому

за единицу времени

из слоя

в слой. Слой / увлекает слой 2, слой

2 — слой 3 и т. д. В

конеч-

 

V

 

 

 

——Ч

'

 

 

1

* - /

2

 

 

 

—)

3

 

 

 

•У'

•••

 

 

 

п-1

 

 

 

 

п

 

 

Рис. 2.9. К

пояснению

природы сопро­

Рис. 2.10. К определению коэффи­

тивления трения

циента трения

 

ном счете силы трения замыкаются на поверхности тела в форме

непрерывной

бомбардировки

молекулами

неровностей зтой

по­

верхности.

 

 

 

 

 

 

 

Используем

введенное в

§ 2.4 представление о

скорости

мо­

лекулы как о результирующей скоростей

хаотического

движе­

ния VT и направленного движения V (VM =

Vr•+ V).

 

 

 

Если за

счет

скорости VT

вероятность

и средняя

сила

ударов

молекул на элементарных участках неровностей поверхности тела, обращенных вперед и назад, совершенно одинаковы, то наличие скорости V увеличивает вероятность и среднюю силу ударов моле­ кул на участках поверхности, обращенных вперед. В результате

суммарная сила воздействия молекул воздуха

на тело

приобре­

тает составляющую, направленную назад (в

сторону

движения

потока). Эта составляющая, численно равная секундной потере количества движения воздуха в пограничном слое, и является со­

противлением

трения.

 

 

 

Выделим на поверхности пластинки, имеющей постоянную ши­

рину / и длину b

(рис. 2.10),

элементарную площадку

dSTp=ldx.

В

соответствии

с

формулой

(1.11) сила трения, действующая на

эту

площадку,

будет

 

 

57

где ( - ^ f ) 0 H р.о — соответственно градиент скорости по нормали и коэффициент динамической вязкости воздуха непосредственно на

поверхности пластинки.

 

 

 

 

 

 

По формуле

(2.3),

полученной

на

основании

теории

подобия,

ту же силу можно представить в виде

 

 

 

 

 

 

dQrp =

c'fqjdx.

 

 

Вошедший в последнюю формулу

коэффициент

 

 

C r ^ = ? : h

 

w .

 

( 2 Л 1 )

называют м е с т н ы м

(для данного

элемента трущейся

поверхно­

сти) к о э ф ф и ц и е н т о м

т р е н и я .

 

 

 

Общее значение коэффициента трения для данного тела в дан­

ном потоке определяется

осреднением

местного

коэффициента c'f

по всей трущейся

о воздух

поверхности

ST p :

 

 

 

 

 

5

 

 

ь

 

 

При малых числах М м коэффициент динамической вязкости воздуха цо= Iх—const и изменение коэффициента c'f вдоль пла-

f

dV\

стинки определяется изменением градиента

y-j^j-

Как было показано, при ламинарном пограничном слое увели­ чение местного числа Re* приводит к более плавному нарастанию скорости по толщине пограничного слоя, т. е. к уменьшению гра-

диента

соответственно

коэффициент

cf с удалением от

передней кромки пластинки уменьшается.

 

При определении общего числа Рейнольдса,

характеризующего

обтекание

данного тела данным

потоком, в

качестве линейной

базы используется какой-либо характерный размер тела. Для пря­ моугольной пластинки таким размером является ее длина Ь в на­ правлении оси Ох. Таким образом,

R e = - = - .

Отношение местного числа Re* к общему числу Re есть отно­ сительная координата данной точки:

Re* х

Не Ь

Очевидно, что при Re<ReK p ламинарная структура погранич-* ного слоя сохраняется на всей пластинке. В этом случае выполне-

58

ние интегрирования в уравнении (2.12) (приводим без доказатель­ ства) дает результат

1.328

Если Re>ReK P , то на пластинке имеются и ламинарный, и тур­ булентный участки пограничного слоя. Турбулизация погранич­ ного слоя сопровождается интенсивным увеличением градиента

("Ж")о'

3

с л е д о в

а т е л ь

н о > и

местного коэффициента трения с'},

С

дальнейшим

увеличением

числа Re* (свыше

ReK p )

коэффи­

циент

c'j постепенно

падает,

но медленнее,

чем

в

ламинарном

слое.

 

 

 

 

 

 

 

 

 

 

 

 

При

больших

и средних

скоростях

полета

Re^Re l t p . В

этом

случае

относительная

протяженность

ламинарного

участка

хт =

 

х

Re Kp

 

 

 

 

 

 

 

 

 

=

— =

—j7^- невелика и можно считать, что весь

пограничный

слой имеет

турбулентную структуру. Тогда

(приводим

тоже без

доказательства) для осредненного коэффициента трения пластинки получается выражение

0.455

 

(2.14)

с = Ц « »

 

(lg Re)2

 

 

При малых скоростях полета, особенно

на больших

высотах,

где коэффициент кинематической вязкости

v = - t - с и л ь н о

возра­

стает за счет падения плотности воздуха, числа Re не намного превышают ReK p . В этих условиях ламинарный участок погранич­ ного слоя по своей протяженности соизмерим с турбулентным и формула (2.14) дает завышенные значения сх

При оценке влияния сжимаемости воздуха на коэффициент трения необходимо иметь в виду, что в реальных условиях полета применительно к основным частям самолета, имеющим значитель­ ную длину в направлении линий тока, большим числам М„ со­ ответствуют и большие числа Re^>ReK p . Кроме того, при боль­ ших числах Мм сопротивление трения обычно невелико по сравне­ нию с сопротивлением давления. Поэтому при оценке влияния сжимаемости воздуха в практических авиационных задачах допу­ стимо считать пограничный слой турбулентным и пользоваться

полученной для

этих условий

приближенной

формулой

 

 

Cf

~ С/нс ^ + о д е ш ^ ) 0 , 6 7 " *

 

Здесь

с / н с — коэффициент

трения

пластинки, определяемый по

формуле

(2.14)

при действительном

числе

Re

без учета сжимае­

мости воздуха (в несжимаемом потоке).

 

 

Поверхности

частей

самолета обычно

имеют небольшую кри­

визну в

направлении

линий

тока.

Поэтому

по коэффициентам

59

Соседние файлы в папке книги из ГПНТБ