Научная статья на тему 'ПОСТРОЕНИЕ КВАЗИОПТИМАЛЬНОЙ ПО БЫСТРОДЕЙСТВИЮ И ЭНЕРГОЗАТРАТАМ ЗАМКНУТОЙ СИСТЕМЫ УПРАВЛЕНИЯ ЭЛЕКТРОУСТАНОВКОЙ'

ПОСТРОЕНИЕ КВАЗИОПТИМАЛЬНОЙ ПО БЫСТРОДЕЙСТВИЮ И ЭНЕРГОЗАТРАТАМ ЗАМКНУТОЙ СИСТЕМЫ УПРАВЛЕНИЯ ЭЛЕКТРОУСТАНОВКОЙ Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
84
8
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЭЛЕКТРОУСТАНОВКА / ОПТИМАЛЬНОЕ УПРАВЛЕНИЕ / БЫСТРОДЕЙСТВИЕ / ЭНЕРГОЗАТРАТЫ / ПРОГРАММНОЕ ДВИЖЕНИЕ / ПРИНЦИП МАКСИМУМА / ОСОБОЕ УПРАВЛЕНИЕ / УСЛОВИЯ ОБЩНОСТИ ПОЛОЖЕНИЯ / УСТОЙЧИВОСТЬ

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Хорошавин Валерий Степанович, Грудинин Виктор Степанович

Введение. Эффективность и работоспособность электротехнологических процессов и установок в динамических режимах их функционирования характеризуются быстродействием, низкими энергозатратами, точностью в переходных процессах и устойчивостью желаемого конечного состояния. Предлагается на базе особого оптимального управления совместить построение систем перехода и стабилизации как в части определения алгоритма управления с минимальными энергозатратами в функции состояний объекта, так и в части определения параметров и условий движения с минимальным отклонением от заданной траектории. Последнее обеспечивает в устойчивой замкнутой системе программного движения эффективные решения по критериям быстродействия и энергосбережения. Материалы и методы. В качестве основного метода нахождения оптимального программного управления используется принцип максимума, дополненный для исследования особых режимов аппаратом условий общности положения для нелинейных объектов с расширением пространства координат, учитывающий вхождение времени и критерий оптимальности. Аппарат условий общности положения использован и для решения задач энергосбережения путем линеаризации в большом исходного объекта. Квазиоптимальность по быстродействию и энергозатратам достигается минимизацией энергии по параметру программного движения, противоречиво влияющему на время переходного процесса и амплитуду управления. Результаты исследования. Для оценки вычислительных трудностей, быстродействия, энергосбережения, точности и устойчивости приведен пример управления инерционным объектом по различным критериям. Получена структура простой в технической реализации замкнутой квазиоптимальной системы со стационарной обратной связью. Обсуждение и заключение. Формализованность подхода к построению квазиоптимальных систем на основе условий общности положения позволяет использовать его в задачах многокритериальной оптимизации и системах автоматизированного проектирования энергоемких промышленных, транспортных, сельскохозяйственных электроустановок.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Хорошавин Валерий Степанович, Грудинин Виктор Степанович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

DEVELOPING A QUASI-OPTIMAL, IN TERMS OF TRANSITION TIME AND ENERGY CONSUMPTION, CLOSED-LOOP CONTROL SYSTEM FOR AN ELECTRICAL INSTALLATION

Introduction. The efficiency and normal work of electrotechnological processes and installations in their operation dynamic modes is characterized by the time of transition from initial to final state, low energy consumption, accuracy in transients and stability of the desired final state. It is proposed, from a single position on the basis of special optimal control, to combine transition and stabilization systems both in terms of determining the control algorithm with minimal energy consumption in the function of the object states and determining the parameters and conditions of movement with minimal deviation from a given trajectory providing the optimal transmission time and energy saving in a stable closed system of an object control. Materials and Methods. The principle of maximum is used as the main method for finding optimal program control, which for the study of special modes was supplemented with the apparatus of the position generality conditions for nonlinear objects with the coordinate space expansion, taking into account the occurrence of time and optimality criterion. The position generality apparatus is also used to solve energy-saving problems through using linearization in a large source object. Quasi-optimality in terms of transition time and energy consumption is achieved through minimizing energy according to the program motion parameter, which has a contradictory effect on the transition time and control amplitude. Results. To assess computational difficulties, transition time, energy saving, accuracy and stability, an example of inertial object control according to various criteria is given. The structure of a closed quasi-optimal system with stationary feedback, which is simple in technical implementation, is obtained. Discussion and Conclusion. The formalization of the approach to the construction of quasi-optimal systems based on the position generality allows it to be used in multi-criteria optimization tasks and computer-aided design systems for energy-intensive industrial, transport, and agricultural electrical installations.

Текст научной работы на тему «ПОСТРОЕНИЕ КВАЗИОПТИМАЛЬНОЙ ПО БЫСТРОДЕЙСТВИЮ И ЭНЕРГОЗАТРАТАМ ЗАМКНУТОЙ СИСТЕМЫ УПРАВЛЕНИЯ ЭЛЕКТРОУСТАНОВКОЙ»



http://vestnik.mrsu.ru

ISSN Print 2658-4123 ISSN Online 2658-6525

ЭЛЕКТРОТЕХНОЛОГИИ И ЭЛЕКТРООБОРУДОВАНИЕ / ELECTRICAL TECHNOLOGIES AND EQUIPMENT

УДК 621.31

doi: 10.15507/2658-4123.032.202202.279-294

Научная статья

Ш Построение квазиоптимальной

по быстродействию и энергозатратам замкнутой системы управления электроустановкой

В. С. Хорошавин н, В. С. Грудинин

Вятский государственный университет (г. Киров, Российская Федерация)

н khoroshavin@vyatsu.ru

Аннотация

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

