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

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

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

Эта задача была решена Н. Е. Жуковским в его теореме о подъ­ емной силе крыла (1906 г.). Докажем эту теорему для простей­ шего частного случая: будем считать толщину и кривизну про­ филя, угол атаки и число М малыми.

Как и в предыдущем случае, представим подъемную силу элементарной полоски участка крыла в виде

dY=(pn—pB)ldx,

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

Рис. 3.7. К выводу теоремы II. Е. Жу­ ковского

ли для струек, проходящих около нижней и верхней дужек про­ филя (рис. 3.7), и вычтем одно из другого:

 

_Р~

+ Я„ =Рп

+ Ян

 

 

 

Р~

+ Я^^Ръ

+ Яв

 

 

 

 

0 = / > „ — Рв +

Яи—Чв,

 

откуда

 

 

 

 

 

PU-PB

= QB- <7Н

= - f (VI

~ VI) =

? ^ 1 Г ^

( у * ~ У * ) -

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

ковы и что, следовательно,

~ — - ж ! / ^ Тогда

 

p*-pB

= ? v j v B - v n )

 

и выражение

элементарной подъемной

силы

принимает вид

 

dY = PVJ(VB

— VH)

dx.

 

Подъемная

сила всего

участка

крыла с размахом / опреде­

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

этого выражения

вдоль

хорды:

 

 

 

/ ь

*

 

 

 

Y =

?Vj\\vBdx-\

Vtidx

 

 

 

 

о

 

 

Представим себе замкнутый контур_К! _ вплотную охватываю­

