# Abstract

It is known that as a result of the device of underground workings on the earth's surface the bowl-shaped hollow caused by movement of soil in the developed country is formed. The cross section of this hollow, carried out in the normal direction with respect to the extension of the production, has a saddle shape. This fact is confirmed by numerous results of theoretical researches and field observations. The article presents an analytical solution to the problem of stress distribution in a homogeneous isotropic soil mass with asymptotically damped movement on the part of its boundary. The resulting solution can be used to determine the stresses arising in the soil mass due to its subsurface mining. Methods of the theory of functions of complex variable (Kolosov-Muskhelishvili methods) were used to construct the analytical solution. Expressions for stress and strain components, depending on the Poisson ratio, are determined based on the solution of the second basic boundary value problem of the plane theory of half-plane, which is very important, because of the sufficiently large difference in the numerical values of this characteristic for different types of rock soils. Graphic images of stress contours are presented. To build the painting from lines of voltages used mathematical shell Maple. Particular cases of the above solutions are the solutions of the problems of the stress-strain state of elastic half-plane with uniform and linear displacements of the boundary of the half-plane. By summing up the stress and strain components given in the article, it is possible to obtain a solution of a similar problem in the asymptotic tendency of displacement to a certain constant value (equal-dimensional displacement).

# Full Text

