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

Разработка математической модели стадии нагрева резиновой смеси и синтез алгоритма управления нагревом с использованием принципа максимума Понтрягина Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
137
16
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
РЕЗИНОСМЕШЕНИЕ / МАТЕМАТИЧЕСКАЯ МОДЕЛЬ / АЛГОРИТМ УПРАВЛЕНИЯ / УПРАВЛЯЮЩИЙ КОНТРОЛЛЕР / ПРИНЦИП МАКСИМУМА ПОНТРЯГИНА / RUBBER MIXING / MATHEMATICAL MODEL / CONTROL ALGORITHM / MAIN CONTROLLER / PONTRYAGIN MAXIMUM PRINCIPLE

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

Статья посвящена разработке алгоритма управления стадией нагрева резиновой смеси для ЗАО «Воронежский шинный завод». Алгоритм предназначен для реализации на базе контроллера Siemens S-300 для управления резиносмесителем РС-270. Для составления алгоритма проведен системный анализ процесса нагрева как объекта управления и разработана математическая модель стадии нагрева на основе уравнений тепловых балансов, описывающая процессы нагрева теплоносителя в теплообменнике и дальнейшего нагрева им смеси в резиносмесителе. Получены динамические характеристики по температуре теплообменного аппарата и резиносмесителя. В виду сложности и нелинейности объекта управления – резиносмесителя, а также наличия методов и большого опыта управления данным аппаратом в промышленных условиях, алгоритм реализован с использованием принципа максимума Понтрягина. Задача оптимизации сводится к определению оптимального управления (подача греющего пара) и оптимальной траектории выходной координаты объекта (температура смеси), обеспечивающих минимальный расход пара при нагреве резиновой смеси за ограниченное время. Для этого математическая модель стадии нагрева записана в матричном виде, сформированы матрицы коэффициентов для каждого состояния управления, векторы управления и возмущения, получена функция Гамильтона и для построения оптимального управления и траектории выхода объекта найдены временные точки переключения. Анализ результатов модельных экспериментов и практических исследований при программировании контроллера, показал снижение расхода греющего пара на 24,4% при проведении стадии нагрева резиновой смеси.

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

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Кудряшов В.С., Алексеев М.В., Юдаков А.А.

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

Development of a mathematical model of the heating phase of rubber mixture and development of the synthesis of the heating control algorithm using the Pontryagin maximum principle

The article is devoted to the development of the algorithm of the heating phase control of a rubber compound for CJSC “Voronezh tyre plant”. The algorithm is designed for implementation on basis of controller Siemens S-300 to control the RS-270 mixer. To compile the algorithm a systematic analysis of the heating process has been performed as a control object, also the mathematical model of the heating phase has been developed on the basis of the heat balance equation, which describes the process of heating of a heat-transfer agent in the heat exchanger and further heating of the mixture in the mixer. The dynamic characteristics of temperature of the heat exchanger and the rubber mixer have been obtained. Taking into account the complexity and nonlinearity of the control object – a rubber mixer, as well as the availability of methods and great experience in managing this machine in an industrial environment, the algorithm has been implemented using the Pontryagin maximum principle. The optimization problem is reduced to determining the optimal control (heating steam supply) and the optimal path of the object’s output coordinates (the temperature of the mixture) which ensure the least flow of steam while heating a rubber compound in a limited time. To do this, the mathematical model of the heating phase has been written in matrix form. Coefficients matrices for each state of the control, control and disturbance vectors have been created, the Hamilton function has been obtained and time switching points have been found for constructing an optimal control and escape path of the object. Analysis of the model experiments and practical research results in the process of programming of the controller have showed a decrease in the heating steam consumption by 24.4% during the heating phase of the rubber compound.

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

BeemHunJBry^T/Proceedings of VSUET, Т. 79, № 2, 207

Оригинальная статья/Original article_

УДК 001.891.573

DOI: http://doi.org/10.20914/2310-1202-2017-2-80-87

Разработка математической модели стадии нагрева резиновой

смеси и синтез алгоритма управления нагревом с _использованием принципа максимума Понтрягина_

Владимир С. Кудряшов 1 kudryashovvs@mail.ru Михаил В. Алексеев 1 mwa1976@mail.ru Андрей А. Юдаков 2 yudakov07@yandex.ru

1 Воронежский государственный университет инженерных технологий, пр-т Революции, 19, г. Воронеж, 394036, Россия

