Рис. 1. График изменения высоты и скорости движения при свободном падении тела в воздухе

 

Результаты построения модели движения каплевидного тела в аналогичных условиях приведены в приложении 2. При проведении расчетов учтено, что коэффициент лобового сопротивления для капли с=0,045 и площадь поперечного сечения уменьшается при вытягивании «хвостика» капли. Чем более вытянутой является капля, тем меньше сопротивление воздуха и выше скорость падения. На рис. 2 представлен график разности скоростей для капель одинаковой массы, но с разным поперечным сечением.

Рис. 2. Зависимость разности скоростей падающих капель разного сечения при свободном падении в воздухе

 

 

3.2. Построение модели падения свинцового шарика в воде с учетом сопротивления среды

 

Решим в MS Excel задачу 2 п. 2.3 на том же рабочем листе, что и задачу 1. Характеристики шарика заданы в таблице 1. Для расчета требуются плотность и динамическая вязкость воды (справочные данные) и значения коэффициентов сопротивления, вычисленные в задаче 2:

Таблица 5.

Плотность и динамическая вязкость воды, коэффициенты сопротивления

 

F

G

H

6

вода

7

μ=

1,002

Нс/м^2

8

ρ_=

1000

кг/м^3

9

k1=

0,1889

Нс/м

10

k2=

0,0628

кг c/м

 

Так как падение в воде происходит медленнее, то начальную высоту зададим 2м и завершим вычисления, когда движение станет равномерным (изменение скорости составляет менее 0,0001 за один шаг вычислений). Проведем расчеты скорости и высоты тела по формулам (9) и (10) аналогично задаче 1 (см. таблицу 6). Для проведения расчетов дальше нужно выделить диапазон ячеек G14:H14 и протянуть выделение вниз для копирования формул (как и в предыдущей задаче). Завершить протягивание нужно при появлении в столбце I надписи «торможение». Приведем в таблице 6 часть расчетных формул, а в таблице 7 – вычисленные значения. Полностью результаты расчетов для этой задачи приведены в приложении 1.

 

Таблица 6.

Формулы для расчета высоты и скорости падения в воде

 

F

G

H

I

12

t

h

v

13

0

=$B3

=G4

14

0,1

=G13-H14*(F14-F13)

=H13+($H$2*$B$2-G$9*H13-G$10*H13*H13)/$H$2*(F14-F13)

=ЕСЛИ(H14-H13<=0,0001; "торможение";""

15

0,2

=G14-H15*(F15-F14)

=H14+($H$2*$B$2-G$9*H14-G$10*H14*H14)/$H$2*(F15-F14)

=ЕСЛИ(H15-H14<=0,0001; "торможение";""

 

Таблица 7.

Расчет высоты и скорости падения в воде

 

F

G

H

I

12

t

h

v

13

0

2

0

14

0,05

1,976

0,490

15

0,1

1,932

0,867

16

0,15

1,875

1,135

17

0,2

1,810

1,314

18

0,25

1,738

1,429

19

0,3

1,663

1,500

20

0,35

1,586

1,543

21

0,4

1,508

1,569

22

0,45

1,429

1,584

23

0,5

1,349

1,593

24

0,55

1,269

1,599

25

0,6

1,189

1,602

26

0,65

1,109

1,604

27

0,7

1,028

1,605

28

0,75

0,948

1,605

29

0,8

0,868

1,606

30

0,85

0,788

1,606

31

0,9

0,707

1,606

32

0,95

0,627

1,606

торможение

33

1

0,547

1,606

торможение

 

Результаты проведения расчетов представим графически. Для этого вставим на рабочий лист диаграмму (график) для рядов данных G12:G33 (высота тела) и H12:H33 (скорость падения), в качестве подписей по оси OX используем значения времени из диапазона F12:F33. Результаты представлены на рис. 3. Из рисунка видно, что скорость погружения сначала возрастает, затем изменение скорости замедляется и движение становится практически равномерным. Изменение высоты тела при этом быстро становится линейным.