砂带磨削静态接触理论建模与有限元仿真
李建勇, 程继发, 樊文刚, 吕洪宾
北京交通大学 机械与电子控制工程学院,北京 100044

第一作者:李建勇(1962—),男,北京市人,教授,博士,博士生导师.研究方向为钢轨打磨技术与装备. email:jyli@bjtu.edu.cn.

摘要

针对砂带表面磨粒与工件微观静态接触行为问题,提取了砂带表面形貌特征并重建,基于赫兹接触理论建立了包含打磨压力、压入深度、磨粒刃端曲率半径、攻角等参量的静态接触理论模型,研究了磨粒与工件表面接触区域和接触应力分布规律.通过有限元仿真对不同打磨压力下的接触应力、接触面积和压入深度进行了分析和模型验证.结果表明:在一定误差范围内仿真值与理论值吻合,验证了理论模型的正确性和有效性;磨粒与工件表面接触应力、接触面积和压入深度随打磨压力的增大而增大并于250 N时趋于稳定,同样接触磨粒数目逐渐增多并趋于稳定.

关键词: 砂带磨削; 磨粒; 微观接触; 有限元仿真
中图分类号:TG74.2 文献标志码:A 文章编号:1673-0291(2018)06-0125-06
Theory modeling and finite element simulation of belt grinding static contact
LI Jianyong, CHENG Jifa, FAN Wengang, LYU Hongbin
School of Mechanical, Electronic and Control Engineering, Beijing Jiaotong University, Beijing 100044,China
Abstract

The problem of the micro-static contact behavior between the abrasive grains and the workpiece, the surface features of the abrasive belt are collected and reconstructed. Based on the Hertzian contact theory, a theoretical model involving parameters of grinding pressure, indentation depth, grains’ lip-end curvature radius and attack angle is established to illustrate the static contact behavior. The distribution rule of the contact area and the contact stress between the abrasive grains and the workpiece surface is obtained. Then, the finite element simulation is conducted to validate the theoretical model through analyzing contact stress, contact area and indentation depth under different grinding pressure. The results show that the simulated values are in good agreement with the theoretical values within a limited error range, which prove the correctness and the validity of the theoretical model. It can be seen that the contact stress, the contact area and the indentation depth increase along with the growth of grinding pressure and become stable when it reaches 250 N. Meanwhile, the number of contact abrasive grains also gradually increases and stabilizes.

Keyword: belt grinding; grain; microscopic contact; finite element simulation

砂带磨削是砂带表面有效接触磨粒对工件表面进行大量重复微量切削的过程, 具有高效磨削、冷态磨削、弹性磨削等优点[1], 砂带表面磨粒形状、磨粒间距、磨粒密度和磨粒出刃高度将直接影响磨粒与工件表面的接触区域和接触应力的分布、有效接触磨粒数以及磨粒切削受力, 进而对砂带磨削力、材料去除率、工件表面质量和砂带磨损产生影响[2].

目前, 国内外对基于砂带真实形貌特征的研究相对较少.王文玺等[3]通过获得砂带表面磨粒密度、磨粒出刃高度及顶锥角等评价参数, 从宏观整体到局部磨粒对表面形貌特征完成了定性和定量评价.张叠等[4]分析了ZrO2工程陶瓷进行砂带磨削加工过程中, 砂带粒度和磨削用量的不同对磨削加工效率、工件表面质量的影响规律.黄云等[5]基于单位面积内每颗磨粒所受正压力相等的假设, 推导了磨粒最大划痕深度与材料硬度、所受正压力、锥顶角大小和有效磨粒粒度的关系.刘伟等[6]选用截角八面体模拟了砂带表面金刚石磨粒, 并且基于Johnson-Holmquist ceramic硬脆材料本构模型, 采用有限元网格法进行单颗磨粒直线磨削仿真.Mezghani等[7]运用分水岭切割算法对砂带表面磨粒信息进行提取, 获得磨粒密度、平均出刃高度和磨粒平均顶锥角3个评价指标, 并以此描述砂带表面形貌特征.Zhao等[8]以单个球形模型为基础推导了3个材料变形阶段对应的接触面积和压力.Zhang等[9]提出一种基于神经网络的砂带磨削接触力分布的计算方法.Khellouki等[10]研究了砂带磨损和磨削工艺参数对工件表面粗糙度的影响规律, 指出平均接触压力和有效接触磨粒数为工艺参数中的关键因子.

