ISSN 1992-6502 (Print)_
2021. Т. 25, № 4 (94). С. 101-110
Вестник УГАТУ
ISSN 2225-2789 (Online) http://journal.ugatu.ac.ru
УДК 621.22-546 doi 10.54708/19926502_2021_25494101
Исследование условий возникновения колебательных процессов
при автоматизированном управлении насосными комплексами
на основе прецедентов
1 2 3
а. г. Лютов , м. б. Новоженин , в. а. Огородов
1 lutov1@mail.ru, 2 novozhenin.maxim@yandex.ru, 3 ogorodov.v@ugatu.su
1 3 ФГБОУ ВО «Уфимский государственный авиационный технический университет» (УГАТУ) 2 ФГБОУ ВО «МИРЭА - Российский технологический университет» (МИРЭА-РТУ)
Поступила в редакцию 07.10.2021
Аннотация. Прецедентный подход, применяемый для автоматизированного управления режимами работы насосных комплексов, позволяет уменьшить или полностью устранить различные нестационарные гидравлические процессы, в том числе кавитацию. Рассмотрены условия возникновения в насосных комплексах колебаний гидродинамического характера, способствующих возникновению или активизации указанных нежелательных процессов, с построением моделей на основе метода электроаналогии и применением критериев устойчивости теории автоматического управления. Полученные результаты позволяют осуществлять формирование и использование базы прецедентов при управлении режимами работы насосных комплексов с учетом колебательности динамических процессов.
Ключевые слова: центробежный насос; прецедентный подход, колебательные процессы; напорно-расходная характеристика; метод электроаналогии.
АКТУАЛЬНОСТЬ
Технологические процессы перекачивания жидкостей широко используются в различных отраслях промышленности - в коммунальном хозяйстве, водоснабжении и водоотведении, при транспортировке нефтепродуктов и т.д., а также в авиационной технике [1]. Для этого обычно применяются насосные комплексы (НК), под которыми подразумеваются один или несколько насосов с прилегающими к ним трубопроводами. Наибольшее распространение получили НК на основе центробежных насосов (ЦН) в силу их значительной производительности, высокого КПД и удобства в эксплуатации.
Среди способов управления режимами работы ЦН можно выделить управление на основе прецедентного подхода [2], который позволяет установить взаимосвязь между текущим режимом работы, реализуемым управляющим воздействием и предполагаемым результатом.
К одному из негативных факторов при работе ЦН относится то, что при определенных условиях в НК могут возникать колебания, в том числе гармонические, которые обуславливают или активизируют различного рода нестационарные гидравлические процессы, в частности такие, как кавитация. Она имеет место при локальном снижении давления в ЦН или трубопроводе до давления насыщенных паров с последующим схлопыванием образующихся полостей, заполненных газом или паром [3, 4]. Это явление приводит к эрозии рабочих органов ЦН, шумам, вибрации и пульсациям давления в системе и, как следствие, уменьшению производительности, к снижению КПД, повышению износа и разрушению элементов, уменьшению срока службы насосного оборудования.
С целью снижения влияния указанных негативных факторов при управлении НК (выборе режимов работы) необходим учет условий возникновения колебательных процессов.
В связи с этим, определение условий возникновения колебательных процессов, выявление соответствующих режимных областей при управлении ЦН и предоставление дополнительной информации для обоснованного выбора прецедентов при автоматизированном управлении НК является важной и актуальной задачей.
ВВОДНАЯ ЧАСТЬ
Основной характеристикой НК, отражающей изменение режима его работы является напорно-расходная характеристика (НРХ). С ее помощью можно определить рабочую точку как пересечение характеристики сети (с увеличением расхода 0 происходит увеличение и напора Н) и характеристики ЦН (с увеличением расхода уменьшается напор). Изменение положения рабочей точки приводит к изменению режима работы НК.
Так как колебательные процессы относятся к категории динамических, то для эффективного исследования динамических процессов в НК методами теории автоматического управления для линейных систем [5] целесообразна линеаризация НРХ в области рабочей точки. При этом для описания данных динамических процессов необходимо воспользоваться зна-
СН
чениями дифференциальных сопротивлений в рабочей точке НРХ: Я =
Ян =■
сса
сси.
> 0,
6с
цн
со
< 0, где Яс - дифференциальное сопротивление сети, Яцн- дифференциальное
с
6с
сопротивление ЦН (рис. 1).
н И
пр
т
О
.7 7 /
чл /?цн
м 1
/ г/ * 1 1 V \ 1 ч
1 1
О Орт 0
Рис. 1. Напорно-расходная характеристика центробежного насоса (1) и сети (2)
Существуют различные схемы включения ЦН в НК. В качестве примера рассмотрим типовую схему включения ЦН, приведенную на рис. 2.
1 насос
Рис. 2. Схема включения ЦН
На рис. 3 приведена схема автоматизированной системы управления НК с применением прецедентного подхода [6].
Рис. 3. Структурная схема системы автоматизированного управления НК с применением прецедентного подхода
Задатчик системы формирует управляющее воздействие, соответствующее требуемой круговой частоте вращения рабочего колеса центробежного насоса и^.
Блок определения кавитационного объема с помощью нейронной сети определяет величину кавитации Ъ\ для используемого режима работы на основе весов н из базы прецедентов. Блок оценки производительности при помощи критерия оценки производит определение отклонения производительности насосного комплекса Л1 от заданной.
Блок поиска решения с помощью алгоритма Хука-Дживса формирует набор подходящих прецедентов Р = [Н\ Q'}, используя значения критерия оценки производительности Л1 и кавитационного объема Ъ\.
Блок выбора прецедента производит проверку предлагаемых прецедентов P на математической модели процесса перекачивания жидкости, формирует заключения об их пригодности Р".
Блок формирования корректирующих сигналов предлагает отобранные воздействия оператору для последующего принятия решения о коррекции управления (Диш] насосным комплексом.
База прецедентов, сформированная в режиме off-line, состоит из результатов моделирования в среде Turbomachinery CFD: возможных значений управляющих воздействий с учетом ограничений, статических и динамических параметров и характеристик насосного комплекса, настройки весов нейронной сети.
Прерывания в системе реализованы по времени и по инициативе оператора.
ПОСТРОЕНИЕ ДИНАМИЧЕСКОЙ МОДЕЛИ
Для построения динамической модели процесса перекачивания жидкости с целью эффективного анализа условий возникновения колебательных процессов целесообразно воспользоваться методом электроаналогии [7], получившим применение при моделировании гидравлических процессов.
Рассмотрим модель НК, сформированную с помощью метода электроаналогии (рис. 4).
Рис. 4. Модель насосного комплекса, построенная на основе метода электроаналогии
На рис. 4: Е - электрический эквивалент напора, развиваемого ЦН; I - электрический эквивалент расхода; Яс, Яцн , Явх - эквиваленты гидравлических дифференциальных сопротивлений сети центробежного насоса, всасывающего (входного) трубопровода (соответственно); С - динамический элемент (электрическая емкость), отражающий упругие свойства жидкости и элементов конструкции насосного комплекса; Ь - динамический элемент (индуктивность), отражающий инерционность потока жидкости [8].
Используя приведенную выше модель, составим уравнение электрической цепи в комплексной области в соответствии со 2-м законом Кирхгофа:
1 • Я
E (Jo) = I (Jo) • (R + Яцн + JoL + ^-), (1)
+ R JoC в
или
(2)
1 + JaCRex
После раскрытия скобок, группировки и перехода, в соответствии с преобразованием Лапласа [5], к операторной форме записи ( p = Ja ), получим:
LCR p2 + (CR -(R + R ) + L)p + (R + R + R ) E(p) = I(p)-(-eXp—(——-)p ( c—^-вх)). (3)
1 + CKxp
Соответствующая передаточная функция будет иметь вид:
W(p) = I(P) =-5--. (4)
^ E(p) LCRexp + (CRex-(Rc + RJ + L)p + (Rc + R^ + Rex) K)
АНАЛИЗ УСТОЙЧИВОСТИ ДИНАМИЧЕСКИХ ПРОЦЕССОВ
Для анализа устойчивости и степени колебательности динамической системы с передаточной функцией (4) в рабочей точке НРХ, необходим анализ характеристического полинома (знаменателя) передаточной функции (4):
LCRxp2 + (CRex -(Rc + R^) + L)p + (Rc + RWi + Rex). (5)
Необходимым условием устойчивости динамических процессов в системах второго порядка является положительность всех коэффициентов характеристического уравнения [5]. Так как коэффициент при p2 заведомо положителен, то анализ сводится к получению условий положительности двух других коэффициентов.
В частности, для среднего коэффициента необходимо, чтобы:
CRex -(Rc + RJ + L > 0, (6)
или
_ т
R + R >-. (7)
с цн Cf>
вх
Учитывая, что значения напоров могут быть определены по известным зависимостям [7]: Нцн = H0 _ aQ2, Hc = bQ2, где #цн - напор, развиваемый центробежным насосом; Hс -
напор гидравлической сети, и, принимая соответственно значения дифференциальных со-
dH dH
противлений: Ri =-= _2aQ, R =-= 2bQ, получим следующее неравенство:
dQ dQ
2Q(a _ b) <-L-, (8)
или
L
а _ b <-. (9)
2QCRex " )
о HoZHe a H с
В то же время а = —~2-, b = —j , отсюда получим выражение
Q Q
е 2есЯ
или
Ь
Ис-Я -Ис Ь -7-1 , (10)
ицн+и > ио——■■ е- (11)
вх
В рабочей точке (пересечении НРХ) напоры равны между собой: + Ис = 2И, тогда из (11) получим условие:
2Я > Ис(12)
вх
или
Т_Г Т
н > Ис--Ь— д. (13)
2 4СЯвХ * ( )
Полученное условие позволяет построить границу колебательной устойчивости и выделить на НРХ области, соответствующие разной степени колебательности гидравлических процессов в насосном комплексе (рис. 5).
Рис. 5. Напорно-расходная характеристика центробежного насоса с выделенными областями различной степени колебательности (неустойчивые процессы, колебательные и апериодические процессы)
Условие положительности значений свободного члена характеристического уравнения передаточной функции (4) исследуемой системы:
Яс + Яцн + Явх > 0. (14)
си си
Учитывая, аналогично вышеизложенному, что =-= —2а6, Я =-= 2Ь6 , полу-
со с со
чим неравенство:
Явх > 20(Ь — а). (15)
И0 ~ ИЦН , _ —
& , Ь=е2 и И ЦН
Принимая, также аналогично, что а =- , Ь = —С и И„н + Ис = 2И, получим
Ь-а = 2—1—1, (16)
или
^ — я е
— + (17)
2 4
Полученное условие, дополнительно к предыдущему позволяет построить границу апериодической устойчивости и выделить на НРХ соответствующие области (рис. 5).
АНАЛИЗ СТЕПЕНИ КОЛЕБАТЕЛЬНОСТИ ДИНАМИЧЕСКИХ ПРОЦЕССОВ
Передаточная функция (4) может быть представлена в обобщенном виде передаточной функции колебательного звена:
Ж(Р) = Т2 2 ^ (18)
Т р + 2Е,Тр + к
где
Т2=ьс^
В . (19)
|2£Т = С^ • (Я + ЯцН) + Ь ( )
Используя данные выражения и учитывая, что условие 0 < Е, < 1 соответствует области колебательных, а 1 < Е области апериодических динамических процессов, а также применяя изложенный выше подход к построению соответствующих областей, можно получить следующее дополнительное условие апериодичности переходных процессов:
М 1
и >—+1 • е 2 2
' ГГ 1 Ь >
СЯвх 2 СДх ,
(20)
Соответствующая граница областей изображена на рис. 5.
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ
Для иллюстрации характера переходных процессов в насосном комплексе, характерных для выделяемых на НРХ областей, осуществлен численный эксперимент (расчетные параметры которого представлены в табл. 1) с использованием динамической модели в виде передаточной функции (4).
Таблица 1
Параметры для расчетов математической модели
Параметр I эксперимент II эксперимент III эксперимент
а 1,6
Ь 1,36 0,383 0,06
Q 4,5 5,5 6
ь 180
с 1,0039
Явх 13,4
Результаты численного моделирования процессов в рассматриваемой динамической системе в математическом пакете Ма^аЬ для различного расположения рабочей точки НРХ представлены на рис. 6.
Рис. 6. Графики переходных процессов при различных условиях моделирования:
а - в области устойчивых колебательных процессов; б - на границе устойчивости; в - в области апериодических (монотонных) процессов
ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ
Схема включения ЦН в НК определяет структуру динамической модели, выстраиваемую на основе метода электроаналогии, а также набор элементов, в том числе динамических, входящих в данную модель.
Полученные в процессе анализа динамической модели результаты позволяют выделить на НРХ НК области, характеризующиеся той или иной степенью устойчивости и колебательности динамических процессов при перекачивании жидкости. Это позволяет предоставить дополнительную информацию для обоснованного выбора прецедентов (с наименьшей степенью колебательности) при автоматизированном управлении НК.
б
а
в
Для определения границ соответствующих областей необходима предварительная идентификация параметров модели (Ь, С, Явх), в частности экспериментально-расчетным путем, а также с применением имитационного моделирования с использованием численных методов.
Доверительный интервал идентифицируемых параметров позволяет уточнить границу каждого из режимов:
н (21)
2 4
H> H _ 1
2 4
( L ^
-±Д
KCRx У
Q, (22)
H > HL + 1 2 2
L
lCR
±>/Д
Л
1Г L }
Q-1 — ±Д • Q, (23)
4
v CR*
где ±5 - поправочный коэффициент Rex для доверительного интервала; ±Д - поправочный коэффициент L для доверительного интервала. В этом случае границы областей будут
CR6X
иметь вид секторов, ширина которых определяется величиной соответствующих поправочных коэффициентов.
ЗАКЛЮЧЕНИЕ
Представленные результаты исследования позволяют определить условия возникновения нежелательных колебательных процессов в системах перекачивания жидкости на основе ЦН в зависимости от статических и динамических параметров НК.
Использованный для построения динамической модели подход на основе метода электроаналогии дал возможность осуществить анализ устойчивости и степени колебательности процесса перекачивания жидкости методами теории автоматического управления. Полученные аналитические соотношения позволили определить границы соответствующих режимных областей на НРХ, соответствующих той или иной степени устойчивости и колебательности динамических процессов.
Предложенный подход рассмотрен на примере одной из схем включения ЦН в НК, но может быть применен к другим возможным схемам с построением соответствующих моделей по методу электроаналогии.
Выявление соответствующих режимных областей на НРХ предоставляет дополнительную информацию для обоснованного выбора прецедентов при автоматизированном управлении НК.
СПИСОК ЛИТЕРАТУРЫ
1. Краев М. В., Лукин В. А., Овсянников Б. В. Малорасходные насосы авиационных и космических систем. М.: Машиностроение, 1985. 128 с. [ M. V. Kraev, V. A. Lukin, B. V. Ovsyannikov, Low-flow pumps of aviation and space systems, (in Russian). Moscow: Mashinostroenie, 1985. ]
2. Джарратано Дж., Райли Г. Экспертные системы: принципы разработки и программирование. 4-е изд.: пер. с англ. М.: ООО «И.Д. Вильямс», 2007. 1152 с.: ил. [ J. Jarratano, G. Riley, Expert systems: principles of development and programming, 4th ed.: translated from English. Moscow: OOO "I. D. Williams", 2007. ]
3. Высокооборотные лопаточные насосы / под ред. Б. В. Овсянникова и В. Ф. Чебаевского. М.: Машиностроение, 1975. 336 с. [ B. V. Ovsyannikov, V. F. Chebaevsky (eds.), High-speed vane pumps, (in Russian). Moscow: Mashinostroenie, 1975. ]
4. Li S. C. Cavitation of Hydraulic Machinery. Imperial College Press, 2000. 492 p. [ S. C. Li, Cavitation of Hydraulic Machinery. Imperial College Press, 2000. ]
5. Певзнер Л. Д. Теория систем управления. СПб.: Лань, 2013. 424 с. [ L. D. Pevzner, Control systems theory, (in Russian). St. Petersburg: Lan, 2013. ]
6. Лютов А. Г., Новоженин М. Б. Методология автоматизированного контроля и управления режимами работы насосного комплекса в условиях возникновения кавитации // Мехатроника, автоматизация, управление. 2021. Т. 22, № 9.
C. 468-474. DOI: 10.17587/mau.22.468-474. [ A. G. Lutov, M. B. Novozhenin, "Methodology of automated monitoring and control of pumping complex operating modes in conditions of cavitation", (in Russian), in Mehatronika, avtomatizacia, upravlenie, vol. 22, no. 9, pp. 468-474, 2021. ]
7. Goppelt F., Hieninger T., Schmidt-Vollus R. Modeling Centrifugal Pump Systems from a System-Theoretical Point of View // IEEE 18th Mechatronika 2018 (Brno, Tschechien, December 5-7 2018), 2018. Pp. 1-8. [ F. Goppelt, T. Hieninger, R. Schmidt-Vollus, "Modeling Centrifugal Pump Systems from a System-Theoretical Point of View", in IEEE 18th Mechatronika 2018, Brno, Tschechien, 2018, pp. 1-8. ]
8. Nollau R. Modellbildung und Simulation technischer Systeme. Berlin Heidelberg: Springer, 2009. Pp. 179-186. [ R. Nollau, Modellbildung und Simulation technischer Systeme. Berlin Heidelberg: Springer, 2009. ]
ОБ АВТОРАХ
ЛЮТОВ Алексей Германович, проф. каф. автоматизации технол. процессов. Дипл. инж. электрон. техн. (УАИ, 1985). Д-р техн. наук по сист. анализу, упр. и обр. информ. (УГАТУ, 2005). Иссл. в обл. управления слож. техн. объектами.
НОВОЖЕНИН Максим Борисович, ст. преп. каф. автоматических систем. Дипл. инженер (УГАТУ, 2011). Канд. техн. наук по автоматизации технол. процессов (ОГУ, 2020). Иссл. в обл. управления слож. техн. объектами.
ОГОРОДОВ Владимир Алексеевич, доц. каф. автоматизации технол. процессов. Дипл. инж.-мех. (УАИ, 1975). Канд. техн. наук по технол. произв. ЛА (КАИ, 1985). Иссл. в обл. автоматиз. техн. процессов.
METADATA
Title: Research of the conditions for the occurrence of vibratory processes in the automated control of pumping complexes based on precedents.
Authors: A. G. Lutov1, M. B. Novozhenin2, V. A. Ogorodov3
Affiliation:
1 3 Ufa State Aviation Technical University (UGATU), Russia.
2 Moscow State Technical University of Radio Engineering, Electronics and Automation - Russian Technological University (MIREA - Russian Technological University), Russia.
Email: 1 lutovl@mail.ru, 2 novozhenin.maxim@yandex.ru, 3 ogorodov.v@ugatu.su
Language: Russian.
Source: Vestnik UGATU (scientific journal of Ufa State Aviation Technical University), vol. 25, no. 4 (94), pp. 101-110, 2021. ISSN 2225-2789 (Online), ISSN 1992-6502 (Print).
Abstract: The precedent approach used for the automated control of operating modes of pumping complexes makes it possible to reduce or completely eliminate various unsteady hydraulic processes, including cavitation. The article discusses the conditions for the occurrence of oscillations of a hydrodynamic nature in pumping complexes, contributing to the emergence or activation of these undesirable processes, with the construction of models based on the method of electrical analogy and the application of stability criteria of the theory of automatic control. The results obtained allow the formation and use of a base of precedents when controlling the operating modes of pumping complexes, taking into account the oscillation factor.
Key words: electric centrifugal pump; oscillatory processes; pressure-flow characteristic; method of electrical analogy.
About authors:
LUTOV, Alexey Germanovich, Prof. Head of Dept. of Automation of Technological Processes. Dipl. Electronics Engineer (USATU, 1985). Dr. of Tech. Sci. (USATU, 2005).
NOVOZHENIN, Maxim Borisovich, Senior Teacher of Dept. of Automated Systems. Dipl. Engineer (USATU, 2011). Cand. of Tech. Sci. (OGU, 2020).
OGORODOV, Vladimir Alekseevich, Prof. Dept. of Aviation Technology. Dipl. engineer-mechanic (UAI, 1975). Cand. of Tech. Sci. (KAI, 1985).