Например, Бобцов

ФОРМАЛИЗАЦИЯ ВЕКТОРА НАБЛЮДЕНИЙ ИЗМЕРИТЕЛЬНОГО КОМПЛЕКСА БЕСПИЛОТНЫХ ЛЕТАТЕЛЬНЫХ АППАРАТОВ

Формализация вектора наблюдений измерительного комплекса беспилотных ЛА

23

УДК 681.5.01

Л. В. АНДРЕЕВ, С. В. БОГОСЛОВСКИЙ, Б. В. ВИДИН, И. О. ЖАРИНОВ, О. О. ЖАРИНОВ, П. П. ПАРАМОНОВ, Ю. И. САБО
ФОРМАЛИЗАЦИЯ ВЕКТОРА НАБЛЮДЕНИЙ ИЗМЕРИТЕЛЬНОГО КОМПЛЕКСА
БЕСПИЛОТНЫХ ЛЕТАТЕЛЬНЫХ АППАРАТОВ

Проанализирована структура измерительного комплекса беспилотного летательного аппарата. Формализован вектор наблюдений.

Ключевые слова: измерительный комплекс авионики, вектор наблюдений.

Введение. В систему управления беспилотного летательного аппарата (ЛА) помимо

сигналов, характеризующих вектор состояния объекта управления, поступают сигналы (воз-

действия), параметры которых система не может воспринимать. Такие сигналы являются

внешними по отношению к системе управления. Внешние воздействия условно разделяют на

полезные, обеспечивающие нормальную работу системы управления, и вредные, снижающие

эффективность работы системы управления, т.е. возмущающие, или помехи.

Рассмотрим принципы управления по обратной связи и управления по возмущению,

учитывающие измерение внешних воздействий и подачу измеренного сигнала через усили-

тель на вход устройства системы управления (УСУ).

Первый шаг математической постановки задачи в общем случае подразумевает конст-

руирование компьютерной модели, т.е. качественное описание исследуемого процесса с по-

мощью тех или иных операторов. Для решения этой задачи исходную блочно-операторную

схему преобразуют в структуру, имеющую один выход.

Преобразование передаточных функций системы управления ЛА. Общий подход к

получению передаточных функций (ПФ), в том числе параметрически заданных, заключается

в нахождении изображения выходной координаты методом преобразования Лапласа в зави-

симости от возмущения, интересующего исследователя.

Правила преобразования структурных схем стационарных систем, в том числе с пере-

крестными связями, подробно изложены в литературе [1]. В качестве примера рассмотрим

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

рирования (рис. 1). На рис. 1 приняты следующие обозначения: x1, x2 , ..., x7 — координаты

вектора состояния системы — выходные координаты отдельных звеньев, стоящих в прямой

цепи (координата x4 на схеме не показана, так как предполагается, что передаточная функ-

ция WC ( p) имеет второй порядок и координата x4 скрыта внутри соответствующего блока

структурной схемы); K — постоянные коэффициенты (индексы относятся к отдельным бло-

кам структурной схемы);

1 τ

— ПФ нестационарного кинематического звена,

τ = t −t0

— ос-

тавшееся время полета ЛА до конечной точки траектории (КТТ), взятое со знаком минус,

t0

=

D0 U

— полное время полета до КТТ, D0

— начальная дальность до КТТ; U — скорость

изменения дальности (скорость сближения с КТТ);

−WC

(

p)=−

Ka p

TC2

T1C p +1 p2 + 2ξCTC

p +1

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 11

24 Л. В. Андреев, С. В. Богословский, Б. В. Видин и др.

— ПФ ЛА в вертикальной плоскости; ξC , TC — коэффициент относительного затухания и посто-

янная времени самолета в горизонтальной плоскости соответственно;

Wkz

(

p)

=

T1C

1 p

+1

— ПФ инерционного кинематического звена, преобразующего угол тангажа в угол наклона

вектора

скорости;

−1 U



стационарный

коэффициент

усиления

нестационарного

кинематиче-

ского

звена;

1 p

— ПФ интегрирующего кинематического звена; V

— модуль вектора скорости

ЛА; Wосϑ ( p) = K1 + K2 p — ПФ стабилизирующей обратной связи по углу тангажа; K1, K2 —

коэффициенты обратной связи по углу и по угловой скорости тангажа соответственно;

WУСУ

(

p)

=

KУСУ p

— ПФ устройства формирования сигналов управления;



K

уWy

(

p

)

=

Ka p

p +

+ Ka1 Ka1



K

y

+



p

— ПФ звена упреждения; K y — коэффициент усиления.

WGRP ( p)

Ka p

x7 KУСУ p

x6

1 x5 p

−KCWC ( p)

x3

Wkz ( p)

x2 V 1 p

x1

WУСУ ( p) –1

Wосϑ ( p) (K1 + K2 p)

1

3

Ka p

p + Ka1 + Ka1



K

y

+



p

2

