Министерство образования Российской Федерации
Московский государственный технический университет имени Н.Э. Баумана
Козлов О.С.
ДИСЦИПЛИНА: УПРАВЛЕНИЕ В ТЕХНИЧЕСКИХ СИСТЕМАХ
ЛАБОРАТОРНАЯ РАБОТА № 5
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ДИНАМИКИ
НЕЛИНЕЙНОЙ САР ЯДЕРНОГО РЕАКТОРА
Москва, 2003 г.
Цель работы:
Типовая система автоматического регулирования ядерного реактора - управляемая
система очень большой размерности. Детальное исследование динамических свойств такой САР с
применением методов численного моделирования - сложная научно-техническая задача, так как
в реакторной установке (РУ) присутствуют технические устройства, для описания динамических
процессов в которых используется информация из большинства фундаментальных и прикладных наук.
До начала 90-х годов прошлого столетия для исследования нестационарных
процессов в сложных управляемых технических системах разрабатывались
специализированные динамические программы (применительно к конкретной установке).
Использование таких программ для исследования нестационарных процессов в случае, например,
значительной модернизации этой же установки требовало серьезной переработки расчетной
программы на уровне исходных кодов (математические модели, алгоритмы и т.п.), что реально
способны были выполнить только программисты, создавшие эту программу.
Значительный прогресс, достигнутый в последнее десятилетие в аппаратных и
программных возможностях современной вычислительной техники, создал необходимую базу
для разработки принципиально новых средств интеллектуального САПР, например,
объектно-ориентированных программных сред для исследования нестационарных процессов в
сложных динамических системах.
В настоящее время в РФ создан ряд современных программных комплексов (ПК) для
исследования нестационарных процессов в реакторных установках.
Так, например, программно-инструментальные комплексы
АИС95 и “ЭНИКОКАД” предназначены для моделирования процессов нейтронной
кинетики, теплогидравлики и автоматического управления применительно к задачам
разработки полномасштабных тренажеров энергоблоков АЭС с реакторами типа РБМК.
Программный комплекс ТЕРМИТ-Д предназначен для проектного обоснования безопасности
ядерных паропроизводящих установок (ЯППУ) транспортных ЯЭУ.
Несмотря на высокий научный уровень вышеуказанных программных средств,
необходимо отметить их два серьезных недостатка:
- во-первых, каждый из них может быть рекомендован для
использования в проектных разработках новых реакторных установок
только аналогичного типа;
- во-вторых, ни один из этих ПКне пригоден для
использования в учебном процессе высшей школы по причине отсутствия соответствующего
методического наполнения.
К программным средствам интеллектуального САПР относится и программный
комплекс “МВТУ”, одним из главных достоинств которого является
инвариантность к предметной области исследуемого объекта или
физического явления. Это позволяет выполнить в среде ПК “МВТУ” численное исследование
рабочих процессов практически в любых сложных технических системах: в электромеханических,
в тепловых, в пневмо- и гидродинамических, в робототехнических, в ряде
других комбинированных динамические систем, в том числе и в реакторных системах.
Для описания динамики нейтронно-физических и теплогидравлических процессов
в ядерных реакторах в ПК “МВТУ” используются более простые, но более
“быстрые” математические модели, чем в вышеупомянутых отраслевых программных средствах.
Программно-технические возможности ПК “МВТУ” качественно превосходят
возможности отечественных отраслевых программных средств в задачах разработки
математических моделей систем автоматического и логического управления, систем защит и
блокировок применительно к проектному обоснованию АСУ ТП для энергоблоков АЭС с
реактором типа ВВЭР и РБМК.
Не менее важным достоинством ПК “МВТУ” является широта области применимости:
от простейших динамических задач учебного назначения до реальных отраслевых
разработок, в том числе и в экспортном исполнении.
Наличие в ПК “МВТУ” подробного учебно-методическое сопровождения
позволяет использовать его в учебном процессе высшей школы по многим инженерным специальностям,
включая и выбранную Вами…
В предыдущих лабораторных работах Вы освоили большинство методов работы в
среде ПК “МВТУ”. В настоящей лабораторной работе Вы освоите еще ряд новых методов
работы в среде ПК “МВТУ”, а также выполните самостоятельное численное исследование
динамических характеристик упрощенной математической модели нелинейной САР ядерного
реактора с регулятором релейного типа.
Перейдем к выполнению заданий настоящей лабораторной работы…
![]()
(к содержанию)
1. МАТЕМАТИЧЕСКИЕ МОДЕЛИ НЕЙТРОННОЙ КИНЕТИКИ
1.1. Описание блоков специализированной библиотеки Кинетика нейтронов
Специализированная библиотека Кинетика
нейтронов ПК “МВТУ” содержит три типовых блока, два из которых
описывают нейтронно-кинетические процессы в ядерном реакторе в точечном односкоростном
приближении, а третий – описывает динамику остаточного энерговыделения с учетом
предыстории реактора (кампании).
Первые два блока позволяют описывать кинетику ядер-предшественников
запаздывающих нейтронов от одно- до n- группового приближений. Математические модели
этих блоков получены на основании известных уравнений кинетики “точечного” ядерного реактора
в односкоростном приближении (т.е. процесс деления ядер осуществляется нейтронами одной
энергетической группы - либо только тепловыми, либо только быстрыми):
|
(1.1) |
|
(1.2) |

