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

Приближенное аналитическое решение задачи о возможности попадания в режим плоской авторотации вращающегося асимметричного тела, входящего в атмосферу Текст научной статьи по специальности «Физика»

CC BY
120
43
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Ученые записки ЦАГИ
ВАК
Область наук

Аннотация научной статьи по физике, автор научной работы — Ларичева В. В.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Приближенное аналитическое решение задачи о возможности попадания в режим плоской авторотации вращающегося асимметричного тела, входящего в атмосферу»

_________УЧЕНЫЕ ЗАПИСКИ ЦАГИ

Том XXXI : 2 000

№1—2

УДК 629.7.015.076.8

ПРИБЛИЖЕННОЕ АНАЛИТИЧЕСКОЕ РЕШЕНИЕ ЗАДАЧИ О ВОЗМОЖНОСТИ ПОПАДАНИЯ В РЕЖИМ ПЛОСКОЙ АВТОРОТАЦИИ ВРАЩАЮЩЕГОСЯ АСИММЕТРИЧНОГО ТЕЛА, ВХОДЯЩЕГО В АТМОСФЕРУ

В. В. Ларичева

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

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

Для входящего в атмосферу, вращающегося около центра масс тела в [1] — [5] выявлено влияние малой асимметрии тела на возможность возникновения явления авторотации — непрекращающегося вращения с нарастающей угловой скоростью.

В данной статье продолжается начатое в [2] — [5] аналитическое исследование случая, когда слабо асимметричное тело входит в атмосферу, вращаясь в продольной плоскости. Аналитически решена краевая задача первого порядка об определении критического решения, являющегося граничным между решениями вращательными: 1) авторотационными и 2) ведущими к захвату в область колебаний около центра масс тела.

В отличие от [2] — [5] в данной статье анализ основывается не на использовании метода усреднения, а на свойствах квазистационарного приближения, достигающего близости к захвату одновременно с точным ре-

шением. В качестве примера использования квазистационарных решений в других задачах сошлемся на [6], [7].

