Динамическая модель отопительных приборов и систем Текст научной статьи по специальности «Строительство и архитектура»

УДК 681.513.6

DOI: 10.14529/ctcr150208

ДИНАМИЧЕСКАЯ МОДЕЛЬ ОТОПИТЕЛЬНЫХ ПРИБОРОВ И СИСТЕМ

В.И. Панферов, С.В. Панферов

Южно-Уральский государственный университет, г. Челябинск

Рассматривается задача построения динамической математической модели отопительных приборов и систем. Система отопления при этом представляется эквивалентным отопительным прибором. Выходной величиной модели в данном случае является мощность (тепловая производительность). Решение выполнено методом «серого ящика», т. е. структура модели определена посредством анализа физики процессов, протекающих в отопительных приборах и системах, а численные значения ее параметров найдены по экспериментальным данным методом параметрической идентификации. Полученное дифференциальное уравнение является линейным относительно искомой величины - мощности отопительного прибора (системы отопления), входные же величины - температура воды на входе системы отопления и температура внутреннего воздуха включены в это уравнение линейным образом, а расход теплоносителя - нелинейным образом. Показано, что при постоянном расходе теплоносителя модель может быть представлена с помощью передаточных функций в изображениях Лапласа: передаточная функция по каналу «температура теплоносителя на входе в отопительный прибор (в систему отопления) - тепловая мощность отопительного прибора (системы отопления)» представляет собой инерционное звено первого порядка, а по каналу «температура внутреннего воздуха - тепловая мощность отопительного прибора (системы отопления)» - инерционно-форсирующее звено с отрицательным коэффициентом передачи. Анализируются качественные особенности параметров модели, в частности, их зависимость от расхода теплоносителя. Приводится алгоритмическая структурная схема объекта управления, который включает систему отопления и собственно тепловой режим здания. Приводятся результаты параметрической идентификации модели по экспериментальным данным. Модель может быть использована при разработке интеллектуальных систем управления отоплением зданий.

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

Введение. Известно, что в нашей стране на теплоснабжение зданий расходуется более одной третьей всего добываемого топлива, причем основными при этом являются затраты на отопление. В связи с этим, вполне естественно, возникает задача разработки современных подходов и решений по рациональному использованию расходуемых на отопление энергетических ресурсов. Причем при решении этой проблемы следует иметь в виду, что наибольшая экономия тепловой энергии достигается за счет автоматизации систем отопления. Поэтому разработка вопросов, непосредственно связанных с автоматизацией систем отопления и их совершенствованием, является одной из актуальнейших проблем настоящего времени.

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

Синтез структуры математической модели. Рассмотрим сначала отдельный отопительный прибор. Известно, что теплота, доставляемая теплоносителем отопительному прибору, может быть вычислена по следующей формуле:

Q = Cm ■ Gm (Îbx "'вЫХ^ (1)

где Ст - удельная теплоемкость теплоносителя, Gm - массовый расход теплоносителя через прибор, tВХ и tВЫХ - соответственно температура теплоносителя на входе и выходе из отопительного прибора.

Часть этого количества теплоты kПР • ^¡^(^ - будет передана отопительным прибором внутреннему воздуху помещения, а оставшаяся часть Ст • Gт - ) - kПР • ^ПР (^ - tВ) пойдёт на изменение температуры теплоносителя в приборе и самого отопительного прибора. Здесь kПР - коэффициент теплопередачи отопительного прибора, ^ПР - его площадь поверхности теплообмена, ^ - средняя температура теплоносителя в отопительном приборе (или, что то, же самое - средняя температура металла отопительного прибора), ^ - температура воздуха в помещении. Поэтому, следуя [1], уравнение теплового баланса запишем следующим образом:

СПР ~Г = Ст ■ Gm (гВХ - ^ЫХ ) - кПР • ^ПР - tВ) , (2)

ат

