Моделирование системы управления химическим реактором

Министерство образования Российской Федерации

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

Кафедра Автоматизации  технологических процессов и  производств

 

.
Пояснительная записка  к курсовой работе по дисциплине «Моделирование систем»
Тема: «Моделирование системы  управления химическим реактором»
Выполнил
студент гр.		АТП-09
	Вислоухова Е.А
(подпись)	(Ф.И.О.)
	Принял доц. каф. АТП:
	А.С. Островский
(подпись)	(Ф.И.О.)

 

 

 

 

Пермь 2012

Оглавление

 

 

 

 

1. Краткие теоретические  сведения

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

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

линейные, когда происходит последовательное превращение одного вещества в другое: А®В®С®D (например, реакции полимеризации);

обратимые, когда происходит одновременное  протекание как прямого, так и обратного превращения: A⇄B (большинство химических реакций, например, окисление сернистого ангидрида в производстве серной кислоты);

реакции синтеза, когда из двух или  более реагентов образуется один целевой продукт реакции: А®В¬С (например, получение нитроглицерина);

реакции разложения, когда  из одного вещества образуются два  или более продукта: А¬В®С (например, процесс кальцинации).

На практике в чистом виде такие реакции, как правило, не встречаются. В промышленных химических процессах обычно протекают реакции смешанных типов, например:

                      

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

, (1)

где   – номер реакции;

 – номер реагирующего компонента;

– константа скорости -й реакции;

– концентрация -го компонента.

В дифференциальной форме  уравнение имеет вид

,          , (2)

где  – стехиометрические коэффициенты;

 – количество реакций;

 – количество реагирующих  компонентов.

Например, для сложной  химической реакции типа

выражения для скоростей  элементарных реакций можно представить в следующем виде:

 (3)

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

. (4)

В матрице (4) 0 обозначает, что реагент не участвует в  реакции; –1 – реагент расходуется  в реакции; +1 – реагент накапливается  в реакции.

Система дифференциальных уравнений (2) с учётом матрицы стехиометрических  коэффициентов (4) примет вид:

 (5)

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

, (6)

где   – температура смеси;

 – универсальная газовая постоянная;

 – энергия активации;

 – предэкспоненциальный множитель.

 

Для нахождения и составляется система уравнений для двух известных температур , :

 (7)

Прологарифмировав левые  и правые части уравнений системы (7), получим

 (8)

Вычитая из первого уравнения  системы (8) второе, получим

. (9)

Из (9) находим 

. (10)

 выражается из первого  уравнения (7):

. (11)

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

 (12)

Отношение вычисляется по формуле (10), а – по формуле (11).

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

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

Схема типового проточного реактора с мешалкой и теплообменным  устройством показана на рис.1.

 

Рис. 1. Схема типового проточного реактора

с мешалкой и теплообменным устройством

 

Основными параметрами, характеризующими материальные и энергетические потоки реакционного процесса в проточном реакторе с мешалкой, являются:

– вектор концентраций реагентов во входном потоке;

– температура входного потока;

– объемный расход хладоагента;

–температура хладоагента на входе;

– объемный расход реакционной массы;

– концентрация реагентов в выходном потоке;

– температура реакционной массы;

– температура хладоагента на выходе.

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

режим идеального смешения реакционной массы;

режим идеального смешения хладоагента в рубашке;

постоянство объема реагирующей  смеси в реакторе;

постоянство расходов реагирующей  смеси и хладоагента.

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

, (13)

где   – объём реактора, м³;

 – интенсивность источников  вещества: .

Учитывая матрицу стехиометрических  коэффициентов (4), составляем систему  дифференциальных уравнений вида (13).

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

 (14)

Подставим в (14) выражения  для скоростей реакций (3):

 (15)

Коэффициенты  находятся по формуле (6).

Химические реакции  отличаются выделением или поглощением тепла, поэтому система уравнений материального баланса (15) дополняется уравнением теплового баланса:

, (16)

где    – теплоемкость смеси, ;

– площадь теплообмена, ;

 – коэффициент теплопередачи;

– интенсивность источника тепла  за счет реакции.

, (17)

где – тепловой эффект i-й реакции, .

Разделим уравнение (16) на .

Для рассматриваемой реакции тепловой баланс (16) с учетом (17) будет выглядеть  следующим образом:

 (18)

где   ;

 – условная теплоёмкость, ;

 – плотность смеси,  .

Итоговая система дифференциальных уравнений, описывающая химический реактор, примет вид:

 (19)

Температура хладоагента в рубашке T_X рассчитывается по формуле

, (20)

где – теплоемкость хладоагента, .

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

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

Рис. 2. Структурная схема САР

 

Для упрощения процесса моделирования  можно предположить, что расходная характеристика регулирующего органа (РО) линейная (см. рис.3).

 

Рис. 3. Расходная характеристика РО

 

Регулирующее воздействие по температуре  производится изменением расхода хладагента. Таким образом, объемный расход хладагента описывается уравнением

, (21)

где – регулирующее (управляющее) воздействие.

Уравнение регулятора записывается в  виде

. (22)

Продифференцировав (21), получим

. (23)

Предполагаем, что  изменяется ступенчато, тогда

 (24)

и

. (25)

Включив в систему уравнений (19) выражения (21) и (25), получаем математическое описание САР химическим реактором.

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

2. Исходные данные

№	Сi0, кмоль/м3					Т1, К	Т2, К	КТ1i, 1/с (при Т1)					КТ2i, 1/с (при Т2)
4	1,5	0	0	0		360	370	0,7	0,5	0,7	1,5	–	KТ2i = 2KТ1i

r_х=1000 кг/м³; С_рх^(у)=1400 Дж/(кг К); v=0,5 м³/мин

Ср(у), Дж/(кг К)	rсм , кг/м3	КТ,  Вт/(м2К)	Тепловой эффект, –DH, Дж/моль					F, м2	V, м3	T0, К	Тх0, К
			–DH1	–DH2	–DH3	–DH4	–DH5
2000	1300	200	12 000	10 000	10 000	20 000		7	1.1	300	275

Схемы превращений 1®3, 2®1, 3®2, 3®4

 

Перевод исходных данных в одну систему  измерения:

v=0,5 м³/мин = 8,333∙10^-3 м³/с

С₁₀=1500 моль/м³

 

Схема превращений имеет следующий  вид

С₁           С₂           С₃           С₄

 

3. Расчеты коэффициентов

 

По формуле (10) и (11) определяем:

 

 

 

 

 

 

4. Система дифференциальных уравнений  химической кинетики

 

Выражения для скоростей элементарных реакций можно представить в следующем виде:

Система дифференциальных уравнений  примет вид:

Скорость химической реакции зависит  от температуры, что отражается на выражении  константы скорости реакции:

С учетом выражения для константы  скорости и исходных данных система  дифференциальных уравнений примет вид:

 

5. Моделирование химической кинетики: графики зависимостей при постоянной и переменной температурах

_{Рис. 4 График зависимости  Сj=f(t)  при постоянной температуре.}

Температура изменяется по закону:

_{Рис. 5 График зависимости Сj=f(t)  при переменной температуре.}

6. Система дифференциальных уравнений,  описывающих реактор

 

 

7. Моделирование химического реактора: графики зависимостей

и

 

                                                                                                                                         

_{Рис. 6 График зависимостей Сj=f(t) и T=f(t).}

 

Установившиеся значения:

8. Переходная характеристика по  каналу температура в реакторе  – расход хладагента

Подадим возмущение по расходу хладагента порядка 5%, т.е.

_{Рис. 7 Кривая разгона.}

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9. Результаты обработки переходной  характеристики

Передаточная функция объекта:

Настройки регулятора, рассчитанные методом расширенных частотных  характеристик: