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

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

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

угловой скорости касательной к траектории и угла атаки, рав­ ные

» у с т = 0 усг = К К ;

(12.63)

a yCT = K te8B.

(12.64)

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

Ѳта* = К К тз*

[(12.65)

и тем выше маневренные свойства летательного аппарата.

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

угловую скорость Ѳ, а нормальное ускорение _/„= РѲ или нор­

мальную перегрузку пу ~ — Ѳ. Тогда рассматривают переда­

точные коэффициенты летательного аппарата по нормальному ус­ корению

( Щ

= К Ѵ

(12.66)

\ 8В /у с т

 

или по нормальной перегрузке

 

 

(ИЯ-\

= К — .

(12.67)

I 6„ /у с т

Г *

 

Очевидно, что максимально возможная установившаяся пере­ грузка, т. е. располагаемая перегрузка, равна

Му расп ~ К ~ ~ ЧВт а х •

( 1 2 . 6 8 )

Располагаемая перегрузка наряду с максимальными значени­ ями нормального ускорения jymax — КѴбвтах И уГЛОВОЙ СКОрОСТИ

касательной к траектории Ѳт ах= 7СбВтах определяет маневрен­ ные возможности летательного аппарата.

Выясним кстати влияние конструктивных и аэродинамических параметров летательного аппарата, высоты и скорости полета на маневренные свойства аппарата. Обычно при достаточной степе­ ни статической устойчивости произведение коэффициентов аца& оказывается в несколько раз меньше коэффициента а12. Поэтому

532

можно приближенно записать, что

 

 

 

 

1

Г X т ь

I

 

/С ~ я 43

а 42

=

Г 5 -------£.(/>_( - Г )

;

(12.69)

mV

 

а п

 

 

К

 

 

 

 

 

 

 

 

 

г

ѵ

1

 

Y s

т* ( Р 4 - Г а)

 

(12.70)

к і

^

 

 

 

<

 

 

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

тойчивости (т. е. с увеличением по модулю коэффициента ml) маневренные свойства летательного аппарата снижаются (про­ тиворечие между статической устойчивостью и маневренностью).

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

/ C t e Ä

- -ms- .-

^

(- 1 2 .=7

ап

т\

 

 

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

аппаратов (m“<C О) схемы «утка» для создания положительного угла атаки требуется отклонить органы управления на положи­

тельный угол, так как в этом случае ml^> 0 и / — ) )>0.

Чтобы сбалансировать статически устойчивый летательный аппарат обычной схемы или схемы «бесхвостка» * на положи­ тельном угле атаки, требуется отклонить органы управления на

отрицательный угол, так как у аппаратов этих схемтг<С0 и по­

этому

< 0 .

\о /у с т

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

* «Бесхвостка» является разновидностью обычной схемы. Между этими схемами нет принципиальных различий. Однако у «бесхвостки» обычно нельзя

пренебрегать подъемной силой органов управления Y 8 6 , тогда как у нормаль­ ной схемы часто возможно допущение, что У86 =о.

18— 3422

533

ш ')> 0 и [ — ) > 0 . В этом случае схема с поворотными

крыльями может рассматриваться как разновидность «утки».

Другие возможные варианты

схемы с

поворотными

крыльями,

когда т\= О (идеальная схема) и

mz<^0 (разновидность

обычной схемы), значительно

уступают первому по

маневрен­

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

4.2. ПЕРЕДАТОЧНЫЕ ФУНКЦИИ ЛЕТАТЕЛЬНОГО АППАРАТА

Общий случай

Будем рассматривать одновременно продольное движение (движение тангажа), описываемое системой уравнений (12.50), и движение рыскания, описываемое аналогичной системой урав­ нений [см. уравнения (13.10)].

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

Передаточные функции

Для получения передаточных функций летательного аппарата запишем систему уравнений (12.50) в операторной форме:

РІР + аи) ^ + {ра'і2+ аи) а = ~ { р а [ 3-\-а^)Ь-,

рѲ — а^а — а^Ь;

(12.72)

б — Ѳ — а = 0.

Чтобы записать соотношения между изображениями б(р), Ѳ(р), а(р), б(р) параметров б(£), Ѳ(Т), a(t), ö(t), разрешим систему (12.72) относительно -б, Ѳ и а. Тогда получим

Ъ(р) = —

' ѳ ( р ) = — . а ір)-

(12.73)

Д

А

 

где А— главный определитель системы (12.72):

 

РІР + ап)

0

ра'п + ап

 

Д =

0

р

— а42

 

' Р

1

- 1

- 1

(12.74)

І Р 2 4~ (а и +

а 12 -J- а 4г) Р-~\~ а і2 а п а 4і]>

534

А®, Аѳ, До — присоединенные определители, полученные из глав­ ного заменой соответствующего столбца столбцом, составленным из правых частей системы (12.72): •

(/^із-Ьа із) 8

О

ра' п - \ - а п

 

Да=

<3:438

Р

— аі2

 

= [а \ ъ p2Jr (аіз +

аі'3«42 — а ' п а ^ )

— 1

- 1

а 12а 43]8;

(12.751

р - \ - а іг а і2 —

Д н [ а 43/?2 (я цЛ^з -j- а і 2а 43 d i 3u 42) р

-j- <212^43

^ 13^ 42] 8>

(12.76)

Ae =

[(^43-f- 0:13) Z7~j- ^ 11^43Ч- ^ 13] 8.

(12.77)

Пользуясь выражениями (12.73) — (12.77), составим переда­ точные функции *. Например, для передаточной функции ®(р)/й(р) получим формулу

а (р)

__ — « 1 3 Р2 — (#13 +

#13#42 — # 12#4з) Р — #13#42 +

# 12#43

(12.78)

8 (Т')

Р

2

+ (# и

+ а[2 +

<342) Р +

а 12 +

#іі#42]

 

или

 

 

 

 

 

 

 

 

 

 

 

» (р )

__“

« 1зР2 — (#із +

# і3#42 # і2#4з) р а 13а4 2 + 0 1 2 0 4 3

(12.79)

8 (Р )

 

^2

( д ц

+

aj2 + #42) Р + #12 "Т #11#42

 

 

 

 

Пусть ßi2+ 011042#

0. Тогда

(12.78)

можно переписать в виде

 

 

 

*

 

 

 

г

г

 

 

 

 

 

 

“ 13

 

 

— #13 — #1з#42 + #12#43

Р + 1

 

8 ( Р ) = к

#12#43 — #13#42

■Р2 +

#12#43 — #13#42

 

 

 

 

 

-, (12.80)

&(Р)

 

 

 

1

 

Р2 +

#11 + # 1 2

+ #42

+

1

 

 

 

 

 

 

 

j

-

 

 

 

 

_ #12 + #И#42

#12 + #11#42

 

 

 

где К —-передаточный коэффициент летательного аппарата, оп­ ределяемый формулой (12.59).

* Эти передаточные функции можно найти либо обычным способом из уравнений (12.50), либо из выражений полных передаточных функций (12.14), (12.36) и (12.39), положив в них коэффициенты с индексами «0» и коэффи­ циент ß44, равными нулю.

18*

535

Аналогично получим

а 43

ß l l a 43 + й 12а 43 — Ö13Ö42

p 2

4 .------------------------ /7 + 1

ѳ » )

К

а -п Ч ъ а \за 42

________ д 12д 43 — #13ö 42

 

 

( 1 2 . 8 1 )

8(Р)

 

 

 

 

 

1

all + ai2 + a42

 

 

 

 

 

 

;

+

#12 +

p

+

1

 

 

 

 

a 12 + a 11^42

a 43 +

a l l a i2

 

 

 

 

 

 

__________

a i3

P + 1

 

 

 

 

а (P)

 

' a 13

 

 

 

 

H P )

Kxa

Q lla 43 +

______________________________

(12.82)

 

 

 

1

p2

a l l ~h a i2 ~i~ a 42

P

1

 

 

 

 

a 12 + a l l a 42

a 12 + ^11^42

 

 

 

где Kr а— передаточный коэффициент летательного аппарата по углу атаки, определяемый формулой (12.61).

Пусть за выходные величины принимаются нормальные пе­ регрузки. Чтобы составить передаточную функцию ny(p)fb(p), найдем связь между изображениями пу(р) и Ѳ(р). Рассматри­ вая V/g как число и переходя в (12.57) к изображениям, по­ лучим

п у ( р ) = — р Ѳ { р ) .

 

 

&.

Отсюда

 

 

Пу (Р) _

V

ѳ (р)

8 ( р )

g

Р 8 ( р )

Характер свободного возмущенного движения

Приравняв нулю знаменатель передаточной функции (12.80) или (12.81), получим характеристическое уравнение звена, у ко­ торого входной величиной является угол отклонения органов уп­ равления б, а выходной — угол # или Ѳ:

р { - - - !- - -

Р 2-}-а и ~~—- +—

2- р - \ - 1) = 0 .

(12.83)

Vо .12 + ЙЦ042

^12 + а 11я 42

/

 

Нулевой корень этого уравнения соответствует медленному возмущенному движению. В самом деле, в результате упроще­ ния системы уравнений продольного движения (11.35) путем от­ брасывания приращения скорости АР порядок системы и ее ха­ рактеристического уравнения (12.3а) понизился на единицу. При этом один из малых корней характеристического уравнения вы­ пал, а другой обратился в нуль. Благодаря нулевому корню при­ ращения § и Ѳ в свободном возмущенном движении, описывае­ мом уравнениями (12.72), не затухают.

536

Корни квадратичного трехчлена

а р 2 - \ - Ь р - \ - \ =

 

1

, ß - r

an . :''~

а,2- + д.42-

р -

|_

1 (12.84)

--------а 12

 

а---------\\а 42

 

а 12

 

 

+

 

+ Лц042

 

 

 

 

 

 

 

 

 

 

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

ной — или + или Ѳ, или ос:

1

Р2 + а и + ^12 4" а 42 /7+1 = 0.

(12.85)

а 12 + а 11а 42

а Ѵ2 + а П а 42

 

Состав сомножителей в выражениях передаточных функций (12.80) — (12.82) зависит от аэродинамической компоновки ле­ тательного аппарата, а именно от его аэродинамической схемы (обычная схема, «бесхвостка», «утка», схема с поворотными крыльями), свойства устойчивости (устойчивый, нейтральный, не­ устойчивый летательный аппарат), а также от степени стати­ ческой устойчивости и демпфирующих свойств аппарата.

Чтобы разложить на элементарные множители знаменатели передаточных функций, необходимо проанализировать квадра­ тичный трехчлен

, N

 

» , а П + а \2 +

а 42

 

/1 О Q С Л

а ( р ) = -------------- р

Ч

--------^----------Р + 1,

(12.85а)

а 12

1а 11а 42

а 12

+

а 11а

12

. .

 

 

+

 

 

 

 

 

 

приняв во внимание, что у любого летательного аппарата коэф­ фициенты ап, й\2 и Ü42 положительны, вследствие чего всегда ап + 0 і2/ + й42>0. Устойчивость летательного аппарата в быстром движении определяется поэтому знаком члена аі2+ аііа42-

При аі2+ ац042>0 летательный аппарат

устойчив, при аіг +

+ апа42=0 он нейтрален, а при 0 і2+ пца42<0 — неустойчив.

Устойчивый летательный аппарат.

 

Обычно

 

а 12 + а 11а42> 0.

(12.86)

Это условие выполняется, если летательный аппарат обладает статической устойчивостью (ßi2> 0), точнее говоря, когда вы­ полняется неравенство

“г

 

ТТІ

А

(12.87)

m ^ < - —

(ср + 57,3с“),

где ц — относительная плотность летательного аппарата:

(12.88)

р ^ 7 ;

537

Cp — коэффициент тяги:

2Р

Cp(12.89) pV2S

При условии (12.86) коэффициенты характеристического уравне­ ния (12.85) положительны и его можно представить в виде

7 Ѵ + 2 5 7 > + 1 = 0 ,

(12.90)

где Т — постоянная времени летательного аппарата:

Т = - --------- с;

~Ѵ а12+ а\\аА2

I— относительный коэффициент демпфирования:

а1\ + а 12 "I" а 42

2 У а-і2 + 0цв42

Корни характеристического уравнения (12.90)

(12.91)

(12.92)

Л,2= у - М

+ У Т ^ )

при выполнении условия (12.86)

всегда имеют отрицательные

вещественные части, и переходный процесс относительно коорди­

нат ö, Ѳ и а является затухающим.

Характер переходного процесса устойчивого летательного ап­ парата (аі2+ аіі<242> 0) определяется относительным коэффици­ ентом демпфирования (12.92). Пренебрегая произведением ацйщ по сравнению с а&можем записать

f

0 1 1

а 42

 

2 ~Ѵа,\ч

«ли

 

Р + Y a

М** + M az

mV

I z

 

 

+ mа

cp + 57,3c“ —

 

z

 

р£6д

 

 

2

m аz

2 m

7' Z2

 

 

 

538

где п — безразмерный радиус инерции:

h

г?

тЪ\

Отсюда следует, что коэффициент относительного демпфиро­ вания g мало зависит от скорости полета и существенно умень­ шается с ростом высоты полета. При возрастании степени стати­ ческой устойчивости коэффициент £ уменьшается.

П ри |< 1 п е р е х о д н ы й п р о ц е с с к о л е б а т е л ь ­ ный. В этом случае корни характеристического уравнения (12.90) являются комплексными сопряженными:

/>1 ,2 = 7 -

Постоянная времени Т пропорциональна периоду колебаний и обратно пропорциональна их частоте. Собственная частота, определяемая формулой

ша = у 1/с

или

/а = 2пТ Гц,

является важной динамической характеристикой летательного аппарата. Влияние конструктивных параметров аппарата, а так­ же параметров движения V и Н на собственную частоту легко проследить с помощью приближенной формулы.

Пренебрегая в (12.91) произведением апа42 по сравнению с ко­ эффициентом йі2, характеризующим статическую устойчивость, получим

V «12

Учитывая, что

maz— тІусУс

Са

 

У '

можно записать

Т = л [ ------ -----------

= і /

----------—^ ----------

. (12,93)

V 57,3mazqSbK

V

57,3 ( х т— х р ) c*qS

 

Отсюда следует, что постоянная времени Т зависит главным об­ разом от момента инерции, т. е. размеров летательного аппарата,

539

степени статической устойчивости, скорости и высоты полета. Конструктор может получить нужную величину Т, меняя цент­ ровку и положение фокуса летательного аппарата.

У летательного аппарата параметры его передаточной функ­ ции Т и I могут сильно изменяться вдоль траектории. В каче­ стве примера на рис. 12.5 и 12.1 приведены величины Т и £ для снаряда «Эрликон».

І

Рис.

12.5. Изменение

относитель­

Рис. 1 2 .6 . Изменение постоянных

ного

коэффициента

демпфирова­

времени

снаряда «Эрликон» в те­

ния снаряда

«Эрликон» в течение

 

чение полета

 

 

полета

 

 

 

При £ =1

п е р е х о д н ы й п р о ц е с с

а п е р и о д и ч е с к и й

(корни квадратичного трехчлена кратные вещественные, отри­ цательные); трехчлен (12.85а) в этом случае имеет вид

а(р) = ( Т р + \ ) \

причем постоянная времени Т определяется по-прежнемѵ форму­ лой (12.91).

При g> 1 п е р е х о д н ы й п р о ц е с с

а п е р и о д и ч е с к и й

(корни трехчлена

(12.85а) вещественные, отрицательные); трех­

член а(р) разлагается на два множителя

 

'

а(р) = (Т 'р+ 1)(Т"р+ \),

где

 

 

г = т { \ - \ - \ г¥=Х)\ Т " =

Т { \ - У ^ Г).

Неустойчивый летательный аппарат.

Рассмотрим случай, когда

 

а 12+ а Па 42 <

О,

540

т. е.

тУ>- -(сР + 57,Зсі).

Теперь коэффициенты а и b квадратичного трехчлена (12.84) от­ рицательны и уравнение (12.85) имеет два вещественных корня

Ь ± V 62 — 4д

Рх, 2

2а

 

из которых один отрицателен, а другой положителен. При этом трехчлен а(р) распадается на множители

а ( р ) = ( Г р + \ ) ( Т " р + \ ) ,

где

Г = ^ ( - Н - у т а ) < 0 , T" = ± { b + y W = T a ) ,

1

а 1Х

+ а \ч + ß 42

a X2 + ß l l ß 42

ß 12 + ß l l ß 42

Нейтральный летательный аппарат: 012+ 011042 = 0 ( у трехчле­ на а(р) один корень нулевой, другой вещественный, отрицатель­ ный); трехчлен а(р) имеет вид

а ( р ) = р ( Т р + 1),

где

 

 

Т =

--------]---------

.

а ХХ + ß 12 +

ß 42

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

ß 12ß 43 — ß 13ß 42

1/с2;

К - -

 

 

а ХХ +

ß 12 4"

ß 42

 

К т .=

ДіЗ +

ß l l ß 43 - 1/c

a ll + ß 12 + ß 42

и другой физический смысл:

 

 

 

К =

уст

 

уст

 

 

*Нт)/ уст

Переходный процесс нейтрального летательного аппарата от­

носительно координат fl, Ѳ и а является апериодическим зату­ хающим.

541