Министерство строительства предприятий тяжелой
индустрии СССР

Проектный и научно-исследовательский институт
«КРАСНОЯРСКИЙ ПРОМСТРОЙНИИПРОЕКТ»

 

ПРАКТИЧЕСКОЕ РУКОВОДСТВО
ПО ПРОИЗВОДСТВУ ТЕПЛОВЫХ РАСЧЕТОВ ОСНОВАНИЙ
В РАЙОНАХ С ВЕЧНОМЕРЗЛЫМИ ГРУНТАМИ

 

 

 

Красноярск

1976

 

В Руководстве освещаются практические методы определения температурных полей оснований, а также обуславливающих их составляющих внешнего теплообмена и теплофизических характеристик. Знание этих величин необходимо для проектирования фундаментов в соответствии с главой СНиП II-Б.6-66 «Основания и фундаменты зданий и сооружений на вечномерзлых грунтах».

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

Руководство написано канд. геогр. наук А.А. Коноваловым. В составления II раздела принимали участие канд. физ.-мат. наук А.Г. Конюшенко и инж. Л.А. Наумова.

ВВЕДЕНИЕ

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

Однако практика и дальнейшее развитие теории ставят перед проектировщиками новые задачи в области тепловых расчетов, прикладное решение которых раньше не являлось первой необходимостью. Например, для расчета несущей способности грунтового основания по существующей главе СНиПа II-Б.6-66 достаточно знать его стационарное (предельное, при стремлении времени к бесконечности) температурное поле. В то же время известно, что предельно-длительные сопротивления мерзлого грунта зависят от температуры, изменяющейся во времени. Но поскольку не было известно, как учитывать при определении сопротивлений изменение температуры во времени, не возникало потребности в расчете нестационарных температурных полей в основаниях зданий и сооружений. В настоящее же время, когда получено [1] уравнение для предельно-длительных сопротивлений, в которое температура входит как функция времени, появилась необходимость в разработке методов расчета нестационарных температурных полей.

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

Руководство состоит из 6 разделов с приложениями. В 1 разделе приводится перечень исходных и искомых величин, их условных обозначений. Для искомых величин указывается номер пункта, в котором излагается способ их определения. Следует отметить, что одни и те же переменные могут быть как исходными, так и искомыми величинами, это зависит от условий конкретной задачи. Поэтому в перечне не проводится грани между исходными и искомыми величинами.

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

Необходимость в определении стационарных полей возникает при установлении расчетных (наихудших) условий проектирования фундаментов на мерзлых грунтах, используемых по принципу 1 (СНиП II-Б.6-66); при вычислении предельной глубины оттаивания; при расчете температур под естественными источниками (водоемами, различными типами местности и т.д.).

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

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

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

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

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

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

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

Исследования института «Красноярский промстройниипроект» показали, что под центрами источников перечисленных конфигураций нестационарные температурные поля удовлетворительно подчиняются уравнению, выведенному Лахенбрухом [2] для центра круглого поверхностного источника.

                                  (1)

где q и q_пр - относительные (избыточные) температуры, соответственно нестационарная и стационарная;

 - число Фурье;

а - коэффициент температуропроводности;

y - глубина, считая от низа источника;

t - время.

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

Таблица 1

Значения q в зависимости от q_пр и F_o для центров полосы (числитель) и круга (знаменатель)

Однотипность нестационарных температурных полей под центрами источников различной формы и сравнительная простота описывающего их уравнения открывают большие возможности для решения задач с более сложными условиями однозначности. В частности, на основе уравнения (1) получено приближенное решение задачи типа Стефана для двух- и трехмерного случаев [4] (см. п.п. 3.6 - 3.7 Руководства), при единственном допущении, заключающемся в следующем: при отыскании производной глубины промерзания (оттаивания) по времени принято, аналогично [5] и по расчетам, что предельное температурное состояние в промерзающем (оттаивающем) грунте наступает практически мгновенно.

Основной исходной величиною в уравнении (1), отражающей влияние формы источника и нуждающейся в определении, является безразмерная стационарная температура q_пр.

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

Решения представлены в окончательном, удобном для практического использования виде (графиков, таблиц), с промежуточными выкладками можно ознакомиться в [2 - 6]. Основные расчетные формулы даны в приложении 1. Графики для стационарных температурных полей в системе «термосвая - грунт» (конвективный теплообмен жидкого теплоносителя в термосвае учитывается методом эквивалентной теплопроводности) и нестационарных температурных полей получены путем обобщения результатов расчетов на электро- и гидроинтеграторе.

Графики для определения положения фазовой изотермы во времени (рис. 7 - 11) построены при допущении а₂/а₁ = 1 и  (l₁, а₁, с₁ - коэффициенты тепло- и температуропроводности и теплоемкость в верхней зоне; l₂, а₂, с₂ - то же в нижней зоне), что близко к действительности. Уравнение для общего случая дано в приложении 1.

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

