Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
M.L.Zocenko_-_Inzh_geol_Mehan_gruntiv_osnovy_i_fund.pdf
Скачиваний:
661
Добавлен:
05.02.2016
Размер:
6.57 Mб
Скачать

стійкості підпірної стінки приймають найбільше з одержаних значень сили Е.

9. ГРАНИЧНИЙ НАПРУЖЕНИЙ СТАН АНІЗОТРОПНИХ ОСНОВ

9.1. УМОВИ ГРАНИЧНОГО НАПРУЖЕНОГО СТАНУ АНІЗОТРОПНОГО ЗА ОПОРОМ ЗРУШЕННЮ ҐРУНТУ І РОЗРАХУНКОВА МОДЕЛЬ

Загальні питання, пов’язані з анізотропією ґрунтів, були детально розглянуті раніше у п. 4.10. Нижче розглядаються основи теорії граничного напруженого стану анізотропних за опором зрушення ґрунтів.

Під характеристикою анізотропії зрушення будемо розуміти функцію τ(β, x, y) граничного опору зрушенню по елементарній площадці, зорієнтованій під кутомβ у довільній точці x, y, що належить ґрунтовій напівплощині. Тоді зручною формою представлення τ(β, x, y) є годографи кута внутрішнього тертя φ(β, x, y) і питомого зчеплення с(β, x, y). Під однорідною анізотропною напівплощиною будемо розуміти напівплощину, яка задовольняє умови:

φ(β,x,y)=φ(β); с(β,x,y)=с(β). (9.1)

Для реальних ґрунтових умов наслідком умови (9.1) є також відповідність залежнос-

тям:

 

φ(β)=φ(π+β); с(β)=с(π+β).

(9.2)

У випадку регулярної площинно-паралельної шаруватості, крім задоволення умов (9.1), (9.2), годографи φ(β) та с(β) як мінімум мають дві ортогональні осі симетрії.

Розглянемо елементарну площадку АВ, що проходить через довільну точку x, y, що належить ґрунтовій напівплощині, анізотропні властивості якої відповідають умовам (9.1), (9.2) і яка зорієнтована під кутом β до вибраної осі відліку. Показники анізотропії за опором ґрунту зрушенню зафіксовані за допомогою φ(β) і с(β). Тоді умову граничного напруженого

стану в головних напругах за аналогією з ізотропним ґрунтом запишемо у вигляді

 

σ1 σ2 =(σ1 +σ2 )sin[ ϕ( β )] ,

(9.3)

де σ1 =σ1 +σc ( β ) і σ2 =σ2 +σc ( β ) – приведені головні напруги, σc ( β ) = c( β )ctg[ ϕ( β )]

– тиск зв’язності.

Повернемо площадку АВ на кут β. Для ізотропного ґрунту ця операція не приведе до необхідності зміни відношення

(σ1 σ2 ) /(σ1 +σ2 ) в умові (9.3). Тобто у цьому випадку, скориставшись поняттям про орієнтова-

ний круг Кулона-Мора, головний діаметр якого суміщений із нормаллю до площадки АВ, при обертанні останньої прийдемо до відомої системи характеристичних кругів (рис. 9.1), котра є графічним аналогом умови (9.3) для ізотропного ґрунту при повороті площадки АВ. Для анізотропного ґрунту виконання умови (9.3) при повороті площадки потребує іншого напруженого стану, тому умова граничного напруженого стану не буде інваріантністю відносно β. Побудова графічного аналога умови (9.3) можлива, якщо виходити з наступного. Неважко побачити, що ліва частина

умови у записі (9.3) σ1 σ2 =σ1 σ2 = d не залежить від характеристик середовища, тому для

О2

 

 

А

β В

φ φ

 

b

(х, у)

R2

α

 

 

О1

β

 

 

 

a

R1

1

2

3

 

0

 

d/2

 

 

 

Рис. 9.1. Графічний аналог умови (9.56)

для ізотропного ґрунту при обертанні площадки АВ: 1-3 – характеристичні площадки

205

а

б

 

75°

60°

 

 

 

φ°=15, 20, 25, 30,

90°

 

 

45°

135°120°105°

φ(β)

φ(β)

30°

 

 

 

150°

 

 

R2(α)

15°

 

 

 

 

165°

a

O1

 

 

 

R1(α)

 

 

 

180°

 

 

 

α=0°, β=0°

d/2

 

b β

 

 

2

 

 

 

195°

3

 

 

 

 

 

 

 

 

 

210°

 

 

 

 

 

225°

255°

270°

 

 

 

240°

 

 

 

 

Рис. 9.2. До прикладу побудови графічного аналога умови (9.56):

а – вихідний годограф кута внутрішнього тертя (див. табл. 9.1);

