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

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

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

Система (2.19) является линейной с постоянными коэффи­ циентами. Если исходное движение иеустамовившееся, то коэф­ фициенты системы уравнений (2.19) являются функциями вре­ мени. Такие системы в общем виде не решаются и требуют специальных методов исследования.

Коэффициенты при неизвестных в системе (2.19) определяются по формулам, приведенным ниже (индекс „нулик" опущен, так как все величины, входящие в формулы, относятся к расчетному движению):

1.

 

P V-

Cr + Сх _

2м pl/S

Производная Р 1' определяется по

характеристикам двигателя,

производная сх — по кривой cx — f { М)

для расчетного угла,

атаки

а.

 

pl/з

 

 

2.

айу =

(Q* — О cos 0) = ■

5 — G cos О

->а I

 

' Производная с“ определяется по кривой cx=f(a) для расчет ного числа М при расчетном угле атаки а.

3.

О*

О*

О1/3

 

 

 

 

 

 

 

4.

а

=

Рп _ су . /

м

.

Производная рГ|

определяется

по кривой

р=/(т|) пли по таб

лице

стандартной атмосферы.

 

 

 

5.a ivp = j j P * p .

Производная Р ор определяется по характеристикам двигателя

р= / ы -

6.4 K = f ( - Q ‘*) = - g - 4 * f- r 5 .

Производная

определяется по полярам, построенным для

различных значений ов при его расчетном значении.

/ .

( Qv cos о +

~

sin 0 — Y

sin 0j =

 

С х + С х - у \ C O S

0 +

(с* — С у )

sin pl/s.

340

8 .

a iT7' = -Q- ^Q!/sin 0 — Qa cos 6 + Y ~ cos О

= ■§- ^ + cj?-yjsinfl + \ { — с’ + су) cos 6 pVS.

9

g y - g Y v

'

b ~ G V ]

2^ cos0

-

g

(r

;

CMM )n S

2g-cos 9

 

l/2

Vfcoc j

~

G

{

y ^

y 2 J p,b

К О С . I

 

Производная

 

определяется

по

кривой

cy= /(M ) для

рас­

четного’ угла

атаки.

 

 

 

 

 

 

 

10.

о

S'

 

■Р —

Ум-G

1 -

V2

 

sin 0

 

 

1/2 .

 

 

0

GV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

К О С .

1

 

 

 

. _ g _

 

P

- c - ^

s + a l i —

^

 

sin 0

 

 

GV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Производная

с* обычно задается как функция числа М в числе

аэродинамических характеристик летательного аппарата.

 

11.

 

g

 

 

m

(

p + c‘Y }r

s

)-

 

 

 

GV l-p + ^ ~

 

12

 

G V

G

y'

V

5

 

 

 

 

 

2

 

 

 

 

.

 

g

у7' = —

c oT|—

 

 

 

 

13

V . =

- £ - r . = X

c *.pJ£s

 

 

 

 

 

о

GV

G

у

2

 

 

 

 

 

Производная

с 8

определяется

по

полярам, построенным для

различных значений 38 при расчетном

числе

М..

 

14.

 

 

 

 

 

 

 

 

 

]_

 

a„Vi = -M ri—Р-тт sin 0—Q -^г sin 0— Y v cos 0 — У '- ^ э т б

 

 

G V \

V

 

 

 

 

 

V

 

 

Я sin 0 4--^- — lcx + c™ -^\cosBу (c; + cv)sin0

p5.

15.f p cos 6 + Q-^r cos 0— Yv sin04-Ya-j^ cosO^ =

G Vs

P cos 0 + -|r I Cy + c f

sin 0 + Y (c* + cv) cos ' ? s .

 

 

G

 

1б. а : =

^ + мГР) = - 1

 

341

Для ракет с крестообразно расположенными крыльями и

баллистических ракет

 

 

 

 

v

1 /

м■М\

,, с т

 

= -J-1

тг + mz

I p V S / .,

где L—полная длина ракеты.

 

то тг 0. Производная

Если расчетный полет прямолинейный,

т? определяется по кривой

 

 

 

тго

- т “а0 +

ml а +

rnze8e= / (М),

где а0— угол атаки, соответствующий су — 0, а и 8в — угол атаки

и угол отклонения руля высоты в расчетном полете, тг0—коэф­ фициент продольного момента при cv— 0.

Производная

М^р определяется

по формуле

*

 

 

 

"

г Р

= P vl

 

 

 

r

‘‘Р

где

плечо

силы тяги

относительно центра тяжести ракеты.

17.

1

i

„ к \

 

1

*pl7

о # ч

 

 

 

 

 

а t-

< = x

t

- ^ ) =

-

x-J- in-

- -

Sba~.

Для ракет с крестообразно расположенными крыльями

