没有合适的资源?快使用搜索试试~ 我知道了~
拟合算法-重物落水后运动过程的动力学分析.pdf
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
5星 · 超过95%的资源 1 下载量 4 浏览量
2022-05-02
06:57:03
上传
评论
收藏 852KB PDF 举报
温馨提示
试读
33页
拟合算法-重物落水后运动过程的动力学分析.pdf
资源推荐
资源详情
资源评论
全
全
国
国
第
第
七
七
届
届
研
研
究
究
生
生
数
数
学
学
建
建
模
模
竞
竞
赛
赛
题 目 重物落水后运动过程的动力学分析
摘 要:
当发生溃坝溃堤时,很难在短时间内将溃口彻底封堵,通过投放重物对尚存的坝体产生一定的保护作
用,可以延缓溃坝过程,为人民群众的撤离争取更多的时间。本文建立了模型用于计算溃坝封堵时重物的
抛投方式,以使重物达到溃坝口,减少无效投放。
本文采用动力学的方法,对物块在流场中的运动进行了受力分析,主要考虑了重力、浮力、拖曳力、
上举力对物体速度和加速度的影响,其中摩擦力、附加阻力和其他流场力都修正在拖曳力和上举力的阻力
系数中。通过动力学微分方程和边界条件分别得到物体静止进入水中和运动进入水中的物体运动轨迹方程。
根据小型试验的四组数据,拟合了各种形状参数下的轨迹方程系数。对于复杂形状,通过工程中常用的体
积等效和截面等效的方法获得当量直径,将物块等效为球形,并且引入了空心率的概念用于表征物块形状
参数对流场中受力的影响。拟合试验数据得到了决定轨迹方程系数的拖曳力系数
D
C
和
L
C
与空心率
关
系。用当量直径和空心率将复杂形状物块简化成球体,使计算方便,使可以推广到各种形状的块石和沙包
中。
本文所述动力学模型,与小型试验结果吻合效果比较好,拟合相关系数接近于 1。当实际溃坝流场进行
重力相似准则转化后,可以通过本模型进行抛投物块的运动轨迹预测,使封堵用的重物落水后能够沉底到、
并保持在预想的位置,尽可能减少无效投放,具有准确、经费省、风险小、时间短、易重复、条件可以改
变等优点。
对于溃口封堵的真实流场,可以采用重力相似准则转化为模型适用的小型流场情况,通过相似的缩尺
比例得到投放重物的水平距离。利用模型计算,得知对于洪水流速分别为
p
4 /v m s
和
p
5 /v m s
,水
深分别为 3m 和4m 时,在离水面 2m 投放重达1t 的重物进行溃口封堵时,应该分别提前0.93086m和2.1534m
投放重物。
关键词 封堵溃口 动力学模型 阻力系数 空心率 小型试验 相似定律
参赛队号 10407002
队员姓名 游森勇、刘源、吴谭
中山大学承办
参赛密码
(由组委会填写)
2
1.问题的提出
我国经常发生洪水,溃坝溃堤进而引发泥石流灾害造成国家和人民生命财产的严重损失。
溃口水流的流量和速度会比较大,在通常情况下很难在短时间之内将溃口彻底封堵,但如果
通过投放重物对尚存的坝体产生一定的保护作用,目前常利用直升飞机投放堵口组件,封住
溃坝缺口。但是投入溃口的重物落水后受到溃口水流的作用会向下游漂移。为了使封堵用的
重物落水后能够沉底到、并保持在预想的位置,尽可能减少无效投放,必须掌握重物落水后
的运动过程,在预定沉底位置的上游一定距离投放达到一定体积和重量的重物。
理论分析和小型试验获取相关数据的方法广泛用于研究溃口封堵中,在此基础上,根据
水力学已经有的方法进行推广,在获得成功并掌握重物在水中运动的规律后才能够最终应用
于实际抢险行动。因此,理论分析并建立数学模型研究不同形状重物在溃坝洪水中的运动规
律是具有经费省、时间短、易重复、易推广的优点,具有重大的实际价值。本文将针对几种
重物形状、四种不同速度的稳定水流、在三种不同的高度多次重复进行的小型试验进行数学
建模分析,分析重物在水中的运动规律。
2. 模型假设及主要变量符号说明
2.1 假设
(1)水流流场等效于无界的均匀剪切流场,重物落水后,不影响流场的整体特性
(2)理想无粘性流体理论适用于小型试验分析
(3)作用于重物上的水流速度,用垂线平均流速代替
(4)真实溃坝流场与模型水槽的均匀流场满足运动相似
2.2 符号说明
D
F
——为水流流速方向拖曳力;
D
f
——为水深方向上举力;
P
——为有效重力;
h
A
——为水平截面等效面积;
P
A
——为竖直方向截面等效面积;
0
u
——为流场垂线平均流速;
u
——为物体的水流方向速度;
w
——为物体水深方向下沉速度;
——为空心率;
——为边界层厚度;
——为水的密度;
g
——为物体的密度;
c
u
——为重物入水时的水平方向初速度;
0
——为物体入水时的垂直方向初速度。
3
3. 模型的建立
在运动的水中抛投重物,重物在下沉过程中受水流的阻力与冲击力的影响,在下落的同
时将出现漂移。根据流体力学知识,水流拖曳力的影响因素主要有:流体的流速
0
u ,流体的
密度,流体的动力粘滞系数
,边界层厚度
,水深
H
,物体所处深度
h
,物体几何尺寸包
括宽度
a
、高度
b
、顺水流方向上长度
L
等。对于空心的物体还得考虑空心率
,对于有一
定角度倾斜的物体,还得考虑倾斜角
。由流体动力力学的知识可以得到,物体在水流中的
受力包括三类
[1][2]
:
(a). 第一类是与流体与块体的相对运动无关的力,有惯性力、重力和压差力等;
(b). 第二类是依赖于流体与物体间相对运动的力,其方向沿着相对运动方向,有拖曳力,
附加质量力,Basset 力,摩擦力等;
(c). 第三类是依赖于流体与物体间相对运动的力,其方向垂直于相对运动方向,有上举
力、Magnus 力和 Saffman 力等。
图 1 物体质心受力图
对于本文的试验和真实溃坝堵截时的流场特性,溃坝时雷诺数
Re
比较大,暂不考虑
Basset 力,Magnus 力和 Saffman 力的影响,在不考虑稳定性的情况下,忽略造成物体翻转的
摩擦力影响。主要考虑压差力(浮力)、重力和水流方向拖曳力
D
F
,水深方向上的上举力
D
f
的影响,物体在水中的受力图见图 1。
图中物体在水中所受有效重力为
P
(指有效重力,重力和浮力的合力)、水深方向上举力
为
D
f
和水流方向拖曳力为
D
F
。根据流体力学,块体在粘性流体中运动或者流体绕过静止块
体时,块体表面将产生表面力
F
。一般来说,作用在单位表面积上的微分力
Fd
可分解为垂
直和平行于该单元表面的压应力和切应力。合力
F
也可以分解成平行和垂直于运动方向上的
分力,平行于运动方向上的力为拖曳力
D
F
,垂直于运动方向上的分力为上举力
D
f
,也就是
说拖曳力和上举力本质上是切应力和压应力的影响。在研究空心块体的实验中块体垂直水流
方向,也就是说拖曳力是由前后两面的压强阻力与四周的摩擦阻力构成即
1 2
0
cos sin
D P f
A A
F F F p dA dA
更多数学建模资料请关注微店店铺“数学建模学习交流”
https://k.weidian.com/RHO6PSpA
4
其中
1
A
包括块体前后面积,
2
A
包括块体四周侧面积,
对
2
A
值为
90
,对
1
A
迎水流面
为
180
,背水流面为
0
D
F
为上述各种因素的函数,工程上一般根据量纲分析的方法结合实际试验值予以给
D
( , , , , , , , , , , )
F f U H h L a b
(1)
根据 Evett 等人总结有计算公式
[3][4]
2
D D P 0
1
C A )
2
F u u
(
(2)
而上举力
2
D h
1
A
2
L
f C
(3)
对于球形物体
2
D D 0
1
C )
8
F u u
2
d(
(4)
2
D L
1
C
8
f w
2
d
(5)
有效重力
gVmgP
,其中
Vm
g
(6)
上述公式中,
D
C
为流动方向阻力系数,即拖曳力系数;
L
C
为下沉方向阻力系数,即上
举力系数;
m
为物体的质量;
g
重力加速度;
为水的密度;
g
为物体的密度;
0
u 为作用
于物体处,流场的水流速度;
u
为物体的水流方向速度;
w
为物体水深方向下沉速度;
h
A
为
水平截面等效面积;
P
A
为竖直方向截面等效面积;
V
为物体的体积。
根据动力学原理,对物体列出力的平衡关系式有:
x
向:
D
du
F m
dt
(7)
h
向:
D
d
P f m
dt
(8)
将
D
F
、
D
f
、
P
、
m
的表达式(1)、(2)、(3)、(4)、(5)、(6)代入(7)、(8),有
2
0
2
2
2
( ) 2( )
=
2 2
P D
g
g g
h L h L
g g g h L
A Cdu
u u
dt V
gV
A C A C
d
dt V V A C
(9)
令
1
2
P D
g
A C
V
其量纲为
1 m
,
2
2
2( )
g
h L
gV
A C
其量纲为
m s
,
3
2
h L
g
A C
V
其量纲为
1 m
,
则有方程:
2
1 0
2 2
3 2
du
u u
dt
d
dt
(10)
5
对时间
t
积分,得
1 1
0
2
2 3
2 2
1
1
ln
2
t
u u
t
(11)
为了求解该微分方程,必须先确定边界条件,边界条件按照重物落水的初始状态来考虑。
现在将考虑物体落水的状态分为两类
3.1 当物体从运动状态中落水
考虑运动状态入水时的边界条件,当
0t
时,
c
u u
、
0
,
c
u
为重物入水时的水平
方向初速度,
0
为物体入水时的垂直方向初速度,由式(11)有
1
0
1
c
u u
,
2 0
2
2 2 0
1
ln
2
将
1
、
2
代入式(11)有
1 0 0
1 0
( )
( ) 1
c c
c
u u u t u
u
u u t
2 2 3 0 2 3
2
2 2 3 0 2 3
( ) ( )
( ) ( )
th t ch t
ch t sh t
(12)
其中
( ) 1
( ) , ( ) , ( ) , ( )
2 2 ( ) ( )
x x x x
e e e e sh x
ch x sh x th x cth x
ch x th x
令
2
0
sh
th
ch
,则式(12)可以简化为
2
2 3
( )
th t
(13)
将式(14)对时间
t
进行积分,且当
0t
时,有
0
x
、
0
h
,即
1 0 0
1 0
2
2 3
( )
( ) 1
(
c c
c
u u u t u
dx
dt u u t
dh
dt th t
)
(14)
可得方程组
:
0 1 0
1
2 3
3
1
ln( ( ) 1)
( )
1
ln
c
x u t u u t
sh t
h
sh
(15)
剩余32页未读,继续阅读
资源评论
- 2201_757616172024-05-23简直是宝藏资源,实用价值很高,支持!
普通网友
- 粉丝: 12w+
- 资源: 9335
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功