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

УДК 519.854

БОТ: 10.20310/1810-0198-2016-21-2-463-470

ТРИЛИНЕЙНАЯ ОКРЕСТНОСТНАЯ МОДЕЛЬ ПРОЦЕССА ФОРМИРОВАНИЯ ТЕМПЕРАТУРЫ СМОТКИ ГОРЯЧЕКАТАНОЙ ПОЛОСЫ

© А. М. Ш^мырин, А. Г. Ярцев, В. В. Правильникова

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

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

1. Введение

Окрестностные модели [1, 2] являются перспективным направлением в моделировании сложных систем, т. к. отличаются гибкостью описания с помощью окрестностей структуры связей между узлами системы. Теория окрестностных систем является универсальным средством моделирования класса дискретных распределенных систем: стационарных и динамических, четких и нечетких, с линейными и нелинейными связями. Данная теория является обобщающей для многих дискретных распределенных систем, использование окрестностного подхода в прикладных задачах охватывает системы разной природы и направленности - от технических до экономических. Окрестностные модели - исключительно эффективный метод моделирования, позволяющий воспроизводить сложные зависимости.

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

Микроструктура горячекатаной стали зависит от температуры конца прокатки и температуры смотки. Для получения мелкозернистой равномерной структуры металла низкоуглеродистой стали ее прокатку необходимо заканчивать при температуре 760 — 900° С, а сматывать рулон при температуре смотки 540 — 720° С [3, 4]. Охлаждение полосы после чистовой группы клетей осуществляется на отводящем рольганге. Отводящий рольганг стана горячей прокатки представляет собой значительный участок пути, через который проходит полоса после прокатки до моталок. Основная функция участка отводящего рольганга - охлаждение полосы от температуры прокатки до температуры, приемлемой для смотки. Сверху и снизу полосы, движущейся на рольганге, установлены бачки с водой для охлаждения. Бачки объединены в полусекции - минимальные управляемые устройства для охлаждения. Система состоит из 40 секций (80 полусекций) охлаждения полосы. Вода душирует поверхность полосы через ряд трубок, установленных по ширине бачка. Для управления последовательным включением и выключением бачков на стане установлена автоматизированная система управления душиру-ющей установкой.

2. Идентификация

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

^ ■шх[а,а]Х[а]+ ^ [а,в]У[в] + ^ [а>ч]уЫ +

a£Ox[a]

ß&Ov [а]

Y&Oy[a]

+ ^ wxv [a,a,ß]X [a]V ß] + ^ wXy [a,a,7]X [a]Y [7] +

a£Ox[a] ß&Ov [а]

a£Ox[a] j£Oy [a]

(1)

+

Wvy[a,ß,7]V[ß]Y[7]+ £ Wxvy[a,a,ß,7]X[a]V[ß]Y[7]

ß&Ov [a] j£Oy [a]

a£Ox[a] ß&Ov [a] Y&Oy [a]

где X[a] € Rn , V[a] € Rm , Y[a] € Rl - состояние, управление и информация в узле системы; Wx[a, а] € Rcxn , Wv[a,ß] € Rcxm , Wy[a,7] € Rcxl, Wxy[a,a,7] € Rcxnxl , Wxv[a,a,ß] € Rcxnxm ,

wvy [a, ß, 7] € R

v

•cxmxl

W

xvy

[a, а, ß, 7] € R

cnml

матрицы-параметры; Ox [a] , Ov [a] , Oy[a]

- окрестности узла а по состоянию, управлению и информации соответственно; а, а, € € А, А = {а1, а2,..., ап} - конечное множество узлов системы. Модель (1) с подробной записью принимает вид:

