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

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

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

М П)

Р и с . 16. Изменение Л(Н) в зависимости от £

Р и с . 17. Определение толщины слоя

Знак производной определяется знаком множителя

вch/ЗН ^

аch пН ( ----

 

chaH

 

Когда И

изменяется

от О до со , функция ch/SH chaH монотонно убы­

вает от

единицы до

нуля. Значит при каждом отношении

существует

а

только значение HQ, при котором выражение в скобках обращается в нуль. При переходе через это значение производная fj-j С^"-Н) меняет знак с отрицательного на положительный и, следовательно, функция f(<f, Н) имеет минимум.

Величина правой части выражения (3 ,3 0 ) всегда больше значений коэффициента при втором члене левой части. Но и сама правая часть не превосходит единицы.

Таким образом, мы получим верхний предел изменения

если при­

равняем единице

правую часть (3 ,3 0 ), т.е.

 

R (f) = О,

 

 

или f = к при Н

=t's •

 

Верхнему значения корня соответствует корень характеристического

уравнения Рэлея для полупространства.

 

В итоге для интервала к<£<к

наблюдается дисперсия волн Рэлея,

скорость которых изменяется от

к до Ь.

,---- —t

Перейдем к интервалу h<f<k.

В этом случае

/З^уТ^-к будет мни­

мым и уравнение частот примет следующую форму:

50

( 2 [ 4 (f2-h2)(f2 к2) +(2f2k2)2]th H y,f2-h 2th H Vk2- f 2 +

. 4 f2\/f2- h V f 2-k2(2 f2-k2) sch H v '^ - h 2 sch H /k2- ,f2 -

- V ^ - h 2 '/ к— ^2\ 4 ^ + (2 f2 - k2)2 ] = 0.

(3 ,3 1 )

 

Это уравнение показывает, что каждому значению f теперь соответ­ ствует не одно, а множество значений Н. Для интервала 0<f<h как а, так и /3 являются мнимыми, поэтому уравнение частот можно предста­ вить в виде

£2|4(f2-h2)(<f2—k2)+(2f2—k2)2] th Н v/h2- f 2 th Н У k2-<f2 +

t 4 f2/h 2-£2 У к2- £2(2 f2~k2)sch Н y h ^ s c h

Н Ук2-<Г2 -

y/h2 -

Ук2 - ^ [ 4е2

I (2 f2 - к2)2 1 - 0.

(3 ,3 2 )

 

 

 

Как и в предыдущем случае

это уравнение для

каждого значения f

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

Для иллюстрации проведем вычисления применительно к двум физи­ чески возможным средам: слоя из несжимаемого материала и слоя рав­

ной прочности сдвигу и растяжению (гипотеза Пуассона).

 

Рассмотрим

первый случай.

Введем новые неизвестные функции

ш

и m

 

 

 

 

 

 

 

 

a il

a),

/3

ИМ,

/ЗН

1Шо

 

 

 

и напишем

 

 

 

 

 

 

 

f2 „2

■A2 i

-

т \

 

 

 

 

 

Л I ц

 

 

 

 

 

Тогда скорость распространения волн будет

 

2

к

ц

(Л : 2

) (1—ш2)

/1

 

^

£ Р A t 2ц т^ц

 

р

 

Когда материал несжимаем

Л

. В этом случае h = О и а =

а

также

 

 

 

 

 

 

 

 

<?Н

■ы. /3-■ mf, 0Н = ты

и У 2 Д 1 - т 2) - -

 

 

 

 

 

 

 

 

Р

 

51

Т а б л и ц а

11

 

 

 

Скорости распространения для различных параметров

 

ш

П

со

Л

V

Н

 

 

 

 

0 ,2 9 5 6

_

 

0

0 ,9 5 5 4

-

 

0,20

5,723

1,098

0,980

0 ,1 5

-

5,316

1,182

0 ,98 9

0 ,1 0

-

5 ,0 9 0

1,234

0 ,99 5

0

0

4,945

1,271

1,0

-

0,2

4,49 1

1,399

1,019

-

0,6

3 ,15 8

1,990

1,166

-

1.0

2,347

2 ,67 7

1,414

-

1,4

1,796

3 ,49 8

1,721

-

1,8

1,348

4,661

2,05 9

-

2,5

0,824

7,625

2,692

-

3,0

0,631

9,957

3,162

-

4,0

0,435

14,444

4,123

-

5,0

0,335

18,756

5,0 9 9

-

 

0

СО

со

П р и м еч а н и е. За единицу скорости

принята скорость поперечных волн

