0 引 言

随着社会的快速发展，研究和开发高效节能的电机来满足电力需求是必然趋势，而电机的温升作为衡量电机运行性能的一个重要因素，直接关系到电机的效率、使用寿命，并成为限制其性能的重要瓶颈。对于TEFC电动机而言，主要是电动机端部外风扇所产生的冷却风量通过散热翅对电动机进行冷却，而风扇损耗在电机机械损耗中占有较大比例，因此找到冷却能力强、风扇自身及电机损耗小的高效外风扇结构具有重要意义。

目前，国内外专家学者已经对电机外风扇的设计及优化进行了大量的研究[1-9]。在风扇设计方面，文[1]在离心式外风扇叶片设计中成功地应用了固定翼飞机机翼结构形状、轴流式风扇叶片和汽轮机叶片形状，从而设计出高效率、低噪声的新型离心式外风扇。文[2]通过理论估算电机实际工作所需的风压和风量,给出了后倾离心式风扇的优化设计方案。文[3]对中小型电机常用的离心式、轴流式及斗式风扇的设计,结合作者多年实践经验作了较完整的介绍。在风扇优化方面[7-9]，文[7]以高压电动机中的径向式离心风扇与改进的后倾式离心风扇为对象，对比分析了流体动力学计算和传统计算方法的结果，证明经验公式在计算风量、风压、风效等性能参数时存在误差较大的问题，应用计算流体动力学进行风扇优化设计更具有适用性。文[8]以一台高压异步电动机外风扇为例，对风扇叶片在不同偏转角下的风扇性能进行了比较分析，得到扇叶偏转角度的最优值。文[9]对某汽轮发电机轴流式通风冷却风扇的原始模型进行了CFD流场数值模拟，通过改变叶片安放角和叶片扭转角，得到冷却风扇计算效率和叶片安放角及叶片扭转角之间的关系。综上所述，国内专家主要是根据电机需要的风量来设计或优化外风扇，针对外风扇自身进行流场分析及性能计算，忽略了风扇结构的改变对电机温度场及流变特性的影响。目前发表的文献中，对不同结构风扇下电机的温升及流变特性如何变化、如何影响风扇冷却能力的研究却很少，尤其是包括电机及外风扇在内的整体温度场及流体场数值研究较少，所以结合电机传热特性研究多结构外风扇的冷却性能具有一定意义及工程实际价值。

应该说，无论哪一种“生理性水肿”都与疾病无太大关联，但是我们依旧可以通过改善生活习惯来减缓甚至避免。大家应当保持乐观情绪，坚持适当锻炼以增强体质，提高适应能力；饮食应坚持低脂肪、低胆固醇，少糖、少盐，多吃瓜果蔬菜和豆制品等；避免久坐久站，经常活动下肢；最后，保证良好的睡眠，起居有律。

本文以55 kW异步电机Y2-250M的风扇原始结构为例，进行流体场数值计算及分析，并对不同结构风扇下电机的温升情况及流变特性进行了详细的研究，并根据电机的温升情况明确不同结构风扇的冷却效果，进而指出了最优的风扇结构方案，为TEFC电机冷却风扇元件选择及优化奠定一定基础。

PA法是由加拿大AP&C公司独有的金属粉末制备技术。原理是通过将金属丝材通过矫直机按照一定的速率加入等离子热源熔化，之后熔融金属液滴通过高压气体进行雾化。使用PA法制备粉体有球形度高、纯度高、含氧量低的特点。但是由于其原料为原料丝，所以提高了原料成本，并且局限了其制备金属粉体的种类。

1 电机及外风扇模型与边界条件

1.1 流体场数学模型

本文仅研究电机稳态运行时电机及外风扇区域的流体流动情况，即流动属于定常流动。由于冷却流体的雷诺数Re>2300，属于湍流运动，采用湍流模型进行求解；马赫数很小，冷却介质不可压缩。流体三维湍流流动遵循以下控制方程：

质量守恒方程为[10]：

 

(1)

6)电机转速为1 475 rpm，旋转方向为顺时针。

