CC BY
17
УПРАВЛЕНИЕ, ИНФОРМАТИКА И ВЫЧИСЛИТЕЛЬНАЯ ТЕХНИКА
УДК 539.3
Д. В. Бережной, М. К.Сагдатуллин
МОДЕЛИРОВАНИЕ ДЕФОРМИРОВАНИЯ ОБДЕЛКИ ТОННЕЛЯ МЕТРОПОЛИТЕНА, РАСПОЛОЖЕННОЙ В ГРУНТЕ, С УЧЕТОМ КОНТАКТНОГО ВЗАИМОДЕЙСТВИЯ
Ключевые слова: метод конечных элементов, контактные задачи, грунты.
Построена конечно-элементная методика расчета деформирования обделки тоннеля метрополитена, взаимодействующей с сухими грунтами. Для моделирования взаимодействия между элементами расчетной области используется специальный «контактный» конечный элемент, позволяющий учесть все случаи взаимодействия контактирующих поверхностей. Решен ряд задач деформирования кольца обделки тоннеля метрополитена, в том числе и состоящего из отдельных контактирующих между собой блоков, расположенного в грунте. Показана необходимость решения контактной задачи для моделирования взаимодействия бетонных блоков кольца обделки.
Keywords: a method offinite elements, contact problems, soils.
The finite and element method of calculation of deforming of undergroundtunnel lining, interacting with dry soil is constructed. For interaction simulation between elements of estimated area the special "contact" terminal element, allowing considering all cases of interaction of contacting surfaces is used. A number of tasks of deforming of a ringlining of a tunnel of the subway, including consisting of the separate units contacting among themselves, located in soil is solved. Need of the solution of the contact task for simulation of interaction of concrete units of a ringlining is shown.
Введение
В процессе строительства и эксплуатации тоннелей метрополитена из-за сложного рельефа грунта иногда может происходить просадка тоннелей (в связи с вымыванием грунтовыми водами слабосвязанных грунтов), в результате чего блоки обделки тоннелей трескаются, причем в нескольких местах. Сквозные трещины проходят вдоль оси тоннеля, и они практически всегда расположены в верхних и боковых блоках. Прочность обделки сохраняется за счет стальной арматуры в блоках и заклинивания частей блоков. Поэтому фактически необходимо провести расчеты для определения влияния появившихся продольных трещин в блоках обделки тоннеля на прочность конструкции (т.е. определить влияние появившихся в результате просадки трещин в блоках обделки тоннеля на возникновение новых трещин) и определить возникающие дополнительные усилия в стальной арматуре.Обделка тоннеля метрополитена составная, при деформировании ее составные части взаимодействуют между собой и с грунтом.
При решении задач подобного типа задача чаще всего применяют методики расчета нелинейного деформирования элементов конструкций [1-4] с учетом их контактного взаимодействия в рамках конечно-элементного анализа [5-8]. Кроме того, в рамках расчетной схемы обычно применяют подкрепленные и специальные конечные элементы, которые учитывают арматуру и геометрическую неоднородность конструкции
[9,10].
В работе дана постановка задачи контактного взаимодействия трехмерных элементов
конструкций, расположенных в грунте, описан контактный конечный элемент со специфическим свойствами, позволяющий реализовывать все возможные ситуации сдвига и отрыва контактирующих поверхностей, решена задача деформирования кольца обделки тоннеля метрополитена, расположенного в сложном многослойном грунте, проведены сравнения с результатами натурных экспериментов.
1. Моделирование механического контакта
Для моделирования механического контакта используется специальный контактный элемент [5,6], позволяющий учесть возможность различные случаи взаимодействия между собой контактирующих поверхностей, в частности, отрыва, проскальзывания с трением и т.д.
Общее разрешающее уравнение деформирования подконструкций с учетом механического контакта в вариационной форме имеет вид
f
\
= Е ЯИ^Т{SV}dn+\\{P}T {SV}dS
(1)
где сумма по m - сумма по объемам блоков, по ^обозначается сумма по накладкам, , ^, -соответственно объемы блоков и накладок; {о-}, , {V} - напряжения, деформации и перемещения элементарных объемов блоков;
{-Я }
Hi
накладках,
- напряжения и деформации в {g} - вектор ускорения свободного
m у и
падения (р{я} - сила тяжести), {Р} - граничная
нагрузка, действующая на части границы Я-.
Будем считать, что первоначальное обжатие контактного конечного элемента всегда существует.
Для решения нелинейной задачи на базе уравнения (1) используется итерационный метод, являющийся комбинацией метода начальных напряжений и метода дополнительной деформации. Базовым для определения к-ой итерации является следующее вариационное уравнение
+гнсндун) аок = xI шр {я}' (зуо + (2)
т V От
л
{рг {д}* + х Нден +ткн5Гн),
к о, / ок
где величины —Н и в (2) являются "начальными напряжениями".
Для реализации описанной в предыдущем параграфе математической модели взаимодействия накладок в рамках МКЭ удобно определить так называемый контактный элемент.
Вводятся аппроксимации лицевых поверхностей
?(-)(#,7)=;£ 4-1N (#,7),
1=1
(3)
где Ni (#,7) - известные функции формы для двумерной аппроксимации. Для аппроксимации вектора перемещений будем использовать представление, аналогичное (3), т.е. введем
I?(-)(#,7) = ХУ2;-1 N.. (#,7) ,
,=1
I?<+)(#,7) = ХУ2,М, (#,7) .
,=1
В процессе деформирования первоначально параллельные лицевые поверхности г1 ' и ' перестают быть таковыми и степень их относительных поворотов в процессе деформирования может достигать большой величины. Поэтому все геометрические, кинематические и силовые характеристики будем определять на обеих лицевых поверхностях самостоятельно. Другими словами, напряженно-деформированное состояние будем определять самостоятельно в каждой накладке (примыкающих,
?(+)), что
соответственно, к поверхностям г и позволит более верно моделировать их состоянии при проскальзывании друг относительно друга.
Используемая итерационная процедура типа "метода начальных напряжений" представляет собой следующую последовательность действий. Первое приближение определяется из решения вариационного уравнения принципа виртуальных
перемещений в предположении справедливости закона Гука:
Ш—Г {де}аУк =
к Ук
= ХЯМ д}аУк +£Ц{г}Т д}^
к Ук к Як
где {—} ,{е} - векторы напряжений и деформаций на к -ом трехмерном фрагменте грунта; {2} ,{Р} -
векторы массовой и поверхностной нагрузок.
В данном уравнении предполагается, что кинематические связи между фрагментами (условия непрерывности перемещений) и кинематические граничные условия выполняются априори. Уравнения равновесия для каждого фрагмента, статические условия сопряжения и статические граничные условия выполняются автоматически в интегральном смысле.
В этом случае разрешающее уравнение линейно и записывается в виде:
к Ук
ХЯМ {ди}Ук +ХЯ{Р}Т
к Ук к Як
Все последующие шаги итераций основаны на линейных уравнениях для приращений
+1}}' [^]{де}аУк = ^Я —Г {де}У ,
к Ук к Ук
из решения которых находятся
{Ди(к+1'}, (Де(к+1'}, {Д—(к+1)} = [Б] {е+1)}. (4)
Истинное деформированное состояние определяется как
{и(к+1' } = {и(к>} + {Ди(к+1'}, {к+1'} = {к)} + (Де(к+1'}.
"Пробные" напряжения для анализа возможного предельного состояния и определения истинных и дополнительных напряжений находятся как
{-к+1' } = {--к > } + {Д-( к+1'}, (5)
где {—'г )} - истинные напряжения, соответствующие
принятой теории прочности. Для их вычисления по напряжениям (5) строится "упругопластическая матрица" и определяются
{—Г } = {—} + [ ]{}. (6)
В правой части уравнения (4) фигурируют так называемые "начальные" или "дополнительные" напряжения, которые вводятся как разности "пробных" напряжений (5) и истинных (6)
{Д—(к+1' } = {—(к>}-{—-к+1'}.
Переход от вариационной задачи к алгебраической производится посредством дискретизации методом конечных элементов.
2. Геометрическое и механическое моделирование задачи
Решается задача упругого деформирования бетонного кольца обделки метрополитена,
,=1
расположенного в толще грунтового массива. Предполагается, что кольцо обделки, поперечное сечение которого показано на рис. 1, состоит из восьми сегментных блоков, соединенных друг с другом посредством резиновых вставок.
В общем случае взаимное расположение сегментных блоков в кольце обделки (по сути дела положение замкового блока А) может быть произвольным. На рис. 1а приведен случай нижнего расположения замкового блока с нумерацией блоков и стыков между ними.
Замковой блок А имеет угол раскрытия 16.363636 градуса, а призамковые блоки В и основные блоки С имеют угол раскрытия 49.090909 градуса. В расчете принимается, что блоки В и С являются идентичными, а стыковка всех блоков производится одинаково.
Рис. 1 - Схема разбиения обделки на блоки и схема зоны стыка блоков
Рис.1Ь иллюстрирует взаимное
расположение резиновой вставки и концов сегментных блоков при их стыковке. Принимается, что в начальном положении блоки расположены друг относительно друга на расстоянии 2 мм, а зазоры между блоками по внешнему периметру обделки заполнены окружающими кольцо обделки грунтовыми массами.
взаимодействии с окружающим грунтом под действием собственного веса (грунта и обделки).
На втором этапе решается задача трансформирования расчетной области при мгновенном появлении трещины (или трещин) в некоторых блоках обделки. Фактически, на этом этапе решения исключаются из рассмотрения элементы трещины (разреза), то есть модифицируется глобальная матрица жесткости конструкции и происходит перераспределение равновесного состояния с учетом геометрической нелинейности. Варианты возможного
растрескивания блоков приведены на рис. 2.
Расчет проводится на основе двумерных квадратичных 8-узловых конечных элементов сплошной среды Серендипова семейства, узловыми неизвестными которого являются проекции вектора перемещений на координатные оси рабочей плоскости.
Считается, что конструкция работает в условиях плоской деформации, а влияние арматуры учитывается путем осреднения ее прочностных свойств по длине тоннеля. На рис.3 приведена схема расположения основных пучков арматуры в обычных (В и С) и замковом (А) блоках обделки соответственно. Расчетная схема предусматривает возможность смыкания бетонных блоков обделки в процессе деформирования.
Рис. 3 - Схемы расположения пучков в бетонных блоках.
арматурных
Таблица 1 - Механические материалов расчетной области
характеристики
Рис. 2 - Схемы растрескивания блоков
Расчет проводится в два этапа. На первом этапе решается задача деформирования блоков обделки тоннеля при их
Материал или грунт Ус [кг/м3] Е, [Мпа] М
Стальная арматура 7800 200000 0.3
Резина 1000 2 0.48
Бетон 2500 35000 0.2
Насыпной грунт 1850 4 0.35
неоднородный (НС)
Глина полутвердая (2А) 1960 20 0.42
Глина мягко-пластичная 1940 7 0.42
(2В)
Суглинок твердый, 1980 6.7 0.35
полутвердый (3А)
Супесь пластичная (4А) 2060 4.5 0.30
Супесь текучая (4Б) 2010 4.5 0.30
Песок пылеватый 2040 22 0.30
водонасыщенный (5 А)
Песок мелкий 2030 33 0.30
водонасыщенный (6А)
Песок средней крупности 1990 33 0.30
водонасыщенный (7 А)
Песок гравелистый 2070 40 0.30
водонасыщенный (8А)
Разрез грунтов на первом аварийном участке приведен на рис. 4.
Рис. 4 -Разрез первого аварийного участка
Кинематические граничные условия определяются следующим образом:
- по боковым границам задается отсутствие горизонтальных смещений.
- по нижней границе определяются нулевые вертикальные смещения.
3. Результаты расчетов
В табл.2 приведены максимальные окружные напряжения в бетоне и в арматуре для трех аварийных участков. Можно отметить, что максимальные сжимающие напряжения в арматуре всегда возникают во внешнем пучке арматуры (ближе к наружной поверхности тоннеля) в верхнем блоке в вертикальном (центральном) разрезе.
Таблица 2 -Окружные напряжения в бетоне и сжимающие напряжения в арматуре в зоне трех аварийных участков
Вариант расчета (растрескивания) _max У в окр в бетоне [кГ/см2] _min ( окр в окр в бетоне [кГ/см2] _max ( окр в окр в арматуре [кГ/см2]
Первый аварийный участок
(рис. 1a) 8.30 -12.1 -324
3-1 (рис.2a) 5.18 -8.76 -618
3-3 (рис.2b) 5.11 -8.75 -610
3-7 (рис.2c) 4.59 -8.64 -598
5-1 (рис.2d) 5.17 -10.36 -623
5-3 (рис^) 5.10 -10.27 -611
5-7 (рис.21) 4.61 -10.09 -600
Второй аварийный участок
(рис. 1a) 7.49 -11.2 -299
3-1 (рис.2a) 5.14 -9.04 -739
3-3 (рис.2b) 4.88 -8.67 -718
3-7 (рис.2c) 4.43 -8.54 -702
5-1 (рис.2а) 5.02 -9.93 -741
5-3 (рис^) 4.87 -9.81 -721
5-7 (рис.21) 4.41 -9.72 -703
Третий аварийный участок
(рис. 1a) 6.56 -10.8 -290
3-1 (рис.2a) 5.12 -8.76 -709
3-3 (рис.2b) 4.71 -8.69 -702
3-7 (рис.2c) 4.38 -8.53 -682
5-1 (рис.2а) 4.94 -10.14 -711
5-3 (рис^) 4.72 -10.03 -694
5-7 (рис.21) 4.43 -9.92 -685
Можно отметить, что рассматривались напряжения в некотором удалении от стыков трещин (разрезов).
Заключение
Для задачи деформирования обделки тоннелей метрополитена, представленной в виде блочного кольца с учетом контактного взаимодействия между блоками, находящейся под действием веса грунтовых массивов и деформирующейся в результате внезапного растрескивания, определены окружные напряжения в блоках обделки и растягивающие напряжения в арматуре. Материал арматуры не переходит в пластическое состояние. Появление трещин в некоторой степени понижает уровень напряжений в бетонных блоках обделки тоннелей. Как показывают расчеты, наиболее опасным с точки зрения прочности является участок тоннеля в зоне первого аварийного участка.
Предложенный метод решения задач механики с конкретными приложениями относится к современной технологии научного сопровождения, проектирования и строительства сложных объектов. Его использование позволяет проследить за изменением напряженно-деформированного состояния и поля перемещений структурно изменяющейся расчетной области от начала и до конца строительства. Это позволяет более точно и технически грамотно принимать проектные решения для различных этапов строительных работ, что зачастую нельзя сделать, опираясь только на существующие СНиПы.
Работа выполнена при частичной финансовой поддержке РФФИ в рамках научных проектов № 12-01-00955, № 12-01-97026, № 12-0131212, № 13-97057, № 13-01-97058.
Литература
1. Golovanov A.I., Sultanov L.U.Postbuckling elastoplastic state analysis of three-dimensional bodies taking into account finite strains. - Russian Aeronautics. - 2008. - Т. 51. - № 4. С. 362-368. Volume 51, Issue 4, 2008, Pages 362-368. doi: 10.3103/S106879980804003X
2. Badriev, I.B., Banderov, V.V., Zadvornov, O.A. On the equilibrium problem of a soft network shell in the presence of several point loads // (2013) Applied Mechanics and Materials, 392, pp. 188-190.
3. Каюмов Р.А., Шакирзянов Ф.Р. Моделирование поведения и оценка Несущей способности системы тонкостенная конструкция-грунт с уче-том ползучести и деградации грунта // Учен. зап. Казан. ун-та. Сер. Физ.-матем. науки. — 2011. — Т.153, №4. — С.67-75.
4. Васидзу К. Вариационные методы в теории упругости и пластичности. - М.: Мир, 1987, 542 с.
5. Бережной Д.В., Голованов А.И., Луканкин С.А., Секаева Л.Р. Моделирование поведения железобетонной обделки тоннеля в деформируемом грунте с учетом одностороннего контактного взаимодействия ее блоков через упругие прокладки // Вестник Казанского государственного технического университета, 2010, №2, С.4-9.
6. Бережной Д.В., Кузнецова И.С., Саченков А. А. Моделирование пластического деформирования многослойного грунта в зоне опоры многопролетного
моста // Учен. зап. Казан. ун-та. Сер. Физ.-матем. науки. - 2010. - Т. 152, кн. 1. - С. 116-125.
7. Бережной Д.В., Карамов А.В., Сагдатуллин М.К. Моделирование деформирования поэтапной выемки грунта при строительстве подземных сооружений // Вестник Казанского государственного технологического университета. - 2012. - №17. -С.137-145.
8. Бережной Д.В., Сагдатуллин М.К., Султанов Л.У. Моделирование деформирования обделки тоннеля метрополитена, расположенной в грунте сложной
физической природы // Вестник Казанского технологического университета. 2013. Т. 16. № 9. С. 250-255.
9. Голованов А.И., Бережной Д.В. Метод конечных элементов в механике деформируемых твердых тел. Казань, 2001, 300 с.
10. Бережной Д.В., Сагдатуллин М.К., Саченков А.А. Универсальный конечный элемент для расчета комбинированных конструкций // Вестник Казанского государственного технологического университета. -2012. - №17. - С.150-157.
© Д. В. Бережной- канд. физ.-мат. наук, доцент каф. теоретической механики, Казанский федеральный университет, [email protected]; М. К. Сагдатуллин - канд. физ.-мат. наук, доцент каф. теоретической механики и сопротивления материалов, Казанский национальный исследовательский технологический университет, [email protected].
© D. V. Berezhnoi-candidate of physico-mathematical sciences, associate professor Department of Theoretical Mechanics, Kazan federal university, [email protected]; M. K. Sagdatullin-candidate of physico-mathematical sciences, associate professor department of theoreticalmechanics andstrength of materials, Kazan National Research Technological University, [email protected].
