Моделирование вспучивания керамзита во вращающейся печи как объекта управления




НазваниеМоделирование вспучивания керамзита во вращающейся печи как объекта управления
страница1/7
Дата19.04.2013
Размер0.54 Mb.
ТипАнализ
  1   2   3   4   5   6   7
ВЕСТН. САМАР. ГОС. ТЕХН. УН-ТА. СЕР. ТЕХНИЧЕСКИЕ НАУКИ. 2011. № 2 (30)


Энергетика


УДК 62.524


МОДЕЛИРОВАНИЕ ВСПУЧИВАНИЯ КЕРАМЗИТА

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

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

1 Самарский государственный архитектурно-строительный университет

443001, г. Самара, ул. Молодогвардейская, 194

2 Самарский государственный технический университет

443100, г. Самара, ул. Молодогвардейская, 244

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

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

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


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

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


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


Рис. 1. Конструкция вращающейся печи и кривая обжига керамзита


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

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

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

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

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

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


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


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


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

µ §, (1)

где T ЁC температура стенки; а ЁC коэффициент температуропроводности, a=л/(cpс), здесь л ЁC коэффициент теплопроводности, ср ЁC удельная теплоемкость, с ЁC плотность; DT/ЃЭф ЁC производная, связанная с движущейся материей,

µ §, (2)

здесь Vr, Vz ЁC скорости движения дымовых газов и материала по координатам Or и Oz соответственно; принимая во внимание, что в расчетной схеме Vr = 0, имеем:

µ §, (3)

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

µ §; (4)

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

µ § (5)

где первое уравнение описывает динамику температуры дымовых газов, второе ЁC материала (керамзита), третье ЁC футеровки и четвертое ЁC корпуса; здесь Tdg, аdg, Vdg; Tm, аm, Vm ЁC температура, коэффициент теплопередачи и скорость движения дымовых газов и материала соответственно; Tf, аf; Tst, аst ЁC температура и коэффициент теплопередачи футеровки корпуса и стальной стенки печи соответственно. Система (5) нелинейна, т. к. коэффициенты дымовых газов и материала изменяют свои значения по длине печи: аdg(z), аm(z) [1, 2].

Уравнения системы (5) связаны между собой граничными условиями (ГУ), выбор которых произведен в соответствии с [4] следующим образом. Рассматриваемые в настоящей работе границы сред разделим (рис. 3) на три группы. К первой отнесем две границы, образуемые холодным (Г1) и горячим (Г2) торцами печи с внешней окружающей средой. Третья граница (Г3) образуется наружной цилиндрической поверхностью печи с окружающей средой. Четвертая граница Г4 ЁC это граница дымовых газов с гранулами керамзита. Ко второй группе отнесем границы: Г5 ЁC образуется между поверхностями материала и футеровки, Г6 ЁC между поверхностью футеровки и внутренней стенкой корпуса. В третьей группе ЁC одна граница Г7 между источником тепла и дымовыми газами. На границах первой группы происходит конвективный теплообмен между поверхностями твердых тел, воздушной средой (Г1, Г2, Г3) и дымовыми газами (Г4). Считаем, что этот процесс происходит при постоянном потоке тепла и имеет квазистационарный режим. Поэтому принимаем, что на этих границах имеет место ГУ 3-го рода. Границы Г1 и Г2 имеют по два участка вдоль оси r ЁC первый определяется отрезком (R2 ЁC R3), второй ЁC (R3 ЁC R4). Границы второй группы соответствуют теплообмену соприкасающихся твердых тел при условии, что температура на границе этих тел одинакова. Поэтому здесь имеет место ГУ 4-го рода [4]. Будем считать, что объемная тепловая мощность Q(ф) горелки передается через поверхность факела. Тогда для седьмой границы можно принять ГУ 2-го рода. Исходя из этого граничные условия T(z, r, ф), приведенные к традиционному виду, в частных производных можно записать в виде системы:

µ § (6)

где Tс(ф) ЁC температура окружающей среды.


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

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