(2)察汗乌苏河地区1∶25 000地球化学测量圈定综合异常60个，显示异常元素种类多，元素组合明显，异常强度高，浓集中心明显，为后期找矿工作提供了较为详实的地球化学资料。

动量守恒方程为[10]：

 

(2)

式中：ρ为流体密度；μ为动力粘度；p为流体压力。

能量守恒方程为[10]：

风影将头凑过去，咬住了她的小耳朵，说话细得像蚊子叫，他告诉了她关于师父的一个小秘密。他说，在山背后的桃花坞村师父有个相好的，那次师父与那个相好的在桃花林中幽会，让他给撞见了。从此，风影偷偷地跑出山门，师父也就开一只眼闭一只眼了。原来和尚也有女人，难怪这世上没有不偷腥的猫。原来，师父有如此隐秘而严重的把柄落到了风影的手中，难怪他胆大包天敢跑出寺院来玩。红琴什么也不说，一个劲地笑，风影让她发誓，不许告诉任何人，红琴就向他保证，这事一定彻头彻尾地烂在肚子里。

协议书模板分为有界桩的界线协议书、无界桩的界线协议书2个模板，分别包括封面、协议主体内容、界线所涉及乡镇政府代表签字3个部分。协议书主体内容包含工作概况、重要问题处理结果、行政区域界线走向和界桩位置说明、行政区域界线的维护和管理、最后条款（补充说明）以及附件6个部分，其中附件包括了附图、界桩登记表、三交点界桩登记表、界桩成果表、界址点成果表、三交点成果表等内容。三交点协议书模板包含标题，协议内容，三交点的位置、坐标，底图所在的图幅以及所涉及的乡镇政府代表签字，所属县民政局代表签字和协议日期。

(3)

其中：cp为流体比热容；T为温度；λ1为流体导热系数。

1.2 电机及外风扇物理模型

如图1为电机及外风扇整体计算域模型。

  

图1 电机及外风扇计算域模型

本文电机中的外风扇为径向离心式风扇，主要由叶片、后盘组成。其中后盘为平直后盘，风扇叶片数量为10个，叶片长度115 mm，高度110 mm。其物理模型如图2所示

  

图2 外风扇物理模型

1.3 基本假设与边界条件

2)认为涡流效应对每根股线的影响相同，定、转子绕组铜耗取其平均值，并认为端部是直的；

1)认为槽内所有绝缘(包括槽楔)的热性能与主绝缘相同；

为合理的简化求解，本文给出的基本假设与边界条件如下：

3)空气自由流入入口，入口采用pressure-inlet边界条件；

4)用空气域包裹出口，空气自由流出出口，出口采用pressure-outlet边界条件；

5)紧贴风扇及紧贴转子的流体域设为旋转流体域，设置成与电机相同的角速度154.46 rad/s，其他区域为静止区域，扇叶及后盘所在的壁面设置为moving wall，速度为0；

式中：u、v和w分别为x、y、z方向的速度分量。

2 电机原始外风扇结构温升及流变特性计算

2.1 电机温度场及流变特性计算

研究对象选取图1所示的计算域模型，为避免回流现象的产生，离心式风扇端部空气域进行了加长，结合电机各部分材料属性，将电机各部分热源损耗及边界条件施加到计算域的相应部分，经过有限体积元计算，得到电机在额定负载时的稳态温度场。其中电机机壳的整体温度分布如图3所示。

为探究更改叶片形状对流经机壳气流的风速影响，将机壳外侧设置三条采样直线，长度从风扇出口到另一侧端盖的外表面，如图8所示,其中直线1位于接线盒上方，直线2位于机壳中部上方，直线3位于机壳下部。经过数值计算，给出了三条采样位置上的风速分布特征，如图9所示。

Discussion on Requirements and Forms of Fire Water Source in Nuclear Power Plant WANG Jinlong,LIU Huanming,ZHANG Longzhou(113)

温度最低的地方出现在机壳端部靠近风扇侧散热翅的地方，是因为此处是外风扇吹动的空气最先经过的地方，空气温度最低，并且流动流畅，故该区域温度最低。

