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

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

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

векторными, либо — для звеньев с распределенными пара­ метрами — функциями пространственных координат.

Реакцию координаты х на воздействие / будем называть устойчивой, если для любого е > 0 существует такое б > О, что для любой функции f (z F, удовлетворяющей условию

 

sup (I/( / ) ! < б,

(4.2)

 

t > о

 

реакция координаты х удовлетворяет условию

 

sup IX (О I < е.

(4.3)

 

t > о

 

(Здесь у

II — норма в пространстве значений, принимае­

мых соответственно х или f; F — класс

функций времени,

для которого определена реакция х;

для / 6 F f (t) = О

при t < 0.)

Если функции х и f принимают скалярные зна­

чения, то F — обычно класс измеримых ограниченных функ­ ций. Если, кроме того, переходная функция от х к f разла­ гается на скачкообразную и абсолютно непрерывную сос­ тавляющие, то необходимым и достаточным условием устой­ чивости реакции х на / является ограниченность вариации переходной функции на (0, оо) [5].

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

для

любого е > 0 существует такое б > 0,

что для лю­

бого

у0 6 У0, удовлетворяющего условию

 

 

Ы | < S ,

(4.4)

реакция х на у0 удовлетворяет (4.3). (Здесь Y0—множество возможных значений у0.)

Для устойчивости реакции х на у0 необходимо и доста­ точно ограниченности реакции х на y0d S yo, равномерной по t и по Sy0 — единичной сфере пространства Е0 [6].

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

117

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

Рис. 4.1. Переходные процессы при возмущении по реактив­ ности.

Д а н н ы е о б у с т о й ч и в о с т и и х а р а к т е р е д и н а м и к и Я Р Д .

Ядерный ракетный двигатель с твердой активной зоной, как показывают аналоговое моделирование и эксперименты [7, 8], обычно оказывается устойчивым. В качестве иллюст­ рации на рис. 4.1, 4.2 приведены осциллограммы переход­ ных процессов для одной из конструкций ЯРД с твердой активной зоной и твердым замедлителем.

На осциллограммах даны процессы, вызванные скачко­ образным изменением реактивности бк (рис. 4.1) и отно­ сительного отклонения коэффициента гидравлического соп­ ротивления основного дросселя рабочего тела р (рис. 4.2). Рассмотрение осциллограмм показывает, что двигатель ус­ тойчив и обладает относительно малой инерционностью: процессы по основным координатам ѵ, -й, я, / (относитель­ ные отклонения соответственно плотности нейтронов в активной зоне реактора и температуры, давления, расхо­ да в камере сопла) заканчиваются, например при возмуще­ нии по bk, в основном за 0,7 сек, имея некоторую колеба­

118

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

Возмущение м=ІО°/о

Статина

0,31 % (J )

г

І2 % П

тг

л

 

 

 

.

н

\

г

 

„-------------------------------

 

0,6 °/о(ж)

,

< r ö J%

;

*

і--------------

І35 °/о N

U

 

У

 

2 і,сек

Рис. 4.2. Переходные процессы при возмущении по коэффициенту гидравлического сопротивления.

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

§ 12. Управление ЯРД на номинальном режиме

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

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

119

мый арсенал средств воздействия на объект (управляющих воздействий) и объем доступной для использования инфор­ мации о текущем состоянии объекта (измеряемые коорди­ наты).

Более традиционным методом синтеза системы управле­ ния является аналоговое моделирование системы объект— регулятор; при этом модель объекта сопрягают с моделью регулятора заранее выбранной структуры, учитывающей располагаемый набор управляющих воздействий и измеряе­ мых координат, и подбором параметров регулятора ста­ раются удовлетворить техническим требованиям к системе. Метод моделирования достаточно оперативен и нагляден, но, как всякий метод подбора, при значительном количест­ ве выбираемых параметров не может дать окончательного ответа на вопрос о разрешимости поставленной задачи син­ теза, описать все множество законов управления, решающих эту задачу, выбрать наилучшее из возможных решений. Поэтому все более широкое распространение получают ана­ литические методы синтеза, дающие возможность (с помощью современного математического аппарата и ЦВМ) получить ответ на эти вопросы.

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

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

 

 

Dx = М[і + Ф/,

(4.5)

где X =

(xj,

j — 1 ~ п )

— вектор-столбец координат объек­

та; р =

(рь

k = 1-рт) — вектор-столбец координат регу­

лирующих

органов;

/ = ( / , ,

I = І-т-s) — вектор-столбец

внешних

воздействий (возмущающих или задающих); Xj,

Рь, fi — скалярные величины; D, М,

Ф — матрицы с раз­

мерами соответственно п Х п ,

п Х т ,

n X s и элементами

из поля

W рациональных дробей от оператора дифферен­

цирования р.

1 20

Уравнение (4.5) задает математическую модель объекта вблизи номинального режима в виде п-мерной линейной сис­ темы с сосредоточенными постоянными во времени па­ раметрами, подверженной действию т регулирующих орга­ нов и s внешних воздействий. Вектор / составлен лишь из тех внешних воздействий, которые будут учитываться при формулировании требований к поведению замкнутой системы; вектор р составлен из тех регулирующих орга­ нов, которые предполагается использовать для подачи на объект управляющих воздействий. В соответствии с физи­ ческими представлениями об определенности поведения объекта при заданных воздействиях и начальных условиях считаем матрицу D невырожденной.

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

J

У = 2 V j { p ) Z j ,

( 4 .6 )

/= 1

 

где у — выход звена; Zj — входы звена; Vj (р) — переда­ точная функция от zj к у (Vj — элемент поля W дробно­ рациональных функций оператора дифференцирования р). В соответствии с этим искомая система уравнений регуля­ тора должна иметь вид (также для нулевых начальных ус­ ловий)

УѴр,==/(хи + І / и + йР;

;

A^ß ß - К$ хи -)-Lß/п +

(X, j

где ß = (ßr, r = 1 — b) — вектор-столбец внутренних координат регулятора (т. е. тех выходов его звеньев, которые

не

подаются на

объект);

хи =

(xj,

j =

nk +

1-няд +

+

пи) (0 sj; nk

 

nh + пя) — вектор-столбец

измеряемых

координат объекта (т. е. тех

координат,

которые

могут

служить

входами

звеньев

регулятора);/и =

{f и

 

I =

1-і-

-t-s„) (0 ^

s„ ^

s) — вектор-столбец

измеряемых

воздейст­

вий; N , К, L,

В — матрицы

с

элементами

из

поля

W;

jVß, Kß,

Lp,

Afß — матрицы

с

элементами

из

поля

W;

элементы главной диагонали у матриц N и N$ должны быть

равны 1 [k-e уравнение 1-й строки записи (4.7)

соответству­

ет звену с выходом [хЛ, г-е уравнение 2-й строки

(4.7) —

звену с выходом ßrl.

 

 

 

 

 

 

 

 

 

121

Вид уравнений (4.7) учитывает ограниченность информа­ ции о состоянии объекта и о воздействиях на него, которая может быть использована при формировании закона управ­ ления. Кроме того, вид уравнений (4.7) показывает, что иные воздействия, кроме /д (/ = 1 s„), которые могут дей­ ствовать на регулятор, здесь не учитываются. (Наличие по­ мех в регуляторе или неточности измерительных устройств во многих случаях могут быть учтены введением соответ­ ствующих звеньев в объект.)

Так как уравнения (4.5) и (4.7) линейны относительно

fi, их неизвестные Xj,

[ih, ßr можно считать линейными фор­

мами ОТ fi'.

 

 

 

 

 

 

 

S

 

 

 

 

 

Х ] =

I= 1Xjifi,

/ = 1 4 - л ;

 

 

s

 

 

 

 

 

 

V

 

 

k = 1

-pm\

(4.8)

 

1= 1М-ы fu

 

 

 

 

І

 

ßrlf„

1

+ b,

 

 

i=1

 

 

 

 

где Xji, \ihi, ßr; — элементы поля W — характеризуют реак­ ции соответствующих координат системы объект—регулятор

на

воздействие / г.

 

 

 

 

 

 

Требования к

поведению замкнутой системы сформули­

руем следующим образом:

 

 

 

 

а) система уравнений ((4.5), (4.7)) должна иметь единст­

венное решение

(х0,

р0, ß0) = £0 относительно вектора

(х,

р,

ß)

= £ линейных

форм от / г;

 

 

 

б) решение £п =

(х0,

р0, ß0) должно удовлетворять соот­

ношениям

 

P x ^ - S ^ + Qf,

 

(4.9)

 

 

 

 

 

 

где

Хр =

(xj, і ==-

1

-р л р) (пк < пр < пк -I- пи) — вектор-

столбец

первых

 

Пр

координат

объекта;

р р =

(рк,

k =

1

-р пгр) (0 г$; тѵ ^

т ); Р, Q,

S — матрицы

с ѵ

стро­

ками

и заданными элементами из поля W.

 

 

 

Выбор

векторов хр,

р р и матриц Р, Q, S должен опре­

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

мента Qu (і —номер строки;

I—номер столбца) матрицы

Q

означает

требование

инвариантности

величины

пр

>

 

 

 

ZPijXj - Z S ih[ih по отношению к воздействию / г.

k = i

k — I

 

 

 

122

Выделив из системы уравнений ((4.5), (4.9)) максималь­ ную линейно независимую подсистему, содержащую (4.5), и соответственно уменьшив ѵ, можно считать далее систему ((4.5), (4.9)) независимой, т. е. имеющей ранг ѵ + п. (Если все соотношения (4.9) зависят от уравнений объекта (4.5), то это означает, что они выполняются при любом виде урав­ нений регулятора, и тогда задача синтеза регулятора не возникает.) Точно так же, исключив из рассмотрения коор­ динаты регулирующих органов и воздействия, которые не входят в систему ((4.5), (4.9)), можем считать матрицу этой системы не имеющей нулевых столбцов.

Отметим непосредственно вытекающие отсюда необ­ ходимые условия разрешимости задачи синтеза: число уравнений в каждой из систем (4.9) и ((4.5), (4.9)) не должно превышать числа входящих в них неизвестных, т. е.

V < пр+ шр;

(4.10)

Первый способ решения задачи синтеза. Решение по­ ставленной задачи синтеза можно было бы проводить пу­ тем составления и исследования уравнений, которым должны удовлетворять неизвестные коэффициенты систе­ мы (4.7) [при выбранных в соответствии с (4.10) множествах входов и выходов звеньев регулятора]. Из требований к поведению замкнутой системы следует, что матрица системы ((4.5), (4.7)) должна быть невырожденной; при выполнении этого условия решение системы имеет вид

l =

Ff,

 

(4.11)

где F — матрица (п + т +

b)

X s с элементами, завися-

 

 

 

 

( х \

щими от коэффициентов системы

(4.11); £ — столбец ц .

Представим матрицу F в блочной записи [9]:

Р/

 

ч

 

\

 

 

F в

 

 

F =

 

 

,*н

(4.12)

 

 

Рр

F w

 

 

Пцн

 

 

 

Fp

/

р

 

(рядом с блоками матрицы F стоят обозначения соответст­ вующих им столбцов — блоков столбца £; хв = (xj,

123

j — я р 'I' 1

як “Ь яп);

лп

{Xj,

j

tiK I

яи

h l .

-7- n); p„ -

(|xK, k =

mp + 1 — tn).

 

 

 

 

Тогда условия

(4.9)

эквивалентны соотношениям

 

 

 

 

 

PFp- S F w ---Q,

 

 

 

(4.13)

задающим

vs

уравнений

относительно

+

пя + s„ +

+ b — 1) (т +

Ь)

неизвестных

коэффициентов

системы

(4.7). Эти уравнения, вообще говоря, нелинейны; их иссле­ дование и решение при значительных размерах матриц систем (4.5) и (4.9) может оказаться затруднительным. Поэ­ тому ниже излагается иной способ построения системы уравнений регулятора, являющийся, по сути дела, обобще­ нием процедуры исключения неизвестных из системы урав­ нений.

Второй способ решения задачи синтеза. Предварительно представим расширенные матрицы систем (4.5), (4.7) и (4.9) в блочной записи:

х к

'Gip

х в

х п

Вр

Мчі

ß

 

 

п=(D „

Ц ,р DB DH -7W p — М я О - Ф и — Фн);

 

/ О - Т С н р - ^ з 0

Np

Nn

- B - L

 

0\

(4.14)

ІО - X ß„ p - X

ßB 0

— Л1Рр

- A f p„

ІѴр - L p

О /’

 

G = ( P K

Рир

О

О

- S

0

 

0 —<2И—-Q„) .

 

(над блоками матриц П, Р, G стоят обозначения соответ­

ствующих им столбцов — блоков столбца ср — /£\;

ч)

 

=

ixj>

І

^ “ ^к)>

-Чір ,(xj> /

я к ~Ь 1

Яр).

 

 

Обозначим Х0, Х р,

Хж множества решений систем урав­

нений соответственно (4.5), (4.7), (4.9) относительно век­ тора £ линейных форм от / г. Тогда сформулированные тре­ бования к замкнутой системе означают, что множество

* 0 П Х р должно состоять из единственного вектора | 0

и этот

вектор должен принадлежать Хж.

 

Пусть существует матрица R, удовлетворяющая постав­

ленным требованиям. Отсюда следует:

 

а)

система

уравнений (4.9) — следствие системы ((4.5),

(4.7)).

Поэтому

[9]

существуют такие матрицы Сд

v X п

и Ср V X +

Ь),

что

 

 

 

 

0 = СоП-(-Ср /?.

(4.15)

124

Из (4.15) следует соотношение, которому должна удов­ летворять матрица С0:

 

 

 

 

 

О„ = С0Пн>

 

 

(4.16)

 

 

 

пн= (Д< Дг —фн);

 

 

 

 

 

Д=(ЯК о

QH).

 

 

а также

включение

Х ѵ с і с,

где обозначено;

Х с— мно­

жество решений системы

(GH—С0Пи)сри = 0,

 

 

 

 

•Чір

 

 

 

 

Ин

ß

 

 

Пи = (Д р Д

- м р - М и 0 и);

 

0И= (ЯИр

0

—S

 

0

0

-Q „);

 

 

 

 

 

 

/

А:ир \

 

 

 

 

 

 

 

 

 

 

 

\

 

 

 

 

 

 

 

Фи

 

Ир

 

 

 

 

 

 

 

 

И-н

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ß

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно,

существует

такая

матрица Л v X (м +

+ Ь),

что

 

 

G „ - C 0n„ = A/?;

 

 

(4.17)

 

 

 

 

 

 

б)

система

уравнений

 

 

 

 

 

 

 

 

 

 

 

П cp = 0;

 

 

(4.18)

 

 

 

 

(Ои- С 0Пи)фи = 0

 

 

 

 

 

 

 

 

 

эквивалентна [с учетом (4.16)]

системе ((4.5), (4.9)), поэто­

му ранг матрицы G„ — С0ПИдолжен быть равен

ѵ.

Из (4.17) вытекает [9], что матрица Л должна иметь ранг

V . Число ее строк ѵ не превышает числа столбцов т + Ь\

дополним

ее

пя = т + Ь — ѵ строками

до

квадратной

невырожденной матрицы

Л0.

Получим соотношение

 

 

 

 

/О и -С о П

 

\

 

 

(4.19)

 

 

 

 

[

 

R ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

*^"ир

 

Ир

Ин

ß

 

 

 

Яд ~ (Дир

Д в

ЛТдр Мяп Бд —Фди) —матрица ПдХ

 

 

 

 

X (л„ + т + b -\г s„).

 

 

125

Так как матрица Л0 невырождена и блок

И"р

!-1п

ß

Nv

NH

- В \

— M ß p

і Щ п

Л /'р /

матрицы R имеет (по условию) диагональные элементы, равные единице, то соответствующий блок

1*р

 

ß

S i- С0 Мр с 0[м„

0

Яідр

Яідн

Л0 R0

Яд

матрицы

/ G„ —С,)ПИ\

1

 

Яд /

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

 

(Ga— C0 Пи)фи =

0;

I

 

 

 

(4.20)

 

Яд Ф„ = 0,

 

1

 

 

 

 

 

 

 

 

 

где матрица С0 удовлетворяет (4.16),

матрица R^

имеет

Пд = т +

b — V строк

и матрица

R '0

не имеет нулевых

столбцов.

 

 

 

 

 

 

 

 

 

 

Обратно, пусть существует матрица С0, удовлетворяю­

щая (4.16),

и пусть система уравнений

 

 

 

 

 

 

 

Пе Фе

:0;

 

 

 

 

 

(4.21)

 

 

Ge Фе

0,

 

 

 

 

 

где

 

 

 

 

 

 

 

■*ир

 

f^p

 

 

 

 

 

 

 

 

 

h.

 

 

 

 

 

 

Агр А Du И р - М и - Ф и - ф „ );

 

Ge- ( P K ЯИр 0 0 - S

 

0

 

 

Qu - А ) ;

 

 

 

Фе

 

 

 

 

 

 

 

 

совместна. Тогда для любого

b

0

совместна и

система

((4.5), (4.9)).

b >

 

 

 

 

 

 

 

 

Выбрав

некоторое

0,

дополним

систему

 

((4.5),

(4.9)) «д =

т + b — V

уравнениями

до

системы

с

квад­

126