導(dǎo)彈燃料相對質(zhì)量因數(shù)計(jì)算

導(dǎo)彈總體設(shè)計(jì)《防空導(dǎo)彈主級燃料相對質(zhì)量因數(shù)計(jì)算》防空導(dǎo)彈主級燃料相對質(zhì)量因數(shù)計(jì)算《導(dǎo)彈總體設(shè)計(jì)》課程大作業(yè)班級姓名學(xué)號

目錄一、實(shí)驗(yàn)?zāi)康?3二、實(shí)驗(yàn)要求 3三、計(jì)算方法 3四、實(shí)驗(yàn)數(shù)據(jù) 61、分離條件 62、氣動(dòng)參數(shù) 63、發(fā)動(dòng)機(jī)參數(shù) 64、導(dǎo)引規(guī)律 75、大氣模型 7五、實(shí)驗(yàn)結(jié)果 71.燃料質(zhì)量因數(shù) 72.彈道圖形 73.攻角變化圖 8附程序源代碼 9

四、實(shí)驗(yàn)數(shù)據(jù)1、分離條件速度V=500m/s，時(shí)間t=3s，x方向位置x0=674m，y方向位置y0=329m，初始攻角α02、氣動(dòng)參數(shù)阻力因數(shù)插值表：（）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升力因數(shù)斜率插值表：（）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、發(fā)動(dòng)機(jī)參數(shù)比沖Is=2156；重力加速度g=9.801；推重比；翼載。4、導(dǎo)引規(guī)律三點(diǎn)法。目標(biāo)勻速直線等高迎頭飛行，，，（彈目斜距）。5、大氣模型ρ0=五、實(shí)驗(yàn)結(jié)果1.燃料質(zhì)量因數(shù)μ=0.3442.彈道圖形圖SEQ圖\*ARABIC2彈道圖形3.攻角變化圖圖SEQ圖\*ARABIC3攻角的變化情況

附程序源代碼%ByJRD%%clear;clc;u=0;%燃料相對質(zhì)量因數(shù)初始值%初始參數(shù)Is=2156;%發(fā)動(dòng)機(jī)參數(shù)g=9.801;P2W=2.2;P0=5880;Vt=420;%目標(biāo)參數(shù)yt=15000;Dt0=34200;q0=asin(yt/Dt0);V0=500;%導(dǎo)彈參數(shù)t0=3;x0=674;y0=329;Attack=1.5;sita0=26;s0=sita0*pi/180;T0=288.15;%大氣參數(shù)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);%音速，馬赫數(shù)計(jì)算Ma=V(i)/V_sound;Attack(1)=1.5;%阻力因數(shù)插值表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');%升力因數(shù)插值表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');%數(shù)值積分%四級四階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(