книги из ГПНТБ / Лебедев А.А. Динамика полета беспилотных летательных аппаратов учеб. пособие
.pdfгде хв — расстояние вдоль оси х от оси вращения передней кон соли до точки схода вихря;
Хц — расстояние вдоль оси х от точки схода вихря до сере дины САХ задней консоли;
уп — смещение по вертикали поверхности 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 а. |
дБ, |
д Б , |
д Б , |
При малых углах а и бі это выражение принимает вид
5і |
~ с |
-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
Подставив сюда выражения 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