论述磨削内圈沟道圆度的几何特征及影响因素

btpit99

# e# P6 `+ N; R- s6 O图1 轴承内圈沟道外圆
滚动轴承内圈沟道圆度是一项重要的质量指标，直接影响轴承的工作精度、平稳性和使用寿命。磨削通常是内圈沟道的半精加工或精加工工序，对成品内圈沟道的圆度起决定性作用，内圈沟道磨削的圆度除了取决于磨削工艺系统的精度和动态特性外，还与磨削的工艺参数密切相关。前者，已为许多学者所重视，并进行了大量的研究，从而大大改善了磨削内圈沟道的圆度。随着磨削工艺系统精度的提高和动态特性的改善，磨削工艺参数成为影响圆度的主要因素，因此研究磨削工艺参数对圆度的影响规律，对合理选择工艺参数，确保轴承沟道的磨削圆度，进而实现磨削轴承沟道的工艺参数优化，有重要意义。图1 是磨削轴承内圈沟道截面实际轮廓与理想轮廓的误差情况，可表示为几何形状误差
5 [" F( {0 X  t. y?Rw（F）=Rw（F）-Rw0
（1）
式中，Rw（F）为实际工件轮廓半径，Rw0为理想工件轮廓半径。图中，O1、O2分别是工件的几何形状中心和测量回转中心。
轮廓误差?Rw（F）是F的周期性函数，可用付立叶级数表示成
* p" d. V6 M' D9 {7 x! z- L?Rw（F）=
& k" [; S, f  N9 R* K' I5 Ax0
: S4 @+ N; i9 n; O+∞Σi=1
) ]7 ~  ~# w: |' }$ |xicos（iF+Fi）
2
7 T" M% H7 ]% N1 {. ?式中:xi（i=0，1，……）为谐波幅值：F、Fi为谐波相位角：x0/2代表加工尺寸误差，也是误差函数Rw（F）的平均值。式（2）中，一次谐波x1cos（F+F1）对应内圈沟道外圆几何中心相对测量回转中心的偏心：二次至十五次谐波xicos（iF+Fi）（2≤i≤15）对应的就是轴承沟道外圆的圆度。
2 圆度的影响因素分析
- `) U) y' O( T) ~# _5 A9 b7 U轴承内圈沟道通常采用变进给速度切入磨削方式，其磨削循环可分为:工件快速趋近、粗进给、细进给和无进给磨削四个阶段。影响磨削轴承沟道圆度的主要因素可归纳为:（1）磨削过程中工艺系统的精度，这取决于磨床的精度和夹具的定位原理、结构参数及精度：（2）工艺系统的动态特性，特别是砂轮的平衡状态：在修整和磨削过程中，砂轮不平衡，会引起强迫振动，由于修整砂轮时修整器和砂轮的相对位置与磨削时工件和砂轮的相对位置的差异，以及这两种不同情况下工艺系统刚度的不同使振动造成磨削的工件表面不圆：而磨削工艺参数会影响砂轮与工件的接触刚度、砂轮的磨损过程以及工艺系统振动的阻尼特性，从而影响磨削过程的振动特性，最终反映在磨削工件的几何形状（即圆度）和其他表面质量上：（3）工艺系统的弹性变形引起的工件原始误差复映：分析磨削循环工件几何形状误差变化规律可知，磨削后的工件几何形状误差复映主要取决于工艺系统刚度、磨削工艺参数、砂轮磨损速度以及工件原始误差?。综合考虑以上影响因素，当工艺系统刚度和动态特性一定时（尤其砂轮必须经过良好的平衡），磨削工件的圆度误差主要取决于磨削工艺参数，因此可通过对磨削圆度试验数据的逐步回归建模，从多工艺参数中挑选重要参数，逐步引入回归方程，从而建立磨削圆度与工艺参数关系的最优回归方程。
3 试验方案
试验因素水平表
9 J5 I. |+ S% q4 a" a水平
因素
  h3 K' H/ e; cVw
/ h0 k5 k: g2 L: R% H( U. {' w' `（m/min）
a1
3 }1 x, f2 C- h+ h6 y" T（?m/r）
a2
（?m/r）
Sd
（mm/r）
td
（mm）
6 c# e  H9 ~5 x; e* yT
（mm3/mm）
/ D/ d9 t# V$ _$ W8 L3 NL2
  W% ]% N% ?& Y9 G5 Z, ]（mm）
