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

книги из ГПНТБ / Пучков, С. В. Закономерности колебаний грунта при землетрясении

.pdf
Скачиваний:
17
Добавлен:
19.10.2023
Размер:
4.9 Mб
Скачать

При этих условиях выражения для смещений будут иметь вид 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

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