Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Механика деформирования и разрушения горных пород

..pdf
Скачиваний:
4
Добавлен:
12.11.2023
Размер:
8.26 Mб
Скачать

Рис. 3.27. Численная схема расчета цели­

У

Р

,<п

 

 

fp

ка

0—0

.

 

0—0

 

 

 

 

 

 

0—0

 

 

 

 

 

 

0—0

 

 

 

 

 

 

о—о

 

 

 

 

с

 

0—0

 

 

 

 

 

0—0

 

 

 

 

 

 

0—0

 

 

 

 

 

 

0—0

 

 

 

 

 

 

0—0

 

 

 

 

I X

 

0 П

П

П

И

П !

 

,

 

 

В

 

 

вие системы целик—кровля—почва. При этом, если жесткость вмещающих пород массива Nn на контакте целик -порода будет намного больше жесткости N целика, то имеет место работа целика в режиме заданных деформаций.

Таким образом, при исследовании упруго-пластичного со­ стояния целиков нужно различать [40] два основных режима их работы: режим заданных нагрузок и режим заданных деформаций.

Заметим, что значения жесткости N и Nn могут быть выражены через деформационные характеристики и геомет­ рические параметры целика и массива пород.

Рассмотрим ‘ленточный целик шириной 2b и высотой 2Л, находящийся в режиме заданной нагрузки. Будем полагать, что горизонтальные перемещения на контакте целик—порода равны нулю (w —0, при г = ±Л), боковые стороны его свободны от напряжений (а„ = тЛУ= 0 при x = ± h ) , а нагрузка, прило­ женная к целику, оу = Р при y =± h . Благодаря геометричес­ кой и силовой симметрии рассмотрим верхнюю правую четверть целика (рис. 3.27). Заметим, что первое допуще­ ние не является единственно возможным, поскольку при деформировании целика в реальных условиях горизонтальные перемещения контактных поверхностей целика изменяются от нуля при полном его сцеплении с массивом пород до максимально возможных при отсутствии трения на контакте целик—порода.

Для решения задачи построим прямоугольную сетку, ко­ ординаты узлов которой: Xi = ib/m; =jh/n (/ = 0, 1......ш; у= 0,

1 , п ) .

Проведем дискретизацию функционала (3.69) конечными разностями и обозначим через г,7, ии вертикальные и горизон­ тальные перемещения в (/, /)-ом узче 1 рапичпые юпня в узлах:

"о ../•■= о : г,-

i ' i м= 0-

(3*75)

где / —0, I. т: у 0. 1. /?.

Тогда полхчим задачу нелинейного программирования по нахождению такого набора величин nij4 vir которые доставляют

I 11

Рис. 3.28. Развитие зон деформационного разупрочнения целика в зависимости

от степени

его

нагружения:

У— упругая

зона;

2

зона деформационного разупрочнения на ниспадающем участке;

3 — зона остаточной

прочности

минимум функционалу (3.69), где полная диаграмма дефор­ мирования пород определяется связями (3.72).

Численный эксперимент нагружения и деформирования целика проводится при следующих значениях исходных па­ раметров:

М= —0,26; v= 0,3;

£ /Г п=Ю 3;

Го = 0,8Г„;

 

п = т= 16;

b = h.

 

Развитие зон

запредельного

деформирования в целике

в зависимости

от степени

нагружения

Р/Тп показано на

рис. 3.28. Анализ результатов численного эксперимента показал, что при небольших значениях нагрузки Р/Тп целик полностью находится в упругом состоянии. При достижении нагрузкой некоторого значения Рг/Тп в центре целика появляется область разупрочнения, в которой материал деформируется по нис­ падающей ветви зависимости интенсивность касательных на­ пряжений Т — интенсивность деформаций сдвига Г

Эта область формировалась в целике при изменении

нагрузки в

диапазоне

1,65Гп< / >,< 1,75ГП. При дальнейшем

увеличении

нагрузки до

значения Р=1,8ГП часть *материала

целика вышла на .горизонтальный участок деформирования диаграммы Т— Г, 'т. е. участок остаточной прочности. Даль­ нейшее увеличение нагрузки привело к значительному увеличе­ нию площадей зон деформационного разупрочнения, при этом зона остаточной прочности оказалась окаймленной зоной деформирования на ниспадающем участке.

Следует отметить, что даже при значительных значениях Р зона деформационного разупрочнения не охватила полностью вертикальное сечение целика. При достижении некоторого критического значения Pt, в данном случае Ркл2,\5Тп, ите-

Рис. 3.29. Изменение

перемещений в

целике упругой (/) и рдзупрочняющей

(2)