上述文献主要是对砂带表面磨粒进行宏观假设描述和接触分析, 未对真实参数特征的砂带表面磨粒微观接触进行分析.本文作者从真实参数特征磨粒的角度出发, 对砂带磨削静态接触理论建模深入研究, 首先对砂带表面形貌特征进行提取并重建, 建立磨粒与工件表面静态接触理论模型, 研究接触区域、接触应力分布规律; 进而对不同打磨压力下的接触应力、接触区域和压入深度进行有限元仿真, 仿真值和理论计算对比分析验证理论模型的正确性和有效性.

1 砂带表面形貌特征提取与重建

砂带采用“ 静电植砂” 工艺, 将破碎机碾压后制得的针状磨粒大体均匀地植在砂带表面, 因植砂过程中磨粒发生尖端放电效应而使绝大部分磨粒刃尖朝上, 从而使外露的磨粒更接近于球顶圆锥模型, 对于单颗磨粒 i, 由刃端曲率半径Ri和攻角αi2个参数描述的简化模型如图1所示.

图1 单颗磨粒形貌示意图Fig.1 Single abrasive morphology

磨粒顶端球形压力工作表面的最大深度为

h=Ri(1-cosαi)(1)

通过基于白光共聚焦原理的三维表面形貌轮廓仪可以测得砂带表面磨粒分布三维点云数据, 三维表面形貌轮廓仪如图2所示, 并通过Matlab软件将三维点云数据生成三维曲面如图3所示.

图2 三维表面形貌轮廓仪Fig.2 Three-dimensional surface topography meter

三维曲面示意图

Fig 3 Three-dimensional surface diagram

由图3可知, 三维坐标点(x, y, z)处的磨粒高度为z , 通过比较相邻八点高度得高度差[1]

Δz(k, l)=zi(k, l)-z(s, r)(2)

式中: s和r表示磨粒i周围的8个磨粒处的横、纵坐标值, z(s, r)为磨粒i周围的8个磨粒的高度; s取k-1, k, k+1; r取l-1, l, l+1; k和l为磨粒i处的横、纵坐标.

如果所有的 Δz(k, l)> 0, zi(k, l)被标记为砂带表面磨粒顶点(“ +” ), 通过Matlab标记砂带表面磨粒顶点如图4所示.

图4 磨粒顶点分布示意图Fig.4 Abrasive grains vertex distribution diagram

对于磨粒球顶表面, 磨粒i分别沿坐标系x、y方向的刃端曲率半径分别为[1]

Rix(k, l)=(180(Δx)2)/(2zk-3, l-27zk-2, l+270zk-1, l-490zk, l+270zk+1, l-27zk+2, l+zk+3, l)(3)

Riy(k, l)=(180(Δy)2)/2zk, l-3-27zk, l-2+270zk, l-1-490zk, l+270zk, l+1-27zk, l+2+zk, l+3)(4)

式中: ΔxΔy分别为沿坐标系xy方向的采样间距, 磨粒 i的刃端曲率半径为

Ri=(Rix+Riy)1/2(5)

磨粒攻角沿 xy方向的分量分别为

tanαix=di/aix(6)

tanαiy=di/aiy(7)

式中: di为磨粒压入深度; aixaiy分别为磨粒沿x、y方向接触半径, 则攻角为

αi=(αix+αiy)/2(8)

由图4可知, 全局坐标系(x, y, z)下磨粒顶点分布, 通过空间坐标变换求得局部坐标系 xi, yi, zi

下磨粒分布模型, 局部坐标系如图5所示.图5中 x0, y0, z0是全局坐标系下磨粒底面中心点处的三维坐标; x1, y1, z1是全局坐标系下磨粒上任一点处的三维坐标; x2, y2, z2是局部坐标系下磨粒上任一点处的三维坐标点.

