Основные уравнения прочности и деформируемости дисперсных пород Текст научной статьи по специальности «Физика»

© В.А. Миронов, О.Е. Софьин, 2007

УДК 622.011

В.А. Миронов, О.Е. Софьин

ОСНОВНЫЕ УРАВНЕНИЯ ПРОЧНОСТИ И ДЕФОРМИРУЕМОСТИ ДИСПЕРСНЫХ ПОРОД

Семинар № 4

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

Условия предельного равновесия Мора-Кулона и Мизеса-Боткина являются наиболее применяемыми в механике сыпучих сред. Первое условие лучше соответствует экспериментальным данным, но крайне затруднительно в использовании из-за наличия на предельной поверхности сингулярных точек. Второе является удоб-

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

В последнее время получили распространение различные обобщения условия предельного равновесия Мора-Кулона, имеющие регулярную поверхность и основанные на учете вида напряженного состояния, например предложенное в работе [1]. В обобщениях такого рода вопрос о физическом объяснении влияния вида напряженного состояния на прочность сыпучих сред остается открытым.

Наиболее приемлемым с физической точки зрения является обобщение условия предельного равновесия Мора-Кулона на пространственный случай напряженного состояния следующего вида [2]:

313

где ст =-

2

I1213 - 9132

(1)

(2)

‘2

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

т =

п =

; (3)

2

12 — ^*1^2 +а\аЪ +^2^3’ 1з — —

инварианты тензора напряжений; сс (/' — 1, 2, 3) - главные напряжения

(сжатию соответствует с,- > 0); ф -эффективный угол внутреннего трения.

Предельная поверхность, соответствующая условию (1), в системе координат с1,с2,сз имеет вид некруговой конической поверхности, описывающей пирамиду Мора-Кулона, с вершиной в начале координат. Таким образом, условие предельного равновесия (1) предсказывает одинаковые по прочности состояния при напряженных состояниях обобщенного сжатия (с2 = с3) и обобщенного растяжения (с1 = с2), и несколько повышенную прочность при промежуточных значениях с2.

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

Для различных состояний дисперсной породы по плотности дила-тансия может быть как отрицательной (разрыхление), так и положительной (уплотнение). Деформация разрыхле-

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

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

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

Тогда уравнение предельного течения имеет вид [5]:

t — tgф0 5, (4)

где

5 — С\11' + сг2тт' + сгъпп t — -^(^т' -с2т1 ')2 + (<с2 тп' -с3пт' )2 +(с3п/ ' -С\1 п' )2 (5)

нормальная и касательная силы, действующие на площадке скольжения V' и численно равные проекциям напряжений, действующих на площадках предельного равновесия, соответственно на нормаль к площадке скольжения и на нее саму (рис. 1); р0 - истинный угол внутреннего трения; I', т', п' - направляющие косинусы нормали V' площадки скольжения.

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

Если к критически напряженной точке приложить приращения напряжений ёст , то вследствие малости этого приращения направления главных напряжений и связанных с ними площадок предельного равновесия не изменится. При этом тензоры напряжений и скоростей пластических деформации будут коаксиальными. Учитывая это и применяя к уравнению предельного равновесия (1) ассоциированный закон течения Мизеса, получим выражения для направляющих косинусов нормали площадки скольжения в виде:

1' =± 73 3ё2РезР - 32 +у1322 - 33/33

1 (-ер)-ер)

1 зе1рез’ - 3 2 +^13 22 - з/1/3

(-ер)-ер) ’

, 1 3е1Ре2Р - 32 +д/322 - 33133

()( еР )

где

31 = + £2Р + £3Р, 32 = е\е\ + £Р£3Р + ,

33 = -

инварианты тензора скоростей пластической деформации;

еР = Ае? (/ = 1,2,3) - главные скоро-

сти пластической деформации.

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