где СПР = Ст • рт • ¥ПР + См • тм - полная теплоемкость отопительного прибора вместе с находящимся в нем теплоносителем, рт - плотность теплоносителя, ¥ПР - объем теплоносителя в отопительном приборе, См - удельная теплоемкость металла отопительного прибора, тм - масса отопительного прибора, т - время.

Известно, что средняя температура теплоносителя в отопительном приборе может быть выражена следующим образом:

^ = + ^ЫХ)/2, (3)

поэтому ^ЫХ = 2tТ - tВХ. Подставляя данное соотношение в (2) и переписав его соответствующим образом, получим, что структура динамической математической модели отопительного прибора будет иметь вид:

_СПР_ ^Т , + _2Ст • Gm_ , ,__кПР • ^Пр_ + ¡л\

----+ ^ =--tвx +---tв . (4)

2Ст • Gm + кПР • ^ПР ат 2Ст • От + кПР • ¥ПР 2Ст • Gm + kПР • ^11Р

Перепишем данное уравнение, используя стандартные для теории автоматического управления обозначения [2]:

^ = к^вх + k2tВ

температура теплоносителя в отопительном приборе», k2 =-———- - коэффициент

Т Т + ^ = к4вх + ^В. (5)

ат

Здесь Т=-ПР- - постоянная времени, к1 =-т-т--коэффициент

2Ст • Gm + кПР • ^ПР 2Ст • Gm + кПР • ^ПР

передачи по каналу «температура теплоносителя на входе в отопительный прибор - средняя

£Пи • К-,

2 2Ст • От + кПР • ^ПР

передачи по каналу «температура внутреннего воздуха - средняя температура теплоносителя в отопительном приборе».

Как видно из приведенных соотношений, с увеличением расхода теплоносителя через прибор Gm его постоянная времени Т уменьшается, следовательно, уменьшается и инерционность переходных процессов в приборе, одновременно уменьшается и коэффициент передачи к2, т. е. уменьшается влияние температуры внутреннего воздуха на среднюю температуру отопи-

дкл 2Ст • к™ • ^р

\2

тельного прибора. Поскольку при этом -—=-т—ПР—ПР-— > 0, то с увеличением

дСт (2Ст • От + кПР • РПР)

расхода теплоносителя через прибор От коэффициент передачи к1 растет, что свидетельствует об увеличении влияния температуры теплоносителя на входе в отопительный прибор на его среднюю температуру.

С увеличением произведения коэффициента теплопередачи отопительного прибора на его площадь поверхности кПР • ^ПР коэффициент передачи к1 однозначно уменьшается, а к2 одно-

дк2 2Ст • Gm

значно растет, так как ---=-—>0, т. е. увеличивается влияние

д(кш • Fпp) (2Ст • Gm + кпр • Fпp)2

температуры внутреннего воздуха на среднюю температуру отопительного прибора, а влияние

ослабевает. Влияние кПР • на постоянную времени Т не так очевидно, так как с увели-

]—• ГГ1 СПР

чением гПР знаменатель выражения Т =- растет, однако и числитель

2Ст • + кПР • -ПР

Спр = Ст -рт • УПР + См • тм тоже растет из-за увеличения массы прибора тм и его объема УПР, поэтому все зависит от того, что растет быстрее - числитель или знаменатель. Если же при постоянных значениях , его массы тм и объема УПР каким-то образом увеличивается кПР, то

Т СПР Т растет только знаменатель выражения Т =-, поэтому постоянная времени Т

2Ст • Gm + кПР • -ПР

будет уменьшаться.

Если режим статический (стационарный), то из формулы (4) следует, что

Т =---Вх +-кПР—-ПР--¿В (6)

2Ст • + кПР • -ПР 2Ст • + кПР • -ПР

Данное соотношение можно переписать следующим образом:

2Ст • °т • tТ + кПР • -ПР • tТ = 2Ст ' °т • tВХ + кПР • -ПР • tВ . (7)

Выполнив несложные преобразования, получим широко известное в литературе соотношение

кПР • -ПР</Т - tВ) = Ст • </вХ _ *ВЫХ ) , (8)

