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

УДК 697.03.01

К ТЕОРИИ МОДЕЛИРОВАНИЯ СИСТЕМ ОТОПЛЕНИЯ

Ю.Н. Денисенко, В,И. Панферов

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

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

По-видимому, исторически первой математической моделью отопительного прибора является широко известное уравнение для определения его мощности:

а =*(0,5(4*+0-0^ (!)

где К и .Г - коэффициент теплопередачи и площадь поверхности отопительного прибора; и {вЫх ~ температура теплоносителя на входе в прибор и выходе из него; ^ - температура воздуха в помещении.

Уравнение содержит как входные переменные tgy.Ts.tg, так и выходную переменную /вых .

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

А і =

ta.

(2)

Формула (2) будет давать точное значение среднего температурного напора только в том случае, если температура теплоносителя вдоль Р будет распределяться по прямой и будет при этом знакопостоянной функцией. В самом деле, хорошо известно, что средняя температура теплоносителя определяется по формуле:

= 0, 5(/ет +Івьа) ;

(3)

где Б1т - площадь прямолинейной трапеции, образованной линейной знакопостоянной зависимостью /(У7) и осью абсцисс на отрезке [0; /<], что

подтверждает справедливость вышеприведенного замечания по поводу вычисления At по формуле (2).

При расчете любого теплообменного аппарата неизвестные составляющие могут быть определены из уравнений теплопередачи или теплового баланса. Не является исключением и отопительный прибор, для которого

Оо=С°(1«х-1вых)’ (4)

где с и в - теплоёмкость и массовый расход теплоносителя.

Таким образом, решая уравнения (1) и (4), найдем первый вариант структуры статической модели отопительного прибора, представляющийся следующими уравнениями:

С°(‘вх ~0.

во!

0,5 +

сЄ

КР

а

ев’

(5)

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

Рис. 1. Структурная схема статической модели отопительного прибора для системы уравнений (5)

На рис. 1 и б - входные переменные;

при этом 1вх и в являются управляющими воз-

действиями, соответственно /в следует отнести к категории возмущений.

Для отыскания другого варианта структуры математической модели выделим бесконечно малый участок площади поверхности отопительного прибора и запишем уравнение теплового баланса для этого участка [2]:

-ОсЛ = К{г-гв)<Ш, (6)

где t - температура теплоносителя в пределах

участка <1Р.

Если считать коэффициент теплопередачи постоянным, то уравнение (6) будет представлять собой уравнение с разделяющимися переменными, которое можно представить в виде:

= ¿Р. (7)

t-tв вс

Решая (7) с учетом того, что ¿(0) = 1вх, получим

( КР

*00 = 1в+((вх-1в)-ехр\-—

Уравнение (8) описывает изменение температуры теплоносителя вдоль поверхности теплообмена при постоянном коэффициенте теплопередачи и если Т7 вся поверхность теплообмена, то

*(р) = гвых.

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

во =

1-е св

сО(!вх~*в)-

О.

Рис. 2. Структурная схема статической модели отопительного прибора для системы уравнений (9)

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

Qo2

_к^ 1-е с°

¿вых ¿в &вх ¿в ) ‘ ®ХР

С^(!вх ^в)>

КР

(9)

вс

Структурная схема модели отопительного прибора в этом случае представлена на рис. 2.

С другой, стороны из теории рекуперативных теплообменных аппаратов (ТО), работающих, например, по схеме «противоток», известно, что мощность теплообменника следует вычислять по формуле:

Ото = £' с<~*\(Цех -¡2вх)’

где ¡1вх и 12вх - температуры первичного и вторичного теплоносителя на входе в теплообменник; сС?1 - эквивалент расхода первичного теплоносителя; £ - коэффициента эффективности ТО, опре-

(10)

(8) КР С°1 1

) 1-ес сС?2

£ =-

1-

сСх

сС?-)

КР

сО\

сС2

где сС2 - эквивалент расхода вторичного теплоносителя.

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

£0=1-есв. (11)

Параметр е0 позволяет определять мощность ОУ через большую разность температур теплооб-менивающихся сред

2о =£о-сС(‘ех-Св)> (12)

откуда совместное решение с уравнением (11), позволяет получить зависимость, полностью совпадающую с уравнением статики (9):

во2

1 е

КУ ^

' св

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

а

I- ^еыл*

а £

Рис. 3. Схема узла смешения перед системой отопления

На рис. 3 смесительным устройством является водоструйный элеватор (Э), перед которым «средствами автоматики» [3] установлен расход сетевой воды Ос с температурой /вхс > . Физический

смыл перехода с температурного режима наружных сетей на температурный график системы отопления определяется через уравнение неразрывности струи посредством вычисления коэффициента смешения

и = -

(13)

Выразим температуру воды на входе в ОУ 1в_ через коэффициент смешения с учетом (5)

.во

t =/ —и -

1вх 1вх.с и ^ ’

си

подставив которую в (5) для £)(Л получим следующую модификацию уравнения статики, позволяющую находить мощность отопительного прибора через разность температур теплообменивающихся сред:

Оо=ео-сС(.*вх.с-*в)> О4)

откуда коэффициент эффективности равен

-Чгг <15>

0,5 + и л--

№

При совпадении установленных температурных графиков в системе отопления и тепловых сетях необходимость устройства узла смешения отпадает, то есть и = 0, тогда коэффициент эффективности примет вид:

1 (16)

0,5 +

св ’ №

откуда ясно, что при подстановке (16) в (12) также прейдем к виду структуры статической модели (5).

Из рис. 3 и уравнения (13) следует, что в = (1 + и)<5с . С учетом этого мощность ОУ можно представить следующим выражением:

Оо=Ео-сСЛ*вх.с-1в)> О7)