б – згідно з результатами розрахунку (див. табл. 9.2)

зручності подальших побудов доцільно зафіксувати d. Як і ран іше будемо використовувати орієнтований круг Кулона-Мора. Тоді при обертанні площадки АВ (рис. 9.2, б) точка О2 буде описувати огинаючу характеристичного круга площадок. На рис. 9.2, б наведений результат цієї побудови для φ(β), що зображений на цьому ж рисунку. Розрахунок R1 та R2, які необхідні для побудови, виходячи з рис. 9.2, б, можна виконати за допомогою елементарних залежностей

R1(α ) = R1

( β +π / 2 ) = 0,5d sin[ ϕ( b )] ;

(9.4)

R2 (α ) = R2 ( β +π / 2 ) = 0,5d / sin[ ϕ( b )] .

 

Отримана система фігур є графічним аналогом умови (9.3) для анізотропного ґрунту. Для того, щоб ним скористатися при довільному куті β1 орієнтації площадки АВ, необхідно виконати побудову у такій послідовності:

з центру О провести лінію ОО2 під кутом α=β+π/2 відносно осі відліку до перетину з огинаючою кута полюсів;

у точці О1 провести хорду аb перпендикулярно до лінії ОО2, яка буде паралельною заданій площадці АВ;

завершують побудову дотичні О2a та О2b, які показують граничний кут відхилення наведеної напруги по АВ.

Круг радіусом d/2, що називають кругом вершин, є геометричною подобою граничного круга Кулона-Мора. У цьому можна впевнитись, виходячи із залежності (9.4).

Можна також запропонувати інші форми графічного аналога умови (9.3), але наведений є найбільш зручним для подальшого аналізу. Відзначимо також, що при φ(β)=const R1(α) та R2(α) перетворюються у постійні величини. Отже, розглянута система фігур вироджується у систему характеристичних кругів для ізотропного ґрунту, чого і варто було чекати.

Приклад побудови графічного аналога умови (9.3). Основа задовольняє умови (9.1), (9.2). Напрямок осі відліку кутів β – горизонтальний. Годограф φ(β) (рис. 9.2, а) заданий значеннями, що наведені у табл. 9.1.

