Анализ динамики линейных и линеаризованных САР, описываемых в переменных «Вход – Выход»

Основные этапы работы в SimInTech

В данном разделе рассматриваются основные процедуры работы, освоение которых является обязательным условием для самостоятельной работы в SimInTech.

Формирование структурной схемы и ее параметров, выбор метода, параметров интегрирования и т.п. целесообразно проводить в следующей последовательности:

• Используя палитру блоков, сформировать схему модели в окне проекта с использованием необходимых блоков.
• Соединить блоки линиями связи.
• Задать свойства блоков на структурной схеме.
• Задать конечное время интегрирования, выбрать необходимый метод интегрирования и настроить другие параметры расчета.
• Сохранить проект.
• Запустить проект на расчёт. Проанализировать результаты моделирования.

Динамика САР ядерного реактора.

Для ознакомления с функциональными возможностями SimInTech рассматривается демонстрационный пример «Релейный регулятор для САР ядерного реактора».

Необходимо запустить SimInTech. Нажатием на кнопку «Открыть» вызвать окно «Загрузка проекта» со списком файлов, имеющих расширение «.prt».

Перейти в папку с каталогом SimInTech (по умолчанию диск «C:»), выбрать «SimInTech\Demo\Автоматика и математика\Регуляторы\Релейные регуляторы\Релейный регулятор для САР ядерного реактора\Релейный регулятор для САР ядерного реактора.prt» и нажать на кнопку «Открыть». Откроется окно проекта со структурной схемой, представленной на рисунке (Рисунок 1), а также откроется график с заголовком «Реактивности в долях β_эфф и относительные отклонения нейтронной мощности».

Рассмотреть САР в которой управляющее воздействие от блока «Ступенька» с подписью «Задатчика мощности» равно нулю, то есть которая выведена из состояния равновесия (начальная реактивность температурной обратной связи ненулевая, а равна «1%» от β_эфф). Левой кнопкой мыши выполнить нажатие на кнопку «Инициализация» в панели инструментов. Таким образом задача будет инициализирована (в нижнем левом углу окна появится сообщение о готовности проекта к старту: Старт). Выполнить нажатие на кнопку «Пуск» левой кнопкой мыши. Начнется процесс моделирования.

Полученные результаты свидетельствуют, что в реакторе установился режим высокочастотных автоколебаний мощности (Рисунок 2), амплитуда которых примерно в два раза превышает уставку по зоне нечувствительности равной одному проценту (блок «Релейное с зоной нечувствительности» с подписью «Управляющее реле»). Причинами автоколебаний являются узкая зона нечувствительности в блоке с подписью «Управляющее реле», а также относительно высокая скоростная эффективность в блоке «Инерционно-интегрирующее звено» с подписью «Привод СУЗ».

Процесс моделирования можно начать без инициализации проекта, выполнив однократной нажатие на кнопку «Пуск».

Двойным нажатием левой кнопкой мыши по блоку с подписью «Управляющее реле» открыть окно «Свойства» этого блока. При помощи шести строк с именами «a1», «a», «b», «b1», «Y1», «Y2» задать в поле «Значение» свойства блока равные «-0.02», «-0.016», «0.016», «0.02», «-1», «1», а в строке «Начальные условия (-1=Y1; 0=0; 1=Y2)» задать значение равное «0» согласно рисунку (Рисунок 3). Настройка этих параметров соответствует зоне нечувствительности по относительному отклонению нейтронной мощности «±2%», коэффициенту возврата «0.8»и значению погрешности при срабатывании реле «±1.0». Снова запустить проект на расчет. Характер поведения графиков свидетельствует, что высокочастотных автоколебаний нет и САР асимптотически возвращается в равновесное состояние (Рисунок 4).

Задать свойства блока с подписью «Управляющее реле» согласно рисунку (Рисунок 5). Двойным нажатием левой кнопкой мыши на блок с подписью «Привод СУЗ» открыть окно «Свойства». Задать свойство «Коэффициент усиления» равным «1e-4».

