книги из ГПНТБ / Применение вычислительной техники на металлургическом заводе
..pdf
|
|
|
Л І 2 |
Л І 3 |
111 2 |
Л Ш 3 |
I i |
III |
Л'11 |
175,5 |
1,12 |
0 |
0 |
0 |
1 |
0 |
|
л 'і 3 |
24,5 |
—0,102 |
I |
0 |
0 |
0 |
0 |
|
•Ѵ1ІІ |
, |
309 |
0 |
0 |
1,12 |
0 |
0 |
1 |
л ' ш |
з |
593 |
0 |
0 |
—0,101 |
1 |
0 |
0 |
L |
|
15 790 |
25,6 |
0 |
11,2 |
0 |
0 |
0 |
Поскольку |
в |
индексной |
строке |
нет отрицательных |
элементов, |
||||
получили |
L m n x |
= 15790 при |
х п = |
1/5,5; . v I 2 = 0 ; .v1 3 =24,5; |
. ѵ І І П =309; |
||||
л '1112= 0; |
х Ш з |
=593. Это решение совпадает с результатом, |
получен |
||||||
ным в примере |
1. |
|
|
|
|
|
|
|
|
Очевидно, |
что |
формализованный |
метод, |
рассмотренный в |
при |
||||
мере 2, |
без особых |
затруднении |
может быть |
запрограммирован |
для |
||||
решения задачи на ЭЦВМ. |
|
|
|
|
|
|
|||
6. Ц И Ф Р О В О Е Д И Н А М И Ч Е С К О Е М О Д Е Л И Р О В А Н И Е |
|
|
|||||||
В О З Д У Х О Н А Г Р Е В А Т Е Л Я Д О М Е Н Н О Й |
ПЕЧИ |
|
|
|
На цифровой динамической модели воздухонагрева теля отрабатывают тепловой режим агрегата, настройки автоматических регуляторов температуры купола п го рячего дутья, оценивают эффективность режимных п конструктивных усовершенствований, намечают пути оп тимизации тепловой работы воздухонагревателя. При моделировании воздухонагревателя приняты следующие исходные предположения:
1. Циклы работы всех воздухонагревателей одинако вы—одни агрегат блока воздухонагревателей стоит на дутье, остальные—на нагреве.
2. Расход горячего дутья и его температура постоян
ны.
3.Тепловая мощность воздухонагревателя постоянна.
4.Для регулирования температуры купола воздухо нагревателя и горячего дутья использованы автоматиче ские изодромпые (ПИ) регуляторы.
5. Потери тепла в окружающую среду постоянны при любых режимах работы и могут быть исключены из рассмотрения.
6. Теплоемкость массы газа в каналах насадки пре небрежимо мала по сравнению с теплоемкостью массьі насадки.
102
7.Теплопроводность насадки в продольном направ лении не учитывается.
8.Теплопроводность насадки в поперечном направле нии учитывается путем коррекции коэффициента тепло
отдачи между газами и стенками каналов насадки.
AI |
ê± ^ |
УZ.L±Z
|
|
J |
|
J-l |
|
|
|
(,0) |
|
|
|
пи |
|
Рис. 21. |
Структурная |
схема |
|
|
|
||||
|
|
воздухонагревателя в |
период |
||
|
|
нагрева |
|
|
|
9. Теплоемкости газов и материалов насадки линейно |
|||||
зависят от температуры: |
|
|
|
|
|
с = с0 |
+ |
с' t |
докІ(кг-град). |
(26) |
|
10. Коэффициент |
теплоотдачи линейно зависит от |
||||
расхода газов М: |
|
|
|
|
|
а = а0 |
- j - |
a'M |
вт |
(м--град). |
(27) |
На рис. 21 изображена структурная схема исследуе мой системы в период нагрева. Как видно из рисунка, распределенная насадка воздухонагревателя разбита на ряд участков высотой Al каждый. Каждый такой доста точно малый участок можно считать звеном с сосредо точенными параметрами.
На . рис. 21 приведены используемые в дальнейшем обозначения: ПИ— изодромный регулятор температуры
купола |
воздухонагревателя, |
изменяющий |
коэффициент |
|
избытка |
воздуха, идущего |
на |
сжигание газа; хвх—вход |
|
ная величина регулятора |
(разность между |
действитель |
ной и заданной температурами купола или верхних ря
дов |
насадки); хвых—выходная |
величина |
регулятора (рас |
||
ход |
воздуха на |
горение); |
Al—интервал |
пространственно |
|
го квантования |
насадки; |
/—номер элемента пространст |
|||
венного квантования; Мв |
— расход дополнительного воз |
духа для снижения температуры продуктов сгорания; Мг , /г —постоянные расход и температура продуктов сгора ния отопительного газа; Мл, t-A—расход и температура
103
продуктов сгорания, поступающих в насадку; /І ( Л —сред
няя температура насадки; / , і |
л . і — регулируемый |
параметр; |
|||
^кл.з — заданная температура |
купола |
(верхних |
рядов на |
||
садки) воздухонагревателя; ^у х —температура |
уходящих |
||||
из воздухонагревателя |
продуктов сгорания; |
г в = 0 ° С — |
|||
температура атмосферного воздуха. |
|
|
|||
Для расчетов на ЭЦВМ примем интервал |
временно |
||||
го квантования Ат. |
|
|
|
|
|
Тогда |
|
|
|
|
|
Мл(к- |
1) |
- Мг |
+ Мл(к- |
1); |
(28) |
*д(й—1,0) = |
^ — . |
(29) |
Здесь и далее первый индекс в скобках означает но мер интервала временного квантования, второй — номер элемента пространственного квантования.
Передаточная функция ПИ-регулятора имеет вид
|
|
w ( р) = |
i ï ï i i i = k p 11 + |
-J - |
) , |
(30) |
|||
|
|
|
|
*вх (P) |
\ |
T„p |
I |
|
|
где |
p — параметр преобразования |
Лапласа; |
|
||||||
|
kp— |
коэффициент передачи; |
|
|
|
||||
Т„—время |
изодрома. |
|
|
|
|
||||
Из выражения |
(30) |
следует |
|
|
|
|
|||
|
|
* в ы х = |
|
|
= kp [tKn |
(1) — /к л .з1 |
+ |
|
|
|
|
|
|
т |
|
|
|
|
|
|
|
+ |
^-\[tKAV-tKnJdT. |
|
|
|
(31) |
||
|
|
|
|
о |
|
|
|
|
|
Если от дифференциалов перейти к |
конечным |
раз |
|||||||
ностям, |
то |
для выходной величины регулятора получим |
|||||||
|
|
|
|
|
|
|
|
1 |
|
Mn{k) |
= k„ stKn(k- |
|
і , і ) - / к |
л . 3 +-£-У |
[t*Ak- |
|
|||
|
|
|
|
|
|
|
А — 1=0 |
|
|
|
|
|
|
- |
1,1) —/кя.э]]. |
|
|
(32) |
|
Рассмотрим тепловой баланс /-того элемента насад |
|||||||||
ки. Тепловой баланс |
газа: |
|
|
|
|
104
/ д ( / - 1 , |
|
1)Мд(А-1)Сд - / „ ( / . Л |
~ 1 ) Л 1 д ( е - 1 ) С д |
||||||||||
а д |
FA/ |
|
М / - І , * - 1 ) |
+ /д ( / , * - ! ) |
- і к л |
|
U,k-i) , |
(33) |
|||||
|
|
|
|
||||||||||
где ад —коэффициент |
теплоотдачи: а д = а 0 д - г - а д M ; |
|
|||||||||||
|
сд —теплоемкость |
продуктов |
сгорания: |
|
|
||||||||
|
|
|
|
|
с д = |
с до ~Ь с д ^ді |
|
|
|
|
|||
|
F—поверхность |
нагрева воздухонагревателя, |
м2£.\ |
||||||||||
|
I — полная |
высота |
насадки, |
м. |
|
|
|
|
|||||
Тепловой |
баланс |
насадки |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
а д |
FДт |
X |
|
|
|
|
|
|
|
|
|
|
|
/пс,( |
|
|
|
|
|
X |
< д ( * - 1 , / - 1 ) |
+ |
<д(*-1,/) |
- , к л ( |
^ 1 , / ) 1 |
(34) |
||||||
где |
скп |
= |
с к л 0 + |
с'кд |
і к п — |
теплоемкость |
материала |
на |
|||||
|
|
|
садки; |
|
|
|
|
|
|
|
|
|
|
|
m— масса насадки, кг. |
|
|
|
|
|
|||||||
Обозначим: |
|
|
|
|
|
|
|
|
|
|
|
||
|
Рд = |
FM |
|
|
1 |
|
«до |
|
|
|
|
||
|
|
Л І д ( А - І ) |
С д о + с ' |
^ ( f t - ! , / - ! ) ' |
|
||||||||
|
|
|
|
|
|||||||||
|
|
|
Уя |
f-Дт |
|
а д о + а д М д (f e |
— ] ) |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Тогда |
математическая |
модель |
воздухонагревателя |
||||||||||
в период нагрева может характеризоваться |
выражения |
||||||||||||
ми: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
[ ( 2 - р ^ д ( А - 1 , / - 1 ) + |
|
||||||
|
|
|
|
2 + |
Рд |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
+ |
2Рд^л (А— I,/")]; |
|
|
(35) |
|||||
^кл (А,/) = |
tM |
(k — |
1,/) |
+ |
Уя |
fr(fe-l,/-i) |
2 |
+ |
fa(ft-i.y) |
_ |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
л (А - 1 , / ) ' . |
|
|
|
(36) |
105"
Уравнения (28), (29), (32), (35) п (36) описывают динамику системы в период нагрева, причем номер / увеличивается по ходу продуктов сгорания — от купола воздухонагревателя. В качестве краевых условии заданы
|
|
|
|
il |
|
|
|
/II |
|
|
âl |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
ri,. |
|
|
|
|
|
|
J j l _ _ |
|
J |
h,0) |
|
||||
|
|
|
|
|
|
|
|
|
|
Ud-o |
||||
|
|
|
|
|
|
|
|
|
l,OJ |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
H* . |
t. |
|
|
|
|
|
|
пи |
|
|
|
|
|
|
|
|
Vue. |
22. |
Структурная |
схема |
но.ілух»жігрспи гол» и |
период пх.чаждспіні |
|||||||||
/д (/г, |
0) |
и /,і Ѵ І (0, |
/) |
для |
любых |
значений |
k |
|
и / |
соответст |
||||
венно. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Структурная схема системы в период охлаждения по |
||||||||||||||
казана на рис. 22. Обозначения па схеме: |
|
|
|
|||||||||||
РО — смесительный |
клапан; |
|
|
|
|
|
|
|||||||
М в 0 — расход дутья па печь; |
|
|
|
|
|
|
||||||||
М'а—расход |
холодного |
дутья |
через |
|
смесительный |
|||||||||
|
|
клапан; |
|
|
|
|
|
|
|
|
|
|
||
М'в — расход дутья |
через |
воздухонагреватель; |
||||||||||||
/ в — температура холодного дутья; |
|
|
|
|
||||||||||
^в.вых—температура |
|
дутья |
после |
воздухонагревателя; |
||||||||||
(ъ — температура |
|
горячего |
дутья, |
поступающего в |
||||||||||
|
|
печь; |
|
|
|
|
|
|
|
|
|
|
|
|
/ц.з — заданная |
температура |
горячего дутья. |
|
|||||||||||
Счет элементов пространственного квантования / ве |
||||||||||||||
дем по ходу |
движения |
воздуха |
( / = |
1 на |
входе воздуха |
|||||||||
в воздухонагреватель, |
/ = т а х |
в районе |
купола |
воздухо- |
||||||||||
пагре вателя). |
|
|
|
|
|
|
|
|
|
|
|
|||
Из балансовых соотношений следует: |
|
|
|
|
||||||||||
|
|
|
Ml(k— |
1) = Мю |
— M'B(k— |
1), |
|
(37) |
||||||
|
|
|
|
|
М„(к— 1)/' + М* (k — Ш |
|
|
|
||||||
|
lB |
( k - |
1) |
= |
— в |
- |
B |
r |
n V |
L s j ü « . . |
(38) |
|||
|
|
|
|
|
|
|
|
Мво |
|
|
|
|
|
ЮГ)
При определении величины lB(k—1), |
так же как ра |
|
нее /д (/г—1,0), полагаем |
(для существенного упрощения |
|
расчетов), что теплоемкости основного |
и подмешиваемо |
|
го потоков одинаковы. |
|
|
Работа автоматического регулятора описывается вы |
||
ражением: |
|
|
M'»(k)~k,\in(k-\)-/h.3 |
+ — y \ ц к - \ ) - І л А , (39) |
A - 1 = 0
где k'p, Tn — настройки регулятора температуры горя чего дутья.
По аналогии с периодом нагрева для периода охлаж дения справедливы уравнения:
'» (Ä — 1,у) = г-тг^- К2 — ßB ) ^в (Ä— 1 , / - 1 ) - Ь
|
^ + |
Рв |
|
|
|
|
+ |
2 ß B ^ K J I ( * - l , / ) l ; |
(40) |
||
/ к л (*;/') = ^кл {k- |
I i ) |
'в ( A — I . / - 1) + / в ( * — ! . / ) |
|||
- И |
|
|
|
||
|
~ / к л |
1,/) |
|
(41) |
|
где |
|
|
|
|
|
FM |
1 |
«ов + |
а в у И в ( / ' ; - 1 ) |
||
Уа |
FAT |
«„в + |
« в |
(* — |
') |
|
|
КЛ |
|<л |
|
|
|
|
|
|
Для моделирования воздухонагревателя на ЭЦВМ составлена программа, содержащая два этапа. Первый этап расчета моделирует процесс выхода воздухонагре вателя на квазистационарный режим из произвольного начального состояния. Этот начальный этап расчета вве ден в программу в связи с тем, что априори неизвестно установившееся распределение температуры по высоте воздухонагревателя в периоды нагрева и охлаждения. Первый этап включает в себя необходимое число циклов «нагрев+охлажденне», когда ЭЦВМ моделирует работу воздухонагревателя без вывода результатов расчета
107
(кроме отметки длительности периода). После установ ления квазпстацноиарпого режима имеется возможность определить установившееся распределение температур кладки и газов по высоте насадки и перевести ЭЦВМ на второй этап расчета—построение термограмм цикла «нагрев-|-охлаждение». При переходе с первого этапа расчета на второй несколько изменяется конечная часть программы расчета.
|
Исходные данные для моделирования получим следующим об |
||||||||||||
разом. |
|
|
|
содержащий 24 г водяных паров |
в 1 ма сухо |
||||||||
|
Влажный |
воздух, |
|||||||||||
го воздуха, имеет теплоемкость |
|
|
|
|
|
|
|
||||||
|
с |
в |
= 0 , 0 2 с н п 1+ !0.98с |
|
ѵ |
= 1018 + |
' |
0,18383/. |
|||||
|
|
|
' |
п 5 и |
в.сух |
|
|
|
|
|
|||
|
Теплоемкость продуктов сгорания |
|
|
|
|
|
|
||||||
с д |
= 0 , 3 0 с с а + |
0 , 0 2 с о + |
О.бЗс^ - I - |
0 , 0 5 с н о |
= 1057,2 + 0,21340/. |
||||||||
|
Теплоемкость |
высокоглиноземистой |
и |
шамотной насадки |
|||||||||
|
|
|
|
|
с к л = |
916,03 + |
0,43618/. |
|
|
|
|||
|
Коэффициент теплоотдачи конвекцией от газов к поверхности |
||||||||||||
насадки |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a* = 0 , 7 4 2 5 t ) ° ' 8 d - 0 ' 3 3 3 |
T 0 ' 2 5 , |
|
|
||||||
где |
аи — скорость газа при 0° С; |
канала |
для прохода |
газов, ,н; |
|||||||||
|
d—гидравлический |
диаметр |
|||||||||||
|
Т—абсолютная |
температура газов. |
|
|
|
|
|||||||
|
Если размер ячеек насадки 45X45 мм и живое сечение 12,4 м2 |
||||||||||||
[32], то |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
о Д 0 = |
0,059/Ид, |
ѵво |
= |
0,064Л/„. |
|
|||||
|
Примем среднюю |
температуру |
продуктов сгорания |
7"д=1273°К |
|||||||||
и воздуха Г В = 9 7 3 ° К - |
Тогда: |
|
|
|
|
|
|
|
« ; = 1 , 5 Ш ° ' 8 ; « ; = 1 , 5 2 < 8 .
Если учесть массивность огнеупорной кладки, то кажущийся коэффициент теплоотдачи можно записать в виде •
|
а * |
где R—половина |
толщины кирпича насадки . (20 мм); |
к—-коэффициент теплопроводности кирпича. Для шамота
X X 1,27.
108
Таким образом: |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
1 + 0,00795Л1Д '8 |
|
|
|||||
|
|
|
|
|
|
|
|
1,52л!"'1 |
|
|
(43) |
|||
|
|
|
|
|
|
|
|
1 + |
0,008Л1°'8 |
|
||||
|
|
|
|
|
|
|
|
|
|
|||||
Линеаризовав |
выражения |
(42) и (43), получим: |
|
|
||||||||||
|
|
|
|
а д |
= |
7,835+ |
0,42546Мд , |
|
(44) |
|||||
|
|
|
|
а „ = |
14,873 + |
0,31232Л1„. |
|
(45) |
||||||
Выражения (44) и (45) справедливы при расходе продуктов |
||||||||||||||
сгорания |
М д = 2 6 н - 6 0 кг/сек, |
|
воздуха — Л!п = 6 4 - ь 9 8 |
кг/сек. |
насад |
|||||||||
Если |
учесть, |
что теплоотдача от продуктов сгорания к |
||||||||||||
ке осуществляется |
также |
изучением, |
то окончательно |
получим |
||||||||||
|
|
|
|
|
|
« |
Д |
А |
28 + |
0,43ЛІд. |
|
(46) |
||
Следовательно: |
с о д |
=1057,2; |
с'а |
=0,2134; |
с о в =1018; |
с'в =0,18383; |
||||||||
с о к л = 9 1 6 , 0 3 ; |
с'кл |
= |
0,43618; |
а о л |
= 2 8 ; |
0 ^ = |
0,43; а о в = |
15; а'в |
=0,31 . |
Рассматриваемый здесь воздухонагреватель имеет параметры: поверхность нагрева / ' = 4 2 470 м2; масса насадки т=2 481 700 кг.
Воздухонагреватель отапливается газом, имеющим теоретическую температуру сгорания 1650°С при номинальном расходе продуктов сгорания Afr =37,5 кгісек. Таким образом,
|
ЛІГ tr =62000. |
|
|
Дутье |
в количестве Л І В 0 = 8 3 , 3 |
кг/сек |
нагревается до 1200° С. |
Максимальная температура верхних |
рядов |
насадки 1300° С, предель |
|
ная температура уходящих продуктов сгорания 400° С. |
|||
Ниже |
приводится программа моделирования воздухонагревателя |
||
на ЭЦВМ |
«Промінь-2». |
|
|
01. |
/ д |
(Л—1,0) |
|
no |
/ |
n, |
у л \ |
02. |
* д ( А - 1 , 1 ) |
||
03. |
t A |
( £ - 1 , 2 ) |
|
1 + л |
tA |
(k~l,n) |
|
51. |
*д |
|
,50) |
52. |
Л!а ( f t - l ) = 0 |
53.ер (1,0)
54.Дт/ТЛ, (0,01)
55.2 = 0
Распределение |
памяти |
ЭЦВМ |
|
|
|
||
ta |
(k—1,0) |
] |
п— |
число |
участ- |
||
t |
іи |
1 n |
|
ков |
разбиения на- |
||
tB(k-l,l) |
\ |
с а д к [ І |
|
( Я т а х = . |
|||
^п |
(А—1,2) |
= 5 0 ) . |
В |
период |
|||
t |
(Ь—і п) |
I |
нагрева здесь хра- |
||||
в |
* |
|
I |
нятся |
^д, |
в пери- |
|
tB |
(k—1,50) |
; |
од охлаждения / в . |
||||
5 6 |
/ к л _ а = 1 2 9 5 ° с |
(с уче |
|
|
|||
том |
работы |
регулятора) |
|
|
|||
57. |
FM/l = ~ ^ = W |
м- |
|
|
|||
5 8 |
а |
о д = 2 8 |
|
|
|
|
|
59. а ' |
=0,43 |
|
|
|
|
|
|
|
д |
|
|
|
|
|
|
109
60.MR(k—1)=37,5
61.fAx/m---2,06
62.M r =37,5
63.Afr /r =620O
64.10000 (адрес)
65.—002 (100—«)=
=—00290 (адрес)
66.00100 (адрес)
67.V = 0
68.—002 (100—л)=
=—00290 (адрес)
69./ у х т а х = 4 0 0 Х
70.(99+rt) 00=10900 (адрес)
71.(99-b«)00=10900 (адрес)
72. |
М'в |
(к—1)=0 |
73- |
k'p |
(0,01) |
74.Дт/7"„ (7,0)
75./ В З =1200°С
Программа
00. |
Чт. |
62. |
21. |
Умн. 152 |
|
01. |
Сл. |
52 |
22. |
Сл. |
151 |
02. |
Зп . |
60 |
23. |
Умн. 60. |
|
03. |
Чт. |
63 |
24. |
Зп . |
172 |
04. |
Д е л . 60. |
25. |
Чт. |
59 |
|
05. |
Зп . |
01 |
26. |
Умн. 60 |
|
06. |
Чт. |
100 |
27. |
Сл. |
58 |
07. |
Выч. 1. 56. |
28. |
З п . |
160 |
|
08. |
Зп . |
164 |
29. |
Умн. 57 |
|
09. |
Сл. |
55 |
30. |
Д е л . 172 |
|
10. |
Зп . |
55 |
31. |
Зп . |
161 |
11. |
Умн. 54 |
32. |
Сл. |
87 |
|
12. |
Сл. |
164 |
33. |
Зп . |
162 |
13. |
У ми. 53 |
34. |
Чт |
I I . 64 |
|
14. |
УП1. 17 |
35. |
Зп . |
167 |
|
15. |
Зп. |
52 |
36. |
Умн. 154 |
|
16. |
БП. |
19 |
37. |
Сл. |
153 |
17. |
Чт. |
00 |
38. Зп. .172 |
||
18. |
Б П . |
15. |
39. |
Чт. |
160 |
19. |
Чт |
I I . 66 . |
40. |
Умн. 61 |
|
20. |
Зп . |
165 |
41. |
Д е л . 172 |
76.a 0 I J =15
77.а ' в = 0 , 3 1
78. М*в |
(к—1)=83,3 |
79.Л 4 В 0 = 8 3 , 3
100.tM
101.ІКЯ (Аг-1,2)
102. |
/ к л |
(ft-1,3) |
149. |
/ к л |
{к—1,50) |
150./' П =150°С
151.с о д = 1 0 5 7 , 2
152.с ' д = 0 , 2 1 3 4
153.со к .,,=916,03
154.с і к л = 0 , 4 3 6 1 8 .
155.с о в = 1 0 1 8
156.с ' в = 0 , 1 8 3 8 3
157.00000
158.00000
моделирования
42. |
Зп. |
163 |
63. |
Чт. |
199 |
|
43. |
Чт. |
167 |
64. |
Зп. |
64 |
|
44. |
Умн. 87 |
65. |
Чт. |
90 |
||
45. |
Умн. |
161 |
66. |
Зп. |
66 |
|
46. |
Зп . |
166 |
67. |
Чт. |
68 |
|
47. |
Чт. |
87 |
68. |
Зп. |
65 |
|
48. |
Выч. |
1.161 |
09. |
Сч. |
П. 157 |
|
49. |
Умн. |
165 |
70. |
Чт. |
( 1 + я ) |
|
50. |
Сл. |
166 |
71. |
Выч. 1.69 |
||
51. |
Д е л . |
162 |
72. |
УП1. 00 |
||
52. |
Зп. |
I I . 65 |
73. |
Чт. |
157 |
|
53. |
Сл. |
165 |
74. |
Ост. 01 |
||
54. |
Умн. 88 |
75. |
Чт. |
150 |
||
55. |
Выч. I. 167 |
76. |
Зп . |
01 |
||
56. |
Умн. |
163 |
77. |
Чт. |
79 |
|
57. |
Сл. |
167 |
78. |
Выч. 1.72 |
||
58. |
Зп. |
I I . 64 |
79. |
Зп. |
78 |
|
59. |
Сч.П. 6-1 |
80. |
Умн. 77 |
|||
60. |
Сч. П. |
66 |
81. |
Сл. |
76 |
|
61. |
Сч. П. |
65 |
82. |
Зп . |
170 |
|
62. |
УПІ . |
19 |
83. |
Чг. |
11.66 |
ПО
84. |
Зп. |
177 |
103. |
Зп . |
171 |
122. |
Выч. Ф. |
90 |
141. |
Сл. |
67 |
|||||
85. |
Умн. |
156 |
104. |
Чт. |
87 |
123. |
Зп . |
70 |
|
142. |
Зп . |
67 |
||||
86. |
Сл . |
155 |
105. |
Выч. |
1.168 |
124. |
Сч. П. 65 |
|
143. |
Умп. |
74 |
|||||
87. |
Умн. |
78 |
106. |
Умп. |
177 |
125. |
УПІ. |
83 |
|
144. |
Сл. |
173 |
||||
88. |
Зп . |
172 |
107. |
Зп . |
176 |
126. |
Чт. |
90 |
|
145. |
Умн. |
73 |
||||
89. |
Чт. |
170 |
108. |
Чт. |
175 |
127. |
Зп . |
66 |
|
146. |
УП1 . |
149 |
||||
90. |
Умн. |
57 |
109. |
Умп. |
87 |
128. |
Чт. |
71 |
|
147. |
Зп . |
72 |
||||
91. |
Д е л . |
172 |
ПО. |
Умп. |
168 |
129. |
Зп . |
70 |
|
148. |
БП.151 |
|||||
92. |
3м. |
168 |
111. |
Сл. |
176 |
130. |
Чт. |
68 |
|
149. |
Чт. |
00 |
||||
93. |
Сл. |
87 |
112. |
Д е л . |
169 |
131. |
Зп . |
65 |
|
150. |
БП. |
147 |
||||
94. |
Зп . 1Б9 |
ИЗ. |
Зп . |
I I . 65 |
132. |
Чт. (1 + н) |
151. |
Сч |
П. |
158 |
||||||
95. Чт. П.70 |
114. |
Сл. |
177 |
133. |
Умн. |
78 |
|
152. Чт. (14-«) |
||||||||
96. |
Зп . |
165 |
115. |
Умп. |
88 |
134. |
Зп . |
174 |
|
153. |
Выч. |
2.75 |
||||
97. |
Умн. |
154 |
116. |
Выч |
|
I. 175 |
135. |
Чт. |
72 |
|
154. |
УП1 . |
77 |
|||
98. |
Сл. |
153 |
117. |
Умн. |
171 |
136. |
Умн. |
150 |
|
155. |
Чт. |
158 |
||||
99. |
Зп . |
172 |
118. |
Сл. |
175 |
137. |
Сл. |
174 |
|
156. |
Ост. 01 |
|||||
100. |
Чт. |
170 |
119. |
Зп . П. 70 |
138. |
Д е л . |
79 |
|
157. |
Чт. |
00 |
|
||||
101. |
Умн. |
61 |
120. |
Сч. П. 66 |
139. |
Выч. |
1.75 |
158. |
Зп . |
157 |
||||||
102. |
Д е л . |
172 |
121. |
Чт. |
70 |
140. |
Зп. |
173 |
|
159. |
Зп. |
158 |
||||
|
|
|
|
|
|
|
|
П е р е х о д к |
команде |
00. Чт. 62 и т.д. |
||||||
|
Приведенная здесь программа предназначена для вывода иссле |
|||||||||||||||
дуемой |
системы |
на |
квазистационарный |
режим. Достижение |
такого |
режима определяется по цифрам, выводимым на печать. Величина в
ячейке 157 |
(команда 74. Ост. |
01) дает |
число интервалов |
времени |
Дт |
|||
в период |
нагрева, |
величина |
в |
ячейке |
158 (команда 156. |
Ост. |
01) |
— |
число интервалов |
времени Дт |
в период охлаждения. Когда эти |
числа |
|||||
начинают |
попарно |
повторяться, |
это |
свидетельствует о |
достижении |
квазистационарного режима. После |
команд |
76. Ост. 00 и 156. Ост. 00 |
|
имеется возможность из ячеек |
100 |
-^-(99+я) получить установившее |
|
ся распределение температур |
по высоте |
насадки соответственно в |
конце периодов нагрева и охлаждения.
Для изучения динамики нагрева и охлаждения воздухонагрева теля внутри периодов иагрева и охлаждения программа расчета дол жна быть несколько изменена:
|
Д л я п е р и о д а |
|
н а г р е в а |
Д л я п е р и о д а |
о х л а ж д е н и я |
||||||
69. |
Чт. |
01 |
74. |
Ост. |
01 |
151. |
Чт. |
( 1 + п ) |
156. |
Ост. |
01 |
70. |
Ост. |
01 |
75. |
Чт. |
( 9 9 + « ) |
152. |
Ост. |
01 |
157. |
Чт. |
( 9 9 + я ) |
71. |
Чт. |
( 1 + н ) |
76. |
Ост. |
01 |
153.. |
Чт. |
173 |
158. |
Ост. |
01 |
72. |
Ост. |
01 |
77. |
БП. |
00 |
154. |
Ост. |
01 |
159. |
БП. |
77 |
73. |
Чт. |
100 |
|
|
|
155. |
Чт. |
100 |
|
|
|
(команды 75.4т. 150, 76.3п.01 и 77.Чт.79 для периода охлаждения вос станавливаются. Пуск осуществляется с команды 75). Остальная программа во время изучения динамики системы остается без изме нения.
111