Моделирование вспучивания керамзита во вращающейся печи как объекта управления Текст научной статьи по специальности «Строительство и архитектура»

Энергетика

УДК 62.524

МОДЕЛИРОВАНИЕ ВСПУЧИВАНИЯ КЕРАМЗИТА ВО ВРАЩАЮЩЕЙСЯ ПЕЧИ КАК ОБЪЕКТА УПРАВЛЕНИЯ

С.Я. Галицков1, А.И. Данилушкин2, А. С. Фадеев1

1 Самарский государственный архитектурно-строительный университет 443001, г. Самара, ул. Молодогвардейская, 194

2 Самарский государственный технический университет 443100, г. Самара, ул. Молодогвардейская, 244

На основании принятых допущений разработана расчетная схема вращающейся печи для производства керамзита в виде четыгрехслойного цилиндра как объекта управления с распределенными параметрами. При моделировании использованы дифференциальные уравнения теплопроводности в цилиндрических координатах с учетом граничных условий, определяемый особенностями вспучивания. На примере печи СМ875Б типоразмера 2,5*40 м синтезирована структура объекта для решения задачи автоматизации процесса вспучивания керамзита по отношению к управляющему воздействию и возмущению (изменение влажности сырца) в виде структуры с сосредоточенными параметрами.

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

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

Конструкция вращающейся печи. Односекционная печь представляет собой [1] наклоненный на 3...5° к линии горизонта стальной барабан 3 (рис. 1, а), внутренняя поверхность которого футерована огнеупорным кирпичом. На корпусе печи имеются два бандажных кольца 4 (ими печь опирается на роликовые опоры 5, смонтированные на тумбах 6) и зубчатый венец 11, в зацепление с которым входит приводная шестерня 10, которая через редуктор 9 соединяется с электродвигателем 7. «Холодный» конец печи входит в загрузочную головку 2, а «горячий» - в разгрузочную откатную головку 12, в которой установлена горелка 13. Сырец керамзита, двигаясь вдоль оси печи, последовательно проходит четыре технологические зоны [2] тепловой обработки - сушку, нагрев, вспучивание и охлаждение (рис. 1, б).

Анализ кривой обжига позволяет осуществить ее линейно-кусочную аппрокси-

Станислав Яковлевич Галицков - д.т.н., профессор. Александр Иванович Данилушкин - д.т.н., профессор. Александр Сергеевич Фадеев - аспирант.

мацию [2]. В середине зоны вспучивания выделим точку С (рис. 1, б), которая определяет температуру вспучивания Тс. Конец зоны нагрева характеризуется градиен-, где ТА и ТВ - температуры в точках А и В. Установлено [1], что один

том

dz

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

Определение объекта управления, основные возмущения. Под объектом управления будем понимать совокупность тепловых и физико-химических процессов во вращающейся печи при вспучивании керамзита. Состояние объекта определяется

значениями температуры керамзита в точках А, В, С и величиной на участке

АВ. Эти процессы можно регулировать [1, 2] изменениями режима работы горелки и величины скорости вращения печи, что определяет температурное поле в печи, протяженность технологических зон (сушка, нагрев, вспучивание и охлаждение). Считаем, что скорость вращения печи постоянна, тогда управляющим воздействием на объект является тепловая мощность Q горелки.

/ 7

3

4

11 4

12 13

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

z

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

При моделировании объекта управления введем ряд допущений.

1. Считаем, что в горелке используется газообразное топливо, состав и температура которого неизменны. Геометрические параметры факела (длина и форма) настраиваются перед началом работы и затем не меняются.

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

[1], эти факторы меняют величину скорости дымовых газов не более чем на 5.. .10%. Поэтому допускаем, что скорость движения дымовых газов в осевом направлении постоянна по всей длине печи. Поток дымовых газов принимаем ламинарным, пренебрегаем движением газов в направлении, перпендикулярном оси печи.

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

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

Расчетная схема. На основании принятых допущений разработана расчетная схема печи (рис. 2) в виде отрезка трубы длиной L с внутренним диаметром футеровки D. Стенка состоит из трех слоев: корпус 1, футеровка 2, материал 3. Внутри цилиндра движутся дымовые газы со скоростью ¥дг, материал перемещается навстречу им со скоростью Ум. Ось z цилиндрической системы координат совпадает с осью цилиндра, начало координат (точка О) расположено на холодном торце. Температуру печи в сечении z на расстоянии г от оси OZ в момент времени т обозначим Т г, т).

Загрузка

маглЕриала

11 і

/

5 6

5 Й

II

Рис. 2. Расчетная схема вращающейся печи