1. ИСХОДНЫЕ И ИСКОМЫЕ ВЕЛИЧИНЫ, ИХ УСЛОВНЫЕ ОБОЗНАЧЕНИЯ

Координаты, параметры (размеры в м, время в час)

х и у - соответственно, горизонтальная и вертикальная координаты пространства, отсчитываемые от центра координат, расположенного под краем теплового источника (или под его центром - в задачах раздела 3 и в решении для горизонтальной трубы в грунте в разделе 2, а также в формулах приложения 1);

r - радиальное расстояние от центра источника (для круга и цилиндра) до точки, в которой определяется температура;

у_ц - глубина под центром источника, отсчитывается от его низа (за вычетом мощности слоя сезонного промерзания и оттаивания - для поверхностных источников);

x - глубина промерзания или для бесконечного цилиндра - радиус промерзания;

b и l - ширина и длина источника;

r_o - радиус круга или цилиндра;

r_тр - радиус горизонтальной трубы;

r_из - внешний радиус теплоизоляционного покрытия;

h_из - толщина теплоизоляции;

h и h_тр - расстояние по оси у, соответственно, до низа заглубленного источника или до центра горизонтальной трубы;

h_тер - общая длина термосваи;

h_п, h_н - длина подземной и надземной частей термосваи;

d₂, d₁ - диаметры большей и меньшей труб, образующих круглую щель;

t - время;

t_г - годовой отрезок времени, 8760 час.

Показатели температуры (°С)

t_o - температура грунта за пределами источника (начальная температура - при расчете нестационарных температурных полей), устанавливаемая по наблюдениям за пределами зоны (или до начала) теплового влияния источника;

t_пр - предельная (стационарная) температура грунта в точке с заданными или искомыми координатами;

t¢_пр - предельная температура с учетом геотермического градиента;

t_t - нестационарная температура грунта в точке с заданными или искомыми координатами;

t_сг - среднегодовая температура в естественных условиях, определяемая разовым замером на глубине 10 м;

t_max - максимальная за год температура грунта;

t_п - температура ограничивающей поверхности (температура на контакте источника c грунтом), устанавливается согласно п.п. 4.1 - 4.6.

Климатические показатели

t_в - температура воздуха, град. (5.2)^х);

^х) Здесь и далее в скобках указан номер пункта, в котором изложена методика определения данного показателя.

А_м - годовая амплитуда колебаний температуры воздуха, град. (5.2);

В - радиационный баланс, ккал/м² мес (5.3);

LE - затраты тепла на испарение, ккал/м² мес (5.5);

Q_c = S + D - суммарная солнечная радиация, ккал/м² меc (5.3);

S и D - прямая и рассеянная солнечная радиация, ккал/м² мес (5.4);

А - альбедо поверхности, доли единицы (прил. III);

α_k - коэффициент конвективного (турбулентного) теплообмена, ккал/м² × час × град. (5.6);

Е - интенсивность испарения, мм/сутки, (5.5);

П_об - общая облачность, доли единицы, по данным метеостанций;

Р - количество осадков, мм/декада, по данный метеостанций;

V_в - скорость ветра на высоте флюгера, м/сек, по данным метеостанций;

h_сн, g_сн - высота и плотность снежного покрова, м, по данным метеостанций с учетом п. 5.7.

Физические и теплофизические константы

l, а, n_к и с - коэффициенты теплопроводности (ккал/м × час × град), температуропроводности (м²/час), кинематической вязкости (м²/час) и объемная теплоемкость; для грунтов, определяемые согласно п.п. 6.1 - 6.4, для прочих веществ и материалов - по справочникам;

С_ус, С_уз, С_ув - удельные теплоемкости частиц грунта, льда и воды (ккал/кг × град) (6.4);

g_ск - объемный вес скелета грунта, кг/м³;

g_у и g - удельный и объемный вес грунта, кг/м³;

W_c и W_н - суммарная влажность и влажность за счет незамерзшей воды, доли единицы (6.3);

q - степень заторфованности (отношение веса торфа к весу всего грунта), доли единицы (6.4);

v - скорость движения хладоносителя, м/час;

G - степень водонасыщения, доли единицы (6.4);

g_в и g_л - плотность воды и льда, кг/м³, (6.4)

W_п - полная влагоемкость, доли единицы (6.4);

l_э - эквивалентный коэффициент теплопроводности хладоносителя в термосвае, ккал/м × час × град, (6.5);

l_из, R_из - коэффициент теплопроводности, ккал/м × час × град и термическое сопротивление изоляции, м² × час × град/ккал, определяемые по справочному материалу;

l_сн - коэффициент теплопроводности снега, ккал/м × час × град (5.7);

R_сн - термическое сопротивление снега, м² · час × град/ккал (5.7);

