基于模拟退火算法的Halbach直线发电机优化设计

刘 娜 谭亦旻 莫伟强 韩欢庆 李 琳

基于模拟退火算法的Halbach直线发电机优化设计

刘 娜1谭亦旻2莫伟强1韩欢庆1李 琳1

（1. 北京理工大学珠海学院工业自动化学院 珠海 519088 2. Sea Electric Energy Incorporation Mississauga L5N 1L2）

在研究海洋波浪能量收集的过程中，需要建立相应的理论模型来降低工程成本，提高研发效率。针对波能转换器在收集波浪能时能量的转换效率不高这一问题，采用无槽Halbach直线发电机作为装置的二级能量转换机构来优化磁场分布，提高能量的转换效率。该文描述Halbach直线发电机的拓扑结构，引入磁矢量势能理论推导出发电机性能的表达式，由此建立直线发电机的综合分析模型。基于此模型，使用模拟退火算法求全局最优解，得到一组直线发电机设计参数，包括永磁体尺寸和绕线线圈的数据。计算结果表明，理论模型解析解的峰值功率与模拟退火算法的结果相差约为3.6%。基于优化后的参数，确定了Halbach直线发电机的设计，为样机实验的验证提供了数据支撑。

波能转换器 Halbach直线发电机 模拟退火算法 优化设计

0 引言

随着世界经济的发展，人口的增加，社会对能源的需求日益增长，各国相继开展了新能源领域的研究。海洋波浪能由于具有能量密度高、分布面广并且取之不尽的特点，受到了越来越多的关注[1]。英国建成了奥克尼波浪发电试验场，完成了功率达2MW的波浪发电设备的研发；苏格兰的一家公司设计了独特的筏式波浪能转换器海蛇，发电量接近1MW；瑞典某公司设计了一种单体式直线发电机垂荡浮子装置，可实现20kW的功率输出；美国某公司研发的双体式浮子装置，可由两个振荡浮子的相对运动产生电能[2-4]。我国对海洋波浪能发电的研究主要集中于固定式和漂浮式波能转换装置[5]。装机容量3kW的珠海市万山岛电站试发电成功后，又相继建成了20kW岸式波浪电站、8kW摆式波浪电站、100kW岸式振荡水柱电站以及30kW摆式波浪电 站[6-8]。目前海洋波浪能源的研发方向主要集中在大型的海洋波浪发电工程，在小型海洋能源设备方面投入的研究较少[9]。本文提出了一种无槽式Halbach直线发电机的理论模型来分析波浪能到电能的转换，用于小型海洋能源设备的研究，致力于为海洋浮标监测系统的建立提供能量。

在进行波浪能量转换器的设计时，一些关键参数需要进行优化。研究发现，以下参数对发电机的性能有直接的影响，即电机中永磁体的长度和厚度、中心轴的半径、空气隙的宽度、绕组导线的直径以及海浪的速度[10-11]。用来解决若干参数的组合优化常用的全局优化算法有遗传算法、蚁群优化算法和模拟退火算法[12]。文献[13]采用遗传算法解决Halbach直线发电机优化设计问题，对电机的关键变量优化后，电机性能得到提升。文献[14]使用蚁群算法解决MW级风力发电机的优化设计，该发电机为Halbach永磁体同步发电机。为了提高功率密度，M. Markovic采用名为ProDesign的软件进行优化设计[15]，但该软件仅返回一个输出参数。S. Kirkpatrick等将退火思想引入组合优化领域，对组合优化问题和统计力学中的热力学平衡问题进行类比，模拟退火过程以获得近似全局最优解来解决大规模组合优化问题[16]。与遗传算法和蚁群算法相比，模拟退火算法能够在优化的精度和效率之间达到一个平衡。在求解组合优化问题时，模拟退火算法能概率性地跳出局部最优解并最终趋于全局最优。所以本设计将采用模拟退火算法对Halbach直线发电机的关键参数进行优化。

1 波能转换器简介