'тх [1,1]х[1] + ых\1, 2]х[2] + ш [1, 1]«[1] + ш [1, 2]« [2]+ + шу[1,1]у[1] + ^ [1, 1, 1]х[1]«[1] + ш™[1,1, 2]х[1Ц2] + + тху [1, 2,1]х[2]«[1] + тху [1, 2, 2]х[2]«[2] + Шху [1,1,1]х[1]у[1]+ + Шху [1, 2,1]х[2]у[1]+ тУу [1,1,1]«[1]у[1] + тУу [1, 2,1Ц2]у[1]+ + тхУу [1, 1,1,1]х[1М1]у[1] + тхУу[1,1, 2, 1]х[1Ц2]у[1] + + шху [1, 2,1,1]х[2]«[1]у[1] + шх^у[1, 2, 2,1]х[2]«[2]у[1] = 0; шх[2,1]х[1] + тх[2, 2]х[2] + и>х[2, 3]х[3] + ш[2,1]«[1] + + ш [2, 2]«[2] + ш [2, 3]«[3] + шу [2, 2]у[2] + ш™ [2,1,1]х[1]«[1] + + тхУ [2,1, 2]х[1]«[2] + шхУ [2,1, 3]х[1]«[3] + ^ [2, 2,1]х[2]«[1] + + тху [2, 2, 2]х[2]«[2] + ^ [2, 2, 3]х[2]«[3] + ^[2, 3,1]х[3]«[1] + + тху [2, 3, 2]х[3]«[2] + ^ [2, 3, 3]х[3]«[3] + ш^у [2,1, 2]х[1]у[2]+ + шху[2, 2, 2]х[2]у[2] + шху [2, 3, 2]х[3]у[2] +

+ тУу [2,1, 2]«[1]у[2] + шУу [2, 2, 2Ц2]у[2] + ю„у [2, 3, 2Ц3]у[2]+ (2)

+ тхуу [2,1,1, 2]х[1Ц1]у[2]+ ^у [2,1, 2, 2]х[1Ц2]у[2] + + тхуу[2,1, 3, 2]х[1Ц3]у[2] + ^у[2, 2,1, 2]х[2]«[1]у[2] + + тхуу [2, 2, 2, 2]х[2Ц2]у[2] + ^у[2, 2, 3, 2]х[2Ц3]у[2] + + шхъу [2, 3,1, 2]х[3Ц1]у[2] + шхъу[2, 3, 2, 2]х[3Ц2]у[2] + + тхуу[2, 3, 3, 2]х[3]«[3]у[2]=0; шх [3, 2]х[2] + шх [3, 3]х[3] + ш [3, 2]«[2] + ш [3, 3]«[3]+ + шу [3, 3]у [3] + тху [3, 2, 2]х[2]«[2] + ^ [3, 2, 3]х[2]«[3] + + шху [3, 3, 2]х[3]«[2] + шху [3, 3, 3]х[3]«[3]+ + шху [3, 2, 3]х[2]у [3] + шху [3, 3, 3]х[3]у [3] + + тУу [3, 2, 3]«[2]у[3] + шУу [3, 3, 3]«[3]у[3]+ + тхУу [3, 2, 2, 3]х[2Ц2]у[3] + ^у[3, 2, 3, 3]х[2Ц3]у[3] + + тхУу [3, 3, 2, 2]х[3Ц2]у[3] + ^у [3, 3, 3, 3]х[3]«[3]у[3] = 0.

Для модели были выделены существенные компоненты состояния х , управления V и третьего параметра, в качестве которого может быть принята информация у (табл. 1).

Таблица 1.

Компоненты состояния, управления и информации.

х[1] Температура смотки для первой зоны участка полосы, °С

х[2] Температура смотки для второй зоны участка полосы, °С

х[3] Температура смотки для третьей зоны участка полосы, °С

v[1] Количество включенных полусекций охлаждения водой для первой зоны участка полосы, шт.

v[2] Количество включенных полусекций охлаждения водой для второй зоны участка полосы, шт.

v[3] Количество включенных полусекций охлаждения водой для третьей зоны участка полосы, шт.

у[1] Температура охлаждающей воды, ° С

у [2] Температура охлаждающей воды, ° С

у[3] Температура охлаждающей воды, ° С

Зададим значения компонентов состояния, управления и информации и проведем идентификацию трилинейной окрестной модели по разработанным ранее алгоритмам [2]. Значения состояния в соответствии с технологическими параметрами (рис. 1): х[1] =587 ° С; х[2] =573 ° С; х[3] =587 ° С.

Рис. 1: Пример распределения температуры смотки по зонам участка полосы

