
книги из ГПНТБ / Пучков, С. В. Закономерности колебаний грунта при землетрясении
.pdf3.5. Некоторые выводы из теории колебаний упругого слоя
на жестком полупространстве под действием
плоской гармонической волны
Как было выяснено, интенсивность сейсмических колебаний на по верхности земли наряду с силой землетрясения определяется различны ми характеристиками грунтов. Существенную роль играет акустическая жесткость среды. Через нее определяется интенсивность упругой вол ны, распространяющейся вдоль сейсмического луча.
Однако в естественных условиях залегания горных пород это ослож няется влиянием угла падения сейсмической волны на основание слоя.
it |
а |
|
л |
|
|
“Т |
. |
Г |
|
|
srи 30-
е»-
!■ |
г |
*
Ри с. 18. Кривые изменения отношения амплитуд колебаний слоя в зазисимости от Н А , Vp /Vpi при Ф - 30° (а) и 70° (б)
6(Г
Вследствие этого сейсмические колебания в эпицентре будут определять ся только вертикальной составляющей, а вне его - горизонтальной со ставляющей.
Другой характеристикой является период резонансных колебаний слоя. Если период вынужденных колебаний падающей волны оказывается бли зок к периоду резонансных колебаний слоя, то на поверхности слоя возникают колебания большей амплитуды. Эти колебания часто сопро вождаются значительными деформациями в поверхностном слое грунта. Таким осложнениям подвержены особенно рыхлые песчано-глинистые отложения.
Интенсивность сейсмических колебаний зависит от мощности слоя в данном месте и длины падающей волны. Каждому отношению этих ве-
v. |
б |
и |
Ч»
V
fи АЛ
|
0,2 |
0,3 |
|
0,0 |
0/Ь |
|
|
0/00-------- |
0/00----------- |
|
• //еззатухания |
||
|
|
0/00 |
—о затуханием |
|
||
Без затухания: |
1 - А /Ап ; 2 - и /ип; |
3 - v/v |
с затуханием - |
|||
4 - А/А0 |
о |
о |
о |
|||
|
|
61
личин отвечает вполне определенное значение амплитуд сейсмических колебаний. Это помогает оценить размеры территории вокруг точки наблюдения, которая может характеризоваться интенсивностью сейсми ческих колебаний на данной сейсмической станции. Очевидно, что пло щадь этой территории определяется тем, в какой мере остается устой чивым отношение Н/л .
Для всех пунктов, где оно постоянно и равно значению под сейсми ческой станцией, интенсивности сейсмических колебаний при всех про
чих равных условь .х будут одинаковы. Там, где оно изменяется, откло няясь от значения в пункте наблюдения, интенсивность сейсмических колебаний будет другая.
Это означает, что для горизонтальных напластований отношение Н/Л будет всюду сохраняться постоянным, и, следовательно, материа лы сейсмической станции могут характеризовать большую территорию вокруг пункта наблюдения. Для пересеченной местности и неустойчиво
го рельефа отношение Н/Л будет резко меняться, и территория, примы кающая к сейсмической станции, будет незначительная.
Следовательно, в том случае, когда территория, подлежащая микро районированию, представляет ровную поверхность с горизонтальными напластованиями, тогда станции можно располагать на больших рассто яниях друг от друга. При сложном рельефе сейсмические станции не обходимо сближать, чтобы полнее охарактеризовать все сейсмические особенности сложного геологического строения района.
Интенсивность сейсмических колебаний на поверхности слоя сущест венно зависит от соотношения между скоростями распространения про дольных волн в слое и полупространстве.
Анализируя решение уравнения о колебаниях слоя на полупространст ве при различном отношении Н/Л, приходим к следующему. Если _1!___».
--- - О, то смешения на поверхности слоя мало отличаются от смещений
на его основании; если |
то колебания на поверхности слоя вы |
|
Л |
рождаются в волну Рэлея для полупространства. Смешения в этой вол не будут определять предельные значения, к которым стремятся ампли туды колебаний на поверхности слоя. Инструментальные исследования показывают, что с увеличением мощности слоя аллювия амплитуды ко лебаний на поверхности возрастают примерно в 6 -8 раз по сравнению со скалой. В промежутке между значениями П/'Л от О до с о мощность слоя мала, но соизмерима с длиной волны, возрастание амплитуд, как было показано, идет по гиперболическому закону.
Полученные данные о закономерностях сейсмических колебаний на поверхности земли (3 ,3 9 ), (3 ,4 0 ) могут быть использованы для рас чета сейсмических воздействий на строительные конструкции и массив ные сооружения.
При расчете зданий и сооружений обычно исходят из того, что зем ля - упругое однородное полупространство. В соответствии с этим дви жение грунта при землетрясениях рассматривают по упрощенной схеме гармонического движения во времени с постоянной амплитудой.
62
На самом деле подавляющее большинство зданий и сооружений воз водится на грунтах, имеющих сложное строение. Характер сейсмическо го движения в этом случае должен задаваться формулами, в которых амплитуды колебаний являются не только функциями времени, но и за висят от И/Л, Vp/Vp и ф- Графики смещений при различных значени ях этих величин представляют резонансные кривые, которые даны на рис. 18, а и б. Как видно на графиках, колебания на поверхности слоя носят резонансный характер, при котором без учета затухания происхо дит значительное увеличение амплитуд смещений. На этих же графиках приведены также кривые с учетом затухания колебаний. Они представ лены сплошной линией.
И
При отношениях — ^ 0 ,6 число резонансных пиков закономерно ра стет с увеличением разницы скоростей распространения волн в слое и
основании: при |
= |
1,5 - |
4 резонансных пика, |
а при----- = 3 - 8 . |
|
Pi |
|
|
|
Таким образом, |
чем |
больше |
различаются скорости |
распространения |
воли в основании и слое, тем чаще должны возникать резонансные яв ления при колебании на поверхности слоя.
3.6. Расчет спектров скоростей по аналитически
заданному закону движения сложного грунта
При расчетах сейсмических нагрузок ускорения или скорости поль зуются записями сильных землетрясений [3 7 ]. Но для конкретных грунтовых условий можно воспользоваться формулами (3 ,4 1 ). Для го ризонтальной и вертикальной составляющих спектра скоростей затухаю щей системы (3 ,4 1 ) перепишем так:
max
(3 ,4 3 )
vZ
шах.
Проинтегрировав в конечном виде, получим
Р
63
v . • - |
|
|
q 2/ (t ’ V |
s k |
|
|
|
|
|
( 3 , 4 4 ) |
|
|
|
|
|
|
|
|
|
|
|
||||
t 2 |
|
|
|
|
|
|
|
|
|
|
|
2n |
|
|
|
|
|
|
|
|
|
2n |
T |
In |
- 8 ■ |
|
|
|
|
|
|
|
|
|
|
|||
Г'де |
|
|
Sens y l |
I (-nr |
i 1) sin — |
•" |
Scos — l - ( — |
-l)sin — t |
|
|||
|
|
|
P |
|
I* |
|
P |
L |
P |
|
P |
|
0(T,T I) . .8) |
’> |
|
|
|
|
277 |
|
|
|
|
|
|
|
|
|
|
|
|
й2 |
I ("7jT h |
1 )' |
|
|
|
|
|
|
|
|
|
|
“t |
|
P |
' |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
2п |
|
2п |
Т |
|
2т7 |
- 8 — t |
|
|
v |
‘ |
- |
|
|
|
е |
Т |
|
||||
- 1) sin — tl |
ЙС.-OS— t 1•(— |
—Dsin— t |
|
|||||||||
|
Р |
Р |
|
т1 |
2 |
р |
|
|
> . |
|
||
|
|
|
277 |
|
|
|
|
|
|
|
||
|
|
|
й2 |
, |
( - 1 . 1 ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
т |
* D |
|
Р и с. 19. Графики функшй 0(T,Tp ,8 ) для различных периодов колебаний зданий и затухания
Тр, Т - преобладающие периоды грунта и здания; 8 - коэффициент за тухания колебаний здания. Функция 0 (Т, Тр, 8 ) _ для двух значений затухания 8 и различных периодов колебаний здания графически изобра жена на рис. 19.
Как видно из приведенных ф орм ул, спектры скброотей зависят не только от затухания и соотношения периодов колебаний здания и грун та, но также могут претерпевать изменения от одного места к друго му. Для того чтобы их использовать, необходимо учесть микрогеологические условия местности.
В некоторых работах I G7 1 поведение земляных плотин при землет рясениях уподобляется поведению рыхлого слоя, лежащего на упругом основании. Поэтому некоторые выводы этой статьи могут иметь значе ние и при расчете земляных плотин на сейсмические воздействия.
Глава 4
КОЛЕБАНИЯ ЗДАНИЙ И СООРУЖЕНИЙ ПРИ ЗЕМЛЕТРЯСЕНИЯХ И ВЗРЫВАХ
4.1. Распространение волн в здании
До сих пор развитие методики определения сейсмических нагрузок Для расчета строительных конструкций базируется на представлении, что поведение зданий можно моделировать невесомым упругим стержнем. При этом принимается, что один конец стержня жестко соединен с грун том, а на другом находится сосредоточенная нагрузка, равная весу здания [10, 20 ]• На самом деле с помощью такого моделирования нельзя описать всю совокупность деформаций реальных зданий при сей смических воздействиях [ 2 2 ].
Реальные здания надо рассматривать как упругие тела конечных размеров, в которых при землетрясениях распространяются сейсмиче ские волны. Вследствие ограниченного объема зданий эти волны отра жаются и преломляются на границах тела и интерферируют между собой. В отдельных точках здания колебания могут усиливаться, а в других ослабляться. В частности, смещения верхних этажей здания могут быть больше в несколько раз из-за отражения приходящей волны. Это приво дит к обрушению потолков внутрь здания, разрывам и выпадению углов в верхних этажах.
Измеряя скорости распространения продольных и поперечных волн в здании, можно определять динамические модули упругости и коэффици
ент |
Пуассона материала здания. Формулы для таких расчетов имеют вид |
|||
|
„ ар ( 1 + г)(1 —2 г) |
(4 ,1 ) |
||
|
Ь = ----------------- -, |
|||
|
|
|
1-г |
|
|
G = Е |
|
1 |
(4 ,2 ) |
|
2 ( 1 + г ) |
|
||
|
|
|
||
|
1 |
- |
2 ( - ) 2 |
(4 ,3 ) |
|
|
|||
|
|
|
а |
|
|
2 — |
2 ( Ь . 2 |
|
|
где |
а и Ь - |
скорости распространения |
продольных и поперечных волн |
|
в здании; |
р - плотность материала. |
|
65
Рассмотрим распространение сейсмических волн в здании. Предста вим здание в виде упругого параллелепипеда, высота которого пример но равна ширине основания, и будем считать, что основание его жестко соединено с земной поверхностью и следует за ней при колебаниях. В таком параллелепипеде могут распространяться продольные, поперечные, крутильные колебания и колебания изгиба, последние будут незначитель ными. Крутильные колебания вокруг вертикальной оси z могут возник нуть вследствие несовпадения центра тяжести с центром жесткости зда
ния.
Распространение продольных, поперечных и крутильных плоских волн в здании, как известно, можно описать дифференциальным уравнением
du |
(4 ,4 ) |
|
где и - отклонение произвольной точки параллелепипеда от положения равновесия; С - скорость распространения волн в его материале.
Примем граничные условия для верхней плоскости
du
= О
dz
и для нижней
uz=0=f(t).
При этом предполагается, что при К О f(t) = О. Примем начальные условия
du
u(z,0 ) = 0 , —----- =i/r(z).
а,1=0
На торцевых поверхностях параллелепипеда будет происходить отраже ние волн, которые будут между собой интерферировать.
Решение дифференциального уравнения, удовлетворяющего граничным и начальным условиям, можно представить в виде [1 5 , 18, 2 7 ].
z |
2 |
2nl |
z |
2п1 |
(4 ,5 ) |
u(z,t) = f (t -----) + X ( - l ) n[ f ( t - F ------- ) - f ( t + --------)]• |
|||||
c |
c |
с |
С |
с |
|
|
П= 1 |
|
|
|
|
Так как при К О |
КО = О, то ряд будет конечным для любого фикси |
||||
рованного момента |
времени. |
|
|
|
|
Рассмотрим неустановившиеся колебания. |
Положим, что движение ос |
||||
нования в горизонтальной плоскости происходит по закону синуса |
|
f(t) = Asinpt,
66
где А - амплитуда и р - круговая частота волны. Тогда движение про извольной точки параллелепипеда с координатой z можно выразить так
|
|
z |
z |
21 |
. . |
z |
21 |
u(z,t)-A[ sinp (t----) + sinp( t +-------)—sinp ( t ----------- |
|||||||
|
|
c |
c |
c |
|
c |
c |
• / |
z |
41 |
z |
41 |
|
|
|
—sinp(t |
\-------- ) + sin p(t-------—) . . . J. |
|
|
||||
|
c |
c |
c |
c |
|
|
(4 ,6 ) |
|
|
|
|
|
|
|
Это выражение представляет вынужденные колебания, происходящие с частотой р.
Найдем выражение для плотности потока энергии в грунте
Fo “ coPovo
и в здании
Fj = clP lv^ .
1
Если t<— т.е. волна не доходит до свободного конца
u(z,t) = A sinp (t |
-----), |
|
|
с |
|
F j = A2 p2 clP l cos2 p( t ----- |
). |
Максимальная плотность потока энергии равна
F 1 = A2 p2 clP l .
Для этого же промежутка времени напряжение в здании, действую щее в направлении оси z, будет
2 |
<?u |
z |
а \ “ Р\ с 1 |
^ ’ “ |
Рi ciA Pc o s p (t - - ) |
и его максимальное значение равно |
||
°1 =~ Арс1Рг. |
21 |
|
_ |
|
J1осле отражения волны от свободного конца при t < —
67
|
F = 4A2 p2 c1 p1 cos2p (----- ) |
cos2 p ( t -----), |
||||
|
1 |
r |
A A |
c |
|
c |
|
|
|
|
1—Z |
|
(4 ,7 ) |
|
= - |
|
|
|
1 . |
|
|
2c1 p1Apsinp ( — -)sin p (t— —). |
|||||
Нели |
частота |
p колебаний грунта достаточно велика, то найдутся точ- |
||||
ки z |
, |
JTC |
|
|
1- z |
п |
= 1—п , для |
которых р ------ = п — . В этом случае максимальные |
|||||
|
|
|
|
|
с |
с |
значения |
потока энергии и напряжения равны |
= 4А2 р2 с 1 р1, = - 2 Арс^рр
Таким образом, после отражения волны от свободного конца здания плотность потока энергии увеличилась в 4 раза, а напряжение и сме щение удвоились. Это подтверждается инструментальными наблюдения ми колебаний зданий [ 1 0 )•
Учтем эффект колебаний за время первого основного периода коле баний параллелепипеда. Для этого необходимо просуммировать первые
четыре слагаемых |
(4 ,6 ). |
В итоге |
получим |
|
||||||
|
16А |
9 9 |
|
2 |
|
1 |
2 |
/ "_ z \ . 2 / 2 \ |
||
|
р c-jpjcos |
|
р — sin |
|
р ( |
----) sinzp( t-------), |
||||
F l - |
|
|
|
|
С |
|
|
С |
С |
|
cti |
|
|
, 1-z |
|
1 |
|
|
, 21. |
(4 ,8 ) |
|
=4Apc1 p1 sinp(-----) sinp —cospU ----->. |
||||||||||
1 |
1 A |
c |
|
|
c |
|
|
c |
|
|
Продолжая процесс суммирования, мы убедимся, что |
||||||||||
F 1 |
= A2 p2 c1 p1 <^(p,k, z ,t ), |
|
|
|
|
crj = А р с ^ ф(p, k, z, t).
Так как A 2 p2 = F / C p и А р = aJ c d>0 » то Для плотности потока энергии и напряжения внутри здания можно записать
clPl |
clPl |
|
|
^\ —Fq------ Ф (p,k,z,t) M a p |
~ |
Ф (Pik,z,t), |
( 4 ^9 ) |
CQpQ |
0р0 |
|
|
где к - собственная частота параллелепипеда; фи ф - функции, отра жающие динамическую восприимчивость здания на колебания земной по верхности.
68
Выражение их легко получается из формулы (4 ,8 ). Для каждого здания они будут принимать свои конкретные выражения.
Формула (4 , 1 ) показывает, что интенсивность колебательного дви жения внутри здания пропорциональна силе землетрясения и функции
(р, к z,t), представляющей спектральные особенности движения в
каждой точке в любой момент времени. Она также зависит от отноше ния акустических жесткостей здания и грунта. При одной и той же акустической жесткости здания, чем прочнее грунт и, следовательно,
больше его акустическая жесткость, тем интенсивность колебательного Движения в здании будет меньше. Наоборот, слабые грунты обладают
^алой величиной |
акустической жесткости, на них сила |
сотрясений выше |
и сейсмическое |
воздействие на здание увеличивается |
[8 0 ] . |
4.2.Оценка скорости колебаний частиц,
относительных деформаций и напряжений в одинаковых зданиях на различных грунтах при прохождении сейсмической волны
Отношение |
F1 |
Fq |
Ф |
(4 ,1 0 ) |
= —— |
= v -i |
|
||
|
СТ1 |
% |
^ |
|
представляет |
собой |
скорость колебаний |
частиц в здании. В каждой точ |
|
ке здания для |
любого момента времени |
ее можно рассчитать и изме |
рить. По величине скорости можно установить относительную деформа цию путем деления ее на скорость распространения волн в здании. Зная модуль сдвига, легко рассчитать величину действующего напряжения.
Мы не знаем ни плотности потока энергии, ни действующего напря
жения в |
здании. Поэтому для определения скорости колебаний частиц |
в здании |
во время землетрясения будем исходить из отношения |
Fo |
ф |
|
|
|
' |
Ф " |
|
|
|
Рассмотрим землетрясение, |
близкое к |
9—балльному, с |
энергией Е = |
|
= 1 0 2 3 |
эрг, глубиной очага |
R = 3 0 км |
и длительностью |
колебаний |
На поверхности земли приблизительно 3 мин. Тогда плотность потока энергии будет иметь значение [4 7 ]
Е
F = --- о— = Ю кг/см • сек.
2п-В2Т
6 9