Реферат: Компьютерное моделирование движения тел
NEXT
END
Данный блок программы используя массив (DIM) случайных чисел ограниченных пространством ( r ) и количеством (n) точек позволяет создать эффектный фейерверк разлета осколков.
Полностью программа представлена в Приложении 2.
Данная задача рассматривалась без учета сопротивления воздуха. Попробуем ввести поправку на сопротивление воздуха.
4. Этап четвертый. Воздействие внешних факторов.
К внешним факторам в решении задачи можно отнести случайные события, которые происходят в природе или при протекании каких либо физических процессов. В нашем случае это может быть сопротивление воздуха, порыв ветра, интенсивность горения пороха, влажность и т.д. Из всех случайных величин рассмотрим только сопротивление воздуха.
Учитывая то, что сопротивление воздуха направленно в противоположную сторону движения тела [4], то силу сопротивления можно представить в виде двух составляющих: вертикальной и горизонтальной, при этом, чем выше скорость, тем больше сопротивление и наоборот.
Выделим блок программы (Приложение 2) с 100-го кадра:
100 q = .0005
t = t + q
s = v * t * COS(f) * 2 (4.1)
h = (v * t * SIN(f) – (g * t ^ 2) / 2) * 5.5 (4.2)
Запишем составляющие горизонтальной и вертикальной скорости:
vs=v*cos(f)
vh=v*sin(f)
Но если пренебречь, тем что с высотой плотность воздуха падает и считать ее в нижних слоях атмосферы неизменной, то упрощенная формула горизонтальной и вертикальной составляющей, в соответствии с законом Аэродинамики [ 4 ], примет вид [ 5 ]:
Тогда выражение 4.1 и 4.2 можно изменить записав блок с 100-го кадра в следующем виде:
100 q = .0005
t = t + q
k = .0002
vs = v * COS(f) - (k - k / (k + v * COS(f)))
vh = v * SIN(f) - (k - k / (k + v * SIN(f)))
s = vs * t * 2
h = (vh * t - (g * t ^ 2) / 2) * 5.5
v = SQR(vs ^ 2 + vh ^ 2) (4.3)
LINE (x2, y2)-(x0 + s / 20, y0 - h / 25), 14
x2 = x0 + s / 20