论文部分内容阅读
自从Iijima发现碳纳米管以来,它已引起了众多科学家的关注和研究。已有的相关研究工作表明碳纳米管可能具有奇特的力学、电学和光学等特性,因此可望具有广阔的应用前景。然而,不论是实验还是理论工作,目前仍然处于探索阶段,且对于碳纳米管的这些特性的预测结果也仍存在着较大的分歧。因此,基于所发展的高阶Cauchy-Born广义连续体本构模型以及范德华力广义参变本构模型,本文分别研究了单壁碳纳米管的一系列相关力学特性、压缩变形机理以及平行或交叉排列的碳纳米管间的黏附行为。所取得的研究进展包括:本文提出了一种基于高阶Cauchy-Born准则建立单壁碳纳米管本构模型的方法。通过引入高阶变形梯度,合理地修正了传统Cauchy-Born准则在描述纳米管变形几何关系时所存在的缺陷。利用原子间相互作用势以及能量等效原理,得到了基于广义连续介质模型的单壁碳纳米管的本构关系。由此得到的本构参数不仅与变形梯度张量F,而且与其梯度▽F相关,同时也包含了碳纳米管的微结构信息,因此是一种考虑微结构效应的广义连续介质模型。基于所发展的高阶Cauchy-Born广义连续体本构模型,首先对石墨片和单壁碳纳米管的一系列无限小变形的力学特性(如:轴向和环向模量、泊松比、剪切模量)以及应变能进行了研究,同时也考察了手性、半径、势能参数以及内位移放松等因素对这些力学特性的影响。研究发现:当单壁碳纳米管的半径较小时,手性和半径对其力学特性的影响十分明显;而当管半径较大时,这种影响不再重要,以至于随着管半径的增大单壁碳纳米管的这些的力学特性最终趋于石墨片相应的结果。从数值结果上看,势能参数和内位移放松对所预测的单壁碳纳米管的力学特性有着很大的影响。进一步,我们也研究了简单的有限变形下单壁碳纳米管的应变能和切向模量的变化。结果表明,单壁碳纳米管的切向模量几乎线性地随着轴向拉伸而减小;且对于半径相近而手性不同的碳纳米管而言,其系统应变能几乎完全一致。继而,手性碳纳米管的轴向变形与其径向和扭转变形间的耦合机理的研究表明,最大的耦合效应发生于手性角为π/12的碳纳米管中。从变形机理方面重新审视了石墨片和碳纳米管的弯曲刚度,并指出了部分文献中关于这两个概念使用和计算上的不恰当性。具体而言,石墨片和碳纳米管的弯曲刚度这两个概念本质上是不同的,前者指的是材料的内禀特性,而后者是结构特性。当采用经典连续体理论模拟碳纳米管时,等效的连续体应该被看作是带有微结构的广义连续介质,因此本文建议把碳纳米管的结构弯曲刚度处理为一个独立的结构刚度参数。进而,基于高阶Cauchy-Born广义本构模型给出了合理的计算石墨片和碳纳米管弯曲刚度的方法和结果。基于高阶Cauchy-Born广义连续体本构模型,结合无网格方法进行了单壁碳纳米管的压缩屈曲变形的分析,从应变能的角度分析了屈曲前后碳纳米管的变形机理。所获得的与分子动力学模拟和原子有限元计算相一致的结果证明了目前所发展的模型的可行性和正确性。本文也提出了处理非成键原子间范德华力的广义参变本构模型以及基于此进行碳纳米管结构力学行为数值模拟的数学规划算法。纳米管中原子间短程力作用采用分子结构力学模型来模拟,而作为长程作用的范德华力用杆单元来模拟,这类杆单元有着特殊的非线性本构关系。对于这种高非线性问题的处理,本文基于所建立的广义参变本构模型与参数二次规划算法求解。与一般的数值方法相比较,本文的方法不需要传统的那种冗长的、反复的迭代,并具有非常好的收敛性,因此为碳纳米管结构力学行为的有效预测提供了保障。基于此我们研究了平行和交叉排列的碳纳米管的稳定构型和接触黏附行为,并发现单壁纳米管间的范德华力直接影响着其最终的稳定构型。所获得的与文献相一致的数值结果也证明了目前所发展方法的正确性和有效性。