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

668

.pdf
Скачиваний:
6
Добавлен:
06.12.2022
Размер:
3.91 Mб
Скачать

Эффективные уплотняющие напряжения в параметрах мгновенной прочности (3.26) для этого случая определим как сумму гидростатических бытовых и дополнительных, рассчитанных по теории линейно деформируемой среды от равномерной полосовой нагрузки (задача Мичелла [32, 37, 22]):

σ*ef = (σ*xef

+ σ*zef ) / 2 = γz + q +

 

+

p*

 

a x

+ arctg

a + x

(3.40)

 

arctg

 

 

 

,

 

 

 

 

π

 

z

 

z

 

 

где p* уплотняющая нагрузка, a = b/2 половина ширины штампа.

Решение достигается численным интегрированием канонических уравнений (3.15) для условия прочности (3.25) и параметров мгновенной прочности (3.26) с учетом (3.40) или численным интегрированием канонических уравнений (3.33) с учетом (3.31), (3.25), (3.26), (3.40) в рамках последовательности краевых задач, показанной на рис. 2.15.

На рис. 3.12 дана сетка линий скольжения в консолидирующемся основании одиночного штампа при следующих исходных данных: ширина штампа b = 1 м, угол внутреннего трения ϕ = 20°, удельное сцепление c = 10 кПа, удельный вес грунта γ = 20 кН/м3, пригрузка q = 3 кПа, коэффициент порового давления β = 0,5, уплотняющая нагрузка p* = 50 кПа. На этом же рисунке приведена искомая эпюра предельного нормального давления pu(x) и касательного τu(x) по подошве шероховатого штампа. Равнодействующая предельного давления на штамп шириной 1 м составила Pu = 137 кН.

Далее рассмотрим результаты серии расчетов, выполненных при следующих исходных данных: ширина штампа b = 1 м, угол внутреннего трения ϕ = 20°, удельное сцепление c = 10 кПа, удельный вес грунта γ = 20 кН/м3, пригрузка q = 0 кПа. В качестве основных факторов, влияющих на несущую способность консолидирующегося основания, примем коэффициент порового давления β и величину уплотняющего давления p*.

На рис. 3.13 показано влияние коэффициента порового давления β на величину среднего нормального предельного давления

233

pu штампа на консолидирующееся основание при различных уплотняющих нагрузках p* = 0, 50, 100, 150 кПа. На рис. 3.14 даны зависимости предельного давления pu от величины уплотняющей нагрузки p* при различных значениях β.

234

Рис. 3.12. Пример сетки линий скольжения в консолидирующемся основании

235

250

 

 

 

 

 

 

 

 

 

 

 

p u , кПа

 

 

 

 

 

 

 

 

 

200

 

 

 

 

 

 

 

 

 

 

150

 

 

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

 

 

50

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

β

 

0

 

 

 

 

 

 

 

 

 

 

0

0,1

0,2

0,3

0,4

0,5

0,6

0,7

0,8

0,9

1

 

p* = 0

p* = 50

 

p* = 100

p* = 150

 

Рис. 3.13. Зависимость предельного давления штампа на консолидирующе-

 

еся основание от коэффициента порового давления

 

250

 

 

 

 

 

 

 

 

 

p u , кПа

 

 

 

 

 

 

 

200

 

 

 

 

 

 

 

 

150

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

 

 

50

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p *, кПа

0

 

 

 

 

 

 

 

 

0

20

40

60

80

100

120

140

160

 

β = 0

β = 0,25

β = 0,5

 

β = 0,75

β = 1

 

Рис. 3.14. Зависимость предельного давления штампа

 

на консолидирующееся основание от величины уплотняющей нагрузки

236

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

3.3.4. Замечание о первой критической нагрузке на консолидирующееся основание

Решения ТМП, как было показано выше, осуществляются при уже известном поле эффективных напряжений в основании. Их определение на сегодняшний день удобнее всего выполнять на базе решений теории фильтрационной консолидации. Если обратиться к практическим методам этой теории, то как в варианте Терцаги–Герсеванова–Флорина, так и в варианте Флорина–Био [39, 36, 9] наибольшее применение находят расчеты, основанные на описании поведения линейно деформируемого скелета грунта. Использование линейно деформируемой модели связано также и с некоторыми проблемами, существующими в упругопластических расчетах. Об этом пойдет речь ниже.

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

§3.4. Об определении предельной нагрузки

вупругопластических расчетах грунтовых оснований методом конечных элементов