Рис. 1.1
Значение времени жизни мгновенных нейтронов (по умолчанию)
соответствует
приблизительно времени жизни в ядерном реакторе типа РБМК.
Значения в последней диалоговой строке (Нормировка) соответствуют
следующим видам выходного сигнала блока: 1 - нормированная мощность
, а 0 - безразмерное отклонение мощности
.
Второй типовой блок соответствует
Модели мгновенного скачка. Входом является изменение реактивности
D
r
(t) = r
(t) - r
(0), а выходом – либо нормированная мощность
, либо безразмерное отклонение
мощности
.
Уравнения кинетики нейтронов после преобразований исходной системы уравнений принимает вид:
|
(1.3) |
|
|
(1.4) |
где
- стационарное значение мощности при t = 0.

Рис. 1.2
Третий типовой блок (Остаточное энерговыделение) в библиотеке Кинетика нейтронов позволяет учесть дополнительный вклад остаточного тепловыделения продуктов деления в тепловую мощность ядерного реактора, что особо актуально при резких снижениях нейтронной мощности. Входной сигнал в блок - относительная нейтронная мощность (нормированная на номинальную нейтронную мощность), а выходной сигнал из блока – относительная тепловая мощность реактора (нормированная на номинальную нейтронную мощность), определяемая выражением:
|
(1.5) |

Рис. 1.3
Данные в диалоговом окне на рис. 1.3 описывают следующую
предысторию работы ядерного реактора: на момент моделирования реактор имеет кампанию 9
×10 6
с, причем: при 0 =< t =< 1×
10 5 c относительная нейтронная мощность реактора
(нормированная на номинальную нейтронную мощность) равнялась 0.5;
при 1×105 < t =< 5
×10 5 c
;
при 5 ×10 5 < t =< 1
×10 6 c
;
при 1×10 6 < t =<
9×10 6 c
.
В прошлом семестре при выполнении лабораторных работ и
домашнего задания по курсу “УТС” Вы использовали простейшую математическую модель
нейтронно-кинетических процессов в ядерном реакторе, а именно: “точечную” модель
с одной эффективной группой запаздывающих нейтронов.
Если рассматриваются “длительные” переходные процессы, то значение
эффективной постоянной распада l
ядер-предшественников запаздывающих нейтронов в одногрупповой модели рекомендуется
получать осреднением “времен жизни” групп запаздывающих нейтронов по соотношению:
, |
(1.6) |
, |
(1.7) |
которое дает значение l
= 0.405 с –1 .
Выполним сравнение частотных и переходных характеристик для одногрупповой
(для обоих вариантов расчета l
) и для “классической” моделей кинетики нейтронов. Сформируйте структурную схему, внешний
вид которой должен быть близок рис. 1.4.

