Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Динамика и управление ядерным ракетным двигателем [Текст] 1974. - 253 с

.pdf
Скачиваний:
38
Добавлен:
23.10.2023
Размер:
8.52 Mб
Скачать

где Qu (i = 1, 2; / = 1 ч- 3) — желаемые передаточные функции замкнутой системы; будем считать их рациональ­

ными дробями от р.

 

 

 

 

М 21,

М 32, D 31, Ф12,

Если передаточные функции D J3,

Ф 23 отличны от нуля, то матрицы

 

 

 

I

1

Dla

D13

0

0

\

/

А 1

1

А 3

^21

0

\

 

£>81

о

1

О

- М 32

\

1

О

о

о

о

 

\

0

1

0

о

о

 

( 1

- ф »

0

\

п „ - Ал - Ф 22 - Ф 2з

 

\Аі о

о

/

 

не вырождены. При этом условие (4.28) разрешимости задачи синтеза выполнено при любом выборе Qu. Общий вид систе­ мы уравнений регулятора можно получить, согласно изло­ женной выше методике, исключением из системы уравне­

ний величин фу, фт, фу, заменяющих

в уравнениях (4.39)

и (4.40) величины соответственно

бР , /т, /у.

 

Фу

Ф і2 Фт

 

= —£>і2 8Р—Di3 бп;

О гі фу

— Ф 22 фт

— Ф гзФ у —^ 2 і

=

б P—D23 б / г;

D3l Фу

 

М 32 Hg*

=

 

б п\

фу

Q12 Фт

— Ql3 фу

 

 

Qu /1! ■

 

---Q22 фт

--- Q23 Фу

=

- б

Р

+Qn h\

 

 

дн,М +0дР_

^ д и р б Р 0 д в б п ' Ф дн/ 1

 

 

 

 

 

 

(4.41)

(4.41) ß — вектор размерности

b,

b ^

0; Мдн, £ д>

£>ДиР> £>дВ. ф Ди — матрицы с 6 строками.]

 

137

При

Q12 Qi3

ф 23Ф12

матрица

 

Q22 Qi'i

 

 

 

1

- Ф 12

0

0

О

 

П2і --Ф22

Ф23

—A4:2i

О

 

D31

0

0

0

'32

 

 

1

Qx2

Qi3

0

о

 

 

0

Q22

Q23

0

о

не вырождена, поэтому можно выбрать b = 0 и построить регулятор по схеме прямой передачи согласно соотноше­ ниям, получаемым путем решения системы (4.41) с Ь = О

относительно і|)у,і|)т, i|v, p-j, pa*:

Нѵ = 6Р + КѴІ бп -h Ln fi’

(4.42)

M'öfc = K22 4 K23 L21 fx,

где

Ки -~(1/Л421 А) [Ф]2 Ф23 + А2з (1 —Dl2 D2i) —Ql2 Ф23 +

+Qi3 22D2l Ф12) —Q22 D12Ф23 -I- Q23 (П12 Ф22—Фхг)1> А^із = (1/Al21 А) [A23 (D23D13 D21)Q22D13 Ф2з +

л Q23 33 Ф22—ПгзФіг)];

Агг - - — (D3J M 32А) (А2з Dn + Ql3Ф12);

К23- (1/М32 А) [А23 (1 D13 D31) Q23 Фі2];

Lu = (\/M2l А) [Aj3 (Ф22—Ö21 Фі2) —(Q21Ф12 + А12) Ф2з] ; L2l = — (Ф32/А4з2) (А13/ А) D31;

Qu

Qij