в слое.

 

Правая часть выражения (3 ,3 0 ),

приравненная единице, в новых пе

ременных будет иметь вид

 

(m - l)2(m3-m 2 + З т -1 ) = 0.

 

Решая это уравнение, найдем т = 0,2956, которое отвечает верхне­ му значению корня. Каждому значению т , взятому между 0,2 956 и 0, будет соответствовать только одно значение <о, полученное из ypaf

нения (3 ,3 0 ),

представленного в новых переменных:

 

4m(m2 + 1)

m(m^ + 2m2 +5)

th a thmo +

----------------. shcush mo> =

+6m2

(3 ,3 3 )

 

m^ + 6m2 + 1

+ 1

Расчеты приведены в табл. 11.

 

 

Когда £<k,

выражение д ля /3 будет мнимым:

m = in, где п имеет

действительное значение. В соответствии с этим будет

/3 - in£, £Н = со, V2 = (1 + п2 ) — ,

Р

5 2

и уравнение

(3 ,3 1 ) примет вид

 

thothno

4п(] п2)

п(п4-2п 2+5)

(3 ,3 4 )

+ -----------

shctjschncn = ---------------- . .

 

 

n4—6n2+l

n4 — 6n + 1

 

Наименьшие значения корней этого уравнения для некоторых величин

п даны в табл. 12. При п =0

уравнение

(3 ,3 4 ) переходит в следующее:

othcj + 4 schw = 5.

 

 

 

 

 

Откуда о = 4 ,9 4 5 .

 

 

 

 

 

 

Рассмотрим случай,

когда

п -»со,о

-*0

и no=const = l.

В этом слу­

чае

уравнение (3 ,3 4 )

переходит в cosl

=

О, откуда 1= (2s +1)-^-, где

•s -

целое число. На основании предыдущих формул V -*и

и £ -* О. Это,

по-видимому, характеризует переход распространяющихся волн в волны стоячие. Период колебаний этих волн выражается следующим образом:

2„

2тг

н

Т

 

 

 

 

 

V 2s

4+1 ^

и определяется в основном отношением моишости слоя к скорости рас­ пространения поперечных волн в данной среде.

Перейдем к среде, которая характеризуется гипотезой Пуассона, т.е.

А = 11 . В этом случае

будем иметь

 

3(1—т 2)'

р

3—т 2

ш2

V2

 

f 2H2

3 - т 2

Р

2

 

и характеристическое уравнение (3 ,3 0 ) перепишется в.следующем виде:

2m(5m4— 18m2 + 45)

8m(m4—7)

sh<ush m at.

th со th тсо --------------------------- -

—-----------------—

(3-m2)(m4 + 22m2+9)

(m2—3)(m4+22m2+9)

(3 ,3 5 )

Для нахождения верхнего предела изменения корня приравниваем правую часть (3 ,3 5 ) единице. В результате приходим к уравнению для определения m

(m -l)2(m4 + 12m3 + 42m2 + 36m - 27) = 0.

Положительный корень этого уравнения будет т = '0,4641, тогда ско­ рость волны

V = 0,9194

Р

5 3

Та бли ц а 12

Наименьшие значения корней уравнения (3 ,3 6 ) для различных п

ш

п

СО

А

V

 

 

 

 

0,4641

_

со

0

0,9194

0,40

-

4,43

1.19

0,942

0,30

_

3,66

1,42

0,968

0,20

3,35

1,54

0,986

0,10

 

3,21

1,60

0,997

0

0

3,145

1,62

1,0

 

0,1

3,13

1,64

1,003

 

0,3

2,87

1,76

1,03

_

0,64

2,27

2,12

1,11

0,8

2,04

2,28

1,16

_

1,0

1,78

2,49

1,22

1,5

1,33

2,92

1,35

 

2,0

1,04

3,23

1,46

 

3,0

0,72

3,56

1,58

_

4,64

0,19

9,43

1,66

СО

0

СО

N/3-

П р и м еч а н и е .

За единицу скорости принята скорость поперечных волн

 

 

в слое.

Этому отвечает со

Л

H"pf = 0.

При т =

о уравнение (3 ,3 5 ) будет иметь такой вид:

со th со +

56

,

10

----sch со

Т

 

27

 

 

 

 

откуда

3,165.

 

Параметры колебаний для различных значений ш даны в табл. 12.

В ней показано изменение отношения длин волн к толщине слоя и ско­ рости распространения волн в зависимости от поведения величины т.

При f <к /3 становится мнимым. Если примем, что