由于风扇端的空气经风扇作用强制吹过机壳，到达机壳传动端时温度有所升高，并且空气由风扇端到达另一端时，风量有所减弱，造成空气冷却能力下降，所以电机机壳传动端的温度比风扇端的温度高。

  

图3 电机机壳温度分布

  

图4 电机机壳外迹线分布

经过电机机壳空气的迹线图如图4所示，迹线图起点为外风扇出口，从图中可以看出，空气经风扇作用后沿着风扇后盘与风扇罩之间的空间吹入电机机壳，一部分空气受到端盖及散热翅的阻碍作用，导致空气经过端盖后速度下降，以较低的速度吹过散热翅风沟到达另一端，而另一部分则在风扇处形成扰流。

气流速度最低的地方出现在机壳中部，这是由于气流在机壳中部相邻散热翅风沟间做横向运动，相邻散热翅对气流的阻碍作用较大，导致此处的气流速度急剧下降。

电机顶部相对于电机侧面机壳来说，风量较小，气流速度较高，是因为此处的气流不在相邻散热翅间做横向运动，而是沿着散热翅风沟做轴向运动，散热翅对其阻碍作用最小所致。由于空气在机壳外表面沿轴向流动时，受到吊装部件、接线盒及基座筋的阻碍，这些位置附近空气流动出现绕流现象，对流换热能力差，因此机壳部分的最高温升出现在近接线盒位置。

2.2 计算结果与实验数据对比

本文中的实验利用PT100对定子绕组及转子结构件进行温升测定，图5为传感器埋设位置示意图。

2.2 两组血清标志物比较 利用3 326例孕妇AFP、β-hCG和uE3的血清学结果对相应孕周血清学标志物中位数值计算出其MOM值，运用独立样本t检验对MOM值进行统计学分析。结果表明，高风险组孕妇血清中β-hCG值大于低风险组孕妇，差异有统计学意义(P<0.05);而低风险组AFP和uE3指标的MOM值高于高风险组孕妇，差异有统计学意义(P<0.05)。见表2。

  

图5 传感器位置分布图

为验证数值计算结果的正确性，将得到的实验温升值与数值计算结果进行对比，如表1所示：

沿周向给定子绕组依次编号，靠近接线盒一侧的定子绕组编号为1，从定子绕组远离电机端部离心式风扇的一侧的远风端开始，取轴向位置“A-F”点定子绕组截面进行分析。其中传感器埋设在1号槽底部，靠近远风端铁心位置30 mm处，即“B”位置，则用“1-B”表示；传感器位于2号槽槽底位置，轴向C点，即“2U-C”，其中U代表靠近槽底绕组，D为靠近槽口绕组，其余位置可依次类推[10]。

 

表1 实验温升值与数值计算结果对比

  

传感器编号实验值/K计算值/K相对误差/%1⁃B8689063562⁃C8385312787⁃A9197106707⁃F92878445215⁃E83841814217⁃D858655182

通过比较可知：测温点当中最大及最小相对误差分别为6.70%和1.42%，由此可见，计算所得数据与实测得到数据基本吻合，证明了计算的准确性。由此可见，此计算模型和计算方法是正确的，能够满足工程需要。

进入二十一世纪，乡愁是一段绘声绘色的视频。语音电话、视频聊天，移动互联让手机功能越来越强大。去年春节，我为母亲购买了一部智能手机，并手把手教会她微信视频聊天。“妹子，你好！”接通视频通话后，大舅和蔼可亲的笑脸出现母亲的手机大屏上。“好，好，我们都很好！”那天，母亲与大舅“面对面”谈了好久，开心的笑容一直挂在母亲的脸上。

3 基于流热特性研究的电机外风扇结构特性优化分析