G_r - геотермический градиент, принимаемый по наблюдениям, или приближенно 0,03 град/м;

r - скрытая теплота разовых превращений воды, 80 ккал/кг,

a - коэффициент теплоотдачи, ккал/м² × час × град (4.6).

Безразмерные обобщенные величины

 - предельная (стационарная) относительная температура, зависящая от формы источника;

 - нестационарная относительная температура;

 - число Фурье для источника любой формы, кроме бесконечного цилиндра;

 - число Фурье для бесконечного цилиндра (рис. 12);

 - число Фурье на графике рис. 7 - 11 (при К_о = 0; );

 - критерий Нуссельта для цилиндра;

 - то же, для кольцевой щели;

 - критерий Прандтля;

 - критерий Рейнольдса;

 то же, для кольцевой щели;

 - критерий Коссовича;

Индексами «М» и «Т» в тексте обозначены теплофизические константы и температуры соответственно мерзлого и талого грунтов.

Индексами «МИН» и «ОР» - удельные теплоемкости соответственно минеральных и органических частиц грунта.

Цифрами 1 и 2 помечены теплофизические константы и температура в мерзлой и талой зонах в случае промерзания и, наоборот, - в случае оттаивания.

2. ПРЕДЕЛЬНЫЕ (СТАЦИОНАРНЫЕ) ТЕМПЕРАТУРНЫЕ ПОЛЯ ГРУНТА

2.1. Предельные температурные поля рассчитываются по формулам:

- для теплофизически однородного (мерзлого или талого) грунта

t¢_пр = t_o + q_пр × (t_п - t_o) + G_r × y_п                                                 (2)

- для двухслойной среды из мерзлого и талого грунта;

в верхней (ближайшей к источнику) зоне

                                        (3)

- в нижней зоне

                                        (3¢)

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

2.2. Величина q_пр в выражениях (1) - (3¢) определяется с помощью формул в приложении 1 графиков на рис. 1 - 7 и в соответствии с указаниями п.п. 2.3 - 2.7.

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

- по формуле (4) вычисляется радиус изотермы r_ит;

                                                        (4)

где d - параметр равный:

 - для полосы и прямоугольника;

 - для круга;

 - для горизонтальной трубы

- описывается окружность радиусом r_ит (искомая изотерма) из точки на оси симметрии, отстоящей от у_ц в сторону источника на величину r_ит.

2.4. При действии нескольких (n) источников с одинаковой температурой ограничивающей поверхности q_пр вычисляется по формуле

q_пр = q_пр(1) + q_пр(2) + ... q_пр(n)                                                    (5)

Величины, q_пр(1) - q_пр(n) в уравнении (5) определяются по формулам приложения или согласно п. 2.3.

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

2.5. Если длина источника в два или более раз превышает ширину, а отношение ширины к глубине, на которой определяется температура, больше некоторой критической величины b/z^кр, зависящей от температурного перепада t_п - t_о (табл. 2) величину q_пр под центром источника с практической точностью можно принимать равной 1.

Таблица 2

Значение величины b/z^кр

Около границ таких источников q_пр определяется по формуле (6)

q_пр = 0,5 ± (1/p)arctg(x/y)                                                    (6)

или графику на рис. 1 как для полосы, но при d = (x/y)² (х - отсчитывается от контакта между участками с отрицательной и положительной температурой). Знаки перед вторым слагаемым подбираются так, чтобы под источником (участком с температурой t_п) q_пр лежала в пределах: 0,5 < q_пр <1; вне его 0 < q_пр < 0,5.

Рис. 1. Стационарные (предельные) температурные поля под центрами (х = 0,5b, z = 0) поверхностных источников.

В углах - величина q_пр рассчитывается по формуле 3 прил. 1 как для сектора круга с центральным углом b = 0,5p.

2.6. Предельные температурные поля заглубленных и приподнятых источников в форме бесконечного бруса, параллелепипеда, вертикальной трубы определяются также как и наземных (см. п.п. 2.2 - 2.5), только вместо X, У и 0,5b в соответствующие формулы и графики подставляются координаты вспомогательной полуплоскости 1 + x, n и 1 + xb, найденные в зависимости от Г и Н с помощью графика на рис. 3, при этом величина xb, определяется при Г = 0,5b/h и H = 0. При расчете температурных полей заглубленных источников используется нижняя часть графика рис. 3 (n > 0), приподнятых - верхняя (n < 0).

Центр координат размещается под низом левого края источника.

Рис. 2. Стационарные (предельные) температурные поля под центром (Х = 0) горизонтальной трубы в области у ³ r_тр + h_тр;

Рис. 3. Сетка для перевода координат X, У в координаты x, n и обратно.

2.7 На рис. 4 - 6 приведены графики для нахождения предельных температурных полей в системе «термосвая - грунт» для трех значений e = l_э/l_м на конец холодного периода.