щий крыло. При малых, значениях, е, [ элемент дтого конту*

ТО

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

контуру,

-охватывающему

крыло:

 

 

 

ft

 

й

ft

о

 

 

\ \\dx~-

\ VJx

=

J' VBdx

+ f VHdx =

ф VdK.

 

о

 

b

b

.

ь

 

В силу непроницаемости крыла скорость в любой точке непо­

средственно

около

его поверхности

направлена

по касательной к

ней. Это значит, чго интеграл

ф VdK = T есть циркуляция скорости

по контуру,

охватывающему

крыло.

 

Таким

образом,

 

 

 

 

 

 

 

 

Y = 9VJT.

(ЗЛО )

Подъемная сила крыла определяется произведением плотности воздуха, скорости невозмущенного потока и размаха крыла на циркуляцию вектора скорости по контуру, охватывающему это крыло. Это и есть теорема Н. Е. Жуковского. Чем сильнее разго­ няется воздух над крылом и тормозится под ним, тем больше создаваемая им циркуляция скорости, тем больше и подъемная сила.

Сравнивая формулу (3.10) с формулой подъемной силы, полу­ ченной на основании теории подобия, можно найти связь между коэффициентом подъемной силы и циркуляцией скорости:

9V2

откуда

г

2/Г

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

§ 3.4. Зависимость коэффициента подъемной силы от угла атаки при малых числах М

При увеличении угла

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

все больше

сужаются, а

под крылом — расширяются (рис. 3.8).

Деформация

воздушных

струек около верхней

дужки

профиля

приводит к увеличению местных скоростей и понижению

давления,

а около

нижней дужки — к уменьшению местных

скоростей и по­

вышению

давления. В результате площадь эпюры распределения

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

по профилю и пропорциональный ей

коэффициент

подъемной

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

атаки

увеличи-

71

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

Угол атаки оо, при

котором ^ = 0,

называют углом атаки ну­

левой подъемной силы.

Очевидно, что

для симметричных профи­

лей ао=0. У несимметричных профилей кривизна верхней дужки

больше,

чем у нижней. При а = 0 струйки над

таким профилем

поджаты

сильнее, чем под ним, и, следовательно,

образуется

поло­

жительная подъемная сила, от которой можно избавиться,

поста-

 

ос=2°

 

 

 

 

 

 

 

—-4ZZZZZZZZ

 

 

 

а= 6°

 

у 2

'

с*= 10°

/^10°

^ < 2 .

 

Рис. 3.8. Преобразование

эпюры р(х)

с увеличе­

нием угла

атаки

 

вив крыло на отрицательный угол атаки. У применяющихся в са­

молетостроении

несимметричных

крыловых профилей

угол

нуле­

вой подъемной силы может достигать —2-.—3°.

 

 

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

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

показывают,

что при углах атаки до I2—18°

зависимость су(а) линейная. Ее

принято записывать в виде

 

 

 

 

с у = с1

( « - « о ) ,

 

(3.12)

где угловой коэффициент с*у =

-~ — частная производная

(при

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

Чем меньше

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

профиля с

и чем острее

его

носок, тем

интенсивнее

деформируются

воздушные

струйки

при изменении угла атаки и, следовательно, больше

величина

с* .

Для

предельно

тонкой пластинки

с* =

2тс =

6,28

а

для

кру­

гового цилиндра, имеющего

с~\,

с* =0. Примерная

зависимость

с [с) показана на рис. 3.9. Применяющиеся в самолетостроении

72

профили с относительной толщиной с = 4 - М6% имеют среднее значение производной

cyzz5,7 1/рад « 0 , 1 1/град.

На первый взгляд может показаться, что коэффициент подъем­ ной силы должен возрастать с увеличением угла атаки до 90°. К сожалению, это не так. При сравнительно небольших значе­ ниях а (12—18°) происходит срыв потока с верхней поверхности крыла и коэффициент су интенсивно уменьшается (см. рис. 3.13). Так как самопроизвольное уменьшение подъемной силы в полете чрезвычайно опасно, летчик должен отчетливо понимать условия образования и характер развития срыва.

о

в

1Z С,

 

 

 

Рис. 3.9.

Зависимость

 

произ­

 

 

 

водной

с* от толщины

про-

Р и с . З Л 0 . к

объяснению

срыва потока

 

филя

 

 

-

с крыла

 

При

увеличении

угла атаки

минимальное сечение

струйки око­

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

занным

сечением

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

струйке возрастают

положи-

тельные

градиенты

dP

характеризующие повышение

давления - — ,

давления на единице длины линий тока. Силы давления

тормозят

поток над крылом. Непосредственно

во внешнем потоке

это тор­

можение не вызывает каких-либо качественных изменений: повы­

шению давления до атмосферного здесь соответствует

уменьшение

скорости до скорости невозмущенного потока.

 

 

 

Принципиально

иная

картина наблюдается

в

пограничном

слое,

где воздух и без того заторможен силами

вязкого

трения.

При

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

угле

атаки положительные

градиенты

давле­

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

большими, что

в некоторой точке S, которую называют точкой

отрыва потока

(рис. 3.10), нижняя часть

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

оста­

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

за этой точкой начинает

течь в обратную

сторо­

ну— навстречу

основному потоку. Теперь к точке S воздух под­

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

73

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

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

Jlo срыва

При срыве '

х

Рис. 3.11. Перераспределение давле­

Рис. 3.12. Распределение дав­

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

ления по верхним дужкам.про-

 

филей

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

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

Угол атаки аК р, при котором коэффициент подъемной силы до­

стигает максимума, называют к р и т и ч е с к и м .

 

Величины <хкр и с у т а х зависят от формы

профиля. Для

профи­

лей с малыми значениями параметров с, хс

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

высо-

74

кие пики разрежения, а следовательно, и большие положительные

градиенты -J~ на заднем скате верхней дужки. Это понятно из

рассмотрения рис. 3.12; существенное сужение струек здесь наблю-' дается лишь на сравнительно небольшом участке вблизи носка. Естественно, что при. таком распределении давления срыв потока

происходит

при

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

небольшом

угле

атаки и

сразу

почти

на

всей

верхней

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

значе­

ния ai ! P

и Сушах

невелики.

 

 

 

 

 

 

 

 

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

ческих

параметров

изменения

сечений струек

и

давления

вдоль

хорды

протекают

более

плавно, эпюра

р(х)

полнее

и, следователь­

но,

те

же значения

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

 

 

 

 

 

су

достигаются

 

при

меньших

пи-

 

 

 

 

 

ках

разрежения

 

и градиентах

dP

 

 

 

 

 

 

-j^-.

 

 

 

 

 

У таких профилей срыв потока про­

 

 

 

 

 

исходит

при

большем

угле

атаки.

 

 

 

 

 

Так

как

точка

отрыва

в

этом

слу­

 

 

 

 

 

чае

расположена

вблизи

задней

 

 

 

 

 

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

 

 

 

 

 

тока охватывает

небольшую

часть

 

 

 

 

 

верхней дужки и обычно вызывает

 

 

 

 

 

не

падение,

а

 

лишь

 

замедление

 

 

 

 

 

увеличения

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

 

 

 

 

 

ной

силы.

По

 

мере

дальнейшего

Рис.

3.13.

Зависимость

коэффи­

увеличения

угла

атаки

точка

от­

циентов су

и сх

профиля

от угла

рыва перемещается

вперед и

коэф­

 

атаки

 

фициент Су падает.

Для

таких

про­

 

 

 

 

филей характерны

существенное различие

углов

атаки ат р ,

при

котором начинается

развитие

срыва и тряска

крыла, и аК р,

боль­

шие значения а к р

и

сутах,

плавное

изменение

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

су на

околокритических углах

атаки.

 

 

 

 

 

Величины а к р

и

Сушах зависят

также от

числа

Re. При

очень

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

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

слоя не

требуются

большие градиенты

, в

связи

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

рт\п

и при

срав­

нительно

небольшом

угле атаки.

В

результате

значения

а к р и

Су max

также

невелики.

 

 

 

 

 

 

 

 

С

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

Рейнольдса

точка

турбулизации

погра­

ничного

слоя

смещается

вверх по

потоку и

при

некотором

значе­

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

75

больше, для его остановки необходима большая

работа сил дав*,

ления. Соответственно величины а к р и сут&%

возрастают,

точка от­

рыва смещается назад, падение коэффициента су

на

закритиче-

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

 

 

 

Необходимо оговориться, что в реальных

условиях полета боль­

шим числам Re соответствуют и большие числа М. В этих усло­ виях проявление сжимаемости воздуха играет обычно решающую роль в изменении величин а к р и cy m ax. Поэтому, рассматривая изо­ лированное влияние числа Re на указанные параметры, нужно считать число М фиксированным, т. е. предполагать, что число Re увеличивается не за счет повышения скорости, а за счет какихлибо других факторов, например за счет увеличения хорды про­ филя.

§ 3.5. Коэффициент профильного сопротивления

при малых числах М

Лобовое сопротивление QP участка крыла бесконечного раз­ маха принято называть п р о ф и л ь н ы м с о п р о т и в л е н и е м . Как и лобовое сопротивление любого тела, профильное сопротив­ ление можно разделить на сопротивление трения и сопротивление давления:

Qp = = Qrp "Т" QuBBn'

Естественно,

что такую

же

структуру имеет

и безразмерный

коэффициент профильного

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

 

 

 

Qp

Qrp

Фдавл

 

 

Cjc р =

~$4~Z "

 

 

"SqZ~ ~ Сд"Гр

°х

ДаВЛ'

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

 

 

 

c x p = 2cfyk,

 

 

(3.13)

где 2cf — удвоенный

коэффициент трения

плоской

пластинки,

определяемый в зависимости от числа Рейнольдса и

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

граничного

слоя

в турбулентную;

 

 

 

 

У]с — эмпирический

коэффициент, учитывающий вихревое со­

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

На угле атаки ао

он

имеет

величину

1,06—1,10

(тем больше, чем больше относительная

тол­

щина профиля)

и с приближением

к углу

атаки

а т р

возрастает

до 1,10—1,15.

 

 

 

 

76

Таким образом, в условиях бэссрывного обтекания сопротив­

ление давления может составлять 6—15% сопротивления трения. Ранее уже указывалось, что коэффициент трения пластинки

метно возрастает при переходе пограничного

слоя

ламинарного

в турбулентное состояние. Из этого следует, что для уменьшенияза­

коэффициента профильного сопротивления

целесообразноиз

иметь

профили, на которых участок ламинарного пограничного слоя был

бы

возможно шире и сохранялся бы до возможно больших чи­

сел

Re.

С началом развития срыва потока интенсивность вихреобразо-

вания

около

профиля быстро возрастает. Поэтому, начиная с угла

атаки

аТ р,

наблюдается крутое увеличение коэффициента сх р

(рис.

3.13).

 

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

задние скаты.

 

 

 

§ 3.6.

Центр

давления

и аэродинамический фокус профиля

Точку

D пересечения полной аэродинамической

силы R с хор-

дой профиля

называют

ц е н т р о м д а в л е н и я

п р о ф и л я

(рис. 3.14). Положение центра давления на хорде принято харак­

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

координатой центра

давления

 

 

 

л-д = ^ .

 

(3.14)

Лобовое сопротивление Q направлено примерно вдоль хорды и

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

на

величину хЛ. Поэтому центр

давления

можно рассматривать

как

точку приложения

подъемной

силы У.

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

При угле атаки ао нулевой подъемной силы точка О полного торможения потока примерно совпадает с носком профиля (рис. 3.15). Обтекая крыло, воздушные струйки сужаются. И под крылом — ближе к носку, и над крылом—дальше от носка струй' ки имеют минимальные сечения, которым соответствуют два пика разрежения. Эпюра коэффициентов давления образует восьмерку, передняя петля которой создает отрицательную (Уоп), а задняя — такую же по величине положительную (Уоз) подъемную силу.

Распределенную воздушную нагрузку профиля при а = ао будем называть начальной нагрузкой. Она сводится к паре сил, момент которой Mzo называют моментом профиля при нулевой подъемной силе. У профилей с положительной кривизной он направлен в сторону уменьшения угла атаки и по общепринятому правилу зна­ ков считается отрицательным.

77

Характер начальной нагрузки и величина момента Мг0

зависят

от формы профиля и прежде всего от его относительной

кривиз­

ны /. Чем больше кривизна профиля, тем больше отрицательный

угол атаки ао, больше и расположены

дальше

друг от

друга

силы Y 0 n

и F 0 3 , больше абсолютная

величина

момента

Мгй.

В случае

симметричного профиля а0 = 0

и начальная нагрузка об­

разует не восьмерку, а единую кривую, общую для верхней и

•нижней дужек, так что момент M z

0 = 0.

 

 

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

жение струек около передней части

верхней

поверхности

крыла и

их расширение около задней части

нижней

поверхности.

Соответ-

Рис.

3.14. Центр

давления

!

Нагрузка профиля

Рис. 3.15.

 

профиля

 

 

при а = ао

ственно

к начальной

нагрузке

прибавляется

дополнительная рас­

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

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

(в основном в задней ее части). Эту

дополни­

тельную нагрузку будем называть несущей.

 

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

деформации струек при обтекании данного профиля на

угле ата­

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

деформа-'

циями струек при

увеличении угла атаки и характер ее

распреде­

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

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

несущей

нагрузки.

Эта

координата

пределах

бессрывного обтекания)