从离心式风扇的设计要点可以看出[11]，根据电机实际结构确定外风扇外径D2、电机所需工作流量qVm及风扇额定工作时的压力p，则风扇内外径、扇叶宽度是唯一确定的，根据设计要点核算此原始结构风扇，各参数均在合理范围内，而通过对机壳温度分布及迹线分布的分析，可以看出在流经机壳中部及顶部的气流中高速气流所占比例不高，特别是机壳中部高速气流较少，而这两处是温升最高的部位，若能提高这两处的风速，则会进一步降低电机温升，提高风扇的冷却能力。那么尝试通过改变风扇的叶片形状或叶片入口角来提高流经机壳表面的风速，进一步降低电机温升及风磨损耗，提高外风扇的冷却能力。

3.1 叶片形状的优化

原始风扇为梯形叶片，叶片切除角度为30°，如图6所示。

  

图6 原始扇叶草图

为了探究合适的风扇叶片形状，本文提出了I、II、III三种叶片，扇叶切除角度分别为10°、20°和35°，并在风罩中部做截面，分别进行数值计算，得到风罩截面的空气迹线图如图7所示。

对比图7不同形状叶片下风罩中的扰流情况，可以看出由于风扇的作用，空气在风罩内呈螺旋状运动，并且在风罩内壁处达到最高速度。各方案扇叶之间均存在扰流现象，但由于扇叶表面积的不同，导致风扇作用下的空气流动情况有所差异，因此扰流情况也有所不同。其中方案I风罩内扰流现象最为严重，方案III扰流现象最为轻微，原始风扇、方案II扰流情况类似，介于二者之间。可见对于梯形叶片，叶片形状同风罩内扰流情况有一定联系，增大扇叶切除角度可以减轻风罩内扰流现象。

2011年建昌县承担国家级玉米高产创建示范区3个，省级玉米高产创建示范区1个，重点推广种植辽单565、丹玉96-22等耐密型玉米品种。全县四个示范区平均产量886 kg/667 m2，对照田平均750.7 kg/667 m2，增产15.3%；全县玉米平均593.5 kg/667 m2，增产49.3%。

  

图7 不同叶片形状下扰流情况

  

图8 机壳外采样直线

由图3及数值计算结果可知，机壳整体轴向温度分布中间较高，逐渐向两侧递减，原因是定子铜耗、定子铁耗及转子通过气隙传递过来的能量传递到定子铁心，而定子铁心与机壳是紧密接触的，所以机壳中部温度较高。而由于电机接线盒对风扇强制吹过空气的阻碍作用，空气冷却能力迅速下降，导致此处的温度最高。

目前，我国相继出台新能源的管理政策，采用创新的服务模式，明确系统的评审职责，对新能源管理给予了政策支持和保护。但是，在实际管理过程中，并没有明确各个环节的职责，导致工程项目漏洞百出，一旦出现问题找不到出错的环节，也无法找到责任人，进而出现推诿的现象，严重影响后续工作的顺利进行。

  

图9 不同叶片形状下风速对比图

图9为在原始结构及三种方案下的三条采样直线风速对比图，气流自风扇流出后，随着轴向长度增加，散热翅顶部附近空气流动速度有升高的趋势，但由于气流受到接线盒和基座筋的阻碍作用及在散热翅风沟内的能量损失，风速均呈现为先增加后减小的趋势。机壳中部的风速总体上要高于机壳顶部及底部，这是因为同机壳顶部及底部不同，气流在机壳中部没有阻碍，气流可在机壳中部沿轴向散热翅风沟径直到达另一端。

由于叶片形状不同，风扇产生的压力有所差异，因此在机壳表面同一位置的风速也因此不同。从图9(a)可以看出，三种方案经过接线盒前的风速均比原始方案大，方案I和方案II的风速分布相似，方案III风速稍低些，但方案I的风速要高于其他方案。而从图9(b)可以看出，方案III的风速最高，原始方案稍低，方案II次之，方案I风速最低。从图9(c)可以看出，三种方案风速均比原始方案低，风速情况很类似，但方案III的风速稍高些。

由图9可以得出结论，适当减小扇叶切除角可以提高流经机壳顶端的风速。

  

图10 不同叶片形状下温升对比图

