欢迎访问一起赢论文辅导网
机械论文
当前位置:首页 > 机械论文
弹塑性粗糙表面实际接触面积演变规律研究_刘检华
来源:一起赢论文网     日期:2021-04-28     浏览数:2336     【 字体:

  57卷第**2021 *月 机  械  工  程  学  报 JOURNAL OF MECHANICAL ENGINEERING Vol.57  No.** ***     2021 弹塑性粗糙表面实际接触面积演变规律研究* 刘检华    张飞凯  丁晓宇(北京理工大学机械与车辆学院    北京    100081) 摘要:粗糙表面的实际接触面积直接影响精密机电设备配合表面的摩擦因数、热导/电导率、接触应力等。然而,针对弹塑性接触行为进行系统性地研究很少,弹塑性粗糙表面实际接触面积的演变规律尚不明确。针对这一问题,通过理论研究获得了弹塑性粗糙表面实际接触面积的影响因素,并引入弹塑性接触力学数值计算方法,对具有不同材料参数和表面形貌特征的弹塑性粗糙表面的接触行为进行数值仿真计算,得到各种情况下的表面实际接触面积随着平均接触压力的变化曲线,总结出了实际接触面积演变规律与这些影响因素之间的映射关系,并最终得到实际接触面积的计算公式,为弹塑性接触力学的相关工程应用奠定了基础。 关键词:接触力学;弹塑性;粗糙表面;实际接触面积;接触数值计算 中图分类号:TB42 Study on the Evolution of Real Contact Area between Elastic-plastic Rough Surfaces LIU Jianhua   ZHANG Feikai  DING Xiaoyu(School of Mechanical Engineering, Beijing Institute of Technology, Beijing 100081) AbstractThe  real contact  area  of  elastic-plastic  rough  surfaces  significantly  affects the  friction  coefficient,  the  thermal/electrical conductivity,  and  the contact stress of the assembly surfaces of precision  electromechanical devices. However,  as far as the authors known,  the  elastic-plastic  contact behavior is rarely  studied  systematically,  and  the  evolution  of  the  real contact  area  has not  been investigated.  To  solve  this  problem,  the  factors affecting  the  real  contact  area  of  elastic-plastic  rough  surfaces are  revealed  by theoretical research.  Based  on  this, the contact behavior of various elastic-plastic rough  surfaces with  different material and  surface topography is numerically studied, and real contact area versus average contact pressure curves are obtained. With these results, the quantitative relationship between  the evolution of the real contact area and the factors is presented,  which  will provide fundamental theory for scientists and engineers working in this field. Key wordscontact mechanicselastic-plasticrough surfacereal contact areanumerical contact analysis 0   前言* 工程中,宏观光滑的机械结构表面在微观尺度上都是粗糙不平的。使用表面轮廓测量仪等设备观测这些表面时,会发现这些表面具有高低不平的波峰波谷。在精密机电设备的装配过程中,在法向装配接触压力的作用下,两个机械加工的微观粗糙配合面的波峰部分首先发生接触,波谷之间仍然会存*  国家自然科学基金资助项目(51935003)20200708收到初稿,20201102收到修改稿在大量的接触间隙。实际接触区域的面积所占整个配合表面宏观名义接触面积的比例很小,这个比例称为接触率(或无量纲实际接触面积,简称实际接触面积)。实际接触面积的大小与法向载荷、表面微观形貌、材料力学性能等参数相关,并直接影响精密机电设备配合表面的摩擦因数、热导/电导率、接触应力等。粗糙表面实际接触面积演变规律研究、尤其是弹塑性粗糙表面实际接触面积演变规律研究,具有重要的学术价值和工程意义。 粗糙表面接触力学研究起始于1881HERTZ提出的Hertz接触理论[1],在该理论的基础上,网络首发时间:2021-03-05 17:16:46网络首发地址:https://kns.cnki.net/kcms/detail/11.2187.TH.20210304.1641.108.html  机  械  工  程  学  报  第57卷第**期期    2 GREENWOODWILLIAMSON[2]1966年率先提出了一种多微凸体接触计算模型(Greenwood- Williamson模型,简称G-W模型)。在此基础上,有学者对G-W模型进行扩展与改进,提出了很多新的微凸体模型[3-12]。微凸体模型具有计算速度快,轻载条件计算精度高的特点,但是它将实际粗糙表面简化为由高度随机或指数分布的球形或其它形状的微凸体构成的表面,因此会为模型引入了误差,使得这一类模型只有在轻载条件下[13]具有较好的计算精度。为了避免这一问题,PERSSON2001年基于分形理论提出了另外一种不对粗糙表面进行近似处理的纯弹性粗糙表面接触计算方法[14-15],这种方法被称为Persson理论。很多学者通过使用微凸体模型和Persson理论等对轻载的无摩擦的纯弹性粗糙表面接触问题进行研究发现[14-16]:纯弹性粗糙表面的实际接触面积与平均接触压力近似成正比例关系,且比例系数与材料的弹性模量、泊松比以及表面梯度的均方根之间的定量关系。最新研究结果表明,实际接触面积不仅与上述参数相关,而且可能与表面的Nayak参数[17]存在一定的联系。对于弹塑性粗糙表面的接触问题,虽然也有部分学者通过研究发现其实际接触面积随着平均接触压力线性增加[7-9],但是尚未有学者对其进行系统性的研究,弹塑性粗糙表面实际接触面积演变规律及其影响因素尚不明确。 随着计算机技术的提高和数值计算方法的发展,数值计算方法逐渐被应用到接触力学领域,成为一种常用的接触力学研究方法。接触力学数值计算方法包括边界元模型(Boundary element methodBEM) [18-22],半解析法(Semi-Analytical) [23-34]和有限元模型(Finite element methodFEM) [35-45]等。使用上述数值计算方法,BORRI-BRUNETTO[36]HYUN[37-38]分别对具有不同泊松比、表面粗糙度等参数的粗糙表面接触问题进行了研究,也得到实际接触面积和接触压力近似成正比例关系的结论。在弹塑性接触力学领域,虽然YAN[46]GAO[47]PEI[39]YASTREBOV[42]NELIAS [28]通过使用数值计算方法研究发现,相对于纯弹性粗糙表面,弹塑性粗糙表面的实际接触面积与平均接触压力更加接近正比例关系;但是同样尚未有学者针对弹塑性粗糙表面的接触行为进行系统性的研究,其实际接触面积演变规律及其影响因素尚不明确。 针对这一问题,本文首先通过理论研究,确定了弹塑性粗糙表面(理想弹塑性)实际接触面积的影响因素。然后对具有不同材料参数和表面形貌特征的弹塑性粗糙表面的接触行为进行数值仿真计算,得到各种情况下的表面实际接触面积随着平均接触压力的变化曲线。最终得出弹塑性粗糙表面实际接触面积演变规律与这些因素之间的映射关系,并总结得到实际接触面积的计算公式。 1   弹塑性粗糙表面实际接触面积演变规律 1.1   弹塑性粗糙表面实际接触面积演变规律的影响因素 当弹塑性粗糙表面受压接触时,在任意平均接触压力下,表面塑性接触区域相对弹性接触区域的比例为定值。GREENWOODWILLIAMSON[2]定义了塑性指数来描述弹塑性粗糙表面的塑性接触区域相对弹性接触区域的比例,而影响塑性指数的参数为材料的等效弹性模量、硬度(或屈服强度)和表面梯度的均方根。使用多微凸体模型[48-49](例如BGT模型)Persson理论[14-15]对纯弹性粗糙表面的实际接触面积演变规律进行研究都能够得出结论:当两个纯弹性粗糙表面受很小的接触压力载荷时,实际接触面积与平均接触压力近似成正比例关系,其比例系数K仅与材料的弹性模量、泊松比以及表面梯度的均方根有关。而根据Tabor的理论[50]可知:当塑性表面(理想刚塑性材料)与刚性平面受压接触时,无论是否考虑材料塑性变形的不可压缩特性,材料的硬度(或屈服强度)都只是唯一影响接触面积与平均接触压力之间的关系的因素。 因此,在弹塑性接触中,综合考虑弹性接触区域、塑性接触区域以及两区域所占的比例,我们认为影响弹塑性粗糙表面实际接触面积演变规律的主要因素有三个:材料的等效弹性模量、硬度(或屈服强度)和表面梯度的均方根。 1.2   弹塑性粗糙表面接触力学数值计算模型 为了探究上述三个因素对实际接触面积演变规律的影响,需要对具有不同参数类型的弹塑性粗糙表面与刚性平面的接触问题进行接触建模和数值计算。当如图1所示的无摩擦理想弹塑性粗糙表面与刚性平面受压接触时,系统总是趋向于最稳定的状态即总余位能最小的状态[51-52]。因此,计算粗糙表面的接触问题等价于求解系统最小势能的问题,即最小势能理论  ( ) ( ) ( )( ) ( )01min  , , , d2, , dV x y p x y u x yp x y g x yWWd*= W +-W éù ëûòò (1) 2021*月  刘检华等:弹塑性粗糙表面实际接触面积演变规律研究  3 式中,V*表示系统的总余位能,W表示粗糙表面的范围,( ) , p x y表示表面接触压力分布,( ) , u x y表示粗糙表面的法向形变量,d表示两个表面的相对刚体位移,可由式(2)求得  ( ) ( ) 0Contact region=Average , + , g x y u x y d éù ëû  (2)  1   弹塑性粗糙表面与刚性平面接触 ( )0, g x y表示粗糙面与平面的初始接触间隙分布,可由式(3)计算得到  ( ) ( )0, = , 0 g x y h x y -  (3) 式中,( ) , h x y为粗糙面的表面高度,假设粗糙面与平面的初始间隙大于等于0,且刚性平面的初始高度为0。同时,在接触过程中,接触压力和接触间隙需要满足以下约束条件  ( )( )( )( ) ( ),,, 0,, 0,, , 0≤≥≥p x y Hp x yg x yp x y g x y ×=  (4) 式中,H为粗糙表面材料的维氏硬度,接触间隙( ) ( ) ( ) 0, , + , g x y g x y u x y d =-。相对纯弹性边界元算法,这里额外增加了约束条件( ) , p x y H,以限制接触压力不能大于材料的维氏硬度。 根据BOUSSINESQ理论[53]可知,表面形变与表面压力分布相关  ( ) ( )( ) ( ), , ( , )d d,,a a a aa a a au x y G x x y y p x y x yG x y p x y= - - =*òò (5) 式中,“*”为卷积运算符,( ) , G x y表示弹性半无限体表面某点上所受的法向集中力载荷对原点法向位移的影响系数,其表达式为  ( )2221,πvG x yE x y-=+  (6) 式中,E为材料的弹性模量;v为材料的泊松比。 综上所述可知,式(1)仅有一个因变量( ) , V x y*和一个自变量( ) , p x y,且自变量受约束条件式(4)限制,所以弹塑性粗糙面与刚性平面的接触问题的实质是二次泛函条件极值问题。为了实现二次泛函条件极值问题的数值求解,需要对粗糙表面接触问题进行离散处理。如图2所示,将粗糙表面离散为MN´点阵,, xy方向的点间距分别为, xy DD,点( ) , ab处的表面高度表示为( ),, h h x yabab =、接触压力表示为( ),, p p x yabab =,, 0,1, 2 1 N ab=-L。式(1)(5)的离散形式分别为式(7)(8)。  ( )11,,0011,,0001min 2MNMNV p upga b a baba b a babd--*==--==éù =+ëûéù-ëûåååå (7)  11,,00MNiju K pa b a z b h zhzh----===åå  (8) 式中,,ijKa z b h --为( ) , G x x y ya x b h--的离散形式,其近似解为[54-56]  2, 2 2222222222222221lg[ ]lg[ ]lg[ ]lg[ ]lg[ ]lg[ ]lg[ ]lg[ ]vK y x x yEx y x yy x x yx y x yy x x yx y x yy x x yx y x ya z b hp--+ + + ++ + + ++ - - +- + - +- + + -+ - + -- - - -- - - --= + + ++ + -+ + -+ + -+ + -+ + ++ + +++ (9) 式中,2 x x x xaz + = - + D,2 x x x xaz - = - - D,2 y y y ybh + = - + D,2 y y y ybh - = - - D。 同时,式(7)需要满足约束条件  ,,,,,,0,0,0≤≥≥pHpgpgabababa b a b= (10)   机  械  工  程  学  报  第57卷第**期期    4  2   粗糙表面的离散化处理 为了更快速地求解式(7)所示的二次泛函在式(10)所示的约束下的条件极值,本文采用共轭梯度法加速迭代收敛速度[55, 57]。同时,在数值计算过程中,采用矩阵运算对式(8)进行求解会更加方便快捷,该式的矩阵形式为  * u = K p  (11) 式中,uKp都是大小为MN´的矩阵,*表示卷积运算。对矩阵进行线性卷积的数值计算需要( )2MN×次乘法运算,如此庞大的计算量极大地限制了数值求解的计算效率。因此,JU[58]首次引入快速傅里叶变换(Fast Fourier transformFFT)求解接触问题的线性卷积求解,基于FFT的卷积计算大幅度降低了接触问题的计算量。使用快速傅里叶变换求解卷积的公式如下所示  ( ) ( ) IFFT FFT FFT × éù ëû u = K p  (12) 为了避免快速傅里叶变换求解卷积为计算物理域的边界引入较大的边界混叠误差,这里采用WANG[57-59]提出的DC-FFT计算方法,提高卷积运算的计算效率,该方法的主要流程如下。 (1)  将大小为M×N的粗糙表面压力分布矩阵p扩充为大小为(2M-1)×(2N-1)的矩阵pnew,扩充区域用0填充。 (2)  (2M-1)×(2N-1)的扩展域内计算影响系数矩阵K,其中矩阵中心元素表示坐标原点。 (3)  分别对pnew K进行傅里叶变换,得到FFT[pnew]FFF[K](4)  FFT[pnew]FFF[K]进行点乘运算,并对其乘积进行逆傅里叶变换可以得到:unew=IFFT[FFT[pnew]·FFF[K]](5)  , 1, 1 MNnewuua b a b + - + -=,即可得到大小为M× N的目标计算区域的弹性变形矩阵u1.3   弹塑性粗糙表面的实际接触面积演变规律 为了探究上述三个因素对实际接触面积演变规律的影响,本节设计了正交试验对具有不同参数类型的弹塑性粗糙表面与刚性平面的接触问题进行数值计算,正交试验的因素和水平如表1所示。对于每一种参数类型,分别建立了5个弹塑性接触对,并使用第1.2节所述的数值计算建模方法对它们的接触行为进行求解。受模型的计算效率和收敛性的限制,只对各个接触对0%30%的实际接触面积范围的接触状态进行求解。如图3所示,由于每种参数类型的5个接触对的实际接触面积的标准差均小于其均值的1%,每一种参数类型的5个接触对的实际接触面积演变规律之间的差异基本可以忽略不计。因此,本节将每种参数类型的5个接触对的实际接触面积均值作为该参数类型的所有接触对的实际接触面积。 表1    弹塑性接触对的参数类型列表 因素 水平 1  2  3  4  5 表面梯度均方根  0.004 08  0.008 07  0.161  0.335  1.21 等效弹性模量/GPa 10  50  100  —  — 屈服强度/MPa  93.33  466.67  933.33  —  —  图3   参数类型15组接触对的实际接触面积变化曲线 对计算结果进行分析发现,当实际接触面积从0%增加到30%时,所有参数类型的弹塑性粗糙表面的实际接触面积都近似随着平均接触压力月2021*月  刘检华等:弹塑性粗糙表面实际接触面积演变规律研究  5 线性增加(实际接触面积与平均接触压力近似成正比例关系)。使用最小二乘法线性拟合实际接触面积-平均接触压力曲线,得到比例系数K如表2所示。通过对比表中数据可以发现,比例系数K随着屈服强度的增加而降低;但是参数类型23的结果似乎不遵循这个规律(它们的比例系数K相等),这是由于这两种参数类型的接触对塑性接触区域所占比例很低所导致的(塑性系数很小)。同时,比例系数K总体上随着弹性模量的增加而降低、随着表面梯度的均方根的增加而降低。定义w为无量纲比例系数,且KH w=×。结合硬度与屈服强度的关系以及塑性指数的定义,计算出每一个接触对的无量纲比例系数w和塑性指数y如表2所示。 表2   不同参数类型弹塑性接触对的计算结果 编号  比例系数K  无量纲比例系数w  塑性指数y  编号  比例系数K  无量纲比例系数w  塑性指数y 1  4.93×10-3  1.38  1.46  24  4.06×10-4  1.14  2.87 2  3.72×10-3  5.21  0.29  25  3.56×10-3  1.00  57.45 3  3.72×10-3  10.41  0.15  26  7.26×10-4  1.02  11.48 4  3.71×10-3  1.04  7.29  27  3.97×10-4  1.11  5.74 5  9.94×10-4  1.39  1.46  28  3.63×10-3  1.02  11.95 6  8.07×10-4  2.26  0.73  29  8.31×10-4  1.16  2.39 7  3.82×10-3  1.07  14.58  30  5.48×10-4  1.53  1.20 8  8.09×10-4  1.13  2.91  31  3.57×10-3  1.00  59.75 9  5.69×10-4  1.59  1.46  32  7.40×10-4  1.04  11.94 10  4.05×10-3  1.13  2.88  33  3.73×10-4  1.04  5.98 11  1.87×10-3  2.62  0.58  34  3.56×10-3  1.00  119.51 12  1.85×10-3  5.18  0.29  35  7.30×10-4  1.02  23.88 13  3.64×10-3  1.02  14.42  36  3.63×10-4  1.02  11.95 14  8.05×10-4  1.13  2.88  37  3.57×10-3  1.00  43.18 15  4.95×10-4  1.39  1.44  38  7.36×10-4  1.03  8.63 16  3.58×10-3  1.00  28.83  39  3.84×10-4  1.07  4.32 17  7.52×10-4  1.05  5.76  40  3.55×10-3  0.99  215.92 18  4.06×10-4  1.14  2.88  41  7.13×10-4  1.00  43.14 19  3.75×10-3  1.05  5.74  42  3.59×10-4  1.00  21.59 20  1.09×10-3  1.53  1.15  43  3.55×10-3  0.99  431.83 21  9.35×10-4  2.62  0.57  44  7.11×10-4  1.00  86.27 22  3.58×10-3  1.00  28.72  45  3.57×10-4  1.00  43.18 23  7.78×10-4  1.09  5.74          无量纲比例系数w随塑性指数y的变化如图4所示(为了便于绘图和观察,该图采用了双对数坐标)。w随y的变化可以近似用图中曲线进行拟合,该曲线的关系式为  图4   无量纲比例系数w随塑性指数y的变化曲线 ( ) ( )( ) ( )( )65 4320.016 9 0.110 80.210 8 0.038 50.555 6lg lg lglg l0.623 7 0.238 2glg lgw y yyyyy- ´ + ´ -´ - ´ +´ - ´ +=  (13) 结合式(13)以及w的定义,即可近似得到弹塑性粗糙表面受压接触时的实际接触面积随着平均接触压力变化的关系式为  10ApaH=   ( ) ( )( ) ( )( )65432lg 0.016 9 0.110 80.210 8 0.03lg lglg l 850.555 6 0.623 7 0.238 2glg lgA yyyyyyw = = - + --+-+  (14)   机  械  工  程  学  报  第57卷第**期期    6 2   试验验证 为了验证本文所总结的弹塑性粗糙表面实际接触面积演变规律,使用式(13)对由如图5所示的各向同性随机粗糙表面所组成的接触对的弹塑性接触问题进行求解,并将计算结果与试验测量[60-61]和有限元仿真结果[61]进行了对比研究。该接触对的等效弹性模量为340 66 MPa,表面硬度为310 MPa,表面梯度的均方根为0.710 9。表3将本文计算的比例系数K与试验测量结果和有限元仿真结果进行了对比,本文计算结果与试验测量和有限元仿真结果的相对误差均可忽略不计,证明了本文所总结的弹塑性粗糙表面实际接触面积演变规律的准确性。  图5   各向同性随机粗糙表面 表3  本文计算结果与试验测量和有限元仿真结果的对比   本文计算  试验测量  有限元仿真 比例系数K  0.003 20  0.003 31  0.003 16 相对误差(%)  —  3.32  1.27 3   结论 (1) 针对弹塑性粗糙表面实际接触面积的演变规律尚不明确的问题,结合多微凸体模型、Persson理论、塑性接触理论和塑性指数的研究得出结论:影响弹塑性粗糙表面接触状态的主要因素为材料的等效弹性模量、硬度(或屈服强度)和表面梯度的均方根。 (2) 在此基础上,通过对多组弹塑性粗糙表面的接触行为进行数值仿真计算,发现当实际接触面积小于30%时,弹塑性粗糙表面的实际接触面积与平均接触压力近似成正比例关系,揭示了该正比例系数与材料的弹性模量、屈服强度和表面梯度的均方根之间的映射关系,并最终得到实际接触面积的计算公式,为弹塑性接触力学的后续工程应用奠定了基础。 参  考  文  献   [1]  HERTZ  H.  On  the  contact  of  elastic  solids[J].  Journal Reine. Angew. Math.188192(1)156-171.   [2]  GREENWOOD  J AWILLIAMSON  J  B  P.  Contact  of nominally  flat  surfaces[J].  Proceedings  of  the  Royal Society of London1966295(1442)300-319.   [3]  GREENWOOD  J  ATRIPP  J  H.  The  elastic  contact  of rough spheres[J]. Journal of Applied Mechanics196734(1)153-160.   [4]  HISAKADO  T.  Effect  of  surface  roughness  on  contact between solid surfaces[J]. Wear197428(2)217-234.   [5]  GREENWOOD  J  A.  The  contact  of  two  nominally  flat rough  surfaces[J].  Proceedings  of  the  Institution  of Mechanical Engineers1970185(1970)625-634.   [6]  BUSH  A  WGIBSON  R  DKEOGH  G  P.  Strongly anisotropic rough surfaces[J]. Journal of Tribology1979101(1)15-21.   [7]  CHANG W RETSION IBOGY D B. An elastic-plastic model  for  the  contact  of  rough  surfaces[J].  Journal  of Tribology1987109(2)257-263.   [8]  ZHAO  YMAIETTA  D  MCHANG  L.  An  asperity microcontact  model  incorporating  the  transition  from elastic  deformation  to  fully  plastic  flow[J].  Journal  of Tribology2000122(1)86-93.   [9]  ZHAO YCHANG L. A model of asperity interactions in elastic-plastic  contact  of  rough  surfaces[J].  Journal  of Tribology2001123(4)857-864. [10]  KOGUT LETSION I. Elastic-plastic contact analysis of a sphere and a rigid flat[J]. Journal of Applied Mechanics200269(5)657-62. [11]  JACKSON  R  LGREEN  I.  A  finite  element  study  of elasto-plastic hemispherical contact against a rigid flat[J]. Journal of Tribology2005127(2)343-54. [12]  WADWALKAR S SJACKSON R LKOGUT L. A study of  the  elasticplastic  deformation  of  heavily  deformed spherical  contacts[J].  ARCHIVE  Proceedings  of  the Institution  of  Mechanical  Engineers  Part  J  Journal  of Engineering Tribology2010208-210(1)1-12. [13]  YASTREBOV V AANCIAUX GMOLINARI J F. From infinitesimal  to  full  contact  between  rough  surfacesEvolution  of  the  contact  area[J].  International  Journal  of Solids & Structures20155283-102. [14]  PERSSON  B  N.  Elastoplastic  contact  between  randomly rough surfaces[J]. Phys. Rev. Lett.200187(11)116101. 2021*月  刘检华等:弹塑性粗糙表面实际接触面积演变规律研究  7 [15]  PERSSON B NBUCHER FCHIAIA B. Elastic contact between randomly rough surfacesComparison of theory with  numerical  results[J].  Physical  Review  B200265(18)184106. [16]  CARBONE  GBOTTIGLIONE  F.  Asperity  contact theoriesDo  they  predict  linearity  between  contact  area and  load?  [J].  Journal  of  the  Mechanics  &  Physics  of Solids200856(8)2555-2572. [17]  YASTREBOV V AANCIAUX GMOLINARI J F. The role of the roughness spectral breadth in elastic contact of rough surfaces[J].  Journal of the Mechanics & Physics of Solids2017107469-493. [18]  HU Y ZBARBER GDONG Z. Numerical analysis for the  elastic  contact  of  real  rough  surfaces[J].  Tribology Transactions199942(3)443-452. [19]  STANLEY  H  MKATO  T.  An  FFT-based  method  for rough  surface  contact[J].  Journal  of  Tribology1997119(3)481-485. [20]  LI  JBERGER  E  J.  A  semi-analytical  approach  to three-dimensional  normal  contact  problems  with friction[J].  Computational  Mechanics200330(4)310-322. [21]  POHRT RLI Q. Complete boundary element formulation for  normal  and  tangential  contact  problems[J].  Physical Mesomechanics201417(4)334-340. [22]  CARBONE GPUTIGNANO C. A novel methodology to predict  sliding  and  rolling  friction  of  viscoelastic materialsTheory  and  experiments[J].  Journal  of  the Mechanics & Physics of Solids201361(8)1822-1834. [23]  CHEN W WLIU SWANG Q J. Fast Fourier transform based  numerical  methods  for  elasto-plastic  contacts  of nominally flat surfaces[J]. Journal of Applied Mechanics200875(1)011022. [24]  WANG W ZWANG Z JWANG W Zet al. A numerical elastic-plastic  contact  model  for  rough  surfaces[J]. Tribology Transactions201053(2)224-238. [25]  WRIGGERS  P.  Computational  contact  mechanics[J]. Computational Mechanics200332(1)141. [26]  JACQ  CNE´  LIAS  DLORMAND  Get  al. Development  of  a  three-dimensional  semi-analytical elastic-plastic  contact  code[J].  Journal  of  Tribology2002124(4)653-667. [27]  BOUCLY VNE´ LIAS DLIU Set al. Contact analyses for  bodies  with  frictional  heating  and  plastic  behavior[J]. Journal of Tribology2005127(2)355-364. [28]  NELIAS  DBOUCLY  VBRUNET  M.  Elastic-plastic contact  between  rough  surfacesProposal  for  a  wear  or running-in model[J]. Journal of Tribology2006128236-244. [29]  NéLIAS  DANTALUCA  EBOUCLY  Vet  al.  A three-dimensional semi-analytical model for elastic-plastic sliding  contacts[J].  Journal  of  Tribology2007129761-771. [30]  CHEN W WWANG Q J. Thermomechanical analysis of elastoplastic  bodies  in  a  sliding  spherical  contact  and  the effects  of  sliding  speedheat  partitionand  thermal softening[J]. Journal of Tribology2008130(4)041402. [31]  WANG Z JWANG W ZHU Y Zet al. A numerical elasticplastic  contact  model  for  rough  surfaces[J]. Tribology Transactions201053(2)224-238. [32]  ZHOU KCHEN W WKEER L Met al. A fast method for solving three-dimensional arbitrarily shaped inclusions in  a  half  space[J].  Computer  Methods  in  Applied Mechanics and Engineering2009198(9-12)885-892. [33]  LIU  SWANG  Q.  Elastic  fields  due  to  eigenstrains  in  a half-space[J]. Journal of Applied Mechanics200572871-878. [34]  LIU SJIN XWANG Zet  al. Analytical  solution for elastic  fields  caused  by  eigenstrains  in  a  half-space  and numerical  implementation  based  on  FFT[J].  International Journal of Plasticity201235135-154. [35]  KOMVOPOULOS  KCHOI D H.  Elastic finite element analysis  of  multi-asperity  contacts[J].  Journal  of Tribology1992114(4)823-831. [36]  BORRI-BRUNETTO MCARPINTERI ACHIAIA B. Lacunarity  of  the  contact  domain  between  elastic  bodies with rough boundaries[M]. NetherlandsSpringer1998. [37]  HYUN SPEI LMOLINARI J Fet al. Finite-element analysis of contact  between elastic self-affine surfaces[J]. Phys.  Rev.  E  Stat.  Nonlin.  Soft.  Matter.  Phys.200470(2)026117. [38]  HYUN SROBBINS M O. Elastic contact between rough surfacesEffect  of  roughness  at  large  and  small wavelengths[J]. Tribology International2007,  40(10)1413-1422. [39]  PEI LHYUN SMOLINARI J Fet al. Finite element modeling  of  elasto-plastic  contact  between  rough surfaces[J].  Journal  of  the  Mechanics  &  Physics  of Solids,  200553(11)2385-2409. [40]  SAHOO PGHOSH N. Finite element contact analysis of fractal surfaces[J]. Journal of Physics D Applied Physics200740(14)4245. [41]  WALTER CMITTERER C. 3D versus 2D finite element simulation  of  the  effect  of  surface  roughness  on nanoindentation  of  hard  coatings[J].  Surface  &  Coatings Technology2009203(20)3286-3290. [42]  YASTREBOV V ADURAND JPROUDHON Het al. Rough  surface  contact  analysis  by  means  of  the  finite   机  械  工  程  学  报  第57卷第**期期    8 element  method  and  of  a  new  reduced  model[J]. Mechanical Reports2011339(7)473-490. [43]  BRYANT  M  JEVANS  H  PSNIDLE  R  W.  Plastic deformation in rough surface line contacts-a finite element study[J]. Tribology International201246(1)269-278. [44]  WRIGGERS  PREINELT  J.  Multi-scale  approach  for frictional contact of elastomers on rough rigid surfaces[J]. Computer Methods in Applied Mechanics & Engineering2010198(21)1996-2008. [45]  ZHANG  FLIU  JDING  Xet  al.  An  approach  to calculate  leak  channels  and  leak  rates  between  metallic sealing surfaces[J]. Journal of Tribology2017139(1)011708-18. [46]  YAN  WKOMVOPOULOS  K.  Contact  analysis  of elastic-plastic  fractal  surfaces[J].  Journal  of  Applied Physics199884(7)3617. [47]  GAO Y FBOWER A F. Elastic-plastic contact of a rough surface  with  Weierstrass  profile[J].  Proceedings Mathematical  Physical  &  Engineering  Sciences2006462(2065)319-348. [48]  BUSH A WGIBSON R DTHOMAS T R. The elastic contact of a rough surface[J]. Wear197535(1)87-111. [49]  GREENWOOD J A. A simplified elliptic model of rough surface contact[J]. Wear2006261(2)191-200. [50]  TABOR D.  The hardness of metals[M]. LondonOxford University Press1951. [51]  RICHARDS T H. Energy methods in stress analysisWith an  introduction  to  finite  element  techniques[M]. New YorkHalsted Press1977. [52]  TIAN  XBHUSHAN  B.  A  numerical  three-dimensional model  for  the  contact  of  rough  surfaces  by  variational principle[J]. Journal of Tribology19961118(1)33-42. [53]  JOHNSON  K  L.  Contact  mechanics[M].  CambridgeCambridge University Press1985. [54]  VENNER C HLUBRECHT A A. Multi-level methods in lubrication[M]. New York:  Elsevier2000. [55]  SAINSOT PLUBRECHT A A. Efficient solution of the dry  contact  of  rough  surfacesA  comparison  of  fast Fourier  transform  and  multigrid  methods[J].  Proc.  Inst. Mech. Eng. Part J-J. Eng. Tribol.2011225(J6)441-448. [56]  王文中.  混合润滑的数值模拟和实验研究[D]. 北京:清华大学,2003.   WANG  Wenzhong.  Numerical  simulation  and  exper- imental study of mixed lubrication [D]. BeijingTsinghua University2003. [57]  王文中,胡元中,王慧.  基于快速傅里叶变换和共轭梯度法求解干接触问题[J].  机械工程学报,200642(7)14-18.   WANG WenzhongHU YuanzhongWANG Hui. Solving dry contact problems based on fast Fourier transform and conjugate  gradient  method  [J].  Journal  of  Mechanical Engineering200642(7)14-18. [58]  JU  YFARRIS  T.  Spectral  analysis  of  two-dimensional contact  problems[J].  Journal  of  Tribology1996118320-328. [59]  WANG WWANG HLIU Yet al. A comparative study of  the  methods  for  calculation  of  surface  elastic deformation[J].  Proceedings  of  the  Institution  of Mechanical  EngineersPart  JJournal  of  Engineering Tribology2003217(2)145-154. [60]  ZHANG FLIU JDING Xet al. A discussion on the capability  of  x-ray  computed  tomography  for  contact mechanics  investigations[J].  Tribology  International2020145106167. [61]  ZHANG FLIU JDING Xet al. Experimental and finite element  analyses  of  contact  behaviors  between non-transparent  rough  surfaces[J].  Journal  of  the Mechanics and Physics of Solids201912687-100.   作者简介:刘检华,男,1977年出生,博士,教授,博士研究生导师。主要研究方向为产品装配与检测技术。 E-mailjeffliu@bit.edu.cn 张飞凯,男,1993年出生,博士。主要研究方向为接触力学、静密封。 E-mailzhangfkbit@163.com 丁晓宇(通信作者),男,1984年出生,博士,副教授,硕士研究生导师。主要研究方向为基于机器视觉的装配精度测量、基于接触力学的装配变形预测、基于数字孪生的装配精度在线跟踪、基于信号处理和机器学习的装配质量控制等。 E-mailxiaoyu.ding@bit.edu.cn   

[返回]
上一篇:多维波形编码SAR俯仰向级联DBF成像性能分析_英文_何峰
下一篇:地铁车辆制动管路动应力分析及结构优化_谢晨希