-航天飞行机设计作业报告-导弹燃料相对质量因数计算.docxVIP

航天飞行器设计作业报告 ——导弹燃料相对质量因数计算 一、问题描述 1、根据已知条件，采用数值积分法求解相对量运动微分方程组，计算燃料相对质量因数。 2、综合运用积分、插值计算等计算方法，采用C、C++，或者Matlab等语言的一种，编制计算程序。 二、模型建立 三、弹道计算 利用python程序可以迭代得到如下的参数变化图像： 最终计算得导弹的相对质量因数为: 0.3453 四、python程序如下 rnok.py文件： import Initialize import Interpolate import ShowResult import numpy from scipy import interpolate from numpy import sin,cos,tan #初始化获得各参数及插值表及迭代初值 (alp1,ma1,cx0,alpha2,ma2,cy0,yh,ph,ah)=Interpolate.interpolation(); (u,t,xt,v,x,y,Drx,cx,cya,alpj,alph,thtj,thth,v0,t0,x0,y0,alp0,tht0,Is,g,P,p0,vt,dt0,yt0,xt0,a,cotq0)=Initialize.initialize(yh,ah); #迭代参量及迭代步长 i=0; h=0.0001; # 龙格库塔迭代 while (y[i]yt0): ma=v[i]/a; cotq=cotq0-u[i]*Is*vt/[P*g*yt0]; p=interpolate.interp1d(yh,ph)(y[i]);#插值得到密度 a=interpolate.interp1d(yh,ah)(y[i]);#插值得到音速 CCX=interpolate.interp2d(alp1,ma1,cx0,kind=linear); cx[i] =CCX(alpj[i],ma); CCYA=interpolate.interp2d(alpha2,ma2,cy0,kind=linear); cya[i] =180/numpy.pi*CCYA(alpj[i],ma); Kv1 = Is / (1 - u[i]) - p * v[i] ** 2 * cx[i] * Is / (2 * P * p0 * (1 - u[i])) - Is / P * sin(thth[i]); Kv2 = Is / (1 - (u[i] + h / 2)) - p * (v[i] + h / 2 * Kv1) ** 2 * cx[i] * Is / ( 2 * P * p0 * (1 - (u[i] + h / 2))) - Is / P * sin(thth[i]); Kv3 = Is / (1 - (u[i] + h / 2)) - p * (v[i] + h / 2 * Kv2) ** 2 * cx[i] * Is / ( 2 * P * p0 * (1 - (u[i] + h / 2))) - Is / P * sin(thth[i]); Kv4 = Is / (1 - (u[i] + h)) - p * (v[i] + h * Kv3) ** 2 * cx[i] * Is / ( 2 * P * p0 * (1 - (u[i] + h))) - Is / P * sin(thth[i]); Kth1 = vt / yt0 * (Is / (P * g) + ( v[i] / (P * g) * Is * v[i] * sin(thth[i]) - y[i] * (Kv1 + 2 * Kv2 + 2 * Kv3 + Kv4) / 6) / ( v[i] ** 2 * sin(thth[i]))) / ( 1 + (cotq0 - u[i] * Is * vt / (P * g * yt0)) / tan(thth[i])); Kth2 = vt / yt0 * (Is / (P * g) + (v[i] / (P * g) * Is * v[i] * sin(thth[i] + Kth1 * h / 2) - y[i] * ( Kv1 + 2 * Kv2 + 2 * Kv3 + Kv4) / 6) / (v[i] ** 2 * si