Запустить проект на расчет. По окончании расчета проанализировать график: высокочастотных автоколебаний не будет, а блок с подписью «Привод СУЗ» изменит своё движение два раза (Рисунок 6).

Рассмотреть неавтономную САР, которая десять секунд находится в равновесии, а после на неё подается управляющее воздействие, которое должно перевести ядерный реактор на повышенный уровень мощности «+10 %».

Двойным нажатием левой кнопкой мыши на блок «Инерционное звено 1-го порядка» с подписью «Обратная связь» вызвать окно «Свойства» и изменить значение «Начальное условие» на «0». Открыть окно «Свойства» блока с подписью «Задатчик мощности» и установить в строки с именами «t», «y0», «yk» значения «10» «0» «0.1» соответственно. Таким образом на вход блока «Сравнивающее устройство» будет подано ступенчатое воздействие равное «0.1» через десять секунд после начала расчёта.

Нажать на кнопку «Параметры расчета» и задать параметр «Конечное время расчёта» равным «100». Запустить проект на расчёт.

Анализ полученных результатов показывает, что ядерный реактор переведен на заданный уровень мощности с точностью до ширины зоны нечувствительности, регулирующий стержень внес дополнительную положительную реактивность «~6...7%» от β_эфф, а реактивность ядерного реактора в переходном процессе в максимуме достигает «~6%» от β_эфф и стремится к нулю обратно пропорционально времени (Рисунок 7).

Для отображения результатов расчета использована векторизованная обработка сигналов: сигналы реактивностей регулирующего стержня и реактора преобразуются в вектор с помощью блока «Мультиплексор», а затем поэлементно умножаются на коэффициент блока «Усилитель». Таким образом сигналы векторно нормированы на значение эффективной доли запаздывающих нейтронов β_эфф.

Блок «Временной график» в данной задаче имеет два входа (первый вход – векторный, а второй – скалярный). С подробным описанием работы других блоков можно ознакомиться в справочной системе SimInTech.

После ознакомления с принципом работы данной схемы необходимо закрыть окно проекта без сохранения изменений. При необходимости можно ознакомиться с другими демонстрационными примерами в папке «Demo» директории установки SimInTech.

Исходные данные для ознакомительной задачи

Цель: ознакомиться с принципами работы в SimInTech, выполнив все этапы построения модели динамики САР, представленной на рисунке (Рисунок 8).

Задача работы: необходимо найти параметр «Коэффициент усиления» в блоке «Интегратор» для формирования САР удовлетворяющей критериям:

• При подаче ступенчатого воздействия перерегулирование отсутствует.
• Время переходного процесса не превышает 20 с.

Формирование структурной схемы

Этап 1 – подготовка структурной схемы.

Нажать на кнопку «Новый проект» и выбрать пункт «Схема модели общего вида», откроется новое окно проекта. Поместить в окно проекта блок «Ступенька» с вкладки «Источники», два блока «Сравнивающее устройство» с вкладки «Операторы», блок «Интегратор» с вкладки «Динамические», «Инерционное звено 1-го порядка», и «Колебательное звено», с вкладки «Вывод данных» добавить блок «Временной график» с вкладки «Вывод данных». Разместить блоки в окне проекта и задать им имена согласно рисунку (Рисунок 8).

Для создания подписи к блокам необходимо выделить нужный блок. При выделении блока появляется дополнительное окно подписи под блоком. Двойным нажатием левой кнопки мыши на окно подписи, курсор устанавливается в окне: задать подпись блока.

Этап 2 – соединение блоков линиями связи.

Нажать на выходной порт блока с подписью «Управляющее воздействие» и, отпустив кнопку, протянуть линию связи к входному порту блока с подписью «Главное сравнивающее устройство» и нажать на него левой кнопкой мыши. Для осуществления поворота линии связи необходимо нажать левой кнопкой мыши по свободному пространству окна проекта, где необходимо выполнить поворот и продолжить построение линии связи.

