基于垂荡运动波浪能发电PTO系统阻尼研究.pdf

SHANDONG ELECTRIC POWER山东电力技术第50卷（总第310期）2023年第9期DOI：10.20097/ki.issn1007-9904.2023.09.004基于垂荡运动波浪能发电PTO系统阻尼研究平力，尚劲光，陈茂祥，凌云凤，刘眉洁（青岛大学物理科学学院，山东青岛266071）摘要：波浪能作为海洋资源中储藏量最丰富的清洁能源，其发展前景十分广阔。因此，如何解决波浪能能量转换，实现波浪能最大层次的规模利用，是当前面临的关键问题。针对求解最优阻尼系数和最大发电功率，建立波浪能装置垂荡运动模型，分别讨论常数阻尼和线性阻尼两种情况下浮子和振子的垂荡运动以及对应的最优阻尼系数和发电功率。通过 ODE45函数、变步长搜索法等精确求解得出：当系统采用常数阻尼时最大发电功率124.20 W，最优阻尼系数为37 940 Ns/m；当系统采用线性阻尼时最大发电功率为124.22 W，最优阻尼系数为46 025 Ns/m。关键词：波浪能；垂荡；阻尼系数；ODE45函数；变步长搜索法中图分类号：P743.2；O17文献标识码：A文章编号：1007-9904（2023）09-0028-07Research on Damping of PTO System Based on Vertical MotionWave Energy Power GenerationPING Li，SHANG Jinguang，CHEN Maoxiang，LING Yunfeng，LIU Meijie（College of Physics，Qingdao University，Qingdao 266071，China）Abstract：As the most abundant clean energy in marine resources，wave energy has a broad development prospect.Therefore，howto solve the wave energy conversion and realize the maximum scale utilization of wave energy is the key problem at present.Inorder to solve the optimal damping coefficient and maximum power generation，the vertical motion model of the wave energydevice is established，and the vertical motion of the float and oscillator under constant damping and linear damping are discussedrespectively，as well as the corresponding optimal damping coefficient and power generation.Through precise solutions such asODE45 function and variable step search method，it can be concluded that the maximum power generation power is 124.20 Wwhen the system adopts constant damping，and the optimal damping coefficient is 37 940 Ns/m；when the system adopts lineardamping，the maximum power generation is 124.22 W，and the optimal damping coefficient is 46 025 Ns/m.Keywords：wave energy；vertical motion；damping coefficient；ODE45 functions；variable step search method0引言在当前社会背景下，人类社会对于可再生能源的需求愈来愈强1，海洋面积约占地球面积的71，海洋中能开发利用的可再生能源也多种多样2-4。近年来，随着海洋能开发技术不断完善，国内外推出了大量政策鼓励社会全面推广使用海洋能源5-8。波浪能作为海洋中一种主要的清洁能源，蕴藏量大，可再生能力强，具有较好的开发前景9-10。但由于波浪能的不稳定性11，开发合适的波浪能装置，提高系统的能量转换率，是波浪能最大化利用的关键所在。为解决该问题，有学者提出利用振荡浮子式波浪发电分段控制策略，以发电效率为基础，结合低速和高速2种发电模式，优化发电量，但该策略下的发电功率不稳定，发电效率低下，无法持续供电12-13；还有学者提出一种可应用于多浮子波浪能发电装置的液压发电系统14，采用蓄能器稳压控制策略，但是实际工作海域的波浪幅度不稳定，蓄能器的选型比较困难，并且成本高，不适合全面推广。综合考虑发电总量、发电效率、系统稳定性、系统成本等因素，针对波浪能装置的垂荡运动情况建立合理的模型，并对模型进行精确求解。该模型可以适应不同的海洋场景，具有普适性，并且成本低，基金项目：国家自然科学基金项目（41976173）；山东省自然科学基金项目（ZR2019MD016）。National Natural Science Foundation of China（41976173）；ShandongProvincial Natural Science Foundation（ZR2019MD016）.新能源 28具有一定的推广价值。1波浪能装置波浪能装置主要由浮子、振子、中轴以及能量输出（power take off，PTO）系统构成15，其中 PTO 系统由弹簧和阻尼器组成，其模型如图1所示。图1波浪能装置Fig.1 Wave energy device图1中，波浪能装置中轴底座固定于隔层的中心位置，弹簧和直线阻尼器一端固定在振子上，一端固定在中轴底座上，振子沿中轴做往复运动。直线阻尼器的阻尼力与浮子和振子的相对速度成正比，比例系数为阻尼系数。在实际中，阻尼系数会直接影响系统的工作效率，而阻尼系数不是定值，其值与材料、结构、阻尼模型有关16。阻尼系数增大，则直线阻尼器的输出阻抗逐渐小于额定负载阻抗，输出阻抗直接影响系统的低频特性。输出阻抗减小，则输出的平均功率不断增大。当增大到一定峰值时，系统转向低频工作，则直线阻尼器的阻尼系数继续增大，影响到工作效率，使输出的平均功率小幅度下滑。选取两种特定情况下的阻尼系数，即常数PTO阻尼系数和线性PTO阻尼系数进行研究。该装置在波浪的冲击下，浮子运动带动振子运动，两者的相对运动会驱使阻尼器做功并输出能量。波浪能的基本形式是微幅波，因此浮子在线性周期微幅波作用下会受到波浪激励力、附加质量力、辐射阻尼力和静水恢复力17，由于海水的黏度和旋转特性、装置内部的摩擦力等因素对能量转化影响较小，在研究中忽略不计。2垂荡运动模型垂荡运动是指浮子在波浪作用下做垂直往复运动。对浮子及振子进行受力分析，得到整个系统中各部件的受力情况如图2所示。图2垂荡运动系统受力分析Fig.2 Force analysis of vertical motion system以浮子在静水中的平衡状态x0为零点，浮子的相对位移为x，振子的相对位移为z，浮子的相对速度为x，振子的相对速度为z，以沿水平面向上为正方向，得到各个力的表达式。2.1波浪激励力波浪激励力，即波浪荷载，由波浪水质点与波浪能装置的相对运动所引起，是波浪对放置在海洋中的波浪能装置产生的推动力。波浪是一种随机性运动，浮子做垂荡运动时，波浪激励力为F0=Acos(t)（1）式中：A为波浪激励力振幅；为波浪频率；t为波浪运动的时间。2.2附加质量力波浪能装置在海水中做垂荡运动时，会引起装置周围的流体运动，因此要使装置在海水中获得一定的加速度，施加在装置上的力等于装置质量与此加速度的乘积和带动装置周围海水运动所需额外力的总和。而该额外力会对应产生一个虚拟质量，称为附加质量18，将此附加质量设为m0，得到附加质量力为F1=m0g（2）式中：g为重力加速度。2.3辐射阻尼力波浪能装置在海水中作垂荡运动时，会带动周围水体流动，而周围水体流动会对波浪能装置的垂荡运动产生一定的阻力作用。此力与垂荡运动的速平力，等：基于垂荡运动波浪能发电PTO系统阻尼研究29山东电力技术第50卷（总第310期）2023年第9期度成正比，方向相反，比例系数为辐射阻尼系数，得到辐射阻尼力为F2=-x（3）2.4弹簧弹力当振子沿中轴做往复运动时，弹簧会给振子趋向平衡位置的弹力，根据胡克定律，弹簧给振子的弹力表达式为F3=kx（4）式中：k为弹簧的弹性系数；x为弹簧的伸长量。2.5静水恢复力静水恢复力是指为浮子做垂荡运动时使浮子回到平衡位置的作用力，与垂荡运动过程中的浮力变化有关。静水恢复力为F4=Ff(m0+m1+m2)g（5）式中：Ff为浮子在水中受到的浮力；m1和m2分别为浮子和振子的质量。2.6直线阻尼力在该波浪能装置中，PTO系统会将浮子在垂荡运动中克服阻尼力所做的功转化为电能。采用线性阻尼模型分析直线阻尼器的阻尼力。直线阻尼器的阻尼力与浮子和振子的相对速度成正比，其比例系数即阻尼系数为，则直线阻尼力为FPTO=(x-z)（6）2.7浮子与振子的运动模型振子受力包括弹簧弹力、重力、直线阻尼力，以垂直向上为正方向，振子在波浪作用下受力情况为Fz=F3-m2g-FPTO（7）浮子受力包括波浪激励力、静水恢复力、辐射阻尼力、弹簧弹力、直线阻尼力、附加质量力。因此，得出浮子和振子整个系统在海洋中的受力情况为F=F0+F2+F4-F3-FPTO（8）进一步可得振子和整个系统的运动微分方程为(m1+m2+m0)x=Fm2z=Fz（9）3模型的求解3.1初始状态求解模型各参数值如表1所示。在初始时刻，浮子和振子在静水中保持平衡。在平衡状态下，浮子的浸水深度大于浮子圆锥部分的高度，则浮子圆锥部分全部浸入水体中，多余的部分为圆柱浸入的体积。求得初始时刻系统的浸水深度为2.8 m，将此时水平面作为零点x0，则平衡时弹簧长度为 0.202 m，振子底部与水平面相距1.798 m。表1模型参数Table 1 Model parameter values参数浮子质量/kg浮子底半径/m浮子圆柱部分高度/m浮子圆锥部分高度/m振子质量/kg振子半径/m振子高度/m海水的密度/（kg/m3）重力加速度/（m/s2）弹性系数/（N/m）弹簧原长/（m）波浪频率/Hz数值4 866130.82 4330.50.51 0259.880 0000.50.7143.2常数PTO阻尼系数海洋波浪是在深海区由风和潮汐力共同作用的结果，潮汐力分为太阳潮汐力和月球潮汐力19，由于太阳对于潮汐力贡献较小，在此忽略不计。所以海波的速度基本与风速和月球引力极大值点在地表的移动速度相同。查阅文献 20 可知，一般海上风速的变化范围为10200 m/s。=V=g2V（10）式中：为深水波频率；V为深水波波速；为深水波波长。由此推出频率变化范围为0.011 Hz。本文波浪频率选取0.714 Hz。当阻尼器的阻尼系数为常量 10 000Ns/m时，利用 MATLAB 软件中ODE45 函数求解 2.7 节建立的运动模型，得到浮子和振子垂荡位移、速度与时间的关系图像，如图 3所示。30（a）垂荡位移（b）垂荡速度图3常数PTO阻尼系数垂荡运动Fig.3 Constant PTO damping coefficient vertical motion由图3可知，当直线阻尼器采用常数阻尼系数时，浮子和振子的位移、速度趋势线吻合度较高，但位移和速度峰值均出现忽大忽小的情况，变化幅度较大；约60 s后，整个系统处于稳定状态，位移和速度均呈周期性变化且振幅近似不变，可近似为简谐运动。系统稳定后，浮子和振子的位移、速度趋势线相似性较高。3.3线性PTO阻尼系数当直线阻尼器采用线性PTO阻尼系数时，取比例系数为 10 000，幂指数为 0.5，由此得到阻尼系数的表达式为=10 000|x-z0.5（11）将式（11）代入式（9），利用 ODE45 函数进行求解，得到此阻尼系数下浮子和振子垂荡位移、速度与时间的关系如图 4 所示。（a）垂荡位移（b）垂荡速度图4线性PTO阻尼系数垂荡运动Fig.4 Linear PTO damping coefficient vertical motion由图4可知，当直线阻尼器采用线性阻尼系数时，浮子和振子的运动曲线几乎重合。浮子系统开始运动时浮子与振子的位移和速度振幅不稳定且无固定周期；约60 s后，浮子和振子的运动近似呈现为简谐振动，有固定的振幅和周期。3.4阻尼系数对发电稳定性的影响对比 3.2 节和 3.3 节计算结果可知，采用线性PTO阻尼时，浮子和振子的垂荡运动在前60 s内更不稳定，位移和速度幅值变化范围更大。经过60 s，系统稳定后，两种阻尼系数下，浮子和振子均近似为简谐运动状态。但当采用线性阻尼系数时，浮子与振子运动周期更短，运动频率更高，且幅值更大。波浪能转换装置的发电稳定性是指装置预热后的发电稳定性21，系统开启后的前60 s可看作装置的预热时间。综合对比可知，当直线阻尼器采用线性阻尼系数时，系统的发电稳定性更好。3.5阻尼系数对功率的影响在波浪的作用下，浮子与振子的相对运动会驱动阻尼器做功，并能将所做的功作为能量输出。功率的表达式为P=FPTO|x z（12）由3.2节和3.3节分析可知，系统稳定后，浮子和振子的运动可以近似看作是简谐振动。则系统平均输出功率为P=12(x-z)2（13）3.5.1常数PTO阻尼对功率的影响直线阻尼器采用常数阻尼时，联立运动模型以及平均输出功率表达式，利用变步长搜索法，绘制出平均输出功率随阻尼系数变化曲线，如图5所示。平力，等：基于垂荡运动波浪能发电PTO系统阻尼研究31山东电力技术第50卷（总第310期）2023年第9期图5常数PTO阻尼下平均输出功率曲线Fig.5 Average output power curve at constant PTOdamping由图 5 可知，随着阻尼系数的增加，平均输出功率呈现出上升趋势且增长速度较快；当阻尼系数达到 2.93104Ns/m 时，增长较为缓慢且图像斜率近似为零；而在阻尼系数达到 3.79104Ns/m时，平均输出功率呈现下降趋势，但下降的速度较为缓慢。平均输出功率最大时对应的阻尼系数范围为 30 00040 000Ns/m。缩小阻尼系数区间后，利用变步长搜索法进行精确求解，最终求得最大平均输出功率为 124.20 W，对应的阻尼系数为37 940Ns/m。将该阻尼系数代入式（13），绘制出每时刻对应的输出功率如图6所示。图6最优常数阻尼系数下的输出功率Fig.6 Output power at optimal constant damping coefficient3.5.2线性PTO阻尼对功率的影响直线阻尼器采用线性阻尼时，得到阻尼系数的表达式为=a|x-zb（14）式中：a为比例系数，a 0，100 000；b为幂指数，b 0，1。将式（14）代入式（9）与式（13），利用变步长搜索法绘制出平均输出功率随a与b的变化曲线，如图7所示。图7线性PTO阻尼下平均输出功率曲线Fig.7 Average output power curve at linear PTO damping由图7可知，随着两个参数的增加，平均输出功率呈现出大幅度上升，而后上升速度减慢并逐渐趋于0，最后缓慢下降的趋势。缩小两个参数取值区间后，利用变步长搜索法进行精确求解，最终求得最大平均输出功率为 124.22 W，对应的比例系数 a和幂指数b分别为97 000和0.4。将这两个参数代入式（13），求解最优阻尼系数为46 025Ns/m，并绘制出每时刻对应的输出功率如图8所示。图8最优线性阻尼系数下的输出功率Fig.8 Output power at optimal linear damping coefficient4结束语建立波浪能发电装置模型，采用ODE45函数求解数值解，研究了该装置在海洋中做垂荡运动的发电性能。在此基础上，进一步求解该装置的最优阻尼系数和最大发电功率，研究了该波浪能装置的最优工作参数，并采用变步长搜索法进行积分。通过代入实际数值对模型进行验证发现，该模型与实际情况拟合较好。因此，所建立的模型具有一定的可行性和研究价值。参考文献1颜越.我国可再生能源发电项目补贴退坡的困境机理解析 J.32能源技术与管理，2022，47（6）：203-205.YAN Yue.Analysis of the difficult mechanism of subsidy decline ofrenewable energy power generation projects in ChinaJ.EnergyTechnology and Management，2022，47（6）：203-205.2彭伟，王芳，王冀.我国海洋可再生能源开发利用现状及发展建议 J.海洋经济，2022，12（3）：70-75.PENG Wei，WANG Fang，WANG Ji.Current situation anddevelopment suggestions of marine renewable energy developmentand utilization in China J.Ocean Economy，2022，12（3）：70-75.3周守为，李清平，朱海山，等.海洋能源勘探开发技术现状与展望 J.中国工程科学，2016，18（2）：19-31.ZHOU Shouwei，LI Qingping，ZHU Haishan，et al.Current statusand prospect of marine energy exploration and developmenttechnology J.Engineering Science，2016，18（2）：19-31.4赵伟国，刘玉田，王伟胜.海洋可再生能源发电现状与发展趋势J.智能电网，2015，3（6）：493-499.ZHAOWeiguo，LIUYutian，WANGWeisheng.Statusanddevelopment trend of marine renewable energy power generationJ.Smart Grid，2015，3（6）：493-499.5吴磊，詹红兵.全球海洋治理视阈下的中国海洋能源国际合作探析 J.太平洋学报，2018，26（11）：56-69.WU Lei，ZHAN Hongbing.Analysis of Chinas marine energyinternational cooperation under the threshold of global marinegovernance J.Pacific Journal，2018，26（11）：56-69.6仲雯雯.我国战略性海洋新兴产业发展政策研究 D.青岛：中国海洋大学，2011.7刘贺青.英国海洋能源产业全球布局背景下的中英海洋能源合作评析与对策 J.太平洋学报，2016，24（10）：39-46.LIU Heqing.Evaluation and countermeasures of Sino-Britishmarine energy cooperation under the background of global layoutof UK marine energy industry J.Pacific Journal，2016，24（10）：39-46.8杨薇，孔昊.基于低碳经济的福建省海洋能源产业发展可行性研究 J.海洋开发与管理，2017，34（11）：61-65.YANG Wei，KONG Hao.Feasibility study on development ofmarine energy industry in Fujian Province based on low-carboneconomyJ.Marine Development and Management，2017，34（11）：61-65.9吴明东，盛松伟，张亚群，等.海洋波浪能浮标发展现状及前景J.新能源进展，2021，9（1）：42-47.WU Mingdong，SHENG Songwei，ZHANG Yaqun，et al.Developmentstatus and prospect of marine wave energy buoyJ.New EnergyProgress，2021，9（1）：42-47.10 肖惠民，于波，蔡维由.世界海洋波浪能发电技术的发展现状与前景 J.水电与新能源，2011（1）：67-69.XIAO Huimin，YU Bo，CAI Weiyou.Development status andprospect of wave energy power generation technology in the worldsocean J.Hydropower and New Energy，2011（1）：67-69.11 刘颖昕.高效稳定的波浪能液压PTO装置设计及控制策略研究D.济南：山东大学，2021.12 姜元，杨金明，谢泽坤，等.振荡浮子式波浪发电系统的分段控制 J.可再生能源，2018，36（12）：1 794-1 799.JIANG Yuan，YANG Jinming，XIE Zekun，et al.Segmentationcontrol of the oscillating floating wave power generation systemJ.Renewable Energy，2018，36（12）：1 794-1 799.13 黄秀秀.振荡浮子式波浪发电系统的功率控制 D.广州：华南理工大学，2019.14 王坤林，盛松伟，叶寅，等.波浪能装置中液压发电系统Boost变换机理及控制策略J.电力系统自动化，2017，41（12）：173-178.WANG Kunlin，SHENG Songwei，YE Yin，et al.The Boosttransformation mechanism and control strategy of the hydraulicpower generation system in the wave energy deviceJ.PowerSystem Automation，2017，41（12）：173-178.15 曹雪玲，高涛，王文胜，等.单浮筒振荡浮子式波浪能装置的动力特性分析 J.太阳能学报，2022，43（11）：364-368.CAO Xueling，GAO Tao，WANG Wensheng，et al.Dynamiccharacteristics analysis of single pontoon oscillating float waveenergy device J.Journal of Solar Energy，2022，43（11）：364-368.16 田坤，李鹏，王元丰.结构等效复阻尼模型对钢梁及钢框架结 构 动 力 响 应 影 响 分 析J.振 动 与 冲 击，2008，27（7）：118-121.TIAN Kun，LI Peng，WANG Yuanfeng.Influence analysis ofstructural equivalent complex damping model on dynamic responseof steel beam and steel frame structure J.Journal of Vibration andShock，2008，27（7）：118-121.17 杨运，张建明，李道朋，等.压力容器自动排气阀工作性能分析及优化 J.液压与气动，2022，46（3）：91-101.YANG Yun，ZHANG Jianming，LI Daopeng，et al.Analysis andoptimization of pressure vessel automatic exhaust valveJ.Hydraulic and Pneumatic，2022，46（3）：91-101.18 马烨，单雪雄.数值计算复杂外形物体附加质量的新方法 J.计算机仿真，2007，24（5）：75-78.MA Ye，SHAN Xuexiong.A new method for numerical calculationof the additional mass of complex shaped objectsJ.ComputerSimulation，2007，24（5）：75-78.19 韦蔚，刘少峰，王凯，等.月岭形成机制及其与潮汐力的相关性分析 J.现代地质，2012，26（1）：198-204.WEI Wei，LIU Shaofeng，WANG Kai，et al.Formation mechanismof lunar ridge and its correlation with tidal forceJ.ModernGeology，2012，26（1）：198-204.20 范开国，徐青，徐东洋，等.星载 SAR 海面风场遥感研究进展平力，等：基于垂荡运动波浪能发电PTO系统阻尼研究33山东电力技术第50卷（总第310期）2023年第9期J.地球物理学进展，2022，37（5）：1 807-1 817.FAN Kaiguo，XU Qing，XU Dongyang，et al.Progress in remotesensing of starborne SAR sea surface wind fieldJ.Progress inGeophysics，2022，37（5）：1 807-1 817.21 盛松伟，游亚戈，张亚群，等.漂浮式波浪能装置能量转换系统研究 J.机械工程学报，2012，48（24）：141-146.SHENG Songwei，YOU Yage，ZHANG Yaqun，et al.Researchon energy conversion system of floating wave energy deviceJ.Journal of Mechanical Engineering，2012，48（24）：141-146.收稿日期：2023-02-21修回日期：2023-05-16作者简介：平力（2002），女，主要研究方向为应用物理学；尚劲光（2003），男，主要研究方向为应用物理学；陈茂祥（2002），男，主要研究方向为应用物理学；凌云凤（2001），女，主要研究方向为应用物理学；刘眉洁（1978），女，通信作者，博士，副教授，主要研究方向为光学和电磁场理论与实验、海洋遥感机理及其应用。（责任编辑车永强）34