Министерство образования Российской Федерации
Московский государственный технический университет имени Н.Э. Баумана
Козлов О.С.
ДИСЦИПЛИНА: УПРАВЛЕНИЕ В ТЕХНИЧЕСКИХ СИСТЕМАХ
ЛАБОРАТОРНАЯ РАБОТА № 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 |
||
× 10 3 , с |
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 дополнительные точки внутри
указанного диапазона).