что подтверждает справедливость проделанной процедуры.

В целом, все вышеизложенное согласуется с простыми физическими соображениями, что свидетельствует, по меньшей мере, о качественной адекватности математической модели.

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

Известно, что мощность отопительного прибора определяется следующим

соотношением:

Жпр = кПР • -ПР^Т - tв). (9)

Поэтому температура теплоносителя может быть вычислена как

Т = tВ + Жпр/(кпр • -ПР). (10)

Подставляя данное выражение в (4) и выполнив некоторые преобразования, получим, что

_СПР__^ +Ж = 2Ст • Ст • кПР • -ПР • ^ _t )_ СПР • кПР • -ПР • ^В (Ц)

2Ст • °т + кПР • -ПР & 2Ст • °т + кПР • -ПР 2Ст • °т + кПР • -ПР ^

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

С

Обозначим выражение -ПР— как постоянную времени ТВ, т. е. будем считать, что

2Ст • Gm

С

Тв =-ПР—, и перепишем формулу (11), используя ранее принятые обозначения, тогда будем

2Ст • Gm

иметь

T+Wnp = k^p^вх - 2CmGmk2 (tB + TB (12)

ат ат

Если полагать, что расход теплоносителя это постоянная величина, т.е. что Gm = const, то все коэффициенты в данном уравнении будут постоянными величинами и поэтому в изображениях Лапласа при нулевых начальных условиях данное уравнение перепишется в следующем виде:

(Tp + 1) Wnp(p)=kiknpFnp tBx(p)-2CmGmk2(1+ Твp)tB(p) . (13)

Здесь Wnp(p), tBX(p) и tB(p) Лапласовы изображения величин Wnp(x) , tBX(x) и tB(x) соответственно.

Bведем следующие обозначения: к3 = k1knpFnp, k4 = 2CmGmk2, тогда уравнение (13) запишется следующим образом:

(Tp + 1) Wnp(p)=k3tBx(p)-k4(1 + Tbp)tB(p), (14)

или же

Wnp(p) = -k-tBX(p)-^ tB(p) . (15)

(Tp+1) (Tp+1)

Из уравнения (15) следует, что передаточная функция по каналу «температура теплоносителя на входе в отопительный прибор - тепловая мощность отопительного прибора» представляет

к - 1 Wnp(p) к3 собой инерционное звено 1-го порядка-=-, а по каналу «температура внутреннего

tBX (p) (Tp + 1)

воздуха - тепловая мощность отопительного прибора» - инерционно-форсирующее звено с

11 Wnp( p) k4(TB p+1) отрицательным коэффициентом передачи -=--. Т ак как

tB( p) (Tp+1)

1в~

TB =-np— >T=-np-, то данное звено по своим свойствам ближе к

2Ст ' <^т 2Ст ' <^т + кПР ' ^ПР

дифференцирующему и инерционно-дифференцирующему звеньям [2]. Причем с увеличением расхода От обе постоянные времени Т и ТВ уменьшаются.

Из уравнения (11) также следует, что статическая модель отопительного прибора будет иметь вид:

Ш - 2Ст ' ^т • кПР • ^ПР (* * ) (16)

ШПР-^—. , , ,----(*ВХ -tВ), (16)

2Ст • ^т + кПР • ^ПР

что полностью соответствует известным в литературе данным, в частности, данным работы [5].

Зависимость мощности отопительного прибора от расхода теплоносителя От весьма сложная, при изменении От меняются как все коэффициенты передачи к1, к2, к3,к4, так и постоянные времени Т и ТВ обоих каналов. В этом случае по принятой в теории автоматического управления классификации [3] отопительный прибор следует отнести к нестационарным объектам с переменными во времени коэффициентами передачи к1, к2, к3,к4 и переменными постоянными времени Т и ТВ . Если при этом расход От является управляемой величиной, то к1, к2, кз, к 4 , Т и ТВ будут детерминированными функциями времени, если же нет - то стохастическими. Тем не менее, в целом, дифференциальное уравнение (11) является линейным дифференциальным уравнением относительно искомой величины - мощности отопительного прибора ШПР, входные же величины *ВХ (т) и (т) включены в это уравнение линейным образом, а расход От - нелинейным образом. Как нам представляется, в данном случае для анализа и расчета динамики предпочтителен метод дифференциальных уравнений.

Заметим также, что, несмотря на то, что все вышеизложенное относится непосредственно к отдельному отопительному прибору, оно применимо также и для всей системы отопления здания в целом. Объясняется это тем, что (см. [4, 5]) систему отопления здания условно можно

представить эквивалентным отопительным прибором. Это весьма удобно и это существенно упрощает расчет, поскольку современные системы отопления зданий являются весьма сложными техническими устройствами, включающими в реальных условиях, как правило, разнотипные отопительные приборы с различными сроками эксплуатации и, возможно, с неодинаковыми схемами подключения к подводящим и обратным теплопроводам. При этом также следует иметь в виду, что структура и состав этих и так достаточно сложных систем довольно часто меняются, в частности, из-за пожеланий собственников квартир и по другим причинам. В связи с этим непосредственный анализ и отслеживание режимов работы таких систем является весьма трудной задачей. Причем для решения задачи эквивалентной замены достаточно разработать процедуры оценки параметров knp и Fnp всей системы отопления по реальным данным. Детально данные вопросы рассмотрены в работе [5].

Алгоритмическая структурная схема отопительного прибора (системы отопления) и теплового режима здания (ТРЗ) для данного случая, когда расход теплоносителя считается постоянной величиной Gm = const, представлена на рис. 1.

На рис. 1 Z (т) - вектор возмущающих воздействий ТРЗ, основным из которых является tH -температура наружного воздуха, сплошными линиями обозначены сигнальные воздействия, а штриховыми - параметрические, т. е. воздействия, определяющие численные значения параметров (коэффициентов) модели.

Рис. 1. Структурная схема отопительного прибора

Из рис. 1 видно, что математическая модель системы отопления вместе с моделью ТРЗ составляют замкнутую систему уравнений, полностью описывающую поведение объекта управления, который в данном случае включает систему отопления и собственно ТРЗ. Входными величинами для данного объекта являются температура теплоносителя на входе системы отопления ^ВХ(т) , расход теплоносителя через систему отопления Gm (т) и вектор возмущений ТРЗ Z(т). Выходные величины объекта - это мощность (тепловой поток) системы отопления (эквивалентного отопительного прибора) WПР(т) , температура внутреннего воздуха ¿В(т) и температура обратной воды (на выходе системы отопления) £ВЫХ(т).

Разработано два варианта модели ТРЗ: 1-й вариант - это модель с сосредоточенными параметрами^], 2-й вариант - модель с распределенными параметрами [7].

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

