
книги из ГПНТБ / Кафаров, В. В. Принципы математического моделирования химико-технологических систем (введение в системотехнику химических производств) учеб. пособие
.pdfфундаментальный цикл, а строки размещать так, чтобы сначала шли строки, отвечающие хордам, а затем — ветвям дерева (в порядке возрастания номеров).
В общем случае дипломатическая матрица для любого цикличе ского потокового графа ХТС может быть записана следующим обра зом:
Г Е
[С] = |
(IV,12) |
|
F |
Здесь [Е] — единичная матрица, порядок которой равен \(е —
— р) X (е — р) ]; [F] — матрица порядка [р X (е — р)].
Для графа, показанного на рис. ІѴ-9, г, цикломатическая мат рица будет:
|
Pi |
P2 |
Рз |
|
Я-2 " +1 |
0 |
0 - |
||
Яз |
0 |
+ 1 |
0 |
|
Яв |
0 |
0 |
+ 1 |
|
Яі |
—1 |
0 |
+1 |
|
Яі |
—1 |
+1 |
+1 |
|
Яъ |
0 |
—1 |
+1 |
- |
3. ХАРАКТЕРИСТИКА II ПРИНЦИПЫ ПОСТРОЕНИЯ ТОПОЛОГИЧЕСКИХ МОДЕЛЕЙ (ГРАФОВ) ХТС
Потоковые графы
Каждой ХТС можно поставить в соответствие потоковый граф, ■гомоморфный рассматриваемой системе и являющийся некоторой топологической моделью одного типа обобщенных или физических
потоков данной системы (рис. ІѴ-11). Потоковые графы строят для установившегося технологического режима ХТС.
Потоковый граф
G = G(A) = (A, Т)
с множеством вершин А, образованным совокупностью элементов, источников и' стоков ХТС, и с множеством дуг Т, элементы которого соответствуют одного типа обобщенным или физическим потокам
128
системы, представляет собой некоторое |
семейство сочетаний или |
пар вида |
Ь£А |
Т — (а, Ъ), где а £ А , |
указывающее, какие вершины графа являются смежными. При пост роении потоковых графов ХТС к элементам множества вершин А отнесем и точки стыковки трубопроводов системы.
Выделим три типа потоковых графов ХТС: материальные, тепло вые (энергетические) и параметрические-
Материальные потоковые графы (МПГ). Эти графы подразделим на графы по общему массовому расходу физических потоков и графы но массовому расходу некоторого химического компонента (химиче ского элемента).
В е р ш и н ы м а т е р и а л ь н о г о п о т о к о в о г о г р а фа н о о б щ е м у м а с с о в о м у р а с х о д у ф и з и ч е с к и х п о т о к о в (МПГО) соответствуют элементам ХТС, кото рые трансформируют общие массовые расходы физических потоков, источникам и стокам веществ физических потоков. Дуги этого графа
отвечают обобщенным материальным потокам типа W 1{ |
[см. выра |
жение (11,5)1. |
г р а ф а |
Вершины м а т е р и а л ь н о г о п о т о к о в о г о |
п о м а с с о в о м у р а с х о д у н е к о т о р о г о х и м и ч е с к о г о к о м п о н е н т а соответствуют элементам ХТС, транс формирующим массовые расходы химического компонента, внешним и внутренним источникам, а также стокам этого компонента в сис теме. Дуги данного графа отвечают обобщенным материальным потокам типа W 2i [см. уравнения (11,6) и (11,7)].
Тепловые потоковые графы (ТПГ). Вершины теплового потокового графа соответствуют элементам системы, которые изменяют расходы тепла физических потоков, внешним и внутренним источникам и стокам тепла ХТС. Дуги теплового потокового графа отвечают обобщенным тепловым потокам W3i [см. выражения (11,8) и (11,9)].
Отметим основные характерные особенности материальных пото ковых графов по общему массовому расходу физических потоков и тепловых потоковых графов ХТС:
1. О р и е н т и р о в а н н о с т ь , так как движение потоков веществ и энергии в системе происходит в строго определенном направлении.
2.А с и м м е т р и ч н о с т ь , потому что не все соседние эле менты системы связаны между собой обратными технологическими потоками.
3.С в я з н о с т ь , так как все элементы в системе взаимосвя заны единой цепью потоков веществ или энергии.
Характерные особенности МПГ по массовому расходу некоторого химического компонента в основном совпадают с особенностями МПГО и ТПГ (ориентированность и асимметричность). Отличие со стоит в том, что эти графы могут быть как связными, так и несвяз ными. Несвязность МПГ по компонентам обусловлена тем, что
9 Заказ 413 |
129 |
образование, химическая убыль или целевое выделение данного химического компонента происходят лишь в некоторых элементах ХТС при наличии необходимых технологических условий.
Пример 1V-1. Для ХТС, операторная схема которой нзображепа на рис. ІѴ-12, а, построить топологические модели в виде материального пото кового графа по общему массовому расходу физических потоков, материального потокового графа по массовому расходу каждого из химических компонентов
и теплового потокового графа.
Каждый г-ый физический поток системы Ьі характеризуется набором сле
дующих параметров: |
qi — общий массовый расход физического потока; |
<?</ = |
= Чісіj — массовый расход /-го химического компонента (/ = А, В, Е) |
в t-ом |
|
физическом потоке; |
с// — концентрация; Ц- — температура; Я, — расход тепла |
|
физического потока. |
|
|
Рис. ІѴ-12. Операторная схема (а) и материальный потоковый граф по общим расходам ХТС (б).
В операторе химического превращения (II) протекает химическая реакция А -j- В Е, сопровождаемая выделением тепла Qi- Компонент А подается на вход оператора химического превращения в избытке по сравнению со сте хиометрическим мольным соотношением компонентов А : В.
Степень разделения компонента А в операторе разделения (III) равна величине ф = сьА/сіА. В операторе охлаждения (IV) происходит охлаждение
нагретого физического потока Ьъ за счет физического потока охлаждающей
воды |
Ь7. |
|
|
|
ХТС |
соответствуют следующие параметры: |
||||||
Физическим потокам данной |
||||||||||||
|
|
L l = { q i , С1Л’ С1В’ |
-^l) > |
|
{92, |
C2A’ c2ß' |
^2, ^ 2}» |
|||||
|
„ |
^3— {Я2, сзАі |
С З В ’ |
*3' |
Нз); |
L i = |
(?4, |
С 4 А |
, C i E |
, |
t i , |
Н4}; |
|
|
^ь~{Яьі |
С5А' |
*5> Ньу, Т6 = |
|д 5, с&А, |
tb, |
Я 5|; |
|
||||
|
|
£ 7= {</7. hi |
Н7}; |
/^8 = {37, |
h> |
Hs) |
|
|
|
|||
МПГО рассматриваемой ХТС представлен на рис. |
ІѴ-12,б. |
Так как опера |
||||||||||
торы |
II и |
IV не преобразуют общие массовые расходы |
физических потоков, |
то отвечающих этим операторам вершин в графе не существует. Порядок построе ния МПГО, МПГ по компонентам А, В и Е, а также ТПГ поясняет таблица соответствия элементов ХТС (технологических операторов) вершинам потоковых
графов (см. табл. ІѴ-1). |
внутренний сток s'A, который соот |
|
В МПГ по компоненту А появляется |
||
ветствует количеству компонента <fA ~ 1 |
(cg/1), вступившего в химическое пре |
|
вращение (рис. ІѴ-13, |
а). В МПГ по компоненту В внутренний сток sB свиде |
|
тельствует о том, что |
компонент В полностью вступает в химическую реакцию |
|
д*в = <72В (рис. ІѴ-13, |
б). В МПГ по компоненту Е внутренний источник дЕ |
130
ТАБЛИЦА V I-1
Таблица соотв етств и я элем ентов ХТС и верш ин потоковы х графов
|
Элемент |
ХТС |
|
|
Внеш ний источник |
питания . . . . |
|||
Оператор смеш ения ( I ) ........................ |
|
|||
Оператор |
химического |
превращ ения |
||
( И ) ................................................................... |
|
|
|
|
Оператор |
разделения |
(III) |
. . . |
|
Оператор |
охлаж дения |
(IV ) . . . . |
||
Внеш ний |
сток п р о д у к т а ........................ |
|
|
Вершины потоковых |
графов |
|
|
|
МПГ по |
компонентам |
|
|
МІІГО |
А |
В |
Е |
ТПГ |
|
|
|||
І |
ІА |
і в |
Не т |
іт |
I |
A I |
ВІ |
Нет |
1 |
Нет |
А П |
в п |
Е ІІ |
2 |
III |
А 111 |
Нет |
Нет |
3 |
Нет |
Н ет |
Нет |
Нет |
4 |
s |
SA |
Нет |
SE |
ST |
отвечает образованию нового химического компонента вследствие протекания химической реакции (рис. ІѴ -13, «). В ТПГ (рис. ІѴ-14) внутренний источник і'т соответствует количеству тепла Qi, которое выделяется в результате экзо
термической реакции.
Ч5Л
Рис. ІѴ-13. |
М атериальные потоковые графы по массовым расходам : |
|
|
||||||||||
|
а — компонента А; б — компонента В, |
в — компонента Е. |
|
|
|||||||||
Пример ІѴ -2. |
Д ля Х ТС , |
операторная схема которой представлена на |
|||||||||||
рис. ІѴ -15, |
а, |
построить материальный потоковый граф по расходам |
хим и |
||||||||||
ческих компонентов В, Е и А. |
В операторах химического превращ ения |
II |
и V |
||||||||||
протекает химическая |
реакция |
А -|- В |
Е. |
Ф изическим |
потокам данной |
ХТС |
|||||||
отвечают следую щ ие |
наборы |
параметров: |
|
|
|
|
|
|
|||||
Li = {qi = qlA}\ |
Т2={<72 = <72В}; |
L3 = {q3, q3A, q3B}; |
|
|
|||||||||
^ 4 = { Я 4 , |
Я4 А 1 Ч4 e } ’ |
|
|
, { ? 5 = 9 б л } ! |
£ e = { 9 e > Я6 А ’ |
?6.e } ' |
|
|
|||||
^7 = |
{<?7=97 B '— ?6e } ’ |
^ 8 = |
\Я8 = |
Я3а ~ Я 6а } ’ |
{?9 = |
99b }; |
|
|
|||||
^ 1 0 = {fflO! |
Ai |
^ i o e } ’ |
^ ' l l = |
{ 9 l l |
? i i A = |
9 i o a }> |
^ 1 2 = {912 = 9 i 2e } |
|
9* |
131 |
Р и с / ІѴ -16. М атериальные потоковые графы по массовым расходам :
а — ко м п о н ен та Е ; б — ком п он ен та А .
132
Потоковые графы GB и GЕ по массовым расходам химических компонентов
В и Е (см. рис. ІѴ -15, б и рис. |
ІѴ-16, а) |
являю тся несвязны ми графами. К аж ды й |
|
из этих графов |
включает два |
связны х |
подграфа. Так, например, несвязны й |
граф GB (см. рис. |
ІѴ -15, б) имеет два компонента связности — связны й подграф |
GlB и связный подграф G2B. В подграф GJB входят вершины i ß , B l l ns'B, в подграф»
G2B — вершины |
і'в , ВУ и sB ; подграф |
GlE образован множ еством верш ин iE f |
|||
E l l , EYV и s’E, |
подграф G2E — множ еством верш ин і ’Е, |
ЕУ и s"E . |
|
||
Поток химического компонента А |
образует единую |
цепь м еж ду |
всеми эле |
||
ментами |
Х ТС , |
поэтом у МПГ по м ассовому р асходу компонента А |
является |
||
связным |
графом (рис. ІѴ -16, б). |
|
|
|
В общем случае материальный или тепловой потоковый граф ХТС (см. рис. ІѴ-11) содержит т вершин-источников, п вершин-стоков, к промежуточных вершин, соответствующих элементам системы, и е дуг, которые отвечают одного типа обобщенным материальным или тепловым потокам системы.
Каждой дуге материального или теплового потокового графа ХТС можно сопоставить неотрицательное число W (е), которое равно значению одного типа обобщенного материального или тепло
вого потока системы. Величину |
W (е) назовем п о т о к о м п о |
д у г е е данного графа. Для |
каждой промежуточной вершины |
материального или теплового потокового графа ХТС на основе зако нов сохранения массы и энергии можно записать « у р а в н е н и е
в е р ш и н » |
для потоков по |
дугам графа: |
|
|
||
|
е |
|
|
|
|
|
|
S |
^ |
* = |
0 |
(IV , 13) |
|
|
^ w iK- |
2 |
WKl=0 |
|
|
|
|
i=1 |
|
/-1 |
|
|
|
где WjK — потоки по /-о й дуге графа, |
инцидентной гг-ой вершине; |
Wік — |
поток |
|||
по г-ои дуге, входящ ей в «--ую верш ину графа; WKJ- — поток по /-о й д у ге, |
выхо |
|||||
дящ ей из к-ой |
вершины графа. |
|
|
|
|
|
Совокупность «уравнений вершин» для потоков по дугам графа (IV,13), которые составлены для всех промежуточных вершин, обра
зует с и с т е м у |
н е з а в и с и м ы х |
у р а в н е н и й в е р ш и н |
||
д а н н о г о г р а ф а : |
|
|
|
|
|
[A*]x[W *] = 0 |
(IV,14) |
||
Здесь [W*] — матрица-столбец |
потоков по |
всем дугам графа, |
||
порядок которой равен (е X 1); [А*] — матрица |
инциденций данного |
|||
графа, порядок |
которой равен |
(* |
X е). |
|
Нетрудно видеть, что матричное уравнение вершин материаль ного или теплового потокового графа (IV, 14) эквивалентно матрич ному уравнению балансов одного типа обобщенных потоков данной ХТС в целом (11,12).
Эквивалентность матричных уравнений (IV, 14) и (11,12) доказы вает, что материальные и тепловые потоковые графы являются топологическими моделями, гомоморфными исследуемой ХТС.
На основе законов сохранения массы и энергии для всей ХТС в целом можно записать общее уравнение вершин-источников и вер шин-стоков материального или теплового потокового графа системы:
2 Wm (Ь) — 2 Wn (с) = 0 |
(IV,15) |
тп
где Wm (Ъ) — поток по è-ой дуге, инцидентной m-ой вершине-источ нику; Wn (с) — поток по с-ой дуге, инцидентной п-ой вершине-стоку; т — число вершин-источников; п — число вершин-стоков данного графа.
Циклическим потоковым графом называют связный граф, полу ченный из материального или теплового потокового графа ХТС путем объединения всех вершин-источников и вершин-стоков в одну общую (нулевую) вершину, для которой справедливо уравнение (IV,15). Таким образом, для каждой вершины указанного графа ХТС можно составить уравнение вершин (IV,13). Циклический потоковый граф G* = (V, В), который соответствует потоковому графу G — (А , Т), имеющему m вершин-источников, п вершин-стоков, к промежуточ
ных |
вершин |
и е дуг, содержит число вершин ѵ = |
к -f- 1 и число |
дуг |
Ъ = е. |
|
|
Циклический потоковый граф ХТС, потоковый граф которой изо |
|||
бражен на |
рис. ІѴ-11, представлен на рис. ІѴ-17. |
Построение и |
исследование топологических особенностей материальных и тепло вых графов позволяет формализовать процесс составления и полу чения оптимальных алгоритмов решения систем уравнений балансов ХТС.
Параметрические потоковые графы (ППГ). Такой граф является топологической моделью, отображающей преобразование элемен тами системы параметров физических потоков ХТС. Вершины ППГ отвечают элементам, представляющим собой технологические опера торы, которые качественно и (или) количественно преобразуют пара метры физических потоков, а также источникам и стокам физиче ских потоков ХТС.
Дуги графа соответствуют физическим потокам системы. Каждой дуге ППГ сопоставим некоторое неотрицательное число яг — п а р а - м е т р и ч н о с т ь этой дуги. Параметричность дуги графа равна параметричности соответствующего физического потока системы. В общем случае все дуги параметрического потокового графа слож ной ХТС разнопараметричны.
Параметрический потоковый граф можно получить совмещением одноименных вершин, т. е. вершин, отвечающих одному и тому же элементу, источнику или стоку ХТС, материальных и тепловых потоковых графов, построенных для данной исследуемой системы. Параметрические потоковые графы ХТС обладают следующими характерными особенностями: они являются конечными, ориентиро ванными, асимметричными, связными как планарными, так и непла нарными графами.
134
Изучение топологических особенностей ППГ позволяет разраба тывать оптимальные алгоритмы анализа функционирования ХТС произвольной технологической структуры.
Пример ІѴ-3. Для ХТС, рассмотренной в примере ІѴ-1, построить пара
метрический потоковый граф.
ППГ исследуемой ХТС (рис. ІУ-18) имеет разнопараметрические дуги, а именно:
р і = р 2 = р з = р 4 = 5; р 5 = р 6 = 4 ; р 7 = р 8 = 3
Все элементы преобразуют параметры физических потоков, поэтому каждому элементу соответствует определенная вершина ППГ.
Чз
Рис. ІѴ-17. Циклический пото- |
Рис. ІѴ-18. Параметрический потоковый |
ковый граф ХТС. |
граф ХТС. |
Построение потоковых графов, для которых известны технологи ческая топология и цель функционирования системы, необходимо осуществлять по следующим этапам:
1.Изучить физико-химическую сущность технологических про цессов ХТС и построить ее операторную схему.
2.Составить таблицы источников и стоков ХТС.
3.Составить матрицы [К] покомпонентного состава физических потоков системы. Общий элемент такой матрицы
1, если /-ый химический компонент входит в состав {г-го физического потока
О—в противном случае
причем г = 1, п — число физических потоков ХТС; = 1, т — чис ло химических компонентов системы.
4.Составить таблицы одного типа материальных и тепловых обобщенных потоков системы.
5.Составить таблицы соответствия элементов, источников и стоков системы вершинам ее потоковых графов.
6.Составить таблицы соответствия одного типа обобщенных
потоков и физических потоков ХТС дугам ее |
потоковых графов. |
7. Графически изобразить потоковые графы |
ХТС. |
135
Структурные графы ХТС
Решение задачи исследования сложных ХТС можно значительно облегчить, используя для получения элементов матриц преобразо вания методику, основанную на построении структурных графов. Топологическая модель ХТС может быть изображена в виде некото рого структурного графа, гомоморфного данной системе. Этот метод иконографического моделирования системы с сосредоточенными параметрами базируется на представлении ХТС как упорядоченной совокупности определенных системных компонентов, коэффициенты
передачи которых |
известны. |
|
В каждой ХТС можно выделить ряд с и с т е м н ы х к о м п о |
||
н е н т о в , или |
п р о с т ы х и д е а л ь н ы х с и с т е м н ы х |
|
э л е м е н т о в , |
характеризуемых определенными свойствами: |
|
1. |
Элементы, создающие потенциальную или кинетическую энер |
|
гию, |
— к о м п о н е н т ы - и с т о ч н и к и . Примерами таких эле |
ментов являются «идеальные» компрессоры и насосы, не обладающие
потерей напора |
или |
утечкой |
потока; |
различные нагреватели |
|
и т. д. |
|
|
|
|
|
2. Элементы, |
рассеивающие |
энергию системы, — р е з и с т и в |
|||
н ы е к о м п о н е н т ы |
( с о п р о т и в л е н и я ) . |
Так, например, |
|||
гидравлическим |
сопротивлением |
являются |
участки |
трубопроводов |
и клапаны; кроме того, гидравлическое сопротивление характери зует потерю напора при истечении жидкости из аппарата.
3. Элементы, обладающие способностью накапливать вещество
или |
энергию системы, а также отражающие свойства упругости |
веществ, — е м к о с т н ы е к о м п о н е н т ы . |
|
4. |
Элементы, характеризующие инерционный эффект массы в по |
токе |
вещества, — и н д у к т и в н ы е к о м п о н е н т ы . |
Для исследования гидродинамических и тепловых процессов любую ХТС можно представить как упорядоченную совокупность рассмотренных системных компонентов, объединенных между собой с помощью т о ч е к с в я з и ( п о л ю с о в с в я з и).
Тип и число системных компонентов, а также способ их соедине ния обусловливают свойства системы и определяют значения эле ментов матриц преобразования гидродинамических и тепловых про цессов. Представим характеристики каждого системного компонента
некоторой ХТС, |
если известны ее параметры, в виде п о л ю с н ы х |
у р а в н е н и й , |
включающих два типа скалярных величин, кото |
рые связаны с двумя классами измерений на «полюсах» изолирован
ного |
системного компонента |
(рис. ІѴ-19, а). |
К |
п е р в о м у т и п у |
с к а л я р н ы х в е л и ч и н отнесем |
сигналы, количественные значения которых могут быть измерены, по крайней мере теоретически, с помощью приборов, включаемых параллельно на «полюсы» компонентов (прибор Ш на рис. ІѴ-19, а) без разрезания системы в точках соединения системных компонентов
(например, манометра или термопары). Такие |
сигналы назовем |
« п а р а л л е л ь н ы м и » п е р е м е н н ы м и , |
а соответствующий |
136
этим переменным класс измерений — « п а р а л л е л ь н ы м и » |
и з |
м е р е н и я м и . |
отне |
К о в т о р о м у т и п у с к а л я р н ы х в е л и ч и н |
сем сигналы, количественные значения которых измеряют последо вательным включением системного компонента и измерительного прибора (прибор П2 на рис. ІѴ-19, а) при разрезании системы в точ
ках соединения компонентов (например, |
расходомера |
или |
калори |
|
В- |
метра). |
Этот тип сигналов назовем |
||
« п о с л е д о в а т е л ь н ы м и» |
||||
п е р е м е н н ы м и , |
а |
отвеча |
||
|
ющий |
им класс |
измерений — |
|
Компонент - № > л |
« п о с л е д о в а т е л ь н ы м и » |
|||
и з м е р е н и я м и . |
|
|
||
Источник |
x!t) |
|
|
x(tj,y(t) |
|
|
|
||
энергии |
|
|
|
О >*і о |
Рис. ІѴ-19. Определение коэффициентов передачи двухполюсного системного компонента (а), линейный граф полюсных измерений для компонента (6) и по люсный граф компонента (в).
Каждому д в у х п о л ю с н о м у к о м п о н е н т у |
соответ |
ствуют полюсное уравнение и п о л ю с н ы й г р а ф , |
который |
представляет собой направленную ветвь, отражающую ориентацию приборов для последовательных и параллельных измерений по отно шению к «полюсам» компонента (рис. ІѴ-19, б и ІѴ-19, в).
В общем случае полюсное уравнение компонента системы с сосре
доточенными параметрами имеет следующий |
вид: |
= |
(IV,16) |
где s — комплексная переменная преобразования Лапласа; К (s) — коэффи циент передачи (комплексная проводимость), характеризующий свойства ком понента; У (s), К (s) — изображения по Лапласу для последовательной у (<) и параллельной х (t) переменных.
Полюсные уравнения компонентов ХТС для исследования гидро
динамических процессов, связанные с |
измерением последовательной |
||
переменной — потока вещества Q (t) |
и параллельной перемен |
||
ной — давления Р (t), |
будут: |
|
|
для |
гидродинамического сопротивления |
||
|
Р (s) = RQ (s) или Q (s) — GP (s) |
||
для |
гидравлической |
емкости |
|
|
|
Q (s) = csP (s) |
|
для |
гидравлической |
индуктивности |
|
|
|
Р (s) = LsQ (s) |
137