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

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

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

Рнс. I ll.fl. Система цилиндрических коор­ динат для осесимметричной задачи теории предельного равнойесия сыпучей среды

Рис. 111.12. Прямоугольные (а ) и полярные

(б)координаты и соотношение между

углами

 

da + 2a tg <р d«>i =

Y(d* -

tg <p dy) =

ydy

C0S*W------(IFI.64)

 

 

 

 

COS (p

COS

( l - f p )

 

будет

 

 

 

 

 

 

 

 

- j g - +

2a tg q> =

yRa

1 ~ JL) .

 

(Ill 65)

 

door

 

cos cp

 

 

 

Второе уравнение будет таким:

 

 

 

 

 

— - 2a tg <р =

-уЛ р C0S(<0| +

^

 

(Ш.66)

 

dooi

 

cos<p-

 

 

 

Если форма линий скольжения задана

заранее,

т. е., если

Rn и

есть известные функции от

a>i,

то

эти

уравнения

могут быть проинтегрированы.

Исследование случая, когда одно из семейств линий скольже­ ния является окружностями, было рассмотрено Како. Второе семейство при этом представляло собой логарифмические спи­ рали. Другой случай, когда одно из семейств линий скольжения представляет собой логарифмические спирали (второе в этом

131

случае также будет представлять логарифмические спирали или в частнолк случае — пучок прямых линий), был рассмотрен Голушкевичем.

Если принять

Ra = Roe-"*',

(ПШ )

то уравнение (Ш .65) может быть проинтегрировано:

а = а0е2“' •«* + ?Яа

sin

- ц + еа);

COS (J

(III.68)

tg во = 2 tg д>+ п.

Аналогичным путём, интегрируется и уравнение (III.66):

о = ао<Г2"'

- Rt cos 60 Sin (wi + ц - Q0).

(Ш.69)

 

COS ф

 

В. В. Соколовский рекомендует применять эти решения для приближенных вычислений, принимая в простейшем случае п = 0.

5. ЗАДАЧА О ПРЕДЕЛЬНО НАПРЯЖЕННОМ СЫПУЧЕМ КЛИНЕ

Рассмотрим задачу о клине, находящемся в предельно напряженном состоянии. Эта задача впервые была рассмотрена Л. Прандтлем [32] и решается в полярных координатах* г, 6. Дифференциальные уравнения равновесия имеют следующий вид:

даг .

1 дтго

,

=

у cos 0;

 

дг

+

г

 

+

 

 

 

 

(III.70)

дт,о

I

дор

2т,о

= —V sin 0.

}

дг

+

г

<?0

+ г

где 0 — угол, отсчитываемый

от

вертикали

против

часовой

стрелки; у — удельный

вес грунта.

 

 

 

 

 

 

 

Этим уравнениям можно удовлетворить, вводя функцию напряжений Эри F(r, 0). Тогда напряжения выразятся следующим образом:

or =

1

dF

1 d2F .

. 2n .

