弹丸刚体外弹道及其引信力学环境数值解算稳定性

针对弹丸刚体外弹道常微分方程组及引信惯性零部件轴向过载求解稳定性问题,借助计算机,以不同数值计算方法和不同积分时间步长对弹丸刚体外弹道常微分方程组及引信惯性零部件轴向过载求解,得到刚体外弹道常微分方程组的数值解,即刚体外弹道含质心运动和绕质心运动的诸元,在此基础上又得到了引信惯性零部件轴向过载的数值解。通过大量仿真发现:应用龙格-库塔法、阿当姆斯预报-校正法和吉尔法解算刚体外弹道常微分方程组,当积分时间步长取为0.05 s时,质心运动诸元解算结果均收敛;在积分时间步长取10^(-4) s时,绕心运动诸元解算结果也均收敛。而在解算引信惯性零部件轴向过载突变峰值时,发现若使峰值稳定,这3种方法均需调整积分时间步长小于或等于10^(-5) s。而求解引信惯性零部件轴向过载,只有吉尔法未失真。