1
9 b; S; _' ~8 t1 |6 c: f30
- y3 _) g$ ^+ ^) n7 a4 b4
' q$ G" J6 H( I4 {- S+ c- A0.5
0.1
0.01
1T0
0.025
; I& }' F# {+ ~( [& b* f2 I2
) c! }; d. t: Y# N50.5
8
1
0.2
0.02
8 B: ?+ w% M1 c2 Q* V% t5T0
: @* U; a4 x. r8 h3
( V& L) T8 J: O' R5 {83
12
1.8
0.3
0.03
& Y: [" S/ C+ O5 J0 \9 q10T0
0.04
注:表中T0为单件单位磨削宽度磨除工件体积，本试验T0=7.07mm3/mm。
$ a$ L  E$ ]1 l- l' a根据前面对磨削工件圆度影响因素的分析，兼顾工厂生产条件和设备性能的限制，把砂轮和磨削液的性能、砂轮线速度Vs、机床刚度和动态特性以及光磨时间作为不变因素，在磨削试验过程中尽量保持稳定。本研究选择的试验因素:修整砂轮的导程Sd，每次行程修整砂轮深度td，工件线速度Vw，粗进给工件每转磨削深度a1，细进给工件每转磨削深度a2，细进给行程L2，单位磨削宽度磨除金属体积T。各工艺参数除细进给行程L2取2水平外，其余均取， 水平，考虑工艺参数间交互效应，按混合型正交表进行试验。同时为了反映不同工件直径引起砂轮等效直径变化对磨削圆度的影响，以分别对208、308 和306 轴承进行磨削试验，并测量磨削工件的圆度。试验因素水平见右上表。
" X4 L- Z) k  b! }/ s6 p( c$ N试验条件:（1）机床为3MZ1310:全自动高速轴承内圈沟道磨床：（2）工件定位方式为双圆弧动支承：（3）砂轮为GB100ZR2A，直径ds=560mm，转速ns=1600r/min：（4）修整工具为单颗金刚石修整器，光修一次：（5）光磨时间2.5s，普通乳化液冷却液：（6）圆度测量仪Taylor-HOBSON：（7）磨削工件为208、306、308轴承内圈外沟道，直径dw=48.1、40、51mm，材料GCr15，硬度60～65HRC。
4 磨削圆度与工艺参数关系的建模
# E7 S. X5 Z0 m( I$ u数学模型假设
磨削内圈沟道圆度的逐步回归建模
" h% i- n0 g6 R利用磨削试验后实测内圈沟道的圆度数据，对上述式（3）线性化的响应函数式（4）进行逐步回归，即引入重要因素，剔除次要因素（注:F检验临界值取Fa=0.4），直到既不能剔除，又无法再引进变量的情况下逐步回归计算结束。定出响应函数式（4）各自变量系数的值，再将线性化后的响应函数代换复原，便可得到磨削工艺参数与圆度关系的数学模型为
R0=0.0623Vw0.442a13.262L20.201De1.08Tb1Sdb2a2b3tdb4
) }' Q) y. s! Q: F! X6 m; X& |
2 r( Z5 ^+ c: R4 L0 Ab1=-0.335lna1-0.368lnSd
" s5 a1 Z1 a1 T) j% k6 tb2=4.833+1.161lna1+2.178lnSd
9 y' V5 f1 W: Q- N! l$ Zb3=-0.233lna1
b4=-0.441lnSd
标准离差s=0.36，相关系数g=0.94，F=26>Fa=0.4
; \; G( s; |+ ^: N9 X. a! L: D从所建立模型的相关系数和方差分析F检验值可知，拟合效果比较满意，模型是可行的。
5 W5 y& {7 j$ h2 s根据前面对磨削圆度影响因素分析，考虑工艺参数间可能存在的交互效应，可假设磨削圆度的数学模型是
R0= KVwa1a1a2L2a3D4a4Twb1Sdb2ab23tdb4
/ v6 v1 R: ?. }
/ i4 D, S* \* t8 D3 P0 a5 {5 {b1=p1+lnVwb1a1b2Tb3Sdb4a2b5Deb6
' I( Y# A" V* G+ y, B* [3 l: Xb2=p2+c1lna1+c2lnSd
6 I" X4 L3 K) Ab3=p3+dlna1
b4=p4+f1lnSd+f2lntd
8 ~2 `9 L* v' Y5 h0 k. [式中R0代表圆度值，其余为工艺参数和待定常数，砂轮等效直径De=dsdw/（d3+dw）。
通过对式（3）两边取对数并进行变量代换后，可线性化为
y=B0+B1x1+B2x2+……+B19x19
（4）
# J9 A* j0 k( D: _1 ]/ z5 结果与讨论

5 R2 Q& ~( S- L  z. r  g（注:a1=8?m/r，a2=1.5?m/r，Vw=50.5m/min，L2=0.03mm，De=45mm，T=58mm3/mm）
图2 R0随Sd、td的变化情况
% U9 m% O4 @5 A2 \/ S2 |  d0 U
$ R! h$ S; z! l8 @. }4 u/ T1 v（注:Sd=0.15mm/r，td=0.02mm，Vw=50.5m/min，L2=0.03mm，De=45mm，T=58mm3/mm）
图3 R0随a1、a2的变化情况
  G$ Z5 u' I, j' E. j
7 v: S8 w+ A! D$ P! Z6 c* C% d! w7 m4 Q（注:Sd=0.15mm/r，td=0.02mm，a1=8?m/r，a2=1.5?m/r，De=45mm，L2=0.03mm）
图4 R0随Vw、T的变化情况
1 q# F! V0 S2 W( G
8 }  }$ p: k0 V8 b1 ?2 h% T（注:Sd=0.15mm/r，td=0.02mm，a1=8?m/r，a2=1.5?m/r，Vw=50.5m/min，T=58mm3/mm）
$ ^6 \6 V" l/ b0 M) o图5 R0随De、L2的变化情况