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

668

.pdf
Скачиваний:
6
Добавлен:
06.12.2022
Размер:
3.91 Mб
Скачать
sin(ε − δ + θ − ϕ)

Из полученной системы находим силу давления E при данном положении линии скольжения:

E = Pcos(θ − ϕ − β) + G cos(θ − ϕ − ω) T cosϕ .

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

Поэтому определяется минимальное значение величины E по углу θ. Этот минимум и будет искомой величиной силы пассивного давления грунта на подпорную стенку Ep.

Отсюда следует еще одно определение пассивного давления:

пассивное давление это минимальное давление стенки на грунт, при котором будет происходить выпор грунта в сторону от стенки.

2.4.8.Формула Понселе

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

(c = 0) с прямолинейной границей, наклоненной к горизонтали на угол α. Единственной массовой силой является удельный вес грунта γ, т.е. ω = 0. Пригрузка на поверхности засыпки отсутствует (P = 0). Расчетная схема к формуле Понселе показана на рис. 2.64.

Линией скольжения BC выделена предполагаемая область обрушения. Вес грунта в пределах этой области составит:

G =

1

γh2

cos(ε − α)sin(ε + θ)

.

(2.91)

2

 

 

 

cos2 εcos(α + θ)

 

При P = 0 и T = 0, с учетом (2.91), величина давления на стенку определится формулой

E =

1

γh2

cos(ϕ + θ)cos(ε − α)sin (ε + θ)

.

(2.92)

2

sin(ε + δ + ϕ + θ)cos2 εcos(α + θ)

 

 

 

 

 

 

 

 

 

193

O

x

 

 

 

C

z

 

 

 

 

 

 

A

α

 

G

π/2

 

 

 

 

 

 

π/2

 

 

 

ϕ

h

Eah

ε

θ

 

T

 

δ

 

R

 

 

 

 

 

E a

Eav

B

 

 

Рис. 2.64. Схема к формуле Понселе

Продифференцировав это выражение по dθ, а затем приравняв нулю, получим

cos(ε + ϕ + 2θ)cos(ε − α)sin (ε + δ + ϕ + θ)cos2 ε×

×cos(α + θ)cos(ε + ϕ + δ + α + 2θ)cos2 εcos(ϕ + θ)× (2.93)

×cos(ε − α)sin (ε + θ) = 0.

Решение этого уравнения дает угол θ = θ0, определяющий положение наиболее опасной линии скольжения, и при котором, очевидно, величина E достигает максимального значения. Полученная величина и является искомой силой активного давления Ea.

Данное решение используется в нормативных документах в следующем виде [24]. Сила активного давления раскладывается на горизонтальную и вертикальную составляющие. Горизонтальная составляющая Eah дается выражением

E

=

1

γh2

 

 

.

(2.94)

λ

 

ah

 

2

 

 

a

 

 

Коэффициент активного давления в рассматриваемом решении λa определяется формулой

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

cos(ϕ − ε)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λa

=

 

 

 

 

 

 

.

(2.95)

 

 

 

 

 

 

 

 

 

sin (ϕ + δ)sin(ϕ − α)

 

 

 

 

 

 

 

 

 

 

cosε 1

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos(ε + δ)cos(ε − α)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

194

Выражения (2.94) и (2.95) называются формулой Понселе. Угол θ0 находится в следующем порядке:

tg θ

 

=

cosα − ηcosϕ

,

η =

cos(ε − α

)

.

(2.96)

0

sin α − ηsin ϕ

 

 

 

 

 

 

 

 

 

 

 

cosε λa

 

 

 

 

 

 

 

Вертикальная составляющая силы активного давления дается формулой

Eav = Eah tg(ε + δ).

(2.97)

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

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

Глава 3. СОВРЕМЕННЫЕ НАПРАВЛЕНИЯ ТПРГ

§ 3.1. Плоская задача статики континуально неоднородной сыпучей среды

3.1.1. Актуальные проблемы ТПРГ

