Комплексное моделирование распределения нагрузок на обделку тоннелей метрополитенов Текст научной статьи по специальности «Строительство и архитектура»

ОСВОЕНИЕ ПОДЗЕМНОГО ПРОСТРАНСТВА

УДК 622.831.3:624.19

М.О.ЛЕБЕДЕВ

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

КОМПЛЕКСНОЕ МОДЕЛИРОВАНИЕ РАСПРЕДЕЛЕНИЯ НАГРУЗОК НА ОБДЕЛКУ ТОННЕЛЕЙ МЕТРОПОЛИТЕНОВ

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

Величины нагрузок зависят от расстояния от забоя, на котором обделка вступит в контакт с массивом и будет выполняться условие непрерывности перемещений контура тоннеля и обделки. Задача по определению нагрузок на крепь решается из условия совместности деформирования системы обделка - грунтовый массив. Рассматривается упругая модель взаимодействия обделки с массивом пород по Г.Н.Савину. При взаимодействии обделки с грунтом обделка вступает в работу с массивом не мгновенно, а на некотором расстоянии от забоя. Для приближения расчетных значений нагрузок к фактическим нужно учитывать влияние забоя. Решение задачи взаимодействия сводится к контактной задаче - задаче механики сплошных сред. При этом модуль деформации заменяется приведенным модулем.

Предложенный метод расчета нагрузок на обделки тоннелей программно реализован в среде DELPHI for Windows. Программа позволяет производить графическую двухмерную и объемную интерпретацию расчетных данных и выводить результаты работы программы на принтер.

Tunnel lining load-forming process generally may be represented by the curve with progress of which three components are to be marked out. They characterize the face and tunnel lining contact conditions with massif depending on lining construction technology.

The load values depend on the distance from the face at which lining gets in contact with massif, therefore tunnel and lining contours replacement continuity condition is met. Support loads determination problem is solved only on lining-massif system straining simultaneity principle. The elastic model of lining-massif interaction according to G.N.Savin estimates is considered,

When interacting with soil the lining starts to contact with massif not instantly but at some distance from the face. For the approximation of loads calculated values to real ones the face effect is to be taken into account. The solving of interaction problem comes to a contact problem - the medium mechanics problem. At that the deformation modulus is substituted by reduced modulus.

The offered tunnel lining loads calculation method is realized as DELPHI software for Windows. The program allows to perform two- and three-dimensional graphic interpretations of calculated data and to out-type the program work results.

60 _

ISSN 0135-3500. Записки Горного института. ТЛ50. Часть 2

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

Аналитическая и статистическая обработка данных наблюдений за показаниями датчиков напряжений и деформаций, полученных сотрудниками лаборатории ЦНИИС и непосредственно автором вместе с институтом «Ленметрогипротранс», показала, что процесс формирования нагрузки на обделку перегонных тоннелей в общем виде можно представить кривой на рис Л. Анализ процесса развития нагрузки на обделку позволяет выделить три составляющих: первая (линейный участок до излома) - изменение нагрузки при влиянии забоя; вторая (линейный участок после излома) - изменение нагрузки в условиях перераспределения напряжений в массиве, связанное с контактными условиями обделки и массива; третья (линейный участок после второго излома) -изменение нагрузки, связанное с реологическими характеристиками массива.

Величины нагрузок зависят от расстояния от забоя, на котором обделка вступит в контакт с массивом и будет выполняться условие непрерывности перемещений контура тоннеля и обделки (рис.2). Для разжатой на грунтовый массив обделки это расстояние соответствует месту установки, т.е. £ = 6,8 м. Для сборной обделки это расстоя-

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

Время

Рис. 1, Схема развития нагрузки на обделку

-•:.....

- У'. Н* . -

о

'V

Р.

Рис.2, Принципиальная схема механизма формирования нагрузки на обделку и смещения

контура тоннеля

делки, применяемых в Санкт-Петербурге, это расстояние составляет 7,5-9,0 м.

Задачу по определению нагрузки на крепь решим из условия совместности деформирования системы крепь - порода

Впервые условие совместности перемещений для расчета нагрузок на крепь стволов было предложено Ф.А.Балаенко в виде

оо Цп

ик (Р),

(1)

где Ц* - радиальное перемещение контура выработки по истечении достаточно большого промежутка времени от момента об-нажения; 17о - радиальное перемещение контура выработки до ввода крепи в с массивом пород; ик(Р) - радиальное перемещение контура крепи под нагрузкой Р.

---—_—-- 61

Санкт-Петербург. 2002

Методы расчета нагрузок по схеме взаимодействия развивались К.В.Руппеней-том, Н.С.Булычевым, А.Г.Протосеней и другими исследователями.

При малых величинах нагрузки глина деформируется в рамках упругой модели. Рассмотрим упругую модель взаимодействия крепи с массивом пород по Г.Н.Савину [4]. Крепь и массив пород считаются выполненными из упругого изотропного материала с различными деформационными характеристиками. Крепь вводится в работу с массивом без зазора, т.е. мгновенно после обнажения массива. Нагрузка на крепь:

р = о = 0,5(1 +

(1 + Х)(с2-1)

А(2-с2) + ВС

+

+ 0,5(Я, -1 )уН

1 + Х

D

х

2 -Зс4 -1) + (2)

1

4 Вс\сА -4с2 4 3).|cos20.

где

D, = с

2

X+

Е

N

Е

А13с4 - + 41+ Вс

+

к J

+

\[А + Вс6)-,

Е

(3)

К

А

Е

