浙江大学学报(工学版), 2020, 54(7): 1355-1361 doi: 10.3785/j.issn.1008-973X.2020.07.014

机械与能源工程

水平轴海流能机组传动链主动阻尼控制技术

李阳健,, 何伟, 刘宏伟,, 李伟, 林勇刚, 顾亚京

Active damping control for drive train of horizontal-axis tidal current turbines

LI Yang-jian,, HE Wei, LIU Hong-wei,, LI Wei, LIN Yong-gang, GU Ya-jing

通讯作者: 刘宏伟,男,副教授. orcid.org/0000-0003-3389-1869. E-mail: zju000@163.com

收稿日期: 2019-05-31  

Received: 2019-05-31  

作者简介 About authors

李阳健(1992—),男,博士生,从事海流能开发研究.orcid.org/0000-0002-9038-7784.E-mail:liyangjian@126.com , E-mail:liyangjian@126.com

摘要

为了减小海流能发电机组在弱阻尼状态下的传动链谐振,借鉴风电领域的动态加阻技术,为海流能传动链设计主动阻尼控制器,通过海流能发电机组的转矩控制系统,给发电机施加和转速波动方向相反的波纹转矩. 主动阻尼控制使传动链系统极点远离虚轴,为传动系统增加等效阻尼,减小谐振. 为了准确得到传动链的低阶模态,利用多体系统法和有限元仿真方法,对传动系统的参数进行辨识,根据传动系统拖动实验的转速波动信号,对理论分析的结果进行验证. 通过Matlab/Simulink的建模仿真以及650 kW海流能机组在低流速工况下的等效海试实验,验证了主动阻尼控制对传动系统的有良好的减振效果.

关键词: 海流能机组 ; 模态分析 ; 阻尼控制 ; 传动系统 ; 载荷控制

Abstract

An active damping controller was designed based on previous research in wind energy in order to minimize the resonant torsional oscillations of the weakly damped transmission system of tidal current turbines. The proposed method of active damping strategy added the compensation torque in opposite direction to the speed deviation. The poles of drivetrain system moved away from imaginary axes which indicated that the equivalent damping increased to reduce the resonant torsional load. Multibody system method and finite element analysis were used to identify the parameters of gearbox system in order to obtain the accurate low order torsional modes of transmission system. The theoretical result of the drivetrain resonant frequency was experimentally verified by an onshore experiment. A simulation model for a 650 kW tidal current turbine was constructed in Matlab/Simulink. A substitution sea trial experiment in low current speed condition was conducted to validate the feasibility of proposed active damping.

Keywords: tidal current turbine ; modal analysis ; damping control ; transmission system ; load control