Значения управления в соответствии с технологическими параметрами: v[1] =15 шт.; v[2] =28 шт.; v[3] =18 шт.

Значения информации в соответствии с технологическими параметрами: у[1] = у[2] = = у[3] =40 ° С. В связи с разным порядком входных данных производим их нормализацию по формуле:

х = х-х,

а

где х - нормализуемое значение, х - среднее арифметическое, а - среднеквадратическое отклонение значений.

После нормализации получаем: х[1] = 0, 70711; х[2] = —1, 41421; х[3]=0, 70711; v[1] = —0, 95962 ; v[2] = 1, 37945 ; v[3] = —0, 41983 ; у[1]= у [2] = у[3] = 1.

Примем, что шу [1,1] = шу [2, 2] = шу [3, 3] = 1 .В результате идентификации получаем:

wx [1,1]= -0,029; wv [1,2]= -0,056; wxv [1,2,1]= -0,055; wxy [1,2,1]= 0,057; wxvy [1,1,1,1]= 0,027; wxvy [1,2,2,1]= 0,079; wx [2,3]= -0,024; wv [2,3]= 0,014; wxv [2,1,3]= 0,009896; wxv [2,2,3]= -0,02; wxv [2,3,3]= 0,009896; wxy [2,3,2]= -0,024; wvy [2,3,2]= 0,014; wxvy [2,1,3,2]= 0,009896; wxvy [2,2,3,2]= -0,02; wxvy [2,3,3,2]= 0,009893; wv [3,2]= -0,071; wxv [3,2,3]= -0,03; wxy [3,2,3]= 0,072; wvy [3,3,3]= 0,021; wxvy [3,3,2,3]= -0,05;

wx [1,2]= 0,057; wxv [1,1,1]= 0,027; wxv [1,2,2]= 0,079; wvy [1,1,1]= 0,039; wxvy [1,1,2,1]= -0,039; wx [2,1]= -0,024; wv [2,1]= 0,032; wxv [2,1,1]= 0,023; wxv [2,2,1]= -0,045; wxv [2,3,1]= 0,023; wxy [2,1,2]= -0,024; wvy [2,1,2]= 0,032; wxvy [2,1,1,2]= 0,023; wxvy [2,2,1,2]= -0,045; wxvy [2,3,1,2]= 0,023; wx [3,2]= 0,072; wv [3,3]= 0,021; wxv [3,3,2]= -0,05; wxy [3,3,3]= -0,036; wxvy [3,2,2,3]= 0,1; wxvy [3,3,3,3]= 0,015.

wv [1,1]= 0,039; wxv [1,1,2]= -0,039; wxy [1,1,1]= -0,029; wvy [1,2,1]= -0,056; wxvy [1,2,1,1]= -0,055; wx [2,2]= 0,047; wv [2,2]= -0,046; wxv [2,1,2]= -0,033; wxv [2,2,2]= -0,065; wxv [2,3,2]= -0,033; wxy [2,2,2]= 0,047; wvy [2,2,2]= -0,046; wxvy [2,1,2,2]= -0,033; wxvy [2,2,2,2]= 0,065; wxvy [2,3,2,2]= -0,033; wx [3,3]= -0,036; wxv [3,2,2]= 0,1; wxv [3,3,3]= 0,015; wvy [3,2,3]= -0,071; wxvy [3,2,3,3]= -0,03;

3. Параметрическое представление

Выразим окрестностные переменные x , v и y через параметрическую переменную U . В качестве примера возьмем представление компонента управления v . Укажем границы допустимых значений компонента управления третьего узла:

vmin — 0, vmax — 80 •

Примем, что значение U изменяется от Ui до U2 , которые взяты константами.

Таблица 2

Пример параметрического представления компонента v

U v v нормализ.

-1000 0 -3,65854

-900 4 -2,93882

-800 8 -2,21911

-700 12 -1,4994

-600 16 -0,77969

-500 20 -0,05998

-400 24 0,659736

-300 28 1,379448

-200 32 2,099161

-100 36 2,818873

0 40 3,538585

100 44 4,258297

200 48 4,978009

300 52 5,697721

400 56 6,417434

