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

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

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

Необходимо отметить, что приведенное решение пространствен­ ной задачи релаксации напряжений под шаровым штампом спра­ ведливо при условии S (T, X 0,0054 так как в этом случае пласти­ ческими деформациями под штампом (Н. А. Цытович, 1973) мож­ но пренебречь.

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

4.4. ИНЖ ЕНЕРНЫ Е МЕТОДЫ РЕШЕНИЯ ПЛОСКОЙ

И ПРОСТРАНСТВЕННОЙ з а д а ч к о н с о л и д а ц и и

ИПОЛЗУЧЕСТИ МНОГОФАЗНЫХ ГРУНТОВ

Впредыдущем параграфе мы рассмотрели ряд задач плоской

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

Вместе с тем предложенный Ц. А. Цытовичем еще в 1934 г.

ив дальнейшем усовершенствованный им метод решения плоской

ипространственной задач консолидации значительно упрощает ре­ шение путем приведения этих задач к эквивалентной одномерной.

Точность определения степени консолидации и пбрового давления в грунтах при этом вполне достаточная для применения в инже­ нерной практике.

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

В настоящем параграфе наряду с изложением разработанных ранее решений рассматривается ряд новых задач с учетом особен­ ностей деформирования грунта и граничных условий.

Мощность активной зоны. Мощность эквивалентного слоя, по Н. А. Цытовичу, определяется путем 1сра1внения осадки поверхно­ сти упругого полупространства под воздействием местной нагруз­ ки и осадки слоя под воздействием оплошной нагрузки, что дает h3 = АыЬ, где А = (1 — р)2:(1 — 2р); Ь — ширина нагруженной поверхности; © —коэффициент формы и жесткости.

Для определения А и ю Н. А. Цытовичем составлены подроб­ ные таблицы (1973), что значительно упрощает расчеты.

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

Для нормально уплотненных грунтов при небольшом диапазо­ не изменений напряжений, когда нелинейными свойствами дефор­ мирования можно пренебречь, активная зона будет максимальной и равной ha=2h3. Однако при наличии структурной прочности, на­ чального градиента напора, существенной нелинейной деформиру­

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

Так, по предложенным Н. А. Цытовичем (1967) формулам ак­ тивная зона (рис. 4.20) при наличии структурной прочности мо­ жет быть определена по формуле

тр 2/гэ,

(4.94)

Р

Рис. 4.20 К определению активной зоны сжатия по методу эквивалентного слоя Н. А. Цытовича

Рис. 4.21. Кривые осадки от действия местной нагрузки для нелинейно дефор­ мируемого полупространства (а) и кри­ вые формоизменения и объемного изме­ нения (б) (схема)

а при наличии к тому же начального градиента напора — по фор­ муле

}ia— 2fi9 1 -

(4.94а)

*0-+

Р ~т~ Рстр

 

2fi3yw

 

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

Метод эквивалентного слоя позволяет устранить этот недоста­ ток инженерными приемами — ограничением активной зоны на ос­ нове учета структурной прочности и начального градиента напора.

Вместе с тем, по-видимому, одной из основных 'причин несоот­ ветствия активной зоны, определяемой по теории линейной дефор-

мируемости, с наблюдаемыми величинами является нелинейная деформируемость грунтовой среды.

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

Приведем здесь результаты приближенного решения задач об осадке жесткой полосы, опирающиеся на нелинейно деформируемое полупространство, в случае, когда физическое уравнение взято в виде (рис. 4.21, б), предложенном А. И. Боткиным (1941):