Рис. 1.4
Внимание:
Набранные структурные схемы, математические выкладки и результаты расчетов необходимо представлять преподавателю для “согласования”.
Рис. 1.5
Переместите в окне Редактора… курсор на
кнопку Вычислить константы и выполните щелчок левой клавишей мыши:
появится таблица с расчетными данными. Используя “прокрутку” таблицы, убедитесь,
что рассчитанные значения эффективных постоянных распада ядер-предшественников для
“коротких” и “длительных” переходных процессов совпадают с приведенными выше значениями...
Закройте окно Редактора… и задайте параметры блоков структурной схемы.
В блоке, описывающем 6-ти групповую модель кинетики нейтронов, оставьте параметры диалоговых
строк по умолчанию. В блоках, описывающих одногрупповую модель кинетики нейтронов, в
диалоговой строке Постоянные распада групп… введите Lam
или Lam_1, соответственно. Что вводить в других диалоговых строках этих
блоков – необходимо определить самостоятельно…
Выполним сначала моделирование “длительного” переходного процесса при скачке
реактивности 0.01×bэфф. Установите в
диалоговом окне Параметры расчета следующие параметры интегрирования:
Время интегрирования – 1000; Минимальный шаг интегрирования –
1е-10; Максимальный шаг интегрирования – 0.1; Шаг вывода результатов
– 0.1. Параметры других диалоговых строк - по умолчанию.
Выполните расчет переходного процесса. Приведите вид линий и параметры осей координат в графическом окне к виду, близкому рис. 1.6.

Рис. 1.6
На рис. 1.6 линией двойной толщины показан график переходного
процесса для “классической” модели кинетики нейтронов, сплошной одинарной линией – для одногрупповой
модели с эффективной постоянной распада Lam, рассчитанной по соотношению (1.6),
а пунктирной линией - для одногрупповой модели с эффективной постоянной распада Lam_1,
рассчитанной по соотношению (1.7).
Результаты рис. 1.6 показывают, что с эффективной постоянной распада Lam
одногрупповая модель кинетики лишь приблизительно соответствует “классической” модели кинетики,
а с эффективной постоянной распада Lam_1 – различие огромно
(при t = 1000 c приблизительно в 40 раз).
Выполним моделирование “короткого” переходного процесса при скачке
реактивности 0.1×
b
эфф . Установите в диалоговом окне Параметры расчета следующие
параметры интегрирования: Время интегрирования – 3; Минимальный шаг интегрирования
– 1е-10; Максимальный шаг интегрирования – 0.001;
Шаг вывода результатов – 0.001. Параметры других диалоговых строк - по умолчанию.
Выполните расчет переходного процесса. Приведите вид линий и параметры осей координат в
графическом окне к виду, близкому рис. 1.7.

Рис. 1.7
Продолжая сравнение, выполним расчет частотных характеристик для
“классической” и одногрупповой моделей кинетики нейтронов (для обоих вариантов вычисления
l).
Выполните щелчок “мышью” по командной кнопке Старт и затем по кнопке Стоп:
расчет будет прерван, так и не начавшись; произойдет инициализация структурной схемы при нулевом
сигнале на входе блоков, описывающих кинетику нейтронов...

Рис. 1.8
Переместите курсор на кнопку Годографов, выполните
щелчок “мышью” и затем заполните появившуюся таблицу…
Далее выполните расчет годографов Найквиста для “классической” и обоих вариантов
одногрупповой модели кинетики нейтронов. Скорректируйте параметры линий годографов и осей координат,
придав графическому окну вид, близкий рис. 1.9.