о1 ; р V е

а„,_ — -J- П12 -g- о L .

Производная ml обычно задается как функция числа М. Для ракет, выполненных по схеме „утка“, „летающее крыло" и для

баллистических ракет /я" = 0 .

18. а \ = Ф ( - Ml - M‘J = - х U Sf>„ + К

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

листических ракет

1 7

a?V*

,

’ о

а =

---- 7-1 mz ‘—г- SL 4-

М ,,

 

Л

 

 

‘ ZP

 

 

 

 

Производная ml обычно задается как-

функция числа М. У ра­

кет, снабженных РДТТ или ЖРД,

М1Р— 0. У ракет, снабженных

ВРД,

М* —

\/у

.

 

 

т г Р

о

' Х р ’

 

342

где Gcei.— секундный

расход воздуха через воздухозаборник,

х р — расстояние

от центра тяжести ракеты до входа в воз­

духозаборник.

19, а к =

Тг

 

m

+

'

[

(

<

*

+ SbЛJ М Щ •

Для ракет с крестообразно

 

расположенными крыльями и бал­

листических

ракет

 

 

 

 

 

 

 

 

 

 

 

а 2 —

1.

r{m;2 +

rrf2 ) p- ^ S L ^ M l

 

 

 

 

 

J,

 

 

 

 

 

 

 

 

Производная

m“z

обычно

задается

как

функция числа М,

 

 

 

СО

 

 

 

 

по формуле

 

 

Производная Mz* определяется

 

 

 

 

 

 

 

m

_ _ @сек Л2

 

 

 

 

 

 

 

z P —

 

о-

 

г.т,

 

 

 

 

 

 

 

 

 

 

to

 

 

 

 

где 1гж — длина

газового

тракта (расстояние

от входа в воздухо­

заборник до

выхода из

реактивного

сопла).

 

 

20.

 

 

 

 

 

 

 

1

PV2

 

 

<

=■ i

К

+

м '- ) =

 

i ( mtf-

f Sb“ + M 4

-

 

 

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

*

1

( u p V * .

 

°ш. ~

j

( mz

2 ^ 1

 

21. < = ~ (л& + М^р) =

-j- ( m / ‘^

+ Л ^,) .

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

а2_— j - ( тУ ^2~ SL + Щ‘zP

Если расчетным движением является прямолинейный полет, то тг = 0. Производная М]р определяется по формуле

M T:p = P 4 L

2 2 .

343

Производная УИ°£ определяется по формуле

 

М ] р = р 1р1р .

23.

- J - Л/:« = у от_> ~ ~ Sba .

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

листических

ракет

 

 

56 Р

 

 

 

 

 

 

 

%

1

С I

 

 

 

 

 

 

а ш®= у

i n /

- S L .

 

 

Производная

от°я обычно

задается как

функция числа М.

24.

1

/

1 .

 

1

к

р

Sb* sin О

= 4 - (

— УИ; 4т sin 0 ] =

---- т- от/

2

-

Jt

\

V’“‘“

I ~

J,

""

 

Для ракет с крестообразно расположенными крыльями <

 

 

 

aj*

---- j- т~

SL2 sin 0 .

 

 

 

 

7

 

 

J*

2

 

 

 

Для ракет, выполненных по схеме

„утка",

„летающее крыло"

и баллистических

ракет,

 

 

=

0.

 

 

 

 

25.

а*} = -

у

(/И* + М;Р) cos 0 +

( М% +

М:р) - 1 sin 6

 

"

Jг

 

 

I/

 

л |

0

^

* Л

 

 

 

 

 

 

 

 

 

 

 

= — аш cos б +

alo -=7- sin 0.

 

 

 

 

 

 

г

 

 

j у

 

 

26.

а*ч = у

Л4/ у

cos 6 =

у

 

 

S b J cos 0 .

Для

ракет с крестообразно расположенными крыльями

 

 

 

a ^

= \ - m ~ 4 r SL1cos0.

 

 

 

 

“г

 

Л

2

 

 

 

 

Для ракет, выполненных по схеме „утка",

„летающее крыло"

и баллистических

ракет,

а ш^ =

0.

 

 

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

27. a® V = -L

 

(ЛГ/ + Мгр) sin 0 +

(/И- +

1

 

Л»“Р) - 4 cos

 

Л

 

 

 

 

 

 

 

 

 

 

 

 

 

: — a

К .

л

0

1

Л

 

 

 

 

“г

sin 0 + а

шл

-ту cos 0 .

 

 

 

 

 

 

 

 

V

'

 

344

28.а% = sin О.

29.a ° =l / c os O.

30.J v = A P e - A Q e .

31../ 0 = AY ,.

32.L z = m zll+AMzPa.

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

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

рА Vn =

а I AV„+ a°v ДО +

а% ДО + а \ А-п

 

 

 

 

pAO — ад A Vn-j- адДО -г ад ДО + ад Дт)

 

 

 

 

ршг =

а-1

Д V„+ (а\ р + аМ ДО + а > г + а®

ДО +

Дт,

(2.20)

 

у?Д0=

2

\

Z

2 1

2

Z

2

 

 

сог

 

 

 

 

 

 

 

 

pAri — a*AVn + о“ ДО

 

 

 

 

 

 

Общее

решение системы (2.20) имеет вид:

 

 

 

Д Vn = A 1ePli + А 2еР*' + Л 3/

3' + л /

'

+ л / 5'

 

 

 

 

ДО: : ВУ 1' + в У ~ + В3еР^ + В .У ^ + Вье ъ*

 

(2.21)

 

 

а», = с / 1*+ с , / 2' + СУ 3' + с У ‘+ с - У 1

 

 

 

д о = d/ 1' + в У 1+ в /*1+ вУ ' +вУ 1

 

 

 

At]= еУ 1 +еУ 3+еУ з>+еУ 4 + Еде ъ

 

 

где A l ,

Bit

С,, В п

Et — постоянные коэффициенты,

pt — корнй

характеристического уравнения системы. Постоянные коэффи­ циенты и корни могут быть как вещественными, так и комплекс­ ными.

343

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

р~ а \

-<

 

.

 

^

1

<3

3

 

м

 

0

V

4

 

0

0

- 4

 

 

a v

 

 

р -

0

0

- 4

- 4

 

ао

 

- ai2P -

■0

а

0

= 0. (2.22)

аШ,

р — яш‘

— Яш

0

6

1

Р

0

 

 

0

0

Р

 

 

Л

 

Раскрывая определитель, получаем характеристическое урав­ нение в окончательной форме:

 

 

р 5 +

atp* +

а2р3 +

а3р2+

я 4д +

я5 = 0,

 

(2.23)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

аЛ

 

 

 

 

 

 

 

 

о2

 

о

 

“ч /

 

 

v

 

 

о \

 

О б ,

 

 

к о

 

е v

 

 

 

С2ш^

^

 

ciy

 

аь)

 

 

о

ш ,

 

о.у

 

 

 

 

 

 

 

 

 

 

 

 

a7)La К

а7]*

