内弹道计算程序

%59nian130

A=0.87; %枪(炮)膛横断面积A dm^2

G=19;%33.4; %弹重 kg

W0=2.04; %药室容积 dm^3

l_g=25.0; %身管行程 dm

P_0 =30000; %起动压力 kpa

fai1=1.02; %次要功系数

K=1.03; %运动阻力系数φ1

theta =0.2; %火药热力系数

%=========================================

f=950000; %火药力 kg*dm/kg

alpha=1; %余容 dm^3/kg

delta=1.6; %火药重度γ

%==================================

ome=2.2;%12.9; %第一种装药量 kg

u1=5.0024*10^-5; %第一种装药烧速系数 dm^3/(s*kg)

n1=0.82; %第一种装药的压力指数n1

lambda=-0.0071; %第一种装药形状特征量λ1

lambda_s=0; %第一种装药分裂点形状特征量λ1s

chi=1.00716; %第一种装药形状特征量χ1

chi_s=0; %第一种装药分裂点形状特征量χ1s

mu=0; %第一种装药形状特征量μ1

et1=1.14*10^-2; %第一种装药药厚δ01

d1=2.5*10^-2; %第一种装药火药内径d1

Ro1=0; %药型系数α1

%=========================================

%常数与初值计算----------------------------------------------------------------- l_0=W0/A;

Delta=ome/W0;

phi=K + ome/(3*G);

v_j=196*f*ome/(phi*theta*G);

v_j=sqrt(v_j);

B = 98*(et1*A)^2/( u1*u1*f*ome*phi*G );

B=B*(f*Delta)^(2-2*n1);

Z_s=1+Ro1*(d1/2+et1)/et1;

p_0=P_0/(f*Delta);

psi_0=(1/Delta - 1/delta)/(f/P_0 + alpha - 1/delta);

Z_0=(sqrt(1+4*psi_0*lambda/chi) - 1)/(2*lambda);

%解算子-----------------------------------------------------------------------

C = zeros(1,12);

C(1)=chi;C(2)=lambda;C(3)=lambda_s;C(4)=chi_s;C(5)=Z_s;%

C(6)=theta;C(7)=B;C(8)=n1;C(9)=Delta;C(10)=delta;C(11)=alpha;C(12)=mu; C;

y0=[Z_0;0;0;psi_0];