со « аН, /3 = ina,

/ЗН = incj,

 

то

 

 

 

 

 

 

 

о

2

3(l+n2)

£

 

f ^2Н2 3+п

со,

V2

р

~2

 

3+п2

 

 

 

54

и характеристическое уравнение частот перепишется в следующем виде:

2(5ш4-18ш2 + 45)

8т(ш4—7)

,

(3 ,3 6 )

thcjthmci)

 

-------------------- — sha>shma>

(3-m2)(m4+22m2+9)

(ш2-3)(ш4+22т2+9)

 

 

Если изменять m от

О до со i

то <н будет изменяться от 3,165 до О.

При стремлении п

>

О и тсо = 0 , уравнение

(3 ,3 6 )

будет

удовлетворяться.

 

 

 

 

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

нения волн

 

 

/ Д '

4

^ ■

V у 3—. £ = 0 и —

Р11

 

Известно,

что каждому значению п

отвечает множество значений

ы.

которые удовлетворяют уравнению

(3 ,3 6 ). Для некоторых значений

п

наименьшие корни уравнений представлены в таблице

12.

 

 

Если положим, что £<h<k, тогда а

и /3 будут иметь мнимые зна­

чения. В соответствии с этим напишем

 

 

 

 

аП ^ \о>, j8

= qa,

|8Н = iqa),

 

 

 

 

 

а характеристическое уравнение (3 ,3 2 )

примет вид

 

 

 

 

 

8q(q4—7)sccjscqco

2q (5q**— 18q2+45)

 

(3 ,3 7 )

 

—Igtilgqa» + ----------------------- ■=

---------- ------- ------

 

 

 

 

(q2-3 )(q 4H22q2+9)

(3 -q)(q4+22q2+9)

 

 

 

Получим выражения скорости распространения волн

 

 

 

