黑匣子数学建模题

1、 我仔细阅读了学院第四届大学生数学建模竞赛参赛须知。我完全明白，在竞赛开始后参赛者不能以任何方式（包括电话、电子邮件、网上咨询等）与任何人研究、讨论与赛题有关的问题。我知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我将受到严肃处理。我选择的题号为： A 题 姓 名 xxxxxx 学 院 xxxxxxxxx 学 号 xxxxxxxx 联系方式 xxxxxxxxxx （以下内容参赛同学不需填写，由竞赛组

2、委会统一编码） 评阅编号： 寻找黑匣子摘要黑匣子具有记录各种飞行参数的功能，常用于飞机事故分析使用。所以在数以万计的飞机灾难事故现象中，通过逻辑分析黑匣子的数据，探究事故发生的原因，以避免更多的灾难。所以在事故中寻找黑匣子是很重要的。本题根据已知的飞机在高空时的速度、地理位置、所处的环境等条件，围绕所给的突发事件，提出方案，建立高准度、低复杂的数学模型，为寻找黑匣子提供理论依据。对于问题一，我们分析飞机降落过程中所受到的力的作用。将运动轨迹分解成水平轨迹和竖直轨迹。考虑到有空气气流的影响，即水平轨迹上有摩擦阻力和压差阻力的阻碍作用，竖直轨迹上的重力和空气阻力的作用。建立起微分方程模型，把竖直轨

3、迹和水平轨迹联系在一起。利用mathematic软件的求解功能，分别得到飞机坠落轨迹和黑匣子的落水点。对于问题二，由于不考虑洋流流动对黑匣子沉降过程的影响，即海水是静止的。飞机在沉降过程中，竖直方向上只受到重力和浮力的作用，水平方向上不受到任何力的作用。根据牛顿第二定律，建立起微分方程模型,建立起水平和竖直位移之间的关系，得出当黑匣子沉到海底时水平方向上的通过的位移。最后据图1的黑匣子在水下沉降过程中的海底剖面图中的区域范围，选出正确的区域。对于问题三，考虑洋流流动对黑匣子在水中沉降的影响，由于洋流流速和方向受温度地形等多因素的影响。我们假设洋流的流速与方向恒定，根据牛顿定律，建立微分方程模型

4、，得到黑匣子在水中沉降轨迹，通过软件mathematic可以绘制出轨迹曲线。分别带入竖直方向上的位移1000m、2000m、3000m到方程中，得到落水方位。关键词：黑匣子 空气气流 洋流流速 mathematic 微分方程模型 一问题重述飞机是远距离航行的交通方式之一，其主要特点是速度快，安全性高。据统计，飞机是汽车、火车、轮船等几种交通方式中事故率最低的交通方式，但是飞机一旦发生事故，乘客的生还几率非常小。黑匣子是飞机专用的电子记录设备之一，它能记录各种飞行参数，供事故分析使用。黑匣子记录的参数包括：飞机停止工作或失事坠毁前半小时的语音对话和两小时的飞行高度、速度、航向、爬升率、下降率、加

5、速情况、耗油量、起落架放收、格林尼治时间、飞机系统工作状况和发动机工作参数等。假设有一架飞机在高空中飞行时突然发生事故，此时飞行高度为10000米，飞行速度是800公里/小时，航向东北方向45°，飞机在地面的投影位置为南纬22.0度，东经88.0度。请建立模型求解以下问题：1. 假定飞机在发生事故时突然失去动力，考虑飞机在降落过程中受到空气气流的影响，建立数学模型，描述飞机坠落轨迹并推测黑匣子的落水点。2. 假设黑匣子落水之后，不考虑洋流流动对黑匣子沉降过程的影响，建立模型描述黑匣子在水中沉降过程轨迹。如图1所示，假设黑匣子落水点所对应的海底位置为1，落水时沿着图1中指定的虚线方向沉

6、海，给出黑匣子沉在海底的位置，并指出在图形中的哪个区域范围。图1 黑匣子在水下沉降过程中的海底剖面图3. 考虑洋流流动对黑匣子在水中沉降的影响，建立模型描述在有洋流流动的情况下黑匣子沉降轨迹方程，并求解出黑匣子沉入水下1000m，2000m和3000m时离落水点的方位。二模型假设1.空气不受温度、压强等因素的影响，为理想气体。2.在空中降落过程中和水中沉降过程中，飞机不受到爆炸、摩擦等因素破坏其表面和损坏黑匣子的影响。3.在空中降落过程中和水中沉降过程中，重力加速度g不随海拔高度和海下深度的变化而变化，取9.8N/kg。4.将海水看作是不可压缩的流体，即海水密度不随深度、水温等因素的影响。三符

7、号说明m:每个小球的质量n:小球运动的平均速度n:单位体积内小球的个数v1:飞机在空中降落过程中竖直方向上的分速度s:飞机竖直投影面积t:在空降过程中到起落点的时间h1：飞机空降过程中离起落点的竖直高度h2：飞机空降过程中离起落点的水平距离z：飞机沉降过程中离入水点的竖直高度Cf：流体粘性的阻力系数L：物体沿流体运动方向的长度B:横宽A：绕流体沿流动法线方向投影面积Cp：压差阻力系数v1：南赤道流流速v0：飞机初始水平速度四模型的建立和求解4.1.1 问题一的分析当飞机失去动力后，飞机空降过程为一抛物线，我们将飞机空降过程分解为水平轨迹与竖直轨迹水平轨迹上，飞机有一个初速度，由于空气气流与飞机

8、的速度是相对的，若在水平方向的飞机无初速度，则飞机是静止的，则不会受到空气气流的影响。在题中，我们已知飞机在水平方向上有一个初速度，则飞机在水平方向会受到空气气流的影响。对水平方向上的飞机所受到空气气流来说，分压差阻力Fp,摩擦阻力Ff两种与速度方向相反的作用力，根据绕流物体受力公式，Fp=Cp··u ²·A/2,Ff=Cf··u ²·B·L/2,中，为空气密度。我们知道，随着海拔高度越高，就越小。根据附表A.3，ICAO标准大气压1，海拔高度为0处的绝对压强为10·10 ²pa,空气

9、密度为1.2250/ m ³,而在海拔高度为10km处的绝对压强为26499pa，密度为0.41351kg/ m ³,即压强变化很大，密度变化小，密度随压强的变化可略去不计，则空去为不可压缩流体，即密度为常数，取0-10000m密度的平均值，而其它参数Cf、U、B、L、Cp、A等可通过计算与查找资料相关数据可知2，求得Ff与Fp。根据牛顿定律，建立水平方向微分方程模型。在竖直轨迹上，飞机同样受到空气气流的影响，受到重力与空气阻力的影响。空气阻力是空气分子在飞机上下表面的压力的差值。而空气阻力与运动速度、物体的形状，体积等因素有关，计算相当复杂。我们假设：1. 空气中所含的理

10、想化分子视为质量和大小相等的弹性小球2. 在空气系统的压强不太强，温度不太低的情况下，弹性小球间的相互作用力忽略不计3. 小球与其它物体的碰撞是完全弹性碰撞，其运动遵守牛顿运动定律4. 所有的小球都处于永不停息的无规则的运动之中查阅资料得到飞机质量M，计算出重力，根据牛顿定律，建立竖直方向微分方程模型，将竖直方向与水平方向的微分方程模型相联系，得到飞机坠落轨迹。求出坠落时间，可得到黑匣子的落水点。4.1.2 建立模型和求解竖直方向：空气阻力的计算飞机上表面的阻力计算：沿竖直向下运动的空气分子以相对速率（u-v1）与飞机上表面发生弹性碰撞，根据动量定理，与飞机上表面碰撞后小球的冲量：I1=m（u

11、-v1）-（-m（u-v1）=2m（u-v1）根据弹性碰撞定律，小球与飞机碰撞后，飞机的冲量大小：I2=I1=2m（u-v1）， 方向指向飞机上表面在单位时间内与飞机上表面任一面积元dS相碰撞的小球个数有N=n（u-v1）dS/6所以飞机任一面积元dS在单位时间内受到空气分子的冲量大小为dI3=NI2=n（u-v1）dS/6×2m(u-v1)=nm(u-v1)dS/3整个飞机上表面在单位时间内受到的总冲量为I3=dI3=nm（u-v1）S/3因为I=Fdt 所以飞机上表面的空气阻力为F1=nm（u-v1）S/(3t)同理可得飞机下表面的空气阻力为F2=nm（u+v1）/(3t)所以空

12、降过程中的飞机所受总空气阻力为f=F2-F1， 方向竖直向上化简整理得到f=2nmv1·S/（3t）把nm带入，根据牛顿第二定律，建立竖直方向上的微分方程模型Mg-2v1·S/（3t）=M·d²h1/dt²水平方向上：摩擦阻力和压差阻力的计算摩擦阻力3：Ff=Cf·u²BL2压差阻力4：Fp=Cp·u²A2根据牛顿第二定律，建立水平方向上的微分方程模型Fp-Ff=M·d ²h2dt ²化简整理得到-u ²（Cf·BLCp·A）2= M·

13、d ²h2dt ²综合水平轨迹和竖直轨迹，得到飞机坠落轨迹模型因为查阅得到Cf=0.0022，Cp=0.45水平方程d ²h2dt ²=-51.3化简，得h2=222.2t-25.6t ²竖直方程，已知M=50000kg，=0.82，S=233.83，化简得到h1=4.9t ²-9.8t在运动的过程中，算出竖直方向上飞机经过的时间t=46.3秒而水平方向上飞机经历了4.3秒，水平速度减小至零，此时经过的水平距离为 h2=482.2米入水时的速度为v=443米/秒飞机在高空的运动轨迹为先做抛物线运动，水平速度减小至零后竖直下落，直到进入

