Файл: Динамика и управление ядерным ракетным двигателем [Текст] 1974. - 253 с.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 23.06.2024
Просмотров: 234
Скачиваний: 0
3. Жуков В. П. О возможности использования уравнений кинетики для распределенного описания нестационарных изменений малых отклонений плотности нейтронов, вызванных неравно мерным возмущением. — В сб. «Автоматическое регулирование двигателей летательных аппаратов», М., ЦИАМ, 1968, № 63, с. 50.
4.Потапенко П. Т. Пространственно-зависимая передаточная функция ядерного реактора. — В сб. «Управление ядерными энергетическими установками». Вып. 4. М., Атомиздат, 1970, с. 5.
5.Бассард Р., Делауэр Р. Ядерные двигатели для самолетов и ракет. Пер. с англ.М., Оборонгиз, 1967, с. 115.
6.Абрамович Г. Н. Прикладная газовая динамика. М. , «Наука», 1969, с. 225.
7.Петров П. А. Ядерные энергетические установки. М. — Л., Госэнергоиздат, 1958, с. 128.
8.Лойцянский Л. Г. Механика жидкости и газа. М., «Наука», 1970, стр. 128.
9.Гухман А. А., Илюхин Н. В. Основы учения о теплопередаче при движении газа с большой скоростью. М., Физматгиз, 1957.
10.Лыков А. В. Теория теплопроводности. М., «Высшая школа», 1967.
11.Березин И. С., Жидков Н. П. Методы вычислений. М., Физмат гиз, 1960.
12.Марчук Г. И. Численные методы расчета ядерных реакторов. М., Атомиздат, 1957.
13.Мак-Кракен Д., Дорн У. Численные методы и программирование на ФОРТРАНе. Пер. с англ. М., «Мир», 1969.
14.Абрамов Н. М. и др. Условия соблюдения устойчивости при чис
ленном счете нестационарных термодинамических уравнений на ЦВМ. — Веб. «Автоматическое регулирование авиадвигате лей», М., ЦИАМ, 1966, № 52, с. 98.
15.Камке Э. Справочник по обыкновенным дифференциальным уравнениям. М., Физматгиз, 1961.
16.Воронов А. А. Основы теории автоматического управления. М. — Л., «Энергия», 1970.
Глава 3
ДИНАМИКА ЭЛЕМЕНТОВ СИСТЕМЫ ПОДАЧИ РАБОЧЕГО ТЕЛА
§ 6. Теплообменные устройства, элементы охлаждения
Общая часть. Тракт рабочего тела в ЯРД содержит участки, на которых происходит теплообмен с другими участками (через стенки) или с элементами конструкции. Рассмотренный в гл. 2 участок прохождения теплоносителя через твэл относится к числу теплообменных, но при его изучении можно было не учитывать влияние на тепловое состояние твэла потока тепла от твэла к теплоносителю охлаждающих участков (поскольку этот поток значительно меньше основного потока тепла от твэла к теплоносителю в каналах последнего). Это упрощение значительно облег чало изучение динамики системы твэл — теплоноситель, ибо передаточные функции от входных воздействий к от клонению температуры теплоносителя подчиняются линей ному дифференциальному уравнению 1-го порядка (по х) и, следовательно, могут быть выражены в квадратурах через параметры теплоносителя и твэла. В этом параграфе рас смотрена динамика (вблизи номинального режима) более сложной теплообменной системы, а именно, системы из двух параллельных, противоположно направленных потоков теп лоносителя 1 и 2 на рис. 3.1, обменивающихся теплом через стенку б, имеющих внутреннее тепловыделение и воспри
нимающих |
тепло через |
стенку а от среды |
тс температурой |
||||
Г т (х, |
t). |
На |
рис. 3.1 |
обозначено: Т} (j |
— 1, 2) — темпе |
||
ратуры |
потоков; |
jTa, Тб — температуры |
стенок соответст |
||||
венно а и б; |
<7б1, |
<7б2 — потоки тепла |
(на единицу длины) |
||||
от стенки |
б к потокам |
соответственно |
1 и 2; qal — поток |
тепла от стенки а к потоку 1; qTa — поток тепла от среды т к стенке а. Схема теплообмена, изображенная на рис. 3.1, может иметь место, например, для системы охлаждения твэла потоками теплоносителя, проходящими по каналам замедлителя. В этом случае роль стенки б играет замедли
72
тель, стенки а — прилегающая к замедлителю часть актив ной зоны твэла, роль <7та — поток тепла изнутри активной зоны твэла к стенке а, роль Т т — температура твэла.
Дифференциальные уравнения теплообмена. Уравнения теплообмена для системы рис. 3.1 составим при тех же упро щающих предположениях, которые использованы в § 2 при составлении одномерного описания динамики системы
твэл + |
теплоноситель; при этих предположениях уравнения |
||||||||||
теплообмена имеют |
вид |
|
|
|
|
|
|||||
|
Сі тх 1(дтJdx) ѵх + (дТJdf)] = qal + qx+ q6x, |
||||||||||
|
c2 tn2 [{— dTjdx) v2 + (dT2/dt)] = q2 + qG2; |
|
|||||||||
|
|
|
ma (dTJdt) = qa + qTa— qal; |
|
|
||||||
|
|
|
Cfj mö (dTftldt) -- q6 |
q61 |
q62, |
|
|
||||
где |
Cj, |
m.j |
( / = 1, 2) —соответственно |
|
i=o |
||||||
удельная теплоемкость |
и |
масса |
на |
|
|
||||||
единицу длины теплоносителя в пото |
|
|
|||||||||
ке I или 2\ Ck, mu (k = |
а, б) — соот |
|
|
||||||||
ветственно |
удельная |
теплоемкость |
и |
|
|
||||||
масса на |
единицу длины материала |
|
|
||||||||
стенки а или б; qt (I = |
1, 2 , а, б) — |
|
|
||||||||
тепловыделения |
на |
единицу |
длины |
|
|
||||||
в соответствующем потоке |
или стен |
|
|
||||||||
ке; |
qkj |
~ akj IT-kj (Тk |
Тj), |
осkj |
|
|
|
||||
коэффициент теплоотдачи; |
II w — пе |
|
|
||||||||
риметр |
теплоотдачи |
для |
соответст |
|
|
||||||
вующей пары (k = |
а, |
б, т; / |
= 1, 2 , а). |
|
|
||||||
Линеаризуем уравнения (3.1) при |
|
|
|||||||||
упрощающих предположениях, кото |
|
|
|||||||||
рые были приняты в § 3; в частно |
|
|
|||||||||
сти, |
коэффициент |
а та будем |
считать |
|
|
||||||
постоянным, а остальные коэффициен |
|
|
|||||||||
ты |
at] |
— зависящими |
от |
расхода |
и |
|
|
||||
температуры соответствующего пото |
|
|
|||||||||
ка |
теплоносителя согласно формуле |
|
|
||||||||
(2.24); |
тепловыделение |
в |
каждом |
из |
Рис. |
3.1. Расчетная |
|||||
потоков или стенок |
и |
при |
каждом |
схема |
теплообмен |
||||||
X £ |
[О, |
Л |
будем |
считать |
пропорцио |
|
ника. |
нальным мощности реактора с коэф фициентом пропорциональности, не зависящим от времени.
При линеаризации уравнений (3.1) используем соотно шения для номинального режима
73
r'io G10(dTxojdx) — (?aio -f- q10+ Ябіо\
^20 ^20 {dTzo/dx) |
УъО“И ^620’ |
|
|
|
(3.2) |
0 |
9 а0 + ^таО |
QalO> |
О = ( 7бо— Ѵбіо |
9б20> |
где Gj0 (/ = 1,2) — расходы теплоносителя в потоке 1 или 2 на номинальном режиме. Линеаризованные уравнения теп
лообмена |
имеют вид |
|
|
|
ах[(dfly/dx) + |
(1/п10) (düjdt)] + |
|
|
+ {da\ ~\~ d61+ |
ex |
laX /бг) 'ö'i = |
-■ = dal + |
d61 6 + <710 V + |
(0,8 qaW+ 0,8 <76 1 0 —gx) yx, |
a2 [(— d b 2jdx) + (l/n2o) (dflydOl +
|
+ |
(^62 |
e2— 4г) 'ö’a — d62,ö6 + <72oV + |
|
||
|
|
|
+ (0.8 <7б20 + ёг) Ѵгі |
(3 .3 ) |
||
|
|
ba (döjdt) + |
(rfal + dT3) |
= |
|
|
|
= |
+ <7a0 V + K l — kl) #1 — 0.8 <7al0 Уі, |
|
|||
|
|
bft (ölig/dt) -)- (rfg1+ d52) ^g — |
|
|||
|
|
= K i — ^6i) ®i + К г —^бг) #2 + |
|
|||
|
|
-!- gG0 V — 0,8 qm yx— 0,8q620 y2, |
|
|||
где |
O’; = |
(АТі/Тв) (I = |
1, 2 , а, б) — отклонения темпера |
|||
тур |
потоков или стенок, отнесенные к некоторому базово |
|||||
му |
значению |
Th, yj = |
(AGj/Gj0) |
( / = 1 , 2), |
v = (Aql/ql0) |
|
( / = 1, і 2 , |
а, б)— относительные |
отклонения |
тепловыделе |
ний в потоках и в стенках (согласно сделанным предположе ниям у} и V не зависят от х); ■Ѳ,т — АТт/Твт, Твт — базо вое значение для Тг,
0-j CjQGjo Tб , bjL— CftQttifoo TBt ë) = Cjo Gjo (dTjo/dx);
Iki = (r(Tj)/f(Tj)) qkj Tb', dhi ~ ahi K ; Tb;
( /= 1 ,2 ; k — a, б, t),
74
d r &та П та Г бт-
Переходя в (3.3) к изображениям по Лапласу при нулевых начальных условиях и исключая d* (k = а, б), получаем
дифференциальные |
уравнения относительно |
|
|
(/ |
= |
1, |
2): |
||||||||||||||||
|
|
(d&Jdx) + 1{р/ѵ10) + Ап ]А = |
|
|
|
|
|
|
|
||||||||||||||
|
— Nt v + |
|
|
|
|
AiYi |
|
A 2Y2A АІ12 А ’ |
|
|
(3.4) |
||||||||||||
|
— (d b2/dx) -f l(p/v20) + Л22] Ä = |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
— A^2 V |
|
|
A l Y l |
|
B 22 У2А A l А> |
|
|
|
|
|
|
||||||||||
где |
|
|
|
|
|
|
|
|
|
|
Jai" ■^ai |
|
^aP + dTa |
|
|
|
|||||||
А і (ж, р) |
|
|
1 |
|
de ю |
|
|
|
|
|
|||||||||||||
|
с10 |
|
dx |
|
|
ai |
|
|
ba p + dTa 4" dai |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
dßi—Ібі |
|
bg P + dö2 |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
ai |
|
^öP + dgi+dßj |
|
|
|
|
|
|
|
||||||
A |
2 |
( xp), |
=d01 (d62—/ |
в2) |
/ [ а |
! |
( 6 |
Öp |
|
!- d |
61 + |
d |
62)J; |
||||||||||
А |
i ( |
x |
, |
p ) |
= |
|
( d |
6 2 |
|
/ f l a ) K |
rf6 i — |
А |
Ѵ |
|
^ |
б Р |
|
+ А і |
|||||
^22 |
|
P) ““ |
|
|
1 |
|
, |
^62 |
^б2 |
^бР“Ь^б1 |
|
|
|
||||||||||
|
|
c2o |
dx |
|
|
&<i |
|
|
bg P + döi + dg2 |
|
|
||||||||||||
Bn (x, p) = |
- ^ |
|
|
|
|
|
|
||||||||||||||||
|
|
dTlt> |
|
0,8^aio |
|
|
ba P+ dTa |
|
|
|
|||||||||||||
|
|
dx |
|
|
аі |
|
|
|
P + dai + ^та |
|
|
||||||||||||
|
|
|
|
|
|
0,8 <7бю |
|
Ö6 P + d62 |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
ai |
|
*6 P+ d6i + d62 |
|
|
|
|
|
|
|
||||||
Я |
2і |
( * |
, |
P |
) |
= |
|
° > |
8 |
962 0 |
d |
6 1 /Рf+f l id a\( |
А+ |
4 |
i |
) ] . |
|
|
|||||
|
ß |
2! |
( x |
, |
p |
) |
= |
760io, 8d62/(b6< |
p |
+ |
A |
|
id62),+ |
|
|
|
|||||||
ß 22 (x, p) ■= |
|
|
|
1 dT2о |
0 ,8 <?б20 |
|
^ ö P + dg i |
|
|
||||||||||||||
|
|
|
|
dx |
|
a2 |
|
P + döi + dß2 |
|
||||||||||||||
|
|
|
1 |
|
|
|
|
|
|
|
|
||||||||||||
Ni(x, p) = |
|
qwj |
+ |
- |
Pao dai |
|
|
|
|
|
Рбо dg1 |
|
|
|
|||||||||
— |
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
ö i |
\ |
|
|
ba p + rfTa 4* ^ a i |
P + d6i + d62 |
|
|
|||||||||||||
|
A (* > |
p) = |
|
— |
(<7so+ |
■ |
|
?6o d62 |
|
|
|
|
|
|
|
||||||||
|
|
|
, , |
, . |
|
|
|
|
|
|
|||||||||||||
|
|
(x, p) = |
|
Ö2 |
V |
|
|
bgP +d6i + dö2 |
|
|
|
|
|
||||||||||
|
А |
dTdal/[a1(ba p+ |
dra + A |
i |
) |
l . |
|
|
|
||||||||||||||
Граничные условия для этих уравнений: |
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
М О , |
0 |
= » п (0 ; |
|
|
|
|
|
|
|
(3.5) |
|||||||
|
|
|
|
|
|
# 2 ( / t), |
= |
A |
i |
|
( 0 > |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
75
где ftа (t) — ATji/Тв (j = 1, 2) — отнесенные к Тъ откло нения входных температур потоков теплоносителя. Урав нения (3.4) при переменных коэффициентах Ajr, BjT, Njt Zx (/, г — 1, 2), как правило, не интегрируются в квадра турах, поэтому определение передаточных функций от входных воздействий ftri, уг (г —1, 2), ѵ и йт к координатам ■бу- (/ = I, 2) проводится путем составления на основании (3.4) дифференциальных уравнений относительно искомых передаточных функций (отдельно для каждого из воздей ствий) и интегрирования этих уравнений с соответствую щими граничными условиями.
Воздействие ^ f. Ясно, |
что реакции координат ftj |
(х, t) |
|||
( / = 1, 2) |
на |
йу* имеют |
чистое запаздывание |
(х) = |
|
X |
|
|
|
|
|
= jd£/[t»10 |
(£)]. |
Поэтому |
изображения этих |
реакций на |
|
о |
|
|
|
|
|
і35і имеет смысл искать в виде |
|
|
|||
öj(x, |
р) = Фп (х,р)ехр( — рт1)Ъи (р), |
j = l , 2 , |
(3.6) |
где функции ФЛ (х, р) удовлетворяют однородной системе дифференциальных уравнений
дФи /дх + Лп ( X , р) Фп = Л12 (х, р) Ф21;
(—дФгіІдх) + [{plv10(х)) + |
(3.7) |
+(рІѵло (х)) + л 22 (х, р)} Ф21 == л 21 (х, р) Фп
сграничными условиями
|
Фп (О, р) = 1; |
Ф*1 (/, р) = 0. |
- |
(3.8) |
|
[Можно показать, что вектор-функция (Фд (х, |
р), j — 1, |
||
2) |
определяемая уравнениями (3.7) с граничными условия |
|||
ми |
(3.8), — обобщенная производная регулярного ориги |
|||
нала, но теряет это свойство |
при умножении |
на ехр (ро), |
где а — любое положительное число. Это означает, что чистое запаздывание реакции (ftj (х, t), / = 1, 2) на Оу,- (t) равно точно тх (х).] Необходимое и достаточное условие одно
значной разрешимости краевой задачи 1(3.7), (3.8)] |
— |
от |
|
личие от нуля величины Ф 'п (/, р), где (Ф 'л (I, р), |
/ = |
1, |
|
2) — решение задачи |
Коши |
|
|
Ф'п (0 , |
р) = 1; Ф 'п (0 , р) = 0 , |
(3.9) |
76