Рис. 1.9
На рис. 1.9 линией двойной толщины показан годограф для
“классической” модели кинетики нейтронов, сплошной одинарной линией – для одногрупповой модели
с эффективной постоянной распада Lam, рассчитанной по соотношению (1.6),
а пунктирной линией - для одногрупповой модели с эффективной постоянной распада
Lam_1, рассчитанной по соотношению (1.7).
Анализ данных рис. 1.9 показывает, что в области высоких частот все 3 годографа почти
совпадают, а в области низких частот - заметно количественное расхождение…
Дополним вышеотмеченные совпадения/расхождения сравнением ЛАХ и ФЧХ для рассматриваемых
математических моделей кинетики нейтронов.
Переместите курсор на кнопку Характеристик,
выполните щелчок “мышью” и далее заполните появившуюся таблицу для расчета ЛАХ.
Выполните расчет ЛАХ для “классической” и обоих вариантов одногрупповой модели кинетики
нейтронов. Скорректируйте параметры линий и осей координат, придав графическому окну вид,
близкий рис. 1.10.
|
|
|
|
(1.8) |

Рис. 1.12

Рис. 1.13
Установите в диалоговом окне Параметры расчета следующие
параметры: Время интегрирования – 200; Минимальный шаг интегрирования–
1е-12; Максимальный шаг интегрирования – 0.1; Шаг вывода результатов
– 0.1. Другие параметры - по умолчанию.
Выполните расчет, скорректируйте параметры линий и осей координат, придав графическому
окну вид, близкий рис. 1.14.
На рис. 1.14 сплошной линией двойной толщины представлено поведение относительной
тепловой мощности, а пунктирной – поведение относительной нейтронной мощности. Данные рис. 1.14
показывают, что примерно к 170 секунде относительная нейтронная мощность достигла уровня 1-го
процента от номинала, в то время как относительная тепловая мощность – более 3 %.
Рис. 1.14 |
Рис.
1.15 |
Структурная схема САР ядерного реактора
Рис. 2.1
Структурная схема Субмодели Привод СУЗ

Рис. 2.2
В настоящей лабораторной работе субструктура Логика аварийной защиты имеет только Входной и Выходной порты (соединенные между собой). Данная Субмодель “отключена” посредством задания в Главном Схемном Окне в блоке Сумматор нулевого весового коэффициента по соответствующему входу (см. рис. 2.1). Набор структурной схемы Логики аварийной защиты будет выполнен в последней (№ 7) лабораторной работе, где данная Субмодель будет “отслеживать” нейтронную мощность, период разгона ядерного реактора и температуру теплоносителя на выходе из активной зоны с выдачей логического сигнала на сброс всех стержней аварийной защиты (режим АЗ-5) при превышении заданных уставок.
![]() |
(2.1) |
Блок Перемещение регулирующего стержня (см. рис. 2.2) описывается следующим динамическим уравнением:
| t пр × x’’(t) + x’(t) = kпр × u(t) , |
(2.2) |
| r ст (х) = 0.5× r ст* × sin(a× x + b) |
(2.3) |
2.3. Задание параметров САР через механизм глобальных параметров
Некоторые параметры математических моделей динамики блоков на рис. 2.1 - 2.2
(включая и ряд параметров блока Температурная обратная связь) необходимо задать в окне
Редактора глобальных параметров Проекта (Субмодели).
Наибольший эффект от использования механизма глобальных параметров будет в тех
случаях, когда один и тот же параметр (константа) используются в качестве параметра
(вводимого в диалоговом окне) во многих блоках структурной схемы. В этом случае при изменении величины
этого параметра можно скорректировать его значение только в окне Редактор глобальных
параметров…, а не корректировать его значение во всех блоках, где он используется.
На рис. 2.3 представлена копия окна Редактор глобальных параметров… с введенными
значениями и комментариями к ним. По аналогии с рис. 2.3 введите в окно Редактор… свои исходные
данные (см. табл. 1 в разделе 3)…