Математическая модель тепловых процессов. Модель объекта управления представим совокупностью уравнений горелки [3], температурного поля и физикохимических процессов, происходящих в материале при его обжиге. Для получения модели четырехслойного цилиндра печи используем уравнение теплопроводности полого цилиндра с одной стенкой в цилиндрической системе координат [4]:

ВТ

дт

(д2Т 1 дТ д2Т ^

- = а

дг2 г дг дz2

(1)

где Т - температура стенки; а - коэффициент температуропроводности, а=Х/(срр), здесь X - коэффициент теплопроводности, ср - удельная теплоемкость, р - плотность; ВТ/дт - производная, связанная с движущейся материей,

дт

дт

дг

дz

(2)

здесь ¥г, V - скорости движения дымовых газов и материала по координатам Ог и Oz соответственно; принимая во внимание, что в расчетной схеме ¥г = 0, имеем:

ВТ^,г,т) дТ(z,г,т) ^ ^ дТ(z, г,т) дт дт 2 дz

поэтому уравнение теплопроводности (1) принимает вид

дТ ^, г,т) дТ (z, г,т)

дт

дz

=а

^ д2Т (z, г,т) 1 дТ (z, г,т) д2Т (z, г,т) ^

дг

дг

дz

(3)

(4)

применяя его к расчетной схеме печи (рис. 2) и учитывая 4-е допущение, получим систему из четырех дифференциальных уравнений:

д^(^ г,т)

дт

дТт ( ^ т)

дт

дТ/ (z, г, т)

дт

дТ5?(z, т)

дт

д Тйг (z, г,т) „ дТс1е (z, г,т)

= ат (z)

д Тт (z, г,т)

дz 2

дz

дТт (z, г,т)

дz '

=а

^ д2Т/ ^, г,т) 1 дТ/ (z, г,т) д2Ту (z, г,т) ^

дг2

дг

дz 2

(5)

=а

^д2Т^(z,г,т) | 1 дТ^(z,г,т) | г,т)^

5?

дг2

дг

дz

где первое уравнение описывает динамику температуры дымовых газов, второе -материала (керамзита), третье - футеровки и четвертое - корпуса; здесь Т^, а^, Р^; Тт, ат, Рт - температура, коэффициент теплопередачи и скорость движения дымовых газов и материала соответственно; Т/, а/, Т^, а51 - температура и коэффициент теплопередачи футеровки корпуса и стальной стенки печи соответственно. Система (5) нелинейна, т. к. коэффициенты дымовых газов и материала изменяют свои значения по длине печи: а,1г^), ат^) [1, 2].

Уравнения системы (5) связаны между собой граничными условиями (ГУ), выбор которых произведен в соответствии с [4] следующим образом. Рассматриваемые в настоящей работе границы сред разделим (рис. 3) на три группы. К первой отнесем две границы, образуемые холодным (Г1) и горячим (Г2) торцами печи с внешней окружающей средой. Третья граница (Г3) образуется наружной цилиндрической поверхностью печи с окружающей средой. Четвертая граница Г4 - это граница дымовых газов с гранулами керамзита. Ко второй группе отнесем границы: Г5 - образуется между поверхностями материала и футеровки, Г6 - между поверхностью футеровки и внутренней стенкой корпуса. В третьей группе - одна граница Г7 между источником тепла и дымовыми газами. На границах первой группы происходит кон-

г

г

г

вективный теплообмен между поверхностями твердых тел, воздушной средой (Г1, Г2, Г3) и дымовыми газами (Г4). Считаем, что этот процесс происходит при постоянном потоке тепла и имеет квазистационарный режим. Поэтому принимаем, что на этих границах имеет место ГУ 3-го рода. Границы и Г2 имеют по два участка вдоль оси г - первый определяется отрезком ^2 - R3), второй - ^3 - R4). Границы второй группы соответствуют теплообмену соприкасающихся твердых тел при условии, что температура на границе этих тел одинакова. Поэтому здесь имеет место ГУ 4-го рода [4]. Будем считать, что объемная тепловая мощность 2(т) горелки передается через поверхность факела. Тогда для седьмой границы можно принять ГУ 2-го рода. Исходя из этого граничные условия Т(г, г, т), приведенные к традиционному виду, в частных производных можно записать в виде системы:

Хг дТ(^,г,т) = « (Тс (г) - Т(0,г,т)),г е [Я 2...Я 3],

^ = «1 (Тс(т) - Т(0,г,т)),г е [Я3...Я4],

дz

Л

бТ (X, г,т)

"б

= «2 (Т (т) - Т (X, г,т)), г є[ Я2...Я3],

Л аТ(1:Г,Т) — «2 (Тс(т) - Т(X,г,т)),г є [Я 3...Я 4],

бг

