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

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

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

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

Когда число Мт е лишь незначительно превышает М 1 ф , сверх­ звуковая зона узка и местный скачок располагается сразу за кри­

тическим сечением. В этих условиях давление повышается

лишь

на небольшом участке верхней поверхности крыла. На остальной

части верхней поверхности давление уменьшается. Кроме того,

происходит повышение давления под крылом. В результате

пло­

щадь эпюры р(х) увеличивается

и коэффициент подъемной силы су

возрастает, обычно

даже круче,

чем при Мт о

< М К р .

При дальнейшем

увеличении

участок

верхней поверхно­

сти, на котором давление повышается, быстро расширяется, а уча­ сток, где возрастают разрежения, сокращается. Это приводит к быстрому замедлению роста, а затем к интенсивному падению производной с*. С приближением Мт е к единице в зависимости от

формы профиля и угла атаки течение функции с* (М„) может быть различным.

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

ном диапазоне

чисел

Мо т

характерно заметное увеличение

про­

изводной с*

Оно объясняется тем, что за счет более заднего

расположения

сечения

[ т \ п

и очень плавного последующего

рас­

ширения струйки нижний скачок раньше, чем верхний, прибли­

жается к хвостику

профиля.

При этом на всей нижней поверх­

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

числа

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

как над крылом между верхним скачком и задней кромкой еще

остается

участок

поверхности, на

котором увеличиваются разре­

жения. При больших углах атаки

воздух под крылом

заторможен

сильнее,

нижний

скачок

образуется

при большем числе

и вы­

ходит на кромку после верхнего

скачка. В этом случае увели­

чения

производной - С,

естественно, нет. С увеличением числа Мт е

свыше

М к р основная часть аэродинамической нагрузки

смещается

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

Зависимость су(а) при смешанном обтекании не является строго линейной, так как каждому углу атаки соответствует свое, частное значение критического числа М. Так, если основным кри­

тическим

числом

М некоторого профиля считать значение МК р при

а = 0, то

при М т е )

несколько меньшем этого значения, обтекание

профиля на малых углах атаки будет дозвуковым. При некотором

угле а > 0

это же

значение

М м окажется

критическим. При не­

сколько большем

угле

атаки

выбранное число Мт е будет

закри-

тическим

и Су за счет

сжимаемости воздуха

дополнительно

возра-

90

стет. При еще больших углах атаки

то же самое

число Мх

 

попа­

дает

в

область

интенсивного падения су, которое

частично

ней­

трализует повышение су

за счет увеличения

угла атаки,

и

наклон

кривой

Су (а)

уменьшится. Однако в

ограниченной

области

углов

атаки,

которая

обычно

используется

 

в полете

при

Моо,

близких

к М к р ,

можно

пренебречь указанным-

обстоятельством

и

считать

зависимость су(а)

при любом фиксированном

числе М м

линейной.

Тогда

производная

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

числа

М м

будет изменяться

так же, как сам

коэффициент су

при

а = const.

 

 

 

 

 

Критический

угол атаки и максимальное значение

 

коэффи­

циента подъемной силы с увеличением числа

 

свыше

М к р

сна­

чала

интенсивно

уменьшаются,

затем

их уменьшение

замедляется,

а при

Мое, близких

к единице,

они

могут даже

несколько возра­

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

положены

вблизи критического сечения. В

этих условиях разви­

тие срыва

практически

сразу

приводит к падению су. С

дальней­

шим

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

Мо т

скачок и совпадающая с его основа­

нием

точка

отрыва постепенно смещаются

к хвостику

профиля,

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

Качественные изменения аэродинамических характеристик про­ филя, обусловленные образованием и развитием местных сверх­

звуковых зон и скачков

уплотнения, называют

в о л н о в ы м

к р и ­

з и с о м

п р о ф и л я ,

а

диапазон чисел Мсо

от

М к р

до М<х>=1,

в котором эти

изменения

происходят,— о б л а с т ь ю

в о л н о ­

в о г о

к р и з и с а .

Наиболее

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

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

явле­

ниями

следует

