Редактирование: КП: Динамика лавины

[[А.М. Кривцов]] > [[Теоретическая механика: физико-механический факультет|Теоретическая механика]] > [[Курсовые проекты ТМ 2015]] > '''Динамика лавины''' <HR>

[[Файл:Lav.jpg|thumb|Сход лавины|450px]]
'''''Курсовой проект по [[Теоретическая механика: физико-механический факультет|Теоретической механике]]'''''

'''Исполнитель:''' [[Сизова Екатерина]]

'''Группа:''' [[Группа 09|09]] (23604)

'''Семестр:''' весна 2015

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


== Формулировка задачи ==
* разработка модели лавины, как набор шариков, обладающих определенной массой и находящихся в неустойчивом состоянии
* расчет параметров отрыва(схода лавины), ее распространения, а также исследование динамических характеристик процесса жизненного цикла лавины


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

Наиболее благоприятными для лавинообразования являются склоны крутизной <math> 25-40^\circ</math>. Сход со склона скопившейся снежной массы обычно провоцируется климатическими причинами: резкой сменой погоды, дождями, а также механическими воздействиями на снежную массу. Иногда, в силу установившегося относительного равновесия между действующей силой трения и силой давления, сход лавины может инициироваться незначительным толчком (например, звуком ружейного выстрела или давлением на снег одного человека — горнолыжника, сноубордиста)

== Решение ==
Моделирование схода лавины происходит в двумерном пространстве. В качестве простейшей модели горы был взят прямоугольный треугольник с углом наклона <math> 30^\circ </math>, а в качестве снежной массы - элементs-круги радиуса <math>r</math>, имитирующие отдельные частицы снежной массы.

Целесообразно разделить движение частицы на 2 части: движение по наклонной плоскости и движение по горизонтальной поверхности.

При движении по наклонной плоскости шары скатываются со склона за счет силы тяжести <math> F_{mg} </math>, однако удерживаются на склоне за счет силы трения <math> F_{тр1}</math>. При данном движении меняются координаты частицы и по оси ординат, и по оси абсцисс. Сила, действующая на частицу, вычисляется по формуле <math> F = mg-Bv</math>, где <math> B </math>  - коэффициент, связывающий скорость частицы и силу.

После того, как шары скатились со склона, они продолжают движение по горизонтальной поверхности, замедляясь за счет силы трения <math> F_{тр2}</math>. Сила, действующая на частицу, считается аналогично <math> F=-B_2v</math>, однако взят другой коэффициент пропорциональности и составляющие, отвечающие за движение по вертикали взаимоуничтожают друг друга, за счёт чего происходит лишь торможение по горизонтальной оси.

Изначально частицы расположены на расстоянии, превышающим радиус частицы в 8 раз.