2 ЗАО «Воронежский шинный завод», ул. Ростовская, 41, г. Воронеж, Россия

Реферат. Статья посвящена разработке алгоритма управления стадией нагрева резиновой смеси для ЗАО «Воронежский шинный завод». Алгоритм предназначен для реализации на базе контроллера Siemens S-300 для управления резиносмесителем РС-270. Для составления алгоритма проведен системный анализ процесса нагрева как объекта управления и разработана математическая модель стадии нагрева на основе уравнений тепловых балансов, описывающая процессы нагрева теплоносителя в теплообменнике и дальнейшего нагрева им смеси в резиносмесителе. Получены динамические характеристики по температуре теплообменного аппарата и резиносмесителя. В виду сложности и нелинейности объекта управления - резиносмесителя, а также наличия методов и большого опыта управления данным аппаратом в промышленных условиях, алгоритм реализован с использованием принципа максимума Понтрягина. Задача оптимизации сводится к определению оптимального управления (подача греющего пара) и оптимальной траектории выходной координаты объекта (температура смеси), обеспечивающих минимальный расход пара при нагреве резиновой смеси за ограниченное время. Для этого математическая модель стадии нагрева записана в матричном виде, сформированы матрицы коэффициентов для каждого состояния управления, векторы управления и возмущения, получена функция Гамильтона и для построения оптимального управления и траектории выхода объекта найдены временные точки переключения. Анализ результатов модельных экспериментов и практических исследований при программировании контроллера, показал снижение расхода греющего пара на 24,4% при проведении стадии нагрева резиновой смеси.

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

Development of a mathematical model of the heating phase of rubber

mixture and development of the synthesis of the heating control _algorithm using the Pontryagin maximum principle_

Vladimir S. Kudryashov 1 kudryashovvs@mail.ru Mikhail V. Alekseev 1 mwa1976@mail.ru Andrei A. Yudakov 2 yudakov07@yandex.ru

1 Voronezh state university of engineering technologies, Revolution Av., 19 Voronezh, 394036, Russia

2 JSC "Voronezh Tire Plant", Rostovskaya str., 41, Voronezh, Russia

Summary. The article is devoted to the development of the algorithm of the heating phase control of a rubber compound for CJSC "Voronezh tyre plant". The algorithm is designed for implementation on basis of controller Siemens S-300 to control the RS-270 mixer. To compile the algorithm a systematic analysis of the heating process has been performed as a control object, also the mathematical model of the heating phase has been developed on the basis of the heat balance equation, which describes the process of heating of a heat-transfer agent in the heat exchanger and further heating of the mixture in the mixer. The dynamic characteristics of temperature of the heat exchanger and the rubber mixer have been obtained. Taking into account the complexity and nonlinearity of the control object - a rubber mixer, as well as the availability of methods and great experience in managing this machine in an industrial environment, the algorithm has been implemented using the Pontryagin maximum principle. The optimization problem is reduced to determining the optimal control (heating steam supply) and the optimal path of the object's output coordinates (the temperature of the mixture) which ensure the least flow of steam while heating a rubber compound in a limited time. To do this, the mathematical model of the heating phase has been written in matrix form. Coefficients matrices for each state of the control, control and disturbance vectors have been created, the Hamilton function has been obtained and time switching points have been found for constructing an optimal control and escape path of the object. Analysis of the model experiments and practical research results in the process of programming of the controller have showed a decrease in the

heating steam consumption by 24.4% during the heating phase of the rubber compound._

Keywords: rubber mixing, mathematical model, control algorithm, main controller, Pontryagin maximum principle

Для цитирования Кудряшов В.С., Алексеев М.В., Юдаков А.А. Разработка математической модели стадии нагрева резиновой смеси и синтез алгоритма управления нагревом с использованием принципа максимума Понтрягина // Вестник ВГУИТ. 2017. Т. 79. № 2. С. 80-87. (1ся: 10.20914/2310-1202-2017-2-80-87

For citation

Kudryashov V.S.,Alekseev M.V.,Yudakov A.A. Development of a mathematical model of the heating phase of rubber mixture and development of the synthesis of the heating control algorithm using the Pontryagin maximum principle. Vestnik VGUIT [Proceedings of VSUET]. 2017. vol. 79. no. 2. pp. 80-87. (in Russian). doi : 10.20914/2310-1202-2017-2-80-87

Введение