На сегодняшний день в инженерных геотехнических расчетах метод конечных элементов (МКЭ) является наиболее часто используемым численным методом [8, 22, 23, 32]. Настоящий факт обусловлен несколькими причинами. Во-первых, это его универсальность и та простота, с которой можно учитывать самые разнообразные свойства грунтов, вести вычисления по очень сложным моделям. Во-вторых, это интенсивное развитие и широкое внедрение в практику проектирования ряда вычислительных комплексов, основанных на МКЭ, таких как PLAXIS, FEM-model, MIDAS

237

и др. То обстоятельство, что расчеты МКЭ становятся повседневным явлением в практике проектирования, заставляет особенно внимательно отнестись к проблемам, возникающим при использовании МКЭ.

Одной из таких проблем, которую можно считать ключевой при проектировании оснований и фундаментов, является проблема оценки предельной нагрузки в упругопластических расчетах МКЭ. В этом параграфе рассмотрим процесс определения предельной нагрузки в упругопластических расчетах и сопоставим эти результаты с решением ТПРГ.

3.4.1. О постановке упругопластической задачи деформирования грунта

Основные механические свойства идеально упругопластической модели были рассмотрены в § 1.2 (см. рис. 1.3, б). Кратко опишем постановку и принципиальную схему решения упругопластических задач деформирования грунта.

В упругопластическом теле существуют области как упругих, так и пластических деформаций. Напомним, что грунт проявляет упругость в очень незначительной степени (см. § 1.2). Поэтому термин «упругий» здесь носит условный характер, указывая лишь на работу грунта в допредельном состоянии и линейную связь между напряжениями и деформациями.

Согласно общим положениям «Механики сплошной среды» [21] в «упругой» области должны выполняться:

уравнения равновесия (1.1) или (1.2); уравнения, связывающие напряжения и деформации, в

частности, закон Гука [22]; уравнения совместности [22].

В пластической области должны выполняться:

уравнения равновесия (1.1) или (1.2);

условие прочности;

уравнения, связывающие напряжения и приращения деформаций [21, 22];

уравнения совместности полных деформаций [21].

Кроме этого, на границе «упругой» и пластической областей должны выполняться условия непрерывности напряжений и де-

238

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

Большая часть из перечисленных уравнений либо знакома читателю из курса «Сопротивления материалов», либо обсуждалась в настоящем пособии, за исключением уравнений совместности и связи между напряжениями и приращениями деформаций. Требование совместности деформаций означает, что если условно разбить тело на элементарные прямоугольники до нагружения, то после приложения нагрузки, сдеформировавшись, эти прямоугольники вновь составят это же тело после деформации, причем будут плотно прилегать друг к другу [7]. О зависимостях между напряжениями и приращениями деформаций речь пойдет чуть ниже.

3.4.2. Основные положения упругопластического расчета МКЭ

Рассмотрим порядок упругопластического расчета в рамках МКЭ. Описание МКЭ можно найти, например, в [8]. Для понимания излагаемого ниже материала не понадобится специальных знаний в области конечно-элементного анализа, а по мере необходимости будут даны соответствующие пояснения.

Нагрузка на основание увеличивается пошаговым методом. После того как сделан очередной шаг нагружения, в каждой точке грунта (в МКЭ это означает для каждого конечного элемента) вычисляется функция (1.18):

f (σij ) = (σx − σz )2 + 4τ2xz (σx + σz + 2c ctgϕ)sin ϕ.

Очевидно, что если ƒ(σij) < 0, то грунт находится в допредельном состоянии (см. рис. 1.13, а) и его деформирование описывается законом Гука. Если ƒ(σij) = 0, то в данной точке грунта достигнуто предельное напряженное состояние (см. рис. 1.13, б) и процесс деформирования будет описываться специальным уравнением состояния. Случай ƒ(σij) > 0 проиллюстрирован рис. 1.13, в такое напряженное состояние не может существовать.

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

239

при ƒ(σij) < 0

 

{dσ} = [D ]{dεe};

(3.41)

e

 

при ƒ(σij) = 0

 

{dσ} = [Dep ]{dε}.

(3.42)

Уравнение (3.41) определяет процесс линейного деформирования грунта в «упругой» области и представляет собой обобщенный закон Гука, записанный в матричной форме. Матричное уравнение (3.42) определяет пластическое деформирование грунта

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

Вформулах (3.41) – (3.42):

{dσ} и {dε} – матрицы-столбцы приращений напряжений и полных деформаций, которые определяют изменение напряжен- но-деформированного состояния грунта в процессе нагружения:

{dσ}T = {dσ

, dσ

, dτ }

и {dε}T = {dε

, dε

, dγ

xz

};

x

z

xz

x

