高速精密主轴轴承热特性的计算及分析

HEATS

高速电主轴是高速机床的核心部件，电主轴单元各零件的精度通常为µm级，如跳动一般为1～3µm，主轴刚度一般为100-300N/µm，因此，电主轴各零件自身的精度误差和变形量很小。然而，高速运转时，轴承发热量很大，导致轴承温升很高，并引起热变形。热变形引起的误差远大于精度误差和变形量。而且，热变形直接改变了轴承的预紧状况，影响轴承的刚度特性和电主轴的加工精度，严重时，甚至导致轴承的热咬合，使电主轴毁坏。电主轴的发热问题降低了其在实际生产应用中的可靠性。故对轴承进行热分析并掌握轴承的热态特性对提高电主轴加工中的稳定性和可靠性都是极其重要的。高速精密角接触球轴承是机床主轴单元中应用最广泛的一类轴承，本文就以此类轴承为研究对象，对高速精密主轴轴承系统作了较有意义的热特性分析与计算。

1 高速精密轴承的摩擦力矩及摩擦热

: h& B7 A) b" R2 \

轴承的发热主要跟摩擦力矩有关，力矩越大，生成的摩擦热也越多。根据Palmgren推导的公式，轴承的摩擦力矩主要有两部分:润滑剂粘性产生的摩擦力矩和与速度无关的载荷作用产生的摩擦力矩。计算公式如下: 8 N% K+ ]% x4 b6 r( l: L- X

2 R( V9 E9 K9 n7 U: Q' j

M0=10-7f0(vn)2/3dm3	(1) 1 S+ H- u3 X3 z
M1=f1P1dm + W* F( L- r. z6 d' s1 e	(2)

式中:f₀为与轴承类型和润滑方式有关的系数;v为工作温度下润滑剂的运动粘度，mm²/r;n为主轴转速，r/min;d_m为轴承节圆，mm;f₁为与轴承类型和所受负荷有关的系数;P₁为确定轴承摩擦力矩的计算负荷,N;M₀,M₁为摩擦力矩，N·mm。 . ~5 Z7 e5 J" H g

轴承接触区热生成还跟滚动体与轴承套圈的自旋运动有关，高速下根据轴承套圈沟道理论，滚动体和内套圈存在自旋运动，自旋摩擦力矩为

- c0 y5 c6 G! F+ q: R k7 l& O, `- Q( h0 _% G6 j7 w2 K8 b0 M0 A8 _5 P1 q9 F

, s6 s; X9 g0 Y% i9 ?. c& }0 {6 l" o$ U; H0 q* o% ?. @ A% u( Z- m Msi= 3µsiQiaiei 0 x) d/ E% i" F. d+ u 0 o; D4 {1 r+ m2 T# z' ^2 b 8 ; C4 {# R: V0 s" A$ o# D

(3)

式中:µ为滚动体与沟道的摩擦系数;Q为滚动体与沟道的法向接触载荷,N;a为赫兹接触椭圆长半轴，m;e为第二类椭圆积分。

因此， 轴承的摩擦热可按如下计算得到

+ F- J$ s2 Z" ]& G5 V2 o! h: P) b! `3 ~/ P; b! J6 ^% Q

( `8 Y% l9 k& B. W# t) u( a* k4 _2 G4 O0 h# E% f; \ \8 J. G6 @( m: r5 I7 F& K7 I* E/ c) I% h Hlot= - b& ~1 | I1 X$ K: G8 b) K 2p $ j( }* z$ [+ o3 n n(M0+M1)+Msiwsi 2 o% K% R6 q t$ @# S1 A) d 0 d/ E( M! J$ ]% n% g* e: [3 u 60 ( N; l* G4 X E1 c D, l1 G

(4)

式中：n为主轴的转速，r/min。w_si为滚动体内圈自旋运动角速度，rad/s。

由式(1)～式(4)可知，要求得摩擦热须计算滚动体的自旋角速度，该角速度计算公式如下：

! p3 d: O ~! a( t& q0 L) r( l, n( _/ h- N

wsi=-wbsin(tgb-ai)+(w-wc)sinai	(5) 9 ^; [, o% E: V& q& h0 [' f6 P
" ?# G3 @- y, J+ E# i A0 z ]9 ?6 A6 Z7 i5 U$ g' r7 K! ^4 t! `: R5 o! V8 H* }0 f wb=w[( % h: c3 h/ k! }8 I, Y0 G' }7 q cosae+tgbsinae + 2 u% E- X5 C" r! \ cosai+tgbsinai )g'cosb]-1 , w, w8 N) T# r, f, N3 S : h. g- c3 b/ h# }3 S5 g8 o 1+g'cosee & |$ _# _3 y) J2 O ^ 1-g'cosei ( T2 y' N) g( N8 c# @: K1 k	(6)
wc= . Q7 M& J7 `8 ~0 f w(1-g'cosei) 0 A5 w% w& ? l 4 K$ W8 X. Z5 Y. \ 1+cos(ai-ae) 1 U( Y$ r; k' D; n7 }( H8 W	(7)

式中：w为主轴的角速度，rad/s；a_i，a_e为轴承内外圈接触角，rad；w_b为自转角速度，rad/s；w_c为滚动体公转角速度，rad/s ! C$ x! g! W0 {1 u4 D

g'=D_b/d_m
式中：D_b为滚动体直径，mm。

" E" Y; Z8 d, |. v

0 Y/ J- u m4 T5 m" w$ `: s1 k6 @- G# I5 Y: w" d7 ~9 D h# i9 n! j$ ? W# N tgb= sinae & n8 Q2 G9 U. ]% y, R) L) x cosae+g' & c% S' P: Z# {- m& O+ x+ C) ] o

(8)

4 U4 I% I q8 v! E2 I

2 “主轴一套圈一滚球一轴承座”系统的热分析

- C: O3 T. r2 H, m5 e, w$ T7 `% y3 C1 f6 N: w' Z; T- W( f1 O% J, g

图1 轴承系统温度节点分布图

将轴承、主轴和轴承座作为一个系统考虑，该系统的摩擦热生成和热传递是复杂的三维热问题，为简化计算，假定系统周向的热特性是一致的，则可以轴向截取一个平面为研究对象作系统的热分析。根据Burton 和steph提出的观点，将摩擦热在轴承滚动体和套圈之间作1:1的分布，即摩擦热有一半进入滚动体，另一半则进入套圈。采用节点网络法求解该系统各关键点的温度，得到整个系统的温度分布状况。这里，一共设立了9个温度节点，满足了一般计算的精度要求。各关键温度节点的分布如图1所示。轴承润滑方式为油气润滑方式，强迫对流系数比较难确定，这里根据经验公式计算强迫对流系数。 3 K$ K* @2 q7 E3 O) w

2 _' Q% O) f1 D+ }+ \0 ~) `& I- p/ Q" M/ u: a$ d+ J+ V6 Z0 o# S$ K6 T# M w

# q! R% n0 K4 H2 \' B3 e) R* C8 R' _0 s h=0.332Kl·Pr1/3·( u * S3 E6 }( k/ [ )0.5 ux

(9)

式中：K_l为润滑油的导热系数，W/(m²• ℃)；P_r为普朗特数冰为对流的速度，m/s；v为润滑油的运动粘度，mm²/s; x为特征长度，m。

对滚动体，对流速度取轴承保持架的表面速度，特征长度取轴承的中径；对内外套圈，对流速度各取保持架表面速度的三分之一，特征长度则取各自的沟道直径。

系统处于热平衡时，各节点的温度不再发生变化。此时，各节点的净热流为零。根据这一点，可以得到9个热流方程，并组成一个包含9个节点温度的方程组。该方程组是个非线性方程组，采用牛顿—拉费逊法求诸计算机可解。

下面给出轴承套圈、主轴径向的导热方程。

v3 p3 a8 ]. U9 A8 P0 C; t: O+ e# Q3 R2 ?4 K( s% O4 C* j. a6 m. i, L# b

( O3 M9 X; e4 Q/ I; D) I( \0 w% X+ A' E7 M8 i& r, k( a1 ~( m# ^ @: o5 o3 Z# \" |1 e) ^6 [4 a3 e$ v& ~& ^$ K/ {. { q= 2pKC 8 C1 B3 m; U1 G7 v* G# b Te-Ti ln(de/di)

(10)

式中：K为导热系数，w/(mK)；C为轴承宽度，m；d_e、d_i为各套圈的外径、内径，m。 & b# O+ g; }* Q, d

主轴径向导热方程为 ; ?3 p& K) A; {4 O% b3 h7 B* k6 j

+ L( H1 A0 E' ?9 |; N, W0 N, b7 _. }9 d8 m6 Y% }$ y9 `7 u' ^& Y R% C4 t" }- r, K

q=KpW(T0-T∞)

(11)

式中：W为选定的主轴宽度，m。 7 H) s4 z6 K6 L) ^* e# H& `/ D* d

3 影响轴承温升的因素分析

9 n. K( `4 m# T' e I

本实例的轴承型号为7013C ，轴承及计算模型参数为：

d=65mm D=100mm , D₀=160mm，B=18mm，L_s=30mm，L_h=20mm。 3 `/ H5 s) k* V6 c( B

已知轴承预紧方式为定压预紧，预紧力为300N 。主轴转速即轴承内圈转速为20,000r/min 。计算得：高速下，轴承内外圈接触角分别为24.3°和7.2°，内外圈法向受力各为729.26N和2391.5N。轴承润滑方式为油气润滑，轴承保持架表面速度为18.9m/s 。计算结果在表1中给出。 3 f G% Y$ j1 u

) a6 f6 S u1 O0 t' ~' Y8 U9 O, \$ [( m) g, u! H: _3 g2 t' E: U5 @1 F) g: ~" ?

表1 节点温度

节点	T1	T3	T4 * {9 [8 F3 w7 Q. J	T5	T6 2 x( `* P. o l/ Z8 X+ l4 C# G	T7 9 }7 ]- b" @1 D! W! m( D' U0 U	T8
温度	50.2 Z+ S) i9 ]) S% ]6 z	56.1	57.2	49.3	37.1	34.9 8 i1 d8 e* X. D; t! W	31.3

# R: Y( ]# A3 K5 _
1. 转速的影响
   轴承的转速对轴承温升影响显著，轴承转速越高，摩擦发热越严重，大量的摩擦热导致轴承的温升也越明显。图2是轴承在不同转速下的温升情况。可以看出，转速越高，温升也越高。 : ~ t/ _# p: s. A) i7 C, V5 q! Q. v) q P# e9 \% s9 ~! B5 D

   图2 转速对轴承节点温度的影响 ' n5 M z% P; H ^& D& p6 a% o

   图3 润滑方式对轴承温度的影响

2. 切削力的影响
   刀具切削工件过程中产生的切削力使轴承承受额外的附加载荷，这些载荷对轴承的摩擦热也有影响，间接影响了轴承的温升。切削力主要影响轴承由外载荷作用产生的摩擦力矩，但是计算表明，该影响值很小，因此，切削前后，轴承的温升变化并不明显。
   ; R7 y) m6 z: h2 `# J# R
3. 润滑方式的影响
   油脂润滑适合转速较低的情况，其d_n值最高一般不超过60万，相比之下，油气润滑能适应更高的转速，其d_n值最高可达200万，是一种可靠的高速润滑方式。图3反映了油气润滑和油脂润滑方式对轴承温度的影响情况，从该图可以看出，油脂润滑在转速升高时，轴承温度攀升很快，说明不适宜高速润滑；油气润滑时，随转速的提高，轴承温度升高比较平缓，润滑冷却效果很好。

4 k: o' V' M' }- i8 B, \& h7 }# B

4 结论

* f( m- m$ C6 s9 y+ h

通过以上分析计算，可以得出如下结论：

1. 转速越高，轴承的摩擦力矩越大，生成的摩擦热也越多，轴承的温升也就越高；
2. 在正常工作范围内，切削力对轴承温升造成的影响很小，考虑更多的因素对轴承摩擦力矩的影响，会使计算结果更符合实际。
3. 对于高速轴承(d_n值在200万的范围内)，采用油气润滑可以很好的抑制轴承的温升。