Корзун А.С., Крайник Д.А., Горбач Н.И.,
Гурвич Ю.А.
Белорусский
национальный технический университет
Дифференциальные уравнения движения снаряда при
квадратичном законе сопротивления решены численным методом. С помощью программы
на языке Turbo Pascal 7.0
определены значения: оптимального угла наклона, при котором дальность
полета будет максимальная; двух углов, когда снаряд попадет примерно в одну и
ту же точку.
Отвлекаясь от
влияния формы снаряда и его вращения, от изменения плотности воздуха с высотой
полета снаряда, от влияния вращения Земли, скорости ветра и многих других
факторов, рассматриваемых во внешней баллистике, примем снаряд за материальную
точку .
Рассмотрим движение снаряда весом Р, которому сообщена начальная скорость под углом
α к горизонту с учётом силы сопротивления , равной по величине и направленной
по касательной к траектории снаряда, противоположно вектору скорости (рис. 1).
Рисунок
1
Составим дифференциальные уравнения движения снаряда в декартовых осях XOY:
Используя обозначения , , , , , , , получим:
Аналитическое
решение уравнений (1) и (2) невозможно, поэтому их решения выполним численным
методом.
Для определения траектории необходимо в
программе задать начальную скорость и шаг
угла М град. Примем значения коэффициента сопротивления воздуха [7].
Построим графические зависимости при следующих данных:
и − рис. 2;
и − рис. 3..
Рисунок 2
Рисунок 3
Для определения таких двух углов, при которых снаряд попадает примерно
в одну и ту же точку с величиной рассеивания 10 м, необходимо в программе
задать начальную скорость и шаг
угла . Результаты вычисления приведены в таблице
2.
Таблица
1
α |
L |
α |
L |
α |
L |
α |
L |
5 |
5387,014 |
25 |
9766,777 |
45 |
9396,144 |
65 |
6530,003 |
6 |
5916,647 |
26 |
9820,47 |
46 |
9307,752 |
66 |
6327,962 |
7 |
6375,204 |
27 |
9865,475 |
47 |
9213,507 |
67 |
6120,309 |
8 |
6777,699 |
28 |
9901,54 |
48 |
9112,959 |
68 |
5907,369 |
9 |
7132,595 |
29 |
9929,864 |
49 |
9006,668 |
69 |
5689,027 |
10 |
7450,127 |
30 |
9949,699 |
50 |
8894,666 |
70 |
5465,317 |
11 |
7734,751 |
31 |
9961,712 |
51 |
8776,543 |
71 |
5236,28 |
12 |
7992,289 |
32 |
9966,077 |
52 |
8652,798 |
72 |
5001,844 |
13 |
8225,346 |
33 |
9962,958 |
53 |
8523,249 |
73 |
4762,182 |
14 |
8436,097 |
34 |
9952,9 |
54 |
8388,14 |
74 |
4517,45 |
15 |
8627,341 |
35 |
9935,245 |
55 |
8247,105 |
75 |
4267,407 |
16 |
8801,439 |
36 |
9910,903 |
56 |
8100,573 |
76 |
4012,428 |
Продолжение таблицы 1
17 |
8959,521 |
37 |
9879,606 |
57 |
7948,38 |
77 |
3752,411 |
18 |
9103,417 |
38 |
9841,819 |
58 |
7790,558 |
78 |
3487,549 |
19 |
9232,274 |
39 |
9797,286 |
59 |
7627,139 |
79 |
3217,883 |
20 |
9349,302 |
40 |
9746,115 |
60 |
7458,148 |
80 |
2943,563 |
21 |
9454,261 |
41 |
9688,716 |
61 |
7283,46 |
81 |
2664,88 |
22 |
9547,694 |
42 |
9624,852 |
62 |
7103,411 |
82 |
2381,839 |
23 |
9630,79 |
43 |
9554,901 |
63 |
6917,722 |
83 |
2094,779 |
24 |
9703,584 |
44 |
9478,639 |
64 |
6726,543 |
84 |
1803,908 |
|
85 |
1509,535 |
Анализ таблицы 1 показал, что существует 10 пар двух различных углов , которым соответствуют 10 пар дальностей полета
снаряда (каждой паре углов – своя пара дальности с величиной рассеивания в паре
приблизительно10 м).
при углах:
, а ;
, а ;
, а .
, а ;
, а ;
, а .
, а ;
, а ;
, а .
, а .
Построим четыре пары траекторий полета
снаряда: 1 и 8; 2 и 7; 3 и 6; 4 и 5
Рисунок 4