книги / Основы прикладной геомеханики в строительстве
..pdf
|
(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 (Ь, е и) |
|