Belgorod, Belgorod, Russian Federation
Belgorod, Belgorod, Russian Federation
Belgorod, Russian Federation
Issleduyutsya dinamicheskie rezhimy i bifurkacii v impul'snoy sisteme upravleniya nagrevatel'noy ustanovkoy, sostoyanie kotoroy opisyvaetsya differencial'nymi uravneniyami s razryvnymi pravymi chastyami. Pokazano, chto rassmatrivaemaya sistema mozhet demonstrirovat' chrezvychayno bol'shoe mnogoobrazie nelineynyh yavleniy i bifurkacionnyh perehodov, takih kak, kvaziperiodichnost', mul'tistabil'noe povedenie, haotizaciya kolebaniy cherez klassicheskiy kaskad bifurkaciy udvoeniya perioda i bifurkacii granichnogo stolknoveniya.
nagrevatel'naya ustanovka, teplovoy ob'ekt, teploprovodnost', tigel', drobnyy poryadok, zakon upravleniya, polevoy tranzistor, shirotno-impul'snaya modulyaciya, bifurkaciya.
Введение. Технология выращивания монокристаллов представляет собой процесс управляемой кристаллизации, при котором качество растущего кристалла определяется точностью управления условиями фазовых переходов [1]. При выращивании кристалла синтетического сапфира необходимо обеспечить закон изменения температуры в тигле от 25 °С до 2050 °С с определённой степенью наращивания и спада температуры, что предполагает применение автоматической системы управления с возможностью программного задания изменения температуры в тигле с заданной точностью.
Теплотехнический объект, нагревательная установка, состоит из следующих зон: внутреннего печного пространства 1, заполненного воздухом или газом; нихромового электронагревателя 2, равномерно распределённого во внутреннем слое футеровки 3, состоящей из магнезитового кирпича и внешнего слоя футеровки 4 из минеральной ваты в цилиндрическом стакане из оцинкованной стали (рис. 1). Геометрическая форма печи – ограниченный цилиндр, сверху и снизу которого располагается футеровка.
Для решения задачи синтеза закона управления классическим методом аппроксимации с использованием свободно распространяемой библиотеки FOMCON по экспериментальной кривой разгона теплового объекта определена передаточная функция нагревательной установки следующего вида:
(1)
где – коэффициент передачи объекта, T1,T2 – постоянные времени объекта [2, 3].
Рис. 1. Нагревательная установка
Используемые в настоящее время регуляторы температуры с тиристорными преобразователями существенно искажают форму кривой тока, потребляемого из сети, приводя к возникновению в питающей сети несинусоидальных режимов.
Для устранения указанных недостатков на базе патента авторов [4] разработана и реализована система управления нагревателем высокой мощности, построенная на основе высокочастотного преобразователя электрической энергии с широтно-импульсным регулированием. Повышение энергетических показателей с упрощением управления технологическим объектом достигается за счет использования в качестве ключевых элементов преобразователя полевых транзисторов с применением дробных законов управления широтно-импульсной модуляцией, улучшающих качество системы [5–7].
Однако в нелинейных импульсных системах при вариации параметров объекта управления, а также воздействии внешних возмущений возможно возникновение сложных нелинейных явлений, включая колебания на пониженных частотах, кратных частоте модуляции, квазипериодические и хаотические режимы [8, 9].
Следствием этого является многократное увеличение амплитуды колебаний температуры нагревательной установки, снижение точности регулирование и нарушение хода технологического процесса.
Целью данной работы является численное исследование бифуркационных явлений в динамике импульсной системы управления нагревательной установкой.
- Постановка и аналитическое решение
задачи
Уравнение движения системы управления нагревательной установкой, непрерывная линейная которой описывается передаточной функцией (1), имеет вид
(2)
где – температура в нагревательной установке; , – сигналы на входе и выходе широтно-импульсного модулятора, соответственно; – коэффициент передачи непрерывной линейной части системы; T1, T2 – постоянные времени.
Введем , и перепишем уравнение движения (2) в нормальной форме Коши:
(3)
Выходной сигнал модулятора
(4)
где T0 – период модуляции, τk – ширина импульса, определяемая видом импульсной модуляции. В работе рассматривается система с широтно-импульсной модуляцией первого рода (ШИМ-1) и пропорциональным корректирующим звеном в цепи обратной связи. Тогда входной сигнал модулятора определяется выражением
Здесь Vref – сигнал задания температуры нагревательной установки, β – коэффициент передачи датчика температуры, α – коэффициент усиления.
При ШИМ-1 величина τk находится как:
где Vs – опорный сигнал модулятора.
Введем обозначения
, ,
. (5)
Уравнения движения (3) примут вид:
(6)
Исследование динамической системы (6) можно свести к изучению свойств двумерного кусочно-гладкого отображения:
, ,
, .
Здесь ширина импульса τk согласно (5) определяется
где – матрица - строка.
В исследованиях были выбраны следующие значения параметров модели:
= 10240 c2; = 352 c; T0 = 10 с; K = 327.8 °C/(B·c); U0= 24 В – напряжение питания; β = 0.01 B/°C; Vs = 5 B; α > 0; Vref = 5 B.
Период T периодического движения динамической системы (6) в общем случае является кратным периоду внешнего воздействия T0: T = mT0, m = 1, 2, . . . . Движение с таким периодом будем называть m-циклом или циклом периода m.
2. Бифуркационный анализ
При проведении бифуркационного анализа в качестве варьируемых параметров были выбраны напряжение питания U0 и коэффициент усиления α. На рис. 2 приведены однопараметрические бифуркационные диаграммы, рассчитанные для разных значений U0 при изменении коэффициент усиления α.
Рис. 2. Бифуркационные диаграммы при различных параметрах U0 и T0, Vs = 5 B, Vref = 5 B
(соответствует уставке 500 °С)
При малых значениях U0 система демонстрирует квазипериодическое поведение с ярко выраженной мультистабильностью. На рис. 2, а изображена бифуркационная диаграмма, иллюстрирующая рождение замкнутой инвариантной кривой, отвечающей двухчастотному квазипериодическому режиму. Как следует из рис. 2, а, при увеличении коэффициент усиления α 1 – цикл теряет устойчивость через бифуркацию Неймарка - Саккера. Потеря устойчивости приводит к возникновению устойчивой замкнутой инвариантной кривой, при этом 1 – цикл продолжает существовать, но становится неустойчивым фокусом. Как известно, характер движения на замкнутой инвариантной кривой определяется числом вращения, когда оно иррационально, точки отображения плотно заполняют инвариантную кривую и динамика становится квазипериодической.
При рациональном числе вращения на инвариантной кривой имеется четное число периодических орбит, половина из которых устойчивые, а половина – седловые, а сама инвариантная кривая образована замыканием неустойчивых многообразий седловых циклов. На рис. 2, а окно с периодической динамикой отвечает области устойчивости резонансного 4 – цикла (области существования замкнутой инвариантной кривой с числом вращения 1:4). При увеличении α резонансный 4 – цикл претерпевает каскад бифуркаций удвоения периода, завершающийся хаотизацией колебаний, при этом замкнутая инвариантная кривая разрушается. Численные эксперименты показали, что с увеличением U0 область устойчивости 1 – цикла уменьшается.
Бифуркационная диаграмма, изображенная на рис. 2, б, иллюстрирует типичный сценарий рождения сосуществующих аттракторов. При изменении коэффициента усиления α жестко возникает устойчивый 3 – цикл. При дальнейшем увеличении α реализуется бесконечный каскад бифуркаций удвоения периода, завершающийся хаотизацией колебаний. По этой причине в широком диапазоне изменения параметров наряду с устойчивым 1 – циклом существуют либо устойчивые периодические колебания, либо хаотические режимы. То есть в зависимости от начальных условий может устанавливаться или периодическое, или хаотическое движение.
На рис. 2, в показан случай субкритического учетверения периода 1 – цикла через так называемую бифуркацию граничного столкновения («border - collision bifurсation», см., например, [8–15]). На рис. 2, г представлен пример рождения 4-х полосного хаотического аттрактора (four-band chaotic attractor) через border-collision flip bifurсation [8, 15].
Заключение. В статье представлены результаты бифуркационного анализа широтно-импульсной системы управления нагревательной установкой.
Выявлено, что при малых значениях напряжения питания системы демонстрирует квазипериодическое поведение с ярко выраженной мультистабильностью [11], но при этом система имеет достаточно большой запас устойчивости по коэффициенту усиления.
Однако, при увеличении напряжения питания область устойчивости 1 – цикла (рабочего режима) сужается, и потеря устойчивости происходит через бифуркацию граничного столкновения, приводящая к внезапной хаотизации колебаний.
*Работа выполнена в рамках Программы развития опорного университета на базе БГТУ им. В.Г. Шухова.
Авторы выражают благодарность проф. Жусубалиеву Ж.Т. за обсуждение результатов исследований и полезные комментарии.
1. Lodiz R., Parker R. Rost monokristallov. M.: Mir, 1974. 540 s.
2. Lykov A.V. Teoriya teploprovodnosti. M.: Vysshaya shkola, 1967. 600 s.
3. Fractional-order Modeling and Control. [Elektronnyy resurs]. URL: http://fomcon.net/ (data obrascheniya: 20.07.2017).
4. Pat. № 2612311 Rossiyskaya Federaciya, MPK G05D 23/22. Ustroystvo regulirovaniya temperatury elektronagreva / Gol'cov Yu. A., Zhusubaliev Zh. T., Kizhuk A. S., Kolenchenko V. V., Rubanov V. G., zayavitel' i patentoobladatel' BGTU im. V.G.Shuhova. - № 2016113209, zayavl. 06.04.2016, opubl. 06.03.2017. Byul. № 7. 5 s.
5. Gol'tsov Yu.A., Kizhuk A.S., Rubanov V.G. Control of high power thermal object in the class of fractional order regulators // International Journal of Pharmacy & Technology. 2016. T. 8. №.4. S. 24790-24800.
6. Rubanov V.G., Kizhuk A.S, Gol'cov Yu.A., Karikov E.B. Realizaciya algoritma approksimacii drobnogo integrodifferencirovaniya s ocenkoy oshibki // Vestnik BGTU im. V.G. Shuhova. 2015. № 2. S. 148-151.
7. Kizhuk A.S., Gol'cov Yu.A. Mikroprocessornaya sistema avtomaticheskogo upravleniya teplovym rezhimom tehnologicheskogo processa vyraschivaniya kristalla sapfira // Pribory i sistemy. Upravlenie, kontrol', diagnostika. 2014. №11. S. 42-49.
8. Moschnyy polevoy tranzistor IRF3205. Tehnicheskaya dokumentaciya. [Elektronnyy resurs]. URL: www.irf.com (data obrascheniya: 20.07.2017).
9. Zhusubaliyev Zh.T., Mosekilde E. Bifurcations and Chaos in Piecewise-Smooth Dynamical Systems. Singapore: World Scientific, 2003. pp: 363.
10. Banerjee S., Verghese G.C. Nonlinear Phenomena in Power Electronics. New York, USA: IEEE Press, 2001.
11. Bernardo Di.M., Feigin M. I., Hogan S. J., Homer M. E. Local Analysis of C -bifurcations in n-dimensional Piecewise-Smooth Dynamical Systems // Chaos, Solitons and Fractals. 1999. 10(1). Pp. 1881-1908.
12. Nusse H.E., Yorke J. A. Border-Collision Bifurcations Including “Period Two to Period Three” for Piecewise Smooth Systems // Physica D. 1992. 57. Pp. 39-57.
13. Banerjee S., Ranjan P., Grebogi C. Bifurcations in Two-Dimensional Piecewise Smooth Maps - Theory and Applications in Switching Circuits // IEEE Trans. Circ. Syst. I. 2000. 47(5). Pp. 633-643.
14. Zhusubaliyev Zh.T., Soukhoterin E.A., Mosekilde E. Border-Collision Bifurcations and Chaotic Oscillations in a Piecewise-Smooth Dynamical System // Int. J. Bifurcation Chaos. 2001. 11(12). Pp. 1193-1231.
15. Bernardo Di.M., Budd C.J., Champneys A.R., Kowalczyk P. Piecewise-smooth Dynamical Systems: Theory and Applications, in: Applied Mathematical Sciences. Springer. 2008. Vol. 163. Pp. 483.
16. Zhusubaliyev Zh.T., Mosekilde E. Multi-stability and Hidden Attractors in a Multilevel DC/DC Converter // Mathematics and Computers in Simulation.