1. Уравнения движения. Рассмотрим входящее в атмосферу (на большой скорости У0 под малым углом 0О < 0 наклона к горизонту) слабо асимметричное тело, вращающееся около своего центра масс в продольной плоскости с постоянной угловой скоростью (О .

Как известно [2] —• [5], в пренебрежении малым в верхних слоях атмосферы аэродинамическим демпфированием, уравнения движения около центра масс спускаемого тела в продольной плоскости имеют вид

^- = у{т2+г1Ат2)д, ^ = (у = ^//4); С1)

здесь а, со _ — угол атаки и угловая скорость тангажа, / — площадь миделя и длина тела, /_ — момент инерции относительно поперечной оси тела, ^=рГ2/2 — скоростной напор.

Тело является почти осесимметричным и имеет небольшое смещение Ду центра масс в продольной плоскости ниже оси геометрической симметрии и малое искажение поверхности. Из-за нарушения симметрии тела возникает дополнительный аэродинамический момент с коэффициентом Дт2(а), малость которого в (1) отмечена малым параметром в1 > 0.

Коэффициенты аэродинамических моментов: основной т2 (а) от строго осесимметричной формы тела и дополнительный Ат. (а) от асимметрии тела — описываются периодическими по а с периодом 2л функциями и могут быть представлены рядами Фурье. Ряд для от. (а) состоит из одних нечетных гармоник 8т(иа), п = 1,...; ряд для Ат!,(а) содержит нулевой член, поэтому

2% 2гс

' т = — Гтг(а)с1а = 0, Ат - — | Ат2(а)с1а ^0. (Г)

• 2л •' 2к •'

о о

Без последнего неравенства (1') авторотация невозможна, см. ниже п. 4.

Как в [2] —[5], предполагаем, во-первых, что на интересующем нас верхнем участке траектории скорость У0 и малый угол наклона |90|«1 постоянны. Во-вторых, плотность атмосферы р(Н) предполагаем изменяющейся по высоте Н по экспоненциальному закону с постоянным логарифмическим градиентом X. В итоге, скоростной напор д(т) медленно изменяется при спуске тела по экспоненциальному закону

9(т) = р(Я)К02/2 = ^0ехр(Я.К0|8те0|-?). (2)

Здесь фактическим малым параметром является безразмерная величина ]зт©01, однако медленность изменения скоростного напора удобно

характеризовать размерной малой величиной еХ = Л,| sin0ОjК0, с-1, вводя безразмерное медленное время как х = sXt.

2. Невозмущенное движение. Сепаратриса. Полагая в(1)б = б, =0, придем к автономной, консервативной невозмущенной системе, для которой первым интегралом служит «энергетическое» соотношение

со2 =Ф2 +2vqF(a), (3)

где

ос .. 2%

F{ol) = Jmz(a)da, F0 = — JV(a)c/a = 0, Fm = max| F(a)|. (3')

2n о

Здесь константа со2 имеет смысл приведенной удвоенной полной энергии невозмущенной системы, по ней можно судить о характере движения около центра масс тела: вращение — при со2 > 2vqFm, колебания — при со2 <2vqFm.

Апериодическому режиму соответствует граничное равенство

со2 =2 vqFm, (4)

которым определяется особое решение уравнений невозмущенного движения — сепаратриса. На фазовой плоскости а,со. она представляет собой самопересекающуюся кривую, проходящую через седловую особую точку. Уравнение сепаратрисы получим из (3) с учетом (4):

<s>l=2vqFm(\ + F(a)/Fm). (4')

3. Возмущенное движение. «Возмущенная» сепаратриса. Следуя теории возмущений, воспользуемся соотношением (3) при переменных q,со2 для замены в (1) переменной coz на со2 и перейдем к эквивалентной (1) системе уравнений возмущенного движения

со2 =G)2(x) + 2vq(x)F(a), ^ = оог, (5)

at

2

= 2v^(T)[s]Amz(a)coz-sA.F(a)l, (eX = ^F0|sin ©о|) (5')

dt \ \ \>

с начальными условиями a0, coq =®z0 ~2v^0F(a0).

Соотношениями вида (4), (4') при q{x) = var опишем «возмущенную» сепаратрису, которая, в отличие от невозмущенной, не удовлетворяет уравнениям возмущенного движения (5), (5'), решения которых могут пересекать такую «сепаратрису».

4. Авторотация и задача о захвате. Если входящее в атмосферу асимметричное тело вращается около центра масс (в продольной плоскости) с достаточно большой по модулю угловой скоростью юг , то далее

вращение спускаемого тела не только продолжится, но из-за асимметрии тела его угловая скорость может нарастать (плоская авторотация [2] — [5]). Разумеется, авторотация невозможна при отсутствии нарушений осевой симметрии тела.

Для безопасного спуска тела в атмосфере необходимо уменьшить ве-

личину

со

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

переход

вращения в колебания около центра масс тела. Поэтому возникает вопрос об определении граничного значения со^ называемого критическим, при котором еще происходит плоская авторотация, но уже при небольшом уменьшении величины ш2коР произойдет захват в область колебаний.

Укажем «энергетический» критерий захвата. Для этого перепишем выражения медленно изменяющейся функции (4) при q = уаг и ее производной как

= 2vq(x)Fm

d(a2s jdt = 2zkvq(x) Fm

/ a \

Fm = шах Jmz(a)da

V J

(6)

Функциями (6) характеризуется энергетический уровень «возмущенной» сепаратрисы и скорость изменения уровня.

Если в некий момент tx характеристика со2 (?) энергетического уровня решения (5), (5') и ее производная dot2/dtне превысят соответствующих сепаратрисных значений (6), то произойдет захват решения (5), (5') в область колебаний; причем в случае точного совпадения с сепаратрисными значениями (6) решение (5), (5') будет критическим — граничным по отношению к захвату.

Отметим, что критерий [2] — [4], основанный на свойствах нижней амплитудной кривой угловой скорости, эквивалентен данному.

5. Амплитудные кривые и линия центров колебаний угловой скорости. В предположении, что на каждом периоде функция mz(a) имеет лишь по одному максимуму и одному минимуму, рассмотренная выше невозмущенная система имеет устойчивое и неустойчивое положение равновесия, в которых балансировочные значения а7 ,а/7 угла атаки являются главными значениями корней уравнения mz{<х) = 0. Поэтому при а = а, или а п функция (3') экстремальна

кк

F(a/,//)= Jwz(a)<*i

-■ + F ---------- J №

1,11

и тогда из (3) следуют выражения для амплитудных кривых (характеристик) угловой скорости

В равенстве (7) слева нижний индекс «в» (или «н») и справа знак «+» (или «-») соответствуют верхней (или нижней) амплитудной кривой. Неравенствами (7), заключенными в скобки, отмечаются пределы изменения точной функции (/).

Пользуясь условием ^(а) = 0, которому удовлетворяет корень с главным значением а»(а/+ад)/2, выделим из (3) срединную характеристику со.?, определяющую линию центров колебаний угловой скорости

На этой линии потенциальная энергия невозмущенной системы обращается в нуль, и полная энергия равняется кинетической, величина которой пропорциональна среднему арифметическому из значений (7). Последние, как видно из (7), (7'), можно определить, если известна величи-

В возмущенном движении характеристики (7), (7') становятся медленными переменными. Из (7), (7') при q = уаг усматривается, что в момент

щенной» сепаратрисы.

В связи с отмеченным свойством синхронности, очевидно, что для выявления условий близости решения (5), (5') к захвату достаточно проследить, например, поведение переменной со2(х).

6. Приближенное решение задачи о захвате. Изложим (пока формально) способ получения бесколебательного приближения к решению (5), (5'), достигающему близости к захвату.

Пользуясь заменой

(7')

на (7').

достижения точной функцией (?) «возмущенной» сепаратрисы ее достигнут также: функция со2 (т) в точке сепаратрисы, где ^(а) = 0, функция

2 2

0)г (т) — в седловой точке, функция со, (т) в вершинной точке «возму-

ю2 =ю2 +2vq(т)E1f(a), £(а) = Р(а) + є,/(а),

\

(8)

а

V

У

перейдем от (5), (5') к эквивалентной системе

О Г^_г О ^

сог = 5 + 2vqF(a), сіа/ск = ю2,

(9)

й(Ь2 |dt = 2zXvq\k■^(дz-F{a'^, (А:] = Є] Ат2о /єА.). (9')

Ряд Фурье периодической функции /г(а) не содержит нулевого члена, поэтому систему (9), (9') можно рассматривать как колеблющуюся около многообразия, определяемого укороченными, не зависящими от а уравнениями:

а = ю, а>г=в>2, с/со2 = 2ъкщк\(Ь. (10)

Последнее уравнение (10) имеет бесколебательную правую часть и с учетом (2) легко интегрируется в виде

<в = <»о+*1'/(9-4,о)> {к\=гхАт2й/ъХ). (10')

Укороченная система (10) выделена из (9), (9') посредством условия

Рт=о, (п)

соответствующего характеристике — линии центров колебаний угловой скорости со2(О- (Аналогичным (11) является условие, приведшее к характеристике (7')).

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Для (10) допустимо любое начальное значение о , но другое условие — частное: а0 = а0, где а0 «(а7 + ап )/2 — главное значение корня уравнения (11), а балансировочные значения а/,а7/ являются главными значениями корней уравнения

тг (а) + е , [ Д/иг (а) - А/яго ] = 0, (11')

в котором учитывается вклад гармоник малого коэффициента е,Ат2(а) аэродинамического момента, возникающего от влияния малой асимметрии тела. Тогда как выше при выводе (7') использовались в качестве а/,ая корни боЛее простого уравнения: т2 (а) = 0, следующего из (11') при £( - 0.

Решение (9), (9'), приближающееся к захвату, будет критическим, если в некий момент компоненты решения м^,(<зй)2/<#) сравняются по величине с соответствующими сепаратрисными, определяемыми соотношениями (6). Поэтому потребуем выполнения при / = равенств

со2-2\д^Рт, 2ekvqlklG)l=2eXvqlFm, из которых определим значения

®1=^Дь 91= ^т/2^, (б^ =1п(?1/?0)), (12)

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

Подставляя (12) в (10') и разрешая последнее относительно <в0,— решая тем самым краевую задачу первого порядка, получим формулу для определения критического начального значения угловой скорости тангажа

“гї =®0=А:1^0+^т/2^ *Рт/2ки (*, = Є,Длїг()/єА.). (13)

На основании (12), (13) устанавливаем, что

®і/®о =(^іЯ/*і)/[(^я./2Лі) + *іУ^0]*2. (13')

Формулы (13), (13') выписаны в учетом того, что величина кх\н]0 пренебрежимо мала уже при Н0 > 110 км (при данных приведенных ниже примеров).

Начальным условием (13) с добавлением к нему значения а0 = а0 « «(а/ +ап)/2, являющегося главным значением корня уравнения (11), определяется критическое решение (9), (9'). Вдоль последнего отношение <7і/<7о может меняться (в зависимости от величины #0) в десятки, сотни и более раз, но величина отношения согласно (13'), остается почти

неизменной.

Как показали расчеты, точная зависимость от а0 такова, что при #0 > 110 км на довольно большом массиве значений а0 величина со^ почти постоянна и к ней близок (с некоторым завышением) результат формулы (13). Практически при тех же условиях близок к точному (с некоторым занижением) результат формулы

*0,41/^/*,, = тах|іиДа)|, *, -.(^Діи^/єА,), (14)

полученной в [5, стр. 163] на базе предложенной автором специальной модификации метода усреднения. Ниже, см. (17) и рис. 1, 3, 4, ограничимся случаем = Ах.

Укажем на связь (13) с результатом [2] — [4] для со^ и г,. При ?0 =• 0 верно ю о «ю0, поскольку, как видно из (8), разность квадратов этих величин при этом пропорциональна <?о * 0 и> значит, пренебрежимо мала. Но тогда, как следствие, (13) совпадает с соответствующей расчетной формулой [2] — [4]. Кроме того, здесь и в [2] — [4] время , потребное для достижения близости к захвату, одинаково определяется формулой вида (12).

Но есть и различие. Вклад от суммы малых гармоник Є[[Дт2(а) -

-Дтго | отсутствует (аннулируется при усреднении) в результатах [2] — [5].

Тогда как выше в балансировочном уравнении (1Г) упомянутая сумма гармоник учитывается как малая добавка, что в итоге приводит к уточнению начального значения угла атаки а0 = а0 для системы (10). Такого рода уточнение наболее значимо в случаях, когда Я0 < 110 км.

Отметим наиболее важное. Уравнениями (10), выделенными из (9), (9') посредством условия (11), описывается образ, относительно которого происходят колебания в системе (9), (9')- Для последней из (10) получены приближенные выражения медленно изменяющихся характеристик ю2,сйо2/<й. Причем, из-за бесколебательности с/<а2/<Й отпала необходимость привлечения метода усреднения для исследования (10).

7. Особенность и геометрический смысл приближенного решения. Убедимся, что приближение (10), (11) можно рассматривать, как начинающееся из бесконечно удаленной особой точки, где # = 0, * = - оо, и совпадающее в этой точке с точным решением (9), (9')- Действительно, при г = - оо (когда, как видно из (2), ^(х) -» 0) из (9), (9') и (10), (11) вытекают одни и те же предельные соотношения

ш2=ю2=со20, йаДй = ю=сого, (15)

которыми описывается стационарный режим кругового вращения около центра масс тела с постоянной угловой скоростью <а2() («кувыркания»).

Сравнивая (10), (11) с (15), заключаем, что приближением (10), (11) описывается квазистационарный режим почти кругового вращения тела с почти постоянной угловой скоростью ю(х). Последняя, как видно из (13'), изменяется в весьма узких пределах на всем интервале 0 < £ < ^ до захвата.

Следовательно, приближение (10), (11) начинается при q& 0 стационарным режимом «кувырканий», совпадает при этом с результатом точного решения (9), (9'), далее (вплоть до захвата) развивается квазистационарно и пересекается с точным решением периодически — в моменты выполнения на последнем условии (11) Р’(а) = 0.

Простые выражения (10), (10') для со2, й?ю2/<Л приближенно характеризуют энергетический уровень и скорость его изменения вдоль решения (9), (9'), поэтому (10), (10') использованы выше при аналитическом решении краевой задачи о близости решения (9), (9') к захвату.

8. Число оборотов ./V, до захвата. Величина ЛГ1 определяется формулой

Л^1 = (а, -а0)/2п,

где а], а0 — концевое и начальное значения угла атаки.

После подстановки (10') в (10) получим приближенное уравнение, описывающее нарастание по времени функции | а(/) |. Поскольку на самом деле вблизи захвата рост | а(0 | замедляется, то неудивительно, что значение а] (и ТУ,) оказывается завышенным. ;

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

а1 -а0 = 27сЛг1 = с50?] ■»

(16)

где /), «о определяются формулами (12), (13).

Соотношение (16) получается в результате интегрирования по времени на интервале [0, ] выражения (10') для <5 при аннулировании к моменту захвата интеграла от асимметричного члена кх V (<? - ), что согласуется

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

Примечание к пунктам 6—8. Отметим, что методическим результатом [2] — [4] являлись также аналогичные (10), (10') укороченные уравнения. Однако в [2] — [4] уравнения вида (9), (9') сначала усреднялись по време-

{ 2к \ ни (с переменным периодом Т = |с?а/со,), затем усреднные уравнения

о

упрощались вдоль линии центров колебаний угловой скорости — в итоге получались уравнения вида (10), (10').

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

Обратим внимание, что, если для (12) предполагать порядок величины д()1 не превышающим 8, то при е -» 0 из (12) следует стремление гх к бесконечности по закону 11п( 1/е) |/в.

9. Сравнение с точными численными решениями. Ограничимся случаем, когда форма спускаемого тела осесимметрична, а центр масс тела смещен ниже оси его геометрической симметрии на малую величину Ду<0: |Ду//|=е,|Ду|«1.

Примем, что аэродинамические коэффициенты данного тела определяются простейшими формулами:

тг=-Ах§1па, сх - сТо + Вх сова. (17)

Тогда

Атг = сг(а)е1Ду = ъхАтпго + гх&уВх сова, (в]А/?г^0 =е!ДуеТо), (18)

а

Р(а) = jmz(a)da, Р(а) = Р(а) + ЕхАуВх$та. (19)

Для (10) в качестве а0 подходящими являются корни а0 =90° и 55° соответственно уравнений F(a)~0 и /’(а) = 0, см. (10), (11'). Далее для простоты ограничимся значением а0 « тс/2, т. е. примем -Р(а0) «Р(а0). Хотя в случае Я0 = 100 км это ухудшает точность приближенного результата.

Для коэффициентов, входящих в (1) или (9), (9'); (17) — (19), зададим численные значения

4=0,05, сТо = 0,2, 5] =1,2, е1Ау = Ау/1 = -0,03, е{АусХ0 =-0,006, у = 57//г =0,0543 м2/кгс-с2.

Пользуясь (17), (18), определим для системы (1) или (9), (9') устойчивое (X/ «-41° и неустойчивое а и =150° балансировочные значения угла атаки.

Для начальных высот Я0 = 100-=-300 км зададим начальные: скорость входа У0 и скоростной напор д0(Н0) в соответствии с формулами

Г0=Л/290(100)/р(100), 90(Я0) = р(Я0)90(100)/р(100),

где р(Я0) — плотность атмосферы на высоте Я0 (по таблицам стандартной атмосферы 1981 г.); в частности, р(100) = 5,65 ■ 10-8 кгс с2/м4, #о(100) = 3 кгс/м2 — плотность атмосферы и начальный скоростной напор при Я0 = 100 км.

Зададимся значениями параметров

еХ, = А.|8т0„|Ко= (0,05 + 0,07), с-1 кх =е1ДусТо/еЛ, = -(0,12 + 0,0858) = -(2,4 + 1,72)4, с-

При перечисленных предпосылках для (1) численно решалась задача захвата — краевая первого порядка — как задача Коши интегрирования уравнений (1), или (9), (9'). Для каждого а0 подбиралось такое значение сого, чтобы в некий, заранее неизвестный момент ^ приближенно выполнялось условие со, = 0; далее проверялось, чтобы при небольшом уменьшении | со2() | происходил захват решения в область колебаний.

Расчеты показали, что для уравнений (1) или (9), (9') решение типа, проанализированного выше (см. пункты 6—8), реализуется в довольно широком диапазоне изменения параметров а0,еХ,Н0.

Так, при Я0 = ЮО км на точном критическом решении на большом массиве значений а0 е[-я, я] достигается вблизи захвата концевое значение угла атаки а1 «-10 рад. При а0 =я/2, еХ = (0,05 + 0,07) с-1 точное решение имеет параметры = <й0 =-0,24 с-1, /,=(47,5 + 38) с,

^ - 1,85; а из приближенных формул (12), (13), (16) при этом следуют данные: а>0 =-0,255 с-1, = (47,2 + 43,2) с, =(1,73 + 2,05), приемле-

мо согласующиеся с точными.

Во многом аналогичные выводы можно сделать по результатам расчетов для Я0 = 120 + 250 км; а при условии гк> 0,0635 с-1 также для Я0 = 300 км.

Н0 =10 0 хм, £к=0,051 с~\ д^Зкгс/м2

КО -1 а ПО формуле (и)

0,и • 1 * - •» » (ш)

-з -г -1 о -0,1 1 2 а.0,рад

~о,г ' и ' ' ' * '

Некоторые результаты точных численных решений даны на рис. 1—4, где темными и светлыми кружочками указаны аналитические значения со%, рассчитанные по

формулам (13), (14). Рис. 1, 2 соответствуют случаю #0 = 100 км, вХ - 0,057 с-1, д0 = 3 кгс/м2 ; рис. 3 — случаю Н0 = 200 км, еХ = 0,0635 с-1, q0 = 1,35-10-3 кгс/м2 (р0 = 2,55-10_и кгс-с2/м4);

рис. 4 — случаю #о=300 км, еЯ, = 0,07 с-1, q0 = 1,04-10-4 кгс/м2 (р0 = = 2-10-12 кгс-с2/м4].

Рис. 1

Ш*Р,С 1 Нкм, е А=0,057с \д0 =3кгс/м г, <х.0=0,95рад

кр -7 Но=200км 0,2 сХ=0,0635ч кр “V* ■ о по формуле(/?) Н0—300км 0,2 • » „ е\—0,07 1 - о по формуле (и) • » >■ {№) 1 1 1

-3 -2 -1 В -0,2 1 2 о.0,ра.д -3 -2 -1 0 1 . -0,2 1 2 <х,0град ~ ’ •

Рис. 3 Рис. 4

На каждом из рис. 1, 3, 4 изображена зависимость со*р от а0, определяющая граничную кривую, ниже которой располагается область авторо-тационных начальных условий, выше — область начальных условий, ведущих к захвату. Влияние д0 при аналитическом определении ю*р проявляется в случае Н0 = 100 км: значение при а0 = л/2 вычисляется по формуле (13), затем пересчетом по «энергетическому» соотношению, входящему в (9), находятся значения при q = qQ, а0 Ф л/2, см. рис. 1. В случаях Н0 = 200 + 300 км величина ц0 пренебрежимо мала.

В правой части каждой из граничных кривых на рис. 1,3,4 возникает небольшой скачкообразный излом, соответствующий наличию небольшой области неединственности решения краевой задачи о захвате [3], [4]. Причем решения, которыми определяется правая (после излома) часть граничных кривых рис. 1, 3, 4, имеют уменьшенное на 2л значение |а,| по срав-

нению с решениями, определяющими левую часть. Здесь (слева) а] « «-10 рад (рис. 1), а, «-41,3 рад (рис. 3), оц »-54 рад (рис. 4), что согласуется с данными аналитических формул (12), (13), (16).

На рис. 2 дан пример точной численной зависимости се>*р от времени t

(для Н0 = 100 км, вЯ, = 0,057 с-1, q0 = 3 кгс/м2 ) при coZ() = -0,26 с-1, а0 = = 0,95 рад. Последнее является средним между балансировочными значениями ai —-41°, аи =150° системы (1).

Для сравнения на рис. 2 нанесена зависимость от времени бесколеба-тельной угловой скорости <в(т) приближения (10') на интервале [0, t\ =46 с]. При t = ty в точном решении угловая скорость близка к нулю. На фазовой плоскости a,coz ему соответствует точка, близкая к седловой особой точке на сепаратрисе, где проходит также частное решение типа зависания в седловой точке. Кроме того, как показано в [3], [4], при изменении а0 может приблизиться к нулю предыдущий min|coz|. Уравнивание предыдущего и последующего значений min|coz(0| происходит в области неединственности решения задачи: при дальнейшем увеличении а0 предыдущий min| (oz(01 становится меньше последующего, что соответствует решению с уменьшенным | (Х|| на 2л.

10. Эмпирические закономерности. Отметим две из них для решений (9), (9').

1) Расчеты показали, что на критическом решении в окрестности вершин нескольких предзахватных нижних витков колебательной функции юz(t) начинает выполняться неравенство со2 < 2vqFm, тогда как признаком области вращения около центра масс тела является неравенство обратного знака со2 >2vqFm. Иными словами, в возмущенном движении на критическом решении наблюдаются некоторые отдельные непродолжительные нарушения условия периодичности ю2 по а.

2) Из расчетов критических решений усматривается, что медленно возрастающая функция со2 становится почти постоянной в течение времени прохождения предзахватного витка колебаний соz(t). (Это объясняется уменьшением в предзахватной ситуации влияния асимметрии тела на параметры решения.)

Значит, можно указать следующие условия близости решения (9), (9') к захвату в некий момент /,

(b\ =2vqxFm', (*йэ2/<#) =0 при F(a) = Fm. (20)

Согласно (20), при t = fj в окрестности верхней части предзахватного витка приближенно выполняется сепаратрисное энергетическое соотношение и становится постоянной величина со2. А далее при t>tx, как под-

тверждают расчеты, вместо (20) начинает выполняться неравенство со -const < 2vqFm, смысл которого пояснен выше.

Из (20) с учетом (9') найдем соответствующие моменту tx концевые, аналогичные (12) значения

®1 = ^fm/2ku v<ll=Fm/4k\’ (21)

После подстановки (21) в (10') и разрешения относительно <в0 получим аналогичные (13) — (14) расчетные формулы

=(b0x°^5Fm/k], =\n[Fm/4kfvq0). (22)

Значения co!J, определяемые первой формулой (22), являются промежуточными по отношению к (13), (14) и хорошо согласуются с точными на весьма большом массиве значений Щ,а0. В течение времени tx, определяемом второй формулой (22), происходит основное приращение угла атаки а, (и числа оборотов N{).

11. В заключение отметим, во-первых, что в случае осесимметричного тела, когда вместо неравенства (1') выполняется равенство

A mZQ = 0, (23)

из (13) следует верный по физическому смыслу результат: со^ = со, т. е.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

невозможность возникновения авторотации, если спускаемое тело строго осесимметрично по форме и распределению масс.

С использованием метода усреднения и понятия линии центров колебаний угловой скорости случай (23) исследован, например, в [8].

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

только непосредственное, описанное здесь выше, но и для упрощения (как

в [2] — [4]) интегралов, входящих в правые части усредненных по времени уравнений вида (5), (5') или (9), (9'), см. выше: примечание к пунктам 6—8.

ЛИТЕРАТУРА

1. Шилов А. А. Влияние массовой и аэродинамической несимметрии тела на характер его пространственного движения//ДАН СССР.— 1968. Т..

183, №5. •

2. Л а р и ч е в а В. В., Ш и л о в А. А. Аналитический метод определения аналога сепаратрис при движении тела около центра масс в атмосфе-ре//Космич. исслед.— 1969. Т. 7, № 1.

3. Ларичева В. В. О разрывном характере аналога сепаратрис при движении асимметричного тела около центра масс в атмосфере//Ученые записки ЦАГИ.— 1972. Т. III, № 3.

4. Л ар и ч е в а В. В. Авторотация при полете в атмосфере асимметричного тела, вращающегося около центра масс в продольной плоскости//Труды ЦАГИ.—1915. Вып. Vm.

5. ЯрошевскийВ. А. Движение неуправляемого тела в атмосфере.— М.: Машиностроение.— 1978, гл. 5, стр. 161—163.

6. Л а р и ч е в а В. В. Эффективное преобразование и асимптотика одного класса систем нелинейных дифференциальных уравнений//ДАН СССР.— 1965, №2.

7. Ларичева В. В., Ефимов Г. Б. Асимптотика несимметричных колебаний маятника и эволюция разгонного движения точки в центральном поле притяжения/ЯТрепринт ИПМ им. М. В. Келдыша.— 1982. № 56.

8. Курьянов А. И., Ларичева В. В. Приближенные параметры перехода вращений в колебания около центра масс осесимметричного тела, входящего в атмосферу//Ученые записки ЦАГИ.— 1979. Т. X, № 6.

Рукопись поступила 31/Х 1997 г.

i Надоели баннеры? Вы всегда можете отключить рекламу.