t • i(\\ n

4 .

*

---- — +

--Г Т 5- + yr (cos'* 0 +

So siir 0) cos 0 — c ctg ф;

 

 

r

or

T . U0

 

 

 

 

OQ=

 

yr (sin2 0 4- io cos2 0) cos 0 c ctg ф;

 

(111.71)

Trt e

^ ( T I ) “ f

(l “ u)eM 9sin 2e-

 

 

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

/ ( о , — <Уд)г + 4 т £ _

■sin ф.

(Ш.72)

о, + «Т + ctg Ф

 

 

132

Уравнения (III.70) и (III.72) составляют замкнутую систему. Поскольку выражения (II 1.71) удовлетворяют полностью уравне­ ниям (III.70), то можно, подставив выражения (III.71) в урав­ нения (III.72), получить одно существенно нелинейное уравнение относительно функции напряжений F\ которую следует найти. Это уравнение будет следующим:

\dF

, \ d2F

d2F

(1 - £ . V cos8cos2e] 2+4[

| ( ^

) +

]_t

2

r dr

V ae2

+

— “

-7fcos0sin20j

 

 

 

 

 

 

 

 

Г 1

dF . 1

d2F , &F , „

, „ ч

Л

Л

12

 

 

\Г 7 Т Л ^ ~ д ¥ + 'д ? '+ (

+ 6о)тг cos 0 +

2c ctg

 

 

 

 

=

sin2 ф

 

 

 

(111.73)

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

Рассмотрим сначала невесомый клин при у = 0. Если вдоль грани клина действует нагрузка, изменяющаяся по степенному закону, то функция напряжений F должна иметь вид

F=*t* + if(9).

(Ш.74)

На основании выражений (III.71) с учетом формула (III.74) получим для невесомой и несвязной среды следующие зависи­ мости для напряжений:

Ог — г* [(л +

2) / + 7" ];

1

 

0о =

(я + 2) (я + 1) г"/;

V

(Ш.75)

ТгО=

—(я +

1) гТ

J

 

Дифференциальное уравнение

(III.73) превратится приу = 0

в обыкновенное однородное уравнение. Положив

f (0) = Сетв и

разрешив уравнение относительно т, получим два корня:

 

m =

±

j / 4 (я +

1)tg2 ф — пг

(III.76)

Решение

(III.26) в

частном

случае при п =

0 переходит

в решение

Прандтля

[32]

(для

которого m = ± 2 t g ср), рас­

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

Компоненты напряжений в исследуемой задаче согласно

выражению (III.71)

имеют вид:

 

 

а, =

Сг" (я +

2 + /я2) в”0 — с ctg <р;

1

 

ае =

Стл (я +

2)(я + 1) е"10 — с ctg q>;

>

(III.77)

Тгэ =

Crnm (я -)- 1) еГ*,

)

 

где m определяется в соответствии с формулой (III.76).

 

 

Угол, составленный направлением нагрузки

(рис.

III. 13) и

нормалью к загруженной грани

 

 

 

6 =

a r c t g | ^ ^ |

 

(III.78)

133

Определенные ло формуле (1J1.76) значения т оказываются действительными при

tg <р > пЦ2 tfn + J).

(111.79)

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

при <р > 19°28'.

 

 

 

 

Компоненты

напряжений (III.77)

можно представить также

в другой записи:

 

 

 

о, =

Сг” £ 1

■+■sin2 ф — ~ cos2

е"19 —- с ctg <р;

'

аа =

