
книги из ГПНТБ / Пучков, С. В. Закономерности колебаний грунта при землетрясении
.pdfМ П)
Р и с . 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тг |
4Н |
н |
Т |
|
|
|
|
|
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 |
|
1Г |
|||||
|
|
|
|
||
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 |
т1н |
Т |
т |
Т |
Т |
н |
Н |
н |
Н |
н |
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