考虑气隙变化的高速电主轴热特性仿真

考虑气隙变化的高速电主轴热特性仿真
李特;芮执元;雷春丽;郭俊锋;胡赤兵
【摘要】为了解决电主轴高速旋转时结构变化引起的自身热特性改变的问题,提出离心膨胀和热膨胀会影响电主轴内定、转子间对流换热的观点.基于弹性力学理论,计算电主轴定、转子受离心力和热载荷而产生的径向膨胀量及由此导致的气隙变化量;根据对流换热理论得出离心力影响下泰勒数随转速和气隙长度的变化规律以及对流换热系数的变化规律.计算发现,气隙长度随定、转子的膨胀而减小,该减小量占设计值的37.7％;泰勒数随转速升高而增大,随气隙长度增大而减小;对流换热系数随气隙长度增加而降低,使得定、转子间传热受到抑制.结果表明,若能提高气隙间对流换热系数,适当提升油水冷却功率,能够大幅度降低主轴转子温升和热位移,提高加工精度.
【期刊名称】《浙江大学学报（工学版）》
【年(卷),期】2016(050)005
【总页数】8页(P941-948)
【关键词】高速电主轴;对流换热系数;气隙;温度分布;热变形
【作者】李特;芮执元;雷春丽;郭俊锋;胡赤兵
【作者单位】兰州理工大学数字制造技术与应用省部共建教育部重点实验室,甘肃兰州730050;兰州理工大学机电工程学院,甘肃兰州730050;兰州理工大学数字制造技术与应用省部共建教育部重点实验室,甘肃兰州730050;兰州理工大学机电工程学院,甘肃兰州730050;兰州理工大学数字制造技术与应用省部共建教育部重点实验室,甘肃兰州730050;兰州理工大学机电工程学院,甘肃兰州730050;兰州理工
大学数字制造技术与应用省部共建教育部重点实验室,甘肃兰州730050;兰州理工大学机电工程学院,甘肃兰州730050;兰州理工大学数字制造技术与应用省部共建教育部重点实验室,甘肃兰州730050;兰州理工大学机电工程学院,甘肃兰州730050
【正文语种】中文
【中图分类】TH122;TH133
电主轴的高速旋转不仅会导致大量生热,还会引发热变形,降低加工精度甚至导致主轴失效.该问题已占主轴总变形的40%～70%,成为影响机床切削性能的主要因素之一.因此对主轴热态特性的精确建模及预测成为设计的关键.
Bossmanns等[1-4]建立了电主轴系统生热、传热、散热的功率流模型及转子-轴承的热力耦合模型,奠定了热问题的研究基础.在热问题中,生热量和传热系数的计算是建立仿真模型的关键.当前电主轴的生热与损耗问题已得到较好解决[5],但因涉及较多因素,与传热、导热有关的边界问题始终难以精确计算,如结合面热阻[6]和气隙间对流换热系数.定、转子间热量主要通过气隙传递,对流换热是影响热力学精确建模和计算的主要因素[7].Yang等[8]基于神经网络法研究了空气对流换热系数与主轴热传递的关系；Li等[9]利用RBF神经网络法计算并优化了空气对流换热系数,得到了与实验数据最接近的最优值；王保民等[10]分析了主轴的产热和散热特性,指出温度场分布的非线性特征,但认为气隙间为自然对流换热；Haitao等[11]则认为空气在气隙中的流态与流过平板时类似,与实际不符.
随着转速升高,与上述边界条件有关的参数会改变,转子的离心膨胀和热膨胀引起的主轴结构的动态变化也会导致显著的动态热传递误差[12].Chen等[13]指出在2种膨胀的综合作用下,气隙的减小会改变对流换热系数,从而改变热传递和热分布；王保民等[14]分析了这2种效应作用下主轴转子的膨胀量及气隙减小程度,发现高速
时气隙减小量占设计值的很大比例；Chang等[15]监测了主轴转子的热膨胀和气隙的变化,指出主轴系统的耦合状况随膨胀量的增加变得更加复杂；Günther等[16]研究了高速转子的离心膨胀和涡动,描绘了膨胀量随转速的变化关系及对实验结果的扰动.上述文献都指出转子膨胀会造成动、热态特性的变化,但未说明这一变化对热态特性能造成何种影响,也未深入分析.
鉴此,本文基于某型号高速电主轴,从转子受离心力和热载荷产生的径向膨胀出发,讨论了定、转子间气隙的减小程度,计算了气隙对流换热系数的变化规律,分析了电主轴的热分布和热变形规律,提出了定、转子间传热效率的改善措施,旨在提高高速电主轴的传热性能与散热的控制水平,为有效的冷却和类似研究提供一定参考.
气隙变化量由2方面因素决定：作用于转子上的离心力和作用于定、转子上的热载荷.定子因静止只具有热膨胀,转子兼有离心膨胀和热膨胀,气隙变化量则取决于初始设计值与定、转子总膨胀量之差.
1.1 气隙变化量
电主轴为轴对称结构,可等效为厚壁圆筒,为计算简便,假设定、转子材料属性相同,二者的热膨胀方程[17]可写为
气隙变化量为
1.2 主轴生热
电主轴生热主要源于电机生热、轴承生热和空气摩擦生热,前两者生热占大部分,后者占比较小.
1.2.1 电机生热电机生热主要源于定子损耗和转子损耗,定子损耗又分定子铜耗和定子铁耗,转子损耗包括转子铜耗和转子铁耗,附加损耗较小可忽略,该部分生热可由电机电磁损耗关系得出[1,5].通常认为定子发热量占功率损耗的2/3,转子发热量占功率损耗的1/3[10].本文电机损耗为3 kW,可知定子和转子的损耗分别为2和
1 kW.
1.2.2 轴承生热[10] 角接触球轴承的生热主要源于滚珠与内、外圈及润滑液的摩擦,其发热量可由下列算式得出
因前轴承尺寸较大,定位预紧时在热诱导预紧力作用下生热大于后轴承,该差距随转速升高进一步增大,同时生热也随润滑液黏度的升高而增加.
1.2.3 风阻损耗因主轴转速很高,定、转子间的空气摩擦损耗较为严重,风阻损耗
[2]可表示为
主轴传热分为热传导、热对流和热辐射,热辐射因数值较小可忽略,热传导与结合部接触状态有关,热对流则取决于间隙中空气的流态.
2.1 结合面接触热阻
主轴内存在大量结合部,主要包括球与轴承内、外套圈,轴承外圈与轴室,轴承内圈与主轴转子,主轴转子与电机转子等.轴承内的热传递通过球与滚道的导热和润滑油气的对流进行.各结合面间的传热系数可由下列关系式确定：
2.1.1 球与轴承内、外圈热阻[10,18]为
2.1.2 轴承外圈与轴室热导[18]为
2.1.3 其他结合部轴承内圈与主轴、电机与主轴之间通常采用过盈配合,接触热阻可表示[6]为
2.2 对流换热系数
对流换热分为强迫对流换热和自由对流换热,在电主轴中,强迫对流换热包括轴承与润滑油气,定子与冷却水套,定、转子与气隙间的热交换等；自由对流换热包括转子端部、主轴壳体与周围空气间的热交换等.对流换热系数为
2.2.1 定转子气隙内的对流换热系数轴承油气润滑系统的气流会使定子和转子温升轻微下降,但因距离较远,影响程度较为微弱,因此认为定、转子气隙中无轴向气流[10].转子在静止的定子中旋转时,气隙内空气流态为泰勒库特流,它会大幅改变两者间的传热特性,当转子转速大于临界转速时,空气将会失稳,并产生二次流.因此对气
隙内的对流换热问题首先需判断流态,继而选取适当的计算方程.通常采用泰勒数(Ta=41.19)判断流态[21],但以临界转速判断更为直观.当气隙长度远小于气隙圆环平均半径时,临界转速为
式中：rm为气隙圆环平均半径,ν为空气的运动黏度.
如图2所示为当气隙δ0=0.3 mm时,泰勒数Ta随电机转子外表面直径d1和转速n的变化关系.由图2可得：泰勒数随转子直径和转速的上升而增大,在同一转速下,转子直径越大,泰勒数上升越快；若直径不变,转速是决定流态的主要因素.当泰勒数大于临界值41.19时,空气将失稳而产生泰勒涡.由图2可知,当转速大于
15 000 r/min时,不论转子直径如何,气隙内空气都将失稳.本文电主轴转子直径为82 mm,从图2中可知只要转速大于7 500 r/min,空气定会失稳,且转速越大失稳越严重.当计及转子的径向膨胀时,该失稳则有进一步增强之势.
式中：rm=(rsotor.e+rstator.i)/2； H=1-(rtator.i-rrotor.e/2rm),ωa为主轴角速度，g为重力加速度.
2.2.2 其他对流换热系数其他努赛尔数Nu及对流换热系数h见表1.表中 Pr为普朗特数,当介质为空气时Pr取0.707,为水时Pr取0.583；Dh=2(rstaor.i-rrotor.e)；Re为雷诺数,Re=v rrotor.e/ν,v为空气流速，l为传热面几何特征长度. 电主轴主要参数如表2所示,建立的三维模型如图3所示.因电主轴为轴对称,因此只需建立1/4模型,同时忽略螺纹孔,倒角等影响不大的几何因素,并将冷却水槽等效为环形槽.精确划分网格后有162 094个单元,290 652个节点,仿真设环境温度为
22 °C.
3.1 主轴膨胀与气隙变化
3.1.1 主轴离心膨胀、热膨胀及气隙减小量根据1.1的分析,离心力和热载荷都将引起主轴膨胀,为准确计算该变形,将离心膨胀与热膨胀单独计算.将主轴模型导入Ansys Workbench软件中,施加绕轴转速即可得到相应的离心膨胀位移,如图4所
示.从图4可见，该膨胀量随转速上升而增大,当转速最大时,膨胀量达到最大值.
与离心膨胀不同,在热载荷作用下主轴转子与壳体都会产生热膨胀.当主轴旋转时,生热与热传递、对流换热及散热等都为动态变化过程,所导致的热膨胀难以测量,由图4分析可知,当转速最大时离心膨胀达到最大值,为简便起见,计算20 000 r/min时所对应的热膨胀量.此时转子内、外表面温度分别为118和113 °C,转子轴内表面温度为115 °C,定子内外表面温度分别为70 和35 °C[14],由此可得壳体内表面和转子外表面的热膨胀量分别为25.2和82.2 μm,总膨胀量为107.4 μm,而该转速时离心膨胀量为5.7 μm,因此气隙减小量为113.1 μm,可见气隙减小程度主要取决于壳体和转子的热膨胀.只有在离心力作用时,气隙变化量占总宽度的3.25%；只有在热载荷时,气隙变化量占总宽度的35.8%；在综合作用下气隙减小量占初始宽度的37.7%；因此精确计算时该变化不可忽略.同时可预测：对于高DmN值高功率电主轴,高转速下发热和离心力引起的综合膨胀更大,气隙将急剧减小,若膨胀量进一步增大,将会引起转子摩擦定子,导致电主轴损坏.
3.1.2 气隙减小对风阻损耗的影响如图5所示为不同条件下空气摩擦损耗的变化,从图5中可见摩擦损耗随转速升高而增大,在同一转速下当考虑气隙变化时,摩擦损耗更大.以n=20 000 r/min时为例,不考虑膨胀时风阻损耗为21.08 W,只考虑离心膨胀时为21.48 W,只考虑热膨胀时为32.97 W,考虑综合膨胀时为33.49 W,与不考虑膨胀时相比,损耗增加了58.87%,这说明气隙变化对风阻损耗有重要影响,而热膨胀影响最大.一方面是因为该损耗与转子直径的三次方成正比,与气隙长度成反比；另一方面由于气隙的减小,速度梯度减小,空气流动更加不稳定,空气剪切力增大,做功增大,摩擦加剧,从而使得损耗增大.
3.1.3 泰勒数与对流换热系数如图6所示为泰勒数变化规律,Ta随转速升高而增大,当转速大于7 500 r/min时即会导致空气失稳,从而产生二次流.同时可看到,考虑气隙变化时泰勒数呈下降趋势,当n=20 000 r/min时,若只考虑离心膨胀,泰勒数
为107.54,只考虑热膨胀时泰勒数为56.95,考虑综合膨胀时则为54.44,泰勒数分别下降了2.8%,48.5%,50.8%.可见热膨胀导致的气隙减小对泰勒数的大小有决定性
影响.
如图7所示为对流换热系数的变化规律,从图7可见ha随转速升高而增大.在同一
转速下,考虑气隙减小时对流换热系数轻微减小,但与转速的作用相比,该变化的影响可以忽略不计.当n=20 000 r/min,不考虑气隙变化、只考虑离心力、只考虑热载
荷及考虑综合因素时对流换热系数分别为202.14、201.76、193.63、
193.10 W/(m2·K),可见对流换热系数随气隙减小而减小.该减小比例依次为
0.18%,4.21%,9.04%,因此热膨胀导致的气隙减小对对流换热系数有重要影响,并可能对温度分布造成一定影响.
3.2 主轴系统热分布
3.2.1 定、转子间的热分布边界条件对热分布的结果有较大影响,如8、9所示
为考虑与未考虑结合面热阻时的热分布图,考虑热阻后各结合部温度梯度增大,以主
轴转子和电机转子结合部的变化最为明显.电机转子温度由112.8 °C上升为
114.84 °C,上升了1.8%；主轴转子外表面温度则由112.8 °C下降至107.57 °C,下降了5.23%.主轴转子和电机转子间温度梯度由6.59 °C增大至7.27 °C,增长了10.3%,说明在热阻作用下,电机转子产生的热量不能顺利传入主轴转子.观察轴承内、外套圈分别与主轴转子和壳体的结合面,也可得出同样的结论,即当考虑结合面热阻后,轴承滚珠、内圈和外圈处的温度均出现热集中,温度值有所上升.这说明在电主轴热分析时不能忽略结合面热阻,否则将引起较大误差.
如图10所示为考虑气隙变化时的热分布,与图9相比可知,转子温度幅值有较大变化,每个温度梯度所对应的温度值有所上升.如电机转子表面最高温度较之前升高
3.34 °C,主轴转子外表面温度则上升了2.58 °C.这也印证了前面的分析,即气隙减小一方面增加了空气摩擦损耗,增加的损耗以生热的形式作用于转子表面,使其温度上
升；另一方面降低了对流换热系数,使转子与定子间的对流换热受阻,转子上的热量难以交换至定子上的水冷系统.两者综合作用的结果使主轴温升增加,这也与文献[22]的实验结果一致,进一步说明了分析的正确性.
3.2.2 转子温升与轴向热变形轴向变形能严重影响电主轴的加工精度,而径向变形对此类主轴影响不大,因此只考虑转子的轴向热变形.为进一步明确转子温度及热变形与对流换热系数和冷却功率之间的关系,绘制如图11、12曲线图.如图13所示为电机转子表面温度与对流换热系数关系,图中Pc为冷却水功率,t为电机转子表面温度,ha为气隙对流换热系数,由图11可知：电机转子温度随对流换热系数的增加而减小,对流换热系数每减小10 W/m2k,温度平均降低2.6 °C.如图12所示为轴向热变形量随对流换热系数的变化规律,图中δrotor.x为主轴轴向热变形量.从图12中可知δrotor.x随ha的增加显著下降,对流换热系数每减小10 W/m2k,热变形量平均下降2.47 μm.如图13所示为不同冷却功率下转子温度变化速率Δt与对流换热系数变化量Δha之间的关系,由图13可知温度并非等比例降低,降低速率Δt非线性变缓.如图14所示为主轴轴向位移变化速率Δδrotor.x与对流换热系数变化量Δha之间的关系,从图14中可知Δδrotor.x随Δha的上升而下降.从图13、14看出：若对流换热系数不变,提高冷却水功率并不能明显降低转子温度,这说明对流换热系数存在一个最优区间,否则若一味的提高该系数,将会造成过高的成本浪费.
温度分布的变化规律进一步说明：定、转子间的热传递主要靠气隙间空气的对流换热,因此如何提高对流换热系数则成为解决该问题的关键,若能提高,则转子温度会大幅度降低.可以效仿中大型电机进行轴向强制风冷,但目前该技术在电主轴中应用较为困难.若能风冷的基础上,一方面提高水冷效率,或者对转子单独进行水冷,则能从根本上改善电主轴转子的热变形特性.
(1)离心膨胀和热膨胀都会使得气隙长度显著减小,热膨胀的影响远大于离心膨胀；中小型电主轴的气隙长度较小,气隙变化明显,若其减小程度过大则会导致转子摩擦
定子,损坏电主轴.
(2)空气摩擦损耗随气隙减小而上升；对流换热系数随转速上升而增大,在同一转速下,该系数随气隙的减小而减小,影响程度不能忽略；其中离心膨胀的影响较小可忽略,但热膨胀对该系数的变化幅度有决定性影响,因此进行精确热分析时必须考虑该因素.本文未考虑齿槽形状、数目及转子表面粗糙度等因素,上述因素能否显著影响定转子间的对流换热还需更进一步研究.
(3)建模时必须考虑结合部热阻,否则会使仿真结果产生较大误差.
(4)提高冷却水功率并不能有效降低转子温升,而提高气隙间对流换热系数则会显著降低转子表面温度、温度梯度及转子轴向热变形量,此时若提高冷却功率,则能显著改善转子轴向热变形特性.。