1. 引言
近年来,能源危机关系到国家安全,地球非再生能源逐渐枯竭以及环境状况日益恶化的形势下,如何高效利用资源充裕、可再生的海洋资源,已经成为各国和地区的普遍共识。《我国海洋能开发利用发展战略研究的基本思路》指出全球大部分海域的波浪能流密度属于可利用范围 [1]。而我国正处于波浪能落后的阶段,在波浪能领域的发电效率远远不如其他发电领域。但波浪能分布广泛、能量密度高,提示着我国必须重视其发展,而波浪能装置的能量转换效率是波浪能规模化利用的关键问题之一。目前的研究主要集中于实现波浪能发电系统的最大功率捕获问题,从PTO系统反作用力的控制角度可以分为无源阻尼控制与无功控制两种方法,其两种控制理论方法都有限定的适用范围,如波浪能发电装置的类型以及机械装置与波浪固有频率大小的限制 [2]。由王亚龙等人提出的波浪能发电装置结构优化显示波浪能收集装置对振荡水柱式波浪发电具有重要的意义和应用价值 [3]。
2. 波浪能装置数学模型建立
2.1. 垂荡纵摇浮子式波浪能装置介绍
如图1所示,为一种波浪能装置示意图,由浮子、振子、中轴以及能量输出系统(PTO,包括弹簧和阻尼器)构成,其中振子、中轴及PTO装置(如图2所示)被密封在浮子内部;浮子由质量均匀分布的圆柱壳体和圆锥壳体组成;中轴底座固定于隔层的中心位置,中轴架通过转轴铰接于中轴底座中心,中轴绕转轴转动,PTO系统连接振子和转轴架,并处于中轴与转轴所在的平面。除了直线阻尼器,在转轴上还安装了旋转阻尼器和扭转弹簧,直线阻尼器和旋转阻尼器共同做功输出能量。在波浪的作用下,浮子进行摇荡运动,并通过转轴及扭转弹簧和旋转阻尼器带动中轴转动。
浮体在波浪的作用下做摇荡运动时,会受到海水的作用,包括附加惯性力(矩)、兴波阻尼力(矩)和静水恢复力(矩) [4]。设定Fe水平面波动造成的波浪激励力,Fr为水平面兴波阻尼力,FPTO为输出系统PTO力,Fhs为静水回复力,设m1为浮子质量,m2为振子质量,A为附加质量,设z1为浮子垂荡位移,z2为振子垂荡位移。
![](//html.hanspub.org/file/1-2820139x8_hanspub.png?20230328082613854)
Figure 1. Heave and pitch float-type wave energy device
图1. 垂荡纵摇浮子式波浪能装置
![](//html.hanspub.org/file/1-2820139x9_hanspub.png?20230328082613854)
Figure 2. Schematic diagram of PTO device
图2. PTO装置示意图
1) 波浪激励力:海平面波浪对装置产生的周期性作用力
(1)
2) 兴波阻尼力:浮体在海水中做摇荡运动时,会兴起波浪,从而产生对浮体摇荡运动的阻力
(2)
3) PTO作用力:PTO装置中弹簧以及阻尼对装置的作用力
(3)
4) 静水回复力:浮体在海水中做垂荡或者纵摇运动时,会受到使浮体回到平衡位置的作用力
(4)
根据振子、浮子运动方式,可以得出静水中浮子的静止水面下的体积随时间的变化
,由后续求解可知浮子振幅小于1 m,静水恢复力在线性区间(圆柱区间)。因此静水恢复力有以下形式:
(5)
将5)带入4)式中可得,
(6)
根据浮子运动状态可知,此外,由于推动浮体做摇荡运动的力不仅要推动浮体运动,还要推动浮体周围的流体运动,因而要使浮体在海水中获得加速度,需要施加额外的力,即为附加惯性力。对应牛顿方程中惯性质量的增加A。
2.2. 浮子、振子运动模型
对振子和浮子分别做受力分析如图3所示:
![](//html.hanspub.org/file/1-2820139x17_hanspub.png?20230328082613854)
Figure 3. Force analysis of vibrator and float
图3. 振子、浮子受力分析
振子受力分析:
(7)
沿中轴运动方向运动方程为:
(8)
垂直中轴运动方向方程为:
(9)
浮子受力分析:对于浮子,只考虑竖直方向的垂荡运动。我们将振子重力项并入PTO装置弹簧恢复力,浮子重力项并入静水恢复力中,综上得竖直方向力为:
(10)
波浪能装置转动惯量计算:
振子对转轴的转动惯量:
(11)
浮子外壳对转轴的转动惯量:圆柱高度H1、圆锥高度H2
(12)
振子转矩受力分析:
(13)
浮子受力分析:
(14)
装置初始状态方程:
(15)
3. 数值计算
3.1. 势流理论模型假设成立条件
据上文中建立的势流理论模型求解波浪能装置整体基于海平面的静止状态特征,假设:
1、将海平面作为参考系;
2、设置垂直于海平面竖直向上为竖直正方向;
3、浮子初始振荡频率与波浪激励振动频率相同,无滞后状态;
3.2. 垂荡/纵摇浮子系统运动状态计算
将式(15)经过拉普拉斯变换求解得到高阶偏微分方程组,运用数值计算软件中计算,解出浮子和振子在波浪激励力
(f为波浪激励力振幅,
为波浪频率)作用下前40个波浪周期内时间间隔为0.2 s的垂荡位移和速度以及纵摇角位移和角速度,再将t = 10 s、20 s、40 s、60 s、100 s带入式子中,得到下列数据,如表1、表2所示。
利用编程软件解出的结果导出为可视化的图像,有利于直观地展示其中的关系,如图4、图5所示。
![](Images/Table_Tmp.jpg)
Table 1. Table of heave displacement and velocity of float and vibrator
表1. 浮子和振子的垂荡位移和速度表
![](Images/Table_Tmp.jpg)
Table 2. Pitch angular displacement and angular velocity of float and vibrator
表2. 浮子和振子的纵摇角位移和角速度表
![](//html.hanspub.org/file/1-2820139x30_hanspub.png?20230328082613854)
![](//html.hanspub.org/file/1-2820139x29_hanspub.png?20230328082613854)
![](//html.hanspub.org/file/1-2820139x32_hanspub.png?20230328082613854)
![](//html.hanspub.org/file/1-2820139x31_hanspub.png?20230328082613854)
Figure 4. Time chart of vibrator and float movement
图4. 振子与浮子运动时间图
![](//html.hanspub.org/file/1-2820139x34_hanspub.png?20230328082613854)
![](//html.hanspub.org/file/1-2820139x33_hanspub.png?20230328082613854)
![](//html.hanspub.org/file/1-2820139x36_hanspub.png?20230328082613854)
![](//html.hanspub.org/file/1-2820139x35_hanspub.png?20230328082613854)
Figure 5. Rotation time diagram of vibrator and float
图5. 振子与浮子转动时间图
3.3. 浮子最大功率计算
在先前的研究中,PTO系统可以是直线电机或液压PTO系统,以液压PTO系统为例,系统中的蓄能器、油液刚度等要素决定了系统刚度kpto,通常kpto的选取主要考虑减少系统内部流量、压力波动,其作为储能元件不参与能量吸收 [4]。而该PTO系统中则主要以阻尼装置作为功率转化输出系统,取浮子为研究对象,建立其运动方程,得到其功率计算方程:
(16)
同时运用局部遍历算法进行计算,考虑浮子在波浪中只做垂荡和纵摇的情形,针对直线阻尼器和旋转阻尼器的阻尼系数均为常量且为相同值的情况,当直线阻尼器和旋转阻尼器的阻尼系数均在区间[0,100000]内取值,计算出最大输出功率及相应的最优阻尼系数为1000,如图6所示。同时还可以看出随着时间变换,其功率会呈现一个震荡的状态,但在高阻尼处的平均功率明显大于低阻尼状态。
![](//html.hanspub.org/file/1-2820139x38_hanspub.png?20230328082613854)
Figure 6. Time relationship diagram of output power-PTO device damping coefficient
图6. 输出功率-PTO装置阻尼系数时间关系图
同时得出最大输出功率时浮子、振子的运动状态,如图7所示:
![](//html.hanspub.org/file/1-2820139x40_hanspub.png?20230328082613854)
![](//html.hanspub.org/file/1-2820139x39_hanspub.png?20230328082613854)
![](//html.hanspub.org/file/1-2820139x42_hanspub.png?20230328082613854)
![](//html.hanspub.org/file/1-2820139x41_hanspub.png?20230328082613854)
Figure 7. Stable motion state of vibrator and float under maximum output power
图7. 最大输出功率状态下振子和浮子的稳定运动状态
4. 总结
本文对于垂荡纵摇振荡体式波浪能装置进行数学动力学模型建立,分析了装置随波浪的运动特性以及PTO装置输出功率的影响因素,并计算出波浪能装置最大输出功率状态下的阻尼系数以及振子和浮子的稳定运动状态,发现此时振子和浮子的运动状态近似呈现“拍振”,即振子和浮子运动频率很相近时PTO装置的输出功率最大。蓸雪玲等人基于流体仿真软件对振荡浮子式波浪能发电设备的运动特性进行研究,分析不同阻尼工况的运动情况并进行数值模拟、实验验证 [5]。相较于该种方法,本文模型分析方法具有更好的适用性,不仅能够分析PTO装置对波浪能的利用,也适合被用来简单地分析所有振子–浮子系统组合形成的水面设施。验证了垂荡纵摇振荡体式波浪能输出装置发电方案的可行性,为未来波浪能装置输出功率提升提供理论依据。