
книги из ГПНТБ / Пучков, С. В. Закономерности колебаний грунта при землетрясении
.pdfПри этих условиях выражения для смещений будут иметь вид i(fx+pt)
u = u0 e |
( 3 ,1 8 ) |
|
v = vQ e' ^ x+Pl\ |
||
|
которые совпадают с выражениями для составляющих в падающей волне. Второй случай. Положим, что мощность слоя стремится к беско
нечности, т.е. Н—*• со.
Тогда определитель системы будет приближаться к
Н
Л(£, Н) - - ( f 2-a/3) [ 4ajSf2 —(2f2—k2)2] e
и предельные формулы для смещений будут иметь следующие выраже ния:
u0a^k2(2^2- k 2) , v0if|8[ (3f2- k 2)2- f 2 ] |
i(fх , pt) |
|
|
-------------------------------------------------- е |
; |
(3 ,1 9 ) |
|
(€2--ajS) [4a/Sf2 (2£2 - k 2)2 ] |
|
|
|
2^a/3k (u()ia * v0f ) |
i(fxrpl) |
|
|
(f2-a/S) [ 4a/З^2 - (2f2 - |
k2)2 ] |
|
|
Эти выражения для смешений в волне Релея для полупространства. Третий случай. Предположим, что мощность слоя Н * Н0 мала, но
соизмерима с длиной волны, и произведения alI и /411 такие, что можно положить
shall |
ч aH0, sh /ЗН |
» /311Q и |
ch all |
ч ch /ЗН |
- 1. |
Тогда определитель системы будет приближаться к
А(£Н) - a/Sk4 + £2а/ЗН[4а2/32 + (2.f2 + k2)2 J.
Функции, выражающие числители компонент смещений, будут стре миться к
f(f,H) - - u0a/Sk4 - v0fia/3ll(4f2k2 + k4),
g(f,H) - - u 0ifa/3k2H(k2 - 2 h 2) - v0 a/3k4
40
Откуда, ограничиваясь первыми степенями Н, найдем для компонент смещения в этом случае выражения вида
u - |
• /о*/1 |
л t |
\ |
e |
i(^xi pt) |
|
u0 I iv0fH (l |
т 4 — |
) |
|
|||
|
, |
h2 |
к1 |
|
|
(3 ,2 0 ) |
v- |
)uQ . vQ]e |
i(£x+pt) |
||||
lifH U -2 — |
|
Найдем величину вектора смещения на поверхности слоя, используя модули предыдущих выражений для компонент смещений и выполняя при этом необходимые операции
Если предположим, что породы слоя удовлетворяют условию а/Ь=3, то получим выражения величины вектора смещения для отдельных част ных случаев.
Вектор смещения в продольной волне. Положим f==h, тогда
Аа ~ Уч20 + v2 + h2H2(0,lu2 , 1,96 v2)'-
Вектор смещения в поперечной волне. Соответствующее этому вол новое число £ = к- И для амплитуды смешений будем иметь выраже ние вида
Аь » У и0(1 ч 0 ,1к2Н2) , vQ(l |
i 9k2U 2)! |
|
|
|
|
|
|
Положим, |
что u2 + v2 = д2 |
и о ,11 |
и() |
= |
А0.чт<£а, |
l t96 |
v() = |
= A„cos<Aa |
для продольной волны и 0 ,1 |
u |
= |
A sin<£b, |
0 ,9 v |
= |
oco$0b для поперечной волны. Вводя эти отношения в вышенаписанные равенства, мы представим их в таком виде
Аа
С 54 <
3
- и 2 if , |
А2
а
- Ат?2 i!f 1.
(3 ,2 1 )
Ао |
Л |
где Аа и Аь .- |
амплитуды |
скальном основании, подстилающем слой рыхлых отложений, А - длина волны.
Геометрический образ полученных выражений представляет гипербо лу (рис. 1 4 ). Как видно из анализа этих выражений, при Н * О отно-
41
Р и с . 14. Зависимость амплитуд сейсмических колебаний на по верхности слоя от его мощности
шения V A0 и А0/А0 стремятся к единице, т.е. амплитуды колебаний на поверхности слоя рыхлых отложений становятся равными амплиту дам колебаний на скальном основании.
Асимптоты гипербол соответственно будут равны
Аа = 1 2гАо 7 |
’ АЬ ^ 2,,Ао Г ' |
Аа |
АЬ |
Если сравнить теоретические представления с экспериментальными данными, которые были получены Гутенбергом [7 7 ], то можно видеть, что они достаточно хорошо совпадают.
Это дает основание использовать приведенные выражения для расче та отношения амплитуд колебаний на поверхности слоя к амплитудам колебаний на его основании в зависимости от изменения мощности рых лых отложений. Последнее обеспечивает возможность надежной интер поляции и экстраполяции интенсивности колебаний на пунктах, располо женных вне и между сейсмическими станциями.
С практической точки зрения важно установить, какое предельное отношение мощности слоя к длине приходящей волны может оказывать влияние на величину смешения при сейсмических колебаниях. Тогда ус ловия третьего случая смешений в продольных и поперечных волнах позволяют дать ориентировку в этом направлении, а именно, заменяя shall и sh/ЗН углами и chall и ch/ЗН единицей, мы тем самым счи таем, что к ним применимо разложение вида
х3 |
+ • • • > |
sh х = х +7jT |
|
х2 |
|
ch х = 1 + --- |
+ . . . , |
2!
которое имеет место для |х| < 1. Следовательно, будем иметь, что для продольных волн
42
а для поперечной волны
Откуда имеем соответственно
11 |
Н |
1 |
Л |
— < |
|
А |
|
Таким образом, формулы (3 ,2 1 ) практически можно использовать, если отношение мощности слоя к длине волны меньше единицы, т.е. JJ_< 1. Вместе с тем, из этих данных видно, что предельная мощность
А
слоя, которая оказывает влияние на интенсивность сейсмических коле баний, зависит от угла падения сейсмической волны на границу слоя и основания и от отношения скоростей распространения в основании и слое. Предельное отношение мощности слоя к длине волны в зависимо
сти |
от отношения скоростей в основании и слое для углов падения Ф = |
||
= О; |
30; |
45; |
6 0 ° приведено в табл. 9. |
На рис. |
15 |
показаны кривые изменения предельных отношений мощ |
ности слоя к длине волны в зависимости от отношения скоростей рас пространения продольных волн в основании и слое для различных углов падения волны на границу слоя и основания. Эти кривые позволяют легко и быстро в каждом конкретном случае оценить влияние того или иного слоя, на котором возводится сооружение, на интенсивность сей смических колебаний.
В качестве примера возьмем среды основания и слоя, у которых отношение скоростей распространения продольных волн равно 4, а угол падения волны на границу между ними равен 60°. Тогда предельное отношение мощности слоя к длине волны по кривой 7 (см. рис. 15) будет равно 0,04.
Предельную мощность слоя можно также получить по формуле (3 ,2 1 ). Если в формулах (3 ,2 1 ) для продольных и поперечных волн примем неравенство
А
43
Т а б л и ц а 9
Отношения мощности слоя к длине волны для различных углов падения и скоростей продольных волн в основании и слое
Pi |
|
72-------- ' |
|
ф° |
|
|
|
2т7|/_1Р—cos2<£ |
3_Р —COS^lfr |
|
■cos^<f> |
9 |
|
|
|
2п |
|
|
—COS ф |
|
|
|
VPl |
|
|
|
|
1.5 |
0 |
0,14 |
0,07 |
30 |
0,13 |
0,07 |
2,0 |
0 |
0,09 |
0,05 |
30 |
0,09 |
0,05 |
2,5 |
0 |
0,07 |
0,04 |
30 |
0,07 |
0,04 |
3,0 |
0 |
0,06 |
0,03 |
30 |
0,06 |
0,03 |
3,5 |
0 |
0,05 |
0,03 |
30 |
0,05 |
0,03 |
4,0 |
0 |
0,04 |
0,02 |
30 |
0,04 |
0,02 |
5,0 |
0 |
0,03 |
0,02 |
30 |
0,03 |
0,02 |
1,5 |
45 |
0,12 |
0,06 |
60 |
0,12 |
0,06 |
2,0 |
45 |
0,08 |
0,05 |
60 |
0,08 |
0,05 |
2,5 |
45 |
0,07 |
0,04 |
60 |
0,07 |
0,04 |
3,0 |
45 |
0,06 |
0,03 |
60 |
0,06 |
0,03 |
3,5 |
45 |
0,05 |
0,03 |
60 |
0,05 |
0,03 |
4,0 |
45 |
0,04 |
0,02 |
60 |
0,04 |
0,02 |
5,0 |
45 |
0,03 |
0,02 |
60 |
0,03 |
0,02 |
Ри с. 15. Графики предельного отношения мощности слоя к длине вол ны для различного отношения скоростей и углов падения волны на грз' ницу слоя О; 30; 45 и 60°
Т а б л и ц а 10
Предельные мощности слоя для различных пород
Порода |
Скорость |
Предельная |
Скорость |
Предельная |
продольных |
МОЩНОСТЬ |
поперечных |
МОЩНОСТЬ |
|
|
волн, м/сек |
слоя Н, м |
волн, м/сек |
слоя Н, м |
Известняки |
3500 |
35 |
2060 |
41,2 |
Песчаники |
1800 |
18 |
1060 |
21,2 |
Галечники |
900 |
9 |
530 |
10,6 |
Глины |
700 |
7 |
410 |
8,2 |
Пески |
200 |
2 |
120 |
2,4 |
то предельное отношение мощности слоя к длине волны, с которого на чинает существенно сказываться его влияние, определяется следующим соотношением:
— >0,05.
А
Для различных горных пород на территории Туркмении, где скоро сти распространения волн достаточно хорошо известны, можно устано вить предельную мощность слоя, с которой начинает сказываться его влияние на величину смещения при сейсмических колебаниях.
Для этого нам потребуются экспериментальные данные о периодах сейсмических колебаний в продольной и поперечной волнах. По материа
лам сейсмических наблюдений известно, что в продольной волне |
Г—- |
= 0,2 сек, а в поперечной волне Т—= 0,4 сек. |
^ |
Предельные значения мощности слоя для различных пород, развитых |
|
на территории Туркмении, приведены в табл. 10. |
|
Можно привести несколько видоизмененные формулы для определе |
|
ния векторов смещений в продольной и поперечной волнах. |
|
Иногда удобнее пользоваться не отношением мощности слоя к длине |
|
волны, а отношением т; Г, где Т — период колебаний, а т=11 V. Тогда |
|
соответствующие выражения для смещений будут иметь вид |
|
Аа » у 0 + 4— ) + v2 (1 + 8*2 Т— 2 )•
А, |
/; |
(1 + 4 — |
) + v“ (1 + .36 |
п |
= V u2 |
||||
Ь |
г\ |
гг,9 |
О |
9 |
45
Зная толщину слоя Н и скорость распространения соответствующей сейсмической волны в нем и измеряя по сейсмограмме период колеба ний Т, можно оценить влияние слоя на амплитуду сейсмических коле баний на поверхности земли.
Предыдущие |
выводы справедливы в том случае, когда знаменатель |
|
выражения (3 ,1 |
5 ) отличен |
от нуля. Но если он обращается в нуль, |
тогда смещения становятся |
бесконечными, что на сейсмограммах обыч |
но соответствует вступлению новой фазы колебаний.
3.3.Исследование корней характеристического уравнения
Проведем анализ корней знаменателя выражения (3, 15), когда он обращается в нуль, т.е.
Д (£,Н ) = 0.
Это - уравнение частот волн Рэлея в слое.
Рассмотрим те корни, которые расположены на вещественной оси. Введем обозначения:
Е = а/3 [ + (2f2 - k2)2 ], |
|
F = £2 [ 4а2/92 + (2£2 —к2) ], |
(3,22) |
Q = 4а/3£2 (2£2 - к2). |
|
Тогда будем иметь |
|
A(f,H) = EchaHch/9H-FshaHsh/3H-Q = 0. |
(3,23) |
Для удобства исследования составим следующие соотношения:
Е - Q = а /3[4£4 + (2f2—k2)(2f2-k2—4f2)] = |
|
|
- ej8[4^4 —(2f2—к2) (2£2 + к2)] = a/3k4 >0, |
(3,24) |
|
E - F = 4 f2a/S + a/3(2£2- k 2)2—4a2j82£2—f 2(2 f2—k2) = |
|
|
= - ( f 2 - a j8 )[(2 f 2- k 2)-4 a | 3 f2]. |
(3 ,2 5 ) |
|
При этих условиях |
(3 ,2 3 ) можно записать |
|
Д(£,Н) = EchaHch/8H - [A+(f2-2a/S)R]shall sh/SH -Q = 0, |
(3 '2 6 > |
|
где |
|
|
R(f) = (2f2- K 2)2 - |
4f2a/3 |
(3 ,2 7 ) |
|
|
46
есть характеристическое уравнение частот для полупространства. Обоз начим его вещественный корень через к .
В силу симметрии Д (£,к) рассмотрим £ в интервале к,со . Значе нием корня уравнения Рэлея для полупространства оно делится на две
части: кк и |
кю ,где |
к > к. |
|
|
|
|
|||
Заметим, |
что для |
значений £ > к, R (f) < 0, |
а при |
f <к, R(£) > 0. |
|||||
Возьмем |
интервал к, |
to |
и исследуем возможность существования |
||||||
корней уравнения |
|
|
|
|
|
|
|
||
Д(£,Н) = Е ch(а—0 )Н - |
(£2-а\3) R shaHsh/ЗН - |
Q = 0. |
(3 ,2 8 ) |
||||||
При этом к<£<со и |
R(£) > 0. |
|
|
|
|
||||
Дадим произвольное значение переменной f |
в рассматриваемом ин |
||||||||
тервале, считая, что Н изменяется от О до |
. При HQ = О на основа |
||||||||
нии (3 ,2 4 ) |
имеем |
|
|
|
|
|
|
|
|
Д(£,0) = Е - |
Q > 0. |
|
|
|
|
|
|
|
|
Полагая R < О, |
уравнение (3 ,2 8 ) будет |
положительно при любом зна |
|||||||
чении Н ф О. Следовательно, |
оно в интервале |
к, со корней не имеет. |
|||||||
Перейдем к интервалу |
к < £ < к. |
|
|
|
|||||
Изучим поведение функции Д (£,Н) |
в зависимости от Н для произ |
||||||||
вольно выбранного значения f |
в заданном интервале. |
Напишем произ |
|||||||
водную этой функции |
|
|
|
|
|
|
|
||
рд |
|
|
|
|
|
|
' |
|
|
— = Ea sh a Hch/SH + E/SchaH sh /ЗН - |
[ E+ (Г-а/9)R ] (a ch aHch/SH + |
+ 0shaIich0H) = [ E (a-/3)-/S(£2_aj8)R j shaHch/3H - [ E(a-/3)
+ a(^2 -aj8)R]chaHsh/SH = [ E(a—/3) +
+ a (f2—a/3)R ] shaHch/3H [ |
- |
l.h^ H ]. |
(3 ,2 9 ) |
|
|
E(a—/3) + a(£ —aj8) R |
thaH |
|
|
При условии k < f< * |
функция R(£) > 0 |
и /3>0, |
следовательно, имеет |
|
место следующее неравенство: |
|
|
|
|
Е(а—/3) - j8(£2-a/3)R |
< 1 |
|
|
|
Е(а—/3) + а(£2 - а/S)R |
|
|
|
|
Положим, что Е(а—/3) —/3(^2—a/3)R <0. В этом случае -^-<0 и
A(f,H ) являются убывающими функциями Н. Но при Н = О
Д(£0) = Е - 0 > 0,
47
а при Н = ^
Д (0 = chaHch/3H i Е - [ Е + ( f 2-a/3)R]thaHth/3H |
Q |
I <0, |
|
chaHchjSH |
|||
|
|
так как второй член в скобках отрицательный, а третий обращается в
нуль. |
|
|
|
|
|
|
|
|
Следовательно, |
при |
к < £< к, |
A(f,H ) |
= 0 относительно Н |
и имеет |
|||
только |
один корень. |
|
|
|
|
|
||
Рассмотрим теперь следующее неравенство: |
|
|||||||
Е(а—/3) -/3(f2 - aj8) R(<f) > 0. |
|
|
|
|||||
В этом случае знак производной |
(3 ,2 9 ) |
будет определяться |
th^H thaH. |
|||||
Найдем производную этого выражения |
|
|
||||||
д_ ( th/ЗН |
_ |
1 |
/3rh2aH - a sh 2/ЗН |
|
|
|||
ЗН |
thaH |
~ |
2 |
sh2a/3ch2i3H |
|
|
|
|
|
|
(2Н )' |
|
(2Н)' |
|
|
|
|
а/32Н [ а |
~зГ~ |
- |
& 3! |
|
0 |
|
||
|
|
|
|
|
|
|
|
|
2 |
s h2 аН с h2 /ЗН |
|
|
|
при
а > 13.
Следовательно, при изменении Н от 0 до со предыдущее отношение изменяется в пределах
< 1.
a thaH
Если
Е (а~13) - р(£2~ар)Н(0 /3
E(a-/3)+a(f2-a/S)R(f)
то производная
ЭА(£ Н)
< 0.
ан
Это подтверждает, что Д(£,Н) в заданном интервале для £ имеет относительно Н только один корень. Полагая
48
Е (а-/8) - |
/3(£2 - g/3) R |
ф |
|
/3 |
|
|
|
|
|
|
|
|
|
|||||
E (а-ф) f a(£2-a/3) R ф |
|
|
a |
|
|
|
|
|
|
|
|
|
||||||
на основании |
(3 ,2 4 ) |
и (3 ,2 9 ) |
получим, |
что ----> |
О. |
|
Таким |
образом, |
||||||||||
на участке изменения мощности слоя от |
№ |
кривая А ф Н)^д воз |
||||||||||||||||
О до Hq |
||||||||||||||||||
растает, |
ибо — > |
0, |
а на участке |
Н0,м |
она убывает, |
так как— |
< 0. |
|||||||||||
Кривая |
|
да |
|
|
|
|
|
|
• |
|
|
|
|
|
да |
|
||
А ф Н) |
в функции Н при заданном значении |
£ в интервале к, |
||||||||||||||||
к показана на рис. |
16. |
|
|
|
<9Д |
|
|
|
|
|
|
|
|
|||||
Толщину слоя Н , для которого |
О, получим из соотношения |
|||||||||||||||||
= |
||||||||||||||||||
E(a-/3) - /3(£2 -a/3)R . |
thjSH |
|
|
|
|
|
|
|
|
|
||||||||
E(a-j8) I a(£2 - a/S)R |
|
thaH |
|
|
|
|
|
|
|
|
|
|||||||
Это будет соответствовать |
точке пересечения |
(рис. |
17) кривой у |
= |
||||||||||||||
= th/ЗН/ thaH |
с |
прямой |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
Е(а -р ) - /3(^2-a/8)R |
|
|
|
|
|
|
|
|
|
|
|
||||||
^ |
Е(а-/9) + a(£^~aj3)R |
|
|
|
|
|
|
|
|
|
|
|
||||||
Для любого произвольного значения £ |
в интервале |
кф<к |
можно |
|||||||||||||||
построить Л ф Н ) |
|
в |
функции |
Н и найти толщину слоя |
Н, |
при которой |
||||||||||||
Д ф Н ) обращается |
в нуль. |
Отсюда следует, что каждое значение £ |
||||||||||||||||
в этом интервале является корнем уравнения |
Л ф Н ) |
|
= О, которому от |
|||||||||||||||
вечает одно определенное значение Н . |
|
|
|
|
|
|
|
|
||||||||||
Найдем верхнее предельное |
значение |
корня уравнения |
А ф Н )= |
О |
||||||||||||||
и соответствующую величину слоя Н. Для этого напишем уравнение |
||||||||||||||||||
(3 ,2 6 ) |
в следующем виде: |
|
|
|
|
|
|
|
|
|
|
|
||||||
thaHth/3H+ |
|
|
О |
|
|
|
|
|
Е |
|
|
|
(3 ,3 0 ) |
|||||
------ ^------- |
schaHsch/9H= --------------- <1. |
|||||||||||||||||
|
|
|
Е + ф -а ,8)R |
|
|
|
Е+ф-а/ЗШ |
|
|
|
|
|||||||
Рассмотрим функцию |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
f(f,H ) = thaHth/3H + schaHsch/3H |
|
|
|
|
|
|
|
|
||||||||||
при Н = О: |
f(f,0) |
= |
1> |
|
|
|
|
|
|
|
|
|
|
|
|
|||
при Н =сл ; |
Рф**) |
= |
1. |
|
|
|
|
|
|
|
|
|
|
|
|
|||
Покажем, |
что |
f(£*H) |
при изменении О < Н |
<»» |
не превосходит |
|||||||||||||
единицы. Иначе говоря, функция f(<f,H) в |
этом |
интервале |
изменения |
|||||||||||||||
Н имеет только минимум. |
|
|
|
|
|
|
|
|
|
|
|
|||||||
Возьмем |
производную рассматриваемой функции |
|
|
|
|
|
||||||||||||
df(f,H) |
|
|
|
|
/ЗсЬаН — ach/311 |
|
|
|
|
|
|
|
ЭИ |
-(sh aH - s h /ЗИ) |
|
ch2«H ch2|8H |
||
|
49