ра частичного излучения OptrisMmisight с диапазоном длин волн 8 ^14мкм. Показания снимались в 8 экспериментальных точках: на подающей и обратной подводках и в центрах 1, 3, 5, 7, 9 и 11-й секций радиатора. Экспериментальные точки представляют собой очищенные от краски участки радиатора и труб. Температура внутреннего воздуха помещения измерялась с помощью термоанемометра АТТ-1004 «Актаком». Расход теплоносителя измерялся накладным ультразвуковым расходомером Portaflow РБ330.

Так как измерение мощности (теплового потока прибора) в данном случае представляет некоторую проблему, то задача параметрической идентификации решалась не для уравнения (8), а для модели (5), поскольку проблема измерения температуры в настоящее время достаточно хорошо отработана.

Задача параметрической идентификации математической модели (5) формулировалась в двух вариантах:

I = Z(т*)" tr (^)

i=i

^ min ;

Г, ki, &2

(17)

значения средней n - количество

I=Ё [T () - T ()]2 ^ min, (18)

i=1 T,k1 ,к2

э

где tT (х,), tT (х,) - соответственно экспериментальные и расчетные температуры отопительного приборав моменты съема данных х, экспериментальных точек (моментов съема данных).

Данная задача решалась методом покоординатного спуска со встроенным методом «золотого сечения», при этом дифференциальное уравнение (5) интегрировалось методом Рунге-Кутта с погрешностью, пропорциональной пятой степени шага по времени. В результате решения были найдены численные значения параметров математической модели, так, в частности, по результатам процесса натопа из решения задачи (18) установили, что постоянная времени T=3126 с , коэффициент передачи k1 = 0,624 °С/°С, коэффициент передачи к2 = 0,829 °С/°С . На рис. 2 приведена расчетная кривая tT(х) , вычисленная по уравнению (5) при указанных значениях его параметров, там же точками показаны экспериментальные значения температуры.