С - [ 1

— sin2 ф + ~ cos2

<pj em0 с ctg ф; ’

(iii.ao)

т о = =F

| / 4 { я + i ) t g 2 9 — я 2 в ” 9

 

На основании полученных зависимостей для компонент напряжений (III.77) или (Ш .80) путем несложных преобразова­ ний можно получить уравнение линий скольжения. Напомним, что вдоль линий скольжения удовлетворяется условие

(III.8I)

Это уравнение будет следующим:

/■= Cie9-C,ew.

(Щ.02)

Угол, составляемый линией скольжения и соответствующим радиусом (см, рис. III.13), может быть определен по формуле

134

— arcsin

~

+ 2) ” 2m(n + 0 tg ip

- M f ’

(111.83)

2

 

=F [ж* + (я + 2Л t8 v

 

 

где т берется по зависимостям

(111,75).

 

 

Для первого семейства линий скольжения выбираем в фор­ муле (III.83) верхние знаки и' k = 0, -для второго — нижние знаки и k = 1. Как следует из выражения (III.82), линии сколь­ жения представляются двумя семействами логарифмических спиралей.

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

Для случая линейного увеличения напряжений вдоль радиуса следует принять п = 1. Тогда получаем вместо зависимостей (III.75) следующие:

 

аг = г (3/ + /");

а0 =

6rf; т,0 ------2rf\

(III.84)

Приняв, как и говорилось,

 

 

 

[ =

/(0) =

Сет0,

(111.85)

будем иметь:

^ ___

 

 

 

т =

 

 

 

± V o tg 2 ср — 1;

 

 

 

Or =

2Сг (1 -Ь 4 tg2 ф) ет0 — с ctg ф;

 

о0 =

6Crem0 — c c tg ф;

 

 

(III.86)

т,о =

—2C<cmemQ;

 

 

 

а д =

(2 tg2 Ф 1 ± т tg ф ) c o s 2 ф ctg ф kn/2. >

 

В этом случае оба семейства линий скольжения будут ло­ гарифмическими спиралями. К сожалению, замкнутые решения этой задачи для весомого грунта получить не удается даже в та­ ком простом случае, когда £ о = 1. При п = 1 получаем сущест­ венное упрощение уравнения (III.73): под радикалом имеется только искомая функция f, уравнение является обыкновенным, однако нелинейным и поэтому не имеет решения в элементарных функциях. Оно имеет следующий вид:

VU" —з/)2+ 16 (f'f

=sin ф.

(111.87)

9/ + /" + 2^ cos 0

 

Это уравнение может быть решено только приближенно. Преобразуем уравнение (III.87) следующим образом:

(/" ^ 3+ i6 (П 2 - sin2 Ф(9f + Г + 2V C O S в)2 = 0.

(111.88)

После этого один из возможных вариантов приближенного решения рассматриваемой задачи будет следующим — уравнение (III.88) заменяется:

135

( Г — З/)5 + 16(/')2 — sin* ф (9/ + /" + 2*1 Ce1"»)2.

(Ш.89)

Применяя, как и ранее, подстановку выражения (Ш.85), мы получим обобщение решения для коэффициента т. Вместо

зависимости

(III.89)

получим биквадратное

уравнение

.отно­

сительно т:

 

 

 

 

 

 

 

 

 

 

(m2 — З)2 + 16m2 — sin2 ф (9 +

от2 + 2*if .

(Ш.90)

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

чительных преобразований в следующем виде:

 

 

т 2 = 4 tg2 ф -

5 +

2*1 tg2 ф +

2(1 + tg2 ф) / 4

+ ( 4 + *i)sin2 <|i.

(III.9I)

Совершенно

естественно,

что при

k\ =

0

формула (111.91)

преобразуется в полученную ранее

 

 

 

 

 

 

 

 

т 2 = 8 tg2 ф— 1 =

т 2.

 

 

(III.92)

Решение это, как указывалось, является

действительным

при /и о > 0,

т.

е. при <р^ 19°28'.

В

табл.

HI. 1 приведены

значения коэффициента т в зависимости от k\ для разных

углов внутреннего трения <р, рассчитанные по формуле (III.91).

 

 

 

Т а б л и ц а Ш.1

Ai

 

Значение ±т при <р, град

 

20

30

40

 

0

0,2445

1,2910

2,1524

0,1

0,3367

1,3421

2,2173

0,2

0,4100

1,3922

2,2812

0,3

0,4732

1,4414

2,3442

0,4

0,5300

1,4896

2,4062

0,5

0,5822

1,5371

2,4673

0,6

0,6309

1,5838

2,5275

1,0

0,8029

1,7635

2,7602

1,5

0,9870

1,9746

3,0340

2,0

1,1510

2,1730

3,2914

3,0

1,4422

2,5384

3,7650

4,0

1,7007

2,8701

4,1944

5,0

1,9361

3,1750

4,5888

В

заключение

вернемся

к

системе

уравнений

(III.49) и

(II 1.50) и ее решению (III.51)

и

(IH.52), которое справедливо,

если

определитель

системы

(III.49) и

(III.50)

не

обращается

в нуль:

 

 

 

 

 

 

 

 

Р _1 sin ф sin 2<02

 

2s sin ф cos 2сог|

 

(III.93)

 

 

1 1 — sin фcos 2m

2s sin ф sin 2шг1

 

 

 

 

 

Раскрывая этот определитель, получим

 

 

 

 

D —2s sin ф(sin ф — cos 2т).

 

(Ш.94)

Равенство определителя нулю, поскольку s ф

0, соответствует

условию sin <р == cos 2(0г, откуда

имеем,

что о>г =

± р с- Следова-

136

тельно, решение (III.51) и (III.52) непригодно для случая

<02 =

± |1 с +

пп

(Ш.95)

и этот случай исключается. Эрто

[42]

исследовал детально

решение системы (III.49) и

(III.50)

именно для этого случая,

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

Рассмотрим сначала случай невесомой среды, когда правые части уравнений (III.47) обращаются в нуль. В этом случае определитель системы оказывается тем же, а уравнения (III.49) и (III.50) претерпевают изменения за счет того, что в правой части первого из них будет отсутствовать cos 0, а в правой части второго уравнения будет отсутствовать sin 0. После разрешения этой системы [42] получено:

^d0+ i .

________ cos2 у________ .

2 sin у (cos 2<I)2 — sin у )r .

 

ds

s sin 2(02

Далее вводится

d0

cos 2(I)2 — sin у

 

 

 

vc =

cos2 y

 

1 — 3 sin2 у

------ -— sin у =

-------------2 sin у -

 

2 sin у

 

и тогда из выражения (II 1.96) можно получить:

d0 ^

sin у +

v*

ds

s sin 2m

dci)2

COS 2(1)2 +

vc *

d(l)2

cos 2O)2 + vc

(III.96)

(III.97)

(111.93)

При sin qp =

1 /3 или

q>=19°28'

имеем: vc=

1 по формуле

(Ш .97). Для остальных

случаев

при

<р>19°28'

имеем vc< 1

и решение [42]

получит вид:

 

 

 

 

 

0 _

0О= ™ 2 ± £ Arth(

 

 

 

tg « ) -

(III.99)

 

У I — %

v w

1

+

Vc

'

S =

so V Icos 2G)2 -f- Vcl ,

 

 

 

 

 

где Oo и so — произвольные постоянные интегрировання.

Зависимость (III.99) имеет смысл только, если

tg 0)2 < /

1 + У с

1 + 2 sin у — 3 sin*<p

(III.100)

1 — Vc

2 sin у + 3 sin2 у — 1

Из сопоставления решений уравнений (III.49) и (111.50)/ имеющих вид (III.51) и (III.52), и зависимостей (III.96) получим, что за счет влияния веса к системе (III.96) в первую формулу добавляется

cos9 — sinycos(2(i)2+Q ) 2s sin у (cos 22 — sin у) ’

и во вторую

sin (2о)2 + 0) cos 2(02 — sin у

137

 

DA

OB

Ins м

1

i

*

2

-Ж/2

JLLC

0

+Л/2 Q

PHC. III.14. К решению задачи о предельно

Рис. 111.15. Зависимость между Ins и

напряженном основании

 

 

углом в

 

Если определитель системы

D = О,

то

при равенстве

нулю

также частных определителей Л = 0 система оказывается неопре­ деленной и наблюдается линейная зависимость уравнений между

собой, т. е. решение получается не единственным.

Если'Д ^ О ,

a D = О, то система несовместна. В первом

случае имеет

место параллельность линий, во втором — линии

скрещиваются

в пространстве и не пересекаются.

 

Рассмотрим теперь граничные условия для задачи о несущей способности оснований, находящихся в предельно напряженном состоянии! Граничные условия нужны для решения задачи с помощью численного интегрирования системы (III.51) и (III.52). Интегрирование ведется для простейшего случая отсутствия пригрузки и сцепления в грунте, т. е. для идеально сыпучего

грунта (рис.

III.14).

Установим граничные условия на линии

ОА

при

0 =

—я /2 +

рс= —я /4 +

ф/2 и на линии ОВ при

0 =

я/2.

В

области

ЛоОЛ имеем

простейшее напряженное

состояние по Ренкину [42]. Положительное направление угла 6 показано на рис. III. 14! Кроме того, вводится обозначение:

(III.101)

Граничные условия сведены в табл. III.2, где знак вопроса означает, что эта величина подлежит определению. На рис. III.15 представлена зависимость In s от 0.

Решение системы уравнений (III.51) и (III.52) выполняется численным методом с учетом граничных условий по табл. Ш.2. Для интегрирования формулы записываются следующим образом:

Si sin 2со2, i — sin (2шг, <+

ft) ^ .

 

cos 2(D2, I — sin ф

 

 

cos 9/ — sin <p cos (2(02, / + 0/) — si cos2

де, (in.Ю2)

Л(1)2, i =

 

2si sin ф (cos 2(I)2, i — sin ф)

 

138

причем

е ,+ ' — 0 i + Д0; <Й2. i+ I = 0)2,, + До>2. 5 ,+ | = Si + b S i . ( I I I .103)

 

 

 

 

 

Т а б л и ц а III.2

Параметр

 

Условия

 

Условия на линии OB

 

на линии OA

при 6

о

при 6 = 0

 

 

 

0

-

— Цс) =

л

 

Л

 

 

 

2

 

2

— ( f + f )

0)2

— Цс

d(i>2

 

~d0~

 

s

1

2 sin fic

 

ds

COS |l.

ill]

2 stif Цс

In 5

—In (2 sin |ic)

я

1 ,

л

“ T “ T (A + e)

2~

—sin <psin (Д -f 6) -f- s cos2 ф

1— 3 sin ф

2s sinф£соб(6+ Д)+ sinф]

2 sin ф

 

>

?

1

cos(Д + б)

1

s cos(Д -Ь 6) И- sin ф

s (1 -|- sin ф)

 

sin (Д + б)

 

cos(Д + б) -f- sin ф

 

 

?

?

Расчет ведется шаговым методом, и интегрирование выпол­

няется от 0 =

+ я /2 до 0 = —(я/4 + ф /2). Задаемся несколькими

значениями

so

и проводим расчет так, чтобы «спустившись»

к 0 = —(я/4 +

ф/2) получить s ( 2 sin цо)“ 1 и т = —цс- Слож­

ность такого расчета заключается в накоплении ошибки с ростом числа итераций. В качестве контрольных значений имеются ве­ личины производных на линии ОЛ (см. рис. III.2). Начальные значения производных на ОВ также даны в этой таблице. В случае необходимости расчет может вестись как для. поло­ жительных, так и для отрицательных значений угла б, причем чем меньше шаг Д0, тем меньше будет ошибка, поскольку ошибка с увеличением 0 накапливается. Естественно, что значения s должны быть не менее получаемых из решения задачи для невесомой среды и вычисляемых по формулам (III.99).

139

6. ПРИБЛИЖЕННОЕ РЕШЕНИЕ ПЛОСКОЙ ЗАДАЧИ ТЕОРИИ ПРЕДЕЛЬНОГО РАВНОВЕСИЯ СЫПУЧЕЙ СРЕДЫ

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

д О х .

(frx y _

дх

ду

д°у

(III.104)

__Q

ду

дх

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

/ (ох — оу)2 + 4т%=

(сг* +

(Jj, + 2с ctg ф) sin ф.

(ШЛ 05)

Система уравнений (III.104) и (III.105) существенно нелиней­

ная, и поэтому в замкнутом

виде

может быть решена

только

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

Эта система имеет две новых переменных и состоит из двух уравнений. Принципиальный подход при такой постановке за­ ключается в следующем. Предлагается.подстановка, с помощью которой удается точно удовлетворить уравнениям (III.105). Од­ нако такая подстановка «портит» систему (111.104), превращая ее из линейной в квазилинейную. Эта система решается численно, таким образом условия равновесия удовлетворяются прибли­ женно, правда всегда имеется возможность повысить степень точности их интегрирования численным путем.

Возможен и другой путь решения, использовавшийся нами [32], когда уравнения равновесия (III.104) удовлетворялись точно, а условие предельного напряженного состояния (III.105) — приближенно. Удовлетворить условиям равновесия можно введе­ нием функции напряжений Эри [49]. Однако в этом случае имеется лишь одно уравнение относительно функции Эри F. При использовании этой функции и после подстановки' в условие

(III.105) зависимостей, тождественно

удовлетворяющих

системе

(III. 104).:

 

 

 

 

<ь = 4 4

- + V*;

°у = 4 з - +

=

(Ш.106)

оу

 

олг

дхду

 

получим это существенно нелинейное уравнение.

140