`

206

Таблиця 9.1. Вихідні дані для побудови годографа

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

β°

0

15

30

45

60

75

90

 

105

120

135

150

165

180

φ°(β)

20

32

30

23

20

17

15

 

14

14

14

15

17

20

 

Розраховуємо коефіцієнт анізотропії μφ

 

 

 

 

 

 

 

 

 

 

 

 

 

µϕ =1/πρmax2

,

 

 

 

 

(9.5)

де = 1 2πρ2 ( β )dβ , а ρ(β) – функція φ(β). 2 0

Для прийнятої залежності φ(β) розраховане за виразом (9.5) значення μφ=0,398, що підтверджує дані про суттєву анізотропію основи.

Приймаємо d=5 (масштаб побудови може бути довільним). Для кожного β із наведених у табл. 9.1 визначаємо R1(α) та R2(α). З умови (9.2) видно, що R1(α+π)=R1(α) й R2(α+π)=R2(α). Це дозволяє скоротити обчислення вдвічі. Результати обчислень наведені у табл. 9.2.

Таблиця 9.2. Результати обчислень для побудови рис. 9.2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

β°

0

15

30

45

60

75

90

105

120

135

150

165

180

α°

90

105

120

135

150

165

180

195

210

225

240

255

270

R1(α)

0,9

1,3

1,2

1,0

0,9

0,7

0,64

0,6

0,6

0,6

0,64

0,7

0,9

R2(α)

7,3

4,7

5,0

6,1

7,3

8,5

9,65

10,4

10,4

10,4

9,65

8,5

7,3

За результатами розрахунків виконуємо побудову (див. рис. 9.2, б). Варто відмітити, що при формулюванні умови граничного напруженого стану у вигляді (9.3) неявним чином використана теорема А. Како відносно врахування зчеплення в основі. Ця ж обставина дала можливість без особливих ускладнень побудувати графічний аналог умови (9.3) для анізотропної напівплощини.

Окрім умови граничного напруженого стану у вигляді (9.3), у подальшому будемо використовувати також його іншу форму

τ( β ) =σβ tg[ ϕ( β )] ,

(9.6)

де σβ =σβ +σc ( β ) – наведена нормальна напруга на площадці ковзання.

Розрахункова модель невагомої ґрунтової напівплощини в теорії граничного напруженого стану, як відомо, містить два рівняння рівноваги у нескінченно малих і умову (9.3), відображену у компонентах напруг:

σ x +

τxy

= 0 ;

 

 

τxy

+

σ y

= 0 ;

 

 

 

 

 

x

y

 

 

x

y

(9.7)

(σ x σ y )2 + 4τxy2

=(σ x +σ y )2 sin2 [ ϕ( β )] ,

 

де σ x =σ

x +σc ( β ) та σ y =σ

y +σc ( β )

– наведені компо-

ненти напруг.

Скористатись рівнянням (9.7) у практичних розрахунках, як правило, дуже складно. Розглянемо дещо іншщий підхід до побудови рівнянь рівноваги, що входять у систему (9.7). Нехай напружений стан у точці заданий компонентами напруг. Навкруги точки побудуємо елементарний прямокутник (рис. 9.3); через точку А проведемо площадку АВ, яка зорієнтована відносно осі y під кутом β, і будемо вважати АВ площадкою ковзання. Для зручності також будемо вважати довжину АВ одиничною. Тоді з умови рівноваги ABD отримаємо три рівняння:

 

σx

 

 

 

A

τxy

 

D

y

cosβ

 

β

 

 

 

 

 

 

1

sinβ

σy

 

 

τβ

 

 

 

 

 

B

τxy

 

 

 

 

E

 

 

C

 

x

 

 

 

 

Рис. 9.3. До аналізу ро-

 

зрахункової моделі

 

207

σβ

=τxy sin 2β σ y sin2 β σ x cos2 β ;

(9.8)

τβ

= 1 (σ y σ x )sin 2β τxy cos 2β ;

(9.9)

 

2

 

τβ sin 2β +σβ cos 2β σ x cos2 β +σ y sin2 β = 0 .

(9.10)

Підставляючи вирази (9.8), (9.9) у рівняння (9.10), у підсумку отримаємо співвідно-

шення

τxy sin 2β σ x +σ y = −(σ x +σ y )cos 2β ,

(9.11)

яке разом з умовою граничного стану, записаного у рівняннях (9.7), при заданих компонентах напруг у точці шляхом послідовних наближень дозволяє знайти напрямок дотичної до кривої ковзання у даній точці. Одержана система рівнянь у дійсних компонентах напруг має такий вигляд:

τxy sin 2β =σ

x

+σ

y

= −[ σ

x +σ

y + 2σc ( β )] cos 2β ;

(9.12)

(σ

x σ

y )2 + 4τxy2 = [ σ

x σ

y + 2σc ( β )] 2 sin2

[ ϕ( β )]

 

Якщо у системі (9.12) прийняти σ

x =σ1 >σ

y =σ

2 , то τxy = 0 , і в результаті отрима є-

мо

= [ σ1

 

 

 

 

+ 2σc ( β )с] cos 2β ;

 

 

 

 

σ1 σ

2

+σ

2

 

(9.13)

 

 

σ1 σ

 

 

= [ σ1

+σ

 

+ 2σc ( β )с] sin[ ϕ( β )] ,

 

 

2

2

 

або в підсумку ± β =π / 4 +ϕ( β ) / 2 . Знаки перед β показують, що реалізуються дві симет-

ричні площадки ковзання. Однак це справедливо лише у тому випадку, коли годографи φ(β) та с(β) мають дві ортогональні осі симетрії, одна з яких збігається з напрямом σ1 або σ2. У протилежному випадку друга площадка, як і перша, відшукується послідовним наближенням та дотична до поверхні ковзання недзеркальна відносно осі y. Таким чином, використання системи (9.12) також приводить до громіздких розрахунків при визначенні кінематичної картини руйнування анізотропної основи.

9.2. ВИРІШЕННЯ ЗАДАЧ ДЛЯ АНІЗОТРОПНОЇ ЗА ОПОРОМ ЗРУШЕННЮ ОСНОВИ

Основна задача

Основною вважають задачу відшукування напрямів характеристик через кінці заданої елементарної площадки, завантаженої наведеною напругою, відхиленою від нормалі до площадки під довільним кутом δ. При вирішенні цієї задачі будемо оперувати анізотропною напівплощиною із заданими φ(β) та с(β), які відповідають умовам (9.1), (9.2).

Розглянемо граничну рівновагу невагомої елементарної призми ABD, вважаючи, що по площадці АВ розподілені наведені напруги, рівнодіюча яких F відхилена від нормалі до площадки під кутом δ. Окрім того, будемо вважати AD i BD площадками ковзання. Відносно осі відліку AB, AD та BD орієнтовані відповідно під кутами β1, β2, β3 (рис. 9.4). Тоді згідно з годографом φ(β) відповідні кути внутрішнього тертя будуть φ(β1), φ(β2) та φ(β3).

Побудуємо навкруги призми коло з центром у точці О, з котрої проведемо перпендикуляр до площадки АВ. На останньому знайдемо точку, де перетинаються дотичні до круга, проведені через кінці площадки АВ. Оскільки масштаб побудови довільний, будемо вважати, що круг 3 відповідає кругу вершин, тобто його радіус дорівнює d/2 відповідно позначенням п.9.1, і тоді виконана побудова збігається з наведеною у попередньому розділі.

Таким чином, якщо призма ABD знаходиться у граничному напруженому стані, кут ρ=φ(β1). У протилежному випадку круг вершин, що є геометричною подобою кола Кулона– Мора, не буде граничним, що суперечить вихідній передумові. Звідси виходить, що точка О1 розташована на огинаючій характеристичних кіл площадок, а О2 – на огинаючій кіл полюсів. Крім того, якщо ABD відповідає мінімальному напруженому станові (а – на рис. 9.4), то кут

`

