降落伞选择问题数学建模
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
降落伞的选择问题
组长:瑜
组员:璐
组员:胡潇
摘要
本文讨论并确定了降落伞的最佳选购方案,在满足空投物资重量的前提下,使购买降落伞的费用最小。该问题是一个优化问题,以购买降落伞的费用最小构造目标函数,以救灾物资2000kg,5种不同半径的降落伞的最大载重量为限制条件,进行线性规划,建立优化模型。通过LinDo软件对模型进行求解,最终得出最佳方案为3m的降落伞数量为6个,其他半径的降落伞不予选购,以及最小费用为4793元。
首先,我们需要计算各规格降落伞的价格,可知其价格由伞面费,绳索费,固定使用费三部分构成,以此进行计算。其次,我们需要计算出阻力系数,我们利用了两种方法确定出阻力系数为2.95747;之后,我们要确定不同半径的降落伞的最大载重量,通过之前计算出的速度与时间的关系式,推出速度与质量的关系,再确定质量与速度的关系,从而通过计算得出不同半径降落伞的最大载重量;最后列出目标函数和约束条件,进行线性规划,利用LinDo软件得出最终结果。
总之,我们的模型在理论分析上提出了选择降落伞最优化,为选择合适的降落伞提供了可行的理论依据。
关键字:优化方案、线性规划、微分方程、MATLAB,LINDO
问题重述
为了向灾区空投救灾物资,需要选择不同类型的降落伞。降落伞根据半径不同分为半径为2m、2.5m、3m、3.5m、4m五种型号,降落伞的造价由伞面费用,绳索费用和固定费用三部分组成。每个降落伞用长为1m的16跟绳索连接重物,重物位于球心正下方的球面处,降落伞在下降过程中除了受到重力的影响外,还受到空气的阻力。并且可以认为阻力的大小与降落伞的速度和伞的面积成正比。其阻力系数可由题中给出的数据确定,问题要求在满足空投物资重量的前提下,使购买降落伞的费用最小。(具体数据见附录中表格1,表格2)
问题的提出
为向灾区空投救灾物资共2000kg,需选购一些降落伞,已知空投高度为500m,要求降落伞落地时的速度不能超过20/
m s。降落伞面为半径r的半球面,用每根长为1m的16根绳索连着载重m的物体位于球心正下方球面处,如图1所示。
图1
每个降落伞的价格由三部分组成。伞面费用由伞的半径r决定;绳索费用由绳索总长度及单价4元/米决定;其他费用为200元。
降落伞在降落的过程中受到了空气的阻力,为了确定阻力的大小,用半径3m、载重为300kg的降落伞从500m高度做降落实验,测得各时刻的高度。
确定降落伞的选购方案,即共需多少个,每个伞半径多大,在满足空投要求的条件下,使费用最低。
模型分析
这是一个优化问题,所求目标函数是降落伞的总费用。针对这个问题我们主要分三部分来分析的。
首先,计算各规格降落伞的价格,由已知其价格的三部分组成:伞面费,绳索费,固定使用费。伞面费为题目所给不同半径决定,绳索由长度决定,固定使用费为常数。
其次,我们分析物资在投放过程中的受力情况。忽略了了其他因素影响,有牛顿第二定律分析可得,是物资受到向下的重力和向上的阻力。接下来的问题就是求出阻力系数。
求阻力系数,我们用了两种方法。
第一种:利用牛顿第二定律,得出速度关于时间的表达式,又由对速度的积分,得出高度与速度的关系,再用MATLAB 作出时间与高度的关系图,分析图像作线性回归,利用MATLAB 软件计算出阻力系数;
第二种:求出的速度的表达式,用MATLAB 软件做出速度与时间的关系图,分析可得出阻力系数的大小。
另外,对于确定不同规格的降落伞最大载重量,利用给速度的关系式,逐步推出速度与质量的关系,进而求得最大载重量。
最后,我们写出了目标函数,并且结合约束条件得出了线性规划,利用LINDO 软件得出结果模型分析
符号说明
1i c (1,2,3,4,5;i =):分别表示购买的半径为2,2.5,3,3.5,4;r =的降落伞的价格,
单位(元)。
2i c (1,2,3,4,5;i =):分别表示购买半径为2,2.5,3,3.5,4r =的降落伞的绳索的价
格,单位(元)。
3i c (1,2,3,4,5;i =):分别表示购买的半径为2,2.5,3,3.5,4;r =的伞面面积,单位
(m )。
i c (1,2,3,4,5;i =): 分别表示购买一个半径为2,2.5,3,3.5,4;r =的降落伞的各
自总费用,单位(元)。
(2,2.5,3,3.5,4)r M r =:指的是半径2,2.5,3,3.5,4,r =的最大载重量,单位(kg )。
()v t : 表示t 时刻降落伞的速度单位(/m s )。 s : 表示降落伞的受力面积,单位(2m )。
t : 表示时间,单位(s )。 k : 表示空气阻力系数
m : 货物的质量,单位(kg )。
g : 重力加速度,单位(2/m s )。
模型假设
1.假设2000kg 物资可以任意分割。
2.假设在降落伞下落过程中只受到重力,和一个可以视为非重力因素共同作用下的合力 空气阻力的影响,不考虑横向受力。
3.假设绳索和伞面的质量忽略不计。
4.假设在受力分析过程中,和下落过程中计算高度时,可将物资看做质点。
5.假设降落伞的阻力与速度和面积的成绩成正比,其系数成为空气阻力系数,为常数。
6.假设绳索的价格每米1元,每个降落伞固定费用是常数为200元。
模型建立
由模型分析可知,这是一个优化问题,要建立费用最小的目标函数,和以不同规格的降落伞总载重量大于等于2000kg 作为约束条件,在LINDO 软件中计算出不同规格降落伞的选择个数,得出一个最优方案。 对此问题分三步进行:
第一步:计算各规格单个降落伞的费用
在建立目标函数时,总费用是各规格降落伞的个数和相应的单个降落伞的费
用,所以首先要计算出各规格单个降落伞的费用。
由题目可知,其价格i c 由三部分组成,伞面费1i c ,绳索费2i c ,固定使用费3i c ,而其中伞面费1i c 为题目中所给的不同半径r 决定,绳索费2i c 由绳索长度和单位长度的价格决定,固定费3i c 为常数
由题目中所给表格1(见附录)及计算可得其费用:
表格 3
r
第二步:计算阻力系数
为表述不同规格总载重量大于等于2000kg 这一条件,并求解降落伞速度满足的微分方程,方正中的重要参数——空气阻力系数是未知的,在此我们需要对已知表格2中的数据进行拟合,从而求出空气阻力系数。 对降落伞进行受力分析见 图1