
книги из ГПНТБ / Динамика и управление ядерным ракетным двигателем [Текст] 1974. - 253 с
.pdfвекторными, либо — для звеньев с распределенными пара метрами — функциями пространственных координат.
Реакцию координаты х на воздействие / будем называть устойчивой, если для любого е > 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