《降落伞的选择——数模优秀获奖论文》由会员分享，可在线阅读，更多相关《降落伞的选择——数模优秀获奖论文（6页珍藏版）》请在金锄头文库上搜索。
1、1降落伞的选择戈亮峰，李昌富，胡 尚（湖南城市学院）摘要：本论文通过列出降落伞空投时的运动方程，分析其下降过程中的速度变化规律，然后运用线性最小二乘法，利用概率论的正态分布函数求出空气的阻力系数 k 为 2.944 。然后根据降落伞在任意时刻的速度 v 是关于载重质量 m 的严格单调递增函数，以及降落伞在接近地面处达到 20 m/s 的最大允许速度 ，计算出各种不同半径大小的降落伞所对应的最大承载质量以及其空投单价。无条件限制时，求解此题可以归结为一个线性规划问题。再运用Mathematica 4解出：= 0 ； = 0 ； = 6 ； = 0 ； = 0。1X23X45X即需要购买半径为 3
2、 m 的降落伞 6 个，其总费用约为 4832 元，当有条件限制时，则求解此题可以归结为一个整数类型的线性规划问题。运用Mathematica 4解出：= 0 ； = 1 ； = 0 ； = 0 ； = 3 。1X23X45X即需要购买半径为 2.5 m 的降落伞 1 个，半径为 4 m 的降落伞 3 个，总费用为 5282 元，从而使得空投费用最低。1 问题重述（略）2 问题的分析根据题意，每种降落伞的价格是可以计算出来的。要确定降落伞的购买方案，即选购降落伞共需多少个，每个降落伞的半径又应为多大（在表 1 中选择） ，在满足空投要求的条件下，使费用最低。首先，根据表 1 中所给的半径 r
3、求出它的最大载重质量 M (r) ，要计算 M (r) 可运用动力学公式和微分方程解出阻力系数 k ，再求出 M (r) 。然后，从降落伞选购方案中运用线性规划选出最优组合模型，得出问题要求的结果。3 合理假设1、忽略气温、压强对空气密度的影响及风向、高空气流对空投物降落（运动）过程的影响，即降落伞在下降的过程中只受重力和空气阻力的共同的作用。22、假设空气密度均匀，即不随高度的变化而改变。3、降落伞落地时的速度不超过 20 m/s 。4、降落伞和绳索的质量问题可以绝对保证以及重量可以忽略。5、假设在降落伞的整个下降过程中不会出现降落伞不能打开和绳索扭曲打结的现象。6、救灾物资 2000kg
4、可以按照 100kg 或 200kg 任意分配组合。7、空气的阻力系数 k 是定值，与外界因素无关。8、M(r)表示半径为 r 的降落伞在满足空投的条件下最大的载重质量。4 变量声明（1）k ：空气阻力系数。（2）t：降落伞下降时间。（3）m：降落伞的载重质量。 （4）g：重力加速度（ ） 。29.8/gms（5）S：降落伞的表面积。（6）H(t) ：降落伞从降落位置到任意时刻 t 下降的位移。（7）Z：最低空投费用。（8） ：某种类型降落伞的个数。nX（9）L：每根绳子的长度（单位：m） 。5 模型的建立与求解（一） 、求解空气系数 k根据题意：降落伞在下降过程中受到重力和空气阻力的共同作用
5、并且初速度 V（0）= 0 FmgfmgkvSkvSa gm合 阻由简单的动力学公式及微分方程，得 ().(1)(0) .(2)kStmdVtkvSgmgekS解 得 ：V(t)=设降落伞从开始下降位置到 t 时刻所下降的位移 H(t) ，即有30()().(3)tHtVdt积分解得： 22() .(4)kStmmgtggt ekS我们进一步得到：t (s)3 3 3 3 3 3 3 3 3 3H(m)3045 53 55 53 49 55 52 5354表 3由物理学知识： 可知，如果在相同的时间间隔内位移相等，则物体作匀速直线运Svt动，也就是说对于此题根据 图（4）的速度与时间的函数关系
6、可知， 降落伞在 t = 9 s 的以后 H(t)与 t 近似呈线性变化，因此我们可以求出在 9 s 后的平均速度 v。求得 v = 17.667 m/s。又因为在此阶段加速度 a 0 ，我们利用 ，易得mgkS出 k=2.944。根据公式 () kStmmggeVtkS=将 r=3m，m=300kg， ， ，k=2.944 代入上式得29.8/s2r0.546()17.6(1)tt e我们运用 Mathematica 4 做出上式的图形如下：5 10 15 20 25 3015.7516.2516.516.751717.2517.5图（4）4由表 3 的数据显示在 t = 9 s 以后，在相
7、等的时间间隔内 H 趋于定值 53，即空投物近似作匀速直线运动。设 , 其中 服从正态分布。()HtVcX=9 12 15 18 21 24 27 30H=128 183 236 285 340 392 445 499调用 Matlab 命令 ；(,1)polyfitXHp=17.5794 -29.2976;所以，V=17.6667(m/s)，空气阻力系数 k=2.944。（二） 、第二步求解当半径为 r 时，伞在满足空投的条件下的最大载重质量 M(r)由 kStmmgekSV()=可知： V(m) 是关于 m 的增函数（证明见附录一） 。因为 V(m) 的反函数为 m(V)，由原函数与反函数
8、单调性一致的性质可知，m(V)是关于速度 V 的递增函数。因此，在所取半径为 r 的伞在满足空投的条件下最大的载重质量 M(r) 就是在速度 V 取最大值时取得。即 V = 20 m/s ，则 M(r) 由以下方程组求出 2 2() kStmkStggetkSmt mgHt kSS由方程组导出（推导过程见附录二）： 2 2ln(1)/().(5)vvggkS由题中所给 H = 500 m , v = 20 m/s , k = 2.944 代入（5）式，得又 2 222.91320ln(1)/(.9513)050.(6)2.953Sg SmgS ，分别代入表 1 中的半径 r = 2 , r =
9、 2.5 ， r = 3 ， r = 3.5 ， r Sr5= 4然后利用 Mathematica 4 解得半径为 r 的伞在满足空投的条件下最大的载重质量 M(r) （详细过程见附录三）数值如下：M(2) Kg M(2.5) KgM(3) Kg M(3.5) KgM(4) Kg151.000 236.067 340.464 462.154 603.998（取整）151236 340 462 604（三） 、计算每种伞的单价如下：单位 /元半径 r 2 2.5 3 3.5 4C1 65 170 350 660 1000C2 181.0193226.2742271.5290316.7838362
10、.0387C3 200 200 200 200 200C 446 596.3 821.5 1176.8 1562备注： C 2=416L=416 r2（四） 、从上表可得出每种伞的单价与其对应的最大的载重质量 M(r) ，接下来就是一个线性规划问题。设 n=(1，2，3，4，5) 为各种半径的降落伞的个数，则有：nXZ=Min 446 +596.3 +821.5 +1176.8 +1562 13X45Xs.t. 23451234523455606202007.(6)XX， ， ， ， 正 整 数6由 Mathematica 4 解得： 123450XX=0， ， 6， =0，此时可以按照 33
11、34+3342 的方式分配，其最低费用为：4832.35 元。当救灾物资以每袋 100kg 或 200kg 等包装空投(每降落伞可多包捆扎空投，但不可将一包分开)时，另设： n=(1，2，3，4，5) 为各种半径的降落伞的个数，则有：nXZ=Min 446 +596.3 +821.5 +1176.8 +1562 13X45Xs.t. 12345234512345000060207.(6)XX， ， ， ， 正 整 数因为 都属于整数，结合 Mathematica 4 算出的解得出：12345， ， ， ， 3X=0， ， ， =0，此时可以按照 6003+2001 的方式分配，其最低费用为：5282.3 元6 模型检验（略）7 模型优缺点分析（略）8 模型的评价及推广（略）参考文献：1 同济大学数学教研室 高等数学（第四版上、下册） 北京 高等教育出版社 19932 邓明成 新编大学物理学（第二版） 北京 科学出版社 19953 秦裕瑗、秦明复 运筹学简明教程 武汉 高等教育出版社 施普林格出版社 2000.84 盛骤 概率论与数理统计（第三版） 浙江 高等教育出版社 20005 降落伞的选择 http:/ 2004 年75 月 3 日
