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

книги из ГПНТБ / Алтухов В.А. Основы аэродинамики летательных аппаратов

.pdf
Скачиваний:
123
Добавлен:
30.10.2023
Размер:
16.73 Mб
Скачать

Ргины хт и силы сопротивления турбулентного участка на длине пластины х т).

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

Х = Х Т — Хт' + Х я\

Здесь Ху — сопротивление пластины на участке хт при турбу­ лентном пограничном слое;

X / — сопротивление того же участка при ламинарном пограничном слое.

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

РУ2

Ху + с'/л

рУ2

-*т-

I

2

2

 

 

'Отсюда найдем коэффициент сопротивления пластины При сме­ шанном пограничном слое

 

Су Суу (^/т

Сул) Ху ,

(4.41)

Ху

 

 

 

где л;т =

—2-.

 

 

 

 

b

 

 

 

Принимая во внимание равенство (4.9), перепишем последнее

выражение в виде

 

Reкр т

 

 

 

Cf — Cf T ~ [ Cfy

(4.42)

 

I

Сул)

 

 

Re

 

140

Так как c'fr и с'^я есть функция только критического числа ReKpTV

то для данного ReKpT коэффициент сопротивления пластины при смешанном пограничном слое будет зависеть только от числа Re пластины.

На фиг. 4.14 показана эта зависимость для двух критических чисел Рейнольдса: ReKpT = 485000 и ReKpT= 2 • 10s. Пунктирные кривые на фиг. 4.14 даны для значений cf ламинарного и турбу­ лентного слоя. Для практического использования часто бывает

2Cf

удобно пользоваться формулой (4.42), графическая интерпрета-.

ция которой представлена на фиг. 4.15. Заметим, что хт=■ О соответствует полностью турбулентному слою на пластине,,

х т= \ — полностью ламинарному.

§8. ВЛИЯНИЕ СЖИМАЕМОСТИ (ЧИСЛА 34) НА ТОЛЩИНУ ПОГРАНИЧНОГО СЛОЯ

ИКОЭФФИЦИЕНТ ТРЕНИЯ ПЛОСКОЙ ПЛАСТИНЫ

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

/

р

V

х \ а

.(4.43)

5Н=■ Ах Rea = Ах I

°°

°° • ) .

\

 

^СО

I

 

Формулы (4.25) и (4.38) для местных и (4.27) и (4.39.) для средник коэффициентов трения для несжимаемого потока в об-, щем виде можно представить .так:

( р

V

х \ а

(4.44.)

c/H = flRe« = £ ( ^

- =

- .

\

Г" оо

/

 

141

Для напряжения трения

на стенке можно записать:

 

=

гр

ОО V

СО2

п Гр«О V 2СО I/ грОО

V СО

х \\ п

Г А А С \

 

 

 

 

 

J .

(4.45)

Значение коэффициентов А, В и а приведены в табл. 4.1.

 

 

 

Т а б л и ц а 4.1

 

А

R *

D *

а

 

-°м

^ср

Ламинарный

4,64

0,664

1,328

- 0 ,5

Турбулентный

0,037

0,0576

0,072

— 0,2

*Вн — для местного коэффициента трения; Вср — для среднего коэффициента трения.

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

ности температура растет, а плотность падает.

и

опытными

Сравнение с имеющимися точными решениями

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

толщины

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

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

с учетом сжимаемости

можно приближенно пользоваться теми

же

формулами,

что и

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

если в эти формулы

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

и коэффи­

циента вязкости

',1.

 

 

р* и р.*.

Тогда для

Обозначим эти средние значения через

газа с учетом сжимаемости

формулы

(4.43) -ь (4Д5)

можно

представить в аналогичном виде:

 

 

 

 

 

■о

= Ах

 

 

 

 

 

(4.46)

 

[А х V,

а

 

 

 

 

(4.47)

 

и.*

 

 

 

 

 

 

 

 

 

 

 

 

 

Р* V-

В Р* v lo ( p*xV,

 

 

(4.48)

 

С, Г оо

 

 

Г

•Сравнивая формулу (4.46) с формулой (4.43) для несжимаемого

потока, получим

 

 

 

=

У Р~°°

Р* Y'

(4.49)

8„

!\ Р-*

Р» /

 

142

При сравнении (4.45) и (4.48) найдем

(4.50)

Средняя плотность газа в пограничном слое вследствие роста температуры меньше, чем плотность невозмущенного потока (см. фиг. 4.6). Это должно повлечь за собой уменьшение напряжения трения. Но, с другой стороны, рост температуры приводит к'уве­ личению вязкости, что влечет за собой увеличение напряжения трения. Однако, как будет показано несколько ниже, уменьше­ ние плотности оказывает более сильное влияние на изменение напряжения трения, чем рост коэффициента вязкости.

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

(4.51)

Отношение плотностей и вязкостей можно выразить через со­ ответствующие отношения температур.

Действительно, имея в виду, что по сечению пограничного слоя р — const, можно записать

=Т\о

(4.52)

Р„ Т*

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

(4.53)

Подставив (4.52) и (4.53) в формулы (4.49) и (4.51), будем иметь

(4.54)

Принимая для воздуха п — 0,76, найдем: а) для ламинарного слоя = —1/2)

(4.55)

143

б) для турбулентного слоя = — 1/5)