Лт бТ(г,1 ^1,т) = «3 (Тс(г,Я1 - 1,т) - Т(г,Я1,т)),г є [0...Х],

бг

Л бТ(гЯ2,т) Л бТ(гЯ2,т) _ є[0

Л/ : — Лт : , 2 є [0...Х],

бг

бг

Л бТ(г.Я3,т) — Л/бТ(Г.Я3,т),2 є [0..Х],

бг

/

бг

Л,, бТ(z,Я4,т) — «4 (Тс(2,Я 4 - 1,т) - Т(г,Я4,т)),г є [0...Х],

бг

2(г^,0,т) — 2(т),г — 0,г є [2...8т],

где Тс(т) - температура окружающей среды.

(6)

\ ГУ 3-го рода ГУ 4-го рода

ГУ 3-го рода /

і ТГ

Рис. 3. Схема граничных условий

На рис. 3: Rl, R2, R3, R4 - радиусы границ раздела сред; а1, а2, а3, а4 - коэффициенты теплоотдачи, Сдг, См, Сф, Сст - теплоемкости, Хгд, Хм, Хф, Хст - коэффициенты теплопроводности, рдг, рм, рф, рст - плотности дымовых газов, обжигаемого материала, футеровки и стального корпуса печи соответственно; Тдгн, Тдгк ,Тмн, Тмк - температуры дымовых газов и обжигаемого материала в начале и конце их движения в печи.

Синтез структуры объекта управления. Исследуем динамику объекта в зоне малых отклонений АТЛ, АТВ, АТС от некоторого установившегося технологического режима печи под действием изменения управляющего воздействия AQ и возмущения Aw. Синтезирована структура объекта управления (рис. 4), где операторы A1, A2 и A3 являются математическими моделями объекта по управлению, и на их выходе получают изменение температур материала в точках А, В и С - ЛТ1, ЛTB, ЛТС. Блоки A4, А5 и A - математические модели по отношению к возмущению, на их выходах получают изменение температур материала ЛГ/, ЛТ®, АТС в тех же точках. В блоке А7 вычисляем производную dT/dz. Для нахождения операторов воспользуемся программной средой SolidWorks, использующей метод конечных объемов [5].

Пример моделирования объекта управления. На примере вращающейся печи СМ875Б типоразмера 2,5*40 м, оборудованной газовой горелкой С 199-08-100 (с расходом газа до 800 нм3/час, тепловой мощностью до 9,3 МВт) и осуществляющей обжиг керамзита из глины Смышляевского месторождения (с влажностью 15...30%, насыпной плотностью после формовки 1100 кг/м3, теплоемкостью 0,92 кДж/(кг К), теплопроводностью 2,56 Вт/(м°С)), что соответствует режиму работы печи, рассмотренному в [2], найдем операторы структуры (рис. 4) в форме передаточных функций и определим их параметры. В программной среде SolidWorks создана трехмерная геометрическая модель вращающейся печи по рис. 2, затем эта модель была модифицирована под приложение Flow Simulation той же программной среды.

Были заданы: параметры газовой среды [3]; материалы (из библиотеки SolidWorks): корпус вращающейся печи - сталь, футеровка - кирпич огнеупорный. Обжигаемый материал создан в редакторе материалов с параметрами: плотность 1000 кг/м3, теплоемкость 880 Дж/(кг-К), теплопроводность 1.74 Вт/(м-К), температура расплавления 1226.85 °C. Была задана на основании экспериментальных данных температура на границе корпуса печи и окружающей среды 300 °C.

Р и с . 4. Структурная схема объекта управления

Задание начальных условий: дымовые газы - температура 20 °C, плотность 1 кг/м3, скорость 30 м/с; загружаемый материал - температура 30 °C, скорость движения 0,016 м/с; мощность объемного теплового источника 35202 Вт/м3. Граничные условия были выбраны программой автоматически, по известным критериям [4, 6], которые использовались при составлении системы (6).

Оценка адекватности построенной вычислительной модели. Было задано время расчета 86400 с (1 сут), что соответствует выводу печи на установившийся технологический режим. В результате получена расчетная кривая обжига керамзита. По результатам ее сравнения с известной экспериментальной кривой, приведенной в

[2], была вычислена ошибка, которая не превысила 7% в интересующих нас зонах (зонах подготовки и вспучивания), что свидетельствует об адекватности модели в установившемся режиме.

Переходные процессы в объекте управления

Переходные процессы в объекте управления по отношению к управляющему воздействию (рис. 5). Для исследования переходных процессов в объекте управления «в малом» по отношению к ступенчатому изменению управления был поставлен эксперимент по следующей методике. Наблюдения температуры материала проводились в точках zA = 35 м, zB = 32 м и z,o = 30 м. Выбор максимальной величины ступенчатого воздействия ±AQ=5 664 Вт/м3 от значения Q0=35 202 Вт/м3 произведен по условию максимальной температуры вспучивания Тс=1200 °C. Этот диапазон разбит на 6 величин: три положительных значения (AQ = 1 888; 3 776; 5 664) и три отрицательных (AQ = -1 888; -3 776; -5 664). При исследовании динамики начальные условия для расчета принимались равными конечным условиям вывода печи на установившийся режим, т. е. моменту т = 86400 с, для исключения погрешностей на начальном участке расчета, связанных с использованием численных методов в SolidWorks; скачкообразные воздействия AQ прикладывалось в момент времени т = 2100 с.

Рис. 5. Переходные процессы по отношению к управляющему воздействию

Переходные процессы в объекте управления по отношению к возмущающему воздействию (рис. 6). Для исследования динамики объекта управления по возмущению эксперимент проводился следующим образом. Была задана начальная влажность подаваемого на обжиг в печь материала w0=15%. Было выбрано 4 значения Aw = -5%; -15%; 5%; 15%. Постановка эксперимента проводилась аналогично тому,

как это выполнялось по отношению к управлению. Наблюдения температуры материала проводились в тех же сечениях печи.

Рис. 6. Переходные процессы по отношению к возмущающему воздействию

Передаточные функции динамических моделей объекта управления

Точка Передаточные функции

по управлению по возмущению

С АТ ( п) Кует’сР т (п) = с (р) = КсЄ ^( р) Туср +1 (_) _Атс(р) _ ^ с Аб( р) Т-р 4-1

В ( р) = АТв (р) = ^ в АQ( п) тур+1 ^ (р) = АТ- (р) _ ^ АQ( р) Т-р + 1

А щ (п)=АОТТ = 0 Аб( р) ( р) = АТл (р) _ А АQ( р) т-р +1

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

Выводы

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

граммной среде SolidWorks ее вычислительная модель. Показано, что модель адекватно описывает формирование установившегося температурного режима печи.

2. Применительно к решению задачи автоматизации процесса вспучивания керамзита синтезирована структура печи в виде многомерной модели с сосредоточенными параметрами, где в качестве выходных координат рассматривается температура печи в трех характерных точках А, В и С, что позволяет установить два основных параметра печи, определяющих плотность керамзита, температуру вспучивания Тс и градиент dT/dz в конце зоны нагрева.

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

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Перегудов В.В., Роговой М.И. Тепловые процессы и установки в технологии строительных изделий и деталей: учебник для вузов. - М.: Стройиздат, 1983. - 416 с., ил.

2. Онацкий С.П. Производство керамзита. - 3-е изд., перераб. и доп. - М.: Стройиздат, 1987. - 333 с.

3. Галицков С.Я., Фадеев А.С. Математическое описание сжигания газа во вращающейся печи для обжига керамзита // Актуальные проблемы в строительстве и архитектуре. Образование. Наука. Практика: материалы 64-й Всероссийской научно-технической конференции по итогам НИР университета за 2006 г. - Самара, 2007. - 564 с.

4. Лыков А.В. Теория теплопроводности. - М.: Высшая школа, 1966. - 560 с.

5. Алямовский А.А. SolidWorks 2007/2008. Компьютерное моделирование в инженерной практике. -СПб, 2008. - 1040 с.

6. Михеев М.А. Основы теплопередачи. Изд. 2-е, стереотипное. - М.: Энергия, 1977. - 344 с.

Статья поступила в редакцию 1 июня 2011 г.

UDC 62.524

MODELING OF HAYDITE SWELLING IN ROTATING FURNACE AS AN CONTROL OBJECT

S. Galickov1, А. Danilushkin2, А. Fadeev1

1 Samara State Architectural-Building University 194, Molodogvardeyskaya st., Samara, 443001

2 Samara State Technical University

244, Molodogvardeyskaya st., Samara, 443100

The paper refers to the design of haydite production rotating furnace in the form of layer cylinder as a control object with distributed parameters. Differential equations in cylindrical coordinates for heat conductivity were used. As an example the furnace SM875B of a standard size 2,5х40м was used for control object for solving the problem of automation of haydite swelling process relatiue to control input and interference in the form of lumped parameters structure/

Keywords: the rotating furnace, expanded clay, a control object, heat transfer, object with distributed parameters, motion equations, control object approximation, adequacy, mathematical modeling.

S. Ya. Galickov - Doctor of Technical Sciences, Professor. A.I. Danilushkin - Doctor of Technical Sciences, Professor. A.S. Fadeev - Postgraduate student.