飞机减速作业

2013高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：B我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：参赛队员(打印并签名)：1.2.3.指导教师或指导教师组负责人(打印并签名)：（论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上内容请仔细核对，提交后将不再允许做任何修改。如填写错误，论文可能被取消评奖资格。）日期：年月日赛区评阅编号（由赛区组委会评阅前进行编号）：2013高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：摘要本文在合理的假设之下，针对本文中的五个问题建立了合适的模型。在求解方面，我们充分利用计算机软件求得结果。问题一：该问题属于物理学问题。本文运用物理学知识结合数学求解方法建立微分方程模型解决。运用物理基本公式解决题目，在求解过程中运用MATLAB数学软件进行绘图求解，既能精确计算出结果也能通过图像清晰看出飞机减速情况。问题二：在问题一得基础上，同样的减速伞装备9T重型轰炸机上，在减速伞的阻力作用下飞机做减速运动直至停下。已知减速伞的作用阻力与飞机的滑跑速度成正比同问题一可建立微分方程模型，并对模型进行求解。将得到的结果与实际跑道长度1500m作对比，当时跑道的长度可以保障飞机安全着陆，反之则不满足要求。问题四：已知减速伞的阻力与飞机的滑跑速度的平方成正比，因此22vkf，其他条件相对问题一二没有变化，可建立微分方程模型。问题五：在增加考虑跑道摩擦力影响的情况下，只需将跑道阻力产生的加速度加在原微分方程模型中。飞机落地以后,最主要的减速手段是以下四种:1引擎反推；2减速板；3机轮刹车；4减速伞。民航喷气式飞机在接地后主要借助反推装置配合传统的主机轮刹车组件来达到最佳减速效果，而反推的原理是利用导流板将发动机排气的方向发生偏转，从原本向后变成向前倾斜以产生拉力，对于军用战机来说由于其速度高吨位大，而且在野战环境下战机跑道受到极大制约。因此，需要减速伞配合减速板与机轮刹车以缩短所需跑道长度。本文建立微分方程模型，将所给的飞机接触跑道时的速度和滑跑的一段距离以及滑跑这段距离未时刻的速度等作为初始条件进行求解。在求解的过程中我们还作出了一些图像，经分析所建模型基本符合实际，能够很好的描述了在减速伞的作用下飞机（轰炸机）在跑道上滑跑的位移的变化规律以及在飞机降落滑跑的不同时刻能够较准确的预测滑跑的距离。该模型可用于飞行员判断飞机（或轰炸机）在减速伞的作用下是否可以在有限的跑道上安全着陆。关键字：MATLAB图像求解微分方程模型一、问题重述当机场跑道长度不足时，常常使用减速伞作为飞机着陆的减速装置。当飞机接触跑道开始着陆时，由飞机尾部张开一幅减速伞，通过增大气动阻力使飞机减少滑跑的距离，保障飞机在较短的跑道上安全着陆，见图1。（图1）问题1.一架重4.5T的歼击机以每小时600km的航速开始着陆，在减速伞的作用下滑跑500m后速度减为每小时100km。设减速伞的阻力与飞机的滑跑速度成正比，并忽略飞机所受的其他外力，试计算飞机的滑跑距离。问题2.将同样的减速伞装备在9T重型轰炸机上，现已知机场跑道的长为1500m，若飞机着陆速度为每小时700m，问该跑道的长度能否保障飞机安全着陆？问题3.就问题1和2的一般情形分别进行讨论。问题4.若实验表明减速伞的阻力与飞机的滑跑速度的平方成正比，试讨论问题1，2和3，并加以比较。问题5.若飞机除受到减速伞的阻力之外还受到跑道摩擦力的影响，模型将做怎样的修改？二、模型假设模型的基本假设：1.忽略飞机所受的其他外力，并把飞机看成物理学上的质点，即不计飞机长度但有质量；2.令飞机所受空气阻力为f，飞机的反向加速度为a,歼击机的质量为mmg,轰炸机的质量为Mmg，设飞机刚接触跑道的速度为0v/kmh,滑行t后的速度为vt/kmh；3.设1k为减速伞的阻力与飞机的滑跑速度成正比的比例系数，即1fkv.设2k为减速伞的阻力与飞机的滑跑速度的平方成正比的比例系数，即22fkv；4.如果考虑摩擦力时，令重力加速度为g,摩擦系数为受到跑道摩擦力的影响。三、符号说明符号说明：f飞机所受空气阻力a飞机的加速度m歼击机的质量M轰炸机的质量0v飞机刚接触跑道的速度1k减速伞的阻力与飞机的滑跑速度成正比的比例系数2k减速伞的阻力与飞机的滑跑速度的平方成正比的比例系数g重力加速度摩擦系数四、问题分析及模型的建立与求解问题一：由于飞机的质量4.54500mtkg，降落初速度0v10006003kmmhs，滑跑距离为1s500m时速度为1v100kmh,假设减速伞的阻力与飞机的滑跑速度成正比，并忽略飞机所受的其他外力，计算飞机的滑跑距离，本文将该问题理解成较为理想化的物理学问题，利用相应的微分知识解决。由相关物理公式：fmadvadtdsvdt和题意得出以下关系式:对其求解，得到相应微分方程：212dsdsmkdtdt，求解得s。具体求解过程如下：是二阶线性齐次微分方程，通过对微分方程求解，得到它的通解为：根据条件:0)0(s，可得到其方程的一个特解：221)(cectstmk；因此可得微分方程组：(2))()(111222tmktmkemkcdtdstvcects；根据初始条件：0)0()()()0(11110svtvstsvv；得到：110110012)(k,svvmvvvsc；求解模型：tsvvtsvvevtvvvvsevvvsts11010010011001)()(；由初始条件：11011()1=(0)0v(0)v()4.5tfkvsvtvmtss阻212dsdsmkdtdt112()ktmstcce45000)0(9250)(500)(3500)0(11mstvtsv解得：121250,c600k当飞机停止运动时，求得所用时间为：-183ln29.60s5500000000t；利用matlab作图(见附录一),01020304050607080901000100200300400500600t轴s轴飞机滑跑位移图像29.429.4529.529.5529.629.6529.729.75599599.5600600.5601X:29.6Y:599.8t轴s轴飞机滑跑位移图像得飞机滑跑距离为599.9sm。问题二：同样的减速伞装备装置在9T重型轰炸机上，在减速伞的阻力作用下飞机做减速运动直至停下。减速伞的作用阻力与飞机的滑跑速度成正比，因此建立微分方程模型进行求解。将得到的结果与实际跑道长度1500m作对比，当1500sm时跑道的长度可以保障飞机安全着陆，反之则不能。根据牛顿第二定律fma可建立下列微分方程：2120112(0)()dsdsmkdtdtvvvtvmM由问题一可知11250k，29000Mkg，又由问题一种所求得的2c，可知14002c于是有微分方程组ttetvets365365625)(14001400)(；由物理知识可知减速伞的阻力与飞机的滑跑速度成正比，当速度越来越小，阻力就会越来越小，这样飞机就会永远运动下去。结合实际的运动情况我们可以假定当速度smv/000001.0时，飞机停止运动。所以当飞机停止运动时，可求得时间t为0ln62500000536t，由此可求得14001500smm。说明该跑道的长度能保障飞机安全着陆。利用matlab软件作出如下图像（具体见附录二）：0501001500200400600800100012001400t轴s轴飞机滑跑位移图像问题三：0501001502002503000100200300400500600t轴v轴飞机滑跑速度图像050100150200250300050010001500t轴s轴飞机滑跑位移图像根据以上图形可知，一般情况下，飞机滑行距离不超过1500m,滑行时间不超过300s。问题四：已知减速伞的阻力与飞机的滑跑速度的平方成正比，因此22vkf，其他条件相对问题一二没有变化，可建立同问题一二的微分方程模型。根据牛顿第二定律maf可得微分方程：2222011111()(0)0(0)()()dsdsmkdtdtsvvvtvstsmM；微分方程：2222dtdsmkdtsd；求得其微分方程的解为：mcktckkmts)(ln)(22122；则建立微分方程组为：(3)/)/(mc)(/)/)(ln()(22112212kctcdtdstvkmctckmts；其初始条件：0)0()()()0(11110smvtvstsvv；求得：1012101201lnklnc,vvsmvvsvc；则求得其微分方程组为：110010100101)(ln)()1lnln(ln)(1stvvvsvtvsvvvvvsts；根据此题给的初始条件：0)0(45009250)(500)(3500)0(11smtvtsv；可求出（具体过程见附录四）：35001c6ln5002c6ln92k；我们知道此问题中说飞机只受减速伞的作用，即阻力只和速度有关系，而且减速伞的阻力与飞机的滑跑速度的平方成正比，当速度smv/1时，速度越小，则速度的平方更小，阻力也就越来越小，这样飞机就会永远的运动下去，因此，结合实际我们假定当速度smv/1时飞机则停止运动。根据以上求出时间t为：277.38s6ln497t；可求得飞机滑跑距离为：1430sm；对于问题二初始条件也没有变化，由牛顿第二定律建立微分方程：22220112()(0)()dsdsmkdtdtvvvtvmM；微分方程组为:11211121()ln(()/)/()mc/()/stmkctcmkdsvtctckdt；根据其初始条件：900091750)0(0)0(mvs；可得方程组：900)6(ln175175000)()1)6(ln367ln(6ln1000)(tdtdstvtts；求出时间t为：6964ln67t；可求得轰炸机的滑跑距离为：29001500smm；因此，该跑道长度不能保障轰炸机的安全着陆。飞机在不同阻力下的飞行降落速度v以及位移s图像如下：050100020406080100120140160180v1=500/3*exp(-5*t./18)0501000100200300400500600s1=-600*