Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Кафаров, В. В. Принципы математического моделирования химико-технологических систем (введение в системотехнику химических производств) учеб. пособие

.pdf
Скачиваний:
86
Добавлен:
22.10.2023
Размер:
12.55 Mб
Скачать

_ х^(і) = <рхг (і/

 

 

а

 

 

x0(s)

x,(s)

1

хг ($) l/-yj

x3(s)

о----

 

£

—•—— о

 

 

-(k+ S)

 

 

 

 

 

X0 (S}

x2(s)

n-fp)

x3(s)

о---

7T~ — 3

 

 

 

 

в

 

 

xg(sj

1

 

 

x3(s)

(i-p+fyi-s Xg(s)

 

о ■

—-

О— ■■

 

о

 

 

г

 

 

 

Xq(S)

!-<P

X3 (S) J'v

 

 

(1-p+k)+s

 

a

Рис. IV-68. Операторная схема XTC (а), соответствующий ей нормализованный сигнальный граф (б) и этапы преобразования этого графа (в — д).

превращения протекает реакция первого порядка с константой скорости к; коэффициент функциональной связи оператора разделения ф = x j x 2.

Составим систему уравнений материального баланса данной ХТС:

 

 

 

(О +

Ф^г (0 —* і ( 0 = 0

 

 

(1>

 

 

 

^ - —x1(t) + Kx2[(t) + x2 (t) = 0

 

 

(2)

 

 

 

x2 (t)—x3 (t) —(px2 (t) — 0

 

 

(3>

Используя преобразования Лапласа и преобразуя исходную систему уравне­

ний, получаем:

 

[

Х1 (і) =

Х0(^)+ фХ2 (*)

 

 

 

(l')

 

 

 

 

 

 

 

I

ЙГ2 (*)= Xi (s) —кХ% (s) —sV2 (s)

 

 

(2')■

 

"

I

X3 (S) = X 2 (S) - f X 2(S)

 

 

 

(3'>

Преобразованной системе уравнений (l')—(З') соответствует нормализо­

ванный сигнальный граф, представленный на рис. ІѴ-68, б.

 

Для определения коэффициентов передач

7\(s) =

Х г (s)/X0 (s) и Т2 (s) =

= Х 3 (s)/X0 (s) необходимо эквивалентно

преобразовать

исходный

сигнальный

граф (см. рис. ІѴ-68,

б). После преобразований (рис. ІѴ-68,

в—б) графа находим:

 

 

 

*2 («)

 

1

 

 

 

 

 

 

1 —ф+ к

 

 

 

 

 

Т1(«)= Хо(*)

1

 

 

 

 

 

 

 

 

1 —ф + К +1

 

1 —ф

 

То (*) =

Х3 (»)

Х 2 (s)

(1 —Ф)=Гх(«) (1 —Ф)= -

1— ф + к

Х 0 (*)

X 0(s)

1

+1

 

 

 

 

 

 

 

 

1 — Ф + К

 

 

 

 

 

 

 

 

 

Для расчета значений Тх =

х2 (t)/x0 (f) и

Т2 = х3 (t)/xn (t) к полученным

выражениям нужно применить обратные преобразования Лапласа.

Пример IV-23.

Определить

динамические

характеристики

колпачковой

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

на

выходе аппарата от составов и расходов жидкости и газа, поступающих

в

колонну.

 

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

функционирования абсорбционной колонны используем следующие допущения: 1) количество газа над тарелкой мало по сравнению с количеством находящейся на ней жидкости; 2) эффективность тарелки 100%; 3) соотношение между равновесными составами газа и жидкости выражается уравнением Y n = mXn + + 6; 4) количества жидкости на всех тарелках одинаковы; 5) тепловым эф­ фектом процесса абсорбции пренебрегаем.

В случае неустановившегося режима материальный баланс переносимого" компонента для n-ой тарелки может быть записан следующим образом:

(ІХ,

(1>

dt

Г (?л-1+^л+1‘—Qn—Rn

 

Все члены, вообще говоря, нелинейны. После линеаризации

они могут

быть представлены соотношениями

 

SQn . , dQn

Q n = - ^ r

 

dyn

Уп

 

 

dg

 

 

( 2>

R

n

f

dRn

 

 

 

Rn- ~ W

l~

dxn n

 

18»

где g a l — отклонения расходов газа и жидкости от установившихся значений;

Q и R — отклонение содержания

переносимого

компонента газа

и жидкости

от установившихся

значений.

 

 

 

 

 

 

Значения частных производных таковы:

 

 

 

dQn = Yn;

d Q n _n.

d R n __

t

dRn

(3 )

 

 

 

dg

 

дУп ~

dl ~

Лп ’

dxn

 

Подставляя их в уравнение (1), находим:

 

 

 

Yn-ig-^r Gyn-\-\- Zn+1Z-f-Lr„+1 Xnl —Lxn— Yng^—Gyn~ h

(4)

Из выведенного уравнения могут быть получены четыре уравнения, связы­ вающие регулируемую величину Yn с одной из независимых переменных g, I, уі и хі, смотря по тому, какая из них фактически изменяется. Когда g = I =

Уп S

х п*/ L

 

Рис. ІѴ-69. Структура потоков

------------------------П

 

id

на тарелках абсорбционной ко­

 

лонны.

хпП

Уп-І

 

= хі = 0, уравнение дает соотношение между составом Y п выходящего из ко­

лонны газа и составом поступающего в нее газа. Уравнение (4) превращается в уравнение (рис. ІѴ-69):

* = Gyn-i~\-LxnJrl Gyn^L x n (5)

Положим, что начальные условия при t = 0 будут: хп = уп = 0, а соотно­ шение между составами газа и жидкости представляет собой линейную зави­ симость. Следовательно

dYn

Уп=—^ - х п=тхп

Преобразуем уравнение (5) и подставим в него соотношение, соответствующее равновесному состоянию

Ѵп+і- (тГ+ОУп-

У п-г-Т ' dt

( 6 )

Gm

h dyn

 

Введем обозначения: Gm/L = X и hjL = а. Затем преобразуем уравнение <6), используя принятые начальные условия:

Ул+і— (1 + X-j-as) уп-\- Хул-і=0

Пусть (1 -f- X + as) = с, тогда справедливо уравнение:

Ул+і = сУп"г Хул-і= О

Для верхней тарелки колонны имеем следующее уравнение материального

баланса:

dxt

вуг Lxi DxB= h —j~-

(7)

190

Подставляя значения xD = у\, уі — тхх и D = G L в уравнение (7), получаем:

Gm

Gm

h

dyt

.

- j - У г - Ѵ і - —

Vi+ «J/i= T

- 5?-

(8>

Яу2=

( і + Я +-|-$ m'j Уі = ( c — m ) y i

(9>

На рис. IV-70, а показан сигнальный граф, соответствующий уравнениям

(7) и (9), а на рис. ІѴ-70, б—г последовательность эквивалентных преобра­

зований исходного сигнального графа (для случая, когда

число тарелок в колонне

п = 3). Для верхней тарелки абсорбера xD =

у г. На

второй тарелке (состав

1

Я

газа равен y t) в вершине у 2 имеем, что у 2 = Уі Н---- у3. После эквивалентных

- Л(с-М) я Л

 

- (с~т1 -

yl

Уз с(с-т)-А Уг (c:mj

У,

Уз cfc-mj-Л Уг

(с-пц Vt

 

 

Кіс-т)

 

 

 

 

c[c(c-mJ-A]

 

 

5

 

 

г

 

Рис. ІѴ-70. Преобразование сигнального графа, соответствующего урав­ нениям материального баланса на верхней тарелке абсорбционной колонны:

а — исходный граф; б—г — этапы преобразования графа.

преобразований

определяем,

например,

значение

передаточной функции W29:

(рис. ІѴ-70, г):

 

 

 

 

 

 

 

 

w

_

У-2, (s)

__

.

У с ________

Я (с — т)

 

23

№ (*)

 

,

Г

Я.

~|

с ( с - т ) - %

 

 

 

 

 

L с ( с — т)

J

 

Теперь выведем уравнения, связывающие состав выходящего газа с изме­ няющимися расходами. Можно предположить, что изменение расхода мгновенно проходит через колонну. Тогда ац, уі и I равны нулю, и уравнение (4) превра­ щается в следующее выражение:

Yn-ig-b Lxn~ Lxn+1— Yng— Gyn= h ^ -

(10)

Это уравнение после введения в него величин а, Я и с и преобразования Лапласа принимает вид:

-JJ (Y/i-i— Yn) éf+2/n+i+ Я?/п-і= суп

Из условия установившегося состояния вытекает соотношение

^ -(У „ -1- У я)= Я » -^ (У / - У 1)

Примем, что m/L (Y і Y x) = ß, тогда

№ng+ уru1 -!-яг/л-і = суп

(И)

19*

Уравнение (И) может быть представлено в виде сигнального графа, изображен­ ного на рис. ІѴ-71, а. Точки уп характеризуют состав газа, уходящего с ге-ой тарелки.

Сигнальный граф (рис. ІѴ-71, а) изображает процесс распространения сигналов для трехтарельчатой колонны. Изучение этого рисунка показывает, что с увеличением числа тарелок представленный граф расширяется простым добавлением подобных контуров.

 

 

 

Уі

 

 

 

 

 

 

 

 

Уз

 

 

 

 

ßA

ßA2fc*AJ

ßAfc'A)

 

 

 

 

Сг-гА+(сг-Л)(сг-2А.)

сг -Л

 

ß(cz-A)+ßA(c+fi

 

 

+ M L

 

 

с(с2-2Л)

 

 

 

 

 

 

 

д

 

 

е

 

 

 

Рис. ІѴ-71. Сигнальный граф,

соответствующий

уравнению для

установившегося состояния

в абсорбционной колонне (а) и этапы

 

его преобразования (б—е).

 

 

 

В результате

эквивалентных

преобразований

сигнального

графа (рис.

ІѴ-71, бд) определяем передаточную

функцию W — у3 (s)/y (s)

для трехта­

рельчатого абсорбера (рис. ІѴ-71,

е).

 

 

 

 

 

Пример 1Ѵ-24.

Найти передаточные функции ХТС однопоточной моноэтанол-

аминовой очистки

синтез-газа от С02

(рис.

ІѴ-72,

а) с

применением правил

эквивалентного преобразования сигнального графа, построенного по структур­ ной блок-схеме системы.

Технологическая схема ХТС (рис. ІѴ-72, а) содержит два основных эле­ мента: абсорбер и регенератор, связанные между собой контуром циркуляции абсорбента, и ряд вспомогательных аппаратов (теплообменники, скруббер-ох­ ладитель и т. д.), которые обеспечивают непрерывное и более экономичное ве­ дение процесса.

В технологической схеме ХТС абсорбер и регенератор через контур цирку­ ляции раствора связаны параметрически по каналам температуры t, расхода

192

X

б

Рис. ІѴ-72. Технологическая схема (а) и структурная блок-схема (б) ХТС одно­ поточной моноэтаноламиновой очистки синтез-газа от двуокиси углерода:

1 — абсорбер; г — сепаратор; 8 — холодильники;

4 — теплообменник,

5 — кипя­

тильник; 6 — регенератор; 7 — скруббер-охладитель;

8 — насосы; ОГ — очищенный газ;

КГ — конвертированный газ; HP — насыщенный раствор; РР — регенерированный раствор; ПГС — парогазовая смесь; П •— пар; К — конденсат; Ф — флегма.

13 Заказ 4 13

Рис. ІѴ-73. Сигнальный граф, соответствующий структурной блоксхеме ХТС моноэтаноламиновой очистки синтез-газа от двуокиси углерода (а), и этапы его преобразования (б — г).

раствора L и содержания двуокиси углерода х в растворе моноэтаноламина (МЭА). На основе связей параметров процесса в технологической схеме и принятых динамических моделей отдельных аппаратов на рис. ІѴ-72, б пред­ ставлена структурная блок-схема ХТС в целом.

Основные обозначения, приведенные на блок-схеме: L — расход раствора МЭА; G — расход неочищенного газа; Q — расход пара в кипятильнике; У — концентрация С02 в неочищенном газе; У0 — концентрация С02 в очищенном газе; X — концентрация С02 в растворе МЭА; t — температура раствора МЭА;

Т — температура пара; [WX, W 2, . . .,

W12]; [W13, Wu , . . .,

W21]; [W{,

W2, . . .,

ІУ4] —передаточные функции

абсорбера, регенератора

и теплооб­

менника

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

 

 

Из анализа технологической схемы и структурной блок-схемы ХТС следует, что: все параметры системы по характеру воздействия можно условно разделить на внутренние (X, t) и внешние (G, У, С, Т, L); в блок-схему входят перекре­ стные обратные связи; все параметры системы взаимосвязаны. Например, изме­ нение любого из внешних воздействий абсорбера даст реакцию на выходной параметр через передаточную функцию колонны и «последействие» через контуры внутренней обратной связи.

На рис. ІѴ-73, а представлен сигнальный граф, соответствующий структур­ ной блок-схеме ХТС. Источниками на исходном графе являются все параметры внешних воздействий (G, У, Q, Г и А), а стоками — переменные внутренних связей процесса и основной выходной параметр У0 (содержание двуокиси угле­

рода в очищенном газе).

 

некоторые преобразования исходного

графа,

На рис. ІѴ-73, б показаны

а на рис. ІѴ-73, в — граф,

преобразованный к виду, удобному для

анализа

взаимосвязей переменных

ХТС.

 

 

По сигнальному графу (рис. ІѴ-73, б) можно определить зависимость

выходной координаты У0 (р) от

внешних возмущающих воздействий

системы

У (р) и Q (р):

+««■)О+

+

«■>-у<')['!

I_____________________________ НМашА___________________________

' 1 — (ez M ? N z + eg) + egez —vw +

(vcoez + vcoM eN z + vcoeg) + vcosgez

где обозначения приняты в соответствии с преобразованиями исходного графа

(см. рис. ІѴ-73).

Т о п о л о г и ч е с к а я ф о р м у л а д л я р е ш е н и я с и г н а л ь н ы х г р а ф о в ХТС

Как уже ранее отмечалось, сигнальные графы можно решать путем составления, а затем решения соответствующих данным гра­ фам систем линейных уравнений ХТС. Одной из привлекательных сторон теории сигнальных графов является возможность решения графов непосредственно, без составления систем уравнений. Эта возможность может быть реализована двумя методами.

Первый метод состоит в том, что исходный сигнальный граф с по­ мощью описанных ранее правил эквивалентных преобразований сводится к конечному графу (рис. ІѴ-74), который состоит из одного источника и одного стока. Однако указанный метод может оказаться

13*

195

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

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

П р я м о й п у т ь Р{ — это путь вдоль ветвей графа от источ­ ника к стоку, при прохождении которого ни один из узлов не встре­ чается более одного раза. П е р е д а ч а прямого пути равна про­ изведению передач ветвей этого пути. Между источником и стоком

сигнального графа может быть несколько прямых путей.

К о н т у р

 

 

( т о ч н е е э л е м е н т а р н ы й к о н т у р )

asi

—о

и п е т л я

сигнального

графа представляют

Хі

собой

замкнутый путь,

вдоль которого ка­

XS

Рис. ІѴ-74.

Конечный

ждый

узел

встречается только

по одному

разу. П е р е д а ч а

к о н т у р а

и л и

сигнальный

граф.

п е т л и сигнального

графа L

равна

про­

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

«касание».

К а с а ю щ и м и с я

называются

контуры (петли),

имеющие хотя бы одну общую вершину.

 

 

 

Графы с одним источником

 

Т о п о л о г и ч е с к а я ф о р м у л а д л я

о п р е д е л е ­

н и я к о э ф ф и ц и е н т а п е р е д а ч и от

и с т о ч н и к а

к с т о к у

с и г н а л ь н о г о

г р а ф а записывается следующим

образом:

П

 

 

 

 

 

 

 

 

2 р і Ср) А/ (р)

 

. . • Рп (Р) Ап (р)

Tis(p)--

)=1

Рі (рі (р)+ Р2(р)А2(р)+

 

А (Р)

А (р)

 

 

 

(IV,39)

 

 

 

 

Здесь р — комплексная частота (оператор Лапласа); Р{ — пере­ дача /-го прямого пути от источника і к стоку s; п — число прямых путей; Д — определитель сигнального графа, выражающийся в об щем виде формулой:

А—1 — L m i + 2 Р т ъ —S ^ т з + • • •

(IV,40)

где 2А?»і — сумма передач всех контуров и петель обратной связи данного сигнального графа; Ьш — произведение передач т-ой воз­ можной комбинации передач из г несоприкасающихся контуров и петель обратной связи (г ^ 2); А/ — минор определителя сигналь­ ного графа, значение которого равно значению определителя А для подграфа, не касающегося /-го прямого пути.

196

Заметим, что предварительные выкладки, связанные с нахожде­ нием значений Р Ау и Д и предшествующие получению формулы (IV, 39), довольно трудоемки. Однако они систематизированы до такой степени, что превращаются в чисто формальные операции, в результате чего в значительной стенени уменьшается возможность субъективных ошибок.

Пример ІѴ-25. Выписать все контуры и пути, начинающиеся в источнике хг

сигнального графа, изображенного на рис. ІѴ-75. Проанализировать касание

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

 

/

Пути, начинающиеся в узле 1:

 

 

7*12—Ъ\ Pi2 z= bd\

P\2 zi— bdf',

Р124= be',

 

P\H3—beg\

P13 = a’,

7*із2= ас;

7*4324= ace;

 

 

P13i= a f

 

 

 

 

Пути, начинающиеся в узле 2:

 

 

 

Р ц —е> Рмз—е8’

7>23 = £^і

Р2 3 4 “

df

 

Пути, начинаюпщеся в узле

3:

 

 

 

Т*3 4 = f j Р 32 = С, Р з 2 4 ~ СС

 

 

Пути, начинаюнщеся

в узле

4:

 

 

 

 

7 * 4 3 = £1

 

Р 4 3 2 = go

 

 

Рис. ІѴ-75. Определение

Контуры графа:

 

 

 

 

 

контуров и путей сиг­

L232=dc;

L ’43 = gf\

L3432—egc;

Z,44= Z0

нального графа.

 

Определяем касание контуров и путей (-*- — касается; \ — не касается):

 

контуры—контуры

пути—контуры

 

7232—* £343

 

 

Р\ 2

—> 72з2

 

L232—* L2432

 

 

7*12 —* 72432

 

7343 —* Lz432

 

 

7*1234—>7 2432

 

7343

Р4 4

 

 

7*1234 "+744

 

72432

* 7,44

 

 

7*1234—>734з

 

7-232

"f 7,44

 

 

7*12

Т 744

 

 

 

 

 

 

7*12

Т 7з4з

Пример ІѴ-26. Определить передачи

Т1 2 = x jx x для сигнальных графов,

изображенных на рис. ІѴ-76, а, б, с применением универсальной топологи­

ческой формулы.

Нахождение передачи

Т1 2 для сигнального графа, представленного на

рис. ІѴ-76, а, поясняют следующие операции.

Определяем передачи путей и контуров сихнального графа:

р 1 = сеЬ,‘

Рз=аЪ\

Li= h;

L3=1g\ L3—de

Отметим касание контуров:

 

 

 

Ls\

L i—>72;

72 >І) 3

197

Соседние файлы в папке книги из ГПНТБ