航天飞行动力学大作业

问题描述

导弹质量$m_0=320kg,$转动惯量$J_z=315kg^.m^2,$发动机推力$P=2000N,$质量秒消耗量$m_s=0.46kg/s,$初始速度$V_0=250m/s,$初始位置$x_0=0m,$初始高度$H_0=7000m,$初始弹道倾角$\theta=0^。,$初始俯仰角$\phi_0=0^。,$初始攻角$\alpha_0=0^。,$初始俯仰角速度$\dot\phi_0=0rad/s,$初始速度$V_0=250m/s,$参考面积$S_{ref}=0.45m^2,$参考长度$L_{ref}=2.5m,$升力系数 $C_y=0.25\alpha+0.05\delta_z,$(角度单位为度).阻力系数$C_x=0.2+0.005\alpha^2,$(角度单位为度).俯仰力矩系数$m_z=-0.1\alpha+0.024\delta_z,$ (角速度单位rad/s).阻尼力矩系数$m_z^{\omega_z}=-2,$ (角速度单位rad/s).升力$Y=C_yq{S_ref},$阻力$X=C_xqS_{ref},$俯仰力矩$M_z=m_zqS_{ref}L_{ref},$动压$q=\frac{1}{2}{\rho}V^2,$大气密度$\rho$计算公式：

$$\rho_0=1.2495$$ $$T_0=288.15$$ $$T=T_0-0.0065H$$ $$\rho=\rho_0(\frac{T}{T_0})^{4.25588}$$
1）飞行方案 $(x<9100m)$
$$H^*=2000cos(0.000314*1.1x)+5000$$ $$\delta_z=k_\phi*(H-H^*)+k_\dot\phi*(\dot H-\dot H^*)$$ $$m_s=0.0kg/s$$
2）飞行方案 $(9100m<x<24000m)$
$$H^*=3050m$$ $$\delta_z=k_\phi*(H-H^*)+k_\dot\phi*\dot H$$ $$m_s=0.46kg/s$$
3）飞行方案 $(x>24000m且y>0m)$
目标位置为 $x_m=30000m$
采用比例导引法

弹道部分具体要求
1） 计算纵向理想弹道，采用瞬时平衡假设
$$m_z^\alpha\alpha+m_z^{\delta_z}\delta_z=0$$
2）不考虑气动力下洗影响，计算飞行器的纵向短周期扰动的动力系数，并分析其在各时间点的自由扰动的稳定性，并计算在各个特征秒弹体传递函数 $W_{\phi\delta}(s),W_{n_y\delta}(s)$

根据理想弹道的调试好的参数按下列要求,进行计算,看看结果变化,分析一下原因
3) 三自由度计算(将瞬时平衡方程更换为旋转动力学方程)

$$J_z\frac{d\omega_z}{dt}=M_z$$
4)考虑舵机和高度表的传递函数
舵偏传递函数：
$$G_{\delta_z}=\frac{\omega_n^2}{s^2+2\xi\omega_ns+\omega_n^2}$$ $$|\delta_z<=15^。|$$
$$\omega_n=150$$ $$\xi=0.7$$
高度表传递函数：
$$G_H=\frac{1}{Ts+1}$$ $$T=0.05s$$

问题解答

详细实验报告及代码见链接：航天飞行动力学大作业