PDF (2059KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

李阳健, 何伟, 刘宏伟, 李伟, 林勇刚, 顾亚京. 水平轴海流能机组传动链主动阻尼控制技术. 浙江大学学报(工学版)[J], 2020, 54(7): 1355-1361 doi:10.3785/j.issn.1008-973X.2020.07.014

LI Yang-jian, HE Wei, LIU Hong-wei, LI Wei, LIN Yong-gang, GU Ya-jing. Active damping control for drive train of horizontal-axis tidal current turbines. Journal of Zhejiang University(Engineering Science)[J], 2020, 54(7): 1355-1361 doi:10.3785/j.issn.1008-973X.2020.07.014

海流能由于高可预测性和高能量密度,具有较高的开发利用价值. 国际上最先进的海流能机组的单机额定功率已经达到2 MW,正在向商业化运营发展. 虽然海流能机组和风力机有着相似的工作原理,但是由于工作介质的不同,海流能机组的发展面临着更大的挑战[1].

近年来,水平轴式海流能机组由于高效稳定的特点,逐渐被更多采用和研究[2]. 传动链作为水平轴海流能机组能量传递的核心部件,具有较大的柔性和较低的阻尼,在复杂的外部激励下,极易发生共振,影响机组的安全和寿命;因此,给传动链增加阻尼对提升机组安全性十分重要[3].

尽管直接为传动链添加物理阻尼是可行的,但是增加了机组的成本;因此在风电领域,利用发电机转矩控制对传动链进行主动阻尼控制是抑制传动链谐振的主要方法[4],在水平轴风力机的相关技术研究中取得了较好的结果[5]. 主动阻尼控制技术使得系统极点远离虚轴,使得振动更快衰减[4, 6]. 此外,White等[7]使用虚拟惯量法,增加系统的转动惯量,减小系统幅频增益. 虚拟惯量法在实际实验中难度大,效果不佳,因此动态加阻法对传动链减振更具有针对性.

尽管在风电领域,主动阻尼技术较成熟,但是目前在海流能发电领域对机组载荷控制的研究甚少. 海流能发电技术尚未完全成熟,我国各高校和研究院所近期的研究重点集中在海流发电机组的机电液装备开发和机组功率等级提高[8-10],对海流能机组的传动链的谐振载荷控制的研究尚处于空白状态.

本文以浙江大学650 kW海流能机组为研究对象,对传动链进行数学建模,结合计算机软件仿真和实验数据,对传动链的低阶扭转模态进行分析;结合海流能机组控制系统,设计主动阻尼加载控制器. 在软件仿真的基础上,开展初步的海试试验验证,结果表明,设计的阻尼控制器能够有效抑制传动链谐振载荷.

1. 传动链模态分析

1.1. 传动链模型

浙江大学650 kW水平轴海流能机组的机械部分主要由叶轮、两级行星齿轮和发电机组成,如图1所示.

图 1

图 1   传动链的结构示意图

Fig.1   Schematic diagram of drive chain


为了研究机组传动链的固有频率,以便对低阶模态进行加阻控制,对该传动系统进行等效建模. 目前,有限元法和多体系统法是对复杂传动系统建模的2种主要方法. 有限元模型虽然能够获得很高的精度,但计算量很大,仿真时间很长;多体系统法是将一个机械系统划分成由若干个柔性弹簧阻尼器连接的刚体,选择合适的刚体数量能够更好地逼近原系统. 本文采用多体系统方法,建立传动链模型.

在风力发电机组的研究中广泛采用两质量弹簧阻尼模型,该方法不能很好地反映传动链的动态特性. 为了获得更加准确的特性,结合650 kW海流能机组齿轮箱的质量分布情况,采用四质量弹簧阻尼模型. 该等效系统的数学模型如下:

$\begin{array}{l} \left[\!\! {\begin{array}{*{20}{c}} {{I_1}}&0&0&0 \\ 0&{{I_2} + i_{_1}^2I_{_2}'}&0&0 \\ 0&0&{{I_3} + i_{_1}^2I_3'}&0 \\ 0&0&0&{{I_4}} \end{array}}\!\! \right]\left[ \!\!{\begin{array}{*{20}{c}} {{{\ddot \theta }_1}} \\ {{{\ddot \theta }_2}} \\ {{{\ddot \theta }_3}} \\ {{{\ddot \theta }_4}} \end{array}} \!\!\right] = \\ \left[\!\! {\begin{array}{*{20}{c}} { - {k_1}}&{{k_1}}&0&0 \\ {{k_1}}&{ - {k_1} - i_{_1}^2{k_2}}&{{i_1}{k_2}}&0 \\ 0&{{i_1}{k_2}}&{ - {i_1}{k_2} - i_{_2}^2{k_3}}&{{i_2}{k_3}} \\ 0&0&{{i_2}{k_3}}&{ - {k_3}} \end{array}} \!\!\right]\left[\!\! {\begin{array}{*{20}{c}} {{\theta _1}} \\ {{\theta _2}} \\ {{\theta _3}} \\ {{\theta _4}} \end{array}}\!\! \right] + \\ \left[\!\! {\begin{array}{*{20}{c}} { - {b_1}}&{{b_1}}&0&0 \\ {{b_1}}&{ - {b_1} - i_{_1}^2{b_2}}&{{i_1}{b_2}}&0 \\ 0&{{i_1}{b_2}}&{ - {i_1}{b_2} - i_{_2}^2{b_3}}&{{i_2}{b_3}} \\ 0&0&{{i_2}{b_3}}&{ - {b_3}} \end{array}} \!\!\right]\left[\!\! {\begin{array}{*{20}{c}} {{{\dot \theta }_1}} \\ {{{\dot \theta }_2}} \\ {{{\dot \theta }_3}} \\ {{{\dot \theta }_4}} \end{array}}\!\! \right] + \left[\!\!{\begin{array}{*{20}{c}} {{T_{\rm{r}}}} \\ 0 \\ 0 \\ { - {T_{\rm{g}}}} \end{array}} \!\!\right] . \\ \end{array} $
(1)

式中:i1i2分别为两级行星齿轮的传动比,I1为叶片和轮毂的转动惯量,I2为一级传动的行星架的转动惯量,I2为一级传动的太阳轮的转动惯量,I3为二级传动的行星架的转动惯量,I3为二级传动的太阳轮的转动惯量,I4为发电机刹车盘等的转动惯量,ki为各轴的刚度,Tr为叶轮转矩,Tg为发电机转矩.

由于传动链处于弱阻尼状态,阻尼系数较小,且阻尼对模态分析结果没有影响,在对传动链进行模态分析的过程中忽略阻尼.

1.2. 传动链模态仿真分析

为了得到传动链系统的模态,需要确定式(1)中的系统参数. 在Pro/E软件中建立650 kW机组传动链的三维模型,如图2所示;分别获取传动链各个部分的转动惯量,如表1所示.

表 1   各个部分转动惯量

Tab.1  Inertia calculated for modeling of system kg·m2

参数 数值
叶片和轮毂 ${I_1}$ 54 151.73
一级传动的行星架 ${I_2}$ 445.68
一级传动的太阳轮 $I_2'$ 404.35
二级传动的行星架 ${I_3}$ 23.96
二级传动的太阳轮 $I_3'$ 19.59
发电机刹车盘等 ${I_4}$ 280.50

新窗口打开| 下载CSV


图 2

图 2   650 kW海流能机组的传动链系统

Fig.2   Drive chain system of 650 kW tidal current energy generating set


传动链系统中每个部分轴的刚度Ksh可以通过下式计算得到:

${K_{\rm{sh}}} = {{G I}}/{L}.$
(2)

式中:G为切变模量,I为截面极惯性矩,L为轴的长度.

由于实际传动链轴的直径不一,通过理论计算的刚度精度较低,无法满足控制需求. 利用有限元仿真软件ANSYS对各轴进行有限元仿真分析,得到施加单位力矩下轴的角位移形变,计算轴的刚度:

${K_{\rm{sh}}} = {{\Delta T}}/{{\Delta \theta }}.$
(3)

式中: $\Delta T$$\Delta \theta $分别为施加的力矩增量和轴的角位移形变.

低速轴的有限元分析如图3所示.

图 3

图 3   低速轴刚度有限元分析

Fig.3   Finite element model analysis for stiffness of LSS


依次对各轴进行有限元分析,得到各轴刚度,如表2所示.

表 2   系统中轴段刚度

Tab.2  Stiffness value for modeling of system 104 N·m/rad

参数 数值
低速轴k1 5.356×103
中间轴k2 4.205×103
高速轴k3 4.601

新窗口打开| 下载CSV


根据式(1),绘制传动系统方块图,如图4所示. 根据系统传递函数,绘制传动链的波特图,如图5所示. 图中,Am为幅值,ϕ为相位,f为频率. 通过波特图可得,该传动链的一阶模态为18.3 rad/s(即2.91 Hz).

图 4

图 4   传动链方块图

Fig.4   Block diagram of transmission system


图 5

图 5   传动系统的Bode图

Fig.5   Mode graph of drive train system


1.3. 传动链模态实验验证

为了验证理论得到的传动链模态的准确性,选取传动链系统拖动实验(见图6)中恒转矩加载的转速数据,对转速数据进行频谱分析,如图7所示,图中,ngA为发电机转速幅值.

图 6

图 6   传动链厂内拖动测试

Fig.6   Onshore test of drive train


图 7

图 7   发电机转速的频谱图

Fig.7   Spectrogram of generator speed


图7可知,得到频谱峰值对应的频率为2.44 Hz,与理论分析结果2.91 Hz相近,但是存在一定的误差. 这表明基于四质量块的传动链等效数学模型能够有效估计传动齿轮箱的低阶模态. 造成误差的可能原因如下:四质量块模型是基于集中质量的假设,忽略了轴的分布质量问题,存在一定的模型误差;忽略了联轴器的影响,使得计算得到的刚度比实际刚度大;忽略了齿轮箱和电机中的润滑油造成的附加质量,使得模型中的质量块转动惯量偏小;忽略了轴系一些小零部件的装配和齿轮的啮合刚度. 在更精准的模型中,可以考虑这些因素,得到更准确的共振频率.

2. 阻尼控制器设计

由于传动链的低阶共振频率较低,阻尼较低,特别是在恒转矩控制过程中,极易发生扭振,造成剧烈的扭矩波动,对传动链造成严重的有害扭振载荷. 本文采用施加电磁转矩的方式,给传动系统人为地增加虚拟阻尼,减小传动链的扭振. 主动阻尼控制原理较简单,通过在给定发电机转矩上附加一个转矩成分 $\Delta {T_{\rm{g}}}$,如下:

$J\ddot \theta + B\dot \theta + k\theta = {T_{\rm{r}}} - {i_1}{i_2}({T_{\rm{g}}} + \Delta {T_{\rm{g}}}).$
(4)

式中:JBk分别为传动链相对于低速轴的等效转动惯量、等效阻尼和等效刚度. 使得

$\Delta {T_{\rm{g}}} = {B_{\rm{g}}}\dot \theta, $
(5)

则有

$J\ddot \theta + (B + {B_{\rm{g}}})\dot \theta + k\theta = {T_{\rm{r}}} - {i_1}{i_2}{T_{\rm{g}}}.$
(6)

即相当于通过电磁转矩控制,给传动链增加了一个额外的阻尼,达到减振的效果. 式中:Bg为电磁阻尼系数.

阻尼控制器须获取有害的转矩成分,即通过检测低速轴和高速轴的转矩完成,但由于安装扭矩传感器不切实际,可以通过对高速轴转速信号进行适当的滤波,获取传动链扭振的信息. 高速轴转速信号经过带通滤波器,得到有害的转速波动信号,并通过海流能机组转矩控制器,施加一个与传动链转速波动方向相反的补偿转矩 $\Delta {T_{\rm{g}}}$,附加到常规控制的最佳转矩上,从而增加传动系统的阻尼,达到减振的效果. 控制系统框图如图8所示.

图 8

图 8   主动阻尼控制系统图

Fig.8   Block diagram of active damping control system


图8中,带通滤波器的传递函数为

$G(s) = \frac{{2\zeta \omega s}}{{{s^2} + 2\zeta \omega s + {\omega ^2}}}.$
(7)

式中: $\zeta $为滤波器的阻尼比, $\omega $为滤波器的中心频率.

由于转速传感器以及控制器的控制信号均为离散信号,且信号通讯、转矩控制及滤波器均会产生一定的延时[11](一般可以等效为一阶惯性环节),直接进行主动阻尼控制,极有可能因为相位滞后、系数选用不当、通讯延迟等原因,不能起到加阻效果,甚至引起振荡. 为了实现良好的控制效果,在滤波器滞后添加相位超前矫正:

$H(s) = \frac{{\alpha Ts + \alpha }}{{\alpha Ts + 1}}.$
(8)

3. 仿真与实验

3.1. 基于MATLAB/Simulink的主动阻尼控制仿真

基于上述控制原理,在Simulink中搭建海流能发电系统的仿真模型,如图9所示. 模型主要由4个部分组成:叶轮、传动系统、控制系统及阻尼控制部分. 叶轮部分主要由Tidal-Bladed中对叶轮的λ-功率系数Cp曲线,采用二维插值法来实现. 传动系统主要由如图4所示的传动链控制框图组成. 海流能机组的控制系统通过发电机转速给定发电机转矩,实现海流能机组的控制[12-13].

图 9

图 9   Simulink仿真框图

Fig.9   Diagram of Simulink simulation


在Simulink仿真模型中,对机组采用恒转矩控制,分别在不同的电磁阻尼系数下对发电机组进行加阻控制,发电机的转速波动仿真结果如图10~12所示,图中,ωg为发电机角速度,t为时间.

图 10

图 10   加阻控制仿真结果(Bg=50 N·m·s/rad)

Fig.10   Simulation result of damping control (Bg=50 N·m·s/rad)


图 11

图 11   加阻控制仿真结果(Bg =100 N·m·s/rad)

Fig.11   Simulation result of damping control(Bg=100 N·m·s/rad)


图 12

图 12   加阻控制仿真结果(Bg=200 N·m·s/rad)

Fig.12   Simulation result of damping control(Bg=200 N·m·s/rad)


仿真结果表明,采用主动阻尼控制能够有效地降低传动链谐振;在一定的范围内,更大的阻尼系数的减振效果更显著. 随后的仿真结果表明,在阻尼系数超过一定的范围之后,随着阻尼系数的增加,减振效果不再增强,有增加振荡的趋势,这可能是带相位补偿的控制器无法完全克服由通讯延迟、执行器件的惯性等带来的相位误差. 这表明主动阻尼控制技术可能会给系统带来潜在的风险,选择合理的阻尼系数至关重要. 仿真结果表明,在阶跃激励下,与风电领域的的相关研究报道[6]相比,海流能机组的谐振更剧烈. 这可能是海流能机组更小的叶片和转动惯量导致的.

由于系统自身的阻尼难以通过理论确定,该仿真结果是在传动系统无阻尼的理想情况下完成的原理性验证. 实际的传动链系统尽管处于弱阻尼状态,但是由于密封件、润滑油、叶轮在水中的水阻力等,存在一定的阻尼,使得实际的传动链谐振振幅比仿真所显示的结果小.

3.2. 海试实验验证

浙江大学650 kW海流能发电机组于舟山摘箬山岛进行海试实验,如图13所示. 由于在额定工况下进行无阻尼控制实验引起传动链振荡,对机组可能造成严重损伤,采用在低流速下的恒转矩控制进行替代性实验. 实验结果如图14所示. 图中,ng为发电机转速. 叶轮在海流的驱动下,开始旋转,并平稳并网. 由于系统的弱阻尼状态,系统逐渐开始振荡,最终触发了机组的停机保护操作.

图 13

图 13   浙江大学650 kW海流能发电机组

Fig.13   650 kW tidal current turbine of ZJU


图 14

图 14   海试未加阻实验转速结果

Fig.14   Speed of sea trail without damping control


与仿真中的无阻尼系统相比,实际系统引起振荡的过程更缓慢,可见系统本身存在一定的阻尼,但是不足以抑制传动链的振荡. 在海流能发电机组的控制系统中添加主动阻尼控制,根据仿真结果调节阻尼,得到的试验结果如图15所示.

图 15

图 15   主动阻尼控制海试结果

Fig.15   Result of sea trail with active damping control


结果表明,通过主动阻尼控制,海流能机组的传动链谐振被显著抑制,发电机转速的振动幅度显得很微小,低频的转速波动主要是由流速变化引起的. 与风电领域的相关研究结果比较,在类似的主动阻尼控制下,海流能机组的谐振抑制更彻底[3, 7],这极有可能是因为海流能机组由于动密封件、配油机构结构,额外给传动链增加了物理阻尼.

4. 结 论

(1)基于多体系统法和有限元仿真,对海流能发电系统的传动系统进行模态分析. 所得的理论结果通过海流能发电机组的厂内拖动试验结果进行验证. 结果表明,采用的多体系统法具有一定的准确性,理论结果与试验结果相比,误差小于20%,对产生误差的原因进行定性分析.

(2)基于海流能发电系统的控制器,设计主动阻尼控制器;根据可能的系统相位滞后,设计超前矫正环节.

(3)在不考虑系统自有阻尼的理想条件下,搭建海流能发电机组主动阻尼控制的Matlab/Simulink仿真模型,研究不同阻尼系数下系统的减振效果. 仿真结果表明,选择合适的阻尼系数,主动阻尼控制可以有效地降低传动链谐振,与未加阻系统相比,传动链转速波动幅值减小69%以上;同时得到更加稳定的叶尖速比,使得机组效率提高约4.3%.

(4)开展650 kW海流能机组在低流速工况的恒转矩控制下的主动阻尼控制实验,实验验证了所设计的主动阻尼控制器对传动链的谐振具有很强的减振效果.

参考文献

WINTER A I. Differences in fundamental design drivers for wind and tidal turbines [C]//Oceans. Santander, Spain: IEEE, 2011: 1-10.

[本文引用: 1]

ZHOU Z, BENBOUZID M, CHARPENTIER J F, et al

Developments in large marine current turbine technologies: a review

[J]. Renewable and Sustainable Energy Reviews, 2017, 71: 852- 858

DOI:10.1016/j.rser.2016.12.113      [本文引用: 1]

GENG H, XU D, WU B, et al

Comparison of oscillation damping capability in three power control strategies for PMSG-based WECS

[J]. Wind Energy, 2011, 14 (3): 389- 406

DOI:10.1002/we.428      [本文引用: 2]

BOSSANYI E A

Wind turbine control for load reduction

[J]. Wind Energy, 2003, 6 (3): 229- 244

DOI:10.1002/we.95      [本文引用: 2]

彭超

风力发电机组传动链调谐惯量阻尼器研究

[J]. 太阳能学报, 2018, 39 (7): 2036- 2043

[本文引用: 1]

PENG Chao

Research of tuned inertia damper for wind turbine drive rain

[J]. Acta Energiae Solaris Sinica, 2018, 39 (7): 2036- 2043

[本文引用: 1]

MANDIC G, NASIRI A, MULJADI E, et al

Active torque control for gearbox load reduction in a variable-speed wind turbine

[J]. IEEE Transactions on Industry Applications, 2012, 48 (6): 2424- 2432

DOI:10.1109/TIA.2012.2227131      [本文引用: 2]

WHITE W N, FATEH F, YU Z. Torsional resonance active damping in grid tied wind turbines with gearbox, DFIG, and power converters [C]//American Control Conference. Chicago: IEEE, 2015.

[本文引用: 2]

LIU H W, MA S, LI W, et al

A review on the development of tidal current energy in China

[J]. Renewable and Sustainable Energy Reviews, 2011, 15 (2): 1141- 1146

DOI:10.1016/j.rser.2010.11.042      [本文引用: 1]

LIU H W, LI W, LIN Y G, et al

Tidal current turbine based on hydraulic transmission system

[J]. Journal of Zhejiang University-Science A: Applied Physics and Engineering, 2011, 12 (7): 511- 518

LI Y J, LIU H W, LIN Y G, et al

Design and test of a 600 kW horizontal-axis tidal current turbine

[J]. Energy, 2019, 182: 177- 186

DOI:10.1016/j.energy.2019.05.154      [本文引用: 1]

顾亚京, 李伟, 刘宏伟, 等

潮流发电系统功率控制及并网运行

[J]. 浙江大学学报:工学版, 2017, 51 (10): 1974- 1980

[本文引用: 1]

GU Ya-jing, LI Wei, LIU Hong-wei, et al

Study on power control and grid-connected operation for tidal current energy conversion system

[J]. Journal of Zhejiang University: Engineering Science, 2017, 51 (10): 1974- 1980

[本文引用: 1]

XU Q K, LIU H W, LIN Y G, et al

Development and experiment of a 60 kW horizontal-axis marine current power system

[J]. Energy, 2015, 88: 149- 156

DOI:10.1016/j.energy.2015.04.018      [本文引用: 1]

GU Y J, LIU H W, LI W, et al

Integrated design and implementation of 120 kW horizontal-axis tidal current energy conversion system

[J]. Ocean Engineering, 2018, 158: 338- 349

DOI:10.1016/j.oceaneng.2018.04.017      [本文引用: 1]

/