208

f=π/2-φ(β1), якщо максимальному (б – на рис. 9.4), то f=π/2+φ(β1), а ν та μ можна визначити за допомогою залежностей

π

±ϕ( β1 ) +δ1 ± arcsin[sinδ1

 

;

ν = 0,5

2

/ sinϕ( β1 )]

 

 

 

(9.14)

π

±ϕ( β1 ) δ1 ± arcsin[sinδ1

 

,

µ = 0,5

2

/ sinϕ( β1 )]

 

 

 

 

де верхні знаки відповідають напруженому станові а, нижні – стану б.

Повертаючись до розгляду граничної рівноваги ABD, відзначимо, що оскільки AD й

AB – площадки ковзання, то F1 i F2 відхилені від нормалей відповідно під кутами φ(β2) та φ(β3) і проходять через середини відрізків AD та BD із причини невагомості основи.

З умови рівноваги ABD виходить, що пряма F повинна проходити через точку е перетину F1 i F2. Оскільки на вид годографу φ(β) обмеження не введені, то в загальному випадку δ≠δ1. Звідси виходить, що для анізотропного ґрунту принцип сполученості площадок ковзання несправедливий. Ця обставина значно ускладнює вирішенння основної задачі у зіставленні з подібним для ізотропного ґрунту. Дійсно, заданому δ може відповідати декілька значень δ1, при яких F1 i F2, відхилені від нормалей під відповідними кутами внутрішнього тертя, будуть перетинатися в одній точці з напрямом F. Отже при розв’язанні практичних задач необ-

хідно використовувати екстремальний принцип механіки.

 

 

 

 

 

 

Побудова, виконана на рис. 9.4, по суті є графічним вирішенням основної задачі для

анізотропного за опором зрушенню ґрунту. Наведемо послідовність процедур:

 

 

у довільному масштабі

 

 

 

 

1 O2

 

 

будують площадку АВ, яка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

зорієнтована під заданим ку-

 

 

 

 

 

F

 

 

том β1

до осі відліку вихідного

 

 

 

ρ δ1

 

 

 

 

 

δ

 

 

годографа φ(β);

 

 

 

 

 

 

 

із точки перетину про-

 

 

 

R2

 

 

 

меня з точки А під кутом φ(β1)

A

 

 

 

 

 

 

 

2

 

 

 

 

із нормаллю до площадки АВ,

ν

a

 

 

B

 

тобто з точки О, будуємо коло

 

β2

 

O1

β1

 

 

μ

(радіус позначаємо d/2), що

 

 

 

 

 

 

 

 

 

β3

 

проходить через точки А і В;

b

 

 

 

 

R1

 

 

на колі вибираємо до-

 

 

 

 

 

 

 

 

 

 

e

 

 

 

вільну точку D і з’єднуємо її з

 

 

 

 

 

 

 

кінцями відрізку АВ;

 

φ(β2)

 

а)

 

O

 

 

і β3

за кутами орієнтації β2

 

 

 

 

3

з

годографа

φ(β),

F1

 

 

F2

 

φ(β3)

d/2

 

з’ясувавши φ(β2) й φ(β3), нано-

 

 

 

 

симо напрямки F1 та F2, які

f

 

 

δ

 

 

 

прикладені у серединах відпо-

 

 

 

 

 

 

 

відних відрізків;

 

D

 

 

 

Fe

 

 

винні

лінії дії F, F1, F2 по-

 

ν

 

O1

B1

 

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

одній

A1

 

β2

б)

2

μ

 

β1

точці, у іншому випадку пере-

 

β3

ходимо до наступної точки D

 

 

 

 

 

f1

 

 

 

 

 

 

 

 

 

 

на колі 3, тобто шляхом пос-

 

F1

 

 

 

D1

F2

 

лідовного

набору точок

D у

 

φ(β2)

 

δ1

 

 

підсумку приходимо до мож-

 

 

φ(β3)

 

 

 

 

 

 

ливих вирішень, серед яких

 

 

 

1

O2

 

 

 

залежно від конкретної задачі

 

 

 

 

 

 

