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

книги / Энергоэффективные конструкции в подземном строительстве

..pdf
Скачиваний:
1
Добавлен:
20.11.2023
Размер:
11.44 Mб
Скачать

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

Говоря об энергоэффективных конструкциях зданий и сооружений, прежде всего следует отметить, что данные конструкции, несомненно, еще молоды: их проектирование и эксплуатация осуществляется не более 30–40 лет.

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

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

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

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

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

31

2.1. Основные положения теплопереноса в грунтах

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

проводимость;

конвективную передачу (конвекцию);

процессы испарения и конденсации (скрытая теплопередача);

теплоизлучение;

ионный обмен;

процессы замерзания – оттаивания.

Теплопередача в незамерзшем грунте происходит, в основном, посредством проводимости (qcond) и, в малой степени, за счет конвекции

(ql,conv, qv,conv).

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

ql,conv = cwρw

υ

w (t t '),

(2.1)

где cw – удельная теплоемкость грунтовой воды; ρw – плотность грунто-

вой воды; υw – вектор скорости воды; t′ – исходная температура. Подобное уравнение получают для конвекции с использованием

пара (поровой газ):

qv,conv = cvρv

υ

v (t t '),

(2.2)

где cv – удельная теплоемкость пара в грунте; ρv – плотность пара

в грунте; υv – вектор скорости пара; t′ – исходная температура. Конвективная теплопередача также происходит, если имеется фа-

зовое изменение воды (скрытая теплота во время испарения и конденсации qlat).

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

32

qlat = L0ρw

υ

v ,

(2.3)

где L0 – скрытое тепло испарения при t′.

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

Процессы замерзания и оттаивания грунтов также могут переда-

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

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

Плотность теплового потока – это тепловой поток, отнесенный к единице площади поверхности. Согласно закону Фурье, плотность теплового потока пропорциональна градиенту температуры, т.е. плотность теплового потока qcond можно записать как

q =

Q

= −λ

t

,

(2.4)

 

 

cond

A

n

 

 

 

 

 

 

где Q – тепловой поток через поверхность площадью А; λ – теплопро-

водность; t – температурный градиент в действительном направле-

n

нии потока n:

t

=

t

e

 

+

t

e

 

+

t

e

 

= gradt.

(2.5)

 

 

x

 

y

 

z

n

x

 

y

 

z

 

 

Таким образом, полную (суммарную) теплопередачу в грунте (qtot) можно определить как

qtot = qcond + qv,conv + qv,conv + qlat ,

(2.6)

где qcond – плотность теплового потока за счет проводимости; ql,conv – плотность теплового потока, образованного конвекцией жидкости;

qv,conv – плотность теплового потока, образованного конвекцией пара; qlat – плотность теплового потока при скрытой теплопередаче.

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

33

Уравнение (2.6) можно записать в прямоугольных координатах как

 

t

 

 

t

 

 

t

 

 

 

q = −λ

 

ex

+

 

ey

+

 

ez

= −λgradt.

(2.7)

x

y

z

 

 

 

 

 

 

 

 

Если теплопроводность и температурный градиент постоянны по всей площади и в нормальном направлении, то уравнение (2.4) для энергоэффективной сваи с радиусом R и длиной l можно записать как

Q = 2πRlλ

dt

,

(2.8)

 

 

dr

 

где Q – тепловой поток через поверхность энгергосваи, Вт.

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

 

 

 

 

 

 

 

−ρc

t

= q +

q + q .

 

 

(2.9)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂τ

x

 

y

z

 

 

 

 

Дифференцирование уравнения (2.4) по отношению к пространст-

венным координатам и соединение его с уравнением (2.9) дает

 

 

t

 

 

2t

 

 

2t

2t

 

 

 

 

 

 

2

 

 

 

 

= а

 

 

 

 

+

 

 

 

 

+

 

 

 

= аdiv(gradt) = а

t,

(2.10)

 

∂τ

 

x

2

 

y

2

z

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где а – коэффициент температуропроводности,

м2/с; 2t

оператор

Лапласа:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а =

 

λ

 

; 2t =

2t

+

 

2t

+

2t

,

 

(2.11)

 

 

 

 

ρс

 

x2

 

 

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

 

 

 

где λ – теплопроводность, Вт/(м·К); c – удельная теплоемкость, Дж/(кг·ºС); ρ – плотность твердой среды, кг/м³.

Если в рассматриваемом объеме грунта есть внутренний источник тепла (внутреннее теплообразование), то основное уравнение теплопроводности принимает следующий вид

t

= а 2t +

qυ

,

(2.12)

∂τ

 

 

ρc

 

где qυ – мощность внутренних источников теплоты.

34

Уравнения (2.10) и (2.12) называют основными уравнениями теплопроводности.

Уравнение (2.10) в прямоугольных координатах можно преобразовать в цилиндрические с радиусом r, азимутом φ и осью z:

t

 

2t

 

1 t

 

1 2t

 

2t

 

 

 

= а

 

 

+

 

 

 

+

 

 

 

 

 

+

 

 

.

(2.13)

∂τ

r

2

 

 

 

 

2

 

∂ϕ

2

z

2

 

 

 

r r r

 

 

 

 

 

 

 

В сферических координатах с радиусом r, углом склонения φ и углом восхождения ψ уравнение (2.10) принимает следующий вид:

t

 

2t

 

2 t

 

1

 

2t

 

 

cos ψ ∂t

 

 

 

1

 

2t

 

 

= а

 

 

+

 

 

 

+

 

 

 

 

 

+

 

 

 

 

 

+

 

 

 

 

 

 

. (2.14)

∂τ

r

2

r r

r

2

 

∂ψ

2

r

2

sin ψ ∂ψ

r

2

sin ψ ∂ϕ

2

 

 

 

 

 

 

 

 

 

 

 

 

 

2.2. Аналитические и численные расчеты

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

Начальные условия определяют распределение температурного поля в начальный период времени и заключаются в том, что в начальный момент времени τ0 должна быть известна функция tc = f (x, y, z, τ0 ) .

Граничные условия определяют особенности протекания процесса на границах тела. Рассмотрим эти условия.

Граничное условие первого рода (граничное условие Дирихле). Тем-

пература поверхности tс постоянна или зависит только от времени τ . Иначе говоря, известна функция распределения температуры по поверхности тела:

tc = f (x, y, z, τ).

(2.15)

Граничное условие второго рода (граничное условие Неймана).

Плотность теплового потока qс на поверхности постоянна или зависит только от времени τ , т.е. задана функция распределения плотности теплового потока по поверхности тела:

qc = f (x, y, z, τ).

(2.16)

Граничное условие третьего рода (граничное условие Коши или смешанное условие Неймана). На поверхности рассматриваемого тела происходит теплообмен с жидким или газообразным окружением этого

35

тела, имеющего температуру tж. По закону охлаждения Ньютона тепловой поток пропорционален разности температур между окружающей средой и температурой поверхности tс. Коэффициент пропорциональности определяется как коэффициент теплоотдачи α (Вт/м²·K). С другой стороны, плотность теплового потока может быть выражена в соответствии с законом Фурье:

 

t

= α (tж tс ).

 

−λ

 

 

(2.17)

 

 

∂τ с

 

 

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

λ1

t

 

1

= λ2

t

 

2 .

(2.18)

 

 

n

 

 

n

 

 

 

 

 

При совершенном тепловом контакте изотермы непрерывно переходят из одного тела в другое, а градиенты температур в точках контакта удовлетворяют условию (2.18).

Из-за математических трудностей аналитические решения этих уравнений возможны только для простых случаев. Можно легко решать одномерные задачи, поскольку уравнение (2.4) зависит только от одной координаты. Двухмерные и трехмерные задачи можно решать в некоторых случаях при помощи комбинирования одномерных решений.

Энергоэффективные конструкции фундаментов имеют различные геометрические формы. Более простые элементы, например, сваи можно смоделировать в виде цилиндров, тогда как сложные для конструирования фундаменты или туннели требуют сложных геометрических расчетов. Зарубежными авторами приводятся три основных случая, которые можно использовать для моделирования большинства проектируемых энергоэффективных фундаментов (ЭЭФ) (табл. 2.1).

Таблица 2 . 1 Основные случаи теплопередачи в грунте

№ п/п

Случай

 

 

Дифференциальное уравнение

1

Полубесконечное тело

 

 

 

 

2

Θ(x,t) =

1

 

Θ(x,t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

α ∂τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

Бесконечное тело с цилиндриче-

2

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

 

ским отверстием

 

 

Θ(r,t) +

 

 

Θ(r,t)

=

 

Θ(r,t)

 

 

r

2

r r

α ∂τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

Бесконечное тело со сферическим

 

2

Θ(r,t) +

 

2

 

 

Θ(r,t)

=

1

 

Θ(r,t)

 

отверстием

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

2

 

r

 

r

 

α ∂τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

36

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Случай 1: полубесконечное тело. Моделирует плоские границы раздела между грунтом и атмосферой или фундаментами (плитами, стенами подвального помещения, подпорными стенами, свайными стенами, «стенами в грунте»).

Случай 2: бесконечное тело с цилиндрическим зазором. Моделиру-

ет длинные вертикальные теплообменники в грунте (сваи).

Случай 3: бесконечное тело со сферическим зазором. Моделиру-

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

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

erfc(ζ) = 1

2

0ζ e− ω2 dω,

(2.19)

π

 

 

 

где ζ – предел интегрирования; ω – переменная интегрирования.

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

Случай 1 был решен для различных граничных условий.

Первое решение случая 1 предполагает гармоничное колебание температуры на поверхности (граничное условие Коши или смешанное условие Неймана) с теплопередачей между грунтом и воздухом. Рассматриваются дневные или годовые колебания температуры, причем учитывается отставание колебания температуры в грунте во времени. Предполагается, что поверхность полубесконечного тела обменивается теплом с воздухом, который производит синусоидальное температурное колебание f (τ) . При этих условиях температура грунта колеблется согласно средней годовой температуре воздуха tн, если не принимать во внимание эффекты излучения и геотермальный температурный градиент. Также учитывается уменьшение с глубиной (d) амплитуды колебания из-за тепловой инерции грунта. С учетом данных граничных условий дифференциальное уравнение (2.10) примет вид

37

 

 

 

 

 

z /d

 

 

 

 

z

 

t (z, τ) = tн.г + ∆tн.днη e

 

 

 

 

 

cos

ω (τ − ε )

 

,

(2.20)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d =

 

2a

=

 

aP

;

 

 

(2.21)

 

ω

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

η =

 

 

 

 

 

 

1

 

 

 

 

 

 

 

;

 

 

(2.22)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1+ 2k + 2k 2

 

 

 

ε = arctan

 

k

;

 

 

 

 

(2.23)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1+ k

 

 

 

 

 

 

 

k =

λ

 

=

 

λ

 

 

π

,

 

 

(2.24)

 

 

 

α

 

 

 

 

 

 

αd

 

 

 

 

 

aP

 

 

 

tн.г – средняя годовая температура воздуха, ˚C; tн.дн

– средняя дневная

температура воздуха, ˚C; tн – амплитуда температуры, ˚C;

ω = 2π / P,

1/с; P – продолжительность периода температурного колебания, с; d – глубина затухания; α – коэффициент теплопередачи из грунта в воздух, Вт/(м2·ºС).

На поверхности (z = 0) решение сводится к следующему виду:

t(0, τ) = tн.г + ∆tнη cos(ωτ − ωε).

(2.25)

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

Глубина затухания d зависит от продолжительности периода и температуропроводности а. В грунтах с высокими значениями а (например, влажный песок) температурная волна проникает глубже. Сравнение отклика почвы на годовые и дневные периоды дает следующее отношение:

dгод

=

Pгод

=

365

19 :1.

(2.26)

 

 

1

dдн

Pдн

 

 

Уравнение (2.26) показывает, что годовая температурная волна проникает в 19 раз глубже в грунт, чем ежедневная волна. Более того, тепловой поток в земле достигает своего максимума раньше, чем температура.

38

Второе решение основного случая 1 достигается путем допущения внезапного падения температуры на поверхности полубесконечного полупространства (граничное условие Дирихле). Это решение позволяет изучить теплопередачу между плоским энергоэффективными эле-

ментом и грунтом.

 

 

 

 

При допущении падения температуры с t0 до tс

во время τ = 0 по-

лучают следующее решение уравнения (2.10):

 

 

 

 

 

z

 

 

 

 

 

 

 

t(z, τ) = tсerfc

 

 

 

,

(2.27)

 

aτ

 

 

2

 

 

 

где t (z, τ) – функция температуры с исходной (начальной) температурой t0 ,

 

t

(z, τ) = t(z, τ) t0 ,

(2.28)

откуда

 

 

 

 

 

(z,0) = 0.

(2.29)

 

 

t

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

Из общего ряда можно выделить следующие программные ком-

плексы: GeoStudio (Канада); COSMOS/M (Россия); ANSIS (США) и др.

Термодинамические расчеты, выполняемые представленными комплексами, сводятся к решению основных дифференциальных уравнений термодинамики (2.10), (2.12)–(2.14) методом конечных элементов (МКЭ).

2.3. Теплофизические характеристики грунтов

При выполнении расчетов ЭЭФ огромную роль играют теплофизические свойства грунтов основания. К основным теплофизическим свойствам грунта принято относить теплопроводность λ, теплоем-

кость с и температуропроводность а.

Теплопроводность λ и теплоемкость с – это зависящие от температуры параметры, которые соединены в уравнении теплопроводности

39

(2.11) с температуропроводностью а. Этот параметр (а, м²/с) описывает глубину и скорость проникновения температурной волны в грунт.

Коэффициент пропорциональности, который связывает плотность теплового потока с температурным градиентом, известен как теплопроводность λ (Дж/(с·м·°С)). На λ сильно влияют содержание воды и плотность грунта. Замораживание значительно увеличивает теплопроводность, потому что для воды λвод = 0,57 Дж/(с·м·°С), а для льда

λл = 2,18 Дж/(с·м·°С).

Теплоемкость с (Дж/(кг·°С)) определяет количество энергии, запасенной в материале на единицу массы, на единицу изменения температуры. Теплоемкость не зависит от микроструктуры. Поэтому, в большинстве случаев, можно вычислить теплоемкость грунта из удельной теплоемкости различных ее составляющих, зная ихобъемное соотношение:

c = cгрxгр + cвозxвоз + cводxвод,

(2.30)

где xгр – удельный объем для твердой фазы, xгр = 1 – е; xвод – удельный объем для поровой воды, xвод = еS; xвоз – удельный объем для порового

воздуха, xвоз = е(1S) , где S – степень влажности, е – пористость грунта. Общая теплоемкость увеличивается от содержания воды и умень-

шается

в случае замерзания. Удельная теплоемкость льда cльда =

= 1884

Дж/(кг·ºС) по сравнению с теплоемкостью воды cвод =

= 4186

Дж/(кг·ºС).

Объемная теплоемкость CV выводится из удельной теплоемкости и объемной плотности грунта и представляет собой взвешенное среднее

арифметическое определенных компонентов грунта. Поэтому

 

 

Сv = pici xi ,

 

 

(2.31)

или

 

 

 

 

 

 

Сv

 

+ свод

 

w

 

= pгр cгр

 

 

,

(2.32)

 

 

 

 

 

100

 

где cгр – удельная теплоемкость минеральных компонентов (большинство минералов имеют cгр = 1000 Дж/(кг·K) при температуре 10˚C); cвод

– удельная теплоемкость воды (cвод= 4186 Дж/(кг·K)); w – влажность, %.

40