z

 

 

{dεe} матрица-столбец «упругих» деформаций:

 

 

 

 

 

{dεe}T = {dεxe, dεze, dγxze};

 

 

 

 

 

1− ν

 

ν

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12ν

12ν

 

 

 

 

E

 

ν

 

1− ν

 

 

 

[D ] =

 

 

 

 

0

 

«упругая» матрица для

1

 

 

12ν

e

+ ν 12ν

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

 

 

 

 

 

 

0

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

условий плоской деформации;

 

D11

D12

D13

 

 

[D ] =

D

D

D

 

упругопластическая матрица.

ep

21

22

23

 

 

 

D

D

D

 

 

 

31

32

33

 

 

Уравнения (3.41) и (3.42) выражают только часть требований, предъявляемых к упругопластическим решениям. Что касается остальных положений упругопластического расчета, то можно считать, что уравнения МКЭ автоматически удовлетворяют условиям совместности и условиям равновесия каждого узла конечноэлементной сетки [8].

240

3.4.3. Уравнение состояния упругопластической среды

Уравнение состояния (3.42) определяется следующими положениями.

Первое. Напряжения, достигнув предельных значений ƒ(σij) = 0, могут на следующих этапах нагружения менять свои значения σij + dσij, но так, чтобы оставаться на поверхности текучести ƒ(σij + dσij) = 0. Сказанное можно записать в виде:

df (σij ) = ∂σfx

или в матричной форме:

где

{F}Ò

dσx

+

f

dσz

+

f

dτxz = 0,

 

 

 

 

 

 

∂σz

 

∂τxz

{F}Ò{dσ} = 0,

(3.43)

 

f

 

 

f

f

 

=

 

 

 

 

 

 

 

 

.

∂σx ∂σz

 

 

 

∂τxz

В качестве условия прочности ƒ(σij) здесь и далее принимаем закон Кулона–Мора (1.18).

Второе. Полагают, что полные приращения относительных деформаций {dε} есть сумма приращений «упругих» {dεe} и пластических деформаций {dεp}:

dε

x

= dεe

+ dεp ,

dε

z

= dεe

+ dεp ,

dγ

xz

= dγe

+ dγ p

,

 

x

x

 

z

z

 

xz

xz

 

или в матричной форме:

 

 

 

 

 

 

 

 

 

 

 

 

{dε} ={dεe}+{dεp},

 

 

 

 

(3.44)

где {dεp} матрица-столбец приращений пластических деформа-

ций:

{dεp}T={dεxp, dεzp, dγxzp}.

Третье. «Упругая» часть приращения деформаций (3.44) связана с приращением напряжений законом Гука:

{dσ} = [D ]{dεe}.

(3.45)

e

 

«Упругая» матрица [De] дана выше.

Следует обратить внимание, что выражение (3.45) определяет «упругую» составляющую приращений полных деформаций в

241

пластической области, а выражение (3.41) – деформирование в «упругой» области.

Четвертое. Пластическая часть приращений полных деформаций определяется законом пластического течения, о котором уже шла речь в § 1.3. Закон пластического течения устанавливает связь между полными напряжениями и приращениями деформаций:

dεp = dλ

P

,

dεp = dλ

P

,

dγ p = dλ

P

,

 

∂σz

 

x

∂σx

z

 

 

 

 

 

xz

∂τxz

 

 

 

 

 

 

 

 

 

 

 

или в матричном виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{dεp} = dλ{P},

 

 

(3.46)

где dλ − некоторый множитель;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

 

P

 

P

 

 

 

{P}Ò =

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂σx

 

∂σz

 

∂τxz

 

 

Функция от напряжений P(σij), напомним, называется пластическим потенциалом. Если она совпадает с функцией пластичности P(σij) ≡ ƒ(σij), то выражения (3.46) представляют собой

ассоциированный закон пластического течения, в противном случае (3.46) представляют неассоциированный закон пластиче-

ского течения.

Исключая из системы уравнений (3.43) – (3.46) приращения компонент «упругих» {dεe} и пластических {dεp} деформаций, а также множитель dλ, определим искомую зависимость (3.42). Указанные преобразования запишем в матричной форме. Из выражения (3.44) получим

{dεe} ={dε}{dεp}

и подставим в (3.45)

{dσ} = [De ]{dε}[De ]{dεp}. Учитывая (3.46), имеем

{dσ} = [De ]{dε}[De ]{P}dλ,

(3.47)

а после подстановки в (3.43):

 

 

{F}Ò[D ]{dε}dλ{F}Ò[D ]{P} = 0.

 

e

e

 

242

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