基于损耗热网络模型的无感电机结温监测系统设计

做过电机驱动的同行应该都有体会绕组温度和功率器件结温这两个参数直接决定了一台机器能压到多少性能、能跑多久不烧。但真正做起来就会发现靠物理传感器去测温度要么测不到热点要么响应慢半拍要么在密闭压缩机壳里压根没地方装。我前阵子在某变频压缩机的项目上做了一整套“基于损耗热网络模型的无感控制电机结温监测系统”本质就是没有编码器也没有额外的温度传感器全靠实时计算损耗、推演热网络把绕组热点温度和逆变器结温给“算”出来。这套东西做下来感触很深今天把设计思路、模型构建、离散化实现和标定流程一次讲透希望对正在做电机控制或变频驱动开发的朋友有帮助。1. 为什么需要一个“无传感器”的结温监测方案1.1 物理温度传感器的三个老问题先说痛点。做过变频驱动的都知道电机绕组里面埋NTC或者PT100是常规操作但实际用起来问题很多。第一是响应慢。常规的热敏电阻探头外面有绝缘层加上热接触本身有热阻它感受到的温度变化往往会滞后几十秒甚至上百秒。可是绕组的热时间常数在小电机里可能只有几分钟堵转工况下几秒钟温度就能冲上去几十度。等你看到NTC报出过温绕组的热点可能已经处于很危险的状态了。对需要快速过温保护的场景传感器响应速度跟不上就非常危险。第二是测不到热点。绕组温度最恨的不是平均温度而是局部热点。尤其在端部绕组、槽内铜线这种地方热点可能比机壳温度和平均绕组温度高出几十K。物理传感器只能装在能装的位置它反映的永远只是“那个点的温度”不代表“最热的那个点的温度”。你要是按这个温度去做保护要么过于保守浪费性能要么过于激进留下烧机隐患。第三是安装和可靠性问题。压缩机、风机、水泵这类封闭式结构你想在绕组里埋一颗探头牵涉到工艺、密封、绝缘寿命成本一下子就上去了。而且传感器本身的引线在振动、高温环境下是故障高发点很多客诉最后查下来就是温度线断了、接触不良了。所以业内的共识越来越清晰能用算法估算的就不要依赖物理传感器。1.2 “无感控制”和“无感温测”是天然一对这就要说到系统里的另一个关键词——无感控制。无感控制的意思是电机转子位置和转速不是靠编码器或旋变反馈而是靠电流电压采样、用观测器推算出来的。这套方案的直接好处是省了编码器省了安装工时提高了整机可靠性。在空调压缩机、冰箱压缩机这种大量、低成本、封闭式场景里“无感化”几乎是必选项。无感控制省了“转速位置传感器”而无感温测省的是“温度传感器”。想一下这个组合一个变频压缩机系统没有编码器没有绕组温度探头但依然能精确知道转子位置和转速依然能实时知道绕组温度和逆变器结温。这不光是成本问题更是整机可靠性的巨大提升。传感器少了故障点就少了器件少了成本就下来了系统还更皮实了。1.3 这套系统的目标定位说白了这个系统要做的事情就是在有高精度电流采样和无感控制算法的基础上增加一层“热状态观测器”。它实时计算各个损耗源的功率把这些损耗当作热源输入一个等效热网络模型再求解一组热阻热容网络的微分方程最终输出绕组热点温度、永磁体温度、逆变器IGBT结温的实时估算值。我在项目上定的设计目标是稳态误差控制在±8℃以内瞬态工况下最大误差不超过15℃保护动作值预留10℃安全余量。实际跑下来稳态误差能收敛到5℃左右整体还是满意的。2. 方案选型为什么是损耗热网络模型2.1 有限元仿真的精度和实时性矛盾开头可以先把几个候选方案放在一起对比方便理解为什么最终选了热网络模型。有限元仿真是最“物理正确”的方法磁场和温度场耦合着算精度可以做得很高。但问题在于计算量太大。一个3D网格模型求解一次瞬态温度场在工控级芯片上几乎是不可能完成的任务在PC级硬件上也要跑几秒到几分钟。这些年也有人尝试用降阶模型ROM把有限元结果压缩成快速查表但建模工作量很大换一个电机规格基本就得重新做一遍。对大批量生产的家电控制器来说这不太划算。另一个方案是纯数据驱动把全工况范围做大量DOE实验测温度、电流、转速、电压然后拟合一个多项式或者神经网络模型。这个方案精度也可以但问题一样在标定成本上。全工况矩阵意味着海量的温升试验每一台电机规格都要重新标开发周期拖得很长。而且基于数据拟合的模型在标定范围之外 extrapolate 的时候很容易“放飞自我”可靠性存疑。2.2 损耗热网络模型的位置算力、精度、标定成本的平衡点损耗热网络模型走的是第三条路先实时计算各类损耗再用等效热阻Rth和热容Cth构成的网络来描述热量从热源到冷却介质的传播路径。它的物理基础是传热学中的集总参数法虽然做了大量简化但在工程精度范围内完全够用。这个方案的好处体现在三个维度算力开销小。实时计算的核心只是一组一阶RC微分方程的离散化迭代在任何一个带FPU的MCU上跑一遍用不了多少微秒。物理意义清晰。热阻决定稳态温升热容决定瞬态响应模型的每个参数都有明确的物理含义调试时有方向感。标定周期可控。做一次堵转温升试验和几组阶跃工况标定就能把主要参数都辨识出来一个规格的电机几天内可以搞定。你可以把热网络想象成一个由管路和水箱组成的水网损耗功率是水源热阻是管径管径越细水流越慢、阻力越大对应温升越高热容是水箱容积水箱越大蓄热能力越强温升变化就越缓慢。电机内部的热量就是从这个水网一路流到机壳表面最后散到空气或冷却水中去的。2.3 典型的电机热网络拓扑谈到具体拓扑以我做的内置式永磁同步电机IPMSM为例整个热网络分成两条主散热路径第一条路径是“定子侧”绕组铜耗→绕组绝缘层→定子齿部→定子轭部→机壳→外部冷却介质。这条路径承载了电机发热的大头也是我们最关心绕组热点温度的原因。绕组到机壳之间的热阻往往占了整个路径的大半部分。第二条路径是“转子侧”永磁体涡流耗和转子铁耗→转子铁芯→转轴→轴承→端盖→机壳。这条路径的热阻通常比定子侧大得多所以在大多数工况下转子温升都比较高但散热很差。在高速或者过载工况下永磁体温度一旦超标就有不可逆退磁的风险所以转子温度估算同样很重要。有些系统还会把“端部绕组到端盖内腔”的散热路径建进去因为在自然冷却条件下端部绕组通过空气对流散掉的热量比例相当可观不能忽略。我建议在模型里把端部绕组单独做一个节点跟槽内绕组分开处理这样对热点温度的估算精度会明显更好。3. 损耗计算与热参数辨识精度从哪来模型框架搭好之后最难的是两件事损耗算得准不准热参数标得准不准。这两件事直接决定了整套系统的估算精度。3.1 损耗实时计算不是简单I²R热网络模型的输入是损耗功率损耗算不准后面全白搭。我把损耗分成三块来算。第一块是铜耗。铜耗是最主要的热源计算公式是P_Cu k × (Id² Iq²) × Rs(T)。注意这里面的Rs不是常温电阻而是要随温度修正的。铜的电阻温度系数大约是0.0039/℃也就是说绕组温度每升高100℃电阻增加将近40%。如果不做温度修正高温工况下铜耗会被低估很多温升估算就会偏低这在保护场景里是很危险的。第二块是铁耗。铁耗的工程计算常用Bertotti模型P_Fe kh×f×B² ke×f²×B²前半部分是磁滞损耗后半部分是涡流损耗。关键在于磁通密度B的获取。在无感控制下可以通过d-q轴电压方程反推磁链再结合电机参数表查出对应的B值。考虑到查表的插值精度和参数误差工程上我通常会在算出的铁耗上乘一个1.1到1.3的修正系数用温升试验去反推。这里有正反馈的问题需要处理温度高了→电阻变大→铜耗增加→温度更高。如果不加任何限制估算模型可能发散或者震荡。所以我在迭代计算里做了阻尼处理——用上一拍的温度估算值来修正当拍电阻同时限制单拍温升变化率不超过物理上限。实际测试下来效果很稳。第三块是机械损耗包括轴承摩擦和风摩损耗。低速工况下这个量很小可以忽略但高速电机或者长时运行时还是要算进去。最简单的办法是测空载功率减去空载铁耗当作机械损耗然后按转速的平方关系做拟合。3.2 热参数辨识堵转试验和直流电阻法热阻和热容这些参数靠看铭牌是拿不到的必须通过专门的试验去标定。我用的主要方法是堵转温升试验。具体做法是锁住转子给绕组通一个稳定的直流电流让电机处于纯粹的电加热状态。此时输入的电功率P基本全部转化为铜耗。用直流电阻法监测绕组温度大家知道铜绕组的电阻和温度呈近似线性关系先测常温下的绕组电阻R0然后在试验过程中不断给绕组通一个很小的测量电流量出当前电阻R通过公式T (R/R0 - 1)/α T0 推算当前平均绕组温度。试验流程大概是让电机在环境温度下充分静置记录环境温度T_amb和冷态电阻R0。锁住转子通电加热每隔一段时间记录电压电流和绕组电阻计算出当前绕组温度T_winding。同时记录机壳表面温度用热电偶贴在机壳外壁。持续加热直到温度趋于稳定判断标准一般是10分钟内温升变化小于1℃。断开电源记录自然冷却曲线一样每隔一段时间记录电阻和温度。有了加热和冷却这两条温升曲线就可以做参数拟合了。一阶RC网络的温升响应是指数形式ΔT(t) P×Rth×(1 - e^(-t/(Rth×Cth)))。从稳态温升可以直接求出Rth ΔT∞ / P从温升曲线的63.2%上升段对应时间能求出时间常数τ Rth×Cth进而算出Cth。这里分享一个容易踩的坑堵转试验加热的是整台电机一套简单的单级RC往往拟合不好因为绕组、铁芯、机壳各自的热惯量差异太大。所以我实际建模时用了两阶甚至三阶RC网络第一阶代表绕组到铁芯的快速热传递第二阶代表铁芯到机壳第三阶代表机壳到环境的慢速散热。拟合的时候分段拟合先拟合长时段的慢速段再回头拟合短时段的快速段这样比一次性拟合整体曲线要准确得多。3.3 一个具体的计算例子拿我们项目里那台2.2kW的IPMSM来说常温下绕组电阻Rs0 0.8Ω额定电流8A。堵转时通8A直流电铜耗大约P_Cu 8²×0.8 51.2W。标定出来的绕组到机壳热阻Rth_cw 1.2 K/W绕组热容Cth_w 300 J/K。那么稳态绕组温升ΔT∞ 51.2 × 1.2 61.4K。如果环境温度40℃绕组稳态温度就约101.4℃对于常见的F级绝缘耐温155℃来说余量已经不算大了。要是这时候不限制电流持续堵转温度继续往上走绝缘老化就会加速。这个计算看着简单实际上很有指导意义它直接告诉我们保护阈值定多少、降频策略在哪里介入合理。比如我设定的过温预警阈值就是95℃在90℃时开始线性降额到110℃强制封锁逆变器给硬件保护留出缓冲。4. 嵌入式端实现从连续方程到控制环代码4.1 控制周期与离散化方法热网络方程本质是连续的一阶微分方程组到了MCU里必须离散化。这里需要谨慎选择离散方法。以最简单的单级RC为例连续域方程是Cth × dT/dt P_total - (T - T_amb)/Rth。我使用的是后向欧拉隐式欧拉离散法。之所以不用前向欧拉是因为热时间常数非常大数百秒但控制周期只有100μs10kHz载波两者相差七个数量级前向欧拉在这种尺度下数值上虽然稳定但计算步长选择稍有不慎就会出现数值刚性。后向欧拉的稳定性是无条件稳定的不存在这个问题。离散化推导过程不复杂T(k1) T(k) Ts × [P(k) - (T(k1) - T_amb)/Rth] / Cth。整理一下得到T(k1) (T(k) Ts×P(k)/Cth Ts×T_amb/(Rth×Cth)) / (1 Ts/(Rth×Cth))因为Ts/(Rth×Cth)这个量级大约是0.0001计算中不用担心溢出用单精度浮点就足够了。系数Ts/Cth和Ts/(Rth×Cth)可以预先算成常量每个中断周期做两次乘法和一次除法。在具体实现中我对每个热节点都维护一个状态变量存储当前温升而不是绝对温度这样在环境温度变化时只需要在输出端加一个偏置避免整个模型状态量受环境温度波动干扰。4.2 无感控制下如何获取高频电压电流热网络模型的输入是损耗损耗计算的输入是电流和电压。这里的难点在于无感控制下没有编码器电压、电流的获取精度直接决定损耗计算的精度。电流是比较容易处理的在FOC控制中每个PWM周期内通过采样电阻或霍尔传感器采集三相电流在特定的PWM开关状态通常是在载波谷底同步采样能获得噪声比较小的电流值。需要注意采样时刻要避开开关管的开关瞬态和续流二极管反向恢复期不然后面的损耗计算会注入大量噪声。电压则有两类做法。一是直接采样母线电压然后利用SVPWM的占空比重构三相电压。这里有个工程细节SVPWM输出的参考电压与实际输出电压之间存在死区时间、管压降、母线电压波动等的误差需要在低电流、低速工况下做死区补偿和管压降补偿。如果不处理损耗计算在轻载工况下可能出现负数导致模型输出的温度估算值明显偏低。二是无感控制本身输出的电压矢量是电流环经过PI调节生成的d/q轴电压命令。这个电压命令在稳定工况下很接近实际电压但在动态过程中会有偏差。所以我的做法是两者结合稳态时以电压命令为准动态过程中做低通滤波过渡既保证响应速度又保证稳态精度。实测下来经过死区补偿后重构电压和真实电压的误差可以控制在2%以内对应的铜耗估算误差约4%热网络最终输出的温度误差在几℃以内完全满足保护需求。4.3 软件架构与计算开销整个系统的软件架构可以分成四层最底层是采样层负责三相电流、母线电压、环境温度可选的采样和滤波。第二层是控制层无感观测器滑模观测器或扩展卡尔曼滤波估算转速转子位置电流环计算d/q轴电压命令。第三层是损耗计算层基于控制层的电流电压输出实时计算铜耗、铁耗和机械损耗。第四层是热网络求解层把损耗功率送入离散化的热阻热容网络迭代求解各个节点的温升。我用的是某款工业级MCU载波频率10kHz热网络求解在这个中断里执行。整个热网络模型如果是三节点两阶RC网络一次迭代运算大概需要30-50次浮点运算在任何一个带硬件FPU的MCU上占用的CPU时间最多几百纳秒对控制环的影响可以忽略不计。另外建议开启一个1Hz的低优先级任务负责把温度估算结果做统计、日志输出和上位机通信这样不占用高频中断的资源。5. 标定流程与实测数据5.1 实验平台与工况设计整个系统搭完后验证环节是重中之重。我在某变频压缩机项目2.2kW IPMSM10kHz载波上做了完整的标定和验证。测试平台包括被测压缩机、测功机或扭矩传感器、直流电源、功率分析仪、预埋热电偶用于对照验证、一个可编程风冷装置。我设计的验证工况覆盖了这些典型场景额定工况转速3000rpm额定负载连续运行。堵转工况锁住转子通额定电流记录温升过程。过载工况1.5倍额定电流持续60s模拟加速和冲击过程。冷却退化工况人为关闭外部风扇模拟换热器堵塞或风机故障。预埋热电偶是验证的重要参照物但不参与控制。一套系统里同时存在“估算路径”和“实测路径”才能公平评价精度。5.2 标定的详细步骤参数标定我一般按下面这个顺序走每一步都会记录并校验防止误差传递第一步常温电阻和温度的基准确认。把电机放在恒温箱里分别在20℃、40℃、60℃、80℃下测绕组直流电阻拟合出准确的R0和α。这个步骤不用求快数据准了后面全准。第二步堵转温升试验。锁住转子通直流额定电流按前面提到的方法记录绕组温度动态和机壳温度动态一直做到热平衡。利用稳态数据拟合热阻Rth利用瞬态曲线拟合热容Cth。这里特别注意一次堵转试验只能提炼出一条散热路径的参数如果电机还有额外的冷却方式比如轴流风扇需要分别在“有风”和“无风”两个状态下各做一轮试验。第三步动态工况验证。用正常FOC驱动电机跑额定、过载、加减速工况把热电偶实测值和估算值做对比记录偏差分布。第四步如果需要更高精度做铁耗修正系数的反推。做法是在额定转速、中低负载下用实测温度和热网络反推总损耗再减去已知的铜耗剩下的部分和理论铁耗做比较求出修正系数。5.3 实测数据对比下面贴一组我印象比较深的实测数据工况是从额定负载阶跃到1.5倍过载持续运行120秒。时间s估算绕组温度℃热电偶实测℃偏差℃052.351.80.53087.679.28.460103.495.18.390112.8107.55.3120119.6114.94.7180124.2121.03.2能看出一个规律阶跃初期估算值跑得比实测快偏差最大到8℃左右但随着时间推进逐步收敛到5℃以内。原因是热网络模型里的绕组热容比真实物理绕组略小导致瞬态响应略微偏快。这在保护场景下其实是偏保守的估算温度比实际温度高安全方向是对的所以我没有强行修正这个偏差。5.4 其他工况的验证结果堵转工况是整个系统最严苛的考验。因为堵转时没有反电动势电流全部转化为铜耗温升速度非常快。实测堵转600s估算绕组温度121℃热电偶实测115℃误差6℃稳态偏差保持在8℃以内。考虑到堵转时冷却条件最恶劣这个精度足够支撑保护决策。冷却退化工况也就是人为关掉风扇让电机自然冷却运行。这个工况下热阻Rth变大如果模型继续用原来强制风冷状态的热阻估算值会明显偏低。我们的解决办法是检测到机壳温度上升速率超过某个阈值时自动切换热网络参数到“自然冷却”模式。切换后估算偏差从15℃以上缩小到7℃左右。6. 常见问题与排查技巧6.1 温度估算滞后导致保护延迟这是大家最容易忽略的一点。热量从绕组传到机壳需要时间而估算模型如果只用慢速路径动态升温时估算温度会明显滞后于实际热点温度。我的处理办法是在模型里加入一个“快速热积累通道”用电流平方的瞬时值做一个带温度系数的积分项在过载开始的瞬间快速抬升估算温度。它表征的是绕组绝缘层附近的热点温升响应速度远快于绕组平均温度。快速通道和慢速主通道做加权混合既保证动态响应又不牺牲稳态精度。6.2 堵转时风机停转冷却条件突变堵转和散热系统失效是关联性很强的复合故障。如果只靠热网络模型内的参数切换还不够建议在策略上加一条简单粗暴的规则当检测到电机处于堵转状态且机壳温度上升速率异常直接按最恶劣冷却条件无风自然冷却的模型参数来计算同时缩短保护时间阈值。宁可误降频不可放过热。6.3 冷启动时初始温度未知系统上电时绕组温度到底是多少在完全没有传感器的情况下是一个未知数。如果你直接初始化成环境温度可能在上一轮停机前电机刚经历过高温一上电估算值就偏低了。解决办法是在允许的场景下做“上电绕组温度辨识”低速零电流运行状态下注入一个短时高频小信号或直流偏置测量绕组电阻反推初始温度。如果应用场景不允许这种注入那么另一个替代方案是做一个“历史温度记忆”断电前把估算温度存入非易失存储下次上电读取并和环境温度取较大者作为初始化值。这个方法简单实用实际效果也不错。6.4 电流传感器偏置误差对温度的影响电流采样偏置哪怕只有几十毫安在轻载时也会导致损耗计算出现明显百分比误差最终反映到稳态温升估算偏差上。我踩过这个坑某次实验发现轻载下估算温度始终比实测高10℃排查到最后是电流传感器零点漂移了30mA。所以建议在每次正常运行前做一次零电流校准或者每几分钟做一次在线偏置校正等效于在零矢量状态下的电流平均值确保损耗计算用的电流是干净的。6.5 从电机模型扩展到逆变器结温估算最后说一个扩展经验。这套热网络思路不光能用在电机本体上还可以直接复制到功率逆变器上。IGBT或MOSFET的损耗模型开关损耗导通损耗可以用电压电流信息和器件手册的热阻参数代入建立从芯片结到散热器底座的热网络。我在这套系统里同步实现了逆变器IGBT结温估算用于过温保护和寿命预测。当你把电机绕组和逆变器结温都纳入同一套“损耗热网络监测”框架后整台变频系统的热状态基本就全在你的掌控之中了。我个人在实际操作中的体会是热网络模型不算新东西但想要在无感控制这种体系里把它做好真正难在“损耗算得准”和“参数标得对”这两个基本功上。很多人一上来就追求模型复杂堆节点堆到十阶八阶结果参数标不出来误差反而更大。先做一个三节点左右的模型用试验把参数磨准再根据实测误差一点点加复杂度这才是稳的路线。另外估算系统永远要和硬件保护做双冗余估算负责性能释放和趋势预警硬件传感器负责安全底线这套组合拳打下来电机系统既跑得狠又烧不了。