не

зависит

и от

угла

атаки,

так

как при

его увеличении

перемещение

пика разрежения

вперед

компенсируется смещением назад области положительных избы­

точных давлений

под крылом.

 

 

Очевидно,

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

несущей

нагрузки — это

подъемная сила

У. Точку F приложения

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

и называют

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

ф о к у с о м

п р о ф и л я .

78

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

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

координата фокуса

хР = - у - = 0,25.

 

 

В

начале

параграфа

было

показано, что подъемная сила про­

филя приложена в центре давления D. Теперь мы установили, что

она приложена в аэродинамическом фокусе F. Никакого противо­

речия здесь нет. Просто речь идет о двух ме­

 

 

ханически

эквивалентных

схемах изображения

 

 

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

нагрузки.

Указанную

на­

 

 

грузку можно

представить

в

виде

аэродина­

 

 

мических

сил, приложенных

в центре

давле­

D

 

ния

(рис. 3.16,а). Если же

выделить

отдель­

 

но

начальную

нагрузку

(при а = ао),

мы

по­

 

 

лучаем вторую схему: аэродинамические силы

 

 

приложены

в фокусе

профиля

и,

кроме

 

того,

а

 

на

профиль

действует

момент

М20

при

нуле­

 

ХдХр

вой

подъемной

силе

(рис. 3.16,6). Так

 

как

обе схемы эквивалентны, то моменты относи­

 

 

тельно произвольной точки, определенные на

 

 

основании

любой

из этих схем, одинаковы. На­

 

 

пример, момент

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

аэродинамиче­

 

 

ского

фокуса

по

первой

схеме — Y (хл

 

xF)

а

 

равен

моменту

Mz0

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

фокуса

по

 

d

 

второй схеме:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

м г0 •

 

 

 

 

 

 

 

 

 

Рис. 3.16.

Две схе­

 

 

 

 

 

 

 

 

 

 

 

 

 

мы изображения

Из

последнего

выражения следует,

что мо­

воздушной

на­

мент

Mz0,

который

ранее

был

назван

момен­

грузки

том

профиля

при

нулевой

подъемной

силе,

 

 

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

Y q~s

Отсюда относительная координата центра давления

х я ~ х Р - ^ .

(3.15)

Имея в виду, что величины mz0 (для данного профиля) и Хр постоянны, причем для применяющихся в самолетостроении про­ филей с положительной кривизной m z 0 < 0 , легко представить себе

характер зависимости хлу).

График такой зависимости показан

на рис. 3.17,

 

79

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