Файл: Лебедев А.А. Динамика полета беспилотных летательных аппаратов учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 09.04.2024

Просмотров: 233

Скачиваний: 16

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Рис. 12.25. Частотные характеристики

передаточных функций b( p ) / ö ( p ) и д(р)/6(р) аппарата с поворотными крыльями (/( = 2,38; Т=0,10; £= 0,15; Гі = 0,084)

Рис. 12.27. Частотные характеристики передаточной функции а(р)/6(р) аппарата с поворотными крыльями (ACT =0,1; 7=0,10; £= 0,15;

Гб= 0,125)

Рис. 12.26. Частотные характеристики

передаточных функций Ѳ(р)/6(р) и Ѳ(р)/6(р) аппарата с поворотными крыльями (/( = 2,38; Т = 0,І0; 1= 0,15; Те =0,09; 1=0,05)

573

новятся аналогичными частотным характеристикам Ѳ(р)/6(р) и Ѳ(р)/Ь(р) аппарата с закрепленными крыльями (сравните рис. 12.25 с рис. 12.18 и 12.20).

Частотные характеристики передаточной функции Ѳ(р)/6(р) строятся суммированием характеристик множителей, входящих в

эту функцию. Характеристики множителей lip

и --------- —---------

^

у

и

Т2р2 + 2ІТр+\

были рассмотрены выше, частотные характеристики множителя 7Ѵр2 + 21нТѳр + 1 приведены на рис. 12.26. Они являются обрат­ ными характеристикам колебательного звена: при частотах, меньших 1/Те, асимптота имеет наклон, равный нулю; при часто­ тах больших 1/Гѳ,наклон равен +40 дБ/дек; при частоте изло­ ма 1/7+ослабление составляет 20 lg2g. При Щз-Ч), Т<$-+Т и частотные характеристики передаточной функции Ѳ(р)/6(р) при­ ближаются к характеристикам множителя КІр (см. рис. 12.26). Различия в динамических свойствах аппаратов с поворотными и закрепленными крыльями выявляются при сравнении рис. 12.26 с рис. 12.18 и 12.20.

Пример частотной характеристики передаточной функции a(p)/ö(p) аппарата с поворотными крыльями приведен на рис. 12.27 (сравните с рис. 12.18).

§ 5. ПЕРВЫЙ ЭТАП ПРОДОЛЬНОГО ВОЗМУЩЕННОГО ДВИЖЕНИЯ НЕМАНЕВРЕННЫХ ЛЕТАТЕЛЬНЫХ АППАРАТОВ

Приближенные уравнения для первого этапа возмущенного движения летательного аппарата получают, отбросив первое уравнение в системе (11.35) или (12.9) и положив А+ = 0 в остав­ шихся уравнениях:

dPда

ЦД9-

, '

db-a.

,

Щ2

« и —dt— Ь«і2—dt— + Я і2Да =

=

— «13 — — -------- «13

+ « 1 5 ^ 2 B I

a 42Д а -f « 44Д Ѳ —

dt

 

> (12.182)

dt

=

Д 8 «-j43-а 45Х в;

— Д & 4 " Да-[- Д Ѳ = 0.

 

 

Составим по известным правилам передаточные функции, со­ ответствующие уравнениям (12.182) *.

Характеристический полином имеет вид

л "(р )= Р3 + Я ’хР1+ А \ р + А \ ,

(12.183)

* Передаточные функции летательного аппарата для первого этапа воз­ мущенного движения можно найти также непосредственно из выражений _ (12.14). (12.36) и (12.39), приняв в них равными нулю все коэффициенты с индексами, содержащими «0», например, дог и а40.

574


где

А\ = d[i -f- Л]2~f“ Ö-42— Ö44;

^ 2 = Ö J 2 - f - аЦЙ 42 — ( a „ 4- f l i 2 ) Д44І

Л з= — а 12«44.

Реакция угла тангажа на отклонение органов характеризуется передаточной функцией

ßiP2 4 -в 2р + в'3

< ( /; ) =

Р2 + А-^р2 А2р 4“ А3

где

#і =

— а\ъ,

#43#12 '— « 1 3;

 

В 2 =

Ö13 (#42 — # 44)

 

Вз = Иіз(112 #j3 (#42'

#44)-

(12.184)

управления

(12.185)

(12.186)

Реакция угла наклона траектории на отклонение органов управления описывается передаточной функцией

 

w t i p y

С \ Р 2+ С 2Р

+ с3

(12.187)

 

Д3+ /^Д24"^2^ 4“Л

 

 

 

где

Cj =

Ц43;

 

 

 

£ 2 =

#43 (й ц -р -Я Н г ) — #13#42>

(12.188)

 

 

 

С 3 =

#43#12— #42# 13’

 

 

 

Реакция угла атаки

 

 

 

 

D\P24- D2P4- D3

(12.189)

 

Wt ( p) =

 

 

 

р 34- А хр 24- А 2р 4- А 3

 

где

P i— — a'i3—ais,

 

 

 

D 2

— Й 13“ Ь #13#44 — # 4 3 # ц !

(12.190)

 

 

 

Z)3 =

CZ13Ö44.

 

 

 

Реакция угла тангажа на возмущающий момент характери-

зуется передаточной функцией

 

 

 

 

£^д 4-е "2

(12.191)

 

W vM{p) = —

 

 

 

Р 3+ A ß 2+ А 2р 4- А з

 

где

Е і ~ Л 15\

I

(12.192)

 

Е 2 = аХ5{ап — аы).

)

 

 

575


Сравнивая передаточную функцию'(12.25) с упрощенной пе­ редаточной функцией (12.185), замечаем, что в результате пре­ небрежения вариациями скорости коэффициенты Л4 и Д4 стано­ вятся равными дулю. В результате этого один из множителей в знаменателе (тзр+1) и множитель в числителе (тір + 1), отве­ чающие медленному движению, обращаются во множитель р и сокращаются.

Коэффициенты А и А2 и Л3 несколько изменяются, так как из

выражений для них выпадают динамические

коэффициенты aik

с индексами, содержащими «О», например, а02

и п40. При этом па­

ра больших корней изменяется незначительно, в то время как значение оставшегося малого корня может заметно отличаться от прежнего значения.

Коэффициенты В2" и В2" теперь немного отличаются от коэф­ фициентов В2' и Bz, приближенно определяющих постоянные времени Т\ и Т2, соответствующие быстрому движению, за счет замены а42 в формулах (12.28) на а42— а44.

Как и прежде, корни характеристического полинома можно рассчитать довольно точно с помощью несложного приема. Старшие члены полинома (12.183) определяют пару больших корней, как корни квадратного уравнения

р* + А Ір + А І = 0 .

(12.193)

Два младших члена полинома (12.183) определяют малый ве­

щественный корень, приближенно равный

 

 

 

Рг =

 

 

(12.194)

Учитывая сказанное, передаточную функцию (12.185)

можно

представить в виде

 

 

 

 

W \ Iр) = к

