![](/user_photo/_userpic.png)
книги из ГПНТБ / Специальные вопросы строительной теплофизики учебное пособие
..pdfГ Л А В А 3
МЕТОДИКА ОПРЕДЕЛЕНИЯ ОТНОСИТЕЛЬНЫХ ОШИБОК ТЕМПЕРАТУР ВОЗДУХА И ПОВЕРХНОСТЕЙ ОГРАЖДАЮЩИХ КОНСТРУКЦИЙ ПОДЗЕМНЫХ СООРУЖЕНИЙ ПРИ НАТОПЕ
Расчет систем отопления подземных сооружений, характерной особенностью которых является нестационарный процесс теплопе редачи через ограждающие конструкции, сводится также к опреде лению мощности этих систем, обеспечивающих заданные темпера туры воздуха и поверхностей ограждающих конструкций сооруже ния к концу выбранного периода теплообмена (натопа). Получе ние заданных температур воздуха и поверхностей ограждающих конструкций необходимо для нормального самочувствия обслу живающего персонала, обеспечения бесперебойной работы техно логического оборудования, электронных приборов и т. п. и исклю чения конденсатообразования на внутренних поверхностях ограж дающих конструкций сооружения в течение периода эксплуатации. Последнее условие является особенно важным, когда требуется повышенная надежность в работе систем автоматики, связи, радио электронной аппаратуры и т. п.
Существующие методы расчета мощности системы отопления подземных сооружений в целях упрощения вычислительных работ не учитывают влияние толщины и разницы в значениях теплофи зических характеристик ограждающих конструкций и грунтового массива на формирование этих температур в процессе теплообме на. Возможность этого упрощения объясняют тем, что теплофизи ческие характеристики железобетонных ограждающих конструкций сооружения по своим величинам больше или соизмеримы с вели чинами теплофизических характеристик грунтов.
Однако ни в методах расчета, ни в специальной литературе нет качественной или количественной оценки принимаемого упроще ния.
В связи с этим возникает вопрос, какова будет величина отно сительной ошибки температур воздуха и внутренних поверхностей ограждающих конструкций подземного сооружения в расчетах, связанных с определением мощности систем отопления, при усло вии, что теплофизические характеристики Ограждающих конструк ций и грунта принимаются равными.
Решению этого вопроса для подземных сооружений прямо угольной формы посвящена настоящая глава.
30
Постановка задачи
В настоящее время существуют точные аналитические решения задач для нестационарной теплопроводности при граничном усло вии второго рода (<7=const) полупространства, как ограниченного, так и неограниченного слоем однородного вещества, толщиной /м, с теплофпзпческими характеристиками, отличными от теплофизи ческих характеристик полупространства. Существующие в практи ке проектирования расчетные формулы для определения мощности систем отопления подземных сооружений в виде параллелепипеда или прямоугольного сечения большой протяженности в своей осно ве содержат эти точные аналитические формулы.
Поэтому исследование изменения величин относительных оши бок температур воздуха и поверхностей ограждающих конструкций сооружения вначале целесообразно произвести путем сравнения результатов общих решений точных аналитических задач для по лупространства, ограниченного слоем однородного вещества (двух слойная задача) и не имеющего этого слоя (однослойная задача).
Решение задачи
1. Для полупространства. Из теории теплопередачи [6] извест-
ккал
но, что при постоянном значении теплового потока q---------- , по-
м г-час
ступающего на поверхность полупространства, повышение темпе
ратуры Д tx слоя, расположенного на |
расстоянии х метров от по |
|||
верхности, будет равно |
х~ |
|
х |
|
Д t X к |
|
|
||
4ат — х erfc |
2 Y a k |
(Ш.1) |
||
где Д tx = t{x,-) — t(o,o) , град; |
|
|
|
|
a = —----- коэффициент температуропроводности, |
— — ; |
|||
у с |
период теплообмена, час\ |
час |
||
т — расчетный |
|
|||
q — удельный |
тепловой |
поток, |
ккал |
|
---------- ; |
|
м2 •час
,, ккал
\— коэффициент теплопроводности,--------------- —;
м■час ■град
ккал
ус —объемная теплоемкость, ---------- .
м ’-град
Изменение температуры на поверхности полупространства с те
чением времени можно получить из выражения |
(III.1), если под |
ставить значение х, равное нулю |
|
Д*(о. ,)= - £ - • 2 | / ^ • |
(III .2) |
31
Изменение температуры как в однородном слое толщиной 1м, имеющем теплофизические характеристики Яь yi и с\ и ограничи
вающем полупространство, так и в самом |
полупространстве, имею |
|||||||||||
щем теплофизические характеристики |
Х2, уд и с2, |
при воздействии |
||||||||||
постоянного удельного |
теплового |
потока |
|
ккал |
|
|||||||
q ---------- . выражается |
||||||||||||
формулами [7]: |
|
|
|
|
|
|
|
|
м2-час |
|
||
|
|
|
|
|
сс |
|
|
2 я / ф * |
||||
t\ (*,7) -- t\ (*,0) |
2 д }/ ах |
|
|
|||||||||
2 |
/4я и |
|||||||||||
|
|
к |
|
|
2 |
|
X |
|||||
|
|
|
|
|
гг—О |
|
|
|||||
“ |
|
/ 2 |
n l — х |
|
О < х < /; |
|
(Ш.З) |
|||||
У Лпо |
|
•------= |
- |
|
|
|||||||
|
|
\ |
2 V а, ' |
|
|
|
|
|
||||
4 <7~|/ а, |
а 2 т |
|
v |
A, |
„ |
i ^ |
|
|
||||
^2 (*,т) — ^2 (х,0) = |
+ |
Х2 ] / flj |
|
2 ] / а ( т |
||||||||
Xj | / а 2 |
"=« |
L2 1/а2 т |
||||||||||
|
|
|
/ < Ф < |
, |
|
|
|
(Ш.4) |
||||
где Xj и а1— коэффициенты |
|
теплопроводности |
и |
температуро |
||||||||
проводности |
|
слоя; |
|
|
|
|
|
|
||||
/ — толщина |
слоя, |
м\ |
|
|
|
|
и |
температуро |
||||
Х2 и а2— коэффициенты |
|
теплопроводности |
||||||||||
проводности |
|
вещества |
полупространства; |
|||||||||
|
|
|
|
|
|
^2"(2^*2 |
|
|
|
|
||
А = |
|
1 |
~ |
1- |
xi Ti <4 . |
|
|
|||||
|
|
|
|
|
V |
Х2Т2Ф |
|
|
|
|||
|
|
|
1 |
+ |
У Ti ci |
|
|
|
|
|||
ф !)= |
— = |
е--2— |
; erfcЕ; |
|
|
|||||||
|
|
|
У * |
|
|
|
|
|
|
|||
_ 2 я/ ф д- |
„ __ 2 nl — x |
|
|
|||||||||
|
2 |
|
|
|
’ |
?2 “ |
~2 У а у У ’ |
|
|
|||
. |
х - 1 |
, |
( 2 я ф 1)ф |
|
|
|||||||
3 |
2 У а 2 т Т 2 У aY- |
|
|
|
||||||||
|
я = |
О, |
1, 2, 3, 4 .... |
|
|
|
|
|||||
Измёнение температуры |
на поверхности |
слоя |
с течением вре |
мени можно получить из выражения (Ш.З), если подставить зна чение х, равное нулю
+ |
. |
4 ql |
г-р,— |
0,28209 |
ос |
я |
|
V A n i |
|||||||
40,0 — Г(0,0) = |
------ V |
г °1 |
(Ш.5) |
УРо,
32
где
Fo1— критерий Фурье; Fol
/2
Относительную ошибку в вычислении температуры поверхности полупространства, ограниченного слоем вещества, при условии ра венства теплофизических характеристик слоя и полупространства можно выразить формулой
|
|
|
|
К = |
( |
1 |
|
п |
100, |
%, |
|
|
(111.6) |
|||
|
|
|
|
М ' |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
где |
|
Д + — изменение температуры |
поверхности полупространст |
|||||||||||||
|
|
|
ва, ограниченного |
слоем |
однородного |
вещества, в |
||||||||||
|
|
|
конце выбранного периода теплообмена т; подсчиты |
|||||||||||||
|
|
|
вается по формуле (III.5); |
|
|
|
|
|
||||||||
|
|
А/п — изменение температуры поверхности полупространст |
||||||||||||||
|
|
|
ва в конце того же периода теплообмена т; подсчи |
|||||||||||||
|
|
|
тывается по формуле (III.2). |
|
|
|
|
|||||||||
Из формулы (III.6) видно, что для определения диапазона из |
||||||||||||||||
менения относительной ошибки |
|
необходимо рассмотреть изме- |
||||||||||||||
|
|
|
|
|
|
|
Д tn |
|
|
|
|
от всех факторов,' |
||||
нение величины отношения ----- в зависимости |
||||||||||||||||
|
|
|
|
|
|
|
Д + |
|
|
|
|
|
|
|
|
|
влияющих на это изменение. |
|
|
|
|
д |
|
|
|
||||||||
Проведем анализ |
величины отношения |
|
|
|
||||||||||||
д t |
’ |
|
|
|||||||||||||
Применяя |
точные |
|
|
|
|
|
|
|
|
|
||||||
аналитические формулы (III.2) и (III.5), |
||||||||||||||||
можно написать следующее выражение для отношения |
Д С |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Д 4 |
|
|
А + |
4 ql VF o, |
0,28209 + |
у, |
А п и |
|
r__ |
А У |
|
|
|||||||
|
|
|
|
|
|
fix |
|
\ VFo, |
|
|
(HI.7.) |
|||||
Д tn |
|
|
|
|
\ q l V |
a |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|||||||
Так |
как в |
выражении |
(III.7) |
\ = |
\ , |
|
а ^ - а 2, т = |
Тг> |
с |
|||||||
V F o l = У |
" |
, то |
оно |
упрощается и принимает вид |
|
|||||||||||
Д + |
2 |
|
V |
и ь |
0,28209 + |
5 |
А ”и |
( у = |
|
. ( 111. 8) |
||||||
|
|
+ |
|
|
|
|
|
|
|
|
||||||
Из |
равенства |
(III.8) |
видно, |
что |
величина |
отношения |
д |
|||||||||
следовательно, |
и of |
зависит |
от |
численных |
значении |
|||||||||||
величин |
3 Зак. 434 |
33 |
1 / |
и 2 |
|
п |
. Основная трудность непосредственно- |
||
A " u h |
||||||
I ^1 Т1С1 |
П |
1 |
\ ’ Fo, |
|
|
|
го вычисления |
|
Д t ' |
заключается в |
|||
отношения — |
по формуле (II 1.8) |
|||||
|
|
|
|
00 |
I |
п |
расчетах, связанных с определением суммы ряда У |
А пч |
- г---- |
||||
|
|
|
|
n = i |
\ |
V Fo1 |
С инженерной точки зрения желательно математические действия, связанные с вычислением этой суммы ряда, заменить более про стыми, которые давали бы требуемую точность для практических расчетов. Эту замену можно произвести только в том случае, ког да принятая простая аналитическая формула будет в допустимых пределах точности эквивалентна выражению суммы ряда.
|
|
|
|
|
|
1 — 'I |
к2Т 2 С2 |
|
||
Необходимо отметить, |
что так как |
А |
V |
\ Ti и |
то |
|||||
1 / |
У |
ТСг2 |
||||||||
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|/ |
I.J "f, С] |
|
||
необходимо проанализировать |
два |
|
/~ \ |
'Y., С% |
1 |
|||||
случая: у |
|
Ti с\ |
||||||||
|
|
^'2 Т2^2 |
|
|
|
|
|
|||
|
При у |
1 |
двухслойная |
задача выро- |
||||||
|
У T l |
|
||||||||
T i |
с \ |
С; |
|
Д tn'_ |
|
|
|
|
||
ждается |
в однослойную. |
Отношение |
принимает |
значение, |
||||||
|
|
|
|
|
Д t„ |
|
|
|
|
|
равное |
единице, a ofn — равное |
нулю |
|
|
|
|
|
|||
В первом случае, когда 1 / |
,'2 12 |
2 |
< 1, |
произведение величин |
||||||
|
|
V |
к \ T i |
с \ |
|
|
|
|
|
теплофизических характеристик слоя больше произведения вели чин теплофизических характеристик вещества полупространства. Коэффициент А в этом случае будет всегда величиной положитель-
°° |
|
|
/ |
а |
\ |
будут так- |
ной и, следовательно, все члены ряда У /4 “ и |
\ у |
---- |
|
|||
|
Д |
К / |
F o J |
|
||
же положительны. Величина отношения |
|
будет меньше еди- |
||||
|
----- |
|
||||
|
Д |
t n |
|
|
|
|
иицы и, таким образом, о, будет иметь отрицательное значение.
Следовательно, в этом случае расчетные температуры поверхности полупространства будут иметь завышенные значения.
/ I */ |
произведение вели- |
. 2 У 2 О 1, |
г'М T сl \
чип теплофизических характеристик слоя будет меньше произведе ния величии теплофизических характеристик вещества полупрост ранства, т. е. слой будет менее теплопроводен по отношению к по лупространству. Коэффициент А в этом случае будет всегда вели-
34
|
00 |
I n |
бу |
|
чиной отрицательной и, следовательно, ряд У |
Ап о I —j=^- |
|||
дет знакопеременным. |
1 |
\ У |
Fol j |
|
м - |
|
|
|
|
Величина отношения |
|
а 8,п |
будет |
|
— будет больше единицы, |
A t,
иметь положительное значение. Следовательно, в этом случае
расчетные |
|
температуры |
поверхности |
полупространства будут |
||||||
иметь заниженные значения. |
формулы, |
заменяющей |
выражение |
|||||||
Для вывода |
расчетной |
|||||||||
0,28209 + |
00 |
/ |
/£ |
\ |
для различных значений |
коэффици- |
||||
у |
Ап и ( |
: |
) , |
|||||||
|
|
о__ |
\ V'Fo, |
|
|
|
|
|
||
еыта А и У Fol |
|
были произведены вычисления сумм рядов |
при |
|||||||
А > 0 и Л -У 0 с точностью до седьмого знака. По результатам |
вы- |
|||||||||
численпй построен график (рис. 11) зависимости s = / |
*(*) |
С2 |
||||||||
У Ti |
с\ |
|||||||||
|
|
|
|
|
|
|
|
Из графика видно, что величина s является линейной функцией в
логарифмических координатах и может быть |
выражена |
зависи |
||||||||||||
мостью |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,28209 |
X, |
, с |
|
|
|
|
|
|
(I1I.9) |
||
где к — угловой коэффициент, зависящий от величины У Fo^ |
||||||||||||||
Зависимость k в пределах от У Fo, = 5 |
до |
У Fol = 100 |
полу |
|||||||||||
чаем из графика (рис. 12) |
зависимости |
k = |
ср ( У Fo^. |
|
|
|||||||||
Эта зависимость |
в аналитической форме имеет вид |
|
|
|||||||||||
|
|
|
k |
|
ypQx |
- |
1 |
|
|
|
|
|
(ШЛО) |
|
|
|
|
|
У Fo, + |
1 |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
||||
Таким |
образом, |
|
сумма s ~=0,28209 |
оо |
Ап ч .I- |
п |
в пре- |
|||||||
|
у |
___ |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
\ |
У Fo, |
|
|
делах от У Fo, = |
5 до У Fo, - |
100 |
может |
быть заменена |
анали- |
|||||||||
тической |
формулой |
вида |
|
|
|
|
|
|
|
|
|
|
|
|
|
s |
= |
0,28209 |
У Тг с 2 |
|
|
|
|
|
|
( 111. 11) |
|||
|
У П И |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
Подставляя значение s |
из |
формулы (111.11) |
в формулу |
(III.8), |
||||||||||
получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A t ' |
|
|
, |
, -L |
|
|
,, |
|
, |
VFF |
|
|
||
= 2 у ~ г |
/ b i i y i |
|
0,28209 |
У Тг сг |
|
|
+ |
|
(III. 12) |
|||||
A |
|
|
У Ti Д |
|
|
|
|
У Ti ci |
|
|
|
|
|
3 * |
35 |
Х2~{2
Рис. 11, График зависимости s от
Ar(i f,
Рис. 12. График зависимости k от Y F ° \
СО-4
Так как произведение величин 2- У я 0,28209, равно единице, а
/ 4 |
. JL , |
ч уТТ~- 1 |
, |
. |
1 |
|
|||
Т2 ^2 \ 2 |
/ 4 Тг С2 \ |
YfbT + 1— / |
Тг ^2 |
jV'/’oi + 1 |
|
||||
\ 4 |
Ti ci ) |
\^i Ti ci I |
|
|
\ 4 |
Ti с\ / |
|
||
то выражение (111.12) упрощается и принимает вид |
|
||||||||
Л.4 |
4 Та ci |
V-Fo, |
+ 1 |
или |
д 4 |
4 |
Ti 4 |
\ Yfo, + 1 |
(111.13) |
Д t„ |
Ti ci |
|
Д 4' |
4 |
Тг с2 |
|
|||
|
|
|
|
||||||
Подставив значение отношения----^ |
из выражения (III.13) в |
||||||||
|
|
|
|
|
Д 4 |
|
|
|
|
формулу (III.6), получим расчетную формулу для определения от носительной ошибки температуры поверхности слоя полупростран ства при равных теплофизических характеристиках слоя и полупро странства. В окончательном виде формула принимает следующее значение:
84 |
А I Ti 4 W/T7+ 1 •100%. |
(III.14) |
|
|
\ 4 Тг С-2 I |
|
|
Для определения относительной ошибки температуры воздуха |
|||
воспользуемся формулой Ньютона |
|
|
|
|
Я = «(4 — 4). |
(III. 15) |
|
, , |
теплоотдачи, |
ккал |
|
где а — коэффициент |
----------------- ; |
|
|
4 — температура |
воздуха, град; |
мг-час-град |
|
|
|
4 — температура поверхности, град.
При нестационарном процессе теплопередачи температура воз духа и поверхности для каждого конкретного случая теплообмена
является функциями времени. |
Исходя из этого |
положения напи |
||||||
шем формулу (III.15) в виде: |
|
|
|
|
|
|||
|
для полупространства |
(однослойная задача) |
|
|
||||
|
(J = <х[(4 |
(т) |
|
t B(0)) |
(4(т> - |
4 (0))]; |
(III. 16) |
|
|
для полупространства, ограниченного слоем |
(двухмерная за |
||||||
дача) |
|
|
|
|
|
|
|
|
|
Ч= « [(*; (т) - |
4 |
(0)) - |
(4 (х) - |
4(о>)], |
(III. 17) |
||
где |
^в(-) и 4 со— температуры |
воздуха в |
конце |
выбранного |
||||
|
периода теплообмена, град; |
|
||||||
|
4<о> и 4<о)— температуры воздуха и поверхности полу |
|||||||
|
пространства, как ограниченного, так и не |
|||||||
|
ограниченного слоем вещества, перед нача |
|||||||
|
лом |
периода теплообмена, град; 4(о>=4(о) |
38
С (т) и *п (-о— температуры поверхностей слоя и полупро
|
|
странства в конце выбранного периода теп |
||||
|
|
лообмена, |
град. |
|
|
|
Приняв следующие обозначения: |
|
|
||||
|
|
|
|
* |
|
|
Д |
^в = |
(-) |
—t B (0), |
M a = |
t n ( - ) — |
7 П(0) , |
Д |
= |
(т ) |
( 0 ) , |
Д tn ~ |
( т ) |
^ П ( О ) , |
напишем формулы (III.16) |
и (III.17) в виде: |
|
||||
|
|
Я =- « |
(Д - |
Д t n), |
|
(111.18) |
|
|
q = а (Д t; — Д tn'). |
|
(III. 19) |
Разрешим уравнения (III.18) и (III.19) относительно Д и Д^':
Д *в - |
— + Д , |
( Ш . 2 0 ) |
|
а |
|
д *„' = |
— + Д С . |
(III.21) |
|
а |
|
Относительная ошибка температуры воздуха, в случае условно го принятия равенства теплофизических характеристик слоя и по лупространства, выражается формулой
Jt |
1 |
|
100%. |
(III.22) |
|
1В |
д |
|
|
|
|
|
|
|
Д^В |
||
Из равенств (III.20) и |
(III.21) |
выразим отношение |
|||
|
|||||
Д tB |
q + аД tn |
(111.23) |
|||
Д t ' |
|
аД t ' |
|||
|
|
||||
Заменяя в равенстве (III.23) |
Д |
его значением из |
формулы |
||
(III.13) и подставляя полученный |
результат в формулу |
(III.22), |
получим расчетную формулу для определения относительной ошиб ки температуры воздуха
аД t„ |
•100%. |
(111.24) |
1 |
q+ «Д tB 2 i2^2 \VFox + 1
2. Для подземных сооружений прямоугольного сечения. Для вывода расчетных формул относительных ошибок температур воз духа и поверхностей ограждающих конструкций подземных соору жений, в случае условного принятия равенства значений теплофи зических характеристик ограждающих конструкций и грунтового массива, воспользуемся следующими формулами для определения мощности систем отопления [8].
39