Приготовление и обработка резиновых смесей - важнейшая операция в шинном и резиновом производствах, от которой зависят пластоэластические и прочностные свойства, структура и гомогенность, а также технологичность смесей резиновых и резинокордных полуфабрикатов [1, 2]. На одном из участков в ЗАО «Воронежский шинный завод» используется резиносмеситель РС-270. Процесс приготовления резиновой смеси начинается с дозирования и загрузки сыпучих материалов из приемных бункеров в резиносмеситель. Заполнение производится системой шнековых питателей и весовых дозаторов. После происходит загрузка жидких ингредиентов, подогретых до необходимых температур. Часть операций производится в ручном режиме. Для отработки рецептов предусмотрен контроллер Siemens S-300 с дискретными и аналоговыми модулями расширения. Контроллер управляет дозированием компонентов, а также выгрузкой резиновой смеси по температуре и времени смешения.

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

Результаты

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

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

Основная часть

Разработано математическое описание стадии нагрева [3-5], представляющее собой уравнения тепловых балансов для малого интервала времени Дт и описывающее процессы нагрева теплоносителя в теплообменнике и дальнейшего нагрева им смеси в резиносмесителе:

для конденсата

Ор-I-Ат = Ок-Тк-ек-Ат +

(1)

+8гр- аг* - (Тк -Ts)- Ат + Мк - ек- АТк;

для воды, нагреваемой в теплообменнике

Оу - Туо - еу- Ат + 8 - ар* - (Т* - Ту) - Ат =

= Оу- Ту- еу- Ат + Му - еу - АТу;

для стенки между конденсатом и водой в теплообменнике

8-аг*-(Тк-Т*)-Ат =

(3)

= 8-ар(Т*-Ту)-Ат + М* -еАТ*;

для воды, выходящей из резиносмесителя (после теплообмена)

Оу- Ту- еу- Ат = Оу - Туо - еу- Ат + 8 2 -

(4)

- а2г* - (Туо -Т2*)- Ат+ Муо - еу - АТуо;

для стенки между водой и резиновой смесью в резиносмесителе

82 - а2& - (Ту - Т2*) - Ат = 82 - а2р* -

- (Т 2* -Туо)- Ат + М 2*- е* - АТ 2*; (5)

для резиновой смеси, нагреваемой в рези-носмесителе

N - Кр-М*т + 82 - а2р* - (Т2* - Туо) -

- Ат = М*т - е*т - АТ*т,

где Ор - массовый расход пара; Оу - массовый расход воды; I - энтальпия пара; Тк - температура конденсата на выходе теплообменника; Ту - температура воды на выходе теплообменника и входе резиносмесителя; Туо - температура воды на выходе резиносмесителя и входе теплообменника; ек, су, ее*т - теплоемкости конденсата, воды, стенки и резиновой смеси; Мк, Му, Муо, Мя, М2 М*т - массы конденсата в теплообменнике, воды в теплообменнике и резиносмесителе, стенки между конденсатом и водой в теплообменнике, стенки между водой и смесью в резиносмесителе,

а также смеси; Ts, T2s - среднее значение температуры стенок между конденсатом и водой, среднее значение температуры стенок между водой и смесью; S1, S2 - площади стенок в резиносмесителе и теплообменнике; ats, aps, a2ts, a2ps - коэффициенты теплоотдачи от конденсата к стенке, от стенки к воде, от воды к стенке и от стенки к резиновой смеси; N - мощность электродвигателя

dTk Gp ■I

рабочего органа смесителя; Кр - коэффициент полезного действия привода резиносмесителя; т - время.

Запишем уравнения (1)-(6) в виде уравнений Коши. Для этого разделим левые и правые части уравнений на Ах, перенесем производную в левую часть уравнений, а все остальное - в правую, и получим систему дифференциальных уравнений:

Ок • Тк • ск Stp • • (Тк - Ts)

dt Mk ■ ck Mk ■ ck Mk ■ ck

dTv Gv ■ Tvo ■ cv + S ■aps ■ (Ts - Tv) Gv ■ Tv ■ cv dt Mv ■ cv Mv ■ cv Mv ■ cv

dTs S ■ ats ■ (Tk - Ts) S ■aps ■ (Ts - Tv)

dTvo dt dT 2s

dz Ms ■ cs

Gv ■ Tv ■ cv Gv ■ Tvo ■ cv

