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

книги из ГПНТБ / Гофман М.Л. Аэродинамика гиперзвуковых скоростей и супераэродинамика

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

Известно, что для конуса в соответствии с теорией Мунка —• Джонса нормальная сила равна

 

 

N = ± p ^ V l S 0 ‘2а .

Как показано

в работе 1124],

даже для сравнительно тонкого

конуса (с полууглом

О порядка

10°) эта формула дает результат,

отличающийся от точного ре­

 

шения. Однако результат для

 

конуса может

быть

улучшен

d ос

(и для больших углов раст­

 

вора), если вместо

скорости

 

1Л» ввести местную скорость

 

на поверхности конуса. Тогда

 

формула примет вид:

 

N =

и cos 0

Sd2a,

 

 

 

 

 

( 4 . 9

)

где и —скорость на поверх­ ности конуса,

Sd — площадь донного се­ чения.

Для расчета подъемной си­ лы тел вращения широкое при­ менение получила теория Нью-

dc,,

тона [125]. Так, величина -т—

0 2 4 6 в ю 12 м la

Рис. 4.29. Зависимость наклона коэф­ фициента нормальной силы от относи­ тельного удлинения цилиндрического тела при различных числах М

при я = 0 для

конусов, вычисленная по формуле Ньютона, очень

близка к рассчитанной по „точной" теории.

Следует, однако,

иметь в виду,

что

при

числах

М = Ю-т-15

для правильного

определения

нормальной

силы

большое значение имеет „тене­

вая" область,

в которой

угол

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

полета и нормалью

к поверхности больше

тг

; формула Ньюто­

2

на в таком случае

неприменима.

 

 

 

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

линдрического тела при различных

числах М. Видно, что при

больших удлинениях наклон кривой

(йсЛ

нормальной силы I

j

стремится к постоянной величине для каждого числа М. Вели­ чину коэффициента нормальной силы можно представить как сумму:

111

Это иллюстрируется на при Моо = 4. Заметим, что
Рис. 4.30 Кривая зависимости коэффициента нормальной силы от угла атаки для тела „конус-ци­ линдр" при числе М = 4

'

dcN

а +

,

с

do.

N ~ \

а=0

 

где Дед,— добавочная сила за .счет угла атаки.

рис. 4.30 для тела „конус-цилиндр" обтекание тела вращения под боль­ шим углом атаки вызывает отрыв в „теневой* области, и этот отрыв является причиной нелинейной за­ висимости кривой коэффициента нормальной силы cN = f(a). Нели­

нейность кривой заметна и на гра­ фике рис. 4.30.

Поясним подробнее характер об­ текания тела под углом атаки.

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

[70].

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

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

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

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

На практике все тела, летающие при гиперзвуковых скоростях, располагаются внутри головного скачка, т. е. они являются «тонки­

112

ми». В работе [67] описан численный метод расчета возмущенного поля течения, вызванного малым углом атаки осесимметричного («тонкого») тела, с помощью линеаризации относительно этого угла атаки.

Опыты показывают [38], что для тел, у которых кривизна поверх­ ности изменяется плавно, формула Ньютона в комбинации с расче­ тами по формулам Прандтля — Мейера дает приемлемую оценку распределения давления и нормальной силы.

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

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

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

обработаны

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

данные [126]

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

см и ст для

комбинации конуса

с цилиндром.

Это дало воз­

можность на основании безразмерных зависимостей определить эти коэффициенты для разных комбинаций конуса с цилиндром при углах полураствора и атаки до 15°. С помощью видоизме-

dc N

In

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

Рис. 4.31. Зависимость производной

от -q

Kt в качестве

параметра

-8- М-.-Л--Гвфм»в — .........

ИЗ

О

0,5

1.0

1.5

2.0

2.5

3.0

3.5

40 l a

 

 

 

 

 

 

 

 

X

Рис. 4.32. Зависимость производной ?Ст от X

при нулевом

угле

атаки

 

 

и при

 

d a

1п

 

 

 

 

 

Kt в качестве

параметра

 

 

Рис. 4.33. Зависимость положения центра давления от _ при нулевом

угле атаки и при Kt в качестве параметра

ненного

закона

гиперзвукового

подобия

параметры

подобия

представляются

в виде

 

 

 

 

 

 

 

 

 

 

 

d Y М2—1

Ka = V М2-

1 а.

 

 

 

 

 

 

 

 

 

 

 

Параметр Kt

зависит

от толщины (удлинения) тела, параметр

К* — от угла

атаки. В расчете установлена следующая область

изменения параметров:

 

 

 

 

 

 

 

О < Ко < 42,5,

 

 

 

 

 

 

 

0,145 <

Kt < 2,07,

 

 

 

 

 

 

0 <

у - < 4,15,

 

 

 

 

 

 

 

 

 

'’П

 

 

 

 

 

 

 

 

где

/„ — длина

конической

 

 

 

 

 

 

 

части,

цилиндра.

 

 

 

 

 

1а — длина

 

 

 

 

 

На рис. 4.31 приведена

 

 

 

 

 

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

 

 

 

 

 

нормальной

 

 

d

c

. i

 

 

 

 

 

силы —— от ~

 

 

 

 

 

 

 

 

 

 

dя

1п

 

 

 

 

 

при К , в качестве

пара­

 

 

 

 

 

метра.

При анализе обнару­

 

 

 

 

 

живается,

что

при

задан­

Рис. 4.34. Изменение величины Сд/j/M 2 —1

ном параметре подобия уве­

от

параметра

 

при Kt =

1,0

личение длины

цилиндриче­

 

 

 

 

 

 

ской

части

сначала* ведет

 

 

 

 

 

к росту

dcN

,

затем к уменьшению и снова к увеличению. Таким

 

образом, для

некоторых длин цилиндрической части нормальная

сила

у

комбинации

конуса с цилиндром может уменьшаться

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

длины

цилиндрической части. Это означает, что

при малых углах атаки имеется

область

на

цилиндрической

части

тела,

 

где нормальная сила отрицательна. На рис. 4.32

представлен

наклон

кривой

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

момента

dc , а на

рис.

 

 

 

 

 

JC

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

нулю (здесь х д—

4.33 — кривая -у-

расстояние от вершины конуса до

центра давления).

 

Рис. 4.31 —4.33

дают аэродинамические коэффициенты при

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

8*

115

мость Сдг]/м*— 1 от параметра у 1 для /(, = 1,0. Для нахожде-

ния коэффициента момента следует , использовать график рис. 4.35, на котором дана зависимость с„ ]/ М2 — 1 от парамет-