Синтез структуры объекта управления. Исследуем динамику объекта в зоне малых отклонений ДТА, ДТВ, ДТС от некоторого установившегося технологического режима печи под действием изменения управляющего воздействия ДQ и возмущения Дw. Синтезирована структура объекта управления (рис. 4), где операторы A1, A2 и A3 являются математическими моделями объекта по управлению, и на их выходе получают изменение температур материала в точках А, В и С ЁC µ §, µ §, µ §. Блоки A4, A5 и A6 ЁC математические модели по отношению к возмущению, на их выходах получают изменение температур материалаµ §, µ §, µ § в тех же точках. В блоке А7 вычисляем производную dT/dz. Для нахождения операторов воспользуемся программной средой SolidWorks, использующей метод конечных объемов [5].


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

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

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


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


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


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

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


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


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


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


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

ТочкаПередаточные функциипо управлениюпо возмущениюСµ §µ §Вµ §µ §Аµ §µ §

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


Выводы

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

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

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

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

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

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

Похожие:

Моделирование вспучивания керамзита во вращающейся печи как объекта управления iconАвтоматизация технологического процесса вспучивания керамзита во вращающейся печи
Работа выполнена на кафедре «Механизация, автоматизация и энергоснабжение строительства» фгбоувпо «Самарский государственный архитектурно-строительный...
Моделирование вспучивания керамзита во вращающейся печи как объекта управления iconОтчет о выполнении ниокр по теме: «Разработка комбинированных методов управления для повышения эффективности работы вращающейся печи»
Вращающаяся цементная печь, Статическая система управления, срез печи, узел, статистические данные
Моделирование вспучивания керамзита во вращающейся печи как объекта управления iconОтчет о выполнении ниокр по теме: «Разработка комбинированных методов управления для повышения эффективности работы вращающейся печи»
Этап № «Декомпозиция печи на узлы и построение нечетких диаграмм поведения этих узлов на базе экспериментальных, статистических данных...
Моделирование вспучивания керамзита во вращающейся печи как объекта управления iconБлок управления управления печи 4,5кВт
Блок управления трехзонной печи (далее «блок») предназначен для точного регулирования температуры в трех зонах печи общей мощностью...
Моделирование вспучивания керамзита во вращающейся печи как объекта управления iconОтчет о выполнении ниокр по теме: «Разработка комбинированных методов управления для повышения эффективности работы вращающейся печи»
Этап № «Синтез и анализ работы комбинированных систем управления печью в статике и динамке»
Моделирование вспучивания керамзита во вращающейся печи как объекта управления iconОтчет о выполнении ниокр по теме: «Разработка комбинированных методов управления для повышения эффективности работы вращающейся печи»
Этап № «Создание нечеткой модели управления печью на основе теории нечеткой логики с использованием экспериментальных и статистических...
Моделирование вспучивания керамзита во вращающейся печи как объекта управления iconОтчет о выполнении ниокр по теме: «Разработка комбинированных методов управления для повышения эффективности работы вращающейся печи»
Вращающаяся цементная печь, комбинированная система управления, статические и динамические режимы, переходные процессы, передаточная...
Моделирование вспучивания керамзита во вращающейся печи как объекта управления iconРгп няц рк
Идентификация химического реактора установки по переработке натрия как технологического объекта управления
Моделирование вспучивания керамзита во вращающейся печи как объекта управления iconИсследование деловой репутации банков 20
Особенности российского коммерческого банковского сектора как объекта управления деловой репутацией 12
Моделирование вспучивания керамзита во вращающейся печи как объекта управления iconJ имитационное моделирование управления рисками инвестиционных проектов ]■
Имитационное моделирование управления рисками инвестиционный проектов [Электронный ресурс]: Дис канд экон наук : 08. 00. 13. И.:...
Разместите кнопку на своём сайте:
Библиотека


База данных защищена авторским правом ©lib.znate.ru 2014
обратиться к администрации
Библиотека
Главная страница