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

книги / Прочность грунтов и устойчивость оснований сооружений

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

Рис. II 1.4. Зависимость между углами

Рис.

III.S. Зависимость Взар ОТ

sin фо/sin ф

at и Q2

для разных значений <р

 

при разных значениях

9

 

 

 

Рис. 111.6. Зависимость Ркр/рм акс

Рис. 111.7. Зависимость Ркр/Рмахс от £о — слева

от Б1П9о/51Пф при разных значе-

н от 1До справа для

несвязного грунта

ннях 9

(с =

0)

аг = const, где р определяется из выражения (III.15):

_________________ — [siny-f-sin(2ai 4-ф)](я4~2д »-hsin 2ai)________________ __

(1 + c o s2 a j){l 4~cos2ai -J- sin Ф[sin 9 +

sin (2a 1-f'(p)'“ (Jl-h 2ai)cos(2ai + 9)]}

_____________________ [sin <p +

sin(2tt2+

ф)](я ~~ 2«2 — sin 2«2)_________________ _

( i + cos 2a2){ 1 + cos 2a2+

sin <p [sin <p 4 - sin (2a24- ф) 4- (я - 2a2)cos(2a24- фШ

(111.29)

121

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

рис.

III.4. Участок /, если

<ц == <Х2=

0зар > —ср,

горизонтален

(си =

const) до пересечения

с линией

максимума,

образующей

участок 2, а если си = аг = 0зар <

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

с той же линией. Участок 3 оказывается всегда горизонтальным.

Таким образом, ход кривой на участках 2 и 3 оказывается опре­

деленным, а положение участка

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

0зар. В решении Федорова (1958, [32]) точки А ц В (см. рис. 111.4) соединялись прямой. В предельном случае, если 0зар = —(лт/4 +

+ ф/2), А и В лежат на одной вертикали,

если Q3ap = п /4 —

ф/2, А и В лежат на одной горизонтали.

 

 

На рис. III.5 приводятся вспомогательные зависимости,

рассчитанные в соответствии с формулами

(III.24)

и (III.25),

0зар от sin фо/sin ф, а на рис. 111.6 — зависимости

ркр/ р иаКс от

sin фо/sin ф. На рис. III.7 даны те же кривые для сыпучего грунта (с — 0), но построенные в зависимости от £о = ?о/рмакс-

Таким образом, оказывается, что нагрузка ркр, при ко­ торой происходит зарождение (начало образования) пластической области, зависит от коэффициента бокового давления £о, причем ркр имеет наибольшее значение при £о = 1. При одном и том же значении угла внутреннего трения ф с увеличением удельного сцепления с уменьшение ркр замедляется. Местоположение за­ рождения пластической области оказалось существенно завися­ щим от угла внутреннего трения ф и коэффициента бокового давления £о. Предельное состояние может возникать как вдоль луча, направленного под фундамент {при £о < 1), так и в сторону от него (при £о > 1).

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

4. ТЕОРИЯ ПРЕДЕЛЬНОГО НАПРЯЖЕННОГО СОСТОЯНИЯ СЫПУЧЕЙ СРЕДЫ

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

Теория предельного равновесия была четко сформулирована применительно к плоской задаче, где система уравнений оказа­ лась определенной. Позднее при наличии допущения о полной сыпучести [1] была сформулирована Пространственная осесим­ метричная задача. Рассмотрим эти peuieHHH.

В плоской задаче имеются два уравнения равновесия:

122

дох |

дхху__

.

дхху__ л

(III.30)

дх

ду

ду

дх

 

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

(а. — 0$ + 4т?, = (о, + о„ + 2 с ctg <p)s sin2 <р.

(HI.31)

Поскольку уравнение (Ш .31) относительно напряжений квад­ ратное, то допускаются два решения задачи, а каким из них воспользоваться — выбирается по физическому смыслу рассмат­ риваемой задачи.

Долгое время имелось лишь простейшее решение приведенной системы уравнений для частного случая тху = О, когда система сразу линеаризуется относительно напряжений. Это решение при­ надлежит Ренкину. В связи с тем, что в исходную систему не входят деформационные условия, ее следует считать статически определимой. Далее изложим вкратце метод интегрирования и педварительных преобразований исходной системы по В. В. Соко­ ловскому [42, 43}, основанный на предложениях С. А. Христановича. Предлагается подстановка, тождествейно удовлетворяющая уравнению (Ш .31) , и тем самым сводящая систему из трех урав­ нений к системе из двух уравнений. Порядок системы (второй) при этом сохраняется. Если положить, что

Ох =

о (1 + sin q>cos 2wi) — с ctg <р;

 

о , =

а (1 — sin ч> cos 2<oi) — с ctg <j>;

(111.32)

тж, =

о -sin <р sin 2o)i.

 

то получим две новые переменные а и а>(, определяющиеся следующим образом:

о = (ох +

а,)/2 +

с ctg <р;

 

 

tg 2(oi =

2тхУ/(ох а Д

 

(111.33)

Далее вводятся подстановки:

 

 

 

ctg q>

..

сг

я

<p

(111.34)

1 ---------2

 

f

~ 7 ~ s

-

 

где <то — произвольная, выбранная заранее величина, имеющая размерность напряже­ ния для того чтобы под знаком логарифма была безразмерная величина <можно, например, положить ао = с или сто = 1 МПа и т. д.).

Затем используется еще одна подстановка — вводятся новые переменные | и rj, причем

i = x + “ i; ч = X^ wi.

(III.35)

Подставляя зависимости (Ш .32) в выражения (IH.30), полу­ чим уравнения равновесия:

(I + sin tpcos 2 m ) ~ + sin у sin 2<a,

_

J

ox

ay

I

i

«

123

,

_

do)|

_

dan4 _

ctg<y...

sin 2CDI --------cos 2<oi

 

— cos qj^sii

 

dx

 

duy '

2oo '

 

 

 

sin Ф sin 2«oi

+

 

 

 

(Ш.36)

(1 sin q>cos 2wi) 4*- +

 

dx

 

 

dy

 

 

 

dm

 

dm i

л

( cos 2ioi —-----h sin 2<й| —-— 1 =

0.

 

 

dx

 

dy s

 

Для упрощения дальнейших выкладок вводятся обозначения:

cti =

02 =

у sin (o)i — ц)

 

sin ф cos (<oi + И)

(Ш.37)

у sin (CDI + Ц)

 

2o sin q>cos(*>i — (i) ■}

 

Умножаяпервое уравнение (III.36)

на sin (coi — p)

или на

sin (ом +

p),

а второе уравнение — на cos (<oi — p) или соответ­

ственно

на

cos (ал + p) и складывая

эти уравнения с

учетом

зависимостей (III. 13), после введения |

и ri по формулам

(III.35)

получим следующие два уравнения:

 

 

g + t g ( « . + n ) g = « . ;

(III.38)

Далее выясняется следующий вопрос. На плоскости х, у пусть задана какая-либо линия с уравнением у = у (х). Следует установить, когда можно определить значения первых производ­ ных от искомых функций I и 1) из этой системы и.когда это сделать нельзя. Пока имеются два уравнения (111.38) и четыре неиз­ вестных производных. Чтобы система была замкнутой, присоеди­ ним еще два уравнения, которые справедливы вдоль указанной линии у(х):

dS- 3 r d' + s ; * ;

(111.39)

dii = ^ d x + ^-dy. dx dy

Таким образом, уравнения (III.38) и (111.39) представляют собой неоднородную замкнутую систему из четырех уравнений

относительно неизвестных производных Ре­

шая эту систему относительно производных, получим следующие их значения:

д% _

flidy — tg((oi 4 -n )d S .

_

dg — aid*

 

dx

dy — tg (<oi + |i)"dx ’

dy

dy — tg (m + ц) dx

 

dr\

fl2dу — tg (coi — |i) dT|

дц _

d4 —^ d x

(III-40)

dx

dy — tg((Dj — ji) dx #

dy

dy tg (coi — |i) dx

 

124

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

Характеристики, а точнее их проекции на плоскость х, у

играют

важную роль в

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

сыпучей

среды. Приравнивая

в выражениях (III.40) отдельно

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

~

= tg(wi + ц); т ^ = « ь

(Ш.41)

ах

ах

 

Для второго семейства характеристик получаются уравне­ ния:

i | = t g ( « , , - ^ ) ; ^ = а2.

(111.42)

Отсюда становится ясным, что характеристики составляют углы с осью х либо (а>1 + |Д либо (coi — р). Поскольку 0)| — угол между направлением наибольшего главного напряжения и осью х (положительные направления отсчета считаются от оси х против часовой стрелки), то ясно, что характеристики совпадают с линиями скольжения, так как линии скольжения составляют углы гЬр с направлением наибольшего главного напряжения. Угол между характеристиками равен 2р. Через каждую точку рас­ сматриваемой области в плоскости х, у проходят две характе­ ристики. Линии разрыва, а их уравнения те же, что и (111.41) и (III.42), представляют собой либо линии скольжения, либо огибающие линий скольжения. Компоненты напряжений претер­ певают на линиях разрыва конечные скачки.

Линиями скольжения называются линии, вдоль которых удов­ летворяется условие-

ITJ = а„ tg Ф + с,

(II 1.43)

где т„/ и агл — соответственно касательная и нормальная компоненты, действующие по площадке с внешней нормалью п н связанные с компонентами напряжений в системе

.координат х. у следующими формулами [42]:

ап=

1/2

(а* + ау) - f 1 /2 (а* — оу) cos 2а + хху sin 2а; 1

 

о< =

1/2

(ох 4- Оу) — 1/2 (ojr — Оу) cos 2а — хху sin 2а;

У

<111.44)

т„| =

1 /2 (ох — ау) sin 2а -f- хху cos 2а,

J

 

где at — компонента напряжения, перпендикулярная направлению п\ а — угол, определяемый по рис. III.8.

125

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

1 — первое

семейство

линий

скольжения;

2 — второе

семейство

линий

скольжения

Рис-. Щ .9. Семейства характеристических, линий

Удобными также являются^ формулы:

=

а [1 ±

sin ф cos 2 (<I)I в)] с ctg ф;

I

Tтnмt =

 

 

 

(1II-45)

a sin ф sin 2 (coi — а),

 

*

где а имеет то же значение,цк что и в-формуле (111III..33),►и

 

 

=

1/2 (ai

+ аз) ± 1/2 (oi — а3) cos 2ам; |

(Ц 1.46)

 

 

 

 

 

Тху =

1/2 (а, — Оз) sin 2<oi.

*

где oi > аз — наибольшее и наименьшее главные напряжения.

 

 

Решение системы уравнений (111.41) и (111.42)

 

вед ется с при

12 6 '

Рис* III. 10. К определению несущей способности, основания по теории предельного равновесия сыпучей среды

а — прямая задача; б — обратная задача; стрелка — направление хода расчета

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

ются

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

I

и

т|

в двух

точках

I и 2 (рис.

III.9), находящихся

на

двух

соседних

линиях

скольжения,, и

координаты х и у

точек 1 и 2, можно определить координаты точки 3, которая будет располагаться на продолжении этих линий скольжения,

иза^ем значение функций £ и т] в этой точке 3. Ход вычисления

исхема достаточно подробно и обстоятельно описаны ы книгах

В. В. Соколовского [42] и В. Г Березанцева [2]. В. Г. Березанцев приводит удобную схему, показывающую порядок вычисления с тем, чтобы удовлетворить требуемым граничным условиям. Порядок этот является определенным и исходит из самой задачи.

Как же ставится задача в теории предельного равновесия? Первая — прямая постановка задачи (рис. ШЛО, а) заключается в том, что вдоль положительного направления оси у задана нагрузкам по величине и по направлению (наклону по отношению к вертикали). На отрицательном участке оси у должна быть пригрузка. Предварительно задается направление этой пригрузки (вертикальная или наклонная) и отыскивается ее величина'. Обратная постановка задачи — по данным величине й направле­ нию пригрузки и заданному направлению нагрузки найти последнюю.

По характеру напряженного состояния всю предельно напря­ женную область основания можно разделить на три зоны — I и II (рис. ШЛ О, б) обобщенные зоны Ренкина и III обобщенная зона Прандтля; В: случае невесомой среды и равномерной: на­ грузки имеются простые решения, полученные Ренкиным и Прандтлем.

В зонах / и II характеристики (линии скольжения) в случае весомой среды — прямые или слабо искривлены. В зоне III они представляют собой логарифмические спирали или слабо отли­ чающиеся от них кривые. Точка 0 представляет собой особую

127

точку, которая на плоскости характеристик £, т) развертывается в вертикальную прямую с уравнением | = const. Поскольку зоны I и II дают характеристики, близкие к прямолинейным, рассматривать напряжение в них уместнее в декартовых ко­ ординатах. Что же касается зоны III — радиальной зоны, то ее рассмотрение удобнее проводить в полярных координатах. Дифференциальные уравнения равновесия в полярных координа­ тах г, 0 имеют вид:

дог .

1 дтго

+

Or — вО

у cos 0;

 

дг

 

г

дВ

Г

 

 

 

(II 1.47)

дтл

,

1

дор

 

2тго = —у sin в.

^

}

дг

+

г

дв

г

 

 

Из уравнений (III.45) получим, полагая нормаль п совпадаю­ щей по направлению с г, t и 0:

 

= o ( l

± sin <pcos2(oj) + cctgf<p;

 

т,о =

о sin <p sin 2о)г;

(111.48)

 

о =

[(af +

CF© )/2] c ctg ф,

 

где ci>2 — угол наклона главной оси к радиусу г, причем *>2= o>i — 0.

Очевидно, что если положить а = yrs(0), то система уравнений (Ш .47) в частных произведениях превращается в систему обыкновенных дифференциальных уравнений:

sin ф sin 2о>2

ds

(s+o-

 

d0

 

— +

2s sin ф cos 2сог|

 

—s (1 +

sin ф cos 2©2) + cos 0;

(III.49)

 

ds

 

 

I — sin ф cos 2<D2 —— + 2s sin ф sin 2(Ог1

 

 

d0

 

 

=

—s sin ф sin 2(02 — sin 0.

(111.50)

Разрешая эту систему относительно

производных

 

cos 0 — sin ф cos (2(02 +

0) — s cos2 ф

5 + 1

2s sin ф (cos 2(02 — sin ф)

ds

__—sin (2co2 + 0) + s sin 2<ог.

d0

cos 2(02 — sin ф

 

получим:

(III.51)

(III.52)

При (02 = ± p c = ± (я /4 — <p/2) знаменатели в правых частях уравнений (III.51) и (III.52) превращаются в нули. Линии скольжения, совпадающие с характеристиками, наклонены к ра­

диусу под углами (02 ± |АС, а к оси х

под углами ©i ±

рс.

Дифференциальные уравнения линий скольжения

 

r 4 ^ s= tg (со2 ±

щ),

(Ш.53)

откуда получаем

 

 

6

 

 

$ctg(u, ± Pc) do

 

r = roe°

 

(Ш.54)

128

Особо простым получается решение при у = 0. В. В. Соко­ ловским дано решение задачи для случая переменного значения

0)2* Однако имеется частное решение

рассмотренной

системы

при (02=

const

(решение этой задачи,

полученное нами,

при­

водится

далее

в п. 5 настоящей главы). Этот случай

как

раз

и относится к решению Прандтля.

При осесимметричном напряженном состоянии уместно приме­ нить цилиндрическую систему координат 70г, представленную на рис. III.11. Дифференциальные уравнения равновесия имеют в данном случае вид:

В осесимметричной задаче В. Г. Березанцев [1] принимает так называемое условие полной сыпучести, согласно которому два главных напряжения (промежуточное и минимальное) оказываются равными. Это лишает задачу неопределенности, так как в пространственном случае имеются шесть неизвестных напряжений, три дифференциальных уравнения равновесия и одно условие предельного равновесия. В силу осевой симметрии одно из уравнений раэновесия выпадает, остаются только два. Два касательных напряжения также оказываются равными нулю. Таким образом остаются два уравнения равновесия. Условие полной сыпучести позволяет сделать систему замкнутой. Усло­ вие предельного равновесия имеет вид

(о, — Ох)2+ 4т?,

= sin ф-

(III.56)

(<Тг 4-

+ 2с ctg ф)2

 

 

Условие полной сыпучести следующее: <7о = 02 = 03,

где <т0 — главное напряжение.

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

с*1 — 03

От + Ог

СС^Ф

(111.57)

2 sin ф

2

 

 

и угол юз, составленный наибольшим

главным напряжением о i

с осью Or (рис. III. 12). Уравнению

предельного равновесия

можно тождественно удовлетворить, если принять:

129

ar — cr (1 + sin q> cos 2юз) — c ctg <p;

 

ox =

a (1 —- §in ф cos 2<оз) — c ctg qp;

(III.58)

Trx =

о sin cp sin 2соз;

 

Oo =

a (1 — sin ф) — c ctg ф.

 

Далее решение ведется с помощью подстановки выражений (111.58) в уравнения (111.55). Получается каноническая система уравнений, которая интегрируется численным методом, анало­ гичным применяемому при решении плоской задачи теории пре­ дельного равновесия. Интегрирование ведется численным методом по схеме, аналогичной решению уравнений для плоской задачи. Если уравнения характеристик

a = о (*, у) = const; Р = 0 (*, у) = const,

(111.59)

то сетка характеристик принимается за криволинейную систему координат на плоскости х, у, причем а и р принимаются за незави­ симые переменные, а х, у у 1, т| считаются функциями а и р . Кано­ ническая система уравнений имеет следующий вид:

(III.60)

! = * < ■ ' + * > $ =

а также

д 1 __

дх дт]__ дх_

др

1 др ’ да

(III.61)

2 да

Решение канонической системы уравнений (III.60) и (III.61), для которого определитель не обращается тождественно в нуль, является решением уравнений исходной основной системы (III.38). Аналогичная система получается и при использовании полярных координат в случае плоской деформации и при решении осесимметричной задачи теории предельного равновесия сыпучей среды. Интегрирование ведется численным способом. Получить результа­

ты

в

замкнутом

виде

представляется возможным

лишь очень

в малом числе случаев.

 

 

 

 

В

заключение

приведем

важные

результаты,

полученные

В.

В.

Соколовским

[42].

Дифференциальные уравнения ха­

рактеристик следующие:

 

 

 

 

 

 

 

d y =

d*tg(coi ±

|i).

(111.62)

 

Вдоль характеристик должно выполняться следующее урав­

нение относительно а:

 

 

 

 

 

 

 

da ±

2а tg ф d<oi = у (dx =F tg ф dу).

(111.63)

Из системы уравнений (II1.63) выводятся известные урав­ нения Кеттера, содержащие радиусы кривизны Ra и /?р характе­ ристик или, что одно и то же, линий скольжения первого и второго семейств. Первое уравнение (III.63), имеющее вид

130