где коэффициент эффективности соответственно равен

1

е° 0,5

и сОс 1+и~+ЯГ

(19)

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

знать только установленные входные параметры и значение коэффициента смешения.

Для получения третьего варианта структуры статической модели отопительного прибора учтем зависимость коэффициента теплопередачи от среднего значения температурного напора и расхода. Как известно [4] эта зависимость имеет вид: К = тАіпЄр, (20)

где т, п и р - экспериментальные числовые показатели, зависящие от типа прибора и характера циркуляции в нем воды; А/ - текущий температурный напор; С/ - относительный расход воды.

Подставив соотношение (20) в уравнение (6), а затем, проинтегрировав, найдем выражение для вычисления текущего температурного напора отопительного прибора:

тп

~вс

.вРр + ({вх-1вГп

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

ґ^.врр+(івх-іву"Л

\ СтС

(21)

позволяющим определить ?вьи при конкретном значении

В работе [2] построены кривые распределения температуры теплоносителя вдоль площади поверхности отопительного прибора для различных значений расхода, из которых видно, что формулы (8) и (21) дают практически одинаковое распределение температуры теплоносителя вдоль поверхности теплообмена. Различие становится существенным только при больших значениях Р и малых значениях расхода теплоносителя.

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

С помощью уравнения (21) можно получить третье уравнение статики для определения О,:

1'

бо =сС

-и-

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

_1

тп —г

ОоЗ ~

(22)

Структурная схема модели отопительного прибора в этом случае представлена на рис. 4.

Рис. 4. Структурная схема статической модели отопительного прибора для системы уравнений (22)

Вместе с тем интересно также оценить и количественное различие моделей. Для примера рассмотрим структуры модели из уравнений (5) и (9). Откуда видно, что отношение мощностей отопительного прибора Q02ІQo\ определяется следующим выражением:

0,2 / ^

бої

1-

св

0,5 +

сО

КР

(23)

Как видно из (23) отношение мощностей является функцией параметра КР/(сО) . Представим последний за X и исследуем поведение функции 602/601 в зависимости от X, для чего проведем прямые вычисления (см. таблицу).

Результаты вычислений 0о2/0О1 по формуле (23)

602/601 0,999992 0,694934 0,599873

ЯР/св 0,01 5,01 10,01

ОогЮох 0,566622 0,549975

ЯР/св 15,01 20,01

0о2І0о1 0,50010 0.50005 0,50003

КР/св 10000,0 20000,0 30000,0

бог/бої 0,50003 0,50002

КР/св 40000,0 50000,0

Как видно из таблицы с увеличением параметра КР/сС отношение мощностей уменьшается.

Для построения графика функции 602/601 =Л^ ) найдем предельные значения этой функции при X стремящемся к бесконечности и к нулю справа:

[(1-е-^){0,5 + 11} = 0,5. (24)

Ііт ^2.

0о\

- Ііт

Х-МОО

X

Для вычисления второго предела функции применим правило Лопиталя

Ііт

А"->+0

1-е-* •= 1ІЛ1 • в"*

X 1

0,5Х + 1 (0,5^ + 1)2

1

=1. (25)

Таким образом, график функции (бя/бл)= f{x) будет иметь вид, приведенный на рис. 5.

0 30 60 1-Ю4 2-Ю4 3-Ю4

Рис. 5. Графики функции Оо2/Оо1

Из рис. 5 видно, при малых значениях параметра КР/(сО), то есть при больших расходах теплоносителя, варианты модели (5) и (9) практически равноценны. При больших значениях этого параметра, имеющих место, очевидно, при малых расходах, мощность отопительного прибора, вычисленная по второму варианту структуры будет в 2 раза меньше, чем мощность найденная по первому варианту структуры. Поэтому при прочих равных условиях температурные графики регулирования, построенные на основе первого варианта, уравнение (5), будут располагаться ниже графиков, построенных на основе второго варианта -уравнение (9). А поскольку физика процесса теплопередачи точнее учитывается вторым вариантом модели чем первым, то теплоснабжение зданий, рассчитанное на основе первого варианта модели (в настоящее время используется именно этот подход) будет приводить к тому, что количество доставленной потребителю тепла не восполнит его тепловую потребность, обусловленную потерями тепла через ограждающие конструкции.

Анализ данных работы [2] говорит о том, что графики распределения температуры теплоносителя вдоль поверхности теплообмена Р, вычисленные по второму, и третьему вариантам структуры математической модели, практически одинаковы и, очевидно, что при использовании этих структур

погрешность в определении теплоотдачи отопительного прибора будет менее значимой.

Вместе с тем это не означает, что одинаковыми будут и вычисленные значения мощностей 0о2 и ()о3, - формулы (9) и (22), поэтому необходимо выяснить какое количественное различие дают эти формулы.

На рис. 6 приведены графики для относительных значений мощностей отопительного прибора, определенные по всем трем вариантам структуры математической модели для различных входных температур теплоносителя (80, 85 и 95 °С). Указанные на рис. 6 кривые позволяют провести анализ по сходимости результатов расчета между структурами уравнений (5), (9) и (22).

За базовое (расчетное) значение мощности 0% отопительного прибора использовалась величина 1536 Вт, найденная стандартным образом [4] для 10-секционного отопительного прибора М 140 АО. Здесь индекс I означает номер варианта структу-

ры. За базовое (расчетное) значение расхода Ор

была взята величина 0,015 кг/с, температурный график 95-70 °С.

Расчет проводился методом последовательных приближений путем увеличения расхода теплоносителя и нахождения ()01 - до достижения

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

Так, сплошные кривые построены по первому уравнению статики (5), пунктирные по второму уравнению статики (9), а сплошные кривые, обозначенные «ромбами», по третьему уравнению статики (22).

Как видно из рис. 6 сплошные и пунктирные кривые практически схожи при данных значениях температуры теплоносителя на входе в отопительный прибор. Такая близость справедлива, так как в

О ¡о*

«¡я» З&э

йЛ 1.0 2.0 ЗД 4,0 5.0 %•

Рис. 6. Кривые изменения мощности отопительного прибора от расхода для различных значений входной температуры теплоносителя

структуру уравнений (5) и (9) входит коэффициент теплопередачи, изменение которого не учитывается.

В то же время, кривые, построенные по уравнению (22) отличаются от остальных и учитывают изменения коэффициента теплопередачи во всем диапазоне регулирования. Наиболее явно различие просматривается при значительном снижении , например до 85 °С, когда для достижения расчетного значения мощности расход необходимо увеличить в 6 раз.

Выводы

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

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

Также видно, что при высоких входных температурах и практически реальных диапазонах измене-

ния расхода мощности <2о1 и ()оЪ практически одинаковы. Следовательно, при управлении процессом теплоснабжения необходимо учитывать и эту особенность, в частности, при качественно-количественном регулировании, и установить целесообразность использования первого варианта структуры статической модели отопительного прибора.

Литература

1. Монахов, ВТ. Моделирование управления режимами тепловых сетей / В.Г. Монахов, Ю.А. Вой-тинская. - М.: Энергоиздат, 1995. - 224 с.

2. Денисенко, Ю.Н. К проблеме построения математических моделей отопительного прибора/ Ю.Н. Денисенко, В.И. Панферов. - Материалы Второй Международной науч.-техн. конф. «Теоретические основы теплогазоснабжения и вентиляции». -М.: МГСУ, 2007. - С. 101-106.

3. Соколов, Е.Я. Теплофикация и тепловые сети: учеб. для вузов /Е.Я. Соколов. - 7-е изд., пе-рераб. - М.: Изд-во МЭИ, 2006. — 472 с.

4. Богословский, В.Н Внутренние санитарно-технические устройства Ч. 1: Отопление/В.Н. Богословский. -М.: Стройиздат. -1990. - 344 с.