Учебное пособие 800429
.pdfПри выбранных параметрах зависимость массы загрузки от времени режима осушки и влагосодержания при заданных параметрах приведена в таблице 2.
Таблица 2 Зависимость массы загрузки от времени режима осушки
t , час |
|
1 |
|
|
3 |
|
|
5 |
|
x , |
0,0135 |
0,026 |
0,04 |
0,0135 |
0,026 |
0,04 |
0,0135 |
0,026 |
0,04 |
|
|
|
|
|
|
|
|
|
|
кг/м3 |
|
|
|
|
|
|
|
|
|
|
q , |
216 |
216 |
|
216 |
|
216 |
|
|
|
216 |
|
216 |
|
216 |
216 |
216 |
||||||
м3/сут |
|
|
|
|
|
|
|
||||||||||||||||
|
Мад , |
2,43 |
4,68 |
|
7,2 |
|
7,29 |
|
|
|
14,01 |
|
21,6 |
|
12,15 |
23,4 |
36 |
||||||
|
кг |
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V , м3 |
0,362 |
0,698 |
|
1,075 |
|
1,08 |
|
|
|
2,091 |
|
3,22 |
|
1,813 |
3,49 |
5,37 |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t , час |
|
10 |
|
|
|
|
|
|
|
|
|
14 |
|
|
|
|
16 |
|
||||
|
x , |
0,00016 |
|
0,00024 |
|
0,0004 |
|
0,00016 |
|
|
|
0,00024 |
|
0,0004 |
|
0,00016 |
|
0,00024 |
0,0004 |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
кг/м3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
q , |
216 |
|
216 |
|
216 |
|
216 |
|
|
|
216 |
|
216 |
|
216 |
|
216 |
216 |
||||
м3/сут |
|
|
|
|
|
|
|
|
|
||||||||||||||
|
Мад , |
0,428 |
|
0,428 |
|
0,72 |
|
0,4 |
|
|
|
0,599 |
|
1 |
|
0,461 |
|
0,685 |
1,152 |
||||
|
кг |
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V , м3 |
0,043 |
|
0,064 |
|
0,107 |
|
0,06 |
|
|
|
0,089 |
|
0,149 |
|
0,0688 |
|
0,102 |
0,172 |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
Проведём гидравлический расчет адсорбера. |
|
|
|
|
||||||||||||||||||
|
Число Рейнольдса Re, определяют по формуле |
|
|
|
|
||||||||||||||||||
|
|
|
|
Re |
0.45 |
|
|
|
|
|
W1ср |
dэ |
, |
|
|
|
(2) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
, |
|
|
|
, |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
где |
dэ – |
средний |
диаметр |
частицы, dэ |
= |
0,0036 |
м [2]; |
|||||||||||||||
, |
– пористость, , = 0,49 [2]; – кинематическая вязкость, м2/с. |
11
|
|
|
|
|
2 Твых |
|
Твх |
, |
|
|
|
|
|
|
(3) |
|||||||||||||||
|
|
f |
|
|
|
|
|
|
Тср |
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
где Тср |
– средняя эффективная температура газа в адсорбере, К; |
|||||||||||||||||||||||||||||
Твх – температура на входе |
|
|
в адсорбер, Твх |
= |
293 К; |
|
Твых – |
|||||||||||||||||||||||
температура на выходе из адсорбера, Твых |
|
= 433 К. |
|
|
|
|||||||||||||||||||||||||
Среднюю температуру в адсорбере Тср.ад ., К определяют по |
||||||||||||||||||||||||||||||
формуле |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Т |
|
|
Твх |
Твых |
. |
|
|
|
|
|
|
|
(4) |
||||||||||||||
|
|
ср |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Среднюю эффективную температуру газа в адсорбере Т , K |
||||||||||||||||||||||||||||||
вычисляют по формуле |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
Т |
|
|
|
|
|
|
Твх Твых |
. |
|
|
|
|
|
|
|
(5) |
|||||||||||
|
|
ср.эф |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
Тср |
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Среднюю эффективную плотность газа в адсорбере |
|
, кг/м3 |
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ср.ад. |
|
|
вычисляют по формуле |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
273 |
|
, |
|
|
|
|
|
|
|
(6) |
||||||||
|
|
ср.ад. |
|
0 |
Tср |
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
где |
0 |
– плотность газа на входе адсорбера, кг/м3. |
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Плотность газа на входе адсорбера |
|
|
|
0 |
, кг/м3 |
определяют по |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
формуле |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
104 |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
(7) |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
R Tвх |
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Принимая время режима |
|
|
|
осушки |
|
t |
|
= |
|
12 |
часов, |
диаметр |
||||||||||||||||||
адсорбера D = 144 мм. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Определяют скорость газа на входе в адсорбер W1 , м/с по |
||||||||||||||||||||||||||||||
формуле |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
W1 |
|
|
|
|
|
|
|
|
m R |
Tвх |
|
|
|
|
|
. |
|
|
|
(8) |
||||||||
|
|
|
|
|
104 0,7854 0,1442 |
|
|
|
вх
Среднюю эффективную скорость газа в адсорбере W1ср , м/с вычисляют по формуле
12
W |
W |
Т |
ср |
W |
Твх |
Твых |
. |
(9) |
|
|
|
|
|
|
|||||
1ср |
1 |
Т |
вх |
1 |
2Т |
вх |
|
|
Подставляя значения в формулы (2) – (9) получаем следующее:
|
|
|
|
W1 |
|
0,00165 52,8 300 |
|
0,082 м/с; |
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
2 104 0,7854 0,1442 |
|
||||||||||||||||||||||||
|
|
|
|
W |
0,082 293 433 0,101 м/с; |
|
|||||||||||||||||||||||
|
|
|
|
|
|
1ср |
|
|
|
|
|
|
|
2 293 |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Т |
|
|
|
293 433 |
363 |
К; |
|
|
|
|
|
2 104 |
1, 29 |
кг/м3; |
|||||||||||||||
ср.ад. |
|
|
2 |
|
|
|
|
0 |
52,8 293 |
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
ср.ад. |
1,29 |
273 |
1,041 |
кг/м3; |
f |
|
|
2 293 433 |
0,771 ; |
||||||||||||||||||||
363 |
363 |
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
Re |
|
|
0.45 |
|
|
|
|
|
|
0,101 0, 0036 |
33, 25 ; |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
1 0, 49 |
|
|
|
|
|
14,71 10 6 |
|||||||||||||||||||||
|
|
0, 49 |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
75 |
|
|
|
|
15 |
|
|
1 |
3,515 ; |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
33,25 |
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
33,25 |
|
|
||||||||||||||||||||
|
|
|
30,6 3,515 |
1, 4614 |
0,771 436629,8 . |
||||||||||||||||||||||||
|
|
|
0,0036 |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, перепад давления на адсорбере при режиме регенерации р , МПа может составить
|
|
|
|
|
|
|
|
|
|
W 2 |
|
||
|
|
|
|
|
|
|
|
ср |
1ср |
, |
(10) |
||
|
|
р |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
2 g |
|
||||
|
|
|
|
|
|
|
|
|
|
||||
|
|
436629,8 |
1, 041 0,1012 |
|
236,32 кгс/м2 = 0,0024 МПа. |
||||||||
р |
2 9,81 |
||||||||||||
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|||||
На |
|
режиме осушки |
|
эта |
|
величина |
составляет примерно |
0,08 МПа.
Литература
1.Коршак А.А. Основы нефтегазового дела / А.А. Коршак, А.М. Шаммазов. – Уфа: ДизайнПолиграфСервис, 2001. – 249 с.
2.Юрьев А.С. Справочник по расчётам гидравлических и вентиляционных систем / Юрьев А.С., Пирогов С.Ю. и др – С.-Пб: АНО НПО «Мир и семья», 2001. – 1154 с.
Воронежский государственный технический университет
13
УДК 620.9
Р.А. Константинов – магистрант, В.Ю. Дубанин, канд. техн. наук, доц. – научный руководитель
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ТЕПЛОЭНЕРГЕТИЧЕСКОЙ СИСТЕМЫ, СОСТОЯЩЕЙ ИЗ ПОСЛЕДОВАТЕЛЬНО СОЕДИНЕННОЙ ГИДРОДИНАМИЧЕСКОЙ СТРУКТУРЫ
Рассмотрена теплоэнергетическая система, состоящая из последовательно соединенных двух зон идеального перемешивания и зоны идеального вытеснения. Получено аналитическое решение математической модели системы, выраженной в виде систем дифференциальных уравнений первого порядка, позволяющее оценить влияние различных параметрических характеристик объектов системы на изменение температуры теплоносителя при выходе его из системы
Ключевые слова: математическая модель, теплоэнергетическая система, теплоноситель, гидродинамическая структура
Теплоэнергетические системы, представляющие собой последовательное соединение зон перемешивания и зон вытеснения, получили широкое распространение в химической и нефтехимической промышленности. Ярким примером таких систем является реакторное оборудование химических производств. Как правило, в целях изучения процессов, протекающих в данных системах, используют их идеальные варианты – аппараты идеального вытеснения и идеального перемешивания [1]. Идеализация протекающих процессов позволяет исключить влияние менее значимых факторов, что в значительной степени облегчает процесс моделирования таких систем.
Рассмотренная система может быть представлена в виде схемы, показанной на рис. 1. Задача моделирования состоит в нахождении температуры теплоносителя на выходе из системы Tвых (t) и
определении влияния технологических параметров (объема V аппарата, расхода U теплоносителя, начальной температуры теплоносителя Tвх (t) и т.п.) на конечный результат.
При математическом моделировании системы, приведенной на рис. 1, будем рассматривать каждый ее объект в отдельности друг от друга. При этом изменение температуры в зоне идеального перемешивания представим в виде системы дифференциальных уравнений первого порядка [2].
14
Рис. 1. Теплоэнергетическая система, состоящая из последовательно соединенных зон идеального перемешивания и идеального вытеснения
Для первого аппарата:
|
|
dT1 (t1 ) |
|
|
1 |
T (t ) Т (t ) ; |
||||||||
|
|
|
|
|
|
|||||||||
|
dt1 |
|
|
|
|
вх 1 |
1 |
1 |
|
|||||
|
|
|
1 |
|
|
|
(1) |
|||||||
|
|
|
|
Т1 (0) T0 . |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|||||||
Для второго аппарата: |
|
|
|
|
|
|
|
|
|
|||||
|
|
dT2 (t2 ) |
|
|
1 |
T (t ) Т |
|
(t |
|
) ; |
||||
|
|
|
|
|
2 |
2 |
||||||||
|
dt2 |
|
1 1 |
|
|
|||||||||
|
|
|
2 |
|
|
|
(2) |
|||||||
|
|
|
|
Т2 (0) T1 (t1 ). |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|||||||
где 1 V1 / U, |
2 |
V2 / U – время пребывания теплоносителя в |
аппаратах 1 и 2 соответственно, с.
Система уравнений, описывающих изменения температуры в зоне идеального вытеснения, выглядит следующим образом
|
T3 |
(z,t3 ) |
u |
T3 |
(z,t3 ) |
; |
|
dt3 |
|
dz |
|||
|
|
|
||||
|
|
|
|
|||
|
|
Т3 (0) T0 ; |
|
(3) |
||
|
|
|
||||
|
|
Т3 (0,t) T2 (t2 ). |
|
|||
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где u U / S3 – скорость движения теплоносителя в аппарате 3,
м/с.
При этом очевидно, что температура на выходе из системы будет равна изменению температуры в третьем аппарате по его длине, то есть
Tвых (t) T3 (l3 ,t),
15
где l3 V3 / S3 – длина аппарата 3.
Для решения систем 1-3 необходимо привести их к безразмерному виду. Для этого введем следующие поправочные
коэффициенты: |
|
t / 1 – |
безразмерное |
время; |
1 / 2 – |
коэффициент, |
учитывающий |
изменение |
времени |
во втором |
|
аппарате; |
u 1 |
коэффициент, учитывающий |
изменение |
||
l3 |
скорости по сечению третьего аппарата.
Сучетом введенных безразмерных коэффициентов системы (1)-
(3)примут вид:
|
d 1 ( ) |
вх ( ) 1 ( ) ; |
|
|||||||||||
|
|
|
|
|
|
|
||||||||
|
|
d |
(4) |
|||||||||||
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
1 (0) 0. |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
||||||
|
d 2 ( ) |
1 ( ) 2 |
( ) ; |
|
||||||||||
|
|
|
|
|
|
|
||||||||
|
|
d |
(5) |
|||||||||||
|
|
|
|
|
|
|
|
|
||||||
|
|
2 (0) 1 ( ). |
|
|
|
|||||||||
|
|
|
|
|
||||||||||
|
3 (Z , ) |
|
3 (Z , ) |
; |
|
|||||||||
|
|
|
|
|
Z |
|
|
|||||||
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
3 (0) 0; |
|
|
|
(6) |
||||||
|
|
|
|
|
|
|
|
|||||||
|
|
3 |
(0, ) 2 ( ). |
|
|
|
||||||||
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
где |
Z z / (u 3 ) ; |
1 |
(T1 T0 ) / T0 ; |
2 (T2 T0 ) / T0 ; |
3 (T3 T0 ) / T0 .
По определению передаточной функции [1] имеем:
вых |
|
вх |
|
|
|
|
W (Р) L( |
) / L( |
) , |
(7) |
|||
где W (Р) – передаточная |
|
функция |
|
всей системы; |
L – |
преобразование Лапласа.
Передаточная функция системы при параллельном соединении ее элементов, определится как произведение передаточных функций отдельных элементов системы, т.е.
W (Р) W1 (Р) W2 (Р) W3 (Р), |
(8) |
где W1 (Р) , W2 (Р) , W3 (Р) – передаточные функции аппаратов 1,
2 и 3 соответственно. При этом имеем:
16
W (Р) |
1 |
|
, W (Р) |
|
, W (Р) е( Р/ ) . |
|
|
|
|||
|
|
|
|||
1 |
Р 1 |
2 |
Р |
3 |
|
|
|
|
Подставляя данные значения в (8), получим выражение передаточной функции для всей системы
W (Р) |
е( Р/ ) |
(9) |
. |
(Р 1)(Р )
Выражая вых ( ) из (7) получим дифференциальное уравнение
для нахождения температуры теплоносителя на выходе из системы (в безразмерном виде)
|
вых |
( ) L 1 W (P) L |
( ) , |
|
|
|
вх |
|
где L – обратное преобразование Лапласа.
Подставляя (9) в полученное выражение, придем к окончательному виду уравнения для нахождения искомой величины
|
|
|
е( Р/ ) |
|
|
|
|
|
вых |
( ) L 1 |
|
|
L вх |
( ) |
, |
(10) |
|
|
||||||||
|
|
|
(Р 1)(Р ) |
|
|
|
|
|
Связь Tвх (t) с |
вх ( ) и |
Tвых (t) с вых ( ) |
выглядит следующим |
образом:
вх ( )
вых ( )
В итоге имеем
Твх (t) Т0 ;
Твых (t) Т0 .
|
е( Р/ ) |
|
|
|
Твых (t) L 1 |
|
L Tвх (t) T0 T0 , |
(11) |
|
|
||||
|
(Р 1)(Р ) |
|
|
|
Зная Tвх (t) , T0 , U , |
V и S с помощью выражения (11) можно |
|||
легко найти Tвых (t) . |
|
|
|
|
Так для системы, изображенной на |
рис.1, при T0 |
= 20 С; |
||
U = 0,015 м3/с; V = 0,1 м3; V = 0,2 м3; |
V = 0,3 м3; S |
= 0,3 м2 и |
||
1 |
2 |
|
3 |
|
изменении температуры теплоносителя на входе в систему по закону Tвх (t) 40(1 e 0,1t ), используя уравнение (11), получим графическое изображение зависимости изменения температуры теплоносителя от
времени на выходе из системы Tвых (t) (рис. 2). Изменяя входные параметры системы, можно получить наилучшее их соотношение с требуемым выходным значением температуры теплоносителя.
17
Рис. 2. Зависимость температуры теплоносителя от времени на выходе из системы, приведенной на рис. 1
Результатом проведенного в данной работе математического моделирования теплоэнергетической системы с последовательным соединением элементов (рис. 1) стала возможность разработки алгоритма наглядной оптимизации теплоэнергетических параметров объектов системы, позволяющая достичь наилучших техникоэкономических показателей системы.
Литература
1.Бондарь А.Г. Математическое моделирование в химической технологии / А.Г. Бондарь. – Киев: «Вита школа», 1973. - 280 с.
2.Цветков Ф.Ф. Тепломассообмен: учебник для вузов / Ф.Ф. Цветков, Б.А. Григорьев. - М: Издательский дом МЭИ, 2011. - 562 с.
3.Берд Р. Явления переноса / Р. Берд, В. Стьюарт, Е. Лайтфут. - М.: «Химия», 1974. - 688 с.
Воронежский государственный технический университет
18
УДК 620.9
С.А. Воробьев – магистрант, Ю.Н. Агапов, д-р техн. наук, проф. – научный руководитель
АНАЛИЗ ВЛИЯНИЯ ВЗАИМНОГО РАСПОЛОЖЕНИЯ АППАРАТОВ ИДЕАЛЬНОГО ПЕРЕМЕШИВАНИЯ ПРИ НАЛИЧИИ ЗАСТОЙНОЙ ЗОНЫ
В данной работе рассмотрена теплоэнергетическая система, состоящая из двух последовательно соединенных зон идеального перемешивания, одна из которых осложнена наличием застойной зоны. Было получено аналитическое решение данной системы, математическая модель которой представлена в виде системы дифференциальных уравнений 1-ого порядка. Вычислительный эксперимент позволяет оценить влияние различных параметрических характеристик объекта на структуру изменения температуры теплоносителя на выходе из системы
Ключевые слова: математическая модель, теплоэнергетическая система, застойная зона, вычислительный эксперимент
Данная система может быть полезна для использования в химической промышленности.
Рассмотрим данный объект, который изображен на рис. 1.
Рис. 1. Теплоэнергетическая система
19
1. Составим математическую модель системы
dТ1 (t) |
|
|
|
1 |
|
Твх |
(t) |
Т1 (t) ; |
||||||
|
|
|
|
|
|
|
|
|
|
|||||
|
|
dt |
1 |
|||||||||||
|
|
|
|
|
|
|
||||||||
Т |
|
(0) T |
|
|
|
|
|
|
|
|||||
|
1 |
|
|
|
0 |
|
|
|
|
|
|
|||
dТ2 (t) |
|
|
|
1 |
|
Т1 |
(t) (1 )Т2 (t) T3 (t) ; |
|||||||
|
|
|
|
|
|
|
|
|
|
|||||
|
|
dt |
|
|
|
|
||||||||
|
|
|
|
2 |
|
|
|
|
||||||
dТ3 (t) |
|
|
|
|
Т2 |
|
|
Т3 (t) |
||||||
|
|
|
|
|
|
|
|
|
|
(t) |
||||
|
|
dt |
|
|
|
|
||||||||
|
|
|
|
3 |
|
|
|
|
||||||
Т |
2 |
(0) Т |
3 |
(0) Т |
0 |
. |
||||||||
|
|
|
|
|
|
|
|
|
|
|
2.Приведём систему к безразмерному виду
t1 , 1 Т1 ( ) Т0 , вх ( ) Твх Т0 ,
1 .
2
Запишем систему, применив введённые параметры
d 1 ( ) |
вх ( ) 1 ( ) ; |
||||||
|
|
d |
|
||||
|
|
|
|
||||
|
|
|
|
|
|
||
1 (0) 0; |
|
||||||
d ( ) |
|
1 ( ) (1 ) 2 ( ) 3 ( ) ; |
|||||
|
|
|
2 |
|
|||
|
|
d |
|
|
|||
2 (0) 0; |
|
||||||
|
d 3 ( ) |
|
|
1 2 ( ) 3 ( ) ; |
|||
|
|
|
|||||
|
|||||||
|
|
d |
|
3 |
|||
|
|
|
(0) 0. |
|
|||
|
|
|
|||||
|
3 |
|
|||||
|
|
|
|
|
(1)
(2)
3. Решение системы:
Для аппарата идеального перемешивания без застойной зоны уравнение передаточной функции запишем в виде:
w1 (р) |
|
, |
(3) |
|
|
||||
р |
||||
|
|
|
а передаточная функция для аппарата идеального перемешивания с застойной зоной имеет вид:
20