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

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

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

где хв — расстояние вдоль оси х от оси вращения передней кон­ соли до точки схода вихря;

Хц — расстояние вдоль оси х от точки схода вихря до сере­ дины САХ задней консоли;

уп — смещение по вертикали поверхности I относительно поверхности II (уи>0, если поверхность I расположе­ на выше поверхности II).

Следует заметить, что знак у в не влияет на величину і, т. е.

Ң - у в) = Ң у в)-

Так как производная е“Р определяется

при а = 0 и 6і=0,

то

У, = Уи-

(3

-37)

В частном случае, когда поверхности I и II расположены в одной плоскости, ув = 0.

Сделаем некоторые замечания по практическому использованию рис. 3.17.

1) Значение

і обычно находят

путем

интерполяции по параметрам

Dn и

1

2zB

 

если

вихрь проходит

вблизи

корпуса,

та-

----- при

—;—

= const. Однако,

Чк II

Л

 

 

 

 

 

 

 

 

2гв

 

кая интерполяция не всегда возможна. Например, при £>п =

0,45

 

и

III

 

 

 

 

 

 

 

 

 

 

 

 

 

:0,55 можно использовать

только график

£>[[ =

0,4,

а график Оц

=

0,6

использовать нельзя, так как

 

2zB

= 0,55

попадет

внутрь корпуса.

точка —;—

 

 

 

 

Iи

 

 

 

2гв/Лі - D И

 

В этом случае рекомендуется вести интерполяцию при

 

1

- D ,

 

 

 

 

 

 

 

 

 

 

 

 

 

= const

I в данном примере

2 г вУи

 

' 0,55 — 0,45

0,182.

На

 

1 — А

 

 

1 — 0,45

=

 

 

 

 

 

 

 

 

 

 

графике Dn =

0,4 находим значение і при

2 z Blln

= 0,4 +

0,182 ( 1 — 0,4) =

= 0,509,

а на

графике Dn = 0 ,6 —при

2zB//n = 0 ,6+0,1-82(1— 0,6) = 0,673.

 

171

2) Если диаметр

корпуса меняется

по

его длине, то координата

вихря

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

 

*ВІІ =

- у ^ Л[0 І + 7 в(/І -

°

І) ]2+ ( ° n - Dl) -

(3'38>

вытекающей из уравнения постоянства расхода.

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

Рис. 3.19. Схема для определения величины ф 6

I

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

f g « >

/ 1 — ,

(3.39)

у

— 1

 

задняя несущая поверхность полностью выходит из зоны влия­ ния вихрей, вследствие чего еСр= 0 (это справедливо при 8і= = */п= 0).

Для учета ограниченности зоны влияния вихрей можно ре­ комендовать следующий приближенный способ: введем в фор­ мулу (3.34) множитель ф£, представляющий собой отношение части площади задней консоли, находящейся внутри конуса Ма­ ха, ко всей площади консоли. Значения фЕ в общем случае за­ висят от М, а, 6і и геометрических параметров летательного аппарата. При расчете производной Вер следует принимать а = ='бі = 0 и находить фЕ путем графического построения. В дозву­ ковом и околозвуковом потоках ф Е= 1.

172

1.5. КОЭФФИЦИЕНТ ТОРМОЖЕНИЯ ПОТОКА

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

ЯI Ян

СтІ : СТІІ:

где <7і, qn — средние скоростные напоры потока в области пе­ редних и задних несущих поверхностей.

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

 

 

 

 

= і ~ м)

 

V,

М,

(3.40)

 

 

 

 

;

*ТІІ = І V

м

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

Мі = Мф &ті;

М ц = М у

kin.

(3.41)

Торможение потока вызывается в основном обтеканием носо­

вой

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

и обтеканием передних

консолей.

Рассмот­

рим,

.например, явления,

 

 

 

 

возникающие при обтека­

 

 

 

 

нии

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

 

 

 

 

ком тела с конической но­

 

 

 

 

совой частью (рис. 3.20).

 

 

 

 

При

прохождении

потока

 

 

 

 

через косой скачок уплот­

 

 

 

 

нения происходит

потеря

 

 

 

 

полного

давления, харак­

 

 

 

 

теризуемая

коэффициен­

 

 

 

 

том

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

 

 

 

 

ления а. Как известно из

 

 

 

 

теории косых скачков, ве­

Рис. 3.20. К определению коэффициента

личина а

зависит

только

торможения потока kT', вызванного об­

от параметра

Msin ß, где

теканием конической носовой части кор­

ß — угол

наклона

скачка.

 

 

пуса

 

В свою

очередь,

угол ß

 

 

 

 

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

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

173

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

Ро =

± № ) к~'