Рис. 4. Стационарные (предельные) температурные поля в системе грунт - термосвая при x = 86.

Рис. 5. Стационарные (предельные) температурные поля в системе грунт - термосвая при x = 43.

Рис. 6 Стационарные (предельные) температурные поля в системе грунт - термосвая при x = 22.

Примечание. Размерные температуры в системе «грунт - термосвая» вычисляются по формуле (2), но вместо t_п в этом случае подставляется t_в - среднезимняя температура воздуха, вместо t_o - t_max.

Величины e могут изменяться от 1 до ¥. При e = 1 величина q_пр = 0,63q_пр.кр. (q_пр.кр. - предельная температура под источником в форме круга, определяемая согласно п.п. 2.2 - 2.3), при e = ¥ величина q_пр определяется как для вертикальной заглубленной трубы (п. 2.6). Для промежуточных значений e предельные температуры ищутся в такой последовательности:

- при заданных значениях переменных Н, m и h определяется q_пр для пяти значений e (0, 22, 43, 86, ¥)

- по этим данным строится график зависимости q_пр от exp(-0,01e)

- c помощью указанного графика находится q_пр для заданного значения e.

2.8. Наличие термического сопротивления между сооружением и грунтом учитывается приближенно - удлинением координат глубины под сооружением (или логарифма относительного радиуса в пределах высоты цилиндрического источника) на величину эквивалентного слоя, т.е. в соответствующие формулы или графики вместо у и h подставляется

                                                           (7)

                                                          (8)

3. НЕСТАЦИОНАРНЫЕ ТЕМПЕРАТУРНЫЕ ПОЛЯ ГРУНТА

3.1. Нестационарные температурные поля в грунтовом полупространстве при отсутствии фазовых превращений влаги (К_о = 0) рассчитываются по формуле

t_t = t_o + q(t_п - t_o) + G_r × y_п                                                   (9)

3.2. Относительная температура q в уравнения (9) для тепловых источников любой формы (за исключением бесконечного цилиндра) рассчитывается в такой последовательности:

- по п.п. 2.2 - 2.4, 2.6 находится q_пр в точке с заданными координатами и затем - глубина под центром источника в которой наблюдается найденное значение q_пр;

- на этой глубине по графикам на рис. 7 при К_о = 0 или, более точно, формуле (1), определяется величина q, которая и будет искомой.

