Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Специальные вопросы строительной теплофизики учебное пособие

..pdf
Скачиваний:
4
Добавлен:
29.10.2023
Размер:
4.32 Mб
Скачать

Г Л А В А 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

Соседние файлы в папке книги из ГПНТБ