УДК 004.942
DOI: 10.14529/mmph240406
МАТЕМАТИЧЕСКОЕ ОБЕСПЕЧЕНИЕ ИЕРАРХИЧЕСКОГО ЦИФРОВОГО УПРАВЛЕНИЯ СЛОЖНЫМ ТЕХНОЛОГИЧЕСКИМ ОБЪЕКТОМ
Е.В. Ларкин1, А.В. Богомолов2, А.Н. Привалов3, Т.А. Акименко1
1 Тульский государственный университет, г. Тула, Российская Федерация
2 Федеральный исследовательский центр «Информатика и управление» Российской академии наук, г. Москва, Российская Федерация
3 Тульский государственный педагогический университет им. Л.Н. Толстого, г. Тула, Российская Федерация
E-mail: [email protected]
Аннотация. Исследуется вопрос синтеза математического обеспечения иерархического цифрового управления сложным технологическим объектом на основе универсальных контроллеров с гибким программированием. Для преодоления ограничений на быстродействие контроллеров, функционирующих в режиме реального времени, при усложнении законов управления и расширении классов управляемых объектов предложено разделение контроллеров на стратегический, тактический и функционально-логический иерархические уровни. Определены функции и построены алгоритмы обработки данных контроллерами каждого иерархического уровня. Показано, что согласованное функционирование исследуемой системы цифрового управления сводится к обмену данными между контроллерами разных уровней, осуществляемому программно, одинаково для контроллеров всех уровней. Получены зависимости для оценки периода перезапуска управляющих программ, работающих по гамильтонову циклу, а также задержки передачи данных для случая, когда на контроллерах тактического и функционально-логического уровней реализован алгоритм пропорционально-интегрально-дифференцирующего регулятора. Показаны особенности управления режимами функционирования контроллеров тактического и функционально-логического уровней за счет реализации на стратегическом уровне нейронной сети.
Ключевые слова: цифровое управление; система управления; иерархические уровни управления; контроллер с гибким программированием; гамильтонов цикл; полумарковский процесс; нейронная сеть.
Одним из ключевых направлений развития современного производства является цифровиза-ция промышленности, в том числе широкое внедрение цифровых систем управления, в качестве аппаратной платформы которых используются универсальные контроллеры с гибким программированием [1-3]. Помимо производства такие контроллеры широко применяются в киберфизи-ческих гетерогенных системах управления робототехническими комплексами специального назначения [4-10].
Реализация алгоритмов управления на универсальных контроллерах с гибким программированием позволяет управлять объектами практически любой структуры, не изменяя аппаратной части системы [11, 12]. Однако усложнение законов управления и расширение класса управляемых объектов ограничивается быстродействием контроллеров, функционирующих в режиме реального времени [13-15]. Это приводит к необходимости разделения задач управления по степени влияния на общую цель управления на стратегические, тактические и функционально-логические задачи [16]. В контроллерах с высоким быстродействием такое разделение может быть виртуальным, а контроллеры с ограниченным быстродействием должны быть организованы в иерархическую структуру с аналогичным разделением. Усложнение структуры приводит к удорожанию аппаратной части системы цифрового управления, поэтому аппаратные средства каждого уровня иерархии должны использоваться максимально эффективно, для чего необходимо четко определить круг задач и задать ограничения на реальное время их решения контроллерами каждого типа [17, 18].
Решение таких задач требует специального математического обеспечения иерархического цифрового управления сложным технологическим объектом, что и определило цель исследования.
Структура системы иерархического цифрового управления сложным технологическим объектом
Структура системы иерархического цифрового управления сложным технологическим объектом приведена на рис. 1.
В состав системы (рис. 1) входят: контроллер стратегического уровня (SLC-контроллер), контроллер тактического уровня (TLC-контроллер) и K контроллеров функционально-логического уровня (FLLCt-контроллеры, ke{1, ..., K} [19]. Разделение контроллеров по уровням может быть реальным или виртуальным. Например, SLC- и TLC-контроллеры могут быть реализованы на одной физической ЭВМ с виртуальным разделением по уровням, а FLLCk-контроллеры могут объединяться в группы, каждая из которых реализуется на одном физическом цифровом контроллере.
При управлении объектом (OUC) на каждом иерархическом уровне реализуется своя часть интегрального алгоритма управления.
На верхнем, стратегическом уровне формируется цель управления W(t) в виде развивающейся в физическом времени последовательности изменений вектора состояния OUC, x(t) = {xk(t)}, которая обеспечивает, например выпуск возможно большего объема продукции в единицу времени, поддержание режима функционирования технологического оборудования, согласованное взаимодействие компонентов робототехнических систем и комплексов специального назначения и т. п. [20-24].
Из глобальной цели W(t) TLC-контроллер формирует вектор u(t) = {uk(t)} текущих желаемых значений вектора x(t) состояний OUC, который передается на FLLCk-контроллеры. Кроме того, на эти контроллеры с выходов SLC-контроллера могут быть непосредственно переданы значения настроек \\ (/) параметров регуляторов, необходимые для реализации оптимальных режимов функционирования приводов Dk с сенсорами SLk в контурах локальных обратных связей.
По вектору желаемых значений состояний u{t) и результатам измерений q(i) = [^¿у, (i) = {qk(t)!
[25, 26] с помощью сенсоров S¡,k физических перемещений q(t)={qk(t)} приводов Dk формируются воздействия на объект управления, которые приводят к изменению его состояния x(t). Результаты измерения в виде вектора x(t) = {xk(t)} возвращаются на уровень TLC, обеспечивая обратную связь по регулируемым параметрам технологической установки. Кроме того, вектор
[u(V), х(7)] = [j/j (7) , формируемый как объединение векторов u(t) и х (/ ). возвращается в SLC
для определения цели управления w(t) на следующем шаге функционирования системы.
Также на уровень SLC передаются элементы векторов z(t) = {zk(t)} и r(t) = {rk(t)}. Элементы вектора z(t) напрямую не используются в процессе управления, но дают дополнительную информацию о состоянии объекта, например, о температуре, давлении и влажности в приобъектной зоне, о срабатывании концевых выключателей, аварийных клапанов и т. п. Элементы вектора r(t) представляют собой логические переменные, которые несут информацию о результатах отработки системы предыдущего шага управления. Совместно элементы векторов \ ( / ), «(/). z(t), r(t) образуют пространство принятия решений в системе управления объектом OUC.
Математическое обеспечение алгоритмов, реализованных в контроллерах
SLC-, TLC- и FLLCk-контроллеры работают по типовому циклическому алгоритму общего вида, показанному на рис. 2, а, состоящему из неисполнимых операторов a0(a) «Начало» и aJ(a)+J «Конец» и исполняемых операторов aj(a), je{1, ..., J}. Циклическим алгоритм делает обратная связь [aJ(a)+1, aj(a)], показанная на рис. 2 пунктирной стрелкой.
Рис. 1. Структура системы иерархического цифрового управления сложным технологическим объектом
Исполнимые операторы реализуют:
в алгоритме SLC-коитроллера - функции ввода данных векторов х(/), «(/). z(t), /*(/). расчет элементов вектора w{t) и передачу рассчитанных значений в TLC, расчет элементов вектора настроек w (/) параметров регуляторов и передачу рассчитанных значений в FLLQ-контроллеры;
в алгоритме TLC-контроллера - функции: ввода векторов х (/) и w(t), расчета и вывода в FLLQfc-контроллеры элементов u(t). расчета и возврата в SLC-контроллер элементов вектора г(1): в алгоритме каждого F LLQ-ко нтро л лс ра - функции ввода значений uk и qk , расчета и вы-
b
Рис. 2. Структура типового алгоритма контроллеров всех уровней: a - общего вида; b - линейно-циклического вида
Вследствие того, что SLC-, TLC- и FLLCt-контроллеры интерпретируют свои алгоритмы в реальном физическом времени t последовательно, оператор за оператором, между обращениями к внешним по отношению к рассматриваемому контроллеру устройствам естественным образом образуются временные интервалы (задержки) [27-29]. Наличие интервалов времени между транзакциями, генерируемыми каждым контроллером, должно учитываться при проектировании программного обеспечения контроллеров системы. Универсальным подходом к оцениванию временных интервалов является построение полумарковской модели типового алгоритма [30-32]:
ц = {A, h (t)},
где A = {a0(a), a1(a), ... aJ(a), aJ(a)+\} - множество состояний, совпадающее с множеством исполнимых и неисполнимых операторов; h (t )= a^(t) = p ® f (t) - полумарковская матрица, описы-
вающая временные и стохастические свойства алгоритма; p =
Pj( a),n( а)
- стохастическая мат-
рица; Г (7) = ^ ) - матрица плотностей распределения времени пребывания процесса в
состояниях множества А до переключения в сопряженные состояния, причем размерности всех матриц [7(а)+2] х [.Да)+2].
Полумарковская модель правильно синтезированного алгоритма обладает следующими свойствами:
для каждого состояния а^еЛ существует хотя бы один путь а0^а^^а](а) и хотя бы один путь ада^^ а^уто есть полумарковский процесс при замкнутой обратной связи является возвратным;
вероятности переключения между состояниями ц-ц удовлетворяют условиям
J ( а)
^ Pj( a ),n( a ) n(a)=0(a)
1, j ( a )* J ( a ) +1,
or j (a) = J (a) +1, feedback is closed; 0, j (a) = J (a) +1, feedback is open; на время пребывания в состояниях полумарковского процесса накладываются ограничения
да
0 < Tj(a),n(a) min - arg fj(a),n(a) (t)-Tj(a),n(a) max < , j fj(a),n(a) (t) dt = 1;
0
в модели алгоритма с замкнутой обратной связью [Да)+1]-я строка полумарковской матрицы h(t) имеет вид fJ(a)+1,0(a)(t),..., 0,..., 0 , где fJ(a)+10(a)(i) - плотность распределения времени
работы операционной системы контроллера по перезапуску исследуемого алгоритма.
Полумарковский процесс является возвратным, плотности распределения времени пребывания в его состояниях имеют область определения, локализованную в положительной полуплоскости, следовательно, процесс является эргодическим [30-32]. Плотность распределения времени возврата в состояние a0(a) полумарковского процесса и и среднее время периода равны соответственно
f (t) = I0(a)L~l
[h'(t )]f
.1=1
IJ ( a )+1* fj ( a )+1,0(a)(t) ; Tc = jVc (t) dt ,
где / ^) - плотность распределения времени от начала текущего цикла до начала следующего цикла работы программы; Ь [...], Ь [...] соответственно прямое и обратное преобразование Лапласа; ^^ - [Да)+2]-элементный вектор-строка, нулевой элемент которого равен единице, а остальные элементы равны нулю; - [Да)+2]-элементный вектор-столбец, последний,
[Да)+1]-й элемент которого равен единице, а остальные элементы равны нулю; ^^) - полумарковский процесс, описывающий алгоритм с разомкнутой обратной связью, у которого [Да)+1]-я строка имеет вид [0, ..., 0]; * - обозначение операции свертки.
На структуру исследуемого алгоритма может быть дополнительно наложено следующее ограничение: в течение одного цикла каждый канал, связывающий один контроллер с другими контроллерами и периферийным оборудованием, опрашивается по одному разу, а, следовательно, исполнимые операторы а](а) образуют гамильтонов цикл на графе (см. рис. 2, а).
Схема линейно-циклического алгоритма, представляющего один из возможных гамильтоно-вых циклов, приведена на рис. 2, Ъ. Для полумарковского процесса с гамильтоновым циклом
/о « ) = /0{а)(< ) */ а)(* ) *...*/,( а)(* ) *..*/Ла)(< ) *Л(а)+1 ^ ) ; То = ТА а) ,
А а)=0(а)
где /да^ (t) и Гд^ - плотность распределения и математическое ожидание времени пребывания
полумарковского процесса в состоянии aj(a) гамильтонова цикла до его переключения в следующее состояние.
Математическое обеспечение обработки данных в ЕЬЬС^-контроллере
Реализация полумарковской модели обработки данных на нижнем (функционально-логическом) уровне иерархии в системах с механизмами, в которых присутствует сухое трение, и естественное ограничение на скорость перемещения терминального элемента (руля, задвижки, клапана и т. п.), а также инерционность сенсора З^, дает следующую систему дифференциальных уравнений третьего порядка в форме Коши:
0
&Як ,1 (г)
= Чк, 2 (г);
&
&Як,2 (') - Г- м
-^-= "«к, 1 • Як,2 - ак,2 • Як,2 (г)] + «к,V (I),
= ~ак,4Як{{) + Як,1 (')>
<%к (О
л
где qk | (/) = с]к (/) - физическое положение терминального элемента; цк - значение положения, измеренное сенсором ; Як 2 - скорость перемещения терминального элемента с учетом ограничения; ак 2 • [Чк 2 (г- сила сухого трения; ак 1 ^ ак 4 - коэффициенты, определяющие динамику и точность работы механизма исполнительного устройства и сенсора.
Величины Чк,2 и Як,2 связаны соотношениями Чк,2 = Як,2 при Як,2ш1п - Як,2 - Як,2тшах ; Як,2 = Як,2ш1п при Як,2 - Як,2шт ; Як,2 = Як,2тах при Як,2 - Як,2тах .
Наиболее часто на практике при обработке данных в FLLCJt-контроллере применяется алгоритм пропорционально-интегрально-дифференцирующего (ПИД) регулятора [33, 34], который имеет вид
Vс,к (г) = «РЬ,Р,к (^Р,к ) • %Ь,к (г) + «РЬ,Б,к (*Ъ,к ) Х ^^ (г) + «РЬ,1п,к (^1п,к ) • | %Ь,к (Г)
ш О
%Ь,к (*)=ик (0~Як(0'>
где ~ ошибка в установлении заданного положения цк (/) терминального элемента;
«рь,р,к (^Р,к ), «рццк () , «рь,1п,к (^1п,к ) - коэффициенты при пPопорционaльной, интегральной, дифференциальной составляющих соответственно, определяющие качество переходных процессов на FLLQrуровне; wc к =[^рк , , ~№]пк J - вектор данных, вводимых в FLLQr
контроллеры с SLC-контроллерa для установления оптимальных режимов управления OUC.
Каждый FLLQrконтроллер является цифровым, поэтому производная и интеграл в алгоритме ПИД-регулятора определяются соответственно по формулам конечной разности и конечной суммы [35, 36], т. е.
(() = %ь,к (п)- %ь,к (п -1); |%ь,к (О-= Е%ь,к 0'),
& о ¿=0
где п - дискретное время, связанное с реальным физическим временем соотношением г = Тс РЬ кп;
Тс РЬ к - период опроса FLLQr контроллеров.
Для вычисления производной каждый FLLCJrконтроллер должен осуществить две транзакции с ^С-контроллером и две транзакции с сенсором Б1:к. Вследствие дополнительных ограничений, которые накладываются на алгоритмы обработки данных по гамильтонову циклу, на это необходимо затратить 2Тс тЬс к единицы физического времени. Поэтому управляющее воздействие на выходе РЬЬС^-контроллера появляется с дополнительной задержкой по времени относительно момента поступления данных ик (/) и с[к (/), т. е.
Ук (г) = Ус, к (г " 2Тс,РЬ,к ) .
Накопление конечной суммы происходит в течение всех периодов работы контроллера и не оказывает существенного влияния на время чистого запаздывания в контуре обратной связи.
Математическое обеспечение обработки данных в ТЬС-контроллере
Объект управления, в котором ^^контроллер обеспечивает заданный режим функционирования, в общем случае описывается системой нелинейных дифференциальных уравнений:
t x(t) dm{x)x(t) dM(x)x(t) q(t) (t) dM(q)q(t)
^ x(t), , m(x) ,.M(x) ^(t),
= 0, 1 < к < K,
х) "'" (х) " ^) "'" (
где q(t) - вектор воздействий, формируемый на уровне FLLCk-контроллеров; х(0 - вектор параметров состояния ОиС, измеряемых сенсорами «к,2 и вводимых в TLC-контроллер.
С целью унификации программного обеспечения при расчете элементов вектора ы(\) применяется алгоритм ПИД-регулятора, который при отсутствии алгоритмически создаваемых перекрестных связей между каналами управления имеет вид
. (I:) = «^С,Р,к ' STLC,k ^) + а1ЪС,Дк & ^ + а^С,1п,к ' |STLC,k (т)
хс,к (t ) = Щ ^) —x ^ (t) ,1 < к < К, где £тьс к (t) - текущая ошибка установления заданного значения хк ^) ; «тьсРк, «тьС Р к, «тьс Р к - коэффициенты передачи ПИД-регулятора по пропорциональной, дифференциальной и интегральной составляющим; Щ ^) - цель управления, определяемая на стратегическом уровне как желаемое значение регулируемого параметра хк ^).
Производная и интеграл в алгоритме ПИД-регулятора также рассчитываются по формулам конечной разности и конечной суммы, поэтому между сигналами х(0 и и(0 возникает задержка по времени, равная 2То тьс к . Эта задержка может быть отнесена к сенсорам к, поскольку перекрестные связи между каналами на уровне алгоритма функционирования ^С-контроллера отсутствуют. С учетом задержек по времени модель преобразования сигналов х(0 сенсорами принимает вид
^к {} ~ 2^с,ТЬС„к) ~ ( ОТ М
аЪ,к, 1 Хк V ~~ ¿*с,Т1С„к ) ~ а$,к,2Хк У) >
где «8 к 1 и «8 к 2 - параметры, определяющие быстродействие и коэффициент передачи к и З2,к соответственно.
Кроме интервала 2То тьс к при работе ТЬС-контроллера возникает такое явление, как перекос данных. Без нарушения общности можно считать, что внутри гамильтонова цикла первой транзакцией является ввод элемента хх (/ ), после чего в порядке, определенном алгоритмом, вводятся оставшиеся элементы указанного вектора и выводятся элементы вектора г(/). Задержки по времени между началом ввода х, и началами ввода х2, ..., хк, ..., хК порождают перекос данных
ТшЛ = Тш,2-> Тт •••> тт,к ПРИ вв°Де- Задержки по времени 7(Ш1| = 0, 7(ш12,..., ТоЫ к,..., ТоЫЛ между началом вывода элемента V] вектора и началами выводов элементов у2 , ..., у* , ..., ук порождают перекос данных при выводе v(t). Временной интервал Ты к 1, алгоритмически формируемый между вводом хк и выводом V;, дает часть времени чистого запаздывания формирования 1-го элемента вектора управления относительно к-го элемента вектора состояния. Полное время чистого запаздывания определяется порядком следования операторов в гамильтоновом цикле.
Вектор г(0 = {гк(0} формируется как вектор подтверждения исполнения команды W(t). Элементы этого вектора обнуляются после передачи в 8ЬС-контроллер и устанавливаются в единицу, если к началу следующей транзакции величина ошибки £т]^с к ^) = Щ ^ ) — хк ^) выполнения команды не превышает заданного порога, то есть
г* =1 , ^Ъеп ©тщк (t) < ©тщк , где ©тьс к - пороговое значение сигнала ошибки по к-му регулируемому параметру ОиС, хк (t) .
Записанные выражения представляют собой полную модель функционирования ТЬС-контроллера в составе иерархической системы управления.
Математическое обеспечение обработки данных в SLC-контроллере
В SLC-контроллере (на стратегическом иерархическом уровне) реализуются алгоритмы, обеспечивающие эффективные режимы работы оборудования. Одна из возможных реализаций подобного алгоритма - это нейронная сеть, которая используется для настройки параметров управления FLLQt-контроллеров [3]. Настройка осуществляется по результатам наблюдения состояния OUC, TLC, Dk. Для принятия решения о настройке параметры, описывающие указанные блоки иерархической системы, в виде векторов u(t), x(t), r(t) передаются с соответствующих уровней в SLC. Дополнительно к перечисленным на этот уровень передаются данные о нерегулируемых параметрах OUC, обозначенные на рис. 1 как вектор z(t). Все переданные данные обра-
yi(у)уj{у)'•••' (y)
принятия реше-
зуют параметрическое пространство y = u ^ x ^ г ^ z =
ния. Решение заключается в установлении значений элементов векторов W(t) и w (t) , передаваемых на уровни TLC- и FLLC^-контроллеров соответственно. Вектор W (t ) = [^(t),..., W (t),..., WK (t)] определяет текущий кластер желаемых значений состояния
OUC. Вектор w(t) определяет векторы wk (t) = [wPk,..., wDk,..., wInk J настроек «^Рк (wPk), ...,
«FL, D,k (wd,к ) , «FL,In,к (WIn,k )] параметров ПИД-регуляторов.
Для выбора векторов W(t) и w(t) в нейронную сеть уровня SLC-контроллеров включаются нейроны, формирующие в пространстве параметров y систему из J(c) гиперплоскостей:
J ( у )
Z bj(0),j(y)yj(y) = Cj(c) j J (c) , j( y )=1 y )
где и - параметры j (c) -го нейрона, задаваемые при настройке нейронной сети.
При подстановке в систему уравнений J(c) гиперплоскостей вектора yt = У (t) = У1(y),t,. ., yj(y),t,. ., yj(y),t текущих оценок параметров формируются значения Сдсу
правых частей линейных уравнений J(c) гиперплоскостей.
J ( У )
Z bj(c),j(y)yj(y),t = Cj(c),t,1(С) ^ j (c) ^ J (С) ; j( y )=1( y)
В том случае если t — > 0, то точка, определяющая состояние иерархической системы управления, находится над плоскостью J(c) гиперплоскостей; если c^ t — < 0, то под плоскостью, а если t — = 0, то принадлежит плоскости. Поэтому в нейронную сеть вво-
дится преобразование вектора yt в вектор о = следующим образом:
°1(c) ,...,aj(c) (c )
логических переменных
_Jl,whencj{c)t -С.(с)>0;
С / \ — s
c) [0 otherwise, где 0, 1 - логический ноль и единица соответственно.
Дизъюнктивные нормальные формы, построенные из элементов вектора с, позволяют разбить пространство y на S областей D1,..., D^,..., D, каждой из которых ставится в соответствие
^-е значения векторов ^ W (t) и ^ w (t), которые дают оптимальные настройки системы управления в этом конкретном режиме эксплуатации OUC. Дизъюнктивные нормальные формы имеют вид:
j(.
Л j(c )(j ))
, 1 <Í<H,
[J
Ас) / \
где /Л - дизъюнктивная нормальная форма, описывающая область 1),: ^ л ^ Л, и Д0",^.) j -
и-я конституента единицы ^-й нормальной формы; Л^ „ дс) ) - логическая функция переменной при ее вхождении в и-ю конституенту ^-й нормальной формы.
Для нормальной работы нейронной сети на дизъюнктивные нормальные формы, описывающие области 1):, накладывается ограничение:
Дг О Ё>£ = 0, если £ Ф £ . Поэтому при принятии управленческого решения из всех областей, на которые разбивается пространство параметров, выбирается та единственная область 1):, для которой /Л = 1 .
Настройка нейронной сети сводится к вариации параметров и Сдс), определяющих
расположение /(с) гиперплоскостей в пространстве параметров у, и выбору таких значений векторов Щ^) и w(t), которые обеспечивают оптимальный режим функционирования ОИС в этом состоянии.
Заключение
Таким образом, в результате исследования сформированы требования к структуре иерархической системы управления сложными многоконтурными объектами. Для контроллеров каждого уровня предложены типовые алгоритмы реализации их функций и определены характеристики контроллеров как физических приборов, выполняющих эти функции. Показано, что реализация оптимальных режимов эксплуатации объекта управления возможна за счет реализации на верхнем иерархическом уровне нейронной сети. Дальнейшие исследования в этой области должны быть направлены на представление контроллера с иерархической структурой как системы массового обслуживания, каждый элемент которой является одновременно и обслуживающим прибором, и генератором потока запросов на обслуживание, что позволит оптимизировать информационные потоки внутри системы иерархического цифрового управления сложным технологическим объектом.
Работа выполнена при поддержке гранта РНФ 24-29-20256.
Литература
1. LofVing, M. Manufacturing Strategy Formulation, Leadership Style and Organizational Culture in Small and Medium-Sized Enterprises / M. LofVing, K. Safsten; M. Winroth // International Journal of Manufacturing Technology and Management. - 2016. - Vol. 30, no. 5. - P. 306-325.
2. Composite Error Learning Robot Control using Discontinuous Lyapunov Analysis / Y. Pan, K. Guo, A. Bobtsov, Ch. Yang, H. Yu // IEEE Transactions on Automatic Control. - 2024. - Vol. 69, Iss. 3. - P. 1705-1712.
3. Ле, В.Н. Обзор интеллектуальных систем управления и робототехнических задач в производстве аквакультуры / В.Н. Ле, А.Л. Ронжин // Морские интеллектуальные технологии. - 2024. -№ 1-1 (63). - С. 171-180.
4. Meshcheryakov, R. Control of hyperlinked cyber-physical systems / R. Meshcheryakov // Proc. 14th International Conference on Electromechanics and Robotics "Zavalishin's Readings". Smart Innovation, Systems and Technologies. - 2020. - Т. 154. - С. 27-33.
5. Эффективное функционирование смешанной неоднородной команды в коллаборативной робототехнической системе / Галин Р.Р., Широкий А.А., Магид Е.А. и др. // Информатика и автоматизация. - 2021. - Т. 20, № 6. - С. 1224-1253.
6. Мещеряков, Р.В. Интеллектуальные робототехнические системы / Р.В. Мещеряков // Датчики и системы. - 2024. - № 2. - С. 12-17.
7. Мещеряков, Р.В. Постановка задачи гетерогенного группового взаимодействия роботов при решении задач умного сада / Р.В. Мещеряков, А.С. Широков // Вестник Южно-Уральского государственного университета. Серия «Математика. Механика. Физика». - 2024. - Т. 16, № 2. -С. 41-49.
8. Осипов, О.Ю. Проектирование цифровых моделей элементов электромашинной части электромехатронных модулей робототехнических систем / О.Ю. Осипов, Р.В. Мещеряков, М.Г. Шепеленко // Экстремальная робототехника. - 2017. - № 1. - С. 160-164.
9. A Probabilistic Approach to Estimating Allowed SNR Values for Automotive LiDARs in "Smart Cities" under Various External Influences / R. Meshcheryakov, A. Iskhakov, M. Mamchenko et al. // Sensors. - 2022. - Vol. 22, no. 2, - P. 609. DOI: 10.3390/s22020609
10. Постановка задачи реализации модели беспилотной авиационной транспортной системы с использованием технологий искусственного интеллекта / А.Л. Смолин, В.П. Кутахов, А.А. Захарова и др. // Искусственный интеллект. Теория и практика. - 2023. - № 1 (1). - С. 41-45.
11. Landau, I.D. Digital Control Systems, Design, Identification and Implementation / I.D. Landau,
G. Zito. - Springer, 2006. - 484 p.
12. Astrom, J. Computer Controlled Systems: Theory and Design / J. Astrom, B. Wittenmark. -Tsinghua University Press, Prentice Hall, 1996. - 557 p.
13. Ivutin A.N., Larkin E.V. Estimation of Latency in Embedded Real-Time Systems // 3-rd Meditteranean Conference on Embedded Computing (MECO-2014). Budva, Montenegro, 2014, P. 236239.
14. Kilian C.T. Modern Control Technology: Components and systems / C.T. Kilian. - Thompson Delmar Learning. - 2000. - 644 p.
15. Ларкин, Е.В. Моделирование надежности бортового оборудования мобильного робота / Е.В. Ларкин, Т.А. Акименко, А.В. Богомолов // Известия Саратовского университета. Новая серия. Серия: Математика. Механика. Информатика. - 2021. - Т. 21, Вып. 3. - С. 390-399.
16. Babishin V., Taghipour S. Optimal maintenance policy for multicomponent systems with periodic and opportunistic inspections and preventive replacements // Applied Mathematical Modelling. -2016. - Vol. 40, Iss. 23-24. - Pp. 10480-10505.
17. Digital Control of Continuous Production with Dry Friction at Actuators / E. Larkin, A. Privalov, A. Bogomolov, T. Akimenko // Smart Innovation, Systems and Technologies. - 2022. -Vol. 232. - P. 427-436.
18. Сенсорные сети контроля состояния авиационной техники при испытаниях и эксплуатации / Е.С. Солдатов, А.В. Богомолов, Е.В. Ларкин, А.С. Солдатов // Авиакосмическое приборостроение. - 2024. - № 2. - С. 61-68.
19. Структура иерархической системы цифрового управления сложным технологическим объектом / Е.В. Ларкин, А.В. Богомолов, А.Н. Привалов, Т.А. Акименко // Математические методы в технологиях и технике. - 2024. - № 6. - C. 63-66.
20. Восстановление динамически искаженного сигнала на основе теории оптимальных динамических измерений / А.Л. Шестаков, А.А. Замышляева, Н.А. Манакова и др. // Автоматика и телемеханика. - 2021. - № 12. - С. 125-137.
21. Ermolov, I.L. Models for Decision Making Support Systems in Robotics / I.L. Ermolov, S.S. Graskin, S.P. Khripunov // Studies in Systems, Decision and Control. - 2021. - Vol. 352. - P. 145149.
22. Развитие теории оптимальных динамических измерений / Е.В. Бычков, С.А. Загребина, А.А. Замышляева и др. // Вестник Южно-Уральского государственного университета. Серия «Математическое моделирование и программирование». - 2022. - Т. 15, № 3. - С. 19-33.
23. Bazhenov, E.I. Software for the Mobile Robot Spatial Orientation System / E.I. Bazhenov, S.A. Mokrushin, S.I. Okhapkin // Вестник ЮУрГУ. Серия «Математическое моделирование и программирование». - 2021. - Т. 14, № 2. - С. 64-69.
24. Алгоритм численного нахождения оптимального измерения, искаженного инерционностью, резонансами и деградацией измерительного устройства / А.Л. Шестаков, С.А. Загребина,
H.А. Манакова и др. // Автоматика и телемеханика. - 2021. - № 1. - С. 55-67.
25. Meyer-Baese, U. Digital Signal Processing with Field Programmable Gate Arrays / U. Meyer-Baese. - Springer-Verlag Berlin, Heidelbrg, 2004. - 527 p.
26. Yeh, Y.-C. Improving the Sampling Resolution of Periodic Signals by using Controlled Sampling Interval Method / Y.-C. Yeh, Y. Chu, C.W. Chiou // Computers & Electrical Engineering. - 2014. - Vol. 40, no. 4. - P. 1064-1071.
27. Discrete Model of Paired Relay-Race / E.V. Larkin, A.V. Bogomolov, A.N. Privalov, N.N. Dobrovolsky // Вестник ЮУрГУ. Серия «Математическое моделирование и программирование». -2018. - Т. 11, № 3. - С. 72-84.
28. Ларкин, Е.В. Дискретный подход к моделированию синхронизированных эстафет / Е.В. Ларкин, А.Н. Привалов, А.В. Богомолов // Научно-техническая информация (НТИ). сер. 2. Информационные процессы и системы. - 2020. - № 2. - C. 17-26.
29. Relay Races along a Pair of Selectable Routes / E.V. Larkin, A.V. Bogomolov, A.N. Privalov, N.N. Dobrovolsky // Bulletin of the South Ural State University. Series: Mathematical Modelling, Programming and Computer Software. - 2018. - Vol. 11, no. 1. - P. 15-26.
30. Bielecki, T.R. Conditional Markov Chains: Properties, Construction and Structured Dependence / T.R. Bielecki, J. Jakubowski, M. Niew^glowski // Stochastic Processes and their Applications. - 2017. - Vol. 127, no. 4. - P. 1125-1170.
31. Markov Chains: Models, Algorithms and Applications / W.K. Ching, X. Huang, M.K. Ng, T.K. Siu // International Series in Operations Research & Management Science. - Vol. 189. - Springer Science + Business Media NY, 2013. - 243 p.
32. Janssen, J. Applied Semi-Markov Processes / J. Janssen, R. Manca. - Springer US, 2006. -310 p.
33. Ang, K.H. PID Control System Analysis, Design and Technology / K.H. Ang, G. Chong, Y. Li // IEEE Transactions of control systems technology. - 2005. Vol. 13, no. 4. - P. 559-576.
34. O'Dwier, A. PID Compensation of Time Delay Processes 1999-2002: a Survey / A. O'Dwier // Proc. 2003 American Control Conference. USA, Denver, Colorado. - 2003. - P. 1494-1499.
35. Pavlov, A.V. About the Equality of the Transform of Laplace to the Transform of Fourier / A.V. Pavlov // Пробл. анал. Issues Anal. - 2016. - Т. 5(23), № 1. - P. 21-30.
36. Li, J. Three Effective Inverse Laplace Transform Algorithms for Computing Time-Domain Electromagnetic Responses / J. Li, C.G. Farquharson, X. Hu // Geophysics. - 2015. - Vol. 81, Iss. 2. -P. 75-90.
Поступила в редакцию 3 сентября 2024 г.
Сведения об авторах
Ларкин Евгений Васильевич - доктор технических наук, профессор, профессор-консультант, Тульский государственный университет, г. Тула, Российская Федерация, e-mail: [email protected].
Богомолов Алексей Валерьевич - доктор технических наук, профессор, ведущий научный сотрудник Федерального исследовательского центра «Информатика и управление» РАН, г. Москва, Российская Федерация, e-mail: [email protected].
Привалов Александр Николаевич - доктор технических наук, профессор, директор института передовых информационных технологий, Тульский государственный педагогический университет им. Л.Н. Толстого, г. Тула, Российская Федерация, e-mail: [email protected].
Акименко Татьяна Алексеевна - кандидат технических наук, доцент, Тульский государственный университет, г. Тула, Российская Федерация, e-mail: [email protected].
Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics" _2024, vol. 16, no. 4, pp. 43-55
DOI: 10.14529/mmph240406
MATHEMATICAL SUPPORT FOR HIERARCHICAL DIGITAL CONTROL OF A COMPLEX TECHNOLOGICAL OBJECT
E.V. Larkin1, A.V. Bogomolov2, A.N. Privalov3, T.A. Akimenko1
1 Tula State University, Tula, Russian Federation
2 Federal Research Center "Informatics and Control" of the Russian Academy of Sciences, Moscow, Russian Federation
Tula State Pedagogical University named after L.N. Tolstoy, Tula, Russian Federation E-mail: [email protected]
Abstract. The paper investigates the synthesis of mathematical support for hierarchical digital control of a complex technological object based on universal controllers with flexible programming. In order to overcome the limitations on the performance of controllers operating in real time, with the complexity of control laws and the expansion of the classes of controlled objects, controllers can be divided into strategic, tactical and functional-logical hierarchical levels. The authors define the functions and construct the algorithms for processing data by the controllers of each hierarchical level. The coordinated functioning of the digital control system under study is reduced to the exchange of data between controllers of different levels, carried out programmatically, equally for controllers of all levels. The authors obtain the dependencies for calculating the estimate of the restart period of control programs operating on the Hamiltonian cycle, and the data transmission delay for the case when the proportional-integral-differentiative regulator algorithm is implemented on the controllers of the tactical and functional-logical levels. The paper highlights the features of controlling the operating modes of the controllers of the tactical and functional-logical levels due to the implementation of a neural network at the strategic level.
Keywords: digital control; control system; hierarchical control levels; controller with flexible programming; Hamiltonian cycle; semi-Markov process; neural network.
References
1. Lofving M., Safsten K.; Winroth M. Manufacturing Strategy Formulation, Leadership Style and Organizational Culture in Small and Medium-Sized Enterprises. International Journal of Manufacturing Technology and Management, 2016, Vol. 30, no. 5, pp. 306-325. DOI: 10.1504/IJMTM.2016.078918
2. Pan Y., Guo K., Bobtsov A., Yang Ch., Yu H. Composite Error Learning Robot Control using Discontinuous Lyapunov Analysis. IEEE Transactions on Automatic Control. 2024. Vol. 69, Iss. 3, pp. 1705-1712. DOI: 10.1109/tac.2023.3326749
3. Le V.N., Ronzhin A.L. Review of Intelligent Control Systems and Robotic Tasks in Aquaculture Production. Marine intelligent technologies, 2024, no. 1-1 (63), pp. 171-180. DOI: 10.37220/MIT.2024.63.1.021 (in Russ.).
4. Meshcheryakov R. Control of Hyperlinked Cyber-Physical Systems. In: Ronzhin A., Shishlakov V. (eds) Proc. 14th International Conference on Electromechanics and Robotics "Zavalishin's Readings". Smart Innovation, Systems and Technologies, 2020, Vol. 154. Springer, Singapore. DOI: 10.1007/978-981-13-9267-2_3
5. Galin R.R., Shirokiy A.A., Magid E.A., Meshcheryakov R.V., Mamchenko M.V. Effective Functioning of a Mixed Heterogeneous Team in a Collaborative Robotic System. Computer Science and Automation, 2021, Vol. 20, no. 6, pp. 1224-1253. DOI: 10.15622/ia.20.6.2 (in Russ.).
6. Meshcheryakov R.V. Intelligent Robotic Systems. Sensors and Systems, 2024, no. 2, pp. 12-17. (in Russ.). DOI: 10.25728/datsys.2024.2.2
7. Meshcheryakov R.V., Shirokov A.S. Statement of the problem of heterogeneous group interaction of robots in solving smart garden problems. Bulletin of the South Ural State University. Series
"Mathematics. Mechanics. Physics", 2024, Vol. 16, no. 2, pp. 41-49. (in Russ.). DOI: 10.14529/mmph240204
8. Osipov O.Yu., Meshcheryakov R.V., Shepelenko M.G. Design of Digital Models of Elements of the Electric Machine Part of Electromechatronic Modules of Robotic Systems. Extreme Robotics, 2017, no. 1, pp. 160-164. (in Russ.).
9. Meshcheryakov R., Iskhakov A., Mamchenko M., Romanova M., Uvaysov S., Amirgaliyev Y., Gromaszek K. A Probabilistic Approach to Estimating Allowed SNR Values for Automotive LiDARs in "Smart Cities" under Various External Influences. Sensors, 2022, Vol. 22, no. 2, pp. 609. (in Russ.). DOI: 10.3390/s22020609
10. Smolin A.L., Kutakhov V.P., Zakharova A.A., Meshcheryakov R.V., Klimakov V.S. Formulation of the Problem of Implementing an Unmanned Aircraft Transport System Model using Artificial Intelligence Technologies. Artificial Intelligence. Theory and Practice, 2023, no. 1(1), pp. 41-45. (in Russ).
11. Landau I.D., Zito G. Digital Control Systems, Design, Identification and Implementation. Springer, 2006, 484 p. DOI: 10.1007/978-1-84628-056-6
12. Astrom J., Wittenmark B. Computer Controlled Systems: Theory and Design. Tsinghua University Press. Prentice Hall, 1996, 557 p.
13. Ivutin A., Larkin E. Estimation of Latency in Embedded Real-Time Systems, 2014 3rd Mediterranean Conference on Embedded Computing (MECO), Budva, Montenegro, 2014, pp. 236-239. DOI: 10.1109/MECO.2014.6862704.
14. Kilian C.T. Modern Control Technology: Components and Systems. Thompson Delmar Learning, 2000, 644 p.
15. Larkin E.V., Akimenko T.A., Bogomolov A.V. Modeling the Reliability of the Onboard Equipment of a Mobile Robot. Izvestiya of Saratov University. Mathematics. Mechanics. Informatics, 2021, Vol. 21, no. 3, pp. 390-399. DOI: 10.18500/1816-9791-2021-21-3-390-399
16. Babishin V., Taghipour S. Optimal Maintenance Policy for Multicomponent Systems with Periodic and Opportunistic Inspections and Preventive Replacements. Applied Mathematical Modelling, 2016, Vol. 40, Iss. 23-24, pp. 10480-10505. (in Russ.). DOI: 10.1016/j.apm.2016.07.019
17. Larkin E., Privalov A., Bogomolov A., Akimenko T. Digital Control of Continuous Production with Dry Friction at Actuators. In: Ronzhin, A., Shishlakov, V. (eds) Electromechanics and Robotics. Smart Innovation, Systems and Technologies, 2022, Vol. 232. Springer, Singapore. DOI: 10.1007/978-981-16-2814-6_37
18. Soldatov E.S., Bogomolov A.V., Larkin E.V., Soldatov A.S. Sensor Networks for Monitoring the Condition of Aviation Equipment during Testing and Operation. Aerospace Instrumentation, 2024, no. 2, pp. 61-68. (in Russ.). DOI: 10.25791/aviakosmos.2.2024.1395
19. Larkin E.V., Bogomolov A.V., Privalov A.N., Akimenko T.A. Structure of a Hierarchical Digital Control System for a Complex Technological Object. Mathematical Methods in Technology and Engineering, 2024, no. 6, pp. 63-66. (in Russ.).
20. Shestakov A.L., Zamyshlyaeva A.A., Manakova N.A., Sviridyuk G.A., Keller A.V. Reconstruction of a Dynamically Distorted Signal Based on the Theory of Optimal Dynamic Measurements. Automation and Remote Control, 2021, Vol. 82, no. 12, pp. 2143-2154. DOI: 10.1134/S0005117921120067
21. Ermolov I.L., Graskin S.S., Khripunov S.P. Models for Decision Making Support Systems in Robotics. In: Gorodetskiy, A.E., Tarasova, I.L. (eds) Smart Electromechanical Systems. Studies in Systems, Decision and Control, 2021, Vol. 352. Springer, Cham. DOI: 10.1007/978-3-030-68172-2_12
22. Bychkov E.V., Zagrebina S.A., Zamyshlyaeva A.A., Keller A.V., Manakova N.A., Sagadeeva M.A., Sviridyuk G.A. Development of the Theory of Optimal Dynamic Measurements. Bulletin of the South Ural State University. Series: Mathematical Modeling and Programming, 2022, Vol. 15, no. 3, pp. 19-33. (in Russ.). DOI: 10.14529/mmp220302
23. Bazhenov E.I., Mokrushin S.A., Okhapkin S.I. Software for the Mobile Robot Spatial Orientation System. Bulletin of the South Ural State University. Series: Mathematical Modelling, Programming and Computer Software, 2021, Vol. 14, no. 2, pp. 64-69. DOI: 10.14529/mmp210206
24. Shestakov A.L., Zagrebina S.A., Manakova N.A., Sagadeeva M.A., Sviridyuk G.A. Numerical Optimal Measurement Algorithm under Distortions Caused by Inertia, Resonances, and Sensor Degra-
dation. Automation and Remote Control, 2021, Vol. 82, no. 1, pp. 41-50. DOI: 10.1134/S0005117921010021
25. Meyer-Baese U. Digital Signal Processing with Field Programmable Gate Arrays. SpringerVerlag Berlin, Heidelbrg, 2004, 527 p. DOI: 10.1007/978-3-662-06728-4
26. Yeh Y.-C., Chu Y., Chiou C.W. Improving the Sampling Resolution of Periodic Signals by using Controlled Sampling Interval Method. Computers & Electrical Engineering, 2014. Vol. 40, no. 4, pp. 1064-1071. DOI: 10.1016/j.compeleceng.2013.12.002
27. Larkin E.V., Bogomolov A.V., Privalov A.N., Dobrovolsky N.N. Discrete Model of Paired Relay-Race. Bulletin of the South Ural State University. Series: Mathematical Modelling, Programming and Computer Software. 2018. Vol. 11, no. 3, pp. 72-84. (in Russ.). DOI: 10.14529/mmp180306
28. Larkin E.V., Privalov A.N., Bogomolov A.V. Diskretnyy podkhod k modelirovaniyu sinkhronizirovannykh estafet (Discrete Approach to Modeling Synchronized Relay Races). Nauchno-tekhnicheskaya informatsiya (NTI). ser. 2. Informatsionnye protsessy i sistemy (Scientific and Technical Information. Series 2: Information Processes and Systems), 2020, no. 2, pp. 17-26. (in Russ.). DOI: 10.36535/0548-0027-2020-02-3
29. Larkin E.V., Bogomolov A.V., Privalov A.N., Dobrovolsky N.N. Relay Races along a Pair of Selectable Routes. Bulletin of the South Ural State University. Series: Mathematical Modelling, Programming and Computer Software, 2018, Vol. 11, no. 1, pp. 15-26. (in Russ.). DOI: 10.14529/mmp180102
30. Bielecki T.R., Jakubowski J., Niewçglowski M. Conditional Markov chains: Properties, Construction and Structured Dependence. Stochastic Processes and their Applications, 2017, Vol. 127, no. 4, pp. 1125-1170. DOI: 10.1016/j.spa.2016.07.010
31. Ching W.K., Huang X., Ng M.K., Siu T.K. Markov Chains: Models, Algorithms and Applications. International Series in Operations Research & Management Science, Vol. 189, Springer Science + Business Media NY, 2013, 243 p. DOI: 10.1007/978-1-4614-6312-2
32. Janssen J., Manca R. Applied Semi-Markov Processes. Springer US, 2006, 310 p. DOI: 10.1007/0-387-29548-8
33. Ang K.H., Chong G., Li Y. PID Control System Analysis, Design and Technology. IEEE Transactions of control systems technology, 2005, Vol. 13, no. 4, pp. 559-576. DOI: 10.1109/TCST.2005.847331
34. O'Dwyer A. PID Compensation of Time Delayed Processes 1998-2002: a Survey. Proc. 2003 American Control Conference, 2003, Denver, CO, USA, 2003, pp. 1494-1499. DOI: 10.1109/ACC.2003.1239802.
35. Pavlov A.V. About the Equality of the Transform of Laplace to the Transform of Fourier. Probl. Anal. Issues Anal., 2016, Vol. 5(23), no. 1, pp. 21-30. (in Russ.). DOI: 10.15393/j3.art.2016.3110
36. Li J., Farquharson C.G., Hu X. Three Effective Inverse Laplace Transform Algorithms for Computing Time-Domain Electromagnetic Responses. Geophysics, 2015, Vol. 81, no. 2, pp. 75-90. (in Russ.). DOI: 10.1190/geo2015-0174.1
Received September 3, 2024
Information about the authors
Larkin Evgeny Vasilyevich is Dr. Sc. (Engineering), Professor, Professor-Consultant of Tula State University, Tula, Russian Federation, e-mail: [email protected].
Bogomolov Aleksey Valerievich is Dr. Sc. (Engineering), Professor, Leading Researcher of thebFederal Research Center «Computer Science and Control» of the Russian Academy of Sciences, Moscow, Russian Federation, e-mail: [email protected].
Privalov Aleksandr Nikolaevich is Dr. Sc. (Engineering), Professor, Director of the Institute of Advanced Information Technologies of Tolstoy Tula State Pedagogical University, Tula, Russian Federation, e-mail: [email protected].
Akimenko Tatyana Alekseevna is Cand. Sc. (Engineering), Associate Professor, Associate Professor of Tula State University, Tula, Russian Federation, e-mail: [email protected].