弹道计算大作业.docxVIP

弹道计算大作业 目录 TOC \o 1-5 \h \z \o Current Document 一、 初始条件和要求 2 \o Current Document 1.1初始条件 2 \o Current Document 1.2仿真要求 2 \o Current Document 二、 模型的建立 2 \o Current Document 2.1升力和阻力模型 2 \o Current Document 2.2大气和重力加速度模型 3 \o Current Document 2.3无控飞行 3 \o Current Document 2.4平衡滑翔 4 \o Current Document 2.5最大升阻比滑翔飞行弹道 4 \o Current Document 三、 仿真结果 5 \o Current Document 3.1无控飞行弹道仿真 5 \o Current Document 3.2平衡滑翔弹道仿真 7 \o Current Document 3.3最大升阻比滑翔弹道仿真 8 \o Current Document 附录 9 一、初始条件和要求 1丄初始条件 已知给定的初始条件如下: 表1初始条件 名称 符号 给定值 单位 参考面积 S 1.7 m2 展弦比 A 0. 86 / 效率因子 e 0.9 / 质量 m 115 炖 重力加速度 g 通过模型计算可得 m/s2 零升阻力系数 0. 02 /rad 密度 P 通过模型计算可得 kg/m3 1.2仿真要求 请使用Simulink或Buildfly完成以下仿真任务： （1） 请完成该导弹的无控飞行弹道仿真； （2） 请完成该导弹的平衡滑翔方案飞行弹道仿真； （3） 请完成该导弹的最大升阻比滑翔飞行弹道仿真; 二、模型的建立 2.1升力和阻力模型 已知展弦比为2的飞行器的升力线斜率为: 根据飞行力学相关知识，飞行器的升力系数和阻力系数为: 其中，升力线斜率由（1）式可得；￡为效率系数：￡ = — o 7ceA 由升力系数和阻力系数，得到导弹的升力和阻力为： * x = c#s 2 ⑶ Y = Cy^pv2S 2.2大气和重力加速度模型 在计算过程屮，大气密度采用如下模型： / \ 4.25588 L 0.0065/7 1 P=Po 1——- \ ‘0丿 其中,po = \225kgl品为海平面的大气密度;7； =288.15K。 重力加速度采用如下模型： （Rd 丫g = g° Ttn E + H丿 其中，g（）=9.8, /?t/ =6371000/??为地球半径；H为飞行器距离地面的高度。 2.3无控飞行 假设导弹的运动始终在铅垂平面，根据飞行力学知识，得到导弹无控飞行时 的运动学和动力学方程为: S S O C g 在上述模型屮，假设俯仰角〃为0。 2.4平衡滑翔 所谓的“平衡可以理解为垂直于速度方向受力平衡，即矿0。因此得到 平衡滑翔时的导弹运动学和动力学方程: dV~dt dV ~dt 兰-gsin m =V cos =Vsi n 由于弹道倾角的变化率为常数，方程组小的第二个方程等于0。这个方程可 以用来求攻角 2.5最大升阻比滑翔飞行弹道 联立(1)式、(2)式可得升阻比的表达式为: 兀九a 1 + 兀九a 1 +J1 + 丿 、2 从上式可以看岀，由于展弦比2、零升阻力系数Cg为常数，因此升阻比只 得到和攻角有关，是关于攻角的函数。因此要使升阻比达到最大，须使 得到 因此，以最大升阻比滑翔时导弹运动学和动力学方程为: dV X . q —= ￡Sin dt m (10)_ Y gcos 0 mV (10) =V cos =V sin 0 三、仿真结果 3.1无控飞行弹道仿真 根据无控弹道模型，写出S函数，搭建的仿真模块如卜?图所示: Stop SimulationScopebCompare To Zero Stop Simulation Scopeb Compare To Zero 图1无控飞行仿真模块 由于初始条件给定，因此模块没有输入；输出有六个，分别为导弹的射程变 化、高度变化、速度变化、弹道倾角变化、攻角变化以及密度变化。模块的仿真 吋间由高度变化决定，当高度降为0（导弹落到地面上）吋仿真结束。导出数据 后画图如下：  图2无控飞行时各参数变化 3.2平衡滑翔弹道仿真 平衡滑翔弹道仿真模块如下图所示: 图3平衡滑翔模块 取仿真时间为150s,无输入，输出分别为：导弹的射程变化、高度变化、速度变化、弹道倾角变化、攻角变化以及密度变化。得到各参量时间变化图如下: \ / 200 / ■ 「\ : ISO / / \ / \ ■ / \ I 140 120 / —