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

Алферов Исследователский ядерный 2012

.pdf
Скачиваний:
28
Добавлен:
12.11.2022
Размер:
3.8 Mб
Скачать

Величина начальной ступеньки по приведенной скорости равна

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) =

γ

(1e

βρ+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

u(t)=Uмакс
u(t)=+Uмакс
линия
переключения управления

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

Процедура синтеза оптимального закона управления для таких систем состоит в следующем [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λ(ttk ) ,

(А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