研究与设计EMCA 迫札与才望制应用2013,40(9) 基于温度场计算的开关磁阻电机冷却系统设计术 宋受俊 , 张 蔓 , 吴晏辰 (1.西北工业大学自动化学院,陕西西安摘710072;2.西北工业大学航海学院,陕西西安710072) 要:针对前期电磁设计得到的开关磁阻电机方案,通过对平均磁密波形的傅里叶分解,计算了电机 的铁耗和铜耗。建立了开关磁阻电机详细的集总参数热网络模型,并得到了电机各主要部位的温升。在此基 础上,设计了一个定、转子双水冷系统。通过对施加冷却系统前后电机主要部位温升的比较,证明了所设计水 冷系统的有效性。 关键词:开关磁阻电机;冷却系统;损耗;温升;热网络 中图分类号:TM 302 文献标志码:A 文章编号:1673—6540(2013)09-0006-06 Cooling System Design for Switched Reluctance Motor Based on Analytic Calculation of Thermal Field SONG Shoujun , ZHANG Man . WU Yanchen (1.School of Automation,Northwestern Polytechnieal University,Xi’an 7 1 0072,China; 2.School of Marine Technology,Northwestern Polytechnical University,Xi’an 7 1 0072,China) Abstract:Aeeording to the initial scheme of switched reluctance machine(SRM)obtained by electromagnetic design,the iron losses and copper losses were calculated by Fourier decomposition of its average flux density waveform.Detailed lumped 1)arameter thermal network model of the considered SRM was built,and temperature rises of some parts were got.On this basis,a stator and rotor dual water cooling system was designed.The validity of the designed cooling system was verified by comparison among temperature rises of some pans in the machine before and after the application of cooling system. Key words:switched reluctance motor(SRM);cooling system;losses;temperature rise;thermal network 算方法主要分为有限元法(Finite Element Method,FEM)和集总参数热网络法(Lumped 开关磁阻电机(Switched Reluctance Motor, Parameter Thermal Network,LPTN)两类。FE 。。 SRM)具有起动转矩大、调速范围宽、控制灵活、适 能够准确得到电机内部的温度场分布及各部分的 温升值,但该方法计算时间长,数据量大,灵活性 差,专业水平要求较高,不便于工程应用 LPTN 能够在保证相当精度的同时,有效提高计算速度, 且具有较强的灵活性。在电机设计阶段,使用 LPTN可大大缩短设计周期,提高研发效率。另 外,LPTN更适用于对冷却系统相关尺寸的迭代 优化。 应恶劣环境及成本较低等优良性能,具备与其他 类型电机相竞争的潜在优势¨j,在多电飞机 、 电动汽车【 J、风力发电 等军民用领域具有广阔 的应片j空间。SRM还具有功率密度高,集中式绕 组散热面积小,热负荷高的特点 J,致使其温升 问题突出。过高的温升会影响电机的运行性能、 可靠性和使用寿命。因此,在SRM设计过程中对 其温度场进行分析很重要。 目前,电机设计阶段的温度场分析及温升计 本文从SRM初始电磁设计方案出发,详细计 算了铁耗及铜耗。利用集总参数热网络模型,得 基金项目:国家自然科学基金(51107100);陕西省自然科学基金(2011GQ7001);教育部博士点基金新教师类 (20116102120033);两北T、I 大学基础研究基金(JC201119) 6—— 迫扎-号}控制应用2013,40(9) 到了电机各关键部位的温升。结果表明,该设计 方案的温升极高。因此,在不影响电机电磁特性 的前提下,设计了一个定、转子双水冷系统。水冷 系统的引人,大大改善了电机的散热性能,各部位 温升均降至合理的数值。 1 SRM初始设计方案 在前期电磁设计过程中,所提出的主要性能 指标要求如下: 结构形式为3相,6/4极;电压Us=270 V DC;额定转速n=27 000 r/min;额定功率P : 30 kW;额定效率”=80%。 针对性能要求和约束,利用基于输出方程的 传统方法,得到了SRM的初始设计方案。图1为 SRM定转子冲片结构示意图,表1给出了初始方 案中各尺寸的数值。 图1 SRM定、转子冲片结构示意图 表1初始设计SRM几何尺寸数值 参数 取值 参数 取值 , 74 Y 12.5 d 24.2 g 0.3 r 37 d 9.7 Y 13.3 r。h 14 t 19.3 t 20.4 性能核算结果表明,上述初始方案的额定功率 为26.3 kW,额定效率为81%,基本满足了设计需求。 初始设计仅考虑了电机的电磁特性,未顾及 其散热性能。电机的热源主要来自于运行过程中 产生的损耗。因此,在对SRM进行热分析之前, 有必要首先对其损耗进行计算。 2 SRM损耗计算 对于SRM而言,其定子中的损耗主要包括由 研究与设计 EMCA 焦耳效应引起的铜耗,以及由磁滞和涡流引起的 铁耗,在转子中,仅存在铁耗。 2.1铁耗计算 电机的铁耗与铁心中的磁密以及电机的换相 频率成正比,在正弦磁通条件下,其近似计算公式 可由斯坦梅茨(Steinmetz)方程得到: Pfe=P +Ph=cfB +ch (1) 式中:P ——涡流损耗; P ——磁滞损耗; C ——涡流损耗系数; C ——磁滞损耗系数; .卜电机的换相频率; B ——铁心磁密的最大值; 凡——指数。 c 、c 和n均可通过所使用硅钢材料M19在 50 Hz及400 Hz下的相关测量数据得到_8 J。 实际上,SRM铁心中的磁通多为非正弦的畸 变磁通波形,因此求解铁耗时必须利用傅里叶分 解将磁通波形转化成不同频率的正弦磁通之和, 再利用式(1)进行计算。 在本文中,首先通过FEM得到定子各极的磁 通波形,如图2所示,然后运用文献[9-11]中给出 的方法得到定、转子其他部分的磁通波形,进而得 到电机各部分的磁密波形。 图2 SRM足f极磁通坡形 对得到的定、转子铁心磁密波形进行傅里叶 分解: B( )=了50+∑[amCOS(m,olt)+ bmsin(//Z/O1t)] (2) 式中:甜 为基波角频率,而待定系数 。…口 ,b。… 6 的求解方法可参见文献[9]。 ——7—— 研究与设计冀 蠡 电札与板书J应用2013,40(9) 在得到上述待定系数后,即可很方便地求解 磁通密度各次谐波的幅值: 3 温度场计算及冷却系统设计 由于表2所示损耗是SRM稳定运行于额定 状态下的平均值,且电机的温升时间常数较大,其 温度场可视为稳态场。因此,本文采用了稳态热 B = ̄/n +6 (3) 对于磁滞损耗,采用J.D.Layers等人提出的 考虑了小滞环的磁滞损耗计算数学模型 。 进行 计算。 网络模型,亦即取消了热容。 3.1初始设计方案温度场计算 本文基于文献[13—14]提出的方法,建立 SRM初始设计方案的稳态热网络模型,建模前需 做出如下假设: 局部磁滞回线因子: K =1+ ・∑AB (4) 式中: ——常数,0.6≤ ≤0.7,本文取0.65; 磁密波形的峰值; (1)轴向与径向热阻相互独立,互不影响; (2)热生成率均匀分布; (3)每一个元件的温度都为平均温度; (4)圆周的热流忽略不计。 3.1.1等效热阻计算 在SRM中,有两类等效热阻,一类是传导热 r AB ——磁密波形的脉动峰值; n——磁密波形一个周期内的脉动次数。 对于涡流损耗的计算,引入如下涡流修正 因子: 耋 P :K P +KhPh (5) 阻,其计算公式为 (6) 最终得到修正后的SRM铁耗计算公式如下: 2.2铜耗计算 R… = i A (8) 式中: ——热传导距离; A——材料的导热系数; A——热传导通道的截面积。 P =m R (7) 铜耗P 的求取非常简单,其计算公式如式 (7)所示: 另一类是对流热阻,其计算式为 1 式巾:m——相数; R… =—L A (9) , ——相电流的有效值; R——电机相电阻。 式中: ——热对流系数。 3.1.2热对流系数计算 需要说明的是,南高频引起的集肤效应会使 相绕组电阻尺增大,进而使P 增大,为了减轻这 一在本文SRM初始设计方案中,热对流由j部 分组成:机壳与外界空气间的自然对流、定转子 不利影响,绕组导线可采用铜绞线。 根据上述公式,对SRM初始设计方案进行了 与气隙空气问的强迫对流,以及端盖内空气与电 机问的强迫对流。 2.3损耗计算结果及分析 (1)机壳与周围空气问对流系数。 当电机向周围空气自然传热时,设机壳壁表 面温度与室外温度始终相同,则表面热对流系数 叭n损耗计算,结果如表2所示。 表2 SRM初始设计方案损耗计算结果 为 厂 h =14(1+ ) √去 (10) 式中: ——气体吹拂效率系数,本文取0.5; 由表2可知,南于额定工作状况下的转速较 高,电机的损耗主要是铁耗,约占94.5%,而铁耗 又主要集中在定子上,定子铁耗约占总铁耗 的74%。 ——外界空气流速; 外界环境温度。 卜(2)定转子与气隙空气间对流系数。 该对流系数可由式(11)计算: 8— 电 乙-号}柱芾J应闭2013,40(9) : × N =0.386 P0 ’ /gapV?'R。:—moor —式中: ,——努赛尔数; ——空气热导率; z ——气隙长度; P ——空气普朗特数,本文取0.7; R ——雷诺数; ——。转子速度; r ——转子外径; ——空气的粘滞度。 (3)端盖内空气与电机间对流系数。该对流 系数可由式(12)计算: 』 c= 5・5(o・29Vsp+ ) (12) L =。Vm。-q 式中: 。——转子等效转速; 77——等效转速系数,本文取0.5。 3.1.3初始设计方案热网络模型 通过以上两类热阻的计算,建立了SRM初始 设计方案的等效热网络模型 ,如图3所示。 表3给出了模型中各热阻的数值。 图3 SRM初始方案等效热模型 研究与设计EMCA 表3 SRlVl初始方案等效热模型热阻值 ℃/W 电阻 阻值 电阻 阻值 电阻 阻值 电阻 阻值 R1 0.733 4 R1O 0.031 Rl9 0.055 3 R28 0.028 5 R2 0.0061 R1l 0.043 2 R20 0.934 7 R29 0.034 9 R 2.317 8 Rl2 5.502 2 R21 19.6191 R3O 0.305 2 R4 0.032 9 R13 0.013 3 R22 0.165 9 R31—0.007 8 R 0.0064 R14 3.027 2 R23 0.749 9 R32 0.0194 R6 0.007 3 R15 2.389 4 R24 0.276 7 R3]0.0301 R7 0.0566 Rl6 0.546 8 R25 10.738 7 R34 0.123 2 R8 0.041 5 Rl7 0.040 5 R26 2,3122 R35 2.】98 5 R。一0.012 R18 0.275 6 R27—0.0104 R36 2.409 2 3.1.4初始设计方案温升计算结果及分析 利用电路中的相关定理,可以很方便得到各 节点的温升数值,如表4所示。 表4 SRM初始方案各部位温升值 R R R ℃ 由表4可看出,由于电机转速较高,损耗较 大,同时空气导热性较差,导致电机内部温升极 高,实际上根本不能使用,因此有必要设计一个冷 却系统。 3.2冷却系统设计 众所周知,热传递有三种方式:传导、对流和 辐射,他们均可用于电机散热系统的设计。电机 所采用的散热方式需要视具体的功率等级和应用 环境而定 J。由于本文中电机功率等级与功率 密度均较大,因此采用强迫水冷散热。 本文设计了一个定、转子双水冷系统,图4给 出了带有冷却系统的SRM截面图,图中阴影部分 为所设计水冷系统的水槽。为了降低水冷系统对 电机电磁特性的影响,定子冷却系统的水槽位于 电机定子外围的铝壳内。转子冷却系统的水槽则 直接开在电机轴的中心,为了兼顾轴的刚性,水槽 直径选为8 mm。 图4带有冷却系统的SRM截面图 一9~ R R 研究与设计EMCA 3.2.1冷却系统热阻计算 以定子水冷系统为例,其等效热阻可通过以 下计算得到: = A :6PL (13) P=2 +(r一 )羔】 式中:hw——水冷系统的热对流系数; A ——铝壳与水的接触面积; P——水槽周长: ——定子轴向长度; r 、r舵、r们——分别为水槽的内、外半径及圆 端半径,且r 1=76 IllrTl,rw2=90 I/lrn,rw3=7 lqlnl。 3.2.2热对流系数h 的计算 h 的求解过程较为繁琐,可分为以下四个 步骤: (1)当量直径的计算。 由图4可见,定子水冷系统的水槽是非圆形 的,所以需要首先计算其当量直径D :Ⅳ l D = 4 一 J (14) R一 7 = +(r ) ] 式中:A ——水槽的截面积。 (2)雷诺数 的计算。 计算雷诺数的目的是判断水冷系统中水流的 形式。液体流动主要有两种形式:层流和湍流。P一王^5一一互 ,一 其中,湍流更有利于散热。在非圆形通道中;当 R <2 300时,液体流动完全为层流;当R。>10 000 时,液体流动完全为湍流;2 300≤R ≤10 000时,液 体流动处于从层流到湍流的过渡状态。 对于本文定子水冷系统而言,R 可由式(15) 得到: R : (15) w 式中:u——水流的平均速度,本文取10 rn/s; ——水的粘滞度。 本文所设计的定子水冷系统,雷诺数计算结果 为R =1 462 000,所以其水流形式完全为湍流。 (3)努赛尔数Ⅳ 的计算。 由于所设计水冷系统的水流形式完全为湍 流,所以其努赛尔数的计算公式可由Gnielinski方 程得到: 迫札再控制应阁2013,40(9) (16) 式中: ——摩擦系数; P ——水的普朗特数,本文取3。 对于光滑的通道而言 可由式(17)得到: = (17) (4)h 的计算。 h 的计算式如下: (18) 式中:kw——水的热导率。 由上述公式可以得到定子水冷系统的热阻。 转子水冷系统热阻的计算与定子相同,区别仅在 于无需计算当量直径。 3.2.3 带冷却系统的SRM热网络模型 图5给出了带冷却系统的SRM完整热网络模 型,表5给出了定、转子冷却系统中主要热阻的数值。 — 享 素 ~ R R R B 凡 L 堂 篓 图5带有冷却系统的SRM完整热网络模型 电札与控制应用2013,40(9) 研究与设计 E eA 表5定、转子冷却系统中主要热阻值 电阻 R41 R“℃/W 阻值 topologies and technologies for electric,hybrid,and 阻值 0.654 3 0.000 327 电阻 R42 R45 阻值 0.000 133 0.001 1 电阻 R43 R39 fuel cell vehicles[c]∥IEEE Vehicle Power and Propulsion Conference,2008:1-6. 0.000 321 0.005 3 [4]CARDENAS R,PENA R,PEREZ M,et a1.Control of a switched reluctance generator for variable—speed 3.2.4带冷却系统SRM温升计算结果及分析 表6给出了带有冷却系统的SRM各部分温 wind energy applications[J].IEEE Transactions on Enery Convergsion,2005,20(4):781—791. 升值。 表6带冷却系统SRM各部位温升值 ℃ 由表6可看出,采用定、转子双水冷系统之 后,SRM各部分的温升得到大幅减小,降至实际 运行允许的范围之内,充分证明了所设计冷却系 统的有效性。需要说明的是,由于采用了转子水 冷,在应用中必须采取一定的动密封措施。 4 结 语 本文针对一个30 kW、3相6/4极SRM的初 始电磁设计方案,通过对其磁密波形的傅里叶分 解计算了铁耗及铜耗,并分析了损耗的分布特点。 建立了该SRM详细的集总参数热网络模型,对其 温度场进行了计算,得到了各关键部位的温升值。 在此基础上,设计了一个定、转子双水冷系统,通 过对采用水冷系统前后电机各关键部位温升值的 比较,验证了所设计水冷系统的有效性。本文的 研究内容和相关结果对SRM温度场计算及冷却 系统设计具有一定的借鉴意义。 【参考文献】 VIJAYAKUMAR K。KARTHIKEYAN R.PARA— MASIVAM S. et a1. Switched reluctance motor modelling,design,simulation,and analysis:A comprehensive review『J].IEEE Transactions on Magnetics,2008,44(12):4605-4617. [2] SONG S J.Detailed design of a 30 kW switched reluctance starter/generator system used in more/all electirc aircraft[M].Germany:Shaker Verlag,2009. [3] ZHU Z Q,CHAN C C. Electircal machine [5]王述成.开关磁阻电机温升研究[D].武汉:华中 科技大学,2004. [6]wU W,DUNLOP J B,COLLOCOTF S J,et a1. Design optimization of a switched reluctance motor by electromagnetic and thermal finite・element analysis [J].IEEE Transactions On Magnetics,2003,39 (5):3334—3336. [7]BOGLIETI'I A,CAVAGNINO A,LAZZARI M,et a1.A simplified thermal model for variable—speed sefr-cooled industry induction motor[J].IEEE Trnasactions on Industry Application,2003,39(4): 945-952. [8] MILLER T J E.SPEED’S PC—SRD version 8.7 user’ S manual[D].University of Glasgow,2007. [9] 丁文,周会军,鱼振民.开关磁阻电机磁通波形与 铁耗计算[J].电机控制与应用,2006,33(6): 11.17. [10]杨红权,辜承林,李承.角度位置控制时开关磁阻 电动机的磁通波形及损耗计算[J].微特电机, 2009,37(5):9—11. [11] 吴建华.开关磁阻电机设计应用[M].北京:机械 工业出版社,2000. [12] MATERU P N,KRISHNAN R,Estimation of switched reluctance motor loss 『J]. IEEE Transactions on Industiral Application,1992,28(3): 668-679. [13] MELLOR P,ROBERTS D,TURNER D.Lumped parameter thermal model for electircal machines of TEFC desing[J].IEE Proceedings B Electric Power Applications,1991,138(5):205—218. [14]JAWAD F,IRANPOUR R,Thermal model for a switched reluctance motor of TEFC design during steady state and transient operation[J].Electric Machines and Power System,1998,26(1):77—91. 收稿时期:2012—11—27