Сход снежных лавин и возникновение воздушных ударных волн.

Жекамухов М. К., Жекамухова И. М. ( Jekamuhova@mail.ru )

Кабардино-Балкарский государственный университет

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

В настоящее время в решении этой проблемы достигнуты определенные успехи, разработаны достаточно эффективные методы профилактического спуска лавин, имеются технические рекомендации по горному строительству и т.д. Однако многие явления, связанные с снежными лавинами, остаются еще недостаточно изученными. Много непонятного таит в себе разрушительное действие лавин и воздушных потоков, возникающих при сходе лавин. Для уменьшения ущерба, наносимого снежными лавинами, и обеспечения безопасности эксплуатации туристических баз и игорно-спортивных сооружений необходимо дальнейшее углубление наших знаний о закономерностях схода лавин и о процессах, связанных с лавинами.

В данной работе рассматриваются отдельные аспекты этой проблемы.

Расчет движения снежной лавины. Существует ряд моделей движения лавин. Наиболее часто пользуются моделями, в которых лавина рассматривается как материальная точка, движущаяся по наклонной плоскости с переменной массой. Эти модели различаются лишь в деталях и дают в основном представление о скоростях схода лавин и о влиянии различных факторов на величину скорости движения лавин. В дальнейшем нам понадобится зависимость скорости движения лавин от времени. Чтобы получить эту зависимость рассмотрим простейшую модель схода лавины.

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

d(mV) /. k ч cf р SV2

-dt- = mg (sinу-kfcosу)-2-, (1)

m ps X

где ps - плотность лавинного снега.

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

Xps0h0HV.

где ps0 - плотность снежного покрова на склоне горы, х - коэффициент захвата снега лавиной, h0 - мощность снежного покрова. При этом отношение

1 dm ps0 h0 V

m dt ps hX Уравнение (2) теперь принимает вид

-dV = g (sm у- fcos у)-a-, (3)

d t X

p ps0h0 где a = cf--1--0-.

В дальнейшем предположим, что отношение а^ остается постоянным при движении лавины.

Приравнивая правую часть уравнения (3) к нулю, найдем установившуюся скорость движения лавины:

Vm = J - (sin у-fcos у) V a

Введем также некоторое характерное время т, определяемое равенством Vmт = X,

X a

g sin у- f cos у

где m - масса лавины, V - скорость, t - время, g - ускорение свободного падения, р -плотность воздуха, S - площадь миделя, у - угол наклона склона горы к горизонту, kf -коэффициент трения скольжения, cf - коэффициент сопротивления. Уравнение (1) перепишем в виде