Если среда пластически несжимаема (£[ +£2’+£3Р = 0), то из выражений (6) следует, что площадка скольжения совпадает с октаэдрической площадкой (I' = т' = п' = 1/л/З).

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

Тангенциальная компонента скорости на площадке скольжения определяется выражением:

Г1=^-32 Ч32 - 3/1 33 . (7)

Условие дилатансии запишем в виде уравнения

(6) ер + Лу' = 0 , (8)

где еР = £1 + еР + £3Р - скорость объемной пластической деформации, Л -скорость дилатансии.

Оказывается, что ассоциированный с условием (1) закон течения соответствует скорости дилатансии Л = - (3/2Ндр.

Предельная поверхность по уравнению (4) также в общем случае не совпадает с предельной поверхностью уравнения (1). Такое совпадение теоретически возможно в случае значительного дилатансионного уплотнения породы, когда Л = (3/2^др (I' = I, т' = т , п' = п) или стремится к ней, если Л ^ - (3/2^др. При Л = -(3/2^др в соответствии с излагаемой моделью истинный угол внутреннего трения р0 равен нулю, что возможно только для нереальной среды с отсутствием трения между частицами. Экспериментальные исследования [7] подтверждают, что «ассоциированная» дилатансия в опыте не наблюдается.

Определяющие связи будем искать в таком виде, чтобы после подстановки их в известные формулы

єІ, = єрї '2 + £2рт '2 + £рп'

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

Проведем сравнительный анализ предлагаемой модели с известными условиями предельного равновесия Мора-Кулона, Треска-Хила, Мизеса-Боткина. Для этого введем обозначения

7 =

п = ■

(11)

2 2 Выражая предельные условия через параметры п/7 = Км и параметр

Ёоде ца = (2<г2 - о! - аъ)/(о! - аъ) ,

получим уравнения [8] Мора-Кулона

1

7

1 + К2

= К2

(12)

УІ, = 2^( -£2р)21'2т'2 + ( -£3р)2т'2п'2 + (3р-£р)п'2!

(9)

результат был тождественен выражениям (5). Такому результату удовлетворяют связи вида

а2т

т

є3р = 2

<73П

(10)

где 2 = d2/d^ - бесконечно малый скалярный множитель (2 > 0).

Треска-Хилла 1 3КТ

7 6-цКт

Мизеса-Боткина

1 (;ЧгКокт +4 2£ + 6 ) Кокт

7“ (2 - ку^/з + 2 ’

(13)

(14)

условия предельного равновесия (1)

Рис. 1. Площадки предельного равновесия и скольжения

условия предельного течения пластически несжимаемой среды (4)

і+7-Н+1

-7+1-1-7

і-7-І+71 = (1«)

=к°°\^ і1+7+М+1 ч1 -7

^ Л

о 1

3 + "7 7

+з

1 -

ґ л2 V 1

1 = К

(15)

где Ям, Кг, Кокт, К - эффективные коэффициенты трения; К0 - истинный коэффициент трения.

Значения теоретически прогнозируемых углов внутреннего трения, по Мору-Кулону, в зависимости от параметра ца в сопоставле-

Рис. 2. Теоретические и экспериментальные зависимости угла внутреннего трения от параметра /ла по Мору-Кулону: а) Опыты Грина и Бишопа (плотный песок из р. Хэм, коэффициент пористости е0 = 0,64); б) Опыты Ёомизе и Крыжановского (кварцевый волжский пылеватый однородный песок, е0 = 0,73): 1 - условие Мора-Кулона; 2 - условие Треска-Хила; 3 - Мизеса-Боткина; 4 - обобщенное условие (1), 5 - условие предельного течения (4) в предположении пластической несжимаемости среды; 6 - условие предельного течения (4) с учетом дилатансии среды, точки - эксперимент

Рис. 3. Зависимость нормальной эн касательной Ь компонент на площадке скольжения (песок из р. Хэм)

3.5' 3' 2.5' 2' 1.5' 1 ■ 0.5' t, МПа

й Па

у

пшу У

\Ф0 S

s>

0 2 4 G 8 10 12

нии с экспериментальными данными по определению пиковой прочности песков, полученными Грином и Бишопом [9], Ломизе и Кры-жановским [10] представлены на рис. 1.

Из графиков видно, что промежуточное главное напряжение ст2 влияет на прочность дисперсных пород. Условия предельного равновесия Треска-Хила и Мизеса-Ботки-на дают завышенные значения прочности, а условие Мора-Кулона вообще исключает влияние на прочность ст2. Значения углов внутреннего трения, полученные по предлагаемой модели, хорошо согласуются с опытными значениями Ломизе и Крыжановского и несколько хуже с данными Грина и Бишопа. На наш взгляд причина этого факта деформационная анизотропия упрочнения, влияние которой значительнее проявляется в опытах Грина и Бишопа. В этих опытах нагружение осуще-

1. Гениев Г. А. К вопросу обобщения условий предельного равновесия сыпучей

ствлялось при ст2 = const и стз = const и возрастающей деформации si (в направлении cti). При этом реализовывалась траектории со значительным ростом среднего напряжения. В опытах Ломизе и Кры-жановкого в основном реализовывались траектории по девиаторной (среднее напряжение постоянно) или близких к ней схемах.

В заключении рассмотрим серию опытов на трехосное сжатие с дренажем насыщенного песка из реки Хэм, выполненную Грином и Скиннером [9]. Опыты проводились при различных значениях начальной пористости песка и величине бокового обжатия, сохраняемого постоянным в процессе всего опыта. Согласно эти опытам разрыхлению соответствует большее значение эффективного угла внутреннего трения, по Мору-Кулону, уплотнению - меньшее значение фм. Обработка опытных данных по уравнению (4) привела к единому значению угла внутреннего трения фо = 16° (рис. 3), что дает возможность принимать его за истинный угол внутреннего трения среды. То же значение ф0 получено и при обработке опытов Ломизе и Крыжановского.

СПИСОК ЛИТЕРАТУРЫ

среды/ Основания, фундаменты и механика грунтов, 1968. № 2. С. 1-2.

2. Nakai T., Matsuoka H. Soils and Foundations, 1983, vol. 23, № 2. p. 26-42.

3. Крыжановский А.Ё., Вильгельм Ю.С., Медведев С. В. Определения угла трения грунтов в приборах трехосного сжатия и срезных приборах/ Основания, фундаменты и механика грунтов, 1985. № 3. С. 2O-23,

4. Ивлев ДМ. О диссипативной функции в теории пластических сред/ В кн.: Механика пластических сред. М.: Физматлит, 2OO1. T. 1. С. 29-33.

5. Миронов В.А, Софьин O.E. Модель предельного течения дилатирующих грунтов при сложном напряженном состоянии/ Будаунщтва. Строительство. Construction. Мн.: БОТУ, 2OO3. № 1 - 2. С. 71-75.

6. Томас T. Пластическое течение и разрушение в твердых телах. М.: Мир, 1964. 3O8 с.

7. Николаевский В.Н. Механические свойства грунтов и теория пластичности/ В кн.: Итоги науки и техники. Механика твердых деформируемых тел. М., 1972. Т. 6. 87 с.

8. Софьин О.Е. Результаты теоретикоэкспериментального анализа моделей прочности грунтов/ Академические чтения Н.А. Цытовича и 2-е Денисовские чтения. М.: МГСУ, 2003. С. 84-88.

9. Бишоп А. У. Параметры прочности при сдвиге ненарушенных и перемятых образцов грунта/ В кн.: Определяющие законы механики грунтов. Механика. Новое зарубежное в науке. М.: Мир, 1975. № 2. С. 7-75.

10. Ломизе Г.М., Крыжановский А.Ё. Основные зависимости напряженного состояния и прочности песчаных грунтов/ Основания, фундаменты и механика грунтов, 1966. № 3. С. 8-11. ЕЕЕ

— Коротко об авторах ------------------------------------------------------------

Миронов В.А., Софьин О.Е. - Тверской государственный технический университет.

Доклад рекомендован к опубликованию семинаром № 4 симпозиума «Неделя горняка-2007». Рецензент д-р техн. наук, проф. С.А. Гончаров.

---------------------------------- ДИССЕРТАЦИИ

ТЕКУЩАЯ ИНФОРМАЦИЯ О ЗАЩИТАХ ДИССЕРТАЦИЙ ПО ГОРНОМУ ДЕЛУ И СМЕЖНЫМ ВОПРОСАМ

Автор Название работы Специальность Ученая степень

ИНСТИТУТ ЦВЕТНЫХ МЕТАЛЛОВ И ЗОЛОТА ФГОУ ВПО СИБИРСКОГО ФЕДЕРАЛЬНОГО УНИВЕРСИТЕТА

НАЗАРОВА Евгения Юрьевна Обоснование технологии разработки месторождений карбонатных пород вблизи селитебных территорий 25.OO.22 к.т.н.