碳纳米管力学行为数值模拟.pdf

摘要 纳米材料被誉为是21世纪的重要材料，并将构成未来智能社会的四大支柱 之一。碳纳米管在纳米材料中最富有代表性，并且是性能最优异的材料之一。 自从日本科学家首次发现碳纳米管，各国科学家对其进行了大量的研究。美、 英、法、德、日及我国均相继成立了纳米材料研究中心，其重要的研究内容就 是碳纳米管。1997年，单壁碳纳米管的研究成果与“克隆羊”和“火星探路者’’ 一起，分别被评为当年的世界十大科学成就之一。碳纳米管具有优秀的力学、 物理、化学及电学性质，在很多方面都有重要的应用。但是，目前在工程应用 方面尚未取得完全的突破性进展。为此，研究人员必须对碳纳米管的结构及其 性能进行全面的研究，这些研究具有重要的理论和工程应用价值。 本文围绕碳纳米管的力学行为，开展了包括建立模型及数值模拟方面的研 究，具有一定的理论意义和实用价值。主要开展了分子结构力学方法、修正的 分子结构力学方法和化学键单元方法在碳纳米管力学行为模拟方面的研究，具 体内容如下。 1．分子结构力学方法在碳纳米管力学行为方面的数值模拟研究 分子结构力学方法是一种基于分子力学的数值计算方法。这种方法的特点 是将原子间的共价键等效成宏观结构的梁，从而，将碳纳米管的微观结构等效 为类似空间刚架的结构，然后用结构力学的方法求解。 迄今为止，学者们利用该方法开展的工作主要集中在理论上，但应用该方 法对复杂的微纳米机电器件进行模拟时会出现一定困难。本文采用有限元方法 求解结构力学问题，使用工程软件实现碳纳米管力学行为的有限元模拟，可以 避免出现这种困难。数值模拟的结果表明，单壁碳纳米管的杨氏模量、剪切模 量与其直径之间存在依赖现象，即碳纳米管的杨氏模量和剪切模量会随着管径 的变化而变化。此外，对单壁碳纳米管系统，在不考虑阻尼时，振动基频与碳 纳米管长度直径比之间也存在尺度依赖关系。最后，分子力学力场弹力常数对 单壁碳纳米管力学行为的影响也在本部分中也做了详细的讨论。 2．修正的分子结构力学方法在碳纳米管力学行为方面的研究 经过仔细地研究分子力学方法在碳纳米管中的应用，本文作者所在研究组 提出了一种基于分子力学的数值计算方法：修正的分子结构力学方法。在工作 组前面工作的基础上，本文将该方法进一步拓展应用于碳纳米管的剪切模量及 动态力学行为的分析。修正的分子结构力学方法不再是简单地将碳纳米管的原 子结构类比成宏观的某些结构，而是用分子力学的力场势能函数表述碳纳米管 系统的势能，这样的处理避免了简单的能量等效可能会产生的误差。本方法将 系统能量按一定的方式离散，从能量最小原理出发，在小变形假设的基础上， 直接建立系统方程。由于选择的力场包含了离面振动能量项，所以该方法有更 广的适用范围及更好的精度。本部分系统地介绍了该方法从理论建模到程序实 现的过程，在课题组前期工作基础上，系统地计算研究了单壁碳纳米管杨氏模 量、剪切模量及振动基频与碳纳米管长度／直径之间的关系。 3．化学键单元方法在碳纳米管力学行为方面的研究 本文建立了一种基于分子力学的三维纳米尺度的有限元方法：化学键单元 方法。化学键单元模拟了碳纳米管原子间碳．碳化学键的力学行为，单元的刚度 矩阵通过联系分子力学与连续介质力学而得到。化学键单元的元素全部是分子 力学力场弹力常数的函数。本方法避免了简单的能量等效方法带来的计算误差， 运算方便，效率更高。应用工程有限元软件提供的扩展单元建立化学键单元， 可以实现碳纳米管力学行为的数值模拟。单壁碳纳米管杨氏模量、剪切模量及 振动固有频率与碳纳米管长度／直径之间的关系在本部分论文中做了详细的讨 论。此外，论文中尝试应用该方法研究了简单的微纳米器件(纳米秤)的力学 行为。最简单的纳米秤结构由一根悬臂的碳纳米管及附着于纳米管的微粒组成。 本部分论文针对如何通过纳米秤系统共振频率的改变得到附加微粒质量的理论 及数值模拟方法进行了讨论。关于单壁碳纳米管直径、长度以及微粒附着位置 对纳米秤系统共振频率的影响及如何有效提高纳米秤敏感度的办法，也通过数 值模拟开展了研究。最后，通过参数拟合的办法对单壁碳纳米管的有效厚度进 行了研究，并采用非线性弹簧模拟了范德华力，对双壁碳纳米管的杨氏模量、 剪切模量以及固定外壁．内壁自由的一阶往复振动进行了分析。 关键词：碳纳米管，分子力学，分子结构力学，修正的分子结构力学(MMSMM)， 化学键单元方法，有限单元法(FEM) ABSTRACT NaIlo—materiaIsarec