论文部分内容阅读
金属氢脆的现象早已为人们所认知,而关于氢脆的机理仍然存在着争议。钛金属是适用于航空等尖端领域的材料,其氢脆的发生易造成较大的经济损失。计算模拟有助于在实验的基础上更深入的研究钛金属氢脆的机理,从而更好的避免氢脆的发生。计算模拟中最关键的因素是原子间作用势的选取,由于钛金属与氢元素的电子排布差异较大,因此其作用势较难描述。改进嵌入原子理论在氢脆相关的计算中得到较为广泛的应用,计算结果能够准确反应金属与氢元素的物理性质。本论文的研究重点即为计算钛金属与氢元素之间的改进嵌入原子势能。使用第一原理方法计算钛氢合金改进嵌入原子势能中的部分参数,计算不借助任何经验值,计算精度较好。使用局域密度近似处理密度泛函理论中的交互-关联能量方程,原子间相互作用势采用Troullier-Martins赝势。电子占据态的处理方式是将钛原子的3p轨道视为价电子。计算参数的设置为：P(?)1m(162号)空间群,240 eV的平面截断能量,0.27 eV的高斯能带展宽计算费米能,k点简化使用4×4×2的Monkhorst-Pack网格。计算模型为氢原子位于钛原子组成的八面体间隙的中心位置。计算得出钛氢合金的体积模量B,形成热ΔE和原子电子密度范围因子ρ0。其中体积模量B和形成热ΔE分别用于计算改进嵌入原子势能参数中的统一方程结构因子α和粘附能Ec,而电子密度范围因子ρ0为合金的势能参数之一。基于Hill模型,本论文计算出钛氢合金的弹性性质。根据PUGH判据,钛氢合金的体积模量与剪切模量之比GHill/Bhill=0.62>0.5,即材料呈现脆性,在内应力的作用下容易产生脆断。钛氢合金的其余势能参数应用改进嵌入原子理论进行了求导。计算得出的势能参数有截断距离,平衡距离,粘附能,统一能量方程的结构因子以及八个屏蔽因子。钛氢合金的势能确定之后,使用分子动力学方法模拟了正则系综下OK时的钛双晶晶界与间隙氢原子相互作用的力学行为,模型中晶界面的垂直方向(Y轴)采用有限边界,而其余两个方向采用周期性边界。使用Verlet速度方法和单独的时间积分共同作用实现温度控制。模拟的时间步长是10-4ps,总计进行106步的运算。其中纯元素的势能参数使用Baskes和Lee等人的计算值,而合金的势能采用本论文的计算值。模拟了氢原子位于∑7(1230)/[0001]、∑13(1120)/[0001]和∑13(1011)/[1210]晶界附近的力学行为,并且计算了钛氢合金的晶格常数,c/a比率,体积模量以及形成热。计算结果基本符合实验值和第一原理计算值,并且通过对比三个晶界受到破坏的程度,发现晶面指数最大的∑7(1230)/[0001]晶界受到的破环最为严重,其次为∑13(1120)/[0001]晶界,而晶面指数最低的∑13(1011)/[1210]晶界受到的破坏最轻。本论文对钛金属氢脆的力学行为进行了模拟,初步研究了钛氢合金的改进嵌入原子势能,使用第一原理方法计算了钛氢合金的势能参数,为钛金属氢脆机理的进一步研究提供了一定的理论基础。