由于在机壳周向上，较高温升出现在机壳顶端及中部，则通过比较采样直线1和2轴向温升情况来确定优化方案。图10为在原始结构及三种方案下的采样直线温度对比曲线图，可以看出机壳温升与机壳表面风速有明显的对应关系，这是因为电机运行产生的热量传递到机壳，而机壳表面风速越高，电机表面与周围空气的热交换速度越快，电机冷却效果越好，电机表面温升越小。从图10(a)可以看出，采用方案III叶片的机壳顶部温升最高，方案I的机壳顶部温升最低，其他形式叶片下的机壳顶端温升基本相同；从图10(b)可以看出，三种方案下机壳中部的温升基本一致，均略低于原始方案。可见扇叶形状的改变可较小程度降低机壳顶部及中部的温升，减少机壳轴向温升梯度，其中采用方案I叶片的机壳轴向温升梯度最小，机壳承受热应力较小，机壳表面温度分布最均匀。

一是“政府工程模式”。贫困地区政府依托“阳光工程”“送教下乡”工程，将新型职业农民培育纳入其自我管理体系中，紧密结合当地农业生产特点和农民实际需求，通过派遣农业专家深入基层授课、推广远程培训、实施“田间课堂”等形式，形成政府主导的新型职业农民培育模式。

不同叶片形状下电机各部分最高温度如表2所示，从表2可以看出，改变叶片形状对电机总体温升影响较小，定子部分和转子部分温升变化趋势一致，叶片形状对定子部分的温升影响比转子部分稍大，这是由于定子部分产生的热量被风扇作用的空气直接带走，而转子部分产生的热量通过热传递传递到定子部分再通过机壳散发出去，此部分温升受外风扇的影响小。

 

表2 不同叶片形状下电机最高温度对比 K

  

电机部位方案I方案II方案III原始方案机壳35650357113585235657端盖34267342933435334237定子压板37413374563755637413定子绕组39613396463972339599定子绝缘39558395913967039544槽楔38642386973877938647定子齿37273373203743137274定子轭36187362443637836191转子压板43760437754381643778转子齿45122451384518045144转子绕组45123451404518145146转子轭45117451334517445139

由于电机外风扇运行过程中产生的风磨损耗转化为散发出风罩的内能，因此在耦合能量方程的情况下，利用有限元软件后处理求出的风罩出入口的内能差即为外风扇的风磨损耗。据此求出各方案下的外风扇的风磨损耗如表3所示。

 

表3 不同叶片形状下外风扇风磨损耗对比 W

  

方案I方案II方案III原始方案18672185431818918565

从表3可以看出，改变风扇叶片形状可以降低风扇的风磨损耗，其中方案III风磨损耗最低，同外风扇原始结构相比下降了2%，这是由于电机风磨损耗由转子旋转引起的转子表面与冷却气体之间的摩擦损耗以及安装在电机转轴上或由电机本身转轴驱动的风扇所需功率组成[12]，而风扇的功率与风扇工作时的压力和流量有关，叶片形状的差异导致风扇产生的压力有所不同，而改变叶片形状可以减少风扇消耗功率，从而减少风磨损耗。

将部署区域T离散化为m×n个目标点集T={T1,T2,T3,…,Tm×n},其中目标点Tl的位置坐标表示为(xTl,yTl)(l=1,2,3,…,m×n)。则目标点Tl与节点Si的欧氏距离为:

3.2 叶片入口角的优化

空气进入风扇时，叶片相对速度的方向与叶片入口角相吻合时，此时对叶片冲击最小，风扇损耗最小。为探究改变叶片入口角对风扇冷却性能的影响，将叶片入口角在20°～30°间变化,分别对四种叶片在不同入口角下进行数值计算，如图11为III型叶片在不同入口角下的风罩截面迹线分布，I型、II型叶片的迹线分布情况与此类似。

  

图11 III型叶片在不同入口角下的扰流情况

从图中可以看出，同种叶片在不同的入口角下风罩内扇叶间的扰流情况有较大差异，其中入口角为20度时扰流现象最为严重，入口角为30度时较为轻微，这是因为叶片入口角不同，空气在扇叶作用下形成的流体场分布也有所差异。可见适当的入口角可以减缓风罩内扰流现象。

 

  

