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

книги / Основы прикладной геомеханики в строительстве

..pdf
Скачиваний:
1
Добавлен:
19.11.2023
Размер:
61.85 Mб
Скачать

консолидационнных процессов, полагая, что эпюра дополнитель­ ных напряжений от собственного веса треугольная:

s =

mV2y*hl

(4.113)

2

Для наглядости полученных решений рассмотрим пример. Оп­ ределим конечную осадку поверхности Земли при следующих па­ раметрах (рис. 4.22): h2= 2000 см; mv2= 0,01 см2/кгс; /г=1000 см; mv= 0,05 см2/кгс; у* —1 т/м3=0,001 кг/см3, тогда

 

 

 

5 = s , + s2=

 

у * h i t r i v i - f У * f t 2m v =

 

 

 

=

0,001.2000 Ц - 2000-0,01 - f 1000-0,05]= 120 CM.

 

 

Для

 

определения

скорости

/■

 

 

осадки

следует

пользоваться

вы­

 

 

шеизложенными

решениями

од­

/ / .'•‘//'Яш

 

1

1 1 1-

 

номерной

задачи

консолидации.

1

1

1

К у с т о в а я о т к а ч к а из во­

1/ 1 4

l/iiL

д о н о с н о г о

п л а с т а .

 

Одним

из наиболее интересных объектов,

 

 

 

на котором детально исследова­

 

 

 

но сжатие горных пород под вли­

 

 

 

янием глубокого

водопонижения,

A A A ■'Y V

A

A A A

является Белозерское железоруд­

 

 

 

ное месторождение (В. А. Миро­

 

 

 

ненко,

В.

М.

Шестаков,

1974). Рис. 4.23. Расчетная схема к прогно­

Снижение

напора

воды

в

водо­ зированию оседания земной

поверх­

носных батских песках на 200 м

ности при откачке подземных вод из

привело

здесь

к

интенсивному

глубоких водоносных горизонтов

сжатию песчано-глинистой и мер­ гельно-меловой толщ, залегающих на глубине около 300 м; в на-

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

Рассмотрим задачу уплотнения грунтов напорного пласта мощ­ ностью ht из которого одновременно производится откачка воды из многочисленных скважин диаметром 2/4 на расстоянии друг от друга 2г2 (рис. 4.23).

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

Особенность этих задач заключается в том, что начальные ус­ ловия заданы в соответствии с условиями установившегося равно­ весия

/?w('ri ) = 2 Yf/Z/— a(xi)»

(4.114)

i=1

 

где yi — объемная масса /-го слоя грунта; hi — мощность /-го слоя; Pv>(%\), a ( x i ) — начальные значения порового давления и эффек­ тивных напряжений в водоносном пласте.

Таким образом, решение задачи об оседании земной поверхно­ сти при кустовой откачке подземных вод из водоносного горизон­ та приводится к решению осесимметричной задачи консолидации, т. е. к уравнению (4.75) и (4.80).

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

Для определения стабилизированной осадки вследствие полно­ го рассеивания начального порового давления [см. формулу (4.114)] следует пользоваться обычным решением одномерной задачи уп­ лотнения:

s[oo)=mmh p J x l)t

(4.114а)

где mv0 — коэффициент относительной сжимаемости и h — мощ­ ность слоя водоносного горизонта.

Для определения изменяющейся во времени осадки следует поль­ зоваться зависимостью вида s(t) — s(oo) U(t) [где U(t) — степень консолидации, определяемая на основании «решений в 4.3].

Так, например, для условий свободных деформаций можно по­ льзоваться решением (4.77а), а для условий равных деформаций— решением (4.81).

В качестве примера рассмотрим следующий случай. Пусть во­

доносный

горизонт мощностью 20 м находится на глубине h =

= 300 м;

первоначальное давление воды в нем составляет Pw(xi) =

=40 кгс/см2, а средняя объемная масса вышележащей толщи рав­ на у = 2 т/м3. Требуется определить величину осадки и скорость ее изменения, если известны следующие характеристики грунта:

е0=1; /пго=0,002

см2/кгс;

т ш=0,001 см2/кгс; /Сф=10~8 см/мин;

fi = 0,1 м; г2 — 1,0 «м. На основании формулы

(4.114а)

s (оо)=0,002-2000-40 =160

см;

 

 

/ = — =10,0;

mV)

0,001

=0,5;

W = < ? 0

0,002

ri

 

 

 