3.3. При расчете температурных полей, возбуждаемых цилиндрическими источниками с m > 70 и (для горизонтального цилиндра расположенных глубже 10 м, в плоскостях, отстоящих от торцов более чем на 4r_о в сторону середины цилиндра, при определении q допускается использовать решение для бесконечного цилиндра (рис. 12).

Рис. 7. Нестационарные температурные поля под центрами источников в форме полосы, прямоугольника, круга, их заглубленных аналогов, а также - горизонтальной трубы в грунте при q_пр = 1 и q_пр = 0,1.

Рис. 8. Нестационарные температурные поля под центрами источников в форме полосы, прямоугольника, круга, их заглубленных аналогов, а также - горизонтальной трубы в грунте при q_пр = 0,9 и q_пр = 0,2

Рис. 9 Нестационарные температурные поля под центрами источников в форме полосы, прямоугольника, круга, их заглубленных аналогов, а также - горизонтальной трубы в грунте при q_пр = 0,8 и q_пр = 0,3.

Рис. 10. Нестационарные температурные поля под центрами источников в форме полосы, прямоугольника, круга, их заглубленных аналогов, а также - горизонтальной трубы в грунте при q_пр = 0,7 и q_пр = 0,4.

Рис. 11. Нестационарные температурные поля под центрами источников в форме полосы, прямоугольника, круга, их заглубленных аналогов, а также - горизонтальной трубы в грунте при q_пр = 0,6 и q_пр = 0,5.

Рис. 12. Нестационарные температурные поля источника в форме бесконечного цилиндра (К_о = 0).

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

Рис. 14. Графики дли определения К в зависимости от Б

3.4. Расчет температурных полей, формируемых несколькими источниками, а также источниками сложной формы, осуществляется по формуле (5), только вместо предельных относительных температур q_пр(1), q_пр(2) ... q_пр(n) подставляются нестационарные температуры q₁, q₂, ... q_n, найденные согласно п. 3.2. - 3.3.

3.5. До глубины около 10 м температуры грунта испытывают сезонные колебания. Расчет температур однородного (талого или мерзлого) грунта в годовом цикле колебаний в естественных условиях, а также под зданиями и наземными сооружениями, основания которых охлаждаются атмосферным воздухом, либо подземными сооружениями, расположенными выше отметки 10 м, производится с помощью графика на рис. 13, при этом в случае расчета температур под зданиями и сооружениями вместо t_сг - подставляется значение t_пр, вычисленное согласно п.п. 2.2. - 2.5.

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

Dt = (t_сг - t_max) × (t - t_o)/3                                                    (10)

здесь: (t - t_о) - время работы теплового источника, час.

3.6. Нестационарные температурные поля в полупространстве при фазовых переходах влаги рассчитываются по формулам:

                                                     (11)

                                                    (12)

q_1y, q_2y и q_1x, q_2x в формулах (11) - (12) определяются согласно п.п. 3.1 - 3.3 при соответствующих значениях теплофизических констант. В выражение F_o при вычислении q_y подставляется заданная глубина, а при вычислении q_x - глубина промерзания, если рассчитываются температуры промерзающего грунта, или оттаивания, если рассчитываются температуры оттаивающего грунта.

3.7. Глубина промерзания (оттаивания) под сооружениями различной формы определяется с помощью графиков на рис. 7 - 11 (а более точно, с помощью формулы в приложении 1).

Радиус промерзания (оттаивания) для источника в форме бесконечного цилиндра определяется в такой последовательности:

из соотношений (13) находятся параметры К и Б.

                                               (13)

- по графикам на рис. 14 определяются значения h = x/r_o и затем x.

В выражение для q при определении глубины или радиуса промерзания (оттаивания) вместо t_t подставляется температура начала замерзания грунтовой влаги t_нз, для песчаных грунтов, а в первом приближении также и для глинистых, принимаемая равной нулю. Более точно величина t_нз определяется по методике, изложенной в приложении IV.

3.8. Учет термической изоляции между тепловым источником и грунтом во всех задачах нестационарной теплопроводности осуществляется приближенно, согласно п. 2.10.

4. ТЕМПЕРАТУРА ОГРАНИЧИВАЮЩЕЙ ПОВЕРХНОСТИ

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

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

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

4.4. Температура поверхности участков естественной поверхности рассчитывается по формулам:

- среднегодовая температура

                                          (14)

- температура у поверхности грунта в теплый период года (с положительными температурами воздуха)

                                                       (15)

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

В правые части формул (14) и (16) подставляются значения входящих величин соответственно средних за год и средних за искомый период (но не менее, чем за декаду).

4.5. Цилиндрической ограничивающей поверхностью обладают различного рода инженерные сети, трубопроводы, а также охлаждающие или нагревающие устройства.

4.6. Температура цилиндрической поверхности вычисляется по формуле

                                                    (16)

где, кроме известных обозначений

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

q¢ - относительная температура, определяемая по табл. 3 в зависимости от критерия

Таблица 3

Зависимость q¢ от В_i

Величина a в выражении для В_i находится из выражений:

- для круглой трубы

                                            (17)

(l_x - коэффициент теплопроводности хладоносителя)

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

                                      (18)

Физические константы некоторых хладоносителей приведены в приложении II.

5. КЛИМАТИЧЕСКИЕ ПОКАЗАТЕЛИ

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

5.2. Температура воздуха и амплитуда ее колебаний устанавливаются по климатическим справочникам и ежегодникам. Среднегодовые температуры воздуха и грунта (левая часть формулы (14) показаны на схематической карте рис. 15.

5.3. Радиационный баланс за месяц вычисляется по формуле

В = q_с(0,9 - 0,7А)² - 0,7 ккал/см² мес.                                  (19)

Суммарная солнечная радиация определяется из выражения

Q_c = Q_o[1 - (1 - K_c)n_об]                                                   (20)

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

Величины Q_o и K_c находятся по табл. 4 и 5, а величина альбедо - по справочному материалу или приложению III.

5.4. Прямая солнечная радиация рассчитывается по формуле

                                                     (21)

5.5. Месячные суммы затрат тепла на испарение (ккал/см² мес) в естественных условиях определяются из выражения

LE = 1,8Е,                                                           (22)

где Е - интенсивность испарения, мм/сут., устанавливаемая с помощью рис. 16.

Таблица 4

Среднеширотные значения суммарной солнечной радиации при безоблачном небе Q_o ккал/см² мес

Рис. 15. Изолинии температур t_п и t_o на территории распространения вечномерзлых грунтов.

Рис. 16. Зависимость испарения от температуры и давления воздуха (е, мб).

Таблица 5

Среднеширотные значения К_с

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

                                               (23)

где Г = Р²/30(t_в + 2)²/256

t_в - среднемесячная температура воздуха

Q_c и D - суммарная и рассеянная солнечная радиация, ккал/см² мес.

5.6. Для перевода месячных сумм солнечной радиации и затрат тепла на испарение с размерностями ккал/см² мес, в систему размерностей ккал/м² час, принятую в технике, результаты, полученные в соответствии с п.п. 5.3 - 5.5, умножаются на коэффициент 1,39.

5.7. Среднегодовые (сг) и среднемесячные (см) значения a_к ккал/(м × час × град) вычисляются по формулам

                                                     (24)

                                                   (25)

Вычисленные по формуле (24) среднегодовые величины a_к показаны на схематической карте (рис. 17).

5.8. Термическое сопротивление снежного покрова м² · час × град/ккал рассчитывается с помощью соотношений

R_сн = h_сн/l_сн                                                         (26)

l_сн = 0,87g_сн + 0,018                                                 (27)

Рис. 17. Изолинии среднегодовые значений конвективного теплообмена.

Рис. 18. Изолинии значений R_сн на территории с вечномерзлыми грунтами

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

h_сн = 0,25(h₀ + h₁ + h₂ + 0,5h₃)                                            (28)

где h₀ - максимальная высота снежного покрова за предшествующую декаду (h₀ = 0, если определяется высота снежного покрова за первый зимний месяц),

h₁, h₂, h₃ - соответственно максимальные высоты снежного покрова за первую, вторую и третью декады месяца.

На рис. 18, 19 представлены схематические карты распределения величин R_сн рассчитанных по формулам (26) - (27) и  в области с вечномерзлыми грунтами.

6. ФИЗИЧЕСКИЕ И ТЕПЛОФИЗИЧЕСКИЕ ХАРАКТЕРИСТИКИ

6.1. Теплофизические свойства грунтов характеризуются коэффициентами тепло- и температуропроводности, а также удельной С_у или объемной С теплоемкостью. Эти величины связаны между собой соотношениями:

а = l/с; С = С_у × g                                                     (29)

6.2. Коэффициенты теплопроводности талого и полностью промерзшего грунта определяются по формуле

l = l^min + (l^max - l^min)G^j = l^min + (l^max - l^min)×e^jlnG                        (30)

где l^min, l^max - коэффициенты теплопроводности при G = 0 и G = 1, а j - безразмерный коэффициент, определяемые по табл. 7 в зависимости от параметра .

Удельная теплоемкость талого грунта С_ут рассчитывается по формуле

                                     (31)

Удельные теплоемкости компонентов грунтового комплекса, входящие в уравнения для теплоемкости грунта, определяется по табл. 6.

6.3 Коэффициент теплопроводности l¢_м и удельная теплоемкость промерзающего (оттаивающего) грунта С_ум рассчитываются по формулам:

l¢_м = l_т[1 + (l_м/l_т - 1) × (1 - W_н/W_c)]                                        (32)

             (33)

Влажность за счет незамерзшей воды и производная этой величины от температуры зависят от температуры (см. приложение IV), но поскольку в уравнения для температурных полей они входят в виде констант, приходится идти на упрощение - брать в расчёт средние значения W_н и  в температурном интервале значительных фазовых превращений. Среднее значение  в этом интервале для всех грунтов составляет минус 0,03 град.

Средние значения W_н принимаются по табл. 6.

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

6.4. Степень заторфованности грунтов и полная влагоемкость W_п в выражении для степени водонасыщения рассчитываются по формулам:

                                                           (34)

                                                  (35)

где А(м³/кг) и В (безразмерный) - параметры, зависящие от состава минеральной части заторфованного грунта, принимаемые соответственно равными 3460 и 1,3 - для песчаных заторфованных грунтов, 3260 и 1,9 - для глинистых заторфованных грунтов.

Удельные веса составляющих грунтового комплекса приведены в табл. 6.

6.5. Эквивалентный коэффициент теплопроводности хладоносителя в автономных охлаждающих устройствах (термосваях) вычисляется по формуле

                                             (36)

здесь, кроме известных обозначений

С_у, l_и - удельная теплоемкость хладоносителя, ккал/кг × град и его истинная теплопроводность, ккал/м × час × град.

g = 9,8 м/сек² - ускорение силы тяжести,

A_p - 1/427 ккал/кг × м - тепловой эквивалент работы,

j - отношение надземной части термосвай к подземной.

Рис. 19. Изолинии значений Dt_сн на территории с вечномерзлыми грунтами

Таблица 6

Удельные веса (кг/м³), теплоёмкости (ккал/кг × град), а также средние значения влажности за счет незамерзшей воды и параметра П (град^0,25) для типичных (незасоленных) грунтов (компонентов грунта)

Таблица 7

Значения l^min, l^max (ккал/м × час × град) и у в зависимости от параметра П

ПРИЛОЖЕНИЯ

I. Формулы для расчетов температурных полей

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

1) бесконечной полосы -

2) прямоугольника -

