CC BY
136
СПИСОК ЛИТЕРАТУРЫ
1. Месчян С.Р. Начальная и длительная прочность глинистых грунтов. - М.: Недра, 1978. - 207 с.
2. Bishop A.W., Green G. The influence of end restrain on the compression strength of a cohesionless soil // Geotechnique. - 1965. -V. 15. - № 3. - P. 243-266.
3. Dounias G.T., Potts D.M. Numerical analyses of drained direct and simple shear tests // Journal of Geotechnical engineering, ASCE. -1993. - V. 119. - № 12. - P. 1870-1891.
Frydman S., Talesnick M. Development of strain during monotonic shear of soft clay // Journal of Geotechnical engineering, ASCE. -1992. - V. 118. - № 5. - P. 704-725.
Lade P. Torsion shear apparatus for soil testing // Laboratory shear strength of soil, ASTM, STP 740. - 1981. - P. 145-163. DIN 18137-3. September 2002. Bestimmung der Scherfestigkeit. -2002. - S. 33.
Поступила 25.06.2007г.
УДК 624.131
МОДЕЛИРОВАНИЕ ОСЕДАНИЯ ПОВЕРХНОСТИ ПРИ ПРОХОДКЕ ТУННЕЛЯ ЩИТОВЫМ СПОСОБОМ
Л.А. Строкова
Технический университет, г. Мюнхен, Германия E-mail: [email protected]
Изложены результаты моделирования оседания поверхности в результате проведения горных работ щитовым способом для строительства линии метрополитена U-8-Nord г. Мюнхена. Для верификации расчетной модели служили замеры оседания поверхности по окончании горных работ. Выбрана модель грунта, наиболее достоверно описывающая оседание земной поверхности при проходке. Определены факторы, контролирующие величину осадки поверхности. Деформации грунта, связанные с осушением массива при выполнении горных работ не учитывались.
Одним из характерных видов деформации грунтов при проходке горных выработок является оседание дневной поверхности (рис. 1). С целью минимизации повреждений существующих зданий и фундаментов на поверхности земли, необходимо прогнозировать эти эффекты и принять соответствующие меры. Р. Пек в 1969 г. предложил для математического описания оседания поверхности [1] использовать формулу Гаусса для нормального распределения
Sv(x) — Sv,
поперечное направление
■+у продольное направление
Рис. 1. Геометрия оседания поверхности, вызванной проходкой туннеля [2]
Типичная поверхность оседания в поперечном направлении показана на рис. 2. Область выше кривой оседания задается уравнением
V — J Sv • dx — yjln • i • Sv]
где V - объем осевшего грунта над туннелем, в результате проведения горных работ, теоретически равен величине сжатия сечения туннеля V, (рис. 2) на единицу длины проходки.
где - максимальная осадка поверхности над осью туннеля, 4 - расстояние до точки перегиба кривой оседания поверхности, соответствующее стандартному отклонению в уравнении Гаусса.
мульда оседания
Рис. 2. Основные характеристики, используемые в расчете [4]
Для оценки оседания поверхности использовался метод конечных элементов, реализованный в программе РЬАХВ которая позволяет моделировать процессы экскавации грунта, возведения обделки туннеля [3]. Расчеты выполнялись для строящейся линии метрополитена г. Мюнхена. Для верификации модели служили данные натурных наблюдений оседания поверхности на уже законченном участке метро, выполненные сотрудниками технического университета под руководством И. Филлибека. Максимальная величина оседания
о ■
Рис. 3. Схематический разрез северной части г. Мюнхена
поверхности составила 7 мм. Оседания зафиксированы на расстоянии до 10 м от оси туннеля.
Геологический разрез г. Мюнхена представлен четвертичными гравийно-щебенистыми грунтами, подстилаемыми третичными песчаными и глинистыми грунтами. Продольный разрез по оси туннеля представлен на рис. 3.
Таблица. Основные характеристики грунтов
Название Обозначение в «PLAXIS» Гравий Песок Глина Ед. изм.
Модель поведения HS: «Hardening soil» Упруго-пластический материал с упрочнением -
Тип поведения Drained Дренированный -
Удельный вес грунта Yunsat 23 21 21 кН/м3
Удельный вес водо-насыщенного грунта Ysat 24 22 21 кН/м3
Модуль Юнга при первичном нагру-жении E50=Eoed 120 95 60 МН/м2
Модуль Юнга при разгружении - повторном нагружении Eur 240 190 120 М/м2
Показатель степени m 0,5 0,4 0,4 -
Коэффициент Пуассона V 0,3 0,3 0,3 -
Сцепление се 1 1.0 25 кН/м2
Угол трения V 37,5 35 25 град.
Угол дилатансии T 0 0 0 град.
Давление Pref 100 100 100 -
Коэффициент бокового давления грунта при нормальной консолидации K0NC=1-sinv 0,391 0,426 0,577 -
Расчетное сечение, по которому мною проводилось моделирование, представлено тремя различ-
J Туннель
А-А
Расчетное сечение отложения
ными по литологии слоями. Верхний слой мощностью 10 м состоит из гравия. Под слоем гравия залегает слой песка мощностью 4 м. Ниже песчаного слоя залегает слой глины мощностью 8,0 м. Показатели физико-механических свойств грунтов приведены в таблице. Модули упругости и прочностные показатели грунтов определялись по данным трехосных и компрессионных испытаний. Модуль сдвига и деформация сдвига назначены по литературным данным. Уровень грунтовых вод расположен на 8 м ниже поверхности земли. Распределение порового давления гидростатическое.
Проходка туннеля осуществлялась щитовым способом в слоях гравия и песка. Туннель имеет диаметр 7,0 м и расположен на глубине 13 м. Основные свойства туннельной обделки следующие: нормальная жесткость 6,0.106 кН/м, изгибная жесткость 2,0.104 кНм2/м, эквивалентная толщина 0,2 м, вес 8,4 кН/м/м, коэффициент Пуассона 0,33. Геометрическая модель для выбранного расчетного сечения представлена на рис. 4.
В качестве основного типа конечных элементов использовался 15-узловый элемент. В области оседания поверхности сетка конечных элементов (рис. 5) более мелкая.
Моделирование процессов экскавации грунта и возведения обделки выполнялось по двум методикам [3]: в-метод (с в: 0,4; 0,5; 0,6 и 0,7 %) и <volume loss> (с величиной сжатия сечения туннеля Vt: 0,2; 0,4 и 0,6 %). После выполнения расчетов в программе PLAXIS, данные по вертикальным перемещениям на верхней границе расчетной модели переносились в Excel, где проводился расчет площади элементарных трапеций поверхности оседания,
Вейпретитр. АмХарт Фрунквуртер Ринг
отложения е ершей терр асы §
" UD
суглинки g- СТЗ граБиино-гапечникоБые
Рис. 4. Геометрическая модель участка
Рис. 5. Сетка конечных элементов
определялась суммарная площадь выше кривои поверхности оседания Ах и, наконец, определялся
объем потери грунта по формуле V для со-
поставления данных объема осевшего грунта V и величины сжатия сечения туннеля V. Результаты моделирования оседания поверхности при экскавации грунта представлены на рис. 6. Наиболее близко к деиствительнои кривои оседания по максимальной величине осадки отвечает расчет экскавации при сжатии сечения туннеля 0,2 %.
Для определения факторов, определяющих величину оседания поверхности при проходке туннеля, были выбраны следующие параметры: тип грунта, отношение мощности пород над туннелем Н и его
Расстояние до оси туннеля, м -30 -20
0,000
-0,005
-10
диаметра D, величина бокового давления в условиях нормальной консолидации Кхс и в переуплотненном состоянии К0ос, тип расчетной модели грунта на основе теории пластически упрочняющейся среды.
Для оценки влияния типа пород на величину осадки проводилось моделирование двух видов: имитировалась проходка туннеля только в гравийных грунтах и только в глинистых, причем для каждого вида, были выполнены расчеты с различным отношением Н/D, равным 1, 2 и 3. Как видно из рис. 7, основной вклад в величину осадки вносит отношение Н/D. При увеличении глубины залегания туннеля резко увеличивается величина сжатия сечения туннеля, что находит отражение в объеме осевшего грунта, сам вид грунта имеет меньшее влияние.
10
20
30
\v\ TU - ti
ч -0,010
ш и О
-0,015 -0,020
Рис. 6. Кривые оседания поверхности по данным натурных наблюдений и расчетов в PLAXIS
VL 0,6 VL 0,4 VL 0,2 Beta 0,7 Beta 0,6 Beta 0,5
0
HS-HSS. Метод <volume loss>. Разрез: гравий, глина. Грунты нормально-консолидированные (K0NC=1-sin<p )
глина Нй=1 К0-№С гравий H:D=1 К0-№С глина H:D=2 К0-№С гравий Н^=2 К0-№С глина H:D=3 К0-№С гравий H:D=3 К0-№С
НЭЭ глина Н^=1 К0-№С НЭЭ гравий Н:С1=1 К0-№С НЭЭ глина H:D=2 К0-№С НЭЭ гравий H:D=2 К0-№С НЭЭ глина H:D=3 К0-№С НЭЭ гравий Н№3 К0-№С
0,8 1,0 Vt [%]
Рис. 7. Сравнение объемов осевшего грунта над туннелем V и величиной сжатия сечения туннеля V при экскавации для нормально-консолидированных грунтов
Расстояние от оси туннеля [м]
уо,7=0,001
-20 -15 -10 -5 0
-0,004
£ -0,006
-0,008
—замеренная -HS VL 0,28 HS VL 0,26 HS VL 0,24 HS VL 0,22 HS VL 0,20 HS VL 0,18 HSS VL 0,28 HSS VL 0,26 HSS VL 0,24 HSS VL 0,22 HSS VL 0,20 HSS VL 0,18 HSS VL 0,28 HSS VL 0,26 HSS VL 0,24
■ HSS VL 0,22
■ HSS VL 0,20 HSS VL 0,18
y 0,001 y 0,001 y 0,001 y 0,001 y 0,001 y 0,001 y 0,0001 y 0,0001 y 0,0001 y 0,0001 y 0,0001 y 0,0001
-0,012
j Jr\ ff
HSS(y0i7=Qcroj= H S HSS(y0i7=Qooi >HS и HSS(yoi7=o,oooi)
Рис. 8. Кривые оседания поверхности: по данным натурных наблюдений и расчетные с использованием моделей НБ и НББ при варьировании деформации сдвига у и величиной сжатия сечения туннеля VL
0
0,002
Для выбора расчетной модели грунта, наиболее достоверно описывающей оседание земной поверхности при проходке туннеля, были выполнены расчеты с двумя моделями грунта: а) HS (hardening soil) - упруго-пластическая модель с изотропным упрочнением, б) HSS (hardening soil with small - strain stiffness) - упруго-пластическая модель с изотро-
пным упрочнением с учетом жесткости малых деформаций. Обе модели учитывают несовпадение модуля упругости ветвей разгружения и повторного нагружения, наблюдаемое при лабораторных испытаниях [3]. Для второй модели параметры грунта модуль сдвига G0 и деформация сдвига у, назначались по литературным источникам [4, 5], поскольку
Расстояние от оси туннеля [м]
-20 -15 -10
0
-0,001
-0,002
-0,003
| -0,004 я и О
-0,005 -0,006 -0,007 -0,008
Рис. 9. Моделирование оседания поверхности при разных коэффициентах бокового давления грунта К0
замеренная Ив VI. 0,2 КО-ЫО Ивв КО-ЫО -&2Еиг, у 0,001 Ив-М. 0,2 К0-00 Ивв VI. 0,2 К0-00
— 0,2
ел >
2 0,15 о
3
0,1
/
_ Л _ ) /
4 О ^ | / < 1/
О ЫО-Ив О ЫО-Ивв О ОО-Ив ■ ОО-Ивв -Мв=М
0,1 0,15 0,2
Величина сжатия сечения тоннеля, VI [%]
Рис. 10. Сравнение объемов осевшего грунта над туннелем и величиной сжатия сечения туннеля Vt при экскавации для нормально-консолидированного (NC) и переуплотненного состояния (OC) массива. Модели поведения материала под нагрузкой: HS и HSS
их определение требует сложного дорогостоящего оборудования, не входящего в стандартные определения показателей механических свойств грунтов. Для корректировки назначенных параметров, проведены расчеты с варьированием модуля сдвига G0: 1,5-Енг; 2-Еиг; 3-Еиг и деформации сдвига у01: 0,01;
0,001; 0,0001. Результаты моделирования представлены на рис. 8. Наилучшее совпадение расчетных кривых оседания с данными натурных наблюдений наблюдается для расчетов с моделью ШБ при значениях у01= 0,001 и при G0=2•E„ которые рекомендуется использовать при дальнейших расчетах.
0,3
0,25
0,05
0
0
0,05
0,25
0,3
Для оценки влияния степени переуплотнения грунтов в недавнем геологическом прошлом при наступлении оледенений были выполнены расчеты с коэффициентом бокового давления К0ОС для переуплотненного состояния массива. По данным лабораторных исследований, специально выполненных инженером Г. Пельцом, в расчетах были приняты следующие значения К0 для: гравия 0,361; песка 0,6; глины 0,8. Результаты моделирования оседания поверхности с учетом К0 отражены на рис. 9. Наилучший результат показали расчеты с учетом переуплотненного состояния массива К0ОС (модель HSS).
Адекватность модели реальному объекту контролировалась совпадением расчетных мульд оседания и мульды оседания, полученной по натурным наблюдениям. Основным параметром сравнения являлся параметр Vs - объем осевшего грунта над туннелем, получаемый в зависимости, от задаваемой в расчете, величины сжатия сечения туннеля Vt при экскавации, рис. 10. Близость точек Vs к теоретической линии V^Vt подтверждала, что модель определена правильно.
Таким образом, результаты моделирования показали, что основными факторами, определяющими величину осадки, являются: глубина залегания туннеля, тип грунта, коэффициент бокового давле-
СПИСОК ЛИТЕРАТУРЫ
1. Peck R.B. Deep excavations and tunneling in soft ground // Proc. of the 7th Intern. Conf. on Soil Mechanics and Foundation Engineering. State of the art volume. Sociedad Mexicnan de Mecanica de Suelos, A. C., 1969. - P. 225-290.
2. Franzius J.N. Behaviour of buildings due to tunnel induced subsidence: Ph.D. thesis. - Imperial College, University of London, 2003. - 360 p.
3. Plaxis V8 (Finite element code for soil and rock analyses). http://www.plaxis.nl
ния грунта. Выполненные исследования позволяют выполнять прогноз оседания поверхности в ходе дальнейшего строительства метрополитена с учетом установленных в этом исследовании положений.
1. Недопустимо изначально считать, что грунты находятся в условиях нормальной консолидации, без проверки этого положения специальными лабораторными исследованиями.
2. Различие между способом определения величины осевшего грунта по ^-методу или <volume loss> незначительное, хотя первый способ более трудоемок.
3. Обе расчетные модели грунта на основе теории пластически упрочняющейся среды хорошо зарекомендовали себя, хотя предпочтение следует отдавать модели В расчетах параметры модуля сдвига и деформации сдвига, до проведения их экспериментального определения, можно ориентировочно принимать равными 70,7=0,001 и при G0=2•Euг
В настоящее время незначительное оседание поверхности, вызванное проходкой туннеля, не представляет угрозы для существующих зданий. Продолжение мониторинга за оседанием поверхности, вероятно, потребует дополнительного учета реологических параметров грунтового массива.
4. Moller S.C. Tunnel induced settlements and forces in linings // Mitteilungen des Instituts fur Geotechnik Universität Stuttgart. - 2006. - № 54. - 174 p.
5. Benz T. Small-Strain Stiffness of Soils and its Numerical Consequences // Mitteilungen des Instituts fur Geotechnik Universitat Stuttgart. - 2007. - № 55. - 209 p.
Поступила 09.10.2007г.
