脱粘和弯曲对炭纳米管复合材料弹性模量的影响.doc

脱粘和弯曲对炭纳米管复合材料弹性模量的影响 邵丽华1，2 罗瑞盈3 白树林1 王建祥1* 1北京大学工学院力学与空天技术系，北京 100871； 2中国科学院力学研究所，北京 100080； 3北京航空航天大学理学院北京 100083jxwang@pku.edu.cn 的等效模量，通过数值分析考察脱粘程度和其他因素对等效模量的影响。 模型和理论 考虑复合材料中一根长度为，半径为的纤维，如图1所示。设脱粘部分长度为，添加的纤维个数为。其中，下标表示完好粘结；下标表示部分脱粘和完全脱粘；为完好粘结的纤维个数；为部分脱粘和完全脱粘的纤维个数之和。设；。 图1 纤维与基体界面发生局部脱粘的模型 在本文中我们提出一种方法分析脱粘纤维对复合材料模量的影响。为此，将一根部分脱粘纤维分成两部分，即脱粘部分与未脱粘的部分，再将这两部分看作是独立的纤维来处理。这样，设基体为第1相，其中含有三种不同的纤维：第2相为与基体完好粘结的纤维；第3相为与基体界面发生脱粘的纤维的未脱粘部分；第4相为与基体界面发生脱粘的纤维的已脱粘部分。 根据上述对纤维的分类可知：第2、3相的弹性性质与原纤维性质相同，假设具有各向同性性质。而第4 相纤维因为已经与基体脱粘，其对基体的增强效应必然会降低，可以理解为增强相的性质发生了变化。本文将采用Zhao and Weng [5]的方法，用一与基体完好粘结并具有横观各向同性性质的等效纤维来取代已脱粘的纤维（即第4相），并且考虑脱粘纤维端部粘结的影响[6]，求得假想的横观各向同性纤维的弹性模量为（采用Walpole记号）[6] (1) 采用Mori-Tanaka方法预测复合材料的模量，得到复合材料的等效体积模量和剪切模量为（和的表达式见文献[6]） ， (3) 1.2 数值分析 设纤维的长径比为[7]。我们分析炭纳米管增强的炭基体复合材料。各向同性炭基体的杨氏模量取为[8]；炭纳米管的杨氏模量取为[9]，泊松比假设为。 图2为纤维体积分数时随纤维脱粘程度的变化情况。由曲线、的变化规律可知，复合材料的等效体积模量随纤维脱粘程度的增加而减小；由曲线可见，当无脱粘纤维时材料体积模量不受影响。复合材料剪切模量也具有相似变化规律。 图2 时随纤维脱粘程度的变化规律 二、纤维弯曲的影响 许多实验观察表明，复合材料中的炭纳米管大多是以弯曲的形式存在的[3]。所以有很多学者认为，炭纳米管的弯曲是导致其实际增强性能比理论值低很多的原因之一[1-3]。Fisher等[10]将弯曲炭纳米管等效为一具有折减模量的直纤维，采用有限元方法计算了它的等效模量，然后用Mori-Tanaka方法预测弯曲炭纳米管增强的聚合物基体复合材料的等效模量，发现纤维的弯曲确实会使复合材料的等效弹性模量降低很多。Shi等人[2]的研究中将弯曲的炭纳米管看作（空间）螺旋线。而根据很多实际观察 ，炭纳米管的弯曲程度没有螺旋线那么大。本节将建立一简单模型，预测纤维弯曲对复合材料等效弹性模量的影响。同时预测当纤维具有不同的几何形状时对复合材料等效模量的影响。 2.1模型和理论分析 从很多文献中可见，炭纳米管在复合材料中一般呈弓状弯曲。Chou and Takahashi [11]和Kuo等人[12]在研究波浪形状的连续纤维增强复合材料的等效模量时，假设弯曲纤维具有正弦曲线形状，并推导出由该弯曲纤维增强的复合材料的横向和纵向等效模量。本节研究含有在基体中随机取向分布的弯曲炭纳米管复合材料的等效弹性模量。如图3所示，将纤维向其水平和竖直方向投影。设水平方向投影所得纤维长度为，直径为，纤维体积分数为；竖直方向投影所得纤维长度为，直径为，纤维体积分数为。 图3 弯曲纤维模型 设所研究材料的体积为，则有 ， (4) 设，且，则由式(4)可求得，。 设图3中水平方向投影的长为的纤维等效模量为，由该单向纤维增强的复合材料的等效模量为。由混合律可知，则 (5) 其中，为基体的弹性模量，、分别为基体和长度为的等效纤维的体积分数。我们引用参考文献[11]的结果计算从而求得。 同理，设图3中竖直方向投影的长为的纤维等效模量为，由该单向纤维增强的复合材料的等效模量为。由混合律可知 (7) 其中， 分别为长度为的等效纤维的体积分数。我们引用参考文献[12]的结果计算从而求得。 此时，可认为复合材料中的