Хорошо известно, что в результате проведения подземных выработок на земной поверхности образуется чашеобразная впадина, вызванная перемещением грунта в выработанное пространство. Поперечное сечение этой впадины, проведенное в нормальном по отношению к простиранию выработки направлении, имеет седлообразную форму. Это подтверждается многочисленными результатами теоретических исследований и натурных наблюдений [1-7] (рис. 1). Поэтому здания и сооружения, построенные на подрабатываемых территориях, испытывают дополнительные вертикальные и горизонтальные перемещения, учет которых необходим при проведении расчетов, связанных с оценкой исходного (существующего на момент начала проектирования) напряженно-деформированного состояния грунтовых оснований. Таким образом, расчет оснований сооружений по деформациям (по второй группе предельных состояний) является в этих условиях обязательным в соотвтетствии со CНиП 2.01.09-2010 «Здания и сооружения на подрабатываемых территориях и просадочных грунтах». Экспериментально установлено [8], что деформации грунтов под фундаментами развиваются преимущественно в верхней зоне основания, поэтому для анализа напряженно-деформированного состояния оснований сооружений могут быть применены расчетные модели, основанные на решениях задач теории упругости [9-13]. Широкое распространение в механике грунтов получила модель линейно-деформируемой среды, в рамках которой проводится анализ напряженно-деформированного состояния грунтовых массивов методами линейной теории упругости [14-16]. а б в г Рис. 1. Ширина зоны влияния строительства коммуникационного тоннеля (а) [3]; мульды оседания земной поверхности поперек простирания пласта (б) [4]; результаты прогноза оседания земной поверхности (в) [6]; графики оседания земной поверхности над квадратной выработкой размером 8´8 м, отработанной на глубине 20 м (г) [7] Fig. 1. The width of the zone of influence of construction of the communication tunnel (a) [3]; of the mould subsidence of the earth's surface across the strike of the reservoir (b) [4]; the results of prediction of subsidence of the earth's surface (c) [6]; the plots of subsidence of a terrestrial surface over a square production 8´8 m of waste at a depth of 20 m (d) [7] Если известна форма поперечника мульды оседания и представляется возможным аппроксимировать граничную линию некоторой функцией, то для отыскания напряженно-деформированного состояния грунтового массива можно воспользоваться одним из наиболее эффективных методов решения задач двумерной теории упругости, которым является метод комплексных потенциалов Г.В. Колосова - Н.И. Мусхелишвили [17, 18], основанный на применении аппарата теории функций комплексного переменного. Как отмечено выше [1-7], деформации земной поверхности, вызываемые горными выработками, представляют собой кривые, концы которых лежат на границе полуплоскости. Это означает, что перемещения на концах участков деформирования отсутствуют. В данной работе полагаем, что перемещения носят асимптотически затухающий характер, что, по нашему мнению, позволит провести более детальный анализ напряженно-деформированного состояния полуплоскости для ряда прикладных задач механики грунтов, связанных с учетом перемещений участков границы полуплоскости. Примеры решения таких задач рассмотрены в наших работах [19-22]. В данной статье в рамках модели линейно-деформируемой среды рассмотрена задача о распределении напряжений в грунтовом массиве при асимптотически затухающем перемещении участка границы на основе решения второй основной граничной задачи плоской теории упругости для полуплоскости, полученного Н.И. Мусхелишвили. Согласно работе [18, с. 353] граничное условие для второй основной задачи плоской теории упругости в случае полуплоскости имеет вид или где - граничные значения функций голоморфных в нижней полуплоскости; и - упругие постоянные, причем , где коэффициент Пуассона, связанный с коэффициентом бокового давления соотношением . Решение задачи дают формулы для функций напряжения и полученные Н.И. Мусхелишвили [18, с. 354, 355], при этом (1) (2) Пусть отрезок оси подвержен асимптотически затухающему перемещению (3) где b и g - заданные положительные действительные числа; h и l - действительные числа, причем при остальных значениях t будем полагать. Определим напряженное состояние упругой полуплоскости при перемещениях, описываемых законом (3) при всех указанных значениях параметра a. Случай α = 0. Формула (3) принимает вид (4) Заметим, что функция имеет максимум в точке с координатами и, следовательно, значение функции в точке максимума зависит лишь от величины параметра На рис. 2 приведены графики функции описывающие перемещения участка оси при различных значениях параметров, входящих в эти функции. Определим напряженное состояние в нижней полуплоскости. Формулы (1) и (2) с учетом (4) дают Рис. 2. Графики функции : а - при ; б - ; в - Fig. 2. Function graphs : а - at ; b - ; c - Следуя [18, с. 352], под выражением на отрезке будем понимать (5) где (6) (7) Тогда, согласно [19, с. 111], имеем откуда (8) где На рис. 3-5 приведены изолинии компонент напряжения, построенные на основании формул (8) при a = 10, h = l = 1 и n = 0,42 (глинистый грунт). Рис. 3. Изолинии горизонтального напряжения: а - при ; б - ; в - Fig. 3. Horizontal stress contours: a - at ; b - ; c - Рис. 4. Изолинии вертикального напряжения: а - при ; б - ; в - Fig. 4. Vertical stress contours at: a - ; b - ; c - Полагая в формулах (8) получаем при а потом при компоненты напряжения для равномерного вертикального, а затем и горизонтального перемещения рассматриваемого в задаче симметричного участка границы, т.е. Рис. 5. Изолинии касательного напряжения: а - при ; б - ; в - Fig. 5. Shear stress contours: a - at ; b - ; c - (9) (10) Формула (9) приведена в нашей работе [19] и в монографии [20], где использована при решении задачи о вычислении полной осадки основания незаглубленного ленточного фундамента с учетом напряжений, возникающих в грунтовом массиве вследствие смещения нагруженного участка границы. Для определения компонент деформации воспользуемся известными формулами [18, с. 95]: (11) которые после подстановки в них соотношений для компонент напряжений (8) дают (12) Случай α = 1. Формула (3) принимает вид (13) Функция имеет максимум в точке с координатами при. Таким образом, значение функции в точке максимума зависит от параметров и . На рис. 6 приведены графики функции , описывающие перемещения участка оси при различных значениях параметров, входящих в . Рис. 6. Графики функции : а - при ; б - ; в - Fig. 6. Function graphs : a - at ; б - ; в - Перейдем к вычислению компонент напряжения. Согласно [19, с. 111], имеем откуда (14) где На рис. 7-9 приведены изолинии компонент напряжения, построенные на основании формул (14) при и (глинистый грунт). Рис. 7. Изолинии горизонтального напряжения: а - при ; б - ; в - Fig. 7. Horizontal stress contours: a - at ; b - ; c - Рис. 8. Изолинии вертикального напряжения: а - при ; б - ; в - Fig. 8. Vertical stress contours: a - at ; b - ; c - Рис. 9. Изолинии касательного напряжения: а - при ; б - ; в - Fig. 9. Shear stress contours: a - at ; b - ; c - Пусть . Тогда, аналогично предыдущему случаю, на основании формул (14) при условии, что а потом при условии, что получаем (15) (16) Формулы (15) и (16) - компоненты напряжения для линейного вертикального и, соответственно, горизонтального перемещения участка границы оси Они приведены в работах [21, 22]. Используя формулы (11) и компоненты напряжения (14), вычислим компоненты деформации. Имеем (17) Случай α = 2. Формула (3) принимает вид (18) Нетрудно видеть, что перемещение (18) можно представить в виде суммы перемещений вида (4), т.е. (19) Решение задачи для перемещения, заданного формулой (19), можно получить путем сложения компонент напряжения, заданных формулами (8) с формулами (9) и (10), предварительно умножив компоненты напряжения на множители из правой части формулы (19). Кроме того, аналогично можно получить решение задачи о напряженно-деформированном состоянии упругой полуплоскости при перемещении вида (20) В этом случае компоненты напряжения могут быть получены сложением выражений (8)-(10) с формулами (14) после их умножения на соответствующие множители. Заметим, что величины перемещений, вычисленные по формулам (19) и (20), в пределе совпадают с равномерным перемещением вида . Выводы 1. Получены в замкнутом виде выражения для компонент напряжения и компонент деформации второй основной граничной задачи плоской теории упругости для полуплоскости при асимптотическом затухании перемещения на участке ее границы. 2. Частными случаями приведенных решений являются решения задач о напряженно-деформированном состоянии упругой полуплоскости при равномерном и линейном перемещениях участка границы полуплоскости. 3. Суммированием приведенных в статье компонент напряжения и деформации можно получить решение аналогичной задачи при асимптотическом стремлении перемещения к некоторому постоянному значению (равномерному перемещению).