О

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

К

7)

 

о

7)

 

 

 

 

 

 

 

а 3—

 

о /

 

V'

 

 

е\

 

 

о о

 

 

 

 

 

о

о к \

 

 

 

 

а к

 

а ь)

 

 

аа а и>^

Ли\ . \ д кйо

a v<xQ)

 

 

 

 

— ^

i

/

 

о к

 

 

 

v

o\

 

o

 

k

,

 

 

 

 

 

 

 

2\ a vaо

 

- a j / fl0j -

a v, a ^ +

 

 

 

Ti Г

К в

 

 

И /

 

в

 

 

Ш V I

,

1 Г

 

О

^

V .

О /

V ,

ш м

+ а Ц - . а в_ап +

ат>(аь +

 

яш*)] + яе' [ - а

 

+

ат>( я к +

яш*)] ,

д 4 —

а шг у 1 уО'ъ

 

а и а 0 J + а шг ( / Ч а е

 

 

а к а е )

 

 

+ 4 г ( 4 4 - 4 4 ) + 4 [ < ( < - а °о4 * +

 

 

+

я^ а !г) +

 

я° (я ^ а “*)] +

 

aj [ят, (я°ка “* —

 

— 4

<

)

+

а\ ( ° 1 , -

4

а “*)]

'*

Ч ( ~

a v an ~

a в 4

) ’

 

 

 

 

 

 

в

 

 

ft

 

,

ft

А

\ .

 

 

 

 

 

 

 

 

 

Ч

< (

-

- ^

 

\

+ Щ < ) + ^ ( 4 4 - 4

4 г ) ] +

 

 

ft

 

 

» 0 \

 

 

 

 

 

V » ,

9 „ Г \ 1 .

 

Г V / в

я

 

0 /

 

 

 

К

( a v

шг

я „

 

я ,

) .

 

 

 

 

 

 

 

 

 

 

 

в

шг /

+ а ч ( — а 1/а ш _ + a v a a . ) ] +

 

 

7]

Г V I »

 

0

 

0 » \ ,

0 ( „ V ft

 

ft К\1

 

 

Ч

Ы

 

d y d f j

 

 

d y d 0 ) + 4

 

