Файл: Пугачев, В. С. Основы статистической теории автоматических систем.pdf

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

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

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

Добавлен: 16.10.2024

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

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

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

0,?

0,6

0,6

0,8

Т,с

Рис. 3.10. Графическое решение уравнений

сечения этих кривых дает искомые коэффициенты.

Система стабилизации лета­ тельного аппарата по крену рас­ сматривалась как линейная. В действительности угол откло­ нения элеронов ограничен. По­ этому необходимо проверить вы­ полненный выше расчет по вели­ чине отклонения элеронов. Ча­ стотная характеристика системы от входного возмущения к углу отклонения элеронов в соответст­ вии с рис. (3.9) имеет вид

Ф (/со) = ___________________ к 4- шт___________________

(3.137)

Т (/со)3 + (1 + а х х Г )

(,W)2 + (аЛЛ- + ядэт) 1(0 + Яд-э*

 

Дисперсия угла отклонения элеронов

 

D6 —

Дсоб,

(3.138)

где эффективная полоса пропускания системы от входного возмуще­

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

элеронов

 

 

 

Acoj = J

 

 

 

 

| к +

Tfсо |2 rfco

(3.139)

| T (ico)3 +

(!-{- axxT) (ico)2 + (axx 4- ахэт) ко 4- ax3k |2

Вычисляя этот

интеграл,

получаем

 

 

 

Дсо. = л —

к (1 +

axxi) + ax3i-

(3.140)

 

°

а..

,[(14

аххТ) (ахх 4- ахэт)—ахзкТ] '

 

Рассмотрим

числовой

пример

при

следующих данных:

ахх —

= 4 с -1; ахз =

40 с-2; Т — 0,05 с;

S.v =

0,2 с-3. По формуле (3.134)

строим зависимость коэффициента усиления от коэффициента демп­ фирования (кривая 1 на рис. 3.10). На этом же графике строим зави­ симость коэффициента усиления от коэффициента демпфирования, вычисляемую по формуле (3.136) при различных значениях £0 — = 0,5, | = 0,7 (кривые 2 и 3 на рис. 3.10). Точки пересечения кривых дают оптимальные значения параметров. Результаты расчетов све­

дены

в таблицу.

 

 

 

 

 

 

£о

*0

То, с

с3

Дсо6, с3

Dy , рад2

D6. рад2

аг

*6

град

град

0,5

3,0

0,17

5,64-Ю-3

6,04Х

1,12- 10-3

1,21 10~2

1,92

6,31

 

 

 

 

X 10'2

 

 

 

 

0,7

5,0

0,40

1,91 -Ю" 3

6,95Х

3,80-10~4

1,39-10--

1,12

6,76

 

 

 

 

XI О’2

 

 

 

 

100


3 .7 . С истем а са м о н а в ед ен и я ракеты

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

11 = ^ —»щ£ц;

 

t = Aa {a — аТ);

(3.141)

а -f- С^а -j- Сасс = С0 С^б;

 

б= — £це — X (t)) -f- koj -}- /ез©г,.

Вэтих уравнениях приняты следующие обозначения: г) — от­ клонение центра массы ракеты Р по нормали от вектора теоретиче­

ской дальности DT\ = v cos (ет — 0T) — проекция вектора ско­ рости ракеты на вектор теоретической дальности; v — модуль век­ тора скорости ракеты; ет — угол ориентации вектора теоретической

дальности;

0Т— угол наклона вектора скорости при теоретическом

движении;

ц1ц = vч cos (ет — 0ЦТ) — проекция вектора

скорости

цели на вектор теоретической дальности;

■— величина

скорости

цели; 0ЦТ—■угол наклона вектора скорости

цели в теоретическом

движении; £, — соответственно вариации углов наклона векторов скорости ракеты и цели; а — угол атаки; а т — угол атаки в теорети­ ческом движении; С0— величина, зависящая от ускорения силы тя­

жести и ее производной; б — угол отклонения руля; е — угловая скорость вектора дальности D\ X (t) — помеха; j — составляющая нормального ускорения, измеряемая акселерометром; со21 — угловая скорость корпуса ракеты относительно центра массы.

На рис. 3.11 показана схема, иллюстрирующая кинематические соотношения. Ракета наводится по методу параллельного сближе-

Р

Рис. 3.11. Кинематика самонаведения ракеты

101


ния, поэтому угловая скорость вектора дальности е является основ­ ным сигналом ошибки в законе управления [последнее уравнение в системе (3.141)]. Закон управления записан в идеализированной форме без учета инерционности измерителей, усилительных устройств

ирулевой машины.

Вуравнениях (3.141) величины Аа, Са, Са, Сб являются пара­

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

чины k lt

