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

Дипломная работа - Компьютеры, программирование

Другие дипломы по предмету Компьютеры, программирование

нив ступенчатое воздействие к вышеописанным каналам управления. Используем следующую программу (Листинг 2):

Исходные данные:

Канал управления: расход разбавителя - концентрация в аппарате:

 

v:=2.1

+0.1

 

теплоёмкость вещества, кдж/кгградК

теплоёмкость хладагента, кдж/кг градК

плотность вещества в аппарате и входных потоках, кг/л

плотность хладагента, кг/л

коэффициент теплопередачи, кдж / м мин градК

поверхность теплообмена, м

тепловой эффект реакции, кдж/моль

универсальная газовая постоянная, дж/моль градК

концентрация компонента А на входе, моль/л

расход на выходе из реактора V, л/мин

расход на входе в реактор V2, л/мин

расход на входе в реактор V1, л/мин

расход хладагента, л/мин

температура на входе в реактор t1, градС

температура на входе в реактор t2, градС

температура хладагента на входе, градС

 

Vp:=550 - объем реактора

Vv:=200 - объем рубашки

 

. Постоянные коэффициенты

 

 

. Константа скорости

 

 

. Вектор-функция правых частей дифференциальных уравнений модели

 

 

. Начальные условия и условия интегрирования

 

Ca

 

6. Решение дифференциальных уравнений модели

 

 

 

 

Рис 6. Реакция на изменение расхода разбавителя v2 (?v2=0.1 литр/мин).

а) изменение концентрации компонента А

б) изменение концентрации компонента В

в) изменение концентрации компонента С

г) изменение концентрации компонента D

д) изменение температуры в реакторе

е) изменение температуры хладагента

 

Рис. 7. Кривая разгона по каналу расход разбавителя - концентрация в аппарате

Канал управления: расход хладагента - температура в аппарате:

Исходные данные:

теплоёмкость вещества, кдж/кгградК

теплоёмкость хладагента, кдж/кг градК

плотность вещества в аппарате и входных потоках, кг/л

плотность хладагента, кг/л

коэффициент теплопередачи, кдж / м мин градК

поверхность теплообмена, м

тепловой эффект реакции, кдж/моль

универсальная газовая постоянная, дж/моль градК

концентрация компонента А на входе, моль/л

расход на выходе из реактора V, л/мин

расход на входе в реактор V2, л/мин

расход на входе в реактор V1, л/мин

расход хладагента, л/мин

температура на входе в реактор t1, градС

температура на входе в реактор t2, градС

температура хладагента на входе, градС

 

+0.1

 

. Постоянные коэффициенты

 

 

. Константа скорости

 

 

. Вектор-функция правых частей дифференциальных уравнений модели

 

 

. Начальные условия и условия интегрирования

 

Ca

 

6. Решение дифференциальных уравнений модели

 

 

 

 

Рис 8. Реакция на изменение расхода хладагента vхл (?vхл=0.1 литр/мин).

а) изменение концентрации компонента А

б) изменение концентрации компонента В

в) изменение концентрации компонента С

г) изменение концентрации компонента D

д) изменение температуры в реакторе

е) изменение температуры хладагента

 

Рис. 9. Кривая разгона по каналу расход хладагента - температура в аппарате

 

Канал управления: расход на выходе из реактора- уровень в реакторе:

Исходные данные:

теплоёмкость вещества, кдж/кгградК

теплоёмкость хладагента, кдж/кг градК

плотность вещества в аппарате и входных потоках, кг/л

плотность хладагента, кг/л

коэффициент теплопередачи, кдж / м мин градК

поверхность теплообмена, м

тепловой эффект реакции, кдж/моль

универсальная газовая постоянная, дж/моль градК

концентрация компонента А на входе, моль/л

расход на выходе из реактора V, л/мин

расход на входе в реактор V2, л/мин

расход на входе в реактор V1, л/мин

расход хладагента, л/мин

температура на входе в реактор t1, градС

температура на входе в реактор t2, градС

температура хладагента на входе, градС

 

-0.1

 

. Постоянные коэффициенты

 

 

. Константа скорости

 

 

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

 

 

. Начальные условия и условия интегрирования

 

 

6. Решение дифференциальных уравнений модели

 

 

 

Рис 10. Реакция на изменение расхода на выходе v (?v=-0.1 литр/мин).

а) изменение концентрации компонента А

б) изменение концентрации компонента В

в) изменение концентрации компонента С

г) изменение концентрации компонента D

д) изменение температуры в реакторе

е) изменение температуры хладагента

 

L0=0.888, L10000=1.254

Рис. 11. Кривая разгона по каналу расход на выходе из реактора - уровень в реакторе