Полный текст

Введение В настоящее время в связи с появлением современных зондов с датчиками порового давления появилась возможность оценить фильтрационные характеристики грунтов при проведении статического зондирования [1-6]. Однако на сегодняшний день отсутствует апробированная методика определения коэффициента фильтрации по данным статического зондирования. В зарубежной практике данные зонды активно применяются, но существующих методик расчета также недостаточно. К наиболее известным и распространенным относятся методика Parez и Fauriel (1988) и методика Robertson (1992) [7-9]. Данные методики расчета основаны главным образом на использовании диссипационных тестов и носят преимущественно эмпиричный характер с привязкой к конкретным разновидностям грунтов и региону их залегания. В связи с этим они имеют весьма ограниченное применение и не могут использоваться на территории России. В данной статье будет рассмотрен способ оценки коэффициента фильтрации грунта по данным статического зондирования с использованием численного моделирования. Оценка коэффициента фильтрации грунта была выполнена на основании результатов статического зондирования, полученных на опытной экспериментальной площадке. 1. Описание опытной площадки Экспериментальная площадка находится в пригороде г. Перми. На данной площадке были выполнены инженерно-геологические изыскания с проведением статического зондирования, а также испытания натурных забивных висячих свай [10, 11]. В настоящей статье представлены результаты диссипационных испытаний основного несущего слоя свайных фундаментов. Данный слой характеризуется водонасыщенным суглинком мягкопластичной консистенции (ИГЭ-1). Для рассматриваемого инженерно-геологического элемента было проведено статическое зондирование с измерением порового давления. Статическое зондирование выполнялось установкой GeoMil LWC100 XS с гидравлической системой задавливания зонда. Зондирование осуществлялось электрическим пьезоконусом с датчиком порового давления и муфтой трения (зонд 2-го типа по ГОСТ 19912). На площадке были выполнены диссипационные тесты (тесты на рассеивание поровой воды). Испытания проводились путем остановки пьезозонда на заданной глубине, после чего фиксировался процесс падения давления воды в грунте с течением времени. Эксперимент проводился до 90 % падения порового давления и длился на каждой отметке около 20 мин. Диссипационные тесты проводились в двух точках на трех разных отметках. Их результаты представлены в конце статьи. Время рассеивания порового давления по результатам диссипационных тестов в различные временные интервалы приведено в табл. 1. Таблица 1 Результаты диссипационных тестов по данным численного моделирования Table 1 The results of dissipation tests from numerical modeling Отметка испытания Номер точки статического зондирования По данным численного моделирования По данным натурных испытаний kф, м/сут t50, с t70, с t90, с kф, м/сут (Parez и Fauriel, 1988) t50, с t70, с t90, с 7,5 м CPTU5 0,00030 155 317 701 0,00018 156 345 1200 CPTU8 0,00018 257 528 1169 0,00007 352 690 1620 8,5 м CPTU5 0,00025 154 347 646 0,00011 226 430 1380 CPTU8 0,00015 339 594 1128 0,00012 217 390 1110 9,5 м CPTU8 0,00020 159 379 539 0,00013 204 390 1260 CPTU5 0,00018 221 431 749 0,00011 231 420 1470 Среднее значение 0,00021 214 432 822 0,00012 231 444 1340 Примечание: kф - коэффициент фильтрации, t50 - время при 50 % рассеивании порового давления; t70 - время при 70 % рассеивании порового давления; t90 - время при 90 % рассеивании порового давления. По времени рассеивания порового давления был произведен расчет коэффициента фильтрации по методике Parez и Fauriel [7], которая является аппроксимацией экспериментальных данных, полученных учеными, и заключается в определении коэффициента фильтрации по времени рассеивания порового давления до 50 % от максимального значения. Результаты расчета по данной методике представлены в табл. 1. Сопротивление грунта по муфте трения и конусу зонда на отметках проведения диссипационных тестов приведены в табл. 2. Таблица 2 Удельное сопротивление грунта по конусу и муфте трения по данным численного моделирования Table 2 Sleeve friction resistance and cone resistance based on the results of numerical simulation Отметка испытания Удельное сопротивление по конусу qc, МПа Удельное сопротивление по муфте трения fsi, кПа Среднее значение по данным зондирования По данным численного моделирования Среднее значение по данным зондирования По данным численного моделирования 7,5 м 1,09 0,81 11.2 12,4 8,5 м 1,18 0,90 13,9 13,5 9,5 м 1,01 1,02 12,6 15,6 2. Лотковые эксперименты Для разработки и калибровки численной модели погружения зонда в слабый водонасыщенный глинистый грунт были проведены лотковые испытания [12] в условиях объемной деформации с использованием испытательного стенда (лотка). Данный стенд представляет собой лоток размером 1´1´1 м, заполненный глинистой пастой, электрический зонд, датчик измерения глубины и скорости погружения зонда, блок регистрации сигналов. Лоток заполнялся глинистой пастой с заданной влажностью и послойным уплотнением грунта. Влажность и степень уплотнения грунта были подобраны таким образом, чтобы по физико-механическим свойствам она была приближена к глинистым грунтами, залегаемым на опытных площадках. В точках проведения зондирования с целью оценки деформации массива грунта глинистая паста укладывалась слоями высотой 10-15 мм с устройством маяков. В качестве маяков использовалась специальная водорастворимая бумага разных цветов, которая укладывалась после уплотнения каждого слоя. Контроль однородности укладки и уплотнения глинистой пасты осуществлялся путем отбора грунта из каждого слоя с определением плотности лабораторным методом. После проведения испытания при откопке грунта были отобраны пробы грунта: для испытания на сдвиг, компрессионных и трехосных испытаний. Всего было выполнено три погружения зонда на глубину от 12 до 13 см. Результаты зондирования в лотке представлены на рис. 1. Среднее удельное сопротивление по конусу зонда на отметке вдавливания 12 см составило qc = 0,21 МПа. После погружения зонда производилась откопка зонда с измерением вертикальной деформации окружающего грунтового массива. Измерение деформаций проводилось с точностью до 1 мм. Результаты полученных деформаций грунтового массива после погружения зонда представлены на рис. 2, а. Полученные деформации глинистой пасты имеют тот же характер, что и у других исследователей, в частности проведенных БашНИИстроем с помощью рентгена и представленных на рис. 2, б [13]. Рис. 1. Графики зависимости удельного сопротивление грунта по наконечнику зонда от глубины погружения зонда в лоток Fig. 1. The graph of cone resistance on the depth in the tray а б Рис. 2. Вертикальные деформации грунта после погружения зонда: а - по данным проведенных экспериментов; б - по данным экспериментов БашНИИстроя (НИИпромстрой) [13] Fig. 2. Vertical deformation of the soil after the penetration of the cone: a - according to the data of the conducted experiments; b - according to the data of the experiments of BashNIIstroi (NIIpromstroy) [13] На основании данных результатов была разработана численная модель погружения зонда в грунт. 3. Численное моделирование Численное моделирование производилось с помощью программного комплекса Plaxis 2D c использованием модели упрочняющегося грунта (Hardening Soil) [14-16]. Параметры грунтовой модели для лотковых испытаний определялись в соответствии с рекомендациями Руководства Plaxis [14-16]. Прочностные характеристики грунта принимались по данным лабораторных испытаний. Деформационные характеристики грунта принимались по данным компрессионных и трехосных испытаний глинистой пасты. Входные параметры глинистой пасты, принимаемые при численном моделировании, представлены в табл. 3. Параметры грунтовой модели для опытной площадки определялись в соответствии с рекомендациями Руководства Plaxis [14-16]. Прочностные характеристики грунта принимались по данным лабораторных испытаний. Деформационные характеристики грунта принимались по данным компрессионных испытаний. Входные параметры для суглинка мягкопластичного, принимаемые при численном моделировании, сведены в табл. 3. Таблица 3 Входные параметры грунтов для модели упрочняющегося грунта Table 3 Input parameters of soils for the model of hardening soil Параметр Глинистая паста Суглинок мягкопластичный Объемный вес грунта γ, кН/м3 17,6 19,4 Удельное сцепление сref, кПа 19 14 Угол внутреннего трении φ, град 7 10 Базовое давление pref, кПа 100 100 Показатель степенной зависимости жесткости от уровня напряжений m 1 1 Секущий модуль деформации грунта E50ref, МПа 2,8 2,8 Одометрический модуль деформации МПа 3,3 3,6 Моделирование производилось в двумерном программном комплексе Plaxis 2D в осесимметричной постановке задачи. В данной постановке рассматривается модель в виде сектора вращения в 1 радиан с предположением, что в любом из радиальных направлений деформация и напряжения являются одинаковыми. При создании модели грунтовый массив разбивался на 15-узловые треугольные конечные элементы с интерполяцией четвертого порядка для перемещений и интегрирования по 12 точкам напряжений. С целью уменьшения времени численного расчета применялось локальное измельчение сетки вокруг конуса и боковой поверхности зонда. Расчетная модель была составлена в натурном масштабе 1:1 со стандартными геометрическими размерами зонда радиусом 18 мм и углом при вершине 60°. Граничные условия в расчетной схеме принимались следующими: по нижней границе модели запрещены перемещения по вертикали и горизонтали, по вертикальной границе запрещены горизонтальные перемещения, по верхней границе перемещения не ограничены. Одна из главных задач численного эксперимента состояла в анализе напряженно-деформированного состояние вокруг боковой поверхности зонда, в частности распределения порового давления в месте установки датчика непосредственно между муфтой трения и конуса зонда. При стандартных подходах моделирования контактной поверхности путем продления интерфейсных элементов в углах зонда возникала концентрация напряжений, которая сильно влияла на распределение порового давления на данном участке и сходимости задачи в целом. С целью уменьшения концентрации напряжений вокруг углов модели зонда, а также минимизации влияния интерфейсных элементов на распределение давления в зоне выше конуса была разработана модель с равномерным сглаживанием сетки конечных элементов в углах зонда (рис. 3). На первом этапе численного моделирования были произведены проверка и калибровка полученной численной модели по данным лотковых экспериментов. По результатам численного моделирования удалось добиться качественной сходимости полученных расчетных сопротивлений грунта по конусу зонда с экспериментальными данными. В конце расчета значение удельного сопротивления по конусу зонда составило qc = 0,23 МПа, что близко к экспериментальному значению qc = 0,21 МПа. На втором этапе был произведен расчет на отметках проведения диссипационных тестов на опытной площадке. В первую очередь была произведена калибровка прочностных параметров грунта. Результаты калибровки представлены на рис. 4 и сведены в табл. 2. Рис. 3. Расчетная конечно-элементная сетка вокруг конуса зонда Fig. 3. Calculated finite element mesh around the cone Рис. 4. График зависимости удельных сопротивлений грунта по конусу и муфте трения вертикальных перемещений моделируемого зонда для опытной площадки: а - на глубине 7,5 м; б - на глубине 8,5 м; в - на глубине 9,5 м Fig. 4. The graph of cone resistance and sleeve friction of the vertical displacements of the simulated cone for the experimental site: a - at a depth of 7.5 m; b - at a depth of 8.5 m; c - at a depth of 9.5 m Как видно из табл. 2, в результате проведенной калибровки результаты численного расчета дают величины qc и fsi, близкие к натурным значениям. На следующем этапе была произведена оценка максимального избыточного порового давления, которое возникает после погружения зонда. Полученные максимальные значения порового давления, возникающего после погружения зонда, представлены в табл. 4. Как видно из табл. 4, максимальное поровое давление, полученное по данным численного моделирования, хорошо сходится с реальными значениями, полученными при зондировании на опытной площадке. Таким образом, можно переходить к расчету консолидации грунта и оценке коэффициента фильтрации. Таблица 4 Избыточное поровое давление, возникающее в процессе погружения зонда Table 4 Excess pore pressure after the penetration of the cone Отметка испытания, м Максимальное избыточное поровое давление в месте расположения датчика зонда (между конусом и муфтой трения) u, кПа Среднее значение по данным зондирования По данным численного моделирования 7,5 201 235 8,5 271 256 9,5 243 270 Расчет консолидации грунта производился до рассеивания избыточного порового давления. Численный эксперимент, как и натурный, проводился до 90 % падения порового давления. В результате проведенных расчетов по времени рассеивания порового давления была откалибрована грунтовая модель и получены коэффициенты фильтрации грунта. На рис. 5 представлены графики зависимости порового давления от времени в откалиброванной модели грунта и их сопоставление с реальными диссипационными тестами. Расчетные коэффициенты фильтрации и время рассеивания порового давления в различные интервалы времени приведены в табл. 1. а б в Рис. 5. Графики зависимости порового давления от времени по результатам численного моделирования: а - на отметке 7,5 м; б - на отметке 8,5 м; в - на отметке 9,5 м Fig. 5. Graphs of pore pressure from time as a result of numerical simulation: a - at a depth of 7.5 m; b - at a depth of 8.5 m; c - at a depth of 9.5 m Как видно из табл. 1 и графиков на рис. 5, удалось получить достаточную сходимость решения задачи в пределах времени, при котором поровое давление рассеивается до 70 % от максимального значения. Коэффициент фильтрации, полученный по численным расчетам, имеет тот же порядок, что и по методике Parez и Fauriel. Таким образом, по результатам проведенных расчетов удалось определить коэффициент фильтрации глинистого грунта, который дает хорошую сходимость с проведенными на опытной площадке диссипационными тестами. Выводы 1. Разработана численная модель погружения зонда в глинистый грунт, которая при соответствующей калибровке обладает хорошей сходимостью с реальным статическим зондированием. 2. С помощью численного моделирования можно определить коэффициент фильтрации путем калибровки параметров грунта при проведении консолидационного расчета. 3. Результаты калибровки параметров грунта по численной модели зонда будут использованы в дальнейших численных расчетах вдавливаемых или забивных свай.

Об авторах

А. Б Пономарев

Пермский национальный исследовательский политехнический университет

М. А Безгодов

Пермский национальный исследовательский политехнический университет

Список литературы

Статистика

Просмотры

Аннотация - 222

Ссылки

• Ссылки не определены.