2022年数模国赛A题思路分析

国赛已经越来越近了，为了让大家更高效的备考，小编今天给大家带来了2022年国赛A题的思路分析，快来看看吧！

原题重现

评阅要点

A题围绕波浪能转换装置输出能量问题，通过波浪作用，使装置中的浮子与振子产生相对运动，驱动阻尼器做功并将其作为能量输出。要求建立浮子与振子的运动模型，并设计最优阻尼系数使得 PTO系统的平均输出功率最大。

问题 1 考虑浮子在波浪中只做垂荡运动，建立在给定波浪激励力下的浮子与振子的振动方程组。振动方程组中应包含直线阻尼器的阻尼力、弹簧力、静水恢复力(或浮力和重力)、附加惯性力和兴波阻尼力等。利用附件3和附件4的参数值，计算并给出线性和非线性阻尼系数下，在给定波浪激励力下各时间点的浮子与振子的垂荡位移和速度。

问题 2 考虑浮子在波浪中只做垂荡运动，针对阻尼系数为线性与非线性两种情况，给出平均输出功率的具体表达式，进而建立使平均输出功率最大的数学模型。利用附件3和附件4的参数值，给出最优阻尼系数及相应的最大平均输出功率。平均功率应该在浮子与振子进入平稳运动状态下进行计算，或选择比较长的时间段进行计算。在线性阻尼系数下，给出振动方程组的解析解、进行稳定性分析的论文应予以鼓励。

问题 3 考虑浮子在波浪中只做垂荡和纵摇运动，建立在给定的波浪激励力和波浪激励力矩下的浮子与振子的振动方程组。除了问题1中的各种作用力外，振动方程组中还应考虑旋转阻尼力矩、扭转弹簧力矩、附加惯性力矩、兴波阻尼力矩和静水恢复力矩等。利用附件3和附件4的参数值，计算并给出直线阻尼器和旋转阻尼器的阻尼系数均为常量时，在给定波浪激励力及力矩下各时间点的浮子与振子的垂荡位移与速度、纵摇角位移与角速度。

问题 4 考虑浮子在波浪中只做垂荡和纵摇运动，给出平均输出功率的具体表达式，进而建立使平均输出功率最大的数学模型。利用附件3和附件4的参数值，给出直线阻尼器和旋转阻尼器的最优阻尼系数及相应的最大平均输出功率。平均功率应该在浮子与振子进入平稳运动状态下进行计算，或选择比较长的时间段进行计算。

问题精简

01 精简问题1

仅考虑垂荡，在波浪激励力为f cos ?t的作用下，要求计算两种情况下浮子和振子在前40个波浪周期内，时间的间隔为0.2s的垂荡位移和速度。

情况1为直线阻尼器的阻尼系数恒为10000Ns/m；

情况2为阻尼系数与浮子和振子的相对速度绝对值的0.5次幂成正比，且比例系数取10000。

02精简问题2

在问题1的基础上，给出阻尼系数范围为 [0,100000]，幂指数区间为 [0,1]，要求阻尼恒定与变化两种情况下，计算最大输出功率及相应的最优阻尼系统。

03 精简问题3

要求考虑垂荡与纵摇，在波浪激励力fcos?t和波浪激励力矩Lcos?t的作用下，要求计算直线阻尼器和旋转阻尼器的阻尼系数分别恒定为 10000 N·s/m和 1000 N·m·s时，浮子与振子在前 40 个波浪周期内，时间间隔为 0.2s 的垂荡位移与速度和纵摇角位移与角速度。

04 精简问题4

在问题3的基础上，给出直线阻尼器和旋转阻尼器的阻尼系数的取值范围均为 [0,100000] ，要求计算最大输出功率及相应的最优阻尼系数。问题分析是很重要的部分，在这里一定要像阅卷老师展现思考的过程，一定要写明想使用的模型和求解方法，同时可以配流程图。

问题分析

01总体分析