Рис. 2.3
Откройте диалоговое окно блока Кинетика нейтронов и в 1-й диалоговой строке вместо числа 6.5е-3 введите beff. Закройте это диалоговое окно (кнопка Да). Откройте диалоговое окно блока с подписью 1/beff (см. рис. 2.1) и введите в диалоговой строке вместо числа 1 выражение 1/beff. Параметр pmax Вам необходимо будет использовать при задании параметров блока Характеристика регулирующего стержня.
Учитывая, что невозможно сформировать абсолютно полную библиотеку
типовых блоков, в ПК “МВТУ” разработаны средства, которые позволяют Пользователю расширить состав
личной библиотеки за счет создания новых типов блоков, например, посредством встроенного
Интерпретатора математических функций, на базе которого функционирует и “Новый” блок,
позволяющий прямо в процессе работы создавать экземпляры блоков со своими оригинальными математическими
моделями.
Использование “Нового” блока в качестве функционального весьма эффективно в
случае наличия в модели сложных функциональных преобразований, когда использование для этих
целей элементарных функциональных типовых блоков приведет к неоправданному усложнению
структурной схемы.
На рис. 2.4 представлена экранная копия окна Редактор интерпретатора математических
функций, где в текстовом виде задана математическая модель динамики блока Температурная обратная
связь.

Рис. 2.4
В табл. 1 приведены исходные данные по параметрам элементов САР ЯР, которую Вам предстоит исследовать. Как видно из табл. 1, параметры САР Вашего варианта не совпадают с данными, которые приведены выше в разделе 2 настоящих указаний…
Таблица 1
|
№ |
Элемент САР |
Параметры элементов САР |
Номер варианта |
||
|
1 |
2 |
3 |
|||
|
1 |
Задатчик мощности |
Время , с |
10 |
10 |
10 |
|
Y0 |
0 |
0 |
0 |
||
|
Y1 |
0,1 |
0,1 |
0,1 |
||
|
2 |
Привод СУЗ |
rст*/b эфф |
0,6 |
0,6 |
0,6 |
|
t пр , с |
0,2 |
0,25 |
0,25 |
||
|
Тхода , с |
5…50 |
5…50 |
5…50 |
||
|
3 |
Управляющее реле |
b |
0,02…0,005 |
0,02…0,005 |
0,02…0,005 |
|
m |
0,4 |
0,6 |
0,8 |
||
|
4 |
Ядерный реактор |
Vтопл , м 3 |
0,1 |
0,2 |
0,4 |
|
g топл , кг/м 3 |
10000 |
9000 |
8000 |
||
|
Стопл , Дж/кг× К |
300 |
350 |
400 |
||
|
|
1 |
0.1 |
0.05 |
||
|
b эфф × 10 3 |
6 |
6,5 |
7 |
||
|
N0 , МВт |
10 |
20 |
50 |
||
|
5 |
Температурная обратная связь |
Т0 , К |
700 |
750 |
800 |
|
Tw , K |
500 |
550 |
550 |
||
|
a × 10 4 , 1/K |
0,7…1,5 |
0,7…1,5 |
0,7…1,5 |
||
|
6 |
Возмущающее воздействие |
D t возм , с |
2…20 |
2…20 |
2…20 |
|
D r возм /b эфф |
0,1…0,3 |
0,1…0,3 |
0,1…0,3 |
||
|
7 |
Запаздывание по каналу измерения |
t зап , с |
0,2…1 |
0,2…1 |
0,2…1 |
Большая часть параметров блоков в структурной схеме
САР не требует особых пояснений и комментариев. Ряд параметров в табл. 1 необходимо прокомментировать.
В субмодели Привод СУЗ параметром Тхода обозначено время перемещения
регулирующим стержнем СУЗ всей активной зоны (сверху-вниз или наоборот). С помощью этого параметра
Вы определите коэффициент скоростной эффективности привода. На рис. 2.2 передаточная функция блока
Перемещение регулирующего стержня неизвестна: ее требуется определить на основании уравнений
динамики…
В блоке Возмущающее воздействие параметр D
t возм задает время, за которое величина возмущающего воздействия линейно
изменяется от нуля до D
r
возм .
Данные в табл. 1 типа 5…50 подразумевают, что Вы должны выполнить какое-то исследование
при варьировании соответствующего параметра в указанном диапазоне (обычно еще 2…3 дополнительные точки внутри
указанного диапазона).