张 静,常 臻,范舒羽
(1.金陵科技学院机电工程学院,江苏 南京 211169;2.盐城工学院机电工程学院,江苏 盐城 224051)
随着海上风能、波浪能等新型能源的不断开发和利用,多种漂浮式能源互补发电平台应用于深海区域。这种能源组合发电系统可以有效提高系统平台输出功率的稳定性,提升海洋能源的空间利用率,减少海洋发电系统平台架构建设成本[1-3]。相较于固定式风机平台,浮式风力机平台要承受风力机的空气动力载荷、海洋波浪和海流的水动力载荷、平台的系锚载荷,在上述环境载荷耦合作用下研究平台的动力学特性是较为复杂的问题。在实际工程应用中,通过设计平台减振装置实现平台的减振功能,将波浪有害的振动能量进行转换,既可以与发电平台能源互补又可以减小结构的振动,提升平台机构的平稳性与安全性[4-5]。
早在20世纪90年代,国外相继出现对筒形直线振荡电机的研究,应用领域主要为冰箱压缩机、斯特林制冷系统和斯特林引擎发电系统。如文献[6]提出了一种无轭部结构的直线振荡电机,其电磁推力最大可达250 N。文献[7]设计了一种复式直线振荡电机,对其齿部和轭部的硅钢片叠制工艺进行研究。在国内,对直线振荡电机的研究主要集中在高校和科研院所,比如浙江大学、华中科技大学、东南大学和沈阳工业大学等对横向磁通双定子永磁直线振荡电机、圆筒形直线振荡电机、单极动磁式直线振荡电机等的模型和电磁性能进行研究[8-12]。
直线振荡电机根据通入电源类型可分为交流振荡电机和直流振荡电机,或分为单相振荡电机和三相振荡电机。根据电机内运动机构的不同分为电枢运动型、磁场运动型等。其中,相比三相振荡电机,单相振荡电机的机电结构简单,运动行程短,响应迅速,控制方法简单可靠;相比电枢运动型振荡电机,单相运动型振荡电机的磁场体积小,电机密封与加工工艺更为简单。
基于上述分析,本文提出了一种用于浮式发电平台调谐减振系统的筒形单相永磁直线振荡电机。电机工作在发电状态时,研究其输出电能特性,实现波浪能与风能组合互补发电;电机工作在电动状态时,研究其振动特性,实现对发电平台的调谐减振功能,提高系统结构的稳定性和电能输出效率。首先,根据波浪载荷、电机的电磁负荷建立电磁耦合模型,分析电机内电磁场的分布情况,计算电机绕组自感、互感参数,以及绕组端输出电动势变化情况和电机内齿槽定力。其次,根据机械振荡原理搭建减振与发电系统振荡时域模型,分析电机振荡频率响应、动子振荡位移及电机电磁推力。最后,通过优化电机动子永磁体结构,提高输出电功率和转换效率。
在浮式平台减振系统中,振荡电机所带的负荷为波浪载荷,即波浪力通过中间传动机构(如浮筒或推板)传导到电机的动子结构部分,此过程的运动物理模型可用带有阻尼的振荡模型来分析,如图1所示。
图1 电机时域分析振荡模型
假定系统振荡过程中机械摩擦力为零,且振荡条件充足,根据牛顿第二定律可得系统的时域振荡方程如式(1)所示。
ma=f(t)-kx-cv
(1)
其中,m为系统质量,a为振荡运动加速度,k为弹性系数,x为系统振荡位移,c为系统运动阻尼系数,v为系统振荡速度。f(t)为系统中随时间变化的振荡物理量,当研究振荡电机的电动状态时,f(t)=0,即此时为自由振荡系统;当研究振荡电机的发电状态时,f(t)为波浪力的垂直分量,计算如式(2)所示。
f(t)=Fsin(ωt+φ)
(2)
其中,F为波浪振荡幅值,ω为波浪力函数角频率,φ为初相位。
将式(1)转换为微分方程,如式(3)所示。
(3)
经过对电机外形特点的分析,结合浮式平台的波浪环境以降低模型的计算复杂性,本文提出一种圆筒形单相永磁直线振荡电机,具体三维结构如图2所示。电机采用磁场运动型结构,即定子内采用半开口槽结构,槽内放置饼式绕组线圈,动子铁心外侧采用表贴式永磁体环结构。
图2 直线振荡电机结构模型
电机定子槽内的4个饼式绕组正反相接形成单相绕组,其振荡电路的振荡频率和机械振荡频率相同,其电动状态的控制方法更简捷。根据虚功法计算电机动子上产生的电磁推力,如式(4)所示。
(4)
其中,W(i)为电机内磁场能。永磁体(PM)表贴式单相振荡电机的电磁推力计算如式(5)所示。
(5)
其中,if为永磁体等效励磁电流,ia为电枢电流,Lf为电机内自感,Laf为电机内互感。
根据上述分析得到电机的基本结构模型,表1为单相永磁振荡电机主要结构的设计参数。
表1 电机的主要设计参数
对本文设计的直线振荡电机采用有限元法R-Z轴模型进行分析,在电机空载运行时间t=0.2 s时,分析得到电机的磁力线分布,如图3所示。电机轴向端部存在漏磁,这是由直线电机的端部效应引起的。研究发现当动子运行到不同位置时,电机内磁场分布随着永磁体位置的变化而发生改变,而交联在定子电枢绕组上的总磁通也会发生周期性的改变,实现电能和机械能的转换。
图3 电机空载时磁力线分布(t=0.2 s时)
根据有限元法可计算出该直线振荡电机定子绕组的自感,计算如式(6)所示。
(6)
其中,L为电感,I为对应的绕组电流,ψ为绕组线圈单位长度内交联的磁链,计算如式(7)所示。
(7)
其中,N为绕组匝数,Φ为绕组对应的磁通。
根据本电机的模型采用有限元法计算出该电机绕组的自感,具体如图4所示,可以得到自感的平均值为23 mH。
图4 电机绕组的自感曲线
根据直线振荡电机工作原理可知,当电机工作在电动状态时,电机电枢绕组中加入不同电流,得到的振荡推力有所不同。本文通过计算得出电机输入电流为i=3sin(10πt) A时的电磁推力变化情况,如图5所示。需要注意的是此时计算的电磁推力为电机电流作用下的电磁力和空载时定位力的合力,最大幅值约为102.70 N。
图5 电机电磁推力变化情况
此外根据分析计算发现,在电流频率不变条件下,随着输入电流幅值的不断增加,电磁推力的幅值也逐渐增大(如表2所示)。
表2 电磁推力与输入电流变化情况
根据电机的设计指标要求,电机工作频率为5 Hz,即振荡的周期为200 ms,在输入电流为i=3sin(10πt) A时,分析得出电机动子随时间的振荡位移情况,如图6所示,振荡的幅值范围为11 mm。
图6 电机振荡位移情况
当电机工作在发电机状态时,由波浪力带动电机动子运动进行发电。发电机绕组端空载电动势是振荡电机的重要指标之一,实现与发电平台其他能源互补发电。分析时,根据海况波浪平均运动速度设定电机动子的振荡速度为v=150sin(10πt) mm·s-1,得出电机空载运行时的电枢电动势,电动势幅值为12 V左右。电机在发电状态时产生空载电动势的波形中含有谐波,波形的总谐波畸变率(THD)是发电机输出电能质量的一个指标,是指各次谐波有效值平方和的方根与基波有效值比值的百分数,计算如式(8)所示。
(8)
其中,U1为基波电压有效值,Un为电压n次谐波有效值。通过傅里叶分解法得出本电机空载电动势谐波分布,如图7所示。通过分析发现波形主要含有5次谐波,这是由电机自身齿槽结构(如槽距和极距比值的不同)产生的,因此可以通过改变齿槽结构和永磁体极距来进行优化。
图7 振荡电机空载电动势波形谐波含量
电机模型的结构参数在最初设定后需要反复优化以使电机获得最优的性能。根据上述分析可知,本文设计的永磁直线振荡电机空载时的磁场是由电机内的永磁体产生的,故永磁体的结构参数将影响电机的空载电动势波形。本电机采用永磁体环阵列结构,在动子运动行程(10 mm)和动子轴向长度(70 mm)范围内分别对不同厚度、宽度的永磁体环作用下的电机空载电动势进行计算,结果如图8所示。由图8可知,当永磁体(PM)厚度和宽度增大时,电动势的幅值也增大,但增大到一定值后增加的幅度趋于缓慢,这是由电机内磁场趋于饱和引起的。考虑到电机的磁场饱和以及永磁体材料的加工工艺与费用,本电机采用厚度为3 mm、宽度为16 mm的永磁体环。
(a)电机动子结构
本文提出一种用于浮式发电平台的减振与互补发电的筒形单相永磁直线振荡电机的设计方法,此种电机在电动状态时进行减振,在发电状态时实现与平台的互补发电;采用有限元分析法对电机空载时的电磁场进行计算,并对电机工作在电动状态和发电状态下的相关特性进行分析,最后对电机内的永磁体环结构进行优化。分析与计算的结果表明:1)电机采用筒形单相绕组结构可以有效提升电机内磁场的磁通密度,减小电机的体积与材料费用,提升电机单位体积的能量密度与发电效率。2)电机在振荡工作状态下,随着输入正弦电流幅值的不断增大,电机的电磁推力也不断增大,考虑电机的温升,额定的输入电流为i=3sin(10πt) A,电磁的振荡位移为11 mm。3)合理优化电机内永磁体环结构,可以在小行程振荡范围内提高电机输出电压,降低电压波形谐波含量,提高输出电能质量。
上述分析表明,该直线振荡电机结构模型的工作原理、设计理念和优化方法适用于小功率直线振荡电机的设计与分析领域。