_ 0 _ = ( 7 \ о _ у

« н — [ т* )

(4.56>

0,648

ч_

£/н

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

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

 

 

 

 

 

 

= 1 +

0,14А^ . (4.57)

О

2

i

6

8

М

Т*оо

 

 

На фиг. 4.16 дано отао-

 

 

Фиг.

4.16

 

 

 

 

 

 

 

 

шение ( -^-)в функции чис-

ла А/,

подсчитанное по формулам (4.56)

\Cfnl

и (4.57) при определя­

ющей средней температуре,

соответствующей формуле (4.57).

Из приведенных графиков видно, что для турбулентного слоя

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

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

чисел

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

интенсивнее,

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

 

 

Заметим, что приводимые в этом параграфе результаты при­ менимы для совершенного газа, т. е. газа подчиняющегося урав­

нению — = gRT, имеющего постоянную теплоемкость.

Р

§ 9. РАСЧЕТ СИЛЫ СОПРОТИВЛЕНИЯ ТРЕНИЯ ЦИЛИНДРИЧЕСКИХ ТОНКИХ ТЕЛ И ТЕЛ ВРАЩЕНИЯ

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

144

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

Ниже 'рассмотрим вопрос об определении силы сопротивле­ ния тонких цилиндрических тел и тел вращения.

А. Дозвуковые скорости

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

Величина давления в рассматриваемой точке А обтекаемого

тела

(фиг. 4.17) характеризуется коэффициентом давления

Р ~ Рсо

р =

-------- —.Примерный график распределения коэффициента

 

р

V 2

 

г оо

со

~~~2

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

Фиг. 4.18

давление невозмущенного потока) отложены вниз, а отрицатель­ ные (давление меньше, чем давление невозмущенного потока) вверх. -

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

ю Изд. № 3831.

145.

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

dp

_

Рсо v l

dp

 

d x

 

2

dx

 

 

При

наличии градиента давления

ф о) расчет погранич-

 

 

 

Vdx

)

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

настоящего учебника.

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

размеры в плане, что и рассматриваемый участок цилиндриче­

ского крыла, но с

введением

близкого к единице

' тон равеч­

ного множителя v

учитывающего влияние

толщины тела на

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

 

 

 

 

Х 1р =

2Х/ г1с.

 

(4.58)

Деля обе части равенства на произведение

Р1/2

S, получим

9

 

сх тр =

2Cf т)с .

 

(4.59)

Поправка на толщину /)с зависит от относительной толщины тела

(фиг. 4.19)

-_ с

~~ 1,

Xj.

и относительной координаты точки перехода х т= —- .

Ь

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

хода х т, либо по формуле (4.41), либо по графику (фиг. 4.15).

При таком расчете необходимо знать положение точки пере-, хода на цилиндрическом теле, которое не совпадает с положе­ нием точки перехода на плоской пластине.

146

Х тр= Ху т]с.

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

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

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

по графику фиг. 4.20,

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

ная зависимость

-

-

---------- ^

.

jct х - =

от числа Re. Здесь

 

 

 

^min

fo

 

х -

(фиг. 4.18)

координата,

где давление достигает мини-

мшп

 

 

 

 

 

мума.

В случае несимметричного обтекания верхней и нижней

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

хт ср

(4.60)

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

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

щины тела с:

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

10*

147

Р

V 2

 

со

и на пло-

мического коэффициента на скоростной напор

2

 

 

щадь наибольшего (миделевого) сечения тела

£■„. Имея это

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

тела 5 б0к, получим для

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

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

тела осевой симметрии выражение

 

 

 

Сх тр =

с/

Чс ■

( 4 . 6 1 )

Поправочный коэффициент

гк

для

тел

осевой симметрии

 

зависит от величины

 

 

 

 

 

1

I

где с — относительная толщина осе­

симметричного тела.

 

Зависимость

^(Х) представлена

 

на фиг. 4.21. Определение точки пере­

 

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

задачей, еще более сложной, чем для цилиндрических тел.

В первом приближении можно на телах вращения положение точки перехода определять так же, как и на тонких цилиндри­ ческих телах.

Б. Сверхзвуковой поток (метод местных пластин)

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

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

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

148

то получим коэффициент сопротивления трения пластины, отне­ сенный к скоростному напору невозмущенного потока

_

Q h 4 ~ Q b

-cos а.

(4.62)

£.1- тр ---

Р

V 2

 

 

 

 

г ОО

со

 

 

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

На фиг. 4.23 приведено отношение коэффициента трения пло­ ской пластины под углом атаки к коэффициенту трения плоской

Фиг. 4.23

149

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