500 60 7,137146

600 64 7,856858

700 68 8,57657

800 72 9,296282

900 76 10,01599

1000 80 10,73571

Линейные зависимости представлены ниже:

х[1] = х[3] = 0.0015и + 1.1617; х[2] = 0.0015и — 1.8688;

v[1] = v[2] = v[3] = 0.0072и + 3.5386;

у[1] = у[2] = у[3] = 0.0005и + 1;

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

^ = —0.00285 • и + 1.11056 • 10"' • и2 + 2.5266 • 10"10 • и3 — 0.0259

4. Условие наличия экстремумов По аналогии с (3) рассмотрим функцию

Уа,ь{х) = х3 + а • х2 + Ь • х

первая и вторая производные которой имеют вид:

<1У 2

—— = 3х + 2ах + Ь <х

<2У <х2

= 6х + 2а

Исключая х из уравнений (6) и (7), получим уравнение параболы (рис. 2):

Ь =3 а2

получаем

(3)

(4)

(5)

(6)

(7)

Рис. 2: Поверхность параболы Ь = 3а2

Определим вид экстремумов выражения (4) в зависимости от условия (7). Для этого найдем критические точки, решив уравнение (5).

Б = (2а)2 — 4 • 3 • Ь = 4а2 — 12Ь;

тт —2а л/4а2 — 12Ь

и1,2 = -—г ±

23

23

Далее возможны два варианта развития событий: 1

3

1) при условии Ь> — а2 , л/4а2 — 126 < 0 - действительных корней нет, т. е. нет критических

точек;

2) при условии b <—a2 , \/ 4a2 — 12b > 0 , т. е. имеются критические точки: 3

Ui,2 = —1 a ± 3\/a2 — 3b Определим характер критических точек, подставив их значения в уравнение (6):

Ui = —31 a + a2 — 3b,

6U + 2a = 6[ —1 a + a2 — 3b] + 2a = —2a + 2\/a2 — 3b + 2a = 2^a2 — 3b > 0, 33

- критическая точка Ui является точкой минимума;

Ui = —1 a — \ a2 — 3b,

ви + 2а = 6[—1 а — —\/а2 — 3Ь] + 2а = —2а — 2л/а2 — 3Ь + 2а = — 2л/а2—3Ь < 0,

3 3

- критическая точка и2 является точкой максимума.

Таким образом, при условии Ь<-а2 (темная зона на рис. 2), функция (4) имеет две крити-

3—

ческие точки типа минимум и максимум, а при Ь> —а2 (светлая зона на рис. 2) критических

3

точек нет.

5. Результаты

Приведем общее параметрическое уравнение (3) к стандартному виду:

2 = и3 + 439.547и2 — 1.13 • 10"7 — 1.025 • 108 (8)

Выполняется условие Ь< 1 а2 , т.е. —1.13 • 10_7 < 1(439.547)2 . Поэтому функция (8) име-

33 ет одну точку минимума и = 1799, 7 и одну точку максимума и = —2092, 8 . При обратном

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

Подстановка найденных значений параметров в уравнения указывает некоторое отклонение от правой части. Минимизация же модуля общего параметрического уравнения (3) приводит к отклонению 6, 4 • 10_9 .

При этом получаем точку минимума и = —9.07, которой соответствуют окрестностные переменные х[1] = х[3] = 590° С, х[2] = 570° С, V = 40 шт., у[1] = у[2] = у[3]=39, 8° С, где V = 0.25v[1] +0^[2] + 0.25v[3] - среднее по длине участка полосы количество включенных полусекций охлаждения водой.

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

соответствии с [5], говорит о потере положения равновесия и возможности перехода системы в новое состояние.

6. Заключение

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

СПИСОК ЛИТЕРАТУРЫ

1. Shmyrin A. M., Sedykh I. A Identification and Control Algorithms of Functioning for Neighborhood Systems Based on Petri Nets // Automation and Remote Control. 2010. V. 71. Issue 6. P. 1265-1274.

2. Блюмин С.Л., Шмырин А.М., Шмырина О.А. Билинейные окрестностные системы: монография. Липецк: ЛГТУ, 2006.

