668
.pdfВыбор одного из двух возможных вариантов граничных условий − σ1 = ƒn, α = 0 или σ3 = ƒn, α = π/2 − осуществляется из дополнительных соображений, изложение которых будет продолжено в следующей главе.
Глава 2.
СТАТИЧЕСКИЕ РЕШЕНИЯ ОСНОВНЫХ ЗАДАЧ ТПРГ
§2.1. Несущая способность основания штампа
Винженерных расчетах предельной стадии работы грунтовых массивов исторически сформировались три больших класса задач: задачи о несущей способности оснований фундаментов, задачи устойчивости откосов и склонов, задачи о предельном давлении грунта на подпорные сооружения. Здесь рассматриваются все перечисленные проблемы с позиций строгого статического метода ТПРГ, и, если это необходимо, даются основные положения приближенных методов расчета, наиболее востребованных в практике проектирования.
Построение решений для конкретных практических схем
начнем с рассмотрения одной из центральных задач в ТПРГ − задачи о внедрении штампа в жесткопластическую среду − задачи Прандтля. В замкнутом виде на сегодняшний день это решение дано для двух частных случаев: невесомой сыпучей среды и иде- ально-связного весомого основания. Сделав несколько предварительных замечаний, дадим оба этих решения, а затем перейдем к построению численного решения для общего случая весомого сыпучего основания.
2.1.1. Минимальное и максимальное напряженное состояние
Рассмотрим предельное равновесие грунтового основания, ограниченного осью Ox, на поверхности которого действует равномерно распределенная нагрузка интенсивностью fn (рис. 2.1). Грунт характеризуется прочностными параметрами ϕ и c, а также удельным весом γ. Запишем граничные условия для этой схемы:
σz = fn , τxz = 0 при z = 0. |
(2.1) |
|
83 |
fn
O
Поскольку по всей границе напряжения не зависят от координаты x, то можно ожидать, что
xони не будут зависеть от нее и внутри основания. Тогда уравнения равновесия (1.2) (или первые два из системы (1.27)) примут вид:
dτxz |
= 0, |
dσz |
= γ . |
|
|||
|
|
||
dz |
dz |
Частные производные здесь заменены на обычные вследствие независимости напряженного состояния от x. Интегрируя полученные выражения с учетом граничных условий (2.1), получим:
|
σz = γz + ƒn, |
τxz = 0. |
|
||
Подставим полученные выражения в формулы (1.33): |
|
||||
γz + fn = σ(1+ sin ϕcos2α) − c ctgϕ, |
0 = σsin ϕsin 2α. |
|
|||
Полученная система имеет два решения [26]: |
|
||||
σ = |
γz + fn + c ctgϕ |
|
α = (1− χ) |
π |
|
|
, |
4 , χ = ±1. |
(2.2) |
||
1+ χsin ϕ |
Поскольку угол α между осью Oz и направлением σ1 всюду постоянен и равен либо 0, либо π/2, то линии скольжения будут представлены двумя семействами прямых, составляющих с вертикалью углы ±µ при χ = +1 или π/2 ± µ при χ = −1.
Соответственно два решения получим и для напряжений:
σx |
= (γz + fn |
+ c ctgϕ) |
1 |
− χsin ϕ |
− c ctgϕ, |
χ = ±1, |
|
+ χsin ϕ |
|||||
|
|
1 |
|
|
||
|
|
|
σz = γz + fn , |
(2.3) |
τxz = 0.
Непосредственно из (2.2) следует, что решение χ = +1 дает меньшее значение среднего напряжения, а решение χ = −1 дает большее. Поэтому напряженное состояние, отвечающее случаю
χ = +1 называется минимальным, а χ = −1 – максимальным.
84
Физический смысл этих двух |
а) |
М |
б) |
М |
решений обычно иллюстрирует- |
||||
ся моделью с гирями и стерж- |
|
A |
|
A |
нем, показанной на рис. 2.2, а, б |
|
B |
|
B |
[26]. Чашка A закреплена на |
|
|
||
стержне, скользящем с большим |
|
|
|
|
трением в направляющих B. Ес- |
|
Fmin |
|
Fmax |
ли масса груза M на чашке мала, |
|
|
||
|
|
|
|
|
то, учитывая наличие трения, си- |
Рис. 2.2. Модель, иллюстрирующая |
|||
стема останется в состоянии по- |
|
понятия минимального (а) |
||
коя. Если же масса груза доста- |
|
|||
|
и максимального (б) |
|
||
точна для преодоления сил тре- |
|
напряженного состояния |
||
ния, то удержать систему в рав- |
|
|
|
|
новесии может только дополнительно приложенная сила F. |
|
|||
Существуют два предельных значения этой силы, при кото- |
||||
рых еще сохраняется равновесие системы: если F станет меньше |
||||
некоторого минимального значения Fmin, то чашка с грузом |
||||
начнет движение вниз, если F станет больше некоторого макси- |
||||
мального значения Fmax, то чашка с грузом начнет движение |
||||
вверх. Положения чашки с грузом на пути ее перемещений пока- |
||||
заны на рис. 2.2 пунктиром. |
|
|
|
|
По аналогии с моделью можно сделать важный вывод о том, |
||||
что минимальное напряженное состояние грунтового массива |
||||
приводит к оседанию поверхности, максимальное − к поднятию |
||||
поверхности или выпиранию грунта. |
|
|
|
|
Итак, минимальное напряженное состояние дается решения- |
||||
ми (2.2) и (2.3) при χ = +1 (рис. 2.3, а): |
|
|
|
|
параметры канонической системы уравнений вдоль линий |
||||
скольжения − |
|
|
|
|
σ = |
γz + p + c ctgϕ |
, α = 0 , x = ±z tgµ + C , |
(2.4) |
|
1+ sin ϕ |
||||
|
|
|
предельные напряжения −
σx = (γz + p + c ctg ϕ)11+− sinsin ϕϕ − c ctgϕ,
σz = γz + p , |
(2.5) |
τxz = 0.
85
Максимальное напряженное состояние дается решениями (2.2) и (2.3) при χ = −1 (рис. 2.3, б):
параметры канонической системы уравнений вдоль линий скольжения −
|
γz + q + c ctg ϕ |
|
π |
|
π |
|
|
|||
σ = |
|
|
, α = |
|
, |
x = z tg |
|
± µ |
+ C , (2.6) |
|
1 |
− sin ϕ |
2 |
2 |
|||||||
|
|
|
|
|
|
предельные напряжения −
σx |
= (γz + q + c ctgϕ) |
1 |
+ sin ϕ |
− c ctgϕ, |
|
− sin ϕ |
|||
|
1 |
|
||
|
σz = γz + q, |
(2.7) |
τxz = 0.
Вформулах (2.4) и (2.6): C − постоянная интегрирования.
а) |
|
p |
|
б) |
q |
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
x |
O |
|
A |
|
O µ |
µ A |
|
|
|
π |
− µ |
2µ |
σ1 |
|
σ1 |
|
2 |
|
|||
|
|
|
|
|
C |
π |
2µ |
C |
α = 0 |
z |
α = |
||
z |
|
|
|
2 |
Рис. 2.3. Сетка линий скольжения в случае минимального (а) и максимального напряженного (б) состояния
Внешнее давление на основание в случае минимального напряженного состояния, когда поверхность основания движется вниз, будем называть нагрузкой, а в случае максимального, когда поверхность основания движется вверх, − пригрузкой.
Отметим важную особенность этих решений, характерную для ТПРГ. Если равномерная нагрузка (пригрузка) задана на ограниченном участке OA горизонтальной поверхности, то решения вида (2.4) или (2.6) будут определены в основании в пределах треугольной области OAC (см. рис. 2.3), образованной самим от-
86
резком OA, на котором действует внешняя нагрузка, и линиями скольжения, выходящими из крайних точек отрезка O и A. В случае минимального напряженного состояния эти линии составляют с горизонталью угол π/2 − µ (см. рис. 2.3, а), а в случае максимального напряженного состояния µ (см. рис. 2.3, б).
Такие области можно рассматривать как результат решения I краевой задачи (см. § 1.4), построенной при известных параметрах x, z, σ и α на OA.
2.1.2. Задача Прандтля. Предельное давление на невесомое основание
Расчетная схема к задаче Прандтля показана на рис. 2.4. На участке AA′ шириной b действует искомая предельная нагрузка pu, на остальной части поверхности − равномерная пригрузка q.
|
|
b |
|
|
|
|
pu |
|
q |
|
q |
D' |
A' O |
A |
D x |
z
Рис. 2.4. Граничные условия к задаче Прандтля
С точки зрения моделирования работы реальных оснований пригрузка представляет собой вертикальное давление от веса грунта, расположенного выше подошвы фундамента, поскольку каждый фундамент обычно имеет некоторую глубину заложения. В расчетах за горизонтальную границу основания принимается, как правило, уровень подошвы фундамента − его нижней плоскости, − а вышерасположенный грунт моделируется как пригрузка, равная q = γ d, где d − глубина заложения фундамента.
Сделаем предположение относительно характера искомой нагрузки pu, которую примем равномерной. В результате решения будет показано, что такое допущение оправдано.
87
Из самой постановки задачи понятно, что поверхность основания в пределах отрезка A′A будет перемещаться вниз, а границы AD и A′D′ будут испытывать поднятие. Таким образом, заключаем, что в область A′AC, примыкающая к искомой предельной нагрузке на отрезке A′A, будет областью минимального напряженного состояния, а области ABD и A′B′D′ − областями максимального напряженного состояния (рис. 2.5).
b
|
|
|
|
pu |
|
|
|
q |
|
q |
|
|
|
µ |
µ |
|
µ |
D' |
µ |
A' |
O A |
|
D x |
µ |
µ σ |
||||
|
σ1 |
µ |
2µ |
1 |
|
|
|
|
|
||
|
α = −π/2 |
C |
α = +π/2 |
||
|
|
|
|||
|
|
B' |
α = 0 |
B |
|
z σ1
Рис. 2.5. Области минимального и максимального состояния в задаче Прандтля
Положив в равенствах (2.4) и (2.6) γ = 0 и p = pu, получим:
в зоне A′AC − |
|
|
|
|
|
|
|||
σ = |
pu + c ctgϕ |
, α = 0 , |
x = ±z tgµ + C ; |
(2.8) |
|||||
1+ sin ϕ |
|||||||||
|
|
|
|
|
|
|
|||
в зонах ABD и A′B′D′ − |
|
|
|
|
|
|
|||
σ = |
q + c ctg ϕ |
, |
α = ± |
π |
, |
x = z tg(α ± µ)+ C . |
(2.9) |
||
1− sin ϕ |
2 |
||||||||
|
|
|
|
|
|
В выражениях (2.9) равенство α = π/2 является решением для области ABD, а α = −π/2 − решением для области A′B′D′ (см. рис. 2.5).
Заметим, что если очертание области A′AC вполне определено, то размеры зон ABD и A′B′D′ зависят от длин отрезков AD и A′D′, которые окончательно определятся только в процессе решения.
88
На этом этапе имеем три области предельного напряженного состояния, симметричные относительно оси Oz (см. рис. 2.5). Необходимо достроить поле предельных напряжений в частях ABC и A′B′C основания и объединить с ними воедино уже построенные области. Замыкать решение в этом смысле будут области радиального веера. Вначале определим их геометрию − получим уравнения линий скольжения, а затем построим непрерывное поле предельных напряжений.
Схема, показанная на рис. 2.5, позволяет предположить, что в переходных областях ABC и A′B′C одно семейство линий скольжения будет представлено пучком прямых, выходящих из точек A и A′. Характеристики другого семейства, соответственно, будут составлять с ними угол 2µ.
Поскольку расчетная схема симметрична (см. рис. 2.5), достаточно построить решение только для правой ее части. Рассмотрим область переходного веера ABC в невесомом основании (рис. 2.6, а). Согласно сказанному ранее из особой точки A выходит пучок прямых − характеристики 2-го семейства. Характеристики 1-го семейства образуют с ними угол 2µ и, следовательно, представляют собой логарифмические спирали с центром в особой точке A, поскольку именно логарифмическая спираль пересекает свои радиусы-векторы под постоянным углом. Ее частным случаем является окружность, когда этот угол равен π/2. Уравнения этих кривых можно получить следующим образом.
На рис. 2.6, б изображены элементарные отрезки линий скольжения 1-го и 2-го семейств. Для удобства будем использовать локальную полярную систему координат z′Ax. В этой системе положение точки задается радиус-вектором r и углом θ (см. рис. 2.6, б). Элементарная дуга спирали ограничена двумя точками − (r, θ) и (r + dr, θ + dθ). Учитывая, что в каждой точке линии скольжения 2-го семейства направлены по радиус-вектору, элементарная дуга составит с ним угол 2µ, и, соответственно, с отрезком rdθ, перпендикулярным радиусу, угол π/2 − 2µ = ϕ. Тогда дифференциальное уравнение характеристики запишется в виде:
dr = tgϕ. rdθ
89
а) |
b/2 |
|
|
|
|
|
б) |
b/2 |
|
π |
|
|
|
|
|
|
|
O |
2 − µ |
A |
|
|
|
O |
A |
|
|
|
|
|
µ |
ЛС-2 |
x |
|
|
|
r0 |
|
|
|
|
|
||
|
|
|
|
|
|
|
||
|
|
|
|
|
ЛС-1 |
|
|
|
|
|
|
|
|
|
|
|
|
C |
θ |
− |
|
θ+ |
|
σ1 |
|
θ |
|
|
|
|
|
||||
|
|
|
|
µ |
|
|
||
|
|
|
|
|
|
|
|
|
z |
σ1 |
|
|
µ |
σ1 |
B |
z |
z' |
z' |
|
|
r |
x |
dθ
|
|
d |
|
|
r |
r |
θ |
|
d |
ϕ |
|
r |
||
|
ЛС-2 2µ ЛС-1
Рис. 2.6. Область радиального веера (ЛС-1 и ЛС-2 − линии скольжения 1-го и 2-го семейства):
а − общий вид, б − схема к выводу уравнений характеристик
Проинтегрировав его, получим:
ln r − ln r = (θ − θ |
) tg ϕ, |
ln |
|
r |
= (θ − θ |
) tgϕ, |
|||||||||||
|
|
||||||||||||||||
0 |
|
|
0 |
|
|
|
|
|
|
|
|
|
0 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r0 |
|
|
где r0 − радиус кривой при значении угла θ = θ0. |
|
||||||||||||||||
Отсюда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r = r e(θ−θ0 )tgϕ . |
(2.10) |
||||||||||||||
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
Постоянные r0 и θ0 определим из рис. 2.6, а: |
|
||||||||||||||||
r |
= |
b |
|
|
1 |
|
|
, |
|
θ |
|
= −µ. |
|
||||
2 sinµ |
0 |
|
|||||||||||||||
0 |
|
|
|
|
|
|
|
|
|||||||||
Запишем окончательный вид уравнения характеристик 1-го |
|||||||||||||||||
семейства в переходной зоне ABC: |
|
|
|
|
|
|
|||||||||||
|
r = |
|
1 |
|
|
b |
|
|
e(θ+µ)tgϕ . |
(2.11) |
|||||||
|
|
2 sin µ |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
Теперь можно найти длину отрезка AB − ему будет согласно схеме на рис. 2.6 соответствовать угол θ = π/2 − µ:
AB = |
1 b |
|
eπ2 tgϕ . |
|||
|
|
|
|
|||
2 sin |
µ |
|||||
|
|
Зная AB, нетрудно рассчитать размеры всей области ABD.
90
Полученная область радиального веера может быть рассмотрена как результат решения II краевой задачи с особой точкой (см. § 1.4).
Таким образом, сетка линий скольжения в задаче Прандтля полностью определена (рис. 2.7).
|
lu |
|
|
b |
|
|
lu |
|
|
|
|
|
|
pu |
|
|
|
|
q |
|
|
|
|
|
q |
|
D' |
|
|
A' |
O |
A |
|
|
D |
x |
hu
B' |
C |
2 |
B |
|
z |
|
|
Рис. 2.7. Сетка линий скольжения в задаче Прандтля (невесомое основание)
Прежде чем переходить к определению предельного давления, отметим еще две важные величины, характеризующие габаритные размеры сетки линий скольжения в целом − это ширина призмы выпирания lu и глубина развития зон пластических деформаций hu (см. рис. 2.7). Ширина призмы выпирания очевидно равна:
π tgϕ |
. |
l = AD = 2ABcosµ = b ctgµ e2 |
|
u |
|
Глубина развития зон достигает максимума в точке, где касательная к логарифмической спирали параллельна оси Ox или, что то же самое, параллельна оси Ox будет характеристика 1-го семейства, т.е. в данной точке выполняется условие (см. рис. 2.6 или рис. 2.7):
θ + 2µ = π2 .
Соответственно радиус-вектор с учетом (2.11) в этой точке равен:
r = 12 sinbµ e(π/2−2µ+µ)tgϕ ,
91
а глубина развития пластических зон |
|
|
|
||||||
h = r cosθ = |
1 |
|
b |
|
e(π/2−µ)tgϕ cos |
π |
− 2µ |
|
= bcosµ e(π/2−µ)tgϕ . |
|
|
|
|
|
|||||
u |
2 sin |
µ |
|
|
|
||||
|
|
2 |
|
|
|
Вернемся к поиску величины предельного давления. Учитывая симметрию расчетной схемы, решение будем строить только с границы AD. Итак, в области ABD имеется полное решение (2.9) при α = +π/2, а в области A′AC − решение (2.8) с неизвестной величиной pu. Свяжем эти зоны. Для этого проинтегрируем в области радиального веера ABC второе из уравнений канонической системы (1.51), выбирая в нем верхний знак, отвечающий линиям скольжения 1-го семейства:
dσ + 2σ tgϕ dα = 0.
Разделяя переменные и интегрируя, получим (см. п. 1.4.5, (1.62)):
σ = σ0e−2tgϕ(α−α0 ) . |
(2.12) |
Здесь α0 и σ0 − постоянные, определяемые граничными условиями.
Определим постоянные α0 и σ0. Как было уже сказано, интегрирование ведем с границы AB до AC. Условия по σ и α на границе AB (область ABD):
σ |
|
= |
q + c ctgϕ |
, |
α |
|
= |
π . |
|
0 |
1− sin ϕ |
0 |
|||||||
|
|
|
|
|
2 |
Условия по σ и α на границе AC (область ACA′):
σ = |
pu + c ctg ϕ |
, α = 0. |
|
1+ sin ϕ |
|||
|
|
Подставляя граничные условия в (2.12), получим
pu + c ctg ϕ |
= q + c ctg ϕ e |
|
0− |
π |
||
|
2 . |
|||||
|
|
|
|
−2tgϕ |
|
|
1+ sin ϕ |
|
1− sin ϕ |
|
|
|
Выразим отсюда искомую предельную нагрузку:
p = (q + c ctgϕ) |
1 |
+ sin ϕ |
eπ tgϕ − c ctg ϕ. |
(2.13) |
|
|
|
||||
u |
1 |
− sin ϕ |
|
||
|
|
92