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

книги из ГПНТБ / Будин, А. Я. Тонкие подпорные стенки

.pdf
Скачиваний:
60
Добавлен:
20.10.2023
Размер:
8.33 Mб
Скачать

Решение для Uy можно искать в виде Uy= ^ n {x)e~ %пУ. После выполнения операций, аналогичных проведенным для

Ux, получено

 

 

е пУ [D1cosXnx Jr D2s'mXnx+(qx/Hi\Xn)cosXх].

(125)

п—0

 

 

На основании условия (д)

следует, что £>i = 0. Тогда

 

ОО

 

 

пУ[As sin К х + (qx/HnK)cosК х\-

(126)

п= 0

 

 

Подстановка соотношений (123) и (126) в уравнение (113)

дает

 

 

2 е пУ [(qlHr\Xn) sin Хпх + (С2 + qx/Hr\Xn) Хпcos Хпх

 

п—0

 

 

D2Xnsin Хпх (qXnx!Н-цХп) cos Хпх\ = 0,

 

откуда

 

 

С2 = 0;

D2 = q/Hr\X2n.

(127)

При учете зависимостей (127) можно найти окончательно:

оо—%пу

Ux = {qx!H'x\) 2 (1А„) е

sinA,„x;

(128)

п—0

 

 

U у — (q/Нц) ^ е %пУ [(*/*.„) cos?i„x+

(l/A,^) sinXnx].

(129)

Для решения задачи о напряженно-деформированном со­ стоянии больверков на ползучих основаниях необходимо знать величины Ux и Uy на напорной поверхности полуполосы.

Полагая в выражении (128) у = 0, можно получить

ОО

Ux | . = (qx/Ht])^ [2Ях/(2п+ 1) л] [sin(rc/2#) (2п + 1)].

Уп = 0

 

ОО

 

 

 

Но так как

2

[1/(2га+ 1)1 sin[(nx/2H)(2ti+ 1)] = п/4,

то

 

п—0

 

 

 

 

 

A ly^o= x q 12ч-

(130)

Аналогично для Uy:

ОО

 

 

 

 

 

и у | =

(q/Нц) {(2Нх/л)

2 [ 1/(2п + 1)] cos [(тсх/2Н) х

 

 

 

!

71=0

 

 

 

ОО

 

X (2п+ 1)] + (4#2/я2) 2

[l/(2n+ I)2] $in[(nx/2H) X

 

 

 

п = 0

 

X (2п + 1)]} = (4<7# / я2т]) {(ях/4Я) lnctg (лх/4Н) -f-

 

 

ОО

 

 

 

+

2

[1/(2я+ I)2] sin [(ях/2Н) (2n + 1)]} .

(131)

 

п= 0

 

 

 

121

Рис.

66.

К

анализу

закономерности релаксации

 

 

 

 

р ( х , 0)

 

 

 

 

а — взаимодействие стенки с ползучим

грунтом;

б — переход

 

 

 

к новой переменной z

 

 

 

Ряд в правой части выражения (131) быстро

сходится. По-

скольку и х =

~

j

и у =

у

(где

ех(х)

и

еу (х) — проек­

ции вектора перемещения точек напорной плоскости соот­ ветственно на оси х я у, a t — время), из уравнений (130) и (131) следует, что

 

ех (х, t) = qxtl2r\\

(132)

еу (х,

t) = (AqHt/n2r\) {(пх/АН) In ctg (пх/АН) +

 

 

ОО

 

+

2 [l/(2n+ I)2] sin [(лх/2Н) (2n+ 1)]} .

/ (133)

 

п

 

Наряду с формулой (133) далее удобно пользоваться упро­ щенным выражением для перемещения точек напорной поверх­ ности полуполосы

х

Н

t

(134)

еу (х, t) = (l/Нц) Цdx j

dx\q{x, t)dt,

o

*

o

 

дающим близкое совпадение с точным решением.

Для вывода зависимостей, характеризующих общие законо­

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

лосой в периоде установившейся ползучести грунта

(рис. 66, а).

В качестве отправных положений приняты заданными: эпюра

реактивного давления грунта р(х, 0), известная

из статиче­

ского расчета для начального момента времени, и внешняя, не­ изменная во времени, активная нагрузка N(x). Исходным усло­

122

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

 

еу (х, t)=y(x, t)— y(x, 0).

(135)

Здесь еу(х, t)

— ползучий прогиб напорной поверхности грунто­

вой полуполосы; у(х,

t ) — прогиб стенки в точке с координатой

х в момент времени /;

у(х, 0) — то же, но в начальный момент

времени (непосредственно после приложения нагрузки).

 

Значения у(х, t) и у(х, 0) определяются из уравнения упру­

гой линии

 

 

 

 

у(х,

t) — (l/EI)^dx^M*(x,

t)dx + C1x + Dx];

(136)

у{х,

0) = (1/EI)[$dx$M*(x,

0)dx + C2x + D2],

(137)

где М * (х, t) — изгибающий момент, действующий в некотором сечении стенки в момент времени t; М*(х, 0) — то же, в началь­ ный момент времени.

Для защемленной консольной стенки при принятом начале координат (см. рис. 66) постоянные интегрирования Си С2, Z?i и D2 равны нулю. Подстановка в соотношение (135) выражений (134), (136) и (137) приводит к интегральному уравнению

х

И

t

 

 

t) M*(x,0)]dx.

(138)

(EI/Hr\) §dx jj dx Цр (х, t)dt = jjdx § [М* (х,

О

я

о

 

 

 

 

 

 

 

Суммарные изгибающие моменты

 

 

 

 

 

 

 

М*(х,

t ) = M a(x, t)— .М(х,

t)\

j

 

 

 

 

 

М*(х,

0 )= M fl(*, 0)—М(*,

0),

j

 

1

где Ma(x,t)

и Ma{x, 0) — составляющие

изгибающих

моментов

от активной

нагрузки ЛДя),

а М(х, t) и М(х,

0 ) — составляю­

щие изгибающих моментов

от реактивной нагрузки

р(х, t) и

р(х, 0). Но,

поскольку при ЛД*) =const

Ма(х,

t)= Ma(x,

0),

из

уравнения

(138) следует:

 

 

 

 

 

 

*н t

(EI/Нц) JdxJ dx§p(x, t) dt = j dx j [M (x, 0)—M (x, t)\dx. (140)

оx о

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

Дважды дифференцируя равенство (140) по х, можно найти

t

(141)

— \р(х, t)dt = (Ht\/EI) [М (х, 0)—M ( x , t ) ].

о

 

123

Выражения для М(х, i) и М(х, 0) удобно записать, введя

новую

переменную

2 , связанную

с х соотношением

z = x—х*

(рис. 66,б):

 

я

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

М{х,

t) — b f

(z— x)p(z,

t)dz\

(142)

 

 

 

X*

 

 

 

 

 

 

 

 

 

Я

(z—x)p(z,

0)dz,

(143)

 

M (x, 0) —б J"

 

 

 

X*

 

 

 

 

 

 

где b — ширина стенки.*

 

 

(142)

и (143) в

формулу

После подстановки

выражений

(141)

получено

И

 

 

 

 

 

 

 

t

 

 

 

И

(*- -х) р (z, t) dz

— [ р (х, t) dt — a2

J

(z— x)p(z,

0)dz-

-I

 

—J р (х, t) dt = a2 j (z—x) \p (z,

t) — p(z, 0)] dz,

(144)

где

0

 

x*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2 = Hr\/EI.

Дифференцируя далее (144) по x, необходимо иметь в виду, что поскольку в интеграле J (z—x)[p(z, t)p(z, 0)]6z, перемен­ ная x является параметром, указанная операция представляет собой дифференцирование по параметру.

В результате дифференцирования

-^-\р(х, t)dt= — а2 Г [p(z, t) — p(z, 0)]dz.

дх о

i

 

Полученное соотношение, будучи продифференцированным

по х, принимает вид

 

 

Я2 t

 

(145)

-jp(x, t ) d t = — а2[р{х, t) р(х, 0)].

Дифференцирование (145) по t дает окончательно

 

д2р (х , t)

dp (х, t)

(146)

дх2

dt

 

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

параболического

типа

используется

в

теории

консолидации

К. Терцаги [76].

 

 

 

 

 

 

Таким образом, отыскание закона изменения во времени ре­

активного давления ползучего грунта

на

упругую

стенку сво­

* В связи с тем,

что

здесь рассматривается

работа

1

пог. м стенки,

в последующих соотношениях величина Ь опущена.

124

дится к решению уравнения (146) при соответствующих гранич­ ных и начальном условиях. Наличие зависимости для р(х, t) позволяет далее найти распределение в стенке изгибающих мо­ ментов М(х, t), а при необходимости и всех их производных; (перерезывающие силы, прогибы).

В случае, если опорное защемление стенки не является жест­ ким, а обладает податливостью, постоянные интегрирования в (136) и (137), определяемые, как обычно, опорными усло­ виями стенки, не равны нулю. Вводя коэффициент податливо­ сти защемления К, можно записать следующие выражения для смещения у и угла поворота 0 сечения стенки с координатой

х = 0:

У(0, 0) = (\IK)flp(x,

0)];

(147)

 

 

У(0,

t) = (\IK)f[p(x,

01;

(148)

 

0(0,

0) = (1IK)F[p(x,

0)];

(149)

 

0(0,

t) = (HK)F(p(x,

01-

(150)

Ввиду того,

что при х=0 у(х, t) = у (0, t) и у(х,

0)=г/(0, 0),

из (136) и (137)

с учетом (147) и (148) можно найти

 

Di(t) = y(0, t)EI = {EIIK)f[p{x, t)]\

D2 (t) = у (0, 0)EI = (EI/K)f[p(x, 0)].

Поскольку — = 0(x), из (136) и (137) на основании dx

и (150) следует:

Ci (0 = 0(0, t)EI = (EI/K)F[p(x, oi; С2 (0 = о (0, 0) El = (EI/К) F [p (x, 0)].

Внесение (151) и (152) в (136) и (137) дает

у(х, t) — (\IEI)\dx\M*(x, t) dx + (1//С) {F [р (х, 01* +

(151)

(149)

(152)

+f[p(x, oil;

(!ЙЗ)

у(х, 0) = (1IEI)$dx$M*{x, 0)dx + (l/K){F[p(x, 0)]* +

+ f[p(x, 0)]).

(154)

Подстановка (153), (154) и (134) в (135) приводит к урав­ нению

хН t

(EI/Hvi) §dx J dx[p(x, t) dt = j dx j [M* (x, t) —M* (x, 0)]dx + 0 x 6

+ (l/K)[x{F[p(x, 0] F [p (x, 0)}}+f[p(x, t)]— f[p(x, 0)]], (155)

заменяющему уравнение (146) в случае податливого опорного защемления стенки.

На основании уравнений (146) и (155) ниже получены зави­ симости для расчета тонких подпорных стенок на ползучих осно­ ваниях на длительную прочность.

125

§ 4. Расчет безанкерных стенок на ползучих основаниях на длительную прочность

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

Здесь первоначально рассмотрен случай, когда податливостью нижнего участка стенки, заглубленного в неползучий грунт, можно пренебречь (рис. 67). Анализ работы жесткозагцемленной стенки позволяет далее перейти к реальной схеме, когда нижний участок стенки, находящийся в подстилающем неползу­ чем грунте, обладает податливостью. Получение расчетных за­ висимостей для безанкерного защемленного больверка сводится к решению уравнения (146) при соответствующих начальных и граничных условиях. Наибольший практический интерес пред­ ставляет нахождение решений, отвечающих начальным эпюрам реактивного давления р(х, 0), используемым в современных ме­ тодах статического расчета больверков [51].

Как известно, в методе Блюма — Ломейера и в последую­ щих его модификациях эпюра р(х, 0) для связного грунта очер­ чивается по трапеции. Поэтому целесообразно получить реше­ ния уравнения (146) для тре­ угольной и прямоугольной на­ чальных эпюр р(х, о), что дает возможность вести расчет при трапецеидальной эпюре путем ее расчленения на прямоуголь­ ную и треугольную составляю­

щие.

При треугольной эпюре

р (х, 0) = рЯ (Я —х), (156)

Рис. 67. Расчетная схема безанкерной стенки на ползучем осно­ вании при неучете ее податливо­ сти в подстилающем грунте

/ — ползучий грунт; II — неползучий грунт

где р, Я и Я — соответственно объемный вес, коэффициент ре­ активного давления и мощность слоя ползучего грунта; началь­ ное и граничные условия запи­ сываются в виде

Р(х, 01 ы = Р(х, 0) = рЯ(Я—х)\

 

( 1 5 7 )

Р ( * . O L o ^ O ;

( 1 5 8 )

P(x,t) \х=н=0.

( 1 5 9 )

Решение уравнения (146) можно получить методом Фурье.

Представляя, как обычно, ис­ комую функцию р(х, t) в виде произведения двух функций

12 6

p(x, t) = ф(х)ф(^), из уравнения (146)

можно найти ф " ( х ) (/) =»

= а2ф(х)1|/(£). Разделение переменных

и введение обозначений

ф" (х)/а2ф(х) = ф/ (^)/ф(/) = —р2 приводят к двум обыкновенным дифференциальным уравнениям

ф" (х) + а2Р2ф (х) = 0; ф'.(0 + Р2ф (0 = 0,

 

решения которых имеют вид:

 

Ф (х) = A cos арх + В sin офх; ф (/) = Се~^н.

 

Тогда

 

р (х, t) = (A cosapx+ В sinaPx) Се~®н,

 

или

(160)

р (х, t) = (М cos aPx + N sin пфх) ё~

где М=АС и N = BC.

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

На основании условия (158) из соотношения (160)

следует'

Ме~$н = 0, откуда

(поскольку е~$н Ф0)

М = 0. Тогда

формула

(160) принимает вид:

 

 

 

р (х, t) = N sinaPxe-|W .

(161)

Использование

условия (159) дает

Я э т а р Я е рг = 0. Но,

поскольку ИфО и ё~^н ф 0, то з т а р Я = 0, откуда а$Н = пя или $п= пл/аН.

Здесь индекс п при найденном коэффициенте р показывает его зависимость от значения п = ± 1, 2, 3 ... Каждому значению

п отвечает свое решение уравнения (146):

2

Р п (х, t) = N n s \ n a $ nx e ^ n t .

Вследствие линейности и однородности уравнения (146) сумма частных решений также является его решением, т. е.

 

°о

 

2

 

 

р(х, t) = 2 Nn sin а$пхе

^nt •

(162)

 

П—1

 

 

 

Используя начальное условие (157), из выражения (162)

можно найти

 

 

 

 

 

00

 

 

 

рХ(Н— х) = 2

Nnsin (пях/Н).

 

 

Я=1

 

 

 

Правая часть

записанного

равенства

представляет

собой

разложение функций рЛ(Я— х)

в ряд Фурье, а Я„ — коэффи­

циент Фурье, который определяется по формуле

 

я

 

 

я

 

Яп = (2/Я )|

р(х, 0)sm(nn/H)dx = (2IH)j рХ(Н— х) X

 

о

 

 

о

 

X sin (пях/Н) dx,

127

из которой после интегрирования

Nn — 2рШ/пл.

Тогда окончательно после подстановки найденного значения Nn в формулу (162)

ОО

Н

р(х, 0 = (2 /Я )2

e-n W H '2sin (rmxlH)j р{х, 0) х

я=1

О

ОО

X sin (ппх/Н) dx = {2рШ/п) 2 (1/л) е_га2п^/а2на sin (ппх/Н) =

П=1

ОО

= (2рХН/п) 2 ( 1 /п) e -n^ B,t!W\ \ n {ппх/Н). (163)

п=1

Полученный ряд быстро сходится.

Для случая, когда начальная реактивная нагрузка распреде­ лена по закону прямоугольника, т. е. р{х, 0) =р(0) = const, гра­ ничные и начальные условия записываются следующим об­ разом:

Р (0, 0 = 0 ;

(164)

др(Н, t) = Q .

( 1 6 5 )

dx

 

Р (0, 0) = р (0).

(166)

Поскольку постоянная М в соотношении (160) по условию (164), как и ранее, равна нулю, остается без изменения и после­

дующее равенство (161).

 

(161)

по х дает

Дифференцирование выражения

■- = Ne~

cos afix,

дх

 

 

 

откуда в соответствии с условием (165)

 

a$Ne~ ^

cos а$Н = 0,

 

или

 

 

 

а$Н — (2п -+- 1) я/2;

(3„ =

(2n +

1) п/2аН.

Тогда, по аналогии с изложенным выше,

_о2 ^

Рп (х, t) = Nпе ” sinap„x;

ОООО

р{х, 0 = 2

Рп(х, 0 = 2 ЛДе_ра' sinа$пх.

(167)

п—0

я=0

 

На основании условия (166)

ОО

р (0) = 2 Nn sin а$пх.

п—0

128

Поэтому коэффициент ряда Фурье равен

н

Nn = (2/Я) J р (0) sin [(2n + 1) пх!2Н] dx=4p (0)/[(2п + 1) я],

о

Подстановка найденного значения Nn в соотношение (167) дает окончательно

р(х, 0 = [4р(0)/я] 2 [l/(2 n + l)]e~ (2n+1)2“WW2

х

п—0

00

 

 

 

X sin [(2я+ 1) пх12Н] =

[4р (0)/я] 2 [l/(2n+ 1)] х

 

п =о

 

X е~(2л'И)а я!£/"4^

sin [(2я+ 1) лх12Н].

(168)

Распространение условия р(х, 0)=р(0) на точку с коорди­ натой х = 0, используемое в методе Блюма — Ломейера и других его модификациях, является допущением. В этой связи пред­ ставляется необходимым произвести проверку устойчивости по­ лученного решения при относительно малом изменении началь­ ного условия (166), или, иными словами, проверить коррект­ ность постановки задачи. Полагая, что интенсивность начального реактивного давления снижается на малом участке длиной | от величины р(0) до нуля в точке х=0, следует принять начальные условия:

р(х 0)= Л р(0)х/^ При

I р (0) при £ < * < # .

Тогда из равенства (167) вытекает

р (0) х/1

=2 sin ар„л\

Р(0) п—О

В таком случае

Nп— (2/Я)

J р (0) (х!\) sin [(2n +

1) nxl2H] dx +

н

 

 

 

 

+ J р (0) sin [(2n -f 1) пх/2Н] dx

8р (0) Н/[(2п + I)2 я 2£] х

 

X sin [(2я-f 1) пЬ,12Н].

Поскольку при g-vO

 

 

4Р(0)

N„ — lim-

(0) Н

sin (2л +

1) я ,

£-0 (2л -f- I)2 я 2|

 

(2л + 1) я

можно констатировать, что полученное решение является ус­ тойчивым.

Наряду с выражениями (163) и (168) полезно получить упрощенные формулы для расчета безанкерного жесткозащемленного больверка на ползучем основании, более удобные для практического использования. Для этой цели при решении

129

исходного уравнения (140) можно принять упрощающее допуще­ ние о том, что

М (х , /) = Ф* [р(х,

0)]/(/);

|

М(х,

0) = ф*[р(х,

0)]/(0),

(169)

J

и поэтому

0 =

1р(х >0)]f(t)

(170)

р(х,

(где/(0) — значение, которое функция f(t) принимает в началь­ ный момент времени).

Тогда для искомой переменной можно получить выражение вида р(х, t) =\x(x)f(t), которое дает результаты, мало отличаю­ щиеся от получаемых по формулам (163) и (168). В соответст­ вии с принятым допущением уравнение (140) принимает вид

х

Н

t

X

X

J dx J

dx J p* [p (x,

0)] f (t) dt = a 2 [/ (0) f (^)] | dx J ф *[р (x , 0)] dx.

0

x

0

o

o

После двукратного интегрирования в левой и правой частях уравнения получено

р[р(*, 0)] j f(t)dt = a2 {Ц[р(х, 0)] /(0)—ф [р (х, 0)]/(*)}■

(171)

о

 

 

Дифференцирование выражения (171) поддает

 

Г (t) =

Р [ р ( х , 0)]

(172)

а 2 -ф [р (х, 0)] V

 

 

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

/(0 = Сехр {—ф [р{х,

0)] /а 2ф [р (х, 0)]).

 

(173)

Поскольку f(t)

|t=0= /(0), то

 

 

0). После внесения в урав­

нение (173) значения а2 и введения обозначения

 

 

можно найти

 

F (х) = ф [р (х,

0)] / р [р (х, 0)],

 

(174)

 

 

 

 

 

 

 

 

 

 

/(*) = /(0)ёхр { -EIt/Hi\F(x)}.

 

(175)

Подстановка соотношения

(175)

в зависимость

(170)

дает

р (х,

t) =

[3* [р (х,

0)] f (0) exp {— E l t / Hf\F (x)}.

 

На основании условия р(х, t)

|*=о= р(х, 0) можно найти, что

f(0)=p(x, 0)/р *[р (я,

0)]. Тогда окончательно

 

 

 

р (х,

t) — p(x,

0)ехр{ Elt I Fl^F (х)}.

 

(176)

Вид функции F (х)

зависит от характера начальной реактив­

ной нагрузки р(х,

0)

и выражается соотношением

 

 

{ X

X

 

 

 

N

/х

^

v

 

jdx fo*[p(x,

0)]dx

: \] d x j р* [р(JC, 0)}dx

(177)

0

0

 

 

 

J

lo

x

j

 

130

Соседние файлы в папке книги из ГПНТБ