л
Научная статья УДК 624.19
Метод прогноза напряженного состояния обделки подземных сооружений квазипрямоугольной и арочной форм
М.А.КАРАСЕВ1, ТИЕН ТАЙ НГУЕН12 Н
1 Санкт-Петербургский горный университет, Санкт-Петербург, Россия 2Ханойский университет горного дела и геологии, Ханой, Вьетнам
Как цитировать эту статью: Карасев М.А., Тиен Тай Нгуен. Метод прогноза напряженного состояния обделки подземных сооружений квазипрямоугольной и арочной форм // Записки Горного института. 2022. Т. 257. С. 807-821. Б01: 10.31897/РМ1.2022.17
Аннотация. Рассмотрен метод прогноза напряженно-деформированного состояния обделок подземных сооружений, форма поперечного сечения которых отлична от кругового очертания. Основной задачей исследования является разработка методики оценки влияния параметров формы поперечного сечения подземных сооружений на напряженное состояние обделки. Для решения поставленной задачи обоснован и развит метод расчета напряженного состояния обделки для тоннелей арочной с обратным сводом и квазипрямоугольной форм. Выполнена апробация методики, которая показала, что точность прогноза напряженного состояния обделки достаточна для выполнения практических расчетов. Предложен алгоритм проведения многовариантного анализа влияния формы поперечного сечения подземных сооружений арочной и квазипрямоугольной форм на напряженное состояние обделки. С использованием разработанного алгоритма выполнены параметрические расчеты и получены закономерности формирования напряженного состояния обделок подземных сооружений для различных инженерно-геологических условий, а также начального поля напряженного состояния. Выполнена количественная оценка влияния геометрических параметров тоннелей на их напряженно-деформированное состояние.
Ключевые слова: тоннель; обделка; напряженное состояние; взаимодействие; нагрузка; метод оптимизации формы; поперечное сечение; грунтовый массив
Поступила: 16.11.2021 Принята: 06.04.2022 Онлайн: 12.10.2022 Опубликована: 10.11.2022
Введение. Разработки в области проектирования тоннелепроходческих комплексов позволяют использовать их для механизированного строительства тоннелей не только кругового, но и арочного, овального и квазипрямоугольного сечений. Выбор формы поперечного сечения определяется рядом факторов, среди которых наиболее значимые: назначение подземного сооружения; геотехнические условия строительства; принятая технология строительства; расчетные экономические затраты. Круговое поперечное сечение является наиболее распространенным при строительстве тоннелей в слабых грунтах, что определяется доступностью и надежностью оборудования для строительства тоннелей таких форм, а также тем, что для большинства горно-геологических условий данная форма наиболее устойчива [1, 2]. Недостатком круговой формы является низкий коэффициент использования сечения - он уменьшается с увеличением площади поперечного сечения. Тоннели некруговых очертаний имеют более высокий коэффициент использования сечения по сравнению с круговым. Однако такое изменение формы ведет к повышению напряжений и снижению устойчивости обделки тоннеля. Снизить напряжения в обделке можно за счет оптимизации формы, в частности, изменения радиуса стенок и свода. С учетом того, что рекомендации для оптимизации геометрии некруговых сечений в общем виде отсутствуют, требуется разработка расчетного решения, которое позволит повысить эффективность использования некруговых форм поперечного сечения подземных сооружений.
Жж © М.А.Карасев, Тиен Тай Нгуен, 2022 ивдЗ™
Исследования в области расчета напряженного состояния обделок подземных сооружений изложены в научных работах, направленных на изучение напряженно-деформированного состояния обделки на основании аналитических решений [3-6], численного моделирования [7-9] и экспериментальных исследований [10-13]. Проведены масштабные экспериментальные исследования напряженного состояния сборных обделок для тоннелей кругового [14-16], прямоугольного [17], квазипрямоугольного [18-20], подковообразного [21] и эллиптического сечений [22]. С учетом сложности проведения масштабных экспериментальных исследований результаты таких работ не носят обобщающий характер, а нацелены на решение задач в конкретных горно-геологических условиях.
Исследователями получены детальные результаты расчета напряженного состояния обделок подземных сооружений, выполненные на основании аналитических решений и численного моделирования [23, 24], предложен подход к параметрическому анализу изменения напряженного состояния обделки тоннелей прямоугольной и эллиптической [25], круговой и прямоугольной форм сечения [26]. В работе [27] оценивается влияние формы поперечного сечения тоннелей на величину оседания земной поверхности. В источнике [28] изучено влияние сдвоенных тоннелей различной формы поперечного сечения на формирование напряженного состояния обделки. В работе [29] выполнен ряд численных расчетов для оценки влияния формы поперечного сечения на напряженное состояние обделки. Результаты выполненных исследований показывают, что напряженное состояние обделки значительно зависит от формы поперечного сечения. Однако вопросы обоснования геометрических параметров обделок некруговых форм в научных работах не рассмотрены.
Основное внимание уделено оценке влияния различных факторов (геологических и геометрических) на формирование напряженного состояния обделок тоннелей некругового очертания. Метод расчета напряженного состояния обделки, основанный на учете породного массива через гиперстатические реакции (Hyperstatic Reaction Method, HRM), более известный как метод Мет-рогипротранса, доказал свою эффективность при расчете обделок тоннелей кругового очертания [13, 30]. Рассмотрено его применение для расчета обделок тоннелей некруговых форм поперечного сечения: арочной с обратным сводом и квазипрямоугольной. Основное преимущество метода HRM - скорость выполнения расчетов - использовано для решения задачи многовариантного анализа влияния формы поперечного сечения тоннеля на напряженное состояние обделки. Разработан алгоритм оценки влияния формы поперечного сечения на напряженное состояние обделки для тоннелей арочной формы с обратным сводом и квазипрямоугольного сечения.
Разработанный метод позволил выявить новые закономерности формирования напряженного состояния обделки тоннелей различной формы поперечного сечения.
Методология. Метод HRM для расчета напряженного состояния обделки подземного сооружения. Метод HRM относится к классу численных методов анализа, в которых напряженное состояние обделки, представляемой в виде одномерных стержневых элементов, может быть определено через интегральные показатели. Взаимодействие между обделкой подземного сооружения и вмещающим породным массивом формируется через набор нормальных и касательных пружин, которые соединяются со стержневыми элементами через смежные узловые точки. Механическое поведение породного массива на допредельной стадии деформирования принимается на основании линейного закона деформирования, в то время как достижение предельного состояния характеризуется идеально-пластическим течением. Узловые перемещения в методе HRM неизвестны, их определение позволяет установить величины интегральных показателей напряженного состояния стержневых элементов, т.е. определить напряженное состояние обделки. Внешнее давление в виде сосредоточенных сил прикладывается к узловым точкам. Детальное описание использования метода HRM для расчета напряженного состояния обделки кругового очертания представлено в [13, 30]. В настоящей работе метод HRM расширен для прогноза напряженного состояния некруговых обделок (арочная форма с обратным сводом, квазипрямоугольная форма).
Взаимосвязь между глобальной матрицей жесткости рассматриваемой системы K, векторами неизвестных перемещений U и узловых сил F можно представить как [31 ]
KU = F.
(1)
Раскрыв уравнение (1), получим:
Ka+ka кс,ь 0 0
Kc ke,d+k2,a k2,b 0
0 k2,c k2,d+k3,a k3,b
0 0 k3,c k3,d+k4
0 0 0 0
k-,b 0 0 0
4, a
0 0 0
-_1,c
0 о о
К-ьь
k--1,d + k-,a
" Fe '
U 2
u3 =
P*.
(2)
где кг,а, кг,ь, к.1,о и кг^ - саб-матрицы матрицы к; иг, Сэ, ... Ц - компоненты вектора перемещений; р1, Рг, Fэ, ... Рп - компоненты вектора сил.
Наличие пружин по периметру тоннеля влияет на жесткость рассматриваемой системы стержневых элементов. Это учитывается за счет изменения жесткости каждого из стержневых элементов, следовательно, и глобальной матрицы жесткости.
Методика расчета нагрузок на обделку тоннелей различного очертания. Существуют различные гипотезы формирования горного давления на обделку подземных сооружений, каждая из которых имеет свою область применения, исходя из инженерно-геологических условий и в зависимости от положения подземного сооружения относительно земной поверхности. Применимость методов расчета нагрузок на обделку исследована многими учеными и широко представлена в литературе [32-35], поэтому не является основным объектом данного исследования. Принята методика расчета нагрузки на крепь, предложенная Терцаги [36], который рассматривал формирование опускающегося столба породы шириной В1, часть нагрузки от которого передается на окружающий породный массив за счет трения, а другая составляющая удерживается обделкой подземного сооружения, что позволяет сохранить равновесие рассматриваемой системы взаимодействия «подземное сооружение -грунтовый массив».
Согласно методике, предложенной Терцаги, при глубине заложения подземного сооружения менее двух его ширин величина вертикального давления на обделку подземного сооружения определяется полным весом столба породы:
öv = yz,
(3)
где у - усредненный объемный вес вышележащих пород; г - расстояние от поверхности земли до расчетной отметки.
Для расчетных случаев, когда величина глубины заложения подземного сооружения составляет более двух его ширин, расчет нагрузки на обделку определяется через величину эффективной мощности пород Ио, вес пород в которой определяет величину давления:
öv = yho.
(4)
При этом принимается, что опускающийся столб породы имеет ширину В1, а в боках выработки формируются призмы сползания, угол наклона которых принимается равным (л/4 + ф/2). Величина эффективной мощности опускающегося столба породы определяется по формуле
h = Вс (С (С / ВсУ)) Л _ e-K,tan(p)(H/Be )) + Ро_ (e-Kotan(<p)(H/В, )) (5)
0 K tan(m) ( ) у ( ),
где В1 - половина ширины опускающегося столба породы (табл.1); с - сцепление; Ко - коэффициент бокового давления; Н - глубина заложения подземного сооружения; Ро - дополнительная нагрузка, приложенная на поверхности земли.
Формулы для определения ширины В1 подземных сооружений различной формы поперечного сечения сведены в табл. 1.
Таблица 1
Определение ширины подземного сооружения
Особенности подготовки расчетной модели в рамках метода HRM. При выполнении расчетов напряженного состояния обделки тоннелей приняты схемы передачи нагрузки, обозначенные на рис.1, а-в. Взаимодействие между обделкой подземного сооружения, описанной стержневыми элементами, и породным массивом осуществляется через пружины, ориентированные в нормальном и касательном направлениях относительно контура подземного сооружения (рис.1, г). Жесткость нормальных и касательных пружин определяется как к« и к соответственно. Значение величин к« и к определяется через величины п« и п«, характеризующие способность вмещающего породного массива сопротивляться деформациям.
Нелинейная зависимость, определяющая связь между деформацией обделки 5 и реакцией породного массивар на эти деформации, представлена в [13, 30]:
Oh
ршшщщс
ks
'/яз ' Rl-^t,
иШШШШШ
Oh
Узел i + 2
Узел I
Рис.1. Расчетные схемы прогноза напряженно-деформированного состояния обделки подземных сооружений различных форм поперечного сечения: а - круговая форма; б - арочная форма с обратным сводом; в - квазипрямоугольная форма; г - схема организации взаимодействия между обделкой подземного сооружения и породным массивом
ау, а, - вертикальная и горизонтальная нагрузки соответственно; ЕА, Е1 - осевая и изгибная жесткости обделки соответственно; Е - модуль деформации материала обделки; А - площадь поперечного сечения обделки; .1 - момент инерции обделки; Я - радиус подземного сооружения (круговая форма поперечного сечения); Я1 - радиус боковой стенки; Я3 - радиус свода подземного сооружения; Я2 - радиус обратного свода (арочная форма поперечного сечения с обратным сводом) или опорной части свода (квазипрямоугольная форма)
б
a
в
о
о
v
v
O
v
v
г
о
x
Р = Рш
1 —
Plim Plim + По5.
(6)
где рнш - предельная величина реакции породного массива; по - начальная величина жесткости породного массива при значении 5, близком к нулю.
Показатель относительной жесткости породного массива можно выразить следующим образом:
П =
Рlim
5
1 --
Рlim
Л
Plim + По5
(7)
Жесткость пружин в нормальном направлении изменяется в зависимости от радиуса подземного сооружения и определяется по формуле
Пп ,0 = 2
1 E
1 + V, R
(8)
в то время, как жесткость пружин в касательном направлении может быть принята равной одной трети от жесткости пружин в нормальном направлении [13]:
_
Л*,0 = ,0(3)'
(9)
где Es - модуль деформации породного массива; Vs - коэффициент поперечной деформации породного массива; Ri - радиус кривизны обделки в каждой рассматриваемой точке ее периметра.
Предельные величины реакции породного массива в радиальном рп,ш и в касательном ps,lim направлениях зависят от прочностных характеристик пород (сцепления и угла внутреннего трения) величины, характеризующей всестороннее обжатие пород Aoоonf:
и,Нш
2c 008 ф 1 + sinф
1 - 8Ш ф 1 - sinф
До
сон/'
(10)
Да , = а + а
сон/
2 1 - V
Жесткость пружин кп,1 и может быть определена по следующим формулам:
к, = лИ,
к,, =
(4-1+ Ц) 1
2
(1Ц) 1
2
?,Иш
1 --
2,Иш
и,Иш
■ л„,05„,г
(А-1).
(11)
^ ,Нш
1 --
^ ,Иш
■ПмА,г
(Ц-1)
где Ь -1 - расстояние между узлами г и г - 1; Ь - расстояние между узлами г и г + 1.
При выполнении расчетов было отмечено, что на формирование напряженного состояния обделки подземного сооружения значительное влияние оказывает ее вес. Так, в работе [37] было показано, что снижение давления в нижней части подземного сооружения позволяет получить лучшую сходимость с результатами прогноза, основанными на численном моделировании. При выполнении расчетов принято, что величина активного давления снижена по отношению к максимальному ее значению на 10 % [38].
Алгоритм определения формы поперечного сечения подземного сооружения по силовому фактору. Алгоритм многовариантного анализа влияния формы поперечного сечения тоннеля квазипрямоугольной формы (рис.2) на напряженное состояние обделки заключается в разбивке сечения тоннеля на отдельные сегменты, параметры (положение центра сегмента и радиус) которых контролируются через опорные точки. В первую очередь определяется площадь сечения тоннеля в соответствии с габаритом приближения строений для двухпутного тоннеля метрополитена [39, 40]. Контур сечения в свету тоннеля квазипрямоугольной формы должен перекрывать эту зону. Принимая во внимание это требование, предложен ряд математических зависимостей для определения геометрических параметров квазипрямоугольной формы сечения тоннеля с точки зрения внутренних сил, возникающих в обделке тоннеля.
Принята следующая последовательность определения геометрических параметров тоннеля квазипрямоугольного сечения (рис.3).
На начальном этапе выбирается геометрический центр О2 и вычисляется радиус R2. Центр О2 (х2, у2) дуги 2 расположен на биссектрисе отрезка прямой А1А2, которая описывается линейным уравнением у = х + (И1 - <^).
Радиус определяется по формуле
^2 =^( й2 - х2 )2 +(К- У2 )2
(12)
где ^2, И1) - координаты точки А2.
2
Рис.2. Общий алгоритм расчета геометрических параметров поперечного сечения (а) и различные конфигурации поперечного сечения (б) для тоннелей некруговой формы хА],, хА2 - горизонтальные координаты точек А1 и А2 соответственно; х2г/ - предельная горизонтальная координата точки 02; Я1, Я2, Я3 - радиусы боковой стенки, опорной части свода и свода; а1Ь а22 - верхний предельный угол и нижний предельный угол углов а1 и а2 соответственно; М,\ - минимальное значение максимального изгибающего момента при положении точки 02 соответствующем случаю , (при изменении положения точек 01 и 03)
У2
О(о,о)
О(о,о)
ЛгеИз
Oi
O
Яз
г
Archi
Оз [
Рис.3. Принципиальная схема расчета геометрических параметров тоннеля квазипрямоугольного поперечного сечения (АгсЫ, АгсЬ2, АгсИз - периметры дуги боковой стенки, опорной части свода и свода соответственно)
б
a
Определяется расположение центров О1, О3. Дуга 2 ограничена точками A1' и А^, расположенными на окружности с центром О2 и радиусом R2. Углы между прямой О2 А, О2 А^ и горизонтальной осью должны удовлетворять условию: а1 = (0, ац) и а2 = (а22, 90).
Центр О1 дуги 1 является пересечением линии, проходящей через точку О2 и имеющей наклон tg(аl), и горизонтальной оси х: О1(х1,0), х = -Уг^(а) + X. Центр О3 дуги 3 является пересечением линии, проходящей через точку О2 и имеющей наклон tg(a2), и вертикальной оси х: О3(0, у3), у3 = tg (а2)(-х2 ) + У 2. Отсюда можно выразить значения радиусов Rl и Rз (м):
ъ =>/( х4- х )2+(у4)2;
__(13)
где ( х^ , уА,) и ( , у^) - координаты точки А и А2 соответственно.
Адаптация метода НИМ для тоннеля некруговой формы поперечного сечения выполнена в ПК Ма1ЬаЬ. Алгоритм расчета параметров поперечного сечения тоннеля также является элементом общего программного кода для выполнения расчета напряженного состояния обделки. В качестве выходных данных выступают интегральные показатели напряженного состояния обделки (перерезывающая сила, изгибающий момент, продольная сила) и положение узловых точек системы после проведения оптимизации.
Схожий алгоритм оценки влияния формы поперечного сечения тоннелей на напряженное состояние обделки разработан и для тоннелей арочной формы поперечного сечения с обратным сводом (см. рис.1, б), но в настоящей работе он не рассматривается.
Обсуждение результатов. Оценка достоверности прогноза напряженного состояния обделки в рамках метода HRM. Рассмотрено формирование напряженного состояния обделки подземных сооружений круговой, квазипрямоугольной и арочной форм поперечного сечения, параметры которых представлены на рис.1 и в табл.2. Расчет напряженного состояния обделки выполнен на основании представленной выше методики. Полученные методом НИМ результаты сравнивались с результатами, полученными на основании численного моделирования взаимодействия обделки с породным массивом методом конечных элементов. Детальная постановка задачи в рамках метода конечных элементов представлена в [38]. При выполнении расчетов толщина обделки принималась равной 0,5 м. Задача решена в плоскодеформационной постановке, влияние призабойной зоны на формирование нагрузки на обделку не учитывалось. Глубина заложения (от поверхности земли до шелыги свода подземного сооружения) - 10 м, коэффициент бокового давления для всех расчетных случаев - 0,6.
Таблица 2
Геометрические параметры сечений тоннелей различной формы
Форма поперечного сечения R, м R1, м R2, м Rэ, м Площадь поперечного сечения, м2
Круговая 4,89 - - - 75,12
Квазипрямоугольная - 5,35 1,00 9,95 59,78
Арочная с обратным сводом - 2,45 8,11 5,24 69,65
Породный массив рассматривался как нелинейно-деформируемая среда с ограничением по величине максимальных касательных напряжений. Физико-механические характеристики породного массива, показывающие механическое поведение пород на рассматриваемом участке, и материала обделки: у - 18/25 кН/м3; Е - 3,6/35000 МПа; V - 0,495/0,15; с - 25,6/- кПа; ф - 16,5°/-соответственно.
Результаты расчета напряженного состояния обделки методом HRM и его сравнение с результатами, полученными на основании решения задачи методом конечных элементов, представлены
ю к
1-е 8
1,4 1
0,6 0,2 -0,2 -0,6 -1 -1,4
Угол, град
§
к
ё
о
«
о &
1,6
1,2 ■■
0,8-;
0,4
■ "'.Л /,-Л »Г к
\\ /Л /1 V
у/
60 120 180 240 Угол, град
---5
........... 6
Рис.4. Напряженное состояние в обделке: а - изгибающий момент; б - продольная сила 1 - круговой (HRM); 2 - сводчатая форма (НИМ); 3 - квазипрямоугольный (НИМ); 4 - круговой (Р1ах1з); 5 - сводчатая форма (Р1амз); 6 - квазипрямоугольный (Р1амз)
300 360
а
б
0
4
в виде эпюр изгибающих моментов и продольных сил по периметру обделки тоннеля (рис.4). Результаты показывают хорошую сходимость между различными схемами взаимодействия породного массива и обделки подземного сооружения. Наибольшее отличие в результатах между двумя методами по величине изгибающего момента наблюдается при рассмотрении тоннеля кругового сечения и составляет 8,3 %, при этом для тоннелей квазипрямоугольного и арочного сечений отличие составляет 1,2 и 2,0 % соответственно. По величине продольных сил отличие между двумя методами составляет 5,8, 5,4 и 6,0 % соответственно. С качественной точки зрения наблюдается хорошая сходимость. Выполненные расчеты при других параметрах породного массива, крепи и напряженного состояния [38] также показывают хорошую сходимость двух расчетных методов.
Можно сделать вывод, что представленный метод позволяет оценить напряженное состояние обделки с достаточной для практических расчетов точностью. С учетом того, что скорость выполнения расчетов в рамках метода НИМ на несколько порядков выше по сравнению с методами численного моделирования, его можно эффективно использовать для решения оптимизационных задач, где количество расчетных итераций может составлять от нескольких сотен до десятков тысяч.
Результаты оценки влияния сечения тоннеля криволинейного очертания на напряженное состояние обделки. Рассмотрим пример практической реализации предложенного алгоритма оценки влияния формы поперечного сечения тоннеля на напряженное состояние обделки. Инженерно-геологические условия строительства тоннеля квазипрямоугольного поперечного сечения приняты для г. Ханой, Вьетнам. Габарит приближения строений тоннеля для рассматриваемых условий представлен на рис.5. Глубина заложения тоннеля - 10 м, толщина обделки - 0,5 м. Расчет геометрических параметров сечения тоннеля выполняется согласно алгоритму, описанному ранее.
Для сокращения времени расчета каждый случай расчета оценивается с использованием шага 1° для углов а1 и а2 от ограниченных углов ап и а22 соответственно; 0,1 м для изменения параметра Х2. При расчете геометрических параметров сечения тоннеля выполнено 12479 расчетов напряженного состояния обделки квазипрямоугольного поперечного сечения, расчетное время для каждого случая составило около двух секунд.
Зависимость внутренних сил и радиальных смещений обделки от величины радиуса обратного свода Яз тоннеля представлена на рис.6. Влияние радиуса Яз на величину максимального изгибающего момента в сечении обделки - значительно. Увеличение радиуса Яз приводит к увеличению максимального изгибающего момента в обделке. Максимальная нормальная сила меняется незначительно при изменении радиуса Яз.
Также увеличение радиуса Яз приводит к тому, что свод и обратный свод тоннеля становятся более плоскими. Тогда нормальная вертикальная нагрузка, действующая на эти сегменты обделки, увеличивается, а вертикальная поперечная нагрузка - снижается. Также снижается вертикальная
Ж © М.А.Карасев, Тиен Тай Нгуен, 2022
\ 1 / г-и
760 | 1680 | 4000 1680
760
шо
Ho
Рис.5. Построение габарита приближения строений тоннеля: а - двухпутный железнодорожный тоннель по стандартам Вьетнама; б - упрощенный габарит приближения строений тоннеля = 8,88 м; Но = 5,50 м; ^ = 0,76 м; С = 3,68 м; С = 4,44 м; Н = 2,75 м; Н = 1,54 м; Н = Ео/2 = Н + ^
2
I 1,8
к ■
gffi 1 6
2 S '
га .
ю f- 1 д
Я И i,4
и и
s I U 1
0,28481п(.г) + »2 = n 0,567 7
15 25 35 45 55 Радиус Яз, м
65
з
сч га
Я ^ S За 2 йч
S я
^ ^ ё ^ I I
^ и я
1,2
0,8
0,6
0,4
15 25 35 45 55 65 Радиус Яз, м
б
33 3 2
1,6
1,4
1,2
и
о &
15
25 35 45 Радиус Яз, м
55
65
0,25
Ё 0,2
% 3»
0,15
ä од
15 25 35 45 55 Радиус Яз, м
65
б
a
h
а
1
5
5
в
г
1
5
5
Рис.6. Внутренние силы и нормальные смещения, возникающие в обделке тоннелей: а - максимальный изгибающий момент; б - максимальная продольная сила; в - максимальная перерезывающая сила; г - максимальное радиальное смещение
нагрузка, передаваемая на боковую стенку тоннеля. Таким образом, чем более плоские свод и обратный свод тоннеля, тем выше максимальный изгибающий момент в обделке тоннеля. Максимальное нормальное смещение обделки тоннеля увеличивается с увеличением радиуса Яз.
Обработка результатов расчетов позволила установить, что соотношение между радиусом Яз и максимальным изгибающим моментом в обделке тоннеля квазипрямоугольной формы поперечного сечения может быть описано уравнением
7 = 0,28481п(Я3) + 0,567 МНм/м. (14)
Связь между максимальным изгибающим моментом в обделке тоннеля и радиусами Я2, Яз, Я1 с достаточной для практических расчетов точностью можно описать уравнением
7 = 0,61038 - 0,00059 ехр (^ ) + 0,287761п (Я3)-0,025051п (Яг ) МНм/м. (15)
Из уравнений (14) и (15) видно, что на максимальный изгибающий момент в обделке в основном влияет радиус обратного свода и незначительно - радиус боковой стенки тоннеля.
Детальный анализ результатов расчетов позволил установить, что величина максимального изгибающего момента в сечении обделки тоннеля для рассматриваемых условий изменяется от 1,05 (наилучшая форма поперечного сечения по силовому фактору) до 1,85 МНм/м (наихудшая форма), т.е. разница в величинах изгибающего момента составляет 75 %. При этом разница в площади между наилучшей и наихудшей формами поперечного сечения составила 21 %. Можно заключить, что анализ данных, полученных на основании многовариантного расчета, позволил в значительной степени снизить напряженное состояние обделки при сохранении минимального размера площади поперечного сечения тоннеля.
Для оценки влияния различных параметров и геометрии тоннеля на напряженное состояние обделки был выполнен многовариантный расчет. В определенном диапазоне изменялись модуль деформации породного массива, глубина заложения тоннеля и коэффициент бокового давления. Модуль деформации грунта является одним из основных параметров, определяющих напряженное состояние обделки тоннеля. Его величина при выполнении расчетов изменялась от 5 до 200 МПа, при этом остальные физико-механические характеристики пород остались прежними. Всего выполнено 26031 расчетов. Результаты представлены в виде зависимостей изменения максимальной величины изгибающего момента и продольной силы в обделке тоннеля (рис.7, а). Как видно из представленных зависимостей, для всех значений модулей деформации породного массива верно следующее соотношение: чем выше значение Я3, тем выше абсолютное значение изгибающего момента; при этом величина продольной силы меняется незначительно. Полученные результаты хорошо соотносятся с данными, представленными в [41].
При изменении Я3 от 10 до 15 м наблюдается резкий рост величин изгибающего момента, в то время как величины максимальных нормальных сил быстро снижаются. Дальнейшее увеличение показателя Я3 не оказывает серьезного влияния на напряженное состояние обделки.
Полученные результаты можно обобщить в виде зависимости, которая характеризует связь между максимальным значением изгибающего момента, геометрическими параметрами и модулем деформации породного массива:
Мтх = 0,76941 - 0,170421п (Е„ ) - 0,004651п Ц - 0,00267ехрЯ2 + 0,335551п Я3 МН м/м. (16)
Оценка влияния начального напряженного состояния породного массива, выраженного через величину коэффициента бокового давления К0 и глубину заложения тоннеля, представлена в виде эпюр изменения интегральных показателей напряженного состояния обделки (рис.7, б). Величина коэффициента бокового давления К0 изменялась от 0,4 до 1,0. Количество выполненных расчетов составило 22778. На основании анализа результатов расчетов можно заключить, что увеличение К0 приводит к уменьшению максимальной величины изгибающего момента, однако величина продольной силы в целом увеличивается. Такая тенденция характерна для любых значений Я3. При достижении Я3 величины 10 м влияние коэффициента К0 на напряженное состояние обделки повышается. С увеличением значения Я3 эффект обжатия снижается, поэтому снижается и влияние показателя Я3 на напряженное состояние обделки. По результатам расчетов получена зависимость для определения максимальной величины изгибающего момента через величину Я3 и коэффициент бокового давления К0:
Мтах = 0,96512-0,7001^-0,02281пЦ -0,0012ехрЯ2 + 0,314811пЯ3 МН м/м. (17)
Коэффициент бокового давления К0 оказывает влияние и на формирование нормальных сил в сечении обделки, однако степень этого влияния значительно ниже влияния на изгибающий момент. Данный результат получен при различных значениях Я3.
Для оценки влияния радиуса Я3 на напряженное состояние обделки при различных глубинах заложения тоннеля выполнены дополнительные расчеты в объеме 19524 (рис.7, в). Величина глубины заложения Е изменялась от 5 до 50 м. На основании выполненных исследований установлено, что влияние параметра Я3 на величину изгибающего момента снижается по мере увеличения глубины заложения тоннеля. Это можно объяснить тем, что с ростом глубины возрастают
2 §
8 н S
2 s
'S о
15
§ 2 s S
Й | a s
2
1,8 1,6 1,4 1,2 1
0,8
0,6 0,4 0,2
1,4
1,3
M 03
3 1,2
к л
H
о
«
о &
1,1
Е =
10 20 30 40 50 60 70 Радиус Ri, м
5 -£=10 — Е = 25 — Е = 50
70
Е = 100
20 30 40 50 6С Радиус Ri, м
-£=150 -£ = 200 МПа
Ii
2 S 9. b
'S о
a s
§ 2 s S
Й | a s
1,8
1,6
1,4
1,2
0,8
0,6
20 30 40 50 60 Радиус Ri, м
70
1,5
1,4
M tf
° 1,3
к л
H
о
«
о &
1,2
1,1
20 30 40 50 Радиус Ri, м
70
K = 0,4
K = 0,5
K = 0,6
K = 0,7
K = 0,8
K = 0,9
K = 1
4,5
s 4
S к 3,5
2 S
и „ -,
3 i 3
« I
=s О 2,5
S 2
| ю 1,5
I e
0,5
3,5
10 20 30 40 50 60 70 Радиус Ri, м
3
§ * S a
2,5
2
l| 1,5
о &
G 1
0,5
10 20 30 40 50 60 70 Радиус Ri, м
D = 5
D = 10
D = 20
D = 30
D = 40
D = 50 м
Рис.7. Влияние параметра Яз на напряженное состояние обделки тоннеля в зависимости от модуля деформации грунта (а), величины коэффициента бокового давления (б) и глубины заложения тоннеля (в)
а
1
0
б
2
1
в
2
1
0
0
и активные нагрузки на обделку, что приводит к возрастанию значений интегральных показателей напряженного состояния обделки. Увеличение показателя R3 приводит к тому, что обратный свод становится более плоским, а соотношение между нормальной и касательной нагрузками возрастает, что уменьшает вертикальную нагрузку в стенках тоннеля. Для расчета величины максимального изгибающего момента в обделке, в зависимости от глубины заложения тоннеля и геометрических параметров тоннеля, используется формула
Mmax = 0,04727H+0,60766lnr + 0,00081expR + 0,00111lnR - 0,68012 МН м/м. (18)
Далее приведена формула, учитывающая влияние комбинации рассмотренных факторов: модуля деформации грунта Es, глубины тоннеля H, коэффициента давления грунта на максимальный изгибающий момент в обделке. Уравнение получено на основании анализа 956675 случаев:
Mmax = 1,10608 - 0,26759lnEs - 1,06703K + 0,03098H - 0,01048lnRi +
+ 0,00395expR2 + 0,50104lnR3 МНм/м. (19)
Как видно из представленных результатов, предложенный алгоритм позволяет эффективно оценивать влияние контролирующих параметров сечения тоннеля на напряженное состояние обделки различной формы поперечного очертания. В то же время автоматизированного способа оптимизации формы тоннеля по заданным критериями (прочность материала, устойчивость конструкции, площадь поперечного сечения и т.д.) не предложено, но будет изучаться в дальнейших исследованиях.
Заключение. Как показывает практика строительства подземных сооружений, их формы поперечного сечения могут быть круговыми, арочными с обратным сводом, квазипрямоугольными и др. Для расчета напряженного состояния обделок некруговых форм может быть использован метод HRM, позволяющий решить поставленную задачу с достаточной для практических расчетов точностью. С учетом того, что скорость выполнения расчетов на несколько порядков выше, чем в методе конечных элементов, где породный массив рассматривается как сплошная среда, описываемая соответствующей моделью ее деформирования, метод HRM может быть использован для решения оптимизационных задач - поиска оптимальной формы поперечного сечения подземного сооружения при заданных горно-геологических условиях.
На основании выполненных исследований:
• развит метод HRM для расчета подземных сооружений некруговой формы поперечного сечения с учетом нелинейного поведения породного массива;
• разработан алгоритм проведения многовариантного расчета напряженного состояния обделки тоннелей арочной и квазипрямоугольной форм поперечного сечения;
• выполнена апробация предложенного метода для различных условий строительства тоннелей, и получена хорошая сходимость с результатами численного моделирования, полученными в программном комплексе PLAXIS 2D;
• выполнен комплекс параметрических расчетов и получены зависимости, связывающие интегральные показатели напряженного состояния обделки (изгибающий момент и продольные силы) с модулем деформации породы, напряженным состоянием породного массива и геометрическими параметрами сечения тоннелей квазипрямоугольной и арочной форм;
• получены аналитические выражения для расчета максимальных значений изгибающих моментов в сечении обделки;
• указано направление дальнейших исследований по разработке алгоритма поиска оптимальной формы поперечного сечения тоннелей по заданным критериям.
ЛИТЕРАТУРА
1. TatiyaR. Civil excavations and tunnelling: a practical guide. London: Thomas Telford, 2005. 322 p.
2. Yingyongrattanakul N., Adachi T., Tateyama K., KurahashiM. Optimal shape of underground structure // Modern Tunneling Science and Technology. London: Taylor & Francis, 2001. 628 p.
3. ГосподариковА.П., ВыходцевЯ.Н., ЗацепинМ.А. Математическое моделирование воздействия сейсмовзрывных волн на горный массив, включающий выработку // Записки Горного института. 2017. Т. 226. С. 405-411. DOI: 10.25515/PMI.2017.4.405
4. Саммаль А.С., Фотиева Н.Н., Булычев Н.С., Хренов С.И. Расчет тоннельных обделок переменной толщины с учетом влияния земной поверхности // Записки Горного института. 2oo4. Т. 15б. С. 24-2б.
5. Protosenya A.G., LebedevM.O. Calculation of the Loads on Linings of Subway Tunnels Constructed in Physically Nonlinear Soil Masses // Journal of Mining Science. 2oo2. Vol. 38. Iss. 5. P. 418-424. DOI: io.io23/A:io23975313128
6. Protosenya A.G., Karasev M.A., Belyakov N.A. Elastoplastic problem for noncircular openings under Coulomb's criterion // Journal of Mining Science. 2o16. Vol. 1. P. 53-б1. DOI: 1o.1134/S1o62739116oio125
7. Belyakov N.A., Smirnova O.M., Alekseev A. V., Tan H. Numerical Simulation of the Mechanical Behavior of Fiber-Reinforced Cement Composites Subjected Dynamic Loading // Applied Sciences. 2021. Vol. 11. Iss. 3. № 1112. DOI: io.339o/app1io31112
8. Nguyen C.T., Gospodarikov A.P. Hyperstatic reaction method for calculations of tunnels with horseshoe-shaped cross-section under the impact of earthquakes // Earthquake Engineering and Engineering Vibration. 2o2o. Vol. 19. P. 179-188. DOI: io.ioo7/s118o3-o2o-o555-o
9. Fan Wang, Biancai Gou, Yawei Qin et al. Modeling tunneling-induced ground surface settlement development using a wavelet smooth relevance vector machine // Computers and Geotechnics. 2oi3. Vol. 54. P. 125-132. DOI: io.ioi6/j.compgeo.2oi3.o7.oo4
10. Безродный К.П., Лебедев М.О. О нагрузках от горного давления на обделки тоннелей закрытого способа работ // Записки Горного института. 2oi7. Т. 228. С. 649-653. DOI: io.25515/PMI.2oi7.6.649
11. Arnau O.,Molins C. Experimental and analytical study of the structural response of segmental tunnel linings based on an in-situ loading test. Part 1: Test configuration and execution // Tunnelling and Underground Space Technology. 2oii. Vol. 26. Iss. 6. P. 764-777. DOI: io.ioi6/j.tust.2oi1.o5.oo2
12. Петраков Д.Г., Пеньков Г.М., Золотухин А.Б. Экспериментальное исследование влияния горного давления на проницаемость песчаника // Записки Горного института. 2o22. Т. 254. С. 244-251. DOI: io.31897/PMI.2o22.24
13. Do N.A., Dias D., Oreste P., Maigre I.D. The behaviour of the segmental tunnel lining studied by the hyperstatic reaction method // European Journal of Environmental and Civil Engineering. 2oi4. Vol. 18. Iss. 4. P. 489-5io. DOI: io.io8o/19648189.2oi3.872583
14. Blom C.B. Design philosophy of concrete linings for tunnel in soft soils. Delft, Netherlands: Delft University, 2oo2.
15. Xiaojun Li, Zhiguo Yan, Zhen Wang et al. Experimental and analytical study on longitudinal joint opening of concrete segmental lining // Tunnelling and Underground Space Technology. 2oi5. Vol. 46. P. 52-63. DOI: io. ioi6/j.tust.2oi4.11.oo2
16. Liang Lu, Xilin Lu, Peifang Fan. Full-Ring Experimental Study of the Lining Structure of Shanghai Changjiang Tunnel // Journal of Civil Engineering and Architecture. 2oii. Vol. 5. Iss. 8. P. 732-739. DOI: io.17265/1934-7359/2oi1.o8.oo7
17. Nakamura H., Kubota T., Furukawa M., Nakao T. Unified construction of running track tunnel and crossover tunnel for Unified construction of running track tunnel and crossover tunnel for // Tunnelling and Underground Space Technology. 2oo3. Vol. 18. P. 253-262. DOI: io.ioi6/So886-7798(o3)ooo34-8
18. Xin Huang, Yeting Zhu, Zixin Zhang et al. Mechanical behaviour of segmental lining of a sub-rectangular shield tunnel under self-weight // Tunnelling and Underground Space Technology. 2oi8. Vol. 74. P. 131-144. DOI: io.ioi6/j.tust.2oi8.oi.oi6
19. Xian Liu, Zhen Liu, Yuhang Ye et al. Mechanical behavior of quasi-rectangular segmental tunnel linings: Further insights from full-scale ring tests // Tunnelling and Underground Space Technology. 2oi8. Vol. 79. P. 3o4-318. DOI: io.ioi6/j.tust.2oi8.o5.oi6
20. Xian Liu, Yuhang Ye, Zhen Liu, Dezhong Huang. Mechanical behavior of Quasi-rectangular segmental tunnel linings: First results from full-scale ring tests // Tunnelling and Underground Space Technology. 2oi8. Vol. 71. P. 44o-453. DOI: io.ioi6/j.tust.2oi7.o9.oi9
21. Trushko V.L., Protosenya A.G., Trushko O.V. Stress-strain behavior of the workings during the rich iron ores development under the confined aquifers // International Journal of Applied Engineering Research. 2oi6. Vol. 11. Iss. 23. P. 11153-11164.
22. Feng Yang, Xinlei Sun, Jiahong Zou, Xiangcou Zheng. Analysis of an elliptical tunnel affected by surcharge loading // Proceedings of the Institution of Civil Engineers - Geotechnical Engineering. 2oi9. Vol. 172. Iss. 4. P. 312-319. DOI: io.168o/jgeen.16.ooi22
23. PalassiM.,MohebbiM. Design of Lining of Tunnels Excavated in Soil and Soft Rock // The Electronic Journal of Geotech-nical Engineering. 2oo8. Vol. 3. P. 1-15.
24. Protosenya A., Vilner M. Assessment of excavation intersections' stability in jointed rock masses using the discontinuum approach // Rudarsko-geolosko-naftni zbornik. 2o22. Vol. 38 (2). P. 137-147. DOI: io.17794/rgn.2o22.2.12
25. Wang H., Zeng G., Jiang M. Analytical stress and displacement around non-circular tunnels in semi-infinite ground // Applied Mathematical Modelling. 2oi8. Vol. 63. P. 3o3-328. DOI: io.ioi6/j.apm.2oi8.o6.o43
26. Abdellah W.R., Ali M.A., Yang H.S. Studying the effect of some parameters on the stability of shallow tunnels // Journal of Sustainable Mining. 2oi8. Vol. 17. Iss. 1. P. 2o-33. DOI: io.ioi6/j.jsm.2oi8.o2.ooi
27. Rostami A., Asghari N., Ziarati M.A. et al. Investigating Effect of Tunnel Gate Investigating Effect of Tunnel Gate Inserted Forces on Its Coverage and Soil Surface Settlement // Journal of Civil Engineering. 2oi6. Vol. 6. Iss. 3. P. 358-369. DOI: io.4236/ojce.2oi6.63o3o
28. Vinod M., Khabbaz H. Comparison of rectangular and circular bored twin tunnels in weak ground // Underground Space. 2oi9. Vol. 4. Iss. 4. P. 328-339. DOI: io.ioi6/j.undsp.2oi9.o3.oo4
29. Do N.A., Dias D., Zhang Z. et al. Study on the behavior of squared and sub-rectangular tunnels using the Hyperstatic Reaction Method // Transportation Geotechnics. 2o2o. Vol. 22. № 100321. DOI: io.ioi6/j.trgeo.2o2o.ioo321
30. Oreste P. A numerical approach to the hyperstatic reaction method for the dimensioning of tunnel supports // Tunnelling and Underground Space Technology. 2oo7. Vol. 22. Iss. 2. P. 185-2o5. DOI: io.ioi6/j.tust.2oo6.o5.oo2
31. Huebner K.H., DewhirstD.L., Smith D.E., Byrom T.G. The finite element method for engineers. New York: John Wiley and Sons, 2ooi. 744 p.
32. БаклашовИ.В., Картозия Б.А. Mеханика подземных сооружений и конструкции крепей. M.: Студент, 2012. 543 c.
33. Булычев Н.С. Mеханика подземных сооружений. M.: Недра, 1994. 382 c.
34. KolymbasD. Tunnelling and Tunnel Mechanics. A rational Approach to Tunneling. Berlin: Springer, 2oo5. 439 p.
35. Széchy K. The Art of Tunnelling. Budapest: Akadémiai Kiadó, 1973. 1097 p.
36. Takano Y.H. Guidelines for the design of shield tunnel lining // Tunneling and Underground Space Technology. 2000. Vol. 15. Iss. 3. P. 303-331. DOI: 10.1016/S0886-7798(00)00058-4
37. Do N.A., Dias D. Tunnel lining design in multi-layered grounds // Tunnelling and Underground Space Technology. 2018. Vol. 81. P. 103-111. DOI: 10.1016/j.tust.2018.07.005
38. Nguyen T.T., Do N.A., Karasev M.A. et al. Influence of tunnel shape on tunnel lining behaviour // Proceedings of the Institution of Civil Engineers - Geotechnical Engineering. 2021. Vol. 174. Iss. 4. P. 355-371. DOI: 10.1680/jgeen.20.00057
39. Нгуен Т.Т., Карасев М.А. Расчет оптимальных геометрических параметров тоннеля квазипрямоугольного поперечного очертания по силовому фактору // Горный информационно-аналитический бюллетень. 2021. Т. 6. С. 59-71. DOI: 10.25018/0236_1493_2021_6_0_59
40. Du D., Dias D., Do N. Lining performance optimization of sub-rectangular tunnels using the Hyperstatic Reaction Method // Computers and Geotechnics. 2020. Vol. 117. № 103279. DOI: 10.1016/j.compgeo.2019.103279
41. Господариков А.П., Зацепин М.А. Математическое моделирование нелинейных краевых задач геомеханики // Горный журнал. 2019. № 12. С. 16-20. DOI: 10.17580/gzh.2019.12.03
Авторы: М.А.Карасев, д-р техн. наук, доцент, https://orcid.org/0000-0001-8939-0803 (Санкт-Петербургский горный университет, Санкт-Петербург, Россия), Тиен Тай Нгуен, аспирант, taitien12@gmail.com, https://orcid.org/0000-0002-5246-9252 (Санкт-Петербургский горный университет, Санкт-Петербург, Россия; Ханойский университет горного дела и геологии, Ханой, Вьетнам).
Авторы заявляют об отсутствии конфликта интересов.