2016全国大学生数学建模大赛国家一等奖论文A题
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
图 1: 传输节点示意图
锚链与锚的链接处的切线方向与海床的夹角不超过16度,否则锚会被拖动。钢 桶竖直时,水声通讯的效果最佳。为了控制钢桶的倾斜角度,链接处可悬挂重物球。
要求通过数学建模来完成以下任务: (1) 某传输节点选用22.05������的II型锚链和1200������������的重物球,布于深18m、海床平 坦、密度为1.025×103������������/������3 的海域。若海水静止,计算海面风速为12������/������和24������/������ 时钢桶和各节钢管的倾斜角度、锚链形状、浮标的吃水深度和游动区域。 (2) 在问题1的假设下,计算海面风速为36������/������时钢桶和各节钢管的倾斜角度、 锚链形状和浮标的游动区域。调节重物球的质量,使钢桶的倾斜角不超过5度,锚链 在锚点与海床的夹角不超过16度。 (3) 因潮汐等影响,布放海域的水深实测在16������ ∼ 20������之间。布放点的海水速 度可达1.5������/������、风速可达36������/������。考虑风力、水流力和水深,请设计系泊系统,分析 不同情况下钢桶、钢管的倾斜角度、锚链形状、浮标的吃水深度和游动区域。
2
二、问题分析
2.1 问题一的分析 问题一要求计算不同海面风速时钢桶和各节钢管的倾斜角度、锚链形状、浮标
的吃水深度和游动区域。要解决此问题,首先对题目中所给出的传输节点示意图建 立直角坐标系,对传输节点进行先整体后局部的受力分析,得出系泊系统各部分的 递推公式,建立递推-校正的反演模型,得到相关相关方程组。
意义 吃水深度 浮标所受浮力 浮标的质量 浮标的重力 浮标直径 浮标高度 锚所受的海床弹力 锚的重力
风力 钢管所受浮力 钢桶及以上各物的浮力和 锚链的浮力和 钢桶及以上各物的浮力和 锚链的重力和
锚链长度 链环长度 重球的重力
4
五、模型的建立与求解
5.1 问题一模型的建立与求解
5.1.1 问题一模型的建立
将其直接代入问题一的模型中,求得此条件下的吃水深度、各钢管倾斜角度、锚链 形状等。此时,发现钢桶倾斜角和锚与锚链链接处的切线方向与海床的夹角均不符 合题目要求,即水声通讯系统设备的工作效果较差,且锚被拖行,因此无法求解游动 区域。
之后,对于第二个子问题,考虑调整重物球的重量,以使倾斜角度满足题目要 求。首先分析重物球对各个角度的影响,发现重物球越重角度越小。据此设立重物 球重力范围,在满足系统各部分递推关系式和钢桶与竖直方向夹角小于5度的条件 下,将多目标优化模型转化为单目标优化模型。然后遍历所有吃水深度的可取值, 寻找可行的情况,再通过分层递进法求出恰好达到要求角度的重物球重力临界值。
对于问题三,考虑到安全因素,将系统设定为处于较极端环境下,即水深、风 速、海水速度分别设为16������,36������/������,1.5������/������。在风速与水流速度同向的情况下,采 用控制变量法对钢桶倾斜角和锚链末端与锚链接处的切线方向与海床的夹角随风 速、海水深度、海水速度影响的变化规律进行分析。从而可以确定使角度满足题意 的变量临界值。然后,以钢桶倾斜角最小和吃水深度最小为目标函数,以水深16������、 系统各部分递推关系式和钢桶倾斜角小于5度为约束条件,以重物球重量、吃水深 度和锚链长度为决策变量,以风速36������/������、海水速度1.5������/������ 为初始条件,得出不同 锚链类型和长度的最优解。最后,选取4组可行解,并对所得系泊系统进行分析。
������2 = ������������������������������������������2
(7)
由此可推导得
������
������
∑︁
∑︁
������������ + ������������������ − ������������ − ������������������
ห้องสมุดไป่ตู้
������������������������������+1 =
基于反演算法的系泊系统设计探究
摘要
本文探究的是近浅海观测网的系泊系统设计问题。通过平衡条件,对系泊系统 的整体和各局部进行受力分析,得到系统各要素之间的递推公式。由此,本文采用 反演、多目标优化等方法分析和解决本问题。
对于问题一,首先对题目中所给出的传输节点示意图建立直角坐标系,进行 受力分析, 得出系统各部分的递推公式, 从而根据递推-验证的反演算法建立模 型。 再将未知量吃水深度ℎ设为已知量,以题目所给水深18������为约束条件,通过 遍历算法运算出结果。接着对异常数据进行作零处理,分析锚链与海床接触的情 况,推导得出钢管倾角、锚链形状等。最后通过系统是否受力平衡对结果进行检 验。结果显示风速为12������/������时吃水深度为0.6870������,第1-4根钢管的倾斜角(度)依次 为1.120、1.127、1.135、1.142,钢桶倾斜角(度)为1.150,游动区域为半径14.34������的旋 转面,24������/������ 时吃水深度为0.7015������/������,第1-4根钢管的倾斜角(度)依次为4.258、4.285、 4.313、4.341,钢桶倾斜角(度)为4.370,游动区域为半径17.46������的旋转面。
(1) 模型的准备
假设吃水深度已知,并令其为ℎ,浮标质量为������������ ,浮标直径为������������ ,浮标高度为������ ,
海水密度为������。
浮标重力������������
������������ = ������������������
(1)
浮标所受浮力������������
(2) 钢管和钢桶的倾斜角度模型 ∙ 浮标的受力分析
图 2: 受力分析图
其中,������������1 为第一根钢管对浮标的拉力。
⎧
⎨ ������������ = ������������ + ������������1������������������������1
(4)
⎩ ������������ = ������������1������������������������1
最后,对模型作出了客观的评价及推广。
关键词: 反演算法;多目标优化;遍历算法;递推公式;系泊系统。
1
一、问题的重述
1.1 问题的背景 为满足海上资源开发利用的需要,海上作业平台、浮游码头及海上浮桥的应用
日渐增加,这些海上漂浮结构物的海上定位需要锚泊系统的约束来实现。 锚泊系统在在海洋海岸工程中的形式多样,应用十分广泛。其中,近浅海观测
此外,可得到:
⎧
⎪ ⎨
������������������������2
=
������������
+
������������1
− ������������ ������������
−
������������1
(6)
⎪ ⎩
������������1
=
������������1 ������
根据反三角关系,得:
为求解此方程组,将题目所给的水深18������作为约束条件,给吃水深度赋以初始 值,遍历所有可能的吃水深度值,记录对应的水下物体竖直方向的高度总和,取 与18������ 最接近时所对应的吃水深度为最优解,推导得相应的倾斜角、游动区域等。
2.2 问题二的分析 问题二的第一个子问题要求在问题一的基础上,将风速设定为36������/������。因此,先
网的传输节点由浮标系统、系泊系统和水声通讯系统组成,浮筒式单点系泊系的设 计与应用对海上科学发展有重要意义。 1.2 问题的提出
已知某近浅海传输节点(如图1所示),将浮标视作底面直径2为������、高为2������、质量 为1000������������的圆柱体,锚的质量为600������������,钢管共4节,每节长度为1������,直径为50������������, 每节钢管的质量为10������������。水声通讯系统安装在一个长为1������、外径为30������������的密封圆 柱形钢桶内,设备和钢桶总质量为100������������。
其中,
������������ = 0.625(������ − ℎ)������������2 ������������ = ������������������
5
根据以上公式可推导出������1,且其与第一根钢管与水平方向的夹角为两平行线间 的内错角,所以相等。
∙ 钢管和钢桶的受力分析
图 3: 受力分析图
对于问题二,首先考虑第一个子问题,将风速36������/������直接代入问题一的模型中, 得出此条件下的吃水深度为0.723������,各钢管倾斜角度(度)依次为8.960、9.014、9.068 、9.123,钢桶倾斜角(度)为9.179,锚链链接处的切线方向与海床的夹角(度)为18.414, 游动区域半径为18.80������。发现此条件下,水声通讯系统设备的工作效果较差,且锚被 拖行。随后考虑第二个子问题,对重物球重量进行调整,以钢桶倾斜角最小和吃水 深度最小为目标函数,以水深18������、系统各部分递推关系式和钢桶倾斜角小于5度为 约束条件,以重物球重量、吃水深度为决策变量,建立多目标优化模型,再将其转 化为单目标优化。之后,通过分层递进法求出达到要求角度的重物球重力临界值, 求得重物球重力最小为2011������������。最后对可能存在的缺失值通过代入法进行了添补。
������������ ������ ������������ ������������ ������������ ������������ ������������������������ ���������′��� ������������������������ ������ ������′ ������ ������′ ������������ * 其余符号在文中详细说明
2.3 问题三的分析 问题三要求考虑海水速度,并给出不同情况下钢桶、钢管的倾斜角度、锚链形
状、浮标的吃水深度和游动区域。 以风速方向和水流速度方向相同时为例,为得出各变量与角度的关系,考虑采
用控制变量法对风速、海水深度、海水速度的变化分别对钢桶的倾斜角度和锚链末 端与锚链接处的切线方向与海床的夹角的影响进行分析,再进一步确定使角度满足 题意的变量临界值。以水深16������、系统各部分递推关系式和钢桶与竖直方向夹角小 于5°为约束条件,将多目标优化转化为单目标优化。通过调节决策变量中锚链的型 号和长度,设定重物球重力的范围,对所有吃水深度的可取值进行遍历,判断对应 的钢管倾斜角度、海水深度是否小于临界值,小于则取其作为可行情况,否则重新 设定锚链的型号和长度,继续运算。最终设计出几种可行的系泊系统方案,并对其 进行分析。再分析在风速与水流速度不在同一方向时的情况,对该模型进行优化。
������������
=
������������������(︀ ������������ 2
)︀2ℎ
(2)
结合近海风荷载公式,得到风力������������
������������ = 0.625������������2 = 0.625(������ − ℎ)������������2
(3)
3
三、模型假设
1. 假设锚链的浮力忽略不计; 2. 假设浮标不倾斜; 3. 假设海床水平; 4. 假设水流方向与海床平行; 5. 假设重力加速度为10������/������������。
四、符号说明
符号 ℎ ������������ ������������ ������������ ������������ ������
其中,������������2 为第一根钢管对第二根钢管的拉力。 对受力分析图进行进一步地分析,建立如下函数关系式:
⎧ ⎨
������������1
=
������������������
′
=
������������
������(
������������ 2
)2ℎ������
(5)
⎩ ������������1 = ������������1������
锚链与锚的链接处的切线方向与海床的夹角不超过16度,否则锚会被拖动。钢 桶竖直时,水声通讯的效果最佳。为了控制钢桶的倾斜角度,链接处可悬挂重物球。
要求通过数学建模来完成以下任务: (1) 某传输节点选用22.05������的II型锚链和1200������������的重物球,布于深18m、海床平 坦、密度为1.025×103������������/������3 的海域。若海水静止,计算海面风速为12������/������和24������/������ 时钢桶和各节钢管的倾斜角度、锚链形状、浮标的吃水深度和游动区域。 (2) 在问题1的假设下,计算海面风速为36������/������时钢桶和各节钢管的倾斜角度、 锚链形状和浮标的游动区域。调节重物球的质量,使钢桶的倾斜角不超过5度,锚链 在锚点与海床的夹角不超过16度。 (3) 因潮汐等影响,布放海域的水深实测在16������ ∼ 20������之间。布放点的海水速 度可达1.5������/������、风速可达36������/������。考虑风力、水流力和水深,请设计系泊系统,分析 不同情况下钢桶、钢管的倾斜角度、锚链形状、浮标的吃水深度和游动区域。
2
二、问题分析
2.1 问题一的分析 问题一要求计算不同海面风速时钢桶和各节钢管的倾斜角度、锚链形状、浮标
的吃水深度和游动区域。要解决此问题,首先对题目中所给出的传输节点示意图建 立直角坐标系,对传输节点进行先整体后局部的受力分析,得出系泊系统各部分的 递推公式,建立递推-校正的反演模型,得到相关相关方程组。
意义 吃水深度 浮标所受浮力 浮标的质量 浮标的重力 浮标直径 浮标高度 锚所受的海床弹力 锚的重力
风力 钢管所受浮力 钢桶及以上各物的浮力和 锚链的浮力和 钢桶及以上各物的浮力和 锚链的重力和
锚链长度 链环长度 重球的重力
4
五、模型的建立与求解
5.1 问题一模型的建立与求解
5.1.1 问题一模型的建立
将其直接代入问题一的模型中,求得此条件下的吃水深度、各钢管倾斜角度、锚链 形状等。此时,发现钢桶倾斜角和锚与锚链链接处的切线方向与海床的夹角均不符 合题目要求,即水声通讯系统设备的工作效果较差,且锚被拖行,因此无法求解游动 区域。
之后,对于第二个子问题,考虑调整重物球的重量,以使倾斜角度满足题目要 求。首先分析重物球对各个角度的影响,发现重物球越重角度越小。据此设立重物 球重力范围,在满足系统各部分递推关系式和钢桶与竖直方向夹角小于5度的条件 下,将多目标优化模型转化为单目标优化模型。然后遍历所有吃水深度的可取值, 寻找可行的情况,再通过分层递进法求出恰好达到要求角度的重物球重力临界值。
对于问题三,考虑到安全因素,将系统设定为处于较极端环境下,即水深、风 速、海水速度分别设为16������,36������/������,1.5������/������。在风速与水流速度同向的情况下,采 用控制变量法对钢桶倾斜角和锚链末端与锚链接处的切线方向与海床的夹角随风 速、海水深度、海水速度影响的变化规律进行分析。从而可以确定使角度满足题意 的变量临界值。然后,以钢桶倾斜角最小和吃水深度最小为目标函数,以水深16������、 系统各部分递推关系式和钢桶倾斜角小于5度为约束条件,以重物球重量、吃水深 度和锚链长度为决策变量,以风速36������/������、海水速度1.5������/������ 为初始条件,得出不同 锚链类型和长度的最优解。最后,选取4组可行解,并对所得系泊系统进行分析。
������2 = ������������������������������������������2
(7)
由此可推导得
������
������
∑︁
∑︁
������������ + ������������������ − ������������ − ������������������
ห้องสมุดไป่ตู้
������������������������������+1 =
基于反演算法的系泊系统设计探究
摘要
本文探究的是近浅海观测网的系泊系统设计问题。通过平衡条件,对系泊系统 的整体和各局部进行受力分析,得到系统各要素之间的递推公式。由此,本文采用 反演、多目标优化等方法分析和解决本问题。
对于问题一,首先对题目中所给出的传输节点示意图建立直角坐标系,进行 受力分析, 得出系统各部分的递推公式, 从而根据递推-验证的反演算法建立模 型。 再将未知量吃水深度ℎ设为已知量,以题目所给水深18������为约束条件,通过 遍历算法运算出结果。接着对异常数据进行作零处理,分析锚链与海床接触的情 况,推导得出钢管倾角、锚链形状等。最后通过系统是否受力平衡对结果进行检 验。结果显示风速为12������/������时吃水深度为0.6870������,第1-4根钢管的倾斜角(度)依次 为1.120、1.127、1.135、1.142,钢桶倾斜角(度)为1.150,游动区域为半径14.34������的旋 转面,24������/������ 时吃水深度为0.7015������/������,第1-4根钢管的倾斜角(度)依次为4.258、4.285、 4.313、4.341,钢桶倾斜角(度)为4.370,游动区域为半径17.46������的旋转面。
(1) 模型的准备
假设吃水深度已知,并令其为ℎ,浮标质量为������������ ,浮标直径为������������ ,浮标高度为������ ,
海水密度为������。
浮标重力������������
������������ = ������������������
(1)
浮标所受浮力������������
(2) 钢管和钢桶的倾斜角度模型 ∙ 浮标的受力分析
图 2: 受力分析图
其中,������������1 为第一根钢管对浮标的拉力。
⎧
⎨ ������������ = ������������ + ������������1������������������������1
(4)
⎩ ������������ = ������������1������������������������1
最后,对模型作出了客观的评价及推广。
关键词: 反演算法;多目标优化;遍历算法;递推公式;系泊系统。
1
一、问题的重述
1.1 问题的背景 为满足海上资源开发利用的需要,海上作业平台、浮游码头及海上浮桥的应用
日渐增加,这些海上漂浮结构物的海上定位需要锚泊系统的约束来实现。 锚泊系统在在海洋海岸工程中的形式多样,应用十分广泛。其中,近浅海观测
此外,可得到:
⎧
⎪ ⎨
������������������������2
=
������������
+
������������1
− ������������ ������������
−
������������1
(6)
⎪ ⎩
������������1
=
������������1 ������
根据反三角关系,得:
为求解此方程组,将题目所给的水深18������作为约束条件,给吃水深度赋以初始 值,遍历所有可能的吃水深度值,记录对应的水下物体竖直方向的高度总和,取 与18������ 最接近时所对应的吃水深度为最优解,推导得相应的倾斜角、游动区域等。
2.2 问题二的分析 问题二的第一个子问题要求在问题一的基础上,将风速设定为36������/������。因此,先
网的传输节点由浮标系统、系泊系统和水声通讯系统组成,浮筒式单点系泊系的设 计与应用对海上科学发展有重要意义。 1.2 问题的提出
已知某近浅海传输节点(如图1所示),将浮标视作底面直径2为������、高为2������、质量 为1000������������的圆柱体,锚的质量为600������������,钢管共4节,每节长度为1������,直径为50������������, 每节钢管的质量为10������������。水声通讯系统安装在一个长为1������、外径为30������������的密封圆 柱形钢桶内,设备和钢桶总质量为100������������。
其中,
������������ = 0.625(������ − ℎ)������������2 ������������ = ������������������
5
根据以上公式可推导出������1,且其与第一根钢管与水平方向的夹角为两平行线间 的内错角,所以相等。
∙ 钢管和钢桶的受力分析
图 3: 受力分析图
对于问题二,首先考虑第一个子问题,将风速36������/������直接代入问题一的模型中, 得出此条件下的吃水深度为0.723������,各钢管倾斜角度(度)依次为8.960、9.014、9.068 、9.123,钢桶倾斜角(度)为9.179,锚链链接处的切线方向与海床的夹角(度)为18.414, 游动区域半径为18.80������。发现此条件下,水声通讯系统设备的工作效果较差,且锚被 拖行。随后考虑第二个子问题,对重物球重量进行调整,以钢桶倾斜角最小和吃水 深度最小为目标函数,以水深18������、系统各部分递推关系式和钢桶倾斜角小于5度为 约束条件,以重物球重量、吃水深度为决策变量,建立多目标优化模型,再将其转 化为单目标优化。之后,通过分层递进法求出达到要求角度的重物球重力临界值, 求得重物球重力最小为2011������������。最后对可能存在的缺失值通过代入法进行了添补。
������������ ������ ������������ ������������ ������������ ������������ ������������������������ ���������′��� ������������������������ ������ ������′ ������ ������′ ������������ * 其余符号在文中详细说明
2.3 问题三的分析 问题三要求考虑海水速度,并给出不同情况下钢桶、钢管的倾斜角度、锚链形
状、浮标的吃水深度和游动区域。 以风速方向和水流速度方向相同时为例,为得出各变量与角度的关系,考虑采
用控制变量法对风速、海水深度、海水速度的变化分别对钢桶的倾斜角度和锚链末 端与锚链接处的切线方向与海床的夹角的影响进行分析,再进一步确定使角度满足 题意的变量临界值。以水深16������、系统各部分递推关系式和钢桶与竖直方向夹角小 于5°为约束条件,将多目标优化转化为单目标优化。通过调节决策变量中锚链的型 号和长度,设定重物球重力的范围,对所有吃水深度的可取值进行遍历,判断对应 的钢管倾斜角度、海水深度是否小于临界值,小于则取其作为可行情况,否则重新 设定锚链的型号和长度,继续运算。最终设计出几种可行的系泊系统方案,并对其 进行分析。再分析在风速与水流速度不在同一方向时的情况,对该模型进行优化。
������������
=
������������������(︀ ������������ 2
)︀2ℎ
(2)
结合近海风荷载公式,得到风力������������
������������ = 0.625������������2 = 0.625(������ − ℎ)������������2
(3)
3
三、模型假设
1. 假设锚链的浮力忽略不计; 2. 假设浮标不倾斜; 3. 假设海床水平; 4. 假设水流方向与海床平行; 5. 假设重力加速度为10������/������������。
四、符号说明
符号 ℎ ������������ ������������ ������������ ������������ ������
其中,������������2 为第一根钢管对第二根钢管的拉力。 对受力分析图进行进一步地分析,建立如下函数关系式:
⎧ ⎨
������������1
=
������������������
′
=
������������
������(
������������ 2
)2ℎ������
(5)
⎩ ������������1 = ������������1������