基于原子势的碳纳米管有限元模拟

【摘要】： 碳纳米管具有非常优异的力学、热学、电学和光学性能,应用前景广阔。自从碳纳米管被发现以来,对其力学性能的理论预测一直是一个热门研究领域。本文基于原子势建立了碳纳米管的有限元模型,基于该模型开展了碳纳米管的若干关键力学问题研究,主要包括以下内容: (1)碳纳米管有限元模型构建。提出一种新的能够完备地替代谐和势描述C-C共价键的等效梁单元,该单元能够区分C-C键的σ-σ键角和σ-π键角的不同弯曲刚度。通过分析石墨层的本构关系确定了等效梁单元的参数。基于该等效梁单元建立了单壁碳纳米管的有限元模型,并用非线性弹簧等效Lennard-Jones势描述多壁碳纳米管层间范德华作用从而建立了多壁碳纳米管的有限元模型。本文模型能够直接在标准商业有限元软件中实现,从而可以利用商业有限元软件的强大功能开展碳纳米管力学问题的模拟。 (2)碳纳米管弹性模量计算。将单壁碳纳米管看作横观各向同性材料,系统地计算了其独立的5个弹性模量随直径和手性的变化。发现其直径越小,弹性模量对直径和手性的依赖性越强,且各向异性越明显;随着直径的增大,弹性模量的各向异性逐渐消失,当直径足够大时候,碳纳米管退化成各向同性材料,弹性模量趋于石墨层的模量。这些结果拓展了前人的文献报道,提供了一个对单壁碳纳米管各向异性性质的完整理解。 (3)碳纳米管屈曲行为模拟。应用特征值屈曲计算方法系统地分析了单壁和多壁碳纳米管临界压缩屈曲模态和临界压缩屈曲应变随其长径比的变化,发现相同长度的单壁碳纳米管承受临界压缩屈曲应变的能力有一个最优直径,处在该直径的单壁碳纳米管承受压缩屈曲应变的能力最强。还用一种附加阻尼的方法对单壁碳纳米管的弯曲屈曲行为进行了初步探讨,结果为相关分子动力学计算所验证。 (4)碳纳米管振动模态分析。应用特征值振动方法计算了单壁和多壁碳纳米管的各种振动模态,着重分析了多壁碳纳米管的“刚性”同轴振动模态(仅由范德华力决定的模态)。本文计算范围内碳纳米管的振动频率均在GHz以上,表明碳纳米管有望成为设计高频纳米器件的理想材料。 (5)碳纳米管拉曼活性模态计算。用特征值振动方法计算了单壁碳纳米管的拉曼模态,计算表明低阶拉曼模态与相关实验和理论报道相当吻合,尤其径向呼吸模态与前人第一原理计算结果高度一致,说明本文模型在拉曼模态计算方面具有一定的适用性。本文还对多壁碳纳米管的径向呼吸模态进行了详细预测,结果表明多壁碳纳米管存在与其层数相同阶数的径向呼吸模态,且其振动频率分别高于各层单壁碳纳米管径向呼吸模态的频率。 (6) Casimir力和静电力驱动下的GHz振动探讨。从多壁碳纳米管得到启发,抽象出双层同轴圆柱管模型,探讨了该模型在Casimir力和静电力作用下的振动特性,结果显示在纳米尺度该模型振动频率即可达到GHz以上,表明在纳米尺度GHz高频振动的普遍存在性。