с дальнейшим построением соответствующей программы контроллера.
Для работы данной БД ЗСУ ПТ потребуется запуск среды МА^АВ. Это, с одной стороны, является недостатком подсистемы, поскольку она не работает без МА^АВ, а с другой - делает БД ЗCУ ПТ абсолютно открытой.
БД ЗСУ ПТ предполагает наличие графического диалогового интерфейса, который позволит вводить исходные данные и отображать результаты расчетов в удобном для пользователя виде (рис. 3). Поскольку направление ветра определяется положением подвижного носителя датчиков измерения, то в дальнейшем рассматривается только составляющая скорости ветра по оси факела.
В таблице приведены примеры отдельных реализаций метеорологических данных, полученных из ССМИ и ЦМИ для построения БД ЗСУ ПТ.
В таблице приведены численные данные, соответствующие графику на рисунке 3, для иллюстрации формирования массивов в БД.
Таким образом, в статье проведен анализ современных подходов к решению задачи минимизации выбросов в режиме реального времени, предложены способы построения БД метеорологической информации в ЗСУ ПТ, приведены требования к построению БД метеорологической информации, предложены схема сбора, обработки и передачи метеорологической информации в ЗСУ ПТ и схема структуры БД метеорологической информации.
Литература
1. Худолей В.В., Ливанов Г.А. Проблемы загрязнения окружающей среды стойкими органическими соединениями, в частности диоксинами // IEHS'04: тр. Междунар. конф. СПб, 2004. С. 39-43.
2. Сольницев Р.И. Построение замкнутой системы «При-рода-Техногеника» // Открытое образование: тр. Междунар. науч.-технич. конф., Украина, Ялта, 2006. С. 404-408.
3. Сольницев Р.И. Вопросы построения замкнутой системы управления «Природа-Техногеника» ЗСУ ПТ // Изв. СПБГЭТУ «ЛЭТИ». 2009. № 7. С. 23-32.
4. Сольницев Р.И., Коршунов Г.И. Системы управления «Природа-Техногеника». СПб: Политехника, 2013. 205 c.
5. Пугачев В.С. Введение в теорию вероятностей. М.: Наука, 1968. 368 с.
6. Сольницев Р.И. Вычислительные машины в судовой гироскопии. Л.: Судостроение, 1977. 312 с.
7. National Hydro-meteorological Service National Center for Hydro-meteorological Forecasting (NCHMF). URL: http://www. thoitietvietnam.gov.vn/web/en-US/62/19/58/map/Default.aspx (дата обращения: 26.11.2012).
8. Solnitsev R. Creation of «Nature-Technogenic» control system of the basis of information technologies, Proc. Intern. Conf. «IEHS'02», 2002, pp. 12-17.
9. Сольницев Р.И., Коршунов Г.И., Шабалов А.А. Моделирование замкнутой системы управления «Природа-Техноге-ника» // Информационно-управляющие системы. 2008. № 2 (33). С. 36-42.
References
1. Khudoley V.V., Livanov G.A., Trudy mezhdunarodnoy konf. IEHS'04 [Proc. of Int.Conf. IEHS'04], ISA, СПб, 2004, pp. 39-43.
2. Solnitsev R.I., Otkrytoe obrazovanie [Open Education], IT+S&E'06, Ukraine, Yalta, 2006, pp. 404-408.
3. Solnitsev R.I., Izvestiya SPBGETU «LETI» [Proc. of St. Petersburg ETU «LETI»], 2009, no. 7, pp. 23-32.
4. Solnitsev R.I., Korshunov G.I., Sistemy upravleniya «Pri-roda-Tekhnogenika» [Managing Systems «Nature-Technogenica»], 2013, 205 p.
5. Pugachyov V.S., Vvedenie v teoriyu veroyatnostey [An introduction to the theory of chances], Moscow, 1968, 368 p.
6. Solnitsev R.I., Vychislitelnye mashiny v sudovoy giroskopii [Calculating machines in marine gyroscopy], Leningrad, 1977.
7. National Hydro-meteorological Service National Center for Hydro-meteorological Forecasting (NCHMF), available at: www.thoitietvietnam.gov.vn/web/en-US/62/19/58/map/Default.aspx (accessed 25 Dec. 2012).
8. Solnitsev R.I., Proc. Intern. Conf. IEHS '02, St. Petersburg, 2002, pp. 12-17.
9. Solnitsev R.I., Korshunov G.I., Shabalov A.A., Information and Control Systems, 2008, no. 2 (33), pp. 36-42.
Сохранение метеорологической информации
День/месяц I II III IV V VI VII VIII
1 83 81 58 90 83 88 87 91
2 84 83 68 81 77 80 90 86
3 82 83 73 81 74 86 84 82
4 81 84 75 80 77 87 81 86
5 82 83 94 80 76 77 86 83
6 67 93 94 76 75 72 81 99
УДК 681.142.2 +621.746
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТЕПЛОФИЗИЧЕСКИХ ПРОЦЕССОВ В СЛОИСТЫХ СТРУКТУРАХ
Г.А. Дубский, к.ф.-м.н., доцент; Л.Г. Егорова, к.т.н.., доцент; Ю.Б. Кухта, к.т.н., доцент (Магнитогорский государственный технический университет им. Г.И. Носова, пр. Ленина, 38, г. Магнитогорск, Челябинская обл., 455000, Россия, egorov-lyudmil@yandex.ru)
Одним из основных факторов, способствующих повышению эффективности функционирования агрегатов непрерывного литья заготовок, является обеспечение работоспособности и оптимальных технологических режимов работы оборудования. Для увеличения срока службы кристаллизатора машин непрерывного литья заготовок в
настоящее время прибегают к определенным технологическим модификациям их рабочих стенок. Это улучшение качества поверхности кристаллизатора за счет использования специальных сплавов либо нанесение покрытий на поверхность медной стенки. В работе представлено физико-математическое моделирование теплофизического процесса системы «расплав металла-покрытие-стенка кристаллизатора» в момент заливки стали. Для математического моделирования разработан программный продукт, использование которого позволяет проследить характер установления температурного поля при тепловом ударе в момент заливки стали, проанализировать температурное поле в зависимости от времени в покрытии и стенке кристаллизатора и определить максимальную величину температуры. Результаты, полученные в ходе использования программного продукта, были использованы при выборе оптимальных режимов нанесения защитно-восстановительного покрытия из никеля на внутреннюю стенку кристаллизатора машины непрерывного литья заготовок, способствующих повышению износоустойчивости оборудования без снижения качества готовой продукции.
Ключевые слова: температурное поле, модифицирующее покрытие, стенка кристаллизатора, программный продукт, математическое моделирование.
MATHEMATICAL MODELLING OF THERMAL PROCESSES IN LAYERED STRUCTURES Dubsky G.A., Ph.D., Associate Professor; Egorova L.G., Ph.D., Associate Professor;
Kukhta Yu.B., Ph.D., Associate Professor (G.I. Nosova Magnitogorsk State Technical University, Lenina Av., 38, Magnitogorsk, Chelyabinsk Reg., 455000, Russia, egorov-lyudmil@yandex.ru)
Abstract. One of the major factors promoting increase of efficiency of functioning of units of continuous molding of preparations, ensuring working capacity and optimum technological operating modes of the equipment is. For increase in service life of crystallizers now resort to certain technological modifications of their working walls. This improvement of quality of a surface of a crystallizer at the expense of use of special alloys or coverings on a surface of a copper wall. Physical and mathematical simulation of thermal process of system «metal melt-a coating-a crystallizer pan wall» at the moment of steel priming is in-process presented. The software product which use allows to track character of an establishment of a temperature field at a heatstroke at the moment of steel pouring is developed for mathematical modelling; to analyse a temperature field from time in a covering and a wall of a crystallizer and to define the maximum size of temperature. The results received during use of the software product, were used at a choice of optimum modes of a covering from nickel on an internal wall of a crystallizer for increase of wear resistance of the equipment without decrease in quality of finished goods.
Keywords: the temperature pattern, an inoculating coating, a crystallizer pan wall, software product, mathematical modelling.
При производстве кристаллизаторов для сортовых, блюмовых и слябовых машин непрерывного литья заготовок (МНЛЗ) используют толстолистовую медь, из которой изготавливают плиты необходимых размеров с последующей доработкой их рабочих поверхностей до требуемых технологических и технических условий. Кристаллизаторы, собранные из таких листов, как правило, показывают низкую производительность из-за быстрого износа рабочих стенок, что сокращает время межремонтных периодов, а значит, увеличивает текущие расходы, повышая тем самым себестоимость единицы продукции.
Во время разливки температура на поверхности меди достигает значений температуры ее размягчения или даже температуры рекристаллизации, что приводит к разрушению поверхности. В результате на поверхности медных кристаллизаторов появляются и развиваются трещины. Прилипая к поверхности сляба, медь вызывает разветвленные паукообразные трещины.
Для увеличения срока службы кристаллизаторов в настоящее время все чаще прибегают к определенным технологическим модификациям их рабочих стенок - улучшение качества поверхности кристаллизатора за счет использования специальных сплавов либо нанесение покрытий на поверхность медной стенки. Физические свойства никеля (твердость, теплоемкость, теплопроводность) позволяют успешно использовать его для восстановления медных плит кристаллизаторов. В
то же время медные плиты, изготовленные из специальных сплавов, очень дороги и требуют значительных затрат на ремонт. Благодаря этим факторам восстановление работоспособности медных плит с использованием никелевого покрытия экономически выгодно. Опыт применения кристаллизаторов с покрытиями показывает увеличение их срока службы до 300 000-800 000 тонн стали с возможностью последующего восстановления покрытия.
Прочность сцепления покрытия со стенкой кристаллизатора во многом зависит от характера поведения температурного поля в покрытии и на границе «покрытие-стенка кристаллизатора». Как показывают экспериментальные исследования [1, 2], температуропроводность, теплопроводность, теплоемкость напыленных покрытий очень сильно зависят от температуры и характера ее применения в пространстве и времени. Целью данной работы являются математическое моделирование теплофизического процесса системы «расплав металла-покрытие-стенка кристаллизатора» в момент заливки стали (тепловой удар) и разработка программного обеспечения для решения данной задачи.
Заполненный расплавом стали кристаллизатор МНЛЗ с введенной в него затравкой представляет собой сложную термодинамическую систему, в которой в пространстве и времени развиваются сложные теплофизические процессы, приводящие к формированию «корки» слитка, термонапряже-
нии в покрытии и медной стенке, на границе «покрытие-стенка». Прочность сцепления покрытия со стенкой кристаллизатора во многом зависит от характера поведения температурного поля в покрытии и на границе «покрытие-стенка кристаллизатора». Схема описываемой термодинамической системы приведена на рисунке 1.
I Расплав стали ©01 II qi —► Покрытие ©1 а1 III q2 Г\ —► О Стенка 0- ©2 , a1 Q ©02, ^2 О --Вода -—-—
0 l L x Рис. 1. Термодинамическая система «расплав металла-покрытие-стенка кристаллизатора»
Для описания этой термодинамической системы воспользуемся нестационарным уравнением теплопроводности, учитывающим сильную зависимость теплопроводности материалов покрытия и стенки от температуры, и выбранными краевыми условиями. Для двухслойной задачи математическая модель может быть представлена в следующем виде:
аэдхт)
с,-р,- —-= к i ©——+
' 5z i dx2
+ dkL fd@i (xt, z) '
a© t ax (i=1, 2), ©i (x, 0) = ©01,
©2( x,0) =©o2, © (I, z) = ©2 (I, z) ,
a©.(/, z) a©2(/, z)
4l - —-- - Л2 -Z- - 42 ,
(1)
(2)
(3)
(4)
(5)
дх дх
\ = [1 - кх © (х, т) + к2 © (х, z)], (6)
Х2 =Хо2(1- х, т)) , (7)
где ©01 - температура расплава; ©02 - температура охлаждающей жидкости; а1 - температуропроводность покрытия; а2 - температуропроводность материала стенки кристаллизатора; ©1(х, т) - переменная температура покрытия; ©2(х, т) - переменная температура стенки кристаллизатора; Х01 -теплопроводность покрытия при 0 °С; Х02 - теплопроводность материала стенки кристаллизатора при 0 °С; к1, к2, к3 - коэффициенты пропорциональности, определяемые из эксперимента. Толщина покрытия намного меньше толщины стенки
I
кристаллизатора, то есть ——- << 1.
Численное решение дифференциальных уравнений с заданными краевыми условиями, соответствующими тепловому удару в момент заливки стали, выполнялось с использованием неявной разностной схемы.
Для задачи моделирования теплового состояния в системе «защитное покрытие-рабочая стенка» кристаллизатора объектом исследования является температурное поле покрытия и рабочей стенки, изменяющееся в пространстве и времени. По результатам математического моделирования получим зависимости температуры покрытия по толщине во времени, температуры границ покрытия с металлом и кристаллизатором, зависимости изменения температуры от времени. По этим зависимостям можно сделать вывод о поведении температурного поля, вынести решение о выборе параметров для расчета. Диаграмма функциональной модели изменения температурного поля в покрытии и рабочей стенке кристаллизатора представлена на рисунке 2.
Для математического моделирования изменения теплофизических процессов в системе «покрытие-стенка кристаллизатора» при различных начальных исходных данных был разработан программный продукт, который позволяет выполнить
- ввод и изменение исходных данных для коррекции модели;
- математический расчет результатов моделирования;
- графическое отображение результатов математического моделирования.
Для реализации программного продукта были выбраны язык программирования C++ и интегрированная среда программирования Borland C++ Builder 6.0.
В программном продукте можно выделить три основных модуля: подготовка исходных данных, реализация математической модели теплового состояния покрытия, графическая визуализация данных.
Результатами работы модуля реализации математической модели в автоматическом режиме являются значения температур, рассчитанных с учетом толщины покрытия и стенки кристаллизатора и количества интервалов разбиения. Для расчета значений температуры в рабочей стенке кристаллизатора входными значениями температуры являются значения температур, получаемые на границе покрытия заданной толщины.
Модуль графической визуализации данных предназначен для отображения результатов реализации математической модели. По результатам математического моделирования получены:
- зависимость температуры покрытия по толщине и времени;
- поверхность температур в координатах (l, ©, z) для точек по толщине покрытия;
Цели и задачи исследования
Алгоритмы и методы моделирования, обработки экспериментальных данных
Требования ГОСТов, нормативы
1
i
Данные об объекте исследования
Данные для расчета
Сбор информации
/
Персонал производственных цехов и научно-исследовательских лабораторий
Математическое моделирование теплового состояния покрытия и рабочей стенки кристаллизатора
Значения
рассчитанных
температур
Обработка результатов моделирования
Информация о новых параметрах
Оборудование производственных цехов и лабораторий
Информация о формировании температурных полей
Решение задачи
оптимизации температурного поля
Корректировка исходных данных для расчета
Решения по поведению температурного поля
Закономерности при тепловом ударе
Средства вычислительной техники и программное обеспечение
Рис. 2. Диаграмма функциональной модели исследования температурного поля в системе «покрытие-рабочая стенка» кристаллизатора
- поверхность линий уровня для температурного поля во времени и пространстве;
- зависимость температуры границы покрытия с металлом и стенкой кристаллизатора;
- зависимость температуры системы «покрытие-подложка» по толщине (рис. 3).
Использование представленного программного продукта позволило проследить характер установления температурного поля в системе «покрытие-подложка» при тепловом ударе в момент заливки стали, проанализировать температурное поле в зависимости от времени в покрытии и стенке кристаллизатора и определить максимальную величину температуры.
Полученные результаты могут использоваться для выбора технологии напыления модифици-
рующих покрытий из никеля и плакированного никеля на рабочую стенку кристаллизатора, способных противостоять их возможному отслаиванию за счет возникающих касательных тепловых напряжений при тепловом ударе. Решение задач моделирования теплофизических процессов в слоистых структурах будет способствовать более интенсивной разработке новых эффективных видов покрытий, позволяющих удовлетворять различные области техники и металлургии.
Разработанные функциональные модели, математическое и информационное обеспечение найдут практическое применение при анализе и прогнозировании производственных ситуаций при непрерывном литье заготовок, решении технологических задач по управлению технологией разливки стали в изложницы, что повысит показатели производительности, качества продукции и экономической эффективности.
Литература
1. Вдовин К.Н., Дубский Г.А., Нефедьев А.А. Теплофи-зические свойства детонационно-напыленного никеля на медные пластины кристаллизаторов МНЛЗ // Вестн. МГТУ им. Г.И. Носова. Магнитогорск. 2008. № 2. С. 57-60.
2. Дубский Г.А., Бутаков С.А., Нефедьев А.А. Теплофи-зические свойства детонационно-напыленной меди при заданных технологических условиях // 66-я науч.-технич. конф.: сб. докл. Магнитогорск: МГТУ, 2008. С. 118-121.
References
1. Vdovin K.N., Dubsky G.A., Nefedev A.A., Vestnik Magnitogorsk State Technical Univ., 2008, no. 2, pp. 57-60.
2. Dubsky G.A., Butakov S.A., Nefedev A.A., Sbornik dokladov 66 konf. Magnitogorskogo Gos. Tekh. Univ. [Proc. 66th Conf.], Magnitogorsk, Magnitogorsk State Technical Univ., 2008, pp. 118-121.
Толщина покрытие-стенка кристаллизатора, мм
Рис. 3. Зависимость температуры системы «покрытие-подложка» кристаллизатора МНЛЗ