|
|
|
|
|
|
|
|
|
11 |
|
|
|
|
|
|
м2; |
– |
коэффициент теплоотдачи |
от |
полимеризатора к |
стенке, |
||||||||||
ккал/(час м2 0С); |
– температура стенки полимеризатора, 0С. |
|
|||||||||||||
|
Температуру стенки полимеризатора можно получить, составив |
||||||||||||||
тепловой баланс для металлических стенок, омываемых полимеризатом: |
|
||||||||||||||
|
|
|
dTст |
F |
1 (TP Tст ) |
|
КЭКВ |
(Tст |
TЭ ) |
, |
(2) |
||||
|
|
|
dt |
|
|
|
|
|
J |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
где |
– |
коэффициент |
теплопередачи; |
– |
температура этилена, |
0С; – |
|||||||||
теплоемкость металла полимеризатора, ккал/ 0С. |
|
|
|
||||||||||||
|
Коэффициент теплопередачи от стенки к этилену определяется |
||||||||||||||
выражением: |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
КЭКВ |
|
|
|
1 |
|
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
1 |
|
|
DL |
|
|
(3) |
|||
|
|
|
|
|
|
( |
|
|
|
|
) |
|
|
|
|
|
|
|
|
|
|
2 |
|
LBD |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где
– коэффициент теплоотдачи от стенки к этилену, ккал/(час
м2
0С);
– толщина стенки, м;
– теплопроводность стали, ккал/(час
м
0С).
Коэффициент теплоотдачи от стенки к этилену определяется по формуле:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(4) |
где |
– давление этилена, |
атм; |
|
– |
тепловая |
нагрузка на поверхность |
|||||||||||
теплопередачи, ккал/(час м2), определяемая по формуле: |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(5) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Скорость образования полимера определяется следующим выражением: |
||||||||||||||||
|
|
|
|
|
|
|
|
GP |
V ( A1 |
MY1 |
A2 MY2 ), |
(6) |
|||||
где |
– молекулярная масса изобутилена, кг/моль; |
– молекулярная масса |
|||||||||||||||
изопрена, кг/моль. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
Материальный баланс реактора можно описать следующей системой |
||||||||||||||||
дифференциальных уравнений: |
|
|
|
|
|
|
|
|
|
||||||||
|
|
dM1 |
|
|
|
M10 |
M1 |
M1 M1 |
C k2 , |
|
|||||||
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
dM 2 |
|
|
|
|
M 20 |
M 2 |
M1 |
M 2 |
C k3 , |
(7) |
|||||
|
|
dt |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
dC |
|
|
|
C0 C |
M1 |
C k1, |
|
|
|||||||
|
|
dt |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
где
– концентрация изобутилена на входе, моль/л;
– концентрация изопрена на входе, моль/л;
– концентрация катализатора на входе, моль/л;
– время пребывания в реакторе, час; 

