MATLAB零维内弹道

实验三固体火箭发动机零维内弹道计算

M文件：

function dy=neidandao(t,y);

dy=zeros(4,1);

rou=y(1);

p=y(2);

Vc=y(3);

e=y(4);

d0=0.016;

h0=0.08;

D0=0.03;

rougr=1750;

k=1.17;

R=300;

Tp=3200;

b=0.002411;

n=0.315;

At=pi*16*10^(-6);

fai=0.95;

ka=0.98;

r=b*(p/1.013/10^5)^n;

gama=(2/(k+1))^((k+1)/(2*(k-1)))*sqrt(k);

c=sqrt(R*Tp)/gama;

if e<=(D0-d0)/2

Ab=pi*(d0+2*e)*h0;

else Ab=0;

end

dy(1)=(1/Vc)*((rougr-rou)*Ab*r-(fai*p*At)/(c*sqrt(ka)));

dy(2)=(1/Vc)*(rougr*Ab*r*k*R*ka*Tp-(fai*p*p*At*k)/(rou*c*sqrt(ka))-p*Ab*r); dy(3)=Ab*r;

dy(4)=r;

end

主程序：

>> [t,y]=ode45('neidandao',[0:0.00001:1.5],[1.29;101300;pi*0.008*0.008*0.08;0]); >> plot(t,y(:,2))

P-t 曲线

00.51 1.5

00.5

1

1.5

2

2.5

3

6

初始段P-t 曲线

00.0050.010.015

02

4

6

8

10

12

14

5

燃烧终了段P-t 曲线

实验总结

这次实验是我们对所学课程即固体火箭发动机零维内弹道计算以及MATLAB 软件的一次练习。通过这次实验，我了解了常微分方程组数值解法的一般过程，掌握了用MATLAB 软件的具体实现方法，得到了零维内弹道压强曲线，完成了实验的要求。

在MATLAB 中，实现常微分方程组数值解法的是ode 函数（在本实验中用的是ode45），它不需要用户自己编程，使用起来比较简单，总体来说这次实验也完成的比较顺利。

1.08 1.09 1.1 1.11 1.12 1.13 1.14 1.15 1.16 1.17

2.8

2.82

2.84

2.86

2.88

2.9

2.92

2.94

2.96

2.98

6