# References

1. Осадки грунта, вызванные проходкой тоннелей (на примере Тегеранского метро) / Х. Чакери, А. Талибинежад, М. Мусави, Б. Юнвер // Основания, фундаменты и механика грунтов. - 2012. - № 4. - С. 12-15.
2. Лапидус Л.С. К вопросу расчета перемещений земной поверхности, вызванных подземными работами // Основания, фундаменты и механика грунтов. - 1961. - № 1. - С. 20-22.
3. Петрухин В.П., Исаев О.Н., Шарафутдинов Р.Ф. Определение зоны влияния строительства коммуникационных тоннелей // Основания, фундаменты и механика грунтов. - 2013. - № 4. - С. 24-27.
4. Спутниковый радарный интерферометрический мониторинг подработанных территорий Карагандинского угольного бассейна / Ф.К. Низаметдинов, Д.В. Мозер, Н.И. Гей, А.С. Туякбай, А.Д. Каранеева // Геоматика. - 2014. - № 4. - С. 70-77.
5. Писаренко М.В., Борисов И.Л. Использование гистехнологий для определения ожидаемых сдвижений и деформаций земной поверхности // Маркшейдерия и недропользование. - 2009. - № 1 (39). - С. 69-71.
6. Карасев М.А. Прогноз оседания земной поверхности при строительстве подземных сооружений глубокого заложения в условиях Санкт-Петербурга // Записки Горного института. - 2013. - Т. 204. - С. 248-254.
7. Численный анализ оседания земной поверхности над горизонтальными выработками / А.Н. Богомолов, Е.А. Степанова, О.А. Богомолова, Е.В. Цветкова, Е.М. Либурацков // Вестник Волгоградского государственного архитектурно-строительного университета. Серия: Строительство и архитектура. - 2016. - Вып. 45(64). - С. 12-26.
8. Далматов Б.И. Механика грунтов, основания и фундаменты. - Л.: Стройиздат, 1988. - 415 с.
9. Цытович Н.А. Механика грунтов. - М.: Госстройиздат, 1963. - 636 с.
10. Кушнер С.Г. Расчет деформаций оснований зданий и сооружений. - Запорожье: ООО «ИПО Запорожье», 2008. - 496 с.
11. Иванов П.Л. Грунты и основания гидротехнических сооружений. - М.: Высшая школа, 1985. - 447 с.
12. Флорин В.А. Основы механики грунтов. Т. 1. - Л.: Госстройиздат, 1959. - 360 с.
13. Тер-Мартиросян З.Г. Механика грунтов. - М.: Изд-во АСВ, 2009. - 551 с.
14. Партон В.З., Перлин П.И. Методы математической теории упругости. - М.: Наука, 1981. - 688 с.
15. Хан Х. Теория упругости. Основы линейной теории и ее применения. - М.: Мир, 1988. - 344 с.
16. Green A.E., Zerna W. Theoretical elasticity. - Oxford: Clareden Press, 1968. - 457 р.
17. Колосов Г.В. Применение комплексных переменных диаграмм и теории функций комплексного переменного к теории упругости. - М.: ОНТИ, 1935. - 224 с.
18. Мусхелишвили Н.И. Некоторые основные задачи математической теории упругости. - М.: Наука, 1966. - 708 с.
19. Богомолов А.Н., Ушаков А.Н. Задача о вычислении осадок ленточного фундамента // Основания, фундаменты и механика грунтов. - 2011. - № 6. - С. 2-7.
20. Богомолов А.Н., Ушаков А.Н. Методы теории функций комплексного переменного в задачах геомеханики. - Волгоград: Перемена, 2014. - 227 с.
21. Богомолов А.Н., Богомолова О.А., Ушаков А.Н. О напряженно-деформированном состоянии упругой полуплоскости при линейном сдвиге участка ее границы // Вестник Волгоградского государственного архитектурно-строительного университета. Серия: Строительство и архитектура. - 2016. - Вып. 46(65). - С. 17-26.
22. Богомолов А.Н., Ушаков А.Н. Напряженно-деформированное состояние упругой полуплоскости при линейном смещении участка ее границы // Вестник МГСУ. - 2017. - Т. 12, Вып. 2 (101). - С. 184-192.

# Statistics

Abstract - 119

### Refbacks

• There are currently no refbacks.

Copyright (c) 2018 Bogomolov A.N., Bogomolova O.A., Ushakov A.N.