Материалы и методы. В качестве основного метода нахождения оптимального программного управления используется принцип максимума, дополненный для исследования особых режимов аппаратом условий общности положения для нелинейных объектов с расширением пространства координат, учитывающий вхождение времени и критерий оптимальности. Аппарат условий общности положения использован и для решения задач энергосбережения путем линеаризации в большом исходного объекта. Квазиоптимальность по быстродействию и энергозатратам достигается минимизацией энергии по параметру программного движения, противоречиво влияющему на время переходного процесса и амплитуду управления. Результаты исследования. Для оценки вычислительных трудностей, быстродействия, энергосбережения, точности и устойчивости приведен пример управления инерционным объектом по различным критериям. Получена структура простой в технической реализации замкнутой квазиоптимальной системы со стационарной обратной связью. Обсуждение и заключение. Формализованное^ подхода к построению квазиоптимальных систем на основе условий общности положения позволяет использовать

© Хорошавин В. С., Грудинин В. С., 2022

© I Контент доступен по лицензии Creative Commons Attribution 4.0 License. ■З^Н This work is licensed under a Creative Commons Attribution 4.0 License.

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

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

Конфликт интересов: авторы заявляют об отсутствии конфликта интересов.

Для цитирования: Хорошавин В. С., Грудинин В. С. Построение квазиоптимальной по быстродействию и энергозатратам замкнутой системы управления электроустановкой // Инженерные технологии и системы. 2022. Т. 32, № 2. С. 279-294. (М: https://doi.org/10.15507/2658-4123.032.202202.279-294

Original article

Developing a Quasi-Optimal, in Terms of Transition Time and Energy Consumption, Closed-Loop Control System for an Electrical Installation

V. S. Khoroshavin V. S. Grudinin

Vyatka State University (Kirov, Russian Federation)

H khoroshavin@vyatsu.ru Abstract

Introduction. The efficiency and normal work of electrotechnological processes and installations in their operation dynamic modes is characterized by the time of transition from initial to final state, low energy consumption, accuracy in transients and stability of the desired final state. It is proposed, from a single position on the basis of special optimal control, to combine transition and stabilization systems both in terms of determining the control algorithm with minimal energy consumption in the function of the object states and determining the parameters and conditions of movement with minimal deviation from a given trajectory providing the optimal transmission time and energy saving in a stable closed system of an object control.

Materials and Methods. The principle of maximum is used as the main method for finding optimal program control, which for the study of special modes was supplemented with the apparatus of the position generality conditions for nonlinear objects with the coordinate space expansion, taking into account the occurrence of time and optimality criterion. The position generality apparatus is also used to solve energy-saving problems through using linearization in a large source object. Quasi-optimality in terms of transition time and energy consumption is achieved through minimizing energy according to the program motion parameter, which has a contradictory effect on the transition time and control amplitude. Results. To assess computational difficulties, transition time, energy saving, accuracy and stability, an example of inertial object control according to various criteria is given. The structure of a closed quasi-optimal system with stationary feedback, which is simple in technical implementation, is obtained.

Discussion and Conclusion. The formalization of the approach to the construction of quasi-optimal systems based on the position generality allows it to be used in multi-criteria optimization tasks and computer-aided design systems for energy-intensive industrial, transport, and agricultural electrical installations.

Keywords: electrical installation, optimal control, time transition, energy consumption, program motion, maximum principle, special control, position generality conditions, stability

Conflict of interest: The authors declare no conflict of interest.

For citation: Khoroshavin VS., Grudinin VS. Developing a Quasi-Optimal, in Terms of Transition Time and Energy Consumption, Closed-Loop Control System for an Electrical Installation. Engineering Technologies and Systems. 2022;32(2):279-294. doi: https://doi. org/10.15507/2658-4123.032.202202.279-294

Введение

Для электротехнических процессов и установок актуальна задача повышения эффективности работы, так как 70 % вырабатываемой электроэнергии потребляют электроприводы, 10 % - электротермические установки1. Установки инфракрасного излучения, светотехнические и насосные установки являются энергоемкими в современных технологиях выращивания сельскохозяйственных культур в защищенном грунте [1]. Многие процессы и установки, особенно тепловые, имеют большую инерционность, то есть длительные динамические режимы перехода в устойчивое требуемое выходное состояние. Динамические режимы электроустановок как объектов управления с достаточной точностью описываются обыкновенными дифференциальными уравнениями, а в качестве управляющих воздействий служат ограниченные электрические переменные: обычно амплитуда, частота или скважность напряжения или тока.

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

60 % потребляемой электроэнергии [2]. Но известные алгоритмы оптимального быстродействия и энергосбережения отыскиваются в виде программного (то есть в функции времени) управления, не обеспечивают устойчивости конечного состояния объекта и в лучшем случае могут быть реализованы в форме нестационарной обратной связи.

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

Обзор литературы

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

1 Браславский И. Я., Ишматов Э. Ш., Поляков В. Н. Энергосберегающий асинхронный электропривод : учеб. пособие для вузов. М. : Издательский центр «Академия», 2004. 256 с. URL: http:// log-in.ru/books/energosberegayushiiy-asinkhronnyiy-elektroprivod-braslavskiiy-i-ya-ishmatov-z-sh-polyakov-v-n-tekhnika/ (дата обращения: 20.01.2022) ; Автоматическое управление электротермическими установками : учебник для вузов / А. Д. Свенчанский [и др.]. М. : Энергоатомиздат, 1990. 416 с. URL: https://search.rsl.ru/ru/record/01001521815 (дата обращения: 20.01.2022).

2 Математическая теория оптимальных процессов / Л. С. Понтрягин [и др.]. М. : Наука, 1969. 384 с. URL: https://search.rsl.ru/ru/record/01006401476 (дата обращения: 20.01.2022).

