航天飞行动力学远程火箭弹道设计大作业.docVIP

航天飞行动力学远程火箭弹道设计大作业 已知火箭纵向运动方程式如公式 GOTOBUTTON ZEqnNum916441 REF ZEqnNum916441 \* Charformat \! (1)所示。 其中，分别为火箭飞行速度、发动机推力、火箭初始质量、弹道倾角、攻角、水平位移和飞行高度；为角度增益系数，为火箭飞行时间，为火箭质量。仿真初始条件如表1和表2所示。 表1初始状态 序号 变量名 变量值 物理意义及单位 0 0 火箭飞行时间，s 1 初始弹道倾角，弧度 2 0 火箭初始速度， 4 0 火箭在地面发射坐标系下的初始水平位置,m 5 0 火箭在地面发射坐标系下的初始高度,m 表2 有关参数 序号 变量名 变量值 物理意义及单位 0 8000 起飞质量 1 28.57 单位时间燃料质量消耗， 2 9.8 重力加速度常数， 3 35 角度增益系数 4 200 发动机推力，KN 5 7000 发动机排气速度， 飞行程序角随火箭飞行时间的关系如公式 GOTOBUTTON ZEqnNum788194 REF ZEqnNum788194 \* Charformat \! (2)： 问题：（1）请根据如上已知条件，完成火箭纵向运动仿真。（2）验证齐奥尔科夫斯基公式。 MATLAB程序： t=0; %火箭飞行时间 theta=pi/2; %初始弹道倾角 v=0; %火箭初始速度 x=0; %火箭初始水平位置 y=0; %火箭初始高度 vl=0; %引力引起的速度损失 m0=8000; %起飞质量 m=m0; md=28.57; %质量秒耗量 g=-9.8; %重力加速度 A=35; %角度增益系数 P=200000; %推力 w=7000; %排气速度 t1=10; t2=130; fig=pi/60; i=1; dt=0.01; while(t(i)=150) if(t(i)t1) fi(i)=pi/2; elseif(t(i)=t1t(i)t2) fi(i)=pi/2+(pi/2-fig)*(((t(i)-t1)/(t2-t1))^2-2*(t(i)-t1)/(t2-t1)); else fi(i)=fig; end a(i)=A*(fi(i)-theta(i)); vv(i)=P/m(i)+g*sin(theta(i)); if(i==1) ththeta(i)=0; else ththeta(i)=(P/m(i)*a(i)+g*cos(theta(i)))/v(i); end xd(i)=v(i)*cos(theta(i)); yd(i)=v(i)*sin(theta(i)); m(i+1)=m(i)-md*dt; t(i+1)=t(i)+dt; v(i+1)=v(i)+vv(i)*dt; theta(i+1)=ththeta(i)*dt+theta(i); x(i+1)=x(i)+xd(i)*dt; y(i+1)=y(i)+yd(i)*dt; vl(i+1)=vl(i)-g*sin(theta(i))*dt; i=i+1; end mk=m(i); vt=v(i); plot(t,y); title(高度随时间的变化); xlabel(t/s); ylabel(H/m); %验证齐奥尔科夫斯基公式 vi1=vt+vl(i) %无速度损失的理想速度 vi2=-w*log(mk/m0) %齐奥尔科夫斯基公式理想速度 dv=vi1-vi2 vi1 = 5.3705e+03 vi2 = 5.3704e+03 dv = 0.1243 因为 Vi1=5370.50655450199 Vi2=5370.38223449984 =0.1243 所以齐奥尔科夫斯基公式正确。