где а,о= (#+cr)tgi|); # it IJV — параметры, характеризующие на­ копление пластических деформаций, определяемые эксперимен­ тально, причем Д г,<0; G, av — модули соответственно сдвига и объемного сжатия; Н — давление связности; tgty коэффициент трения.

Задача об осадке поверхности нелинейно деформируемого по­ лупространства при действии местной (нагрузки, как правило, ре­ шается численными методами механики грунтов (М. В. Малышев, Ю. К. Зарецкий и др., 1973), что связано с большими математиче­ скими трудностями.

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

(4.96)

где р — интенсивность приложенной нагрузки; b — ширина полосо­ образной нагрузки или диаметр загруженной площади; © — коэф­ фициент, зависящий от формы и жесткости конструкции; П — коэф­ фициент, учитывающий накопление пластических деформаций в массиве грунта; р0 — интенсивность предельной критической на­ грузки, определяемая по решениям В. Г. Березанцева (1952) и В. В. Соколовского (1960):

Ро=^хУЬ-\- N 2yh -f N zc,

(4.97)

где h — глубина приложения местной нагрузки; у — объемная мас­ са грунта; N u N2) — коэффициенты, зависящие от угла внутрен­ него трения грунта; значения их даются ,в табличной форме (Н. А. Цытович, 1963); с — сцепление.

Очевидно, что яри небольшом диапазоне изменений приложен­ ной местной нагрузки р эта зависимость носит линейный характер, а с возрастанием ее — существенно нелинейный и имеет асимптоту, параллельную оси осадок (рис. 4.21, а). На этом же рисунке пока­ заны (схематично) зависимости между осадкой поверхности осно­ вания для случаев чисто упругий деформации (пунктирная линия) и для пластической деформации (точка-пунктир). Отметим, что коэффициент Я отличается от коэффициента n v, поскольку в ос­ новании имеют место деформации формоизменения и объемного изменения. Значение коэффициента Я может быть определено по результатам полевых штамповых испытаний по ветви нагрузки и разгрузки, причем по ветви разгрузки определяются упругие ха­ рактеристики. В случае пластичных глин, когда внутреннее трение незначительное и в нестабилизированном состоянии уплотнение не реализуется, зависимость (4.97) упрощается и принимает вид

P o = (2 -fjt)c-fy h .

(4.98)

Для определения мощности уплотняемой толщи воспользуемся методом эквивалентного слоя Н. А. Цытовича. Определим сначала осадку слоя грунта толщиной Я при действии уплотняющей на­ грузки, распределенной по глубине слоя по закону треугольника, полагая, что в нем выполняется условие компрессионного сжатия. Для удобства направим ось z снизу вверх, а начало координат—па уровне основания слоя. Тогда p(z)= pz/H .

Из условия отсутствия бокового расширения имеем

(*)=/>(*)“ а2 И ; V = Р (*)+2*2 (*);

(4.99)

У з " е, (z) = 2e1(z),

 

где 02 (-г) — боковое давление на уровне z; еДг) — осевая дефор­ мация. Для удобства в дальнейшем опустим в скобках z} полагая, что все величину зависят от z.

Совместное рассмотрение выражений (4.95), (4.96) и (4.99) да­

ет

/ 3