Изменить расположение портов блока с подписью «W₃(s)». Двойным нажатием левой кнопкой мыши по блоку вызвать окно «Свойства». В пункте «Порты» выбрать необходимой порт и установить его расположение на блоке согласно рисунку (Рисунок 8).

Повторить соединение линиями связи для остальных блоков структурной схемы.

При необходимости подключения входа блока к существующей линии связи необходимо провести линию связи от входа блока до необходимой линии связи и нажать на нее левой кнопкой мыши для подключения. Для удаления линии связи необходимо выделить данную линию и нажать кнопку «Delete».

Сохранить проект. Для этого открыть меню «Файл» в главном окне и выбрать пункт «Сохранить проект как...». В появившемся окне задать название проекта (например «lesson1»). Нажать на кнопку «Сохранить».

Для масштабирования структурной схемы по окну проекта необходимо зайти в меню «Масштаб» и выбрать пункт «Показать все»

Сохранить проект, нажатием на кнопку «Сохранить».

Этап 3 – настройка свойств блоков.

Двойным нажатием левой кнопкой мыши на блок с подписью «Управляющее воздействие» открыть окно «Свойства». Во вкладке «Свойства» задать значения параметрам блока согласно рисунку (Рисунок 9).

Ввести соответствующие параметры для блоков с подписью «W₂(s)»и «W₃(s)» согласно рисункам (Рисунок 10 и Рисунок 11).

Двойным нажатием по блоку «Временной график» открыть окно с заголовком «График». Задать необходимый размер окна.

Сохранить проект.

Примечание: окно «Свойства» можно открыть нажатием правой кнопки мыши по блоку и выбором пункта «Свойства объекта» (Рисунок 12).

Открыть окно «Свойства» блока с подписью «W₁(s)» и в строке «Коэффициенты усиления» задать значение равным «1».

Сохранить проект.

Этап 4 – установка параметров расчета.

Нажатием на кнопку «Параметры расчета» открыть окно «Параметры проекта» (Рисунок 13).

Вкладки окна «Параметры проекта» предназначены для:

• «Параметры расчёта» – параметры математического ядра, используемые при расчете, зависят от выбранного решателя и типа проекта.
• «Синхронизация» – настройка синхронизации времени расчета с реальным временем.
• «Рестарт» – параметры сохранения и запуска исходных состояний расчетной модели (параметры проекта сохраняются в виде файла с расширением «.rst»).
• «База данных» –структурированная файловая база данных, содержащая в себе сигналы, используемые в одном или нескольких проектах (представляет из себя файл с расширением «.db»).
• «Вид» – параметры внешнего вида расчетных схем при моделировании и редактировании.
• «Настройки» – дополнительные настройки внешних библиотек.

Перейти во вкладку «Параметры расчёта». В пункте «Метод интегрирования» выбрать из списка «Адаптивный 1». Установить значение в пункте «Конечное время расчёта» равным «40», в пункте «Минимальный шаг» задать значение равным «0.001», а в пункте «Максимальный шаг» значение равным «0.1».

Сохранить проект.

Моделирование переходных процессов

Запустить проект на расчёт нажатием на кнопку «Пуск». По окончании расчета появится сообщения в нижней части она проекта: Ошибка: Заданная точность не обеспечивается.

Открыть окно «Параметры расчета» и задать параметр «Минимальный шаг» равный «1e-10» и повторите процесс моделирования.

При изменении минимального шага расчёта внешний вид переходного процесса не изменился, так как при первоначальном минимальном шаге интегрирования заданная точность не обеспечивалась только на первом шаге моделирования. Сообщение о точности можно проигнорировать.

