CC BY
31
ESTIMATE OF PROJECTILE INITIAL VELOCITY AS A SOLUTION OF A TWO-POINT BOUNDARY VALUE PROBLEM
FIELD: Mechanics - Ballistics ARTICLE TYPE: Original Scientific Paper ARTICLE LANGUAGE: English
Abbreviations
ADR - Artillery Doppler Radar IPV - instantaneous projectile velocity
BC - ballistic coefficient ODE - ordinary differential equation
BVP - boundary value problem PIV - projectile initial velocity
IC - initial condition TOF - time of flight.
C\l CD !± cp
E"
JD .a o
<u
Vadim L. Khaikov 5
independent researcher, Krasnodar, Russian Federation, ^
e-mail: [email protected],
d
ORCID iD: http://orcid.org/0000-0003-1433-3562 |
http://dx.doi.org/10.5937/vojtehg66-15097
b
o p
a o
Summary:
The estimating of a projectile initial velocity is formulated as a two-point boundary value problem. To solve it, the data of a Doppler Radar or the results of solving the Cauchy problem can be used. The projectile initial velocity v0 estimation process is based on the numerical solution of a o system of ordinary differential equations and the bisection method. The iterative calculating process is interrupted when a predetermined accuracy of a projectile initial velocity and a predetermined value of the width of velocity's search interval is reached. In the article, the block diagram of the algorithm for the projectile initial velocity process is developed. The Mathcad program code for mathematical modeling and a computer 2 simulation of the projectile initial velocity estimation process for a 57mm
o
armor-piercing projectile of ZIS-2 anti-tank gun 1943 model is given. a>
a
Key words: applied external ballistics ofbarreled weapon, projectile E initial velocity estimation process, mathematical modeling, computer simulation, Mathcad software.
v o ik ai
h K
Introduction
g A large number of practical problems of external ballistics of barreled
artillery systems require information about a projectile initial velocity (PIV), briefly denoted as v0. This type of assessment is based on experimental
data containing information about a projectile velocity1 at trajectory points oc with known time or at trajectory points with known horizontal ranges.
LU
dc Typically, methods and devices that implement v0 estimates use empirical
ZD
CO mathematical dependencies based on assumptions about the nature of a projectile deceleration process (Alon, 1989), (Zaytsev et al, 2004) or o general information about a law of frontal air resistance (Semenov, 2013).
X
o
The aim of this paper is to substantiate the PIV estimating process for EC tubed artillery based on the numerical solution of a two-point BVP. The £ PIV2 estimating is based on:
^ - multiple repeated solution of the system of ordinary differential
equations (ODEs) of the longitudinal movement of an artillery projectile using a method of numerical integration; - comparison of the received result with the right boundary condition. Unlike previously known solutions, the technique assumes the use of !5 the actual law of air drag, ballistic coefficient (BC), atmospheric parameters ^ and actual angle of departure. If necessary, this data can be corrected o before the start of computational or physical experiments.
X LU I—
o z
o >
An instantaneous projectile velocity (IPV) as a function of time
It is known that the longitudinal motion of an artillery projectile (stabilized by rotation) in the Earth's atmosphere can be described by the system of ODEs with an independent argument time (t) (Burlov et al, 2006, p.249). Systems of ODEs (1) are obtained by assuming the absence of the influence of lifting and lateral forces. There is no wind effect.
1 The article considers only the longitudinal motion of the artillery projectile. Projectile movements such as rotational, nutational and precessional movements are not taken into account.
2 Under PIV or v0 is meant a fictitious velocity for TOF t=0. Muzzle velocity is the actual velocity of the projectile that left a gun barrel. Information about PIV is used to calculate tabular (graphical) firing tables or for developing firing control equipment.
(1)
— = v cos(û) dt
dy = vsin(^) dt
d9 g (a)
dV = -(( + g-sin(»))
dt
here
x - the abscissa of the trajectory, m; y - the ordinate of the trajectory, m;
6 - angle of velocity vector relative to the horizon plane, degree or
radian; v - IPV, m/s; t - TOF, s;
g - acceleration of gravity, m/s2; J - acceleration of the air resistance force, m/s2:
J = CH (y )F (v ), (2)
where
C-BCis estimated by the expression of Siacci (Germershausen,
1982, p.159), m2/kg; H(y) - a function characterizing the dependence of air density vs altitude;
F(v) - a function characterizing the dependence of the air resistance vs an IPV.
The ballistic coefficient is calculated as
d2
C = 1000 • i • —
q
where
i - value of the coefficient taking into account the shape of the projectile (dimensionless); from a formal point of view, this is a fitting coefficient;
d - caliber of projectile (main diameter), m; q - projectile weight, kg.
ф (Л
CD CD
"5
>
О
см
ОС Ш ОС ZD О
О _|
< с;
X О ш
I—
>-
ОС <
(Л <
-J
О >о
X ш I—
о
о >
On the basis of the foregoing, the Cauchy problem, namely the problem with initial conditions (IC), is formulated. The totality of its elements is shown on Figure 1.
Figure 1 - The main task of external ballistics as the Cauchy problem Рис. 1 - Основная задача внешней баллистики в виде задачи Коши Слика 1 - Основни задатак спо^ше балистике у облику Коши]евог проблема
The acronym ADR means an Artillery Doppler Radar, which is intended for a projectile velocity measurement.
Since the system of ODEs (Figure 1, box 6) does not have "stiff" properties, it is therefore recommended to implement a numerical solution by dint of the fourth-order Runge-Kutta method or Bulirsch-Stoer method. As a result of the solution, there are five trajectory parameters: TOF (t), horizontal range (x), height of trajectory (y), angle of velocity vector (see (1)) (в), IPV (v).
From a mathematical point of view, the difference between PIV and IPV is as follows: PIV is a constant; IVP is the interval from the PIV value to the so-called collision velocity (that is the velocity of the projectile, which has an impact on a target).
As an example, we can estimate the external ballistics of a 57mm 53-BR-271SP Armor-Piercing Tracer Solid projectile with a weight of 3.14 kg,
fired from an anti-tank gun ZIS-2 M1943, whose gun barrel had a length of 73 caliber (Burlov et al, 2006, p.425). This artillery cannon was widely used in the World War II on the Soviet-German front. The serial production of the ZIS-2 was carried out in the period from 1941 to 1949.
The left graphical window (Figure 2) shows the dynamics of the deceleration process as a change of the IPV value vs TOF3.
Figure 2 - IPV vs TOF graphs and trajectories of a gun projectile Рис. 2 - Графики изменения скорости и траектории снаряда Слика 2 - Графици тренутне брзине проjектила (ТБП) и времена лета и путане артиъерщског проjектила
The angle of departure is assumed equal to19 angular min. « 4.945 10-3 radians, which corresponds to the flat trajectory of a projectile flight. The PIV is 990 m/s. The three curves (Figure 2) have different BCs C=0.4, 0.5, 0.6 m2/kg. For the computer simulation of the projectile flight path, Siacci air drag law was used. The effect of the atmosphere influence was imitated by Eberling formula:
46
1000000
(y)
H (y )= 10"
The time interval for integration is taken as 0-1.1 s.
3 This kind of IPV vs TOF dependence is characteristic for the flat trajectory weapon (for example, tank or anti-tank artillery). The summit of flight path for a distance ~900 m is a little more than 1 m.
The values of the variables t, x, y, 6, v were determined every 1.1 milliseconds of the projectile's flight, which corresponds to 1000 points inside the time interval. For each curve, the value of BC was constant for the whole domain of functions V=f(t). > The curves of Figure 2 are monotonically decreasing functions of
TOF. In (Semenov, 2013) they are called "trajectories" of velocities. In the time interval (0-1.1 s), the graphs in form are close to linear functions. The p( dashed line t = 0.5 s (time after leaving the gun muzzle), crossing the
LU
E velocity-time curves, forms three points of intersection. Their position is o determined by the time interval t and the value of the BC. The higher the ° ballistic coefficient, the more significant the reduction of IPV in the chosen time interval by the constant angle of departure. Thus, an advanced x knowledge of the BC value and time delay will allow to restore v0.
uj In the Cartesian coordinates Y-X in (the right graphic window of
>_ Figure 2), the shape of the trajectories is shown. Numerals 0.6, 0.5, 0.4 < indicate a discrete decreasing of the BC value.
The growth of the BC causes a decrease in the horizontal range and a decrease in the value of the summit of a projectile flight path. The displacement of these parameters causes significant transformations of SN the descending branch of the trajectory in comparison with the ascending branch.
In order to simplify the calculations and reduce the computation time,
>o Siacci air drag law was used in the form of the F-type curve
X LU
(Germershausen, 1982, p.160)
g Fsj (v ) = 0.2002v - 48.05^ (0.1648v - 47.95)2 + 9.6 + 00442v(v ^ § 371 + ( v 1 (3)
200
with the velocity range 0-1200 m/s replaced by a quadratic dependence, which has the form:
F, (v) = B<v2, (4)
where the value of the coefficient Bsi is 0.0002696 (dimensionless).
The maximum degree of the power by variable v by Siacci air drag law corresponds to 10, but formula (4) allows its reduction to a power of 2. The result of this approximation is shown in Figure 3.
400
300
200
100
F(V)
/ Fsiap
-
- (1),/ -
- ,y {2) 1 1 -
v, mfc
500
i> Ю3
L 5x
103
C\l CO !±
s=
JD О
<u
J3
со >
CO T3
О .Q
О
CO
4—
О
Figure 3 - Approximation of Siacci air drag law using the curve (4) Рис. 3 - Аппроксимация закона Сиаччи кривой (4) Слика 3 - Апроксимаци^а Siacci-jевог закона отпора ваздуха помоПу криве (4)
о
<Л
(Л СО
It should be noted that the comparison of the trajectories in the Y-X coordinates calculated according to Siacci air drag law and using formula
(3) (only the law of the frontal air resistance was changed, the remaining parameters retained the previous values) leads to a decrease of a horizontal range about 10% regarding the result for Siacci. However, the replacement of (3) by (4) does not affect an accuracy of the solution of the BVP, since the conditions on the boundary are also calculated using
(4). The gain from replacement consists of mathematical expressions for compactness and saving of calculation time, which is important for repeatedly recurring calculations at the testing stage and for checking the algorithm correctness. Time measurement for the numerical solution for the system of ODEs (1) using the fourth-order Runge-Kutta method fulfilled on the PC Lenovo Ideapad Z510 and Mathcad 15 software gave the following results: 17 ms for Siacci air drag law contra 14 ms for the quadratic law. The gain in calculating time was 3 ms, that is 17%.
If the BC is constant and PIV is a variable value, then the family of curves will have the form shown in Figure 4.
о о
<D >
О <D
сГ
О <U CO
E +-<
(Л
Ш
>
о
'со .с
ф (Л
CD CD
"о
>
О
см
ОС Ш ОС
о
о _|
< с;
X О ш
н
>-
ОС <
(Л <
о >о
X ш н
о
о >
Figure 4 - Velocity vs. time graphs and trajectories of a projectile Рис. 4 - Графики изменения скорости и траектории снаряда Слика 4 - Графички прикази брзине и путане проjектила
The invariable ballistic coefficient (C = 0.5 m2/kg) and a discrete increasing of PIV (v0 = 980; 990; 1000 m/s) lead to the velocity change in
the velocity-time coordinates in the form of three monotonically decreasing curves (the left graphic window of Figure 4). As before, the three curves characterize the trajectories of a projectile flight by the angle of departure of 19 angular minutes. Just like before, the presented graphs are close to linear dependences. However, the difference between Figure 4 and Figure 2 is the parallelism of the graphs of functions relative to each other, which predetermines the absence of intersection points. That is why the approximation of the functions (Figure 4) by the polynomial of the first degree will give the same tangents of the line slope (the condition of parallelism).
The right graphic window of Figure 4 shows the changes in the trajectories of the projectile in the Cartesian Y-X coordinates. The PIV increases discretely (v0= 980, 990, 1000 m/s) but the value of BC and the angle of the projectile departure remain constant.
If a velocity value at a given time is determined through the ADR measurement or by means of the Cauchy problem solution and if the angle of a projectile departure and the coordinates of the departure point are known, then in this case one can formulate a two-point BVP.
The two-point boundary value problem
The essence of the two-point BVP is that, in contrast to the Cauchy problem, the boundary conditions are not strictly defined on the left (as ICs), but also on the right boundary. For example, the value of velocity condition (IPVt=k) restricts the result on the right boundary. Thus, the search for a solution to the systems of ODEs (1) that satisfies the boundary conditions for the time interval t=k (see Table 1) is organized as a procedure consisting of the following items:
1. Solution of the Cauchy problem described in the previous section;
2. Comparison of the result with a given condition on the right boundary;
3. Making a decision about changing the ICs in the ODEs with the subsequent return to point 1 or interrupting the calculations.
Table 1 - Boundary conditions on the right and left boundaries Таблица 1 - Краевые условия на правой и левой границах Табела 1 - Гранични услови на десноj и лево] граници
№/№ Left border: t=0 s. Right border: t=0.5 s.
1 x(0)=0 -
2 y(0)=0 -
3 phi(0)=k -
4 sought value (PIV) v(0.5)=a2
c\i
CD d Ci
S=
JD _Q О Ci Ф J3
ro
>
ro
TD
О .Q
О Ci
ГО
4—
О
О <Л
(Л
го
о
о ф
>
о ф
сГ
о ф
го Е
(Л
Ш
>
о
'го .с
The value of the projectile velocity at the right boundary of the interval is determined by a theoretical calculation or a physical experiment. In the first case, the Cauchy problem (1) is solved according to the beforehand known ICs. The estimation of v0 by an experimental determination is
carried out by using a ballistic chronograph or ADR.
The result of solving the boundary problem is obtaining the velocity of the projectile on the left border of the interval, herewith the fourth element of the matrix of ICs, v0 (see Figure 1, box 5) is determined.
It should be noted that sometimes the boundary problem is called "bringing the projectile velocity to the muzzle". This idiomatic expression is valid only from a mathematical point of view. In view of the fact that the physics of the projectile motion process in a gun barrel and in the air by the shot is very complex, the muzzle velocity (vm) and PIV (v0) are not equal
to each other due to the presence of the aftereffect. "PIV is a fictitious velocity that a projectile would have at the point of departure, so that in the absence of a period of aftereffect, its velocity outside the aftereffect line
segment would be equal to the actual velocity. Such a velocity is 1-2% greater than the projectile muzzle velocity" (Burlov et al, 2006, p.115). ço Thus, the fictitious value of v0 exists as a phenomenon of a simplified
CD
>
O CM
point t = 0.5 s, then the solution will be found in the iteration at which
V0 )" v0.5 = 0 . (5)
ol. mathematical description of a shot as a physical process, but plays an independent role in the preparation of tabular or graphical firing tables. Consequently, it is necessary to know it when the two-point boundary problem needs to be solved on the basis of experimental data.
yj If the condition on the right boundary is regulated by the IPV at the
oc
o
o <
In this expression, the decrement vc (t05, v0) correspond the IPV at a uj point with a time t=0.5 s by PIV= v0. The velocity v0 is the solution of the
> Cauchy problem. This is reminiscent of the subscript "c". The subtrahend
< is the user-defined right boundary condition. Expression (5) is a nonlinear equation of the argument v0. In fact, the right-hand side of (5) is an
admissible error of the difference denoted as sf. Under these conditions,
ot (5) is transformed to the form
<
^ t \
2 vc (t0. 5> v0 )-v0 . 5 =Sf . (6)
0
1 Since solution (6) is carried out in the interval of PIV v01 - v02, it is
o rational to impose an additional condition on its limiting width (after which ^ further reduction of the interval is not rational). This restriction makes the
> iteration interrupt mechanism more flexible. Nonlinear equation (6) is solved numerically by the method of
bisection (division of a line segment in half). The algorithm of action is shown in Figure 5.
Box "1" is intended for entering the necessary constants. The block diagram assumes that the right boundary condition for velocity is determined from outside, for example, using the ADR or another method. Box "2" sets the limits of the search interval for initial speeds. It is assumed here that condition (5) has different signs on the left and right selected boundaries, and the left-hand side of function (6) is continuous. The task of box "3" is to divide the line segment of the search in half.
C\l CO !±
E"
JD .a (3
<u
J3
ro
>
ro
T3
о
.Q
О
ГО
4—
О
О <Л
(Л
го
о
о ф
>
о ф
(Г
о ф
го Е
(Л
ш
>
о
'го
h
Figure 5 - Block diagram of the algorithm for solving the BVP Рис. 5 - Блок-схема алгоритма решения граничной задачи Слика 5 - Шема алгоритма за решаваше проблема граничне вредности
<u tn
CO CD
O
O CM
The solution of the Cauchy problem is carried out in box "4". A comparison of the calculated value and the data of the ADR is performed in box "5". Boxes "6" and "7" serve to interrupt the iterative process of estimating the value of the PIV. The calculations are terminated upon > reaching the predetermined accuracy sf of determining the initial speed
and the pre-set value of the width of the search interval sv. In boxes "8-
oc 11", the process of controlling the change in the boundaries of the interval oE of the search for the PIV when it is divided in half is carried out. The o number is 0.5 s present in the description of blocks "1.1" and "5" shows ^ that the solution is sought for the right boundary condition of 0.5 s (the < velocity is fixed at 0.5 s from the start of the flight).
The computer simulation of a BVP solution will be illustrated due to o the Mathcad engineering math software and the sbval function (Solodov & m Ochkov, 2004, p.101). Such commented program code is presented dc below.
^ At the first stage, it is necessary to determine the angle of departure
Q0 (zero in subscript is the symbolical index) in radians (u=19 angular
minutes), the ballistic coefficient and the acceleration of gravity as a S constant
^ u := 19 O0:= —■— C := 0.5 g := 9.81.
£ 0 180 60
Next, we will enter the starting value to find the required PIV (m/s), the
w starting and ending values of the time interval (s)
° z := 983 t1:= 0 t2 := 0.1.
The starting value z is entered to the program-code as a matrix that contains one element z0.
The matrix-column D is the right-hand parts of the system of ODEs (1). It includes variables with the following notation:
- the abscissa of the trajectory x - y0;
- ordinate of the trajectory y - y1;
- angle of inclination of the tangent Q - y2,
- IPV v - y3.
All four y0.. y3 elements described variables from the virtual matrix Y
D((, y ):=
y3 • cos(y2 )
y3 • sin(y2) - g
• cos
y3
(C • (l00000046(yi))• (0.00002694y3)+ gsin(y2))_
(y2 )
The fourth element of the matrix D stores information about the BC, about the function H(y) and about the law of an air drag. When using an air drag law other than quadratic, it is necessary to modify this element of the D-matrix.
The matrix-column L contains the ICs in the form of their known numerical values and unknown variables. In our case, these are d0 and z0
L((1, z ):=
r 0 ^
00
V z0 J
The parameter of the left side of expression t1 corresponds to the conditions on the left boundary of the interval. Here z0 is specified as an
element of the matrix z.
The matrix-column score stores the difference between the initial condition given at time-point t2 and the value of the sought solution at that point. The dimension of the matrix L is equal to the dimension of the matrix z. The parameter t2 denotes an employment of conditions from the right boundary
score((2, y) := y3 - 976.045 .
(7)
C\l cn
S±
cp
E"
JD _Q O
<U J3
ro
>
ro
T3
o
.Q
O Ci
CO
M—
O
O tn
ro
o
Oo >
o 'cT
o <u ro E
tn LU
The value of the velocity is 976.945 m/s and it is the speed of the projectile at the right boundary of the time interval 0-0.1 s. and corresponds to the solution of the Cauchy problem with the following parameters and ICs:
the parameters: C=0.5 m2/kg; law of air resistance - formula (3);
change in air density - Eberling's formula; time interval of integration 0-1.1 s. with 1000 calculation points; ICs: x0 = 0; y0 = 0; d0 = 19 min.; v0 = 990 m/s.
>
o
'ro .c
0
<u tn
CO CD
15 >
o
CM
QC UJ QC
O
o <
o 2:
X
o
LU
I—
>-
CC <
< -J
CD >Q
x
LU I—
O z
—3
O >
Calling the solver function sbval for solution of the two-point BVP is performed by the command:
S:= sbval(z, t1, t2, D, L, score).
The matrix-column of the solutions S has the same dimensions with the score and z matrices. As a result, we obtain
S = 989.987 m/s.
The error in the numerical solution of the BVP was 0.0013% (the actual value of the PIV was 990 m/s).
If instead of the instantaneous velocity v in the ODEs, the velocity projections on the OX, OY axis are used, then the system of ODEs (1) is transformed into expression (8)4 (Chujko, 1976, p.38, 42)
(Q)
dx
— = vcos dt
dy = v sin(Q) dt
dQ g (q)
— = — COS(Q) dt v
— = - J cos(Q)
dt
da dt
= -(J sin(Q) + g)
(8)
where u is the projection of IPV on the abscissa axis, m/s; w is the projection of IPV on the ordinate axis, m/s.
Formulating the instantaneous velocity v through its projections on the axis OX, OY, we modify the first three equations (8) and obtain (9):
4 the first - third equations of (8) describe the kinematic relationships, the fourth and fifth -aeroballistic
dx dt dy_ dt
■ = u
= a
_ = _ g
dt v du dt da dt
au g 2(a\ — = _ —cos (0)
= _ J cos(<9) = _(( sin(0) + g)
(9)
C\l CD
s=
JD .a (3
<1J J3
<0 >
ro
T3
o
.Q
O
Similar actions can be performed with respect to the quadratic law of air resistance (4):
Fsla (v ) = Bs
(
u
cos
0
(10)
As in the previous example, we find the values of the parameters u, a in 0.1 s after a 57mm projectile has left the gun barrel. The following result is a solution of the Cauchy problem for the system of ODEs (9): u = 976.938 m/s; a = 3.856 m/s. Since the horizontal and vertical projections of the velocity v in the Cartesian coordinate system are related to the velocity by the relation
v = vu
then by substituting the values we have
lu2 +a2
v = V976.9382 + 3.8562 = 976.946 m/s .
The obtained value of the IPV agrees with the value which was found in the previous example (the condition of expression (7)). The value of the relative error is 0.0001%.
Solving the BVP for the system of ODEs (9) at u = 976.938 m/s; a =3.856 m/s (right boundary conditions) we obtain the velocity projections on the OX, OY axis at time t = 0 (unknown IC) equal to u0 = 989.975 m/s;
is equal to:
= V989.9752 + 4.9032 = 989.987 m /s.
a0 = 4.903 m/s, from which it follows that the PIV v0
V
22 u + a0
ro
M—
O
O tn
ro
o
Oo
<D >
o <D
'¡IT
o
<1J
ro E
tn LU
>
o
<0 -C
2
ф (Л
CD CD
"о
>
О
см
ОС Ш ОС
о
о _|
< с;
X О ш
н
>-
ОС <
(Л <
о >о
X ш н
о
о >
The errors in the numerical solution of the BVP for the system of ODEs (1) and the system of ODEs (9) are the same and equal to 0.0013%, so the uniqueness of the solution of the two-point boundary problem takes place.
From the point of view of practice, it is of interest to determine the dependence of the error in estimating the PIV SPIV as a function of error in
determining the boundary condition SPIV_RB. The results of such an assessment are shown in Figure 6.
Figure 6 - Change in the PIV error of a projectile (SPIV) as a function of the assessment
of the IPV on the right boundary SIPV_RB Рис. 6 - Изменение погрешности определения начальной скорости снаряда (SPIV ) с учётом погрешности задания скорости на правой границе Sipv-rb Слика 6 - Промена грешке ТБП проjекmила (SPIV) у функции процене ТБП-а на
десноj граници SIPV RB
As an upshot, the graph of the dependence (Figure 6) on the error in determining the PIV (SPIV) is obtained taking into account the error of the
IPV specification on the right boundary (SIPV-RB), which can be approximated by the formula
tS (45° К
SPIV
y IPV-RB
IPV-RB
Such a dependence is symmetric about the axis OY passing through the origin. In the case of a velocity error at the right boundary, SIPV_RB = 2%, the corresponding error in estimating the PIV is also equal to SPIV = 2%.
Conclusion
The paper proposes to estimate the projectile initial velocity v0 of the
artillery projectile on the basis of the system of ODEs (1) and (8, 9) taking into account the known law of air resistance, ballistic coefficient, atmospheric parameters and angle of departure. In the process of implementing this approach, the problem of estimating the PIV in the form of a two-point boundary problem is formulated. An algorithm for estimating the PIV is based on the numerical solution of the system of ODEs of the main problem of external ballistics and the method of bisection (division of a line segment in half-and-half).
The paper gives the Mathcad code of mathematical modeling of the problem of determining the initial velocity of an armor-piercing projectile of a 57-mm ZIS-2 anti-tank gun of the 1943 model, which can be modified for the benefit of the user.
The form of the dependence of the error in estimating the PIV on the error in determining the boundary conditions is established. There is a direct proportionality between SPIV and öIPV_RB. This allows us to specify
the requirements for the ADR for the IPV estimating by working in the artillery firing control system.
References
Alon, Y. 1989. Doppler radar method and apparatus for measuring a projectile's muzzle velocity. US4837718.
Burlov, V.V. et al., 2006. Ballistika stvol'nyh sistem. Moskva: Mashinostroenie (in Russian).
Chujko, V.S., 1976. Vneshnjaja ballistika. Moskva: Voenno-vozdushnaja inzhenernaja akademija (in Russian).
Germershausen, R., 1982. Handbook on weaponry. Düsseldorf: Rheinmetall GmbH.
Semenov, V.L., 2013. RLS izmerenija nachal'noj skorosti snarjada ustanovlennaja na stvole orudija. RU2501033 C1 (in Russian).
Solodov, A.P. & Ochkov, V.F., 2004. Differential Models: An Introduction with MathCAD. Springer.
Zaytsev, N.A., Iljuha, A.A. & Pyraev, V.V., 2004. Sposob izmerenija nachal'noj skorosti snarjada i ustrojstvo dlja ego osushhestvlenija. RU2250476 C2 (in Russian).
C\l CD !± ü
s=
.Q O
<U J3
ro
>
ro
T3
o
.Q
o ü
O
J
ro
M—
O c
o +-<
J3
o tn ro </) ro
'o o
<D >
<D
ü M—
O <U
ro E
tn LU
> >
o
CO sz
ф (Л
О <
ОЦЕНКА НАЧАЛЬНОЙ СКОРОСТИ СНАРЯДА КАК РЕШЕНИЕ ДВУХТОЧЕЧНОЙ ГРАНИЧНОЙ ЗАДАЧИ
Вадим Л. Хайков независимый исследователь, Краснодар, Российская Федерация
ОБЛАСТЬ: механика - баллистика
ВИД СТАТЬИ: оригинальная научная статья
о см
ЯЗЫК СТАТЬИ: английский
Е з
о Резюме:
Задача оценки начальной скорости снаряда сформулирована в виде
0 двухточечной граничной задачи. Для её решения могут быть использованы данные доплеровского измерителя скорости или
^ результаты решения задачи Коши. Приведен алгоритм оценки ч0>
базирующийся на совокупности численного решения системы ос дифференциальных уравнений (СДУ) полёта снаряда и метода
бисекции. Итерационный процесс оценки начальной скорости прерывается при достижении заранее назначенной величины погрешности и заблаговременно установленного значения ширины интервала поиска. В статье представлена блок-схема ^ алгоритма начального процесса скорости снаряда. Приведен
щ пример математического моделирования задачи определения
2 начальной скорости 57-мм бронебойного снаряда
'У противотанковой пушки ЗИС-2 обр. 1943 г. с использованием
1 программной среды МэМсэсС.
^ Ключевые слова: прикладная внешняя баллистика ствольных
г систем, начальная скорость снаряда, математическое и
о компьютерное моделирование, среда МэМсэС.
ПРОЦЕНА ПОЧЕТНЕ БРЗИНЕ ПРОJЕКТИЛА КАО РЕШЕТЕ ПРОБЛЕМА ГРАНИЧНЕ ВРЕДНОСТИ ЗА ДВЕ ТАЧКЕ
Вадим Л. Ха]ков
независни истраживач, Краснодар, Руска Федераци]а
ОБЛАСТ: механика - балистика
ВРСТА ЧЛАНКА: оригинални научни чланак
иЕЗИК ЧЛАНКА: енглески
Сажетак:
Процена почетне брзине про]ектила формулисана }е у облику проблема граничне вредности у две тачке. За ъегово решаваъе могу да се користе подаци доплер радара или резултати решеног
со !±
S= <u
тз
Коши]евог проблема. Процес процене почетне брзине про]ектила v0 гм заснива се на нумеричком решаваъу система обичних диференцщалних ]едначина и методи бисекцце. Итеративни рачунски процес се прекида када се постигну унапред дата тачност почетне брзине про]ектила, као и унапред одре^ена вредност ширине -§ интервала претраживаъа брзине. У чланку jе представлен блок- ^ ди}аграм алгоритма процеса почетне брзине про]ектила. Приказан }е код математичког моделираъа у Маткеду, као и компуутерска симулаци]а процеса процене почетне брзине 57-милиметарског панцирног про]ектила противоклопног топа ЗИС-2 из 1943. године.
Къучне речи: применена спорна балистика артиъерцског ору^а, процес процене почетне брзине про]ектила, математичко моделираъе, комп}утерска симулаци}а, Mathcad софтвер.
го
Paper received on / Дата получения работы / Датум приема чланка: 20.09.2017. о
Manuscript corrections submitted on / Дата получения исправленной версии работы / Датум достав^а^а исправки рукописа: 05.10.2017. Paper accepted for publishing on / Дата окончательного согласования работы / Датум ° коначног прихвата^а чланка за об]ав^ива^е: 07.10.2017.
© 2018 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/). ^
© 2018 Автор. Опубликовано в «Военно-технический вестник / Vojnotehnicki glasnik / Military Technical Courier» (www.vtg.mod.gov.rs, втг.мо.упр.срб). Данная статья в открытом доступе и распространяется в соответствии с лицензией «Creative Commons»
(http://creativecommons.org/licenses/by/3.0/rs/). ъ
ф
© 2018 Аутор. Обjавио Воjнотехнички гласник / Vojnotehnicki glasnik / Military Technical Courier 2°
(www.vtg.mod.gov.rs, втг.мо.упр.срб). Ово jе чланак отвореног приступа и дистрибуира се у складу са Creative Commons licencom (http://creativecommons.org/licenses/by/3.0/rs/). °
e t a
im ti
s
ш
v o ai h
s a
