Синтез контура командного теленаведения летательного аппарата в программном комплексе «Моделирование в технических устройствах» (МВТУ 3.6) (96

Исходные данные

Исходными данными для синтеза командной системы теленаведения (КСТН) являются результаты синтеза системы угловой стабилизации летательного аппарата (ЛА), разработанного в рамках курса «Управление движением летательных аппаратов». В настоящих методических указаниях в качестве примера рассматривается ЛА со следующими характеристиками:

T 0,072 — постоянная времени углового движения ЛА, с;

T1 0,572 — постоянная времени поступательного движения ЛА, с;

ξ 0,063 — коэффициент демпфирования колебаний; K 0,608 — коэффициент усиления ЛА, с 1;

Kу 140,1 — статический коэффициент усилителя;

Kрп 0,005 — статический коэффициент усиления рулевого привода, рад В;

Kдг 115 — статический коэффициент усиления демпфирующего гироскопа, (В с) рад;

V 480 — скорость зенитной управляемой ракеты (ЗУР) в рассматриваемой точке, м с;

jn 124,5 — нормальное ускорение ЗУР в рассматриваемой точке, м с2;

r 1858 — расстояние до рассматриваемой точки, м; ϕц 0,5 — угол положения цели, рад.

На начальном этапе синтеза сформируем расчетную структурную схему КСТН при наведении без упреждения в виде, представленном на рис. 1.

Здесь Kh и Kh — коэффициенты двучленного закона управления, величины которых находятся в процессе синтеза КСТН; Kлс — коэффициент усиления линии связи, который в нулевом приближении принимаем равным единице, e τлсs — звено с чистым сдвигом, которое приводит к повороту аплитудно-фазовой частотной характеристики (АФЧХ) системы в разомкнутом состоянии по часовой стрелке, поэтому можно проводить синтез контура наведения без учета влияния этого фактора, но ввести в расчет дополнительный

θ 0 1θ ˜

запас устойчивости по фазе лс , 0; WСС s — передаточная

3

Рис. 1. Схема КСТН

Рис. 2. Схема угловой стабилизации ЛА

функция системы стабилизации (СС) ЛА, расчетная структурная схема которой для нормального режима полета jn jn max представлена на рис. 2.

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

4

Kрп 0,005 — статический коэффициент усиления рулевого привода, рад В;

Tрп 0,014 — постоянная времени рулевого привода, с; ξрп 0,5 — параметр демпфирования рулевого привода;

Kу 140,1 — статический коэффициент усиления входного усилителя;

Tу 0,002 — постоянная времени входного усилителя, с;

Kдг 0,005 — статический коэффициент усиления блока демпфирующего гироскопа, (В с) рад;

Tдг 0,05 — постоянная времени блока демпфирующего гироскопа, с;

ξдг 0,5 — параметр демпфирования блока демпфирующего гироскопа.

Порядок работы при синтезе КСТН

Этап 1. Упрощение математической модели СС

Первым этапом синтеза является работа по упрощению математической модели СС для понижения порядка синтезируемой системы.

Выходом СС, используемым в контуре наведения, является уг-

ловая скорость θ. Этот сигнал является низкочастотным для СС

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

Передаточная функция системы стабилизации при этом будет иметь вид

5

Обозначим Kw KуKрпK, рад В с , тогда

2ξw Tw 2ξT KwKдгT1 , с;

1 KwKдг

тогда характеристическое уравнение передаточной функции будет иметь вид 5,184 10 3s2 28,03s 50 0, а его корни равны: λ1 5404 и λ2 1,784. Передаточную функцию (1) можно представить в виде

Так как T2 T1 , то можно считать СС по рассматриваемому выходу устойчивым апериодическим звеном с коэффициентами

6

Полученный результат аналитического анализа проверяем моделированием в программном комплексе (ПК) МВТУ.

Для того чтобы провести анализ переходных характеристик

вПК МВТУ, необходимо выполнить следующие действия.

1.Переместите курсор на кнопку «Новый» и сделайте одно-

кратный щелчок левой клавишей мыши: откроется чистое «Схемное Окно». Переместите курсор на закладку «Источники входных воздействий» и щелкните один раз левой клавишей мыши: Вы ини-

циализировали соответствующую библиотеку типовых блоков. Переместите курсор на блок «Постоянное воздействие» (под-

пись «Константа») и сделайте однократный щелчок левой клавишей мыши: фон блока в «Линейке» изменился. Это означает, что блок можно переносить в «Схемное Окно». Переместите курсор в верхний левый угол «Схемного Окна» и щелкните один раз левой клавишей мыши: в поле «Схемного Окна» появился переносимый

блок.

Переместите курсор на закладку «Операции математические» и выполните однократный щелчок левой клавишей мыши: Вы ини-

циализировали соответствующую типовую библиотеку. Переместите курсор на блок «Усилитель» и сделайте однократный щелчок:

фон блока изменился. Переместите курсор в поле «Схемного Окна» на место, где Вы хотите расположить «Усилитель» и щелкните один раз левой клавишей мыши: перенос блока «Усилитель» в «Схемное

Окно» выполнен.

Переместите курсор на закладку «Динамические звенья», инициализируйте ее, перенесите Апериодическое звено в «Схемное Ок-

но» по вышеописанной процедуре приблизительно на желаемые места.

Переместите курсор на закладку «Данные», инициализируйте данную библиотеку типовых блоков, перенесите блок «Временной

график» в «Схемное Окно» примерно на желаемое место.

2. Далее проведите линии связи на структурной схеме. Для этого переместите курсор на выход блока «Константа», сделайте щел-

чок левой клавишей мыши и, отпустив клавишу, протяните горизонтальную линию связи к верхнему входному порту «Усилителя».

Снова сделайте однократный щелчок левой клавишей: на верхнем

7

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

Если требуется сделать поворот на 90 в линии связи, выполните щелчок левой клавишей мыши и продолжайте проведение линии связи в новом направлении.

Если Вы желаете прервать процедуру проведения линии связи (например, в случае ошибочного соединения блоков), нажмите клавишу Shift и, не отпуская ее, сделайте щелчок левой клавишей мыши: линия оборвется. Далее можно удалить эту линию: выдели-

те ее (щелчок левой клавишей мыши по линии) и затем удалите линию с помощью командной кнопки «Вырезать» (пиктограмма

«ножницы»).

Для задания узла (ответвления) переместите курсор на линию связи, нажмите на клавиатуре клавишу Ctrl и, не отпуская ее, сделайте щелчок левой клавишей мыши: появится темная точка. От-

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

Описанные процедуры пояснены схемами на рис. 3.

Рис. 3. Набор структурных схем в пакете МВТУ

3.Сохраните введенную часть задачи.

4.Переместите курсор на блок «Усилитель» и сделайте двойной

щелчок левой клавишей мыши: откроется диалоговое окно этого блока с активной закладкой «Параметры». Инициализируйте

диалоговую строку, введите требуемый коэффициент усиления и нажмите на кнопку «Да». Повторите аналогичные процедуры для апериодических блоков и введите соответствующие значения T .

8

5. Переместите курсор на командную кнопку «Параметры рас-

чета» и сделайте однократный щелчок левой клавишей мыши: откроется диалоговое окно «Параметры интегрирования» с активной закладкой «Основные». При активной закладке «Основные» выберите численный метод, например, Адаптивный 1. Далее введи-

те: время интегрирования 10 с; минимальный шаг интегрирования 0,001 с; максимальный шаг интегрирования 0,1 с; шаг вывода результатов 0,1. Параметр точности 0,001. Закройте диалоговое окно, щелкнув левой клавишей мыши по кнопке «Да».

6.Переместите курсор на командную кнопку «Продолжить»

ищелкните левой клавишей мыши: Вы запустили созданную задачу

на счет. По окончании расчета появится специальное окно «Сообщения», например, с информацией: «Ошибка: Заданная точность не обеспечивается».

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

ните минимальный шаг интегрирования на новое значение (1e-7) и повторите процесс моделирования.

После завершения моделирования для получения графика пе-

реходной функции исследуемой системы дважды щелкните левой клавишей мыши по блоку «График» на выходе структурной схемы

исследуемой системы.

В рассматриваемом примере процесс апериодический. Моделирование второй системы показывает, что малоинерционная составляющая не оказывает влияния на вид переходного процесса, и ею можно пренебречь и использовать эту математическую модель на начальном этапе синтеза КСТН.

Этап 2. Анализ КСТН

Вторым этапом работы является анализ КСТН, описанной моделью третьего порядка.

Расчетная структурная схема контура теленаведения с упрощенной системой стабилизации приведена на рис. 4.

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

9

Рис. 4. Расчетная структурная схема КСТН третьего порядка

формулировке при синтезе системы, будем рассматривать вход

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

Полученная система представляет собой систему третьего порядка с астатизмом второго порядка. Поэтому использовать ее для определения мгновенного промаха неудобно, так как надо на вход подавать сигнал второго порядка. Для определения желаемой величины коэффициента усиления по ошибке Kh используем другую систему, где входом является требуемое нормальное ускорение, а выходом — промах (рис. 5). Здесь:

Минимально допустимое значение коэффициента Kh, В (рад м) рассчитывается для установившего предельного состояния КСТН по формуле

Рис. 5. Расчетная структурная схема для определения Kh

10