/га, к3 являются параметрами системы управления.

Первое

уравнение в

системе

(3.141) описывает

кинематику движения

центра

массы

ракеты,

второе— динамику

центра массы,

третье

уравнение характеризует кинематику и динамику движения корпуса ракеты относительно вектора скорости (по углу атаки), последнее выражение описывает закон управления.

Закон управления в системе уравнений (3.141) записан в полных выражениях. Запишем его в вариациях. Измеряемое акселерометром

ускорение

 

 

j =

vAaa.

(3.142)

Угловая скорость корпуса ракеты есть сумма угловых скоростей

угла атаки и угла наклона вектора скорости:

 

coZi = a -j- 0.

(3.143)

В свою очередь, величина

 

 

0 = Аар. +

,

(3.144)

где gy — составляющая ускорения силы тяжести на ось у.

 

Следовательно, угловая скорость ракеты

 

= ос -|- Ааа -j——.

(3.145)

Угловую скорость вектора дальности представим в виде разности теоретического значения и вариации:

е = ет — бе.

(3.146)

Подставляя соотношения (3.142), (3.145), (3.146) в закон управ­ ления, запишем его в следующем виде:

б == —■kx(ет — бе — X (t)) -j- Аа (k2v -)- к3) -}- k3a -f- k3 • (3.147)

С целью получения аналитических выражений для вероятностных характеристик промаха примем два допущения: вторая производная

угла атаки а = 0; скорость сближения

| = vr = const.

(3.148)

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

D = 0; D (t) = D0— vrt\ D0 = vrT,

(3.149)

102


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

D 0— начальная дальность стрельбы.

С учетом сделанных допущений систему уравнений (3.141) можно

представить в следующем виде:

 

 

£ =

(сс ССТ),

( (3.150)

та -)- а — — /г0 (бе + X)

k0&T-j- /г0

ф ) >j

где т — эквивалентная постоянная времени; /г„ — коэффициент уси­ ления, характеризующий динамику движения ракеты относительно центра массы:

Cd +

Сц 4-

{к-р +

/;з) ’

(3.151)

k0= ______ Сб^1______

 

Со +

С(,Аа (k2v

/г3)

 

Промах самонаводящейся ракеты описывается формулой

^ = 1/т + 'Ч + 'Пд*.

(3.152)

где г/т — промах в теоретическом движении; At

— время полета

до встречи с целью

после выключения

координатора на

расстоянии

DB от цели.

 

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

= D5e (см. рис. 3.11), представим формулу

(3.152) в следующем виде:

 

У =

ут-j- D5e—губе At -f- D6e At.

(3.153)

При постоянной скорости сближения время полета после выклю­ чения системы управления

At = D/vr.

(3.154)

Подставляя At в формулу (3.153), получаем

У — Ут уут бе.

(3.155)

Определим отсюда вариацию угловой'скорости вектора дальности, подставив которую в третье уравнение системы (3.150), получим

та -j- а = —- /е0 ( ^ У + * ) +

+ /?08т+^о1

Со

(3.156)

Cfi/ei

крг

103


Из второго уравнения системы (3.150) определяем угол атаки и подставляем его в уравнение (3.156), тогда

Я + i = - М а

г + х ) + м Я Н-

<зл57)

Представим значение промаха в следующем виде:

У = Ут

(3.158)

 

vr

и продифференцируем уравнение (3.158) по времени:

У = Уx + +

+

(3.159)

Дифференцируем первое уравнение системы (3.150) в предполо­ жении постоянства скоростей:

“Л= °iS — ищ£ц-

(3.160)

Подставляя это соотношение в формулу (3.159), получаем

 

Y =

 

 

vr

 

(3.161)

 

 

 

 

 

 

Отсюда находим

 

 

 

 

 

 

 

£ =

 

 

и1ц^ц~) >

(3.162)

у- _

vr I

 

Dvm V

и1цСц D +

 

S ~ v^D- I У — Ут —

- Vr

£ц +

 

 

 

У Ут

 

D

 

(3.163)

 

Ь

 

Vr ^ 1ц ?ц \ .

Подставляя Си

С в уравнение (3.157)

и преобразовывая,

полу­

чаем уравнение для текущего промаха самонаводящейся ракеты:

тY + ( l + ^ f ) y + k ^ Y =

= ~ (Т—5Г" Я ^Vr^ ) Я /т (^)> (3.164)

где приняты следующие обозначения: k — обобщенный коэффициент

усиления,

 

 

 

(3.165)

k = k 0Aav J v r\

 

У т— проекция нормального

ускорения

цели

на вектор

теоретиче­

ской дальности,

 

 

 

 

jlu, --- /ц COS (бт

9цт)

(®т

®цт) >

(3.166)

104