3) круга

где   

П(p/2, n, k) - полный эллиптический интеграл третьего рода (табулированная функция).

Для центра круга (b = 2p, r = 0) это уравнение значительно упрощается -

4) горизонтальной трубы в грунте -

б) нестационарные температурные поля грунта под центрами сооружений любой формы при отсутствии фазовых переходов -

5)

в) то же, при наличии фазовых переходов -

6) уравнение связи времени и глубины промерзания (оттаивания) под центром сооружений

где:

В выражение s¢ в отличие от s величина F₀ входит c множителем К_a.

Примечание: В формулах (5 и 6) глубина отсчитывается от низа сооружения, а его форма учитывается через величину q_пр

II. Физические и теплофизические свойства некоторых хладоносителей

Таблица 8

Примечание: В таблице 8 приведены средние значения в интервале температур от -10 °С до -50 °С.

III. Значения альбедо (%) некоторых поверхностей

IV. Определение содержания незамерзшей воды в мерзлом грунте в зависимости от температуры

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

                                                                  (*)

где Т = t/t_нз

t - температура грунта,

t_нз - температура начала замерзания свободной воды в грунте, определяемая по формуле:

                                                         (**)

где Р - концентрация солей в грунтовой (поровой) влаге,

Р_э и t_э - эвтектические значения концентрации и температуры, определяемые по табл. 9.