режима, когда принцип максимума не устанавливает однозначной связи между оптимальным управлением и вспомогательными переменными3. Известные способы вычисления особого управления используют вспомогательные переменные непосредственно или с помощью скобок Пуассона [3; 4]. Найти решение для вспомогательных переменных, даже в функции времени, сложно, а для построения замкнутой системы в функции координат, позволяющей компенсировать действующие на объект возмущения, невозможно из-за неразрешимости двухточечной граничной задачи, к которой сводится решение оптимальной задачи. В задачах энергосбережения с квадратичным по управлению критерием особых режимов не возникает, но даже если оптимальное управление в функции времени будет получено, то синтез замкнутой системы осуществляется с помощью зависящих от времени и граничных условий обратных связей [2]. Тогда как для простоты технической реализации системы требуются стационарные обратные связи.

Для достижения предельного быстродействия или минимума ресурсов в задачах с линейным управлением необходимо кусочно-постоянное управление с максимально возможными амплитудами воздействия, что ведет к перерегулированию координат, а моменты переключения определяются нелинейными нестационарными поверхностями переключения или стыковки особых и неособых траекторий [5; 6]. Поэтому алгоритмы оптимального быстродействия неприменимы в задачах программного движения

или оптимального по точности управления по достижению минимального отклонения от заданной траектории, особенно для объектов с длительными и частыми динамическими режимами перехода [7]. К тому же при оптимальном быстродействии не обеспечивается устойчивость конечного состояния в отличие от программного движения, гарантирующего ее на полубесконечном интервале времени.

Что касается связи решений задач управления с критериями быстродействия и энергосбережения, то известны противоположные точки зрения: от их полного противопоставления, что объясняется необходимостью больших энергозатрат для увеличения быстродействия, до их полного совпадения, что справедливо при управлении в большом, когда при больших диапазонах задания начального и конечного состояний выхода минимальные энергозатраты получаются при максимально возможном управлении, то есть как в задаче быстродействия4 [8]. Сравнительные анализы оценки критериев быстродействия, энергозатрат и точности (программного движения) выполнены при управлении реальными тепловыми процессами в металлургии и теплоснабжении [9-11]. На основе анализа особых, в смысле принципа максимума, режимов в задачах с линейным вхождением управления показано, что алгоритмы оптимальных управлений по быстродействию и на минимум ресурсов совпадают, если критерий на минимум ресурсов физически адекватно отражает поступление энергии в управляемый процесс [12].

3 Габасов Р., Кириллова Ф. М. Особые оптимальные управления. М. : Наука, 1973. 314 с. URL: https://litgu.ru/knigi/tehnicheskie_nauki/308008-osobye-optimalnye-upravleniya.html (дата обращения: 20.01.2022).

4 Макеева А. В. Оптимальное по быстродействию управление одной линейной системой второго порядка // Сборник докладов всерос. науч.-практ. конф. с междунар. участием «Фундаментальные и прикладные проблемы механики, математики, информатики». Пермь : ФГБОУ ВО «Пермский государственный национальный исследовательский университет», 2015. С. 68-71. URL: https://elibrary.ru/item.asp?id=24231084&pff=1 (дата обращения: 24.01.2022).

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

Материалы и методы Пусть математическое описание объекта управления (электроустановки) задается системой обыкновенных дифференциальных уравнений, нелинейных по координатам с линейным скалярным управлением в векторно-матричной форме:

X = A (x ) + B (x )U,

(1)

где х е Rn — вектор координат объекта; и - управление, 0 < и < 1; элементы вектор-столбцов А(х), В(х) непрерывны и дифференцируемы по х.

Задачами системы управления являются, во-первых, выработка ограниченного управления в функции времени или координат, переводящего объект из начального состояния х(0) в конечное х(Т) за незаданное или заранее заданное время и обеспечивающего минимум критерию оптимальности в виде интегрального функционала, учитывающего временные или энергетические затраты или отклонение текущей траектории от заданной (задача оптимального управления); во-вторых, последующая стабилизация конечного состояния в замкнутой системе со стационарными обратными связями (задача устойчивости или, в конечном счете, обеспечения работоспособности установки).

Для решения первой задачи оптимального управления с определением

существования особого режима и вычисления особого управления в явном виде от координат и параметров нелинейного объекта, что необходимо на практике для синтеза замкнутой системы с обратной связью, применим аппарат условий общности положения (УОП) для нелинейных объектов. УОП был выведен в работе В. А. Олейникова для задач быстродействия и расширен для задач на минимум ресурсов в одном из наших исследований5.

Покажем применение УОП для нелинейных объектов в расширенном пространстве координат в задаче программного движения с критерием

7=ИЫо-х»с))2*> (2)

учитывающим квадратичные отклонения текущей выходной координаты объекта от требуемой [13]. Задача (1), (2) является нелинейной по координатам, поэтому в ней при использовании принципа максимума, даже для линейных объектов, возможно возникновение особых ситуаций, если на некотором интервале времени I е [¿2]:

dk (у,B (X))

dtk

= 0, k = 0,1,2,

(3)

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

dy dt

дЛ (x) dB (x )U

\T

die

die

y.

У „„ „„ у

Условие существования особого режима (3) требует анализа вспомогательных переменных Чтобы его избежать, в аппарате УОП вычисляются

5 Олейников В. А. Оптимальное управление технологическими процессами в нефтяной и газовой промышленности. Л. : Недра, 1982. 216 с. URL: https://search.rsl.ru/ru/record/01001089712 (дата обращения: 24.01.2022) ; Хорошавин В. С., Зотов А. В. Особое оптимальное управление нелинейными объектами : моногр. Киров : Науч. изд-во ВятГУ, 2019. 219 с. URL: https://www.elibrary.ru/ item.asp?id=37523319 (дата обращения: 24.01.2022).

векторы Bj, j = 2,..., n + 2, по рекуррентному соотношению:

B (x ) = в (x),

8BJ_1 dU

Bj =

8A (X ) + 8B (X )U

8U dt

\

8je

8je

BJ i

8Bj_i dx

8x dt

, (4)

где x =

B (x ) =

A (x ) =

v x«+i J

()- (t)) A (x) 1

2^

B (*) 0

D =

Vb3i

u Dn+i =

u32j

yki

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

Интегральный критерий энергозатрат для большинства электротехнических устройств имеет вид

I =

f TU 2dt. J 0

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(5)

После образования матрицы Dn+2 размером (п + 2) • (п + 2) из векторов (В1 ... Ви+2) и приравнивания к нулю определителя det Dn+2 = 0 находятся уравнения особых траекторий и особых управлений, множество которых можно расширить путем приравнивания к нулю функциональных элементов матрицы Оп+2, так как при этом тождественно выполняются условия особого режима (3).

При исследовании функциональных элементов матрицы Dn+2 учитываются не только особые режимы задачи программного движения, но и задач быстродействия и на минимум ресурсов с линейным управлением, так как в матрицу Оп+2 вложены матрицы

1333 У

Как и М. Атанс и П. Фалб, предположим, что «и(0 - скаляр, пропорциональный напряжению или току, тогда величина и2(0 пропорциональна мощ-

¡•Т 2

ности, а I и dt пропорционален энер-

■> to

гии, израсходованной на интервале [¿0, Т ]»6. Такой критерий энергозатрат используется в большинстве работ, но иногда, как заметил Ю. П. Петров, «в задачах частотного управления рекомендуется использовать критерий !т 15

I и ' dt, так как потери на гистерезис

■> to

и вихревые токи пропорциональны соответственно первой и второй степени частоты электрического сигнала»7.

Задача энергосбережения (1), (5) нелинейна по управлению, поэтому особых режимов в ней не возникает, и при ее решении можно однозначно выразить оптимальное программное управление, как отмечает В. И. Ловчаков, через вспомогательные переменные с помощью принципа максимума или множители Лагранжа в классическом вариационном исчислении или методом динамического программирования [2].

6 Атанс М., Фалб П. Оптимальное управление. М. : Машиностроение, 1968. 734 с. URL: https:// search.rsl.ru/ru/record/01005951981 (дата обращения: 24.01.2022).

7 Петров Ю. П. Оптимальное управление электроприводом с учетом ограничений по нагреву. Л. : Энергия, 1971. 144 с. URL: https://search.rsl.ru/ru/record/01007323479 (дата обращения: 24.01.2022).

x

о

x

ENGINEERING TECHNOLOGIES AND SYSTEMS

Vol. 32, no. 2. 2022

Чтобы исключить трудоемкую операцию нахождения дополнительных переменных, предложен способ поиска дифференциального уравнения для оптимального управления в функции координат объекта, основанный на применении аппарата УОП для нелинейных объектов, для чего производится приведение исходной задачи (1), (5) к редуцированной задаче с линейным вхождением управления, в которой возможно существование особого режима. Как мы писали в одной из предыдущих работ, «УОП для нелинейных объектов в расширенном пространстве координат Rn+2 применяются к редуцированной задаче энергосбережения (1), (5):

< х = А (х) + В (х) хй, Х«+1 =(х° )2.

Из выражения для детерминанта и элементов матрицы 0„+2 = (В1 В2 ... Вп Вп+1 Вп+2), вычисленных с учетом замены и на у по рекуррентному соотношению (4), в общем случае получается множество особых управлений у и его производных по времени в функции координат редуцированного объекта

¿й Dn+1 = F0 (х, х0, у, у,...) = 0,

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

X0 = и, —рт у = —ти, к = 1,2,3,... —к —к

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

управления в функции координат исходного объекта»8. А именно:

F [х,и ,иД...) = 0.

Постоянные интегрирования в общем решении системы дифференциальных уравнений объекта (1) и дифференциального уравнения для оптимального программного управления (6) определяются из решения двухточечной граничной задачи для заданных состояний выхода объекта и его (п - 1) производных в начальный и конечный фиксированный момент времени Т. Получаемое управление и соответствующие фазовые траектории являются трансцендентными функциями от времени, и замкнутая система может быть реализована с помощью нелинейных функциональных преобразователей, нестационарной обратной связи или путем аппроксимации траекторий и управления [2; 14-16].

Что касается второй задачи, стабилизации, то алгоритмы быстродействия и энергосбережения не обеспечивают устойчивость конечного состояния. Алгоритм программного движения может гарантировать устойчивость, поэтому для его эффективности по быстродействию и энергосбережению предлагается ввести в желаемый закон движения выходной координаты объекта такой параметр т, который противоречиво влияет на время переходного процесса Т и амплитуду управляющего воздействия и. Таким параметром, исходя из физических соображений по второму закону Ньютона, может быть мера инерционности процесса. Если взять апериодическое звено первого порядка с дифференциальной связью хх = и - х, то Т прямо пропорционально т и составляет (3^4)т, а амплитуда и обратно пропорциональна т. Тогда

8 Хорошавин В. С., Зотов А. В. Особое оптимальное управление нелинейными объектами. Electrical technologies and equipment 285

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

дЕ ((ч, xкон

,r,U ,T)

дт

= 0. (7)

x„

(t) = xH

+ x„

1 - e>

(8)

(9)

Если при вычисленном по (7) топт выполняются ограничения по управлению и координатам, то в задаче программного движения с критерием (2) при

где коэффициент усиления и постоянная времени равны 1, на управление и состояние наложены ограничения 0 < и < 1, 0 < х1 < 1. Граничные условия х1нач = 0, х1кон = хк. Выбор объекта объясняется меньшей громоздкостью вычислений и большей прозрачностью результатов для оценки вычислительных, временных и энергетических затрат и оценки устойчивости получаемых решений.

Задача оптимального энергосберегающего управления

УОП в расширенном пространстве координат R3 применяются к редуцированной задаче

обеспечивается квазиоптимальное (в смысле локальной оптимальности, в отличие от глобальной оптимальности при оптимальном программном управлении) управление по быстродействию и энергозатратам.

Устойчивость квазиоптимальной системы можно проверить по функции Ляпунова методами качественной теории дифференциальных уравнений или критериев устойчивости [17]. Структура замкнутой квазиоптимальной системы для линейных объектов управления содержит линейные стационарные обратные связи с коэффициентами передачи, зависящими от параметра топт.

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

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

х = .у,

Вычисляются векторы B, j = 1, 2, 3, по соотношению (4), и образуется матрица D3:

D3 = (B B2 B3 ) =

0 -1 -2 х0

0 -1 -2 7

Из det Б3 = 0 определяется особое управление в редуцированной задаче у - х0 = 0, после перехода к оптимальному управлению в исходной задаче решается система дифференциальных уравнений

\ ии — и.

Общим решением системы дифференциальных уравнений объекта и управления являются

9 Евланов Л. Г. Теория и практика принятия решений. М. : Экономика, 1984. 176 с. URL: https:// hram-sveta.ru/9189-teoriya-i-praktika-prinyatiya-resheniy.html (дата обращения: 24.01.2022).

\ х () = С1е' + С2е ', \ и (г) = 2С1е<.

Постоянные интегрирования определяются из граничных условий. Пусть х^0) = 0, х1кон = 0,5 при Т = 2. Графики переходных процессов для координаты х1*(^) = 0,069(е' + еи оптимального программного управления и *(0 = 0Д38е 'приведены на рисунке 1. При этом затраты энергии

Е = |2(и*^)}2А = 0,52, что на 30,7 %

меньше, чем энергозатраты в разомкнутой системе с управлением Ц (0 = 0,5 на интервале времени ' е [0, 3]. Но аналитически выразить Ц* через х* из-за трансцендентности невозможно, к тому же нужно будет дополнительно решать задачу стабилизации конечного состояния.

Задача квазиоптимального управления

Рассмотрим задачу программного движения для объекта (9) с критерием

1,2

1,0

_ T ( f t 2

1=f J 0 x K 1 - e 7 dt

V V У

учитывающим переход из х1(0) = 0 в х1(Т ^ да) = хк с минимальными квадратичными отклонениями от переходного процесса апериодического звена тХ — и - х с Ц = х .

к

Применим УОП в расширенном пространстве координат R3 к поставленной задаче

х0 = I

X1 = U - x1,

xo \

1 - e

-x

После вычисления векторов B, j = = 1, 2, 3, по соотношению (4) образуем матрицу D3 = (B1 B2 B3):

D =

0 ^

1

b33 j

0,8

0,6

0,4

0,2

U, x x1'{t), U*(t), x1(U = 0,5)

t

0,5

1,0

1,5 U*

2,0

2,5

3,0

-x1*(U*) -U* -x1(U)

Р и с. 1. Траектории x* и U* оптимального программного управления F i g. 1. Trajectories x1* and U of optimal programmed control Electrical technologies and equipment

2

0

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

0

из элементов которой

Ь32 = 2*к

1 - e

- 2xi,

Ъ33 = 2

X

v v

1 - e

+ 2 - 2U + 2xx

определяются в функции времени особая траектория, совпадающая с желаемой по критерию (10), и особое управление, являющееся оптимальным в задаче программного движения

U (t, хк,т) = хк

1 + е

. (11)

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

ихк,т) = Хк " Х(1 "Т). (12)

На рисунке 2 приведена графическая иллюстрация связей между переменными в уравнении (12). Из уравнений (11) и (12) следует, что начальная амплитуда управления и(х1(0) = 0) обратно пропорциональна параметру т. При т < тмин (рис. 2) начальный участок управления в программном движении, как в задаче быстродействия, принадлежит ограничению на управление и = 1 и чем больше х тем при больших тмин выходим на это ограничение, что следует и из связи

и (х1 (0) = 0) < 1 = . Это означает,

^мин

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

Величина энергии, затрачиваемой на переход из х1(0) = 0 в хк за время Т = = 3т при управлении (13), равна

x0 1

X

X

-Tau=oi -Tau=2 -Tau=1 -Tau=opt -Tau=min -Tau

Р и с. 2. Связи переменных в программном движении F i g. 2. Relationships of variables in the program motion

E = f > ( -к^))2

dt = x

(3т2 + 2т + 1)

Минимум энергии Е по параметру т определится из условия

8E 2

er-*

(2 - О

2r

= 0,

откуда топт = 0,58. Отметим, что полученный параметр топт не зависит от хк. Зависимость энергии Е от параметра т при х = 0,5 показана на рисунке 3. Спадающий в положительном направлении т участок функции Е(т) характеризует эффективные по быстродействию и энергозатратам режимы [18]. Оптимальное в программном движении решение с топт = 0,58 дает 9,3 % относительной экономии энергии, по сравнению с управлением в разомкнутой системе с т =1, что существенно ниже 30,7 % экономии энергии при оптимальном программном управлении, но при программном движении обеспечивается устойчивость конечного состояния.

При выборе хк > тмин, то есть, когда Тмин > Топт (рис. 2, 3), переходим на восходящую ветвь функции Е(т), быстродействие и энергозатраты непротиворечивы, их эффекты хотя и снижаются, но сохраняются по сравнению с т = 1. Отметим, что выбор меньшего т дает больший эффект по быстродействию или производительности установки, чем эффект по энергозатратам, что следует из слабой вогнутости кривой Е(т). При больших хк наблюдается более выраженный экстремум в зависимости Е(т). Графики переходных процессов в программном движении с топт = 0,58 и т = 1 для хк = 0,5 приведены на рисунке 4.

Структура позиционной системы управления, квазиоптимальной по быстродействию и энергозатратам, полученная из (12) с т и техниче ски про сто

опт

реализуемая, приведена на рисунке 5. Если граничные условия хк > топт или

тмин > топт, то хк > топт или тмин > топт, то для ограничения амплитуды управления введено звено насыщения с линейным участком единичного наклона в диапазоне [0, 1].

1,8 1,6 1,4 1,2 1,0 0,8 0,6 0,4 0,2 0

E

0,5

1,0

1,5

2,0

2,5

т

0

Р и с. 3. Зависимость энергии E от параметра т в программном движении F i g. 3. Dependence of the energy E on the parameter т in the program motion Electrical technologies and equipment

-^(TAUorn) -U(TAUonT) -*i(TAU=1)

Р и с. 4. Графики переходных процессов в программном движении F i g. 4. Graphs of transients in the program motion

Р и с. 5. Структура позиционной системы управления F i g. 5. Structure of the positional control system

Непосредственное применение результатов данной работы можно показать на примере оптимального энергосберегающего управления сушильной камерой, приведенного в работе В. И. Ловчакова [2]. Хотя поведение объекта описывается дифференциальным уравнением второго порядка, для конкретных граничных условий можно построить адекватную модель первого порядка с помощью метода наименьших квадратов по кривой переходного процесса или частотных характеристик. Далее, применяя пошагово методику данной работы, можно в замкнутой устойчивой системе получить квазиоптимальное по энергозатратам

и быстродействию управление, обеспечивающее 10-30 % экономии электроэнергии и в 1,5-2,0 раза большее быстродействие при корректном сравнении с управлением в разомкнутой системе, соответствующим конечному значению выходной координаты.

Обсуждение и заключение

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

Vol. 32, no. 2. 2022

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

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

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

Результаты работы могут быть использованы при исследовании динамических режимов электроустановок в промышленных и сельскохозяйственных тепловых процессах (жилые и нежилые помещения, теплицы, печи, сушильные камеры, автоклавы); в светотехнических установках; на транспорте в автономных электроприводах; в меха-тронных и робототехнических устройствах, например для обеспечения мягкого пуска электроприводов, и в других процессах и установках [19; 20].

СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ

1. Доманов В. И., Певчева Е. В. Анализ основных узлов энергосистемы тепличного комбината и способов снижения энергозатрат // Промышленные АСУ и контроллеры. 2017. № 3. С. 3-10. URL: https://www.elibrary.ru/item.asp?id=29115500 (дата обращения: 20.01.2022).

2. Ловчаков В. И. Аппроксимационный подход к синтезу систем регулирования на основе оптимального программного управления // Известия Тульского государственного университета. Технические науки. 2017. № 3. С. 225-236. URL: https://clck.ru/h7ejP (дата обращения: 20.01.2022).

3. Nikol'skii M. S. Singular Sets of Extremal Controls in Optimal Control Problems // Proceedings of the Steklov Institute of Mathematics. 2019. Vol. 304. P. 236-240. doi: https://doi.org/10.4213/tm3970

4. Цинцадзе З. Вычисление особого оптимального управления в квазилинейных управляемых системах со смешанными ограничениями // Компьютерные науки и телекоммуникации. 2005. № 2. С. 71-73. URL: http://gesj.internet-academy.org.ge/ru/list_aut_artic_ru.php?b_sec=&list_aut=1248 (дата обращения: 24.01.2022).

5. Gao Z. On Discrete Time Optimal Control: A Closed-Form Solution // Proceeding of the 2004 American Control Conference (30 June - 2 July 2004). Boston, 2004. P. 52-58. URL: https://folk.ntnu.no/ skoge/prost/proceedings/acc04/Papers/0009_WeA02.6.pdf (дата обращения: 24.01.2022).

6. A Simple Discrete-Time Tracking Differentiator and Its Application to Speed and Position Detection System for a Maglev Train / H. Zhang [et al.] // IEEE Transactions on Control Systems Technology. 2018. Vol. 27, Issue 4. P. 1728-1734. doi: https://doi.org/10.1109/TCST.2018.2832139

7. Филимонов Н. Б. Проблема качества процессов управления: смена оптимизационной парадигмы // Мехатроника, автоматизация, управление. 2010. № 12. С. 2-10. URL: https://www.elibrary. ru/item.asp?id=15510887 (дата обращения: 20.01.2022).

8. Табунщиков Ю. А., Бродач М. М. Экспериментальное исследование оптимального управления расходом энергии // АВОК. 2006. № 1. С. 32-36. URL: https://www.abok.ru/for_spec/articles. php?nid=3132 (дата обращения: 20.01.2022).

9. Плешивцева Ю. Э., Попов А. В., Дьяконов А. И. Двумерная задача оптимального по типовым критериям качества управления процессом сквозного индукционного нагрева // Вестник Самарского государственного технического университета. Серия «Технические науки». 2014. № 2.

С. 148-163. URL: https://joumals.eco-vector.com/1991-8542/article/view/19977/16230 (дата обращения: 25.01.2022).

10. Панферов В. И., Анисимова Е. Ю., Нагорная А. Н. Об оптимальном управлении тепловым режимом зданий // Вестник Южно-Уральского государственного университета. Серия: Энергетика. 2007. № 20. С. 3-9. URL: https://cyberleninka.ru/article/n/ob-optimalnom-upravlenii-teplovym-rezhi-mom-zdaniy (дата обращения: 20.01.2022).

11. Biyik E., Kahraman A. A Predictive Control Strategy for Optimal Management of Peak Load, Thermal Comfort, Energy Storage and Renewables in Multi-Zone Buildings [Электронный ресурс] // Journal of Building Engineering. 2019. Vol. 25. URL: https://app.dimensions.ai/details/publication/ pub.1117015634 (дата обращения: 20.01.2022).

12. Хорошавин В. С. Сравнение алгоритмов управления тепловым процессом по быстродействию и на минимум ресурсов // Известия Тульского государственного университета. Технические науки. 2020. № 7. С. 211-216. URL: https://www.elibrary.ru/item.asp?id=43895260 (дата обращения: 20.01.2022).

13. Хорошавин В. С., Грудинин В. С. Синтез программного движения на основе особого оптимального управления // Мехатроника, автоматизация, управление. 2021. Т. 22, № 8. С. 395-403. doi: https://doi.org/10.17587/mau.22.395-403

14. Дубровин В. С., Никулин В. В. Способ построения управляемых функциональных генераторов // Наукоемкие технологии в космических исследованиях Земли. 2013. Т. 5, № 2. С. 16-23. URL: https://www.elibrary.ru/item.asp?id=22897398 (дата обращения: 24.01.2022).

15. Moreau L., Aeyels D. Periodic Output Feedback Stabilization of Single-Input Single-Output Continuous-Time Systems with Odd Relative Degree // Systems & Control Letters. 2004. Vol. 51, Issue 5. P. 395-406. doi: https://doi.org/10.1016/j.sysconle.2003.10.001

16. Шумафов М. М. Стабилизация линейных систем управления. Проблема назначения полюсов. Обзор // Вестник СПбГУ. Математика. Механика. Астрономия. 2019. Т. 6, № 4. C. 564-591. doi: https://doi.org/10.21638/11701/spbu01.2019.404

17. Борковская И. М., Пыжкова О. Н. Задачи управления и стабилизации для гибридных динамических систем // Труды БГТУ Серия 3: Физико-математические науки и информатика. 2018. № 2. С. 5-9. URL: https://elibrary.ru/item.asp?id=36367417 (дата обращения: 24.01.2022).

18. Колесников В. Л., Бракович А. И., Жук Я. А. Решение многокритериальных задач, оптимальных по Парето // Труды БГТУ. Серия 3: Физико-математические науки и информатика. 2014. № 6. С. 128-130. URL: https://elibrary.ru/item.asp?id=27707178 (дата обращения: 24.01.2022).

19. Mahmoud M. S., AL-Sunni F. M. Control and Optimization of Distributed Generation Systems. Cham : Springer, 2015. 578 p. doi: https://doi.org/10.1007/978-3-319-16910-1_1 (дата обращения: 20.01.2022).

20. Review of Optimum Temperature, Humidity, and Vapour Pressure Deficit for Microclimate Evaluation and Control in Greenhouse Cultivation of Tomato: a Review / R. R. Shamshiri [et al.] // International Agrophysics. 2018. Issue 32. P. 287-302. doi: https://doi.org/10.1515/intag-2017-0005

Поступила 25.01.2022; одобрена после рецензирования 11.02.2022; принята к публикации 03.03.2022

Об авторах:

Хорошавин Валерий Степанович, профессор кафедры электропривода и автоматизации промышленных установок Вятского государственного университета (610000, Российская Федерация, г. Киров, ул. Московская, д. 36), доктор технических наук, ORCID: https://orcid.org/0000-0002-4355-3866, Researcher ID: G-5298-2018, khoroshavin@vyatsu.ru

Грудинин Виктор Степанович, доцент кафедры электропривода и автоматизации промышленных установок Вятского государственного университета (610000, Российская Федерация, г. Киров, ул. Московская, д. 36), кандидат технических наук, ORCID: https://orcid.org/0000-0002-1615-6195, Researcher ID: G-5550-2018, grudinin@vyatsu.ru

Заявленный вклад авторов:

В. С. Хорошавин - постановка задачи и выбор методов решения.

В. С. Грудинин - анализ материалов и моделирование процессов.

Все авторы прочитали и одобрили окончательный вариант рукописи.

REFERENCES

1. Domanov V.I., Pevcheva E.V. Analysis of the Greenhouse Complex Main Power Grid Nodes and the Ways to Reduce Energy Consumption. Industrial Automatic Control Systems and Controllers. 2017;(3):3-10. Available at: https://www.elibrary.ru/item.asp?id=29115500 (accessed 20.01.2022). (In Russ., abstract in Eng.)

2. Lovchakov V.I. Approximation Approach to Synthesis of Control Systems Based on Optimal Program Control. News of the Tula State University. Technical Sciences. 2017;(3):225-236. Available at: https://clck.ru/h7ejP (accessed 20.01.2022). (In Russ., abstract in Eng.)

3. Nikol'skii M.S. Singular Sets of Extremal Controls in Optimal Control Problems. Proceedings of the Steklov Institute ofMathematics. 2019;304:236-240. doi: https://doi.org/10.4213/tm3970

4. Tsintsadze Z. [Computation of Particularly Optimal Control in Quasi-Linear Controlled Systems with Mixed Constraints]. Kompyuternye nauki i telekommunikatsii. 2005;(2):71-73. Available at: http://gesj. internet-academy.org.ge/ru/list_aut_artic_ru.php?b_sec=&list_aut=1248 (accessed 24.01.2022). (In Russ.)

5. Gao Z. On Discrete Time Optimal Control: A Closed-Form Solution. In: Proceeding of the 2004 American Control Conference (30 June - 2 July 2004). Boston; 2004. p. 52-58. Available at: https://folk. ntnu.no/skoge/prost/proceedings/acc04/Papers/0009_WeA02.6.pdf (accessed 24.01.2022).

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

6. Zhang H., Xie Y., Xiao G., et al. A Simple Discrete-Time Tracking Differentiator and Its Application to Speed and Position Detection System for a Maglev Train. IEEE Transactions on Control Systems Technology. 2018;27(4):1728-1734. doi: https://doi.org/10.1109/TCST.2018.2832139

7. Filimonov N.B. The Problem of Quality of Control Processes: Change of an Optimizing Paradigm. Mechatronics, Automation, Management. 2010;(12):2-10. Available at: https://www.elibrary.ru/item. asp?id=15510887 (accessed 20.01.2022). (In Russ., abstract in Eng.)

8. Tabunshchikov Yu.A., Brodach M.M. [Experimental Study of Optimal Energy Management]. Ven -tilation, Heating, Air Conditioning, Heat Supply and Building Thermal Physics. 2006;(1):32-36. Available at: https://www.abok.ru/for_spec/articles.php?nid=3132 (accessed 20.01.2022). (In Russ.)

9. Pleshivtseva Yu.E., Popov A.V., Dyakonov A.I. Two-Dimensional Problem of Optimal with Respect to Typical Quality Criteria Control of through Induction Heating Processes. Samara State Technical University Bulletin. Technical Sciences Series. 2014;(2):148-163. Available at: https://journals.eco-vector. com/1991-8542/article/view/19977/16230 (accessed 25.01.2022). (In Russ., abstract in Eng.)

10. Panferov V.I., Anisimova Ye.Yu., Nagornaya A.N. [On the Optimal Management of the Thermal Regime of Buildings]. Bulletin of South Ural State University. Power Engineering Series. 2007;(20):3-9. Available at: https://cyberleninka.ru/article/n/ob-optimalnom-upravlenii-teplovym-rezhimom-zdaniy (accessed 20.01.2022). (In Russ., abstract in Eng.)

11. Biyik E., Kahraman A. A Predictive Control Strategy for Optimal Management of Peak Load, Thermal Comfort, Energy Storage and Renewables in Multi-Zone Buildings. Journal of Building Engineering. 2019;25. Available at: https://app.dimensions.ai/details/publication/pub.1117015634 (accessed 20.01.2022).

12. Khoroshavin V.S. Comparing Optimal Control Algorithms by Time and Minimum Resources Heat. News of the Tula State University. Technical Sciences. 2020;(7):211-216. Available at: https://www. elibrary.ru/item.asp?id=43895260 (accessed 20.01.2022). (In Russ., abstract in Eng.)

13. Khoroshavin V.S., Grudinin V.S. Synthesis of Programmed Motion Based on Special Optimal Control. Mechatronics, Automation, Control. 2021;22(8):395-403. (In Russ., abstract in Eng.) doi: https:// doi.org/10.17587/mau.22.395-403

14. Dubrovin V.S., Nikulin V.V. [A Way of Constructing Controlled Function Generators]. High Tech in Earth Space Research. 2013;5(2):16-23. Available at: https://www.elibrary.ru/item.asp?id=22897398 (accessed 24.01.2022). (In Russ.)

15. Moreau L., Aeyels D. Periodic Output Feedback Stabilization of Single-Input Single-Output Continuous-Time Systems with Odd Relative Degree. Systems & Control Letters. 2004;51(5):395-406. doi: https://doi.org/10.1016/j.sysconle.2003.10.001

16. Shumafov M.M. Stabilization of Linear Control Systems. Pole Assignment Problem. A Survey. Vestnik SPbGU. Matematika. Mekhanika. Astronomiya. 2019;6(4):564-591. (In Russ., abstract in Eng.) doi: https://doi.org/10.21638/11701/spbu01.2019.404

17. Borkovskaya I.M., Pyzhkova O.N. The Problems of Control and Stabilization for Hybrid Dynamic Systems. Trudy BGTU. Seriya 3: Fiziko-matematicheskie nauki i informatika. 2018;(2):5-9. Available at: https://elibrary.ru/item.asp?id=36367417 (accessed 24.01.2022). (In Russ., abstract in Eng.)

18. Kolesnikov V.L., Brakovich A.I., Zhuk Ya.A. [Solving Pareto-Optimal Multicriteria Problems]. Trudy BGTU. Seriya 3: Fiziko-matematicheskie nauki i informatika. 2014;(6):128-130. Available at: https://elibrary.ru/item.asp?id=27707178 (accessed 24.01.2022). (In Russ.)

19. Mahmoud M.S., AL-Sunni F.M. Control and Optimization of Distributed Generation Systems. Cham: Springer; 2015. 578 p. doi: https://doi.org/10.1007/978-3-319-16910-1_1 (accessed 20.01.2022).

20. Shamshiri R.R., Jones J.W., Thorp K.R., et al. Review of Optimum Temperature, Humidity, and Vapour Pressure Deficit for Microclimate Evaluation and Control in Greenhouse Cultivation of Tomato: a Review. InternationalAgrophysics. 2018;(32):287-302. doi: https://doi.org/10.1515/intag-2017-0005

Submitted 25.01.2022; approved after reviewing 11.02.2022; accepted for publication 03.03.2022

About the authors:

Valeriy S. Khoroshavin, Professor of the Chair of Electric Drive and Industrial Equipment Automation, Vyatka State University (36 Moskovskaya St., Kirov 610000, Russian Federation), Dr.Sci. (Engr.), ORCID: https://orcid.org/0000-0002-4355-3866, Researcher ID: G-5298-2018, khoroshavin@vyatsu.ru

Viktor S. Grudinin, Associate Professor of the Chair of Electric Drive and Industrial Equipment Automation, Vyatka State University (36 Moskovskaya St., Kirov 610000, Russian Federation), Cand.Sci. (Engr.), ORCID: https://orcid.org/0000-0002-1615-6195, Researcher ID: G-5550-2018, grudinin@vyatsu.ru

Contribution of the authors:

V. S. Khoroshavin - problem statement and choice of solution methods.

V. S. Grudinin - material analysis and process modeling.

All authors have read and approved the final manuscript.

i Надоели баннеры? Вы всегда можете отключить рекламу.