(3.42)

Р

 

 

Аналогичная зависимость имеет место и в произвольной точке / в окрестности цилиндрической части тела:

 

 

 

к

 

Z0L = (

1+ _ ^ ± м Г ) * _І .

(3.43)

Р\

\

2

/

Л

Если расстояние от начала цилиндрической части тела до точки 1 достаточно велико (порядка двух диаметров тела), то, как показывают расчеты, р \ ~ р . Тогда, разделив (3.43) на (3.42), получим

откуда

(3.44)

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

Зависимость &/(М, ЯНОс), рассчитанная по приведенным выше

 

 

 

 

 

 

формулам,

изображе­

кт

 

 

 

 

 

на на

рис. 3.21.

если

 

 

 

 

 

 

В том

случае,

 

 

 

 

 

 

форма

носовой

части

 

 

 

 

 

 

корпуса

отличается

от

 

 

 

 

 

 

конической,

можно

ре­

0,8

 

 

 

 

 

комендовать

следую­

 

 

 

 

 

 

щий

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

0,7

 

 

 

 

 

способ определения k?':

 

 

 

 

 

1)

подсчитать

схвос

 

 

 

 

 

 

 

 

 

 

 

 

для заданной носовой

0, 6 :

 

 

1,8

 

2,6 Д дас

части,

как

указано в

0,6

1

0

2,2

гл. IV;

 

 

 

 

Рис. 3.21. График для

расчета

коэффициен­

2)

пользуясь графи-

 

та торможения потока к?'

- К О М Сх н о с

(М, К о с ) для

174

конических носовых частей (см. рис. 4.11), найти удлинение экви­ валентного конуса, имеющего тот же сх1Юс при данном числе М;

3) по рис. 3.21 найти &/.

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

Рис. 3.22. График для расчета коэффициента тормо­ жения потока kT*

лом Маха. Влияние последних двух факторов отражено на эмпи­ рическом графике рис. 3.22, относящемся к летательным аппа­ ратам обычной схемы, у которых S i^ S n . Чтобы использовать эту зависимость для аппаратов любой схемы с любым соотно­ шением Si и Sn, можно рекомендовать приближенную формулу

+ S n jS l

(3.45)

1 + Sjj/Sj

 

где Л?* определяется по рис. 3.22.

что при обычной

Из выражения (3.45) следует, в частности,

схеме летательного аппарата (Si^>Sn) k r

с , а для аппара­

тов схемы «утка» (Si<cSn) k!r ~ 1.

Таким образом, расчет коэффициента торможенш потока должен производиться по формулам:

— в области передней несущей поверхности

 

£ТІ = £т;

(3.46)

175

— в области задней несущей поверхности

k т11--kikf'

(3.47)

§2. ПРОИЗВОДНЫЕ КОЭФФИЦИЕНТА ПОДЪЕМНОЙ СИЛЫ ЛЕТАТЕЛЬНОГО АППАРАТА ПО УГЛАМ ОТКЛОНЕНИЯ

ОРГАНОВ УПРАВЛЕНИЯ с*' и с*11

Продифференцируем выражение (3.1) по углу 6і:

дсу

дсу1 cos а —

дсхі sin а.

дБ,

д Б ,

д Б ,

При малых углах а и бі это выражение принимает вид

~ с

-8І

(3.48)

сУ

У1 '

Поделив равенство (3.6) почленно на qS и взяв производную по 6і, получим

I __

т / І Г

(3.49)

 

 

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

4 і = Кі„з.кРtfwAjr

(3.50)

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

ац = а — sscp8b

откуда следует

4 п = “ (^1из.кр^-)іІ£ср-

(3-51)

Аналогичным способом можно получить выражение частной производной коэффициента подъемной силы летательного аппа­ рата по углу бц:

с " =

= (cl i u3.KpKionSkT)u.

(3.52)

Таким образом, для расчета производных коэффициента подъемной (или нормальной) силы летательного аппарата по

углам бі и 6ц необходимо определить величины Као, п и еср.

176

2.1.КОЭФФИЦИЕНТЫ ИНТЕРФЕРЕНЦИИ £s0 и Кьо

В§ 1 была рассмотрена картина интерференции корпуса и несущей поверхности для случая «сш», когда их геометрические углы атаки были равны между собой. Рассмотрим второй слу­ чай, когда угол атаки корпуса равен нулю, а консоли отклоне­ ны относительно его оси на угол б (случай «60»).

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

