Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,441

СИНТЕЗ СИСТЕМ РЕГУЛИРОВАНИЯ ПЕРВИЧНЫХ ДВИГАТЕЛЕЙ СИНХРОННЫХ ГЕНЕРАТОРОВ

Мельников Д.В. 1 Чжо Ту Аунг 1 Окар Мин 1 Мин Чжо Ту 1
1 Калужский филиал ФГБОУ ВО «Московский государственный технический университет имени Н.Э. Баумана (национальный исследовательский университет)»
В статье на базе современных проекционно-матричных методов представлен алгоритм параметрического синтеза системы регулирования первичного двигателя электрических синхронных генераторов, которым служит паровая турбина. Алгоритм основан на представлении процессов в форме разложения по базисной системе функций, использует аппарат структурных преобразований и методы математического программирования. Отличительной особенностью алгоритма является его ориентация на класс нелинейных систем, причем количество нелинейных элементов не имеет принципиальных ограничений и существенно не отражается на быстродействии вычислительных алгоритмов. Это стало возможно благодаря разработке «быстрых» алгоритмов замены нелинейных элементов соответствующими эквивалентными матричными операторами. В случае детерминированной постановки нелинейный элемент заменяется эквивалентным линейным элементом с переменным коэффициентом передачи, что эквивалентно матричному оператору умножения. Критерий эквивалентности – равенство выходных сигналов нелинейного и эквивалентного ему линейного элемента, заданного матричным оператором, при отработке конкретного входного сигнала. Алгоритм приводится на примере конкретной энергетической турбины, математическая модель которой описывается структурной схемой.
энергетическая турбина
генератор
матричный оператор
ортонормированный базис
система регулирования
синтез регулятора
алгоритм
математическая модель
1. Аунг Ч.С., Макаренков А.М., Мьо П.С. Идентификация случайных параметров математической модели электрогидравлического следящего привода // Фундаментальные исследования. – 2016. – № 2–2. – С. 231–235.
2. Гусев А.С., Свечкарев С.В., Плодистый И.Л. Математическая модель первичных двигателей синхронных генераторов // Известия Томского политехнического университета. – 2005. – Т. 308. – № 7. – С. 216–221.
3. Корнюшин Ю.П, Мельников Д.В., Егупов Н.Д., Корнюшин П.Ю. Исследование и расчет параметров элементов системы регулирования частоты вращения ротора турбины с учетом параметрической неопределенности математической модели // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. – № 1. – С. 78–93.
4. Мельников Д.В., Егупов Н.Д. Синтез систем регулирования энергетических турбин в условиях параметрической неопределенности // Известия Тульского государственного университета. Технические науки. – 2011. – № 5–1. – С. 108–113.
5. Мельников Д.В., Корнюшин Ю.П., Мин Чжо Ту, Чжо Ту Аунг, Окар Мин. Алгоритм синтеза системы регулирования частоты вращения ротора энергетической турбины // Научное обозрение. – 2015. – № 20. – С. 138–143.
6. Мельников Д.В., Корнюшин П.Ю., Мин Ч.Т., Чжо Т.А., Окар М. Проекционно-матричный подход к анализу и синтезу систем управления электроэнергетических систем // Научное обозрение. – 2015. – № 2. – С. 88–97.
7. Мельников Д.В. Метод автоматизированного исследования систем регулирования энергетических турбин при случайных возмущениях: дис. ... канд. техн. наук. – Калуга, 2002.
8. Мельников Д.В., Мин Чжо Ту. Математическая модель систем регулирования энергетических турбин с регулируемыми отборами пара // Инженерный журнал: наука и инновации. – 2015. – № 2 (38). – С. 13. – http://engjournal.ru/catalog/pmce/tctp/1366.html.
9. Мельников Д.В., Окар Мин, Чжо Ту Аунг. Математическое моделирование систем регулирования энергетических турбин // Научное обозрение. – 2015. – № 20. – С. 144–149.
10. Мельников Д.В. Проекционно-матричный метод синтеза контура регулирования частоты вращения ротора паровой турбины // Вестник Московского государственного технического университета им. Н.Э. Баумана. Серия: Машиностроение. – 2013. – № 4 (93). – С. 43–53.
11. Мельников Д.В., Фишер М.Р. Математическая модель контура регулирования частоты вращения ротора паровой турбины К-800-130 // Вестник Московского государственного технического университета им. Н.Э. Баумана. Серия: Машиностроение. – 2011. – № 7. – С. 197–215.
12. Нормы участия энергоблоков тепловых электростанций в нормированном первичном регулировании частоты и автоматическом вторичном регулировании частоты и перетоков активной мощности // Стандарт организации. Открытое акционерное общество «Системный оператор единой энергетической системы». Издание официальное. – М., 2013. – http://so-ups.ru/fileadmin/files/laws/standards/sto_002-2013_freq_regulation.pdf.
13. Пупков К.А., Егупов Н.Д., Колесников Л.В., Лукашенко Ю.Л., Мельников Д.В., В.М. Рыбин, Трофимов А.И. Матричные методы расчета и проектирования сложных систем автоматического управления для инженеров / под ред. К.А. Пупкова и Н.Д. Егупова. – М.: Издательство МГТУ им. Н.Э. Баумана, 2007. – 664 с.
14. Окар Мин, Мельников Д.В. Алгоритм расчета нелинейных систем управления проекционно-матричным методом // Инженерный журнал: наука и инновации. – 2014. – № 12 (36). – С. 17. – http://engjournal.ru/catalog/it/asu/1268.html.
15. Чжо Ту Аунг, Мельников Д.В. Алгоритм исследования нелинейных систем автоматического управления в стохастических режимах // Инженерный журнал: наука и инновации. – 2014. – № 4 (28). – С. 3. – http://engjournal.ru/catalog/it/asu/1270.html.

