2009302026王勋..docVIP

航空外弹道学 课程设计 姓 名 ： 王勋 班 级 ： 学 号 ： 2009302026 时 间 ： 2012.7.12 一、已知条件及题目要求 查表可知，标准下落时间 2．气象： 温度梯度 弹丸参数： 弹长 空气阻力系数： 初始条件： 要求：列出弹道参数，并画出炸弹弹道曲线。 二、题目分析 对于所给题目，取直角坐标系Oxy，坐标原点取在投弹高度为H的投弹点O上，x轴取在飞机投弹瞬间速度的铅垂面（投弹面）内的水平方向，y轴铅直向下，如图1所示。转角方向规定顺时针为正。设弹道上任一点M速度向量在x轴与y轴上的分量分别为u与w。 图1 水平轰炸图 由题目给出的条件，P=0，W=0可以知道，所要设计的题目类型为无外力的水平轰炸，由于飞机速度V在X轴方向上，因而初始条件为： 当时，建立方程组如下： 三、相关参数计算 (1)弹道系数C的计算 其中d为弹直径，q为炸弹的重量。 (2)， 其中 (3)查表得，其中： 可作如下近似 = （4）的计算 其中由（3）可求得，τ由（2）计算可求得。 （5）的计算 (6)气压的计算： 其中 (7)音速的计算： (8)密度的计算： 其中为空气比重，为地面空气比重标准值。 四、方程变形 由如下方程： 令 把上述方程组改变为： 积分的初始条件：当时， 带入上述转换化简得： 五、运用RungeKutta法求解 RungeKutta法的主要优点是计算精度较高，能满足通常的计算要求；每次计算，只用到前一步的计算结果；计算过程中改变步长不受限制。 求解C代码如下： #include math.h #include stdio.h #includestdlib.h #define H 2000 #define G 5.682e-3 #define R 29.27 #define tao0 288.4 #define C 0.10165 #define U 400 #define g 9.806 #define dt 0.01 //求相关参数； double tao(double y2) //求τ { double tao_; tao_=tao0-G*(H-y2); return(tao_); } double Hg(double y2) //比重函数Hg { double Hg_; Hg_=pow((1-G/tao0*(H-y2)),1/(R*G))*sqrt(tao0/tao(y2)); return(Hg_); } double V(double y0,double y1) //速度 { double V_; V_=sqrt(y0*y0+y1*y1); return(V_); } double Cx0(double y0,double y1,double y2) //弹丸阻力系数Cx0 { double a; double Cx0_; a=20.074*sqrt(tao(y2));//音速为a if(V(y0,y1)/a0.8) Cx0_=0.6; else Cx0_=V(y0,y1)/a-0.1; return(Cx0_); } double Gt(double y0,double y1,double y2) //空气阻力函数 { double Gt_; Gt_=4.737e-4*V(y0,y1)*Cx0(y0,y1,y2)*sqrt(tao0/tao(y2)); return(Gt_); } //x速度u，y为P，z为高度 double f1(double x, double y, double z) //方程组（1）式 { return (-1)*C*Hg(z)*Gt(x,x*y,z)*x; } double f2( double x) //方程组（2）式 { return g/x; } double f3( double x, double y) //方程组（3）式 { return x*y; } double f4(double x) //方程组（4）式 { return 1.0*x; } double