用于C2C 共价键力能关系表征的等效C2C 键合单元3 - PSU MNE.PDF

第 28 卷 第 4 期 固 体 力 学 学 报 Vol . 28 No . 4 2007 年 12 月 A C TA M EC HAN ICA SOL IDA SIN ICA December 2007 用于 CC 共价键力能关系表征的等效 CC 键合单元 曾　攀 　杜 　婧 　杨学贵 (清华大学机械工程系 ,先进成形制造教育部重点实验室 ,北京 , 100084) σ π 摘 　要 　基于 BornOpp enheimer 近似的二次函数形式的势能关系 ,研究了 CC 键沿键轴伸缩 、在 平面和 平面内旋转的力学表征. 在原子的平衡位置附近和线弹性的范围内,得到了基于原子核位移表达的平衡关系. 在定 义一个 CC 共价键上两个原子的 6 个位移 自由度的基础上 ,构建了完整的等效 CC 键合单元. 当对构成晶形碳的 CC 共价键的六环结构进行建模时 ,等效 CC 键合单元方法能够反映出共价键中的键长 、键角和键能这三大特征. 通过对石墨片和 C60 分子振动的计算分析和实验比较 ,得到等效 CC 键单元中的力常数 ,并讨论了力常数的变化状 况. 该文还对单壁碳纳米管的拉伸弹性模量进行了计算分析 ,取得较好的效果. 等效 CC 键合单元可以为碳材料的 大规模分子力学计算提供一种实用方法. 关键词 　等效 CC 键合单元 ,分子结构表征 ,力常数 ,石墨 ,碳纳米管 关系 ,完全以结构力学的方式进行分析 ,这些方法的 0 　引言 特点是直接采用宏观力学的有限单元进行等效 ,而不 是直接从共价键的受力分析来推导相应的力能表征 , 石墨及碳纳米管中的碳原子之间是通过 CC 因而所描述的 CC 共价键的力能关系不够准确 ,甚至 共价键进行结合的[ 1 ,2 ] , 它是 自然界中已知的最强 会出现错误 ,如文献[ 16 ]在采用梁单元进行等效时 , 的化学键之一 ,单壁碳纳米管有极高的轴向强度 、韧 在原子节点上出现了像梁单元中的节点转角 ,而共价 ( ) 性和弹性模量 ,因此被称为超级纤维材料 ,在航天 、 键的转动是在键上 ,不是在原子 节点 上. 航空领域有广泛的应用前景. 因此 ,石墨及碳纳米管 本文在小位移和线性关系的基础上来直接表征 的力学性能已引起了众多学者的研究兴趣 ,有关碳 CC 共价键力能关系 ,并建立较为准确的等效 CC 纳米管的力学建模和分析也是多种多样[3 ] ,主要有 键合单元 , 具体地 , 基于 Bor nOpp enheimer 假设 , ( ) [35 ] 基于势函数的分子动力学方法 MD 、基于结合 应用分子力学方法给出 CC 共价键的伸缩能和旋 ( ) [ 67 ] 键变形能分析的分子力学方法 MM 、基于连续 转能的表达 ,然后以一个 CC 共价键的两个原子位 介质力学的等效方法[ 810 ]