纳米单晶铜杆压缩时的稳定性研究.docVIP

纳米单晶铜杆压缩时的稳定性研究 张任良 张文飞 燕山大学建筑工程与力学学院，秦皇岛，066004 摘要：应用嵌入势函数和分子动力学方法研究了纳米单晶铜杆压缩时的稳定性。通过一端固定，另一端施加轴向作用力驱使原子运动，模拟了纳米单晶铜杆失稳弯曲破坏的过程，详细分析了弛豫阶段表面效应的影响。画出了材料压缩时的应力应变关系，得到其压缩时的弹性模量以及失稳时的临界载荷，求出其临界应力。 关键词：分子动力学；单晶铜；表面效应；动力屈曲 引言 随着微电机系统技术的迅速发展[1]，人们不断追求制造出尺度越来越小、性能越来越完善的微型装置。特别是在生物工程、医疗器械、航空航天、数字通信、传感技术以及灵巧武器等领域中机械装置微型化的要求日益增长。了解并控制这些原子尺度的力学现象显得十分必要。所以在纳米尺度下的材料的力学及物理性能的研究显得日益迫切。从统计物理学衍生出来的分子动力学模拟，能够确定一群有相互作用的粒子在时空中的演化规律，这对澄清纳米材料微观结构有重要作用，现在被公认为是一种描述原子和分子级微观现象的有效方法[2~3]。在低温条件下有些学者已经研究了金属纳米杆受压稳定性[4~5][6] （1） （2） （3） （4） （5） （6） （7） 这里，表示系统得总能量，为原子所具有的势能，是其他原子对原子的电子云密度，代表吸引项。 图3 初始二维横截面 所模拟的纳米杆尺寸为6a×6a×30a(2.166nm×2.166nm×10.83nm)，原子个数为5155．初始三维构型如图1，原子按照面心立方(fcc)理想晶格点阵排列，纵向二维y-z截面初始构形如图2，为了看清楚原子标号，放大横向x-y截面如图3。在模拟过程中，动力学方程积分格式采用Beeman算法[7]，初始构型作为粒子的初始位置，而粒子的初始速度服从Maxwell-Boltzmann分布，为了加快计算效率采用verlet列表法[8]，用直接速度标定方法[9~10]进行等温调节，温度控制在293K。时间步长取为1fs。 弛豫过程 先对理想初始构型进行弛豫60000步，势能随时间变化如图4，试件长度随时间变化如图5。从图4和图5中可以看出在30000步以后，构型进入稳定阶段，从30000步开始，采用时间平均[11]，得到势能趋于稳定值如图6，同样试件的长度趋于稳定值10.66nm如图7，初始构型达到能量的稳定状态总势能为-17040.9eV，平均每个原子的势能为－3.3eV，与初始构型相比弛豫后试件缩短了1.57％。 由势函数公式得到的原子应力[12]表达式为： （8） 根据对称原则，可以把图3中的标号原子15，16，17，18作为第一组，标号原子1，7，8，14作为第二组是对称的，标号原子2，6，9，13为第三组，标号原子3，5，10，12为第四组，4号原子单独为第5组，从30000步以后，在各组内取时间平均得到5条曲线如图8。 图4 弛豫过程中势能随时间的变化 图5 试件长度随时间的变化 图6 时间平均得到势能 图7 时间平均得到试件长度 图8 五组原子应力的平均值的时间平均 从图8中可以看出，第三组，第四组，第五组的原子（认为是体内原子）应力大小相近约为-3.2Gpa。第二组原子应力也就是横截面正方形四边中点原子的应力次之约为7.01Gpa是体内原子应力的2倍。第一组原子应力即四个角点原子的应力最大约为10.6Gpa是体内原子的3.3倍。这是因为表面原子相对于体内原子，表面上的原子缺少一个或数个近邻原子，配位不足，有大量的悬键或不饱和键存在。从影响范围看表面效应影响到向内的两层原子。而从应力的正负号中可看出整个杆件呈现出表面受拉内部受压的应力状态。 加载过程 初始构型