книги / Эксергетические расчеты технических систем
..pdfнечного числа уравнений
- Ц - = 0 (i = 1,2, |
k; / = 0, 1, . . . . t). |
(9.78) |
Решение системы уравнений (9.78) позволяет вычислить значения экстремалей в заданном числе точек и на их основании построить кусочно-линейные функции, аппроксимирующие экстремали функционала (9.72). При этом считается, что в пределах каждого промежутка времени экстремали остаются постоянными, а их изменения происходят скачкообразно на границах участков. Чем больше вы бранных участков, тем меньше погрешность полученных результатов, и поэтому при t оо задачу можно решить точно.
С учетом изложенного выражение функционала, определяющего годовые приведенные затраты для рассматриваемой одноцелевой холодильной уста новки (рис. 9.8), представляется в виде
П = [Zu (е?', |
©к11’. Д О |
+ |
Z12 (е2, ек, ATJ + |
Z13 (e f\ |
0<13\ AT™) + |
||||||
+ |
Z21 (e3, 0 И> AT5) + |
Z22(e?>, AT?') + Z31 (eQoxji, |
0 B, ДГВ) + |
|
|||||||
+ 2 |
32( ^ |
(32k |
|
+ ц ЭЛ V |
(Еп (еУ, ©!/’, АТ1») Ат] + |
|
|||||
AT?')]л , |
т г |
|
|||||||||
|
|
|
|
|
|
|
/ ± 1 |
|
|
|
|
|
+ |
£ |
[Е13(еУ, ©if, АТш') Дт] + |
£ [Егг( ^ , А?/') Дт] + |
|
||||||
|
|
/=1 |
|
|
|
|
|
/=1 |
|
|
|
+ |
S |
[^ з2 ( С |
л . AT?)Ax]\ + |
nwi |
|
©к71, AT'w) Д т] + |
|
||||
|
/=1 |
|
|
|
|
) |
/—1 |
|
|
|
|
|
|
|
+ |
S |
|
|
©!/’, ДП71) - |
е^1 Дт} + |
|
||
|
|
|
|
/=1 |
|
|
|
|
|
|
|
+ |
Б |
|
[Е3 (4 0ХЛ, ©в\ |
АПЛ) - |
е?) Дт} |
( / = 1 . 2 , . . . , о, |
(9.79) |
||||
|
/=1 |
|
|
|
|
|
|
|
|
|
|
где — функциональный множитель Лагранжа; |
надстрочный индекс / соот |
ветствует текущему значению номера участка, а индексы 11, 13, 22, 32 — режи му максимальной производительности соответствующего (по номеру) элемента оборудования холодильной установки.
Используя выражение функционала (9.79) и соотношения, описываемые зависимостями (9.76) и (9.78), можно получить следующую систему уравнений:
|
|
+ |
. dZ\\l) |
dZ\f |
i \ - щ п [цэл (£</‘ + |
^ |
Н - ^ Г Г + |
= 0; |
|
|
сэе(11> |
d0il3) |
||
/ = 1 |
|
|
|
к |
/ |
- (гТ е У + Ztf) |
— 0; |
|
|
V |
|
|||
/=1 |
дв? |
|
|
|
/=21 |
dQ(j' |
+ Z'i) |
= 0; |
|
|
|
|
|
|
|
|
|
z\$] + |
|
/5 7 (H) |
|
л 7(13) |
|
||
[цЭЛ(£$? + |
£$) + ц.И§ + |
^ |
т |
г |
+ 1 ^ |
| г = 0; |
(9.80) |
|||||
d&lW |
|
|
|
|
|
|
|
|
|
|
|
|
, |
|
|
|
|
|
|
л 7(22) |
|
|
|
||
^ 0 U » + r f W + « > + * £ r - f c |
|
|
||||||||||
|
|
|
|
|
^7(32) |
|
|
|
|
|
|
|
•j- (Цэл-^32 |
+ 4 |
Л^ |
+ 4 ? ) |
+ |
0 Д 7 ,« 2) |
= |
® |
( |
/ = 1 |
. 2 |
...............О» |
|
d&T(J> |
|
|
|
|
|
|
|
|
|
|
|
|
г д е |
|
|
|
|
|
|
|
|
|
|
|
|
г? = ^ т г [Цэл (М? + |
Е% + |
nwv?2 + № |
|
|
/57(11) |
|
/57ОЗ) |
|
||||
+ |
|
- ^ 5 - |
+ |
; |
(9.81) |
|||||||
дер |
|
|
|
|
|
|
|
|
|
|
|
|
ГФ |
= |
л |
|
|
+ |
4?) + |
|
<97^2’ |
|
|
||
ГЗ |
(Цэл^М + |
|
|
22 |
‘ |
|
|
|||||
|
|
|
|
|
|
|
|
|
* Г |
|
|
Система уравнений (9.80) с учетом выражений (9.81) представляет собой запи санное в общем виде решение задачи динамической оптимизации одноступенча той одноцелевой холодильной установки, позволяющее получить набор зна чений экстремалей, обеспечивающих с учетом сделанных допущений минимум функционала годовых приведенных затрат (9.79).
Решение системы уравнений (9.80) с учетм (9.81) дает наборы оптими
зирующих |
переменных |
©Д |
АTw\ ©Д ДГД ©Д |
AT{J' ( / = 1 , 2, |
t), |
объемных |
расходов |
сред |
у}#, иД 0Д V{b и |
соответствующих |
этим |
объемным расходам мощностей электродвигателей насосов охлаждающей сре ды е(/з, промежуточного хладоносителя djl, электродвигателя вентилятора £32
и электродвигателя компрессора |
£п, обеспечивающих |
минимальное значение |
||
функционала (9.79). Набор значений i//2\ vKJ\ |
И/’ (j |
= 1, 2, |
..., t) отражает |
|
закон оптимального управления |
рассматриваемой |
холодильной |
установкой с |
учетом использованных допущений. В то же время значения V\$\ vf2\ г432), ИД служат для подбора соответствующего оборудования (насоса охлаждающей воды, насоса хладоносителя, вентилятора воздухоохладителя, компрессора) по их максимальной в течение рассматриваемого периода производительности. Поверхности теплопередающих аппаратов (конденсатора FK, испарителя Fuy воздухоохладителя FBO) согласно принципу их полного использования (9.75) могут определяться по оптимальному значению соответствующего среднелога
рифмического перепада температур (©[/*, ©Д ©Д для любого рассматриваемого промежутка времени при условии, что тепловой поток и коэффициент тепло передачи определены для того же промежутка времени.
Аппроксимация полученного набора значений для каждой функции дает выражение для экстремалей (минималей) функционала, т. е. закон изменения оптимизирующих переменных в течение года, обеспечивающий минимальную себестоимость производимого холода.
Переход к динамической оптимизации приводит к резкому возрастанию числа уравнений в задаче. Так, если при статической оптимизации для рас сматриваемой установки необходимо решить систему из шести уравнений (9.35), то в динамической оптимизации, даже при рассмотрении только среднемесяч
ных значений температур (t = 12), количество этих уравнений увеличивается до тридцати девяти.
Решение системы такого количества уравнений, которое для более сложных схем установок еще более возрастет, даже при наличии современной вычисли тельной техники затруднительно. Однако в Ленинградском технологическом институте холодильной промышленности на базе теории динамического про граммирования разработана методика решения, позволяющая свести решение /^-мерной задачи к решению tk одномерных задач вида у (х) = 0 [126].
Сущность метода сводится к решению задачи динамической оптимизации в две стадии. На первой стадии при некоторых значениях параметров системы Ft (9.75) за счет выбора оптимального управления ДТц находится локальный минимум; на второй стадии определяется оптимальное значение параметров Fi, обусловливающих оптимальные значения 0 t/, и осуществляется выбор оп тимального управления АТ1*/, обеспечивающего минимум рассматриваемого функционала (9.79).
Учитывая, что применение математического аппарата, в частности решение задачи Лагранжа, справедливо для непрерывных функций, решение оптими зационной задачи и выбор оборудования осуществляются в три этапа. На пер вом решается задача динамической оптимизации холодильной установки, ис ходя из непрерывного ряда используемого основного оборудования (тепло передающих аппаратов, компрессоров, насосов, вентиляторов). В результате получается строгое решение оптимизационной задачи, являющейся предельным идеализированным вариантом.
На втором этапе на основе методов комбинаторного анализа подбирается реальное, выпускаемое промышленностью оборудование. Подбор осуществля ется на основе сопоставления характеристик реального оборудования с харак теристиками оборудования, полученными на первом, идеализированном, эта пе решения. Возможны три варианта подбора оборудования: полной заводской готовности (конденсаторы, испарители, воздухоохладители), комплектуемое из стандартных модулей (охлаждающие батареи) и с регулируемой производитель ностью (компрессоры, насосы, вентиляторы). Критерии для выбора наибо/iee рационального варианта, которыми могут служить, например, максимальная степень приближения по производительности к идеализированному решению первого этапа, минимальная масса оборудования или занимаемая площадь, задаются ЭВМ либо варианты выбираются проектировщиком по комплексу показателей в процессе диалога человек — машина.
На третьем этапе, используя основную программу динамической оптими зации, для выбранного оборудования решается первая стадия задачи — опре деляется закон оптимального управления подобранным оборудованием. Сопо ставление значений функционала, полеченных на третьем и первом этапах, по казывает степень удаления от оптимального решения, обусловленную реальной
номенклатурой основного оборудования.
Рассмотренная модель одноцелевой одноступенчатой холодильной установ ки принимается в качестве базовой. Моделирование более сложных систем осу ществляется путем некоторого дополнения базовой модели.
[Цэл (Еп + Е[>] + |
Е\<}) + |
и У £ + |
+ 2 Й 1 + |
dzft1’ |
d Z \f |
dZ\14' |
|
++ дЬТ™ + адтй4» ~~ ’
|
|
[Дэл (£# + |
Е(й + |
Ей) + п У й + У |
+ z$ ) + |
|||||
|
|
|
dZ\\l) |
dZfj' |
|
dZ\14' |
_ 0 . |
|
||
|
|
|
д Д Г ^ |
дД7<п13) |
^ Э Д Г ^ |
|
|
|
||
|
|
_ |
|
|
|
|
|
dZK ’ |
|
|
|
|
- £ ц г |
(Дзл i t f |
+ ^ |
+ |
ZSP) |
+ - |
^ |
= |
0; |
|
|
|
|
|
Q7(32) |
|
|
|
|
|
|
d ^ jW (Цэл-^32 + |
Гз^Ез' + |
Z3i’) + д^ |
(з2) |
— ® |
(/ = |
1. |
2, . . . , / ) , (9.83) |
||
г д е |
Bjj |
|
|
|
|
|
|
|
|
|
|
Г(2) = |
[Дэл (■£'п) + |
£ $ + |
£ Й ) |
+ Ц .У & |
+ Zi2 |
+ |
Z15] + |
||
|
|
|
dZ*,1.11 |
dZ(,33) |
|
dZ^4’ |
|
|
||
|
|
|
# - + ^ l r + - ^ r - |
|
|
|||||
|
r f = |
д |
|
|
|
Л7*22*dz!21 |
|
|
|
|
|
( u „ £ g + / f |
+ Z ^ ) + |
^ - J l p |
( / = |
1 , 2 , . . . , / ) . |
|||||
|
|
|
|
|
|
Ц |
|
|
|
Выражения (9.80) и (9.81), относящиеся к зонам II и / / / модели, не пре терпевают никаких изменений по сравнению с предыдущим вариантом одноце левой одноступенчатой установки, что подтверждает универсальность полу чаемых решений.
Решением системы уравнений (9.83) получаются наборы значений, от ражающие оптимальные законы изменения расхода сред (рабочего тела в ступе ни высокого и низкого давления, охлаждающей воды, промежуточного хладоносителя, воздуха), оптимальные значения площади теплопередающей поверх ности аппаратов (конденсатора, змеевика промежуточного сосуда, испарителя, воздухоохладителя) и максимальные значения производительностей компрес соров, насосов и вентиляторов, обеспечивающие минимум переменной части приведенных затрат для рассматриваемой одноцелевой двухступенчатой холо дильной установки.
Моделирование и оптимизация холодильных установок с оборотной системой водоснабжения
В условиях эксплуатации для отвода теплоты конден сации от холодильных установок широко применяются системы оборотного водоснабжения, в которых вода направляется в атмосферный охладитель (гра дирню), после чего возвращается в конденсатор холодильной установки. По-
Для установки, схема которой изображена на рис. 9.12, а модель — на рис. 9.13, функционал, определяющий переменную часть приведенных затрат,
имеет вид
П = [Zu (eT, А Т?\ Т ? ) + Z12(e2, 0К, АТт Тщ) +
+ Z13{e?\ 0<13\ ДГ"3), T™) + Z21(e3, 0„, ATS) + Z22( e f A 7 f 2’) +
+ Z ai(eQ<JXJI, ©в, АГВ) + Z32( 4 1 . АП32))]тг+
+ |
Цэл(Б [Д и (4л, 0 к \ АГш’, TSl) А т ] + |
S [^ 1з(е2л, 0 ^ , АГ<?, 7 ^ ) Дт] + |
|
|
+ S |
АТ?) Ат] + £ |
[Е32(е?охл, АТ?) Ат]) + |
|
/=1 |
/=1 |
) |
+ |
£ { М ? , |
0 к \ AT?, g?, T?J][V12(e?, в ? , AT?, g?, Т?)]Ах) + |
|
|
/=1 |
|
|
+£ {г? [Е2(е?, 9?, АТ?) — е?] Ат} +
/=1
+ £ [г? [Е3 (е?охл, в ? , АТ?) - е?] Ат} (/= 1, 2, |
t). (9.84) |
Выражение (9.84) написано с учетом определения тарифа на оборотную воду по выражению
/ _ “с.вУс.в + W l7 + г!66’ + г<™
Hw --
”l2
что позволяет использовать одну и ту же унифицированную оптимизационную программу как при прямоточном, так и при оборотном водоснабжении.
На основе функционала (9.84) составляется система уравнений:
£ ( - ^ [ « . и * + m |
+ ш р v e + z g j ) + , ( - Щ |
. + i j g . ) - « |
||
/=1 ' К |
/ |
\ |
к |
к J |
2 |
д0(Л (r?E? + Z?) |
= |
0; |
|
/=1 |
г |
|
|
|
/ |
0; |
|
|
|
2 |
= |
|
|
|
|
|
|
|
|
/=1 |
|
|
|
|
|
■(r?E? + Z?) |
= |
0; |
|
/=1 L д в ? |
ns |
^7(13) |
||
|
|
|||
|
|
dz?> |
|
|
гг P u . ( £ i f + £ 8 ) + Ц « ? + 7 1 Й + — З п г + — i s - _ 0 ; |
||||
дАТ<? |
|
дАТ?> |
<?д г <‘ 3 ' |
|
д г [им(£!f + £!Й + ц“Т8 + z|{j + -Ц£- + |
аг"3’ |
|||
Wt |
|
О*ш, |
||
|
10, |