基于Ansys Workbench雅阁ISG温度场仿真分析

李新华，杨国威，李哲然

（1.湖北工业大学电气与电子工程学院，湖北 武汉 430068；2.华中科技大学控制科学与工程系）

轻度混合动力汽车集成式起动-发电机ISG（ISG:Integrated Starter Generator）功率和转矩密度高、运行工况多变，特别是工作环境温度高、散热条件差，这些都给电机设计带来了新的挑战，仅按常规电磁设计是不够的，还需要对其进行温度场的仿真分析与设计。

本田汽车公司的雅阁 （Accord）混合动力汽车采用中度混合 （mild hybrid），即并联混合动力方案，发动机与电动机同轴，传动结构简单，与普通汽车发动机室差别不大。据本田厂方数据，雅阁混合动力汽车在城市路况下，百公里综合油耗仅为8.1L，这对于一台3.0L 6缸发动机来说已经相当不易了。雅阁ISG是一款在业界具有广泛影响的电机。该电机为16极／24槽配合永磁同步电机，采用组合式定子铁心结构，q=0.5分数槽集中绕组，见图1；转子为内置式V形磁钢 （接近一字形），每极有3个磁桥， 两极之间有V形沟， 见图2［1］。

目前电机温升计算方法有3种，即简化公式法［2］、等效热路法［3］和温度场法［4，5］。 简化公式方法比较简单，只能计算电机的平均温度，计算结果不太精确；等效热路法计算精度比简化公式法高，但要提高计算精度，需要增加网络节点和热阻数，计算工作量大大增加；温度场法采用现代数值方法来求解热传导方程，将求解区域离散成许多小单元后，在每个单元中建立方程，再对总体方程组进行求解。温度场法是一种快速和准确的数值计算方法，是现代电机温升计算的主流方法。

1 温度场仿真方法

电机内部存在损耗导致电机发热，使其温升增高。电机内部损耗的组成比较复杂，但主要是铜耗和铁耗。本文研究铁耗和铜耗所引起的电机发热。铜耗可按路的方法计算，而铁耗必须用有限元方法来计算。

电机铁耗由定子铁耗和转子铁耗两部分构成。对于同步运行ISG电机，由于转子铁耗较小，从简化ISG温度场仿真的角度出发，可以认为定子铁心内各处铁耗密度相等，即视定子铁心内为均匀磁场，且只考虑定子铁心内的铁耗。

然而由于电机定子铁心各处磁密并不相同，定子铁心各处的铁耗密度也会不一样；另一方面，谐波磁场也会在转子铁心中产生一定铁耗。如果均值铁耗密度代入，进行电机温度场仿真会有一定误差。为了考虑上述问题，可以采用Ansoft二维有限元与Ansys Workbench联合仿真方法分析电机中的温度场，基本步骤如下。

第1步，计算电机给定工况下的内热源，铁耗通过二维有限元动态仿真计算，即根据已知数据建立二维有限元铁耗模型并进行给定工况下ISG的铁耗仿真。

第2步，建立ISG三维温度场有限元模型，将该模型和第1步铁耗仿真结果导入Ansys Workbench软件。

第3步，给定边界条件，在Workbench环境下进行三维有限元温度场仿真。

根据文献［6，7］提供的数据，建立雅阁ISG的1／8二维有限元铁耗模型 （图3）并仿真；运用Solidworks软件建立ISG三维温度场模型 （图4），并将该模型和第1步给定工况下ISG铁耗仿真结果导入Ansys Workbench软件中，图5为最高转速工况时导入后的铁耗分布。图6a为该工况下铁耗不均匀分布时定子温度场的仿真结果。如果视铁耗为均匀分布，可将内热源计算结果直接施加于ISG三维温度场模型并仿真，仿真结果见图6b。

从图6可以看出，铁耗均匀分布时，定子最高温度出现在定子齿中心，达195℃，且每个齿的温度分布相同，但这只是一种近似算法。铁耗不均匀分布时，定子最高温度出现在定子齿边缘，达197.7℃，且每个齿的温度分布并不相同。显然，铁耗不均匀分布时温度磁场的仿真结果更接近于实际情况，但数据导入比较费时。以上2种方法温度场最高温度计算结果比较接近，定子最高温度误差只有1.4%。

2 温度场仿真结果及分析

按前面所述方法对雅阁ISG不同工况下的温度场进行仿真分析。仿真时，ISG的环境温度设定为70℃，ISG模型与空气对流换热系数取5×10-6W／mm2·℃。表1给出了雅阁ISG不同工况下铜耗的计算结果。

表1 雅阁ISG不同工况下铜耗的计算结果

将不同工况时的铁耗仿真结果直接导入Ansys Workbench软件，然后使用Ansys Workbench对雅阁ISG进行温度场仿真。考虑到最大功率为短时工作制，这里采用暂态温度场仿真，仿真时间取100 s，其它工况采用稳态温度场仿真，仿真时间取1h。额定工作点、最大功率工作点和最高转速工作点3种工况下ISG温度场仿真结果如图7所示。表2给出了ISG内各部件的最高温度及所在的位置。