图12 不同入口角下风速对比图

由于空气流入风扇时的角度不同，导致在风罩内相应的流体场分布有所差异，气流在流出风扇时周向位置上的速度分布也因此不同。如图12所示为不同入口角下三处采样直线的风速对比图，在从图12(a)可以看出，随着入口角的增加，机壳顶端的风速也随之缓慢增加，增量比较小,但改变叶片入口角后的风速均比原始结构小；图12(b)可以看出，改变叶片入口角并不能提高在机壳中部的风速，风速在此处的最大值随着入口角的增大而减小；而从图12(c)可以看出，机壳底部的风速情况与叶片入口角并无明显的线性关系，改变叶片入口角不能提高此处的风速。由此可见，改变扇叶入口角并不能有效提高流经机壳表面的风速。

  

图13 不同入口角下温升对比图

从图13可以看出，改变叶片入口角不能降低机壳的轴向温升，这是因为改变叶片入口角不能有效提高流经机壳表面的风速，与之相对应的机壳表面也不能与周围空气有更快的热交换速率所致。从图13(a)可以看出，叶片入口角为25度时的机壳轴向温升梯度最高，20度次之，而其他入口角时的温度梯度轴向分布基本一致，比前两种情况温升梯度小；从图13(b)可以看出，原始结构风扇作用下机壳轴向温升梯度最大，其余入口角下的温升梯度均较低，其中叶片入口角为30度时梯度最小，入口角为20度的次之。由此可见，改变扇叶入口角并不能有效降低机壳表面的温升，但可以较小程度的的减小机壳表面的温升梯度，温度分布更均匀。

 

表4 不同叶片入口角下电机各部件最高温度对比 K

  

入口角I型叶片II型叶片III型叶片20°25°30°32°20°25°30°32°20°25°30°32°原始方案0°机壳35515355233559135621355393555535604356653569335738357353574935657端盖34204342143423934266341973421534237342753427134305343023430934237定子压板37315373243737837407373243734037381374333743437473374733748837413定子绕组39540395483959539624395393955339594396403961839653396563967339599定子绝缘39483394913953839567394833949739537395843956439599396013961839544槽楔38605386113862338645385693858038606386693869038692386983870638647定子齿37166371753723137260371783719437237372913730137340373423735437274定子轭36059360693612936159360793609636141362003622736268362703628136191转子压板43722437334377043803437034372143757438004514443769437804380143778转子齿45089450994513745170450714508745123451664510945134451454516845144转子绕组45090451014513845171450724508345124451674511145135451464516945146转子轭45084450944513245165450654508245118451614510445128451404516245139

从表4可以看出，更改扇叶入口角对电机各部件温升影响不大，在相同叶片入口角下改变叶片形状时电机各部分温升变化趋势不一致。

 

表5 不同入口角下外风扇风磨损耗对比

  

叶片入口角/°风磨损耗/WI型II型III型201880925187973018460321830920189422518876301858732182902018721251861230184953218295

通过表5可以看出，三种叶片下的风磨损耗均随着入口角的增加而降低，对于同种叶片而言，改变风扇叶片入口角也可以降低风扇的风磨损耗，这是由于空气进入风扇时，叶片入口角越接近叶片与空气相对速度的方向，空气对叶片冲击越小，风扇对应的风磨损耗就越小。

4 结 论

本文通过对TEFC电机不同结构特征风扇条件下的电机内流体场及温度场性能的分析研究，综合考虑了风扇的出风能力及产生损耗，明确了电机风扇设置的最优方案，具体得到如下结论：

1)减小扇叶面积或者采用合适的叶片入口角可以减轻风罩内的扰流现象，提高风扇的利用效率，更充分的发挥风扇的冷却能力；

2)改变扇叶切除角对电机总体温升影响较小，但会大幅减小风磨损耗，适当减少扇叶切除角可提高机壳顶部的风速；

3)单独改变扇叶入口角并不能有效提高流经机壳表面的风速，对降低电机各部件温升及风磨损耗作用很小；