Ms ■ cs S2 ■ a2ts ■ (Tvo - T2s) Mvo ■ cv Mvo ■ cv Mvo ■ cv

S2 ■ a2ts ■ (Tv - T2s) S2 ■ a2ps ■ (T2s - Tvo)

(7)

dx

dTsm dt

M 2 s ■ cs N ■ Kp ■ Msm Msm ■ csm

Рассчитаем динамические характеристики стадии нагрева (рисунки 1, 2).

Т, °С

¿S* Конденсат " Condeiisale

// » * * t : * / • ' с t Ст *нка de

l* * / \ в Т1 ода на е water выходе at the exi t

М 2 5 • cs S2 • а2р5 • (Т25 - Туо)

М5Ш • С5Ш

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

Приведем уравнения системы (7) к матричной форме:

" 5 57.5 75 112 5 150 1S7J ::5 1Ы.5 300

t, С

Рисунок 1. Динамические характеристики по температуре теплообменного аппарата

Figure 1. The dynamic characteristics of the temperature of heat exchanger

T oC

dX = Л ■ X.

B ■ U + C ■ Z.

(8)

Tvo, -T:sl Tsni^

Вода The w на бых ater at t \ эде из ie exit о месите f the mb ЛЯ ter q генка с Mixer «ecKTej side [я

^ * * * • * * * : t * r t t \Cw B1 есь aid

** ** it *t .t :t /

300 t, С

Рисунок 2. Динамические характеристики по температуре резиносмесителя

Figure 2. Dynamic characteristics by Rubber mixer temperature

где X, и, Ъ - векторы переменных состояния, управления и возмущения; А, В, С - матрицы коэффициентов.

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

В уравнениях вектор переменных состояния

имеет вид {Тк Ту Туо T2s Tsm}T ; вектор управления - {Ор 0 0 0 0 }Т; вектор возмущений -(К • Кр 0 0 0 0}Т. Так как заведомо предполагая, что расход поступающего греющего пара Ор равен расходу конденсата, выходящего из теп-лообменного аппарата, то можно заметить, что система уравнений нелинейна. Но если организовать управление подачей греющего пара в двух состояниях: греющий пар подается в максимальном расходе или не подается совсем, то система уравнений (7) в каждом

<Ветшк&ГУИТ/Фгоиг£^ о/Т. 79, № 2, 2017,

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

81-Б + Ск-ек

А1 :=

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

Мк-ек 81-Б Му-еу

0 0

А2 :=

81^ Мк-ек 81-арз Му-еу

0 0

Му2- еу 0

0

81^ Мк-ек 81-аps + Су - еу Му-еу Су - еу

Му2- еу 0

В соответствии с векторами переменных состояний запишем матрицы коэффициентов для каждого состояния управления:

81-Б Мк-ек 81-Б + Су-еу Му-еу Су-еу

0

Су-еу Му-еу S2-a2ts + Су -еу Му2- еу S2-a2ts

Ms-cs 0

0

Су - еу

Му- еу S2-a2ts + Су-еу Му2- еу 82- a2ts

0

Ms-еs 0

0 0

S2-a2ts Му2- еу 82- (a2ts + a2ps) Ms-еs 82-a2ps Msm- еsm

0 0

82-a2ts Му2- еу S2-(a2ts + a2ps) Ms-еs 82-a2ps Msm- еsm

0 0 0

82-a2ps Ms- еs 82- a2ps Msm-еsm

0 0

82-a2ps Ms- cs 82- a2ps Msm-еsm

(9)

Также запишем и возмущения:

( ]раг ^ Мк-ек 0 0 0

V 0 .

В :=

векторы

Г

С :=

МЛ/

0 0 0 0

К-КР

управления

Л

(10)

V

У.

Чтобы получить математическое описание в виде зависимости выхода объекта Tsm от управления Gp и возмущения N - ^ или в виде соответствующих передаточных функций, можно исключить из уравнений (7) промежуточные переменные Tk, Tvo, T2s и получить одно уравнение пятого порядка для Tsm. Это проще всего выполнить, переходя к преобразованиям Лапласа:

X (s,щ z) = (E - s - Al)-1 - B - u + +(E - s - A2)-1 - C - z,

(11)

где Е - единичная матрица; s - комплексная переменная Лапласа.

Так как уравнение записано для приращений, все начальные значения переменных равны нулю, что удовлетворяет требованиям для передаточной функций [7-9].

Из уравнения (11) достаточно найти вектор -Ж^), используя операции с матрицами в системе Mаthcаd. После нахождения передаточных функций для выхода объекта, находим корни характеристических уравнений для двух состояний управления:

( -3.3092024252193001718 ^

Х2 :=

-.50525582941656373523 -.19599241914230279279 -.24858249769714896775е-1 ^.20955764521184033564е-2 у ( 0 > -3.3092023859043633231 -.50357336558886473926 -.16758754463218023878 ^.23707838501281881098е-1 у

(12)

Функция Гамильтона для данной системы будет выглядеть как:

0

0

0

h = ть

Gp • I Gk • ck + Stp • ats

Mk- ck Y3

Mk- ck

• Tk

+ Y2<

Gv • cv Gv • cv + S • aps --Tvo---Tv

Mv cv

Mv cv

Gv^ cv Gv^ cv + S2^ a2ts S2^ a2ts

• Tv---Tvo +---T 2s

Mvo cv

Mvo cv

Y4

S 2^ a2ts S 2^ a2ts S 2^ a2ps• (T 2s -Tvo) --Tv---T 2s--—---

. M2s • cs

M 2s • cs

M2s cs

+T5

Mvo cv

N Kp • Msm + S2^ a2ps • (T2s-Tvo)

(13)

В соответствии с принципом максимума Понтрягина оптимальное управление должно обеспечивать максимум функции Гамильтона в каждый момент времени в явном виде, т. е.

u = <

fa, если ц/1 • Kim > 0 [6, если ц1 • Kim < 0

(14)

Мш - с'ш М'ш - с'ш

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

ф (4, С, п ) = С1 -s1n-es14 +

(16)

-е"44 + С -'5п -е"5-4

+С2 • s2n • es

+ С3 • s3" • e

+C4 • s4n • es4 •t + C5 •s5n

t

s 5 4

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

Решение для стадии без подачи греющего пара будет иметь вид:

y (t, C,uc) = C1 + C2 • e122 •t +

y (t, C,uc) = C1 •esl •t + C2 e21 + +C3 •es3 •t + C4 •es4 •t + C5 •e +Ko •uc + Q.

+C3 •e +C5 •e

s33 • t s 551

+ C4 • e

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

(17)

Q.

(15)

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

dy (t, C, n ) =

|C2 • s22n • e ' + C3 • s33n • < C, •s 22"^ 22J + a•s33n

+ C4 • s44n • es44'' + C5 • s55n • es55 ' + Q, если n = 1;

331 , /-1 л л n s44^t , - -n s55.t ^ Л

+ C. • s44 •e + C. •s55 •e , если n > 1.

(18)

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

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

Примем в начальный момент управление положительным и равным а. Рассчитаем выход объекта при управлении и=ис для первого участка постоянства управления (рисунок 3).

Рисунок 3. Выход объекта управления для первого участка Figure 3. The output of the control object for the first section

Задавая последовательно точки переключения в пределах положительных значений у, определим начальные приближения констант интегрирования по участкам постоянства управления. Для второго участка в качестве начального приближения для констант интегрирования С2 возьмем С1, рассчитанные для первого участка. Т. °с

Выражение для выхода объекта при управлении и=0 для второго участка будет иметь вид (рисунок 4):

, ч Г у (/, С1, а ), если t < t1 У1( /, /1)4 (19)

[у1 (/ - /1, С2, Ь ), если t > t1

120 100

Yl(t,tl) 80 60 40

0 100 Рисунок 4. Выход объекта управления для второго участка Figure 4. The output of the control object for the second section

200 t,,

Для третьего участка принимаем вектор констант С3=С2. Точку переключения также задаем приближенно по графику: /2=80.

Выражение для выхода объекта при управлении и=ис для третьего участка будет иметь вид (рисунок 5):

ч ГУ1(/,/1), если t < 12 У 2 (/, /1, /2 Н (20)

[ у (/ -/2, С3, а), если 1 >12

Т,°С 200

150

Для четвертого участка принимаем вектор констант С4=С3. Точку переключения также задаем приближенно по графику: /3=85.

Выражение для выхода объекта при управлении и=0 для четвертого участка будет иметь вид (рисунок 6):

ч ГУ 2 (/,/1, /2), если 1 < 13 У 3 (/,/1, /2, /3 Н (21)

[ у1(/ -/3, С4, Ь), если 1 >13

Y2(t,tl,t2)l00

0 100 Рисунок 5. Выход объекта управления для третьего участка

Figure 5. The output of the control object for the third section

T,°C

200

t, с

100

Y3(t,tl,t2,t3)

50

0'

0 100 Рисунок 6. Выход объекта управления для четвертого участка

Figure 6. The output of the control object for the fourth section

200

t, с

Для пятого участка вновь принимаем вектор констант С5=С4. Точку переключения также задаем приближенно по графику: ^=100.

Выражение для выхода объекта при управлении и=ис для пятого участка будет иметь вид (рисунок 7):

7 (t)=

Y4(t,t1,t2,t3,t4) =

[Y3 (t, t1, t2, t3), если t < t4 [y(t -t4, C5, a), если t >t4

(22)

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

т, °с

100

Описания координаты у и функций управления и (рисунки 8, 9):

у С1, а), если 1 < 11; у1 (t-t1, С2, Ь), если 1 < 12 л 1 > 11; у^ -12, С3, а), если 1 < 13 л 1 > 12; (23) у1^ -13, С4, Ь), если 1 < 14 л 1 > 13; у( t -4, С5, а), если 1 >14; а, если 1 < 11; Ь, если 1 < 12 л 1 > 11;

a, если 1 < 13 л 1 > 12; (24)

b, если 1 < 14 л 1 > 13; а, если 1 > 14.

UU (t ) =

Y4(t,tl ,t2,t3,t4)

Рисунок 7. Выход объекта управления для пятого участка Figure 7. The output of the control object for the fifth section

T,°Ci 20 110 inn 90 SO

Y(l) 70 60 50 -10 30 20

0 50 100

Рисунок 8. Оптимальная траектория выхода объекта управления (температура смеси) Figure 8. The optimal path of the control object's output (temperature of the compound)

t, с

O.Sr

0.6-

UU(t)0.4

0.2"

0 50

Рисунок 9. Оптимальное управление объектом (подача греющего пара) Figure 9. Optimal control of the object (supply of the heating steam)

100

t, с

Сравнивая результаты (при постоянной подаче греющего пара в теплообменный аппарат и при реализации оптимального управления), установлено, что время необходимое для разогрева резиновой смеси до 110 °С при постоянной подаче пара составляет 110 с (с общим количеством поданного пара 55 кг), а при оптимальном управлении - 130 с (с общим количеством поданного пара 41,5 кг). Таким образом, экономия греющего пара составила 24,4%.

Выводы и рекомендации

Анализ результатов моделирования и проведения испытаний на предприятии при

ЛИТЕРАТУРА

1 Резниченко С.В., Морозов Ю.Л. Большой справочник резинщика. Том 1. Каучуки и ингредиенты. М.: Техинформ, 2012. 744 с.

2 Резниченко С.В., Морозов Ю.Л. Большой справочник резинщика. Том 2. Резины и резинотехнические изделия. М.: Техинформ, 2012. 648 с.

3 Красных В.К)., Королев В.Н. Тепломассообмен. Основные формулы, задачи и способы их решения. Екатеринбург: УрФУ, 2012. 64 с.

4 Алексеев MB., Юдаков A.A. Разработка математического описания и алгоритмов управления процессом резиносмешения // Сборник трудов XXVIII международной научной конференции ММТТ 28. 2015. С. 34-37.

5 КудряшовВ.С., Алексеев М.В., Юдаков A.A., Попов А.П., Чертов Е.Д. Разработка математической модели процесса резиносмешения и синтез корректирующего алгоритма управления температурным режимом процесса с использованием искусственной нейронной сети // Вестник ВГУИТ. 2016. № 2. С. 52-59.

6 Программирование с помощью STEP 7 V5.3: Руководство. Siemens AG, 2004. 602 с.

7 Candido J. J., Justino P. A. Modelling, control and pontryagin maximum principle for a two-body wave energy device//Renewable Energy. 2011. T. 36. №. 5. C. 1545-1557."

8 Laschov D., Margaliot M A maximum principle for single-input Boolean control networks //IEEE Transactions on Automatic Control. 2011. T. 56. №. 4. C. 913-917.

9 Kim N., Cha S., Peng H. Optimal control of hybrid electric vehicles based on Pontryagin's minimum principle //IEEE Transactions on Control Systems Technology. 2011. Т. 19. №. 5. С. 1279-1287.

СВЕДЕНИЯ ОБ АВТОРАХ

Владимир С. Кудряшов д.т.н., профессор, кафедра информационных и управляющих систем, Воронежский государственный университет инженерных технологий, пр-т Революции, 19, г. Воронеж, 394036, Россия, kudryashovvs@mail.ru Михаил В. Алексеев к.т.н., доцент, кафедра информационных и управляющих систем, Воронежский государственный университет инженерных технологий, пр-т Революции, 19, г. Воронеж, 394036, Россия, mwa1976@mail.ru

Андрей А. Юдаков инженер, ЗАО "Воронежский шинный завод", ул. Ростовская, 41, г. Воронеж, Россия, yudakov07@yandex.ru

КРИТЕРИЙ АВТОРСТВА

Все авторы в одинаковой степени принимали участие в подготовке материала для рукописи.

КОНФЛИКТ ИНТЕРЕСОВ

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

ПОСТУПИЛА 02.02.2017 ПРИНЯТА В ПЕЧАТЬ 12.04.2017

программировании управляющего контроллера Siemens S-300 [6], показал эффективность синтезированного алгоритма управления стадией нагрева.

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

Синтезированный алгоритм можно использовать для настройки контроллера в системе управления стадией нагрева в резиносмесителе РС-270.

REFERENCES

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

1 Reznichenko S.V., Morozov Yu.L., ingredienty Bol'shoi spravochnik rezinshchika. Tom 1. Kauchuki i ingredient [A large directory of reinsta. Volume 1. Rubber and ingredient]. Moscow. Tekhinform 2012. 744 p. (in Russian).

2 Reznichenko S.V., Morozov Yu.L. Bol'shoi spravochnik rezinshchika. Tom 2. Reziny i rezinotekhnicheskie izdeliya [ A large directory of reinsta. Volume 2. Rubber and rubber products]. Moscow: Tekhinform 2012. 648 p. (in Russian).

3 Krasnykh V.Yu., Korolev V.N. Teplomassoobmen. Osnovnye formuly, zadachi i sposoby ikh resheniya [Heat and mass transfer. Basic formulas, tasks and ways of their solution]. Ekaterinburg. UrFU 2012. 64 p. (in Russian).

4 Alekseev M.V., Yudakov A.A. The development of the mathematical description and control algorithms of the rubber mixing process. Proceedings of the XXVIII International Scientific Conference MMTT 28. Ryazan, 2015. pp. 34-37. (in Russian).

5 Kudryashov V.S., Alekseev MV., Yudakov A.A., Popov A.P., Chertov E.D. Development of rubber mixing process mathematical model and synthesis of control correction algorithm by process temperature mode using an artificial neural network. Vestnik VGUT [Proceeding of VSUET]. 2016. no. 2 pp. 52-59. (in Russian).

6 Programming with STEP 7 V5.3: Guide. Siemens AG, 2004. 602 p. (in Russian).

7 Cándido J. J., Justino P. A. Modelling, control and pontryagin maximum principle for a two-body wave energy device. Renewable Energy. 2011. vol. 36. no. 5. pp. 1545-1557.

8 Laschov D., Margaliot M A maximum principle for single-input Boolean control networks. IEEE Transactions on Automatic Control. 2011. vol. 56. no. 4. pp. 913-917.

9 Kim N., Cha S., Peng H. Optimal control of hybrid electric vehicles based on Pontryagin's minimum principle. IEEE Transactions on Control Systems Technology. 2011. vol. 19. no. 5. pp. 1279-1287.

INFORMATION ABOUT AUTHORS Vladimir S. Kudryashov doctor of technical sciences, professor, information and control systems department, Voronezh state university of engineering technologies, Revolution Av., 19 Voronezh, 394036, Russia, kudryashovvs@mail.ru Mikhail V. Alekseev candidate of technical sciences, engineer, information and control systems department, Voronezh state university of engineering technologies, Revolution Av., 19 Voronezh, 394036, Russia, mwa1976@mail.ru Andrei A. Yudakov engineer, JSC "Voronezh Tire Plant", Rostovskaya str., 41, Voronezh, Russia, yudakov07@yandex.ru

CONTRIBUTION

All authors equally participated in the preparation of material for manuscript.

CONFLICT OF INTEREST The authors declare no conflict of interest.

RECEIVED 2.2.2017 ACCEPTED 4.12.2017

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