Одним из наиболее сложных элементов современных энергетических систем является первичный двигатель электрических синхронных генераторов, которым обычно служит паровая турбина со своими системами регулирования. Основным способом получения информации, необходимой для проектирования, эксплуатации и развития электроэнергетических систем, является математическое моделирование [2]. Динамика первичных двигателей существенно влияет на нормальные и аварийные процессы производства, распределения и потребления электроэнергии. Поэтому требование к уровню адекватности модели первичного двигателя постоянно возрастает. А это означает, что на сегодняшний день описывать модели систем регулирования приходится в классе нелинейных систем.

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

Решение поставленной задачи осложнено тем, что в настоящее время в электроэнергетических системах находится в эксплуатации множество различных типов турбин и котлоагрегатов с разнообразными системами и законами регулирования. Хотя передаточные функции отдельных элементов и звеньев известны, разработка полной математической модели для конкретного типа оборудования и компоновки первичного двигателя, а также расчета системы регулирования является уникальной [2, 8, 9, 11]. Поэтому возникает задача разработки универсального алгоритма синтеза, который может использоваться для стационарных и нестационарных систем, как для линейных, так и нелинейных систем управления. Такие алгоритмы могут быть построены на использовании метода матричных операторов расчета и проектирования сложных систем [6, 13]. Для иллюстрации таких алгоритмов в работе приводится синтез системы регулирования турбины производства Калужского турбинного завода. Хотя модель не является громоздкой, но на ее примере довольно подробно иллюстрируются все особенности алгоритма, что не снижает его общности по отношению к системам высокого порядка [10].

Общий алгоритм синтеза системы регулирования

Рассмотрим систему регулирования турбины ПТ-25/30-90/10М производства Калужского турбинного завода в классе нелинейных систем (рис. 1) [7]. Эта система предназначена для поддержания частоты вращения ротора турбины на требуемом уровне.

На рис. 1 обозначено: φ – относительное отклонение угловой частоты вращения ротора от номинального значения ωНОМ, Melnikov01.wmf; λ – относительное отклонение электрической нагрузки генератора; ξ – относительное изменение расхода пара в турбине; μ – относительное отклонение поршня сервомотора; z – относительное отклонение поршня золотника; δ – степень неравномерности системы регулирования; Tz – постоянная времени золотника; Tμ – постоянная времени сервомотора; Tξ – постоянная времени камеры эквивалентного объема (время полного заполнения камеры при определенном расходе пара); Tφ – постоянная времени ротора; θ – эквивалентный коэффициент самовыравнивания.

pic_22.wmf

Рис. 1. Структурная схема нелинейной системы регулирования энергетической турбины ПТ-25/30-90/10М (ПТ-12/15-35/10М)

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

