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

На правах рукописи

ВЕСЕЛОВ ЮРИЙ ВЛАДИМИРОВИЧ

МЕТОД И АЛГОРИТМЫ ОБРАБОТКИ ИНФОРМАЦИИ В СИСТЕМЕ ПРОГНОЗИРОВАНИЯ ОХЛАЖДЕНИЯ АГЛОМЕРАТА НА ПРЯМОЛИНЕЙНОМ ОХЛАДИТЕЛЕ

Специальность 05.13.01- Системный анализ, управление и обработка информации (в металлургии)

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

Череповец - 2009

003492601

Работа выполнена в Государственном образовательном учреждении высшего профессионального образования «Череповецкий государственный университет»

Научный руководитель: доктор технических наук, профессор

Кабаков Зотей Константинович

Официальные оппоненты: доктор технических наук, профессор

Аншелес Валерий Рудольфович

кандидат технических наук Бирюков Евгений Николаевич

Ведущая организация: ЗАО «Фирма «СТОИК», г. Череповец

Защита состоится 18 декабря 2009 г. в 16 часов на заседании диссертационного совета Д 212.297.02 в Государственном образовательном учреждении высшего профессионального образования «Череповецкий государственный университет» по адресу 162600. г. Череповец, Вологодская обл., ул. Луначарского, д.5, конференц-зал.

С диссертацией можно ознакомиться в библиотеке Череповецкого государственного университета

Автореферат разослан «17» ноября 2009 г.

Ученый секретарь диссертационного совета

кандидат технических наук, доцент ----К.А. Харахнин

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

После окончания процесса спекания шихты среднемассовая температура агломерата обычно находится в пределах 700^800°С. Высокая температура агломерата создает трудности при последующей его обработке на аглофабрике, транспортировке и загрузке в доменную печь. В связи с этим на современных аглофабриках агломерат охлаждают на охладителях различного типа, в том числе на прямолинейных, что облегчает условия работы обслуживающего персонала и оборудования. На практике известны случаи возгорания транспортерных лент при транспортировке агломерата в доменный цех.

Благодаря фундаментальным работам большого круга ученых, достигнуты значительные успехи в области совершенствования технологии производства агломерата. Наибольший вклад в развитие теории и технологии агломерационного процесса внесли работы ученых Бабушкина Н.М., Тимофеева В.Н., Коро-тича В.И., Пузанова В.П., Фролова Ю.А., Майзеля Г.М., Братчикова С.Г., Бази-левича C.B., Вегмана Е.Ф., Белоцерковского Я.Л., Бабошина В.М., Шкляра Ф.Р., Сигова A.A., Шурхала В.А. и др. Однако, практически все известные методы ориентированы на прогнозирование процесса спекания шихты. Поэтому в настоящее время разработка метода, позволяющего прогнозировать процесс охлаждения агломерата на прямолинейном охладителе в системе управления агломерационным производством, представляется актуальной научно-технической задачей.

Объект исследования: процесс охлаждения агломерата на прямолинейном охладителе.

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

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

Для достижения поставленной цели в работе решаются следующие основные задачи:

- анализ существующих математических моделей и методов прогнозирования хода технологических процессов агломерационного производства;

- разработка метода прогнозирования, основанного на математической модели процесса охлаждения агломерата на прямолинейном охладителе;

- разработка алгоритмов обработки информационных сигналов, позволяющих прогнозировать технологические параметры;

- проведение экспериментальных исследований работоспособности и эффективности предложенных метода и алгоритмов.

Методы исследований. Для решения поставленных задач в работе использовались теоретические основы газодинамики и теплотехники агломерационного процесса; методы статистической обработки информации; методы математического моделирования. Математическое моделирование производилось путем разработки численной модели для реализации в среде разработки Delphi 7. Достоверность полученных результатов подтверждена путем сравнения с экспериментальными данными.

Научная новизна и основные положения, выносимые на защиту:

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

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

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

Практическая ценность работы заключается в следующем:

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

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

Реализация результатов работы. Разработанные метод и алгоритмы обработки информационных сигналов, реализованные в системе прогнозирования процесса охлаждения агломерата, прошли опытно-промышленную эксплуатацию в агломерационном производстве на Череповецком металлургическом комбинате (ЧерМК) ОАО «Северсталь».

Теоретические результаты внедрены в учебный процесс на кафедре программного обеспечения ЭВМ ГОУ ВПО «Череповецкий государственный университет» по дисциплине «Основы теории управления».

Апробация работы. Основные положения диссертационной работы докладывались и получили положительную оценку на 4-ой международной научно-технической конференции «Автоматизация и энергосбережение машиностроительного и металлургического производств, технология и надежность машин, приборов и оборудования» (Вологда, 2008), Всероссийской научной конференции «Молодые исследователи - регионам» (Вологда, 2009) и на постоянно действующем научно-техническом семинаре кафедры программное обеспечение ЭВМ Череповецкого государственного университета.

Публикации. По материалам диссертации опубликовано 5 печатных работ, в том числе 2 статьи опубликованы в рецензируемых научных журналах,

входящих в перечень ВАК Министерства образования и науки Российской Федерации.

Структура н объем работы. Диссертационная работа состоит из введения, четырех глав, заключения, списка литературы, включающего 88 наименований и приложений. Работа содержит 127 страниц, 50 рисунков и 17 таблиц. Приложения включают 6 страниц.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

1. Математические модели спекания шихты достаточно хорошо разработаны.

2. Процессы спекания и охлавдения не рассматриваются совместно.

3. В известной модели Бабушкина Н.М., Братчикова С.Г. и др. охлаждения агломерата не учитываются газодинамические процессы.

Показана необходимость создания системы прогнозирования процесса охлаждения агломерата на прямолинейном охладителе.