– константы скорости реакции.
Константы скорости реакции определяются из уравнения Аррениуса:
|
|
|
Ei |
|
(8) |
k |
k |
e R T , |
|||
i |
i0 |
|
|
|
|
12
где
– энергия активации i-ой константы, Дж/моль; 
– предэкспонента i-ой константы;
– температура в полимеризаторе, К; R – универсальная газовая постоянная, R=8,28 Дж/моль.
Зависимость давления этилена от температуры кипения описывается следующим уравнением с четырьмя постоянными:
lg P |
1243, 766 |
11, 213927 lgT 0, 01102331T 30, 470741. |
(9) |
|
|
|
|||
|
T |
|||
|
|
|
|
|
Обратная зависимость, полученная в результате аппроксимации, имеет следующий вид:
|
|
|
(10) |
T 84,9 5 1,15 P 190,7. |
|||
Помимо описанных выше зависимостей, в модели объекта учтены начальные концентрации шихты и катализатора в соответствующих растворах.
В результате проведенного над моделью эксперимента были получены статические характеристики, наглядно иллюстрирующие зависимость температуры полимеризации от входных параметров (расхода катализатора, расхода шихты, давления этилена) (рис. 3) и позволяющие формализовать метод управления исследуемой системой.
Рис. 3. Статические характеристики объекта управления
Следующим этапом построения СППР процесса пуска является разработка метода вывода объекта в рабочую точку (т. А рис. 3) пространства параметров. При этом данная задача разбита на два этапа: выход на рабочую температуру по регламентной траектории (что обеспечит безопасность производства) и выход в рабочую точку (что обеспечит заданное качество готовой продукции в связи с верным поддержанием молекулярно-массового распределения).
Первая подзадача решена с помощью применения системы нечеткого вывода, разработка которой включает в себя ряд этапов:
13
1)формирование базы правил систем нечеткого вывода;
2)фаззификация входных переменных;
3)агрегирование подусловий в нечетких правилах продукций;
4)активизация или композиция подзаключений в нечетких правилах продукций;
5)аккумулирование заключений нечетких правил продукций;
6)дефаззификация выходных переменных.
При формировании базы правил системы нечеткого вывода были определены множества входных и выходных лингвистических переменных, для каждой из которых были сформулированы соответствующие терм-множества (рис. 4), а также было составлено множество правил нечетких продукций.
Рис. 4. Терм-множества входных (а, в) и выходных (б, г) лингвистических переменных СППР этапа пуска
Для агрегирования подусловий был использован наиболее распространенный метод min-конъюнкции:
|
|
А А |
(x) min |
А (xi ), |
А |
(xm ) , |
(11) |
|
|
ij im |
ij |
im |
|
|
|
где А (xi ), |
А |
(xm ) – |
численные |
значения |
функций |
принадлежности по |
|
ij |
im |
|
|
|
|
|
|
каждой из предпосылок.
В качестве операции для активизации подзаключений использовалась min-активизация:
( y) min |
i |
, |
( y) , |
|
(12) |
Bi |
|
Bi |
|
|
|
где i – степень истинности предпосылки i-го правила; |
В |
( y) – функция |
|||
|
|
|
|
i |
|
принадлежности его заключения.
Для аккумулирования заключений использовался метод maxдизъюнкции, а для проведения процедуры дефаззификации – метод центра тяжести:
14
Ymax
y В ( y)dy
Ymin
y , (13)
Ymax
В ( y)dy
Ymin
где Ymin , Ymax – границы интервала носителя нечеткого множества выходной
переменной у.
В результате моделирования была синтезирована система нечеткого вывода с применением выбранных ранее методов для реализации различных этапов ее построения. Таким образом были получены поверхности вывода для каждой выходной лингвистической переменной (рисунок 5), отражающие поведение системы при тех или иных входных управляющих воздействиях.
Рис. 5. Поверхности вывода выходных лингвистических переменных СППР этапа пуска (скорость изменения расхода катализатора – а, скорость изменения давления этилена - б)
Из статических характеристик объекта управления видно, что, если при выходе на рабочую температуру реакции использовался только канал расхода катализатора, то вывод в рабочую точку происходит без дополнительных воздействий. Если же в совокупности с ним использовался канал давления этилена, то для выхода в рабочую точку необходимо перейти на поверхность, соответствующую номинальному давлению. Для этого в первую очередь необходимо определить момент завершения I-го этапа, что осуществляется следующим образом: при достижении -60 0С включается таймер; если в течение определенного интервала времени температура в полимеризаторе остается в заданном диапазоне, что соответствует установившемуся состоянию, оператору подается сигнал о возможности подключения контура вывода давления этилена на номинальное значение. После выхода объекта в рабочую точку пространства параметров системы процесс пуска считается завершенным, оператор передает управление локальной системе регулирования.
В общем виде метод вывода объекта управления в рабочую точку пространства параметров иллюстрируется рисунком 6.
15
Рис. 6. Метод вывода объекта управления в рабочую точку пространства параметров