14、水中此时黑匣子的落水点离南纬22度、东经88度482.2米处，方向为东北方向45度。4.2.1 问题二的分析题二中不考虑洋流流动对黑匣子沉降过程的影响，水流流速为零，即海水是静止的。当飞机以一个倾斜的速度进入水中时，水平方向不受力的作用，飞机做匀速直线运动。而竖直方向上飞机受到水对飞机的浮力作用、飞机自身的重力以及水对飞机的水平压力的作用。在液体中，水分子朝任何方向都有压力的作用，且压强的分布一般是不均匀的【3】。在飞机沉降过程中任意一个位置，同一深度，不考虑温度的影响，压强为定值。水平面上，左右前后的投影面积对应相等。根据流体平面面积上的压力公式：F=PdA，可知水平面上压力合力为零，即飞机

15、在水平面上不受水分子压力的作用。浮力为物体所受液面向上向下的压力之差，即受到各方向流体静压力的向上合力。浮力与水中的密度密切相关。通常情况下，海水是不可压缩流体，密度为定值。通常海水的密度取1.025×10³kg/m³。查阅资料可知飞机的体积【4】，从而得到浮力，根据牛顿第二定律，建立微分方程模型，将此模型与水平方向的运动方程相结合，即得到黑匣子在水中的沉降过程轨迹。根据已给的黑匣子在水下沉降过程中的海底剖面图得出，在不同水平区域范围内，海底离海面竖直方向上的深度不同。将各个区域的不同深度的数值带入模型中，分别算出水平方向上的位移，判断是否在此区域中，选出正确的区

16、域范围。4.2.2 建立模型和求解浮力的计算：海水为不可压缩流体，即海水密度为常数【5】，一般取值1.025×10³kg/m³。因为浮力F=gV，V为飞机排开液体的体积，g=9.8N/kg。已知M=50000kg根据牛顿力学运动定律，水平方向X=v2·t竖直方向，Mg-gV=M·d²z/dt²整理，v用dz/dt代替化简，z=v3·t+（Mg-gV）t²/2将两者相结合，可得黑匣子在沉降过程中轨迹的模型从问题一中，我们已经知道v2=0，即入水后水平方向上没有速度，所以进入水后，飞机沿竖直方向运动，没有水平

17、位移。由于已知条件中，黑匣子落水点所对应的海底位置为区域I，所以据数据分析结果来看，黑匣子的沉在海底的位置为区域I，且在东北方向45度，离南纬22度，东经88度地区的482.2米的水下4200米左右的位置。4.3.1 问题三的分析洋流又叫海流，是指大洋表层海水常年大规模的沿一定方向进行的较为稳定的流动。洋流是有着流速和方向的海流，随季节而不断发生变化【6】。而南纬22度、东经88度地区为南赤道洋流所流经的，其平均流速为2530厘米/秒，冬季流速最大，约为50厘米/秒【7】。我们假设南赤道洋流不随季节、水温、地形等因素的变化而变化，方向恒为自东向西，平均流速恒为30厘米/秒。竖直方向上，仍旧受重

18、力和浮力的影响，东北45度方向上，飞机做匀速直线运动。正西方向上，飞机做以洋流流速为速度的匀速直线运动。将三者不同方向上的运动方程相结合，建立微分方程模型，表示了有洋流流动的情况下黑匣子沉降轨迹方程。将z=1000m、2000m、3000m带入模型，分别求出对应的数值，得到入水点的方位。4.3.2建立模型和求解考虑洋流流动对黑匣子在水中的影响，即海水以恒为0.3米/秒的速度自东向西流动，水平方向上，水分子对飞机的压力差依旧为零，所以水平方向上不受力的作用。在竖直方向的运动情况，与问题二的情况一样，z=-v3·t-（Mg-gV）t ²/2以入水点为坐标原点建立空间三维坐标系。

19、令x轴为东北方向45度，y轴为西北方向45度，z轴为水下深度。x轴线上的位移：x =（v2-2/2·v1）ty轴线上的位移y=2/2·v1·t联立三者的轨迹关系，建立数学模型来描述在有洋流流动的情况下黑匣子沉降轨迹方程：通过mathematic软件的图形描绘出大致轨迹：当z=1000米，t=0.12秒，x=0.03米，y=0.03米；当z=2000米，t=0.16秒，x=0.034米，y=0.034米；当z=3000米，t=0.20秒，x=0.043米，y=0.043米；所以，黑匣子沉入水底1000米处离落水点西北方向45度0.03米远；沉入水底2000米处离落水点西北方向45度0.034米远；沉入水底3000米处离落水点西北方向45度0.043米远；五模型的分析优化由