.

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

Язык: русский
Формат: реферат
Тип документа: Word Doc
142 3424
Скачать документ

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

Санкт-Петербургский государственный технологический институт

(технический университет)

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

Г.В. Иванова

«Автоматизация технологических процессов

основных химических производств»

Методические материалы по курсу лекций

(в двух частях)

Часть 2.

2003г.

УДК 66-52:66(075)

Иванова Г.В. Автоматизация технологических процессов основных химических
производств: Методическое пособие. Часть 2 / СПбГТИ(ТУ).-СПб., 2003.-
70с.

Методическое пособие предназначено для курса лекций по учебной
дисциплине «Автоматизация технологических процессов основных химических
производств», являющейся дисциплиной специализации 210201 –
«Автоматизация технологических процессов химической промышленности»
учебного плана по специальности 210200.

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

Часть 2 методического пособия включает в себя физико-химические основы
технологических процессов, технологические схемы рассматриваемых
объектов управления, математические описания объектов управления,
постановку задачи автоматизации, типовые схемы автоматизации, типовые
решения автоматизации для массообменных и реакторных процессов.

Утверждено на заседании методической комиссии факультета Информатики и
управления 23 июня 2003г., протокол № 6.

Материалы к лекции №9

Общая характеристика массообменных процессов

Вид диаграммы равновесия

для системы с 3 степенями свободы

c1* =f(c2) при ?=const и P=const;

Рис.2.

На диаграмме обозначено:

c1 – концентрация компонента в газовой фазе, cy.

c2 – концентрация компонента в жидкой фазе, cх.

при cх =c2, равновесное значение концентрации в газовой фазе будет cy* =
c1*.

Вид диаграммы равновесия

для системы с 2 степенями свободы

c1* =f(c2) при Р=const.

Рис.3.

Каждая точка кривой рис.3 соответствует равновесному состоянию при
различных температурах.

Отношение концентраций фаз при равновесии называют коэффициентом
распределения:

(3).

,

т. е.как тангенс угла наклона касательной к линии равновесия, если она
нелинейна,

или как тангенс угла наклона самой линии равновесия, если она линейна.

Уравнение рабочей линии процесс массопередачи

при противотоке распределяющих веществ.

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

противотоком:

G1 c1

с1н с1к

с2к с2н

G2 c2

Рис.1.

Материальный баланс процесса по целевому компоненту можно записать в
виде:

(1),

или

(2).

Для произвольного сечения аппарата с концентрациями с1 и с2 можно
записать:

(3а).

или

(3б).

Выразим из (3б) зависимость с1 = f(с2):

(4).

Выражение (4) – уравнение рабочей линии (рабочих концентраций)
массопередачи.

.

Рабочая линия всего аппарата ограничена точками с координатами с1н , с1к
, с2н , с2н .

Вид рабочей линии

при противотоке распределяющих веществ.

Рис.2.

Уравнение рабочей линии процесса при прямотоке

распределяющих веществ.

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

прямотоком:

G1 c1

с1н с1к

с2н с2к

G2 c2

Рис.3.

Материальный баланс процесса по целевому компоненту можно записать в
виде:

(5),

или

(6).

Для произвольного сечения аппарата с концентрациями с1 и с2 можно
записать:

(7а)

или

(7б).

Выразим из (7б) зависимость с1 = f(с2):

(8).

Выражение (8) – уравнение рабочей линии (рабочих концентраций
массопередачи) при прямотоке распределяющих веществ.

.

Рабочая линия всего аппарата ограничена точками с координатами с1н , с1к
, с2н , с2н .

Вид рабочей линии

при прямотоке распределяющих веществ.

Рис.4.

Диаграмма с1-с2 при расположении рабочей линии

выше равновесной.

Рис.1.

Из диаграммы рис.1 следует:

Это означает, что целевой компонент будет переходить из фазы G1 в фазу
G2.

Диаграмма с1-с2 при расположении рабочей линии

ниже равновесной.

Рис.2

Из диаграммы рис.2. следует:

Это означает, что целевой компонент будет переходить из фазы G2 в фазу
G1.

Средняя движущая сила процессов массопередачи.

:

(1).

:

(2).

Основное уравнение массопередачи.

(1),

где М – масса переносимого вещества, кг/с;

?ср – средняя движущая сила процесса массопередачи.

Определение движущей силы по каждой из фаз:

(2),

(3),

где

?срy и ?срx – средние движущие силы, в фазах y и x;

Ky и Kx – коэффициенты массопередачи для фаз y и x.

Материалы к лекции №10

Автоматизация процесса кристаллизации

Система с резким возрастанием растворимости.

Рис.1.