у _ Crt _____ {____. Q ________ _______

г2

250-10000

Yw

О"Ь /z/Ищ)

_________ HHJ________ __

10-8

1

0,001 (0,002-1-0,5.0,001) ~ 0,001-0,0025 = "25o' См2/мин-

На основании кривых рис. 4.15

и

0,3

0,4

0,5

0,6

0,7

0,8

Тг

10-4

10“ 2

зхю -1

5хЮ -1

1,0

1,5

t мин

250

2,5'ХЮ5

7.5ХЮ5

1,25Х106 2,5X10® 3,75X10®

Таким образом, полная осадка произойдет примерно через 10 лет. О б р а з о в а н и е д е п р е с с и о н н о й воронки. При откачке подземных вод из одной скважины большого диаметра или из близ­ ко расположенных друг от друга

нескольких

скважин

(кустовая

 

 

 

откачка)

в

подземном

простран­

 

 

 

стве

образуется

Депрессионная

 

 

 

воронка

большого диаметра (рис.

 

 

 

4. 24), вследствие чего в окружа­

 

 

 

ющем массиве изменяется напря­

 

 

 

женное состояние, что приводит в

 

 

 

конечном итоге к оседанию зем­

 

 

 

ной

поверхности

и образованию

 

 

 

мульд оседания.

 

 

 

 

 

 

Для

прогнозирования

напря­

 

 

 

женно-деформированного

состоя­

Рис. 4.24. Расчетная схема к прогно­

ния массива горных пород вслед­

зированию оседания

земной

поверх­

ствие образования депрессионной

ности при образовании депрессионной

воронки вследствие

откачки

подзем­

воронки

может быть использова­

ных вод

 

 

но решение задачи о силе, прило­ женной внутри упругого полупространства (Р. Миндлин, 1956), в

предположении, что по всему объему внутри депрессионной ворон­ ки в каждом единичном объеме действует сила интенсивностью у* (разность объемных масс пород до и после образования депресси­ онной воронки).

По Р. Миндлину, осадка поверхностных точек упругого полу­ пространства от действия сосредоточенной силы Р, приложенной на глубине с от поверхности, определяется по формуле

+

( 4 Л , 5 )

где

R =='Vс~ -\-х* -f-у".

На основании этого решения легко определить осадку в произ­ вольной точке поверхности полупространства от действия равно-

меньше рассматриваемых величин с*, можно пренебречь значения­ ми Рг, тогда получим

i ) + i -

V 1 +a|

Для значений A t составлены таблицы (табл. 4.4), что в значи­ тельной степени облегчает расчеты.

 

Значения коэффициентов Ai=f(cii, ц0)

Таблица 4.4

 

 

 

 

 

^0

 

ai.

0.2

0,3

0,4

0,5

 

0,1

0,0135

0,0125

0,0113

0,0099

0,5

0,3110

0,2887

0,2616

0,2298

1,0

1,0296

0,9588

0,8716

0,7677

5,0

8,5132

8,0229

7,3686

6,5504

10,0

18,0961

17,1011

15,7440

14,0250

20,0

37,2880

35,2905

32,5320

29,0125

50,0

94,8832

89,8842

82,9248

74,0050

100,0

190,8816

180,8821

166,9224

149,0025

200,0

382,8808

362,8810

334,9212

299,0012

300,0

574,8805

544,8807

502,9208

449,0008

400,0

766,8804

726,8805

670,9206

599,0006

500,0

958,8803

908,8804

838,9204

749,0005

1000,0

1918,8801

1818,8802

1678,9202

1499,0002

10000,0

19198,8802

18198,8800

16798,9200

14999,0000

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

р

E \h\ -J- E<ihi 4- •••+ Enhn

<

_

+ V-2^2 + •••+ V-nhn

nP“

h

 

h

где h — суммарная мощность рассматриваемого массива; hn—мощ­ ность /г-го слоя.

В качестве примера рассмотрим случай оседания поверхности Земли вследствие образования депрессионной воронки при откач­ ке подземных вод из скважины диаметром 1 м на глубину 100 м, при этом диаметр депрессионной воронки на поверхности равен 1000 м, а депрессионная кривая описывается логарифмическим за­ коном в виде а = а0 (In 100 — lnz), где а0 — коэффициент пропор­ циональности, равный 500 м, z—глубина от поверхности, м. Тогда для элементарных слоев грунтов, обладающих модулем деформа­ ции EQ=200,0 кгс/см2 и коэффициентом Пуассона р=0,3, получим следующие значения коэффициентов Ли полагая толщину эле­ ментарного слоя Ас = 10 iM, пористость п — 0,39, объемную массу у» == 1 т/м3, а ys — 2,65 т/м3:

zi = 10 м

гг=20 м

23=30

м

24 = 40 м

2б = 50 м

Г! =500,0 м

г2= 349,5 м

Гз=261,45 м

г4= 198,95 м

г 5 — 155,5 м

Ci=5 м

Сг= 15 м

с3=25 м

с4=35 м

£5= о5 м

си=200

аг=22,7

а3= 10,458

а4=5,684

а5= 2,827

/11=362, 87

Аг= 40,2

А3= 17,92

Д4=9,25

Д5 = 4,18

 

гв=60 м

г7= 70 м

20=8О м

29= 90 М

 

г6= 110,9 м

г7= 77,45 м

г8=48,45 м

гэ=22,9 м

 

с6=65 м

с

м

 

Сз=95 м

 

7=75

 

с8= 8 5 м

 

 

а6= 1,706

а7= 1,032

а8=0,57

ад= 0,241

 

Лв=2,10

А7=0,99

48=0,364

Лд=0,075

По формуле (4.110а),

Y*=Y-(YJ- Y « )(l - » )= 2 —(2,65— 1) (1-0,39) =

 

= 1 т/м3=0,001

кг/см3.

 

 

На основании выражения (4.115г)

s0 = 254 см.

 

Таким образом, в однородной грунтовой толще с модулем £ 0 =

= 200 кгс/см2 при

образовании депрессионной

воронки

глубиной

100 и диаметром

1000 м максимальная осадка

в центре будет

равна 254 см, что совершенно недопустимо в условиях

городской

территории. Для построения чаши прогиба поверхности земли мож­ но пользоваться подобием кривых осадки и депрессионной ворон­ ки (Н. М. Герсеванов, Д. Е. Полышин, 1948):

s(r)= s 0exp ( —

(4.116)

где а0 = 500 м; г — рассматриваемый радиус чаши

прогиба, м;

а -параметр кривой чаши прогиба.

 

Для точного определения чаши прогиба поверхности Земли мо­

жно решить

 

интеграл (4.115а)

с учетом его упрощения в виде

s(r) =

 

2 лЕ 0

2(1 — р.0) (l^г?—2 /у сое о + с2+ г2 —

 

 

Р О + fa)

 

 

 

 

 

 

 

УГо — 2гго cos ср —J—с2 -j- г2)

2г cos ср(1 — р.0)

X

Х1п-

У г\

2 Г\Г COS <Р +

с2 +

г2 +

Г\ — Г C0S ?

_j_

1

Уr\

 

 

 

 

 

 

 

 

 

 

 

Г COS <р

 

 

2 r r o COS <? +

С2 +

г2 +

Го —

^1 + Г S^n У j

X

с2 +

г2 — ГГо co s <р

 

 