відшукують екстремальне.

Рис. 9.4. Графічне вирішення основної задачі для анізо-

 

За

необхідності

вирі-

тропного за опором зрушенню ґрунту

 

 

209

шення основної задачі для площадок АВ, зорієнтованих під різними кутами до осі відліку, спочатку доцільно побудувати графічний аналог умови граничного напруженого стану відповідно до послідовності процедур, наведених у п. 9.1. Тоді послідовний перебір точок для задоволення умов рівноваги можна проводити на кругу 3 рис. 9.2.

Графічне вирішення, як видно з наведеного алгоритму, достатньо громіздке і для досягнення необхідної точності потребує спеціальних навичок. Тому розглянемо аналітичне вирішення основної задачі, зберігаючи наведені вище позначення.

Із геометричних співвідношень, наведених на рис. 9.4, можна визначити, що відрізок О1е дорівнює у відносних одиницях

 

sinν sin[ ϕ( β1 ) ϕ( β2 )] sinη = sin µ sin[ ϕ( β1 ) ϕ( β3 )] sinη1 ,

(9.15)

де

η =ν ±ϕ( β ) δ ;

 

η= µ ±ϕ( β3 ) δ , (9.16)

аν і μ відповідають залежностям (9.14). Залежність (9.15) зв’язує δ1 та δ2, тобто сумісно із залежностями (9.14) дає розв’язання основної задачі. Розглянемо окремі випадки отриманого вирішення.

1.φ(β)=const. Вираз (9.15) зліва і справа перетворюється в нуль, тобто відрізок О1е вироджується у точку, розташовану в точці прикладання F. Таким чином, принцип суміщеності площадок ковзання буде справедливим, і вирішення, що визначається залежностями (9.14), вироджується, як і належало чекати, у вирішення основної задачі для ізотропного ґрунту.

2.δ1=φ(β1). Із залежностей (9.14) отримаємо μ=0. У цьому випадку вираз (9.15) задовольняється при η=0. Тоді з рівнянь (9.16), оскільки β31, виходить, що δ=δ1, тобто по площадці АВ реалізується площинне зрушення, що відповідає фізичному розумінню. До аналогічного висновку приходимо при δ1=-φ(β1).

3.φ(β2)=φ(β1)≠φ(β3) або φ(β3)=φ(β1)≠η(β2). Тоді з виразів (9.15), (9.16) виходить, що у першому випадку δ=ν±φ(β2) і δ=μ±φ(β3) – у другому.

4.φ(β2)=φ(β3)≠φ(β1). Вираз (9.15) набуває вигляду1 2

sinν / sin µ = sin[ν ±ϕ( β2 ) δ ] / sin[ µ ±ϕ( β2 ) +δ ] .

Звідси при δ=0 отримаємо ν=μ, що відповідає фізичному розумінню. Таким чином, окремі випадки свідчать, що отримане розв’язання основної задачі узагальнює відоме вирішення для ізотропного ґрунту і містить його в собі як окремий випадок. Із виразів (9.15), (9.16) шляхом елементарних перетворень можна перейти до форми вирішення, більш зручної для практичних розрахунків.

tgδ = (sinε k sinψ ) /(cosε + k cosψ ),

(9.17)

де

 

 

sinν sin[ ϕ( β1 ) ϕ( β2 )]

ε =ν ±ϕ( β2

); ψ = µ ±ϕ( β3

); k =

sin µ sin[ ϕ( β1 ) ϕ( β3

)]

 

 

 

При вирішенні конкретних задач окрім кінематичної частини, необхідно вирішувати статичну частину, найважливішою ланкою якої є визначення напрямку дії і значень активних сил. При вирішенні кінематичної частини розглядалась невагома призма АВD. Згідно із загальноприйнятим, впливом власної ваги на кінематичну картину можна знехтувати. Статична частина задачі відрізняється від вирішення для ізотропного ґрунту урахуванням зчеплення. Покажемо це.

При визначенні напрямку F (див. рис. 9.4) ураховувався тиск зв’язності σс( β1 ) = c( β1 ) / tg[ ϕ( β1 )] . Водночас прийнято припущення, що по площадках AD та BD

діють також σc ( β1 ). Найбільш зручно у практичних задачах вибирати у якості базового зна-

чення тиск зв’язності, тобто відносно цього значення вести врахування зчеплення по усіх площадках ковзання. Тоді, розглядаючи ∆σс( βi ) =σс( βi ) σс( β1 ) як навантаження на і-й

поверхні ковзання, що коригує прийняте раніше рівним σс( β1 ), залишається її рівнодіючу

