降落伞模型

数学建模大赛

论文题目：降落伞在下降过程中安全性问题

姓名1：马颖涛学号********专业：土木工程

姓名2：刘雷学号：********专业：土木工程

姓名3：崔磊学号：********专业：土木工程

2012 年5月3日

目录

一.摘要： (3)

二.问题的提出 (4)

三．问题的分析 (4)

四.建模过程 (5)

1模型假设： (5)

2.定义符号说明： (5)

3.模型建立 (5)

4模型求解： (8)

五.模型的评价与改进 (9)

六.参考文献以及附录代码 (10)

摘要：

“降落伞在下降过程中的安全问题”数学模型是通过研究人体的重力、伞

的空气阻力(与受力面积成正比)、弹性绳的拉力之间的关系，建立人在竖直方向上的运动模型，进而给出运动方程。通过查阅资料我们可得一般人落地速度不得大于5m/s ，空气阻力系数为2.9378，重力加速度9.82/m s 。因此通过数据模拟拟合最终的外出最优值。首先考虑最简单的情况，即不考虑绳子的强度，忽略水平方向的风速影响，忽略绳子和伞衣的重量，把人和伞衣看成整体，运用物理学中力与运动的关系和微分方程给出速度和下落时间的微分关系，用matlab 软件给出解析关系。然后用该软件求出人体质量m 和伞衣面积的对应关系，并用表格表示。使不同的人可以根据自己的体重选择降落伞，也可以统计人的平均体重，确定降落伞的一般尺寸。使人们根据自己的体重可以选择适合自己的降落伞。计算过程中，把伞衣视为半圆柱面，并且设定半圆柱面的长度和直径的关系。伞衣面积234

S d π=。但是，这种情

况只能粗略估计体重与伞衣面积的关系，实际中应考虑绳子的强度，即人和伞衣的运动不同步。由图4可知，十秒之后速度趋向恒定，加速度近似为零。此时绳子拉力最大。

关键词：安全问题 运动方程 拟合 Matlab

一、问题的提出

在我们的生活环境中，任何事物都不可能一帆风顺的，尤其是在高空中，此时降落伞成了必不可少的救生设备。降落伞的使用越来越广泛，大到战争中士兵的空降，小到救灾中物资的空投，降落伞都起到了非常重要的作用，降落伞设计制作的优劣对空降人员安全和物资完整变得尤为重要。因此为避免这种现象，必须找出一个合理的方案，一般，通过改变降落伞的表面积增大与空气的接触面积、改变弹性绳与人之间的夹角来控制，使人下落时速度控制在一定的范围内，以便安全着陆。

二、问题的分析

由题意可知，目的就是为了建立一种模型，解决降落伞伞衣面积，伞衣的弯曲弧度，绳长和绳与伞的夹角，从而使降落伞的安全性达到标准的目的。在我们所研究的问题中，首先建立人体竖直方向的运动模型，求解出降落伞在菜蔬一定的情况下额外呢的下落速度并且求出最优解。

图1 图2 图3

三、建模过程

1模型假设：

①降落伞的各个组成部分都完好无损，没有任何质量问题

②降落过程中，弹性绳的伸缩变量在弹性范围内

③在绳与伞，人或者物资与背带系统等打结处不会松

④忽略伞衣因空气阻力而产生的变形

⑤降落伞与绳子质量忽略不计

⑥空气阻力的阻力系数看做定值与其他因素无关

2 定义符号说明

g ——重力加速度 k ——空气阻力系数 m ——人体质量 s ——伞衣的面积

a ——下落某一时刻人的加速度 H ——初始下落高度

v ——下落某一时刻人体速度 L ——每根弹性绳的长度 T ——每根绳子拉力

D ——每根绳子水平投影长度

3.模型的建立

为了方便对人进行受力分析，我们将人和降落伞看成一个整体，由图2可知整体只受到竖直向下的重力和竖直向上的空气阻力，题上给出空气阻力与伞衣面积成正比，而根据查阅资料我们可以得出空气阻力f ksv =这样的关系式，因此根据牛顿第二定律得出：

F =mg-f

合

由运动学公式：

F =ma 合

由速度v 与时间t 的一次积分建立方程（如下）

dv

a dt

mg f ma f ksv

v

=

-===

得出()ks

t m

mg

mg v t e ks

ks

-

=-+

由位移H 与时间t 的二次积分建立方程

22

d H mg ksv

a m dt -==

根据ＭＡＴＬＡＢ（见附录1）得出

2

2

2222

m H t =kst m

g

mgt m g e

ks k s

k s

-

+-（）

如果假设时间t 为常量的话可以得出速度 V 质量m 函数关系为

()sk t m

gm

mg v m e

ks

ks

-

=-+

速度

v 与伞衣面积s 的函数关系为

()sk t m

gm

mg v s e

ks ks -=-+

因为在摘要中我们已经得知max 5v ≤

根据如下方程（运用MATLAB 见附录2）得出一系列数据以及S-M 图表，取其中部分数据如下表：

()2

22222max 2m H t =H t =600V 34

sk t m

kst m gm

mg v t e ks

ks g mgt m g

e ks k s

k s S d

π-

-=-

++-≤=（）

（）5

d=[1.5,1.7，……6.3,6.5]