[ и -

^ ----- 1

2 { Р ± 2 ' г >

 

( Р - Ъ )

 

 

L

«го— \Р— °г) J

av

(4.100)

где

 

 

 

 

 

 

 

 

 

,.0=(X±2ZL-f-//jtg<J);

G= 2(1 -Мх) ; aD= - 1 —2ц

. (4.101)

Решая уравнение

(4.100),

можно получить выражение

для 02

в зависимости от р, Е, р, Я*, IJV. Тогда на основании формулы (4.95) осевую деформацию слоя можно определить следующим об­ разом:

£ ± ^ - [ 1 + ^ + 2 * , ) ] ,

(4.102)

а осадку — путем интегрирования от 0 до Н:

 

s = \ - ^ ^ . [ \ + n ,(p + 2 4 )]dz.

(4.103)

о

Сравнивая выражение (4.102) с (4.96), получим для Я

н

|и ( 1 - ^ ) + _ _ Д _ ] = Г(/>+2.8)[1+ Л ,(/,+ 2*,)](<2.

 

 

 

(4.104)

Полагая Я* = Я* = 0 и °2= Р — -— ,

получим

 

1

—{4

 

 

t f = 2 б э=о)£

 

...

(4.105)

 

1 —2(х

 

Проведенный анализ изложенного выше решения показал, что глубина сжимаемой толщи существенно зависит от ширины фунда­ мента и величины передаваемой нагрузки на основание. С увели­ чением нагрузки сжимаемая толща увеличивается. Этот вывод соответствует выводам, сделанным Ю. К. Зарецким (1973) на ос­ новании анализа аналогичной задачи, решенной численным мето­ дом.

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

Определение мощности активной зоны имеет важное значение для прогнозирования скорости осадок оснований сооружений во времени по нелинейной теории консолидации грунтов.

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

быстрее.

Решения на базе метода эквивалентного слоя Н. А. Цытовича.

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

Чтобы учесть ползучесть скелета и сжимаемость поровой жид­ кости при односторонней фильтрации, можно пользоваться реше­ нием (4.46), для которого составлены графики и таблицы (Н. А. Цытович, Ю. К. Зарецкий и др., 1967); в случае необходимости прогнозирования вторичной консолидации следует пользоваться зависимостью (4.47).

Д ля учета старения скелета и сжимаемости поровой жидкости при двусторонней фильтрации может быть использовано решение

(4.51а). Здесь следует отметить, что с увеличением мощности эк­ вивалентного слоя, т. е. ширины фундамента, относительная осад­ ка, по этой теории, будет уменьшаться.

При учете ползучести и старения скелета и сжимаемости норо­ вой жидкости при двусторонней фильтрации можно пользоваться решениями вида (4.56), где вместо h следует подставлять 21г.

При учете линейной деформируемости скелета, союимаемости поровой жидкости и ограниченной фильтрации на нижнем уровне активной зоны сжатия граничное нулевое условие для порового давления ставится не на линии z = 2пэ, а на бесконечности. В этом случае решение основного дифференциального уравнения (4.37) существенно меняется. Оно может быть получено с помощью функ­ ции источника на полубесконечной прямой (А. Н. Тихонов, А. А. Самарский, 1966; Н. Г. Араманович, В. И. Левин, 1979), которая

при начальном pw{z, ti)= A 0/ ^ l -----и

граничных

/?„(0, t) =

= pw (оо, /) = 0 условиях имеет вид

 

 

 

А*(*. 0 = 2 У ftcvt

 

 

4Cjjt

J

—exp[

Gfr + CFI)

dZ.

(4.106)

4cvt J)

Интегрирование этого уравнения дает

где Ф (х) — интеграл вероятности.

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

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

U{t)

Подставляя сюда решение (4.107), с учетом уравнения равновесия

б1(*.

j~ A v ( z >t) получим

2 Y cvt

£ /(/)= ! 1 — exp (

Уп

 

(4.108)

Окончательно имеем

 

s(t)= U (t)m v/i3p.

(4.109)

Сравнивая структуру формул (4.107) и (4.108) с решениями пространственных задач консолидации при действии местной наг­ рузки, изложенными в предыдущем разделе, видим, что они во многом совпадают, а именно во всех решениях присутствует ин­ теграл вероятности, зависящий от геометрического параметра за­ дачи, коэффициента консолидации и времени. Если в нашем ре­ шении мощность активной зоны заменить на h0 = Aab, то полу­ чим еще большее сходство.

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

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

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

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

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

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

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

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

Вследствие понижения уровня воды грунт становится тяжелей

и объемная масса его увеличивается на величину

 

у * = у —у',

(4.110)

где Y —объемная масса грунта выше уровня грунтовых вод; у' — объемная масса скелета грунта, взвешенная в воде, равная

Y'=(Y,—Y»)(l—*). (4.110а)

Рис. 4.22. Расчетная схема к прогнозированию осадки земной поверхности при равномерном понижении грунтовых вод

Тогда изменяющееся во времени уплотняющее давление на ни­ жележащую толщу слабых водонасыщенных грунтов будет опре­ деляться выражением

P = y*[h2- H ( i ) ] t

(4.111)

где hz —первоначальный уровень грунтовых вод над поверхностью уплотняемой толщи мощностью h; H (t) — изменяющийся во вре­ мени уровень грунтовых вод над поверхностью уплотняемой тол­ щи.

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

P w i t ) = M x(t)\ pw(z, *!)= ?(«, tj);

/\Л*» /)= Ж 2(^); p w{z, ti)=<l>(z, тх).

Если принять, что скелет грунта обладает ползучестью, описы­ ваемой уравнением вида (4.36), и что в процессе уплотнения име-

ет место уравнение равновесия вида

p(t) -{-р„,(7) = pw(z, t) -f

a(z, t), то дифференциальное уравнение (4.37) примет вид

№pw + Л

 

(4.112)

дР

 

 

где F{t) = AQP (/)- f B p (/) + A0pw(t)+ B pw{t)\

^ _ 4 ( m v t + n m w +

m V i)

/Сф

m Vi + t i m w

Cv

Vw(mv l+ n m w)

•f|(mvl + mv2)

(4.112a)

mvl + nmw

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

Не останавливаясь на подробном изложении решения (4.112), которое дается в нашей -работе (3. Г. Тер-Мартиросян, 1973), при­ ведем окончательный результат

00

р „ ( г ,

0 = Л .(° . <)+-§■ Л .И +

2

v * {t) sin Н

Г ’ (4Л126)

 

 

 

 

л=1,3,...

 

где

 

 

 

 

 

 

,,

V^foi) —X2Vii(Ti)_cx.f I

 

* -----

 

 

 

Xj —Хг

 

 

h ~ h

 

 

СГ

M<-*>

 

 

 

 

 

e*1^-

X2— Xi 1

 

(.*)№;

(4Л12в)

 

+ n — X2

 

A,,2 = ^ -

[ — (a„+bn) ± V a l+ b ly ,

a„—/i+ c „ ( h )2;

 

> ' - ч * ( т р F ^

 

F ( i y

 

/„ (< )= —

p - ( ° . t)+ -\p „ {t)-\-A [p w{h. t)+

,

 

ttn

V

 

 

 

 

 

 

 

}’

 

 

V'e (T1)= - ^ |[/K * i)+ /> « (‘fi)Mo —[а Л°» TI)+

(4.112r)

199