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

книги из ГПНТБ / Динамика полета и конструкция крылатых летательных аппаратов

..pdf
Скачиваний:
52
Добавлен:
30.10.2023
Размер:
21.98 Mб
Скачать

ветра будем пренебрегать. С учетом указанных допущений из

(1.46)

второго

раздела имеем

Vncos 0

Р cos (ft + Та) Р 8 + ^sin — a) sin Р з

§

 

 

 

 

— Q cos (0 — a) cos p sin p3 -f Z xcos p3.

Из фиг. 3.2 видно, что сумму сил

Y sin (& — a) sin рз — Q cos (ft — a) cos p sin (33

Ф и г . 3 . 2

можно заменить суммой

 

{ Y sin 0 — Q cos6) sin(P„3 — p3).

 

Тогда

 

 

 

 

 

 

 

 

 

Vncos 0„

at

= — P cos (9-+ Та) sin Р з + (Y sin 0 —

 

g

 

 

 

 

 

 

 

— Q cos 6) sin (ря3 — p3) + Z xcos p3.

(3.5)

Для

линеаризации (3.5)

будем

полагать,

что

Р = const, . Y =

const,

Q = const,

 

0 =

const,

0„ = const,

& =

const, Vn= V 0.

Тогда, имея в виду, 4fo

 

 

 

 

 

 

 

 

Фл =

ФлО +

Д Фн > ■

 

 

 

 

 

Рлз

РлзО "t" Д Рлз >

 

 

 

 

 

Рз —

Рло +

Д Рз 1

 

 

 

 

 

«1 =

2Г10 + Д21,

 

 

370

после преобразований уравнение (3.5) в отклонениях может бытьпредставлено в следующем виде:

Цг Vt ^

= - Р„— COS 9

cos Р|»оДрлв - (Qo *os 0„ + F0 sin 0о) х

g

rtO

 

X-C°S(P-

Рло) (Дряд - Д|У -

Z10sin M P« + c o s

iW

 

cost)

 

 

 

 

 

 

 

В общем случае

 

 

 

 

 

Поэтому

 

■Zi=/(P, s*.-.-)-

 

AZX=

Zj0Ap + Z°* Дон + AZ,C.

 

 

 

 

р и

Из формул (1.38), (1.39)

и (1.41) второго раздела

при

следует,

что

 

(cos g)o

 

 

 

 

 

ДРЛ

ДР„>

 

 

 

 

COS^o

 

 

 

 

ЛР я ~ ДР + - ^ -

 

Откуда

 

. .

c o s & 9

д а

 

 

 

 

др -

 

 

 

 

"'"(cos а)е

Рпз

Ущ

 

(3.6)

(3.7)

малых

(3.8)

Следовательно,

(3.7)

с учетом (3.8) принимает вид

 

 

iZ l =

^ « - ( ^ 7 4Р"

- 4> Щ

+ Z ‘f

i8« + 4Z i.

(S-9)

Кррме того, легко определить разность

Дрял —Др,, если иметь

в виду,

что

 

 

 

 

 

 

 

 

sin Рдз =

cos60sin р3+

 

cos &0 1/0

 

Откуда,

принимая

 

 

 

 

 

 

 

 

 

Р п л — Р л л О + Д Р п з >

 

 

 

 

 

 

 

Р л —

РлО + Д Р л

 

 

 

 

 

 

 

(cos.a)0 cos 0О— 1

О,

 

 

имеем

 

 

cos &,

0

 

 

 

 

 

 

_ (cosa)0

 

ПК;

 

 

 

 

 

 

 

(ЗЛО)

 

 

ДРял-ДРл

cos 8-0 cos P„j0

V0

 

 

 

 

 

24*

37]

Наконец, пользуясь

фиг.

3.3,

найдем

связь между

Др„.,, Дфя

и .16. Из фиг.

3.3 следует,

что

 

 

 

 

откуда

 

Ф в ^

Ф

+ Р я в .

 

 

 

ЛРпл =

АФ,, — ДФ-

 

(3.11)

 

 

 

Теперь, подставляя

значения

Ъ-Z^ из (3.9), (дрпя — ДА,)

из (3.10)

и Др„3 из (3.11) в (3.6),

имеем

 

 

 

 

d ^ n

 

 

 

 

 

WC W, + L

 

dt

= а1пдФп+ aj. ДФ+^1"Д5 + а"

 

ГП

тп

 

у п

V

 

где

а*п = -М—

р

c o s (&„ + ?,?)

c o s в _

z io

i

Ч

GV0

р 0

eos»„0

c°sP*>

cosOn0

 

 

 

1

Z?o

 

COS O p

~

 

 

 

'

cos 0nO

^,M0 cos(a)0J ’

-аф" 4 ’

Л11

cosP" ''’

 

 

 

rnЙ - o l y

 

 

 

g

(Qo cos 0O

Y Cin ft 4

C0S (РплО - P * > )

COSp

GVj

 

0

°'

. COS0„o

COS&pCO

+COS0„o C0S |3",°.

AZ'«COS4o-

372

Для решения задачи еще необходимо иметь уравнение мо­ ментов в проекции на связанную ось Оуу х и две кинематические связи:

Л -fif = Му+ муР,

 

 

-Ц- =

Ш,С08»,

(3.13)

 

 

dl

„ • ,

.

 

 

 

sin

Тогда имея в виду, что

 

 

 

т у= м

др + му шу + Мь"Доя+ Ш ув,

J

 

>ч»

J

Уо .

дм,уР~ м;родр + м у0и)у ,

уравнения (3.12) и (3.13) в отклонения с учетом (3.8) в оконча­ тельной форме имеют вид:

рЦп =

aln%

+

а* ДФ+

Д8К+ ajfc W : 4-Д

 

 

 

p w y а*п Дф„ + а “ у Шу + а ' ^ Дф + а °« Д5„ +

 

+ / Иу

(3.14)

=

а “ушу

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

р&(, =

а*'1Дф„

 

 

 

 

 

 

 

 

 

 

принять

cos (0О+

уд)

 

1

 

.

cos &0 cos рПл0

Если

----- ------------~ 1 ,

cosOn0~

l , ------ й— — ---- ~ 1

 

v

 

 

COS 0ло

 

 

 

>

j

 

> cos Qno

cos a o

COS P»30 — РзО

cos a

 

то

 

коэфф,щиенты при

не"3'

--------- —

со8 8,совр „.;= !1'

 

вестных в системе (3.14) определяются

следующими

формулами

(индекс

„нулик1 опущен):

 

 

 

 

 

 

 

 

1. а5” = сд 7 P c o sP+ Z i - Z j s i n l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

PV

 

 

=

-

-

g y P c o s

P + -

§ -

(

c l - ^

sinP )-2~5;

 

2■ a* =

7 ^ 7

(P co sP -Z j*

+

Zj sin p) =

-

aj" ;

 

 

 

G V

 

 

 

 

 

 

 

 

 

373

3.

V

 

g

7*h- 8

c 6„p—

S '

 

 

 

 

Q V Z '

G C*‘ 2

^

 

4.

Wr

 

[z\ — Q cos 0 — Y sin 6) =

 

a,

с

GV2

 

 

 

 

 

 

 

 

 

 

 

 

=

-C-(c1!

c^cos 0 — cv sin 0) -4-5.

 

 

 

 

O

'

 

2,

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

ct и

zi

обычно задаются как функции числа М.

 

 

 

 

 

5'

 

 

.m*y ? v L s i - M l P ] ,

 

 

 

 

V

 

 

 

 

«• < ’ = т, ( « ? + = i ( и> * 7 - SP +

где Л1"ур = Л ^ р ,

 

 

7.

 

 

 

 

8 -

ak = 7 ; M >"= i m>‘ ^ r S!-

 

 

9‘

 

= Ту ( “ м у ~ ЖУ ~v

а “у ’ ■

 

10. . a£y=cos&,

11.а\п — — Vcos^n,

12.Д п = - A Z Iecosp,

13.Д у —&М.у>а -f- ДуИу,р« •

Производные m?, т “у т°н задаются как функции числа М.

'Для ракет с крестообр'азно расположенными крыльями и аэро­ динамически симметричных баллистических ракет

Р

a

т\, — т, ,

О)

ш_

туутгг >

тун = т1*»

374

При расчете коэффициентов уравнений (3.14) для этих объектов необходимо в формулы подставлять L вместо I и 2L2 вместо I2.

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

Следует отметить, что уравнения (3.14) дают точное решение для переходного процесса в боковом движении при угле крена, рав­ ном нулю. Это условие близко к реальному, так как; автоматически управляемые летательные аппараты, у которых изменение угла кре­ на заметно сказывается на характере переходного процесса, стаби­ лизированы по крепу. Поэтому решение системы (3.14) вполне точно 'воспроизводит переходные процессы, типичные для таких ле­ тательных аппаратов.

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

§ 3 .2 . П Е Р Е Д А Т О Ч Н Ы Е Ф У Н К Ц И И В Б О К О В О М Д В И Ж Е Н

В качестве примера входных воздействий для движения кре­ на рассмотрим угол отклонения элеронов 5Э и угол скольже­ ния р. Тогда линеаризованное уравнение движения крена при­ мет вид

(3.15)

чт<^ можно записать также в виде

Л~( == W j (р) Д р + И ^ ( / 7 ) А о э ,

(3.16.)

где W?(р) — передаточная функция угла крена по углу скольже­

ния, Wx>{p) — передаточная функция угла крена по углу откло­ нения элеронов, и

(3.17)

375

Здесь

 

*? =

 

 

 

 

 

 

 

 

 

Т7 =

— А.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tf' ^ =

F ( 7 ^ T T y

 

 

 

 

(3.18)

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

а э

 

 

 

 

 

 

 

 

 

м.г

 

 

 

 

 

 

 

 

 

ш

 

 

 

 

 

 

 

 

 

а *

 

 

 

 

 

 

Звено

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

 

 

 

функцией

вида

(3.18)

 

 

 

имеет

логарифмическую

 

 

(D

амплитудную

характери-

 

 

стику,

показанную

на

 

 

 

фиг. 3.4.

 

 

 

 

 

 

 

го

В уравнениях плоско­

 

 

 

бокового

движения

будем считать угол отклонения руля

входными

воздействиями

направления

8К и

боко­

вой ветер W(,.

 

 

 

случае

принимает

вид

Система уравнений движения в этом

рЦ п = а%»Ьй„ +

а* Д«1>+ с°«До„ -f

W

 

 

 

 

*п

' П

■п

 

п

 

 

 

(3.19)

р « у = а!» Дфл + а“ш у +

а*уЦ +

 

Зк +

 

W ,

 

рД'!)=а“уо)у

 

 

 

 

 

 

 

 

 

Откуда, перенося все члены, кроме возмущений,

в левую часть

и исключая (Oj,, получаем

 

 

 

 

 

 

 

 

(3.20)

376

Главный определитель системы (3.20) определяется из очевиД' ного соотношения

 

ф

 

 

 

 

Р ~ Ч п

 

Ч

 

дО) =

•п

 

(3.21)

ф (1)

 

 

агпа

у

р 2 - аУр

 

 

ы»уыф

 

 

 

 

У

 

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

ак

 

й\Ф

 

 

 

 

 

К

 

 

 

а « а , у

Р 2 — a j p

 

 

 

у

 

У

 

Д8„ +

 

 

 

 

Ч

р )

 

 

а.Wr

 

 

 

 

 

~

Ф

 

 

 

 

 

 

а \

6

ш.

 

a " V

 

у

 

 

 

*

 

у

*“уа<ь5

 

юу

 

 

 

w.

(3.22)

+

 

 

Ч р)

 

или

 

 

 

 

 

 

 

 

 

 

 

(3.23)

 

 

 

 

 

 

 

а г н

 

 

 

 

«7-

 

р - ч п и д>

 

 

 

р - ч

 

ч

 

 

 

 

 

 

 

Дф =

или

ф

ш

а « а , у

а;п а,, у

Wr

<»„

— а

л а , у

а , ‘■ а,,1У

у

ш у т

у 1

“ у

W Q, (3.24)

 

*(р)

д 2 „ +

 

 

Чр)

 

 

 

Дф =

^ : « Ы д5«+

W p ( p ) W b

(3.25)

Найдем расчетные выражения для

передаточных

функций,

а) П е р е д а т о ч н а я ф у н к ц и я п у т е в о г о у г л а по у г л у

о т к л о н е н и я

р у л я

н а п р а в л е н и я

W°H(р).

 

Раскрывая

определители

в

формуле

*П

 

(3.22) и имея в виду,

что

di

ю

 

ф

Ф

to

г\

 

 

-

 

 

a+na^ V

^

fl»yfl+y =

0'

 

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

К \ р 2+ К \ р

& п

(3.26)

 

 

■w *n{p)- Т у , р { ; г ^ + ‘х 9Т л + \) '

371

где

* 5

__

 

 

 

 

“ *ф

 

 

 

 

*

4

2 ~

a *

'

 

 

ф

со

,

ф

со \

a “ W, - aW

 

+ al : s )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

о

со

 

 

 

ь N

— ■

 

 

— а ф« а со

у

 

 

 

 

 

 

 

уп

 

у

 

 

 

4 , 1

<

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а + а

6

 

5

Ф

 

 

 

» 0

 

 

 

* — а , кДш

 

 

 

 

 

N

“ у

 

фл

 

 

(3.27)

«ф

— "

 

 

Ф

0)

1

ф

ш \

*

гп

 

а-1 ( <

 

v

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т И' =

лЪ „<■>.

 

 

 

 

 

 

 

 

 

ф

и)

 

ф

о>

 

 

 

 

 

aZ а.у + а,»а

у

 

 

 

 

 

 

 

У

 

 

“ у

 

 

 

 

 

— а*»

 

а

у

 

 

 

'*{5"

2]Л - а£ а“у + at«u:y

