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

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

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

Компоненты распределены в теле равномерно. Отличие данной структуры от предыдущей заключается в том, что здесь элементы 2, содержание которых в 2 раза выше, имеют непосредственный контакт друг с другом, в то время как элементы 1 друг от друга разрознены и непосредственного контакта не имеют.

Число возможных плоскостей со по элементам 2 в 2 ра­ за больше, чем по элементам 1. При этом угол а ориенти­ ровки плоскостей по обоим элементам одинаков, что озна­ чает равенство коэффициентов поперечных остаточных дефор­ маций.

a

По элементам 2 материал одинаков в направлениях стх и а 2,

по элементам

1

имеется

анизотропия. По элементам

1

в направлении

ст2

плоскости

со могут проходить только

под углом 45°, что в условиях, близких к одноосному сжатию, в большинстве случаев не может реализоваться. В условиях, близких к одноосному сжатию, сдвиг по элементам 1 в этом случае может произойти по плоскостям, угол наклона которых относительно направления а 2 близок к 90°, коэффициент поперечной деформации здесь будет велик, а пластичность

мала. Таким образом, хрупкость

материала

может

меняться

в зависимости от ориентировки

плоскостей

сдвига

со.

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

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

Картина изменяется, если размер структурных элементов компонент будет разный.

Пример. Представим среду, состоящую из двух компонент с одинаковым содержанием их (по массе) при одинаковой плотности. Линейный размер минеральных зерен первой компоненты в 10 раз меньше, чем второй,

например, Я!=0,1 мм

и а2 = 1 мм.

площадью 1 см2, расположенной

Для образования

одной плоскости со

под углом 45°, требуется число элементов

/? = (1/а)2. Для

компонентжбудем

соответственно иметь:

/i' = (10/0,1)2 = 104; л" = (10/1)2= 10.

В единице объема

полное число N структурных элементов компонент:

Ni =(10/0,1)3= 106;

уу2=(10/1)3=103

со в единице объема:

 

Число плоскостей

 

/ii =Ni jn = 102; n2 = N2/n "= 10.

Плоскость, расположенная под углом 45 , как это следует из модели, полностью «выложена» из элементов без зазоров. Наиболее вероятное расположение плоскостей при равномерном заполнении под углом 30 В этом случае для построения плоскости œ требуется примерно в 2 раза меньше элементов, т. е. л'%0.5 104; 0,5 10“