С2 — /*2 — ГГ\ COS <р

фр.

 

 

 

 

 

 

УГj — 2rri COS + с2 +

|/ г\— 2rro cos <f +

с2 4- ,-2

г2

 

 

 

 

 

 

 

 

 

 

(4.115Д)

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

опд

В данном случае iz = 0; ix — sin р. Интересно отметить, что при угле р = 30° в водонасыщенном пласте вследствие установивше­ гося фильтрационного потока вдоль пласта возникают фильтраци­ онные силы интенсивностью 0,5уш, что соответствует интенсивности силы от действия собственного веса скелета, взвешенного в воде.

На стационарное напряженное состояние, описываемое уравне­ нием (4.117), будет накладываться нестационарное напряженное состояние, вызванное действием граничных нагрузки p(t) и напора Н(0, t). Решение этой задачи нами было рассмотрено выше [см. (4.1126) и (4.112д)].

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

Суммируя эти решения с (4.117), получим общее решение по­ ставленной задачи в виде

aziz. 0 = [ Y , ( A - 2:) + / ;(/)]c o e ?-jP w(z,

t)\

0=(Y /+ Yw)(/f— zbiop; ax (z,

t)= to 2(z, *)•

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

 

°/

|

— в/о

(4.118)

 

•'ll (а)

т,2 ( о )

 

 

гдет]1 (а)

и г|2 (а) — коэффициенты вязкого

и вязкопластнческого

течения;

шо— порог ползучести

при

формоизменении, определяе­

мый выражением

 

 

 

 

cio=(a - l- ^ )tg ,l»,

(4.118а)

где Н — связность; tgop — коэффициент трения.

Очевидно, что при напряжениях Ог>Пго первым членом в выра­

жении (4.118) можно пренебречь,

так как т|| (a) ^>ri2 (ст)

и ско­

рость вязкопластического течения

значительно больше

скорости

вязкого течения, а при Oi<Oio следует пренебречь вторым членом, так как деформации пластического течения отсутствуют.

Аналогичное уравнение течения необходимо рассматривать для описания относительного смещения грунтовых блоков по плоско­ сти скольжения. В данном случае наиболее вероятной плоскостью скольжения является плоскость z — const. Тогда скорость смеще­ ния по этой плоскости может .быть описана уравнением вида

*

Тагz

| Чгz *0

(4.1186)

\

 

ri2( a z )

Ю

 

где т)!*(г) и rj2* (-г) — коэффициенты вязкого и вязкопластическо­

го скольжения по образовавшейся плоскости; То — порог ползуче­ сти при плоскостном -скольжении, определяемый уравнением

 

*о=в*(*£?)+<?»

(4.118в)

где с — сцепление; ср — угол внутреннего трения.

 

Уравнение (4.118)

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

переходе к декартовой системе

координат записывается

в виде

з =

2 ^ 2(з) +

^ а) tg ‘ j(sz-e);

(4.118г)

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

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

зависит от х, а

е* = еу = у*v — yyz — 0, а также полагая, что е2<

jCyxz, получим

еi = yXz, т. е. вместо выражения

(4.118г) получим

 

2 % (з ) - ^ = 2 т г ,г- ( з 2+ Я )‘е(|.

(4.118д)

Подставляя сюда значения напряжений для установившегося режима напряженного состояния из формул (4.117), получим

[(Y'+ YJ (A — z) si-n р — ^ [ Н

у ' { h - z ) cos р]J.

Полагая, что Tj2 (cr) = т]2 = const, получим

после интегрирова­

ния

 

1+ Y”) ------2----- 5

Х

х Г Я г - Y '

(4.119)

Постоянную интегрирования С можем определить, полагая в общем случае, что на плоскости z = 0 имеет место вязкопластиче­ ское скольжение, причем rj2*(a) = r\2* — const, тогда окончатель­ но имеем

Vj,=

|(V' + УдГ)~

~ ^ ~

г)~

sin J3-----j-tg + X

X f v ' V - t o - * * .. coeP +

A r z ] } +

- i - [ ( Y ' + V . ) A 3 i n ? - c -

 

— у'Л сое p tg <р].

(4.119а)

Из этой формулы видно, что скорость смещения vx зависит от z нелинейно и что вязкопластическое течение начинается с неко­ торой глубины hu, выше которой могут иметь место деформации вязкого течения, описываемого уравнением вида

= ~[(У '+ V J Л2~ ~ 2)2 sin р j - f «*(*..)» (4.1196)

где vx (hn) — скорость смещения на характерной глубине hn, ниже которой имеет место вязкопластическое течение, определяемое урав­ нением (4.119а).

Если на всех глубинах касательные напряжения не превышают порог ползучести, т. е. hu> h, то зависимость (4.1196) будет спра­ ведлива для всей глубины слоя; при этом следует учитывать гра­ ничные условия обобщенного вида:

Ч с= — [ ( v '+ v j

sin р]

+-^г(У / + V j h sin р.

“П L

J

тн

 

 

(4.119B)

Н е у с т а н о в и в ш е е с я с о с т о я н и е . В случае неустановившегося напряженного состояния ,в пространстве и во времени под воздействием переменной нагрузки p(t) и напора H{t) в уравне­ ние (4.118д) следует подставлять напряжения, определяемые урав­ нениями (4.117а), а также изменяющееся во времени поровое дав­ ление, определяемое в общем случае уравнением вида (4.1126). Тогда с учетом граничного условия при 2 = 0 , vx = vQполучим

ч>х(z> t) =

|(у' +

уJ

sin p----- Y

tg <bX

X |jv'

~ Z)~■+P {t) zj cos p-f H z pw(0,

t)z +

- - - [pw{hi

 

 

 

"b 1) ] -

 

 

/1=1,3,...

 

 

 

~ H z V+

{(?'+Y W) A sin p- c —

 

 

)

rl2

 

 

-1 Ь 'А + М * )]« * Р -Л .(0 ,

<))tg<p).

(4.120)

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

Приведем результаты комплексных полевых и лабораторных исследований динамики развития оползневых смещений Ахангаранского оползневого склона, проведенных аспирантом кафедры