4)调整扇叶切除角或叶片入口角均可以降低电机轴向温升梯度，叶片入口角为30度时轴向温升梯度最小，温度分布最均匀。

医疗体制逐步改革，法律法规逐渐完善，尊重和保护病人权益的呼声越来越大，导致医患矛盾加剧，病人往往不配合床边教学。病人周转快、住院天数缩短等也使得可供选择的床边教学病人减少。

参 考 文 献:

[1] 李贤明，刘文超，张永义．大功率高速电动机高效低噪声离心风扇的研究与设计[J]．防爆电机，2005，40(1)：9-12．

[2] 汪书苹，赵争鸣，冯垚径．YKK355-630 系列高压三相异步电动机高效风扇设计[J]．清华大学学报：自然科学版，2009，49(1)：9-12．

[3] 魏书慈，李士菊．浅谈中小型电机风扇的设计[J] ．电机技术，1996(2)：29-35．

[4] 罗次申．动力机械测试技术[M]．上海：上海交通大学出版社，2001．

[5] 王成．电机冷却风扇出口流量分配研究[D]．上海：上海交通大学，2008．

[6] 刘立峰．大型电机冷却风扇的实验测量与模拟计算[D]．上海：上海交通大学，2007．

[7] 李光宇，袁立强，赵争鸣，等．基于计算流体动力学的高效交流电机风扇优化[J]．电工电能新技术，2014，33(11)：24-25．

[8] 温嘉斌，唐跃，何金泽，等．高压异步电机外风扇流场数值计算及优化[J]．电机与控制学报，2013，17 (9)：79-86．

[9] 刘全忠，宫汝志，王洪杰，等．汽轮发电机冷却风扇的数值模拟及优化[J]．哈尔滨工业大学学报，2010，42 (3)：442-445．

[10] 王海涛．不同驱动条件下感应电机通风与传热分析[D]．哈尔滨：哈尔滨理工大学，2014．

[11] 陈世坤．电机设计[M]．北京：机械工业出版社，1997．

[12] 丁树业，孙兆琼，徐殿国，等．3MW双馈风力发电机传热特性数值研究[J]．中国电机工程学报,2012，32(3)：137-143．

[13] 魏永田，孟大伟，温嘉斌．电机内热交换[M]．北京：机械工业出版社，1998．

[14] 丁舜年．大型电机的发热与冷却[M]．北京：科学出版社，1992．

[15] 王福军．计算流体动力学分析-CFD软件管理与应用[M]．北京：清华大学出版社，2004．

[16] POTTER Wiggert．Mechanics of Fluids[M]．北京：机械工业出版社，2003．

[17] STATON D A，CAVAGNINO A．Convection Heat Transfer and Flow Calculations Suitable for Electric Machines Thermal Models[J]．IEEE Transactions on Industrial Electronics，2008，55(10)：3509-3516．

[18] TRAXLER-SAMEK G，ZICKERMANN R，SCHWERY A．Cooling Airflow， Losses，and Temperatures in Large Air-cooled Synchronous Machines[J]．IEEE Transactions on Industrial Electronics，2010，57(1)：172-180．

[19] LIU Y J，LEE Y，JUNG H K，et al．3D Thermal Stress Analysis of the Rotor of an Induction Motor[J]．IEEE Transactions on Magnetic,2000，36(1/4)：1394-1397．

[20] XYPTERAS J，HATZIATHANASSIOU V．Thermal Analysis of an Electrical Machine Taking into account the Iron Losses and the Deep-bar Effect[J]．IEEE Transaction on Energy Conversion，1999，14(4)：996-1003．

[21] SUHAS V Patankar．Numerical Heat Transfer and Fluid Flow[M]．Minnesota University Minneapolis Department of Mechanical Engineering，2009：130-156．

[22] HUANG Y，SRINATH V E．Detailed Heat Transfer Distributions in Two-pass Square Channels with Rib Tabulators and Bleed Holes[J]．International Journal of Heat and Mass Transfer，1998，41(23)：3781-3791．