3. Коновалов Ю.В. Справочник прокатчика. Книга 1. Производство горячекатанных листов и полос. М.: Теплотехник, 2008.

4. Mukhin U.A., Soloviev V.N., Makarov E.V. Improving of Low-Carbon Hot-Rolled Strips Production at the High Acceleration Rolling // New Technologies and Achievements in Metallurgy and Materials Engineering: a collective monograph. Czestochowa. 2012. P. 423-427.

5. Гилмор Р. Прикладная теория катастроф: В 2-х книгах. Кн. 1. Пер. с англ. М.: Мир, 1984.

БЛАГОДАРНОСТИ: Работа поддержана грантом РФФИ (код проекта 16-07-00854 а).

Поступила в редакцию 21 марта 2016 г.

Шмырин Анатолий Михайлович, Липецкий государственный технический университет, г. Липецк, Российская Федерация, доктор технических наук, профессор, заведующий кафедрой высшей математики, e-mail: amsh@lipetsk.ru

Ярцев Алексей Геннадьевич, Липецкий государственный технический университет, г. Липецк, Российская Федерация, магистрант физико-технологического факультета, e-mail: yartsevlekha@mail.ru

Правильникова Валерия Владимировна, Липецкий государственный технический университет», г. Липецк, Российская Федерация, ассистент кафедры высшей математики, e-mail: pravilnik@mail.ru

UDC 519.854

DOI: 10.20310/1810-0198-2016-21-2-463-470

TRILINEAR NEIGHBORHOOD MODEL OF THE PROCESS OF FORMING THE TEMPERATURE OF HOT-ROLLED STRIP COILING

© A.M. Shmyrin, A. G. Yartsev, V. V. Pravilnikova

The trilinear neighborhood model of the process of forming the temperature of hot-rolled strip coiling where the parameters are a state, control and information is considered. The purpose of the work is to find the values of the components of the model ensuring steady functioning of the system. The technique of defining the structure of the extrema is presented. The existence condition for extrema is received and checked on a concrete example. The suggestion about the area in which it is impossible to speak with definiteness about stability of the system is made. The hypothesis about a condition of the stable balance loss and of the transition of system to a new state is proposed.

Key words: neighborhood system, hot-rolled strip, run-off table, coiling temperature, general parametric equation.

ACKNOWLEDGEMENTS: This research was supported by the Russian Fund for Basic Research, project code 16-07-00854 a.

REFERENCES

1. Shmyrin A. M., Sedykh I. A Identification and Control Algorithms of Functioning for Neighborhood Systems Based on Petri Nets // Automation and Remote Control. 2010. V. 71. Issue 6. P. 1265-1274.

2. Blyumin S.L., Shmyrin A.M., Shmyrina O.A. Bilineynyie okrestnostnyie sistemyi: monografiya. Lipetsk: LGTU, 2006.

3. Konovalov Yu.V. Spravochnik prokatchika. Kniga 1. Proizvodstvo goryachekatannyih listov i polos. M.: Teplotehnik, 2008.

4. Mukhin U.A., Soloviev V.N., Makarov E.V. Improving of Low-Carbon Hot-Rolled Strips Production at the High Acceleration Rolling // New Technologies and Achievements in Metallurgy and Materials Engineering: a collective monograph. Czestochowa. 2012. P. 423-427.

5. Gilmor R. Prikladnaya teoriya katastrof: V 2-h knigah. Kn. 1. Per. s angl. M.: Mir, 1984.

Received 21 March 2016.

Shmyirin Anatoliy Mihaylovich, Lipetsk State Technical University, Lipetsk, the Russian Federation, Doctor of Techniques, Professor, the Head of the Higher Mathematics Department, e-mail: amsh@lipetsk.ru

Yartsev Aleksey Gennadevich, Lipetsk State Technical University, Lipetsk, the Russian Federation, Master program student, Faculty of Physics and Technology, e-mail: yartsevlekha@mail.ru

Pravilnikova Valeriya Vladimirovna, Lipetsk State Technical University, Lipetsk, the Russian Federation, Lecturer of the Higher Mathematics Department, e-mail: pravilnik@mail.ru