轴承沟道磨削表面微观形貌逐步回归建模
<FONT face=宋体>轴承沟道是滚动轴承的关键工作表面，其微观形貌直接影响整个滚动轴承的使用性能。磨削通常被作为轴承沟道的精密或半精密加工工序，研究其磨削工艺参数与加工表面微观形貌的关系十分必要。本文根据轴承沟道磨削表面形成机理，确定磨削表面微观形貌建模的试验因素，同时考虑因素间的交互效应。为使试验次数尽可能少，并且能获得可靠的有代表性的试验结果，结合混合型正交表设计多因素试验方案。采用逐步回归法，从多因素中挑选重要变量，逐步引入回归方程，从而建立试验数据的最佳回归方程。<BR><BR><STRONG>1 试验方案拟定</STRONG></FONT><P align=left><FONT face=宋体><STRONG>1.1 试验因素安排<BR> </STRONG>磨削是砂轮表面大量磨刃接连不断作用的结果。砂轮的磨粒大小，磨刃间距是随机分布，再加上磨削过程中存在磨粒严重的耕犁现象，支承砂轮磨粒的粘接剂弹性变形以及磨粒分布不均等现象，决定了磨削过程的复杂性<SUP>［1］</SUP>。尽管如此，仍可通过分析磨削表面粗糙度形成的几何关系和影响砂轮形貌的因素，同时考虑轴承沟道磨削是变切入进给速度，有光磨阶段磨削循环的特点，确定影响轴承内圈外沟道磨削表面粗糙度的因素，兼顾实际生产条件的限制以及简化试验过程，不可能把所有的影响因素都作为试验因素。把砂轮和磨削液的性能，光磨时间作为不变因素，在试验过程中尽量保持稳定。此外，通常砂轮均经过良好的平衡。可以认为它对磨削表面粗糙度的影响是不变的因素。本研究选择的试验因素：修整砂轮的导程和深度S<SUB>d</SUB>，t<SUB>d</SUB>，砂轮和工件的线速度V<SUB>s</SUB>，V<SUB>w</SUB>，粗进给和细进给时工件每转磨削深度a<SUB>1</SUB>，a<SUB>2</SUB>，细进给行程L<SUB>2</SUB>，单位磨削宽度磨除金属体积T，等效砂轮直径D<SUB>e</SUB>。结合混合型正交表L<SUB>54</SUB>（2<SUP>1</SUP>×3<SUP>25</SUP>)安排90次不同工艺参数组合的试验。测量的表面微观形貌指标包括表面粗糙度轮廓算术平均偏差R<SUB>a</SUB>，微观不平度平均间距S<SUB>m</SUB>。</FONT></P><P align=left><FONT face=宋体><STRONG>1.2 试验设备及试验条件<BR> </STRONG>(1)机床：3MZ1310A全自动高速轴承内圈沟道磨床；<BR> (2)工件定位方式：双圆弧动支承；<BR> (3)砂轮规格：GB100ZR<SUB>2</SUB>A，直径d<SUB>s</SUB>=512～600 mm，转速n<SUB>s</SUB>=1 600 r/min；<BR> (4)磨削对象：208，306，308轴承内圈沟道，直径d<SUB>w</SUB>=48.1 mm，40 mm，51 mm；材料GCr15,硬度HRC60～65。<BR><BR><STRONG>2 磨削表面微观形貌特征参数模型</STRONG></FONT></P><P align=left><FONT face=宋体><STRONG>2.1 模型的假设<BR> </STRONG>为了便于分析对比，本文假设了两种形式的数学模型，第一种是幂函数形式，不考虑各因素间的交互效应；第二种形式是根据机理分析，把可能存在交互效应的项引入数学模型，即考虑因素间可能存在的交互效应。<BR>2.1.1 不考虑因素间的交互效应</FONT></P><P align=left><FONT face=宋体>R=K<SUB>r</SUB> V<SUP>α<SUB>1</SUP>w</SUB> α<SUP>α<SUB>2</SUP>1</SUB> L<SUP>α<SUB>3</SUP>2</SUB> D<SUP>α<SUB>4</SUP>e</SUB> V<SUP>α<SUB>5</SUP>s</SUB> T<SUP>β<SUB>1</SUB></SUP> S<SUP>β<SUB>2</SUP>d</SUB> α<SUP>β<SUB>3</SUP>2</SUB> t<SUP>β<SUB>4</SUP>d </SUB>(1)</FONT></P><P align=left><FONT face=宋体>式中，R代表表面微观形貌特征参数R<SUB>a</SUB>或S<SUB>m</SUB>；K<SUB>r</SUB>,α<SUB>1</SUB>，α<SUB>2</SUB>，α<SUB>3</SUB>，α<SUB>4</SUB>，α<SUB>5</SUB>，β<SUB>1</SUB>，β<SUB>2</SUB>，β<SUB>3</SUB>，β<SUB>4</SUB>为待定常数。对式(1)两边取对数并进行变量代换后可线性化为</FONT></P><P align=left><FONT face=宋体>Y=B<SUB>0</SUB>+B<SUB>1</SUB>*X<SUB>1</SUB>+B<SUB>2</SUB>*X<SUB>2</SUB>+…+B<SUB>9</SUB>*X<SUB>9 </SUB>(2)</FONT></P><P align=left><FONT face=宋体>2.1.2 考虑因素间交互效应</FONT></P><P align=left><FONT face=宋体>R=K<SUB>r</SUB> V<SUP>α<SUB>1</SUP>w</SUB> α<SUP>α<SUB>2</SUP>1</SUB> L<SUP>α<SUB>3</SUP>2</SUB> D<SUP>α<SUB>4</SUP>e</SUB> V<SUP>α<SUB>5</SUP>s</SUB> T<SUP>β<SUB>1</SUB></SUP> S<SUP>β<SUB>2</SUP>d</SUB> α<SUP>β<SUB>3</SUP>2</SUB> t<SUP>β<SUB>4</SUP>d </SUB>(3)</FONT></P><P align=left><FONT face=宋体>式中，K<SUB>r</SUB>，α<SUB>1</SUB>，α<SUB>2</SUB>，α<SUB>3</SUB>，α<SUB>4</SUB>，α<SUB>5</SUB>为常数，β<SUB>i</SUB>(i=1,2,…，4)为</FONT></P><P><IMG src="http://www.chmcw.com/upload_files/article/20/1_bhafop2008319152955.gif"></P><P align=left><FONT face=宋体>同样通过对式(3)两边取对数并进行变量代换后可线性化为<BR>Y=B<SUB>0</SUB>+B<SUB>1</SUB>X<SUB>1</SUB>+B<SUB>2</SUB>X<SUB>2</SUB>+…+B<SUB>21</SUB>X<SUB>?</SUB>21 (4)</FONT></P><P align=left><FONT face=宋体><STRONG>2.2 磨削表面微观形貌特征参数R<SUB>a</SUB>、S<SUB>m</SUB>逐步回归建模<BR> </STRONG>利用式(1)、式(3)各线性化的响应函数，分别对实测的磨削表面微观形貌轮廓算术平均偏差R<SUB>a</SUB>和微观不平度间距S<SUB>m</SUB>试验数据进行逐步回归，即引入重要因素，剔除不重要的因素(F检验临界值取F<SUB>α</SUB>=0.4)。直至既不能剔除，又无法再引进变量，最后定出式(2)、式(4)自变量各系数的值，再将其代换复原，可得R<SUB>a</SUB>，S<SUB>m</SUB>的数学模型。<BR>2.2.1 不考虑交互效应的模型</FONT></P><P align=left><FONT face=宋体>R<SUB>a</SUB>=0.928×10<SUP>-2</SUP>V<SUP>0.282</SUP><SUB>w</SUB>α<SUP>0.088</SUP><SUB>1</SUB>T<SUP>0.039</SUP>S<SUP>0.763</SUP><SUB>d</SUB>α<SUP>0.081</SUP><SUB>2</SUB>t<SUP>0.105</SUP><SUB>d</SUB>D<SUP>1.761</SUP><SUB>e</SUB>V<SUP>-0.620</SUP><SUB>s </SUB>(5)</FONT></P><P align=left><FONT face=宋体>标准离差σ<SUB>a</SUB>=0.14，相关系数γ=0.84；<BR>S<SUB>m</SUB>=4.376×10<SUP>3</SUP>V<SUP>-0.060</SUP><SUB>w</SUB>a<SUP>-0.125</SUP><SUB>1</SUB>L<SUP>-0.163</SUP><SUB>2</SUB>S<SUP>0.586</SUP><SUB>d</SUB>a<SUP>-0.131</SUP><SUB>2</SUB>t<SUP>0.123</SUP><SUB>d</SUB>V<SUP>-0.783</SUP><SUB>s </SUB>(6)</FONT></P><P align=left><FONT face=宋体>标准离差σ<SUB>m</SUB>=12.7,相关系数γ=0.83。<BR>2.2.2 考虑因素间交互效应的模型</FONT></P><P align=left><FONT face=宋体>R<SUB>a</SUB>=5.89×10<SUP>-2</SUP>V<SUP>0.317</SUP><SUB>w</SUB>D<SUP>1.609</SUP><SUB>e</SUB>V<SUP>-0.890</SUP><SUB>s</SUB>S<SUP>β<SUB>2</SUP>d</SUB> α<SUP>β<SUB>3</SUP>2</SUB> t<SUP>β<SUB>4</SUP>d </SUB>(7)</FONT></P><P><IMG src="http://www.chmcw.com/upload_files/article/20/1_r9uxqq2008319153017.gif"></P><P align=left><FONT face=宋体>标准离差σ<SUB>a</SUB>=0.11,相关系数γ=0.90；</FONT></P><P align=left><FONT face=宋体>S<SUB>m</SUB>=4.268×10<SUP>5</SUP>a<SUP>0.150</SUP><SUB>1</SUB>L<SUP>-0.176</SUP><SUB>2</SUB>V<SUP>-0.938</SUP><SUB>s</SUB>T<SUP>β<SUB>1</SUB></SUP>S<SUP>β<SUB>2</SUP>d</SUB>t<SUP>β<SUB>4</SUP>d </SUB>(8)</FONT></P><P><IMG src="http://www.chmcw.com/upload_files/article/20/1_pjan372008319153033.gif"></P><P align=left><FONT face=宋体>标准离差σ<SUB>m</SUB>=9.98，相关系数γ=0.90<BR> 从以上各模型的相关系数γ和方差分析(本文略)可知，拟合效果比较满意；在生产现场三种不同工艺参数条件下，磨削6件208、306、308轴承内圈沟道，其R<SUB>a</SUB>、S<SUB>m</SUB>的实测值与式(7)、式(8)的计算值偏差分别小于1.96σ<SUB>a</SUB>，1.96σ<SUB>m</SUB>(即可信度97.5%)，说明所拟合的模型是可行的。在R<SUB>a</SUB>，S<SUB>m</SUB>模型的两种形式中，第一种形式比较简单、直观，容易看出各工艺参数对表面微观形貌的影响情况，但回归效果不如第二种形式的模型；第二种形式比较复杂，但回归效果较好，能较为精确地拟合出试验数据的变化规律，作为磨削机理分析以及工艺参数优化的约束函数较为合适，由式(7)，式(8)可作出图1～图4，表示磨削工艺参数对表面微观形貌特征参数R<SUB>a</SUB>，S<SUB>m</SUB>的影响情况。<BR><BR></FONT><FONT face=宋体><STRONG>3 结果讨论</STRONG></FONT></P><P align=left><FONT face=宋体><STRONG>3.1 修整砂轮的导程和深度对磨削表面微观形貌的影响</STRONG></FONT></P><P align=center><IMG src="http://www.chmcw.com/upload_files/article/20/1_dwkcyg2008319153055.gif"></P><P align=center><FONT face=宋体><STRONG>图1 磨削表面微观形貌特征参数R<SUB>a</SUB>,S<SUB>m</SUB>随S<SUB>d</SUB>，t<SUB>d</SUB>的变化情况<BR></STRONG>a<SUB>1</SUB>=1μm/r，a<SUB>2</SUB>=1.5μm/r，<BR>L<SUB>2</SUB>=0.03mm，T=58mm<SUP>3</SUP>/mm,<BR>v<SUB>w</SUB>=50m/min，v<SUB>s</SUB>=47m/s，<BR>D<SUB>e</SUB>=44.21mm</FONT></P><P align=left><FONT face=宋体> 图1a,lb分别表示修整砂轮的导程S<SUB>d</SUB>和深度t<SUB>d</SUB>对磨削表面微观形貌特征参数R<SUB>a</SUB>，S<SUB>m</SUB>的影响情况。分析图1a可知，当t<SUB>d</SUB>≤0.02 mm，S<SUB>d</SUB>的取值范围S<SUB>d</SUB><1.1 mm/r，R<SUB>a</SUB>随S<SUB>d</SUB>的减小而增大，这是因为过于精细修整，砂轮表面磨刃数显著增加，以致容屑空间太小，相对于磨削深度而言，实际磨刃变钝，耕犁现象和磨刃的破裂加剧；当S<SUB>d</SUB>增大到一定值（S<SUB>d</SUB>≥1.1 mm/r)或t<SUB>d</SUB>>0.02 mm，R<SUB>a</SUB>随S<SUB>d</SUB>的增大而增大。其主要原因是修整导程小，修整的砂轮其表面有效磨刃数多，磨刃的等高性好，因此磨削的工件表面微观形貌R<SUB>a</SUB>小，反之，R<SUB>a</SUB>则增大。而R<SUB>a</SUB>随t<SUB>d</SUB>的变化情况是：S<SUB>d</SUB>≥1.1 mm/r时，t<SUB>d</SUB>的增大，R<SUB>a</SUB>随之增大。这说明增大修整深度t<SUB>d</SUB>，修整时砂轮表面破碎的微粒增大，修整的砂轮表面更加不平。同样从图1b可知，S<SUB>m</SUB>随修整导程S<SUB>d</SUB>和深t<SUB>d</SUB>的变化规律与R<SUB>a</SUB>类同。</FONT></P><P align=left><FONT face=宋体><STRONG>3.2 工件和砂轮线速度对磨削表面微观形貌的影响</STRONG></FONT></P><P align=center><IMG src="http://www.chmcw.com/upload_files/article/20/1_65n66a2008319153120.gif"></P><P align=center><FONT face=宋体><STRONG>图2 磨削表面微观形貌特征参数R<SUB>a</SUB>,S<SUB>m</SUB>随v<SUB>w</SUB>，v<SUB>s</SUB>的变化情况<BR></STRONG>a<SUB>1</SUB>=8μm/r，a<SUB>2</SUB>=1.5μm/r，<BR>L<SUB>2</SUB>=0.03mm，T=58mm<SUP>3</SUP>/mm,<BR>D<SUB>e</SUB>=44.21mm，t<SUB>d</SUB>=0.02mm，<BR>s<SUB>d</SUB>=0.15mm/r</FONT></P><P align=left><FONT face=宋体>图2a,2b分别表示磨削表面微观形貌特征参数R<SUB>a</SUB>，S<SUB>m</SUB>随工件和砂轮线速度v<SUB>w</SUB>、v<SUB>s</SUB>的变化情况。分析图2a可知，磨削表面微观形貌的轮廓算术平均偏差R<SUB>a</SUB>随工件速度v<SUB>w</SUB>的增大而增大。这是由于随着工件速度v<SUB>w</SUB>的提高，单颗磨粒切深度变大，磨刃的切痕相应增大。这虽使砂轮表面参加磨削的有效磨刃数增加，但同时也使砂轮表面参加工作磨刃的不平度增加。此外，单颗磨刃切深增大所增加的砂轮表面工作磨刃切深较小，在磨削形成工件表面过程中，其作用相应形成工件表面微观不平的峰顶，这些磨刃耕犁引起的隆起使磨削的工件表面微观不平的峰顶增高，从而使R<SUB>a</SUB>增大；另一方面，工件速度增大，磨粒切深增大，砂轮的磨损形式更趋于磨粒局部破碎和脱落，也使R<SUB>a</SUB>增大。R<SUB>a</SUB>随砂轮线速度v<SUB>s</SUB>的提高而减小；v<SUB>s</SUB>提高，单颗磨刃的切深减小，相应磨刃的切痕减小；磨刃的切痕减小，虽使砂轮表面有效磨刃数有所减少，但v<SUB>s</SUB>提高，单位时间内通过磨削区磨刃数增加，每个磨刃对工件表面的作用次数增加，因此实际工作磨刃的等高性好。同时，v<SUB>s</SUB>提高使形成工件表面的包络线变多，最后导致R<SUB>a</SUB>减小。从图2b可以看出，工件速度v<SUB>w</SUB>对磨削表面微观不平度平均间距S<SUB>m</SUB>的影响不明显。原因在于：一方面，v<SUB>w</SUB>提高，磨刃的切深增大，砂轮表面有效磨刃数增多，产生的工件表面不平度平均间距s<SUB>m</SUB>减小；另一方面，v<SUB>w</SUB>提高，砂轮的磨损更趋于破碎或脱落，使不平度平均间距S<SUB>m</SUB>增大。对于本研究的砂轮和工件性能，两者的作用效果相互抵消，综合作用结果表现为v<SUB>w</SUB>和S<SUB>m</SUB>影响不明显。砂轮线速度v<SUB>s</SUB>对S<SUB>m</SUB>的影响情况与R<SUB>a</SUB>类同。</FONT></P><P align=left><FONT face=宋体><STRONG>3.3 粗进给和细进给磨削深度对磨削表面微观形貌的影响</STRONG></FONT></P><P align=center><IMG src="http://www.chmcw.com/upload_files/article/20/1_wozmhc2008319153142.gif"></P><P align=center><FONT face=宋体><STRONG>图3 磨削表面微观形貌特征参数R<SUB>a</SUB>,S<SUB>m</SUB>随a<SUB>1</SUB>，a<SUB>2</SUB>的变化情况<BR></STRONG>v<SUB>w</SUB>=50m/min，v<SUB>s</SUB>=47m/s，<BR>L<SUB>2</SUB>=0.03mm，T=58mm<SUP>3</SUP>/mm,<BR>D<SUB>e</SUB>=44.21mm，t<SUB>d</SUB>=0.02mm，<BR>s<SUB>d</SUB>=0.15mm/r</FONT></P><P align=center><IMG src="http://www.chmcw.com/upload_files/article/20/1_ezyhif200831915324.gif"></P><P align=center><FONT face=宋体><STRONG>图4 磨削表面微观形貌特征参数R<SUB>a</SUB>,S<SUB>m</SUB>随D<SUB>e</SUB>，L<SUB>2</SUB>的变化情况<BR></STRONG>v<SUB>w</SUB>=50m/min，v<SUB>s</SUB>=47m/s，<BR>a<SUB>1</SUB>=8μm/r，a<SUB>2</SUB>=1.5μm/r，<BR>T=58mm<SUP>3</SUP>/mm,t<SUB>d</SUB>=0.02mm，<BR>s<SUB>d</SUB>=0.15mm/r</FONT></P><P align=left><FONT face=宋体>分析图3a可知，当粗进给工件每转磨削深度a<SUB>1</SUB>≥5 μm/r，磨削表面微观形貌特征参数R<SUB>a</SUB>随细进给磨削深度a<SUB>2</SUB>的增大而增大。主要原因：a<SUB>2</SUB>的增大使单颗磨粒切深增大，从而使R<SUB>a</SUB>增大；当a<SUB>1</SUB><5 μm/r时，R<SUB>a</SUB>随a<SUB>2</SUB>的增大反而小幅度减小。R<SUB>a</SUB>随a<SUB>1</SUB>的变化趋势与之类似：当a<SUB>2</SUB><0.8 μm/r，R<SUB>a</SUB>随a<SUB>1</SUB>增大而减小；当a<SUB>2</SUB>≥0.8 μm/r，R<SUB>a</SUB>随a<SUB>1</SUB>的增大而增大，这些情况主要与粗进给和细进给磨削阶段砂轮表面磨刃的磨损形式和自励性有关。分析图3b可知，磨削表面微观形貌特征参数S<SUB>m</SUB>随a<SUB>1</SUB>、a<SUB>2</SUB>的增大而减小。这是由于同时参加磨削的砂轮表面磨刃数随a<SUB>1</SUB>、a<SUB>2</SUB>的增大而增多。其中粗进给工件每转磨削深度a<SUB>1</SUB>通过影响砂轮表面磨刃磨损形式对s<SUB>m</SUB>起作用，a<SUB>1</SUB>的增大，单颗磨刃的受力增大，磨粒更容易局部破碎和磨损，部分磨粒局部破碎后，单颗磨粒上出现了多个磨刃，这就使磨削后工件表面微观不平度平均间距S<SUB>m</SUB>减小。</FONT></P><P align=left><FONT face=宋体><STRONG>3.4 砂轮等效直径和细进给行程对磨削表面微观形貌的影响<BR> </STRONG>图4a、图4b分别表示磨削表面微观形貌特征参数R<SUB>a</SUB>、S<SUB>m</SUB>随砂轮等效直径D<SUB>e</SUB>和细进给行程L<SUB>2</SUB>的变化情况。分析图4a可知，微观形貌特征参数R<SUB>a</SUB>随D<SUB>e</SUB>的增大而增大。主要原因是：随D<SUB>e</SUB>的增大，单颗磨刃的切除体积增加，磨刃耕犁引起的隆起高度增加；细进给行程L<SUB>2</SUB>对R<SUB>a</SUB>无明显影响。分析图4b可知，微观形貌特征参数S<SUB>m</SUB>随D<SUB>e</SUB>和L<SUB>2</SUB>的增大而减小。</FONT><STRONG></P></STRONG>
页:
[1]