Алферов Исследователский ядерный 2012
.pdf
Величина начальной ступеньки по приведенной скорости равна
a(t0 ) = |
β a↑ |
|
β + S / n0 . |
(А2.31) |
Важной особенностью процесса является то, что в момент достижения скорости a↑ реактивность отрицательна.
А2.5 Увеличение мощности в энергетическом диапазоне
Основная особенность управления мощностью ИЯР в энергетическом диапазоне – существенное влияние обратной связи по реактивности на вид процессов увеличения физической мощности. Эта обратная связь, приведенная к физической мощности, в рамках точечной модели кинетики описывается дифференциальным уравнением (А2.15).
Полная реактивность системы получается суммированием реактивностей управления и обратной связи.
Подстановка в (А2.15) требуемого закона изменения n(t) (А2.16, А2.17) дает возможность получить зависимость реактивности обратной связи от времени.
Возможный вид переходных процессов в энергетическом диапазоне для различных значений τρ показан на рисунках А2.5, А2.6 (величина kρ предполагается отрицательной).
τρ << (tк-t0) |
τρ ≈ (tк-t0) |
|
Рис. А2.5 Рис. А2.6
101
Основной результат проведенного рассмотрения состоит в том, что для реакторов с отрицательной внутренней обратной связью по реактивности стабилизация на достигнутом уровне мощности nк происходит мгновенно.
А2.6 Снижение мощности
Для режима снижения мощности с постоянной приведённой скоростью справедливы соотношения (А2.21) и (А2.26), в которых величину a↑ следует заменить на (- a↓).
Отличительной особенностью процесса снижения мощности является наличие в показателях экспонент разностей (λi - a↓),
которые могут быть отрицательными при a↓ > λi , что означает расходимость соответствующих экспонент. Графики этих процессов для различных значений a↓ показаны на рисунках А2.7, А2.8.
Рис. А2.7 |
Рис. А2.8 |
102
А2.7 Линейное изменение реактивности
Случай линейного изменения реактивности представляет особый интерес при анализе динамики ядерного реактора. Это обусловлено тем, что такое изменение реактивности имеет место при пусках реактора и при отказах в целях управления двигателями приводов СУЗ.
В общем виде уравнения кинетики для случая ρ(t)=ρ0+γ·t не разрешимы.
Характер процесса может быть оценен из рассмотрения модели мгновенного скачка с одной группой запаздывающих нейтронов.
Полагая в уравнении (А2.13) ρос = 0 и подставляя выражения для С(t) и dC(t)/dt в уравнение (А2.15), получаем
d n(t) |
|
n(t) dρ(t) |
λS |
. |
(А2.32) |
||||
|
|
|
|
|
|
|
|||
dt |
= β − ρ ( dt +λρ(t)) + |
β − ρ |
|||||||
|
|
||||||||
При стационарном начальном состоянии выполняется условие
S = −n0 ρ0 = n0 | ρ0 | . |
|
|
|
|
|
|
(А2.33) |
|||||
С учетом этого из (А2.28) следует |
|
n0 λ| ρ0 | |
|
|||||||||
1 |
|
d n(t) |
1 |
|
dρ(t) |
|
|
. (А2.34) |
||||
a(t) = |
|
|
dt |
= |
|
( |
dt |
+ |
λρ(t)) + |
|
β − ρ |
|
n(t) |
β −ρ |
n(t) |
||||||||||
Уравнение (А2.30) дает возможность оценить сверху значение приведенной скорости a(t) в момент достижения критического состояния, т.е. при ρ = 0
|
|
|
γ |
|
λ| ρ0 |
| |
. |
(А2.35) |
|
|
|||||||
a |
|
ρ=0 ≈ β +λρ(t)) + |
10β |
|
||||
|
|
|
|
|||||
В [1] приведены менее грубые оценки мощности и приведенной скорости в момент достижения критичности
n |
|
|
ρ=0 |
≈ n0 | ρ0 | πλ / 2γβ |
, |
(А2.36) |
||
|
||||||||
|
||||||||
a |
|
|
|
≈ |
γ |
+ 2λγ / πβ |
. |
(А2.37) |
|
|
|
||||||
|
|
ρ=0 |
β |
|||||
|
|
|
|
|
|
|
||
Следует отметить, что при отсутствии источника (т.е. S = 0) уравнение (А2.32) разрешается относительно приведенной скорости
|
1 d n(t) |
1 |
|
dρ(t) |
. |
(А2.38) |
|||||
|
|
|
|
|
|
|
|
|
|||
a(t) = n(t) dt |
= β −ρ ( |
dt +λρ(t)) |
|||||||||
|
|
||||||||||
Отсюда легко получается закон изменения мощности при линейном изменении реактивности
103
|
β |
)(1+ |
λβ |
) |
|
|
n(t) = n0e−λt ( |
.γ |
(А2.39) |
||||
|
|
|||||
|
β −γt |
|
|
|
||
График, описываемый |
|
|
|
|||
выражением (А2.35) и |
|
|
|
|||
точное решение приведены на рисунке А2.9.
Рис. А2.9
А2.8 Ступенчатые приращения скорости изменения реактивности
Полезно оценить характер изменения мощности при ступенчатых
приращениях γ, т.е. скорости изменения реактивности. Учитывая, что
S = 0, и полагая
ρос = 0, упростим структурную модель
рисунка А2.2. Новая структура показана на рисунке А2.10.
Рис. А2.10
Реакция интегратора с обратной связью (β+la-ρ) на ступеньку dρ(t)/dt=Δγ представляет собой экспоненту (аналогично так называемому апериодическому звену)
104
a(t) = |
γ |
(1−e |
− |
β−ρ+a l |
) . |
(А2.40) |
β −ρ +a l |
|
l |
Таким образом, с постоянной времени τ = l/(β -ρ+a·l) устанав-
ливается значение |
a = |
γ |
. |
|
β − ρ +a l |
|
При а < 0,1 с-1 (периоды больше 10 с) и l < 10-3 выполняется
соотношение |
|
(β -ρ) ≈ 30a·l . |
(А2.41) |
Поэтому в (А2.40) можно пренебречь членом a·l |
(структуры на |
рисунке А2.11). |
|
Приращения приве- |
|
денной скорости |
|
a = Δγ /(β - ρ) будут |
|
устанавливаться за |
|
времена порядка |
|
t = l /(β - ρ) ≤ 0,2 c. |
|
Рис. А2.11
С точки зрения задач управления ИЯР это означает практически мгновенную реакцию
приведенной скорости a(t) на конечное изменение скорости введения реактивности.
105
А3 АНАЛИЗ И СИНТЕЗ СИСТЕМ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ (НАЧАЛЬНЫЕ СВЕДЕНИЯ)
А3.1 Анализ управляемых процессов на фазовой плоскости
Поскольку процессы управления ядерным реактором являются существенно нелинейными, необходимо для их исследования использовать соответствующий аналитический инструмент. Удобство и наглядность выделяют в качестве такового метод фазовой плоскости [3].
В общем случае уравнения, описывающие поведение нелинейной динамической системы, могут быть представлены в виде
y j = Φ j (x1...xm ,u1...us ,t),
dxi = Ψi (x1...xm ,u1...us ,t), dt
где y j ( j =1...n) - выходныевеличины, xi (i =1...m) - переменныесостояния, uk (k =1...s) - управления, t - время.
(А3.1)
(А3.2)
Решение системы уравнений (А3.1, А3.2) при заданных начальных условиях позволяет для каждой координаты xi определить описывающую функцию времени yj(t). Процесс, наблюдаемый в системе, можно изобразить в виде линии в п-мерном пространстве состояний.
В том случае, когда управляющие воздействия u1…us не зависят от времени, а зависят только от значений переменных состояния (или выходных величин), систему (А3.1, А3.2) можно, исключая время, свести к (n-1) уравнению. Это дает возможность рассматривать процесс, протекающий во времени неограниченно долго, в ограниченной области (п-1)-мерного пространства, т.к. каждая переменная состояния физически ограничена.
106
Пример 1. Рассмотрим движение объекта, описываемого уравнением
d 2 y = 1 u(t), dt2 m
под действием управлений u1(t) = +Uмакс,
u2(t) = − Uмакс.
Переменными состояния в этом случае являются x1 (t) = y(t),
x2 |
(t) = |
dy(t) |
. |
|
|||
|
|
dt |
|
С учётом этого, из (А3.3) получаем dxdt1 = x2 ,
dxdt2 = m1 u(t).
(А3.3)
(А3.4) (А3.5)
(А3.6)
(А3.7)
(А3.8)
(А3.9)
Исключая из (А3.8), (А3.9) переменную dt, приходим к уравнению
x |
2 |
dx |
2 |
= |
u(t) |
dx . |
(А3.10) |
|
|||||||
|
|
|
m |
1 |
|
||
|
|
|
|
|
|
|
В случае, когда управление не зависит от времени и принимает
значение u(t) = = ±Uмакс, непосредственным интегрированием (А3.10) получаем
|
1 |
(x2 |
− x2 ) = dx |
2 |
= |
|
(x − x ). |
(А3.11) |
||
|
|
|
||||||||
2 |
2 |
20 |
|
|
m |
1 |
10 |
|
||
|
|
|
|
|
|
|
|
|||
где x10 и x20 - значения перемен- |
|
|||||||||
ных x1 и x2 в нулевой момент |
|
|||||||||
времени. |
Траектории |
движения |
|
|||||||
на плоскости |
показаны |
на |
|
|||||||
рисунке А3.1. |
|
|
|
|
|
|
|
|||
Рис. А3.1
107
Метод фазовой плоскости оказывается очень полезным для анализа и синтеза оптимальных по быстродействию систем управления. Обусловлено это тем, что для объектов, описываемых линейными дифференциальными уравнениями, и для некоторых других оптимальное управление состоит в последовательном переключении управляющего воздействия с максимального положительного на максимальное по модулю отрицательное.
Процедура синтеза оптимального закона управления для таких систем состоит в следующем [4]:
-определяют цель управления;
-исследуют условия, позволяющие удержатьобъектвобласти цели;
-исследуют движения объекта под воздействием управлений
u(t)=±Uмакс, т.к. эти управления являются элементами оптимальной последовательности;
- определяют геометрические места точек переключения управления, разделяющие фазовую плоскость на области отрицательного и положительного управления, т.е. закон управления, учитывая при этом, что количество интервалов постоянства управления для неколебательных объектов на единицу меньше порядка дифференциального уравнения, описывающего этот объект.
Пример 2: Рассмотрим задачу оптимального по быстродействию перевода объекта примера 1 в точку с координатами x1=0, x2=0.
Область цепи в этом случае совпадает с началом координат. При этом нулевое управление обеспечивает нахождение объекта в области цепи, после её достижения. Поскольку на завершающем этапе движение всегда будет происходить по траекториям, проходящим через начало координат, то линия переключения
представляет собой комбинацию отрезков фазовых траекторий при u(t)=−Uмакс и u(t)=+Uмакс, проходящих
через ноль. Линия переключения и фазовые траектории переходных процессов показаны на рисунке А3.2.
Рис. А3.2
108
Пример 3: Необходимо синтезировать оптимальный по быстродействию закон управления выходной координатой объекта, описываемого уравнениями
dy(t) |
|
= ( |
dρ(t) |
+λ ρ(t)), |
(А3.12) |
|
dt |
dt |
|||||
|
|
|
||||
dρ(t) |
= u(t). |
|
(А3.13) |
|||
dt |
|
|||||
|
|
|||||
Отметим, что (А3.12) аналогично линейному одногрупповому уравнению кинетики ядерного реактора.
Запишем для полноты ограничение области допустимых управлений
−Uмакс ≤ u(t) ≤ +Uмакс. |
(А3.14) |
Обозначим x1(t) = y(t) и x2(t) = ρ(t).
Поскольку система линейна, решение должно получиться инвариантным относительно начального и конечного значения Поэтому в качестве цели управления рассмотрим состояние
y(t) = 0 и dy(t)/dt = 0. (А3.15)
Из уравнения (А3.12) следует, что в области цели (А3.15)
|
dx2 (t) |
+λ x2 (t) = 0 |
(А3.16) |
|||
|
|
|||||
|
dt |
|
||||
Поэтому в области цели |
|
|||||
−λ x2 |
(t) = |
dx2 (t) |
= u(t), |
(А3.17) |
||
|
||||||
|
|
|
|
dt |
|
|
и, следовательно |
|
|||||
−Uмакс |
≤ x2 (t) ≤ +Uмакс . |
(А3.18) |
||||
λ |
|
|
|
λ |
|
|
Итак, областью цели на фазовой плоскости x1, x2 является отрезок оси x2, удовлетворяющий (А3.18). Управление, удерживающее объект в области цели, получается из (А3.16)
u(t) = −x2 (tк )λe−λ(t−tk ) , |
(А3.19) |
где -tк - момент достижения области цели. |
|
Траектории движения на фазовой плоскости получаются из (А3.12), (А3.13) исключением времени
109
|
dx1 |
=1+ |
|
λ |
|
x |
. |
|
|
|
|
|
|
|
|
|
|
|
(А3.20) |
|||||||
|
|
u(t) |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
dx2 |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Отсюда при u(t) = ±uмакс |
|
|
|
|
|
|
|
|
||||||||||||||||||
(x − x |
) = x |
|
− x |
|
+ |
|
λ |
(x2 |
− x2 |
). |
|
|
(А3.21) |
|||||||||||||
2 |
20 |
|
|
|||||||||||||||||||||||
1 |
10 |
|
|
|
|
|
|
|
|
|
|
|
2u |
2 |
20 |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Семейства траекторий (А3.21) показаны на рисунке А3.3а. |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
−Uмакс |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u(t)=+Uмакс |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u(t)=−Uмакс |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
линия |
|
|
|
|
|
|
|
+Uмакс |
|
|
|
|
|
|
|
|
|
|
|
|
|
переключения |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
управления |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. А3.3 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Вершины парабол (А3.21) лежат на прямых |
|
|
||||||||||||||||||||||||
|
x2 |
= ± |
uмакс |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(А3.22) |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
λ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Линия переключения управления образуется в данном случае комбинацией отрезков фазовых траекторий и области цели.
Линия переключения и фазовые траектории переходных процессов показаны на рисунке А3.3б.
Отличительным свойством рассматриваемого объекта управления является то, что существует область начальных значений (при x2 =0), из которых объект переводится в область цели управлением одного знака. Эта область описывается соотношением
y |
|
≤ |
uмакс |
. |
(А3.23) |
|
|||||
|
|
||||
|
|
|
2λ |
|
|
|
|
|
|||
А3.2 Уравнения систем в изображениях по Лапласу
Во многих случаях поведение систем управления исследовательским реактором исследуется в предположении малых отклонений от некоторых средних значений. В этих случаях системы часто могут
110