== Визуализация задачи==
{{#widget:Iframe |url=http://tm.spbstu.ru/htmlets/Sizova/dl_avalanche_finish.html | справа |width=900|height=500 |border=0 }}

<div class="mw-collapsible mw-collapsed">
'''Текст программы на языке JavaScript:''' <div class="mw-collapsible-content">
Файл '''"dl_avalanche_finish.js"'''
<syntaxhighlight lang="javascript" line start="1" enclose="div">
function main_dl_avalanche_finish() {

    var ctx = canvas_avalanche_finish.getContext("2d");
    var width = canvas_avalanche_finish.width;
    var height = canvas_avalanche_finish.height;

    ctx.fillStyle="#8888ff";
                           
        var vx = 2, vy = 2;                         // скорость шара
        var r = 5;                                 // радиус шара

	var grass_height = 5;
	var mountian_length = 400;

	var sin_alpha_1 = 0.5;
	var cos_alpha_1 = 0.866;

	var sun_r = 50;
	var rx = 0, ry =  mountian_length * sin_alpha_1 / cos_alpha_1 - r;
	
	var g = 2;
	var dt = 0.1;
	var m = 1;
	var B = 0.06;
	var B2 = 0.09;
	var fy = 0;
	var fx = 0;

    function step() {                           
        tick();
        draw();
    }

    function tick() {  
		if (rx + r > 0 && rx - r < mountian_length && ry < (height - grass_height - mountian_length * sin_alpha_1 / cos_alpha_1)
                                                                    * rx + mountian_length * sin_alpha_1 / cos_alpha_1) {
			fy = m * g - B * vy;
			fx = m * g - B * vx;

			vy = vy + fy / m * dt * sin_alpha_1;
			vx = vx + fx / m * dt * cos_alpha_1;

			ry = ry + vy * dt * sin_alpha_1;
			rx = rx + vx * dt * cos_alpha_1;
		}
		if (rx - r > mountian_length){

			fx = - B2 * vx;
			vx = vx + fx / m * dt;
			rx = rx + vx * dt;
		}
        
    }

    function draw() { 
		
        ctx.clearRect(0, 0, width, height);     //clear
		ctx.beginPath();
		
		ctx.fillStyle="#49b9ed"; //sky
		ctx.beginPath();
		ctx.fillRect(0, 0, width, height);
		
		ctx.fillStyle="#16a600"; //grass
		ctx.beginPath();
		ctx.fillRect(0, height - grass_height, width, height);
		
		ctx.fillStyle="#57564c"; //mountain
		ctx.beginPath();
		ctx.moveTo(0, height - grass_height);
		ctx.lineTo(mountian_length, height - grass_height);
		ctx.lineTo(0, mountian_length * sin_alpha_1 / cos_alpha_1);
		ctx.lineTo(0, height - grass_height);
		ctx.fill();
		
		
		ctx.fillStyle="#fff200"; //sun
                ctx.beginPath();
                ctx.arc(width, 0, sun_r, 0, 2 * Math.PI);
                ctx.fill();
		
		ctx.fillStyle="#e6e5d1"; //snow
		ctx.beginPath();
                ctx.arc(rx, ry, r, 0, 2 * Math.PI);
                ctx.fill();
		        
				
    }

    setInterval(step, 1000 / 60);               // функция step будет запускаться 60 раз в секунду (60 раз / 1000 мс)

}

</syntaxhighlight>

Файл '''"dl_avalanche_finish.html"'''
<syntaxhighlight lang="javascript" line start="1" enclose="div">
<!DOCTYPE html>
<html>
<head>
    <title>Avalanche</title>
    <script src="dl_avalanche_finish.js"></script>
</head>
<body onload="main_dl_avalanche_finish();">
    <canvas id="canvas_avalanche_finish" width="800" height="400" style="border: 1px solid #000000"></canvas>
</body>
</html>
</syntaxhighlight>
</div>

== Обсуждение результатов и выводы ==
В ходе работы над курсовым проектом была разработана простейшая модель, визуализирующая процесс схода лавины. Было установлено, что при увеличении коэффициента пропорциональности между силой <math> F </math>, действующей на частицу и скоростью частицы <math> v </math>  уменьшается координата конечного положения частицы. При увеличении начальной скорости <math> v </math> увеличивается координата конечного положения частицы.

Из визуализации программы видно, что после схода расстояние между частицами уменьшается по сравнению с расстоянием между начальными положениями частиц примерно в <math> 8 </math> раз, т.е. происходит "кучкование" частиц. Следовательно можно сделать вывод о том, что линейные размеры района, охваченного снежной массой, значительно меньше линейных размеров снежной массы, из которой образуется лавина.  

Дальнейшие предложения по проекту:
* Увеличить число частиц, дав возможность пользователю выбирать количество с помощью движка
* Дать возможность пользователю менять коэффициент пропорциональности с помощью движка
* Добавить движок, позволяющий менять угол наклона горы


<br>
[[Медиа:avalanche.doc|Скачать отчет]]
<br>
Скачать презентацию:

== Ссылки по теме ==
*[https://ru.wikipedia.org/wiki/%D0%9B%D0%B0%D0%B2%D0%B8%D0%BD%D0%B0  Лавины]
*[http://cyberleninka.ru/article/n/modelirovanie-shoda-snezhnoy-laviny-pri-izmenenii-temperatury-okruzhayuschey-sredy Моделирование схода снежной лавины]

== См. также ==

* [[Справка:Содержание|Справочная информация]]
* [[Курсовые проекты ТМ]]


[[Category: Студенческие проекты]]