dV (. f ч рV2 S V 1 dm

- = g (sm у-fcos у)-с^ ----V---. (2)

dt 2 m m dt

Пусть X - длина лавинного снега, h - его средняя высота, Н - ширина захвата лавины. Тогда отношение

S 1

где u - скорость течения, c0 - скорость звука в воздухе, f(u) - произвольная функция, зависящая от скорости, у = 1,4 - показатель адиабаты воздуха.

В случае, когда поршень движется с постоянным ускорением, время и место образования ударной волны определяются формулами

a (у +1) a (у +1)

Тогда уравнение (3) примет вид

= 1-V2, (4)

где V = -- безразмерная скорость, t = - - безразмерное время. Vm т

Решение уравнения (4), удовлетворяющее начальному условию V(0) = 0, имеет вид

V = Vmth -. (5)

т

Отсюда видно, что скорость V стремится асимптотически к установившемуся значению Vm.

Для закона движения снежной лавины получаем формулу

x = Vm т Inch-. (6)

т

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

Как известно [1], впереди поршня, толкающего воздух в трубе, профиль волны в каждый момент времени описывается формулой

Y +1

mthJ t + f (VmthT

t Z Обозначим Vm т th - = z, тогда t = т arcth-.

T Vm

Подставляя значение t в равенство (9) и учитывая соотношение

chl:

1 - th2

z2 V 2

lnch t = -1 ln

будем иметь

c0 +-Y-A1 z т arcth-- + f(z). 0 2 J Vm W

Отсюда получим

f(u) =-VmI ln

y +1

c0 н--u I т arc th

Равенство (7) теперь принимает вид

( u2

+1 1

Vm т,

x = 1 c0 + J-u 11--- ln

2 J 2

У+1 1 , u c0 н--u тarcth-

0 2 J Vm

Время образования разрыва (ударной волны) может быть найдено из условия [1]

(dx 1 V5u J t

Оно приводит к равенству

У+It - = 0.

2 Vm

Таким образом, получаем формулы, определяющие время и место образования ударной волны

x = c0t =

2c02

(10)

т

Рассмотрим теперь случай, когда скорость поршня задается формулой (5). Поскольку равенство (7) удовлетворяется для частиц газа, непосредственно примыкающих к поршню, то, подставляя (5) и (6) в (7), получим

Vmт lnch~ = 1 c0 Vmth-

где a1

т

На рис.1 схематически показано распределение скорости воздуха между поршнем и фронтом ударной волны в произвольный момент времени t.

Vm т Inch-

Рис. 1 Распределение скорости газа перед поршнем в момент времени t.

Приведенные выше количественные формулы относятся к движению поршня внутри бесконечно длинной трубы. При движении лавин ситуация совершенно другая -лавина обтекается воздушным потоком. Следовательно, условия возникновения лавин будут более жесткими. Но если даже допустить, что лавина толкает весь воздух, находящийся впереди нее, как это имеет место в случае движения поршня в трубе, то подставляя в формулу (10) значение т = 5,5 с, получающееся при Vm = 18 м/с и Х= 100 м, получаем t = 84 с, х = 280 км; при максимально возможной скорости движения лавин, равной 50 м/с, t = 29 c, х =9 км.

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

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

т

p = p pi -ps0 lpi -ps1

(12)

где pi - плотность льда, ps1 - плотность снега и р1 - давление в ударной волне, параметр k играет роль показателя адиабаты снега. Хорошее совпадение с экспериментальными данными [2] получается, если показатель k принять равным 3,5.

Таким образом, ударная адиабата снежно-воздушной массы представлена в виде достаточно простой аналитической формулы; она справедлива до давлений, сравнимых с пределом прочности кристаллов льда, т. е. до 5 - 10 атм. При более высоких давлениях показатель k будет возрастать.

необходимость дальнейшего усовершенствования методов расчета силы удара лавин. Одним из основных моментов при этом является аналитическое представление ударной адиабаты лавинного снега. При этом могут быть построены интерполяционные формулы, пригодные лишь для отдельных типов снега и при определенных интервалах давлений. Так, например, в выше процитированной работе [2] для представления ударной адиабаты снега в аналитическом виде используется интерполяционная формула вида

p = Apk, (11)

где А - своя постоянная для каждого типа снега, а показатель k меняется в пределах от 7 до 14, в зависимости от типа снега.

При проведении конкретных расчетов с использованием формулы (11) возникают большие неудобства из-за неопределенности в выборе значений параметров А и k. Кроме того, эта формула не отражает физики процесса сжатия снега. Поэтому желательно иметь физически более обоснованное представление ударной адиабаты лавинного снега в виде достаточно простой аналитической формулы, удобной при проведении практических расчетов.

В дальнейшем будем рассматривать лавинный снег, представляющий собой снего-воздушную массу. Лавины из такого снега обычно обладают наибольшей разрушительной силой. При ударном сжатии снега, в силу быстротечности процесса, можно пренебрегать выдавливанием воздуха из снежной массы. В зоне сжатия снега можно также пренебрегать сжимаемостью кристаллов льда, образующих скелет снега. Таким образом, за фронтом ударной волны снег можно уподобить продуктам детонации ВВ, т. е. снег можно наделять свойствами сильно сжатого газа. Уравнение состояния такой среды можно получить лишь экспериментально.

Будем задавать это уравнение в виде

Рис.2. Кривые зависимости давления в ударной волне от плотности снега.

Пусть снежная лавина длиной X ударяется о неподвижное препятствие, расположенное перпендикулярно к направлению ее движения. После удара по лавинному снегу бежит отраженная ударная волна; между препятствием и ударной волной возникает застойная зона.

Законы сохранения массы и импульса на фронте ударной волны можно записать в

Ps0Ds =Ps1vDs u1 p0 +ps0Ds2 = p1 +ps1(Ds

(13)

-s U1 )2, (14)

где Ds - скорость ударной волны в снеге, u1 - массовая скорость за фронтом ударной волны.

Из теоремы сложения скоростей при сложном движении тел следует, что массовая скорость равна скорости удара: u1 = V.

Присоединяя к системе (13) - (14) уравнение ударной адиабаты (12), получим замкнутую систему уравнений относительно трех неизвестных величин : р1, p s1 и Ds.

Из равенств (13) и (14) легко получается известная формула С. А. Христиановича, выражающая зависимость между давлением и плотностью за фронтом ударной волны

ps0 ps1 v2 ps0 -ps1

На рис.2 приводятся ударные адиабаты лавинного снега при различных значениях начальной плотности снега, построенные по формуле (12); сравнение кривых на этом рисунке с аналогичными кривыми [2], полученными экспериментально, показывает, что формула (12) достаточно хорошо описывает ударную адиабату снежно-воздушной массы и ею можно пользоваться при расчете удара лавины о препятствие.

Эта формула проста и удобна при оценочных расчетах. Из равенства (13) следует

1 -ps0

(15)

а из равенств (13) и (14) получим

p1 = p0 +ps0VDs Подставляя (15) в (16), будем иметь

PL = 1,ps0V2 (

p0 p0

Из равенства (12) выразим p s1 через р1:

ps0 V ps1 J

(16)

(17)

ps1 =pi -(pi -ps0 f -

(18)

Подставляя (18) в (17), получим уравнение, связывающее V и р1

-V.=( pl -1

Vs Vps0 ,

-1 + ps0

(19)

где $ = -p1;Vs = p0 \

£0

постоянная величина, имеющая размерность скорости,

Vs 10,5 м/с.

Подставляя значения p s1 и V в (15), будем иметь

= pl

Vs Vps0

ps0 Vpi -ps0

-1 + ps0

(20)

Скорость распространения звука в снежной массе, сжатой ударной волной, определяется формулой

1 k+1

dp1 12 с 2k

V d ps1 J

cs0 S

(21)

где cs0

Vpi -ps0

где

На рис.3 приводятся кривые зависимости относительного давления

A p = p1 - p0, от скорости удара лавины, построенные по формуле (19) при различных значениях плотности снега до удара.

Рис. 3 Зависимость давления в ударной волне от скорости удара лавины о препятствие.

Кривые 1, 2, 3, 4 - для ps = 0.1, 0.2, 0.3, 0.4 кг/м3 соответственно.

Из расчетов следует, что давления в зоне застоя, рассчитанные по формуле (19) при различных скоростях удара, превосходят динамические напоры, рассчитанные по

известной формуле p = pS0V , которой часто пользуются при расчетах противолавинных

сооружений, примерно в 2 - 4 раза в интервале плотностей снега от ps0 = 200 кг/м3 до ps0 = 400 кг/м3.

На рис.4 приводятся кривые зависимости Ds от скорости удара лавины,

построенные по формуле

(20) при

значениях ps0.

различных

Рис. 4 Зависимость скорости ударной волны в снеге от скорости удара лавины о препятствие.

Заметим, что формулу (12) нельзя экстраполировать к малым амплитудам р1 - 0, при которых скорость звука в снеге определяется в основном жесткостью снежного каркаса, так, что cs0 не есть скорость звука в снеге при обычном его состоянии.

Таким образом, все параметры снежно-воздушной массы за фронтом ударной волны выражаются через давление.

Электронный журнал ИССЛЕДОВАНО в РОССИИ 52 1 http: zhurnal.ape.relarn.ru/articles/2003/044.pdf

Кривые 1, 2, 3, 4 - для ps = 0.1, 0.2, 0.3, 0.4 кг/м3 соответственно.

Эффект действия лавины на неподвижное препятствие определяется величиной ударного импульса K = p-At, где At - длительность удара лавины о препятствие. Величина

At состоит из двух слагаемых: At = At1 + At2, где A ti = yD - промежуток времени, в

течение которого ударная волна, отраженная от препятствия, достигает хвостовой части

лавины, At2 = у - промежуток времени, в течение которого волна разгрузки бежит по /cs1

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

1 + J Ds cs1.

(22)

Зная скорость удара V, по графикам 4 могут быть определены соответствующие значения Ds, а затем по формуле (22) и ударный импульс.

Воздушные ударные волны, возникающие при ударе снежной лавины о препятствие. Обширный материал наблюдений за сходом снежных лавин показывает, что сход мощных лавин часто сопровождается возникновением воздушных ударных волн. При этом имеются случаи, когда лесоповалы и другие разрушения происходят не впереди лавины, а в стороне от нее [3]. Эти разрушения уже не связаны с воздушными потоками, образующимися вокруг лавины, а являются следствием действия ударных волн. Однако до настоящего времени механизм образования этих волн остается неясным. Ниже предпринята попытка установления механизма возникновения воздушных ударных волн при ударе лавин о препятствие.

Будем исходить из упрощенной схемы удара лавины о неподвижное препятствие, рассмотренной в предыдущем пункте. После выхода ударной волны на тыльную поверхность лавины начинается процесс разгрузки сжатой снежно-воздушной массы, происходит разгон поверхностного слоя снега. При этом слой воздуха, непосредственно примыкающий к поверхности снега, получает резкий толчок (импульс), который передается окружающей воздушной среде в виде ударной волны. Движущаяся поверхность снега играет роль поршня, который толкает воздух впереди себя.

На рис.5 схематически показана разгрузка снега в ударной волне.

т