моделей:

 

 

 

а)

Р/Т„ от г//г К)3: 6

и/Л-Ю'3 от y/h:

v/lr 103

x/h

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

Изменение наибольшего вертикального перемещения от­ носительно середины плоскости целика в зависимости от заданных нагрузок показано на рис. 3.29, а. При нагрузке Р<РЪ перемещения увеличиваются по линейному закону, поскольку материал целика деформируется упруго. С даль­ нейшим увеличением нагрузки перемещения возрастают, при­ ближаясь к горизонтальной асимптоте. Эпюры распределения горизонтальных для х = Ь и вертикальных для y = h перемеще­ ний для упругой (/) и разупрочняющейся (2) моделей целика при нагрузке Р=1,85ГП показаны на рис. 3.29, б, в.

Анализ данных расчета показывает, что максимальное горизонтальное перемещение для упрочняющейся модели в 3 раза больше соответствующего перемещения для упругой модели. Сравнение эпюр вертикальных перемещений показы­ вает, что разупрочнение материала привело к их увеличению

в1,91 раза в центре целика, которое снижается с приближением

ккраю целика. Этот результат объясняется формированием значительной зоны разупрочнения в центре целика и наличием упругой зоны вблизи края.

Влияние отношения TQ/Tn на развитие зон разупрочнения при фиксированной нагрузке Р=1,85ГП показано на рис. 3.30. Численный эксперимент проводился при следующих прочност­ ных и деформационных характеристиках: М = 0,2G; v= 0,3;

Е/Тп= 103 Геометрические размеры целика составляли б = /?, а относительные значения остаточной прочности Т0/Т при­ нимались равными 0,9; 0,75; 0,6.

Анализ результатов численного моделирования показывает,

что конфигурация

и

размеры зон

разупрочнения изменяются

с изменением Т0.

С

уменьшением

Т0 зона деформирования

на ниспадающем участке диаграммы «Т— Г» становится более

0,9

0,6

Рис. 3.30. Развитие зон упругого деформирования, разупрочнения на нис­ падающем участке и остаточной прочности в зависимости от Т0, Т„

протяженной и узкой, при этом она окаймляет зону дефор­ мирования на участке остаточной прочности.

Деформирование однородных целиков в режиме заданных деформаций. Режим заданной нагрузки трудно реализовать экспериментально при испытании образцов, поскольку кон­ тактная поверхность образца через жесткие плиты испытывает одинаковое вертикальное перемещение. Данный режим дефор­ мирования, имеющий место при испытаниях образцов, назовем [7, 40] режимом заданных деформаций.

Постановка задачи. Будем полагать, что целик состоит из двух слоев различных по прочности материалов, при этом более прочный слой породы (аналог плиты нагружающего устройства) расположен по краям целика (рис. 3.31). Оба материала способны деформироваться за пределом прочности и имеют аналогичную диаграмму деформирования Т — Г Однако для более прочного слоя значения Т и Т0 будут значительно больше.

Пф описанной выше схеме построим с узлами (/, у) равномерную сетку в целике. Значение модуля объемного сжатия К и функцию g (Г) будем выбирать в зависимости от того, к какому слою принадлежит /, у-й узел. К ограничениям

Рис. 3.31 Расчетная схема:

/ прочный слой: 2 целик

Рис. 3.32. Развитие зон упругого деформирования, разупрочнения на нис­ падающем участке и остаточной прочности в целике в зависимости от нагружения (вертикального перемещения)

также следует добавить условие непрерывности вертикальных перемещений на контактной поверхности целика.

 

 

 

Vi,„ = v,

 

(3.76)

где V— вертикальное

перемещение верхнего

слоя

на контакте

с целиком; (/ = 0; 1;

; т).

нагружения целика

выполняется

Численный

эксперимент

в режиме заданных деформаций при ограничениях

(3.75), (3.76)

и следующих

исходных данных: М = 0,26;

v= 0,3; Е/Тп= 103;

Го = 0,8Гп; т = п= 16;

b = h\

T n/Tnl = E J E = 102; Л1/Л= 0,25.

Здесь параметры

£ \, hi относятся

к жесткому слою,

модуль упругости и предел прочности которого принимаются на два порядка больше.

Для удобства дальнейшего описания процесса деформирова­ ния введем безразмерное перемещение V = \ 0 0 v / h .

Развитие зон деформационного разупрочнения в зависимо­ сти от вертикальных перемещений V показано на рис. 3.32. Принципиальное отличие этого режима от режима заданной нагрузки проявляется в образовании и развитии зон разуп­ рочнения. В режиме заданных деформаций образование зон разупрочнения начинается в углу целика и с ростом нагружения развивается под углом 45е к серединной линии целика. При дальнейшем увеличении V зона разупрочнения образовывается в центре целика и в последующем сливается с зоной, развивающейся от угла целика. В дальнейшем развитие зон приводит к значительному увеличению их площадей, особенно зоны остаточной прочности. •