本文主要解决的是一种波浪能装置的能量转化效率最优化问题。问题一和问题二要解决的是在该装置只做垂荡运动时浮子与振子的运动状况以及如何确定最优的阻层系数使得平均输出功率最高。而问题三和问题四要解决的是在该设备只做垂荡和纵摇运动时浮子与振子的运动状况以及通过调节阻尼系统实现平均输出功率的最优化。

02具体分析

问题一波浪能装置垂荡运动模型分析

建立浮子与振子的垂荡运动模型需要准确反映浮子和振子振动的物理过程。对此，基于牛顿第二定律对装置进行受力分析，建立浮子和振子垂荡运动的二阶非齐次线性微分方程，并采用四阶龙格库塔法进行数值求解。

问题二波浪能装置最大输出功率设计分析

波浪能装置输出功率优化设计，需要在垂荡运动模型的基础上，求出给定情景下装置的最大输出功率以及对应的最优阻尼系数．通过推导平均功率函数，得到最大输出功率的设计方法，利用龙格库塔法粗略寻优后，采用遗传算法，分别求出两种给定情景下的最大平均输出功率，以及对应的最优阻尼系数。

问题三装置垂荡与纵摇运动模型分析

浮子在水中只做垂直和纵摇运动，中轴和底座铰接，并在转轴处增加了旋转阻尼器与旋转弹簧，其余条件与问题1相同。假定浮子和振子的角位移很小，因此考虑垂荡时，可以忽略纵摇影响。

在考虑纵摇时，由于计算得到的浮子和振子相对垂荡位移很小，可将二者之视为一个整体计算转动，此后再构建以转轴为原点的平动非惯性参考第，分析振子的转动。分别列出两条动力学方程和2条转动方程，利用MATLAB求解。

问题四利用数值积分法求解PTO

需要对阻尼系数在 [0,100000]内直线阻尼器和旋转阻尼器取值，得到浮子同时进行垂荡和纵摇下最大输出功率及相应的最优阻尼系数。类似于问题二，利用数值积分法求解出系统PTO，再对直线阻尼器和旋转阻尼器最优阻尼系数从0-100000进行遍历，计算出不同条件下PTO。通过多次细分目标区间，找出最优参数。

基本假设

1.假设海水是无粘且无旋的;

2.底座、中轴、隔层、PTO的质量与所有的摩擦不计;

3.忽略附件未给出的有关底座、中轴、中轴架、转轴以及其他结构的大小、厚度、高度等;

4.海平面足够大且浮子的运动对海平面高度不产生影响;

5.浮子顶端不会浸没且圆锥部分不露出水面;

6.在解决问题二时，假设模型将在有限时间内趋于稳定，最终所有广义坐标的圆频率与激励力圆频率相等;

7.在解决问题三时，相比于浮子的垂荡运动，振子与浮子的垂荡运动的相对位移可视为小量;

8.在解决问题三时，浮子与振子的纵摇运动的角位移为小量;以上为2022年国赛A题的基本思路，供同学们参考，并根据自身情况建立模型并完成问题的求解。

此外，一篇合格的数模竞赛作品还需要模型的检验与评价，很多同学总是只给出模型的评价未给出模型的检验，这里给出针对该思路的模型检验方向。

模型检验

01方差检验

该题波浪能装置的运动属于简谐激励下线性系统的受迫振动，此类物体振动稳定后的振动周期应与波浪激励的周期接近，故据此进行方差检验，检验模型的正确性。

02稳定性分析

对于不同大小的波浪激励力矩，计算振子相对于浮子的角位移。若波浪激励力矩大小在一定范围内变动，振子相对于浮子位移均满足小量的假设。此外，还可以计算该模型对扭转弹簧刚度的稳定性分析，计算输出功率，检查是否符合题目中的数据区域。

注意事项

①历年国赛A题经常考察与物理结合的数学建模问题，在审题的时候需要源于机理分析，然后再拓宽至具体数学模型的建立。

②国赛A题常考察微积分问题，建议同学们至少掌握线性微分方程和积分求解方法。CORPORATE CULTURE