纳米级粗糙峰法向接触特性研究.pdf
《纳米级粗糙峰法向接触特性研究.pdf》由会员分享,可在线阅读,更多相关《纳米级粗糙峰法向接触特性研究.pdf(8页珍藏版)》请在咨信网上搜索。
1、为了突破有限元法在微观尺度下界面接触特性分析中的局限性,采用分子动力学方法研究纳米级粗糙峰的法向接触行为。利用大规模原子/分子并行模拟器(L AMMP S)建立3种不同半径的C u单粗糙峰与金刚石刚性平面法向接触的分子动力学模型,从原子尺度模拟纳米级粗糙峰法向接触动态过程,分析法向接触力、实际接触面积、接触应力等参量的动态变化特性。结果表明:粗糙峰半径一定时,法向接触力、接触应力、实际接触面积与接触变形率正相关;接触变形率一定时,法向接触力、接触应力、实际接触面积均随粗糙峰半径的增大而递增;不同于宏观接触行为,由于原子间相互作用力,纳米级粗糙峰在法向初始接触阶段出现负的法向接触力以及非零的实际
2、接触面积;粗糙峰的接触变形过程伴随着原子位错以及应力集中。关键词:纳米级粗糙峰;单粗糙峰;分子动力学;法向接触特性中图分类号:O 3 1 3.5;TH 1 1 7.1 文献标志码:A 文章编号:1 6 7 4-3 6 4 4(2 0 2 3)0 4-0 2 8 1-0 8 机械零部件的相互连接形成接触界面,其接触特性对于力和扭矩的传递至关重要1-2。基于连续介质力学的接触理论可以解释宏观尺度上的界面接触问题,但无法用于微观尺度的界面接触行为分析。微观尺度下的接触界面不是光滑的,而是 由具有不同 高度和曲率 半径的粗糙 峰构成3。在法向负载作用下,粗糙峰之间会先接触,导致实际的接触面积远小于接触
3、界面的面积。研究粗糙峰的法向接触特性能够揭示微观与宏观接触行为的差异性,从而为机械连接的可靠性分析提供理论依据。单粗糙峰的界面接触行为分析是揭示整个粗糙界面接触特性的基础,相关研究成果有不少。J o h n s o n等4对经典H e r t z接触理论进行了修正,得出 两 弹 性 球 体 的 黏 着 接 触 模 型(J K R模型)。C h a n g等5将单微凸体接触变形区域分为弹性区和塑性区,将接触过程定义为弹性接触和塑性接触两个阶段,并提出了C E B接 触模型。K o g u t等6基于有限元方法提出了粗糙表面的弹塑性接触模型(K E模型)。A b d o7考虑弹性、塑性等多因素改进
4、了单粗糙峰的弹塑性接触模型,并基于极限应力概念对粗糙表面进行有限元分析。J a c k s o n等8-9考虑了塑性硬化,并依此建立了 单一粗糙表面接触模型(J G模型)。M e g a l i n-g a m等1 0基于单粗糙峰的高斯分布接触模型推导出弹性、弹塑性和全塑性状态下接触载荷与接触面积的经验表达式。分子动力学方法突破了有限元法的尺度局限性,被广泛用于微纳尺度的接触特性研究。L u a n等1 1建立了球凸面接触的分子动力学模型,研究法向负载和黏着度对粗糙峰接触特性的影响规律,探讨了连续介质力学理论在纳米尺度的适用性。S i等1 2-1 3模拟多粗糙峰的滑动接触过程,分析了粗糙峰数量
5、、曲率对黏附力及摩擦力的影响规律。G u1 4等基于单粗糙峰的法向加载、卸载过程,揭示出纳米尺度下C u-C u引线键合的微焊形成和断裂机理。C h e o n g等1 5分析了粗糙峰的位置变化对纳米级表面接触磨损的影响。C h a等1 6研究了单粗糙峰及其在两个表面接触时的变形。J a c o b s等1 7采用分子动力学方法,结合仿真和实验研究了原子级粗糙度对碳基材料与类金刚石(D L C)或超纳米晶金刚石的纳米级微凸体之间黏附的影响。目前,对纳米尺度下粗糙峰接触参量变化特性的研究较少。本文基于分子动力学原理建立纳米级单粗糙峰的法向接触模型,模拟不同半径的单粗糙峰与刚性平面的法向接触变形过
6、程,研究法向接触力、接触应力、实际接触面积、接触变形武汉科技大学学报2 0 2 3年第4期量等接触参量的动态变化特性,揭示纳米级粗糙峰的法向接触变形机理,以期为微纳尺度下的界面接触行为分析提供理论参考。1 单粗糙峰法向接触的分子动力学模型采 用 大 规 模 原 子/分 子 并 行 模 拟 器(L AMMP S)进行编程,建立单粗糙峰与刚性平面的法向接触力学模型,如图1所示。单粗糙峰法向接触模型的基本尺寸为1 0n m1 0n m8.5n m。模型由上下两部分构成:上半部分为刚性平面,即刚性层,由晶格间距为0.3 5 7n m的金刚石立方晶胞C原子组成;下半部分为单粗糙峰,包含牛顿 层、恒 温
7、层 和 边 界 层,均 由 晶 格 间 距 为0.3 6 1 5n m的面心立方晶胞C u原子构成。刚性层、恒温层以及边界层均为规则基体,且z向高度均为1n m,牛顿层则是半径为R的半球体。通过改变R值建立不同半径的单粗糙峰法向接触模型,如图2所示,其中从左至右,粗糙峰半径R依次为2、3、4n m,相关参数在表1中列出。图1 单粗糙峰法向接触的分子动力学模型F i g.1M o l e c u l a rd y n a m i c sm o d e l o f s i n g l e-a s p e r i t yn o r m a l c o n t a c t图2 不同半径的单粗糙峰法向接
8、触模型F i g.2S i n g l e-a s p e r i t yn o r m a l c o n t a c tm o d e l sw i t hd i f f e r e n t r a d i i表1 不同半径的单粗糙峰法向接触模型相关参数T a b l e1R e l e v a n tp a r a m e t e r so f s i n g l e-a s p e r i t yn o r m a l c o n-t a c tm o d e l sw i t hd i f f e r e n t r a d i iR/n m牛顿层原子数/个总原子数/个213 2 8
9、1 75 3 0346 0 02 34 1 641 10 1 82 98 3 4 图3所示为模拟体系的条件设置。为消除因计算规模引起的边界效应,将x、y向设置为周期性边界条件,将z向设置为非周期性边界条件;图3 模拟体系的条件设置F i g.3C o n d i t i o ns e t t i n g f o r s i m u l a t i o ns y s t e m采用速度标定法将恒温层温度控制为3 0 0K;通过正则系综(NV E)约束模拟体系以保证系统的稳态;模拟时间步长为0.0 0 5p s。刚性层被视为刚性体,刚性层内部C原子间的相互作用力可忽略;Em b e d d e d
10、A t o m M e t h o d(E AM)势通常用于描述金属原子间的相互作用力,因而下半部分C u原子间相互作用力选用C u_u 3.e a m势描述;法向接触过 程 中,接 触 界 面 原 子 间 相 互 作 用 力 为M o r s e势,其表达式为(d)=De(-2(d-d0)-e(-(d-d0)(1)式中:D为结合能参数,为势能曲线梯度系数,d为原子间瞬时距离,d0为截止距离。根据洛伦 兹-贝塞洛定律 确定接触界 面的M o r s e势参数:DC-C u=(DCDC u)1/2(2)C-C u=(CC u)1/2(3)dC-C u=(dCdC u)1/2(4)由式(2)式(4
11、)得到M o r s e势函数参数值:DC-C u=0.0 8 7e V,C-C u=0.5 1 4n m,dC-C u=0.2 0 5n m。2822 0 2 3年第4期涂福泉,等:纳米级粗糙峰法向接触特性研究 法向接触前,单粗糙峰法向接触的分子动力学模型需进行充分弛豫,以确保模拟系统处于平衡状态,整个弛豫阶段持续2 0p s。在充分弛豫阶段,刚性层以V=1 0m/s的速度沿z轴负向移动。随着刚性层位移Xz的增大,刚性层与牛顿层开始接触,即单粗糙峰与刚性平面间发生接触,如图4所示。图4 单粗糙峰法向接触过程F i g.4N o r m a l c o n t a c tp r o c e s
12、 so f s i n g l ea s p e r i t y 用接触变形率表示接触过程中粗糙峰的变形程度:=u/R(5)式中:u为等效接触变形量。用面积接触率S表示粗糙峰与刚性基板的接触程度:S=Sr/S0(6)式中:Sr为实际接触面积;S0为名义几何面积,这里定义为粗糙峰底面积。法向接触力FN实质为接触界面原子间相互作用力的法向分量,故将刚性层与牛顿层间相互作用力F的z向分量定义为单粗糙峰与刚性平面的法向接触力FN。2 纳米级单粗糙峰法向接触特性分析2.1 法向接触力通过模拟不同半径的单粗糙峰法向接触过程,得到法向接触力FN随接触变形率的变化曲线,如图5所示。可以看出,粗糙峰半径R相同时
13、,法向接触力FN随接触变形率的增加总体呈增大趋势,其中,粗糙峰半径R=4n m时,FN随的变化率最大。由于粗糙峰与刚性平面间的相互接触表现为接触原子间相互作用力F,接触变形程度加大时,粗糙峰与刚性平面间的接触原子数量随之增多,故而法向接触力FN也增大。另外,接触变形程度相同时,法向接触力FN与粗糙峰半径R正相关。这是因为相同时,接触原子数量与粗糙峰半径R正相关,致使法向接触力FN也随粗糙峰半径R的增大而增大。图6所示为不同半径粗糙峰的法向接触力FN随接触变形量u的变化情况,可以看出,FN图5 法向接触力FN随接触变形率的变化F i g.5V a r i a t i o no fn o r m
14、a l c o n t a c t f o r c eFNw i t hc o n t a c t d e-f o r m a t i o nr a t e(a)R=2n m(b)R=3n m(c)R=4n m图6 法向接触力FN随接触变形量u的变化F i g.6V a r i a t i o no fn o r m a l c o n t a c t f o r c eFNw i t hc o n t a c t d e-f o r m a t i o nu382武汉科技大学学报2 0 2 3年第4期与u总体呈非线性正相关。经典H e r t z接触理论揭示了弹性接触阶段法向接触力与接触变形
15、量的变化关系,根据该理论有如下公式成立:P=ERr1/2(7)式中:P为接触载荷,r为接触半径,E为有效弹性模量。结合H e r t z接触理论构建曲线y=k x3/2,然后与通过分子动力学模拟得到的FN-u变化曲线进行对比,来验证本文建立的纳米级单粗糙峰法向接触模型的有效性。基于最小二乘法进行曲线拟合(结果见图6),拟合参数如表2所示。可以看到,随着粗糙峰半径R由2n m增至4n m,系数k值也在增大,并且3个拟合曲线的决定系数较大(0.5 50.9 3),表明分子动力学模拟结果与H e r t z接触理论公式的重合度较高,证明所建立的纳米级单粗糙峰法向接触模型是有效的。表2 FN-u曲线拟
16、合参数T a b l e2F i t t i n gp a r a m e t e r so fFN-uc u r v e sR/n mk决定系数27 7.4 10.5 531 0 1.9 00.7 841 4 2.1 40.9 3 不同于宏观接触行为,纳米尺度下初始法向接触力FN出现负值(见图5、图6小方框处),例如在00.0 5区间变化时,法向接触力FN在零点附近波动。原子间相互作用力为原子斥力和原子引力的共同作用,原子间距直接影响原子间相互作用力,原子斥力与原子引力均随原子间距减小而增大,且原子斥力变化速率更大。单粗糙峰在初始接触状态下,接触原子引力作用更加明显,这时表现为负的法向接触力
17、;接触变形程度较大时,接触原子斥力作用更加明显,致使法向接触力FN为正;粗糙峰半径较大时,法向接触力的变化愈加明显,而接触原子的引力相对斥力变化较小,致使初始法向接触力为负值的区域较少。此外,原子热运动直接影响原子间距,导致了法向接触力FN的波动。2.2 接触变形下面通过分析单粗糙峰的应力和晶体结构变化来研究其接触变形机理。法向接触变形过程中,不同半径单粗糙峰的法向接触应力H分布如图7所示。随着接触变形程度的加深,接触区域应力值不断增大,应力集中主要发生在接触体内部的中心区域,且应力值由内向外递减。原子间距与接触变形程度负相关,这就意味着,随着接触变形率的增大,原子间距会减小,原子间作用力会增
18、大,由于原子斥力作用变化速率相对更大,因此表现出正的接触应力。中心区域的原子间相互作用相对于边缘区域表现得更加明显,因而中图7 不同半径单粗糙峰的法向接触应力分布F i g.7N o r m a l c o n t a c t s t r e s sd i s t r i b u t i o no f s i n g l ea s p e r i t i e sw i t hd i f f e r e n t r a d i i4822 0 2 3年第4期涂福泉,等:纳米级粗糙峰法向接触特性研究心区域接触应力相对较高。粗糙峰半径R较大时,粗糙峰内部原子相互作用力也较大,表现为更高的接触应力。法
19、向接触过程中的原子结构演变如图8所示,可以看出,在接触变形过程中,粗糙峰的原子结构由面心立方晶格结构(F C C)向体心立方晶格结构(B C C)、六方密堆积晶体结构(HC P)、非晶结构(AMO)演变,其中面心立方结构原子所占比例如图9所示。由图9可知,面心立方结构原子比例减小速率与粗糙峰半径R正相关,与图7中接触应力分布的变化相一致,这是因为原子结构受原子间距影响,即原子发生了位错。图8 不同半径单粗糙峰在法向接触过程中的原子结构演变F i g.8A t o m i c s t r u c t u r ee v o l u t i o no f s i n g l ea s p e r i
20、 t i e sw i t hd i f f e r e n t r a d i i d u r i n gn o r m a l c o n t a c tp r o c e s s图9 法向接触过程中面心立方晶格结构原子占比的变化F i g.9V a r i a t i o no f t h ep r o p o r t i o no fF C Ca t o m sd u r i n gn o r-m a l c o n t a c tp r o c e s s 图1 0所示为接触变形过程中不同半径单粗糙峰的原子结构变化情况。随着变形程度的增加,F C C结 构 的 原 子 数 量 非 线
21、 性 减 少,B C C、HC P、AMO结构的原子数量非线性增加。这是因为接触变形过程中,原子间相互作用力改变了原子的初始排布,致使其结构发生位错,同时,原子热运动效应导致原子结构变化出现波动。2.3 实际接触面积接触面积是反映粗糙峰法向接触特性变化的重要参量,单粗糙峰在法向接触过程中接触面积的变化情况如图1 1所示。图1 1(a)揭示了实际接触面积Sr随接触变形率的变化规律:粗糙峰半径R不变时,实际接触面积Sr随接触变形率近似线性增长;接触变形率不变时,实际接触面积Sr与粗糙峰半径R正相关。面积接触率S随接触变形率的变化规律如图1 1(b)所示,粗糙峰半径R不变时,面积接触率S与实际接触面
22、积Sr具有相同的变化趋势,即随接触变形率近似线性增长;粗糙峰半径R由2n m增至4n m时,面积接触率S的变化曲线下移。因为实际接触面积(a)R=2n m582武汉科技大学学报2 0 2 3年第4期(b)R=3n m(c)R=4n m图1 0 不同半径单粗糙峰在法向接触过程中的原子结构变化F i g.1 0A t o m i cs t r u c t u r ev a r i a t i o no fs i n g l ea s p e r i t i e sw i t hd i f f e r e n t r a d i i d u r i n gn o r m a l c o n t a
23、c tp r o c e s sSr、面积接触率S与接触原子数量正相关,接触变形程度加深时,接触原子数量增大,实际接触面积Sr、面积接触率S也随之增大,而接触变形程度相同时,接触原子数量随粗糙峰半径R的增大而增多,法向接触过程中的实际接触面积Sr也越来越大。面积接触率S受粗糙峰半径R的影响,这是因为粗糙峰内部原子间相互作用力与粗糙峰半径R正相关,即粗糙峰半径R越大时,原子间扩散作用愈加微弱,粗糙峰抵抗形变能力越强,面积接触率S随之减小。(a)Sr-变化曲线(b)S-变化曲线图1 1 不同半径单粗糙峰在法向接触过程中的接触面积变化F i g.1 1C o n t a c ta r e av a
24、r i a t i o no fs i n g l ea s p e r i t i e sw i t hd i f-f e r e n t r a d i i d u r i n gn o r m a l c o n t a c tp r o c e s s 与连续介质力学理论方法不同,采用分子动力学方法进行分析时,单粗糙峰在法向接触的初始阶段,其接触面积变化曲线的起始点并非零点,见图1 1小方框处。参照文献1 8,通过原子间距阈值0.5n m来定义原子接触,即原子间距小于0.5n m时发生接触,反之则未发生接触。当接触变形率=0时,虽然未发生如宏观物体那样的直接接触,但牛顿层原子与刚性层原
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 纳米 粗糙 接触 特性 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。