数学建模-安全跳伞的研究文档格式.doc
- 文档编号:6873423
- 上传时间:2023-05-07
- 格式:DOC
- 页数:13
- 大小:371.56KB
数学建模-安全跳伞的研究文档格式.doc
《数学建模-安全跳伞的研究文档格式.doc》由会员分享,可在线阅读,更多相关《数学建模-安全跳伞的研究文档格式.doc(13页珍藏版)》请在冰点文库上搜索。
评
阅
人
分
备
注
全国统一编号(由赛区组委会送交全国前编号):
全国评阅编号(由全国组委会评阅前进行编号):
11
安全跳伞的研究
摘要
本文从建立跳伞安全的数学模型开始,从跳伞运动员在下落过程中各个时刻的速度和达到第一收尾速度的时刻出发,分别通过对这两个方面的深入研究从而制定出跳伞运动员打开降落伞的最佳时机,最后再综合考虑这两个主要因素,进一步深入并细化,从而求得最优解。
模块Ⅰ中,我们将焦点锁定运动的独立性上。
我们通过建立数学模型,并利用MATLAB软件编程求得的v—t图中可比较直观地了解到速度的变化特点。
我们可以发现发现跳伞运动员在空气中下落时,由于受到的摩擦力正比于速度的一次方或二次方,故当经过一段时间后,竖直方向所受的力会达到平衡,之后跳伞运动员的速度将通过一个极小值,随后开始增加,逐渐趋于速度,我们称之为第一收尾速度。
跳伞运动员必须等待这个速度极小值以减小开伞时的震动。
开伞后,经过一段时候后,竖直方向所受的力会达到第二次平衡,之后跳伞运动员的速度将通过另一个极小值,随后也会逐渐增加,直到趋于第二收尾速度。
跳伞时应该有足够的高度,以确保张伞后能有充分时间达到第二收尾速度。
最后,我们通过题中所给出的数据和利用公式。
求出和。
但我们通过计算发现,当时,本方程无解,说明此模型不太合理,我们采用模块Ⅱ。
在模块Ⅱ中,我们从运动的迭加性出发,利用跳伞运动员受到与v的平方成正比的阻力和运动的合成与分解建立相关数学模型,我们发现跳伞者在一个方向的运动状态必然影响到另一个方向的运动。
因此,这两个分运动彼此交叉关联,不能独立,即跳伞运动员离机后的动力学方程在x、y轴方向的投影不是模块Ⅰ中的(3)、(4)式,而是本模块的(3)、(4)式。
由此可见运动的合成和分解时普遍的。
无条件的,运动的独立性只有在特定条件下才成立。
所以,所谓“运动的独立性”不是运动本身的特征。
利用MATLAB软件编程求得的a—v图中可比较直观地了解到加速度的变化特点,我们发现加速度逐渐减小到零。
说明在此模块中,当其加速度将为零时,即为运动员开伞的最佳时刻,此外,我们从图中也可以看出:
当初速度大于125m/s时,就不会出现极小值。
我们发现当跳伞运动员速度降为44.72m.s-1时,就可以打开降落伞,以实现安全降落。
最后,我们从本论文研究方向考虑,为优化安全跳伞指出了一些参考性意见,如:
适当增加飞机高度,把握好打开降落伞的时间等。
【关键词】第一收尾速度运动的独立性雷诺数
一、问题提出
安全跳伞事关跳伞运动员的生命安全,因此组织方、跳伞运动员及社会各方面的广泛关注。
安全跳伞的一个核心指标是确定打开伞的时间后,能安全地完成软着陆。
跳伞属于一项惊险、刺激、挑战自我的运动,优秀的跳伞运动员能在空中做各种惊险动作,从而给人以深刻的印象。
除此之外,自从苏联1930年建立世界上第一支正式伞兵部队以来,各国相继建立了自己的伞兵部队。
作为国防中的重要力量,我国对伞兵部队的训练和投入也在加大,伞兵的安全面临着新的挑战。
随着我国经济的发展,跳伞运动逐渐在民间兴起,跳伞的安全也面临着严重的矛盾。
跳伞安全涉及每一个跳伞运动员和伞兵的生命安全以及其引起的一系列社会问题。
由于跳伞在对很多国人还比较陌生,缺乏经验,民间并没有完全掌握其中的规律,且缺乏一些紧急处理方案。
因此,跳伞安全问题将来可能会逐渐出现在人们的视野中。
从跳伞安全的探索出发,通过建立数学模型,就跳伞者打开降落伞的时间的标准进行定量分析,并从中得出明确、有说服力的结论。
二、问题分析
(一)关于跳伞运动员离机后的运动的讨论
跳伞运动员在下落过程中,究竟做什么运动,需要我们进行讨论。
人教版物理必修1的第1章“运动描述”第1节中,为了说明选择不同的参考系观察同一物体的运动结果会不同,绘制了一幅图片。
图片内同为一跳伞运动员跳伞的过程,飞机做匀速直线运动,运动员从飞机上跳下做平抛运动。
但是这幅图片是不妥的,因为实际生活中永远不可能看到这样的景象。
看过跳伞运动的人都应该注意到,人从飞机中跳出后,很快就落到了飞机的后面,而不是像图中画的那样。
以前这类图片是以炸弹为例面的,将炸弹换成人,就从正确走向了错误,因为炸弹都是流线型的,在空中运动阻力较小,炸弹的运动可看成平抛运动。
但由于人的身体构造复杂,加上背上的背包和自身密度较小,所以阻力对运动的影响很大。
这就导致了炸弹和人的运动规律完全不同[1]。
本论文为了说明问题的简便性,我们假设跳伞员动员在降落过程中做平抛运动。
(二)跳伞运动员离机后的下降过程
在高空中从飞机上跳伞后,下降的第一阶段是在阻尼介质中做平抛运动,持续一段时间后,跳伞者打开降落伞完成软着陆。
对于跳伞者来说至关重要的是打开降落伞的时机,显然不能打开太晚,否则着陆时速度过大,会造成伤害乃至致命。
另一方面,由于高空空气稀薄,下落速度迅猛,能大大增强跳伞的乐趣,有时可能由于跳伞者与飞机或其他跳伞者太近,以致无法打开降落伞。
我们需要考虑跳伞运动员的重量,准确估计空气阻力及打开伞后空气阻力的明显变化,着陆速度,讨论打开伞的最佳位置。
因受到与速度的平方成正比的阻力,其竖直方向所受的力会逐渐达到平衡,速度将随时间减小,再通过一个极小值后,逐渐达到稳定的收尾速度。
跳伞者须在该速度极小值到来之时张开降落,以减小伞的震动,我们称之为第一收尾速度[2]。
打开降落伞后,会逐渐达到另一个稳定的速度值,我们称为,第二收尾速度。
我们采用MATLAB编程画出的图形得的v—t和a—v图中可比较直观地了解到速度和加速度的变化特点.
(三)雷诺数
雷诺数就是表征流体流动特性的一个重要参数,雷诺数的大小取决于三个参数,即流体的速度、流束的定型尺寸以及工作状态下的粘度。
实验表明,外部条件几何相似时(几何相似的管子,流体流过几何相似的物体等),若它们的雷诺数相等,则流体流动状态也是几何相似的(流体动力学相似)。
这一相似规律正是流量测量节流装置标准化的基础。
可见,雷诺数确切地反映了流体的流动特性是流量测量中常用的参数。
在本次跳伞安全研究中,我们使用公式,来估算雷诺数的大小。
Error!
Referencesourcenotfound.Error!
Referencesourcenotfound.
三、模型假设
1、假设本次模型结果不受跳伞运动员主观因素影响;
2、假设跳伞运动员下降过程中不受其他外力作用;
3、假设不考虑题中假设数据所带来的影响。
四、定义与符号说明
——雷诺数
——介质密度
——物体垂直运动方向的横截面积
——在工作状态,流体的动力粘度
——常数
——速度极小值
——第一收尾速度
——开伞最佳时刻
——跳伞运动员离开飞机时的初速度
a——加速度
五、模型的建立与求解
(一)基于运动的独立性的跳伞安全标准——模块Ⅰ
1.动力学特征
跳伞员的跳伞过程可视作在阻尼介质空气中的平抛运动,在此过程中,动力学方程为:
公式1-1
为了确定空气阻力与速度的关系,我们先估算雷诺数的大小由:
取=1.25kg.m-1=1.87ⅹ10-5Pa.s=125m.s-1=2m
可得到:
≈1.67ⅹ107
根据这一结果,可将阻力近似地取为:
公式1-2
式中为介质密度,为物体垂直运动方向的横截面积,为阻力系数,与物体的形状大小有关
2.张伞前的速度变化
为简化计算,我们令,我们就可以将
(2)式写为。
我们选取脱离飞机开始跳伞处为坐标原点,竖直向上为y轴正向,沿飞机飞行方向为x轴正向。
将
(1)式在x、y轴上投影得:
近似地取作,则
公式1-3
公式1-4
初始条件:
t=0时,,,将(3)、(4)两式分离变量后积分可得到、随t的变化规律为:
公式1-5
公式1-6
由此可知
公式1-7
以=75kg,=1.20,=1.25kg.m-3,=0.5m2,=125m.s-1,g=10N.kg-1作为一般值,将以上数据代入公式后,并化简得:
用v作为纵坐标,t为横坐标,利用MATLAB作图得到图5.1,如下所示(相关程序代码见附录1):
图表五1速度v随时间t的变化曲线
我们发现,当跳伞员降落后,其速度减小,经过一段时间后,在w时刻达到极小值,随后开始增加,逐渐趋于速度,我们称之为第一收尾速度。
显然,根据曲线所反映的规律,我们可以看出w便是开伞的最佳时期。
3.关于最佳跳伞时刻t
由(7)式可得
故:
当时,应有
即公式1-8
将(5)、(6)式代入(8)式得
公式1-9
便是方程(9)的解,利用利用MATLAB编程求解(相关程序代码见附录2),并且准确地确定及相应的值。
将数据代入公式,并化简得:
通过解方程,我们发此方程无解,说明当初速度时,不会出现极小值,说明此数学模型的建立不太合理,我们通过采用模块Ⅱ来进行计算并求解。
(二)基于运动迭加跳伞安全模型的建立——模块Ⅱ
1.运动的叠加
跳伞者在跳伞过程中,其动力学方程为:
公式2-1
空气阻力的大小,式中k为常数。
选取跳伞者跳离飞机处为坐标原点,竖直向上为y轴正向,飞机飞行方向为x轴正向。
则
(1)式可写为:
公式2-2
(2)式在x、y轴投影为:
公式2-3
公式2-4
(3)式和(4)式中都包含着另一方向的运动参量,也就是说,跳伞者在一个方向的运动状态必然影响到另一个方向的运动。
因此,这两个分运动彼此交叉关联,不能独立。
2.张伞前的a和v的变化规律
由(3)式和(4)式组成的方程组是有交叉项的一阶非齐次微分方程组,数学上没有解析,但由方程组(3)、(4)式可见,知道了跳伞者任一时刻的速度、,就可算出该时刻的加速度、。
又根据加速度的定义[3]:
公式2-5
公式2-6
并做如下近似计算
公式2-7
公式2-8
进而可知,若知道了跳伞者前一时刻的速度、和加速度、,由(7)、(8)式可以算得下一时刻速度、。
这样,只要给出跳伞者离机时的初始速度,就可以计算出跳伞者在任一时刻的速度。
参照模型Ⅰ中的雷诺数,我们令k=0.375,此外,以=75kg,=125m.s-1,g=10N.kg-1作为一般值,将以上数据代入公式
(2)后,并化简得:
公式2-9
用a作为纵坐标,v为横坐标,利用MATLAB作图得到图5.2,如下所示(相关程序代码见附录3)
图表五2加速度a随速度v的变化曲线
3.关于最佳跳伞时刻
我们发现,当跳伞员降落后,其加速度逐渐减小,直至减少为0。
显然,根据曲线所反映的规律,我们可以看出当a=0,速度变为恒定,这时候便是开伞的最佳时期。
当a=0时,。
利用MATLAB编程(相关程序代码见附录4)对公式(9)的计算得。
所以,当跳伞运动员速度降为44.72m.s-1时,既可以打开降落伞,以实现安全降落。
六、模型的评价
模块Ⅰ:
在论证运动独立性上,我们通过建立相关数学模型,并利用MATLAB软件编程作图,我们可以看出,跳伞运动员的速度会经过一个极小值,最后速度趋于。
这时,运动员就可以打开降落伞了。
但“运动的独立性”不是运动本身的特征,我们通过计算发现此模型建立的方程无解,说明当初速度大于一定值时,跳伞运动员在下降过程中就不会出现极小值了。
因此,模块Ⅰ的结果的误差可能会比较大,可能不适于跳伞运动员采用。
在模块Ⅱ:
首先,我们考虑跳伞运动员下降过程中的两个分运动;
接着我们根据动力学方程在x、y轴方向的投影分别计列出出、的计算公式;
最后利用MATLAB软件编程作出a—v的曲线图,我们根据图形就可以看出加速度一直在减小,说明速度也一直在减小,也就是说跳伞运动员下降过程中不出现极小值。
我们通过图形可以看出,当a=0时,跳伞运动员速度达到第一收尾速度。
因此,当速度降为44.72m.s-1,就是跳伞运动员打开降落伞的最佳时刻了。
七、参考文献
[1]徐峥嵘,徐州高等师范学院,跳伞运动员是做平抛运动吗?
—对“运动的描述”中插图的质疑,物理教师,第29卷第四期,3,5,2008年
[2]韩振海,魏瑛源,徐州高等师范学院,浅析跳伞过程中的速度变化规律,物理通报,第1期,46,1998年
[3]田杨萌,周学麒,贾金萍,河北科技大学基础部,也谈“跳伞过程过程中的速度变化规律”,物理通报,第11期,24,1998年
八、附录
附录1:
clear
clc
ezplot('
v.^2-(375/(3+1.875*t)).^2-2000*((1-exp(-2*t*sqrt(0.05)))/(1+exp(-2*t*sqrt(0.05)))).^2'
[0,140,10])
axis([01000140])
arrow([00],[0140])
arrow([00],[1000])
gtext('
v-t曲线'
)
t/s'
v/m.s-1'
w'
saveas(gcf,'
jianmo1.jpg'
附录2:
t=solve('
32*sqrt(5)*(1-exp(-20*t*sqrt(5)))*exp(-20*t*sqrt(5))/(1+exp(-20*t*sqrt(5)))^3-390.625/(0.625*t+1)^3'
t=double(t);
A=t==real(t);
tt=t(A)
附录3:
v=linspace(125,0,1000);
a=(0.375*v.^2-750)/75;
plot(v,a)
axis([0125070])
arrow([400],[1300])
arrow([400],[4070])
a-v曲线'
a/m.s-2'
jianmo3.jpg'
附录4:
v=solve('
v^2-2000'
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数学 建模 安全 跳伞 研究