内弹道计算程序(国外英文资料).docVIP

内弹道计算程序(国外英文资料) % 59nian130 = 0.87；%枪（炮）膛横断面积DM ^ 2 G = 19；% 33.4；%弹重公斤 W0 = 2.04；%药室容积DM ^ 3 l_g = 25；%身管行程DM p_0 = 30000；%起动压力kPa fai1 = 1.02；%次要功系数 k = 1.03；%运动阻力系数φ θ= 0.2；%火药热力系 % ========================================= F = 950000；%火药力公斤的DM /公斤 α= 1；%余容DM ^ 3 /公 δ= 1.6；%火药重度 % ================================== F = 2.2；% 12.9；%第一种装药量公斤 U1＝5.0024×10 ^ - 5；%第一种装药烧速系数DM ^ 3 /（S *公斤） N1 = 0.82；%第一种装药的压力指数N1 λ= -0.0071；%第一种装药形状特征量 lambda_s = 0；%第一种装药分裂点形状特征量λ1 χ= 1.00716；%第一种装药形状特征量 chi_s = 0；%第一种装药分裂点形状特征量χ1 亩= 0；%第一种装药形状特征量μ ET1 = 1.14 * 10 ^ - 2；%第一种装药药厚δ0 D1 = 2.5 * 10 ^ - 2；%第一种装药火药内径D1 RO1 = 0；%药型系数α % ========================================= 常数与初值计算----------------------------------------------------------------- % l_0 = W0 /； 三角洲=青梅/ w0； φ= k（3×g）； v_j = 196 * F *人/（φ×θ* G） v_j = SQRT（v_j）； B = 98 *（ET1 *）^ 2 /（U1 * * *人*披U1 F * G）； b =（* *δ） z_s = 1 + RO1 *（D1 / 2 + ET1）/ ET1； p_0 = p_0 /（F *三角）； psi_0 =（1 / 1 /三角洲）/（F / p_0 +α- 1 /三角洲） z_0 =（sqrt（1 + 4 * psi_0×λ/卡）1）/（2×λ）； 解算子----------------------------------------------------------------------- % C =零（1,12）； C（1）=智；C（2）=λ；C（3）= lambda_s；C（4）= chi_s；C（5）= z_s； c（6）＝θ；c（7）＝b；c（8）＝N1；c（9）＝δ；C（10）＝δ；C（11）＝α；C（12）= C； y0 = [ z_0；0；0；psi_0 ]； 选项= odeset（outputfcn ，odeplot”）； [译文]，Y = ode45（@ ndd_fun，0:100，[ 0 ]，z_0；0；选项C）； l＝y（：，2）； L = L * l_0； FL =找到（L＞＝l_g）； min（FL）； [译文]，Y = ode45（@ ndd_fun，0:0.005：FL，[ 0 ]，z_0；0；选项C）； Z = Y（：，1）；LX = Y（：，2）；VX = Y（：，3）； ψ=（Z = 0和z＜1）。*（智* Z. *（1 +λ* Z +μ×Z））+…%%%%%%% （Z = 1和z z_s）。*（chi_s * Z *（1 + lambda_s * Z））+… （Z = z_s）* 1； l_psi = 1 -（Δ/Δ）*（1-psi）-α×δ×PSI PX =（PSI，VX，VX）/（LX + l_psi）； p = PX * *δ/ 100 V = Vx * v_j / 10； L = LX * l_0； T = TT * l_0 * 1000 / v_j； FL =找到（L＞＝l_g）； FL = min（FL）+ 1； p（FL；结尾）= []；v（FL；结束）= []；l（FL；结束）= []；T（FL：结束）]； PD = PX * F *三角洲/ 100 /（1 +一/ 3 / fai1／g）； Pt = PD *（1 +一/ 2 / fai1／g）； 最大值（PX）； m =查找（PX = aa）； PM = [ TT（M）* l_0 * 1000 / v_j LX（M）* l_0 VX（M）* v_j / 10像素（M）* F *三角/ 100角（M）PD（M）PSI（m）Z（M）]； % =长度（TT）； 查找（z = 1）； RAN