捕获海洋波浪能有多种方式，常见的有激荡水柱式、液压传动式和直线发电机直驱式[17-18]。三种波能转换器的工作原理如图1所示，其中，直线发电机直驱式利用发电装置不同部件间的相对运动实现发电。这种直接驱动的能量转换方式减少了能量转换的中间环节，降低了过程中的能量损耗，能够提高能量转换效率。本设计将采用Halbach直线发电机作为二级能量转换机构。

图1 波能转换器的工作原理

2 Halbach直线发电机模型

2.1 拓扑结构

图2 磁场分布的有限元计算结果

Halbach直线发电机的拓扑结构如图3所示。图中，s为直线电机中心轴半径，p为中心轴半径与永磁体宽度之和，i为直线电机的半径（不包含硅钢壳的厚度），为空气隙的宽度，为单相线圈的高度。

图3 Halbach直线发电机的拓扑结构

2.2 综合分析模型

根据法拉第电磁感应定律，系统产生的感应电动势应为

式中，M为被磁化的磁感应强度；下标1、2、3分别为所对应的空气隙区域、永磁体区域和电枢区域。

通过傅里叶变换可求得Halbach阵列区域被磁化的磁感应强度为

在区域的边界处，连续性应满足

假设电枢的磁导率趋于无穷，则空气隙区域和电枢区域的磁通密度为0。联立边界条件和区域方程组求解，可得对应区域的磁通密度为

其中

式中，I和K分别为第一类和第二类修正贝塞尔函数；为贝塞尔函数的系数，则发电机的每相线圈的磁链和反作用力分别为

式中，为海浪的速度。假设线圈的感应电动势有效值为，那么发电机的输出电压和输出功率分别为

3 Halbach直线发电机模型的研究

3.1 线圈的铜线直径对Halbach直线发电机输出功率的影响

通常认为线圈中铜线的直径将对直线发电机的输出电压和功率有直接影响。原因在于在相同的绕组空间中，选择的线规不同，对应的线圈匝数将不同。线圈匝数的增加不仅会增大系统内部阻抗，也会产生更大的感应电动势，因此有必要找到一个最佳的铜线直径，通常使用美国线规（American Wire Gauge, AWG）使电机功率输出最大化。

由式（13）可知，电机的输出功率与铜线的直径无关，即与AWG的选择无关。但不可否认的是，绕线线圈的线径对直线电机所能承受的电流有着重要的影响，这对于样机的验证实验至关重要。

3.2 永磁体厚度对Halbach直线发电机输出功率的影响

图4 Halbach直线发电机的输出功率与永磁体厚度之间的关系

3.3 中心轴半径对Halbach直线发电机输出功率的影响

Halbach直线发电机的输出功率与中心轴半径之间的关系如图5所示。由图5可知，当t、m、air和AWG为指定常数时，即t=20mm、m=20mm、air=26mm、AWG=22，随着s的增加，直线电机的峰值功率亦呈单调递增趋势。

图5 Halbach直线发电机的输出功率与中心轴半径之间的关系

3.4 Halbach直线发电机理论模型的磁通密度分布

为了验证直线电机的理论模型，对永磁体阵列的磁通密度分布分别进行有限元的仿真计算和解析计算，如图6～图8所示。

图6 PM阵列在轴上的磁通密度分布

图7 PM阵列在轴上的磁通密度分布

图8 Halbach直线发电机PM阵列的磁场分布

4 模拟退火算法的应用

4.1 模拟退火算法流程的确定

图9 SAA流程

4.2 模型的验证

Halbach直线发电机的峰值功率如图10所示。由图10可知，Halbach直线发电机的峰值功率为508W。

图10 Halbach直线发电机的峰值功率

