Научный журнал
Современные наукоемкие технологии
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,899

ЗАДАЧА О ВЫЧИСЛЕНИИ НАПРЯЖЕНИЙ В СЛОИСТОМ МАССИВЕ С ВЕРТИКАЛЬНОЙ ЦИЛИНДРИЧЕСКОЙ ШАХТОЙ

Бобылева Т.Н. 1 Голобоков О.С. 1
1 Национальный исследовательский Московский государственный строительный университет (НИУ МГСУ)
В работе рассмотрена задача о напряженном состоянии неоднородного полубесконечного горного массива, состоящего из попарно чередующихся изотропных слоев, параллельных горизонтальной плоскости. В массиве имеется вертикальная цилиндрическая полость кругового поперечного сечения. С помощью метода усреднения исходная задача для неоднородного массива, уравнения которой содержат быстро осциллирующие коэффициенты, заменяется аналогичной задачей для однородной среды с постоянными (эффективными) коэффициентами. В то время как исходные слои являлись изотропными, полученная однородная среда представляет собой трансверсально-изотропный материал. Для определения напряжений было использовано имеющееся в литературе решение о напряженно-деформированном состоянии трансверсально-изотропного массива с вертикальной цилиндрической шахтой кругового поперечного сечения. В статье рассмотрены два случая: слои, составляющие массив, являются чисто упругими и упруго-ползучими. Упруго-ползучие слои описываются линейной интегральной моделью Больцмана – Вольтерра, в качестве ядер которой выбраны экспоненциальные функции, содержащие две константы. Применение преобразования Лапласа сводит данную задачу к задаче теории упругости с коэффициентами, зависящими от параметра. Получено решение поставленной задачи в аналитическом виде для нескольких типов граничных условий на поверхности шахты.
слоистый массив
усредненная модель
гомогенизация
трансверсально-изотропная среда
ползучесть материала
экспоненциальное ядро ползучести
1. Олейник О.А. Математические задачи теории сильно неоднородных сред / О.А. Олейник, Г.А. Иосифьян, А.С. Шамаев. – М.: Изд-во МГУ, 1990. – 311 с.
2. Ильюшин А.А. Основы математической теории термовязко-упругости / А.А. Ильюшин, Б.Е. Победря. – М.: Наука, 1975. – 281 с.
3. Андреев В.И. Механика неоднородных тел / В.И. Андреев. – М.: Юрайт, 2015. – 255 с.
4. Шамаев А.С. О спектре одномерных колебаний в периодической комбинированной слоистой среде / А.С. Шамаев, В.В. Шумилова // Вестник Нижегородского университета им. Н.И. Лобачевского. – 2011. – № 4 (4). – С. 1882–1883.
5. Шамаев А.С. Прохождение плоской звуковой волны через слоистый композит с компонентами из упругого и вязкоупругого материалов / А.С. Шамаев, В.В. Шумилова // Акустический журнал. – 2015. – Т. 61, № 1. – С. 10–20.
6. Савенкова М.И. Применение метода осреднения в задаче упругопластического изгиба пластины / М.И. Савенкова, С.В. Шешенин, И.М. Закалюкина // Вестник МГСУ. – 2012. – № 9. – С. 156–164.
7. Савенкова М.И. Сравнение результатов конечно-элементного анализа с результатами асимптотического метода осреднения в задаче упругопластического изгиба пластины / М.И. Савенкова, С.В. Шешенин, И.М. Закалюкина // Вестник МГСУ. – 2013. – № 8. – С. 42–50.
8. Шамаев А.С. Асимптотическое поведение спектра одномерных колебаний в среде из слоев упругого материала и вязкоупругого материала Кельвина – Фойгта / А.С. Шамаев, В.В. Шумилова // Современные проблемы механики: сборник статей. Труды МИАН. – М.: МАИК «Наука/Интерпериодика», 2016. – Т. 295. – С. 218–228.
9. Шамаев А.С. Усреднение уравнений состояния для гетерогенной среды, состоящей из слоев двух ползучих материалов / А.С. Шамаев, В.В. Шумилова // Современные проблемы механики: сборник статей. Труды МИАН. – М.: МАИК «Наука/Интерпериодика», 2016. – Т. 295. – С. 229–240.
10. Bobyleva T.N. An Efficient Algorithm for Calculating Rheological Parameters of Layered Soil Media Composed from Elastic-Creeping Materials / T.N. Bobyleva, A.S. Shamaev // Soil Mechanics and Foundation Engineering. – 2017. – V. 54(4). – P. 224–230.
11. Лехницкий С.Г. Теория упругости анизотропного тела / С.Г. Лехницкий. – М.: Наука, 1977. – 416 с.
12. Bobyleva T. Approximate Method of Calculating Stresses in Layered Array / T. Bobyleva // Procedia Engineering. – 2016. – № 153. – P. 103–106.