Полученная картина образования зон деформационного разупрочнения в качественном отношении согласуется с на­ блюдаемыми в шахтных условиях конусами разрушения в кра­ евой части целика.

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

Зависимость среднего давления Р/Тп на целик от вертикаль­ ного смещения V контактных поверхностей показана на рис. 3.33, а, из которого видно, что она качественно повторяет диаграмму Т — Г для образца горной породы. При этом наибольшая несущая способность целика Рп/Тп наблюдается на пределе прочности. При дальнейшем увеличении V происходит разупрочнение целика и снижение его несущей способности до Р0/Т п при выходе на горизонтальный участок, соответствующий

остаточной прочности. Далее, с ростом смещения

несущая

способность целика не изменяется: Р= Р0 = const.

 

Следует отметить, что отношения между значениями ин­

тенсивности касательных напряжений Т на пределах

прочности

иостаточной прочности для материала целика (горной породы)

идля всего цедика не сохраняются. Так, в данном численном

эксперименте для горной

породы при отличии

Тп и Т0 на

20% (Го = 0,8Гп) отличие

между предельной Рп и

остаточной

Р0 несущей способностью целика составило 10%, т. е. оказалось в 2 раза меньше. Это объясняется влиянием упругих зон в целике, а также наличием касательных напряжений по контакту целик—порода, поскольку при эксперименте рас­ сматривался случай их полного сцепления.

Распределение вертикальных напряжений вдоль горизон­ тального сечения y = h при Г=0,6ГПв зависимости от смещения V приведено на рис. 3.33, 6. Анализ эпюр напряжений Gy показывает, что снижение несущей способности целика при

изменении

V от 0,1 до 0,15 было вызвано его разупрочнением

в

краевой

части

целика.

При

изменении

V от 0,15 до 0,2

происходило

снижение

прочности

по всей длине

целика.

С

дальнейшим ростом

V до 0,3 происходило некоторое

перераспределение

напряжений,

однако,

среднее

давление

 

1 ь

^

на

целик

практически не

отличалось

от его

Р= - Jxg(.Y)d.v

 

" о

 

V—0,2. На

диаграмме

Р/Тп— V этому

режиму

значения при

деформирования целика соответствует горизонтальный участок. Следует отметить, что необходимое условие (3.71) сущест­ вования экстремума функционала полной энергии системы при всех значениях V выполняется, т. е. целик обладает устой­

чивостью в физическом смысле.

Выполненный выше анализ поведения целика в режиме заданных деформаций и полученные при этом результаты могут быть также использованы для установления закономер­ ностей деформирования и разрушения образцов горной породы при испытаниях их на «жестких» прессах. При этом удается установить с ростом V кинетику образования и развития в образце зон разупрочнения на ниспадающем участке диаграм­ мы Т — Г и на участке остаточной прочности.

éy/Tn

а

 

1

0,1

0,2

0,3

V

 

О

0,25

0,5

0,75 х/6

 

О

 

Рис. 3.33. Закономерности

распределения

средних

нагрузок

на целик (а)

и

вертикальных

напряжений

в целике

при

y=h(ô):

 

 

/ -

Г= 0,1;

2 — Г=0,125: 3 -

I = 0.150; 4 -

У=()Л15: 5 - 1 = 0 ,2 0 0

 

 

Отметим, что закономерности деформирования образцов породы (целиков), полученные с помощью численных экс­ периментов для режимов заданных нагрузок и заданных деформаций, согласуются с результатами лабораторных ис­ следований. Так, при лабораторных испытаниях образцов с использованием щеточных плит наблюдается формирование начала разрушения в центре образца, а при использовании традиционных монолитных плит, с помощью которых удается реализовать режим заданных деформаций, разрушение начина­ ется с опорных плоскостей. Полученные при этом конусы или клинья разрушения при обычном испытательном обору­ довании согласуются с результатами (см. рис. 3.33) численных экспериментов.

Деформирование неоднородных ленточных целиков. Разработ­ ке методов расчета несущей способности неоднородных целиков посвящен ряд работ [40], которые носят экспериментальный характер. В этих работах исследуется влияние толщины, взаимного расположения и несущей способности каждого из слоев. Анализ и теоретическое обобщение закономерностей деформирования неоднородных целиков выполнено в работах А. К. Черникова [40]. При этом в качестве основной исполь­ зована модель жесткопластического тела.

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

