Исследование и построение причинно-следственной математической модели системы управления частотой вращения ротора турбины парового турбомашинного комплекса. Реализация нелинейной модели в среде Matlab/Simulink. Выбор параметров регулятора процессов.
При низкой оригинальности работы "Моделирование системы управления частотой вращения ротора парового турбоагрегата", Вы можете повысить уникальность этой работы до 80-100%
Схема реализации нелинейной модели системы управления в среде Matlab/Simulink представлена на рисунке 1.Определим состояние равновесия системы для номинального режима (давление пара, поступающего из внешней паровой емкости, , коэффициент нагрузки ) с использованием команды пакета Matlab trim: " [sys,x0,xstr]=static " u0=40.79 " y0=1.30 " ix= [] " iu= [1] " iy= [] " [x,u,y,dx]=trim("static",x0,u0,y0,ix,iu,iy) x =1.0e 03 * 0.0010 0.0010 0.0010 0.0010 2.5000 0.0700 u = 1.0000 y = 1.0000 Теперь, используя в качестве начального вектора равновесия значения, отвечающие номинальному режиму, найдем состояние равновесия для заданного режима: " [x,u,y,dx]=trim("static",x0,u0,y0,ix,iu,iy) x = 1.0e 03 * 0.0010 0.0005 0.0005 0.0005 1.2542 0.0351 u = 1 y = 1.0000 Осуществим поиск других состояний равновесия при заданном фиксированном воздействии: в цикле будем изменять значения начального вектора и, фиксируя результат, выполнять команду trim. Анализ полученного результата позволяет сделать вывод о том, что в диапазоне 0,01-1 для каждой из координат вектора состояния другого состояния равновесия для заданного режима не найдено.Для моделирования системы при переходе с номинального на заданный режим выберем метод интегрирования ode4 (Рунге-Кутта) с шагом интегрирования h=0,001, t = 10000.Максимально допустимый шаг для каждого из методов численного интегрирования представлен в таблице 2. В таблице 2 для каждого метода приведена длительность моделирования поведения системы на временном интервале [0 c; 15000 c]. По рисункам видно, что при величине шага h из таблицы 2 методы остаются устойчивыми, а при незначительном увеличении шага устойчивость теряется. Euler (ode1) 0.0279 226.69. 0,5 24,7 % / 165,9 21,2 % / 110,3 18,9 % / 91,7 19,4 % / 91,2 21,1 % / 75,2.По известной передаточной функции регулятора найдем его дискретную передаточную функцию .Заменим использовавшийся ранее непрерывный регулятор дискретным. Структурная схема системы в среде Matlab/Simulink представлена на рисунке 15. Рисунок 15. Структурная схема системы с дискретным регулятором Проведем анализ поведения системы с дискретным регулятором при переходе с номинального на заданный режим (рисунки 21-23). Рисунок 16. Переходный процесс на выходе системы
Вы можете ЗАГРУЗИТЬ и ПОВЫСИТЬ уникальность своей работы