SAA算法在Matlab环境下进行，目标函数的执行次数大约为两万次，最优解被接受的次数约为8 500次，恶化解被接受的次数约为5 440次。随着程序的运行，接受恶化解的概率逐步下降，恶化解被接受的次数随之减少。SAA算法结果如图11所示。由图11可知，Halbach直线发电机最大输出功率为527.2W。因此，理论模型的解析解与SAA算法的数值相差19.2W，约为3.6%。二者的误差在可接受的范围内。图11中，运行次数为19 600，最大功率为527.2W，空气隙宽度为6mm，AWG=22，t=390mm，m=24mm，s=20mm。

图11 SAA算法结果

对应于全局最优解，发电机参数的设计尺寸见表1。

表1 Halbach直线发电机的设计参数

Tab.1 Design parameters of Halbach linear generator （单位: mm）

5 实验验证

Halbach直线发电机实验平台如图12所示。永磁体阵列由绕在铝制中心轴上的环形磁铁组成，其内径为40mm，外径为80mm。三相线圈由硅钢包裹，每相线圈的内径为92mm，外径为132mm，高度为67mm。通过线性电动推杆模拟海浪的作用力，其速度可达1.2m/s。样机中位移传感器用于测量相对运动并向电动推杆提供反馈信号，反作用力的大小通过力传感器测得。

图12 Halbach直线发电机实验平台

当负载与线圈的阻抗匹配，系统速度为0.53m/s时，输出的平均峰值功率为142W；系统速度为1m/s时，预计总输出功率为570W。实验平台的等效阻尼系数为1 102.8N/m，阻尼系数的计算结果为1 103N/m，这表明所建立的理论模型有效地描述了动力输出系统的性能，可用于结构优化。

6 结论

本文描述了Halbach直线发电机的拓扑结构，构建了其理论模型，基于该模型计算出电机的最大输出功率为508W。选择了几个对发电机的磁通密度分布产生影响的参数进行分析，得到了这些参数和电机输出功率的关系。使用SAA得到了一组电机的设计参数及其最大输出功率，功率可达527.2W。利用得到的电机参数制作了样机，考虑到安全因素进行了部分标称速度的验证。本设计中Halbach直线发电机气隙磁通密度的理论值为1.4T，有限元计算的结果为1.39T，二者基本匹配。值得注意的是，过大的磁场密度会导致电机的过度饱和，因此在样机的制作过程中要注意永磁体的用量和尺寸。

[1] 游亚戈, 李伟, 刘伟民, 等. 海洋能发电技术的发展现状与前景[J]. 电力系统自动化, 2010, 34(14): 1-12.

You Yage, Li Wei, Liu Weimin, et al. Development status and perspective of marine energy conversion systems[J]. Automation of Electric Power Systems, 2010, 34(14): 1-12.

[2] Shami E A, Zhang Ran, Wang Xu, et al. Point absorber wave energy harvesters: a review of recent developments[J]. Energies, 2018, 12(1): 47.

[3] 洪岳, 潘剑飞, 刘云, 等. 直驱波浪能发电系统综述[J]. 中国电机工程学报, 2019, 39(7): 1886-1900.

Hong Yue, Pan Jianfei, Liu Yun, et al. A review on linear generator based wave energy conversion systems[J]. Proceedings of the CSEE, 2019, 39(7): 1886-1900.

[4] Henderson R. Design, simulation, and testing of a novel hydraulic power take-off system for the Pelamis wave energy converter[J]. Renewable Energy, 2006, 31(2): 271-283.

[5] 方红伟, 宋如楠, 姜茹, 等. 振荡浮子式波浪能转换装置的全电气化模拟研究[J]. 电工技术学报, 2019, 34(14): 3059-3065.

Fang Hongwei, Song Runan, Jiang Ru, et al. An oscillating buoy type wave energy converter with all- electrical analogue method[J]. Transactions of China Electrotechnical Society, 2019, 34(14): 3059-3065.

[6] 王鑫, 李大鸣, 李雪临, 等. 大万山波浪能试验场泊位间水动力影响分析[J]. 海洋技术学报, 2015, 34(4): 69-73.