Melnikov02.wmf

где Δz0 – относительные высоты треугольного профиля кромок поршня отсечного золотника; ε0 – величина, характеризующая степень нечувствительности сервомотора.

Требуется определить параметры системы регулирования δ, Tz, Tμ таким образом, чтобы φЭ(t) и φ(t) (реальный переходной процесс) были близки, например, в среднеквадратическом смысле:

Melnikov03.wmf

Представим модель системы регулирования турбины с нелинейностью в операторной форме (в матричном виде) (рис. 2).

На схеме обозначено:

Melnikov04.wmf

Melnikov05.wmf Melnikov06.wmf

Melnikov07.wmf Melnikov08.wmf

где AИ – матричный оператор интегрирования. Из рис. 2 имеем следующие векторно-матричные соотношения:

Melnikov09.wmf Melnikov10.wmf

Melnikov11.wmf Melnikov12.wmf Melnikov13.wmf

Melnikov14.wmf

или

Melnikov15.wmf Melnikov16.wmf

pic_23.wmf

Рис. 2. Структурная схема нелинейной системы регулирования энергетической турбины в операторной форме

Последнее выражение представляет собой матричный оператор всей системы, который связывает спектральные характеристики входного и выходного процессов. Матричный оператор нелинейного элемента вычисляется по формуле [6]:

Melnikov17.wmf (1)

Для вычисления спектральных характеристик остальных процессов системы нужно воспользоваться структурными преобразованиями. В результате можно получить следующие зависимости:

Melnikov18.wmf Melnikov19.wmf Melnikov20.wmf Melnikov21.wmf Melnikov22.wmf

Алгоритм, определяющий процедуру решения задачи расчета параметров системы управления турбины, представлен на рис. 3.

pic_24.wmf

Рис. 3. Структурная схема алгоритма синтеза системы регулирования турбины в классе нелинейных систем

Вычислительный алгоритм синтеза системы регулирования первичного двигателя

Рассмотрим подробно алгоритм параметрического синтеза системы регулирования турбины ПТ-25/30-90, представленного на рис. 3.

Этап 1. При заданной структуре регулятора формируем вектор искомых параметров: p = [δ, Tz, Tμ]. Здесь же инициализируем вектор параметров начальными приближениями p0 = [1, 1, 1].

Этап 2. Накладываем ограничения на параметры, связанные с их физической реализуемостью. Эти ограничения будут пересланы в функцию оптимизации в качестве опции. Данная операция позволяет значительно сократить время поиска оптимальных параметров. Например, все полученные параметры регуляторов должны быть положительными.

Этап 3. Задаем эталонный переходный процесс при резком сбросе нагрузки (λ = –1):

Melnikov23.wmf

Стоит обратить особое внимание на то, чтобы эталонный переходный процесс соответствовал требованиям, предъявляемым в [12].

Этап 4. Выбор ортонормированного базиса на промежутке исследования [0, T]. В дальнейшем всегда будем использовать ортонормированный базис функций Уолша, упорядоченных по Адамару, и обозначать его H – квадратная симметричная матрица размера l×l. Вычисление спектральных характеристик входного и выходного воздействий.

Melnikov24.wmf Melnikov25.wmf

где λ, φЭ – вектор-столбцы отсчетов функций λ(t), φЭ(t) соответственно.

Этап 5. Определяем матричные операторы тех звеньев, параметры которых известны (как правило, это параметры объекта и часть параметров системы управления):

Melnikov26.wmf

Melnikov27.wmf

где AИ – матричный оператор интегрирования; I – единичная матрица; θ = 0,05 (коэффициент саморегулирования); для ПТ-25/30-90/10М: Tφ = 12 c, Tξ = 0,5 c. Эти матричные операторы, так же как и спектральные характеристики, рассчитываются один раз и хранятся в памяти ЭВМ.

Этап 6. Определяем матричные операторы звеньев, параметры которых подлежат определению (матричные операторы вычисляются при текущих значениях вектора искомых параметров p = [δ, Tz, Tμ]).

Melnikov28.wmf Melnikov29.wmf

Melnikov30.wmf

Здесь же можно задать начальное приближение для выходного сигнала (или его спектральной характеристики Melnikov31.wmf).

