永磁调速器(permanent magnet coupler, PMC)是通过调节气隙长度控制输出转矩的一种调速节能装置, 采用机械式结构, 实现了非接触能量传递, 与传统的接触式传动机构相比, 高效永磁调速器具有维护方便、容忍对中误差、消除振动传递、避免电磁干扰、实现软启动和过载保护等优势, 已经成功应用于电力、石油化工和军事等诸多领域中[1].
高效永磁调速器主要靠铜盘上的涡流产生磁场进而实现转矩的传递, 传递过程中各个部件都会产生或多或少的损耗, 而损耗是衡量高效永磁调速器性能的主要指标, 损耗也直接影响着高效永磁调速器的工作效率.
目前针对高效永磁调速器各部件的损耗问题, 国内外的研究成果还较少, 且主要是对铜盘涡流损耗的研究.文献[2]建立了涡流损耗模型, 分析了永磁联轴器涡流损耗及其影响, 并进行了相关的理论计算, 但并未对本体各部分进行损耗精确计算; 文献[3]中建立了轴向永磁联轴器的三维解析模型, 考虑了设备的几何尺寸和材料特性, 利用三维有限元方法进行了验证; 文献[4]提出了一种开槽结构的二维解析模型, 分析了磁场分布情况, 与三维有限元方法进行了比较, 但并未对损耗进行详细的计算; 文献[5]主要是基于有限元模型分析永磁体和铜盘对损耗的影响因素并进行了优化设计, 但并未对其他部分进行损耗计算与分析, 缺乏系统的理论计算与研究.目前对永磁调速器的各部分损耗研究主要还是参考永磁同步电机相关文献, 例如永磁电机的损耗计算[6-8]、电机的损耗模型分析[9-11]、电机的电磁场与温度场分布[12].
针对以上研究现状, 本文对永磁调速器上各类电磁损耗进行理论分析, 建立了运动涡流场的三维有限元分析模型, 对1 MW永磁调速器的铜盘损耗、永磁体涡流损耗以及永磁体背衬钢盘对磁场及转矩大小影响进行了详细分析并与实验数据进行对比, 最后深入分析了分段永磁体对永磁体涡流损耗的影响.
1 机械结构与时步有限元模型本文以一台1 MW高效永磁调速器为损耗研究对象, 高效盘式永磁调速器的机械结构示意图如图 1所示.1 MW高效永磁调速器样机参数如表 1所示.各部分材料参数如表 2所示.
在ANSYS有限元软件的电磁场模块中, 利用时步有限元法对1 MW高效永磁调速器的损耗进行计算, 按表 1参数建立的1 MW高效永磁调速器全域时步三维有限元模型示意图如图 2所示.为了得到更精确的结果, 仿真中设有500个计算点.
为保证计算精度, 铜盘和永磁体进行手动剖分, 并考虑集肤效应.永磁调速器有限元模型共有剖分四面体202 535个, 其中永磁体分段四面体73 582个, 占总数的36.3%.
2 各部件损耗分析与计算 2.1 铜盘涡流损耗分析与计算铜盘的作用是传导感应的涡流, 当永磁调速器高速运行时, 涡流密度随时间不断变化, 涡流的频率同步变高, 此时铜盘上涡流引起的集肤效应不能忽略.导体集肤深度为
(1) |
式中:ω为激励源的角频率, 此处为铜盘中感应电动势变化的频率, 且fc=sfPM; fPM为永磁体磁场的交变频率; μ为磁导率; γ为电导率; s为转差率.
由于位移电流密度远小于传导电流密度, 可以忽略不计, 电磁场所满足方程组的微分形式为
(2) |
式中:H为磁场强度; J为电流面密度; E为电场强度; B为磁感应强度; D为电通密度; ρ为自由电荷体密度.
对于各向同性的媒质, 有
(3) |
磁通密度B可由磁矢位A表示:
(4) |
在永磁调速器中, 涡流密度表达式为
(5) |
式中:Je为涡流密度; φ为标量电位; v为导体运动速度, 常数.
采用A, φ-A法[13]建立三维动态电磁场数学模型, 得到涡流区域的控制方程为
(6) |
考虑集肤效应在体积V的铜盘中消耗的平均功率为[12]
(7) |
在相对速度一定的情况下, 整个铜盘的涡流总损耗在任意时刻都相同, 因此暂态电流面密度也可以当作有效值来直接计算铜盘涡流值[6].所以永磁调速器的输出功率P及传输扭矩T表达式为
(8) |
(9) |
式中n为原电动机转速.
对整个求解区域的外边界施加狄利克莱边界条件, 可计算出铜盘上的平均涡流密度.由有限元计算可得出, 铜盘上的涡流损耗为40 kW.当转差率为0.04时, 可以求出输出功率为950 kW, 输出转矩为6 500 N·m, 轴向力为-5 500 N.
永磁调速器实验平台如图 3所示, 在基座上从右往左依次为交流电动机、永磁调速器、扭矩/转速仪、交流电动机和直流电动机, 设备相互之间通过弹性联轴节进行连接.通过有限元分析软件得到的数据与实验数据相比较, 如表 3所示.可以看出, 仿真计算得到的数据与实验数据很接近, 进而验证了理论推导的正确性与有效性.
由于永磁体背衬钢盘的厚度对输出转矩有一定的影响, 所以本文重点研究永磁体的背衬钢盘对磁场的影响情况.对永磁调速器进行动态电磁仿真分析, 二维有限元仿真模型如图 4所示, 永磁体背衬钢盘厚度与转矩的关系数据如表 4所示.
从表 4中可以看出, 如果钢盘厚度过薄, 气隙磁通密度变小, 传递转矩的大小也会减小; 如果磁盘厚度太厚, 转动惯量变大, 会增加永磁调速器的起动转矩.本文选用了钢盘厚度为9.5 mm, 磁力线通过永磁体背衬钢盘的二维瞬态路径如图 5所示.
由图 5可知, 磁力线从永磁体的N极发出, 经过气隙、铜盘以及其后面的背衬钢盘, 到达相邻永磁体的S极, 最后通过永磁体的背衬钢盘构成闭合的主磁路[14].
2.3 永磁体涡流损耗分析与计算永磁调速器铝盘上的永磁体块在实际的损耗计算与分析中一般都是忽略不计的, 然而在后续的实际调查中发现, 工程上使用永磁调速器工作往往是在高速、大功率和一些工作环境非常恶劣的情况下, 这样就使得永磁体由于自身涡流损耗发热的影响, 加大了高温退磁的风险系数.
永磁体涡流损耗在数值上的计算可由式(7)计算得到.为了方便分析永磁调速器永磁体涡流损耗的分布情况, 本文对永磁调速器永磁体上的涡流分布进行了形象化分析, 建立了图 6所示的永磁体涡流计算模型.该模型中, 永磁体的长为l, 宽为w, 厚为h.本文研究坐标在(x, y)处的某回路(见图 6中的阴影).此回路x方向的单位宽度为dx, y方向的单位宽度为dy.
根据法拉第电磁感应定律[15]:
(10) |
由式(4), 得
(11) |
式中:Φ为涡流路径所包围的磁通量; S为涡流路径所包围的面积.
当永磁体涡流路径为图中阴影回路时, 由式(11)可得
(12) |
该计算模型中,
因此, 由式(11)可得
(13) |
由式(13)可知, 永磁体涡流密度的分布情况与永磁体的电导率、永磁体的形状、永磁体上磁通密度的变化率以及所处永磁体上的位置有关.
1 MW永磁调速器永磁体涡流密度分布如图 7所示.从图中可以看出, 永磁调速器上的永磁体块虽然直接决定着铜盘上生成涡流的大小, 但是永磁体块本身的磁通密度和涡流并不大, 且整体处于相对平均的状态, 每块永磁体的状态及场量值大小几乎接近.由有限元软件分析可得, 系统达到稳定状态时永磁体涡流总损耗为2.63 kW.
永磁体分段是降低永磁体涡流损耗的一个有效的途径, 具有较高的研究价值[9].在1 MW高效永磁调速器的永磁体盘上, 相邻永磁体之间存在些许的横向漏磁, 分段永磁体同样可以解决漏磁通的问题.结合前面的理论研究, 对分段处理永磁体如何能降低永磁体涡流损耗的原理进行分析.
3.1 等效电路由式(13)可知, 在其他条件不变的情况下, 当永磁体的电导率减小时, 永磁体内的涡流密度会相应变小.对永磁体进行分段处理后, 可以降低永磁体电导率[15], 分段永磁体径向分段等效电路如图 8所示.由本文的1 MW高效永磁调速器的永磁体涡流密度分布图可以看出, 涡流在径向的路径较长且涡流密度更高, 所以可以采用径向分段的方式增大径向方向的电阻值.对永磁体块进行分段处理实际上就是增大径向回路中的电阻值, 从而减少了永磁体中感生出的涡流面密度[9].
沿径向永磁体分段前后, 永磁体内涡流损耗关系可近似为
(14) |
式中:G为永磁体的段数; Ψ为垂直穿过永磁体的磁链值;W和L分别为单块永磁体的宽度和径向长度; RW, RL为单块永磁体轴向与径向的等效电阻值.
由式(14)可以得出如下结论:当永磁体分段数G大于1时, 永磁体内涡流损耗的情况与分段数成平方关系衰减.与不分段时相比, 伴随着分段数的增加, 永磁体内涡流损耗被抑制得相当明显, 且永磁调速器各部分的涡流损耗仍保持着各自相应规律不变.
3.2 结果分析为了更加直观地分析永磁体分段的作用, 本文对1 MW高效永磁调速器采用三维瞬态电磁场求解, 永磁体平均涡流损耗与永磁体分段关系如图 9所示.
从图 9可以看出, 径向永磁体分段可以有效降低永磁体内的涡流损耗.当永磁体分段数大于两段后, 永磁体内部涡流损耗减少量较为明显; 当多于一定段数以后, 涡流损耗减少量越来越小, 最终趋于稳定.很显然, 这一变化规律与式(14)是一致的.把1 MW高效永磁调速器铝盘上的每块永磁体沿径向平均分成几小段时, 实际上是增加了永磁体内径向路径的电阻, 有效减弱了永磁体内产生的涡流面密度, 全部永磁体的涡流损耗可以明显减少, 而且这并不会影响其输出扭矩的大小.
4 结语本文对1 MW高效永磁调速器的损耗机理进行较为详细的研究, 使用三维电磁场对永磁调速器各部件损耗的分布进行了仿真分析.对铜盘涡流和永磁体涡流进行了损耗计算, 对钢盘厚度对转矩和磁场的影响进行了分析, 进一步分析了分段永磁体对永磁体涡流损耗的影响情况.永磁体沿径向由一整块分段成几小块时, 可以显著减少永磁体上的涡流密度, 永磁体内的涡流损耗也会降低.本文的研究内容对永磁调速器的效率分析、散热设计以及损耗优化有很好的借鉴价值.
[1] | Lubin T, Mezani S, Rezzoug A. Experimental and theoretical analyses of axial magnetic coupling under steady-state and transient operation[J]. IEEE Transactions on Industrial Electronics, 2014, 61(8): 4356–4365. DOI:10.1109/TIE.2013.2266087 |
[2] | Zhu Z N, Meng Z. 3D analysis of eddy current loss in the permanent magnet coupling[J]. Review of Scientific Instruments, 2016, 87(7): 074701. DOI:10.1063/1.4954970 |
[3] | Mohammadi S, Mirsalim M, Vaez-Zadeh S, et al. Analytical modeling and analysis of axial-flux interior permanent-magnet couplers[J]. IEEE Transactions on Industrial Electronics, 2014, 61(11): 5940–5947. DOI:10.1109/TIE.2014.2311391 |
[4] | Dai X, Liang Q H, Cao J Y, et al. Analytical modeling of axial-flux permanent magnet eddy current couplings with a slotted conductor topology[J]. IEEE Transactions on Magnetics, 2016, 52(2): 8000315. |
[5] |
李桃, 林鹤云, 黄允凯, 等.
基于三维运动涡流场分析的永磁涡流联轴器特性[J]. 东南大学学报(自然科学版), 2010, 40(2): 301–305.
( Li Tao, Lin He-yun, Huang Yun-kai, et al. Characteristics study of permanent magnet eddy current coupling based on 3D moving eddy current field analysis[J]. Journal of Southeast University (Natural Science), 2010, 40(2): 301–305. ) |
[6] |
孔晓光, 王凤翔, 邢军强.
高速永磁电机的损耗计算与温度场分析[J]. 电工技术学报, 2012, 27(9): 166–172.
( Kong Xiao-guang, Wang Feng-xiang, Xing Jun-qiang. Losses calculation and temperature field analysis of high speed permanent magnet machines[J]. Transactions of China Electrotechnical Society, 2012, 27(9): 166–172. ) |
[7] | Yamazaki K, Abo A. Loss investigation of interiorpermanent-magnet motors considering carrier harmonics and magnet eddy currents[J]. IEEE Transactions on Industry Applications, 2009, 45(2): 659–665. DOI:10.1109/TIA.2009.2013550 |
[8] | Martin F, Zaim M, Tounzi A, et al. Improved analytical determination of eddy current losses in surface mounted permanent magnets of synchronous machine[J]. IEEE Transactions on Magnetics, 2014, 50(6): 1–10. |
[9] |
张磊, 高春侠.
永磁同步电机磁钢涡流损耗模型及其衡量指标[J]. 电机与控制学报, 2013, 17(7): 46–53.
( Zhang Lei, Gao Chun-xia. Eddy-current loss model and index for magnets of permanent magnet synchronous motors[J]. Electric Machines and Control, 2013, 17(7): 46–53. ) |
[10] | Wu L J, Zhu Z Q, Staton D, et al. Analytical modeling and analysis of open-circuit magnet loss in surface-mounted permanent magnet machines[J]. IEEE Transactions on Magnetics, 2012, 48(3): 1234–1247. DOI:10.1109/TMAG.2011.2157351 |
[11] | Wang J, Lin H Y, Fang S H. Analytical prediction of torque characteristics of eddy current couplings having a quasi-halbach magnet structure[J]. IEEE Transactions on Magnetics, 2016, 52(6): 1–9. |
[12] | Qiu H B, Yi R, Li W L, et al. Influence of rectifiers on high speed permanent magnet generator electromagnetic and temperature fields in distributed power generation systems[J]. IEEE Transactions on Energy Conversion, 2015, 30(2): 655–662. DOI:10.1109/TEC.2014.2366194 |
[13] | Bryant C F, Emson C R I, Trowbridge C W. A general purpose 3-D formulation for eddy currents using lorentzgauge[J]. IEEE Transactions on Magnetics, 1990, 26(5): 2373–2375. DOI:10.1109/20.104735 |
[14] |
王旭. 永磁调速器磁路设计与建模分析技术的研究[D]. 沈阳: 东北大学, 2013.
( Wang Xu.Research on adjustable permanent magnetic coupler magnetic circuit design and modeling analysis technique[D]. Shenyang:Northeastern University, 2013. http://cdmd.cnki.com.cn/Article/CDMD-10145-1016009121.htm ) |
[15] |
陈萍, 唐任远, 佟文明, 等.
高功率密度永磁同步电机永磁体涡流损耗分布规律及其影响[J]. 电工技术学报, 2015, 30(6): 1–9.
( Chen Ping, Tang Ren-yuan, Tong Wen-ming, et al. Permanent magnet eddy current loss and its influence of high power density permanent magnet synchronous motor[J]. Transactions of China Electrotechnical Society, 2015, 30(6): 1–9. ) |