Р и с . 4.35. И з м е н е н и е в е л и ч и н ы cm Y М 2— 1

от п а р а м е т р а ^ п р и

= 1 , 0

ра ~ при Kt ~ 1,0. Аналогичные графики имеются и для дру-

!-п

гих значений К,.

§ 4.6. ТЕЛА ВРАЩЕНИЯ МИНИМАЛЬНОГО СОПРОТИВЛЕНИЯ

Вопросу выбора формы носовой части ракеты, обеспечиваю­ щей минимальное сопротивление, посвящен целый ряд работ [127, 128]. Так, используя закон Ньютона, можно найти формы тел вращения, обладающих минимальным сопротивлением при гииерзвуковых скоростях и различных исходных данных, таких, как, например, длина тела, объём и боковая поверхность тела, диаметр его основания и др.

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

тивления Ньютона р —2 sin2 6. Если известно давление на по­

116

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

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

легко определяется

как

 

 

 

 

 

 

 

cxqx~d2

 

 

 

 

(4Л0)

 

А----—

 

 

 

где / — длина тела.

 

 

в другой форме:

 

Это выражение можно представить

 

 

А? =

= Jip y /d x .

 

(4.11)

 

 

 

о

 

 

 

Подставив в это равенство

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

из фор­

мулы Ньютона и имея в виду, что

 

 

 

 

2

 

 

sin3 0 “

получим функционал в виде

 

 

 

 

 

/ «

- J t + t »*'-

 

 

(4Л2)

 

 

о

 

 

 

 

В дальнейшем

задача

сводится к

тому,

чтобы отыскать та­

кую функцию у = . у ( х ) ,

при

которой

функционал I Q имеет ми­

нимальное значение при различных заданных условиях; зная функцию у = у { х ) , можно по значениям у и х построить тело.

Результаты определения формы тел с минимальным сопротивле­ нием [129] указывают на две общие особенности очертаний этих тел: если длина тела задана, то оптимальное тело получается с тупым но­ ском, и, наоборот, если длина тела не задана, то оптимальное тело имеет острый носок. Это объясняется тем, что при фиксированной длине тела его сопротивление уменьшается при перераспределении давления таким образом, что большое давление действует на срав­ нительно малый участок поверхности вблизи носка, а на остальную, относительно большую часть тела с малым наклоном поверхности действует малое давление.

С другой стороны, если длина тела не фиксирована, то, очевид­ но, минимальное сопротивление будет иметь не тупое, а остроносое тело.

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

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

1J7

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

У

 

 

 

 

 

 

 

 

 

 

0,8

Дано:

 

 

 

 

 

""d и 5

 

0.6

 

i и S \

 

 

 

 

 

 

1и й

 

 

 

 

 

 

du V

 

0.4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.2

 

 

 

 

 

 

 

 

 

 

О

0.1

0.2

0.3

0.4

0.5

0.6

0,7

0,8

0.9

г

Рис. 4.36. Контуры тел минимального сопроти­ вления при различных заданных условиях

На рис. 4.37 показана зависимость коэффициента сопротив­ ления сХо (при су — 0) испытанных тел вращения от числа М.

PhCi 4.37. Зависимость коэффициента сопроти­ вления при нулевой подъемной силе испытан­ ных тел вращения от числа М