Число плоскостей (о в этом случае также удвоится:

ni = 10е/(0,5 • 104)= 2 • 102;

/,2= 103/(0,5 • 102)= 20.

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

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

Рассмотренные выше результаты по влиянию скорости деформирования на прочность и предел упругости при раз­ личных видах напряженного состояния поддаются количествен­ ному описанию с помощью кинетического уравнения. Это относится как к результатам, полученным при заданной скорости деформирования, так и к результатам, полученным при ползучести [32].

Однако наряду с этим необходимо отметить и ограничения в непосредственной применимости кинетической концепции. Главная причина этого обстоятельства заключается в неоп­ ределенности применения к твердым телам понятия грамммоля, а отсюда и неопределенности понятия входящих в урав­ нение величины £/0, Y и è0. Это вызвано тем, что неизвестна доля (в процентах), участвующих в деформировании и раз­ рушении элементарных частиц, заключенных в единице объема исследуемого твердого тела. Кроме того, в процессе дефор­ мации и разрушения твердого тела относительные перемещения частиц имеют различный характер: полный разрыв связей, перемещения типа замещения и внедрения в пределах кристал­ лической решетки, диффузия и т. д.

Каждое из этих перемещений сопровождается преодолением характерного для этого случая потенциального барьера, пред­ ставляющего энергию активации данного «процесса. Количест­ венно упомянутые барьеры могут различаться в несколько раз, что приводит к дополнительной неопределеннее™ Уп0~ мянутых выше величин.

В координатах т и \gt{ по чисто формальным признакам кинетическое уравнение может описывать только те случаи, когда напряжение возрастает с увеличением скорости дефор­

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

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

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

Рассмотрим смысл входящих в кинетическое уравнение величин.

-—эн ер ги я актив ац ии, характеризующая энергию свя­ зи элементарных частиц в твердом теле. Для каждого тела энергия связи имеет свое значение и, таким образом, является константой. Параметр у— с т р у к т у р н о ч у в с т в и т е л ь н ы й ко эф фициент, который интерпретируется в кинетической теории как произведение активационного объема на коэф­ фициент локального перенапряжения. Значение у зависит от термической обработки, режима отжига и размера зерен компонент. С ростом размера зерен у возрастает. Под активационным объемом подразумевается часть объема кри­ сталлической решетки размером 10"23 см, в которой проис­ ходит процесс ее перестройки и разрушения, вызванный действием приложенного внешнего напряжения т, усиленного коэффициентом локального перенапряжения.

Величина

({У—ут) называется эф ф е к т и в н о й

энергией

а к т и в а ц и и

и обозначается С/Эф. Прикладываемое

напряжение

т уменьшает энергию £/эф и тем самым снижает потенциальный

барьер, препятствующий атому покинуть свое положение

равновесия Чем

меньше

{Уэф, тем

меньше атому

требуется

времени ожидания для

выхода

из

потенциальной

ямы.

В уравнении

для

долговечности

величина

 

t0

означает

пе ри о д с о б с т в е н н ы х к о л е б а н и й

а т о м а около положе­

ния равновесия. Еще эту величину,

равную 10“ 12 с,

называют

ч а с т о т н ы м

мн о ж и т е л е м .

В

уравнении

для

скорости

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

[32] величина

е0

представляет

п р е д е л ь н о

в о з м о ж н у ю

с к о р о с т ь

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

в

твердом

теле и также играет роль частотного множителя. Об этом более подробно будет сказано ниже.

Несколько ранее работ С. Н. Журкова и его школы были начаты исследования пластических деформаций металлов при разных значениях скорости деформирования, включая рас­ пространение в них волн пластичности под руководством H. Н. Давиденкова (1950 г.). В этих исследованиях развивалась кинетическая концепция пластической деформации. Уравнение, связывающее время релаксации / с приложенным внешним напряжением а, имеет вид

t = t0e\p{[U0—U1th(ya)]/(AT)},

где U0— энергия активации, имеющая тот же смысл, что и в уравнении С. Н. Журкова; Ux— некоторая постоянная, харак­ теризующая коллективность процесса диффузии; tQ— период со­ бственных колебаний атома около положения равновесия; у — структурно чувствительный коэффициент, имеющий тот же смысл,

что и в формуле С. Н. Журкова; К — постоянная

Больцмана;

Т — абсолютная температура; КТ= 2463,7 Дж/моль

при 20° С.

Как видно, приведенное уравнение отличается от уравнения

С. Н. Журкова для долговечности наличием гиперболического тангенса при произведении уа. По физическому смыслу рассматриваемые уравнения очень близки. В обоих случаях имеются структурно нечувствительные постоянные £/0, умень­ шающиеся под действием внешнего приложенного напряжения. Разница между U0 и работой внешних напряжений, названная эффективной энергией активации £/эф, преодолевается атомом за счет энергии тепловых флуктуаций, что и приводит в одном случае к разрушению тела, в другом— к образованию в теле остаточной деформации. Гиперболический тангенс введен с це­ лью наиболее точно отразить форму потенциальной кривой, описывающей энергию связи между двумя атомами.

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

Исследуя зависимость скорости движения дислокации от температуры, авторы доказали справедливость следующего

уравнения:

 

»=/(<*) ехр [ - и0/[КТ)].

(2.21)

Оказалось, что U0 от приложенного напряжения не зависит, а по значению близка к энергии активации диффузии от­ рицательного иона в решетке фтористого лития и равна 0,6 эВ.

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

игорные породы.

Вдальнейшем будет использовано кинетическое уравнение

С.Н. Журкова для скорости деформирования, поэтому подробнее рассмотрим входящую в это уравнение величину è0. В отличие от универсальной константы t0 в уравнении долговечности величина

е0 не имеет такого универсального смысла.

В названной выше работе Джойстона и Гилмана получена эмпирическая зависимость (2.21) скорости движения дислокации vs от напряжения. Это выражение можно заменить близким

ему по смыслу:

 

vs= v0exp [ —t/эф/ (КТ )],

(2.22)

где va— предельная скорость распространения дислокации. Выражение для скорости относительной деформации è имеет

вид (М. И. Бессонов, 1964 г.) è = bmvs или, учитывая выражение

(2,22)

è = bmv0exp [ - С/Эф/ (К Т )] .

 

При U3ф = 0

 

è = bmv0.

(2.23)

Выражение (2.23) позволяет подсчитать е0, физический

смысл которой — м а к с и м а л ь н о в о з м о ж н а я

с к о р о с т ь

п р о т е к а н и я п р о ц е с с а . Максимально возможная скорость передвижения дислокации v0 близка к скорости распространения поперечной волны и составляет около 105 см/с. Так как 80 отражает процесс разрушения твердого тела, то плотность дислокаций т должна достигать критической величины. В ка­

честве таковой

можно

принять т = 109 [32].

Подставляя это

выражение в

формулу

(2.23)

и

принимая вектор

Бюргерса

b = 10~8, получаем: è0 = 10“8

109

105 = 106c -1

другим путем,

Оценку величины s0 можно получить и

если использовать

эмпирическое

соотношение

[32]

 

 

 

 

/е$ = е&.

 

(2.24)

Постоянная

г'0

по порядку

величины составляет

10_3 (см.

табл. 2.1). Показатель степени к (см. табл. 2.1), как было установлено специальными расчетами, очень близок к единице. Подставляя значение 10 " 3 в выражение (2.24) и принимая во внимание / = 10 “ 12, получаем

é0= 10“3 1012 • 1= 109с -1,

т. e. большее значение, чем при первой оценке.

Причина такой разницы объясняется отсутствием точных

данных о предельной (критической) плотности дислокации

т в первой

оценке.

Величина

т в зависимости от различных факторов [32]

может принимать значения от 109 до 1012. Внесение последнего приблизило бы первый результат ко второму. Во второй оценке использование г=10~*2 с также достаточно произволь­ но. Более обоснованный выбор указанных величин может быть сделан только на основании непосредственных экс­ периментальных определений и, в частности, при введении в эксперимент фактора температуры. Не располагая такими данными, целесообразно провести расчеты для всего диапазона возможных значений в0 от 106 до 109с -1, что позволит оценить соответствующие вариации U0. Результаты таких расчетов приведены в табл. 2.1.

Как видно, изменение U0 при вариации è0 на три по­ рядка не превышает 10%, что объясняется структурой ки­ нетического уравнения. Для мрамора, галита и сильвина рассчитанные значения U0 сравниваются со значениями энергии активации диффузии Q (см. табл. 2.1), заимствованными из других источников, где они были получены экспери­ ментально совершенно другими методами. Так, Q для мрамора определялась методом меченных атомов (В. А. Вейбулл, 1939 г.) и оказалась близкой к энергии диффузии молекул СаСОэ в решетке кальцита. Как видно из табл. 2.2, рассчитанные значения U0 на основании механических испытаний и Q ока­ зались весьма близкими. Значения Q энергии активации в галите и сильвине заимствованы из справочного материала. Значение Q и в этом случае также оказалось весьма близким к рассчитанным по результатам механических ис­ пытаний значениям U0.

2.4. ОПРЕДЕЛЯЮЩИЕ УРАВНЕНИЯ МЕХАНИКИ РАЗРУШЕНИЯ, ОСНОВАННОЙ НА КИНЕТИЧЕСКОЙ И СТАТИСТИЧЕСКОЙ ТЕОРИЯХ ПРОЧНОСТИ

Рассмотрим несколько моделей предельных состояний в ко­ ординатах lge и т. Простейший вариант такой модели показан на рис. 2.20 для однокомпонентной среды с умеренной дис­ персией свойств соответствующей функции распределения 2 (см. рис. 2.1). Луч 1 соответствует прочности пород при боковом давлении: ст2 =pi =0, луч 2 — а 2=/?2, луч 3 — а 2=/?3, луч 4 — ст2 а- При этом р2<Ръ<Р4- Лучи, изображенные штри­ ховыми линиями, дают значения соответствующих пределов упругости. Все лучи пересекаются ô одном полюсе с абсциссой

 

U0 103,

Константа

Yn»

U " 103.,

YP ,

Q 103,

Yy.

 

 

 

 

 

 

 

 

 

 

 

к

Дж

Рп

Ру

Дж • см

Дж

Д ж •см 2

Д ж

Д ж •см 2

 

моль

моль•к

М О Л Ь

м о л ь •кг

м о л ь

м о л ь •кг

Мрамор

106

и

222,1

19

22,7

377,1

3.2-

10"

ю 7

1,1

221л

19

22,7

377,1

3.2-

10"

10е

1,1

232,5

19

22,7

377,1

3.2-

10"

ю 9

1,1

238,8

19

22,7

377,1

3.2-

10"

Диабаз

3 89 ,7

1 1 5 2 2

2 50

 

3 9 7 ,2

И 5 22

2 50

528

4 0 2 ,2

11 5 22

2 50

4 0 8 ,5

11 522

2 50

 

W6

1

204

8,4

12,9

105,6

4 1 0 "

2 5 2 ,6

3 16 3

 

130

 

130

107

1

209,5

8,4

12,9

105,6

4-10"

2 5 7 ,7

3 163

 

 

10е

1

215,8

8,4

12,9

105,6

4-10"

2 63 , Г

3163

130

ю 9

1

221,2

8,4

12,9

105,6

4-10"

2 6 9 ,4

3163

130

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

_

387,6

Галит

I

 

 

 

 

106

1,2

127,8

 

 

 

_

 

150,8

1022

 

 

 

 

 

150,8

1022

ю 7

1,2

133,2

387,6

 

 

 

 

_

 

150,8

1022

10е

1,2

139,1

387,6

 

 

 

 

 

 

150,8

1022

ю 9

1,2

145

387,6

 

 

Галит

II

 

 

 

 

 

 

_

358,2

 

 

 

 

106

1

127,8

 

 

 

_

 

 

8 46

ю 7

1

133,2

358,2'

 

 

8 46

358,2

 

 

_

_

8 46

ю 8

1

139,1

 

 

 

 

 

8 46

109

1

145

358,2

 

 

 

 

 

 

 

 

 

 

 

Сильвинит

 

 

 

 

106

0,1

99,7

 

 

240,9

 

 

_

 

125,7

4 27

 

 

 

 

 

125,7

4 2 7

ю 7

0,1

105,1

240,9

 

 

 

 

_

_

125,7

4 2 7 '

ю 8

0,1

111,0

240,9

 

 

ю 9

0,1

116,5

240,9

 

 

125,7

4 27

 

 

 

 

 

 

 

 

I/o Ю3,

Константа

Y,M

 

 

и „ ю 3,

 

е ю 3,

 

 

к

 

 

 

 

V

V

 

Дж

 

 

Дж • см 2

 

 

Дж

Дж • см 2

Дж

Дж • см 2

 

 

Р».

Ру

 

 

 

 

моль

моль • кг

 

 

моль

моль • кг

моль

моль • кг

 

 

 

 

 

 

 

 

 

 

 

Квар\цевый

песчаник

 

 

 

 

106

0,1

144,5

6

9,1

192,7

6,2 - 10 ~24

310

10014

2

248

107

0,1

149,6

6

9,1

192,7

6,3

Ю"24

318,4

10014

248

ю 8

0,1

155

6

9,1

192,7

6,3-10 ~24

324,7

10014

248

ю 9

0,1

161

6

9,1

192,7

6,3 10 24

329,3

10014

248

 

 

 

 

 

 

Песчаник

( ВО )

 

 

 

 

106

0,8

105,6

 

_

_

_

_

_

ю 7

0,8

155

 

_

_

10е

0,8

158

 

_

ю 9

0,8

166,8

 

Рис. 2.20. Зависимость пределов прочности от скорости деформирования е для однокомпонентной среды с умеренной дисперсией свойств структурных элемен­ тов

èt = 0. Такой пучок лучей описывается единым кинетическим уравнением (2.1) для скоростей деформирования.

Входящие в это уравнение величины были объяснены ранее. При описании разных лучей в пучке в уравнении меняется значение у, которое и определяет наклон луча. Произведение ут при какой-либо одной скорости ё сохраняется постоянным для всех лучей. Например, при ét величина ут в точках a, b, с, d, е, /, g, h имеет одинаковое значение как для пределов прочности, так и для пределов упругости независимо от бокового давления а 2. С увеличением а 2 величина ' у уменьшается. Вследствие постоянства произведения ут при заданной скорости с уменьше­ нием у соответственно увеличиваются предел прочности и пре­ дел упругости. По мере повышения скорости деформирования б1<ё2<ёз<ё4 и т. д. ут также растет за счет т.

Значения пределов прочности и упругости при какой-то одной скорости деформирования и при разных боковых давлениях ст2 описываются уравнениями предельных состояний (1.2) и (1.7).

С изменением скорости в этих уравнениях меняются т," , т", А, В и С, которые являются функциями скорости дефор­ мирования.

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

Величины т„(ё), т"(ё) получаются из кинетических ура­ внений:

Рис. 2.21. Паспорт прочности (сплош­ ная линия) и упругости (штриховая линия) для однокомпонентной среды при разных скоростях деформирова­ ния

Tn(è) = [ln(è/è<>) ATT"+ f/0] 1/Тп!

(2.25)

TÏ(è) = [\n(è/è0)KT+U0] 1/уу°,

(2.26)

где у" и у"— коэффициенты, относящиеся

к лучам

при <т2 = 0.

Пользуясь уравнениями (2.25) и (2.26),

можно

подсчитать

значения т" и ту для любой скорости и для любого луча, относящегося к разным ст2. Для этого достаточно в уравнения (2.25) и (2.26) подставить соответствующие этим лучам значения коэффициентов у° и у® и задаться нужной скоростью. Значения у" и у" получают по экспериментальным данным. Методы

нахождения

этих величин

рассмотрены ниже.

скорос­

Проведя

такие

расчеты

для нескольких значений

ти è1<è2<è3<è4,

можно

построить в координатах

lg т и с

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

тистической

моделью,

где одновременно

происходят сдвиг

и отрыв.

сдвиг без

отрыва наступает

на горизонтальном

Чистый

участке паспорта прочности при с„^спг. Эта область начина­

ется при

гидростатическом давлении

ст2=дг = а 2г (луч 4) (см.

рис. 2.20),

т. е. а 2=р4= а 2г.

имеет значение уп г для

Коэффициент у для этого давления

пределов прочности и уп у для пределов упругости. С воз­ растанием скорости наклон прямых линий 1— 4 (см. рис. 2.21) увеличивается, что объясняется соответствующим изменением величин А и В в условиях (1.2) и (1.7). Предельные напряжения, входящие в паспорт на горизонтальном участке при разной скорости, определяются из уравнений:

tn.r(è) = [ln(è/è(,)A:7’+ t/0] 1/уп. г ;

(2.27)

ту.г(ё) = [1п(е/б0)Л Т + [/ 0]1/уу.г.

(2.28)