−K yWy ( p)

−1 ξ

1



Рис. 1

Для преобразования схемы рис. 1 к схеме с одним входом и одним выходом необходи-

мо последовательно найти передаточные функции W1( p) , W2 ( p) , W3 ( p) контуров, обозна-

ченных на рис. 1 цифрами (1—3). При этом передаточная функция стационарной части сис-

темы (от выхода ξ нестационарного звена

1 τ

до выхода по координате

x1) будет (передаточ-

ная функция W1( p) включается в состав W2 ( p) )

Φ( p) =

x1( p) ξ( p)

=



1 U

W3

(

p)W2

(

p)Wkz

(

p)

VС p

.

(1)

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 11

Формализация вектора наблюдений измерительного комплекса беспилотных ЛА

25

При Ka1 =1 получим

W1 (

p)

=



(T1С

p +1)( K1

−KС + K2

(T1С p)+

p +1) p(TС2

p2

+

2ξС TС

p +1)

=

P1( p) Q1( p)

;

P1( p) = −KС (T1С p +1) ;

Q1( p) =TС2 p3 +(2ξСTС + KСT1С K2 ) p2 +(1+ KСT1С K1 + KС K2 ) p + KС K1;

W2

(

p)

=

P2 Q2

( (

p) p)

;

P2 ( p) = −KУСУ KС (T1С

p+1) ;

Q2 ( p) =TС2 p4 +(2ξСTС + KСT1С K2 ) p3 +(1+ KСT1С K1 + KС K2 + KωK y KСT1С ) p2 + + ⎣⎡KС K1 +(1− K y )KУСУ KСT1С + KωKУСУ KС ⎦⎤ p+(1− K y )KУСУ KС ;

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

W3

(

p)

=

P3 ( p) Q3 ( p)

;

P3 ( p) = Ka ;

Q3 ( p) =

p+Ka .

Φ( p) =

a1

p6

+ a2

p5

+

a3

b1 p4 +

a4

p3

+

a5

p2

+

a6

p

,

b1

=

VС U

K y KС Ka ;

a1 = TС2 ;

a2

= TС2 Ka + 2ξcTС + KСT1С K2 ;

(2)

a3 = Ka (2ξСTС + KСT1С K2 )+1+ KСT1С K1 + KС K2 + KωKУСУ KСT1С ; a4 = (1+ KСT1С K1 + KС K2 + KωKУСУ KСT1С )Ka + KС K1 +(1− K y )KУСУ KСT1С + Kω KУСУ KС ;

a5 = ⎡⎣KС K1 + (1− K y )KУСУ KСT1С + Kω KУСУ KС ⎦⎤ Ka + (1− K y )KУСУ KС ;

a6 = (1− K y )KУСУ KС Ka .
Приведение системы к нулевым начальным условиям. Формула (2) представляет собой искомую передаточную функцию стационарной части нестационарной системы управления с одним входом и одним выходом. Подобные передаточные функции используются в теории стационарных систем для построения переходных процессов по вещественной частотной характеристике, при исследовании систем с одной нелинейностью методом гармонического баланса и в некоторых других случаях, когда на первый план выступают требования наглядности получаемых результатов. В системах нестационарного управления с особой точкой такое представление передаточных функций необходимо для получения интегрального представления решения по первой координате. Рассмотрим подробнее внешние и начальные возмущения, действующие в астатической системе нестационарного управления.
Под возмущением будем понимать факторы, влияющие на параметры движения динамического объекта и не зависящие от системы управления.
Начальные возмущения, или начальные условия, определяются предшествующим маневрированию участком движения ЛА. Чаще всего маневрированию предшествует участок установившегося движения, например, горизонтальный полет ЛА на заданной высоте или прямолинейное движение в плоскости маневрирования. В этом случае, если уравнения движения ЛА и рулевого тракта записать в вариациях параметров относительно балансировочных значений, начальные условия будут нулевыми.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 11

26 Л. В. Андреев, С. В. Богословский, Б. В. Видин и др.