тл у

Так как А?«, и А°« по абсолютной величине значительно меньше,

•Л“ VII

чем AN, то можно пользоваться приближенным значением UPN(p),

определяемым по формуле

К

н(р)

 

п

где

А N:

 

 

А к

 

 

 

 

4

 

 

(3.28)

7'v^ ( 7 ' ^ 3H - 2 ^ V + l ) ’

 

 

a

4

 

 

 

ф

%

I

ф

to

(3.29.)

a>

a ■ a . у 4 - a N а у

 

“ у

Ф

1

Ф „

“ у

 

При малых углах наклона траектории у ракет с крестооб­ разно расположенными крыльями и баллистических ракет

Звено, имеющее передаточную функцию вида (3-28), обладает логарифмической амплитудной характеристикой, изображенной

на фиг. 3.5. Частота ш = -!г- есть собственная частота летатель-

ного аппарата в боковом колебательном движении.

378

б) Пе ’р е д а т о ч н а я ф у н к ц и я п у т е в о г о у г л а по б о ­ к о в о м у в е т р у W ^ (р). После раскрытия определителей в

формуле (3.22) передаточная функция W ^ ( p ) приводится к виду

 

w

\ Р ~ + k l ' \ P

 

где

 

 

 

 

 

< <

 

 

Ч 2~

 

 

Ф„1'

 

 

 

 

ф

4

 

 

«I йсоу V

 

k w(-=

 

По

аналогии с

преды-

щ,))\

дущей

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

функ­

 

цией можо записать

 

 

* 7' ( Р ) «

УП

Ч

Т у . р [ Т № + Х 9Т,р+1)'

(3.32)

(3.30)

(3.31)

но это приближение яв­

 

 

ляется

менее

точным

по

Фиг. 3.5

 

сравнению с

формулой

 

(3.28).

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

в)

лу о т к л о н е н и я р у л я

направления W^h(p). Из формулы (3.24)

находим

 

5

пК?1.рг+~ лКн

 

 

 

(3.33)

 

 

W*H{p) =

TV,P (ТрУ +

2 ^ + Т )

 

 

 

Соседние файлы в папке книги из ГПНТБ