CC BY
31
УДК 625.711.84+625.31 А.С. Миляев
С.-Петербургский государственный лесотехнический университет имени С.М. Кирова
Миляев Александр Сергеевич родился в 1936 г., окончил в 1959 г. Ленинградское высшее военное инженерно-техническое училище ВМФ, доктор технических наук, профессор, заслуженный работник высшей школы РФ, зав. кафедрой теоретической и строительной механики С.-Петербургского государственного лесотехнического университета. Имеет более 150 научных работ в области механики деформируемого твердого тела, в том числе механики силового взаимодействия конструкций и сооружений с грунтом при статических и динамических нагрузках. E-mail: [email protected].
ПРОГНОЗИРОВАНИЕ ПАРАМЕТРОВ ПРОМЕРЗАНИЯ СЛОИСТЫХ ОСНОВАНИЙ ЗИМНИХ ЛЕСОВОЗНЫХ ДОРОГ ПО ТЕКУЩИМ ЗНАЧЕНИЯМ МЕТЕОДАННЫХ
Приведена методика расчета глубины и скорости промерзания слоистых оснований зимних лесовозных дорог на базе фактических метеоданных. Методика построена на численном решении нестационарной задачи о теплопередаче в термодинамической системе «воздух - грунт» с учетом фазовых превращений поровой воды в лед.
Ключевые слова: зимние лесовозные дороги, теплопередача, промерзание оснований, фазовые превращения воды в лед.
При проектировании и эксплуатации зимних лесовозных дорог необходимо располагать информацией о распределении температуры по глубине естественного основания. Эта информация позволяет определить толщину мерзлого слоя грунта и выполнить прочностные расчеты дорожной конструкции.
Для расчета глубины и скорости промерзания оснований зимних дорог требуются сведения о состоянии атмосферы (метеоданные), о геологическом строении основания и теплофизических характеристиках грунтов в талом и мерзлом состояниях, температуре на поверхности основания в начале расчета, а также об условиях теплообмена на подошве основания.
К настоящему времени создана научная база таких расчетов [1, 7, 8, 22], разработаны руководства и рекомендации [16, 18], стандарты [2-5, 20], методики [6, 13, 14, 21], однако нормы по расчету глубины и скорости промерзания слоистых грунтовых оснований зимних лесовозных дорог на основе текущих значений метеоданных отсутствуют.
В практических расчетах используют методики, построенные на полуэмпирических формулах решения задачи о теплопередаче от атмосферы к грунту, которые базируются на статистических метеоданных, осредненных по многолетним наблюдениям [6, 16].
© Миляев А.С., 2012
Й vjl
В работе [16] глубину сезонного промерзания неосушенных торфяников рекомендуется определять по формуле, построенной на базе аналитического решения известной задачи Стефана [10], в которой реальные аргументы заменены на среднестатистические, что в конечном итоге приводит к отрыву математической модели [16] от конкретных условий протекания процесса промерзания. В частности, по методике [16] нельзя установить зависимость изменяющегося во времени распределения температуры по глубине грунта при изменении температуры воздуха, теплопроводности и теплоемкости в зависимости от температуры слоев грунта.
Кроме того, эта методика не учитывает теплоприток от слоев грунта ниже подошвы основания в зависимости от их термодинамического состояния.
В работе [6] время, необходимое для промерзания слоя грунта заданной толщины, рекомендуется определять по формуле, базирующейся на широко известном решении одномерной стационарной задачи о передаче тепла через стенку [9], причем температура воздуха на поверхности дороги зимой отсчи-тывается от среднестатистической нулевой температуры за 10 лет. Это означает, что независимо от реальных текущих значений метеоданных время промерзания грунта должно отсчитываться от среднестатистической даты наступления холодов. Количество холода, необходимого для промерзания слоя грунта заданной толщины, определяется в [6] по теплоемкости массы грунта в слое при его температуре - 0,1 °С , что не соответствует физическим условиям. Так, согласно литературным данным [22], поровая вода в слоях торфа, расположенных на некоторой глубине от поверхности, замерзает при температуре -0,3...-0,4 °С, а переход в лед связанной воды, заполняющей внутренние поры растительных волокон и их остатков, происходит при температуре ниже -1,0 °С.
Методика [6] не учитывает суточные изменения температуры воздуха в расчетном отрезке времени, теплоту фазовых переходов воды в лед, зависимость теплопроводности и теплоемкости грунта от температуры, а также теп-лоприток от слоев грунта ниже подошвы основания в зависимости от их термодинамического состояния.
В учебном пособии [13] для теплотехнических расчетов оснований зимних лесовозных дорог используется методика [6].
Реальными суточными метеоданными (в противоположность среднестатистическим) будем называть текущие значения температуры и влажности воздуха, скорость ветра, толщину и температуру снежного покрова или поверхности грунта, измеренные в данном месте в течение расчетного интервала времени.
Методику прогнозирования глубины и скорости промерзания оснований сезонных зимних лесовозных дорог разработаем, используя реальные суточные метеоданные и численное решение нестационарной задачи о теплопередаче в термодинамической системе «воздух - грунт» с учетом фазовых превращений поровой воды в лед и зависимости теплофизических свойств грунта от температуры.
Объекты и методика исследования
Объектами исследования служат слоистые основания зимних лесовозных дорог. Методика исследования носит теоретический характер и состоит в следующем. Рассматривается замкнутая термодинамическая система (рис. 1), состоящая из слоев воздуха и грунта с известными физико-механическими и теплофизическими свойствами; начальная температура и фазовое состояние системы заданы; температура воздуха изменяется по заданному закону. Поверхности у = 0, у = Ъ, х = Н4 теплоизолированы. Требуется определить изменение температуры в каждой точке системы в зависимости от изменения температуры воздуха и заданных условий на ее границе.
Распространение теплоты в слоях грунта подчиняется уравнению теплопроводности [9]:
8 t дх Iх дх
где Н - энтальпия, Дж/м3
^ - время, с; Т - температура, °С;
_8_
8у
Я.
8Т ' 8у
8 (у 8 Т Л
+ & Iк * I + ^
(1)
\х ,
, - теплопроводность в направлении осей х, у, г, Вт/м ; - мощность внутренних источников энергии, Вт/м3. В правой части уравнения (1) теплопроводность в направлении осей х, у, г зависит от температуры. Тепловые потоки через поверхности у = 0, у = Ъ, х = И4 отсутствуют: Q(x, у, г, ^ = 0. Приращение энтальпии
Т
АН = | р С(Т)ёТ,
Ть
где ТЪ, Т - начальное и текущее значение температуры, °С; р - плотность грунта, кг/м3; С - удельная теплоемкость, Дж/(кг-0С).
(2)
Рис. 1. Общая расчетная схема: 1 - 4 - слои грунтового основания; Н1 - Н4 - отметки подошв слоев; Та1г - температура воздуха,
Тагт Тагг(^')
Внутренних источников энергии в системе нет (qV = 0), однако в слоях грунта при замерзании поровой воды выделяется теплота кристаллизации (Qgr, Дж/м3), которую необходимо учесть в расчете. В соответствии со СНиП 2.02.04 - 88 [20] Qgr принимается равной количеству теплоты, необходимой для замерзания воды в единице объема грунта, и определяется по следующей формуле:
Qgr = Qw(Wgr - wm)ps (3)
где Qw - теплота кристаллизации воды, Qw = 3,35 10 Дж/кг;
Wgr - суммарная влажность грунта, доли единицы;
Wm - влажность мерзлого грунта между включениями льда, доли единицы;
ps - плотность скелета грунта, кг/м3.
Примем, что Qgr выделяется в интервале температур от —1 до 0 °С; Wm = 0.
На границе x = 0 происходит конвективный теплообмен:
Q = aA(Tair - Tgr), (4)
где Q — тепловой поток через границу x = 0;
а — коэффициент теплоотдачи;
A — площадь теплообмена;
Tgr — температура грунта на его поверхности.
Коэффициент а на границе «воздух-снег» (snw) и «воздух-лед» (ice) можно вычислить по следующим эмпирическим формулам СНиП [19]:
а snw = 23^/7+0,3 ; (5)
^ snw =
6^/v + 0,3, (6)
где v — скорость ветра над поверхностью теплообмена, м/с.
Уравнение (1) — нелинейное, нестационарное. Его решение возможно только численным методом. Для обеспечения устойчивости численной процедуры необходимо соблюдать определенное соотношение между «шагами» по времени At и пространству Ah [17]:
^Aí 1
-7 ^-, (7)
pC (Ah)2 2
где X — теплопроводность в направлении распространения тепловой волны;
Ah - приращение пространственной координаты в направлении распространения тепловой волны.
Сформулированная задача может быть решена методом конечных элементов, например, с помощью программного комплекса «ЗЕНИТ».
Пример теплотехнического расчета
Основанием зимней дороги служит слоистый грунт, состоящий из пластов торфяного грунта (1), суглинка (2), супеси (3) и песка (4).
Ограничимся определением температуры на одной вертикали основания вдоль трассы дороги. С этой целью на оси трассы «вырежем» грунтовую колонку с размерами в плане b*b и глубиной h4: b = 1,0 м; h4 = 2,4 м.
Используем прямоугольную систему координат, как показано на рис. 2.
Физико-механические и теплофизические свойства грунта принимаем по СНиП 2.02.04-88 [20].
Верхний слой (1) - торфяной грунт с плотностью скелета = 400 кг/м3; толщина слоя - 0,5 м; относительная влажность торфяного грунта = 2,00 (в долях плотности скелета грунта); теплопроводность талого торфяного грунта Хй = 0,93, мерзлого - 1,39 Вт/(м°С); объемная теплоемкость талого торфяного грунта Сл = 3,78 106, мерзлого - 2,73 • 106 Дж/(м3°С).
Торфяной грунт подстилается слоем суглинка (2) с плотностью скелета = 1400 кг/м3; толщина слоя - 0,5 м; относительная влажность суглинка = 0,05; теплопроводность талого суглинка Хй = 0,46, мерзлого - 1,39 Вт/(м°С); объемная теплоемкость талого суглинка =1,6 106, мерзлого - 1,47106 Дж/(м3°С).
Третий слой (3) - супесь с плотностью скелета = 1400 кг/м3; толщина слоя - 0,7 м; относительная влажность = 0,10; теплопроводность талой супеси Хй = 0,93, мерзлой - 1,05 Вт/(м°С); объемная теплоемкость талой супеси Са = 1,89-106, мерзлой - 1,74^ 106 Дж/(м3°С).
Четвертый слой (4) - песок с плотностью скелета = 2000 кг/м3; толщина слоя - 0,7 м; относительная влажность = 0,1; теплопроводность талого песка Х^ = 2,73, мерзлого - 2,90 Вт/(м°С); объемная теплоемкость талого песка Са = 2,68 106, мерзлого - 2,26-106 Дж/(м3°С).
Плотность грунта в слоях считаем постоянной величиной.
Обозначим теплопроводность, теплоемкость и энтальпию при температуре 0 °С символами Х0 , С0 и Н0. Тогда зависимости теплопроводности Х, теплоемкости С и энтальпии Н от температуры можно представить в виде
Рис. 2. Грунтовая колонка к примеру расчета
Х = ^/(Т); С = С0 ¿(Г); Н = Н0/(?).
Введем безразмерные функции
X = £ = т, С = С = /2(Г); Н = Н = А(П
(8)
(9)
Безразмерные функции (9) представлены на рис. 3.
Приращение энтальпии, определяемое выражением (2), для каждого слоя грунта вычисляем при начальном (Тъ = -15 °С) и текущих значениях температуры (Т_1 = -1 °С; Т0 = 0 °С; Т10 = 10 °С), промежуточные - по программе линейной интерполяции.
В табл. 1 приведены промежуточные значения приращения энтальпии грунтовых слоев по рис. 3.
Таблица 1
Грунт АН10-6, Дж/м3, при Т, °С
-15 - 1 0 10
Торф 0 38,20 306,20 344,00
Суглинок 0 20,58 44,00 60,00
Супесь 0 24,30 71,20 90,10
Песок 0 31,64 98,64 125,44
Рис. 3. Зависимости теплопроводности X (1), теплоемкости С (2) и энтальпии Н (3) от температуры Т
Граничные условия: температура воздуха Tair = -16 °С (формула (4)); на прямых у = 0 и у = 1 м тепловые потоки Q = 0. В подошве основания, на прямой х = 2,4 м, тепловой поток Q Ф 0, поэтому поддерживается постоянная температура Т^п = 5 °С.
Скорость ветра у поверхности торфа принимаем равной 5 м/с. Полагаем, что на границе «воздух-торфяной грунт» при температуре Т = -0,5 °С образуется ледяная пленка, поэтому коэффициент теплоотдачи аСе можно вычислить по формуле (6) СНиП 2.06. 04-82 [19]:
асе = 675 + 03 = 13,8 Вт/(м2 -°С). (10)
Начальные условия: при ^ = 0 на поверхности основания, т.е. на прямой х = 0 м, температура Т = -1,5 °С, во всех остальных точках основания, кроме подошвы, Т = -0,5 °С.
Ставится задача: найти распределение температуры Т по толщине мерзлого слоя грунта в момент времени ¿120 = 120 ч.
Расчет выполняется в одномерной постановке, так что в правой части уравнения теплопроводности (1) второе и третье слагаемые равны нулю. Приращение пространственной координаты в направлении распространения тепловой волны Ах полагаем 0,01 м, поэтому в направлении оси ох геометрическая модель имеет 241 узел.
Используя условие устойчивости численной процедуры (7), вычисляем допустимые шаги по времени А( для каждого слоя грунта, результаты вычислений сводим в табл. 2.
6*
83
Таблица 2
Грунт At, с, при Т, °С
-15...-1 0.10
Торф 49,1 101,6
Суглинок 70,7 87,0
Супесь 41,4 50,8
Песок 19,5 24,5
В расчете использованы: среднее значение шага по времени At = 15 с; минимальное — Atmm = 10 с; максимальное - Atmax = 30 с.
Для решения сформулированной задачи применяется метод конечных элементов, шаговое по времени решение и итерационная процедура приближения на каждом шаге. Точность решения на каждом шаге по времени контролируется в программном комплексе «ЗЕНИТ» с помощью итерационной процедуры по двум показателям: по стабилизации количества тепла, протекающего через конечный элемент, и по стабилизации температуры.
Аналогичная схема решения нелинейной нестандартной задачи теплопроводности была использована автором в статье [12].
Машинное время счета одного варианта задачи на персональном компьютере типа Pentium (R) с двумя центральными процессорами несущей частоты 2,6 ГГц и оперативной памятью 2 ГБ составляет около 20 мин. В результате решения задачи на каждом шаге по времени получаются значения температуры в узлах конечных элементов.
Поскольку в данном случае необходимо иметь распределение температуры в грунте по прошествии 120 ч с момента установления температуры воздуха Tair = —16 °C, в файл выходных данных записаны результаты расчета только на последнем шаге по времени tend = 120 ч.
По этим данным построены соответствующие графики зависимости температуры T120 от координаты х в узлах на прямой y = 0 м вблизи от поверхности грунтового основания дороги и у подошвы этого основания (рис. 4).
0 0.083 0,167 0,250 0.333 0.417 л. м 7",,
-3,75
-7.50
-11.25
/
/ /
hfr = 0,24 м
3,6«
2.40
1,20
-1.20
/
/ 7
/
/
40 1.60 1,80 2.00 2.20 .v. м Л,„ " 0,5 м к--—--?
Рис. 4. Распределение температуры по глубине грунта у поверхности (а) при
Tair = -16 °C; tend = 120 ч и у подошвы (б) при Tdwn = 5 °C; tend = 120 ч
Символом кг на рис. 4, а отмечена нижняя граница промерзания грунта. При этом принято, что торфяной грунт замерзает при температуре -0,3 °С. Толщина мерзлого торфяного грунта к^ , полученная в расчете, составляет около 29 см. Скорость промерзания торфяного грунта равна 29/5 = 5,8 см/сут.
Символом кг на рис. 4, б отмечена толщина незамерзшего песка вблизи от подошвы грунтового основания. Она составляет 0,5 м.
Заключение и выводы
Разработана методика автоматизированного расчета процесса промерзания грунта естественного слоистого основания зимних лесовозных дорог, в которой учитываются фактические (не усредненные за предыдущие годы) метеорологические данные в заданном интервале времени, фазовые превращения поровой воды в лед и зависимость теплофизических свойств грунта от температуры.
Методика базируется на численном решении нестационарного нелинейного уравнения теплопроводности, в котором теплопроводность, теплоемкость и энтальпия слоев грунта изменяются в зависимости от температуры этих слоев.
Эту методику можно использовать для решения теплотехнических задач при оценке несущей способности оснований сезонных зимних лесовозных дорог.
Приведенный пример расчета промерзания слоистого грунта естественного основания зимней дороги показывает, что скорость промерзания верхнего слоя основания из торфяного грунта при достаточно низкой стабильной температуре воздуха (-16 °С в течение 5 сут) может достигать 5,8 см/сут, что обеспечивает значительную толщину (29 см) мерзлого слоя.
В примере расчета влажность мерзлого грунта между включениями льда Жт принята равной нулю, поскольку, с одной стороны, ее можно определить только опытным путем после промерзания грунта, с другой - отличная от нуля Жт может только увеличивать теплоту замерзания грунта Qgr в формуле (3), что может уменьшить скорость промерзания. Следовательно, в приведенном примере расчета получена нижняя граница скорости промерзания грунта.
Для оценки достоверности рассчитанной в примере скорости промерзания верхнего слоя основания из торфяного грунта при температуре воздуха -16 °С в течение 5 сут, равной 5,8 см/сут, воспользуемся результатами расчетов толщины промерзания торфа по методике работы [6]. В табл. 15 этой работы [6] приведены следующие результаты расчетов промерзания зимней дороги на болоте, покрытой снегом: влажность торфа - 1000 %; средняя температура воздуха -12° С в течении 9 сут; средняя толщина снежного покрова -6 см; толщина промерзания торфа - 30 см. Таким образом, средняя скорость промерзания торфа 30/9= 3,33 см/сут.
Учитывая, что в нашем примере расчета отсутствует снежный покров и температура воздуха принята ниже, чем в работе [6], степень достоверности приведенной в статье методики расчета можно оценить как удовлетворительную.
При известных механических характеристиках и толщине мерзлого слоя напряженно-деформированное состояние конструкции сезонных зимних лесовозных дорог можно определить по методике [11].
СПИСОК ЛИТЕРАТУРЫ
1. Вялов С.С. Реология мерзлых грунтов / Под ред. В.Н. Разбегина М.: Строй-издат, 2000. 464 с.
2. ГОСТ 24847-81. Грунты. Методы определения глубины сезонного промерзания. М.: Госстандарт, 1981. 10 с.
3. ГОСТ 25358-82. Грунты. Метод полевого определения температуры. М.: Госстандарт, 1982. 10 с.
4. ГОСТ 26263-84. Грунты. Методы лабораторного определения теплопроводности мерзлых грунтов. М.: Госстандарт, 1984. 9 с.
5. ГОСТ 5180-84. Грунты. Методы лабораторного определения физических характеристик. М.: Госстандарт, 1985. 24 с.
6. Дорофеев А.Г., Дорофеев В.Н. Зимние лесовозные дороги с увеличенным сроком действия // Лесоэксплуатация и лесосплав. Обзор. информ. Вып. 10. М.: ВНИПИЭИлеспром, 1987. 52 с.
7. Ершов Э.Д. Общая геокриология. М.: Недра, 1990. 559 с.
8. Инженерная геокриология: справ. пособие / Под ред. Э.Д. Ершова М.: Недра, 1991. - 439 с.
9. Исаченко В.П., Осипова В.А., Сукомел А.С. Теплопередача М.: Энергоиз-дат, 1981. 416 с.
10. ЛыковА.В. Теория теплопроводности. М.: Высш. шк., 1967. 600 а
11. Миляев А.С. Автоматизированный расчет конструкций зимних лесовозных дорог: учеб. пособие. СПб: Изд-во СПбГЛТУ, 2006. 303 с.
12. Миляев А.С. Влияние тепловыделений биомассы подстилающего слоя на промерзание оснований зимних лесовозных дорог // Лесн. журн. 2010. №1. С. 53-58. (Изв. высш. учеб. заведений).
13. Морозов В.С. Теплотехнический расчет оснований сезонных зимних дорог на болотах: учеб. пособие. Архангельск: Изд-во АГТУ, 2004. 79 а
14. Невзоров А.Л. Фундаменты на сезоннопромерзающих грунтах: учеб. пособие. Архангельск: Изд-во АГТУ, 1999. 156 а
15. Павлов А.В. Расчет и регулирование мерзлотного режима почвы. Новосибирск: Наука, 1980. 240 а
16. Рекомендации по определению теплофизических характеристик торфяных грунтов и расчетам их промерзания и оттаивания. М.: НИИОПС Госстроя СССР, 1978. 54 с.
17. Рихтмайер Р., Мортон К. Разностные методы решения краевых задач. М.: Мир, 1972. 418 с.
18. Руководство по определению физических, теплофизических и механических характеристик мерзлых грунтов / ПНИИИС Госстроя СССР. М.: Стройиздат, 1973. 191 с.
19. СНиП 2.06.04-82. Нагрузки и воздействия на гидротехнические сооружения (волновые, ледовые и от судов) / Госстрой России. М.: ГУП ЦПП, 1982. 37 с.
20. СНиП 2.02.04-88. Основания и фундаменты на вечномерзлых грунтах / Госстрой России. М.: ГУП ЦПП, 2000. 52 с.
21. Строительство промысловых сооружений на мерзлом торфе / С.С. Вялов, Г.Л. Каган, А.Н. Воевода, В.И. Муравленко. М.: Недра, 1980. 144 с.
22. Цытович Н.А. Механика мерзлых грунтов (общая и прикладная). М.: Высш. шк., 1973. 448 с.
Поступила 11.01.11
A.S. Milyaev
St. Petersburg State Forest Technical University
Forecasting of Freezing Parameters of Layered Bases of Winter Wood-Roads According to the Current Meteorological Data
The method of calculating the depth and speed of freezing of layered bases of winter wood-roads according to the current meteorological data is introduced in the article. The method is based on a numerical solution of a non-stationary problem of heat transfer in a thermodynamic "air-ground" system with an allowance for phase transitions of interstitial water into ice.
Key words: winter wood-roads, heat transfer, freezing of bases, phase transitions of water into ice.