Здесь представлено пять моделей с удлинением ^- = 3 и три

I г „

модели с удлинением — =г 5. Контуры всех моделей в мери­

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

118

при значениях n = 1; 0,75; 0,50 и 0,25.

При п — 0,75 форма тела близка к экстремальным телам при заданных значениях длины и диаметра основания. При п = 1 получаем конус. Одна модель была выполнена в форме тан­ генциального оживала, т. е. в форме

тела, контур которого

в меридиональ­

 

 

/

ной плоскости образован дугой окруж­

 

 

ности.

 

Ньютона

была

применена

 

 

Теория

 

 

*/

для выбора оптимальной формы снаря­

0,3

Конус

да

при

других

заданных

условиях —

 

 

толщине

и

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

[130].

По­

 

 

 

казано, что форма оптимального

тела

0,2

 

 

не зависит от подъемной силы. Однако

 

 

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

связан

с

формой

и с

 

 

 

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

сх0 лобового

сопротив-

0.1

 

 

ления при су = 0-

 

 

 

 

 

 

 

 

Некоторое уточнение при выборе оп­

 

•С Оптималь-

тимальной формы тела может быть про­

 

мое тело

изведено применением формулы Бузема-

 

'

на,

учитывающей

центробежные

силы.

0

0.2

0.4

На рис. 4.38 приведена кривая коэффи­

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

циента лобового сопротивления для ко­

 

тела

 

нуса и оптимального

тела.

Из

графика

Рис. 4.38. Кривые

коэффи­

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

циента

лобового сопротив­

ния дают существенное уменьшение со­

ления

для конуса

и опти­

противления

по сравнению

с

конусами

 

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

того

же

удлинения. Так,

переход

от ко­

 

 

 

нуса к телу оптимальной формы с той же относительной толщиной уменьшает сопротивление более чем в 2 раза.

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

Расчет параметров течений при k = 1,4 представляет прак­ тический интерес, jaK как при М->-оо и &->1,0 тело вращения

степенной формы (г — сх0'15) является при заданном удлинении оптимальным по волновому сопротивлению.

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

119

канию

тела

заданной толщины

при меньших

числах М.

По­

этому результаты расчета обтекания тонких

 

тел

вращения

сте­

пенной формы при М -оо

могут быть

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

в широком

диапазоне сверхзвуковых

скоростей.

С этой

целью

по

такому

методу

было_ рассчитано

обтекание

 

тел

вращения

степенной

формы

(r =

cxn+l). Приведены результаты

 

расчета

для п = 0;

 

 

 

-0 ,1 0 ;

 

-0 ,1 5 ;

— 0,20;

-0 ,2 5 ;

 

 

 

— 0,30;

— 0,35.

 

Представленные в

 

 

 

работе таблицы позволяют опре­

 

 

 

делить

 

коэффициент

волнового

 

 

 

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

 

 

 

 

 

 

 

 

 

 

Расчеты показывают, что для

 

 

 

гиперзвуковых

 

скоростей

при

 

 

 

k =

1,4

наименьшее

волновое соп­

 

 

 

ротивление

из

тел

вращения

сте­

 

 

 

пенной формы имеет т^ело, уравне­

 

 

 

ние

которого

 

г ~ сх0,70,

а

не

 

 

 

г = ел?0,75, как

 

по

Ньютону

(при

 

 

 

k -> 1,0). Следовательно, при k = 1,4

 

 

 

оптимальная форма тела вращения

 

 

 

при

больших

числах М будет бо­

 

 

 

лее

полной,_

 

чем

ньютоновская

 

 

 

форма

г — еле0,75. На рис. 4.39

по­

 

 

 

казаны результаты расчета обтека­

 

 

 

ния

конуса

по этому методу.

 

§4.7. ВЛИЯНИЕ ФОРМЫ ГИПЕРЗВУКОВОГО АППАРАТА НА ЕГО НАГРЕВ.

АЭРОДИНАМИЧЕСКИЕ МЕТОДЫ УМЕНЬШЕНИЯ ТЕПЛОПЕРЕДАЧИ И СОПРОТИВЛЕНИЯ

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

Дело в том, что при спуске нужно погасить ту огромную кинети­ ческую энергию, которой обладает летательный аппарат, двигаясь с очень большой скоростью. Так, например, кинетическая энергия те­ ла, летящего со скоростью, близкой к первой космической скорости, настолько велика, что только одной трети ее, превращенной в тепло, достаточно для испарения тела. Однако существенное повышение температуры уже начинается при значительно меньших скоростях, соответствующих числам М порядка 4—5; нагреву прежде всего под­ вергаются передние кромки крыльев и рулей и носовые части кор­ пусов (рис. 4.40). Нагрев может привести к постепенному оплавле­ нию, а затем и к полному разрушению летательного аппарата. Проб­ леме теплопередачи при гиперзвуковых скоростях и методам умень­

120

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