
3.2 磁流变液制动器的温度场仿真分析
3.2.1 瞬态温度场仿真数学模型
磁流变液制动器在制动过程中内部热量集聚会导致其温度急剧上升,降低制动效果[2-4],基于ANSYS Workbench仿真平台进行磁流变液制动器瞬态温度场仿真分析。在磁流变液制动器中取一微元六面体,根据傅里叶定律可知在某时刻微元体六个面上的热量变化可表示为[5-7]

式中,dx、dy、dz分别为所选取微元六面体的长、宽、高;λx、λy、λz分别为微元六面体在x、y、z三个方向上的导热率;Qx、Qy、Qz分别为某时刻沿着x、y、z三个方向流入微元六面体的热量;Qx+dx、Qy+dy、Qz+dz分别为该时刻沿着x、y、z三个方向流出微元六面体的热量。
根据能量守恒定律,单位时间内该微元六面体的热量变化等于其内部产生的热量加上外部流入的热量,再减去流出的热量。结合式(3.2),微元六面体的瞬态温度场微分方程可表示为[8]

式中,ρw为微元六面体材料的密度;cw为微元六面体材料的比热容;Tw为微元六面体的瞬时温度;为微元六面体单位体积的发热功率。
磁流变液制动器处于制动工况下,其内部产生的热量主要来源于磁流变液的剪切发热和线圈的工作发热。制动盘与磁流变液之间剪切摩擦产生的热量可以等效于车辆减少的行驶动能,其是磁流变液制动器的主要热源。在车辆制动t时刻、Δt时间段内,车辆行驶动能的减少量ΔE(t)可表示为

式中,v0为车辆制动初速度;a为车辆制动减速度。
结合式(3.4),在车辆制动t时刻、Δt时间段内,车辆的动能损失功率Pm可表示为

根据欧姆定律可知,磁流变液制动器工作时线圈的发热功率Pc可表示为

式中,I为线圈电流;Rc为线圈电阻。
根据所设计线圈的结构特点,其电阻Rc可表示为

式中,N为线圈匝数;为线圈导线的平均直径;Sc为线圈导线的横截面积;ρc为导线材料的电阻率。
磁流变液制动器工作时,由于温差的存在制动器壳体表面与周围环境间存在对流换热。根据牛顿冷却定律,对流换热功率Pa可表示为

式中,ha为对流换热系数;Su为制动器壳体表面的换热面积;Tu为制动器壳体表面的温度;Te为周围环境的温度。
结合式(3.5)~式(3.8),根据能量守恒定律可得其功率表示形式为

式中,ci、mi、dTi/dt分别为制动器各部分材料的比热容、质量和温度变化率。