飞行力学弹道设计与分析

飞行力学弹道设计与分析一．题目重述..........................................................................................................................................1二．第一次作业......................................................................................................................................21.求解思路...........................................................................................................................................22.计算结果：（在MATLAB环境下编程求解）.................................................................................3三．第二次作业......................................................................................................................................41.要求...................................................................................................................................................42.计算结果展示及问题分析...............................................................................................................53.精度分析...........................................................................................................................................7四．第三次作业......................................................................................................................................81.作业要求...........................................................................................................................................82.简单思路分析...................................................................................................................................83.计算结果及简单问题分析...............................................................................................................84.下面我们讨论步长DT对结果的影响...........................................................................................10............................................................................................................................................................10五．第四次作业....................................................................................................................................111.作业要求.........................................................................................................................................112.简要分析与计算结果.....................................................................................................................12一．题目重述某一型号导弹，给定参数如下0m=320kg;P=2000N;sm=0.46kg/s;0v=250m/s;0x=0m;0H=7000m0=0o;0=0o;0=0o;Sref=0.452m;Lref=2.5m;0.250.05Cyz;20.20.005Cx;0.10.024zmz;(气动计算里的角度单位都为度)yYCqSref;xXCqSref;zMzmqSrefLref;动压20.5qv;其中01.24950288.15T00.0065TTH;00=TT4.25588（）;现要求设计,KK两个参数在瞬时平衡假设下满足一下飞行方案(飞行中舵偏角015z)1.9100xm0.0/smkgs*2000cos(0.0003141.1)5000Hx**()()zKHHKHH2.910024000mxm*3000Hm(原题3050m感觉不太合理，这里改为3000m)*()zKHHKH3.24000mx,0y目标位置30000,0mmxmym**()zKK（）采用比例导引法攻击目标二．第一次作业1.求解思路选取一个微小的时间步长dt从初始的飞行状态（可看作初始点）的0v=250m/s;0=0o;可以计算出dt时刻后下一点的x,y.从而可以计算出下一点的H和H，带入式*()zKHHKH中求出z，又根据瞬时平衡假设求出攻角，再计算出升力，阻力带入导弹纵向平面的运动方程利用欧拉法求解,v的微分方程组就可以求出第二个点的,v。然后这样迭代把所有点连线得导弹飞行轨迹。我们可以用图表示11111,,,,nnnnnnnvxyHzYfXfv通过此迭代过程我们可以得最后的轨迹2.计算结果：（在matlab环境下编程求解）K=-0.08，K=-0.032；其他参数曲线如下我们可以看到，计算的结果满足要求三．第二次作业1.要求在前面2段方案飞行的基础上，这次我们加上了第三段的比例导引飞行（攻击目标段）即要设计*bddqKdtdt中的比例系bK(这里bK为了与前面一二段的方案飞行的K区别)和**()zKK（）中的K，K三个参数2.计算结果展示及问题分析首先我们选取了==-0.055=-0.782.5bKKK，，（K，K两个参数取变量对结果影响不大，故这里取常值）得到如下结果我们不难发现在舵偏角的变化曲线中在方案飞行与导引飞行衔接处（24000,3000xmym）会发生4度左右的突变，这样的阶跃信号会引起调整时的振动。是我们不希望看到的。而这样的突变靠改变后面的K，K两个参数（即使使它们随位移或时间变化）是很难实现的。这里我们经过matlab环境中反复的实验给出了一种解决方案，我认为这里的根本问题在于方案飞行和导引飞行的舵偏角控制方程不一样所引发的突变。所以我们给出了以下的平滑过渡方案再第三段飞行中给方案飞行的的舵偏角控制方程和导引飞行的舵偏角乘上相应的系数12和，即：**1()()zKHHKHH（方案飞行的控制方程）**2()zKK（）（导引飞行本来的控制方程）修正后的导引的舵偏方程为121+2zzz其中sin(2*pi/(4*6000)*X(i1=1-))+sin(2*pi/(4*6000)*X2=(i))12和趋势如图所示经过以上过度处理后我们可以得到如下结果我们可以清楚地看到这次弹道倾角的变化趋于平滑。但在方案飞行中的一段和二段之间舵偏角因发动机突然施加推力且控制方程相同还尚未找到合适的办法使之平滑过渡。3.精度分析对攻击目标时做精度分析取弹道中后面的几个点分析X=[2.99962.99993.00013.0003]*410mY=[1.74860.6614-0.4246-1.5096]m取(X(2)，Y(2)),(X(3),Y(3))两个点做线化分析可以得到距目标点（30000，0）的误差为0.12m。四．第三次作业1.作业要求在原有的基础上去掉瞬时平衡假设，考虑导弹的动态特性。进一步得到新的弹道，与原有的理想弹道做对比。2.简单思路分析在去掉瞬时平衡假设后，攻角的计算由原来的通过力矩平衡由舵偏角计算z变为=-其中的计算和原来的一样，而可由z得出3.计算结果及简单问题分析经过稍微的调整系数后K=-0.03，K=-0.032dt=0.01s可以看到考虑了动态特性后，导弹在比例导引段与原有的理想弹道相差较大而方案飞行段与理想弹道吻合度很高。分析其他参数我们可以看到在比例导引阶段有三次明显的震荡。4.下面我们讨论步长dt对结果的影响当dt0.024s时弹道发散不能形成完整弹道图。当dt=0.0224sdt=0.0215sdt=0.0212我们可以看到当dt取得越小时弹道拟合的越好（dt0.024s），这也一定程度上说明了欧拉法的选择收敛性，而非无条件收敛，所以最终我们选取dt=0.01s（在与瞬时平衡弹道对比图中dt=0.01s）就能很好的拟合。五．第四次作业1.作业要求在前面的基础上，求出5个动力系数，选取合适的特征点，并分别求出舵偏角对弹道倾角，及高度的传递函数（分别几位,HGG）。2.简要分析与计算结果1）特征点选取我们共选取了4个特征点，主要的依据是画出5个动力系数的图后，在图中选择一些拐点，驻点，最值点。他们所对应的为特征点。下图为5个动力系数的图我们选取如下特征点特征点Xa1a2a3a4a5注19101.2m109.405.47-1.310.0640.0059发动机开始工作230000m298.3414.92-3.580.0950.0133攻击目标时323999m156.327.91-1.880.0720.0079转为导引飞行时46349.7m150.907.55-1.810.0330.0067第一段飞行拐点（a1,a2,a3,a4,a5的单位在上图中有说明）2）传递函数的求解在忽略长周期干扰的影响下我们利用公式2515342532123()=-=()asaasaaaasGzssAsAsA（1）其中114Aaa;2214Aaaa;30A可以求出对的传函再利用公式00(sin)(