( а ^ е

a v ( l i JJ •

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

346

Когда расчетное движение устойчиво, в общем случае корни характеристического уравнения' (2.23), если нумеровать их в по­ рядке убывания модулей, получаются такими:

P i , Рз

большие по модулю комплексные сопряженные корни;

Рз — малый по абсолютной величине вещественный корень;

Pi, Р ь ~ наименьшие по модулю комплексные сопряженные корни.

Или

такими:

Р\ >Р'1 — большие по модулю комплексные сопряженные корни;

Рз> Pi ~ малые по модулю

комплексные сопряженные корни;

р ь — наименьший по

абсолютной

величине

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

корень.

 

 

 

Такое распределение корней соответствует физическому харак­

теру переходного процесса в продольном

движении.

Переходный

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

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

На фиг. 2.2 изображена осциллограмма переходного процесса малой ракеты с дозвуковой скоростью полета, вызванного скачко­ образным отклонением руля высоты. Собственное демпфирование у ракеты мало, поэтому время затухания корошопериодичеокогр движения довольно велико. Частота короткопериодических колеба­ ний около одного герца. Частота длиннонерио-дических колебаний около 0,02 гц. .

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

347

■ния скорости и высоты .равны нулю. При этом для получения реше­ ния, соответствующего коротко1периоднческо'Му движению, доета-

Фнг. 2.2

точно в системе (2.20) отбросить первое н последнее уравнения, а в остальных уравнениях положить

ДV „ = Дт] = о.

>

§ 2.3. ПЕРЕДАТОЧНЫЕ ФУНКЦИИ В ПРОДОЛЬНОМ ДВИЖЕНИИ

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

ского регулирования. Динамические

свойства летательного

аппа­

рата при этом задаются с помощью

передаточных функций.

 

Определим передаточные функции

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

в продольном движении, пользуясь

тем,

что продольное

возму­

щенное движение можно разделить

на две фазы. В первой фазе

переходных процессов в продольном движении наиболее важным является выяснение характера изменения таких параметров, как угол наклона траектории 0„, определяющий направление вектора путевой скорости V„, угол тангажа 0, определяющий направле­ ние связанной оси OxX i, и нормальная перегрузка пу .

При определении передаточных функций для этих параметров, как это было указано выше, необходимо учесть, что отклонение руля высоты или действие вертикального порыва ветра обусло­ вливают существенное изменение углов 0Л и 4, а также нормаль­ ной перегрузки, причем изменение этих параметров относительно больше, чем изменение скорости п высоты полета. Поэтому для вывода интересующих нас передаточных функций воспользуемся уравнениями короткопернодичеекого движения. Необходимо иметь в виду, что полученные таким способом передаточные функции отражают с необходимой точностью динамические свойства ле­ тательных аппаратов только в первые секунды переходных про­ цессов. Для получения точных результатов следует, пользуясь общим методом, найти передаточные функции из уравнений (2.19). Однако передаточные функции, полученные из усеченной систе-

348

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

С. учетом введенных входных воздействий уравнения коротко­

периодического

движения,

т. е.

когда

Д-/| =

0 и AVn =

0,

имеют

вид:

 

 

 

 

 

 

 

 

 

 

рЛ6„ =

а,0Д0„ + <4 ДО +

a;* ДО* +

UP,

 

 

рч>г — аш_рАОп+ а<0 Д0„ -J-

 

юг -|-

ДО -Ь

 

 

 

+ а \

А8« +

 

Р Wr>+

 

 

дда =

шг

 

 

 

 

 

 

 

 

 

Подставляя третье

уравнение

во второе,

перенося все

члены,

кроме входных

воздействий,

в левые

части

уравнений,

находим

-

а';) Д0л - аI ДО =

а*« До, + a ^ W n

 

 

[a!L_p + а1г) ДО,, + [р '—а ^ р

Д0=

 

(2.24)

 

— a il До, +

[ а ^ р +

а«гт‘) WTi

 

 

Главный определитель системы (2.24) можно представить в сле­

дующем виде:

 

 

 

 

 

 

 

 

 

 

Д (р) =

P -

а

е

 

а ь

 

 

(2.25)

О

 

О

О

“г

о

 

 

аш р — аш р- аш‘ р —аш

 

 

 

Z

 

 

2

 

2

2

 

 

Тогда, решая

систему (2.24) относительно Д0Я, получим

 

а!в AZe - f -

a ^ v. W r,

 

 

- а *

 

 

aj> До„ + (асГ1р + а„Гч) WT, р- a2z р —а!г

 

 

ДО.

 

 

 

 

Ч р )

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ао°

 

 

 

 

 

 

 

 

 

 

а */г — а. —аш

 

 

 

Д0.=

(D~

 

(.0*

 

г

Дов

 

 

 

 

Ч / >

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

\\7

 

\Г>'

«>

«>-

 

 

 

 

+

а„,^ /? + а„,’i/?"— а„,;д-

 

 

(2.26)

 

 

 

 

 

 

 

 

 

349

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