弹道计算大作业(11页).doc
《弹道计算大作业(11页).doc》由会员分享,可在线阅读,更多相关《弹道计算大作业(11页).doc(11页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、-弹道计算大作业-第 9 页弹道计算大作业目录一、初始条件和要求21.1 初始条件21.2 仿真要求2二、模型的建立32.1 升力和阻力模型32.2 大气和重力加速度模型32.3 无控飞行42.4 平衡滑翔42.5 最大升阻比滑翔飞行弹道5三、仿真结果63.1 无控飞行弹道仿真63.2 平衡滑翔弹道仿真73.3 最大升阻比滑翔弹道仿真8附录9一、初始条件和要求1.1 初始条件已知给定的初始条件如下:表1 初始条件名称符号给定值单位参考面积m2展弦比/效率因子/质量115kg重力加速度通过模型计算可得m/s2零升阻力系数/rad密度通过模型计算可得kg/m31.2 仿真要求请使用Simulink
2、或Buildfly完成以下仿真任务:(1)请完成该导弹的无控飞行弹道仿真;(2)请完成该导弹的平衡滑翔方案飞行弹道仿真;(3)请完成该导弹的最大升阻比滑翔飞行弹道仿真;二、模型的建立2.1 升力和阻力模型已知展弦比为的飞行器的升力线斜率为: (1)根据飞行力学相关知识,飞行器的升力系数和阻力系数为: (2)其中,升力线斜率由(1)式可得;为效率系数:。由升力系数和阻力系数,得到导弹的升力和阻力为: (3)2.2 大气和重力加速度模型在计算过程中,大气密度采用如下模型: (4)其中,为海平面的大气密度;。重力加速度采用如下模型: (5)其中,为地球半径;为飞行器距离地面的高度。 无控飞行假设导弹
3、的运动始终在铅垂平面,根据飞行力学知识,得到导弹无控飞行时的运动学和动力学方程为: (6)在上述模型中,假设俯仰角为0。2.4 平衡滑翔所谓的“平衡”可以理解为垂直于速度方向受力平衡,即。因此得到平衡滑翔时的导弹运动学和动力学方程: (7)由于弹道倾角的变化率为常数,方程组中的第二个方程等于0。这个方程可以用来求攻角。2.5 最大升阻比滑翔飞行弹道联立(1)式、(2)式可得升阻比的表达式为: (8)从上式可以看出,由于展弦比、零升阻力系数为常数,因此升阻比只和攻角有关,是关于攻角的函数。因此要使升阻比达到最大,须使得到 (9)因此,以最大升阻比滑翔时导弹运动学和动力学方程为: (10)三、仿真
4、结果3.1 无控飞行弹道仿真根据无控弹道模型,写出s函数,搭建的仿真模块如下图所示:图1 无控飞行仿真模块由于初始条件给定,因此模块没有输入;输出有六个,分别为导弹的射程变化、高度变化、速度变化、弹道倾角变化、攻角变化以及密度变化。模块的仿真时间由高度变化决定,当高度降为0(导弹落到地面上)时仿真结束。导出数据后画图如下:图2 无控飞行时各参数变化3.2 平衡滑翔弹道仿真平衡滑翔弹道仿真模块如下图所示:图3 平衡滑翔模块取仿真时间为150s,无输入,输出分别为:导弹的射程变化、高度变化、速度变化、弹道倾角变化、攻角变化以及密度变化。得到各参量时间变化图如下:图4 平衡滑翔飞行时各参数变化3.3
5、 最大升阻比滑翔弹道仿真按最大升阻比飞行时弹道仿真模块如下图所示:图5 最大升阻比飞行模块取仿真时间为180s,无输入,输出分别为:导弹的射程变化、高度变化、速度变化、弹道倾角变化、攻角变化以及密度变化。得到各参量时间变化图如下:图4 最大升阻比飞行时各参数变化附录附表1 无控弹道飞行时完整的s函数无控弹道function sys,x0,str,ts,simStateCompliance = trace2(t,x,u,flag)switch flag, case 0, sys,x0,str,ts,simStateCompliance=mdlInitializeSizes; case 1, sy
6、s=mdlDerivatives(t,x,u); case 2, sys=mdlUpdate(t,x,u); case 3, sys=mdlOutputs(t,x,u); case 4, sys=mdlGetTimeOfNextVarHit(t,x,u); case 9,sys=mdlTerminate(t,x,u); otherwise DAStudio.error(Simulink:blocks:unhandledFlag, num2str(flag);endfunction sys,x0,str,ts,simStateCompliance=mdlInitializeSizessizes
7、= simsizes;sizes.NumContStates = 4;sizes.NumDiscStates = 0;sizes.NumOutputs = 5;sizes.NumInputs = 0;sizes.DirFeedthrough = 0;sizes.NumSampleTimes = 1;sys = simsizes(sizes);x0 = 0; 2000; 100; -5/180*pi;str = ;ts = 0 0;simStateCompliance = UnknownSimState;function sys=mdlDerivatives(t,x,u) S=1.7; %参考面
8、积,m2 AR=0.86; %展弦比 e=0.9; %效率因子; m=115; %质量,kg g0=9.8; %海平面重力加速度, m/s2 Rd=6371000; %地球半径 r=Rd/(Rd+x(2); g=g0*r2; %飞行器所在高度的重力加速度 rho0=1.225; %海平面大气密度,kg/m3 T0=288.15; rho=rho0*(1-0.0065*x(2)/T0)4.2288;%飞行器所在高度的大气密度 alpha=-x(4); %无控飞行时 CLa=3.141592* AR/(1+sqrt(1+(AR/2)2); %升力线斜率, /rad CDo=0.02; % 零升阻力
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 弹道 计算 作业 11
限制150内