– c*=f(() – кривая растворимости, характеризующая равновесие
концентрированного раствора при изменении (.

– сп= f(() – линия условной границы метастабильной области.

А – неустойчивая, лабильная область массового образования центров
кристаллизации;

Б – относительно устойчивая метастабильная область образования и
растворения кристаллов;

В – область ненасыщенных растворов.

Система с плавным изменением растворимости.

Рис.2.

Переход в область пересыщенных растворов происходит только при
значительном охлаждении.

При этом выделяется незначительное количество твердой фазы.

Рекомендуемый способ кристаллизации – получение сп удалением части
растворителя из раствора.

Рекомендуемый метод кристаллизации – вакуум-кристаллизация.

Система с незначительным изменением

растворимости.

Рис.3.

Рекомендуемый способ кристаллизации – получение сп путем выпаривания
растворителя из раствора.

Рекомендуемый метод кристаллизации – изотермическая кристаллизация.

Изотермическая кристаллизация – это кристаллизация с удалением части
растворителя испарением или вымораживанием.

Система с незначительным изменением

растворимости.

Рис.3.

Рекомендуемый способ кристаллизации – получение сп путем выпаривания
растворителя из раствора.

Рекомендуемый метод кристаллизации – изотермическая кристаллизация.

Изотермическая кристаллизация – это кристаллизация с удалением части
растворителя испарением или вымораживанием.

Кинетика процесса кристаллизации.

Скорость образования центров кристаллизации:

,

где

– число частиц, образующихся в единице объема в единицу времени;

,

EN – энергия активации зародышеобразования, (кдж/кг);

Сп и С* – концентрации пересыщенного и насыщенного растворов, (кг/м3);

m=2 – 4 – кинетический коэффициент, зависящий от типа
кристаллизующегося вещества.

Качественная характеристика скорости роста кристаллов.

Зависимости скорости кристаллизации от времени.

Рис.4.

;

.

t0 – t1 – индукционный период, т.е. период подвижного равновесия
зародышей с раствором.

Кривая 1 – при большой степени пересыщения имеет резкий максимум
скорости процесса в момент tmax.

Кривая 2 – при малой степени пересыщения имеет пологий максимум в
течение времени t2 – t3.

Количественные оценки скорости роста кристаллов на основании
диффузионной теории.

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

(1),

где

– коэффициент массоотдачи, кг/м2*с;

Сп-С*=(, (кг/кг);

F – поверхность кристалла, (м2).

Для аппаратов с мешалками коэффициент массоотдачи ( зависит от следующих
параметров:

,

где

а – характеристический размер кристалла;

n – число оборотов мешалки, об/мин;

dм – диаметр мешалки, м.

Процесс подвода вещества к поверхности кристалла идет с большой
скоростью. Кинетика процесса определяется скоростью встраивания молекул
в кристалл:

(2),

где

КB – константа скорости встраивания молекул в кристалл;

n – эмпирическая постоянная.

Оба процесса протекают с соизмеримыми скоростями:

(3),

где

К – общий коэффициент скорости процесса, определяемый из соотношения:

.

Учитывая, что К=f((, KB), а (=f(n), в целом можно считать:

.

Таким образом, скорость роста кристаллов определяется поверхностью
кристалла, движущей силой процесса и скоростью мешалки.

Объект управления

Изогидрический кристаллизатор непрерывного действия с мешалкой.

Рис.5.

В схеме принято:

Gс=Gмр+Gкр; (мр = (кр = (с =(;

Скр = 1, т.е. кристаллы чистые.

Работа объекта

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

В результате пересыщения раствора и при интенсивном перемешивании
происходит кристаллизация целевого компонента из раствора с образованием
кристаллов (Мкр(Gкр).

При этом концентрация раствора понижается и оставшаяся жидкая фаза Gмр в
смеси с Gкр в виде потока суспензии Gc выводится из процесса.

Показатель эффективности процесса – диаметр кристаллов, dкр.

Цель управления процессом – обеспечение dкр = dкрзд.

Материальный баланс по всему веществу

Уравнение динамики:

(1).

:

(2).

Материальный баланс по кристаллизуемому веществу.

Уравнение динамики:

(3).

:

(4).

из (2) и полагаем Скр=1:

(5).

Из (5) выразим Gкр в явном виде:

(6).

Выражение (6) представляет Gкр на основе материального баланса процесса
кристаллизации.

, которое определяется на основе физики процесса массопередачи:

(7),

где

– число кристаллов, которое образуется

за время (пр в объеме Vс;

– изменение массы одного кристалла в ед. времени, кг/с.

, то на основе физики массопередачи можно считать:

.

В целом, на основании (6) и (7) можно записать:

.

Математическое описание для размера частиц.

На основании диффузионной теории и правила Мак-Бена скорость роста
кристаллов можно представить через радиус частиц:

(9),

где

(10);

r – радиус кристалла, м; t – время, с;

;

Сп, С* – концентрации пересыщенного и насыщенного растворов, кг/м3;

Еа – энергия активации, кдж/кг;

– температура, К;

R – универсальная газовая постоянная, кдж/кг*К.

На основании (9) и (10) – диаметр кристалла можно представить:

.

Если процесс кристаллизации вести при ( = const = (зд и обеспечить Сн=
const, то Сп и С* будут предопределены, т.к. система имеет 2 степени
свободы (s=2).

Таким образом, dкр=dкрзд можно обеспечить стабилизацией ( при условии
Сн= const.

Тепловой баланс процесса кристаллизации.

Уравнение динамики:

(11).

Можно принять ( = (мр = (кр = (с .

:

(12).

На основании (11) и (12) можно считать:

Предпочтительное управляющее воздействие Gхл .

Информационная схема кристаллизатора

Рис.6.

;

Возможные неконтролируемые возмущения:

.

В целом, кристаллизатор является сложным многосвязным объектом.

Типовая схема автоматизации процесса кристаллизации

Рис.7.

Регулирование.

Регулирование ( в аппарате по подаче хладоагента Gхл – обеспечивает
косвенное регулирование показателя эффективности процесса: (= f (dкр).

Регулирование h по отбору маточного раствора Gмр – для обеспечения
материального баланса по жидкой фазе.

Стабилизация расхода исходного раствора Gр – для обеспечения заданной
производительности установки.

Контроль.

.

.

Уровень: h.

Сигнализация.

Значительные отклонения температуры ( от задания.

3. Материалы к лекции №11

Автоматизация процесса абсорбции

Равновесие в процессе абсорбции.

Число степеней свободы для системы бинарный газ+жидкость:

S = k – f + 2=3-2+2=3.

Переменные для данной системы: температура (, давление Р; концентрации
С.

Равновесие такой системы при постоянных ( и Р описывается законом Генри:

(1),

где m – коэффициент распределения:

(2),

где Е – константа Генри:

(3),

где q – дифференциальная теплота растворения; R – универсальная
газовая постоянная; С – константа.

На основании (2) и (3) коэффициент распределения m зависит от P и (
следующим образом: при Р(, m( ; при ((, Е (( . m(.

, увеличивается с увеличением давления Р( и уменьшением температуры ((.

Влияние Р и ( на среднюю движущую силу процесса абсорбции.

(фазовые диаграммы при противотоке распределяющих веществ)

Рис.1.

При Р1 и (1 , ?ср1; 2 – При Р2 > Р1 , ?ср2 ; 3 – При (3 > (1 , ?ср3

Результаты анализа диаграмм:

?ср =f ((, Р, сгн , сгк , сан , сак );

?ср2 > ?ср1 ; при Р(? ?ср (;

?ср3 Влияние направления движения потоков на средние движущие силы процесса абсорбции. Рис.2а. - рабочая линия процесса абсорбции при противотоке распределяющих веществ; - рабочая линия процесса абсорбции при прямотоке распределяющих веществ; - равновесная линия процесса абсорбции. Движение распределяющих веществ противотоком. Рис.2б. ). . Движение распределяющих веществ прямотоком. Рис.2в. ). . Выводы по характеристикам схем противотока и прямотока: . Кинетика процесса абсорбции. Уравнения массопередачи в процессе абсорбции: (4а), или (4б), где Мга - масса распределяемого компонента, переходящая из газа в абсорбент в единицу времени, кг/ч; F - поверхность массопередачи м2; ; ; ; , где (г - коэффициент массоотдачи от потока газа к поверхности контакта фаз, кг/(м2 *ч); (а - коэффициент массоотдачи от поверхности контакта фаз к потоку абсорбента, кг/(м2 *ч). Схема насадочного абсорбера. Рис.6. Объект управления Схема абсорбционной установки. Рис.5. Работа схемы. Исходная газовая смесь Gг и абсорбент Gа в холодильниках 1 и 2 охлаждаются до заданных температур (г0 и (а0 и противотоком подаются в колонну 3. В колонне 3 происходит извлечение целевого (распределяемого) компонента из исходной газовой смеси с помощью жидкого абсорбента. В результате массообменного процесса между газовой и жидкой фазами получают: в низу колонны - насыщенный абсорбент Gна с концентрацией целевого (распределяемого) компонента сна; в верху колонны - обедненную газовую смесь Gог с концентрацией целевого (распределяемого) компонента сог . Показатель эффективности процесса - концентрация распределяемого компонента в обедненной газовой смеси сог. Цель управления - обеспечение сог = согзд на минимально возможном для данной установки значении. Материальный баланс по целевому компоненту. Материальный баланс по целевому компоненту в газовой фазе. Уравнение динамики: (1), где Мгна - масса целевого компонента, переходящая из газовой фазы в жидкую в единицу времени, кг/ч. : (2). (3), где Мгна - определяется уравнением массопередачи. Материальный баланс по целевому компоненту в насыщенном абсорбенте. Уравнение динамики: (4). : (5). (6), где Mгна - определяется уравнением массопередачи. Материальный баланс по общему количеству целевого компонента в процессе абсорбции. Уравнение динамики: (7), : (8). (9). (10). Материальный баланс по жидкой фазе. Уравнение динамики: , (11), Уравнение статики: (12) . (13). Материальный баланс по газовой фазе. Уравнение динамики: (14), где Мог - мольная масса обедненной газовой смеси, кг/моль; Рог - давление в колонне, Па; (ог - температура в колонне (по газовой фазе), К, Vог - объем газовой фазы в колонне, м3 . Уравнение статики: (15). На основании (14) и (15) можно считать: (16), Предпочтительное управляющее воздействие Gог. Тепловой баланс в абсорбере. Уравнение динамики для холодильника1: (17). : (18). На основании (17) и (18) можно считать: (19). Предпочтительное управляющее воздействие Gхл1. Уравнение динамики для холодильника 2. (20). На основании (20) можно считать: (21). Предпочтительное управляющее воздействие Gхл2. Информационная схема для установки с показателем эффективности сог . Рис.7 . . . . Схема абсорбционной колонны как многосвязного объекта при показателе эффективности сог. Рис.8 Информационная схема для установки с показателем эффективности сна . Рис.9 Возможные управляющие воздействия, контролируемые и неконтролируемые возмущения те же, что и в системе с показателем эффективности сог. . Схема абсорбционной колонны как многосвязного объекта при показателе эффективности сна. Рис.10. Типовая схема автоматизации процесса абсорбции. Рис. 11. Регулирование. Регулирование сог по подаче абсорбента Gа - как показателя эффективности процесса абсорбции. Регулирование давления верха колонны Рв = Рог по отбору обедненной газовой смеси Gог - для обеспечения материального баланса по газовой фазе. Регулирование уровня hна по отбору насыщенного абсорбента Gна - для обеспечения материального баланса по жидкой фазе. Регулирование температуры исходных материальных потоков газа (г0 и абсорбента (а0 по подаче хладоагентов Gхл1 и Gхл2 соответственно - для обеспечения теплового баланса установки. Стабилизация расхода исходной газовой смеси Gг - для обеспечения заданной производительности установки. Контроль. расходы - Gг, Gа, Gог, Gна, Gхл1, Gхл2; ; давление - Рв, Рн, (Р; уровень насыщенного абсорбента - hна; концентрация - сог . Сигнализация. существенные отклонения сог от согзд ; значительное повышение Рв( > Рпред , при этом формируется сигнал «В
схему защиты».

Система защиты.

По сигналу «В схему защиты» – открывается магистраль Gог, закрываются
все остальные магистрали.

4. Материалы к лекции №12

Автоматизация процесса сушки

Основные параметры сушильного агента и материала, как влагоносителей.

Относительная влажность сушильного агента ( :

(1).

На основании уравнения Менделеева-Клапейрона можно получить :

(3).

Относительная влажность материала ( – это

отношение массы влаги Мвл к общей массе влажного материала М=Мсм + Мвл ,
или к массе абсолютно сухого материала Мсм:

Влажность, отнесенная ко всему веществу:

, где М=var

Влажность, отнесенная к массе абсолютно сухого материала:

, где Мсм= const .

Диаграммы равновесия при сушке.

при ( = (*.

На диаграмме принято (1* >(2*>(3*.

Рис.1.

Из диаграммы следует:

( ( (*( – эффективность сушки повышается;

При ((((*( – эффективность сушки повышается.

Диаграмма (*= f(() при (=const.

Рис.2.

На рис.2. приведены: 1 – диаграмма при (1; 2- диаграмма при (2;

-диаграмма при (3.

При условии, что (1 (1.

Из диаграммы рис.4 следует:

Кинетика процесса сушки.

Уравнение массопередачи 1-го периода сушки с постоянной скоростью (Т1 ):

(1),

где

К=(са – коэффициент массоотдачи от поверхности контакта в газовую
фазу;

.

.

наступает 2-ой период сушки.

Уравнение массопередачи 2-го периода сушки с уменьшающейся скоростью
(Т2 ):

(2),

где

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

.

.

Кривая скорости сушки.

Рис.5.

Возможны следующие соотношения периодов сушки:

Тсушки = Т1

Тсушки = Т2

Тсушки = Т1+ Т2.

При анализе равновесия мы рассмотрели процесс сушки, характеризующийся
первым периодом и уравнением массопередачи (1).

Объект управления.

Схема барабанной сушилки прямоточного действия

– топка; 2 – смесительная камера; 3 – сушильный барабан;

– циклон; 6 – вентилятор

Рис.6.

Газообразное топливо Gт подается с первичным воздухом Gпв через горелки
в топку 1, где сжигается для получения сушильного агента.

Формирование сушильного агента осуществляется в смесительной камере 2,
куда подается вторичный воздух Gвв.

Влажный материал подается с помощью автоматического дозатора 7 в
сушильный барабан 3. Барабан наклонно расположен и вращается со
скоростью 4-5 об/мин, так что материал перемещается вдоль барабана и
высушивается к моменту попадания в бункер 4 до определенной влажности
(см.

Сухой материал Gсм отгружается из бункера 4 автоматическим дозатором 7.

Отработанный сушильный агент Gса в циклоне 5 очищается от пыли и
вентилятором 6 выводится из процесса.

Показатель эффективности процесса – влажность сухого материала (см .

Цель управления процессом – поддержание (см = (смзд.

Структурная схема топки и смесительной камеры.

Рис.7.

Материальный баланс по газовой фазе.

Уравнение динамики:

(1),

Уравнение статики:

.
(2)

Тепловой баланс.

Уравнение динамики:

(3).

Уравнение статики:

(4)

Информационная схема топки

и смесительной камеры.

Рис.8.

– возможные неконтролируемые возмущения;

– возможные контролируемые возмущения;

– возможные управляющие воздействия.

Выводы по анализу топки и смесительной камеры как объекта управления:

;

.

Структурная схема сушильного барабана.

Рис.9.

Материальный баланс по твердой фазе.

Уравнение динамики:

,
(5),

где hм – высота слоя материала в барабане;

Sапп – сечение аппарата;

(м – плотность материала;

Wмса – количество влаги, удаляемой в процессе сушки из материала в
единицу времени, кг/с.

Уравнение статики:

(6).

На основании (5) и (6) можно принять:

.

В типовой схеме автоматизации стабилизируют Gвм и Gсм, что обеспечивает
также стабилизацию hм

Материальный баланс по газовой фазе.

Уравнение динамики:

(7),

где

Мсаб – мольная масса сушильного агента в барабане,

кг/моль;

Рсаб – разрежение в барабане, Па;

(саб – температура в барабане (по газовой фазе), К,

Vсаб – объем газовой фазы в барабане, м3.

Уравнение статики:

(8).

На основании (7) и (8) можно считать:

.

Предпочтительным управляющим воздействием является Gсавых.

Тепловой баланс в сушильном барабане.

Уравнение динамики:

(9).

:

(10).

В выражениях (9) и (10) принимаем:

;

r – удельная теплота испарения влаги , дж/кг.

На основании (9) и (10) можно принять:

(11),

где

(12).

Материальный баланс по количеству влаги в материале.

Уравнение динамики:

(13),

:

(14).

Из выражений (13) и (14) следует, что:

(15),

где Wмса – определяется уравнением массопередачи.

Материальный баланс по количеству влаги в сушильном агенте.

Уравнение динамики:

(16).

:

(17).

Из выражений (16) и (17) следует, что:

(18).

Материальный баланс по общему количеству влаги в процессе сушки.

Уравнение динамики:

(19),

.

:

(20).

На основании выражений (19) и (20) можно считать:

(21).

При анализе топки мы получили:

(22).

Из всех возможных управляющих воздействий, перечисленных в выражениях
(21) и (22), для регулирования конечной влажности материала наиболее
информативным является расход топлива Gт.

Информационная схема сушильного барабана.

Рис.10.

Возможные управляющие воздействия:

.

Возможные контролируемые возмущения:

.

Возможные неконтролируемые возмущения:

.

Возможные управляемые переменные:

.

Сушильный барабан является сложным многосвязным объектом.

Типовая схема автоматизации процесса сушки.

Рис.11.

Регулирование.

Регулирование (см по подаче топлива Gт – как показателя эффективности
процесса сушки.

– для обеспечения эффективности сгорания топлива.

Регулирование температуры сушильного агента на входе в барабан (савх по
подаче вторичного воздуха Gвв.

Регулирование разрежения в барабане Рсаб по отбору отработанного
сушильного агента Gсавых – для для обеспечения материального баланса
по газовой фазе.

Стабилизация расходов влажного и сухого материала Gвм и Gсм
автоматическими дозаторами – для обеспечения материального баланса по
твердой фазе.

Контроль.

расходы – Gт, Gпв, Gвв, Gвм, Gсм;

;

разрежение – Рб;

влажность – (см ((к).

Сигнализация.

существенные отклонения (см от (зд;

значительное повышение (савх;

Незапланированное отключение привода, при этом формируется сигнал «В
схему защиты».

Система защиты.

По сигналу «В схему защиты» – прекратить подачу материала и топлива в
сушильный агрегат.

5. Материалы к лекции №13

Автоматизация процесса ректификации (ч.1)

Диаграмма равновесия «? – x – y » при Р= const.

Рис.1.

Смесь двухкомпонентная: компонент А – низкокипящий; компонент В –
высококипящий.

Диаграмма позволяет определять:

?x и ?y – для любого состояния смеси при заданном Р;

x и y – для любой температуры смеси.

По данным диаграммы «? – x – y » может быть построена диаграмма
равновесия « x – y » при Р= const.

Диаграмма равновесия « x – y » при Р= const.

Рис.2.

На осях y и x откладывают значения yА и xА соответствующие 100%
содержанию компонента А в паровой и жидкой фазах соответственно;

Диагональ ОА соответствует составам смеси yА = xА.

Расположение равновесной линии выше диагонали означает, что пары
обогащены низкокипящим компонентом.

Чем ближе линия равновесия к диагонали, тем меньше разница составов
жидкой и паровой фаз и тем труднее разделяется смесь при ректификации.

Схема колонны ректификации для математического

описания рабочих линий процесса.

Рис.3.

Обозначения параметров процесса:

Gf, Gфл, Gд, G к – расходы питания, флегмы, дистиллята, куба;

xf, xд, xк – молярные составы соответствующих потоков;

R=Gфл/Gд – флегмовое число.

Расчет рабочего флегмового числа.

,

где xf – заданный состав жидкой фазы потока питания;

yf – соответствующий ему равновесный состав паровой фазы потока
питания.

Рабочее флегмовое число находят из соотношений:

R=1,3*Rmin+0,3

или

R=?* Rmin,

где ?=1,2-2,5 – коэффициент избытка флегмы.

Уравнение рабочей линии верхней (укрепляющей) части колонны.

(1),

где y, x – составы паровой и жидкой фаз верха колонны, молярные доли;
xд – состав дистиллята.

.

Уравнение рабочей линии нижней (исчерпывающей) части колонны.

(2),

.

Построение рабочих линий процесса ректификации.

Рис.4.

1 – Равновесная линия процесса ректификации;

2 – Линия равных концентраций паровой и жидкой фаз компонента А в
колонне.

Составы паровой и жидкой фаз в верхней и нижней частях колонны равны,
т.е соответствуют линии 2.

По xд находят точку С.

По xк находят точку А.

По xf проводят вертикаль.

из соотношения (1) и получают точку D.

Получают точку М на пересечении СD с вертикалью от xf.

Получают СМ – рабочую линию верха колонны.

Проводят МА – рабочую линию низа колонны.

Диаграмма «Рабочая линия-линия равновесия»

процесса ректификации.

Рис.7.

1 – равновесная линия процесса бинарной ректификации;

2 – Линия равных концентраций в паровой и жидкой фазах

Линия МА – рабочая линия низа колонны с параметрами:

.

Анализ диаграммы «Рабочая линия – линия равновесия».

По диаграмме определяют следующие параметры:

, где ? – к.п.д. тарелки.

Движущую силу процесса в любых точках колонны:

.

зависят:

от режимных параметров процесса, определяющих положение равновесной
кривой (Р и ?);

от xf , xд, xк – определяющих положение рабочей линии процесса;

, т.е. от Gфл, Gд, Gf;

от тепловых параметров Gf.

Требуемую движущую силу процесса ректификации можно обеспечить:

стабилизацией режимных параметров Р или ?;

стабилизацией параметров потока питания Gf и ? f;

стабилизацией или изменением флегмового числа R=G фл /Gд .

Объект управления

Схема ректификационной установки.

– ректификационная колонна; 2 – подогреватель потока питания;

– кипятильник; 4 – конденсатор (дефлегматор); 5 – флегмовая
емкость.

Рис. 8.

Описание установки.

Объект управления – ректификационная установка для выделения из
исходной жидкой смеси целевого компонента в составе дистиллята.

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

соответственно.

Работа установки.

Исходная смесь Gп (Gxf) нагревается в подогревателе потока питания 2 до
температуры кипения ?п0 и подается в колонну 1 на тарелку питания (i=f).

Исходная смесь стекает по тарелкам нижней части колонны в виде
жидкостного потока Gx в куб колонны, участвуя в массообменном процессе с
паровым потоком Gy.

Из куба колонны выводится кубовый продукт Gкуб. Часть кубового продукта
подается в кипятильник 3, где испаряется с образованием парового потока
Gy0 , который подается в низ колонны.

Паровой поток поднимается вверх колонны, контактируя с жидким потоком
и обогащаясь целевым компонентом.

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

Конденсат собирается во флегмовой емкости 5. Из сборника флегмы
отбирается два потока:

поток дистиллята Gд – целевой продукт;

поток флегмы Gфл – жидкая фаза, используемая для орошения верха
колонны.

Показатель эффективности процесса сд – концентрация дистиллята.

Цель управления процессом – обеспечение сд.=сдзд.

6. Материалы к лекции №14

Автоматизация процесса ректификации (ч.2).

Структурная схема ректификационной установки.

Рис.1.

Математическое описание низа колонны

Структурная схема куба и кипятильника.

Рис.2.

Тепловой баланс низа колонны (? н = ? 0 ).

Уравнение динамики:

(5).

Уравнение статики:

Gгр rгр + Gх1 Cрх1?х1 = Gy0 rk + GkCpk??н
(6).

На основании (5) и (6) можно считать:

?н = f (Gгр, Gк ).

Предпочтительное управляющее воздействие Gгр .

Материальный баланс по всему веществу.

Уравнение динамики :

(1).

Уравнение статики :

G x1 = G k + G y0
(2),

где ? k – плотность кубовой жидкости , кг/м3;

S k – сечение куба колонны, м2 ;

h k – уровень кубовой жидкости, м;

G x1 , G k , G y0 – массовые расходы потоков в кубе колонны.

На основании (1) и (2) можно считать:

h k = f(G k ,G y0 ).

Предпочтительное управляющее воздействие Gk .

Материальный баланс по легколетучему компоненту.

Уравнение динамики:

(3).

Уравнение статики:

G x1 C x1 = G k C к + G y0 C y0
(4)

Основные допущения:

Кипятильник с полным испарением, т.е. C y0=C x0;

Тепловой баланс кипятильника:

Обозначения:

М0 – масса жидкости в нижней части колонны, кг;

r гр – удельная теплота конденсации пара, дж/кг;

r k – удельная теплота испарения кубовой жидкости, дж/кг.

На основании (3) и (4) можно записать:

.

Предпочтительное управляющее воздействие Gгр .

Информационная схема низа колонны.

Рис.3.

Информационная схема низа колонны

или hk, Ck .

Рис.4

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

Структурная схема дефлегматора с флегмовой емкостью.

Рис.5.

Материальный баланс по всему веществу.

Уравнение динамики:

(7)

где ?фл – плотность флегмы , кг/м3 ;

Sфл – сечение флегмовой емкости , м2 ;

hфл – уровень флегмы , м;

Gyn, Gфл, Gдист – массовые расходы, кг/с.

Уравнение статики:

Gyn = Gфл + Gдист
(8).

На основании (7) и (8) можно считать:

hфл=f(Gyn, Gфл, Gд )

Предпочтительное управляющее воздействие Gдист .

Материальный баланс по целевому компоненту.

Уравнение динамики:

(9).

Уравнение статики:

Gyn Cyn = Gдист Cx n+1 + Gфл Cx n+1
(10).

На основании (9) и (10) можно считать:

Cдист=f(Gyn, Gфл, Gд )

Предпочтительное управляющее воздействие Gфл .

Тепловой баланс верха колонны ( ?в = ?н ).

Структурная схема n–ой тарелки

Рис.6.

Уравнение динамики:

(11).

Уравнение статики:

Gyn-1*Cpyn-1 *?yn-1 + Gфл*Cрфл *?фл =

Gyn *Cpyn *?в + Gxn *Cpxn *?в
(12).

Обозначения:

Мxn – масса паровой фазы наверху колонны;

Cpyn, Cpy,n-1, Cрфл, Cpxn – удельные теплоемкости паровой и жидкой фазы
на n-ой тарелке;

Gyn-1, Gyn, Gxn – расходы паровой и жидкой фазы на n-ой тарелке.

Предпочтительное управляющее воздействие Gфл .

Баланс по паровой фазе.

Структурная схема конденсатора без флегмовой емкости.

Рис.7

Уравнение динамики:

(13).

Уравнение статики:

(14).

Особенности:

Решение уравнения динамики для pв дает выражение для интегрального
звена.

Если учесть выражение Gyn = f (pв ), то звено получается апериодическим
1 порядка.

Gyк = f (Gхл ), можно получить на основании теплового баланса
конденсатора:

.
(15).

На основании (13), (14) и (15) можно принять:Pв =f(Gхл).

Информационная схема верха колонны.

Рис.8.

Информационная схема верха колонны как многосвязного объекта по (в и pв

Рис.9а.

Информационная схема верха колонны как многосвязного объекта по hфл и (в

Рис.9б.

Информационная схема колонны как многосвязного объекта по (в и (н.

Рис.9в.

Математическое описание подогревателя потока питания.

Тепловой баланс

Уравнение динамики.

(16).

Уравнение статики.

Gт Cрт ?твх – Gт Cрт ?твых = Gп Cpп ?n0 – Gп Cpп ?n ,
(17)

Обозначения:

?т вх ,??т вых , ?n 0 ,?n – температуры потоков теплоносителя и питания
на входе и выходе из теплообменника;

Vn – объем потока питания в трубах теплообменника;

Cpп, Срт – удельные теплоемкости потоков питания и теплоносителя;

Gт, Gп – массовые расходы теплоносителя и питания, кг/ч.

На основании (16) и (17) можно считать:

.

Предпочтительное управляющее воздействие Gт.

Информационная схема подогревателя потока питания

как объекта управления ?п0

Рис.10.

Типовая схема автоматизации процесса ректификации.

Рис.11.

Типовое решение автоматизации

процесса ректификации.

Регулирование.

Регулирование ?в=f(cд) по подаче флегмы – косвенное регулирование
показателя эффективности процесса cд.

Регулирование Pв по подаче хладоагента Gхл – обеспечивает материальный
баланс по паровой фазе.

Регулирование hфл по отбору флегмы Gфл – обеспечивает материальный
баланс по жидкой фазе верха колонны.

Регулирование hк по отбору кубового продукта Gк – обеспечивает
материальный баланс по жидкой фазе низа колонны.

Стабилизация расхода питания Gп – обеспечивает:

материальный баланс по всему веществу,

снятие наиболее существенных возмущений,

заданное положение рабочей области колонны;

стабилизацию производительности установки.

Стабилизация расхода греющего пара Gгр – обеспечивает:

тепловой баланс установки;

стабилизацию Gy0 .

Регулирование ?п0 по подаче Gт обеспечивает:

заданное положение рабочей линии;

эффективность процесса разделения;

тепловой баланс

Контроль.

Температуры и расходы всех исходных потоков.

Температуры – ?в, ?н, ?кв, ?кн, ?п0.

Давление – Рв, Рн.

Уровень – hфл, hк.

Концентрации – сд или ск .

Сигнализация.

существенные отклонения hфл, hк, ?в от заданий:

;

.

7. Материалы к лекции №15

Автоматизация реакторных процессов (ч.1)

Упрощенная структурная схема химического реактора.

Рис.1.1.

Химический процесс (3) определяется:

уравнениями кинетики

взаимодействием гидродинамических, массообменных и тепловых процессов в
аппарате,

от которых зависят концентрации реагентов, температура ( и давление Р
реакторного процесса.

Химические превращения (3) приводят к изменению тепловых (2) и
гидродинамических процессов (1) в реакторе.

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

:

процесс идет в кинетической области;

скорость процесса определяет химическое взаимодействие;

массообмен не влияет на скорость химической реакции.

:

процесс идет в диффузионной области;

процесс характеризуется массообменом;

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

:

процесс идет в переходной области;

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

Показатели эффективности реакторного процесса.

Степень превращения.

Степень превращения Un , представленная через мольные доли:

(1а),

где n0 – число молей компонента в исходном потоке;

n – число молей компонента в реакционной смеси.

Степень превращения Um , представленная через массовые доли:

(1б),

где m0 – масса компонента в исходном потоке;

m – масса компонента в реакционной смеси.

Факторы, влияющие на степень превращения:

( и Р – температура и давление, влияют на смещение химического
равновесия в реакции;

tр – продолжительность химической реакции;

С0 – концентрации исходных веществ;

подбор катализатора;

величина потока рециркуляции.

Выход продукта.

Определение выхода продукта Х:

(2),

где

mф – масса фактически полученного продукта;

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

Факторы, влияющие на выход продукта Х:

? – температура;

Р – давление;

Сi – состав реакционной смеси;

tпр – время пребывания реакционной смеси в аппарате.

Выход продукта характеризует:

степень совершенства технологического процесса: чем ближе Х ( 1, тем
ближе расходные коэффициенты к стехиометрическим;

экономические показатели технологического процесса: чем ближе Х ( 1,
тем лучше экономические показатели реакторного процесса.

Избирательность химического процесса.

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

(3).

где

– количество молей исходного продукта, превращенных в целевой
продукт;

– количество молей химически превращенных исходных веществ.

Избирательность влияет на экономические показатели процесса.

Скорость химического процесса.

Скорость химического процесса – это количество вещества, которое
реагирует или образуется в единицу времени в единице объема (или на
единице поверхности):

(4),

где

– движущая сила процесса,

определяемая для химических реакций как произведение концентраций
компонентов в степенях, равных их стехиометрическим коэффициентам;

К – константа скорости реакции,

определяемая на основании уравнения Аррениуса.

Скорость процесса характеризует производительность химического реактора.

Основные факторы, влияющие на скорость реакции:

С0 – состав исходных реагентов;

? – температура;

Р – давление.

Схема реактора непрерывного действия с мешалкой

с экзотермической реакцией 1-го порядка А ( В.

Рис.4.1.

Показатель эффективности реакторного процесса в общем случае –
концентрация целевого продукта в реакционной смеси СВ .

Цель управления в общем случае:

.

Цель управления для данного процесса:

.

Обозначения на рис.4.1:

G0 , G , Gхл – массовые расходы исходного реагента, реакционной смеси и
хладоагента, кг/с;

cp0 , cp , cpхл – удельные теплоемкости соответствующих потоков,
дж/(кг*град);

( – плотность реакционной смеси, кг/м3 ;

(0 ,( , (хл – температуры исходного реагента, реакционной смеси и
хладоагента;

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

V – объем реакционной смеси, м3 ;

h – уровень реакционной смеси, м;

-(H – тепловой эффект экзотермической реакции, дж/кг;

r – скорость химической реакции, кг/(м3*с).

Математическое описание реактора

на основе физико-химических особенностей процесса.

Стехиометрическое уравнение химической реакции.

Стехиометрическое уравнение химической реакции, характеризующее ее
материальный баланс:

(1а).

Стехиометрическое уравнение для данного реактора:

(1б),

где (i , (i , ( , ( – стехиометрические коэффициенты, числа
реагирующих молей компонентов процесса.

Уравнение кинетики химической реакции.

Результирующая скорость обратимой реакции:

(2),

где

– общий порядок реакции;

– скорость прямой реакции;

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

– константы скорости,

ZА и ZВ – эмпирические коэффициенты;

ЕА и ЕВ – энергии активации.

Кинетика для реактора типа рис.4.1:

скорость прямого процесса:

(3а);

скорость обратного процесса:

(3б).

Уравнение гидродинамики процесса.

Уравнение гидродинамики процесса характеризует давление в потоке как
сложную функцию:

(4)

где

x, z – координаты потока;

( – ламинарная вязкость, характеризующая силы внутреннего трения;

k – турбулентная вязкость;

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

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

t – текущее время.

Материальный баланс реакторного процесса

по всему веществу.

Уравнение динамики:

(1),

где

:

(2).

На основе(1) и (2):

(3).

Материальный баланс реакторного процесса

по расходуемому веществу А.

Уравнение динамики:

(4),

где

I – изменение количества вещества А в реакторе в единицу времени;

II – количество вещества А, поступающего в реактор в единицу времени;

III – количество вещества А, отводимого из реактора в единицу времени;

IV – количество вещества А, расходующегося в реакторе на химическую
реакцию в единицу времени, где V=S*h.

:

(5).

На основе (4) и (5):

(6).

Тепловой баланс реакторного процесса.

Уравнение динамики:

(7),

где

I – изменение количества тепла в реакторе в единицу времени;

II – количество тепла, поступающего в реактор с исходным реагентом в
единицу времени;

III – количество тепла, отводимого из реактора с реакционной массой в
единицу времени;

IV – количество тепла, выделяющегося в реакторе в результате химической
реакции в единицу времени;

V – количество тепла, отводимого из реактора с хладоагентом в единицу
времени.

:

(8).

На основании (7) и (8):

(9).

Информационная схема реактора на основе балансовой модели.

Рис.6.1.

.

.

.

.

8. Материалы к лекции №16

Автоматизация реакторных процессов (ч.2)

реакторного процесса в адиабатическом режиме.

Рис.1.1.

Влияние изменения температуры (0 на температуру реакторного процесса (.

При изменении входной температуры (0 в пределах от (02 до (04 работа
реактора осуществляется по петле гистерезиса D-A-B-C-D.

Участок статической характеристики А-3-С – путем изменения входной
температуры не реализуется.

Возможные стационарные состояния реактора.

При (01 – одно устойчивое состояние в (()1, но температура в реакторе
низкая и реакция идет слабо.

При (02 – два стационарных состояния в (()D и (()С:

(()D – устойчивое состояние, но при низкой температуре;

(()С – неустойчивое состояние с возвратом в (()D при уменьшении
начальной температуры.

При (03-три стационарных состояния в (()2, (()3, (()4:

(()2 – устойчивое состояние, но с низкой скоростью реакции;

(()4 – устойчивое состояние, почти полное превращение реагентов, но
реакция идет с высокой температурой – вне рабочего диапазона;

(()3 – изменением начальной температуры не реализуется.

При (04 – два стационарных состояния в (()А и (()В:

(()А – неустойчивое состояние, с проскоком в (()В при увеличении
начальной температуры;

(()В – устойчивое состояние, почти полное превращение реагентов, но
реакция идет с высокой температурой – вне рабочего диапазона.

При (05 – одно устойчивое состояние в (()5, но температура в реакторе
вне рабочего диапазона.

от концентрации С0 исходного реагента.

Рис.1.2.

При больших значениях С0 (С01 ) – статическая характеристика
неоднозначна даже при низких температурах.

Единственность стационарного состояния возможна только при высоких
значениях (0 .

Область отрицательных значений (0 практического смысла не имеет, поэтому
показана пунктиром.

Уменьшение входных концентраций (С02 и С03) вызывает смещение
неоднозначности статической характеристики в область более высоких
температур.

от концентрации С0 исходного реагента.

Рис.1.2.

При больших значениях С0 (С01) – статическая характеристика
неоднозначна даже при низких температурах.

Единственность стационарного состояния возможна только при высоких
значениях (0 .

Область отрицательных значений (0 практического смысла не имеет, поэтому
показана пунктиром.

Уменьшение входных концентраций (С02 и С03) вызывает смещение
неоднозначности статической характеристики в область более высоких
температур.

Определение стационарных состояний

по диаграмме «выделения – отвода тепла».

Рис.2.1.

1, 2, 3 – кривые отвода тепла Q при изменении температуры в реакторе (
c различными значениями коэффициента теплопередачи (: (1 = (3 и (2 Зависимость характеристики выделения тепла от времени пребывания реакционной массы в реакторе. Рис.2.2. 1, 2, 3 - характеристики выделения тепла при Тср1 > Тср2 > Тср3 .

4 – характеристика отвода тепла.

, где Vрм и Qрм – объем и объемный расход реакционной массы.

При увеличении времени пребывания характеристика выделения тепла
смещается влево.

Система 1-4:

одно устойчивое стационарное состояние в (()Е;

практически полное превращение реагента;

но очень высокая температура, которая может быть вне рабочей зоны
реактора;

Система 3-4.

одно устойчивое стационарное состояние в (()А;

очень низкая температура в реакторе;

низкая скорость реакции.

Система 2-4.

три стационарных состояния в (()В , D и С;

(()В – очень низкая температура в реакторе; низкая скорость реакции,
стационарное состояние устойчивое;

(()D – практически полное превращение реагента; но очень высокая
температура, которая может быть вне рабочей зоны реактора; стационарное
состояние устойчивое;

(()С – температура в рабочей зоне реактора; но стационарное состояние –
неустойчивое.

Оценка устойчивости стационарных состояний

по диаграмме «выделения – отвода тепла».

Qр = f(() – характеристика выделения тепла реакции;

Qт = f(() – характеристика отвода тепла.

Рис.2.3.

Устойчивость реактора в стационарном состоянии В:

;

;

стационарное состояние в (()В – устойчивое.

Устойчивость реактора в стационарном состоянии D:

;

;

стационарное состояние в (()D – устойчивое.

Устойчивость реактора в стационарном состоянии C:

переход в(()D;

переход в (()В;

стационарное состояние в (()С – неустойчивое.

9. Материалы к лекции №17

Автоматизация реакторных процессов (ч.3).

Схема для реакции типа А(ж) + В (г)( D (г-ц.пр)+С (ж).

Рис.1.1.

Схема для реакции типа А(ж) + В (г)( D (г)+С (ж-ц.пр).

Рис.1.2.

Схема автоматизации для реакции типа

А(ж) + В{Всв(ж) + Врец(ж)} ( С(ж- ц.пр)+В(ж).

1- реактор; 2 – сепаратор; 3 – сборник; 4, 5 – насосы;

А, В – исходные реагенты; С – целевой продукт; Врц – возвратный
(рецикловый) поток реагента В; Всв – свежий поток реагента В.

Рис.1.3.

Система регулирования концентрации

из двух одноконтурных АСР.

Рис.1.5.

Каскадная АСР концентрации со вспомогательным контуром

стабилизации соотношения расходов исходных реагентов.

Рис.1.6.

Каскадная АСР концентрации со вспомогательным контуром

стабилизации температуры исходного реагента.

Рис.1.7.

Каскадная АСР концентрации со вспомогательным контуром

стабилизации температуры хладоагента.

Рис.1.8.

Одноконтурная АСР температуры в реакторе

по подаче исходного реагента.

Рис.2.1.

Каскадная АСР температуры в реакторе со вспомогательным контуром

стабилизации температуры исходного реагента.

Рис.2.2.

Каскадная АСР температуры в реакторе со вспомогательным контуром

стабилизации температуры хладоагента.

Рис.2.3.

Каскадная АСР температуры в реакторе со вспомогательным контуром

стабилизации давления хладоагента.

Рис.2.4.

10. Рекомендуемая литература

Основная литература

1. Дудников Е.Г. Автоматическое управление в химической
промышленности. – М.: Химия, 1987.- 368 с.

2. Полоцкий Л.М., Лапшенков Г.И. Автоматизация химических производств.
– М.: Химия, 1982.- 295 с.

3. Автоматизация технологических процессов легкой промышленности: Учеб
пособие для вузов по спец. «Автоматизация технологических процессов и
производств» / Под ред. Л.Н. Плужникова. – 2-е изд., перераб. и доп. –
М.: Легпромбытиздат, 1984.- 366с.

4. Мамиконов А.Г. Проектирование АСУ.- М.: Высшая школа, 1987.- 303 с.

5. Стефани Е.П. Основы построения АСУ ТП.- М.: Энергоиздат, 1982.- 352с.

6. Пиггот С.Г. Интегрированные АСУ химических производств. – М.: Химия,
1985.- 410 с.

7. Кафаров В.В., Макаров В.В. Гибкие автоматизированные системы в
химической промышленности: Учебник для вузов. – М.: Химия, 1990.- 320с.

8. Плютто В.П. Управление химико-технологическими процессами. Процессы
массообмена: [Учеб. пособие].- М.: МХТИ, 1984.-48с.

9. Плютто В.П. и др. Автоматизированные системы управления
периодическими пролцессами химической технологии. – М.: МХТИ, 1985.-48с.

10. Ицкович Э.Л., Сорокин Л.Р. Оперативное управление непрерывным
производством. – М..: Наука, 1989.-155с.

Дополнительная литература

11. Автоматизация технологических процессов легкой промышленности: Учеб
пособие для вузов по спец. «Автоматизация технологических процессов и
производств» / Под ред. Л.Н. Плужникова. – 2-е изд., перераб. и доп. –
М.: Легпромбытиздат, 1993 .

12. Уланов Г.М. и др. Методы разработки интегрированных АСУ
промышленными предприятиями. – М.: Энергоатомиздат, 1983.- 320 с.

13. Практикум по автоматике и системам управления производственными
процессами: Учеб. пособие для вузов / Под ред. И.М. Масленникова – М.:
Химия,1986.- 336 с.

14. Плютто В.П., Путинцев В.А., Глумов В.М. Практикум по теории
автоматического управления химико-технологическими процессами: Цифровые
системы – М.: Химия,1989.- 168 с.

15. Голубятников В.А., Шувалов В.В. Автоматизация производственных
процессов в химической промышленности. – М.: Химия, 1985.-352с.

16. Шувалов В.В. Огаджанов Г.А., Голубятников В.А. Автоматизация
производственных процессов в химической промышленности. – М.: Химия,
1991.-480с.

17. Процессы и аппараты химической технологии. Учебник для вузов./ Н.Н.
Смирнов, М.И. Курочкина, А.И. Волжинский, В.А. Плессовских. – СПб.:
Химия, 1996.-400с.

18. Математическое моделирование основных процессов химических
производств. Учеб. пособие для вузов. / В.В. Кафаров, М.Б. Глебов. –
М.:Высш.шк., 1991.-399с.

19. Математическое моделирование химико-технологических систем: Учеб.
пособие в 3ч./ под ред Л.С. Гордеева. – М.:РХТУ, 1999- 48с.( ч1); 47с.(
ч2), 67с.( ч3)

.

СОДЕРЖАНИЕ

TOC \o “1-1” \h \z HYPERLINK \l “_Toc48800337” 1. Материалы к
лекции №9 PAGEREF _Toc48800337 \h 3

HYPERLINK \l “_Toc48800338” 2. Материалы к лекции №10 PAGEREF
_Toc48800338 \h 8

HYPERLINK \l “_Toc48800339” 3. Материалы к лекции №11 PAGEREF
_Toc48800339 \h 17

HYPERLINK \l “_Toc48800340” 4. Материалы к лекции №12 PAGEREF
_Toc48800340 \h 25

HYPERLINK \l “_Toc48800341” 5. Материалы к лекции №13 PAGEREF
_Toc48800341 \h 34

HYPERLINK \l “_Toc48800342” 6. Материалы к лекции №14 PAGEREF
_Toc48800342 \h 41

HYPERLINK \l “_Toc48800343” 7. Материалы к лекции №15 PAGEREF
_Toc48800343 \h 50

HYPERLINK \l “_Toc48800344” 8. Материалы к лекции №16 PAGEREF
_Toc48800344 \h 57

HYPERLINK \l “_Toc48800345” 9. Материалы к лекции №17 PAGEREF
_Toc48800345 \h 62

HYPERLINK \l “_Toc48800346” 10. Рекомендуемая литература PAGEREF
_Toc48800346 \h 68

PAGE

PAGE 69

Нашли опечатку? Выделите и нажмите CTRL+Enter

Похожие документы
Обсуждение

Ответить

Курсовые, Дипломы, Рефераты на заказ в кратчайшие сроки
Заказать реферат!
UkrReferat.com. Всі права захищені. 2000-2020