считать

интенсивное увеличение

коэффициента сх,

падение коэффициента су и сдвиг аэродинамического фокуса про­ филя назад.

§ ЗЛО. Аэродинамические характеристики профиля при сверхзвуковом обтекании

i

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

91

тем, как профиль деформирует линии тока в обтекающем его воз­

душном потоке.

 

 

В сверхзвуковом

потоке изменения

скоростей и давлений мож­

но связать с углами

|3 поворота потока

(см. § 1.15).

В силу непроницаемости крыла вектор скорости воздуха около

произвольной точки

А его поверхности (рис. 3.28) направлен по

касательной к дужке профиля. Следовательно, угол 8 поворота линий тока относительно направления невозмущенного потока в

Рис. 3.28. К определению давления в произвольной точке профиля

этой точке складывается из угла атаки а и угла 8, заключенного между касательной к дужке профиля и хордой:

Р == я - j - 8.

Считая углы а и 8 достаточно малыми, а число Moo достаточно большим, определим избыточное давление в произвольной точке профиля по формуле линейной теории (1.42):

д

2 (а + I)

 

Р

^ ос Г ""'о

"

 

К м! - 1

 

Соответственно коэффициент давления в произвольной точке

будет

 

 

 

 

 

-

_

кр_

2S

(3.22)

 

 

 

 

 

Как видно,

распределенная

нагрузка

профиля в сверхзвуке-"

вом потоке делится

на две части. Первый

член уравнения (3.22)

выражает несущую нагрузку, обусловленную поворотом потока на угол а, а второй член — начальную нагрузку и зависит от формы профиля (при а = 0). В рамках линейной теории никакой взаимной связи между этими нагрузками нет, так как несущая нагрузка не зависит от формы профиля и ее можно рассмотреть на про­ стейшем профиле — на тонкой плоской пластинке АВ бесконечного размаха, установленной под небольшим углом а в сверхзвуковом потоке (рис. 3.29).

Передняя кромка пластинки разделяет поток на две части, об­ текающие ее сверху и снизу. Верхняя часть потока, входя на пла­ стинку, поворачивается на угол а, расширяется и разгоняется в секторе между волнами разрежения' АК\ и AKj. За волной АК2

92

воздух движется параллельно пластинке. Коэффициент давления здесь не меняется вдоль хорды и во всех точках верхней поверх-

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

Км ! _ i

Узадней кромки поток поворачивается на угол —а (во вну­ треннюю сторону) на хвостовом скачке уплотнения BL\. При этом (по линейной теории, пренебрегающей потерями механической энергии на скачках) восстанавли­ ваются параметры невозмущенного потока.

Нижняя часть потока около пе­

 

 

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

поворачивается

на

 

 

внутренний

угол — а,

сужается

и

 

 

тормозится

на

головном

скачке

 

 

уплотнения

АЬ2.

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

ко­

 

 

эффициент давления на нижней по­

 

ВП

верхности

пластинки

будет

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

о

 

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

исходного

на­

 

нп

правления

и исходных

значений

па­

 

раметров нижнего

потока

осущест­

Рис. 3.29. Картина сверхзвукового

обтекания и

аэродинамические

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

силы

пластинки

в секторе

КзВКл-

 

 

 

 

 

 

Таким образом, эпюра несущей нагрузки пластинки и любого тонкого профиля при сверхзвуковом обтекании имеет вид прямо­

угольника, симметричного относительно оси Ох.

Коэффициент

подъемной силы профиля

в сверхзвуковом потоке

будет

су=* j

(P*-PB)dx--

(3.23)

6

 

 

Поскольку несущая нагрузка распределена вдоль хорды рав­ номерно, то аэродинамический фокус F любого профиля (точка приложения равнодействующей несущей нагрузки) расположен на средине хорды:

3^ = 0,5.

(3.24)

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

93

Обратим внимание на тот факт, что это сопротивление обуслов­ лено несущей нагрузкой- и непосредственно связано с подъемной силой. При дозвуковом обтекании крыла бесконечного размаха

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

нет. Сопротивление QB j

называют в о л н о ­

в ы м и н д у к т и в н ы м

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

Его энергетиче­

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

Коэффициент волнового индуктивного сопротивления профиля

найдем, разделив

выражение

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

Sqx:

 

 

 

 

 

 

 

 

 

Суа-.

 

4 а 2

 

 

 

 

 

 

 

(3.25)

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Начальной нагрузки

(при а = ао) у тонкой

пластинки нет, так

как ее поверхности совпадают с хордой

и для

любой

точки

угол

8 = 0.

Но

такая нагрузка обязательно

имеется

у

любого профиля,

 

 

^

 

 

 

имеющего конечную

толщину.

 

 

 

 

 

 

 

Рассмотрим,

 

например,

клино­

-

3

^

 

 

 

видный

профиль, образованный тре­

 

 

1

b

Г |

 

мя

плоскими

гранями

(рис.

3.30).

 

 

 

 

 

 

При

а = 0

нижняя

 

поверхность

кли­

 

 

 

 

 

 

на,

как

и

пластинка,

не создает в

 

 

 

 

 

 

потоке

возмущений

и на

ней р н = 0.

 

 

 

 

 

I

При

входе

на передний

скат

 

верх-

 

 

 

 

 

ней

поверхности

поток

поворачи­

 

 

 

 

вп

 

вается

на угол §п . Соответственно

 

 

 

 

на этом

скате

 

 

 

 

 

 

 

О НП нп

 

 

X

 

 

Ра=о-

Ум*й

- 1

 

 

 

 

 

вп

 

 

 

При переходе на задний скат по­

Рис. 3.30.

Клиновидный профиль

в

сверхзвуковом

потоке

 

ток

поворачивается

на

угол

8П во

 

 

 

 

 

 

внешнюю сторону

(при этом

восста­

навливаются параметры невозмущенного потока) и еще в

ту же

сторону на угол Ь3. Поэтому на заднем

скате

 

 

 

 

 

 

 

 

 

 

 

Р*=я

=

28,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выражая

максимальную

толщину, с профиля

через углы

накло­

на переднего

и заднего

скатов:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с = х1Ьв = (Ь —

х1)Ь„,

 

 

 

 

 

 

 

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

1л. — Ь — хх

Это означает, что положительная и отрицательная площади эпюры р(х) равновелики и что, следовательно, начальная нагруз-

94

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

щем виде этот вывод сводится к тому; что в сверхзвуковом

потоке

у любого профиля угол

нулевой подъемной силы 'а0

= 0. Как и при

дозвуковом обтекании

несимметричного профиля,

эпюра

началь­

ной нагрузки образует

восьмерку: передняя часть

клина

создает

отрицательную подъемную силу, а задняя — положительную. По­ этому момент при нулевой подъемной силе направлен на пикиро­ вание и при а > 0 центр давления расположен позади аэродинами­ ческого фокуса.

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

ния переднего

и

заднего

участков

приложены

на

срединах

этих

 

 

 

 

 

 

 

 

 

 

 

 

 

х, ,

Ь — х,

 

b

участков и, следовательно, плечо между ними — -\

^=

" у •

Сама

равнодействующая

сила

для

переднего ската

(для

заднего

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

только

знаком)

будет

 

 

 

 

 

 

 

 

 

 

Г 0 п

= -

, j g L _

qJXl.

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

м

_

v

b

 

 

 

ZKqJxib

 

 

 

 

 

 

 

 

^zO

* On " у — —2 "J/MI

 

 

 

 

 

Для перехода к безразмерному

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

разделим

мо­

мент

Mz0

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

q^Sb

=

qjb2

и одновременно

учтем,

что

8 „ =

с

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

л

х

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тг0

= хЛ

cqjxxb

 

 

 

 

 

 

 

 

 

 

 

Wt„-\qjb*

 

 

I

М2

 

 

 

 

Поскольку

при

а = 0

давление

на

переднем

скате больше,

чем

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

Разделив Q„o на произведение qтс5 = q Jb, найдем коэффи­ циент волнового сопротивления ромбовидного профиля при нуле­ вой подъемной силе:

95

Коэффициенты mzQ и с х в 0 зависят от распределения давления по профилю. Поэтому для других' профилей их выражения, есте­ ственно, отличаются от формул, полученных для клина. Но в лю­

бом случае избыточные

давления

остаются пропорциональными

углам поворота потока и обратно

пропорциональными величине

УМ.2^ — 1. Поскольку для

каждого

типа профиля средние углы

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

носительную толщину

с и ее относительную координату х и то коэф­

фициенты т г 0 и с х в 0

в общем виде выражаются следующим об­

разом:

 

 

 

C " ° = F § = T

 

( 3 ' 2 7 )

Здесь

(i(c, Х\) и

f2(с, *i) — функции с

и Х\,

выражения

кото­

рых зависят от типа

профиля.

~

 

 

 

Рассмотрим более подробно картину сверхзвукового обтекания

профиля

с криволинейными

обводами при

а > 0

(рис. 3.31).

На

верхней поверхности у передней кромки углы а и 3 имеют проти­ воположные знаки. Если' угол атаки меньше, чем модуль угла на­ клона касательной непосредственно у носка, то здесь образуется слабый верхний головной скачок уплотнения, на котором поток по­ ворачивается параллельно этой касательной и несколько тормо­ зится. Далее по всей верхней дужке профиля происходит плавное расширение струек, увеличение скорости и падение давления. В не­

которой промежуточной точке N касательная к профилю

парал­

лельна вектору

У». Здесь р = 0, скорость и давление имеют

значе­

ния Voo и рж. На

задней кромке углы а и 8 складываются. Поэтому

верхний хвостовой скачок за криволинейным профилем всегда мощ­

нее, чем за пластинкой при том же угле атаки.

 

Если угол атаки равен модулю угла Ь у передней кромки, то

верхнего головного

скачка нет.

При

входе на верхнюю

поверх­

ность крыла поток

не изменяет

своих

параметров — точка

N сов­

пала с носком профиля. Если угол атаки еще больше, то при входе

на крыло происходит поворот потока во

внешнюю

сторону.

На нижней поверхности около носка

углы а

и

5 всегда скла­

дываются. Поэтому нижний головной

скачок

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

всегда интенсивнее, чем перед пластинкой. За скачком поток не­

прерывно расширяется и разгоняется. При малых углах

атаки

здесь имеется точка N, в которой параметры равны своим

зна­

чениям в невозмущенном потоке. При таких углах а на задней кромке образуется слабый нижний хвостовой скачок уплотнения. С увеличением угла атаки точка Л/ смещается к задней кромке. Когда она совпадет с задней кромкой, нижний хвостовой скачок исчезнет — воздух плавно выходит из-под крыла, имея параметры

96

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

Если профиль имеет закругленный носок (рис. 3.32), то, как указывалось в § 3.8, его передняя часть обтекается с дозвуковой скоростью. Непосредственно около носка имеется точка О полного торможения. Однако давление рог в ней будет меньше, чем оно было бы в случае плавного торможения воздуха от М м до М = 0, за счет потерь механической энергии на головном скачке уплот­ нения:

А>2 = ° г Л . О + 0.2М1)3 '5 ,

где d — коэффициент восстановления давления на головном скачке.

Рис. 3.31. Влияние угла атаки

Рис. 3.32. Профиль с закругленным

на картину сверхзвукового об-

носком в сверхзвуковом потоке

текания профиля

 

Коэффициент давления

в точке О соответственно

будет

 

_

 

сг (1 + 0 , 2 M J ) 3 - 5 - l

 

 

/'max—

Q J J ^

 

Например, при

Мт е = 2

аг =0,72 и pm ax=l,66.

 

За точкой О струйка, обтекающая верхнюю дужку профиля,

быстро сужается и

в

некоторой

точке Л, которая с

увеличением

M M приближается

к

носку, ее

сечение /в .кр становится критиче­

ским. Давление в этой

точке

 

 

 

 

АсР =

0,528/>0 2

 

4-831

 

 

 

 

97

и коэффициент давления

О.528ог (1 + 0 , 2 М ^ ) 3 ' 5 — 1

Например, при Моо = 2 р1,07. За критическим сечением струй­ ка непрерывно расширяется, скорость возрастает и давление па­ дает вплоть до хвостового скачка уплотнения.

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

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

может

протекать либо

одинаково круто

(тогда эпюры Рп(х) и

рв(х)

эквидистантны),

либо в несколько

различном темпе. За счет

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

закругленным

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

рями механической энергии потока и приводит к увеличению коэф<

фициента с х в

0 .

Что касается применимости формул линейной теории для ре­ шения практических задач, то необходимо отметить, что при опре­ делении суммарных характеристик профиля и, в частности, коэф­ фициентов cv, сХВг, с х ъ й , эти формулы дают значительно мень« шие ошибки, чем при определении местных скоростей и коэффи­ циентов давления в конкретных точках. Дело в том, что неучет изменений числа М в процессе поворота потока во внешнюю сто­ рону и потерь механической энергии на скачках уплотнения дает

ошибки в одну и ту же

сторону — занижает

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

дав­

ления. При

этом эпюры

р(х) смещаются

в сторону

больших

раз­

режений, но их площади искажаются сравнительно мало.

 

 

Анализируя

все сказанное,

приходим

к

следующим

выводам.

1. Коэффициент подъемной силы тонкого

профиля не

зависит

от его формы. Зависимость су(а)

при любом

фиксированном числе

М е т > 1

остается

линейной. Производная

с*

в

сверхзвуковом

по­

токе плавно уменьшается с увеличением

числа

М м .

Как

следует

из формулы

(3.23):

 

 

 

 

 

 

 

 

2. Аэродинамический

фокус

любого

тонкого

профиля

при

М о о > 1

расположен на

средине

хорды

(.v> = 0,5).

Угол

нулевой

98

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

3. Коэффициент лобового сопротивления профиля при сверх­ звуковом обтекании складывается из трёх частей:

сх — Схтр + СхвО + Схв I-

(3.29)

Коэффициент волнового сопротивления профиля при нулевой подъемной силе отражает потери механической энергии на скачке

уплотнения

при а = 0. Для

каждого данного типа профиля он тем

больше,

чем

больше

относитель-

C i

 

ная толщина, так как с ее увели­

 

 

чением

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

 

средние

 

 

углы наклона скатов профиля, а

 

 

следовательно,

и

интенсивность

 

 

скачков

уплотнения.

Для

разно­

 

 

типных профилей

при

одинаковой

 

 

относительной

толщине

коэффи­

 

 

циент с* в о тем

больше, чем

боль­

 

 

ше неравномерность

распределе­

 

 

ния углов наклона дужек вдоль

 

 

хорды. Эта неравномерность

при­

 

 

водит к

увеличению

углов

8 не­

2 5

ма

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

около

передней

и

 

 

задней кромок,

а

следовательно,

 

 

и к усилению

скачков

уплотне­

 

 

ния. Закругление

передней кром­

 

 

ки вызывает

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

увели­

 

 

чение коэффициента

с х в 0 .

С

уве­

10 1.5 2.0 2,5

Ма

личением

числа Мм

он

 

плавно

уменьшается, так как уменьшают­

Рис. 3.33. Аэродинамические характе­

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

ристики профиля

 

ханической

энергии

потока

на

 

 

единицу скоростного

напора.

 

 

 

 

Коэффициент волнового индуктивного сопротивления отражает дополнительные потери энергии, обусловленные повышением ин­ тенсивности скачков уплотнения с увеличением угла атаки. Он не зависит от формы профиля и быстро (пропорционально а2 ) воз­

растает при

увеличении а.

С

увеличением

числа Мт е при a = const

коэффициент

сХъх

плавно

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

 

Кривые изменения основных аэродинамических характеристик

профиля

во

всем

диапазоне

практически

возможных чисел М м

показаны

на

рис.

3.33.

 

 

 

4*

 

 

 

 

 

99

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