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

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

论文摘要

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

论文目录

  • 摘要
  • ABSTRACT
  • 目录
  • 插图
  • 第1章 绪论
  • 1.1 引言
  • 1.2 碳纳米管的结构
  • 1.3 碳纳米管的力学特性与应用
  • 1.4 碳纳米管力学性能的实验检测
  • 1.5 碳纳米管力学性能的理论研究及模拟
  • 1.5.1 理论预测方法
  • 1.5.2 数值模拟方法
  • 1.6 本文的研究目标和内容
  • 第2章 分子模拟方法的基本介绍
  • 2.1 分子模拟的概况
  • 2.2 力场
  • 2.2.1 力场简述
  • 2.2.2 力场作用项的一般式
  • 2.2.3 常见的力场
  • 2.3 分子力学方法
  • 2.3.1 基本原理
  • 2.3.2 分子稳定结构的确定
  • 2.4 分子动力学模拟的原理
  • 2.4.1 分子动力学发展史
  • 2.4.2 分子动力学求解基本步骤
  • 2.4.3 作用势与动力学计算
  • 2.4.4 时间步长与约束动力学
  • 2.4.5 分子动力学模拟在碳纳米管力学行为模拟中的应用
  • 2.5 特卡罗计算方法
  • 2.6 本章小结
  • 第3章 分子结构力学方法
  • 3.1 分子结构力学方法
  • 3.1.1 分子结构力学方法基本原理
  • 3.1.2 分子结构力学方法进展
  • 3.2 数值算例
  • 3.2.1 有限元模型及参数选取
  • 3.2.2 弹性模量
  • 3.2.3 动力学特性
  • 3.3 本章小结
  • 第4章 修正的分子结构力学方法
  • 4.1 修正的分子结构力学方法简述
  • 4.2 修正的分子结构力学的基本思想
  • 4.3 系统能量的等效离散处理
  • 4.4 内禀坐标与原子位移间的关系
  • 4.4.1 符号规定
  • 4.4.2 键长的伸缩变化
  • 4.4.3 键角的弯曲变化
  • 4.4.4 离面角的表示
  • 4.5 修正的分子结构力学方法的建立
  • 4.6 应用修正的分子结构力学方法研究单壁碳纳米管的静力学特性
  • 4.6.1 纳米管的均匀拉伸
  • 4.6.2 碳纳米管的弯曲
  • 4.6.3 碳纳米管的扭转
  • 4.7 应用修正的分子结构力学方法研究单壁碳纳米管动力学特性
  • 4.7.1 多自由度系统的振动
  • 4.7.2 单壁碳纳米管的固有频率计算
  • 4.7.3 连续模型的固有频率
  • 4.8 本章小结
  • 第5章 化学键单元方法
  • 5.1 化学键单元方法基本原理及方程
  • 5.1.1 平面刚度矩阵的推导
  • 5.1.2 原子抗扭刚度与键角弯曲弹力常数的关系
  • 5.1.3 空间刚度矩阵的推导
  • 5.2 化学键单元方法的特点
  • 5.2.1 与宏观梁模型的不同
  • 5.2.2 化学键单元方法的优势
  • 5.3 工程软件中化学键单元方法的实现
  • 5.3.1 ANSYS软件基本单元介绍
  • 5.3.2 单元刚度矩阵“病态”的消除
  • 5.3.3 单元刚度矩阵从局部坐标到整体坐标的转换
  • 5.3.4 转换矩阵可信度的检验
  • 5.4 碳纳米管基本力学特性的数值算例
  • 5.4.1 碳纳米管的均匀拉伸
  • 5.4.2 碳纳米管的扭转
  • 5.4.3 力场弹力常数选取对弹性模量计算结果的影响
  • 5.4.4 有效厚度选取对弹性模量计算结果的影响
  • 5.4.5 泊松比
  • 5.4.6 碳纳米管的动力学特性
  • 5.5 本章小结
  • 第6章 纳米秤及多壁碳纳米管的数值模拟
  • 6.1 纳米秤模拟
  • 6.1.1 纳米秤的制造
  • 6.1.2 纳米秤的理论及数值模拟研究
  • 6.1.3 化学键单元方法在纳米秤数值模拟中的应用
  • 6.2 单壁碳纳米管等效厚度的推导
  • 6.3 多壁碳纳米管力学行为的数值模拟
  • 6.3.1 范德华力
  • 6.3.2 双壁碳纳米管的杨氏模量及剪切模量
  • 6.3.3 双壁碳纳米管的振动
  • 6.4 本章小结
  • 第7章 全文总结和展望
  • 7.1 工作总结
  • 7.2 工作展望
  • 在读期间发表的学术论文与取得的其他研究成果
  • 致谢
  • 相关论文文献

    标签:;  ;  ;  ;  ;  ;  

    碳纳米管力学行为的数值模拟
    下载Doc文档

    猜你喜欢