УДК 681.511
АНАЛИЗ ДИНАМИЧЕСКИХ ПРОЦЕССОВ В СИСТЕМЕ СТАБИЛИЗАЦИИ ТЕМПЕРАТУРЫ В ИНКУБАЦИОННОЙ КАМЕРЕ
А.А.Варенов, В.П. Зенцов, Н.А. Малёв
Казанский государственный энергетический университет, Казань, Россия
Резюме: Исследованы тепловые процессы в инкубационной камере, разработаны рекомендации по управлению температурой в этой прецизионной системе на основе синтеза корректирующих устройств, позволяющих повысить качество её функционирования.
Ключевые слова: функциональная схема, передаточная функция, процесс регулирования, моделирование линейной непрерывной системы.
DOI:10.30724/1998-9903-2018-20-9-10-101-108
Для цитирования: А.А.Варенов, В.П. Зенцов, Н.А. Малёв Анализ динамических процессов в системе стабилизации температуры в инкубационной камере // Известия высших учебных заведений. ПРОБЛЕМЫ ЭНЕРГЕТИКИ. 2018. Т. 20. № 9-10. С. 101-108. DOI:10.30724/1998-9903-2018-20-9-10-101-108.
DYNAMIC PROCESSES ANALYSIS IN THE TEMPERATURES STABILIZATION SYSTEM, IN THE INCUBATIOM CHAMBER
A.A. Varenov, V.P. Zentsov, N.A. Malyov
Kazan State Power Engineering University
Abstract: Thermal processes in incubation chamber were investigated. From the investigations, the recommendations were developed for controlling the temperature in this precision system based on the synthesis of correcting devices, allowing one to improve its functional quality.
Keywords: Functional scheme, transfer function, process regulation, modeling of linear continuous system.
For citation: A.A. Varenov, V.P. Zentsov, N.A. Malyov Dynamic processes analysis in the temperatures stabilization system, in the incubatiom chamber // Proceedings of the higher educational institutions. ENERGY SECTOR PROBLEMS 2018. vol. 20. № 9-10. pp. 101-108. DOI: 10.30724/1998-9903-2018-20-9-10-101-108.
Поддержание постоянной температуры в замкнутом объеме является важной задачей во многих производственных процессах. Теплообмен при этом характеризуется значительной инерционностью, т.к. для мгновенного изменения температуры необходимо затрачивать бесконечно большую энергию [1-3]. Распространенным источником тепловой энергии является электронагреватель, преобразующий электрическую энергию в тепловую.
На рис. 1 показана функциональная схема системы автоматического регулирования (САР) температуры в инкубационной камере (ИК) [6].
Исследуемая САР относится к разряду прецизионных систем. Динамические процессы в этой системе описываются уравнениями: сравнивающее устройство А—(/) = —0 - Е(/);
усилительное устройство и(?) = К АU (t); тиристорный блок питания итр (?) = Ктриу (^); нагревательный элемент Гн@н (/) + ©н (/) = Л',,/',, (/) + ©к (/); инкубационная камера Гк©к (/) + ©к (/) = А',,©,, (/) + А0©о (/); термопара '/',/•'(/) + /■'(') = А"т@к (/), где ©н , ©к, ©0 - температуры нагревателя, камеры и окружающей среды соответственно;
и 2
Е - эдс термопары; рн = —н - мощность, потребляемая нагревателем (ин - напряжение,
К
прикладываемое к нагревателю, а К - его сопротивление); к., , К™, Кн =—1—,
У тр н ТТ О
К тн °И
К — КБ
К =-IИ-И-, К =-то_о-' К - коэффициенты передачи усилителя, блока
КтнБн + Кто —о Ктн —н + Кто —о
управления тиристорами, нагревателя, инкубационной камеры по управляющему воздействию —н, инкубационной камеры по возмущающему воздействию ©о и термопары
С С т С т соответственно; т„ = —н ^ , т =--, т = —I—т. - постоянные времени
н ^ О К у- о , тт о т ТТ о
Ктн—н Ктн—н ^ Кто—о Ктт—т
нагревателя, камеры и термопары соответственно; Сн, С , С - удельные теплоемкости; тИ, т, тт - массы; ктн, Кто, Ктт- коэффициенты теплоотдачи нагревателя, внутренних стен камеры и термопары; £ , 5 , 5 - площади поверхности этих элементов.
Рис. 1. Функциональная схема САР температуры в инкубационной камере: 1 - ИК; 2 - термопара; 3 - лоток для яиц; 4 - нагреватель; 5, 6 - тиристорный блок питания; 7 - блок управления тиристорами; 8 - усилитель; 9 - задающее устройство
После записи приведенных уравнений в изображении по Лапласу с использованием источников [3, 4, 5, 8] запишем передаточные функции элементов САР:
Жу (5 ) = —Щ = К у- усилитель;
- блок управления тиристорами;
тр( ) иу (5) *
ш (*) = ®н(5) = Кн - нагреватель по управлению; ну ( ) Рн (5) +1
Ш (5) = —н (51 =_1_ - нагреватель по возмущению;
Ш (5) = Е(5) = Кт - термопара;
А ' ©к (5) Тт 5 +1
ш (5) = —к ) =_К_- ИК по управлению;
Шку)=©н (5) = Тк5 +1 у р
© ($) К ш (5) = —=-0— - ИК по возмущению.
кв( ) —(5) Тк5 +1
и 2
Мощность нагревателя Рн = является нелинейной функцией напряжения.
К
Учитывая, что изменения напряжения ин в процессе регулирования небольшие, проведена линеаризация этой функции разложением в ряд:
чСин Г н но/ 2!
С 2 Рн
V Си2 ,
1!
откуда
СР СР
Рн (ин )-Рн (и0 ) = ^ Аин или ЛРН = ^ .
С
Учитывая СРн
и 2 ^ К
син син К
2ин , при и = Цп получим
СРн _ 2ин0 _ тт -—-— К „ .
син К р
2и
Таким образом, АР =-— Аин = АР А^н является уравнением в отклонениях
К р
переменных, в котором удержаны только составляющие первого порядка малости и опущены составляющие более высокого порядка, что является допустимым для инженерных расчетов. Передаточные функции элементов САР позволяют построить структурную схему динамической модели (ССДМ) [2], показанную на рис. 2. Свернув внутренний замкнутый контур с положительной обратной связью [9], получим передаточную функцию разомкнутой некоррелированной системы:
Ку Ктр Кр Кн (тн * +1) Кк Кт
Ш ( 5 ) =
(Тн * +1)[[ * + 1)Т * +1)-Кк ] К _ В ( 5 )
[(Тк * +1 )(ТН* +1)-Кк ] П ( 5)'
где К = КуКтрКрКнКкКт - коэффициент передачи разомкнутой некоррелированной системы.
где
Рис. 2. Структурная схема динамической системы модели Характеристический многочлен замкнутой системы 4(5) = В(5) + В(5) или
А(5) = ТнТкТт 53 + (ТнТк + ТнТт + ТкТт 2 +
+ (тн + Тк - ТтТк )5 + (1 + К) = а353 + а252 + ах5 + а0, а = ТТТ ^ 0,
н к т
a-
= (ТИТк+ТИТт+ТкТт)>~0
«1 =тп+тк-тптк>0
а0 = 1 + К ^ 0.
Оценка устойчивости САР по критерию Гурвица (А2 = а^а2 — а0а3 ^ 0) для
физически реализуемых параметров показывает, что система является устойчивой и обладает большими запасами устойчивости, что подтверждается графиками логарифмических частотных характеристик (ЛЧХ), показанными на рис. 3. Из полученных ЛЧХ следует, что запас устойчивости по амплитуде составляет 95,3 дБ, тогда как запас устойчивости по фазе сложно оценить, в силу сравнительно небольшого значения общего коэффициента передачи системы на порядок меньшего единицы и, как следствие, отсутствия точки пересечения ЛАЧХ с осью 0 дБ.
Frequency (rad/s)
Рис. 3. ЛЧХ САР температуры
ТI
Рис. 4. Схема моделирования САР температуры в среде МайаЪБтиЫпк
В связи с этим научная задача, решаемая в работе, сводится к анализу временных характеристик системы регулирования температуры в инкубационной камере с последующей выработкой рекомендаций по применению соответствующих способов коррекции, позволяющих улучшить показатели качества переходных процессов.
Результат моделирования некоррелированной системы в соответствии со схемой моделирования в среде МаЛаЪБШиНпк [7] (рис. 4) при заданных параметрах показан на рис. 5. Из него следует, что температура в ИК достигает своего установившегося значения 38,2 °С за 3-104с.
е„ 5с
Рис. 5. Переходная характеристика системы по задающему воздействию
Из графика видно, что в системе имеет место ошибка по возмущению и при температуре окружающей среды ©0 = -40 °С абсолютное значение ошибки составляет 1,7 °С. Кроме того, САР температуры обладает существенной инерционностью и малым быстродействием, что также является весомым недостатком некоррелированной системы.
Моделирование нелинейной САР температуры в инкубаторе показало, что графики переходных процессов не имеют принципиальных отличий по сравнению с
линеаризованной системой, что говорит о корректности принятого выше допущения. На рис. 6 показана зависимость мощности нагревателя от времени.
Р(Вт)
2*1 I I I I I I-1 ■
Рис.6. График зависимости мощности нагревателя от времени
Из графика следует, что в момент запуска системы потребляемая на нагрев мощность составляет 27,83 Вт, постепенно устанавливаясь к значению 25,5 Вт. При воздействии возмущения в виде изменения температуры окружающей среды наблюдается некоторое увеличение мощности до величины 25,62 Вт.
Таким образом, исследуемая система регулирования температуры в инкубаторе является чувствительной к возмущающим воздействиям, а время достижения установившегося состояния составляет около 8 часов. С учетом перечисленных особенностей некоррелированной САР возникает задача синтеза коррелирующих устройств, позволяющих повысить качество функционирования системы.
Литература
1. Бахмуров А.В., Гильфанов К.Х., Гольбрайх Л.Я. Совершенствование алгоритма управления объектом теплоэнергетики на основе интеллектуального контроллера // Энергетика Татарстана. 2013. С.42-46.
2. Варенов А.А., Якимов Н.Д., Ерашова Ю.Н. Инженерная модель теплового режима в индивидуальном доме // Известия вузов. Проблемы энергетики. 2016. № 3-4. С.22-28.
3. Варенов А.А. Структура и динамическая модель электропривода с импульсным управлением для беспилотного летательного аппарата мониторинга энергообъектов // Известия вузов. Проблемы энергетики. 2013. № 5-6. С.78-83.
4. Варенов А.А., Малев Н.А. Гибридный регулятор напряжения релейного типа в системе бесперебойного электроснабжения // Известия вузов. Проблемы энергетики. 2013. № 7-8. С.47-53.
5. Варенов А.А., Малев Н.А., Погодицкий О.В. Аналоговый и цифровой регуляторы автоматической системы стабилизации частоты // Известия вузов. Проблемы энергетики. 2014. № 5-6. С. 120-129.
6. Математическое моделирование систем автоматического регулирования: Практикум по автоматике / Под ред. Карташова Б.А. М.: "Колосс", 2006.
7. Гайдук А.Р., Беляев В.Е., Пьявченко Т.А. ТАУ в примерах и задачах с решениями в MATLAB. Изд. Лань, 2017. 464с.
8. WONGH.Y. Heat transfer for engineer. LONGMAN, 1977, p.166.WONG, H.Y., 1977. Heat transfer for engineer. LONDON.: LONGMAN.
9. Mehta B.R.,Y. 26 November 2014. Jaganmohan Reddy. Industrial Process Automation Sistemsist Edition. Design and Implementation..Butterwrth-Heinemann.
10. Mohd Badli Ramli, Mohd Faiz, Mohd Zin. Egg Hatching Incubator Using Conveyor Rotation System// www.sciencedirect.com URL:https://doi.org/10.1016/j.promfg.2015/07/091.
Авторы публикации:
Вареное Александр Андреевич - профессор кафедры «Теоретические основы электротехники» (ТОЭ) Казанского государственного энергетического университета (КГЭУ). E-mail: [email protected].
Зенцое Василий Павлович - канд. физ.-мат. наук, доцент кафедры «Теоретические основы электротехники» (ТОЭ) Казанского государственного энергетического университета (КГЭУ). E-mail: [email protected].
Малёе Николай Анатольевич - канд. техн. наук, доцент кафедры «Теоретические основы электротехники» (ТОЭ) Казанского государственного энергетического университета (КГЭУ). E-mail: [email protected].
References
1. Bakhmurov A.V., Gilfanov K.H., Golbraikh L.Ya. Perfection of the algorithm of control of the object of heat power engineering on the basis of an intelligent controller. Power engineering of Tatarstan. 2013. P.42-46.
2. Varenov A.A, Yakimov N.D, Yerashova Yu.N. Engineering model of the thermal regime in the individual house. //Известия HIGH SCHOOL. Problems of energy. 2016. № 3-4. P.22-28.
3. Varenov A.A. Structure and dynamic model of an electric drive with pulse control for an unmanned aircraft for monitoring energy facilities.// Известия HIGH SCHOOL. Problems of energy, 2013. № 5-6. P.78-83.
4. Varenov A.A, Malev N.A. Hybrid voltage regulator of relay type in the uninterruptible power supply system // Известия HIGH SCHOOL. Problems of Power Engineering, 2013. No. 7-8. P.47-53.
5. Varenov A.A, Malev N.A, Pogoditsky O.V. Analog and digital regulators of automatic frequency stabilization system. //Известия HIGH SCHOOL. Problems of Power Engineering, 2014. № 5-6. P.120-129.
6. Ed. Kartashov B.A. Workshop on automation. Mathematical modeling of automatic control systems. M .: "The Colossus", 2006.
7. Gaiduk A.R, Belyaev V.E, Piavchenko T.A. TAU in examples and problems with solutions in MATLAB. Ed. Lan. 2017. 464p.
8. WONG H.Y. Heat transfer for engineer. LONGMAN, 1977, p.166. WONG, H.Y., 1977. Heat transfer for engineer. LONDON .: LONGMAN.
9. Mehta B.R., Y. 26 November 2014. Jaganmohan Reddy. Industrial Process Automation Sistems ist Edition. Design and Implementation .. Butterwrth-Heinemann.
8. WONGH.Y. Heat transfer for engineer. LONGMAN, 1977, p.166.WONG, H.Y., 1977. Heat transfer for engineer. LONDON.: LONGMAN.
9. Mehta B.R.,Y. 26 November 2014. Jaganmohan Reddy. Industrial Process Automation Sistemsist Edition. Design and Implementation..Butterwrth-Heinemann.
10. Mohd Badli Ramli, Mohd Faiz, Mohd Zin. Egg Hatching Incubator Using Conveyor Rotation System// www.sciencedirect.com URL:https://doi.org/10.1016/j.promfg.2015/07/091.
Authors of the publication
Alexander A. Varenov - Professor, Department of Theoretical Foundations of Electrical Engineering, Kazan State Power Engineering University (KSPU). E-mail: [email protected].
Vasiliy P. Zentsov - Ph.D. (physical&mathematical), Associate Professor, Department of Theoretical Foundations of Electrical Engineering, Kazan State Power Engineering University (KSPU), E-mail: [email protected].
Nikolay A Malev - Ph.D. (Techn.), Associate Professor, Instrument Making and Mechatronics Department, Kazan State Power Engineering University (KSPU), E-mail: [email protected].
Поступила в редакцию 24 апреля 2018 г.