1; В = 1 + Xi

Е

К

Е

К

% = 3-4У; XI =3-4Ук;

X - коэффициент бокового распора; у -удельный вес пород; Н - глубина заложения тоннеля; с = Я\/Яо] Я\ ~ внешний радиус

крепи; 7?о - внутренний радиус крепи; 0 -угол между осью х и рассматриваемым сечением обделки; Е и Ек - соответственно модуль упругости пород и крепи; V и ук -соответственно коэффициент Пуассона пород и крепи.

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

При взаимодействии обделки с грунтом крепь вступает в работу с массивом не мгновенно, а на некотором расстоянии от 62 _

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

В работе [3] показано, что решение задачи взаимодействия (1) сводится к контактной задаче - задаче механики сплошных сред. При этом модуль упругости нужно заменить приведенным модулем^ определяемым по формуле

£* =exp(-aj// R)EKX

(5)

где В. - радиус выработки; он - константа, определяемая на основе натурных наблюдений; / - расстояние от места возведения крепи до забоя.

Анализ натурных данных взаимодействия различных типов обделок с массивом дает основание преобразовать величину / (рис.2):

(6)

где /д - расстояние от забоя до места возведения обделки; 1К - дополнительное расстояние от места возведения обделки, на

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

Для нахождения он используем аппроксимацию [3]:

и0 = и[ж)/{1\

где /(/) - функция, учитывающая влияние

забоя на развитие перемещений в призабой-ной зоне выработки.

Функцию /(/) определим из экспоненциальной зависимости [1,2]:

/(/)=l-exp(-ä1/1/*X

(8)

где /1 - расстояние, на котором влияние забоя на смещение контура в сторону выработки прекращается (на основании натурных наблюдений).

Для определения константы он внесем функцию (8) в (7):

(9)

Показатель он находится из условия, что при / = 4^1, С/о= 0,01С/(ао) •

ISSN 0135-3500. Записки Горного института. ТЛ50. Часть 2

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

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

(10)

где 8 - характеристики ползучести, полученные экспериментально.

Для решения задач теории ползучести вместо упругих постоянных Работновым предложено использовать временные функции:

г Е 0,5 - V Е< =--; у = 0,5--

1 + Ф

где Ф - функция ползучести,

1 + Ф

(П)

Ф

81

1-а

1 -а

Внося соотношения (11) в формулы по определению нагрузки на обделку тоннелей (2) - (4), получим окончательные формулы по определению нагрузок при вязкоупругой модели.

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

Объемный вид К*1 В*3

1 ? в -108 46 1 Г

| Показать оси I Тоже ль

: 20 проекция ! Показ проекции

183.08 60.76

.21

155.36 5211 148.31

43.70

32.13

по гориз.

и _} ±1

Поворот по верт.

« I I ►

Зарыть

1.;___

Рис.3. Объемный вид решаемой задачи

Сравнительные результаты величин расчетных и фактических (поданным натурных исследований)

нагрузок на обделку перегонных тоннелей

Тип обделки Я, м /?св, м м Ек, МП а ЕГр, МПа % 1, м сут, , кПа Яр, кПа

Обжатая 60 2,665 2,815 30000 762 0,15 0,25 5,6 750 490 450

Сборная 40 2,5 2,75 35000 762 0,15 0,25 8,0 420 220 232

_ 63

Санкт-Петербург. 2002

массива и обделки, а также от технологии строительства.

Программный модуль «Геомеханика» выполнен в среде DELPHI for Windows. Объем рабочей программы составляет 267 кБ. Программа написана с использованием стандартных объектов библиотеки и функций DELPHI, имеет главное и несколько вспомогательных окон. Во вспомогательных окнах реализованы диалоги и панели управления для удобства работы с программой, которые управляются главным меню. Программа позволяет производить графическую двухмерную и объемную интерпретацию расчетных данных и выводить результаты работы программы на принтер, работать с объемным видом формирования напряжений на контакте обделки с массивом (рис.3), основными зависимостями напряжений от влияющих факторов.

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

Рассмотрим численное решение определения нагрузок на обделку по предлагаемой методике для двух технологических схем проходки. Первая схема - щитовая проходка с обжатой на породу обделкой. Вторая - эректорная проходка со сборной железобетонной обделкой.

Исходные данные и средние нагрузки, полученные по данным натурных исследований <7ф и расчетным путем др сведены в

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

чем на 9-10 %.

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

ЛИТЕРАТУРА

1. Протосеня А.Г Расчет нагрузок на крепь вертикальных стволов при больших глубинах / А.Г.Протосеня, Б.В.Бокий. Ю.С.Обручев V Шахтное строительство. 1974. № 1С.4-6.

2. Протосеня А.Г. Расчет нагрузок на крепь глубоких стволов, сооружаемых в сложных горно-геологических условиях / А.Г.Протосеня, А.М.Козел, В.А.Борисовец // Шахтное строительство. 1984. №6. СЛЗ-15.

Протосеня А.Г. О постановке задач по расчету нагрузок на крепь капитальных выработок и тоннелей. // Устойчивость и крепление горных выработок. Крепле-

■'V . '■;<•

ние и поддержание горных выработок / Санкт-Петербургский горный ин-т. СПб, 1992. С.4-8.

4. Савин Г.И. Концентрация напряжений около отверстий / ГИТТЛ. М., Л., 1951. 494 с.

Научный руководитель профессор, д.т.н. А.Г.Протосеня

64 _-

ISSN 0135-3500. Записки Горного института. ТЛ50. Часть 2