Файл: Динамика и управление ядерным ракетным двигателем [Текст] 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б , bjLCftQttifoo 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

 

 

 

 

 

 

 

Я

( *

,

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