Wang Xin, Li Daming, Li Xuelin, et al. Analysis on the hydrodynamic impact among the berths in the dawanshan wave energy test site[J]. Journal of Ocean Technology, 2015, 34(4): 69-73.

[7] 鄂世举, 金建华, 曹建波, 等. 波浪能捕获及发电装置研究进展与技术分析[J]. 机电工程, 2016, 33(12): 1421-1427.

E Shiju, Jin Jianhua, Cao Jianbo, et al. Development of permanent magnet linear synchronous motors structure and research[J]. Journal of Mechanical & Electrical Engineering, 2016, 33(12): 1421-1427.

[8] Engstrom J, Kurupath V, Isberg J, et al. A resonant two body system for a point absorbing wave energy converter with direct-driven linear generator[J]. Journal of Applied Physics, 2011, 110(12): 124904.

[9] 张静, 余海涛, 施振川. 一种波浪发电装置用低速双动子永磁直线电机运行机理研究[J]. 电工技术学报, 2018, 33(19): 4553-4562.

Zhang Jing, Yu Haitao, Shi Zhenchuan. Research on a tubularl linear permanent magnet machines with dual translators for low speed wave energy con- version[J]. Transactions of China Electrotechnical Society, 2018, 33(19): 4553-4562.

[10] 张春雷, 张辉, 叶佩青, 等. 两相圆筒型永磁同步直线电机无传感算法[J]. 电工技术学报, 2019, 34(23): 4901-4908.

Zhang Chunlei, Zhang Hui, Ye Peiqing, et al. Research on sensorless algorithm of two-phase tubular permanent magnet synchronous linear motor[J]. Transactions of China Electrotechnical Society, 2019, 34(23): 4901-4908.

[11] Tan Yimin, Lin Kejian, Zu J W. Analytical modelling of Halbach linear generator incorporating pole shifting and piece-wise spring for ocean wave energy harvesting[J]. AIP Advances, 2018, 8(5): 056615.

[12] 倪有源, 刘跃斌, 王群京. 组合磁极结构的永磁电机解析法建模与分析[J]. 电机与控制学报, 2019, 23(1): 52-62.

Ni Youyuan, Liu Yuebin, Wang Qunjing. Analytical modeling and analysis of permanent magnet machine with combined magnetic poles[J]. Electric Machines and Control, 2019, 23(1): 52-62.

[13] Eigebaly A E, Elnemr M K. Optimal design of slotless PM Halbach array linear generator for wave energy converters at maximum power transfer condition[C]//2015 IEEE 8th GCC Conference and Exhibition, Beijing, 2015: 49-56.

[14] Alshibani S. Application of particle swarm optimi- zation in the design of Halbach permanent magnet synchronous generators for megawatt level wind turbines[C]//2018 7th International Conference on Renewable Energy Research and Applications (ICRERA), Beijing, 2018: 456-463.

[15] Markovic M, Perriard Y. Optimization design of a segmented Halbach permanent-magnet motor using an analytical model[J]. IEEE Transactions on Magnetics, 2009, 45(7): 2955-2960.

[16] Kirkpatrick S, Gelatt C D, Vecchi M P. Optimization by simulated annealing[J]. Science, 1983, 220(4598): 671-680.

[17] 肖曦, 摆念宗, 康庆, 等. 波浪发电系统发展及直驱式波浪发电系统研究综述[J]. 电工技术学报, 2014, 29(3): 1-11.

Xiao Xi, Bai Nianzong, Kang Qing, et al. A review of the development of wave power system and the research on direct-drive wave power system[J]. Transactions of China Electrotechnical Society, 2014, 29(3): 1-11.

[18] 卢琴芬, 沈燚明, 叶云岳. 永磁直线电动机结构及研究发展综述[J]. 中国电机工程学报, 2019, 39(9): 2575-2588.

Lu Qinfen, Shen Yiming, Ye Yunyue. Development of permanent magnet linear synchronous motors structure and research[J]. Proceedings of the CSEE, 2019, 39(9): 2575-2588.

