Scientific journal
Modern high technologies
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,940

INVESTIGATION OF MULTI-CONNECTED SYSTEM OF AUTOMATIC CONTROL WITH DELAY DECOMPOSITION METHOD

Ilyasov B.G. 1 Saitova G.A. 1 Elizarova A.V. 1
1 Ufa state aviation technical university
The problem of controlling objects with time delay is quite complicated. The presence of a delay in the control loop leads to an increase in the phase shift, which can provoke instability of the closed system, including in the presence of small gain factors of the regulator. The problem is get worse if the control object is multidimensional. In article investigates a multi-connected automatic control system separately and concurrently with delay in direct and cross communication channels. The estimation of the critical value of the delay on the basis of a systematic approach to the description of a multi-connected system in combination with frequency methods of stability analysis is proposed. The technique of estimation of the critical value of delay in direct and cross-links of a multi-connected automatic control system using the decomposition method is developed. Another problem solved in this article is to determine the stability of a multi-connected system with delay in subsystems and connections using a system approach to describe a multi-connected automatic control system through the characteristics of connections and characteristics of subsystems. In article presents examples that demonstrate the use of the proposed techniques and to find the critical value of the delay, and to assess the stability of the MCSAC. The results of the study were confirmed by simulation using MATLAB SIMULINK package.
delay
stability
multi-connected automatic control system
cross-links
direct channels
decomposition method
critical delay
characteristic equation

В работе исследуется многосвязная система автоматического управления (МСАУ), где одновременно осуществляется регулирование нескольких взаимосвязанных координат и наблюдается явление запаздывания. Далее будем рассматривать однотипную МСАУ – многосвязная система с одинаковыми подсистемами.

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

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

ily01.wmf (1)

где ily02.wmf, X(s), U(s) – векторы задающих, регулируемых, управляющих воздействий соответственно; ily03.wmf – матричная передаточная функция (МПФ) многомерного объекта по управляющим воздействиям, с запаздыванием; ily04.wmf – МПФ сепаратных регуляторов.

С целью изучения МСАУ с запаздыванием используется подход, представленный в труде [1], при котором линейная МСАУ описывается на уровне физических подсистем и многомерных элементов связи между ними.

Цель исследования: разработка методов оценки критического значения запаздывания в прямых и перекрестных связях МСАУ с помощью метода декомпозиции, на основе описания системы через характеристики связей и характеристики подсистем [1].

Материалы и методы исследования

Рассмотрим однородную МСАУ с запаздыванием в подсистемах. Поскольку подсистемы однотипные, то их передаточные функции (ПФ) можно представить как

R1W11e-ts = R2W22e-ts = ... = RnWnne-ts = RWe-ts, (2)

а ПФ индивидуальных характеристик подсистем примут вид [1]:

ily05.wmf ily06.wmf ily07.wmf (3)

Для полной МСАУ, состоящей из n подсистем и соответствующей системы уравнений (1), характеристика связи (ХС) в общем виде между k подсистемами имеет вид

ily08.wmf, ily09.wmf ily10.wmf ily11.wmf, ily12.wmf.

Характеристическое уравнение однотипной МСАУ имеет вид

ily13.wmf. (4)

Перепишем равенство (4) сравнительно разомкнутых подсистем W(t, s) = RWe-ts в общем виде для удобства вычисления критического значения запаздывания tк [2]:

ily14.wmf (5)

Проанализируем уравнение связей

ily15.wmf (6)

относительно переменной x, что выходит из (5) с помощью подстановки W(t, s) = x, где запаздывание в подсистемах.

Перенеся на комплексную плоскость годограф функции W(jω) c запаздыванием, а также корни характеристического уравнения связи (6) ily16.wmf, можно найти критическое значение запаздывания τiкр [3].

Критическое запаздывание находится в точке пересечения корня уравнения (6) ily17.wmf и годографа подсистемы W(jω). В МСАУ для нахождения критического значения запаздывания tкр необходимо, чтобы годограф W(jω) построенный на одной комплексной плоскости с корнями уравнения (6), проходил через ближайший из них и не охватывал при этом другие корни, то есть МСАУ оказывается на границе устойчивости. Значения ωiср (частота среза) и τiкр определяются из выражения: ily18.wmf. Из этого условия получаем систему из двух уравнений относительно t и ω0 [4]:

|W(t,ω0)| = |ily19.wmf|, arg W(t,ω0) = arg ily20.wmf. (7)

Минимальным из найденных значений ti является критическое значение запаздывания ily21.wmf, ily22.wmf [4].

Результаты исследования и их обсуждение

Пример 1. Многосвязная система автоматического управления с запаздыванием в подсистемах.

В качестве примера рассмотрим дву- связную систему. Объектом исследования является многомерный объект с запаздыванием в подсистемах. ПФ разомкнутой подсистемы равна

ily23.wmf,