Однако следует иметь в виду, что в момент перехода системы автоматического управления в режим маневрирования к заданному положению, нули некоторых измерителей и звеньев могут не соответствовать балансировочным параметрам движения. Типичным в этом смысле возмущением для астатической системы автоматического управления (САУ) является возмущение fУСУ0 = ϑGRP0 − K yϑy0 , возникающее в том случае, когда не производится „обну-
ление“ по сигналу ϑGRP привода антенны радиоприемника и сигналу ϑy упреждения.
Учет обнуления заключается в том, что при моделировании полагают fУСУ0 = 0 и вво-
дят возмущение fε = −ε0 , компенсирующее начальное значение угла места. К наиболее характерным возмущениям, действующим в САУ в вертикальной плоско-
сти, относятся: — начальное значение высоты полета H (0) = H0 ;
— установочный угол fУСУ (t) = ϑGRP , вводимый для исключения возможности занижения САУ высоты полета ЛА;
— случайные возмущения, вызванные наличием пассивных и активных помех в антенном тракте, эти возмущения могут быть учтены заданием функции fx1 (t) ;
— движение посадочной платформы (может быть учтено, если положить U =V −Vпx , где
Vпx — проекция скорости движения посадочной платформы); — изменение состояния атмосферы, приводящее к движению воздушных масс; движе-
ние воздушных масс принято делить на крупномасштабное (ветер) и мелкомасштабное (турбулентность), при этом ветер учитывается возмущением fW (t) , непосредственно приводящим к изменению скорости движения центра масс, а турбулентность приводит к изменению углового движения вокруг центра масс и может быть учтена возмущением рулевого тракта fPT (t) [2].
Наиболее характерными возмущениями САУ в горизонтальной плоскости являются: — движение посадочной платформы с постоянной скоростью, движение платформы в этом случае может быть учтено заданием функции fx1 =VпZ ; — случайные возмущения, вызванные наличием пассивных и активных помех в антенном тракте; — начальное значение бокового отклонения Z (0) = Z0 , которое может быть вызвано помехами в момент начала маневрирования; — переменный ветер, влияние скорости которого можно учесть согласно рис. 2.

fW − pWkz ( p) fPT



KC p

WC

( p)

fW

Wkz ( p)

Z V p

Рис. 2

На рис. 2 приняты следующие обозначения: fW — возмущение, обусловленное воздей-

ствием ветра — отношение скорости ветра к скорости ЛА;



KC p

WC ( p) = −

KC p

T1C p +1 TC2 p2 + 2ξCTC

p +1

— передаточная функция ЛА в горизонтальной плоскости, T1C — постоянная времени кине-

матического звена.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 11

Формализация вектора наблюдений измерительного комплекса беспилотных ЛА

27

Все указанные выше детерминированные возмущения могут быть описаны функциями времени, для которых существуют и известны преобразования по Лапласу, а случайные возмущения могут быть описаны известными случайными функциями. В частности, предполагается, что заменой переменных все начальные возмущения могут быть учтены в виде внешних возмущений таким образом, чтобы в новых переменных система имела нулевые начальные условия. Приводимость системы к нулевым начальным условиям необходима для получения передаточных функций отдельных звеньев и последующего их объединения в общую структурную блочно-операторную схему.
Заключение. Анализ структуры САУ позволил путем последовательного преобразования перейти от многоконтурной схемы к одноконтурной. Такой переход позволяет определить параметрически заданные передаточные функции стационарной части нестационарной системы, применяемые при исследовании систем с одной нелинейностью. Классификация внешних возмущений позволила формализовать вектор измерений и привести систему к нулевым начальным условиям для получения передаточных функций отдельных звеньев и последовательного их объединения в общую структурную блочно-операторную схему.

СПИСОК ЛИТЕРАТУРЫ

1. Герасимов А. Н., Григорьева Н. Н., Жаринов И. О., Жаринов О. О., Исаков В. И., Орлов А. П., Кулин А. Н., Шепета А. П. Линейные системы автоматического управления: Учеб. пособие. СПб: ГУАП, 2009. 232 с.

2. Богословский С. В., Богословский В. С. Теория нестационарного управления: Учеб. пособие. СПб: ГУАП, 2005. 380 c.

Леонид Владимирович Андреев Сергей Владимирович Богословский
Борис Викторович Видин Игорь Олегович Жаринов Олег Олегович Жаринов
Павел Павлович Парамонов Юрий Иванович Сабо

Сведения об авторах — аспирант; Санкт-Петербургский государственный университет
аэрокосмического приборостроения, кафедра стабилизации, навигации и управления; E-mail: Lio8300@mail.ru — д-р техн. наук, профессор; Санкт-Петербургский государственный университет аэрокосмического приборостроения, кафедра стабилизации, навигации и управления; E-mail: postmaster@elavt.spb.ru — канд. техн. наук, профессор; ОКБ „Электроавтоматика“ им. П. А. Ефимова, Санкт-Петербург; E-mail: postmaster@elavt.spb.ru — канд. техн. наук, доцент; ОКБ „Электроавтоматика“ им. П. А. Ефимова, Санкт-Петербург; E-mail: igor_rabota@pisem.net — канд. техн. наук; Санкт-Петербургский государственный университет аэрокосмического приборостроения, кафедра моделирования вычислительных и электронных систем; E-mail: zharinov@hotbox.ru — д-р техн. наук, профессор; ОКБ „Электроавтоматика“ им. П. А. Ефимова, Санкт-Петербург; E-mail: postmaster@elavt.spb.ru — д-р техн. наук, профессор; ОКБ „Электроавтоматика“ им. П. А. Ефимова, Санкт-Петербург; E-mail: post-master@elavt.spb.ru

Рекомендована кафедрой систем управления и информатики

Поступила в редакцию 01.07.09 г.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 11