Идеально-жесткопластическая модель среды, принятая в теории предельного равновесия грунтов, давно завоевала себе прочное место как в инженерных расчетах, так и в научных исследованиях. И хотя эта модель накладывает ряд серьезных ограничений на поведение реальных пластических сред, существенно идеализирует их, она остается крайне востребованной, ибо, вопервых, основанные на решениях ТПРГ методы расчета получили широкую апробацию и подтвердили свою надежность, а вовторых, в теоретическом отношении решения ТПРГ зачастую становятся базовыми и для других моделей. На сегодняшнем этапе развития «Механики грунтов» совершенствование расчетных схем идет не только по пути создания более универсальных моделей (упругопластических, вязкоупругопластических), одновременно описывающих несколько механических свойств материала

195

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

ТПРГ развивается сегодня по нескольким направлениям. Вопервых, это получение новых решений в классической постановке [33, 34, 18], во-вторых, различные обобщения классической ТПРГ [1, 40, 17], в-третьих, развитие осесимметричной задачи [12] и многое другое. Также здесь можно выделить в отдельное, строго говоря, смежное направление создание принципиально новых по своей надежности методов расчета в рамках существующих приближенных схем, о чем шла речь выше [13].

Здесь будут рассматриваться решения ТПРГ для континуально неоднородной среды с условием прочности, обобщающим закон Кулона–Мора в части учета нелинейности графика сдвига [26, 3, 40, 11]. Термин континуально неоднородная среда означает, что параметры грунта непрерывно меняются от точки к точке, например, по некоторой экспериментально установленной непрерывной функции. Полученные при таком условии прочности уравнения ТПРГ позволяют рассмотреть целый ряд очень важных практических задач, в том числе перейти к задачам о несущей способности консолидирующихся оснований [27]. В завершение настоящей главы рассмотрим одну частную задачу, имеющую очень большое значение для современных геотехнических расчетов и проектирования оснований и фундаментов.

3.1.2.Прочность континуально неоднородной сыпучей среды

снелинейным графиком сдвига

Как уже говорилось в первой главе, прочность грунтов исследуют в сдвиговых испытаниях. В результате таких испытаний получают опытную зависимость предельных касательных напряжений τn от нормальных σn по площадке скольжения:

τn

 

= f (σn ).

(3.1)

 

График этой функции, напомним, называется графиком сдвига или графиком среза грунта. Ясно, что функция ƒ(σn) является огибающей предельных кругов Мора (рис. 3.1).

Наибольшее распространение получила простейшая аппроксимация экспериментальной зависимости (3.1) прямая Кулона

196

(1.12) и ее частные случаи (1.13) и (1.14). Решениям ТПРГ с линейным законом Кулона посвящены первые две главы настоящего пособия.

 

τ n

 

τ n = f

(σn )

 

B'''

 

 

 

 

 

 

 

 

 

 

 

B''

 

 

 

 

 

H

B'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

µ''

 

 

 

 

σ'

µ'

σ'

σ'''

σ''

µ'''

σ'''

 

 

3

1

3

1

 

1

A

O

σ''3

C'

C''

 