ле изолированных крыльев: Уік~ Киз.крКроме нее, по тем же причинам, что и в случае «<ш», возникает индуцированная нор­ мальная сила корпуса Унф-

По аналогии с коэффициентами k aa и Ки введем коэффи­

циенты интерференции &ао и Къо!

 

YІк

Пс + Уи

Кьо--

 

к5

1 из.кр

л 1 из.кр

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

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

Каа. П рИ ВѲДеМ ИХ В ОКОНЧЭТеЛЬНОМ ВИДв!

 

&ао =

^£0*ц.сХм;

(3.53)

/Сзо =

[^50+ (ÄTSO-- ^50) (IJ1 ѵ с*м..

(3.54)

Коэффициенты k*o

и Кго

связаны с коэффициентами kla

и Кааследующими соотношениями, вытекающими из теории тон­ кого тела:

(С У .

(3.55)

и*

'

 

^аа

 

 

КІ0 = С

 

(3.56)

(см. рис. 3.11).

Функция F(LXB), учитывающая влияние длины хвостовой масти корпуса, при М >1 по-прежнему определяется выражени­ ем (3.30). При M s^l F( LXB) = 1. Поправка на влияние числа Маха им берется по рис. 3.13.

Поправочный множитель кп.с, учитывающий наличие погра­ ничного слоя корпуса, несколько отличается от аналогичного множителя Хп.с, полученного ранее для случая «аа» [см. вывод формулы (3.16)]. В случае «60»

— __ fe80S K ^ S K

^80^к

177

Рис. 3.23. Рули типа поворот­ ного оперения

Подставив сюда выражения 5 Ки S / , найдем:

[ 1- D (1+ в*)]

— D (1 + В * )

 

 

2D

(3.57)

 

 

(1 - 0 ) 1

D

 

 

+ 1 -2D

 

Относительная толщина вытеснения б* подсчитывается по фор­ муле (3.17).

2.2. ОТНОСИТЕЛЬНАЯ ЭФФЕКТИВНОСТЬ ОРГАНОВ УПРАВЛЕНИЯ

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

с8

(3.58)

п = - £ ~ .

Г9

 

СУ1

 

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

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

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

иметь угол стреловидности %р¥=0. Тогда угол отклонения рулей отсчитывается ,в плоскости, перпендикулярной оси вращения.

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

178

1. Изменение угла установки стабилизатора вызывает точно такое же изменение его угла атаки, в то время как поворот руля на некоторый угол б в общем случае приводит к изменению угла атаки на меньшую величину, равную б cos %р- В пределе, когда 7р= 90°, при повороте руля угол его атаки совершенно не ме­ няется *.

Отсюда следует, что поворот рулей вокруг их оси вращения в l/cos5Cp раз менее эффективен, чем изменение угла установки стабилизатора.

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

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

 

 

 

 

 

cs

(3.59)

 

 

 

 

n = —^— = kщсоэхр.

 

 

 

 

 

eh

 

Ориентировочные значения &щ:

 

при

М <М Кр

£щ = 0,8-у0,85;

 

при М>1,4

= 0,95-М.

 

Концевые рули. В настоящее время некоторое распростране­

ние

получили

концевые

рули

 

(рис. 3.24), крепящиеся к не­

 

подвижному центроплану, ко­

 

торый называется

иногда так­

 

же пилоном. Согласно теории

 

тонкого

тела

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

 

эффективность

таких рулей не.

 

зависит от диаметра корпуса, а

 

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

только отноше­

 

нием размаха рулей к размаху

 

консолей

оперения

/р//к. Хотя

 

этот

вывод приближенный, он

 

дает основание воспользовать­

 

ся

для

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

более

_

точными

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

дан-

ными, полученными для

изолированных крыльев

(.0 = 0). Отно­

сительная эффективность концевых рулей при 0 = 0 определяет­ ся выражением

n = n1kmcos Хр.

(3-60)

Входящий сюда множитель cos%p, как и ранее, учитывает тот факт, что при отклонении рулей на угол б их угол атаки меняет­ ся на величину б cos хР-

* Строго говоря, это справедливо только при а = 0.

179

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

Рули, расположенные вдоль задней кромки стабилизаторов.

На летательных аппаратах, предназначенных для полета с до-

относительную эффективность конце-

торов

вых рулей

 

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

 

 

Рис.

3.27.

Характер

зависимости-:

 

 

n = f( М) для рулей,

расположен­

 

 

ных

вдоль

задней кромки стаби­

о

1

г м

 

лизаторов

 

 

 

 

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

180