Этап 7. Вычисление матричных операторов нелинейных звеньев. Для этого необходимо предварительно найти спектральные характеристики процессов на входе нелинейных элементов:

Melnikov32.wmf Melnikov33.wmf

Melnikov34.wmf Melnikov35.wmf z = HCz;

Melnikov36.wmf

здесь Melnikov37.wmf – матричный оператор умножения функции представленной в скобках; z – вектор-столбец отсчетов функции z(t). Непосредственное применение формулы (1) для вычисления матричного оператора умножения является далеко не эффективным с точки зрения быстродействия его расчета. Так как система Matlab ориентирована на операции с матрицами, то для расчета матричного оператора умножения эффективно использовать следующую функцию:

function Ay=m_ymn1(H,f)

%Матричный оператор умножения на некоторую функцию %f – вектор-столбец

%H – матрица-базис функций Уолша

Ay=H*(H.*(ones(size(f))*f’)’);

Этап 8. Вычисление матричного оператора всей системы и спектральной характеристики выходного сигнала

Melnikov38.wmf

Melnikov39.wmf

Этап 9. Уточнение матричных операторов нелинейных элементов. Для этого проверяем решения задачи анализа при текущих значениях p: степень близости спектральных характеристик выходного сигнала, полученных на текущей и предыдущей итерациях.

Melnikov40.wmf (2)

Если (2) не выполняется, то Melnikov41.wmf, и переходим к этапу 7, если выполняется, то переходим к этапу 10.

Этап 10. Вычисление искомых параметров:

Melnikov42.wmf (3)

Реализация (3) представляет собой итерационный процесс, в котором параметры меняются в соответствии с выбранным алгоритмом нахождения минимума функции многих переменных. При очередном вычислении целевой функции (3) при конкретных значениях p, следует определять Cφ(p), т.е. возвращаться к этапу 6. Если минимум найден, то переходим к анализу системы при найденных параметрах регулятора.

Этап 11. Анализ системы регулирования турбины при найденных параметрах регулятора:

φ = HCφ; z = HCz; Melnikov43.wmf

μ = HCμ; ξ = HCξ,

где φ, Melnikov44.wmf z, μ , ξ  – вектор-столбцы отсчетов функций φ(t), Melnikov45.wmf z(t), μ(t), ξ(t) соответственно.

pic_25.wmf

Рис. 4. Эталонный и реальный переходные процессы относительного изменения частоты вращения ротора

Результаты синтеза системы регулирования турбины ПТ-25/30-90/10М представлены на рис. 4 при λ = –1, ε0 = 0,0167. Размер базиса составлял 512 функций Уолша. В качестве меры близости спектральных характеристик в выражениях (2) и (3) использовалась евклидова норма. Минимизация целевой функции (3) осуществлялась методом Гаусса – Ньютона, который реализован в среде Matlab с помощью функции lsqnonlin. Эта функция позволяет осуществлять условную минимизацию целевых функций, представляющих собой свертку вектора.

Для системы регулирования турбины ПТ-25/30-90/10М были найдены следующие параметры (ограничения на параметры не накладывались):

δ = 0,03; Tz = 0,0630 с; Tμ = 0,1014 с; (4)

Melnikov46.wmf

Для получения параметров (4) понадобилось 8 итераций. Степень нечувствительности сервомотора составляла ε0 = 0,0167, что соответствует 0,05 % нечувствительности системы регулирования относительно регулируемого процесса φ(t). В работах [1, 3, 4, 5, 14, 15] отражены детальные алгоритмы построения законов регулирования энергетических систем, в том числе и в робастной постановке, а также алгоритмы решения задач идентификации.

Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований и Правительства Калужской области (гранты № 14-48-03013, № 16-41-400701).


Библиографическая ссылка

Мельников Д.В., Чжо Ту Аунг, Окар Мин, Мин Чжо Ту СИНТЕЗ СИСТЕМ РЕГУЛИРОВАНИЯ ПЕРВИЧНЫХ ДВИГАТЕЛЕЙ СИНХРОННЫХ ГЕНЕРАТОРОВ // Фундаментальные исследования. – 2016. – № 10-3. – С. 509-515;
URL: http://www.fundamental-research.ru/ru/article/view?id=40886 (дата обращения: 08.05.2021).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074