C'''

 

σn

Рис. 3.1. Нелинейный график сдвига τn = ƒ(σn)

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

В качестве примера приведем несколько таких функций, описывающих нелинейные графики сдвига [40]:

− параболическая

 

 

 

σ

n

+ H 1/n

 

 

τn

= k

λ

 

 

 

;

 

 

 

 

 

 

 

 

k

 

− экспоненциальная

τ

n

 

= k

1− exp

 

−λ σn

+ H

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

− условие М.М. Протодьяконова (младшего)

 

 

 

σ

n

b 3/8

 

τn

= 0,73a 1

+

 

 

.

 

 

 

 

 

 

 

 

a

В этих формулах: k, H, n, λ, a, b − параметры, определяемые при обработке опытных сдвиговых данных. Графики этих функ-

197

ций в плоскости τnOσn даны на рис. 3.2, а, для следующих значений параметров: k = 0,1; H = 0,03; n = 2; λ = 1; a = 0,05; b = 0,01.

Экспериментально полученную нелинейную зависимость (3.1) можно обобщить в виде условия прочности (1.7), представляющего собой функцию главных напряжений. Предположим, что далее это условие прочности можно записать в виде:

τ = τ(σ),

(3.2)

где τ и σ − полуразность и полусумма главных напряжений:

τ =

σ1 − σ3

,

σ =

σ1 + σ3

.

2

2

 

 

 

 

Так, если условие прочности Кулона–Мора в главных напряжениях (1.15) переписать в виде (3.2), то получим

τ = σsin ϕ + ccosϕ .

(3.3)

Условие прочности вида (3.2) оказалось весьма удобным как для теоретического анализа, так и для практических расчетов, и именно к этому виду условия прочности грунта мы будем в дальнейшем обращаться.

Приведем несколько условий прочности, записанных в виде (3.2) [40]:

степенное

 

 

λ

σ + H 1/n

;

 

τ = k

 

k

 

 

 

 

 

 

 

 

 

экспоненциальное

 

 

 

 

 

 

 

 

 

 

 

−λ

σ + H

;

τ = k 1

exp

k

 

 

 

 

 

 

 

 

условие В.В. Соколовского

 

 

 

 

 

 

σ + H

при

σ + H

π

.

τ = k sin

 

 

2

 

k

 

k

 

 

Графики указанных функций пластичности в плоскости τOσ даны на рис. 3.2, б для следующих значений параметров: k = 0,1; H = 0,03; n = 2; λ = 1.

198

Рис. 3.2. Нелинейные аппроксимации графика сдвига (а) и условия прочности (б): 1 параболическая функция; 2 экспоненциальная;

3 уравнение М.М. Протодьяконова; 4 степенное условие прочности; 5 экспоненциальное; 6 условие В.В. Соколовского

199

3.1.3. Постановка плоской задачи статики континуально неоднородной сыпучей среды

Рассмотрим плоскую задачу ТПРГ с условием прочности (3.2). Для континуально неоднородной среды параметры прочности являются функциями координат. Тогда плоская задача предельного равновесия континуально неоднородной сыпучей среды с нелинейным графиком сдвига примет вид:

∂σx +

∂τxz

= X ,

∂τxz

+

∂σz

= Z,

 

x

z

 

x

z

 

(3.4)

τ = τ(σ, x, z).

Преобразования системы (3.4) выполняются в целом аналогично тому, как это делалось в главе 1 для случая линейного закона Кулона, поэтому здесь будем следовать описанной ранее последовательности действий, по возможности избегая повторений.

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

σ

 

 

τxz = τ sin2α ,

 

 

x

= σ m τ cos2α,

(3.5)

σz

 

 

 

где α − угол между направлением σ1 и осью Oz (см. рис. 1.19, б). Заметим, что, приняв закон прочности Кулона–Мора (3.3),

уравнения (3.5) можно без труда привести к виду (1.33), полученному ранее для этого случая.

Подставим выражения для компонент напряжений (3.5) в уравнения равновесия. Для этого предварительно выполним дифференцирование функции τ = τ(σ, x, z), определяющей прочность неоднородной среды в каждой точке, по координатам:

∂τ = ∂τ ∂σ + ∂τ , x ∂σ ∂x x

∂τ = ∂τ ∂σ + ∂τ . z ∂σ ∂z z

С учетом полученных выражений возьмем частные производные от функций напряжений (3.5):

∂σx

=

(σ − τ cos 2α) =

∂σ

∂τ

cos 2α + 2τsin 2α

∂α

=

x

 

 

 

 

 

x

 

 

 

x

x

 

 

x

 

 

=

∂σ

∂τ ∂σ

+

∂τ

 

 

∂α

=

 

 

 

 

x

 

cos2α + 2τsin 2α

x

 

 

 

 

 

∂σ ∂x

 

x

 

 

 

 

 

200

=

∂σ

∂τ ∂σ cos2α −

∂τ cos2α + 2τsin 2α

∂α =

 

x

∂σ ∂x

 

 

x

 

x

 

=

∂σ

 

∂τ

 

∂τ

cos2α + 2τsin 2α

∂α

.

1

∂σ

cos2α

x

x

 

x

 

 

 

 

 

 

 

Выполнив сходные операции дифференцирования по x и z остальных уравнений (3.5), по аналогии с (1.34) сможем записать:

∂σx

 

= A

∂σ

+ A

∂α

+ B ;

∂τxz

= A

∂σ

+ A

 

∂α

+ B ;

x

 

 

 

 

 

 

11

x

12

x

11

z

13

z

14

 

z

12

∂τxz

 

∂σ + A

∂α + B ;

 

∂σz

 

∂σ + A

∂α

(3.6)

= A

 

= A

+ B .

 

z

z

x

21

x

22

x

21

 

23 z

 

24

22

В формулах (3.6) для краткости введены обозначения:

A11 =1− τ′σ cos 2α, A13 = τ′σ sin 2α, A21 = τ′σ sin 2α, A23 =1+ τ′σ cos 2α,

τ′σ ∂σ∂τ ,

A12

 

= 2τsin 2α,

B11 = −τ′x cos 2α;

A14

 

= 2τcos2α,

B12

= τ′z sin 2α;

A22

= 2τcos 2α,

B21 = τ′x sin 2α; (3.7)

A24

= −2τsin 2α,

B22

= τ′z cos 2α;

τ′

∂τ

,

τ′

∂τ

.

 

 

x

 

 

x

z

 

z

 

 

 

 

 

Подставим полученные после дифференцирования зависимости (3.6) в уравнения равновесия, а затем приведем подобные слагаемые:

A

∂σ + A

∂α + A

∂σ + A

∂α = X B B ,

11

x

12

x

13

z

14

z

11

12

 

∂σ + A

∂α + A

∂σ + A

∂α

 

(3.8)

A

= Z B B .

21

x

22

x

23

z

24

z

21

22

Присовокупим к уравнениям (3.8) выражения для частных дифференциалов искомых функций:

∂σ dx = dσ −

∂σ dz,

∂α dx = dα −

∂α dz .

x

z

x

z

Теперь умножим оба уравнения (3.8) на dx и подставим в них выражения для частных дифференциалов. Результат запишем в виде:

201

a ∂σ + a

∂α = b ,

a

∂σ + a

∂α = b .

11 z

12 z

1

 

21 z

22 z

2

В уравнениях (3.9):

a11

= −A11dz + A13dx,

a21

= −A21dz + A23dx ,

a12

= −A12dz + A14dx,

a22

= −A22dz + A24dx;

 

b1

= (X B11 B12 )dx A11dσ − A12dα ,

b2

= (Z B21 B22 )dx A21dσ − A22dα .

(3.9)

(3.10)

Решим систему (3.9) методом Крамера относительно производных искомых функций:

∂σ =

b1a22

b2a12

,

∂α =

a11b2

a21b1

.

(3.11)

 

 

 

 

z a a

22

a

21

a

z a a

22

a a

 

11

 

12

 

11

21

12

 

 

Как и в случае с линейным законом прочности, решение будем искать, приравнивая нулю числитель и знаменатель одновременно. Из равенства нулю знаменателя получим уравнения характеристик, а приравняв нулю любой из числителей (3.11), получим уравнения, связывающие искомые функции α и σ вдоль этих характеристик.

Итак, положив равным нулю знаменатель:

a11a22 a21a12 = 0,

а также выполнив необходимые алгебраические преобразования с учетом (3.7) и (3.10), получим:

(dx)2 (2τ)(τ′σ + cos 2α)+ 4τsin 2α dx dz +

+2τ(cos 2α − τ′σ )(dz)2 = 0,

или

dx 2

dx

 

 

 

 

(τ′σ + cos 2α)2sin 2α

 

(cos 2α − τ′σ ) = 0.

 

dz

dz

 

 

Корни данного квадратного уравнения определят уравнения характеристик:

dx

 

sin 2α ± 1− τ′

2

 

 

 

=

σ

 

.

(3.12)

dz

τ′σ + cos2α

 

 

 

 

 

202

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