характеристика связи h2 = 0,2 [5].

Характеристическое уравнение системы:

ily24.wmf.

Корни уравнения связи ily25.wmf равны:

x1 = –0,56 + 0,5j; x2 = –0,56 – 0,5j.

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

Так как годограф с запаздыванием W(jω) пересекает один из корней уравнения связи, следовательно, система на грани устойчивости. С помощью графоаналитического способа найдено критическое значение запаздывания τкр = 0,15 сек.

iljs1.tif

Рис. 1. Годограф МСАУ с критическим значением запаздывания

iljs2.tif

Рис. 2. Годограф МСАУ с запаздыванием в подсистемах, построенный на плоскости корней уравнения связи

Чтобы определить устойчивость МСАУ с запаздыванием в подсистемах, воспользуемся частотным критерием. По характеристическому уравнению (4) получаем уравнение связи

ily26.wmf (8)

На одной комплексной плоскости строим годограф Ф*(jw,t) (подсистема с запаздыванием) для всех w∈(0, +∞) и корни уравнения связи (8). Согласно критерию устойчивости для многомерных систем необходимо и достаточно, чтобы годограф амплитудно-фазовой характеристики подсистем Ф*(jw,t), для всех w∈(0, +∞), построенный на плоскости корней уравнения связи, не охватывал ни один из его корней (рис. 2). На рис. 2 показана устойчивая МСАУ, где Ф(jw) – годограф без запаздывания.

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

Рассмотрим однородную МСАУ с запаздыванием в перекрестных связях. Для данного вида системы характеристическое уравнение однотипной МСАУ:

ily27.wmf

где

ily28.wmf, ily29.wmf,

ily30.wmf.

Введем переменную Wp, которая обозначает разомкнутую систему, и перепишем относительно нее характеристическое уравнение. Так как Wp = R(s)W(s), следовательно характеристическое уравнение МСАУ примет вид [6]:

ily31.wmf (9)

Уравнение связи относительно переменной Wp = x (10):

ily32.wmf (10)

Решением уравнения связи (10) будут корневые годографы xi(jω) [3, 4].

Пример 2. Рассмотрим МСАУ с двумя подсистемами и запаздыванием в перекрестных связях. ПФ подсистемы равна ily33.wmf, характеристика связи ily34.wmf.

Корни уравнения связи (10) при t = 0 в перекрестных связях и w = 0: ily35.wmf (10) равны x1 = –0,5611 + 0,5442i; x2 = –0,5509 – 0,4537i.

iljs3.tif

Рис. 3. Годограф МСАУ с запаздыванием в перекрестных связях

При постепенном увеличении значения запаздывания t видим, что корневые годографы xi(jω) не пересекают годографы подсистем. Это означает, что запаздывание не влияет на устойчивость системы (рис. 3). На устойчивость разомкнутой системы можно повлиять, изменив коэффициенты характеристик связей h2.

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

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

В общем виде характеристическое уравнение МСАУ с однотипными подсистемами и запаздыванием одновременно в подсистемах и перекрестных связях имеет вид

ily36.wmf

где

ily37.wmf, ily38.wmf,

ily39.wmf.

Запишем по аналогии характеристическое уравнение МСАУ относительно разомкнутых подcистем:

ily40.wmf (11)

и уравнение связи

ily41.wmf. (12)

Далее, если построить на комплексной плоскости годограф разомкнутых подсистем Wр(jω,t) и корневой годограф уравнения связи (12) , то можно найти критическое значение запаздывания τiкр и исследовать систему на устойчивость (рис. 4).

iljs4.tif

Рис. 4. Годограф МСАУ с запаздыванием в прямых и перекрестных связях

Значения τiкр определяются из выражения

y1 = φо(ω) – ω1срτiкр, τiкр = (φо(ω) – y1)/ω1ср,

где y1 – argx1, φо(ω1ср) – argWp(jw) на частоте ωiср = 1.

Минимальным из найденных значений ti – критическое значение запаздывания τкр = min{τiкр}, i = 1,n. В примере 2 критическое значение запаздывания равно τкр = 0,087 сек.

Заключение

В ходе выполненной работы применен способ оценки критического значения запаздывания в прямых и перекрестных связях МСАУ с помощью метода декомпозиции.

С помощью известного критерия устойчивости для МСАУ рассмотрена ее устойчивость с запаздыванием в подсистеме, перекрестных связях. Также рассмотрен алгоритм определения критического значения запаздывания для МСАУ по разомкнутым и замкнутым подсистемам. Проверена эффективность разработанной методики нахождения критического значения запаздывания tкр. Результаты и проверка системы на устойчивость подтверждены моделированием МСАУ по результатам переходных процессов.

Работа выполнена при поддержке Российского фонда фундаментальных исследований (Гранты РФФИ № 18-08-00702 А, 18-08-01299 А).