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

книги / Основы прикладной геомеханики в строительстве

..pdf
Скачиваний:
1
Добавлен:
19.11.2023
Размер:
61.85 Mб
Скачать

 

(3.19а)

а сумма главных напряжений

 

а»=(1 -ЬнО (адг+°у)-

(3.196)

Полагая, что линии фильтрационного потока параллельны оси х, получцм iy= 0, ix—sin а. Интересно отметить, что при наклоне поверхности полуплоскости к горизонту под углом а =30° /*=0,5, а интенсивность объемной фильтрационной силы будет равна 0,5ую, что составляет 100% от интенсивности объемной силы тяжести взвешенного в воде скелета.

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

(3.19в)

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

З а д а ч а о клине . Решение задач о напряженно-деформиро­ ванном состоянии в клиновидной области с учетом сил тяжести и граничных нагрузок рассмотрено в работах А. Лява (1935),

Б. Г. Галеркина (1952).

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

©(*,

y )= c lx^-\-c2x 1y‘ -\-czxy2-\-c4iyz,

(3.20)

а граничные условия

(3.146) примут вид:

 

при х —f/lg'6

 

 

 

ах cos 8 -{-хХу sin В = q xy

cos В+т^г/ sin В;

-r^cosB

sin b — - q xy

sin 8—x { y c o s

8;

при х = у tg6

ох cos 8-\-Х ху sin 8= q 2 у cos 8 — х2у sin 8;

хху cos 8 sin b = ~ - q 2 y sin 8—x2y cos 8.

Полагая, что на клин действуют также объемные силы тяжести параллельно оси у , получим компоненты напряжений на основании формул (3.15а) и (3.19):

=

[у

+■-j-

tg 8 (tj -f х2)

q2) (tg2 8 - 2)jу -

— 7

ITl _ т2+ (л —Я2) ctg 8]л:;

 

 

.-Н

+ Q 2

Т1 Ч~ ^2

+ т [ - ^ 1~

^

ctg 8(ctg2 8 “ 2)+

 

 

 

+ 1

4

4 tg Б

tg2 5 (fi—t 2)J х;

 

 

(3.20а)

 

 

 

*ху-------— [ ( ^ 1 — ^ 2) ctg 8 - J - T j —

Х2]у-\-^-~;-\-

*1 + Т2

 

 

 

 

 

 

4 tg 5

Ч\ +Я2

 

 

 

 

 

 

' )

х ■

 

 

 

 

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

 

вов

горных

пород, имею­

 

щих

треугольный

про­

 

филь, включая дамбы, на­

 

сыпи, плотины и т. д. Сле­

 

дует

отметить,

что это

 

решение не

может

быть

 

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

при б= я/2.

 

Не

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

трудно­

 

сти получить также реше­

 

ние

при несимметричном

 

относительно

оси

у

кли­

 

не

(Б. Г. Галеркин,

1952).

 

 

З а д а ч а д л я п а р а ­

 

б о л и ч е с к о г о р е л ь е ­

 

фа. В этом случае, пола­

 

гая, что уравнение конту­

 

ра

записывается

в

виде

 

у —ах2 и что на

массив

Рис. 3.8. Расчетные схемы к задаче о кли-

действуют силы

тяжести

не (а), о параболическом рельефе (о, в) и

п

квазистатические инео-

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

и-иаажлсиичеише

инер

вырезом (г)

ционные силы сеисмики,

или

ar ~ уг sin 0 ( l -----Y " ]; О0= у г sin 0; тгг0=О ,

где 0 — угол между радиусом г и осью х.

Интересно отметить, что при действии сил тяжести мы имеем во всей области растягивающие напряжения в направлении оси ху которые не зависят от координаты, в то время как остальные компоненты G y и % Х у зави­ сят линейно от коорди­

нат.

Рис. 3.9. Расчетные схемы для определения осесимметричного напряженного состояния массивов горных пород:

а — конус; б — парабола вращения; в — полу­ пространство с вертикальной выработкой; г — полупространство с полусферическим вырезом

О с е с и м м е т р и ч ­ н ы е з а д а ч и д л я п р о с т ы х к о н т у р о в (рис. 3.9). В некоторых частных случаях осе­ симметричного напря­ женного состояния уда­ ется получить простые решения. Так, напри­ мер, для напряженного состояния горных по­ род вокруг вертикаль­ ной выработки радиу­ сом г0 с учетом сил гра­ витации и гидростати­ ческого давления жид­ кости, заполняющей выработку, С. Г. Лехницким (1938) получе­ но простейшее решение*.

ar= y z

4wz

(3.22)

Исследования Д. М. Ахпателова показали; что напряженное со­ стояние осесимметричных массивов можно определить, исходя из решения соответствующей, плоской задачи. Так, например, для ко­ нуса и параболы вращения можно исходить из решений плоской задачи для клина [формула (3.20а)] и параболы [формула (3.216)]:

для конуса с углом при вершине 26 (рис. 3.9, а)

Jr = ~ - tg2S, s9 = Y Y z tg 28; oz= ± * y z ; *r2= \ yr;

(3.22a)

для параболы вращения (рис. 3.9, б)

аг

У

yz; тrz

j yr. (3.226)

8/*0

 

 

 

Аналогичным образом получено решение Д. М. Ахпателовым и И. М. Юдиной (1978) для полупространства с полусферической вы­ емкой радиусом R на поверхности:

/D

- vr sin 0 1

где 0 — угол между радиусом г и плоскостью 2 = 0.

З а д а ч и с л о ж н ы х р е л ь е ф о в . Решение задач в этом-слу­ чае может быть получено методом комплексных потенциалов, раз­ работанным Г. В. Колосовым и значительно развитым Н. И. Мусхелишвили (1935, 1966). Оказывается, что функцию напряжений <р(х, У)> удовлетворяющую бигармоническому уравнению (3.15), можно представить в виде комбинации из двух функций комплекс­ ной переменной ФДг) и <pi(2 ), называемых комплексными потен­ циалами. При этом компоненты напряжений определяются из со­ отношений

3-' + °»==4Re ® '(z) j

1

(3.23)

Компоненты напряжений в криволинейных координатах £, rj на плоскости связаны с соответствующими напряжениями в декарто­ вых координатах х, у\

°5 + ат1= 0 дг+ву5

1

(3.23а)

cri—

=

у ах ~т2/t^)e2'“,J

 

тде а — угол между

осями |

и х

(рис. 3.10), a exp (2iа) =ю'(£)

©'(£).

Если ввести обозначения <р(£) =<pi[<o(£)]; ЧЧ£) = 4 fi[o>(£)];

Ф а ) = Ф 1[а>(С)]-ф'(С) :со/ (а ;Ч гШ = № а ) ] = х1г,а ) :<**>'(£), то со­ отношения (3.23) и (3.23а) дают:

«5+ о л = 2 [Ф (С)+ ф (С)]= 4 Re Ф (0;

(3.24)

а, - «£+ 2rtn, =

№ Ф' (С)+ «' (С) V (С)].

•ЧО Сложив две строки этой системы, получим

+ Л и = ФГО + ® (С)+

НГ) Ф' (С)+ »' (С) V (С)].

(3.24а)

*Re означает действительную часть комплексного потенциала.

**Здесь и далее черта над величинами означает знак с о п р я ж е н и я .

Компоненты перемещений щ и в криволинейных координатах связаны с компонентами перемещений и и v в декартовых коорди­ натах следующим образом:

20 (в + lv) = (3 - 4р) <КС)-

~

(3'24б)

где G — модуль сдвига; р, — коэффициент Пуассона.

£ на / —

Граничные условия можно получить, если

заменить

точки на действительной оси в плоскости Z:

 

 

Л Г+гГ=Ф (*) + Ф (/) + -==[<°

 

(3.24в)

Рис* ЗЛО. Расчетная схема к применению метода комплекс­ ных потенциалов при решении плоской задачи теории упру­ гости для полубесконечных областей с криволинейной грани­ цей:

а — для плоскости Z; б — для плоскости £

где N и Т — заданные на границе значения нормального и каса­ тельного компонентов напряжений, т. е. ап и •%,. Функции Ф(£) и (О находят из граничного условия (3.24в) и ему сопряженного путем умножения на ядро Коши [2лi(t—I)]-1 и интегрирования в

пределах от —оо до + оо.

Однако такой метод нахождения комплексных потенциалов предполагает отсутствие объемных сил, т. е. однородность уравне­ ний равновесия (3.14). Это условие легко удовлетворяется, если найти какое-либо частное решение неоднородных уравнений (3.14) и представить общее решение в виде суммы частного и дополни­ тельного решений:

ау = ау + °2 ’

д

(3.25)

ху-

 

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

Если переписать граничные условия

(3.146) в виде

 

 

 

N =

ах Чу

Су

 

^ху sin 2а;

 

 

 

 

2 cos 2и

 

 

 

 

 

— ву

 

 

(3.14B )

 

 

Т =

 

 

 

 

 

~2---- sin 2а-\-хху cos 2а,

 

то с учетом формулы (3.25) получим

 

 

 

N--

аЛ +

ал

ал — ад

 

 

 

х

у

л:

^-cos2a —r*„ sin 2 a 4- - * + °y—

 

 

 

 

2

^

 

1

2

„Ч

 

u

 

 

 

 

 

 

oл:— о

 

 

 

sin 2 a,

 

 

 

 

—cos 2a —

 

 

(3.25a)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

o4 — a4

 

T = -

 

 

sin2a-]-T^cos2a

JIT

sin 2 a -f-

 

 

 

-f-t'j.

cos 2a.

 

 

 

 

 

Это дает возможность при известных значениях граничных на­ грузок N и Т определить граничные значения для дополнительных

компонентов напряжений, пользуясь

соотношениями Эйлера:

2 cos 2 a= е2г'а+ e~2t‘g; 2i sin 2a= e2i“—e-2<a,

а также зависимостью

e2ia= (o(0 :<о'(0» т. е. получим

 

 

о4 -4- о4

«(О

(3.256)

N A + iT * = N + iT ----- М

2^(0

 

 

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

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

Разработанная нами (Н. А. Цытович, 3. Г. Тер-Мартиросян и др.> 1977) обобщенная рациональная функция, позволяющая осущест­ вить конформное отображение симметричных и несимметричных полубесконечных областей с криволинейными границами, имеет вид

z==e>(C)=A j^C-f

К

Bi

з т з г ] •

<3-26>

t + a— i

(С-О2

 

 

 

где z —x+iy; £=£+^п;

 

С, а, Ь— постоянные коэффици­

енты; h — коэффициент

пропорциональности.

Разделяя

действи­

тельную и мнимую части уравнения (3.26), полагая g= f и т)=0, получим параметрические уравнения:

{A- b t){ t + a)

Ш

,

С ( < 2 - 3 ) <

(< + а)2 + 1

' (<2 + 1)2

(<2 + 1)3

Л— Ы

Д(<2—1)

 

 

. (3.26а)

■ С(<2— 1)

(< + а)2 + 1

( < 2 + 1 ) 2

"Г ( < 2 + 1 ) 3

На рис. 3.11 и 3.12 показаны симметричные и несимметричные относительно оси у полубесконечные области с криволинейной

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

.включая дамбы, котлованы, каньоны, склоны и т. п.

•Рис. 3.11. Симметричные границы полубесконечных областей, описываемые функ­ цией (3.26а) при а=Ь— 0 и различных значениях Л, В и С в относительных ко­ ординатах:

1 — А = —0.696;

В= —0,87; С=0,522; 2

Л= —0,4;

£ = —0,75;

С=0,5;

3

— А = — 0,083;

В = —0,83; С=0,5; 4

Л= —0,5;

£ = —0,3;

С=0,3;

5

Л =0;

£ = —0,6;

С=0,35; 6

Л=0,8; £=1,0;

С = —1.0;

7 —

Л=0,05;

£=0,83; С = -0 ,5 ;

8

Л= - 0 ,1 ;

£=1,8; С = —1,0;

•9 — Л= 1,6; £=2,0; С = —1,6;

10

Л = 5,0; £ = 0 ; С=1,0;

И — Л=50;

£ = 0 ; С=10

Решения

ряда

задач

при

различных

значениях

параметров

уравнения

(3.26)

были

рассмотрены

Д. М. Ахпателовым (1972),

Р. Г. Манвеляном

(1976)

и Г. Е. Шалимовым

(1977)

под руковод­

ством и при участии авторов настоящей книги. В качестве частного решения при этом были приняты напряжения, определяемые урав­ нениями (3.19). Из-за громоздкости мы здесь не приводим выра­ жения для комплексных потенциалов. Приведем лишь некоторые результаты вычислений на ЭВМ, полученные на основе этих реше­ ний; на наш взгляд, они представляют интерес для прикладной геомеханики. Результаты решения даны в относительных координа­ тах х}Н и у/Н, а напряжения — в относительных величинах.

В л и я н и е к р и в о л и н е й н о с т и р е л ь е ф а , к о э ф ф и ц и ­ е н т а б о к о в о г о д а в л е н и я и и н е р ц и о н н о й г о р и з о н ­ т а л ь н о й с и л ы (рис. 3.13—3.18). Очевидно, что характер напря­

Рис. 3.13. Изолинии

максимальных

касательных напряжений Стах в массиве горных пород при действии

сил гравитации

у (а, б, в), а также

инерционных горизонтальных сил ксу интенсивностью 0,1 у(г д, е) и 0,3у(ж, з, и)

при различных

значениях коэффициента бокового

давленния £о=0,5 (а. г, ж), |о=0,8 (б, д, з)и |о=1,0 , е и)