图5 磨粒三维坐标变换示意图Fig.5 Three-dimensional coordinate transformation of abrasive grains

通过Matlab软件对重建后的砂带表面磨粒三维点云数据进行处理, 生成系统能够识别的数据文件, 并进行砂带形貌曲面建模, 如图6所示.

图6 提取特征后砂带表面形貌Fig.6 Surface morphology of the abrasive belt after its feature extraction

由图6可知, 在保留磨粒高度、刃端曲率半径、攻角等主要特征参数后, 重建后的磨粒为球顶圆锥, 符合简化模型.基于Pro/Engineer三维建模软件, 将曲面转化为三维实体模型如图7所示.

图7 真实参数特征的砂带表面形貌Fig.7 Real parametric characteristics of the abrasive belt surface morphology

2 磨粒与工件接触理论模型

图8 砂带表面磨粒与工件横截面接触示意图Fig.8 Cross-sectional contact diagram of abrasive grains on the belt and workpiece

假设以砂带表面作为 z轴原点, 距离工件表面高度

h0, 施加打磨压力Fn 见图8, 所有高度z> h0的磨粒都与工件表面接触, 高度为zi的磨粒i压入工件表面深度为

di=zi-h0(9)

当磨粒压入工件表面深度较小时, 即满足条件: dih, 这时仅磨粒顶端球形与工件表面接触; 如果 di> h, 这时磨粒锥形表面与工件表面接触.

2.1 磨粒顶端球形接触模型

如果接触应力 PY(屈服极限), 接触表面发生弹性变形, 根据赫兹接触理论, 磨粒弹性接触作用力、面积分别为[11]

Fs=43E* Ri1/2d3/2(10)As=πai2=πdRi(11)

式中:d是磨粒弹性压入深度; ai为磨粒接触半径; E* 为等效弹性模量, 可得

1E* =1-v12E1+1-v12E2(12)

式中: E1、E2分别为磨粒和工件的弹性模量; v1、v2为对应泊松比.则弹性接触应力为

Ps=Fs/As=0.4E* d/Ri(13)

如果 Y< PH(材料硬度), 接触表面开始发生弹塑性接触[12], 磨粒弹塑性压入深度为

d1=(0.92HE* )2Ri(14)

式中: H为塑性接触应力.

弹塑性接触应力、面积分别为[13]:

Pe=H(1-0.6lnd2-lndlnd2-lnd1)(15)Ae=πRid1[1-2(d-d1d2-d1)3+3(d-d1d2-d1)2](16)

式中: d2为磨粒完全塑性压入深度.

磨粒弹塑性接触作用力为

Fe=AePe(17)

如果P> H, 接触表面发生完全塑性接触, 此时 , d2=54d1, 塑性接触面积、作用力分别为[14]:

Ap=πai2=2πRid2(18)

Fp=2πd2RiH(19)

2.2 磨粒圆锥表面接触模型

磨粒与工件表面发生弹性接触时, 弹性接触应力、面积为[14]

σc=0.2E* tanαi(20)

As=πd2cot2αi(21)

此时, 磨粒弹性接触作用力为

Fc=σcAc(22)

如果 Y< PcH, 接触表面开始发生弹塑性接触, 弹塑性压入深度 d1=d, 此时, 弹塑性接触面积为 Ac, 弹塑性接触应力为[14]

Pc=23(1+lnE* tanαi3Y)(23)

如果 Pc> H时, 接触表面发生完全塑性接触, 塑性压入深度 d2=54d1, 具体实施流程见图9.

图9 磨粒与工件接触理论计算流程图Fig.9 Flow chart for calculating the contact theory of abrasive grains and workpiece

3 有限元仿真验证

采用Abaqus有限元软件建立砂带与工件之间的三维仿真模型, 陶瓷磨粒与工件材料物理特性参数如表1所示, 失效准则参数 D1~D5分别取0.05, 3.44, -2.12, 0.002, 0.61, 工件材料本构模型参数中屈服极限A为792 Mpa、硬化模量 B为510 Mpa、硬化系数 n为0.26、应变率敏感指数 C0.014、热软化系数m为1.03 .