[19] Younesian D, Alam M. Multi-stable mechanisms for high-efficiency and broadband ocean wave energy harvesting[J]. Applied Energy, 2017, 197: 292-302.

[20] 许孝卓, 孙震, 汪旭东, 等. Halbach交替极永磁同步直线电机特性分析[J]. 电工技术学报, 2019, 34(9): 1825-1833.

Xu Xiaozhuo, Sun Zhen, Wang Xudong, et al. Characteristic of a novel permanent magnet linear synchronous motor with Halbach array consequent- pole[J]. Transactions of China Electrotechnical Society, 2019, 34(9): 1825-1833.

[21] 倪有源, 崔征山. 组合磁化无槽永磁电机磁场分析与优化研究[J]. 电机与控制学报, 2020, 24(3): 79-87.

Ni Youyuan, Cui Zhengshan. Magnetic field analysis and optimization of slotless permanent magnet machines with combined magnetization[J]. Electric Machines and Control, 2020, 24(3): 79-87.

[22] Johnson D S, Aragon C, Mcgeoch L A, et al. Optimization by simulated annealing: an experi- mental evaluation. part I, graph partitioning[J]. Operations Research, 1989, 37(6): 865-892.

[23] Chams M, Hertz A, De Werra D, et al. Some experiments with simulated annealing for coloring graphs[J]. European Journal of Operational Research, 1987, 32(2): 260-266.

[24] Otten R H, Van Ginneken L P. Stop criteria in simulated annealing[C]//IEEE International Conference on Computer Design: VlSI in Computers Society, Los Alamitos, 1988: 549-552.

[25] 卢宇婷, 林禹攸, 彭乔姿, 等. 模拟退火算法改进综述及参数探究[J]. 大学数学, 2015, 31(6): 96-103.

Lu Yuting, Lin Yuyou, Peng Qiaozi, et al. A review of improvement and research on parameters of simulated annealing algorithm[J]. College Mathematics, 2015, 31(6): 96-103.

[26] 康立山, 谢云, 尤矢勇, 等. 非数值并行算法[M]. 北京: 科学出版社, 2003.

Optimization Design of Halbach Linear Generator with Simulated Annealing Algorithm

12111

（1. School of Industrial Automation Beijing Institute of Technology of Zhuhai Zhuhai 519088 China 2. Sea Electric Energy Incorporation Mississauga L5N 1L2 Canada）

It is necessary to establish a corresponding theoretical model to save engineering cost and improve research efficiency during studying energy harvesting of ocean wave. Because of low energy conversion efficiency existed in capturing ocean wave energy by wave energy converter, the mechanism of slotless Halbach linear generator was adopted as the secondary energy conversion device to optimize magnetic field distribution of the generator in this paper. The topology of Halbach linear generator was defined, then the magnetic vector potential theory was adopted to deduce the expressions of Halbach linear generator’s performance. Finally, a comprehensive analysis model of the generator has been established. With the derived model, a global optimal solution was gained by simulated annealing algorithm (SAA). Meanwhile, the linear generator design parameters have been obtained, including dimensions of permanent magnet and winding coil. The difference between the linear generator peak power of the analytical solution and SAA result is about 3.6%. Based on the optimized parameters, the design of Halbach linear generator is determined, which provides data support for the verification of the prototype experiment.

Wave energy converter, Halbach linear generator, simulated annealing algorithm, design optimization

TM359.4

10.19595/j.cnki.1000-6753.tces.200442

广东高校省级重点平台和重大科研项目（2018KQNCX344，2018KQNCX345）、广东高校省级重点平台和重大科研项目（2019KQNCX196）资助。

2020-05-05

2020-08-01

刘 娜 女，1984年生，硕士，讲师，研究方向为振动能量的收集、直线电机。E-mail: 43736176@qq.com（通信作者）

谭亦旻 男，1989年生，博士，研究方向为非线性振动、智能材料应用等。E-mail: yimin.tan@hotmail.com

（编辑 陈 诚）