На основе проведенного анализа сформулированы задачи работы и определены параметры прогнозирования: температура материала и газа, скорость газа и давление. На рис.1 представлена структурная схема процесса агломерации железных руд, которая включает структурные схемы для системы спекания шихты (газодинамики слоя, процессов теплообмена, процессов массообмена при сушке, горения твердого топлива в слое, диссоциации известняка, формирование газовой смеси в слое) и системы охлаждения агломерата (газодинамики слоя и процессов теплообмена).

Сформулированы функциональные требования к математическому обеспечению системы прогнозирования охлаждения агломерата. На рис. 2 представлена функциональная схема системы прогнозирования процесса охлаждения агломерата на прямолинейном охладителе.

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

Схема расчетной области представлена на рис. 3 для нахождения скорости и давления в охлаждаемом слое.

Модель процессов газодинамики в слое представлена системой уравнений:

дР _

- уравнением движения газа в слое дисперсного материала: т~ - а(11'г )"г; (1)

- уравнением неразрывности газового потока: — (рг№г) = О,

ду

где а(0 = А"+ к,р .п'г, Кг вязкостный коэффициент газодинамического сопротивления, м"2; К2 - инерционный коэффициент газодинамического сопротивления, м"1; Р - давление, Па; у - координата высоты слоя, м; IV;.- скорость газа на свободное сечение, м3/(м2-с); ц - динамическая вязкость газа, кг/(м с); рг - плотность газа, кг/м3.

Начально« состояние системы

Парвшары входа

Рис. 1. Структурная схема процесса агломерации железных руд

Из (1) и (2) получено уравнение для определения давления в каждой точке по высоте слоя дисперсного материала:

4-^1 = 0, (3)

- начальное условие т = 0: Р = Рат.,,, где т - время, с;

Рис. 2. Функциональная схема системы прогнозирования процесса охлаждения агломерата на

прямолинейном охладителе

Рис. 3. Схема расчетной области для определения температуры агломерата (Тм=(2, у)), газа (Тг=(2,у)), давления (Р=(г,у)) и вертикальной скорости движения газа по толщине слоя агломерата (1Уг=(г,у)): 1 - газ, 2 - отходящий газ, Н-высота слоя, Ратм и Рк - атмосферное давление и давление в камере.

{Р \ ~ Н

о '

где Р„= Рк- ДРкр\ Рк - давление в камере, Па; ДР,р = /2- сопротивление

колосниковой решетки, Па, рЛ -плотность газа на выходе из слоя при нормальных условиях, кг/м3; 1Угк - скорость газа на выходе из слоя при нормальных условиях, м/с; - коэффициент сопротивления колосниковой решетки.

Уравнения (1)^(3) решаются на каждом временном шаге общего алгоритма решения систем уравнений для сечения г, которое перемещается со скоростью движения охлаждаемого слоя. При этом положение сечения г на прямолинейном охладителе определяется по формуле ъ = и 1, где и - линейная скорость движения охлаждаемого слоя, м/мин. После расчета давления в каждой точке слоя определяется скорость газа по формуле:

1Г,=_1_* (4)

2К2р, ^

Модель процессов конвективного теплообмена представлена следующими уравнениями:

- уравнение баланса тепла газов:

-К-Ргсг ^=«,(/,-О; (5)

ду

- уравнение баланса тепла агломерата:

- начальное условие: гм = ;

- граничное условие =(Д

где !г и /„ - температуры газа-теплоносителя и материала, °С; т - время, с; координата по высоте слоя, м; с,„ -теплоемкость материала, Дж/(К-кг); IVг -

скорость газа, м/с; р„ - кажущаяся плотность материала, кг/м ; рг - плотность кгЛ Вт

газа, кг/м3; т - порозность слоя; а г - объемный коэффициент теплоотда-

чи,

м ■ и

Объемный коэффициент теплоотдачи а у определяется по формуле Китаева Б.И.:

.„ 0.9 0.3

(7)

где (1 - диаметр куска агломерата, м; М' = 1,0 для идеальных условий, М' = 0,5 при наличии мелочи 20% (Л/'зависит в основном от равномерности распределения газов по поперечному сечению слоя). Алгоритмы решения системы уравнений (3) (7) разработаны с помощью метода конечных разностей.

Тестирование алгоритма расчета процессов конвективного теплообмена в неподвижном охлаждаемом слое проводилось путем сравнения результатов моделирования с точным решением задачи Шумана, полученным В.Н. Тимофеевым в интегральной форме. В результате тестирования установлено, что для того, чтобы относительная погрешность определения температуры не превышала 1 %, необходимо выбирать количество узлов в расчетной области не менее 30 для заданной высоты слоя 0,5м.

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

Расстояние от загрузки агломерата, м

Рис. 4. Зависимость температуры отходящего газа от расстояния загрузки агломерата: " - экспериментальные данные; — - статистическая обработка экспериментальных данных; ""- результаты моделирования.

1,6

0 1.4

1 1,2

§ 0,6 | 0,4 ° 0,2 0

0 10 20 30 40 50 60 70 80 90 Расстояние от загрузки агломерата, м

Рис. 5. Зависимость скорости газа от расстояния загрузки агломерата: - экспериментальные данные; _в— - результаты моделирования.

Для температуры отходящего газа погрешность составила 7,1%, а для скорости газа - 4,0%. Установлено, что разработанное математическое обеспечение отображает основные закономерности конвективного теплообмена и газодинамики. Из рис. 4, 5 следует, что результаты моделирования удовлетворительно согласуются с экспериментальными данными.

Разработанная модель составила основу математического и программного обеспечения.

В третьей главе на основе математической модели охлаждения агломерата выполнен анализ входных параметров для выявления доминирующих. С этой целью для условий прямолинейного охладителя ЧерМК ОАО «Северсталь» изучали влияние начальной температуры агломерата, скорости движения охлаждаемого слоя, высоты охлаждаемого слоя и перепада давления газа на изменение температуры отходящего газа на контрольной отметке 1, в которой установлена термопара для измерения температуры отходящего газа на промышленном охладителе. Эта отметка находится на 45 метре от загрузки агломерата на прямолинейный охладитель.

Предварительно определили диапазон изменения входных параметров. Для получения диапазона изменения начальной температуры агломерата применили математическую модель спекания шихты на агломерационной машине для условий ЧерМК ОАО «Северсталь». Для определения диапазона остальных параметров использовали данные, полученные на промышленном охладителе.

Результаты определения диапазона изменения входных параметров приведены в табл. 1. На рис. 6^9 представлены результаты изучения влияния входных параметров на температуру отходящего газа в указанной контрольной отметке 1.

Результаты анализа влияния входных параметров на температуру отходящего газа представлены в табл. 1.

рлИ

o J--------

600 650 700 750 800 850 900

Начальная температура агломерата, С

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

160 -Ч 140 -

ю

s 120 -S.ioo-

I 80 -0>

= 60 • н 4020 * 0 -

1 1,2 1,4 1,6 1,8 2 72 2,4 2,6 Скорость движения охлаждаемого слоя, м/мин

Рис. 7. Зависимость температуры отходящего газа на отметке 1

от скорости движения полотна: - - результаты моделирования; ♦ - экспериментальные данные.

200 180 160 Ч 140

I 120

I юо g 80 | 60 40 20 0

0,3 0,35 0,4 0,45 0,5 0,55 0,6 0,65 0,7 Высота слоя, м

i

------'

Рис. 8. Зависимость температуры отходящего газа на отметке 1 от высоты охлаждаемого слоя: - - температура газа; -"- температура материала.

4,5

Рис. 9. Зависимость температуры отходящего газа на отметке 1 от перепада давления: - - температура газа; температура материала.

Таблица 1

Входные параметры процесса охлаждения

№ Наименование параметра Ед.изм. Пределы изменения параметра (Д/7 ) Пределы изменения температуры отходящего газа (АТ) на от-1 о,-, метке 1, С Д Г ЛП

Нижний Верхний Нижний Верхний

1. Начальная температура агломерата °С 600 900 39 56 0,05

2. Скорость движения охлаждаемого слоя м/мин 1 2,6 14 150 85

3, Высота охлаждаемого слоя М 0,3 0,7 18 118 250

4. Перепад давления кПа 0,5 4,5 34 280 61,5

Анализ показал, что значительное влияние на процесс охлаждения оказывает высота охлаждаемого слоя (^ =280), которая зависит от способа загрузки

агломерата на прямолинейный охладитель. Остальные параметры в меньшей мере влияют на температуру отходящего газа: начальная температура агломерата - 0,05; перепад давления - 61,5 и скорость движения охлаждаемого слоя -85. Однако температура агломерата и высота слоя являются нерегулируемыми параметрами. Поэтому в системе прогнозирования в первую очередь будем

\ \

\ \

к

—^н

~~~ -—•■

0,5 1 1,5 2 2,5 3 3,5 4

Перепад давления, кПа

учитывать регулируемые параметры: скорость движения охлаждаемого слоя и перепад давления.

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

Скорость движения охлаждаемого слоя зависит от скорости движения паллет на агломерационной машине:

■9«,=^..., (8)

где к - эмпирический коэффициент, определяется опытным путем; За ч - скорость движения паллет на агломерационной машине, м/мин; Зоа - скорость движения охлаждаемого слоя на прямолинейном охладителе, м/мин.

В работе выведена формула для вычисления эмпирического коэффициента:

Игт^-0-/»). (9)

Н о'оРо

где р* - насыпная плотность готового агломерата на агломерационной машине, кг/мЗ; Ро - насыпная плотность агломерата на прямолинейном охладителе, кг/мЗ; нл - высота агломерируемого слоя, м; Н0 - высота охлаждаемого слоя, м; 1Л - ширина агломерируемого слоя, м; 10 - ширина охлаждаемого слоя,

м; & - доля горячего возврата агломерата с агломерационной машины, кг/с.

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

Шаг 1. Ввод технологических параметров по агломерационной машине.

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

Шаг 3. Расчет коэффициента к по формуле (9).

Шаг 4. Расчет прогнозируемой скорости движения охлаждаемого слоя по формуле (8).

Алгоритм прогнозирования давления газа в дутьевой камере реализован на основе математической модели процесса охлаждения агломерата путем прогнозирования перепада давления в контрольной отметке 1.На рис. 10 представлена типичная зависимость температуры отходящего газа от перепада давления газа.

Алгоритм прогнозирования давления газа в дутьевой камере представлен в виде нескольких этапов.

Этап 1. Построение функции / температуры отходящего газа от перепада давления газа_/(А/>)=Л

Зависимость температуры отходящего газа от перепада давления (рис. 9) аппроксимировали с помощью квадратного уравнения:

¡=/{АР) = аАР+ ЬД/>+ с.

Рис. 10. Типичная зависимость температуры отходящего газа от перепада давления газа

Для построения функцииЛ АР) найдены коэффициенты а, Ь и с путем вычисления температуры отходящего газа I с помощью математического обеспечения в трех точках (АРЫп, 1т(п), ( ДРтах, 1тах), (ДР.„, /„), где АР^- минимальное значение перепада давления газа, Па; АР шах - максимальное значение перепада давления газа, Па; АРМ- моделируемое значение перепада давления газа, Па; 1тт, ¡тах и (,, - температура отходящего газа минимальная, максимальная и моделируемая, соответственно,°С.

Решая систему уравнений (10), можно найти значения коэффициентов а, Ь и с в конкретном случае:

•1ав=аЬРав2+Ь&Ряв+с. (Ю)

/„ =аАР„2 +ЬАРМ+с

Этап 2. Расчет прогнозируемого перепада давления газа по заданной температуре отходящего газа 13 на отметке 1.

Для вычисления прогнозируемого перепада давления газа необходимо использовать квадратное уравнение (11):

13 = аАР* + ЬАР7 + с, (11)

где а, Ь и с - коэффициенты; ДР3 - прогнозируемый перепад давления газа по заданной температуре, Па; Г, - заданная температура отходящего газа на контрольной отметке 1, °С.

Решением уравнения (11) будет только одно положительное значение.

Этап 3. Расчет давления газа в дутьевой камере.

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

Блок схема алгоритма прогнозирования давления газа в дутьевой камере на

прямолинейном охладителе представлена на рис. 11.

-

Г

Расчет АРтп

Расчет

ДР™.

Вычислить а,Ь,с

Конец

Вычислить

Р

Растет АРМ

Вычислить

АР,

Растет давлений с помощью математического обеспечения

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

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

Рис. 11. Блок схема алгоритма прогнозирования давления газа в дутьевой камере на прямолинейном охладителе

В четвертой главе приведены результаты экспериментальных исследований метода и алгоритмов прогнозирования скорости движения охлаждаемого слоя и давления газа в дутьевой камере на ЧерМК ОАО «Северсталь».

Экспериментальные исследования выполнены для агломерационной машины АКМ 312 и прямолинейного охладителя ОП 315. На рис. 12 представлены экспериментальные данные по скорости движения паллет на агломерационной машине и скорости движения охлаждаемого слоя на прямолинейном охладителе.

На рис. 13 -И 4 представлены результаты прогнозирования коэффициента к (9) и скорости движения охлаждаемого слоя (8) во времени. Величина к в реальных условиях принимают постоянной в интервале 0,5 0,6, что способствует неравномерной высоте охлаждаемого слоя по длине прямолинейного охладителя и колебанием температуры отходящего газа и материала на контрольной отметке 1. Полученное изменение скорости движения охлаждаемого слоя

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

3,5 -

1,5

О 1

0,5 0

10

20

40

50

60

30

Время, с

Рис. 12. Изменение скорости движения во времени: — скорость движения паллет на агломерационной машине; скорость движения охлаждаемого слоя на прямолинейном охладителе.

0,66

0,64

^ 0,62

£ 0,6 х

I 0,58 вв- 0,56 «

* 0,54 0,52 0,5

щ и * ■ у 1 к

«.......... (\ш ■■

Р. м ! ш ' > • \ *

■ *

0 10 20 30 40 50 60

Время, м

Рис. 13. Изменением коэффициента к во времени. - - фактические данные, прогнозируемые.

I 2,1 у— X

: а

5 2,05 -Г

0 10 20 30 40 50 60

Время, м

Рис. 14. Изменение скорости движения охлаждаемого слоя: - - фактические данные, прогнозируемые.

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

Рис. 15. Изменение скорости движения охлаждаемого слоя

Время, ч

Рис. 16. Изменение температуры отходящего газа во времени: —' - фактическая температура отходящего газа, °С; - заданная температура отходящего газа, °С.

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

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

Время, ч

Рис. 17. Изменение давления газа в дутьевой камере во времени — - фактическое давление газа, Па; -»- - прогнозируемое давление газа, Па.

Алгоритм прогнозирования давления газа в дутьевой камере имеет погрешность 7,35%.

Ожидаемый годовой экономический эффект от внедрения системы прогнозирования, обеспечивающей устранения возгорания транспортерной ленты, составит 1 млн. 44 тыс. руб.

Заключение содержит перечень основных результатов диссертационной работы.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

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

2. Разработана математическая модель, которая прогнозирует в зависимости от начальной температуры агломерата, скорости движения и высоты охлаждаемого слоя и давления в дутьевой камере, температуру отходящего газа с точностью 7,1%, а скорость газа - 4,0%.

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

4. Разработаны алгоритмы прогнозирования скорости движения охлаждаемого слоя и давления газа в дутьевой камере прямолинейного охладителя, обеспечивающие такие значения технологических параметров, при котором температура отходящего газа поддерживается на заданном уровне. Погрешность прогнозирования давление газа в дутьевой камере - 7,35%.

Разработанная система прогнозирования прошла опытно-промышленную проверку. Результаты переданы в агломерационное производство ЧерМК ОАО «Северсталь». Ожидаемый годовой экономический эффект от внедрения системы прогнозирования составит 1 млн. 44 тыс. руб.

Список публикаций по теме диссертации по перечню ВАК:

1. З.К. Кабаков, А.А.Елисеев, Ю.В. Веселов. Математическая модель тепло- и массообменных процессов при агломерации шихты // Изв. вуз. Черная металлургия. 2008. № 1. С. 19 - 22.

2. З.К. Кабаков, A.A. Елисеев, Ю.В. Веселов. О косвенном критерии качества агломерата// Изв. вуз. Черная металлургия. 2008. № 5. С. 10-12.

в прочих изданиях:

3. З.К. Кабаков, A.A. Елисеев, Ю.В. Веселов. Математические модели спекания шихты на агломерационных машинах //Автоматизация и энергосбережение машиностроительного и металлургического производств, технология и надежность машин, приборов и оборудования. Материалы четвертой международной научно-технической конференции. I том. г.Вологда. ВГТУ. 2008. с. 178-181

4. З.К. Кабаков, A.A. Елисеев. Ю.В. Веселов. Математическая модель тепловых и газодинамических процессов при охлаждении агломерата// Современные наукоемкие технологии. 2009. № 8. С. 28 - 30.

5. З.К. Кабаков, A.A. Елисеев, Ю.В. Веселов. Исследование процесса охлаждения агломерата // Молодые исследователи - регионам. Всероссийская научная конференция, г. Вологда. ВГТУ.2009. С. 12 - 13.

Лицензия А № 165724 от 11 апреля 2006 г.

Подписано к печати 16.11.09г. Тир. 100. Усл. печ. л. 1. Формат 60x84. Зак. 656. ГОУ ВПО «Череповецкий государственный университет» 162600 г.Череповец, М.Горького, 14.

ВВЕДЕНИЕ.

1 Общая характеристика проблемы прогнозирования охлаждения агломерата.

1.1 Состояние проблемы прогнозирования охлаждения агломерата.

1.2 Характеристика охлаждения агломерата как объекта прогнозирования.

1.3 Функциональные требования к математическому обеспечению 34 системы прогнозирования охлаждения агломерата.

1.4 Выводы.

2 Математическое обеспечение метода прогнозирования охлаждения агломерата.

2.1 Математическое описание тепло-массообмена в агломерируемом слое.

2.1.1 Подмодель процессов газодинамики в слое.

2.1.2 Подмодель процессов конвективного теплообмена в неподвижном плотном слое.

2.1.3 Подсистема горения твердого топлива.

2.1.4 Подсистема диссоциации известняка.

2.1.5 Учет тепловых эффектов от фазовых превращений при плавлении и кристаллизации.

2.1.6 Подмодель массообмена при испарении и конденсации влаги.

2.1.7 Подмодель формирования газовой смеси.

2.2 Математическое описание тепловых и газодинамических процессов в охлаждаемом слое.

2.3 Разработка дискретной модели тепловых процессов в охлаждаемом слое.

2.4 Тестирование модели расчета процессов конвективного теплообмена в неподвижном плотном охлаждаемом слое.

2.5 Проверка адекватности модели.

2.6 Выводы.

3 Алгоритмы обработки информации и прогнозирования параметров охлаждения.

3.1 Формирование входных параметров процесса охлаждения.

3.1.1 Начальная температура агломерата.

3.1.2 Скорость движения полотна.

3.1.3 Высота охлаждаемого слоя.

3.1.4 Давление в дутьевой камере.

3.2 Разработка алгоритмов прогнозирования.

3.2.1 Алгоритм прогнозирования скорости движения охлаждаемого слоя на прямолинейном охладителе.

3.2.2 Алгоритм прогнозирования давления в дутьевой камере прямолинейного охладителя.

3.3 Особенности применения алгоритма прогнозирования давления в дутьевой камере в промышленных условиях.

3.4 Выводы.

4 Экспериментальные исследования метода и алгоритмов обработки информации в системе прогнозирования охлаждения.

4.1 Описание экспериментального оборудования.

4.2 Прогноз технологических параметров.

4.3 Исследование эффективности алгоритмов прогнозирования.

4.3.1 Алгоритм прогнозирования скорости движения охлаждаемого слоя на прямолинейном охладителе.

4.3.2 Алгоритм прогнозирования давления в дутьевой камере прямолинейного охладителя.

4.4 Перспективы применения системы прогнозирования в системе управления прямолинейным охладителем.

4.5 Выводы.

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

После окончания процесса спекания шихты среднемассовая температура агломерата обычно находится в пределах 70(К800°С. Высокая температура агломерата создает трудности при последующей его обработке на аглофабрике, транспортировке и загрузке в доменную печь. В связи с этим на современных аглофабриках агломерат охлаждают на охладителях различного типа, в том числе на прямолинейных, что облегчает условия работы обслуживающего персонала и оборудования. На практике известны случаи возгорания транспортерных лент при транспортировке агломерата в доменный цех.

Благодаря фундаментальным работам большого круга ученых, достигнуты значительные успехи в области совершенствования технологии производства агломерата. Наибольший вклад в развитие теории и технологии агломерационного процесса внесли работы ученых Бабушкина Н.М., Тимофеева В.Н., Коротича В.И., Пузанова В.П., Фролова Ю.А., МайзеляГ.М., Братчикова С.Г., Базилевича С.В., Вегмана Е.Ф., Белоцерковского Я.Л., Бабошина В.М., Шкляра Ф.Р., Сигова А. А., ШурхалаВ.А. и др. Однако, практически все известные методы ориентированы на прогнозирование процесса спекания шихты. Поэтому в настоящее время разработка метода, позволяющего прогнозировать процесс охлаждения агломерата на прямолинейном охладителе в системе управления агломерационным производством, представляется актуальной научно-технической задачей.

Объект исследования: процесс охлаждения агломерата на прямолинейном охладителе.

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

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

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

- разработка метода прогнозирования, основанного на математической модели процесса охлаждения агломерата на прямолинейном охладителе;

- разработка алгоритмов обработки информационных сигналов, позволяющих прогнозировать технологические параметры;

- проведение экспериментальных исследований работоспособности и эффективности предложенных метода и алгоритмов.

Методы исследований. Для решения поставленных задач в работе использовались теоретические основы газодинамики и теплотехники агломерационного процесса; методы статистической обработки информации; методы математического моделирования. Математическое моделирование производилось путем разработки численной модели для реализации в среде разработки Delphi 7. Достоверность полученных результатов подтверждена путем сравнения с экспериментальными данными.

Научная новизна и основные положения, выносимые на защиту:

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

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

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

Практическая ценность работы заключается в следующем:

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

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

Реализация результатов работы. Разработанные метод и алгоритмы обработки информационных сигналов, реализованные в системе прогнозирования процесса охлаждения агломерата, прошли опытно-промышленную эксплуатацию в агломерационном производстве на Череповецком металлургическом комбинате (ЧерМК) ОАО «Северсталь».

Теоретические результаты внедрены в учебный процесс на кафедре программного обеспечения ЭВМ ГОУ ВПО «Череповецкий государственный университет» по дисциплине «Основы теории управления».

Апробация работы. Основные положения диссертационной работы докладывались и получили положительную оценку на 4-ой международной научно-технической конференции «Автоматизация и энергосбережение машиностроительного и металлургического производств, технология и надежность машин, приборов и оборудования» (Вологда, 2008), Всероссийской научной конференции «Молодые исследователи - регионам»

Вологда, 2009) и на постоянно действующем научно-техническом семинаре кафедры программное обеспечение ЭВМ Череповецкого государственного университета.

Публикации. По материалам диссертации опубликовано 5 печатных работ, в том числе 2 статьи опубликованы в рецензируемых научных журналах, входящих в перечень ВАК Министерства образования и науки Российской Федерации.

Структура и объем работы. Диссертационная работа состоит из введения, четырех глав, заключения, списка литературы, включающего 88 наименований и приложений. Работа содержит 127 страниц, 50 рисунков и 17 таблиц. Приложения включают 6 страниц.

4.5 Выводы

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

2. Разработанные алгоритмы прогнозирования позволяют улучшить экономические показатели путем экономии затрат по замене транспортерных лент.

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

108

ЗАКЛЮЧЕНИЕ

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

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

2. Разработана математическая модель, которая прогнозирует в зависимости от начальной температуры агломерата, скорости движения и высоты охлаждаемого слоя и давления в дутьевой камере, температуру отходящего газа с точностью 7,1%, а скорость газа - 4,0%.

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

4. Разработаны алгоритмы прогнозирования скорости движения охлаждаемого слоя и давления газа в дутьевой камере прямолинейного охладителя, обеспечивающие такие значения технологических параметров, при котором температура отходящего газа поддерживается на заданном уровне. Погрешность прогнозирования давление газа в дутьевой камере - 7,35%.

1. Авдеенко, А.А. Методика расчета горения топлива и окислительно-восстановительных процессов при агломерации Текст./ Авдеенко А.А., Боковиков Б.А., Исаенко Г.Е., Клейн В.И., Ярошенко Ю.Г. // Сталь. 2002. № 4, с. 34 - 36.

2. Ашихмин, А.А. Закономерности выгорания углерода в псевдоожиженном слое применительно к котлам малой мощности Текст.: дисс. канд. Тех.наук:05:14.04, защищена 14.03.87/Ашихмин А.А. Свердловск: УПИ.-1987.-24с.

3. Бабушкин, Н.М. Охлаждение агломерата и окатышей Текст. / Бабушкин Н.М., Братчиков С.Г., Намятов Г.Н., Швыдкий B.C., Шкляр Ф.Р., Ярошенко Ю.Г. //М.Металлургия, 1975 г. 208с.

4. Бабушкин, Н.М. Теплообмен в слое агломерационной шихты (математическая формулировка задачи) Текст./ Бабушкин Н.М., Тимофеев В.Н.//Теплотехника доменного и агломерационного процессов: Сб. науч. тр. ВНИИМТ. -М.: Металлургия, 1966. № 14. С. 160 171.

5. Бабушкин, Н.М. Горение топлива в слое агломерационной шихты Текст./ Бабушкин Н.М., Тимофеев В.Н.// Теплотехника доменного и агломерационного процессов: Сб. науч. тр. ВНИИМТ. М.: Металлургия, 1966. № 14. С. 160-171.

6. Базилевич, С.В. Теплотехнические расчеты агрегатов для окускования железорудных материалов Текст./ Базилевич С.В., Бабошин В.М., Белоцерковский Я.Л. и др. М., «Металлургия», 1979. 208 с.

7. Базилевич, С.В. Агломерация Текст./ Базилевич С.В., Вегман Е.Ф. М.: Металлургия, 1967. 368с.

8. Бондаренко, В.Д. Пути снижения энергетических затрат при агломерации Текст./ Бондаренко В.Д., Мусиенко К.А., Руденко Н.Р.// Металлургическая и горнорудная промышленность, 2006 №4.

9. Братчиков, С.Г., Худорожков И.П. «Изв. вузов. Черная металлургия», 1970, №4, с.32-36 с ил.

10. Вегман, Е.Ф. Металлургия чугуна Текст./ Вегман Е.Ф., Жеребин Б.Н., Похвистнев А.Н., Юсфин Ю.С., Клемперт В.М.//М.: Металлургия, 1989. -512 с.

11. Викулов, Г.С. Производство агломерата заданного химического состава Текст./Вилков А.Е., Кабанов Ю.А., Татаркин Н.Л.// Металлург.-2002.-№6.-С.47-49.

12. Дмитриева, Е.Г. Совершенствование технологии производства агломерата на основе анализа закономерностей горения твердого топлива Текст.:диссертация. кандидата технических наук : 05.16.02 Екатеринбург, 2007 185 с. РГБ ОД, 61:07-5/2917

13. Елисеев, А.А. Исследование тепло-массообменных процессов при агломерации шихтыТекст.: дис. канд.техн.наук : 05.14.04 Череповец, 2006 196 с. РГБ ОД, 61:06-5/3822

14. Журавлева, А.Г. Инженерная модель теплообмена в слое агломерируемой шихты Текст./ Журавлева А.Г., Мысик А.Ф., Фролов Ю.А.//

15. Теплотехническое обеспечение основных технологических процессов черной металлургии: Сб. науч. тр. ВНИИМТ. М.: Металлургия, 1988. С. 5 - 9.

16. Кабаков, З.К. Математическая модель тепло-массообменных процессов при агломерации шихты Текст./ Кабаков З.К., Елисеев А.А., Веселов Ю.В. // Изв. вуз. Черная металлургия. 2008. № 1. С. 19 22.

17. Кабаков, З.К. Математическая модель тепловых и газодинамических процессов при охлаждении агломерата Текст./ Кабаков З.К., Елисеев А.А., Веселов Ю.В.// Современные наукоемкие технологии, «Академия естествознания». 2009. № 8. С. 27 30.

18. Кабаков, З.К. Исследование процесса охлаждения агломерата Текст./Кабаков З.К., Елисеев А.А., Веселов Ю.В.// Молодые исследователи -регионам. Всероссийская научная конференция студентов и аспирантов. 1618 апреля 2009г, г. Вологда, № 1. С. 12 13.

19. Кабаков, З.К. О косвенном критерии качества агломерата Текст./ Кабаков З.К., Елисеев А.А., Веселов Ю.В.// Изв. вуз. Черная металлургия. 2008. №5. С. 10-12.

20. Карабасов, Ю.С. Использование топлива в агломерации Текст./ Карабасов Ю.С., Валавин B.C. М.: Металлургия, 1976. 263 с.

21. Ксендзык, Г.В. Изменение гидравлического сопротивления агломерационной шихты по ходу процесса спекания Текст./ Ксендзык Г.В.// Известия вузов Чер. мет. — №7, 1958, стр. 3 — 16.

22. Киссин, Д.А. Механизм минералообразования офлюсованного агломерата Текст./ Киссин Д.А., Литвинова Г.И.// Сталь.-1960.-№5.-с.9-14.

23. Китаев, Б.И. Теплообмен в доменной печи Текст./ Китаев Б.И., Зобнин Б.Ф., Ратников В.Ф., Телегин А.С., Лисиенко В.Г., Братчиков С.Г., Казяев М.Д., Маркин В.П., Суханов Е.Л., Сучков В.Д. Изд-во «Металлургия», 1966г.

24. Коротич, В.И. Агломерация рудных материалов Текст./ Коротич В.И., Фролов Ю.А., Бездежский Г.Н.Научное издание. Екатиринбург: ГОУ ВПО "УГТУ УПИ", 2003,. 400 с.

25. Коротич, В.И. Влияние зоны переувлажнения на газопроницаемость слоя материалов в процессе агломерации Текст./ Коротич В.И., Пузанов В.П.// Известия вузов Чер. мет. 1965. № 4, с. 53 - 58.

26. Коротич, В.И. Газодинамика агломерационного процессаТекст./ Коротич В.И., Пузанов В.П. М., «Металлургия», 1969. 208 с. с ил.

27. Коротич, В. И. Горение топлива и окислительно-восстановительные процессы при агломерации железорудных материалов Текст./ Коротич В. И.: Текст лекций. УГТУ, 1996. 64 с.

28. Коротич, В.И. Основы теории и технологии подготовки сырья к доменной плавке Текст./ Коротич В.И. М., «Металлургия», 1978. 208 с. с ил.

29. Коротич, В.И. Образование зоны переувлажнения при агломерации методом просасывания Текст./ Коротич В.И., Пузанов В.П.// Известия вузов Чер. мет. 1964. № 10, с. 28-33.

30. Коротич, В.И. Поведение влаги и газодинамика слоя в начальный период агломерации Текст./ Коротич В.И., Пузанов В.П., Фролов Ю.А. Сообщение 1 // Известия вузов Чер. мет. 1968. № 10, с. 26 - 30.

31. Коротич, В.И. Поведение влаги и газодинамика слоя в начальный период агломерации Текст./ Коротич В.И., Пузанов В.П., Фролов Ю.А. Сообщение 2 // Известия вузов Чер. мет. 1968. № 12, с. 37 - 41.

32. Лавров, Н.В. Введение в теорию горения и газификации топлива Текст./Лавров Н.В., Шурыгин А.П.-Москва: Изд. АН ССР.-1968.-217с.

33. Ладыгичев, М.Г. Сырье для черной металлургии Текст./ Ладыгичев М.Г., Чижикова В.М., Лобанов В.И., Винтовкин А.А., Буткарев А.П.,

34. КокоринЛ.К., Жилкин В.П., Доронин Д.Н.: Справочное издание: В 2-х т. Т. 1. Сырьевая база и производство окускованного сырья. М. Машиностроение-1, 2001 .-896с.

35. Лисиенко, В.Г. Теплофизика металлургических процессов Текст./ Лисиенко В.Г., Лобанов В.И., Китаев Б.И.Учебник для вузов. М.: Металлургия, 1982. 240 с.

36. Лыков, А.В. Теория тепло- и массопереноса Текст./ Лыков А.В., Михайлов Ю.А. М. Л. Госэнергоиздат. 1963. 536.

37. Маковский, В.А. Цифровая динамическая модель агломерационного процесса Текст./ Маковский В.А., Власюк Ю.Н.// Известия вузов Чер. мет. -1976. №8, с. 136- 139.

38. Махорин, К.Е. Сжигание топлива в псевдоожиженом слоеТекст./ Махорин К.Е., Хинкс П.А./ АН УССР. Ин-т газа. Киев: Наукова думка. 1989.-200 с.

39. Миллер, В.Я. Теоретические основы агломерации железных руд Текст./ Миллер В.Я.// Труды НТО ЧМ.-Т.УП.-М.:Металлургиздат, 1956.-е. 152-177.

40. Миллер, В.Я. Требования к качеству агломерата и мероприятия по его улучшению Текст./ Миллер В.Я.// Труды научно-технического общества черной металлургии.-T.XVI.-M.: Металлургиздат, 1958.-С.285-302.

41. Миллер, В.Я. Исследование прочности агломерата Текст./ Миллер В.Я., Базилевич С.В., Худорожков И.П., Майзель Г.М.// Сталь,- 1961.- с.769-777.

42. Натансон, И.П. Краткий курс высшей математики, издание четвертое стереотипное Текст./ Натансон И.П. Издательство «Лань», Санкт-Петербург, 2001г., 727с.

43. Нурмаганбетов, Ж.О. Удельный расход воздуха на агломерацию Текст./ Нурмаганбетов Ж.О., Коротич В.И.// Известия вузов Чер. мет. 1992. № 4, с. 10, 11.

44. Нурмаганбетов, Ж.О. Удельный расход воздуха на агломерацию Текст./ Нурмаганбетов Ж.О., Коротич В.И.// Известия вузов Чер. мет. 1992. № 6, с. 1 -3.

45. Попов, Г.Н. Новое в исследовании механизма переувлажнения шихты в агломерируемом слое Текст./ Попов Г.Н. // Известия вузов Чер. мет. 1987. № 1, с. 23-26.

46. Основы практической теории горения. Под ред. В.В. Померанцева. Учебное пособие для студентов вузов. Л.: «Энергия», 1973. 264 с.

47. Радованович, М. Сжигание топлива в псевдоожиженом слое / Пер. с англ. -М.: Энероатомиздат, 1990. -248 с.

48. Раева, М.В. Модель тепло- и массообмена при сушке дисперсного слоя Текст./ Раева М.В., Шкляр Ф.Р., Фролов Ю.А.// Металлургическая теплотехника: Сб. научн. тр. ВНИИМТ. М., Металлургия, 1974. №2. С. 154-162.

49. Ростовцев, С.Г. Агломерация криворожских железных руд Текст./ Ростовцев С.Г., Мееров С.М.//Домез.-1934.-№11-12,-с. 7-39.

50. Ростовцев, С.Г. Физико-химические основы процесса агломерации криворожских руд Текст./ Ростовцев С.Г.//Теория и практика металлургии.-1938.-№6.-c.3-9.

51. Сигов, А.А. Диссоциация карбоната кальция при агломерации Текст./Сигов А.А.// Известия вузов Чер. мет. 1958. № 3, с. 3 - 12.

52. Сигов, А.А. Перераспределение влаги при агломерации железных руд. Текст./ Сигов А.А. // Известия вузов Чер. мет. 1958. № 8, с. 7 - 12.

53. Телегин, А.С. Тепло-массоперенос Текст./ Телегин А.С., Швыдкий B.C., Ярошенко Ю.Г.: Учебник для вузов. М., «Металлургия», 1995. 400с.

54. Телегин, А.С. Термодинамика и тепло-массоперенос Текст./ Телегин А.С., Швыдкий B.C., Ярошенко Ю.Г. М., «Металлургия», 1980. 264 с.

55. Фастовский, М.Х. Механическое и транспортное оборудование агломерационных фабрикТекст./ Фастовский М.Х., Дакалов Г.В., Носовский А.А. М.: Металлургия, 1983.-265с.

56. Фролов, Ю.А. Трехмерная математическая модель для исследований и управления процессом агломерации Текст./ Фролов Ю.А., Полоцкий Л.И., Кобелев В.А., Конопляник В.В.// ОАО «Черметинформация». Бюллетень «Черная металлургия» ноябрь 2005 г., С. 29-30.

57. Хзмалян, Д.М. Теория топочных прочесовТекст./ Хзмалян Д.М.: Учеб. пособие для вузов. М.: Энергоатомиздат, 1990. - 352 с.

58. Хохлов, Д.Г. Физико-химические особенности процесса получения офлюсованного агломерата и пути улучшения его качества Текст./Хохлов Д.Г. // Окускование железорудного сырья: Тем.сб. науч.тр.-Свердловск: Уралмеханобр. Вып. 12.-1965.-c. 112-129.

59. Хохлов, Д.Г. Исследования по формированию агломерата и разрушаемости его при восстановлении Текст./ // Хохлов Д.Г.,

60. Фофонова Е.Г., Климова Н.С. Окускование железных руд и концентратов: Тем.отрасл.сб.-Свердловск: Уралмеханобр.-1973.-№1 .-с.5-17.

61. Шапран, А.А. Модель переувлажнения и скорости спекания агломерируемого слоя Текст./ Шапран А.А.// Известия вузов Чер. мет. -1985. №4, с. 104- 108.

62. Шкляр, Ф.Р. Прогрев слоя газом с движущимся фронтом постоянной температуры Текст./ Шкляр Ф.Р., Бабушкин Н.М., Тимофеев В.Н. //Теплотехника доменного и агломерационного процессов: Сб. науч. тр. ВНИИМТ. М.: Металлургия, 1966. № 14. С. 172 - 189.

63. Шкляр, Ф.Р. Закономерности нагрева неподвижного слоя Текст./ Шкляр Ф.Р., Тимофеев В.Н., Раева М.В.// Нагрев и охлаждение стали. Теплотехника слоевых процессов: Сб. науч. тр. ВНИИМТ. М.: Металлургия, 1970. № 23. С. 180 - 194.

64. Шурхал, В.А. Расчет состава газовой фазы, образующейся при горении углерода в процессе агломерации и обжига Текст./ Шурхал В.А. // Известия вузов Чер. мет. 1969. № 8, с. 27 - 31.

65. Шурхал, В.А. Расчет расхода воздуха и количества отходящих газов при агломерации железорудной шихты // Известия вузов Чер. мет. 1969. № 10, с. 26-29.

66. Fukagama, Т. Автоматизированная система управления агломерационным процессом с искусственным интеллектом. Текст./ T.Fukagama // Кавасаки сэйтэцу тихо — Kawasaki steel gino. — 1991.-№3.-С.203-209

67. Hamada, Katsushige. Прогрессивная АСУ ТП агломерационным процессом Текст./ Hamada Katsushige, Murai Tatsunori, Jyoko Tadasugu, Nakamura Yuji, Morioki Keiji // Сумитомо киндзоку Sumitomo Metals. -1992. №1.- C.151-160

68. Hideyuki, Yamaoka Development of a 3-D Sinter Process Mathematical Simulation Model Текст./ Hideyuki Yamaoka and Takazo Kawaguchi // ISIJ International, Vol. 45 (2005), No. 4, pp. 522-531.

69. Iwada, Kunihiro. Применение экспертной системы на аглофабрике №3 в Тобата Текст./ Iwada Kunihiro // Дзайре то пуросэсу Curr.Adv.Mater. and Proc. - 1991.-№1.-С.123.

70. Yang, W. Modeling of Combustion and Heat Transfer in an Iron Ore Sintering Bed with Considerations of Multiple Solid Phases Текст./ Yang W., Ryu Ch., Choi S., Choi E., Lee D., Huh W.// ISIJ International, Vol. 44 (2004), No. 3, pp. 492 499.

71. Young, R.W. Dynamic mathematical model of sintering processTeKCT./ Young R.W./ Ironmaking and Steelmaking, 1977, № 6, p. 321 328.

72. Young, R.W. // Ironmaking and Steelmaking, 1977, № 6, P. 321 328.