Масса снаряда

Зависимость от массы очень похожа на зависимость от баллистического коэффициента(БК), например увеличение массы вдвое эквивалентно уменьшению вдвое БК. Однако при этом необходимо помнить что для придания снаряду большей массы той же начальной скорости что и снаряду стандартной массы необходимо увеличение количества взрывчатого вещества (пороха) в снаряде (здесь под снарядом подразумевается не только «болванка», полет которой мы изучаем, но и гильза, ВВ и т.п.). Кроме того скорее всего потребуется упрочнение орудия, что может быть признано нерациональным для орудия данной категории и калибра.

Листинг (m-файл)

function polet

 

clear all

 

global Cx S Pv m a V0 g y

 

V0 = 680; % начальная скорость

a = 45; % угол бросания (возвышения орудия)

g = 9.8; % ускорение свободного падения

m = 6.3; % масса снаряда

cal = 76; % калибр (в мм)

S = pi * ( cal / 1000 ) ^ 2 % площадь поперечного сечения снаряда (мидель)

Cx = 0.2; % аэродинамический(баллистический) коэффициент

Pv = 1.225; % плотность воздуха (на уровне моря)

NU = [0; 0]; % координаты в начальный момент времени

t = 100; % время, в течение которого ведется расчет

 

colour = ['r' 'g' 'b' 'c' 'm'] % цвета графиков

kV0 = [340 580 680 780 1360] % начальная скорость

ka = [25 35 45 55 65] % угол бросания (возвышения орудия)

kCx = [0.1 0.15 0.2 0.25 0.3] % аэродинамический(баллистический) коэффициент

km = [3.15 5.3 6.3 7.3 12.6] % масса снаряда

 

 

% Зависимость от Vo

figure

hold

for i = 1:5

V0 = kV0(i)

[T2,Y2]=ode45(@Complex, [0 t], NU);

plot(Y2(:, 1), Y2(:, 2), colour(i))

end

title 'Зависимость от Vo'

xlabel 'Дальность'

ylabel 'Высота'

legend 'Vo = 340' 'Vo = 580' 'Vo = 680' 'Vo = 780' 'Vo = 1360'

axis([0 10000 0 10000])

grid

 

V0 = kV0(3)

 

 

% Зависимость от a

figure

hold

for i = 1:5

a = ka(i)

[T2,Y2]=ode45(@Complex, [0 t], NU);

plot(Y2(:, 1), Y2(:, 2), colour(i))

end

title 'Зависимость от a'

xlabel 'Дальность'

ylabel 'Высота'

legend 'a = 25' 'a = 35' 'a = 45' 'a = 55' 'a = 65'

axis([0 10000 0 10000])

grid

 

a = ka(3)

 

 

% Зависимость от Cx

figure

hold

for i = 1:5

Cx = kCx(i)

[T2,Y2]=ode45(@Complex, [0 t], NU);

plot(Y2(:, 1), Y2(:, 2), colour(i))

end

title 'Зависимость от Cx'

xlabel 'Дальность'

ylabel 'Высота'

legend 'Cx = 0.1' 'Cx = 0.15' 'Cx = 0.2' 'Cx = 0.25' 'Cx = 0.3'

axis([0 10000 0 10000])

grid

 

Cx = kCx(3)

 

 

% Зависимость от m

figure

hold

for i = 1:5

m = km(i)

[T2,Y2]=ode45(@Complex, [0 t], NU);

plot(Y2(:, 1), Y2(:, 2), colour(i))

end

title 'Зависимость от m'

xlabel 'Дальность'

ylabel 'Высота'

legend 'm = 3.15' 'm = 5.3' ' m = 6.3' 'm = 7.3' 'm = 12.6'

axis([0 10000 0 10000])

grid

 

 

function Complex = Complex(t, x)

global Cx S Pv m a V0 g

Complex = [ ( 2*m*V0*cosd(a) ) / ( Cx*S*Pv*V0*cosd(a).*t + 2*m ) ; ( 2*m*V0*sind(a) ) / ( Cx*S*Pv*V0*sind(a).*t + 2*m ) - g.*t ];