CC BY
37
2. далее агент-координатор передает данные поисковому агенту;
3. поисковый агент считывая данные, ищет в базе прецедентов сочетание данных, если такое сочетание находится, то выдается готовое решение;
4. если не находится, то идет обращение к базе знаний. В базе знаний хранятся правила для создания опытных решений. Поисковый агент находит нужные правила, и передает их в блок вычисления, для получения ответа, затем полученное опытное решение передается в базу прецедентов.
Алгоритм работы многоагентной системы показан на рис.2.
Список литературы:
1. Туманов В.Е. Применение многоагентных технологий при разработке встроенных в портал экспертных систем // Программные продукты и системы, №3, 2009.
2. Гречишников В.А. Моделирование систем инструментального обеспечения автоматизированных производств. М: ВНИЦТЭМР, 1988, 450 с.
3. FIPA Agent Management Specification. 1996-2002 Foundation for Intelligent Physical Agents http://www.fipa.org/
List of references:
1. Tumanov V.E. Use of multi agents technologies for development of the built in portal expert systems // Software and Systems №3, 2009.
2. Grechishnikov V.A. Simulation systems provide automated production tool. M: VNICTEMP, 1988, 450 p.
3. FIPA Agent Management Specification. 1996-2002 Foundation for Intelligent Physical Agents http://www.fipa.org/
А.В. Польщиков
аспирант, ФГБОУ ВПО «Магнитогорский государственный технический
университет им. Г. И. Носова»
В.Д. Тутарова
к. техн. н., доцент, ФГБОУ ВПО «<Магнитогорский государственный технический
университет им. Г. И. Носова»
ПОСТАНОВКА ЗАДАЧИ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ ЗАТВЕРДЕВАНИЯ НЕПРЕРЫВНОЛИТОЙ КРУГЛОЙ ЗАГОТОВКИ
Аннотация: в статье рассматривается математическая модель затвердевания непрерывнолитой заготовки круглого сечения в зоне вторичного охлаждения с учетом конструкционных особенностей машин непрерывного литья заготовок №1 ОАО «Уральская сталь» (г. Новотроицк).
Ключевые слова: математическое моделирование, машины непрерывного литья заготовок, зона вторичного охлаждения, уравнение теплопроводности.
A.V. Polschikov, Magnitogorsk state technical university
V.D. Tutarova, Magnitogorsk state technical university
FORMULATION OF THE PROBLEM OF MATHEMATICAL MODELING OF SOLIDIFICATION
OF CONTINUOUS CASTING ROUND BILLETS
Abstract: in the article the mathematical model of solidification billets of circular cross section i n the zone of secondary cooling taking into account the structural features of CCM № 1 of «Ural Stee» (Novotroitsk).
Keywords: mathematical modeling, information technology, the zone of secondary cooling, the mold, heat equation.
В последнее десятилетие производство непрерывнолитых круглых заготовок получило большое развитие в промышленно развитых странах и к 2010 году составило около 10,1 млн. тонн, что на 50,3% выше, чем в 2006 г. В связи с этим идет непрерывный поиск лучших конструкций машин непрерывного литья заготовок (МНЛЗ) и совершенствуются отдельные узлы существующих установок для повышения производительности на уже действующих МНЛЗ для производства заготовок круглого сечения. Решение данных задач весьма затруднительно без использования инструментов математического моделирования.
Каждая МНЛЗ имеет свои конструкционные особенности, что необходимо учитывать при разработке математической модели затвердевания непрерывнолитой заготовки. В качестве объекта моделирования была выбрана четырехручьевая комбинированная МНЛЗ №1, установленная в электросталеплавильном цехе ОАО «Уральская сталь» (г. Новотроицк). Это уникальный агрегат, позволяющий разливать трубную заготовку от 300 до 600 мм в диаметре, которую в России делают только на прокатных станах.
Математическая модель затвердевания непрерывнолитой круглой заготовки была построена на основе следующих общих положений.
При нестационарных условиях разливки наблюдаются переходные режимы в температурном поле заготовки, которые статические модели расчета температуры поперечного сечения заготовок во времени не позволяют отследить. Для этой цели предлагается использовать динамическую модель охлаждения заготовок, которая позволяет определять температуру по всему объему заготовок во времени.
Процессы теплообмена во всех фазах заготовки описываются нестационарным уравнением теплопроводности с внутренними источниками тепла в цилиндрических координатах вида:
дТ „ Гд2Т 1 дТ 1 д2Т
0 < г < R, 0 <ç< 2п.
(1)
pc — = Л —г- +-----+ -г---г-
и дт дг2 г дг г2 dç2
где
Т = Т(г,ф,т) - функция распределения температур по сечению заготовки во
времени, 0С;
т - время, с;
с(т) -теплоемкость металла, Дж/(кг0С); Л(т) -теплопроводность металла, Вт /(м0С); Ф - угол по радиусу заготовки, рад. р(т) - плотность металла, кг/м3.
Известно, что при затвердевании бинарный сплав железа с углеродом претерпевает не только фазовый, но и структурный переход. Наибольшее тепловыделение происходит при фазовом переходе, но при точных расчетах следует учитывать и тепловыделение, возникающее при охлаждении за счет структурных переходов. Для учета выделяющейся теплоты структурных переходов вводится величина относительного количества вещества с определенной структурой у.
Относительное количество твердой фазы находящейся в равновесии с жидкостью при Т, может быть определено из диаграммы состояния сплава Ре-С. Если принять, что линии ликвидуса и солидуса - параллельные прямые, то
у = (тл - т)/(тл - Тс), (2)
где Тл - температура ликвидуса, 0С; Тс - температура солидуса, 0С.
Температура ликвидуса и солидуса для углеродистых и низколегированных сталей определяется с учетом влияния содержания в металле химических элементов согласно зависимостям:
Тл =1536-(78 [С]+ 7,6•[3'] + 4,9■[№]+ 4,7•[0/]+1,3[Сг] + 3,1[М] + 3,6[А/] + 38^[в] + 34,4[Р]) (3)
Тс = 1536 - (200 [С] +16 • [в/] + 6 • [Мп] +1,7 • [Сг] + 3,9 • [М ] +1100 [в] + 93 • [Р]), для низколегированных -
Тл = 1539 - (80 • [С] +1,09 • [Сг] +1,80 • [М] + 2,39 • [Мо] + 2,13 [V] + 45 • [в] + 27,8 • [Р]) (4)
ТС = 1539 - (80 • [С] +1,54 • [Сг] + 2,06 • [М] + 2,70 • [Мо] + 2,37 [V] + 870 [в] + 97,8 • [Р])
где [С], [в/], [Мп]. и др.- процентное содержание соответствующих химических элементов в стали.
Ю.А. Самойловичем был предложен следующий способ учета тепла кристаллизации. При постоянном темпе кристаллизации сплава внутри интервала температур Тл - Тс величина - ду/дТ = -[1/(Тл - Тс)].
При наличии теплообмена тела в различных точках имеют различную температуру. Для упрощения методики решения задач затвердевания теплоту кристаллизации и структурных превращений учитывают при помощи введения эффективных теплоемкости и теплопроводности металла сэ (Т) и Яэ (Т).
Коэффициенты сэ (Т) и Яэ (Т) определяются в зависимости от температуры исследования фазового состояния металла по формулам:
с
при Т > Тл
сэ(Т) = \с(Тс) + ^ при Тс <Т <Тл ,
при Т < Тс
при Т > Тл
Я (Т) = \Ят9 + Яж (1-р) при Тс < Т < Тл ,
при Т < Тс
где ^ - скрытая теплота кристаллизации, кДж/кг;
В результате преобразований уравнение теплопроводности примет вид:
с(Тс) + L ^ V с> дТ
Ст
Яж
ЯтР+Яж (1 -р)
ят
(5)
(6)
дТ д
рсэ ~эТ =дг
(
0 дТ 1 д2Т
гЯэ— +---2-
э бг г др
\
0 < г < Г1,0 <р< 2п.
(7)
Разность расходов охладителя по большому и малому радиусу машины можно учитывать несимметричностью граничных условий третьего рода:
дг
= аг ( - Тпов ),
(8)
где - радиус заготовки, м,
аг - суммарный коэффициент теплоотдачи с поверхности к заготовке
а = Г (г,Т,У, Т, ^), Вт/м 0С;
Тпов - температура охладителя, подаваемого на грани заготовки, 0С.
Особенностями существующей системы вторичного струйного охлаждения на МНЛЗ №1 ОАО «Уральская сталь» при затвердевании круглой заготовки является несимметричность охлаждения относительно большого и малого радиуса не только при подаче охладителя, но и при расположении заготовки относительно коллекторов сегментов охлаждения, что необходимо учесть в граничных условиях 3-его рода. На каждом коллекторе первого и второго сегментов располагается по 6 круглофакельных форсунок в шахматном порядке, относительно рядом находящихся коллекторов.
При решении задачи предлагается рассматривать половину круглой заготовки (рис.1). Расходы охладителей, подаваемые на коллектора верхней части сегмента I и нижней II в соответствии с технологической инструкцией отличаются. Для расчета
коэффициента теплоотдачи необходимо учитывать локальную площадь орошения заготовки от каждой форсунки.
Нечетный коллектор Четный коллектор
Рис. 1. Расположение форсунок на четном и нечетном коллекторах
Введем обозначения: R1 - радиус коллектора, мм, R2 - радиус заготовки, мм, с центром в точке О, симметрично расположенной относительно коллектора, r - радиус заготовки произвольного размера, мм с центром в точке О1, смещенной относительно центра О. В точке A нижняя граница радиусов заготовок совпадает.
Локальная площадь орошения зависит от расхода воды, подаваемой на форсунку foxn, расстояния от сопла форсунки до поверхности заготовки s, типа
форсунки (плоскофакельная, круглофакельная, конусная и др.), а также от ее геометричеких параметров (^ - угол раскрытия факела, dc - диаметр сопла форсунки)
Lk _ op = f ((, ^ dc ).
На рис. 1 и 2 жирным выделены расстояния от сопла форсунки до поверхности заготовки, которые необходимо рассчитать.
В соответствии с теоремой косинусов определяем:
|ОС| = |001 • cos ZC001 ^|00112 • cos2 ZC001 -100112 + |01С|2 |СС2| = |0С2| - \0С\.
Введем обозначения: |00^ = R2 - r, |0С| = r, |0С2| = R1, |0С.,| = R2 Следовательно, расстояния от сопла форсунки до поверхности заготовки вычисляются по формуле:
s = R1 -((R2 - r) • cos в + ^(R2 - r )2 • cos2 в- (R2 - r )2 + r2). (9)
L«_орх =п-1 S • tan£ I. (10)
Рис. 2. Геометрия расчета расстояний от форсунки до заготовки
Зная угол раскрытия факела форсунки, расстояние от сопла форсунки до поверхности заготовки и учитывая геометрию факела форсунки, можно определить площадь орошаемой поверхности.
Геометрия «пятна» орошаемой поверхности (рис. 3) для круглофакельной форсунки для заготовки круглого сечения принимает форму эллипса.
Локальный фронт орошения вдоль линии заготовки по направлению х составит: (
Следовательно, размер «пятна» вдоль оси х составит \K1K2 \ = 2 • fn0K орх.
На рис. 4 C2O1 - расстояние от сопла форсунки до центра заготовки, C2C-расстояние от сопла до поверхности заготовки, O1N1- радиус заготовки, известен угол раскрытия факела р.
Определим размер «пятна» вдоль оси у (рис.4).
Рассмотрим треугольник O1N1C2.
|N1C2|2 - 2 • |N1C2| • |N1O1| • cos ZN1C2O1 - |N1O1|2 = 0.
Решаем квадратное уравнение, относительно N,C2.
В соответствии с теоремой косинусов определяем:
_орУ = NC = VNC|2 + |CC212 -2 • |CC21 • \NC\ • cosZNAC . (11)
Распределение охладителя по площади орошения неодинаково. Основная часть его приходится на центральный участок «пятна». Для количественного отношения вводим следующее допущение: распределение охладителя происходит в соответствии с правилами 6а.
Следовательно, размер «пятна» вдоль оси у составит |М1М21 = 2 • Глок_ору.
Используемые конструкции форсунок характеризуются углом раскрытия 900 и малым размером проходных отверстий (диаметр меньше 1 мм), что в свою очередь может приводить к засорению форсунок примесями воды и выходу форсунок из строя, и, как следствие, к появлению дефектов в МНЛЗ. В виду последнего, очень часто их производительность снижается. Так лабораторные испытания, проведенные на лабораторном стенде в лаборатории ЗАО «Корад» продемонстрировали, что угол раскрытия используемых форсунок нестабилен в заявленном рабочем диапазоне расхода воды и заметно меняется.
Таким образом, коэффициент теплоотдачи в ЗВО, учитывающий несимметричность охлаждения определяем по формулам: в зоне плотного орошения
1570,0 • (( • ГлоК_ор ) [1 - 0,0075 • ( - 273,5)]
Рис.3. Геометрия пятна орошаемой поверхности
аг = -
(12)
между зонами плотного орошения и контакта с роликами
а (*4 - - 1С),
где е - приведенный коэффициент черноты тела, Ф - постоянная Стефана-Больцмана.
(13)
Рис. 4. Геометрия расчета размера «пятна» вдоль оси у
В начальный момент времени температурное поле принимается однородным: она во всем объеме постоянна и равна температуре металла в промежуточном ковше.
а
t = 0: t = ^ 0 < г < r . (14)
Данная модель позволит рассчитать температурное поле по сечению заготовки с учетом конструкционных особенностей МНЛЗ №1 ОАО «Уральская сталь» (г. Новотроицк).
Список литературы:
1. Процессы непрерывной разливки. Монография/ Смирнов А.Н., Пилюшенко В.Л., Минаев А.А. и др., Донецк: ДонНТУ, 2002 г.-536 с.
2. Тепловые процессы при непрерывном литье стали, под редакцией Ю.А.Самойловича.- М.: Металлургия, 1982 г.-152 с.
3. Соболев В.В, Трефилов П.М., Теплофизика затвердевания металла при непрерывном литье - М.: Металлургия, 1988 г.-160 с.
List of references:
1. Сontinuous casting process. Monograph / Smirnov, AN, Pilyushenko VL, Minaev, AA and others, Donetsk: Donetsk National Technical University, 2002 - 536.
2. Thermal processes in the continuous casting of steel, edited by A. Samoilovich .- Moscow: Metallurgiya, 1982 - 152.
3. Sobolev, Trefilov, Thermophysics of metal solidification in continuous casting -M.:Metallurgy, 1988 - 160.
В.П. Божко
д. экон. н., профессор, ГОУ ВПО "Московский государственный университет
экономики, статистики и информатики"
А.М. Батьковский
к. экон. н., ОАО «ЦНИИ «Электроника»
А.П. Мерзлякова ОАО «ЦНИИ «Электроника»
МОДЕЛЬ ФОРМИРОВАНИЯ ОПТИМАЛЬНОЙ ПРОГРАММЫ ИННОВАЦИОННОГО РАЗВИТИЯ ЭКОНОМИЧЕСКОЙ СИСТЕМЫ
Аннотация. Разработанная и представленная в статье модель обеспечивает возможность проводить многовариантные расчеты и формировать наиболее эффективную программу инновационного развития экономических систем.
Ключевые слова: инновационное развитие, оптимизация, модель, программа.
V.P. Bogko, Moscow State University of Economics, Statistics and Informatics A.M. Batkovsky, Central Research Institute of Electronics А.P. Merzlyakova, Central Research Institute of Electronics