`

210

Q( βi ) = ∆σс( βi ) i , де i – довжина площадки і, додати до власної ваги G призми ABD. При цьому слід користуватися прави-

лом: від’ємні результуючі Q( βi ) спрямо-

вуються по нормалі від площадки ковзання в середину ABD, позитивні – у протилежному напрямкові. На рис. 9.5 показано приклад урахування зчеплення для анізотропного ґрунту і його вплив на величину й напрям активної сили G (рис. 9.5, б).

Вирішення основної задачі є найважливішою складовою ланкою, необхідною для розгляду прикладних задач теорії граничного напруженого стану анізотропних ґрунтів, які будуть наведені у п. 9.3.

а

F σg

б

δσ

 

A

 

 

B

σc2)

G

G1

 

 

 

Q3

 

Δσc3

 

 

 

 

 

 

 

 

 

 

 

 

 

φ(β2)

G

 

 

 

 

 

 

 

 

φ(β3)

 

 

F1

 

 

Q2

 

Δσc2 D

 

F2

 

 

Q3

Рис. 9.5. До врахування зчеплення, заданого

σс(β) при вирішенні основної задачі

Загальна задача Ренкіна

Припустимо, що до плоскої зорієнтованої до горизонту під кутом β1 поверхні напівнескінченного ґрунтового масиву прикладене навантаження інтенсивністю σ0 із кутом нахилу δ′ відносно нормалі до поверхні. Ґрунтове середовище характеризується φ(β) і с(β), осі відліку β у яких збігаються з горизонтом, а також питомою вагою γ. Поставимо задачу про відшукування поверхонь ковзання у масиві.

Оскільки у даному випадкові напруги не залежать від координати x (рис. 9.6), рівняння рівноваги у позначеннях (9.7) набувають вигляду:

dτxz

= γ sin β1 ;

dσ z

= γ cos β1 ,

(9.18)

dz

 

dz

 

 

Після інтегрування і врахування граничних умов отримаємо

б

а δ0

σc(β)

σ0

 

δ

γzcosβ1

π/2

x

z

β1=165°

σ

δ0 σc1)

 

 

δ′σ0 x

z

z=zi

Рис. 9.6. Задача Ренкіна:

а – визначення; б – сітка ліній ковзання, φ(β) і с(β) за рис. 4. , d0=20°, γ=20 кН/м3

211

 

 

 

 

 

 

τxy =γ z sin β1 +σ0 sinδ ;

 

 

 

 

 

 

 

 

 

 

 

 

 

σ z =γ z cos β1 +σ0 cosδ ′+σc ( β1 ) .

 

 

 

 

 

 

 

У підсумку для довільної точки z одержимо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tgδ

 

sinδ ′+ λ1 sin β1

,

 

 

 

 

(9.19)

 

 

 

 

 

 

= cosδ ′+ λ1 cos β1 + λ2

 

 

 

 

де λ1 =γ z / σ0

та λ2 =σc ( β1 ) /σ0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Прирівнявши вирази (9.19) і (9.17),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sinδ ′+ λ1 sin β1

=

sinε k sinψ

,

 

 

 

(9.20)

 

 

 

 

 

cosδ ′+ λ1 cos β1 +

λ2

cosε + k cosψ

 

 

 

 

 

 

одержимо залежність, за допомогою якої можна відшукати напрямки поверхонь ковзання у

будь-якій точці ґрунтового напівпростору. На рис. 9.6, б наведений приклад побудови ліній

ковзання у загальній задачі Ренкіна за допомогою залежності (9.20).

 

 

 

 

 

 

Особливістю вирішення у зіставленні з ізотропним ґрунтом є наявність зони неодно-

значних вирішень, котра легко встановлюється за допомогою залежності δ÷δ1

для площадки

із заданою орієнтацією.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Розглянемо визначення положення площин zi, zj, які відповідають зоні неоднозначних

вирішень δi≤δ≤δj. З виразу (9.19) видно, що

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z = σ

 

 

λ2 + cosδ

)sinδ

 

 

 

 

(9.21)

 

 

 

 

 

 

0 tgδ(

 

.

 

 

 

 

 

 

 

 

 

 

γ

 

sin β1 tgδ cos β1

 

 

 

 

 

 

 

 

Підставляючи у вираз (9.21) замість δ послідовно δi

та δj, у підсумку отримаємо відпо-

відно zi й zj, що визначають межі зони.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Зрозуміло, що для невагомого напівпростору права частина залежності (9.20) зали-

шиться постійною, оскільки ліва частина

tgδλ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= const

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1+ cos2δ

 

 

 

 

 

 

 

 

 

 

 

незмінна для довільної площини z. Тому поверхні ковзання вироджуються у площини, на-

прямки яких визначені аргументами μ і ν в ε і ψ, k із пояснень до виразу (9.17).

 

 

 

β°

 

а

 

 

 

 

 

 

 

 

Аналогічний

висновок

можна

зробити

 

 

 

 

 

 

 

 

при λ12=0. У цьому випадку δ´=δ. Крім того,

200

 

β

β-x

 

 

 

 

оскільки

 

 

залежність

(9.20)

виходить

з

 

 

 

 

 

розв’язання основної задачі, розглянуті окремі

 

 

 

А

 

 

 

 

 

150

 

 

 

 

 

 

 

випадки необхідно доповнити наведеними при

 

βn

βn-nа1

 

 

 

 

 

аналізі залежності (9.15).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

Е

 

 

 

 

 

 

 

 

Звідси можна зробити висновок, що

100

φ(βn)

 

 

 

 

 

отримане вирішення задачі Ренкіна для анізот-

 

 

 

 

 

 

 

 

 

ропного за опором зрушенню ґрунту узагальнює

50

 

 

 

 

 

 

 

 

аналогічне

для ізотропного

масиву

або

при

 

б

 

 

 

 

 

 

φ(β)=const залежність (9.20) вироджується у за-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

лежності (9.14), де φ(β1)=φ і δ1.

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

βn=200°

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Зона Прандтля

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-50

 

 

 

 

 

 

 

 

 

 

Вирішення задачі для проміжної зони є

 

 

 

 

 

 

 

 

складовою частиною діагностики несучої здат-

 

 

 

 

 

 

 

 

 

-100

 

 

 

 

 

 

 

 

ності основ і бічного тиску для анізотропних за

 

45

90

125

180

225

β°n

 

опором зрушенню ґрунтів. Тому її розгляд необ-

 

 

 

хідний

для

отримання

загального розв’язання

Рис. 9.7. До побудови поверхні ковзання

вказаних задач.

 

 

 

 

 

 

другої сім’ї в зоні Прандтля:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Не визначаючи спочатку граничні умови,

а – схема до виводу основного співвідно-

 

 

шення; б

– допоміжна залежність β-βn

 

розглянемо

ґрунтовий

невагомий

напівпростір,

`

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