Рис. 2. Расчетные и экспериментальные значения средней температуры отопительного прибора в процессе нагрева

Как видно из рис. 2, качество настройки математической модели на экспериментальные данные вполне удовлетворительное.

Практически такие же результаты были получены и из решения задачи (17).

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

Литература

1. Сканави, А.Н. Переходные тепловые процессы в отопительных приборах / А.Н. Сканави, Л.М. Махов, В.Э. Сварич //Изв. вузов. Строительство и архитектура. - 1986. - № 4. - С. 86-88.

2. Теория автоматического регулирования / под ред. А.В. Нетушила. - М.: Высш. школа, 1976. - 400 с.

3. Теория автоматического управления. Ч. 1. Теория линейных систем автоматического управления / под ред. А.А. Воронова. - М. : Высш. школа, 1977. - 304 с.

4. Панферов, В.И. К теории управления режимами централизованного теплоснабжения /

B.И. Панферов, С.В. Панферов // Вестник ЮУрГУ. Серия «Строительство и архитектура». -2011. - Вып. 12. - № 16 (233). - С. 41-45.

5. Панферов, С.В. Адаптивная идентификация параметров систем отопления зданий./

C.В. Панферов, В.И. Панферов // Вестник ЮУрГУ. Серия «Строительство и архитектура». -2014. - Т. 14. - № 2. - С. 33-37.

6. Панферов, В.И. Идентификация и управление тепловым режимом зданий /В.И. Панферов, А.Н. Нагорная, Е.Ю. Пашнина //Вестник УГТУ-УПИ. Серия «Строительство и образование». -2005. - Вып. 14. - № 14 (66). - С. 351-353.

7. Панферов, В.И. Об оптимальном управлении тепловым режимом зданий /В.И. Панферов, Е.Ю. Анисимова, А.Н. Нагорная // Вестник ЮУрГУ. Серия «Энергетика». - 2007. - Вып. 8. -№ 20 (92). - С. 3-9.

8. Панферов, В.И. Об одном решении задачи построения динамической математической модели отопительных приборов и систем / В.И. Панферов, А.Н. Нагорная, Ю.В. Кунгурцева // Вестник ЮУрГУ. Серия «Строительство и архитектура». - 2012. - Вып. 15. - № 38 (297). -С. 46-49.

Панферов Владимир Иванович, д-р техн. наук, профессор, заведующий кафедрой теплога-зоснабжения и вентиляции, Южно-Уральский государственный университет, г. Челябинск; tgsiv@mail.ru.

Панферов Сергей Владимирович, канд. техн. наук, доцент кафедры теплогазоснабжения и вентиляции, Южно-Уральский государственный университет, г. Челябинск; tgsiv@mail.ru.

Поступила в редакцию 10 февраля 2015 г.

DOI: 10.14529/ctcr150208

DYNAMIC MODELOF HEATING DEVICES AND HEAT SUPPLY SYSTEMS

V.I. Panferov, South Ural State University, Chelyabinsk, Russian Federation; tgsiv@mail.ru, S.V. Panferov, South Ural State University, Chelyabinsk, Russian Federation; tgsiv@mail.ru

In order to develop previously performed work we consider the problem of constructing a dynamic mathematical model of heating devices and systems. The heating system at the same time

