水平轴海流能机组传动链主动阻尼控制技术
Active damping control for drive train of horizontal-axis tidal current turbines
通讯作者:
收稿日期: 2019-05-31
Received: 2019-05-31
作者简介 About authors
李阳健(1992—),男,博士生,从事海流能开发研究.orcid.org/0000-0002-9038-7784.E-mail:
为了减小海流能发电机组在弱阻尼状态下的传动链谐振,借鉴风电领域的动态加阻技术,为海流能传动链设计主动阻尼控制器,通过海流能发电机组的转矩控制系统,给发电机施加和转速波动方向相反的波纹转矩. 主动阻尼控制使传动链系统极点远离虚轴,为传动系统增加等效阻尼,减小谐振. 为了准确得到传动链的低阶模态,利用多体系统法和有限元仿真方法,对传动系统的参数进行辨识,根据传动系统拖动实验的转速波动信号,对理论分析的结果进行验证. 通过Matlab/Simulink的建模仿真以及650 kW海流能机组在低流速工况下的等效海试实验,验证了主动阻尼控制对传动系统的有良好的减振效果.
关键词:
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:
本文引用格式
李阳健, 何伟, 刘宏伟, 李伟, 林勇刚, 顾亚京.
LI Yang-jian, HE Wei, LIU Hong-wei, LI Wei, LIN Yong-gang, GU Ya-jing.
海流能由于高可预测性和高能量密度,具有较高的开发利用价值. 国际上最先进的海流能机组的单机额定功率已经达到2 MW,正在向商业化运营发展. 虽然海流能机组和风力机有着相似的工作原理,但是由于工作介质的不同,海流能机组的发展面临着更大的挑战[1].
本文以浙江大学650 kW海流能机组为研究对象,对传动链进行数学建模,结合计算机软件仿真和实验数据,对传动链的低阶扭转模态进行分析;结合海流能机组控制系统,设计主动阻尼加载控制器. 在软件仿真的基础上,开展初步的海试试验验证,结果表明,设计的阻尼控制器能够有效抑制传动链谐振载荷.
1. 传动链模态分析
1.1. 传动链模型
浙江大学650 kW水平轴海流能机组的机械部分主要由叶轮、两级行星齿轮和发电机组成,如图1所示.
图 1
为了研究机组传动链的固有频率,以便对低阶模态进行加阻控制,对该传动系统进行等效建模. 目前,有限元法和多体系统法是对复杂传动系统建模的2种主要方法. 有限元模型虽然能够获得很高的精度,但计算量很大,仿真时间很长;多体系统法是将一个机械系统划分成由若干个柔性弹簧阻尼器连接的刚体,选择合适的刚体数量能够更好地逼近原系统. 本文采用多体系统方法,建立传动链模型.
在风力发电机组的研究中广泛采用两质量弹簧阻尼模型,该方法不能很好地反映传动链的动态特性. 为了获得更加准确的特性,结合650 kW海流能机组齿轮箱的质量分布情况,采用四质量弹簧阻尼模型. 该等效系统的数学模型如下:
式中:i1和i2分别为两级行星齿轮的传动比,I1为叶片和轮毂的转动惯量,I2为一级传动的行星架的转动惯量,I′2为一级传动的太阳轮的转动惯量,I3为二级传动的行星架的转动惯量,I′3为二级传动的太阳轮的转动惯量,I4为发电机刹车盘等的转动惯量,ki为各轴的刚度,Tr为叶轮转矩,Tg为发电机转矩.
由于传动链处于弱阻尼状态,阻尼系数较小,且阻尼对模态分析结果没有影响,在对传动链进行模态分析的过程中忽略阻尼.
1.2. 传动链模态仿真分析
表 1 各个部分转动惯量
Tab.1
参数 | 数值 |
叶片和轮毂 | 54 151.73 |
一级传动的行星架 | 445.68 |
一级传动的太阳轮 | 404.35 |
二级传动的行星架 | 23.96 |
二级传动的太阳轮 | 19.59 |
发电机刹车盘等 | 280.50 |
图 2
图 2 650 kW海流能机组的传动链系统
Fig.2 Drive chain system of 650 kW tidal current energy generating set
传动链系统中每个部分轴的刚度Ksh可以通过下式计算得到:
式中:G为切变模量,I为截面极惯性矩,L为轴的长度.
由于实际传动链轴的直径不一,通过理论计算的刚度精度较低,无法满足控制需求. 利用有限元仿真软件ANSYS对各轴进行有限元仿真分析,得到施加单位力矩下轴的角位移形变,计算轴的刚度:
式中:
低速轴的有限元分析如图3所示.
图 3
依次对各轴进行有限元分析,得到各轴刚度,如表2所示.
表 2 系统中轴段刚度
Tab.2
参数 | 数值 |
低速轴k1 | 5.356×103 |
中间轴k2 | 4.205×103 |
高速轴k3 | 4.601 |
图 4
图 5
1.3. 传动链模态实验验证
图 6
图 7
从图7可知,得到频谱峰值对应的频率为2.44 Hz,与理论分析结果2.91 Hz相近,但是存在一定的误差. 这表明基于四质量块的传动链等效数学模型能够有效估计传动齿轮箱的低阶模态. 造成误差的可能原因如下:四质量块模型是基于集中质量的假设,忽略了轴的分布质量问题,存在一定的模型误差;忽略了联轴器的影响,使得计算得到的刚度比实际刚度大;忽略了齿轮箱和电机中的润滑油造成的附加质量,使得模型中的质量块转动惯量偏小;忽略了轴系一些小零部件的装配和齿轮的啮合刚度. 在更精准的模型中,可以考虑这些因素,得到更准确的共振频率.
2. 阻尼控制器设计
由于传动链的低阶共振频率较低,阻尼较低,特别是在恒转矩控制过程中,极易发生扭振,造成剧烈的扭矩波动,对传动链造成严重的有害扭振载荷. 本文采用施加电磁转矩的方式,给传动系统人为地增加虚拟阻尼,减小传动链的扭振. 主动阻尼控制原理较简单,通过在给定发电机转矩上附加一个转矩成分
式中:J、B、k分别为传动链相对于低速轴的等效转动惯量、等效阻尼和等效刚度. 使得
则有
即相当于通过电磁转矩控制,给传动链增加了一个额外的阻尼,达到减振的效果. 式中:Bg为电磁阻尼系数.
阻尼控制器须获取有害的转矩成分,即通过检测低速轴和高速轴的转矩完成,但由于安装扭矩传感器不切实际,可以通过对高速轴转速信号进行适当的滤波,获取传动链扭振的信息. 高速轴转速信号经过带通滤波器,得到有害的转速波动信号,并通过海流能机组转矩控制器,施加一个与传动链转速波动方向相反的补偿转矩
图 8
图8中,带通滤波器的传递函数为
式中:
由于转速传感器以及控制器的控制信号均为离散信号,且信号通讯、转矩控制及滤波器均会产生一定的延时[11](一般可以等效为一阶惯性环节),直接进行主动阻尼控制,极有可能因为相位滞后、系数选用不当、通讯延迟等原因,不能起到加阻效果,甚至引起振荡. 为了实现良好的控制效果,在滤波器滞后添加相位超前矫正:
3. 仿真与实验
3.1. 基于MATLAB/Simulink的主动阻尼控制仿真
图 9
图 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. 海试实验验证
图 13
图 14
与仿真中的无阻尼系统相比,实际系统引起振荡的过程更缓慢,可见系统本身存在一定的阻尼,但是不足以抑制传动链的振荡. 在海流能发电机组的控制系统中添加主动阻尼控制,根据仿真结果调节阻尼,得到的试验结果如图15所示.
图 15
4. 结 论
(1)基于多体系统法和有限元仿真,对海流能发电系统的传动系统进行模态分析. 所得的理论结果通过海流能发电机组的厂内拖动试验结果进行验证. 结果表明,采用的多体系统法具有一定的准确性,理论结果与试验结果相比,误差小于20%,对产生误差的原因进行定性分析.
(2)基于海流能发电系统的控制器,设计主动阻尼控制器;根据可能的系统相位滞后,设计超前矫正环节.
(3)在不考虑系统自有阻尼的理想条件下,搭建海流能发电机组主动阻尼控制的Matlab/Simulink仿真模型,研究不同阻尼系数下系统的减振效果. 仿真结果表明,选择合适的阻尼系数,主动阻尼控制可以有效地降低传动链谐振,与未加阻系统相比,传动链转速波动幅值减小69%以上;同时得到更加稳定的叶尖速比,使得机组效率提高约4.3%.
(4)开展650 kW海流能机组在低流速工况的恒转矩控制下的主动阻尼控制实验,实验验证了所设计的主动阻尼控制器对传动链的谐振具有很强的减振效果.
参考文献
Developments in large marine current turbine technologies: a review
[J].DOI:10.1016/j.rser.2016.12.113 [本文引用: 1]
Comparison of oscillation damping capability in three power control strategies for PMSG-based WECS
[J].
Wind turbine control for load reduction
[J].
风力发电机组传动链调谐惯量阻尼器研究
[J].
Research of tuned inertia damper for wind turbine drive rain
[J].
Active torque control for gearbox load reduction in a variable-speed wind turbine
[J].DOI:10.1109/TIA.2012.2227131 [本文引用: 2]
A review on the development of tidal current energy in China
[J].DOI:10.1016/j.rser.2010.11.042 [本文引用: 1]
Tidal current turbine based on hydraulic transmission system
[J].
Design and test of a 600 kW horizontal-axis tidal current turbine
[J].DOI:10.1016/j.energy.2019.05.154 [本文引用: 1]
潮流发电系统功率控制及并网运行
[J].
Study on power control and grid-connected operation for tidal current energy conversion system
[J].
Development and experiment of a 60 kW horizontal-axis marine current power system
[J].DOI:10.1016/j.energy.2015.04.018 [本文引用: 1]
Integrated design and implementation of 120 kW horizontal-axis tidal current energy conversion system
[J].DOI:10.1016/j.oceaneng.2018.04.017 [本文引用: 1]
/
〈 |
|
〉 |
