没有合适的资源?快使用搜索试试~ 我知道了~
因子分析-探讨重物落水后运动过程的数学模型.pdf
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 7 浏览量
2022-05-01
20:03:03
上传
评论
收藏 555KB PDF 举报
温馨提示
试读
24页
因子分析-探讨重物落水后运动过程的数学模型.pdf
资源推荐
资源详情
资源评论
2
1 问题重述
我 国 经 常 发 生 洪 水 、溃坝溃堤,进而引发泥石流灾害对国家和人民生命财
产 造 成 严 重 的 损 失 。因 此 ,在 汛 期 河 道 堤 防 一 旦 出 现 决 口 ,必 须 迅 速 堵 口 以减
轻 洪 水 破 坝 造 成 的 危 害 。其中 ,利用 直 升 飞 机 投 放 堵 口 组 件 是 开 展 封 堵 溃 口 的
有 效 方 法 。 由 于 投 入 溃 口 的 重 物 在 落 水 后 会 受 到 溃 口 水 流 的 作 用 而 向 下 游 偏
移 ,因此 我 们 必 须 掌 握 重 物 落 水 后 的 运 动 过 程 ,在预 定 沉 底 位 置 的 上 游 一 定 距
离 投 放 达 到 一 定 体 积 和 重 量 的 重 物 ,使其 落 水 后 能 够 沉 底 、并保 持 在 预 想 的 位
置 , 以 减 少 无 效 投 放 。
由 于 溃 坝 溃 堤 的 高 度 危 害 性 、不可 重 复 性 和 经 济 损 失 过 大 ,肯 定 无 法 通 过
相 关 实 物 试 验 去 研 究 封 堵 用 重 物 落 在 溃 口 后 的 运 动 过 程 ,而只 能 先 通 过 理 论 分
析 和 小 型 试 验 获 取 相 关 数 据 的 方 法 进 行 研 究 。由 于 具 体 情 况 不 同( 如 溃 口 的 纵 、
横 断 面 形 状 ;溃 口 底 面 的 坡 度 、粗 糙 度 ;流速 等 )都 会 影 响 封 堵 重 物 落 在 溃 口
后 的 运 动 过 程 和 沉 底 后 的 状 况 ,会使 问 题 复 杂 化 。因 此 ,本 文 只 在 简 单 情 况 下 ,
通 过 固 定 大 多 数 因 素 、条 件 ,采 用 单 因 子 分 析 法 ,分 别 考 虑 在 不 同 的 影 响 因素
下 物 体 落 入 水 中 的 运 动 状 态 模 型 。
本 次 试 验 选 取 几 种 重 物 形 状 、四种 不 同 速 度 的 稳 定 水 流 、在 三 种 不 同 的 高
度 多 次 重 复 进 行 投 放 试 验 。 通 过 试 验 试 图 解 决 :
(1)分析 影 响 重 物 在 水 中 运 动 过 程 的 因 素 ,并建 立 大 实 心 方 砖 落 水 后 运 动
过 程 的 数 学 模 型 。
(2) 在 试 验 基 础 上 , 建 立 包 含 上 述 各 种 因 素 , 从 而 能 够 适 应 不 同 情 况 的、
描 述 重 物 水 中 运 动 过 程 的 数 学 模 型 。
(3)对所 建 立 的 数 学 模 型 进 行 误 差 分 析 ,并利 用 已 有 的 资 料 验 证 所 建 立 的
数 学 模 型 的 合 理 性,并 根 据 数 学 模 型 提 出 让 堵 口 重 物 恰 好 落 在 最 有 效 位 置 触 底
的 猜 想 和 应 该 进 行 哪 些 新 试 验 来 验 证 这 个 猜 想 。
(4)根据 水 力 学 理 论 模 型 ,选择 对 过 程 影 响 较 大 的 相 似 准 则 对 本 试 验 及 所
建 立 的 新 模 型 的 成 果 加 以 推 广 , 并 对 未 来 的 试 验 和 研 究 工 作 提 出 建 议 。
2 问题分析
2.1 问题 分 析
对 于 问 题 一 ,要建 立 大 实 心 方 砖 落 水 后 的 运 动 方 程 ,需 分 析 大 实 心 方 砖 在
水 中 的 受 力 情 况 ,根据 牛 顿 第 二 定 律 建 立 微 分 方 程 表 达 式 ,可分 别 求 出 大 实 心
方 砖 在 水 平 方 向 和 竖 直 方 向 上 位 移 随 时 间 的 变 化 关 系 ,也 就得 到 了 大 实 心 方 砖
落 水 后 对 时 间
t
的 参 数 方 程 。
3
对 于 问 题 二 ,问题 一 的 分 析 已 考 虑 不 同 重 物 放 置 方 式 、不 同 水 流 速 度 和 不
同 高 度 下 对 重 物 落 水 后 的 运 动 方 程 的 影 响 ,但 依 据 题 意 本 题在 问 题 一 的 基 础 上
还 应 考 虑 重 物 形 状 参 数 的 影 响 ,从 而 建 立 出 适 应 不 同 形 状 、不 同 重 物 放 置 方 式 、
不 同 水 流 速 度 和 不 同 高 度 下 重 物 落 水 的 运 动 方 程 。
对 于 问 题 三 ,对所 建 立 的 模 型 进 行 误 差 分 析 ,首 先 对 小 型 试 验 数 据 利 用 最
小 二 乘 原 理 求 出 模 拟 曲 线 ,分 别 计 算 试 验 数 据 和 模 拟 曲 线 、理 论 曲 线 的 均 方 差 ,
分 析 理 论 模 型 的 合 理 性 。要求 堵 口 重 物 在 最 有 效 的 位 置 触 底 ,即考 虑 重 物 在 水
平 方 向
x
上 的 位 移 , 同 时 还 需 保 证 重 物 触 底 时 的 稳 定 性 。
对 于 问 题 四 ,依 据 流 体 力 学 知 识 ,选 取 对 过 程 影 响 较 大 的 相 似 准 则 ,将 小
型 试 验 推 广 到 真 实 情 况 。利用 选 取 的 相 似 准 则 ,将实 际 数 据 反 馈 到 所 建 的 模 型
中 , 求 得 重 物 沉 底 时 在 水 平 方 向 上 的 移 动 距 离 。
2.2 数据 分 析 及 校 正
由 于 数 据 来 源 于 摄 像 机 隔 着 玻 璃 的 录 像 截 图( 如 图 1),将 不 可 避 免 地 引入
误 差 。 因 此 , 首 先 需 要 对 数 据 进 行 预 处 理 。
消
能
池
转
页
式
尾
门
设置流速仪
试验区域
摄像机
40m
15m 5m
0.5m
电脑
水流方向
1.2m
图 1 试 验 平 面 布 置 图
依 题 意 可 知 ,附 件 数 据 中 的 NaN 表示某 时 刻 重 物 已 经 触 底 或 实 验 录 像 模 糊
无 法 读 取 的 数 据 。显 然 ,对 于 已 经 触 底 的 情 况 ,可 直 接 做 舍 弃 处 理 ;而 对 由 于
录 像 模 糊 无 法 读 取 的 数 据 ,则 需 要 通 过 一 定 插 值 方 法( 如 多 项 式 拟 合)对 其 进
行 插 值 。
附 件 中 所 给 的 数 据 为 显 示 在 玻 璃 网 格 上 的 重 物 的 重 心 纵 横 坐 标 值 ,其 值并
不 完 全 代 表 重 物 的 真 实 位 置 。 由 于 重 物 在 水 中 的 运 动 轨 迹 与 玻 璃 具 有 一 定 距
离 ,使得 所 读 出 的 数 据 与 重 物 运 动 的 真 实 位 置 存 在 一 定 的 偏 差 。根据 折 射 定 理
有
4
sin*tan*20'
cos*tan*20'
'25
20'
tan
120/)'20()'25(tan
sin/sin
22
ryy
rxx
x
y
yxi
nri
(1)
式 中
i
为 入 射 角 ,
r
为 折 射 角 ,入 射 角 与 折 射 角 的 关 系 如 图 2,
为 入 射 线 在
xy
平 面 上 的 投 影 与
x
轴 的 夹 角 。在 等 式
sin*tan*20' ryy
中 ,当
20
y
中 ,取 正
号 ,否则 取 负 号 。
n
为 水 相 对 于 空 气 的 折 射 率(
n
=1.33),
'
x
为 实 验 数 据 ,
x
为
修 正 后 的 数 据 。修 正 后 的 数 据 会 出 现 负 数 ,显 然 这 是 不 合 理 的 ,判为 误 读 ,需
要 采 取 一 定 的 方 法 再 对 其 进 行 校 正 。
图 2 坐标 轴 的 建立
图 3 数 据校 正 原 理 图
以 图 2 所示 , 取 实 验 区 域 的 左 下 角 为 坐 标 原 点 O, 水 流 方 向 为
x
坐 标 轴 ,
垂 直 于 河 面 的 方 向 为 y 坐标轴建立坐标 系(若 无 特 殊 说 明 ,下文均采用此坐标
系 )。 A( 0,27.5,0) 为 重 物 重 心 在 水 面 投 放 时 的 初 始 位 置 ,
x
为 理 想 情 形 下 在
玻 璃 网 格 上 读 取 的 A 点 坐 标 值 。 经 过 matlab 迭 代 计 算 , A 在 无 误 读 情 形 的 在
玻 璃 网 格 上 的 数 值 为( 2.770,26.669,0)。一般说来,当 入 射 角
i
愈 大 ,所 读 数 据
与 重 物 真 实 位 置 的 误 差 越 大 。定 义 平均 角 度 误 差
_ ( ) /
c b
delta avg e e i
,其 中
c
e
为 重 物 在 A 点 时 所 读 的 实 验 值 与 真 实 值 之 间 的 误 差 ;
b
e
为 重 物 运 动 到 横 坐 标 值
与 摄 像 机 的 位 置 横 坐 标 相 等 时 , 实 验 值 与 真 实 值 的 误 差 , 通 常 可 视 为 零 。 表 1
列 出 大 实 心 方 砖 在 距 水 面 0.12m 以 竖 放 状 态 投 入 水 面(取 速 度
U
=0.55
/m s
)所
更多数学建模资料请关注微店店铺“数学建模学习交流”
https://k.weidian.com/RHO6PSpA
5
测 得 的 原 始 数 据 与 校 正 后 的 数 据 。同 理 ,可 得 其 他 测 量 数 据 的 校 正 值( 未 附其
他 校 正 值 )。
表 1: 试 验 数 据 与 校 正 后 的 数 据 对 比
试 验 数 据
纵 坐 标(cm)
25.5
18.5
14.5
9.5
5.5
横 坐 标(cm)
2.5
3.2
4.5
7
10.5
校 正 后 数
据
横 坐 标
0.0008
0.7529 2.2048 4.9979 8.9021
纵 坐 标
26.189
18.312 13.811 8.186 3.6887
3 问题假设及基本符合说明
3.1 问 题 假 设
1) 假 设重 物 在 水 中 运 动 时 不会 和 水 槽 边 壁 发 生 碰 撞 , 且 重 物 触 底 的 瞬 间
即 将 重 物 认 为 静 止 不 动 , 不 考 虑 重 物 在 底 面 发 生 的 偏 转 和 滑 移 ;
2)假设水流为速度不变的稳定水流,且重物运动过程中不发生偏转;
3) 由 于本 实 验 中 投 掷 高 度 较小 , 假 设 重 物 从 某 一 高 度 下 落 至 水 面 时 无 能
量 损 失 , 即 竖 直 方 向 初 速 度 不 变 ;
4) 假 设重 物 在 水 面 放 置 , 其运 动 过 程 为 一 质 点 的 运 动 , 不 考 虑 形 状 的 变
化 , 也 就 是 不 考 虑 物 体 在 水 表 面 的 运 动 过 程 ;
5)假 设 对 于 蜂 巢 类 和 方 砖 类 试 件 ,平放 时 试 件 有 两 条 边 平 行 于 水 流 方 向 ,
立 方 和 竖 放 时 试 件 有 两 条 边 平 行 于 竖 直 方 向 。
3.2 基 本 符 号 说 明
名 称
符
号
单 位 名 称 符 号
单 位
质 量
m
kg
拖 曳 力 系 数
D
C
密 度
3
/kg m
上 举 力 系 数
L
C
重 物 的 体 积
V
3
m
重 物 的 投 影 面 积
A
2
m
溃 口 平 均 水 流 速 度
U
/m s
重 物 的 高 度
s
m
水 深
H
m
重 物 的 长 度
l
m
重 物 离 地 面 的 高 度
h
m
重 物 的 宽 度
b
m
6
4 问题一求解
4.1 模 型 建 立
影 响 水 流 作 用 力 因 素
[1]
主 要 有:流 体 的 流 速 ,以
U
表示 ;流 体 的 密 度
(常
在 研 究 分 层 流 时 考 虑 );流 体 的 动 力 粘 滞 系 数
k
;边 界 层 厚 度
;水 深
H
,重 物
投 放 的 高 度
h
,
即 初 始 摆 放 位 置 ;重物在水平方向的初始速度
xo
;
重 物 最 低 部
( 底 边 ) 与 壁 面 间 的 距 离
y
; 重 物 宽 度 b、 高 度 s、 顺 水 流 方 向 上 的 长 度
l
,即
重 物 的 形 状 、大 小 的 影 响 。另 外 ,在 考 虑 物 体 在 水 下 翻 转 中 还 常 考 虑 物 体 迎流
面 的 法 线 方 向 与 流 体 方 向 夹 角
, 同 时 如 果 重 物 不 是 实 心 块 体 , 在 研 究 空 心 块
体 的 时 候 需 要 考 虑 空 心 率
。 因 此 , 水 流 作 用 力 可 表 达 为 :
( , , , , , , , , , , , , )
xo
F f U H h y l s b
本 文 不 考 虑 重 物 和 边 界 层 以 及 壁 面 间 的 影 响 ,重 物 在 水 中 发生 的 偏 转 的 现 象 也
不 考 虑 , 再 加 上 流 体 水 的 动 力 粘 滞 系 数 已 知 , 所 以 , 水 流 作 用 力 可 简 化 为
( , , , , , , , , )
xo
F f U H h y l s b
取 单 个 块 体 进 行 研 究 ,块体 在 水 中 的 受 力 如 图 4 所 示 ,大实心方块受到的
作 用 力 有 水 流 拖 曳 力 、有 效 重 度、上 举 力 。即 水 平 方 向 上 的 拖 曳 力
D
F
,竖 直 方
向 的 有 效 重 力
G
有
效
和 上 举 力
L
F
, 各 力 的 计 算 表 达 式 如 下 :
图 4 块 体在 水 中 运 动的 受 力 示 意图
根 据 牛 顿 运 动 定 律 , 块 体 在 水 平 方 向 上 的 拖 曳 力 提 供 其 水 平 流 向 上 加 速
度 , 其 方 程 如 下 :
剩余23页未读,继续阅读
资源评论
普通网友
- 粉丝: 12w+
- 资源: 9335
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功