Горные массивы в большинстве своем являются неоднородными по составу. Как показывает практика, свойства таких массивов могут значительно отличаться от свойств составляющих их компонент. Задачи о напряженно-деформированном состоянии неоднородной среды содержат дифференциальные уравнения с быстро меняющимися коэффициентами, которые характеризуют свойства отдельных компонент материала. Это серьезно усложняет решение данных краевых задач даже при использовании современных компьютеров, так как требует больших затрат машинного времени. Поэтому создаются модели, содержащие более простые уравнения с усредненными (эффективными) коэффициентами. Общая теория усреднения и примеры решения задач с применением этой теории даны, например, в работе [1]. Кроме того, подобные среды обладают способностью с течением времени медленно деформироваться при постоянных напряжениях, то есть обладают свойством ползучести [2, 3]. С помощью метода асимптотического усреднения в [4] построена математическая модель малых перемещений для комбинированной среды, состоящей из взаимно чередующихся слоев вязкоупругого материала и вязкой сжимаемой жидкости. Задача о прохождении плоской звуковой волны через композит конечной толщины со слоями упругого и вязкоупругого изотропных материалов решается в [5] также с использованием усредненной модели.

Работы [6, 7] иллюстрируют применение метода усреднения в задачах упругопластического изгиба пластин. Для гетерогенных сред, состоящих из упругого и вязкоупругого материалов [8], а также из двух ползучих материалов [9], построены соответствующие усредненные модели, описывающие совместное движение слоев. В работе [10] получены эффективные модули упругости слоистой упруго-ползучей среды.

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

Методы решения. Использованы механические модели теории упругости и теории ползучести. В качестве основных реологических соотношений принимаются соотношения теории линейной наследственности Больцмана – Вольтерра

bob01.wmf,

где ε – деформация, t – время, E – модуль Юнга, σ – напряжение, K(t) – ядро ползучести, которое является убывающей функцией времени [2]. Параметры реологической модели определяются на основе лабораторных испытаний. Применение преобразования Лапласа сводит задачу для упруго-ползучего массива к задаче теории упругости с коэффициентами, зависящими от параметра. Кроме того, использовано имеющееся в литературе решение о напряженно-деформированном состоянии трансверсально-изотропного массива с вертикальной цилиндрической шахтой кругового поперечного сечения.

В статье рассматривается слоистый горный массив, ограниченный горизонтальной плоскостью. Изотропные слои, из которых он состоит, взаимно чередуются и расположены также горизонтально. В данном массиве имеется вертикальная полость в виде цилиндра кругового поперечного сечения радиуса R, идущая от граничной плоскости. Требуется определить напряжения в данном массиве в случае, когда объемная сила – это его вес.

Выберем цилиндрическую систему координат bob02.wmf с началом в центре верхнего кругового сечения полости и осью z, направленной вертикально вниз (рисунок).

bob1.tif

Сечение полуплоскостью Orz слоистого горного массива, содержащего вертикальную цилиндрическую полость кругового поперечного сечения радиуса R (ε – относительный период ячейки периодичности)

Массив, состоящий из чисто упругих слоев

Уравнения равновесия (1) и уравнения состояния (2) для каждого слоя, составляющего массив, имеют вид

bob03.wmf

bob04.wmf (1)

bob05.wmf

bob06.wmf

bob07.wmf

bob08.wmf (2)

bob09.wmf bob10.wmf bob11.wmf,

где bob12.wmf компоненты напряжений, bob13.wmf – компоненты деформаций сплошной среды, λ, μ – постоянные Ламе, γ – величина удельного веса породы. Нагрузки на поверхности полости отсутствуют, то есть σr = 0 и σrz = 0 при r = R. На граничной поверхности z = 0 выполняются условия: σz = 0 и σrz = 0.

На горизонтальных граничных плоскостях слоев заданы условия идеального контакта: непрерывны все компоненты перемещений bob14.wmf и нормальная компонента напряжений σz, параллельная оси oz, то есть bob15.wmf и [σz] = 0.

Метод усреднения неоднородной упругой среды базируется на построении асимптотического решения по отношению к периоду составной среды [1]. Все модули упругости и удельный вес являются периодическими функциями координаты bob16.wmf (ε – относительный период ячейки) и кусочно-постоянными функциями этой переменной:

bob17.wmf bob18.wmf (3)

После применения метода усреднения получаем однородный трансверсально-изотропный массив. Уравнения состояния такой среды можно записать в следующем виде [11]:

bob19.wmf bob20.wmf bob21.wmf

bob22.wmf bob23.wmf bob24.wmf, (4)

где bob25.wmf – упругие податливости, являющиеся постоянными коэффициентами.

