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

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

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

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

Добавлен: 09.04.2024

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

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

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

ние на множители квадратичных трехчленов зависит от распре­ деления их корней на комплексной плоскости, что зависит уже от динамических свойств данного летательного аппарата. Разло­ жение на множители квадратичного трехчлена вида ap2 + bp+l приведено в табл. 12.2.

 

 

 

 

 

Таблица 12.2

 

Разложение на множители квадратичного трехчлена äp2 + bp+ 1*

Корни трехчлена

Знак —Ь

Дискриминант

г

Множители трехчлена

 

 

а

 

 

 

Пара

комплексных со­

 

 

 

 

пряженных с отрица­

 

/>2 <[ 0

г < і

Т2р2 + 2-Тр 1

тельной

вещественной

 

частью

1

> 0

Кратные

вещественные

а

 

 

è2 4а = 0

$== 1

СГр + 1)2

отрицательные

 

 

 

 

 

 

 

Вещественные

отрица­

 

 

 

62 > 0

г > і

(Т'р 4- 1 ) ( Т " р + Ц

тельные

 

 

 

 

 

 

 

 

Нулевой

и веществен­

1

=

0

Р (Гр + 1)

ный отрицательный

а

 

 

 

 

 

Вещественные

отрица­

1

 

 

( Т ' р + 1 ) ( Т " р + 1)

тельный и

положитель­

<

0

 

 

7 '< 0

ный

 

 

а

 

 

 

 

 

 

 

 

 

 

 

 

 

* В отличие от

книги [14],

где все постоянные времени считались положительными,

и множитель,

соответствующий

положительному корню,

записывался в виде — {Т'р 1),

здесь постоянные времени рассматриваются как алгебраические величины и положитель­ ному корню отвечает множитель Т'р-И , где Т ' <0.

В дальнейшем для сокращения записи трехчлен, характери­

зующий

быстрое

движение, будем

записывать

 

в виде

Т2р2 + 2^Тр

+ \,

имея

в виду, что он может

распадаться

на два

множителя

в

соответствии с табл.12.2. Квадратичный

трехчлен,

соответствующий малым корням характеристического уравнения, условимся также записывать в общем виде т2р2 + 2 |т т р + 1, не за­ нимаясь детально его структурой. Таким образом, характерис­

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

 

А{р)=,А, (Т2р 2 + 2\Тр + 1) (т2р2+ 2Ъ%р+ 1).

(12.18)

Приближенный способ определения корней

Благодаря тому, что две пары корней характеристического по­ линома (12.17) сильно отличаются друг от друга по модулю, для их определения могут быть использованы различные приближен­

17

515


ные способы. Простейший из них состоит в решении дйух квад­ ратных уравнений:

р2-{- А гр-\- Л2= 0 ;

(12.19)

Л2/?2-|-Л3/?-{-Л4 = 0.

(12.20)

Первое из них определяет большие корни, второе — малые. Полученные такие путем значения корней можно при желании

уточнить [19].

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

Л (р) = /'4 + ^і/?3-Ь^2/с’2_г Л 3/ ? + Л 4 = 0

(12.21)

путем отбрасывания членов А 3р и Л4. Поскольку, положив Л3 = = Л 4 = 0, мы пошли на весьма существенное упрощение задачи, то чтобы быть последовательными, целесообразно соответствен­ но упростить и выражения для коэффициентов Лі и А%. Нетруд­

но видеть, что А3 и Л4 обращаются

в нуль,

если положить, что

<% = <*w= ° и

<*44 = 0,

(12.22)

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

 

 

Рэ+ А '1Р+ а ;= 0,

 

 

(12.23)

в котором

 

 

 

 

 

 

 

 

 

 

Лі -—

-j-

 

**42)

 

 

(12.24)

 

A<i= а.п -\- ana.ft.

 

 

 

 

 

 

 

Для иллюстрации сказанного приведем пример. Рассмотрим

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

на

высоте 12 000 м со

скоростью

800 км/ч. Коэффициенты уравнений (11.35)

имеют

следующие

значения:

аоо= 0,0134

1/с; а02= 9,35

м/с2; аоз=0;

а04= 9,81

м/с2;

аю= 0; ац = 0,49 1/с; а12 = 2,98

1/с2;

аа' = 0,22

1/с;

йіз= 2,32

1/с2;

аіз' = 0; а40 = 0,0007 1/м; а42 = 0,577

1/с; а43 = 0;

а44 = 0. В этом

при­

мере характеристический полином

 

 

 

 

 

 

 

р4-\- 1,299/73 -j- 3,278/?2 -|- 0,0455/? -ф0,0201

 

 

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

 

 

 

р і

2 = - 0 ,6 4 4 ±

1 ,687/; р ъ л =

- 0 ,0 0 5 6

± 0,0783/.

 

Корни квадратного уравнения (12.19)

 

 

 

 

 

р

- \ - 1,299/7+ 3,278 = 0

 

 

 

 

516


равны

р х 2 — — 0,650 ± 1,690/, а корни уравнения (12.23)

р2 + 1,287/7 + 3,263 = 0

равны

р12 = -0,6435 ;г 1,688/,

т.е. практически не отличаются от «точных» значений корней. Корни уравнения (12.20)

3,278/72+ 0,0455/? + 0,0201 = 0

равны

. />3)4= —0,0069. ± 0,0780/.

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

2.3. РЕАКЦИЯ УГЛА ТАНГАЖА НА ОТКЛОНЕНИЕ ОРГАНОВ УПРАВЛЕНИЯ

Определение реакции угла тангажа на отклонение органов управления необходимо при проектировании и исследовании сис­ темы стабилизации, в которой при помощи свободного гироскопа измеряется отклонение угла* тангажа от некоторого заданного значения. При помощи такой системы может быть реализован полет по программе (см. пример в гл. II, § 7).

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

го гироскопа. Такие системы

стабилизации

получили широкое

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

(см. [14]).

 

 

Реакция угла тангажа на отклонение органов управления опи­

сывается передаточной функцией (12.14)

=

W

?

(

В'РІ±o .^РІ) ± -В*Р= + —

 

 

P4

+ A ip s +

A2p 2 + АъР + A4

А (p)

в которой коэффициенты A lt ..., Л4 и В\, ..., ß 4 определяются фор­ мулами (12.11) и (12.15).

* Эта система уравнений приведена ниже, в § 3, см. (12.49).

517


Для приближенного разложения на множители числителя пе­ редаточной функции (12.25)

В ( р } = В іР3+ В 2рч- + В 3р і - В 4

(12.26)

воспользуемся свойствами продольного возмущенного движения. Корни полинома В(р), соответствующие быстрому движению, найдем, положив в выражениях (12.15) коэффициенты аОІ, ам и Ö44 равными нулю. Это будут корни квадратичного трехчлена *

В'(р) = В іР^ В 2 р - { - В І

(12.27)

так как В4 обращается в нуль. Здесь

 

 

+ = — an -,

 

 

В ч —

— # + < 7 42 - J - а 43а і2 —

а 13;

(12.28)

 

 

Вз — — 0 - 1 3 ^ 4 2 “ |~ ^ 1 2 ^ 4 3 -

Корни полинома (12.26) сильно отличаются друг от друга по модулю: он имеет два больших и один малый корень. Поэтому старшие члены этого полинома

ВіР2-\- В 2р-\- В 3~ В гр2-\- Вчр-\- Вз

(12.29)

определяют два больших корня, соответствующих быстрому дви­ жению, а младшие члены полинома (12.26) Взр + В^ определяют малый вещественный корень, соответствующий медленному дви­ жению. Этот корень приближенно равен

Р г = ~ ^ т -

(12.30)

 

Рассматривая пример с тем же реактивным самолетом, най­ дем, что полином в числителе передаточной функции (12.25) име­ ет второй порядок, так как ßi= 0 :

-(2,32/-+1,370/7 + 0,01868),

причем его корни равны рі = —0,576 и р3= —0,0140. Приближенное значение малого корня

В±

0,01868

0,0136.

Рз~-

1,370

В з

 

Приближенное значение большого корня находим из урав­ нения

2,32/7+ 1,370 = 0,

* Квадратичный трехчлен (12.27) совпадает с числителем передаточной функции (12.78), полученной в предположении, что ЛУ и а44 АѲ равны нулю. Разложение этого трехчлена на множители представлено в табл. 12.4.

518


т. е.

Рі =

1,370

0,590.

2,32

 

 

Учитывая сказанное выше, полином в числителе передаточ­ ной функции (12.25) можем разложить на множители следую­ щим образом:

в (р)= В 4(7 > + 1)( 7 > + 1) (тіР + 1),

(12.31)

где Ті и Т2—’Постоянные времени, соответствующие быстрому движению, а ті — постоянная времени, отвечающая медленному движению *.

Таким образом, передаточная функция летательного аппара­ та (12.25) может быть представлена в виде

 

W ѵ(р)

К &____ (^іР + 1)(Т2р + 1) (тцр +

1)

(12.32)

 

5 { Р ) ~

S

( Т 2 р 2

+ 2$7> +

1) { Х 2 р 2

+

26ХТр

+ I)

 

если

В г= — а'пфО,

или в

виде

 

 

 

 

 

 

wl{p) = K

__________________ ( Т 1р +

1) ( У \ р

+

1)__________________ _

(12.33)

 

( Т 2 р 2 + 2 $ 7 > + 1) { х 2 р 2 + 2 5 тТ р

+ 1) ’

 

 

 

 

если

В х= — аіз =

0.

 

 

 

 

 

 

 

В этих формулах

 

 

 

 

 

 

 

К \ = - ^ - — передаточный коэффициент

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

аппара-

 

А4

 

 

 

 

 

 

 

 

та; Ті,т, |т — параметры, соответствующие медленному движению;

Ті, Тг, Т, § — параметры, отвечающие быстрому движению. Для рассматриваемого реактивного самолета:

К І = - 0,929; Т = 0,544 с и £=0,357; т = 12,73 с и $т = 0,0713;

7 ^ = 1 ,7 62 с ; 7'2= 0 ; т1= 71,43 с .

В качестве второго примера приведем результаты разложе­ ния на элементарные множители передаточной функции (12.25) для баллистической ракеты типа Ѵ-2 на 50-й секунде активного полета:

 

w\{p)=Kl-----(ТД’+іНѵ+В----_

(12.34)

 

(Г2р2 +

2 £ 7 > + 1 )(т 2/ > + 1 )(тз/> + 1 )

 

 

Здесь

/Сг— —4,125;

Г =

0,23 с; $ = 0,06; 7’1=

3,4 с;

т2 =

= - 210

с; т3= 100 с;

т1 = 75 с.

 

 

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

тательного аппарата характеризуется постоянными

времени, из-

* Постоянные времени Тх и Т2 приближенно определяются выражениями (12.94а) и (12.97а), а т — формулой т = —Л/р3.

519