表1 砂带与工件材料特性参数 Tab.1 Parameters of abrasive belt and workpiece material characteristics
3.1 不同打磨压力下接触应力仿真分析

选取10个不同的打磨压力 Fn, 其值分别为50, 100, …, 500 N, 这样压力的分布范围广, 趋势观察更加明显, 然后分别利用理论模型和仿真模型计算接触区域内最大点接触应力, 选取100、300、500 N 3个打磨压力值查看有限元应力仿真结果, 如图10所示.

图10 不同打磨压力下的接触应力仿真结果Fig.10 Simulation results of contact stress under different grinding pressure

由图10知, 随着打磨压力的增大, 接触应力逐渐增大, 同时, 接触磨粒数目逐渐增多, 接触区域逐渐增大, 理论计算与仿真获得的接触应力曲线见图11.由图可知, 随着打磨压力的增大, 接触应力趋于稳定在800 MPa左右.在打磨压力 Fn≤ 250 N时, 随打磨压力Fn的增大, 接触压力P呈现增大趋势, 同时有效接触磨粒数目逐渐增多; 当Fn≥ 250 N时, 接触压力P逐渐趋于稳定. Fn在50~250 N范围内时, 理论与仿真获得的 P值出现较大误差为15%, 且误差值在允许的范围0%~20%内, 故仿真数据与理论计算数据基本吻合.

图11 不同打磨压力的接触应力理论曲线和仿真值Fig.11 Theoretical curves and simulation values of contact stress under different grinding pressures

3.2 不同打磨压力下压入深度仿真分析

通过仿真求得磨粒最大压入深度与打磨压力的关系, 理论计算与仿真获得的压入深度变化曲线见图12, 可以看出, 随着打磨压力的增大, 压入深度趋于稳定在70 μ m左右.当 Fn≤ 200 N时, 压入深度随打磨压力的增大逐渐增大, 同时接触磨粒数目逐渐增多; 当200 N≤ Fn≤ 500 N时, 随打磨压力的增大, 压入深度趋于稳定, 同时接触磨粒数目趋于稳定.当 Fn在50~200 N范围内时, 理论计算与仿真获得的压入深度值出现较大误差为13.3%, 且误差值在允许的范围0%~20%内, 仿真值和理论计算结果基本吻合.

图12 不同打磨压力的压入深度理论曲线与仿真值Fig.12 Theoretical curves and simulation values of indentation depth under different grinding pressures

3.3 不同打磨压力下接触面积仿真分析

依据理论模型计算求得砂带表面磨粒与工件接触面积, 并通过仿真求得相应打磨压力下磨粒与工件之间的接触面积, 理论计算与仿真获得的接触面积变化曲线如图13所示.可以看出, 随着打磨压力的增大, 接触面积趋于稳定在0.65 mm2左右, 当打磨压力 Fn≤ 400 N时, 接触面积随打磨压力的增大而增大且增加趋势较稳定; 当打磨压力400 N≤ Fn≤ 500 N时, 接触面积随打磨压力的增大趋于稳定, 同时接触磨粒数目也趋于稳定.当 Fn在200~350 N范围内时, 理论与仿真获得的接触面积出现较大误差为9.3%, 且误差值在允许的范围0%~20%内, 仿真值和理论计算结果基本吻合.

图13 不同打磨压力的接触面积理论曲线和仿真值Fig.13 Theoretical curves and simulation values of contact area under different grinding pressures

4 结论

1)对砂带表面特征进行提取并重建, 建立了砂带表面磨粒与工件表面接触理论模型, 实现了接触区域和接触应力分布的理论计算。获得了打磨压力、压入深度、磨粒刃端曲率半径、攻角等参量对接触面积和接触应力分布的影响规律.