В работе [11] решена задача о распределении напряжений в массиве из трансверсально-изотропного материала с вертикальной цилиндрической полостью, имеющего плоскость изотропии, перпендикулярную оси z, получены следующие выражения для компонент напряжений (при условии, что uθ = 0 и, следовательно, τrθ = 0, τθz = 0):

bob26.wmf bob27.wmf bob28.wmf bob29.wmf. (5)

Подставим в (5) усредненные значения упругих податливостей а11, а12, а13, а33, а44, [12]. В результате получим

bob31.wmf

bob32.wmf (6)

bob33.wmf bob34.wmf

Напряжение σθ в радиальных сечениях у поверхности шахты:

bob35.wmf

Это напряжение будет в два раза больше, чем в сплошном массиве без полости на таком же расстоянии z от горизонтальной плоскости, ограничивающей массив.

Если перемещения на поверхности шахты в радиальном направлении невозможны (например, когда боковая поверхность шахты закреплена трубой), то граничные условия на ней будут ur = 0 и σrz = 0 при r = R. В этом случае радиальное напряжение, оказываемое массивом на стенку трубы:

bob36.wmf

Если на поверхности шахты действует давление qz, меняющееся пропорционально расстоянию [11], то к напряжениям σr и σθ в формулах (6) добавятся соответственно

bob37.wmf и bob38.wmf.

Массив, состоящий из упруго-ползучих слоев

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

bob39.wmf, (i, j = 1, 2, 3). (7)

Здесь символ * обозначает оператор свертки:

bob40.wmf. (8)

То, что ядра релаксации зависят от разности t – τ, как доказал Вольтерра, вытекает из условия инвариантности величин компонент тензора напряжений σij относительно изменений начала отсчета времени t.

Так как исходные слои упруго-ползучих материалов предполагаются изотропными, тензоры bob41.wmf и bob42.wmf (i, j, k, h = 1,2,3) в (7) соответственно имеют вид [2]:

bob43.wmf (9)

bob44.wmf. (10)

Здесь λ, μ – постоянные Ламе, δij – символ Кронекера, Ds(t) и Dν(t) – регулярные части ядер сдвиговой и объемной релаксаций соответственно. Предположим, что Dν(t) = kDs(t), где k = 0 или k = const > 1/3. Далее Ds(t) обозначим через D.

Пусть для каждого слоя ядро ползучести имеет вид: bob45.wmf где αi, di – постоянные, bob46.wmf, i – номер слоя, t – переменная, задающая время.

К уравнениям (7) с учетом (9), (10) применим преобразование Лапласа по времени:

bob47.wmf.

В результате получим задачу теории упругости, в уравнениях которой все коэффициенты зависят от параметра p. К такой задаче применим метод усреднения [1]. Используя формулы (5) для компонент напряжений трансверсально-изотропного тела и выражения для усредненных (эффективных) модулей упругости, полученные в работе [12], запишем решение поставленной задачи в изображениях Лапласа:

bob48.wmf;

bob49.wmf; (11)

bob50.wmf

После применения к (11) обратного преобразования Лапласа формулы для компонент напряжений будут иметь вид

bob51.wmf

bob52.wmf (12)

bob53.wmf

bob54.wmf (13)

bob55.wmf bob56.wmf

где δ(t) – дельта-функция Дирака. В этих формулах введены следующие обозначения

(i – номер слоя, i = 1,2):

bob57.wmf bob58.wmf

ki – коэффициенты пропорциональности ядер объемной релаксации ядрам сдвиговой для каждого слоя. Кроме того, p1, p2 – корни следующего квадратного уравнения:

bob59.wmf, (14)

имеющего положительный дискриминант

bob60.wmf.

Отсюда следует, что p1, p2 – действительные числа.

Напряжение σθ рядом с поверхностью круговой полости (r = R) для любых углов θ будет

bob61.wmf

bob62.wmf

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

bob63.wmf

bob64.wmf

Если на поверхности шахты действует давление qz, меняющееся пропорционально расстоянию z (например, внутри находится жидкость [11]), то к напряжению σr в формуле (12) добавится bob65.wmf, а к напряжению σθ в формуле (13) bob66.wmf.

Выводы

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

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


Библиографическая ссылка

Бобылева Т.Н., Голобоков О.С. ЗАДАЧА О ВЫЧИСЛЕНИИ НАПРЯЖЕНИЙ В СЛОИСТОМ МАССИВЕ С ВЕРТИКАЛЬНОЙ ЦИЛИНДРИЧЕСКОЙ ШАХТОЙ // Современные наукоемкие технологии. – 2018. – № 6. – С. 34-39;
URL: https://top-technologies.ru/ru/article/view?id=37029 (дата обращения: 28.11.2021).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074