CC BY
79
ASSESSMENT OF THE SINGLE SHOT HIT PROBABILITY AS A FUNCTION OF THE HORIZONTAL RANGE TAKING INTO ACCOUNT DIFFERENT TARGET TYPES AND POINTS OF AIM
s>
Vadim L. Khaikov
independent researcher, Krasnodar, Russian Federation, e-mail: wadimhaikow@inbox.ru, ORCID iD: http://orcid.org/0000-0003-1433-3562
DOI: 10.5937/vojtehg67-18522; https://doi.org/10.5937/vojtehg67-18522
FIELD: Applied Mathematics, Operations Research ARTICLE TYPE: Original Scientific Paper <o
ARTICLE LANGUAGE: English
Abstract:
.s
o
SS
This article presents a method for assessing the hit probability for stationary shooting targets as a function of the projectile horizontal range | (PHR), the target type and the point of aim (POA) location. The proposed l procedure consists of three blocks. Block I is devoted to the double integral formula taking into account changes in the point of impact (POI) location as a function of the PHR. The characteristics of bullet accuracy and precision versus the PHR are determined in block II. The basis of blocks I and II is the bi-variate uncorrelated Laplace-Gauss probability distribution. The functions of the POI location (ordinate/abscissa) versus the PHR and the functions of the characteristics of the bullet dispersion versus the PHR are represented in the form of polynomials. The description of the target silhouette contour is given in block III. Mathematically, the target contour is a piecewise function which defines the upper and lower edges of the shooting target and it also represents the limit for the double integral formula of block I. The proposed method is built on a modular basis and allows a user to change types of weapons and shooting targets. For demo calculations, the accuracy and precision characteristics of the 5.45*39 Kalashnikov assault rifle model MPi AK-74N were selected. Five types of Swiss military targets imitating an OPFOR combatant's silhouette were used as shooting targets. For illustrating the em operability of the proposed method, the hit probabilities for the K, H, G, F, and E shooting targets were evaluated for the PHR from 50 to 400 m. All | computations are implemented in the PTC Mathcad v. 15. >
Key words: hit probability, numerical solution, shooting target, double integrals, bullet dispersion, horizontal range, Mathcad.
¡8
<D
CD
O >
CD
0
01
of
UJ Cd ZD O o
_J <
o
X
o
LU
I— >-
a: <
I—
CO <
-J
CD >o
X LU I—
o
o >
Introduction
Estimation of the hit probability is an important part of the shooting effectiveness assessment (Rodney, 2012), (Svateev, 2013), (Svateev, 2014), (Peelen, 2017), (Wollschlager, 2017). In (Khaikov, 2018), we show how to estimate the single shot hit probability by shooting to a fixed target with an arbitrary contour form owing to a numerical solution of the double integral. During the calculation process, the following parameters are set: the position of the mean point of impact (POI) and the characteristics of the bullet dispersion ellipse. However, during firing, there is uncertainty about the horizontal range to the selected target; it is also necessary to have a priori information about the change in the hit probability for the interval from the minimum to the maximum distances.
Since the position of the dispersion center and the dispersion characteristics are the functions of the projectile horizontal range (PHR), the hit probability can be expressed as a function of the «target-shooter» distance1.
The aim of this contribution is to show how to calculate the hit probability (Phit) as a function of the PHR with respect to different target types and the location of the POA. As a result, we get not only one (or a point) estimation of the Phit for a single (selected) value of the horizontal range, but also a change of this parameter within a certain distance interval. As in (Khaikov, 2018), all calculations will be implemented in the computer algebra system (CAS) Mathcad v.15.
In order to make the estimations original and more diverse, we will use here Swiss military targets and the external ballistics of the 5.45*39 Kalashnikov assault rifle model MPi AK-74N2. The information about this firearm can be found in (NVA der DDR, 1985), (TRADOC, 1975, 2015), (http://weaponland.ru, nd).
Creating a mathematical model
Let us assume that the attacking opposing force (OPFOR) is represented by running or lying combatants whose silhouettes can be replaced by plane (or 2D) shooting targets with a certain equivalent area3
1 In this article, the expressions "the target-shooter distance" and "the projectile horizontal range" are synonymous.
2 MPi AK-74N - (in German) Mashinenpistole-Automat Kalaschnikow AK-74 für Nachtsichtgerat. Assault rifles MPi AK-74N were manufactured on "VEB Fahrzeug und Waffenfabrik Ernst Thälmann", Suhl (Thuringia).
3 As a rule, shooting targets are a set of geometric primitives: rectangles, trapezoids, triangles, and squares.
(Fig. 1a). Human silhouette targets of different types are used for military and police firearms training. More schematized targets are applied for firearm sports competitions.
Let it be necessary to estimate the hit probability for several targets at the distances xi, xn, and xm (Fig. 1b). All shooting targets at different distances are the same and in Fig. 1b they are rotated for 90 degrees counter-clockwise. The targets are rotated only for the simplicity of visualization.
<D Cp
Sy
S
e .s
a)
b)
Figure 1 - Replacement of combatant silhouettes with shooting targets (a) and the change in the POI location and the characteristics of the dispersion ellipse as a
function of the PHR (b) Рис. 1 - Замена силуэтов комбатантов в виде стрелковых мишеней (a) и изменение местоположения СТП и характеристик эллипса рассеивания в зависимости от горизонтальной дальности «стрелок - мишень» (b) Слика 1 - Замена силуета во]ника у облику мета за гамаше (а) и промена локаци^а тачака погодака (ТП) и карактеристика елипсе дисперзи^е у зависности од хоризонталног расто]аша измену стрелца и мете (ХРМ) (б)
When firing at a shooting target, depending on the PHR, the mean bullet (ballistic) trajectory4 will pass through the target or above (below) it (Fig. 1). The difference in the PHR between the shooter and the shooting target (xI, xII, xIII) determines the difference in the POI location and the constantly increasing parameters of the dispersion ellipse for each target. The location of the POA in Fig. 1 is not indicated.
4 This trajectory is a projectile path obtained by theoretical calculation and agreed upon with experimental shooting tests. In this way, the mean ballistic trajectory (MBT) is the space curve near which the remaining trajectories are grouped taking into account the regularity of bullet dispersion.
.q ra
-Q
p
E
Я
I >
15
Let us define that the shooting mode (single shots) and the aiming parameters are not changed by shooting. The ABCD curve (Fig. 1b) is the mean ballistic trajectory (MBT) with respect to the specified angle of £ departure (or number/clicks of sight) and the previously specified zero ° range. Uncertainty in the PHR magnitude is transformed into a vertical projectile dispersion at a chosen shooting target. Naturally, the hit o probability value (Phtt i, Phit ii, Phtt iii) will change with the variations of the 2,R POI location and with the change in the dispersion ellipse ^ characteristics5.
UR The mathematical model of the hit-probability vs the PHR will
o consist of three blocks. Block I will be devoted to the double integral
< formula taking into account the changes in the POI location and bullet ° dispersion. Block II will include the bullet dispersion characteristics vs the g horizontal range. Finally, block III will describe the functions of a target TE silhouette contour.
a:
< Block I. This block of the mathematical model is based on the results of (Khaikov, 2018). If the T region (letter "T means a target) is a plane figure, and the hit probability of random variables Z and Y is determined by the 2D uncorrelated Laplace-Gauss distribution f (z, y),
< then the target hit probability is defined by formula (1) (Venttsel', 2006,
o p.196)
P[(Z, 7)c r] = jj f(z, y)dzdy (1)
>
X LLI I—
o
g where: T - the region of integration in the zy-plane; z, y - the Cartesian coordinates; a , a - the standard deviation of the continuous random
z 5 y
variables Z, Y; mz, m - the coordinates of the center of dispersion
(the mean of Z and Y) or the POI, since it is the same one; f (z, y) - the integrand and the uncorrelated bi-variate Laplace-Gauss distribution (2):
5 The first important factor is the type of weapon and the type of cartridge used. In the second place, the POI location depends on the target type. For the pre-selected weapon, the cartridge and the shooting target, the POI location and the characteristics of the dispersion ellipse are the functions of the PHR. In the case under consideration, skills of the shooter and combat stress as psychophysical conditions are not taken into account.
f (z, y ) =
2n<jzay
exp
(z - mz )2 (y - my )2 " 2<2
2<
(2)
In (Khaikov, 2018, pp.739-756), it is shown that if the target silhouette is described by two functions (the function of the bottom contour line of the target $1(z) and the function of the upper line of the target's contour) ^(z) and two constants a and b (a<b), then the hit probability can be expressed as
b %(z)
P[(Z )] = J J
1
7 % 2/K<J< „
% (z) z y
exp
(z - mz )2 (y - my } " 2<2
2<
2
dzdy.
(3)
If in formula (3) the constants mz, my, az, ay are transformed into functions that depend on the distance x "shooter - target" (mz (x), ,(x), &z(x), ay(x)), we get a formula that realizes the transition from
m
p[(z, y)] to p[(z, y, x)]
b %(z)
P[(Z,Y,X)]=J J
a % (z) 2
2< (x)<y (x)
exp
(z - mz (x))2 (y - my (x))2
+ ■
2<z (x)2 2<y (x)
dzdy.
(4)
The method of determining the functionsmz(x), my(x), &z(x), and ay(x) for the assault rifle MPi AK-74N will be described in block II.
Block II. Characteristics of mz(x), my(x), az(x), ay(x) functions.
The trajectories of the steel core bullet of the 5.45*39 Kalashnikov cartridge (bullet weight - 3.4 g; muzzle velocity - 960 m/s) under the line of departure and over/under the line of sight for mechanical sight from 1 to 10 are described in (NVA der DDR, 1985, pp.90-93) as tabulated data. Let us visualize this table and obtain the approximation formula (Fig. 2).
The data corresponding to sight number 3 (NVA der DDR, 1985, pp.90-93) can be approximated using a polynomial. As a result, we obtain the function
mv =-0.029 +1.969• 103 x - 2.645• 10-6 x2 -1.192• 10-8 x3,
(5)
1
a
1
CO
О >
0
01
of
Ш
a:
ZD О О
_J
<
О 2:
X
о ш
I—
>-
Ql <
(Л
с
-J
CD >0
x
Ш I—
О
О >
where x is the PHR («shooter-target» distance) in meters (50 < x < 400 m) and my is the ballistic trajectory profile in meters for the sight setting 3 (position 3 for standard mechanical sight of the Kalashnikov assault rifle). The red curve in Fig. 2 shows the change of the POI height. At 50 m, shots will be about 6 centimeters=0.06 m high. The vertex of the flight path is 27.66 cm = 0.267 m, observed for the PHR 168.3 m (Fig. 2a). The zero distance is 300 m. Beyond 300 m (PHR>300 m), the bullet trajectory will drop rapidly to about 43 centimeters « 0.427 m low at 400 m.
a)
b)
Figure 2 - Characteristics of the location of the center of dispersion (a) and the
characteristics of the dispersion ellipse (b) as a function of the PHR Рис. 2 - Характеристики положения центра рассеивания (a) и эллипса рассеивания (b) как функция дистанции до цели Слика 2 - Карактеристике положа^а центра дисперзи^е (а) и карактеристике елипсе дисперзи^е (б) у функщи расто]а^а до мете
The right-hand bullet drift or the bias of the POI to the right side due to the bullet gyroscopic stabilization can be expressed by the polynomial model
m.
= -1.667 • 105 x +1.667 • 107 x2
(6)
where x is the PHR in meters (50 < x < 400 m) and mz is the right-hand bullet drift in meters (the blue dotted curve in Fig. 2a).
The bullets dispersion in the vertical (ct ) and the lateral directions
(ctz ) in the form of standard deviation are expressed as two polynomial functions PHR-argument (Fig. 2b)
ctv = 5.935 • 104 x,
18
<7г = 2.967 • 10"4 л,
(8)
where x is the PHR in meters (50 < x < 400 m), and a and az represent the standard deviation in meters.
Thus, two functions (a (x) and az(x)) determine the shape of the
dispersion ellipse and two more, i.e. m (x), mz(x), determine the
locations of the dispersion center versus the PHR from 50 to 400 m.
In order to use formula (4), we need to determine the functions $1(z), $2(z) and the constants a and b.
Block III. Functions of the target silhouette contour $1(z) and $2(z). Depending on an OPFOR combatant silhouette seen by a shooter, there are «head», «head-shoulders», «head-burst», «head-hip», and «man sized» military targets (Fig. 3)6.
Ш CP
.S
Figure 3 - Military shooting targets (cibles de campagne): K, H, G, F, and E Рис. 3 - Армейские мишени K, H, G, F, E для стрельб Слика 3 - Во]не мете за гамаше K, H, G, F и E
Е
.s
-Q
га
.о
р
Among the five shooting targets, the K and H targets have the minimum height of 33 centimeters (0.33 m); conversely, the E target has the maximum height of 1.65 m. The target area is increased in the
6 in comparison with Fig. 1, the number of types of targets in Fig. 3 was increased from 3 to 5. In different countries, shooting targets of the same type (for example, the K target «head in a protective helmet» or the E target «man sized target») may have different shapes and different geometric areas. These characteristics are determined in individual countries by their armed forces regulations.
< >
19
direction from K to E. The K, H, G, F, and E human silhouette targets7 are used in the armed forces of the Swiss Confederation for shooting training (http://www.wikiwand.com, nd).
The target contour can be designated in two ways (Khaikov, 2018, pp.739-756): in option I, the upper (lower) contour should be identified as the y(z) functions or, in option II, the right (left) contour should be identified as the z(y) curves. For all five targets in Fig. 3, the upper contour lines (UCL) are red, and the bottom counter lines (BCL) are yellow. If we consider the UCL of the target K, then we use the notation UCLk. This procedure for the BCL of the target E is similar: BCLe.
The geometry analysis and the contour lines identification for the G target are shown in Fig. 4. It is a flat shooting target which can be described as the ABCDEFGHIJ polygon. The G target has an axis of symmetry that passes through the EJ (Fig. 4a). The left side of the G target is the ABCDEJ polygon, and the right-hand side is the EFGHIJ polygon. The areas of these polygons are the same. The presence of the right-hand bullet drift causes the POI to be located on the right side (the red ellipse in Fig. 4b ).
a)
b)
Figure 4 - Analysis of the geometry of the G shooting target (a). A location of the
dispersion ellipse (b) Рис. 4 - Анализ геометрии стрелковой мишени G (a). Положение эллипса
рассеивания (b)
Слика 4 - Анализа геометри^е мете за гамаше G (а). Локацща елипсе дисперзще
(б)
7 The following names are used for target designation (in German): Feldscheiben, Ordonanzscheiben, K-Sheibe, G-Sheibe, H-Sheibe, and so on.
Not taking into account the symmetry of the G target, its UCL can be described as the BCDEFGH line. When symmetry is taken into account, the geometric formula of the left-hand side is bCdE. Similar reasoning applies to the description of the BCL: AJI (without symmetry) and AJ (the left side of the figure with symmetry taken into account).
It is necessary to determine the coordinates of the vertices of the ABCDEFGHIJ polygon and to obtain an analytical expression of the sides of the G target. The coordinates of the points for the G target contour and its analytic functions are collected in Table 1.
Table 1 - Coordinates of the points for the G target contour line Таблица 1 - Координаты точек контура мишени типа G Табела 1 - Координате тачака контуре мете типа G
Points Z, m Y, m Boundary points of line Analytic function
A -0.225 0.0 A-B z(y)=-0.225
B -0.225 0.302 B-C y(z)= 0.384*z + 0.388
C -0.1 0.35 C-D z(y)=-0.1
D -0.1 0.55 D-E y(z)=0.55
E 0.0 0.55 E-F y(z)=0.55
F 0.1 0.55 F-G z(y)=0.1
G 0.1 0.35 G-H y(z)= -0.384*z + 0.388
H 0.225 0.302 H-I z(y)= 0.225
I 0.225 0.0 I-J y(z)= 0
J 0.0 0.0 J-A y(z)= 0
Notes: The coordinates of the points are given according to Fig. 3. The symmetry of points is shown with the orange and blue background. The locations of the vertices for the G target can be used in the analysis of the H target and the F target.
The contour line (CL) CLg(z) of the G shooting target (the UCL (BCDEFGH) UCLg(z) and the BCL (AJI) BCLg(z) ) can be described by a piecewise function:
CLg (z ) =
( 0.384z + 0.388
- 0.225 < z <-0.1^
UCLg(z) = 0.55 - 0.1 < z < 0.1
- 0.384z + 0.388 0.1 < z < 0.225
ч
BCLg (z) = (0 - 0.225 < z < 0.225)
Since the geometry of the G target has an axis of symmetry (EJ axis) and if the POI lies on the EJ axis, then the hit probabilities for the
right and left parts are the same. If the POI location is not on the EJ axis, then the hit probabilities for the right and left parts of the G target are different. The side of the target where the POI is located has a higher hit probability. In view of the fact that the bullet gyroscopic stabilization causes a shift of the POI to the right, the right-hand side of the G target has a larger hit probability than the left one.
The coordinates of the points for the H target contour line are collected in Table 2.
Table 2 - Coordinates of the points for the H target contour line Таблица 2 - Координаты точек контура мишени типа H Табела 2 - Координате тачака контуре мете типа H
Points Z, m Y, m Boundary points of line Analytic function
A -0.225 0.0 A-B z(y)=-0.225
B -0.225 0.082 B-C y(z)= 0.384*z + 0.1684
C -0.1 0.13 C-D z(y)= -0.1
D -0.1 0.33 D-E y(z)=0.33
E 0.0 0.33 E-F y(z)=0.33
F 0.1 0.33 F-G z(y)=0.1
G 0.1 0.13 G-H y(z)= -0.384*z + 0.168
H 0.225 0.082 H-I z(y)= 0.225
I 0.225 0.0 I-J y(z)= 0
J 0.0 0.0 J-A y(z)= 0
Notes: The coordinates of the points are given according to Fig. 3. The symmetry of the points is shown with the orange and blue background.
The analysis of Table 2 allows us to describe the contour needed for integrating the double integral (4) and determining the hit probability. The contour of the H shooting target (the UCL (BCDEFGh) and the BCL (AJI)) can be described by a piecewise function:
CLH (z ) =
( 0.384z + 0.1684
- 0.225 < z <-0.Л
UCLr(z)= 0.33 -0.1 < z < 0.1
- 0.384z + 0.1684 0.1 < z < 0.225
ч
BCLr (z)=(0 - 0.225 < z < 0.225)
The coordinates of the points for the F target contour line are presented in Table 3.
Table 3 - Coordinates of the points for the F target contour line Таблица 3 - Координаты точек контура мишени типа F Табела 3 - Координате тачака контуре мете типа F
Points Z, m Y, m Boundary points of line Analytic function
A -0.225 0.0 A-B z(y)=-0.225
B -0.225 0.752 B-C y(z)= 0.384*z + 0.838
C -0.1 0.8 C-D z(y)=-0.1
D -0.1 1.0 D-E y(z)=1.0
E 0.0 1.0 E-F y(z)=1.0
F 0.1 1.0 F-G z(y)=0.1
G 0.1 0.8 G-H y(z)= -0.384*z + 0.838
H 0.225 0.752 H-I z(y)= 0.225
I 0.225 0.0 I-J y(z)= 0
J 0.0 0.0 J-A y(z)= 0
Notes: The coordinates of the points are given according to Fig. 3. The symmetry of the points is shown with the orange and blue background.
The contour of the F shooting target (the UCL (BCDEFGH) and the BCL line (AJI)) can be described by a piecewise function:
CLf (z ) = <
UCLf (z ) =
( 0.384z + 0.838 -0.225< z <-0.1^ 1.0 - 0.1 < z < 0.1
- 0.384z + 0.838 0.1 < z < 0.225
ч
БСЦ(z) = (0 - 0.225 < z < 0.225)
The coordinates of the points for the E target contour line are listed in Table 4.
Table 4 - E shooting target: coordinates of the points Таблица 4 - Координаты точек контура мишени типа E Табела 4 - Координате тачака контуре мете типа E
Points Z, m Y, m Boundary points of line Analytic function
A -0.15 0.0 A-B z(y)=-0.15
B -0.15 0.4973 B-C y(z)= -5.817*z - 0.375
C -0.17625 0.65 C-D y(z)=0.65
D -0.225 0.65 D-E z(y)=-0.225
E -0.225 1.402 E-F y(z)= 0.384*z + 1.488
F -0.1 1.45 F-G z(y)=-0.1
G -0.1 1.65 G-H y(z)=1.65
Points Z, m Y, m Boundary points of line Analytic function
H 0.0 1.65 H-I y(z)=1.65
I 0.1 1.65 I-J z(y)= 0.1
J 0.1 1.45 J-K y(z)= -0.384*z + 1.488
K 0.225 1.402 K-L z(y)=0.225
L 0.225 0.65 L-M y(z)=0.65
M 0.17625 0.65 M-N y(z)=5.817*z - 0.375
N 0.15 0.4973 N-O z(y)=0.15
O 0.15 0.0 O-P y(z)=0
P 0.0 0.0 P-A y(z)=0
Notes: The coordinates of the points are given according to Fig. 3. The symmetry of the points is shown with the orange and blue background.
Figure 5 - Geometric analysis of the E target (turn 90° counterclockwise) Рис. 5 - Геометрический анализ фигуры E (повёрнута на угол 90° против часовой
стрелки)
Слика 5 - Геометри]ска анализа мете Е (окренуте супротно од правца казачке
на сату за 90°)
The contour of the E shooting target (the UCL (EFGHIJK) and the BCL (DCBAPONML) (Fig. 5)) can be described by a piecewise function:
CLe (z ) =
UCLe (z) =
BCLe (z) =
( 0.384z +1.488 1.65
- 0.384z +1.488
0.65
- 5.817z - 0.375
0
5.817z - 0.375 0.65
- 0.225 < z <-0.11
- 0.1 < z < 0.1 0.1 < z < 0.225
- 0.225 < z <-0.1761
- 0.176<z <-0.15
- 0.15 < z < 0.15 0.15 < z < 0.176 0.176< z < 0.225
The K target has the smallest area. The silhouette and the geometry of the target are shown in Fig. 6.
в
G
В
□ J
с к
N
M
Figure 6 - Geometric constituent parts of the K target (see Fig. 3) Фиг. 6 - Геометрические составляющие мишени K (см. рис. 3) Слика 6 - Геометр^ски елементи мете К (видети слику 3)
The hit probability of this target can be found as an algebraic subtraction of the hit probabilities of its three domains
PhitK= Phit(ABEFGHILM) - Phit(BEDC) - Phit(KJIL) .
The coordinates of the vertices of the K target as a geometric figure and the analytical functions of its sides are collected in Table 5.
Table 5 - K shooting target: coordinates of the points Таблица 5 - Координаты точек контура мишени типа K Табела 5 - Координате тачака контуре мете типа K
Points Z, m Y, m Boundary points of line Analytic function
A -0.11 0.0 A-B z(y)=-0.11
B -0.11 0.0652 B-C y(z)= 0.795*z + 0.153
C -0.0714 0.0959 C-D z(y)=-0.0714
D -0.0714 0.1458 D-E y(z)= -0.795*z + 0.089
E -0.11 0.1765 E-F y(z)= 2.883*z + 0.494
F -0.0714 0.2878 F-G y(z)= S1(z)
G 0.0 0.33 G-H y(z)= S2(z)
H 0.0714 0.2878 H-I y(z)= -2.883*z + 0.494
I 0.11 0.1765 I-J y(z)= 0.795*z + 0.089
J 0.0714 0.1458 J-K z(y)=0.0714
K 0.0714 0.0959 K-L y(z)= -0.795*z + 0.153
L 0.11 0.0652 L-M z(y)=0.11
M 0.11 0.0 M-N y(z)=0
N 0.0 0.0 N-A y(z)=0
S1(z) = -57.9679*zA3 - 12.4167*zA2 + 0.33 z [-0.0714; 0 ] S2(z) = 57.9679*zA3 - 12.4167*zA2 + 0.33 z [ 0; 0.0714]
The piecewise functions CL for the K target are determined below.
CLu (z ) =
(2.883z + 0.494
UCLaU (z ) =
-0.11< z <-0.0714
A
- 57.9679z3 -12.4167z2 + 0.33 -0.0714< z < 0
57.9679z3 -12.4167z2 + 0.33 - 2.883z + 0.494
0 < z < 0.0714 0.0714< z < 0.11
clbedc (z
BCLall(z) = (0 - 0.225 < z < 0.225)
UCLbedc(z) = (-0.795- z + 0.089 -0.11 < z < -0.0714)
(z ) = ■
BCLbedc(z) = (0.795- z + 0.153 - 0.11 < z < -0.0714)
clkjil (z
(z ) = ■
UCLkjil (z) = (0.795- z + 0.089 0.0714< z < 0.11) BCL(z) = (-0.795-z + 0.153 0.0714< z < 0.11)
The sequence steps - the development of the calculation program in the CAS Mathcad will be described in the next section.
Calculation in the CAS Mathcad
The projectile horizontal range (in meters) is defined in the Mathcad as a range of variable x with the initial value of 50, the end value of 400 and the increment of 51 (the step size is 1 meter)
x:= 50, 51 .. 400 .
The text below shows the usage of the variable range to define the values of the vectors ay and az (eqs. 7 and 8) which characterize the bullet dispersion in the vertical and lateral directions
ayx:= 5.935 10-4x , azx:= 2.967 10-4 x ,
and the values of the vectors my and mz (eqs. 5 and 6) (which assign the POI location according to the MBT with the zero range at 300 meters)
myx:=a -0.029 +1.969 10-3x -2.645 10-6x2 -1.192 10-8x3
a:=0 ;
mzx:= -1.667 10-5 x + 1.667 10-7x2
Using the constant a, the my vector can be displaced vertically in the up-down direction. If a>0, then the zero range increases (>300 m). Usually, a is equal to 0 (zero range = 300 m). Thus, the parameter a makes it possible to simulate the position of the POA along the height of the shooting target.
For each value of x, set the xth element of the vector variables oy and oz, and my and mz with a corresponding value of x. The tables of numbers for the vector variables oy, oz and my, mz are of the same size. Each of them consists of 351 numbers. Using the variables sYi and sZi, we determine the center of the shooting target. For the K target, its center will be defined as (sYK, sZK). For the H target, as (sYH, sZH), and so on. The data on the position of the centers of the K, H, G, F, and E targets are shown in Table 6 (in meters).
Table 6 - Coordinates of the shooting target centers Таблица 6 - Координаты центра стрелковых мишеней Табела 6 - Координате центара мета за гамаше
Number Type (kind) of targets Center height coordinates, m Center coordinates in width, m
1 K sYK= 0.165 sZK= 0
2 H sYH= 0.165 sZH= 0
3 G syG= 0.275 sZG= 0
4 F syF= 0.5 sZF= 0
5 E syE= 0.825 sZE= 0
For the K target, the following expression can be proposed:
PhitK = Phit(ABEFGHILM) - Phi(BEDC) - Phit(KJIL) = PhitA - PhitL - PhitR, here Phit(ABEFGHILM), PhitA - probability of hitting the ABEFGHILM domain as a function of the PHR; Phit(BEDC), PhitAL - probability of hitting the BEDC domain as a function of the PHR; Phit(KJIL), PhitR -probability of hitting the KJIL domain as a function of the PHR.
Due to the conditional operator if, the UCL and the BCL of the constituent elements for the K shooting target are defined as:
UCLM (z ) =
f2.883z + 0.494 f(- 0.11 < z < -0.0714)
-57.9679z3 -12.4167z2 + 0.33 if(-0.0714< z < 0)
57.9679z3 - 12.4167z2 + 0.33 if (0 < z < 0.0714)
- 2.883z + 0.494 f (0.0714< z < 0.11)
Л
BCLdl (z) = (0 if (- 0.225 < z < 0.225)),
UCLbedc(z)=(-0.795-z + 0.089 if(-0.11 < z <-0.0714)), BCLbedc(z)= (0.795-z + 0.153 if(-0.11 < z <-0.0714)),
UCL^(z)= (0.795-z + 0.089 if(0.0714< z < 0.11)), BCL^(z) = (- 0.795-z + 0.153 if (0.0714< z < 0.11)).
The assessment of the hit probabilities of the elements of the K target is realized by the following formulas
0.225 UCLai (z)
Phit_ Ax J J
-0.225 BCL,,, (z)
1
2nazxoy}
- exp
[z - mzx ]2 [y - rnyx ]
+ -
2(azx )2 2(ayx )
dydz
(-0.0714) (-0.795-z+0.089)
Phit _ Lx := J J
(-0. 11) (0. 795-z+0 . 153)
1
2n<azxayx
- exp
[z - mzx ] , [y - myx 1
2(azx)2 2(avx)
dydz
(0.0714) (-0.795-z+0.089)
Phit _ Rx := J J
(0.11) (0.795- z+0.153)
1
2nazxOyx
- exp
z - mzx]2 + [y - myx]2 2(azx )2 2(ayx )2
dydz
For the UCL of the H, G and F target silhouettes, there are the
expressions
UCLH (z) :=
UCLa (z) :=
(0.384z + 0.1684) if (-0.225 < z < -0.1) (0.33) if (-0.1 < z < 0.1) (-0.384z + 0.1684) if (0.1 < z < 0.225) (0.384z + 0.388) if (-0.225< z <-0.1) (0.55) if (-0.1 < z < 0.1) (-0.384z + 0.388) if (0.1 < z < 0.225)
2
UCLF (z) :=
(0.384z + 0.838) if(-0.225< z <-0.1)
(1.0) if (-0.1 < z < 0.1)
(-0.384z + 0.838) if (0.1 < z < 0.225)
Finally, the UCL (UCLE(z)) and the BCL (BCLE(z)) of the E target silhouettes are designed by the formula:
UCLe (z) :=
(0.384z +1.488) if (-0.225< z <-0.1)
(1.65) if(-0.1 < z < 0.1)
(-0.384z +1.488) if (0.1 < z < 0.225)
BCLe (z) :=
(0.65) if (-0.225 < z <-0.176) (-5.817z - 0.375) if (-0.176< z < -0.15) (0) if (-0.15 < z < 0.15) (5.817z - 0.375) if (0.15 < z < 0.176) (0.65) if (0.176< z < 0.225)
The hit probability assessment for the H target is
0.225 UCLH (z)
Phu _ Hx := J J
-0.225 0
1
2nazxoyx
■• exp
+
[y - my, ]
2(ayx )2
[z - mzx ]2
2(ozx )2 dydz .
+
The formulas for the G and F targets are given in the appendix of this paper. The hit probability assessment for the E target is
0.225 UCLE (z)
Put _ Ex := J J
-0.225 BCLe (z)
1
2^ozxoy}
• exp
z - mzx]2 +
2K )2
CO
О >
О CM
of
Ш
a:
ZD О О
_J
<
С)
2: x о ш
I-
^
а: С
<л С
о >о
X ш н
О
О >
+
1у - тУх ] 2(<У )2
dydz
Results of the calculation
The curves characterizing the hit probabilities functions of five different types of targets (K, H, G, F, and E) that were determined using the MBT (5) are shown in Fig. 7. The space under the curve in the interval 50-400 m is highlighted in yellow. For all five targets, the POA (red square) was the same location - the middle of the BCL. The shooting conditions were considered to be ideal.
a)
b)
Figure 7 - Hit probability vs the PHR for the K and H targets (a) and the G, F, and E targets (b) for the range of 50-400 m Рис. 7 - Функции вероятности попадания в мишени K, H и мишени G, F, E для дистанций 50-400 м Слика 7 - Функци]е вероватноПе погодака мета K, H и мета G, F, E на pacmoja^y
од 50 до 400 м
The hit probability change vs the PHR for all five targets is a classical decreasing curve (Shereshevskiy et al, 1979). The dashed lines in Fig. 7a show the ranges that correspond to the hit probabilities of 0.9, 0.8, and 0.7 for the K target.
Changing the shape for the hit probability vs the PHR with respect to the F target and the location of the POA (red square) is shown in Fig. 8.
a)
b)
c)
Figure 8 - Hit probability vs the location of the POA for the F target for the range of 50400 m
Рис. 8 - Вероятности попадания в мишени F для дистанций 50-400 м и положения
точки прицеливания Слика 8 - ВероватноПе погодака мете F на растоjашу од 50 до 400 м и положаjа
тачке нишашеша
Fig. 8 shows that if the range to the shooting target is 200 meters and if the POA is positioned on the BCL (see the red square as the graphical symbol), then the hit probability is 0.92 (Fig 8a). If the POA is 0.5 m higher, then the hit probability for the PHR=200 m is increased to be 0.98 (Fig. 8b). If the POA is increased to 1.0 meters, then the hit probability decreases to 0.08. Thus, if the distance to the shooting (real) target is known as well as the target geometric dimensions and the ballistic characteristics of the bullet flight path (MBT and bullet dispersion), the POA determining algorithm for a remote weapon station of a ground combat vehicle can be proposed.
Conclusions and Implications
This paper represents a method for assessing the single shot hit probability for a stationary (fixed) shooting target as a function of the PHR, the selected target type and the location of the POA. The proposed procedure has three blocks. The proposed method is built on a modular basis and allows a user to change types of weapons and types of shooting targets.
The single shot hit probability depends on five variables: height of the POI; lateral displacement of the POI relative to the target; bullet dispersion in the height (vertical direction); lateral bullet dispersion (horizontal direction); and the forms of the shooting target contour. To order these variables, three auxiliary quantities can be introduced. These are the projectile horizontal range (PHR), the altitude of the point of aim
and the lateral displacement of the point of aim. The last two variables can be combined into the «POA location».
For demo calculations, the external ballistics and the bullet dispersion characteristics of the 5.45*39 Kalashnikov assault rifle model MPi AK-74N were selected. Five types of Swiss military targets imitating OPFOR combatant silhouettes were used as shooting targets. For illustrating the operability of the proposed method, the hit probability for the K, H, G, F, and E shooting targets was evaluated as well as for the PHR from 50 to 400 m. The calculations have shown that the PHR and the form of the target cause a complex dependence between the PMt and
the shooter-target distance.
The developed mathematical model can be used in tactical computer interactive simulators that emulate combat situations of a tactical strike. Likewise, the proposed mathematical models of shooting targets can be used for improving digital target range simulators or in fire-control systems of remote weapon stations.
Appendix
Formulas for the hit probability assessment of some shooting targets
The hit probability vs the PHR for the G and F targets can be expressed as
0.225 UCLG (z)
Put _G := J J
-0.225 0
1
2Kozxoyx
■ exp
+
[y - myx ] 2{oyx)2
dydz ;
[z - mzx ]2
2K )2
+
0.225 UCLf (z)
Phu _ F := J J
-0.225 0
1
2nazxayx
■ exp
+
[y - myx ] 2(oyx )2
dydz .
[z - mzx ]2
2K )2
+
2
2
References
http://weaponland.ru/load/avtomat_shturmovaja_vintovka_serii_mpi_k_mpi _ak/10-1-0-31. (in Russian). Accessed: 20.06.2018.
http://www.wikiwand.com/en/Shooting_target. Accessed: 20.06.2018.
Khaikov, V.L. 2018. Single shot hit probability estimation as a result of a numerical solution of double integrals using Mathcad. Vojnotehnicki glasnik/Military Technical Courier, 66(4), pp.739-756. Available at: https://doi.org/10.5937/vojtehg66-17433.
-NVA der DDR. 5.45-mm-Maschinenpistole AK74 und leichtes Maschinengewehr RPK74. Beschreibung und Nutzung1985 (in German).
Peelen, J. 2017. There is more to short barrels than just velocity. [Internet]. Available at: www.quarryhs.co.uk/shortbarrel2.pdf. Accessed: 20.06.2018.
Rodney, S. 2012. RPG encounter modeling. SURVIAC Bulletin, 27(1), pp.1-5.
Svateev, V.A. 2013. Avtomatchik dolzhen i mozhet porazhat' golovnuyu figuru. Bulletin of the Academy of Military Sciences, 43(2), pp.127-131 (in Russian). (In the original: Сватеев В.А. 2013. Автоматчик должен и может поражать головную фигуру, Вестник академии военных наук, №.2(43). С.127-131).
Svateev, V.A. 2014. Tochniy sposob raschota veroyatnosti popadaniya v figurnuyu tsel'. Bulletin of the Academy of Military Sciences, 49(4), pp.61-66 (in Russian). (In the original: Сватеев В.А. 2014. Точный способ расчёта вероятности попадания в фигурную цель, Вестник Академии военных наук. 2014. T.49, № 4. C.61-66.)
-TRADOC. 1975. Range and lethality of US and Soviet anti-armor weapons.
-TRADOC. 2015. Worldwide Equipment Guide, Ground Systems. Vol. 1.
Venttsel' Ye.S. 2006. Teoriya veroyatnostey. Moscow: Vysshaya shkola (in Russian). (In the original: Вентцель Е.С. 2006. Теория вероятностей. Москва: Высшая школа).
Wollschlager, D. 2017. Analyzing shape, accuracy, and precision of shooting results with shotGroups. [Internet] Available at: https://cran.r-project.org/web/packages/shotGroups/vignettes/shotGroups.pdf. Accessed: 20.06.2018.
МЕТОДИКА ОЦЕНКИ ВЕРОЯТНОСТИ ПАПАДАНИЯ ОДИНОЧНЫМ ВЫСТРЕЛОМ КАК ФУНКЦИЯ ГОРИЗОНТАЛЬНОЙ ДАЛЬНОСТИ С УЧЁТОМ ВИДА СТРЕЛКОВОЙ МИШЕНИ И ПОЛОЖЕНИЯ ТОЧКИ ПРИЦЕЛИВАНИЯ
Вадим Л. Хайков,
независимый исследователь, г. Краснодар, Российская Федерация
ОБЛАСТЬ: прикладная математика, исследование операций ВИД СТАТЬИ: оригинальная научная статья ЯЗЫК СТАТЬИ: английский
Резюме:
Оценка вероятности попадания является важной частью процедуры анализа эффективности стрельбы. В статье предложен метод оценки вероятности попадания в статическую стрелковую мишень как функция горизонтальной дальности до мишении (ГДМ), её типа и положения точки прицеливания. Предлагаемая процедура оценки состоит из трех частей. Часть I посвящена двойной интегральной формуле с учётом изменений в местоположении средней точки попадания (СТП) на поверхности мишени в зависимости от ГДМ. Характеристики дисперсии пуль в зависимости от ГДМ определены в части II. Основой части I и части II является двумерное некоррелированное распределение вероятностей Лапласа-Гаусса. Как правило, функции ординат (абсцисс) СТП от аргумента ГДМ и функции ординат (абсцисс) рассеивания пули как функции ГДМ представлены в виде классических полиномов. В части III даётся описание контура мишени. В математическом плане контур цели является кусочной функцией, которая описывает верхний и нижний края стрелковой мишени, а во-вторых, контур цели в формуле двойного интеграла выступает как предел интегральной формулы части I. Все вычисления реализованы в PTC Mathcad. Предлагаемый метод построен на модульной основе и позволяет пользователю изменять типы оружия и виды мишеней. Для демонстрационных расчётов использованы данные характеризующие внешнюю баллистику и характеристики рассеивания 5.45*39 автомата Калашникова. В качестве рассматриваемой модели автомата была выбрана MPi AK-74N. В качестве стрелковых мишеней использовались пять видов швейцарских армейских мишеней, которые имитируют различные части силуэта комбатанта. Для иллюстрации работоспособности предложенного метода оценена вероятность попадания в мишени K, H, G, F, E и для диапазона ГДМ от 50 до 400 м.
Ключевые слова: вероятность попадания, численное решение, стрелковая мишень, двойной интеграл, рассеивание пуль, горизонтальная дальность, Mathcad.
ОДРЕЪИВА^Е ВЕРОВАТНОЪЕ ПОГОТКА JЕДНИМ ХИЦЕМ У ЗАВИСНОСТИ ОД ХОРИЗОНТАЛНЕ РАЗДА^ИНЕ, РАЗЛИЧИТИХ ВРСТА МЕТА И ПОЛОЖКА ТАЧКЕ НИША^Е^А
Вадим Л. XajKOB
независни истраживач, Краснодар, Руска Федерац^а
ОБЛАСТ: применена математика, операциона истражива^а ВРСТА ЧЛАНКА: оригинални научни чланак иЕЗИК ЧЛАНКА: енглески
Сажетак:
У чланку je представлена метода за одре^иваше вероватноЯе поготка стационарних мета у зависности од хоризонталног растоjаша до мете (ХРМ), врсте мете и локац^е тачке цила. Предложени поступак састоjи се од три блока. Први блок се бави формулом двоструког интеграла ща узима у обзир промене у локации тачке поготка у зависности од ХРМ. Карактеристике дисперз^е метка, у зависности од ХРМ, одре^у се у другом блоку. Бивар^антна некорелисана Гаус- Лапласова расподела вероватноЯе узета je за основ прва два блока. Функц^е локац^е тачке цила (ордината/апциса) у зависности од ХРМ и функц^е карактеристика дисперз^е погодака у зависности од ХРМ приказане су у облику полинома. У треЯем блоку наведен je опис контура силуета мета. У математичком смислу, контура мете je хибридна функц^а ща дефинише горше и доше ивице мете и, тако^е, представка границу формуле двоструког интеграла првог блока. Сви прорачуни су се вршили у PTC Mathcad. Предложени метод заснован je на модуларном принципу и омогуЯава кориснику да меша врсте оруж1а и мета. За илустрац^у je изабран приказ израчунаваша карактеристика дисперз^е поготка ко\и je постигнут из аутоматске пушке „калашшиков", калибра 5.45*39, модел MPi AK-74N. За мете je изабрано пет типова шваjцарских воних мета ще имитираjу различите делове силуете воjника. Операбилност предложеног метода илустрована je одре^ивашем вероватноЯа пога^аша мета К, Н, G, F и Е на удалености од 50 до 400 м. Клучне речи: вероватноЯа поготка, нумеричко решаваше, мета за гамаше, двоструки интеграли, дисперз^а погодака, хоризонтално растоjашe, Mathcad.
Paper received on / Дата получения работы / Датум приема чланка: 06.08.2018. Manuscript corrections submitted on / Дата получения исправленной версии работы / Датум достав^а^а исправки рукописа: 26.11.2018.
Paper accepted for publishing on / Дата окончательного согласования работы / Датум коначног прихвата^а чланка за об]ав^ива^е: 28.11.2018. © 2019 The Author. Published by Vojnotehnicki glasnik / Military Technical Courier (www.vtg.mod.gov.rs, втг.мо.упр.срб). This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/rs/).
© 2019 Автор. Опубликовано в «Военно-технический вестник / Vojnotehnicki glasnik / Military Technical Courier» (www.vtg.mod.gov.rs, втг.мо.упр.срб). Данная статья в открытом доступе и распространяется в соответствии с лицензией «Creative Commons» (http://creativecommons.org/licenses/by/3.0/rs/).
© 2019 Аутор. Обjавио Воjнотехнички гласник / Vojnotehnicki glasnik / Military Technical Courier (www.vtg.mod.gov.rs, втг.мо.упр.срб). Ово jе чланак отвореног приступа и дистрибуира се у складу са Creative Commons лиценцом (http://creativecommons.org/licenses/by/3.0/rs/).