V2 „ 3(.л! д 12

J L , ^2Н2 = 4 ^ 3

^

 

 

 

 

 

 

q2 — 3

Р

 

2

 

 

 

 

 

При этом

очевидно,

что q

меньше

у 3

быть не может.

 

 

 

При q

»V

3

уравнение

(3 ,3 7 )

имеет такой вид:

 

 

 

со.чысоку >о> -

2

 

 

 

 

 

(3 ,3 8 )

 

»

 

 

 

 

 

 

откуда наименьший корень ^ = 2,856. В этом случае

V

и £ =0,

и распространяющиеся волны вырождаются в стоя ие. Период этих волн равен

Т - -----

2»гН

-

V

- у (я2 ~ D-*>

55

Т а б л и ц а

13

 

 

 

Значение скоростей,

см /сек, при

различных

q,&>,—

 

q

 

OJ

A

V

 

 

H

 

 

 

 

 

 

 

 

 

 

 

0

CO

V3"

 

 

 

 

 

10

 

0 ,2 1 9

4,11

1,75

 

5

 

0,411

4 ,6 0 9

1,81

 

2

 

0,785

11,32

3,0

 

V T

 

2,856

со

oo

П р и м еч а н и е. За единицу скорости

принять скорость поперечных волн

 

 

в слое.

 

 

 

а для q2 = 3

 

 

 

 

Т =

2п\\

,

 

 

 

---------

 

 

 

 

v s№

 

 

 

 

где со

- один из корней

(3 ,3 8 ).

 

 

Каждому численному значению q

получим из

(3 ,3 8 ) соответствую­

щее значение со (табл. 13).

 

 

Из этой таблицы видно, что с увеличением q величина скорости V

возрастает и при q = 3

она равна бесконечности.

3.4.Траектория движения частиц на поверхности слоя

Для выявления характера движения частиц на поверхности слоя об­ ратимся к уравнениям (3 ,1 5 ). Положим, что

2п

С= — cos ф,

ивозьмем только действительную часть выражения. Тогда на поверх­ ности слоя при х = О и у ■=-Н для смещений получим следующие вы­ ражения:

V V P

и = - 6 —В (cosaH—cos/ЗН) + 9 ---- cosaH Scospt — Pi

56

V2

"о'-ЯФ

9 ( 1 ~

_ l

V"

<

^^ —

sina П I 6 — 1 /з —— —cos2<Д sin^H ■sinpt,

Л

 

Pi

 

 

 

- £

cos20

 

Pl

 

4 / V2

 

 

 

>1

 

(3 ,3 9 )

 

 

 

 

f

'

v *

v - __2- ^ 6 —

A V2

Pl

(cosaH — cos/ЗН) ч 9 —P

cosaH1 W:os>cc pt —

v *

J

“1

 

о

vp > hi

Atg<£

6 —— \ /— — — c o s ^ sinaH —

V2 \/ V2

 

V2 / V2

3 i

- 2И "№ 1

Pl \

Pi

sinpt,

V2

3 v 2 “ COs2<^

(3 ,4 0 )

где характеристическое выражение А имеет виц

при этом

 

 

 

all = 2т

Н

, /3Н = 2тт U 3 ----

н

— cosлф

—соsz<£ — ,

 

'

j 2

А

 

Pl

Pl

 

57

V

где ф - угол падения

продольной волны на основание слоя,

Vni и

-

скорости распространения продольных волн в слое и полупростран­

v p

стве,

Н -

мощность слоя

и Л - длина волны.

 

 

 

Обозначим через Р и

Q выражения, стоящие перед cospt,

и через

 

и Qj

соответствующие выражения, стоящие перед sinpt,

в

форму­

лах

(3 ,3 9 ) и (3 ,4 0 ). Используя эти данные, можно компоненты сме­

щений по осям представить в следующем виде:

 

 

 

u = uQ

J Р 2 + О2 • sin(pt t у),

 

 

 

 

 

/— ----- Г '

 

 

(3 ,4 1 )

 

v = vQ

V Pi + Oj sin(pt + У ]).

 

 

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

Нетрудно заметить резонансный характер этих выражений. Наличие Д(Н/Л,ф,Ур/Vp ) в знаменателе показывает, что при некоторых зна­ чениях переменных он обращается в нуль. На сейсмических записях это явление характеризуется колебаниями большой амплитуды с преоб­ ладающим резонансным периодом. Зная этот период Ни ^'p/Vp] , можно рассчитать отношение амплитуд колебаний на поверхности и ос­ новании слоя. В этом случае, как показывает опыт, при учете затуха­ ния оно может быть в 3 -4 раза больше, чем при обычных колебаниях нерезонансного характера.

Может быть решена и другая задача. Рассчитывая отношение амп­ литуд смещений на поверхности слоя и на его основании для различных величин Н/Л при заданных ф и Vp,' , можно найти резонансные пе­ риоды колебаний. Они, как известно, должны соответствовать макси­ мальным значениям отношений u/u0 и v. vQ.

В табл.

14 приведены резонансные периоды горизонтальных Т „ и

вертикальных

Т колебаний слоя при разных углах падения продольной

волны для

V

V = 3.

Рг j

Как видно из таблицы, резонансные периоды вертикальных и гори­ зонтальных колебаний различаются. Периоды горизонтальных колебаний несколько больше, чем периоды вертикальных. Подобное явление наблю­ дается на записях сейсмических колебаний с большой разверткой во времени. Это показывает, что жесткость слоя в вертикальном направ­ лении больше, чем в горизонтальном. Резонансные периоды медленно уменьшаются при уменьшении угла падения волны. Однако если срав­ ним периоды 1-го тона для углов падения 15 и 70°, то убедимся, что периоды горизонтальных и вертикальных колебаний уменьшаются соот­ ветственно лишь на 4,3% и 2,5%. Такую же картину можно наблюдать и для периодов 2-го тона колебаний.

58

Т а б л и ц а 14

Резонансные периоды колебаний слоя для различных углов

Угол,

град

Тон

30

15

- Т

1-й тон колеба­ ний

2-й тон колеба­ ний

Тон

- Т

1-й тон колеба­ ний

2-й тон

колеба­

ний

Тн

Т

Тн

 

Т

Н

И

н

 

Н

0,049 Vp

0,082 Vp

0,049 Vp

0 ,0 8 lV p

Н

И

И

 

И

0,149Vp

0,25Vp

0,148 Vp

0,248 Vp

 

Угол,

град

 

 

 

45

70

 

90

т

Т

т

Т

Т

н

Н

н

Н

н

0,048 Vp

0,08 Vp

0,047 Vp

0,08 Vp

0,078 V p

Н

И

И

 

II

0 ,145 V р

0,245 Vp

0,142 Vр

0,242 Vp 0.24 Vр

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

Исключая время t из выражения (3 ,4 1 ), найдем траекторию движе­ ния частиц на поверхности слоя

 

2cos (

о2)

А , (3,42)

vo (Р? Q2i uoV/p 2 + Q2> ? + Q

где f = У1 ~ У-

 

Геометрическая картина этого выражения есть эллипс, наклон кото­ рого к горизонту является функцией тех же аргументов, как и смеще­ ний и и V.

59

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