2)通过Aaqus有限元仿真和理论模型计算求得不同打磨压力下磨粒与工件之间接触应力、压入深度与接触面积随打磨压力的变化趋势, 当打磨压力达到250 N时, 接触应力、压入深度、接触面积都趋于稳定, 接触应力稳定在800 Mpa, 压入深度稳定在70 μ m, 接触面积稳定在0.65 mm2.

参考文献
[1] JOURANI A, HAGEGE B, BOUVIER S, et al. Influence of abrasive grain geometry on friction coefficient and wear rate in belt finishing[J]. Tribology International, 2013, 59: 30-37. [本文引用:3]
[2] NGUYEN T A, BUTLER D L. Simulation of precision grinding process, part 1: generation of the grinding wheel surface[J]. International Journal of Machine Tools & Manufacture, 2005, 45(11): 1321-1328. [本文引用:1]
[3] 王文玺, 李建勇, 樊文刚, . 砂带三维表面形貌特征量化评价方法[J]. 华南理工大学学报 (自然科学版), 2016, 44(12): 14-22.
WANG Wenxi, LI Jianyong, FAN Wengang, et al. Quantitative evaluation method for 3D surface topography of abrasive belt[J]. Journal of South China University of Technology(Natural Science Edition), 2016, 44(12): 14-22. [本文引用:1]
4 (in Chinese) [本文引用:1]
[4] 张叠, 黄云. ZrO2工程陶瓷砂带磨削实验及工艺研究[J]. 机械科学与技术, 2015, 34(12): 1966-1970.
ZHANG Die, HUANG Yun. Experiment and technology of abrasive belt grinding for ZrO2 engineering ceramics[J]. Mechanical Science and Technology for Aerospace Engineering, 2015, 34(12): 1966-1970. (in Chinese) [本文引用:1]
[5] 黄云. 核电高压容器强力高效砂带磨削方法及应用研究[D]. 重庆: 重庆大学, 2009.
HUANG Yun. Research on heavy abrasive belt grinding and its application in machining of nuclear power high pressure vessel[D]. Chongqing: Chongqing University, 2009. (in Chinese) [本文引用:1]
[6] 刘伟, 邓朝辉, 万林林, . 单颗金刚石磨粒切削氮化硅陶瓷仿真与试验研究[J]. 机械工程学报, 2015, 51(21): 191-198.
LIU Wei, DENG Chaohui, Wan Linlin, et al. Simulation and experiment study for silicon nitride cutting with single diamond grain[J]. Journal of Mechanical Engineering, 2015, 51(21): 191-198. (in Chinese) [本文引用:1]
[7] MEZGHANI S, MANSORI M E. Abrasiveness properties assessment of coated abrasives for precision belt grinding[J]. Surface & Coatings Technology, 2008, 203(5): 786-789. [本文引用:1]
[8] ZHAO Y, MAIETTA D M, CHANG L. An asperity microcontact model incorporating the transition from elastic deformation to fully plastic flow[J]. Journal of Tribology, 2000, 122(1): 86-93. [本文引用:1]
[9] ZHANG X, KNEUPNER K, KUHLENK B. A new force distribution calculation model for high-quality production processes[J]. International Journal of Advanced Manufacturing Technology, 2006, 27(7): 726-732. [本文引用:1]
[10] KHELLOUKI A, RECH J, ZAHOUANI H. The effect of abrasive grain’s wear and contact conditions on surface texture in belt finishing[J]. Wear, 2007, 263(1/6): 81-87. [本文引用:1]
[11] JOHNSON K L, KENDALL K, ROBERTS A D. Surface energy and the contact of elastic solids[J]. Proceedings of the Royal Society Y, 1971, 324(1558): 301-313. [本文引用:1]
[12] QUARRELL A G. The hardness of metals[J]. Nature, 1952, 170(4333): 818. [本文引用:1]
[13] FRANCIS H A. Phenomenological analysis of plastic spherical indentation[J], Journal of Engineering Materials and Technology, 1976, 98(3): 272-281. [本文引用:1]
[14] BRISCOE B J. Contact mechanics[J]. Tribology International, 1985, 19(2): 109-110. [本文引用:3]