seems to be equivalent heater. The output value of the model in this case is the power (heat capacity). Decision is executed by the "gray box", that is the structure of the model is determined by analyzing the physics of the processes occurring in heating devices and systems, and the numerical values of its parameters are found from the experimental data by parametric identification. The resulting differential equation is linear with respect to the unknown quantity-power heater(heating), input the same size-water inlet temperature of the heating system and the temperature of the indoor air are included in this equation in a linear fashion, and the flow is in non-linear manner. It is shown that a constant flow of coolant model can be represented by the transfer functions in the iso-mapping Laplace transfer function of the channel "coolant temperature at the inlet to the heater(heating system) -thermal capacity of the heater(heating)" is a inertia component of the first order, and the channel "indoor air temperature-heat capacity of the heater(heating)" - inertia-forcing component with a negative coefficient of transmission. The qualitative features of the model parameters, in particular, their dependence on the coolant flow are analyzed. Algorithmic block diagram of the control object, which includes the heating system and the thermal regime of the building itself is provided. The results of the parametric model identification from experimental data is described. The model can be used in the development of intelligent control systems heating buildings.

Keywords: dynamic mathematical model, the thermal regime of buildings, heating devices and systems, structural synthesis, parametric identification, algorithmic block diagram of the control system.

1. Skanavi A.N., Mahov L.M., Svarich V.Je. [Transient Thermal Processes in Heaters]. Proceedings of Higher Education Construction and Architecture, 1986, no. 4, pp. 86-88. (in Russ.)

2. Netushil A.V. Teoriya avtomaticheskogo regulirovaniya [The Theory of Automatic Control], Moscow, Higher School, 1976, 400 p.

3. Voronov A.A. Teoriya avtomaticheskogo upravleniya [Automatic Control Theory], Moscow, Higher School, 1977, 304 p.

4. Panferov V.I., Panferov S.V. [To the Theory of Control Modes for District Heating]. Bulletin of South Ural State University. Ser. Construction Engineering and Architecture, 2011, vol. 12, no. 16 (233), pp. 41-45. (in Russ.)

5. Panferov V.I., Panferov S.V. [Adaptive Identification of Parameters for Heating Systems in Buildings]. Bulletin of South Ural State University. Ser. Construction Engineering and Architecture, 2014, vol. 14, no. 2, pp. 33-37. (in Russ.)

6. Panferov V.I., Nagornaya A.N., Pashnina E.Ju. [Identification and Thermal Management of Buildings]. Bulletin of Ural State Technical University. Ser. Construction and Education, 2005, vol. 14, no. 14 (66), pp. 351-357. (in Russ.)

7. Panferov V.I., Anisimova E.Ju., Nagornaya A.N. [About the Optimal Thermal Control of Buildings]. Bulletin of South Ural State University. Ser. Power Engineering], 2007, vol. 8, no. 20 (92), pp. 39. (in Russ.)

8. Panferov V.I., Nagornaya A.N., Kungurtseva Ju.V. [Solution to the Problem of Dynamic Mathematical Model of Heating Appliances and Systems Development]. Bulletin of South Ural State University. Ser. Construction Engineering and Architecture, 2012, vol. 15, no. 38 (297), pp. 46-49. (in Russ.)

References

Received 10 February 2015

БИБЛИОГРАФИЧЕСКОЕ ОПИСАНИЕ СТАТЬИ

REFERENCE TO ARTICLE

Панферов, В.И. Динамическая модель отопительных приборов и систем / В.И. Панферов, С.В. Панферов // Вестник ЮУрГУ. Серия «Компьютерные технологии, управление, радиоэлектроника». - 2015. - Т. 15, № 2. - С. 75-82. DOI: 10.14529Мсг150208

Panferov V.I., Panferov S.V. Dynamic Model of Heating Devices and Heat Supply Systems. Bulletin of the South Ural State University. Ser. Computer Technologies, Automatic Control, Radio Electronics, 2015, vol. 15, no. 2, pp. 75-82. (in Russ.) DOI: 10.14529/ctcr150208