导弹燃料相对质量因数计算

导弹总体设计《防空导弹主级燃料相对质量因数计算》防空导弹主级燃料相对质量因数计算《导弹总体设计》课程大作业班级姓名学号

目录一、实验目的 3二、实验要求 3三、计算方法 3四、实验数据 61、分离条件 62、气动参数 63、发动机参数 64、导引规律 75、大气模型 7五、实验结果 71.燃料质量因数 72.弹道图形 73.攻角变化图 8附程序源代码 9

四、实验数据1、分离条件速度V=500m/s，时间t=3s，x方向位置x0=674m，y方向位置y0=329m，初始攻角α02、气动参数阻力因数插值表：（）2468101.50.04300.05110.06510.08470.11202.10.03600.04360.05580.07360.09732.70.03080.03720.04810.06410.08493.30.02650.03230.04190.05600.07464.00.02220.02720.03560.04780.0644升力因数斜率插值表：（）12468101.50.03020.03040.03060.03090.03110.03132.00.02790.02800.02840.02860.02880.02902.50.02610.02640.02670.02690.02720.02743.00.02470.02480.02510.02540.02570.02593.50.02260.02270.02310.02330.02360.02384.00.02090.02100.02130.02160.02190.0221注：以上攻角单位均为度。3、发动机参数比冲Is=2156；重力加速度g=9.801；推重比；翼载。4、导引规律三点法。目标匀速直线等高迎头飞行，，，（弹目斜距）。5、大气模型ρ0=五、实验结果1.燃料质量因数μ=0.3442.弹道图形图SEQ图\*ARABIC2弹道图形3.攻角变化图图SEQ图\*ARABIC3攻角的变化情况

附程序源代码%ByJRD%%clear;clc;u=0;%燃料相对质量因数初始值%初始参数Is=2156;%发动机参数g=9.801;P2W=2.2;P0=5880;Vt=420;%目标参数yt=15000;Dt0=34200;q0=asin(yt/Dt0);V0=500;%导弹参数t0=3;x0=674;y0=329;Attack=1.5;sita0=26;s0=sita0*pi/180;T0=288.15;%大气参数r0=1.2495;H=y0;i=1;whileH<=ytx(1)=x0;y(1)=y0;s(1)=s0;V(1)=V0;ify(i)<=11000%不同高度的温度和密度T=288.15-0.0065*y(i);r=r0*(T/T0)^4.25588;elseT=216.65;r=0.36392/exp((y(i)-11000)/6341.62);endV_sound=20.0468*sqrt(T);%音速，马赫数计算Ma=V(i)/V_sound;Attack(1)=1.5;%阻力因数插值表attack1=[246810];Mach1=[1.5,2.1,2.7,3.3,4.0];cx=[0.0430,0.0511,0.0651,0.0847,0.1120;0.0360,0.0436,0.0558,0.0736,0.0973;0.0308,0.0372,0.0481,0.0641,0.0849;0.0265,0.0323,0.0419,0.0560,0.0746;0.0222,0.0272,0.0356,0.0478,0.0644];Cx=interp2(attack1,Mach1,cx,Attack(i),Ma,'spline');%升力因数插值表attack2=[1,2,4,6,8,10];Mach2=[1.5,2.0,2.5,3.0,3.5,4.0];cy=[0.0302,0.0304,0.0306,0.0309,0.0311,0.0313;0.0279,0.0280,0.0284,0.0286,0.0288,0.0290;0.0261,0.0264,0.0267,0.0269,0.0272,0.0274;0.0247,0.0248,0.0251,0.0254,0.0257,0.0259;0.0226,0.0227,0.0231,0.0233,0.0236,0.0238;0.0209,0.0210,0.0213,0.0216,0.0219,0.0221];Cy=interp2(attack2,Mach2,cy,Attack(i),Ma,'spline');%数值积分%四级四阶Runge-Kutta法h=0.001;%步长K1=Is/(1-u)-r*V(i)^2*Cx*Is/(2*P2W*P0*(1-u))-Is*sin(s(i))/P2W;K2=Is/(1-(u+h/2))-r*(V(i)+h/2*K1)^2*Cx*Is/(2*P2W*P0*(1-(u+h/2)))-Is*sin(s(i))/P2W;K3=Is/(1-(u+h/2))-r*(V(i)+h/2*K2)^2*Cx*Is/(2*P2W*P0*(1-(u+h/2)))-Is*sin(s(i))/P2W;K4=Is/(1-(u+h))-r*(V(i)+h*K3)^2*Cx*Is/(2*P2W*P0*(1-(u+h)))-Is*sin(s(i))/P2W;V(i+1)=V(i)+(K1+2*K2+2*K3+K4)*h/6;y(i+1)=y(i)+Is/(P2W*g)*V(i)*sin(s(i))*h;H=y(i+1);x(i+1)=x(i)+Is/(P2W*g)*V(i)*cos(s(i))*h;p=P2W;p0=P0;Ks1=Vt/yt*(Is/(P2W*g)+(V(i)*Is/(P2W*g)*V(i)*sin(s(i))-y(i)*(K1+2*K2+2*K3+K4).../6)/V(i)^2/sin(s(i)))/(1+(cot(q0)-u*Is*Vt/P2W/g/yt)*cot(s(i)));Ks2=Vt/yt*(Is/(P2W*g)+(V(i)*Is/(P2W*g)*V(i)*sin(s(i)+Ks1*h/2)-y(i)*(K1+2*K2+2*K3+K4).../6)/V(i)^2/sin(s(i)+Ks1*h/2))/(1+(cot(q0)-(u+h/2)*Is*Vt/P2W/g/yt)*cot(s(i)+Ks1*h/2));Ks3=Vt/yt*(Is/(P2W*g)+(V(i)*Is/(P2W*g)*V(i)*sin(s(i)+Ks2*h/2)-y(i)*(K1+2*K2+2*K3+K4).../6)/V(i)^2/sin(s(i)+Ks2*h/2))/(1+(cot(q0)-(u+h/2)*Is*Vt/P2W/g/yt)*cot(s(i)+Ks2*h/2));Ks4=Vt/yt*(Is/(P2W*g)+(V(i)*Is/(P2W*g)*V(i)*sin(s(i)+Ks3*h)-y(i)*(K1+2*K2+2*K3+K4).../6)/V(i)^2/sin(s(i)+Ks3*h))/(1+(cot(q0)-(u+h)*Is*Vt/P2W/g/yt)*cot(s(i)+Ks3*h));s(i+1)=s(i)+(Ks1+2*Ks2+2*Ks3+Ks4)*h/6;Attac(i)=(V(i)*(Ks1+2*Ks2+2*Ks3+Ks4)/6+Is*cos(s(i))/P2W)/(r*V(i)^2*Cy*180/pi*Is.../(2*P2W*P0*(1-u))+Is/(1-u));Attack(i+1)=Attac(i)*180/pi;u=u+h;i=i+1;endu=u-h;Attack(1)=Attack(2);%图像绘制figure(