Двойным нажатием по блоку «Временной график» вызвать окно «График». Нажать правой кнопкой мыши по окну «График» и выбрать пункт «Свойства». На вкладке «Общие» и введите название графика «График переходного процесса при К=1» в поле «Заголовок».

На вкладке «Графики и оси» задать подписи осям графика в полях «Название оси» согласно рисунку (Рисунок 14).

При установке в блоке «Интегратор» значения «Коэффициент усиления» равным «1», данные расчета свидетельствуют, что САР неустойчива, а переходной процесс становится расходящимся (Рисунок 15).

При необходимости в окне «Свойства графика» можно задать: тип графика, цвет и стиль линии, цвет фона, масштабы и нормировочные значения по обеим осям.

Что бы закрыть окно «Свойства графика» необходимо нажать на кнопку «Ok», при этом все изменения будут сохранены. Сохранить проект.

В блоке «Интегратор» задать свойство «Коэффициент усиления» равным «0.2». Запустить проект на расчёт.

Данные расчёта свидетельствуют, что перерегулирование отсутствует, но время переходного процесса значительно превышает необходимые критерии.

Изменить название графика согласно рисунку (Рисунок 16).

В блоке «Интегратор» задать свойство «Коэффициент усиления» равным «0.35». Запустить проект на расчёт.

Данные расчёта свидетельствуют, что перерегулирование отсутствует и время переходного процесса соответствует необходимым критериям.

Изменить название графика согласно рисунку (Рисунок 17).

ПРЯМОЕ МОДЕЛИРОВАНИЕ ПЕРЕХОДНЫХ ПРОЦЕССОВ

Сформировать структурную схему простейшей математической модели динамики САР ядерного реактора (САР ЯР) и определить динамические свойства САР прямым моделированием переходных процессов при подаче управляющих и возмущающих воздействий.

Структурная схема САР ЯР имеет примерно следующий вид (Рисунок 18), где:

• u^~(t)=(u(t)-u₀)/u₀ - нормированное отклонение управляющего напряжения в обмотке возбуждения электродвигателя;
• n^~(t)=(n(t)-n₀)/n₀ - нормированное отклонение плотности нейтронов (мощности или нейтронного потока);
• p_ст(t), p_возм(t), p_ос(t) - реактивность, вносимая в реактор регулирующим стержнем, внешним возмущающим воздействием, местной (внутренней) обратной связью, соответственно;
• ε(t)=u^~(t)-n^~(t) - рассогласование (ошибка).

Для изображения двух блоков в структурной схеме (Рисунок 18) использована пиктограмма в виде вопросительного знака. Это означает, что для описания динамики этих блоков необходимо определить вид каждой передаточной функции и решить, какое типовое звено использовать.

Для изображения блока с подписью «Кинетика нейтронов» использована пиктограмма блока «Передаточное звено общего вида» из библиотеки «Динамические», что можно рассматривать как подсказку о том, какой типовой блок описывает кинетику нейтронов в этой задаче.

На схеме не изображен ряд блоков, например, блоков графического отображения результатов расчета и блок преобразования сигналов. Каких блоков из библиотек «Вывод данных», «Операторы», «Функции»не хватает необходимо определить самостоятельно.

Необходимо построить на разных графиках зависимости величин n(t), T(t), d[p_ст(t)/β_эфф]/dt, а также на одном графике реактивности p_ст(t), p_возм(t), p_ос(t) в долях β_эфф в зависимости от модельного времени. Самостоятельно найти способы формирования и отображения этих динамических переменных.

Известно, что переход к нормированным отклонениям плотности нейтронов (что эквивалентно нормированным отклонениям мощности или нейтронного потока) и последующая линеаризация дифференциального уравнения для плотности нейтронов, позволяют представить математическую модель точечной кинетики нейтронов с одной эффективной группой запаздывающих нейтронов в следующем виде:

Где β_эфф – эффективная доля запаздывающих нейтронов; l – время жизни мгновенных нейтронов; λ – постоянная распада ядер-предшественников запаздывающих нейтронов; c^~(t)=[c(t)-c₀]/c₀ – нормированные отклонения концентрации ядер-предшественников запаздывающих нейтронов.

Блок с подписью «Местная обратная связь», определяется отрицательным температурным эффектом реактивности и описывается следующими уравнениями:

где: α – температурный коэффициент реактивности;

T0, T^~(t) – стационарная температура топлива в активной зоне и нормированное отклонение температуры топлива в активной зоне от стационара, соответственно;

τос – постоянная времени (инерционность) топлива в активной зоне;

A – безразмерный коэффициент.

Блок с подписью «Привод регулирующего стержня» состоит из электродвигателя постоянного тока, редуктора, муфт, преобразователя движения, непосредственно регулирующего стержня и т.п., однако для упрощения задачи все эти элементы объединены в одно звено.

Нестационарные процессы в блоке с подписью «Привод регулирующего стержня» описываются следующим дифференциальным уравнением:

где K_пр – коэффициент скоростной эффективности; τ_пр – постоянная времени.

Для выполнения самостоятельной части лабораторной работы необходимо выполнить следующие этапы:

• Представить динамику САР в переменных «вход-выход» (в передаточных функциях), определив вид передаточных функций для всех блоков структурной схемы.
• Подставить численные значения всех параметров и вычислить значения всех коэффициентов в полиномах числителей и знаменателей передаточных функций.
• Используя освоенные процедуры работы в SimInTech, сформировать структурную схему динамики САР ядерного реактора.
• При моделировании переходных процессов необходимо построить зависимости: n(t)/n₀, T(t), d(p_ст(t)/β_эфф)/dt,а также на одном графике реактивности в долях β_эфф в зависимости от времени. Найти способы формирования и отображения динамических переменных.
• Представить сформированную структурную схему (с введенными значениями свойств звеньев, блоков формирования и отображения сигналов) преподавателю для проверки ее корректности.
• Выполнить пробное моделирование переходного процесса в исходной САР при отсутствии возмущающего воздействия и подаче управляющего воздействия u(t) = 0.05×1(t) для конечного времени моделирования равного «40», и на основании вида переходных процессов сделать вывод об устойчивости исходной САР и найти способ последующей коррекции САР посредством изменения ее параметров.
• Методом последовательных приближений достичь устойчивости САР.
• Выполнить моделирование переходного процесса в скорректированной САР при отсутствии возмущающего воздействия и подаче управляющего воздействия u(t) = 0.05×1(t) для конечного времени моделирования равного «400», и на основании результатов расчета (по графикам) выполнить анализ поведения отображаемых динамических переменных.
• Выполнить моделирование переходных процессов в скорректированной САР при отсутствии управляющего воздействия и подаче возмущающего воздействия pвозм(t)=0.1·β_эфф·1(t) и pвозм(t)=0.7·β_эфф·1(t) (поочередно) для конечного времени моделирования равным «100», и на основании результатов расчета выполнить анализ поведения отображаемых динамических переменных.
• Выполнить моделирование переходных процессов при подаче вышеуказанных возмущающих воздействий для нового значения температурного коэффициента реактивности α = 5·10-5 K^-1, и на основании результатов расчета выполнить анализ поведения отображаемых динамических переменных.

Исходные данные: β_эфф=0.006, l=0.072 c^-1, τ_пр=0.2 с^-1, T₀=700 K, A=0.75, α=1·10-4 K^-1.

Внимание: ряд параметров (β, λ, τ_пр, T₀, A, α) одинаковы для всех подгрупп, а значения параметров l, τ_ос, K_пр будут заданы преподавателем индивидуально для каждой подгруппы.

Заключение

Демонстрационно-ознакомительная задача на этом завершена. Необходимо сохранить проект и сформировать отчёт о проделанной работе.