仿真结果表明，对于额定工况，电枢绕组端部温度最高，定子铁心 （齿部）次之，转子磁钢处最低；对于最高转速工况，定子铁心 （齿部）温度最高，电枢绕组 （槽内导体）次之，转子磁钢处最低；最大功率工况由于工作时间短，电机内各部件的温度都较低。由于铁耗急剧增加，3种工况中最高转速工况电机内各部件的温度最高，特别是磁钢温度接近200℃，应引起高度关注。

表2 ISG内各部件的最高温度及所在的位置

文献 ［1］公布了雅阁ISG电枢绕组的温升试验结果，见图8。当负载转矩为70Nm时，绕组 （端部）最高温度由起始的160℃直线上升；持续一段时间后负载转矩下降至65 Nm左右，此时温度停止上升，维持在160℃左右；再经过一段时间后负载转矩升至75 Nm左右，绕组端部温度再次上升，最后稳定在200℃左右。

为了验证本文温度场仿真方法的正确性，对上述工况雅阁ISG进行了温度场仿真。绕组端部最高温度仿真结果见图9，仿真时没有考虑负载转矩的波动。比较两图可知，绕组端部最高温度曲线略低于试验曲线，但2条曲线基本吻合。由于仿真时只考虑了铁耗和铜耗，而没有计及其它损耗对温度的影响，使得绕组端部最高温度曲线略低于试验曲线。

3 温度对转子的影响

温度上升对ISG的性能以及结构安全性等都会带来一系列影响，下面重点讨论对转子稀土磁钢和磁桥结构安全性的影响。

3.1 温度对磁钢的影响

众所周知，稀土磁钢是一种耐温能力较差的永磁材料。温升的较大上升，不仅导致稀土磁钢磁性能的下降，还可能会造成磁钢的不可逆退磁。图10、图11分别为雅阁ISG额定和最高转速两个工况下转子磁钢的温度场仿真结果。

从图10、图11可以看出，磁钢高温区始终出现在轴向端部，最高温度点在共顶角，因此必须加强对磁钢端面的散热；另一方面，额定工况下磁钢最高温度为146.7℃，如选用耐温150℃的稀土磁钢尚可安全运行，而最高转速工况磁钢最高温度达194.2℃，耐温150°C的稀土磁钢则有失磁风险。

图12为最高转速工况下磁钢最高温度与工作时间关系的温度场仿真曲线。曲线表明，最高转速工况下ISG工作时间小于28 min，磁钢最高温度不超过150℃，选用耐温150℃的稀土磁钢磁钢是安全的，否则会存在失磁风险。

3.2 温度对转子磁桥结构的影响

表2仿真结果表明，ISG转子最高温度一般出现在磁桥处。出现这一现象的原因在于磁桥磁密高，铁耗大。由于转子高速运行，磁桥要承受很大应力，它的结构安全引人关注。为了考察温度对转子磁桥应力的影响，分别在不同温度条件下对雅阁ISG转子进行了应力仿真，图13、图14分别给出了22℃和197℃时的应力仿真结果。

由图13、图14可见，22℃时，雅阁ISG中间磁桥应力为47.2 MPa，热形变为0.00214 mm；197℃时，应力为1181.2MPa，形变为0.04576mm。温度上升了5.5倍，应力和形变分别上升了25.1倍和20倍。图15为雅阁ISG转子磁桥应力和形变随温度变化的关系曲线。可见，随着温度的增加，转子磁桥应力和形变会直线上升，威胁到转子结构的安全。因此，从转子结构安全的角度出发，必须限制最高转速工况下ISG的工作时间，使磁桥应力和形变在允许范围内。

4 结语

按本文所述方法对雅阁ISG进行了温度场仿真，仿真结果与文献［7］中的相关试验结果基本吻合；不同工况下ISG的温升情况存在较大差异，最高转速工况电机温升最高；从转子结构安全的角度出发，必须限制最高转速工况下ISG的工作时间。

［1］R.H.Staunton， T.A.Burress， L.D.Marlino.EVALUATION OF 2005 HONDA雅阁［R］.U.S.DEPARTMENT OF ENERGY， 2005.

［2］李文字.同步电机热计算与通风计算［J］.船电技术，2007， 27 （2）:72-75.

［3］J.Pyrhonen， T.Jokinen， V.Hrabovcová， Design of Rotating Electrical Machines［M］.John Wiley&Sons， Ltd， 2008.

［4］陈志刚.等效热网络法和有限元法在电机三维温度场计算中的应用与比较［J］.中小型电机，1995，22（1）:3-6， 35.

［5］李广德，何刚，何文秀.大型水轮发电机定子三维温度场计算［J］.大电机技术， 2000， 4 （2）:1-5.

［6］王水发，陈德为.基于ANSYS的异步电动机二维稳态温度场分析［J］.电气传动自动化， 2011， 33 （2）:23-26.

［7］李科，胡欲立.永磁直流电机温度场分析［J］.鱼雷技术，2010， 18 （3）:228-230， 235.