CC BY
23
УДК 622.28+550.34
Деев Петр Вячеславович
доктор технических наук, доцент, Тульский государственный университет, 300012, г. Тула, пр. Ленина, 92 e-mail: dodvsva@vandex.ru
DOI: 10.18454/2313-1586.2017.04.020 Deev Petr V.
doctor of technical sciences, assistant professor, Tula State University, 300012 Tula, 92 Lenin av. e-mail: dodysya@yandex.ru
Тормышева Ольга Александровна
кандидат технических наук, доцент, Тульский государственный университет e-mail: knopka.85@bk.ru
Петрухин Максим Андреевич
аспирант,
Тульский государственный университет e-mail: mpetruhin@mail.ru
ВЛИЯНИЕ ГРАНИЦЫ РАЗДЕЛА ПОРОД НА НАПРЯЖЕННОЕ СОСТОЯНИЕ КРЕПИ ГОРНОЙ ВЫРАБОТКИ ПРИ СЕЙСМИЧЕСКОМ ВОЗДЕЙСТВИИ*
Tormisheva Olga A.
candidate of technical sciences, assistant professor, Tula State University e-mail: knopka.85@bk.ru
Petruhin Maksim A.
post-graduate,
Tula State University
e-mail: mpetruhin@mail.ru
INFLUENCE OF ROCK INTERFACE ON STRESSED STATE OF MINE SUPPORT SUBJECTED TO SEISMIC ACTION
2 s
X
re S a <u EC О
U
Аннотация:
Предложена методика оценки напряженного состояния подземных конструкций при сейсмических воздействиях землетрясений, позволяющая учитывать отражение и преломление сейсмических волн на границе раздела пород с разными деформационными характеристиками. В основе методики лежит аналитическое решение плоской квазистатической задачи теории упругости о напряженном состоянии кольца, подкрепляющего отверстие вблизи линии контакта двух полуплоскостей. Среда, имеющая отверстие, нагружена на бесконечности и прямолинейной границе напряжениями, взятыми из известного решения динамической задачи об отражении и преломлении плоских гармонических волн на границе раздела. Решение рассматриваемой квазистатической задачи получено с использованием теории аналитических функций комплексного переменного, интегралов типа Коши и комплексных рядов. Решение сведено к итерационному процессу, на каждом шаге которого уточняется влияние границы раздела на напряженное состояние кольца, моделирующего крепь выработки. Выявлены условия, при которых наблюдается возрастание амплитуд горизонтальных напряжений на границе раздела и сейсмических напряжений в подземной конструкции. Получены зависимости амплитуд напряжений на линии контакта от угла падения сейсмических волн при разных характеристиках контактирующих сред. Выполнено сравнение результатов расчета крепи круговой выработки, полученных с использованием предлагаемого подхода и стандартной методики. Использование разработанной методики позволило выявить возможность значительного увеличения сейсмической нагрузки на крепь горной выработки при наличии вблизи нее границы раздела разных пород. Сделаны предложения по распространению используемого подхода на расчет крепи выработок некругового поперечного сечения и подземные сооружения, расположенные на небольшой глубине.
Ключевые слова: крепь, напряженное состояние, землетрясение, граница раздела, квазистатическая задача, теория упругости
Abstract:
The methods estimating the stressed state of underground structures under seismic action of earthquakes and taking into account reflection and refraction of seismic waves at the interface of rocks with different deformation properties are proposed. The methods are based on an analytical solution of the plane pseudo-static elasticity theory problem on stressed state of the ring supporting an opening near the interface of two half-planes. The half-plane with the opening is loaded at infinity and along straight boundary by stresses found from known solution of the dynamic problem on reflection and refraction of plane harmonic waves at interface. Solution of the problem has been obtained with the help of complex analytical functions theory, cauchy type integrals and complex series. The solution is reduced to an iteration process specifying interface influence on every step of approximation. The favourable conditions for growth of stress amplitudes at the interface and stresses in an underground structure are discovered. The dependencies of stress amplitudes on wave incident angle have been obtained. The design results obtained for a circle opening support are compared with the data received by the standard technique. The developed methods permit to find the possibility of significant i seismic load increasing on the underground structure located near the rock interface. It is proposed to apply the approach used for the design of non-circular shape supports and shallow underground structures.
Key words: support, stress state, earthquake, interface, pseudo-static problem, elasticity theory.
' Исследования выполнены в рамках гранта Президента РФ МД-1280.2017.5
Наличие в горном массиве протяженных границ, разделяющих породы с существенно отличающимися деформационными характеристиками, приводит к отражению и преломлению сейсмических волн, распространяющихся от очага землетрясения. Возникновение отраженных волн существенным образом влияет на напряженное состояние массива пород, а также подземных конструкций, расположенных вблизи границы раздела. На напряженное состояние подземных сооружений мелкого заложения большое влияние могут оказывать также сейсмические волны, отраженные от земной поверхности.
Существующие в настоящее время отечественные методы расчета подземных сооружений на сейсмические воздействия землетрясений [1, 2] не позволяют учесть влияние преломленных и отраженных волн на напряженное состояние подземных конструкций. Зарубежные подходы к сейсмическому расчету подземных сооружений [3, 4] принципиально отличаются от методов, используемых в России, однако также не рассматривают воздействие отраженных и преломленных волн на подземные сооружения.
В настоящей работе предлагается методика расчета подземных конструкций на сейсмические воздействия землетрясений, позволяющая учесть действие отраженных и преломленных волн на основе использования решений динамической и квазистатической задач теории упругости.
Напряженное состояние массива пород вблизи границы раздела предлагается определять с помощью известного решения динамической задачи об отражении и преломлении сейсмических волн на границе раздела двух упругих сред (рис. 1).
Рис. 1 - Отражение и преломление сейсмических волн на границе раздела
Здесь плоская продольная (Р) или поперечная (£К) волна падает под углом, соответственно, а или Р на границу раздела сред £о/ (/ = 1, 2), свойства которых характеризуются модулями упругости Ео/, коэффициентами Пуассона Уо/ и объемным весом уо/. В общем случае одна волна, падая на границу Ь, приводит к образованию двух отраженных и двух преломленных волн.
Перемещения точек сред £о/ выражаются через скалярный и векторный потенциалы, представляемые в виде экспоненциальных функций ф/ и определяемых формулами:
2 5 X
(С
£ а ш
ф _ ^ eik1(/1X-ЩУ-C\t). ^ _ A eik2(l2X-n2У)•
ф^ _ Aeik1(/1^+П1У+C1t); - Aeik2(/2^+Щ2У+c2¿).
_ J (/3.X-П3y-C3¿ ). _ j (/4.X-П4y-C4t)
(1)
где Ai, A2, ... , Aó - амплитуды волн; ш - круговая частота; ci, ci и ci7, c2 - скорости продольной и поперечной волн в средах So,i и So,2;
k - ш / c;k2 -ш / c2;k - ш / q;k4 -ш / c2; /х - sin а; /2 - sin р; щ - cos а; щ - cos р;
/3 - sin а/; /4 - sin р/; щ - cos а/; щ - cos р/. (2)
Граничные условия задачи отражают равенство нормальных и касательных напряжений на границе L в любой момент времени, при этом отыскиваемое решение сводится к выражению амплитуд отраженных и преломленных волн через амплитуды волновых потенциалов в падающей волне (Ai или A2). Полное решение задачи приведено в работе [5].
В результате решения определяется напряженно-деформированное состояние контактирующих сред при падении на границу раздела продольной или поперечной волны. Амплитуды возникающих напряжений зависят от угла падения волн и расстояния от рассматриваемой точки до границы раздела.
7Г/12 7Г/6 71/4 7Г/3 571/12 р
Рис. 2 - Зависимости амплитуд горизонтальных (а), вертикальных (б) и касательных (в) напряжений на границе раздела от угла падения волн
На рис. 2 показаны зависимости амплитуд напряжений на линии раздела от угла падения продольной и поперечной волн, полученные при следующих исходных данных:
2 X
ге *
а
ф
ЕС
о и
Ео,2 / Ео,1 = 0,1; Уо,1 = Уо,2 = 0,25; уод = уо,2. Цифрами «1» и «2» обозначены кривые, соответствующие значениям горизонтальных напряжений в нижней и верхней полуплоскостях. Пунктирные линии соответствуют напряжениям в сплошной среде.
Из зависимостей, представленных на рис. 2, видно, что наличие границы раздела оказывает существенное влияние на напряженное состояние контактирующих сред, причем в некоторых случаях отмечается значительное увеличение напряжений вблизи границы. Так, в рассматриваемом случае можно заметить почти двукратный рост горизонтальных напряжений в нижней полуплоскости при углах падения поперечной волны, близких к критическому углу, т. е. углу падения волны, при котором угол преломления равен 900 (см. рис. 2 а).
Для обеспечения безопасности подземных сооружений важно определить условия, при которых возможно усиление сейсмического воздействия. Рассмотрение различных комбинаций исходных данных позволило выявить факторы, существенно влияющие на величину амплитуд напряжений вблизи границы раздела двух сред при прохождении через границу продольных и поперечных сейсмических волн.
а, / РГ
Рис. 3 - Зависимости амплитуд горизонтальных напряжений от угла падения волн, иллюстрирующие сделанные выводы
2 г х
а *
а
ш
Ниже перечислены условия, при которых отмечается возрастание амплитуд напряжений вблизи границы раздела. Зависимости, иллюстрирующие сделанные выводы, представлены на рис. 3, применяемые обозначения аналогичны использованным на рис. 2.
1. Модуль упругости среды, в которую преломляются волны, больше модуля упругости второй среды £0,2 > £о,1 (рис. 3 а, уо,2 = уод; £0,2/ £0,1 = 2,0; Уо,2 = Уо,1 = 0,25).
2. Волны распространяются из среды с большим объемным весом в менее плотную среду уо,2 < уо,1 (рис. 3 б, уо,2 = 0,67уо,1; £0,2 = £0,1; Уо,2 = Уо,1 = 0,25).
3. Вторая среда имеет больший коэффициент Пуассона Уо,2 > Уод (рис. 3 в, уо,2 = уо,1; £0,2 = £0,1; Уо,1 = 0,20; Уо,2 = 0,35). Возрастание напряжений наблюдается только в случае распространения продольных волн.
4. Модуль упругости второй среды значительно ниже модуля упругости первой £0,2 << £0,1. Увеличение амплитуды напряжений наблюдается только при распространении поперечных волн. Например, в случае, когда поперечная волна падает на свободную поверхность под углом, близким к критическому, отмечается шестикратный рост амплитуды горизонтальных напряжений (рис. 3 г, £0,2 = 0; Уо,1 = 0,25).
В целом, как показали выполненные исследования, влияние границы раздела пород может привести к значительному увеличению горизонтальных напряжений вблизи линии контакта по сравнению с сейсмическими напряжениями, возникающими в однородном массиве пород.
Напряжения в падающей волне могут быть определены по формулам [6]:
Р = ^ ЛК1е1Т0; 5 = ^ ЛК1уе2Т0, (2)
2л 2л;
где А - коэффициент, зависящий от интенсивности землетрясения; К1 - коэффициент, учитывающий допускаемые повреждения подземных сооружений; у- объемный вес горных пород; То - преобладающий период колебаний частиц грунта; С1, С2 - скорости продольных и поперечных волн.
Решение рассматриваемой динамической задачи позволяет определить напряжения, которые могут быть использованы для сейсмического расчета подземных конструкций, расположенных вблизи линии контакта разных пород. С этой целью рассматривается плоская квазистатическая задача теории упругости, схема которой представлена на рис. 4. Использование квазистатического подхода в рассматриваемом случае оправдано тем, что наибольшую опасность для подземного сооружения представляют длинные сейсмические волны, длина которых намного превосходит поперечные размеры подземного сооружения. Сравнение решений квазистатической и динамической задач, подтверждающих указанное предположение, приводится в работах [7 - 9].
Рис. 4 - Расчетная схема крепи выработки, расположенной вблизи границы раздела пород
Здесь прямолинейная граница L разделяет упругие полуплоскости £0,1 и 50,2. В нижней полуплоскости £0,1 имеется круговое отверстие, подкрепленное кольцом 51. Деформационные свойства каждой из рассматриваемых сред характеризуются модулем деформации и коэффициентом Пуассона. Нижняя полуплоскость нагружена на бесконеч-
Т ^(0,1) _(0Д) „.(0,1)
ности и на границе L напряжениями ст^. , СТ^ , Т , полученными из решения динамической задачи, рассмотренной ранее. Указанные напряжения определяются в точке, расположенной на расстоянии Н от границы раздела сред при угле падения и фазе волны, соответствующих наиболее невыгодному напряженному состоянию, которое может возникнуть в кольце £1, моделирующем обделку тоннеля. На линиях контакта L и Lo выполняются условия непрерывности смещений и напряжений, внутренний контур кольца Ll свободен от действия внешних сил.
Решение рассматриваемой задачи получено с использованием теории аналитических функций комплексного переменного [10], интегралов типа Коши и комплексных рядов. Решение сводится к итерационному процессу, в каждом приближении которого уточняется влияние границы раздела на напряженное состояние кругового кольца, моделирующего крепь выработки. Ниже приводятся результаты расчета крепи круговой выработки радиусом Я0, центр которой расположен на расстоянии 1,5Я0 от границы раздела пород. Толщина крепи принята равной 0,1Яо. Отношения модулей упругости материала обделки и горных пород приняты равными £1; £0,1 : £0,2 = 30 : 10 : 1, коэффициенты Пуассона и = 0,2, ю,1 = 0,25, ю,з = 0,35, отношение объемного веса пород 70,2/70,1 = 2/3. Углы падения продольной и поперечной волн приняты равными а= Р = 35,30. При заданном угле падения волн на границе раздела возникают максимальные горизонтальные напряжения.
Эпюры напряжений сте, возникающих на внутреннем и наружном контурах поперечного сечения крепи выработки при распространении в массиве продольных и поперечных волн, представлены на рис. 5. Для сравнения пунктирными линиями даны напряжения в крепи выработки, пройденной в однородном массиве. Расчетные значения напряжений, моделирующие действие сейсмических волн, и прикладываемых к среде £0,1 на границе L и на бесконечности, приняты следующими:
- при действии продольной волны
ст(0Д) = 0,63Р; ст(0Д) = 0,37Р; т^ = 0,14Р;
- при действии поперечной волны
ст(01 = 2,365; ст^1 = -0,325 т(°Д) = 0,375.
Следует принимать во внимание, что при рассмотрении сейсмического воздействия в противоположной фазе расчетные напряжения, эпюры которых приведены на рис. 5, изменят свой знак на противоположный.
Из результатов расчета, представленных на рис. 5, видно, что в рассматриваемом случае влияние границы раздела пород приводит к существенному изменению напряженного состояния крепи выработки. При этом максимальные напряжения, возникающие в крепи при действии продольных волн, несколько меньше аналогичных напряжений, возникающих в случае, когда выработка расположена в однородном массиве. При действии поперечной волны максимальные напряжения в крепи в два с половиной раза превышают аналогичные напряжения, возникающие в случае, когда выработка расположена в однородном массиве пород.
Предложенная методика дает возможность учитывать возрастание сейсмической нагрузки на подземные сооружения, обусловленное влиянием отражения длинных сейсмических волн от границ раздела пород или земной поверхности, также может учитываться преломление сейсмических волн на границе раздела. Методика вводится в рамках
о.
ф
ЕС
квазистатического подхода к расчету подземных сооружений на сейсмические воздействия землетрясений, принятого в Российской Федерации, и снимает некоторые противоречия, касающиеся применения рассматриваемого подхода к расчету объектов мелкого заложения. Использование методики при практическом проектировании должно способствовать повышению сейсмостойкости объектов, сооружаемых в сейсмически опасных районах.
Рис. 5 - Напряжения а0 на внутреннем (а) и наружном (б) контурах поперечного сечения крепи
Найденное аналитическое решение рассматриваемой плоской квазистатической задачи теории упругости позволяет определить напряженное состояние только крепи круговой выработки, однако получение решения задачи о напряженном состоянии кольца произвольной формы, моделирующем крепь выработки некругового сечения, не должно вызвать принципиальных затруднений.
Целесообразным представляется использование предложенного подхода для определения сейсмонапряженного состояния крепи выработки, расположенной в слое породы или в приповерхностном слое относительно слабых грунтов.
2 s х
ГС *
а
ш
Литература
1. Саммаль А.С. Расчет обделок комплексов тоннелей мелкого заложения на сейсмические воздействия землетрясений / А.С. Саммаль, С.В. Анциферов, П.В. Деев // Современные проблемы гидрогеологии, инженерной гидрогеологии и гидроэкологии: матер. всеросс. конф. / Томский политехнический университет. - Томск, 2015. -С. 366 - 370.
2. Деев П.В. Расчет обделок параллельных транспортных тоннелей на действие длинных сейсмических волн заданного направления / П.В. Деев // Сб. докл. IX междунар. конф. ПГУПС, Санкт-Петербург 27 - 28 мая 2014 г. - СПб.: ПГУПС, 2014. - С. 59 - 65.
3. Hashash Yu.M.A. Seismic design and analysis of underground structures / Ha-shash Yu.M.A. [et. al]] // Tunneling and underground space technology. - Vol. 16. - 2001. -P. 247 - 293.
4. Rostami A. Evaluation of seismic behavior and earth's surface acceleration, by interaction of tunnels with different shapes and different types of soils / A. Rostami [et al.] // Open journal of civil engineering - No 6. - 2016. - P. 242 - 253.
5. Miklowitz J. The theory of elastic waves and waveguides. Noth-Holland series in applied mathematics and mechanics. - Vol. 22. - 1978. - 618 p.
6. СП 268.1325800.2016. Транспортные сооружения в сейсмических районах. Правила проектирования. - М., 2017. - 107 с. [Электронный ресурс] - Режим доступа: http://docs.cntd.ru/document/456044283
7. Pao, J.H. Dynamic stress concentration in elastic plate with rigid circular inclusion / J.H. Pao, C.C. Mow // Proc. of IV US national congress of applied mechanics. - Berkeley, USA, 1962. - P. 335 - 345.
8. Lemos, J.V. Numerical and analytical modelling of the seismic behaviour of shallow tunnels / J.V. Lemos [et al.] // Proc. of the int. symp. on modern tunneling science and technology. Kyoto, Japan. 30.10-01.11.2001. - Vol.2. - P. 675 - 680.
9. Саммаль А.С. О применимости решений квазистатических задач для определения напряженного состояния крепи некруговых выработок при действии длинных сейсмических волн / А.С. Саммаль и др. // Матер. 13-го междунар. симп. 25 - 29.09.2015 Белгород, Россия. - С. 247 - 253.
10. Мусхелишвили Н.И. Некоторые основные задачи математической теории упругости / Н.И. Мусхелишвили. - М.: Наука, 1966.