Соотношения между прочностными и деформационными характеристиками пород целиков при численных экспериментах

а

6

в

Рис. 3.34. Схемы строения не­

 

 

 

однородных слоистых целиков:

1— слабый слой; 2 — материал цели­ ка

ФФФФФAФФФФФФ

были приняты

следующими: Е = Е 1;

7^ = 0,75 Г;

Го = 0,8Г;

Го=0,8Гп,

где

7^,

— пределы

прочности

и модуль

упругости

слоя.

 

 

 

 

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

Последовательное развитие зон деформационного разуп­ рочнения в 1/4 области целика в зависимости от нагрузки Р/Тп можно проследить на рис. 3.35, которые соответствуют моделям строения целика на рис. 3.34. Отметим, что харак­ терной чертой для всех моделей является локализация зон деформационного разупрочнения в слабом слое, но происходит это при разных значениях сжимающей нагрузки. Так, для модели (а) разупрочнение началось при нагрузке, равной 7>=1,2Гп,г до этого целик находился в упругом состоянии.

При росте нагрузки до значения Р=1,5ГП в слабом слое впервые появляется зона остаточной прочности, т. е. материал деформируется на горизонтальном участке диаграммы Т — Г С увеличением нагрузки до 1,7Т в слабом слое происходит рост зоны остаточной прочности, при этом в более прочном основном материале целика начинается развитие зоны разуп­

рочнения на ниспадающем участке зависимости Т

— Г

При дальнейшем увеличении нагрузки рост зон

разупроч­

нения в основном материале целика продолжается, направление их развития происходит от уже разрушенного и вышедшего на остаточную прочность материала ,слабого слоя к углам целика. При значениях нагрузки Р/Тп^ 2 итеративный процесс минимизации функционала завершить не удается. Это соот­ ветствует физическому разрушению целика, а математически — нарушению необходимого условия (3.71) существования ми­ нимума потенциальной энергии, которое без труда переносится на многослойную среду.

Для модели, представленной на рис. 3.35, <5, до значения нагрузки Р=1,5Тп зона разупрочнения локализуется в слабом слое, однако при Р=\,6Тп в центре целика уже образовывается достаточно обширная зона разупрочнения в более прочном

материале. При P = \ J T n

в слабом слое впервые возникает

зона деформирования на

участке остаточной прочности зави-

Рис. 3.35. Раш ите юн упругого деформирования (У). разупрочнения на ниспадающем участке (_) и ocrai очной прочности {3) в 1,4 части неоднородных

педиков

в

зависимости от

нагрузки

Р Т„:

</, о,

в

схемы

строения педиков, соответствующие рис. 3.34

симости

Т — Г

При

нагрузке

Р= 1,8; 1,9Тп рост зон разуп­

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

Рис. 3.36. Зависимость средней нагрузки на целик PjTn от наибольшего вертикаль­ ного перемещения г по контакту целик — порода:

I однородная модель среды; 2, 3. 4 ветсгвенно схемы строения целика в.

Для модели (см. рис. 3.35, в), в которой слабые слои располагаются на контактных поверхностях, разупрочнение начинается при более высоком; по сравнению с другими моделями, значении нагрузки Р=1,5ГПдо значения Р=1,65ГП, разупрочнение локализуется в углу слабого слоя. При нагрузке Р=\,1Тп зона разупрочнения «прорастает» от угла к центру целика. Увеличение нагрузки приводит к интенсивному раз­ витию зоны остаточной прочности.

Для двух последних моделей значения нагрузки, при ко­ торой итеративный процесс минимизации завершить не уда­

ется, оказались одинаковыми и

приблизительно равными

Р= 1,95ГП.

экспериментов показывает,

Анализ результатов численных

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

Проанализируем влияние слоев на зависимость между нагрузкой на целик Р/Тп и наибольшим вертикальным пере­ мещением V по контакту целик — порода. Соответствующие зависимости для различного расположения слоев, а также однородного целика показаны на рис. 3.36.

Анализ численных эксперйментов показывает, что до не­ которого значения нагрузки Рп, своего для каждой модели, между Р и г имеет место линейная пропорциональная зави­ симость. Это объясняется тем, что большая часть материала целика находится в упругом состоянии. При дальнейшем увеличении нагрузки Р>РП зависимость сильно зависит от нагрузки. В первом приближении на этом участке зависимость можно аппроксимировать линейной функцией с небольшим углом наклона графика к горизонтальной оси. При этом небольшое увеличение нагрузки приводит к существенному

росту вертикальных перемещений

и проявляется различие

в деформировании неоднородных

целиков.