Таблица 9

Производная величины W_н по температуре, которая входит в формулу для расчета эффективной теплоёмкости, равна

                                                (***)

V. Примеры расчета

а) к разделу 2

1. Определить стационарную температуру грунта под длинным холодильником в форме заглубленного бруса при следующих исходных данных: t_п = -5, t₀ = -1, х = 3 м, у = 4 м, h = 6 м, b = 12 м (помнить, что начало системы координат в левом углу заглубления).

Находим безразмерные переменные

Г = 3/6 = 0,5; H = 4/6 = 0,67; 0,5В/h = 1

С помощью рис. 3 по найденным значениям Г, H определяем значения

x = 3,8; xb = 5,7; n = 2,5

В формулу (1) приложения 1 вместо х подставляем 1 + x = 4,8 вместо 0,5b 1 + xb = 6,7, вместо у n = 2,5 тогда имеем

По формуле (2) основного текста вычисляем t_пр

t_пр = -1 + 0,76 · (-4) + 10 · 0,03 = -3,74

2. Определить стационарную температуру в системе «грунт - термосвая» при следующих исходных данных: h_п = 4 м, r₀ = 0,1 м, y = 1 м, r = 0,5 м, t_в = -21, t₀ = -1, e = 86.

Определяем безразмерные переменные

m = 40; h = 5, Н = 0,25.

По графику на рис. 4 находим искомую предельную температуру q_пр = 0,12.

По формуле (2) вычисляем

t_пр = -1 + 0,12 × (-20) + 1 × 0,03 = -3,37

б) к разделу 3.

3. Определить нестационарную температуру грунта при исходных данных примера 1 и а_м = 0,004 м²/час; t = 1200 час.

По графику на рис. 1 определяем относительную глубину, на которой q_пр = 0,76 наблюдается под центром источника на вспомогательной плоскости с координатами x и n.

С учетом того, что xb = 5,7 (см. пример 1), величина n равна 6,7 × 0,35 = 2,35.

С помощью рис. 3 переводим величину n = 2,35 (при xb = 5,7) в у_ц; у_ц/h = 0,67, у_ц = 4.

Определяем F₀:

По графику на рис. 9 при К₀ = 0, q_пр = 0,76 ≈ 0,8 и F₀ = 0,3 находим: q = 0,23.

По формуле (9) вычисляем t(t):

t(t) = -1 + 0,23 × (-4) + 10 × 0,03 = -1,62

4. Определить время, за которое нулевая изотерма достигнет глубины у = 4 м под центром здания (b = 12 м, L = 19 м) c холодным подпольем (t_г = -3,5); необходимые для расчета исходные характеристики грунта: t₀ = +0,5; а_м = 0,004 м²/час; W_c = 0,1; g_ск = 1600 кг/м³; С_м = 400 ккал/м³ град.

С помощью графика на рис. 1 в зависимости от у/0,5b = 0,67 находим величину q_пр: q_пр = 0,6.

Вычисляем К₀ = 80 × 0,1 × 1600/(400 × 4) = 8 и q = (0 - 0,5)/(-3,5 - 0,5) = 0,125.

По графику на рис. 11 находим P₀ = 6, откуда

t = 6 × 16/0,004 = 24000 час ≈ 2,7 года

в) к разделу 4.

5. Определить температуру поверхности скважины (d₂ = 0,4), охлаждаемой холодным воздухом через трубу меньшего диаметра (d₁ = 0,1 м), коаксиально расположенную в скважине; l_м = 1 ккал/м × час × град.

V = 5 м/сек, t_ср = -10, t₀ = -1.

Из формулы (18), подставив в нее указанные значения V, d₂ и d₁, а также значения n и а для воздуха, по приложению II находим величину коэффициента теплоотдачи a = 13,5 ккал/м² × час × град.

Вычисляем критерий Био: Bi = 13,5 × 0,2/1 = 2,6