- —

{ТіР +1} <Т2Р+ 1-}—

(12.195)

Ѵ

{ Т 2 р 2 + 2 ^ Т р + \ ) { х 2 р

+ \ )

 

Здесь К — передаточный

коэффициент

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

аппа­

рата, равный

 

 

 

 

К

 

 

(12.196)

т2 — постоянная времени медленного движения, рав

ная

 

 

 

 

*2

1

 

(12.197)

 

 

Рг

Т, I, Ті и Т2— параметры, характеризующие быстрое движение.

576


Пример использования передаточной функции (12.195) мож­ но видеть ниже.

Рассмотрим теперь передаточную функцию баллистической

ракеты Wb {р), составленную без учета вариаций скорости [14]. Выполнив вычисления, получим

W\{p):

' а 13Р + Ö43Ö12 + ^13 (ß44 — a ¥l)

+ a\\ aii) P2 +

+ a\2 al\^AA) P а 12а 4і

P 3

так как здесь+Ві(^42= Аі = 0, или в стандартной(a llÖ42

форме

w l ( p ) = K l -

TiP+l

(12.198)

 

 

(Т2р2 +

21Тр + 1) (т2/>+ 1)

В этом случае параметры Ти Т и £, характеризующие быст­ рое движение, по величине не отличаются от соответствующих параметров передаточной функции (12.47), а постоянная време­ ни т2 принимает другие значения, указанные на рис. 12.4 пунк­ тирной линией. На этом примере можно еще раз убедиться в том, что динамические коэффициенты с индексами, содержащими «О», практически не влияют на характеристики быстрого движения и при исследовании свойств последнего можно не принимать во внимание вариации скорости.

Если построить частотные характеристики для двух переда­ точных функций (12.47) и (12.198), то легко обнаружить, что при достаточно больших частотах эти характеристики, отвечающие одному и тому же моменту времени, совпадают, тогда как при очень малых частотах они могут существенно различаться. Так, в нашем примере эти частотные характеристики практически сов­ падают при значениях частоты ©>0,03 1/с. Если учесть, что час­ тота среза системы стабилизации в нашем примере близка к 10 1/с, то отсюда следует, что при анализе устойчивости замкну­ той системы стабилизации можно не учитывать вариаций скоро­ сти и использовать более простое выражение передаточной функ­ ции (12.198). Естественно, этот вывод может быть основан толь­ ко на анализе конкретных числовых характеристик.

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

Можно также отметить, что быстрое движение характеризует­ ся весьма малым затуханием (на большей части активного участ­ ка траектории коэффициент относительного демпфирования £<0,1), что вызывается отсутствием развитого оперения и боль­ шой высотой полета. Поэтому система стабилизации должна обес­ печивать хорошее демпфирование этого движения.

В заключение составим выражение передаточной функции W t (р ), характеризующей изменение угла Ѳ при отклонении ру-

577


лей высоты. Это выражение без учета вариаций скорости имеет вид .

Wt(p)

_____________ Й43І2 + ДцД43Р +

«12^43 — «42^13____________

Р3 + (Ö42 + «11 — а0о) Р2 + («Па 42 + «12 — 0ЦЙ44) Р ~

012^44

ИЛИ

 

 

Wl{p) = Kb

( П р

 

ті Р

 

 

 

(12.199)

 

 

(7-2р2 +

26Гр +

1) (т2р +

1)

 

 

 

+ 1) (

 

+

1)

 

 

 

Зависимость параметров К& , Т3 и Т4 от времени приведена на

 

 

 

 

рис. 12.28.

 

динамические

 

 

 

 

Основные

 

 

 

 

особенности

баллистической

 

 

 

 

ракеты

как

объекта

управле­

 

 

 

 

ния наглядно выявляются и при

 

 

 

 

рассмотрении ее частотных ха­

 

 

 

 

рактеристик.

 

 

 

 

 

 

 

Для

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

 

 

 

 

частотные

характеристики, от­

 

 

 

 

вечающие

передаточной

функ­

 

 

 

 

ции

(12.198), значения

пара­

 

 

 

 

метров

которой соответствуют

Рис. 12.28. Зависимость парамет-

рис. 12.4.

 

 

 

 

ров

передаточной функции W j

На рис. 12.29 показано се-

продольного движения от време-

мейство частотных характерис-

 

 

ни

 

тик для

 

различных

моментов

 

 

 

 

времени

полета. Обращают на

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

Рис. 12.29. Амплитудные и фазовые частотные характеристики ракеты для продольного движения

578