(і> / —1 "г* 3), А — А2з

Q23 Фіг*

*u = (?2j

Qzj

 

 

 

Если А =

0, т. е. матрица

вырождена,

то регулятор

не может быть построен по схеме прямой передачи. При этом, если каждый из столбцов матрицы Мп линейно незави-

сим от столбцов

матрицы II Пн\I , т.

е. если в каждой из пар

(<2із> Q 2 3 ) и

[ Ф 1 2 Ф 2 3

С іг Ф г з

“ Ь

<3із ( Ф 22

П 21Ф12),

Q23 2 2— П 21Ф12) — ФггФгзІ хотя

бы

одна из составляю­

щих отлична от нуля, регулятор может быть построен без внутренних координат (но с перекрестными связями между выходами регулирующих органов). В противном случае (например, при требовании абсолютной инвариантности

138

обеих регулируемых

координат относительно воздействия

/,/) необходимо ввести внутренние координаты.

 

С и н тез р е г у л я т о р а

д л я

Я Р Д

по з а д а н н ы м

ст а т и ч еск и м

со о т н о ш е н и я м . Отдельно

рассмотрим задачу

обеспечения

в ЯРД заданных статических

соотношений

вида (4.40),

в частности, потребуем астатизма регулируемых координат бРу и ÖP по отношению к /т и fj, т. е. обращения в нуль при

р = 0 передаточных функций Qtj (і =

1, 2; j = 2, 3). В ста­

тике уравнения объекта (4.39) принимают вид

 

 

6РУ + Р і2

(0) ЬР + D 13(0) бп «= Ф12(0)/т,

 

 

 

D2i(0)6Py-|-

бР +О2з(0)бгс =Ф22(0)(т +Фгз(0)(ур4^2і (0)[Ху;

Р зі бРу

 

=

 

 

 

(4.43)

 

 

 

 

 

 

 

где Du (0), Фп (0){і =

1, 2; / = 1 -

3;

/

= 2,

3), Л421 (0),

Dзі,

ЛІ32 — числовые

коэффициенты; будем считать, что

D13

(0),D§1( Af21 (0), M l2, Ф12 (0) и Ф 23 (0) отличны от нуля.

Система (4.41)

в статике принимает вид (при 6 = 0)

 

\|)у —Ф12 (б)

—Di2(0) 6Р

—£>із(0) 6n;j

D2i(0) i|5y —Ф22 (0) 4'т —Фгз(0)'Фу —442і (0) fly

=

—бР —02з(0)б/і;|

D31 il’y

 

Л4"2

 

=

0;

j

 

4 y

 

 

 

=

Q11 (0) fi -,

 

 

 

 

0

=

- 6 P

ч Qal

 

 

 

 

 

 

 

(4.44)

где Qn (0), Q21 (0) — заданные статические коэффициенты передачи от / х соответственно к бРу и бР.

 

о

 

 

Л421(0)

 

Столбец

0

здесь зависим,

 

о

 

 

о

 

0

0

а столбец

независим от столбцов матрицы

о

о

139

 

' 1

--Фі2:(0)

о

¥ п

Dt l (0)

--Ф22 (0)

Ф 23(0)

Д°,і

0

о

 

1

0

о

 

0

0

о

поэтому регулятор должен содержать по крайней мере одно интегрирующее звено. Из системы (4.44) видно, что для обес-

печения астатизма хр =

/ 0 Р у \

I / т \

I gp ) по отношению

к / н == I* j

и заданных значений Qjx (0)(/ ==1,2) регулятор должен иметь уравнения следующего общего вида:

 

 

N

 

бР

= о,

(4.45)

 

 

М'бі

І ап ) - ^

 

 

 

 

 

где N , К,

L — матрицы с двумя строками и элементами из

поля W, удовлетворяющие условиям: N (0)

/0 0ѵ

~

\1 0'

(0) =

 

 

 

Dt г Qu (0)

 

0

0

, L (0) =

невырожденная

— 1

0

Q21 (0)

 

матрица 2x2 .

Один из простейших по структуре законов управления, удовлетворяющих этим условиям, имеет вид

І*/ = К і2 (Р) ЬР + Рі Л р) fl,

(4.46)

M'öé ~~ ^21 (Р) fl.

где /\12 (0)— оо, — I р=0 _

Qn (0), L21 (0) — jj1 Qn (0);

Ді2

M32

функция К12 (р) выбирается такой, чтобы обеспечить устой­ чивость замкнутой системы объект — регулятор.

С и н тез р е г у л я т о р а с у ч ет о м в о зм у щ ен и я п о р еа к т и в н о с т и .

Рассмотрим задачу синтеза регулятора для того же объек­ та, но с иным набором учитываемых внешних воздействий; уравнения объекта теперь будем считать имеющими вид

ÖPy

+ Di2 8P

\-Dl3 6n -= 0;

 

\I

D2 1 bPy

+

бP

4 - D23 8n -

Ф28 fy

+ M2lpy;

(4.47)

D3l bPy

+

 

бn ;

Ф22 fbk

+ Мз2р6Ь

1

140

гдеП 31, Л432, Ф 32 пропорциональны Wvn;fök — возмущение по реактивности реактора, которое считаем неизмеряемым. Выбор управляющих воздействий, регулируемых и изме­ ряемых координат оставим тем же, что и раньше.

Соотношения (4.40) для замкнутой системы теперь имеют вид

Ру — Q u f i^ Q iz f f t k

Q u fj',

 

&öp — Q l / l + Q22/6ft

'IIQ23/./, II

(4.48)

где Qn(i = 1, 2; / 2

 

синтезе

= 1 -у- 3) — задаваемые при

рациональные дроби от р.

Будем считать отличными от нуля передаточные функции Dis, Л421, М 32, Ф 23, Ф 32, что обеспечивает невырожденность

матриц (По

I

1

0

0

\

и Пн =

Ап

0

—Фаз ). При этом условие

 

VA.1 — Ф 32

О

/

(4.28) разрешимости задачи синтеза выполняется при

любом

выборе

Qi 1-

 

При

Q12

Q13

Ф 0 матрица

 

 

1

 

0

 

0

0

0 \

Q22

Q23

 

 

/

 

 

 

 

 

 

Ö21

 

0

—Ф23 -М 21 о

 

 

 

 

 

 

 

 

 

Яп =

 

Пзі

—Ф32

0

0

—М32

не вырождена; поэто­

 

 

1

--Q12 — Q13

0

0

 

 

 

 

V

0

Q22 —Q23

0

0 /

 

по схеме прямой пе­

му регулятор

может быть

построен

редачи согласно соотношениям (4.42),

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

(1 /Л 4 2 і)

[1 —

D la

D 21----( Ф 23/ ^ 2 з )

(Q l2 + Q 22 П 12) ];

/C is = (1 /Л 4 21) [ D 23

П 13 D 21

( 1 / А 23)

Q22 Ф 23 D 13];

Бц —

 

Ф23 ^ 1 2 / ^ 2 1

 

 

 

 

 

(4.49)

^ 2 2 =

( 1 / 4 4 32) [

D 12 D 31 г

( Ф 32^^2з) ( П 12 С ;* H“ Q is)l>

^ 2 3 ~

( 1 / / И 32) [1 —

D 13 D 31 +

( 1 / Aga) Q23 D 18 Ф з2];

A21

^13 Фза/^23 M32!

А-;:

Qu Q u!

(С /

1 - 3 ) .

 

 

Qzj Q2J I

При А23 = 0 регулятор нельзя построить по схеме пря­ мой передачи. Схема регулятора может не содержать внут­ ренних координат, если каждый из столбцов матрицы Мп

независим от столбцов матрицы ^ Hj , т. е. если в каждой из

нар (Q12, Q22) и (Q13, Q23) хотя бы одна из составляющих от­

141

лична от нуля, ß противном случае, когда требуется аб­ солютная инвариантность обеих регулируемых координат по отношению хотя бы к одному из воздействий /т и не­ обходимо ввести в схему регулятора внутренние коорди­ наты.

Пусть

теперь

для

объекта

(4.47)

требуется

обеспечить

астатизм

хр = ^

I

 

 

 

г

//т\

"

 

по отношению к / н =

 

j

и заданные

статические коэффициенты передачи

QJX (0)(/

=

1, 2) соот­

ветственно от бРу и ЬР к /j.

 

 

 

 

 

 

Уравнения (4,"47) в статике имеют вид

 

 

 

 

 

6Py+ D12(0)ÖP +

Dl8(0)6/i = 0;

 

 

 

)

ö 2i(0)6ßy+

 

 

D23 (0) ön —- Ф23 (0) fj -J- УИ2І (0) [Ху; I

D ° t l

6Py

 

 

 

 

 

Н-7И°2рбь

J

 

 

 

 

 

 

 

 

 

 

(4.50)

где

Dij (0)(t = 1 ,

2;

j

1

3),

M 21(0),

Ф23 (0),

'31.

Фз2 — числовые коэффициенты, из которых D13

(0), М 21

(0),

М %2 отличны от нуля,

что обеспечивает выполнение для ста­

тики

условия (4.28).

Считая, кроме того, отличными

от

 

 

 

 

получаем, что каждый из столбцов

 

 

 

 

зависит

от столбцов матрицы

 

При этих условиях регулятор

должен содержать два интегрирующих звена.

Уравнения регулятора, обеспечивающего выполнение

сформулированных

выше требований

к поведению системы

в статике,

имеют

общий

вид

(4.45),

где N (0) — 0,

12 (0) D13 (0)\

 

т

=

/Qn(0)\

~

невырожденная

А (°) = (

1

п

),

 

L

 

, I —

матрица 2 x 2 .

0

1'

 

 

\Qn(fl)J '

законов управления.,

Один

из

простейших

удовлетворяющих этим условиям, имеет вид

 

 

 

Рб к “ ТС22 б Р Ф ^С23

Ф

f i t і

142

где

K n (0) = оо;

LU

Р=О -

Qzi (0); ^ 2 3 (°) = °°;

Кгг

 

 

■^22

__ ^12 (0) .

£21

Q11

(0) ■

^23

р=0 £>13 (0)

^23

Di3

(0) ’

передаточные функции КЛ2 и /(22 выбираются так, чтобы обеспечить устойчивость замкнутой системы.

Сравнение законов управления, получаемых при двух разных наборах учитываемых воздействий [(Д, /т, fj) и (/ъ /б/г. /,/)1 показывает, что учет возмущения /в* обусловли­ вает необходимость использования в регуляторе информа­ ции о 6«; если же ожидаемый уровень возмущений по реак­ тивности значительно ниже уровня возмущений со сторо­ ны тракта рабочего тела, то закон управления может быть взят в виде (4.46), использующем информацию лишь о (6Р, /і) и требующем лишь одного интегрирующего звена (для получения астатизма 6Р по отношению к неизмеряемым возмущениям).

§13. Динамика ЯРД на пусковых режимах

ирежимах останова

Пусковые режимы. Пусковые режимы довольно сложны во многих отношениях и предъявляют более жесткие тре­ бования к конструкции двигателя, чем режим номинальный. Во время пуска двигателя должны быть осуществлены сле­ дующие процессы.

Физический пуск реактора. Задача этого процесса — вывод реактора из подкритического состояния в крити­ ческое и подъем мощности до нескольких десятков ватт, после чего начинается следующий процесс — разгон дви­ гателя. Во время физического пуска плотность нейтронов в реакторе возрастает во много раз—от стационарного уров­ ня, который создается в реакторе источником, до уровня, при котором тепловая мощность достигает нескольких ватт.

Разгон двигателя (или вывод на рабочий режим). Во вре­ мя этого процесса происходит подъем мощности до номи­ нального уровня. Одновременно организуется отвод тепла из реактора путем подачи в него рабочего тела, для чего осу­ ществляется, как и в ЖРД, раскрутка турбонасосного агрегата или включение газобаллонной системы подачи рабочего тела.

143

Установление номинального режима. В конце разгона параметры двигателя приближаются к номинальным. При этом важно сделать их на этом уровне стабильными. Этот ответственный процесс обычно осуществляется совместной работой пусковой системы управления и системы управле­ ния номинального режима. Рассмотрим этапы пуска под­ робнее.

Физический пуск реактора. Плотность нейтронов в реак­ торе перед началом физического пуска определяется источ­ ником нейтронов и величиной подкритичности [13] и равна

п = sl*/(Ak),

где п — плотность нейтронов, нейтрон!сма\ s — мощность источника, нейтрон!(см3 -сек); I* — среднее время жизни мгновенных нейтронов. На время физического пуска не накладывается никаких ограничений со стороны тепловых процессов. Поэтому увеличение плотности нейтронов про­ изводится в таком темпе, какой диктуется оператором либо системой управления. Обычное ограничение во время физи­ ческого пуска — ограничение на период реактора

Т = nj(dnldt) « (ß— Sk)!(k6k)

в том смысле, что период реактора не должен быть меньше некоторой величины Т мин. Эта величина определяется вре­ менем, необходимым для остановки реактора в случае пус­ ковой аварии [13].

Разгон двигателя. В этот период действуют следующие ограничения.

Ограничение на период реактора.

Ограничение на скорость нарастания температуры реактора. Это ограничение при разогреве определяется воз­ никающими в конструкции реактора механическими напря­ жениями. Причина их — температурные градиенты. Смысл данного ограничения состоит в том, чтобы напряжения не превысили опасных значений.

Ограничение на максимальную температуру. В процессе разгона температура реактора, твэлов, газовых трактов и других элементов двигателя не должна превышать пре­ дельных значений.

Ограничение на амплитуду и частоту колебаний в газо­ вых трактах двигателя. Колебания в газовых трактах, обя­ занные своим возникновением волновым явлениям, не долж­ ны достигать опасных значений, а частота их не должна ле­

144

жать в резонансных областях для элементов конструкции двигателя.

Ограничение на время разгона. Время разгона не должно быть слишком велико, так как это приводит к чрезмерному непродуктивному расходу рабочего тела.

Кроме этих существует еще ряд технических ограниче­ ний, ограничений по работе с активными материалами, огра­ ничений по технике безопасности и т. п.

В этом параграфе рассмотрены некоторые особенности периода разгона двигателя, связанные с распределенностью

Рис. 4.3. Распределение тепловыделения в твэле.

параметров. Во время разгона в реакторе, во-первых, про­ исходит интенсивное повышение температуры, определяе­ мое подъемом мощности и, во-вторых, перемещение макси­ мума температуры к выходному концу реактора. Остановим­ ся подробнее на последнем явлении. Максимум тепловыделе­ ния имеет место в центральной части реактора, как это пред­ ставлено на рис. 4.3. В том случае, когда разогрев реактора идет без подачи рабочего тела, профиль температур прибли­ зительно повторяет профиль тепловыделения. Отклонение профиля температур от профиля тепловыделения происходит вследствие теплопроводности по длине реактора, которая весьма мала. Поэтому, чем быстрее протекает разогрев реак­ тора, тем ближе профиль температур повторяет профиль тепловыделения. Если по достижении определенной темпе­ ратуры прогрева включить расход рабочего тела, то в пер­ вый момент профиль температур рабочего тела повторит (с естественными отклонениями) профиль температур тепло­ выделяющего элемента. Отклонения температур будут оп­ ределяться динамикой теплопередачи: в передней части твэла, где имеется нарастающий фронт температуры, темпера­

1 4 5

тура газа ниже температуры стенки, и из-за разницы темпе­ ратур происходит разогрев газа. В выходной части, где имеется падающий фронт температуры, температура газа вы-

Рис. 4.4. Перераспределение температур в твэле после пуска теплоносителя.

ше температуры стенки, и стенка твэла прогревается от газа. Постепенно передняя часть твэла охлаждается, хвостовая прогревается и максимум температуры смещается к выход­ ной части. Описанный процесс иллюстрируется рис. 4.4,

146