212

який знаходиться у граничному напруженому стані Прандтля. Це означає, що існує два сімейства поверхонь ковзання, перше з яких являє собою пучок прямих (площин), що проходять через точку А, котру звичайно називають полюсом. Як і раніше, будемо вважати, що анізотропія задана годографами φ(β) та с(β) з осями відліку, фіксованими відносно ґрунтового масиву. Поставимо задачу відшукування поверхонь ковзання другого сімейства. Для цього звернемося до рис. 9.7, а.

Нехай довільна точка Е розташована на довільній лінії АЕ (яка належить першому сімейству поверхонь ковзання), що зорієнтована відносно осі відліку φ(β) під деяким кутом β. Будемо вважати, що аа1 є дотичною до поверхні ковзання другого сімейства, яка проходить через точку Е. При цьому кут орієнтації дотичної аа1 до осі відліку складає βn. Тоді, оскільки наведена напруга у точці Е з умови граничного напруженого стану відхилена від нормалі до аа1 під кутом φ(βn), а з умови належності до зони Прандтля спрямована вздовж лінії АЕ (див. рис. 9.7, а), можна записати таке рівняння:

βr β4 β2

A

 

 

r+dr

r

 

 

F1 φ(β2)

k

φ(β4) F3

βn

D φ(βn)

E F4

Рис. 9.8. Елемент проміжної зони

(до побудови поверхні ковзання)

β = βn ϕ( βn ) π / 2 .

За допомогою рівняння (9.22) можна відшукати напрям дотичної до поверхні ковзання другого сімейства у будь-якій точці, що належить до проміжної зони. Для практичних розрахунків перед побудовою поверхонь ковзання другого сімейства доцільно для заданого годографа φ(β) накреслити залежність β-βn. На рис. 9.7, б зображено цю залежність для годографа, наведеного на рис.4.???. Пунктиром показано порядок визначення βn за заданим β.

У загальному випадку рівняння (9.22) не дає однозначного вирішення у проміжній зоні. При цьому інтервали β, що відповідають неоднозначним вирішенням, легко фіксуються на графіках, аналогічних рис. 9.7, б. На цих інтервалах із збільшенням β значення βn зменшуються або залишаються постійними. Виходячи з екстремальних принципів механіки, в цих випадках із ряду можливих реалізацій кінематичної картини вибирається найбільш невигідна.

Маючи рівняння (9.22), можна перейти до побудови загальної поверхні ковзання, що має важливе значення у практичних задачах. Для цього розглянемо невагомий елемент (рис. 9.8) проміжної зони з нескінченно малим кутом розкриття , сторона AD якого обумовлена і дорівнює r. Дотична до поверхні ковзання у точці D зорієнтована під кутом βn. Будемо вважати, що поверхня ED є поверхнею ковзання, а AE=r+dr. З’єднаємо точки D та Е площиною. У трикутнику ADE кут при вершині Е залишається невідомим, але його можна обумовити через нескінченно малий приріст таким чином:

AED = 90 ϕ( βn ) dϕ . Тоді з геометричних співвідношень виходить:

 

r + dr

 

 

 

=

 

r

 

 

.

sin[ π

dβ +ϕ(

β

n

) + dϕ ]

sin[ π

ϕ( β

n

) dϕ ]

2

 

 

 

2

 

 

 

Після елементарних перетворень рівняння (9.23) запишемо у вигляді

 

r + dr

=

cos[ ϕ( βn ) dβ ] cos dϕ sin[ ϕ( βn ) dβ ] sin dϕ

 

r

 

=

 

.

 

cos[ ϕ( βn )] cos dϕ sin[ ϕ( βn )] sin dϕ

 

(9.23)

(9.24)

Нехтуючи у знаменниках рівняння (9.24) лівими частинами як нескінченно малими порівняно з правими і скорочуючи на cosdφ, в підсумку отримаємо:

213

r + dr

 

r

 

 

=

 

.

(9.25)

cos[ ϕ( βn )] cos dβ sin[ ϕ( βn )] sin dβ

cos[ ϕ( βn )]

Приймаючи у рівнянні (9.25) cosdβ=1 та sindβ=dβ, після інтегрування і знаходження постійної інтегрування отримаємо вираз

 

 

 

 

 

βm

 

 

 

 

 

 

 

 

 

 

 

,

(9.26)

r( β

 

) = r( β

 

 

tg[ ϕ( β

 

m

k

)exp

β

n

)] dβ

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

який дозволяє виконати побудову загальної поверхні ковзання другого сімейства. Неважко побачити, що при φ(βn)=const рівняння (9.26) вироджується у відоме для ізотропного ґрунту в зоні Прандтля. Оскільки на вид годографа φ(β) обмежень не накладають і він задається табличною функцією, визначення інтеграла у рівнянні (9.26) виконується чисельно. У загальному випадку вирішення на окремих елементах зони може бути неоднозначним, тому інтегрування слід виконувати для найбільш невигідної кінематичної реалізації. При цьому попередньо за заданими граничними кутами β розкриття зони з’ясовують значення βk, βm, які і будуть межами інтегрування у рівнянні (9.26). На рис. 9.9 наведений приклад побудови загальної поверхні ковзання, який виконаний на основі годографа, зображеного на рис. 4.??. За допомогою побудови, наведеної на рис. 9.7, б, замість рівняння (9.26) у практичних розрахунках можна використати наближений прийом, що полягає у заміні поверхні ковзання ламаною, напрямки елементів якої виконуються для пучка першого сімейства поверхонь ковзання за допомогою побудов, аналогічних рис. 9.7, б. При цьому похибка не перевищує 3÷5%. Зі с- тавлення із загальними поверхнями ковзання для ізотропного ґрунту, поданими на рис. 9.9, приводить до висновку про суттєву відмінність результатів.

Для практичних задач необхідно мати напрям реакції F4, результуючої реактивні напруги по поверхні ковзання ED. Для розв’язання цього питання розглянемо рівновагу невагомого елемента, зображеного на рис. 9.8, замінивши кінцевим значенням кут dβ=β24. Крім того, обумовимо r(β2) та r(β4) за допомогою виразу (9.26) або графічного засобу, викладеного вище. За умовами задачі елемент ADE знаходиться у граничному напруженому стані Прандтля. Тому наведені напруги по площадках ковзання AD й AE і їх результуючі відхилення від нормалей відповідно під кутами φ(β2) та φ(β4). Крім того, F1 і F3, з причин невагомості ADE, прикладені у точках, розташованих від полюса А відповідно на відстанях r(β2)/2 та r(β4)/2. Продовжуючи лінії дії F1 і F3, отримаємо точку їх перетину К. З умови рівноваги ADE лінія дії F4 повинна включити точку К і полюс А, виходячи з виду напруженого стану. Ця обставина дозволяє з геометричних співвідношень отримати напрямок F4, що визначається кутом βr, значення якого можна визначити за формулою

 

 

 

2

 

 

170°

1

 

 

 

 

 

3

 

r(β)

 

 

 

 

 

 

 

 

160°

 

 

 

 

 

150°

 

 

 

 

 

140°

130°

120°

110°

1090°

Рис. 9.9. Поверхня ковзання:

1 – для анізотропного за опором зрушенню грунту в гранично напруженому стані Прандтля; 2 – відповідає φ=minφ(β)=const; 3 – φ=maxφ(β)=const

`

214

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]