简介:本资源是一篇聚焦电动汽车轮毂电机热管理关键技术的学术论文PDF,面向新能源汽车研发工程师、电机设计与仿真技术人员及高校相关专业研究生。针对轮毂电机因空间受限、功率密度高导致温升快、散热难、易失效的工程痛点,论文提出一套基于Maxwell与ANSYS协同的磁热耦合仿真方法,并创新性引入部件等效简化策略,在保证物理准确性的同时显著降低计算复杂度;全文含电磁损耗分析、温度场分布模拟、边界条件建模及试验验证全流程,具备强实践指导价值。资源为单个12.56MB PDF文件,内容完整覆盖摘要、引言、建模方法、仿真结果、实验对比与结论,含中英文摘要、参考文献及作者单位信息,结构规范,可直接用于技术复现与方案参考。目前已有187人学习下载,是开展轮毂电机热设计、电磁-热多物理场联合仿真不可或缺的专业参考资料。
1. 轮毂电机磁热耦合仿真不是“加个温度场就完事”:它决定电动车续航衰减快慢、爬坡力矩突降是否真实、甚至影响制动能量回收效率
一辆搭载轮毂电机的纯电轻卡,在-10℃环境连续爬坡15分钟后,驱动轮扭矩骤降18%,仪表未报任何故障码;另一台同型号乘用车在夏季高速工况下,单侧轮毂电机温升比对侧高23℃,导致ESC系统误判为车轮打滑。这些现象背后,都不是简单的“电机过热停机”,而是磁场分布畸变→铁芯涡流损耗激增→局部热点形成→磁路饱和加剧→反电势失真→控制环路震荡的级联过程。轮毂电机因直接集成于车轮,散热路径长、风冷效率低、空间约束严,其电磁场与温度场存在强非线性双向耦合:温度升高不仅降低永磁体剩磁,还会改变铜绕组电阻率和硅钢片磁导率,进而反向重塑磁通路径与损耗分布。本文聚焦“电动汽车用轮毂电机磁热耦合仿真分析”这一工程刚需,拆解从Maxwell瞬态电磁求解到Fluent稳态热场映射的完整闭环链路,给出ANSYS Workbench中实现双向耦合的参数配置关键点、网格敏感度验证方法、以及实测温升误差控制在±4.2℃内的典型收敛判据——所有步骤均基于2023版Workbench R2及后续版本验证,不依赖第三方插件。
2. 为什么必须用双向耦合而非单向热分析:磁性能参数随温度变化的非线性曲线是仿真可信度的生死线
2.1 永磁体与硅钢片的温度-磁特性必须分段建模,而非简单线性外推
轮毂电机常用钕铁硼永磁体(如N42SH)在80℃以上时,剩磁Br每升高1℃下降约0.11%,但该系数在120℃附近会加速至0.15%/℃;而无取向硅钢片(如35W440)的相对磁导率μr在20℃时约为1500,当温度升至150℃时,μr可能跌至900以下,且B-H曲线膝点明显右移。若在仿真中仅将材料设为“常温属性+固定温升”,会导致铁芯磁密计算偏差达35%,进而使涡流损耗预测值偏离实测值2.3倍以上。常见错误做法是直接导入供应商提供的20℃B-H曲线,忽略温度映射表。
提示:必须获取材料供应商提供的多温度点B-H曲线数据(至少含20℃、80℃、120℃、150℃四组),在Maxwell材料库中创建“Temperature-Dependent B-H Curve”类型材料。对于永磁体,需同时定义Br(T)与Hcj(T)函数——后者影响退磁风险判断,不可省略。
2.1.1 在Maxwell中构建温度依赖型永磁体材料的具体操作
在Maxwell 2023R2中,右键Materials → Add Material → 选择“Permanent Magnet”类型 → 在“BH Curve”选项卡点击“Import” → 选择包含温度列的CSV文件(格式:Temp(℃), Br(T), Hcj(T))。关键参数设置如下:
Br Temperature Coefficient:设为“User Defined Table”,禁用默认线性系数Demagnetization Curve:必须勾选“Use Temperature-Dependent Demag Curve”Operating Temperature:在Simulation Setup中设为变量,而非固定值
# 示例:用于生成温度依赖B-H曲线CSV的Python脚本片段(供材料测试数据拟合) import numpy as np temp_points = [20, 80, 120, 150] br_data = [1.32, 1.18, 1.05, 0.97] # 单位:T hcj_data = [-950, -820, -680, -590] # 单位:kA/m np.savetxt('n42sh_temp_brcurve.csv', np.column_stack([temp_points, br_data, hcj_data]), delimiter=',', header='Temp,Br,Hcj', comments='')此脚本输出CSV可直接导入Maxwell。注意:Br与Hcj必须同为负温度系数,若实测Hcj在高温下衰减更剧烈,则曲线斜率需更陡——这直接影响电机在峰值负载下的抗退磁能力评估。
2.2 绕组铜线电阻率必须按IACS标准动态更新,而非采用20℃标称值
轮毂电机绕组在持续大电流下温升可达130℃,此时铜电阻率ρ比20℃时高约42%。若仿真中仍用ρ=1.724e-8 Ω·m(20℃值),将导致焦耳热计算偏低,进而使温度场迭代发散。ANSYS Maxwell支持通过Conductivity vs Temperature表格定义铜的电阻率温度函数,但需注意:该函数必须与实际铜纯度匹配(如IACS 100%对应ρ₂₀=1.724e-8,IACS 98%则为1.759e-8)。
2.2.1 铜电阻率温度函数的工程化校准方法
实测某轮毂电机绕组在不同温升下的直流电阻,拟合得到ρ(T)=ρ₂₀[1+α(T−20)],其中α为温度系数。但标准α=0.00393/℃仅适用于20℃附近,高温区需修正。推荐采用双段线性模型:
- 20–100℃区间:α₁=0.00393/℃
- 100–150℃区间:α₂=0.00412/℃(实测铜线老化后α增大)
在Maxwell材料属性中,选择Copper → Edit → “Electrical Conductivity” → “vs Temperature” → 输入两段折线坐标点:(20, 5.8e7), (100, 4.1e7), (150, 3.5e7) S/m。此处电导率σ=1/ρ,单位S/m,数值需换算准确。
| 温度(℃) | 电阻率ρ(Ω·m) | 电导率σ(S/m) | 对应工况 |
|---|---|---|---|
| 20 | 1.724e-8 | 5.80e7 | 冷态启动 |
| 100 | 2.44e-8 | 4.10e7 | 连续满载30分钟 |
| 150 | 2.86e-8 | 3.50e7 | 极限爬坡或制动回收 |
该表格直接决定焦耳热源项Q=J²/σ的精度,是磁热耦合收敛稳定性的基石。
3. 在ANSYS Workbench中实现磁-热双向耦合的最小可行配置:从Maxwell瞬态求解到Fluent热场映射的6步闭环
3.1 创建双向耦合项目结构:System Coupling模块是强制必选项
单靠Maxwell导出热源再导入Fluent属于单向耦合,无法反映温度升高导致磁路变化的反馈效应。Workbench R2起必须启用System Coupling模块构建闭环。正确流程为:Maxwell(电磁瞬态)→ System Coupling(数据交换引擎)→ Fluent(热传导求解)→ System Coupling → Maxwell(更新材料属性)。项目树中需包含三个独立系统:Maxwell 3D、Fluid Flow (Fluent)、System Coupling,三者通过Setup链接。
注意:Maxwell必须运行在Transient模式(非Magnetostatic),因轮毂电机在加速/制动时电流波形含高频谐波,稳态假设会漏掉趋肤效应引起的附加损耗。时间步长需满足奈奎斯特采样定理——对基波频率f₀=200Hz的PWM驱动信号,建议dt≤1/(10×f₀)=0.5ms。
3.1.1 Maxwell瞬态求解设置的关键参数表
| 参数项 | 推荐值 | 说明 |
|---|---|---|
Time Step | 0.2 ms | 覆盖PWM开关频率(如10kHz对应0.1ms),此处取2倍余量 |
Total Time | 2.0 s | 覆盖一个完整工况循环(如NEDC工况中单次加速+匀速+减速) |
Adaptive Passes | 8 | 确保磁场收敛,尤其关注气隙区域网格 |
Save Fields | Every 5th step | 平衡存储空间与数据完整性,热耦合只需关键时间点 |
Solver Type | Time Domain Solver | 必须启用,否则无法输出瞬态损耗密度 |
3.2 定义热源映射面:轮毂电机中定子铁芯背轭与端部绕组是热源加载主区域
Maxwell输出的损耗密度(W/m³)需映射到Fluent网格上。但轮毂电机结构特殊:定子嵌入轮辋内腔,转子为车轮本体,二者间气隙极小(通常<1mm)。若将整个电机域作为热源加载,会因网格不匹配导致热量“泄漏”到空气域。正确做法是仅提取固体部件表面的损耗通量——即在Maxwell中对定子铁芯、绕组、永磁体分别设置Loss Density输出,再在Workbench中通过Geometry→Named Selections创建对应面组。
3.2.1 在Maxwell中精确提取绕组焦耳热与铁芯涡流热的分离策略
轮毂电机绕组常采用分布式绕组,槽内导线与端部引线损耗占比不同。在Maxwell中需:
- 对槽内部分:设置
Solid对象属性为Conductor,指定Current Density和Conductivity - 对端部引线:单独建模为
Stranded Conductor,启用AC Resistance计算 - 对铁芯:划分
Laminated Core区域,设置叠片厚度(如0.35mm)和叠压系数(0.95)
运行后,在Results→Fields中分别导出:
OhmicLoss(绕组焦耳热)CoreLoss(铁芯涡流+磁滞损耗)PermanentMagnetLoss(永磁体涡流损耗,常被忽略但高温下不可忽视)
导出格式必须为ANSYS CFD兼容的.cfd文件,而非通用.csv——后者无法保留几何拓扑关系。
3.3 Fluent热场求解器设置:轮毂电机特有的强制对流边界条件建模
轮毂电机无轴向风扇,主要靠车辆行驶时的迎面气流冷却。Fluent中需设置Moving Wall边界模拟车轮旋转带来的离心通风效应。典型设置:
- 轮辋外表面:
Wall类型,Heat Flux设为0(绝热),Shear Stress启用Moving Wall,速度按v=ω×r计算(ω为车轮角速度,r为半径) - 制动盘通风孔:
Pressure Outlet,静压设为-50Pa(模拟气流抽吸效应) - 轮毂罩内腔:
Interior类型,启用Porosity模型模拟多孔介质(等效散热片)
# Fluent命令行快速设置旋转壁面(在Text User Interface中执行) define/boundary-conditions/velocity-inlet/rotating-wall name wheel_outer_surface rotation-axis 0 0 1 # Z轴旋转 angular-velocity 31.4 # 300rpm对应31.4 rad/s该命令将轮辋外表面设为旋转壁,生成的切向速度场能显著提升近壁面换热系数,使仿真风速与实车30km/h行驶状态匹配。
4. 磁热耦合收敛性验证与误差溯源:用3组关键监测点温度对比锁定模型缺陷
4.1 在Fluent中布置3类温度监测点,覆盖热传递瓶颈区域
轮毂电机热失效常始于局部——非整体温升过高。必须在以下位置设置Surface Monitor:
- 绕组端部热点:定子绕组伸出轮辋端面的最外层导线表面(监测绝缘漆耐受极限)
- 永磁体背铁界面:永磁体与转子铁芯接触面中心(此处热阻最大,易引发退磁)
- 轴承座安装孔:轮毂轴承支撑座内壁(机械强度临界点,温升超120℃将加速润滑脂老化)
监测点需绑定到几何面而非节点,避免网格变形导致位置漂移。在Fluent中:Solution→Monitors→Surface→Create→选择对应面→Field Variable选Static Temperature。
4.1.1 收敛判据:温度残差与磁密残差必须同步满足
双向耦合求解中,Maxwell与Fluent各自有残差,但真正收敛需满足:
- Maxwell磁密残差
< 1e-4(相对值) - Fluent能量方程残差
< 1e-6 - 跨软件温度偏差:同一监测点在连续3次耦合迭代中,Fluent输出温度与Maxwell反馈的材料温度差异
< 0.5℃
若仅满足前两项而第三项波动>2℃,说明材料温度-磁特性曲线拟合不准或网格不匹配。此时应检查Maxwell中Temperature变量是否已关联到材料属性——常见疏漏是忘记在Setup→Model中勾选Enable Temperature-Dependent Materials。
4.2 实测对比验证:用红外热像仪定位热点,反向修正网格敏感度
某15kW轮毂电机实测发现:Fluent预测绕组端部最高温为112℃,红外热像仪实测为108℃,但热点位置偏移15mm。根源在于Maxwell中绕组端部网格过粗(平均尺寸1.2mm),未能解析导线束间的微尺度涡流。解决方案:
- 在Maxwell中对绕组端部区域施加
Inflation Layer网格,首层高度0.1mm,增长率为1.2 - 在Fluent中对应区域启用
Enhanced Wall Treatment湍流模型 - 重新运行耦合后,热点位置误差缩至2mm内,温升绝对误差降至±3.7℃
提示:网格无关性验证必须针对“热源区域”而非整个电机。对定子铁芯背轭进行网格加密(尺寸从2mm减至0.8mm)仅使温升预测变化0.3℃,但绕组端部网格细化可带来8.2℃偏差——证明热源定位精度远比全域网格密度重要。
5. 工程落地技巧:用参数化扫描快速定位轮毂电机热设计瓶颈,3个关键变量决定散热方案选型
5.1 基于DesignXplorer的3变量联合扫描:轮毂电机散热优化的黄金三角
轮毂电机热设计受制于空间、重量、成本三重约束,需在有限自由度内寻找最优解。在Workbench中启动DesignXplorer,对以下3个参数做拉丁超立方采样(LHS):
Wheel Rim Thickness(轮辋厚度):范围12–22mm,影响散热面积与结构刚度Ventilation Hole Diameter(轮毂通风孔直径):范围4–10mm,决定强制对流换热系数Thermal Interface Material Thickness(TIM厚度):范围0.1–0.5mm,控制定子与轮辋间接触热阻
扫描样本数设为45组(3参数×15),每组自动触发磁热耦合求解并提取3个监测点温度。结果生成响应面图,可直观识别:
- 当轮辋厚度<15mm时,绕组端部温升呈指数上升(热容不足)
- 通风孔直径>8mm后,永磁体界面温升改善趋缓(气流短路效应)
- TIM厚度>0.3mm时,轴承座温度反而升高(导热路径变长)
5.1.1 从响应面提取设计规则:用温度梯度比锁定主导传热路径
定义Thermal Gradient Ratio(TGR)=(绕组端部温−永磁体界面温)/(永磁体界面温−轴承座温)。TGR>3.0表明热量积聚在绕组端部,需强化端部散热;TGR<0.8表明热量滞留在永磁体区域,应优化背铁导热路径。某款商用车轮毂电机TGR实测为2.1,响应面显示:增大通风孔直径对TGR影响微弱(ΔTGR=0.05),而减薄TIM至0.15mm可使TGR升至2.8——证实界面热阻是当前瓶颈,优先选用导热系数≥6W/m·K的相变TIM材料。
5.2 快速验证散热方案:用Fluent瞬态热分析替代全耦合,节省85%计算时间
对已确定的散热改进方案(如增加轮辋散热筋),无需重复耗时的磁热双向耦合。可采用简化流程:
- 在Maxwell中运行一次稳态电磁场,导出各部件稳态损耗密度分布
- 将损耗作为恒定热源加载到Fluent瞬态热模型中
- 设置初始温度为环境温度,求解300s内的温升过程
该方法将单次仿真时间从12小时压缩至1.8小时,且对温升趋势预测误差<5%(因稳态损耗已涵盖95%以上热源)。关键在于:稳态电磁求解必须启用Advanced Core Loss Model(如Bertotti模型),而非默认的Standard Core Loss,否则铁芯损耗低估20%以上。
# Fluent中快速启动瞬态热分析的TUI命令序列 /solve/initialize/initialize-flow /define/models/unsteady? yes /solve/iterate 200 /report/fluxes/surface-integrals temperature wheel_outer_surface执行后,report/fluxes命令直接输出轮辋外表面平均温度,可与红外实测值比对——这是产线工程师每日快速验证散热改模效果的标准动作。
本文还有配套的精品资源,点击获取