По табл. 3 находим значение q¢: q¢ = 0,86

По формуле (16) вычисляем температуру поверхности скважины

t_п = -1 + 0,86(-10 + 1) = -8,7

г) к разделу 5.

6. Определить величину радиационного баланса бетонного (светло-серого) покрытия за июль на 70° широте при общей облачности n_об = 0,6.

По таблицам 4 и 5 определяем Q₀ и К_с, а по приложению А:

Q₀ = 20,3 ккал/см² × мес; К_с = 0,5; А = 0,23.

По формуле (20) вычисляем Q_c:

Q_c = 20,3[1 - 0,5 × 0,6] = 14,2 ккал/см² мес.

По формуле (19) находим величину радиационного баланса

В = 14,2(0,9 - 0,16)² - 0,7 = 7,1 ккал/см² мес

7. Определить величину Dt_cн в районе Якутска

По схематической карте (рис. 19) находим Dt_cн = 4,8.

д) к разделу 6.

Определить теплофизические характеристики заторфованной глины со следующими характеристиками W_c = 0,43; g_узг = 2,320 кг/м³, g = 1620 кг/м³; t = -2, t_нз = -0,5.

С помощью формул (34), (35) и по приложению IV определяем исходные обобщенные характеристики:

g = 3260 : 2320 - 1,19 = 0,21

G = 0,43/0,48 = 0,89;

По формуле (*) приложения IV рассчитываем влажность за счет незамерзшей воды: W_п = 0,27.

Определяем параметр  и по табл. 6, в зависимости от величины этого параметра, находим l^min = 0,19; l_т^max = 1,24; l_м^max = 1,9; у = 1,27.

По формуле (30) вычисляем l_т и l_м.

l_т = 1,1 ккал/м × час × град; l_м = 1,6 ккал/м × час × град.

По формуле (32) рассчитываем коэффициент теплопроводности при t = -2.

l = 1,11 + (1,45 - 1) × (1 - 0,63) = 1,28 ккал/м × час × град.

По формулам (30), (33) с учетом (***), табл. 6 и выражения для объемного веса скелета грунта g_ск = g/(1 + W_c) определяем эффективную объемную теплоемкость

С = 1140[0,22 × 0,79 + 0,38 × 0,21 + 0,5 × 0,43 - 0,27 × (-0,5) + 0,02 × 80] = 2508 ккал/м³ · град.

ЛИТЕРАТУРА

1. Вялов С.С. Длительная прочность мерзлых грунтов при переменной нагрузке и температуре. Труды V Всесоюзного совещания-семинара по обмену опытом строительства в суровых климатических условиях. Т. 2, вып. 5. Тюмень, октябрь, 1968. Красноярск, 1968. (Институт «Красноярский промстройниипроект»).

2. Lachenbruch A.H. Three-dimensional heat conduction in permafrost beneath heated buildings. - Jeol. Surv. Bull. № 1052-В, 1957.

3. Лыков А.В. Теория теплопроводности. М., «Высшая школа», 1967.

4. Коновалов А.А. Метод определения температурных полей промерзающих или оттаивающих оснований зданий и сооружений. В сб.: Строительство в районах Восточной Сибири и Крайнего Севера, вып. 26, Красноярск, 1973. (Институт «Красноярский промстройниипроект»).

5. Порхаев Г.В. Тепловое взаимодействие зданий и сооружений с вечномерзлыми грунтами. М., «Наука», 1970.

6. Конюшенко А.Г., Коновалов А.А. К расчету стационарных температурных полей заглубленных сооружений. В сб.: Строительство в районах Восточной Сибири и Крайнего Севера, вып. 26, Красноярск, 1973. (Институт «Красноярский промстройниипроект»).

7. Иванов Н.С. Тепло- и массоперенос в мерзлых горных породах. М., «Наука», 1969.

8. Пособие по проектированию оснований и фундаментов зданий и сооружений на вечномерзлых грунтах. М., Стройиздат, 1969.

9. Руководство по расчетным способам определения прочностных деформативных и теплофизических характеристик мерзлых, промерзающих и оттаивающих грунтов. Красноярск, 1972. (Институт «Красноярский промстройниипроект»).

10. Будыко М.И. Тепловой баланс поверхности земли. Гидрометеоиздат. Л., 1956.

11. Константинов А.С. Испарение в природе. Гидрометеоиздат. М., 1963.

12. Кутателадзе С.С., Боришанский В.М. Справочник по теплопередаче. Госэнергоиздат. Л., - М. 1959.

13. К расчету температурных полей промерзающих и оттаивающих грунтовых оснований зданий и сооружений. В сб. «Строительство в районах Восточной Сибири и Крайнего Севера», вып. 29, Красноярск, 1974. (Институт «Красноярский промстройниипроект»).

 

СОДЕРЖАНИЕ