表面效应对碳纳米管中弯曲波波动特性的影响

第51卷第8期2020年8月
中南大学学报(自然科学版)
Journal of Central South University(Science and Technology)
V ol.51No.8
Aug.2020表面效应对碳纳米管中弯曲波波动特性的影响
黄彬1，武井祥1，金花2，周强2
(1.湘潭大学土木工程与力学学院，湖南湘潭，411105；
2.湘潭大学物理与光电工程学院，湖南湘潭，411105)
摘要：基于广义梯度弹性梁理论研究表面效应对自由空间和弹性介质中碳纳米管(CNTs)波动性能的影响。

碳纳米管采用同时考虑弯曲变形和剪切变形的剪切梁进行描述，弹性介质采用双参数Pasternak-type弹性基模拟。

建立考虑表面效应的广义梯度弹性剪切梁控制方程，推导碳纳米管中弯曲波的色散关系式，并通过分子动力学模拟结果进行验证。

探讨表面效应、尺度因子、弹性介质对碳纳米管中弯曲波相速度的影响。

对于多壁碳纳米管(MWCNTs)，研究波速与范德华力的相关性，讨论MWCNTs层数、表面效应、尺度因子和弹性参数对MWCNTs相速度的影响。

研究结果表明：由所推导的色散关系式所得理论结果与分子动力学模拟结果较吻合，表现所推导的色散关系式理论模型能较好地表征碳纳米管弯曲特性。

关键词：碳纳米管；表面效应；尺度因子；弹性介质；广义梯度弹性理论
中图分类号：TB121文献标志码：A
文章编号：1672-7207（2020）08-2289-10
Influence of surface effect on wave behavior of flexural waves in
carbon nanotubes
HUANG Bin1,WU Jingxiang1,JIN Hua2,ZHOU Qiang2
(1.School of Civil Engineering and Mechanics,Xiangtan University,Xiangtan411105,China;
2.School of Physics and Optoelectronic Engineering,Xiangtan University,Xiangtan411105,China)
Abstract:Based on the theory of generalized gradient elastic beam,the influence of surface effect on the wave behavior of carbon nanotubes(CNTs)in free space and elastic medium was Ts were described by shear beam that considered both bending and shear deformation.The elastic medium was simulated with a two-parameter Pasternak-type elastic foundation.The governing equation of generalized gradient elastic shear beam considering the surface effect was established.The dispersion relation of bending wave in CNTs was derived, which was verified by the results of molecular dynamics simulation.The influence of surface effects,scale factors and elastic medium on the phase velocity of bending waves in CNTs was investigated.For multi-walled carbon nanotubes(MWCNTs),the dependence of wave velocity on van der Weals forces was studied,and the effects of MWCNTs layer number,surface effect,scale factor and elastic parameters on phase velocity of MWCNTs were DOI:10.11817/j.issn.1672-7207.2020.08.023
收稿日期：2019−11−04；修回日期：2020−01−23
基金项目(Foundation item)：湖南省教育厅−湘潭大学配套科研项目(KZ03030)(Project(KZ03030)supported by the Scientific Research Program of Department of Education of Hunan Province−Xiangtan University)
通信作者：武井祥，博士，讲师，从事微纳米力学研究；E-mail:*******************.cn
第51卷
中南大学学报(自然科学版)
discussed.The results show that the theoretical rusults obtained by the diperson reletion of bending wave in CNTs
are agreement with those by molculer dynamies simulation,which incicates that the theoretical model can characterize the wave behavior of flexural waves in carbon nanotubes.
Key words:carbon nanotubes;surface effect;scale factor;elastic medium;generalized gradient elastic beam theory
1991年，IIJIMA[1]研制了碳纳米管，掀起了人们对碳纳米管的研究热潮。

碳纳米管具有很多优异性能，其中，其具有良好的柔韧性，在工业上常用作增强型纤维，是一种理想的高强度纤维材料，在工业生产中，通常以其他工程材料为基体与碳纳米管制成复合材料，以改善复合材料的强度及韧性等性能。

连续介质力学理论常用于研究碳纳米管力学性能，但其模型不涉及尺寸依赖性。

而大量实验结果表明，碳纳米管的力学性能具有较强的尺寸依赖性。

在此情况下，ERINGEN[2]提出了非局部弹性理论作为研究纳米材料尺度效应的理论方法。

基于此理论，LIEW等[3]结合非局部Timoshenko梁理论与分子动力学这2种方法探讨了碳纳米管中弯曲波的传播特性。

PANG等[4]结合非局部弹性理论探讨了尺度效应下黏弹性单壁碳纳米管中横波的传播特性，发现表面尺度效应对碳纳米管中横波的特性影响取决于波数与管的直径这2个因素。

余阳等[5]基于非局部应变梯度Euler梁模型，研究了尺度效应下充流碳纳米管的波动情况，发现应变梯度对低频波动起促进作用，对高频波动起阻尼作用。

王碧蓉等[6]探讨了非局部效应对碳纳米管中弯曲波频散特性的影响，发现非局部效应对弯曲波频散特性的影响主要体现在高波数阶段。

SHEN等[7]基于非局部Timoshenko梁模型研究了被弹性介质包围时碳纳米管中弯曲波的特性，发现弹性介质对低频波速有重要影响。

GAFOUR等[8]基于非局部Euler-Bernoulli梁理论，探讨了弹性介质中“之”字形双壁碳纳米管中弯曲波的特性。

碳纳米管具有较大的比表面积，研究其性能时也需要考虑表面效应的影响。

CAMMARATA[9]研究了界面材料与纳米结构材料的表面效应。

LI等[10]基于非局部应变梯度理论，研究了表面效应对黏弹性单壁碳纳米管中波特性的影响。

ZHEN[11]分析了表面效应对黏弹性单壁碳纳米管中波特性的影响，发现表面效应的影响在小波数及小管径下显得尤为显著。

NARENDAR等[12]探讨了表面效应下纳米管的非局域波特性，发现考虑表面效应的弯曲波数比不考虑表面效应的弯曲波数多，并且考虑表面效应时，弯曲波表现出压缩性质。

综上可知，考虑表面效应的Euler-Bernoulli梁模型未考虑转动惯量和剪切变形的影响。

而考虑表面效应的Timoshenko梁模型，其表面效应模型中曲率采用∂2ω/∂x2，显然忽略了剪切变形的影响。

人们对基于非局部弹性梁理论下碳纳米管的波动性能进行了大量研究，但大多只考虑了应力梯度的影响，未考虑应变梯度的影响。

为了更准确地表征表面效应下碳纳米管中弯曲波的波动特性，本文作者建立考虑表面效应的广义梯度剪切梁模型。

与其他模型不同，表面效应理论模型中曲率采用-∂φ/∂x的形式，并且同时兼顾了弯曲变形和剪切变形与表面效应的相关性。

此外，广义梯度综合考虑了应力和应变双梯度的影响。

与碳纳米管极大的弹性模量相比，其周遭介质可被视为弹性介质，弹性介质采用双参数Pasternak-type弹性基描述。

与Winkler-type基相比，Pasternak-type弹性基同时考虑了轴向和径向的弹簧刚度。

为此，推导碳纳米管中弯曲波的色散关系式，并通过分子动力学模拟对理论结果进行验证，分别探讨表面效应、尺度因子、弹性介质对单壁碳纳米管和多壁碳纳米管中弯曲波相速度的影响。

1考虑表面效应的广义梯度剪切梁模型
基于广义梯度弹性梁理论，本构方程可以表述为[13]
(1-l2
1∇2)σij=(1-l22∇2)(λδijεkk+2μεij)(1)
式中：σ
ij
和ε
ij
分别为应力张量与应变张量；λ和μ
均为Lamé常数；l
1
=e
1
a，l
2
=e
2
a，均为尺度系
数；a为C—C键长度；e
1
和e
2
均为量纲一材料常
2290
第8期黄彬，等：表面效应对碳纳米管中弯曲波波动特性的影响数，其值可通过实验或分子动力学模拟得到，∇2
为拉普拉斯算符。

当l
1
=0时，上述理论变为应变
梯度弹性理论；当l
2
=0时，上述理论变为非局部
弹性理论。

对于一维纳米结构，其应力张量σ
yy =σ
zz
=
σyz =σ
xy
=0，此时，式(1)简化为
(1-l21∂2∂x2)σxx=(1-l22∂2∂x2)Eεxx(2) (1-l21∂2∂x2)σxz=(1-l22∂2∂x2)Gγxz(3)
式中：E和G分别为弹性模量与剪切模量；x为纵
坐标；z为垂直于未变形梁中性轴的坐标；σ
xx 和ε
xx
分别为正应力与正应变；σ
xz 和γ
xz
分别为剪应力与
剪应变。

对于剪切梁模型，正应变ε
xx 与剪应变γ
xz
可以
表示为
ε
xx
=z∂ϕ∂x(4)
γ
xz
=ϕ+∂w∂x(5)其中：w为横向挠度，记作w(x,t)；ϕ为横截面的
转角，记作ϕ(x,t)。

将式(4)与(5)代入式(2)与(3)，得
(1-l21∂2∂x2)σxx=E(1-l22∂2∂x2)z∂ϕ∂x(6) (1-l21∂2∂x2)σxz=G(1-l22∂2∂x2)(ϕ+∂w∂x)(7)
对方程(6)进行变换处理，在方程两边同时乘
以z，再对其在横截面上积分，得
(1-l21∂2∂x2)M=EI(1-l22∂2∂x2)∂ϕ∂x(8)其中：M为弯矩，M=∫A zσxx d A；I为横截面惯性矩，I=∫A z2d A；A为横截面面积。

对方程(7)两边积分，得
(1-l21∂2∂x2)Q=G s A(1-l22∂2∂x2)(ϕ+∂w∂x)(9)其中：Q为剪力，Q=∫Aσxz d A；G s为有效剪切模
量，G
s
=κG；κ为剪切修系数，由所选横截面形状决定。

基于剪切梁理论，运动方程可以写作
∂M
∂x=Q(10)
∂Q
∂x+q=ρA
∂2w
∂t2(11)其中：ρ为碳纳米管密度；q为作用于横向表面上的分布载荷。

对式(10)与式(11)进行变形处理，得
∂2M
∂x2=
∂Q
∂x(12)
∂2Q
∂x2=ρA
∂3w
∂t2∂x-
∂q
∂x(13)当考虑表面效应时，有如下新增变量：
(EI)*=EI+E s J(14)
P=-H∂ϕ∂x(15)
其中：J=πR3；H=4τ
R；E s为表面弹性模量；τ
0为残余表面张力；R为碳纳米管横截面半径；()
EI*为梁体和表面等效刚度。

将式(12)和式(13)代入式(8)和(9)，得
M=l2
1(ρA
∂2w
∂t2-q
)+EI(1-l22∂2∂x2)∂ϕ∂x(16)
Q=l2
1(ρA
∂3w
∂t2∂x-
∂q
∂x
)+G s A(1-l22∂2∂x2)(ϕ+∂w∂x)
(17)
将式(14)，(15)，(16)和(17)代入式(10)和(11)，可得到考虑表面效应的广义梯度剪切梁控制方程：(EI)*(1-l22∂2∂x2)∂2ϕ∂x2=G s A(1-l22∂2∂x2)(ϕ+∂w∂x)
(18)
G
s
A∂∂x(1-l22∂2∂x2)(ϕ+∂w∂x)=
()
1-l2
1∂
2
∂x2
()
ρA∂2w∂t2+H∂ϕ∂x-q(19)上述控制方程是关于转角ϕ和挠度w的耦合方程。

为了便于下一步计算，用1个新的函数F进行简化。

在此情况下，转角ϕ和挠度w为
ϕ=-(1-l22∂2∂x2)∂F∂x(20) w=(1-l22∂2∂x2)F-()EI*G s A(1-l22∂2∂x2)∂2F∂x2(21)
因此，考虑表面效应的广义梯度剪切梁控制方程可以用函数F为
(1-l22∂2∂x2)2∂4F∂x4+
2291
第51卷
中南大学学报(自然科学版)(
)1-l
21
∂2
∂x 2
(
)
1-l 2
2
∂2∂x 2∂2∂t 2éëêêù
û
úúρA ()EI *F -ρG s ∂2F ∂x 2-H
()
EI *
(
)1-l
21
∂2
∂x 2
(
)
1-l 22
∂2∂x 2∂2F ∂x 2
=
1
()
EI *
(
)
1-l 21
∂2∂x 2
q
(22)
2单壁碳纳米管中的弯曲波
在研究过程中，碳纳米管通常被用作增强材
料而嵌于复合材料中，相比碳纳米管极大的弹性模量，周围介质可视为弹性介质。

本研究中用Pasternak-type 弹性基模型模拟，其表达式如下：
q =-β0w +β1a
2
∂2w
∂x 2
(23)
式中：β0与β1分别为径向约束的弹簧刚度与切向约束的弹簧刚度。

当弯曲波在碳纳米管中传播时，函数F 的解析式如下：
F (x ,t )=-F e i ()
kx -ωt (24)
式中：-F 为振幅；i =
-1；k 为波数；ω为圆频
率。

此时，式(23)可变化为
q =-(β0+β1a 2k 2)(1+l 22k 2)-F e i ()
kx -ωt (25)
将式(24)和(25)代入式(22)，得
b 0ω4-2b 1ω2+b 2=0
(26)
其中：b 0=0;b 1=12éëêêρA ()EI *+ρG S k 2ù
û
úú(1+l 21k 2)(1+l 22k 2)；b 2=
1
()
EI *(β0
+β1a 2k 2+Hk 2)(1+l 21k 2)(1+l 22k 2)+
(1+l 22k 2)
2
k 4。

由于b 0为0，式(26)可以被简化成2b 1ω2-b 2=0，并可求得ω为
ω
=
(27)
进一步可得到弯曲波在碳纳米管中传播时的相速度c
=
b 1及b 2，可得到单壁碳纳米管中传播的弯曲波的相速度
c ：
c =
{G s ()β0+β1a 2k 2+Hk 2[
]G s A +()EI *
k
2
ρk
2
+
}G s ()EI *
k 2
[]
G s
A +()EI *
k 2
ρ1+l 22k 2
1+l 21k
2
1/2
(28)
2.1考虑表面效应的非局部剪切梁理论（NSBT ）
基于非局部剪切梁理论(l 2=0)，弯曲波的相速度变为
c =
{G s ()β0+β1a 2k 2+Hk 2[
]G s A +()EI *
k
2
ρk
2
+
}
G s ()EI *
k
2[]
G s
A +()EI *
k 2
ρ11+l 21k
2
1/2
(29)
2.2考虑表面效应的应变梯度剪切梁理论
（SGSBT ）
令l 1=0，可以得到基于应变梯度剪切梁理论，弯曲波的相速度为
c =
{G s ()β0+β1a 2k 2+Hk 2[
]G s A +()EI *
k
2
ρk
2
+
}
G s ()EI *
k 2
[]G s
A +()EI *
k 2
ρ
()
1+l 22k 21/2
(30)
2.3经典剪切梁理论（CSBT ）
当不考虑尺度参数和弹性介质的影响即碳纳米管处于自由空间，且l 1=l 2=0时，可得到经典剪切梁理论下弯曲波在单壁碳纳米管中传播的相速度为c =
(31)
2.4单壁碳纳米管的结果分析与讨论
本研究中，选择(5，5)扶手型单壁碳纳米管作为研究对象，得出弯曲波基于不同梁模型下在单壁碳纳米管中传播时的色散关系图，并与分子动
力学模拟结果进行对比，结果如图1所示。

密度ρ=2.237g/cm 3，弹性模量E =0.39TPa ，泊松比ν=0.28，剪切模量G =0.5E /(1+ν)，厚度t =0.342nm [14]；尺度因子l 1=0.8a ，l 2=0.1a 。

对于
横截面为圆形的碳纳米管，剪切修正系数κ=0.8[15]。

表面弹性模量E S =5.1882N/m ，表面
2292
第8期
黄彬，等：表面效应对碳纳米管中弯曲波波动特性的影响残余张力τ0=0.9108N/m [16−17]。

图1中，GSBT ，GEBT 和MD 分别代表考虑表面效应的广义梯度剪切梁理论、考虑表面效应的广义梯度Euler-Bernoulli 梁理论和分子动力学模拟结果。

从图1可见：GEBT 预测的波速只在小波数区间内与MD 预测的波速较接近；随着波数增大，与MD 模拟结果偏差越大；而GSBT 预测的波速与MD 预测的波速[14]在整个波数区间都较吻合。

因此，考虑表面效应的广义梯度剪切梁模型能很好地表征CNTs 的波动特性。

图2所示为考虑与不考虑表面效应这2种情况下单壁碳纳米管中弯曲波相速度的对比，以便考察表面效应对CNTs 波速的影响。

从图2可知：在小波数范围内，2条相速度曲线存在明显偏差；随着波数增大，这种偏差逐渐减小。

这说明表面效应对单壁碳纳米管在小波数区间内的弯曲波特性影响显著，而随着波数增大，这种影响逐渐减小，直至消失。

图3所示为基于不同剪切梁理论下，单壁碳纳米管中弯曲波的相速度曲线，GSBT ，NSBT ，SGSBT 和CSBT 分别代表考虑表面效应的广义梯度剪切梁理论、考虑表面效应的非局部剪切梁理论、考虑表面效应的应变梯度剪切梁理论和经典剪切梁理论。

由图3可知：在小波数范围内，弯曲波在4种不同剪切梁理论下的相速度曲线非常接近；随着
波数增大，SGSBT 模型的相速度曲线、NSBT 模型的相速度曲线以及CSBT 模型的相速度曲线都逐渐偏离GSBT 模型的相速度曲线；当波数足够大时，NSBT 模型的相速度曲线无限趋近于0；SGSBT 模型的相速度和CSBT 模型的相速度则随着波数增加
而增加。

SGSBT ，CSBT 和GSBT 这3个模型不符合“弹性波可认为是长波的极限”这一事实。

因此，当考虑尺度系数的影响时，选择广义梯度剪切梁理论研究单壁碳纳米管中弯曲波的特性是合
理的。

图4和图5所示分别为尺度因子e 1和e 2
对弯曲
图1
不同梁模型下弯曲波相速度与分子动力学结果的
比较
Fig.1
Comparison of phase velocity and molecular
dynamics results of flexural wave in different beam
models
图2
考虑表面效应与不考虑表面效应情况下单壁
碳纳米管的弯曲波相速度曲线
Fig.2
Phase velocity curves of flexural waves in SWCNTs with or without surface
effect
图3不同剪切梁理论下单壁碳纳米管中弯曲波的
相速度
Fig.3
Phase velocities of bending wave in SWCNTs with
different shear beam theories
2293
第51卷
中南大学学报(自然科学版)波相速度的影响。

分析图4和图5可知：在小波数
范围内，尺度因子e 1和e 2对弯曲波相速度的影响可以忽略；随着波数增加，e 1和e 2对弯曲波相速度的影响逐渐增强；不同e 1和e 2对相速度的影响程度也不同，对于大波数范围，e 1增大使相速度减小，而在波数较大时，相速度随e 2增大而增大。

选取参数径向约束的弹簧刚度β0与切向约束的弹簧刚度β1，弹性介质对弯曲波相速度的影响如图6和图7所示。

从图6和图7可见：弹性参数β0和β1对CNT 中弯曲波相速度有较大影响；在小波数范围内，随着β0和β1增大，对应的相速度均呈现增大趋势，
说明弹性介质的存在使弯曲波相速度的预测结果
偏大；波数越小，弹性介质对单壁碳纳米管中弯曲波相速度的影响越显著；当波数足够大时，弹性介质对弯曲波相速度的影响可以忽略，也就是说，波数越小，单壁碳纳米管中弯曲波的相速度对弹性介质越敏感。

从图6和图7还可看出：β1对弯曲波相速度的影响范围比β0的影响范围要广，这是因为β0与β1分别为径向约束的弹簧刚度与切向约束的弹簧刚度，所以，可以认为弹性介质的切应力对波数的影响范围要比正应力对波数的影
响范围更广。

图4
尺度因子e 1对单壁碳纳米管中弯曲波相速度
的影响(e 2=0.1)
Fig.4Influence of scale factor e 1on phase velocity of
flexural wave in SWCNTs when e 2
=0.1
图5
尺度因子e 2对单壁碳纳米管中弯曲波相速度
的影响(e 1=0.8)
Fig.5Influence of scale factor e 2on phase velocity of
flexural wave in SWCNTs when e 1
=0.8
图6弹簧刚度β0对单壁碳纳米管中弯曲波相速度
的影响(β1=0)
Fig.6
Influence of spring stiffness β0on phase velocity of flexural wave in SWCNTs when β1
=0
图7弹簧刚度β1对单壁碳纳米管中弯曲波相速度
的影响(β0=0)
Fig.7
Influence of spring stiffness β1on phase velocity of flexural wave in SWCNTs when β0=0
2294
第8期
黄彬，等：表面效应对碳纳米管中弯曲波波动特性的影响
3
多壁碳纳米管中的弯曲波
3.1
多壁碳纳米管中相速度的推导
不同于单壁碳纳米管，多壁碳纳米管相邻管
之间存在范德华力，且范德华力可由Lennard-Jones 势推导得到[18−21]。

第j 层碳纳米管的范德华力q j 为q j =-αj (w j -w j +1)-αj -1(w j -w j -1)；
j =1,2,⋯,N -1
(32)
其中：αj 为范德华作用系数(单位为10−7J)，
αj ={
0;j =0
320(2R j )
0.16a 2
;j =1,2,⋯,N -1
(33)
R j 为相应内管的半径(单位为nm)。

对于最外层碳纳米管，范德华力的表达式可以写成
q N =-αN -1(w N -w N -1)-β0w N +β1a 2
∂2w N ∂x 2
(34)
最外层碳纳米管会同时受到范德华力与周围介质的影响。

当弯曲波在多壁碳纳米管中传播时，函数F 的表达式为
F j (x ,t )=-F j e i ()
kx -ωt (35)
其中：j 为相应碳纳米管的层数；-F j 为振幅；k 为
波数；ω为圆频率。

将方程(35)代入式(32)与(34)得
q 1=-α1(-F 1--F 2)(1+l 22k 2)e i ()kx -ωt (36)
q j =-[]
()αj +αj -1-F j -αj -1-F j -1-αj -F j +1⋅
()1+l 22k 2e i ()kx -ωt ;j =2,3,⋯,N
-1(37)
q N =[αN -1-F N -1-
]
()α
N -1
+β0+β1a 2k 2-F N ()1+l 22k 2e i ()kx -ωt (38)
再将式(36)，(37)和(38)代入式(22)，得
éë
êêêù
û
úú
ú1+l 22
k 2
1+l 21
k 2k 4+α1()EI 1*+Hk 2()EI 1*
-(
)
ρA 1
()
EI 1
*
+ρG s k 2ω2⋅-F 1-
α1
()
EI 1
*
-F 2=0
(39)
éë
êêêù
û
úú
ú1+l 22k 21+l 21k 2k 4+αj +αj -1()EI j *+Hk 2
()
EI j *
-(
)
ρA j
()
EI j
*
+ρG s k 2ω2⋅-F j -αj -1()EI j *
-F j -1-αj ()
EI j
*
-F j +1=0
(40)
éë
êêêêùûú
úúú1+l 22k 21+l 21k 2k 4+αN -1+β0+β1a 2k 2()EI N *+Hk 2()EI N
*-æèçççöø÷÷÷ρA N ()
EI N *+ρG s k 2ω2⋅-F N -αN -1()
EI N
*
-F N -1=0(41)
为了使上述关于-F j 的方程组具有非平凡解，
方程组的系数矩阵所构成的行列式值应为0，即
det f (ω,k )=0
(42)
以双壁碳纳米管为例，有
f (ω,k )=
(
)
D 11D 12D 21D 22
(43)
式中：
D 11=1+l 22k 21+l 21k 2k 4+α1()EI 1*+Hk 2()EI 1*-é
ëêêêêρA 1()EI 1*
+ρG s k 2ù
ûúúúúω2
；D 12=-α1()EI 1*；D 21=-α1()EI 2*
；D 22=1+l 22k 21+l 21k 2
k 4+α1+β0+β1a 2k 2+Hk 2
()EI 2*-(
ρA 2()EI 2*
+ρG s
k 2)
ω2。

对上述行列式进行求解，得
(m
1
-m 2ω2)-F 1+m 3-F 2=0；
n 3-F 1+()n 1-n 2ω2-F 2=0。

其中：
m 1=1+l 22k 21+l 21k 2
k 4+α1+Hk 2
()
EI 1*
；m 2=
ρA 1
()
EI 1
*
+ρG s k 2；m 3=-α1()
EI 1*；n 1=1+l 22k 21+l 21
k 2
k 4
+2295
第51卷中南大学学报(自然科学版)
α1+β
+β
1
a2k2+Hk2
()
EI
2
*
；n
2
=
ρA
2
()
EI
2
*
+
ρ
G
s
k2；n
3
=
-
α
1
()
EI
2
*。

弯曲波在双壁碳纳米管中传播时的相速度为
c
1,2
=
3.2多壁碳纳米管的结果分析与讨论
为了进一步研究范德华力对相速度的影响，以双壁碳纳米管为例，针对范德华力设置以下3种情况：1)考虑范德华力对相速度的影响；2)忽略范德华力对相速度的影响；3)假设范德华力趋于
无穷大。

范德华作用系数α
1由式(33)得到，α
1
=0
或α
1→∞。

当α1=0时，双壁碳纳米管可以看作是由2个独立的单壁碳纳米管组成[22]；当α
1→∞时，双壁碳纳米管相当于2个碳纳米管固结在一起，此时，双壁碳纳米管也可看作是1个单壁碳纳
米管，且横截面面积A=A
1+A
2
，惯性矩I=I
1
+
I
2
[23−24]。

验证结果如图8所示。

事实上，α
1=0和α
1→∞是2种理论上的极
限状态。

α
1
取由式(33)确定的值时，对应的波速位于上述2种情况下的波速之间。

同时，图8也说明范德华作用只对较小波数范围内的相速度有一定影响，而对大波数的相速度影响较小。

图9与图10所示分别为弹性介质及自由空间中几种碳纳米管相速度的对比结果。

尺度因子e
1
=
0.8，e
2
=0.1，弹性介质中弹性参数β
/E=0.01，β
1
/E=20，自由空间中弹性参数设为0。

研究对象
包含单壁碳纳米管(SWCNT)、双壁碳纳米管
(DWCNT)、三壁碳纳米管(TWCNT)及四壁碳纳米
管(QWCNT)。

由图9与图10可知：当波数
ka>1时，弹性介
质对多壁碳纳米管相速度的影响可以忽略不计；
当ka<1时，弹性介质的存在会使相速度的预测结
果偏大[25]，且波数越小，弹性介质对多壁碳纳米管
相速度的影响越明显；当多壁碳纳米管处于弹性
介质中时，随着碳纳米管管层数增加，相速度反图8自由空间中，范德华作用系数α
1
对(5，5)@(10，10)
扶手型双壁碳纳米管相速度的比较
Fig.8Comparison of phase speed of(5,5)@(10,10)
armchair DWCNT in free space duo to van der Waals
interaction coefficientα
1
图10自由空间中多壁碳纳米管相速度的比较
Fig.10Comparison of phase velocity of MWCNTs in
free space
图9弹性介质中多壁碳纳米管相速度的比较
Fig.9Comparison of phase velocity of MWCNTs in
elastic medium
2296
第8期
黄彬，等：表面效应对碳纳米管中弯曲波波动特性的影响而逐渐减小；当多壁碳纳米管处于自由空间中时，相速度随着碳纳米管管层数增加而变大[26]。

因此，碳纳米管管层数越少，相速度对弹性介质的反应越敏感。

为了研究表面效应对多壁碳纳米管相速度的影响，给出考虑表面效应与不考虑表面效应这2种情况下双壁碳纳米管中弯曲波相速度的对比结果，如图11所示。

从图11可知：在小波数范围内，2条相速度曲线存在明显偏差；随着波数增大，这种偏差逐渐减小直至消失。

这也说明表面效应影响多壁碳纳米管中弯曲波的特性，且影响范围主要集中在小波数范围内。

4结论
1)表面效应影响碳纳米管中弯曲波的相速度，
且这种影响主要体现在小波数范围内。

2)对于小波数范围，尺度系数e 1和e 2对碳纳米管相速度的影响可以忽略。

随着波数增大，尺度系数e 1和e 2对碳纳米管相速度的影响逐渐增强。

3)弹性介质对碳纳米管中弯曲波相速度的影响主要体现在小波数范围内，且刚度系数β0与β1对碳纳米管中弯曲波相速度的影响范围不同。

波数越小，弹性介质对多壁碳纳米管相速度的影响越明显。

参考文献：
[1]IIJIMA S.Helical microtubules of graphitic carbon[J].Nature,1991,354(6348):56−58.
[2]
ERINGEN A C.On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves[J].Journal of Applied Physics,1983,54(9):4703−4710.
[3]LIEW K M,HU Yangao,HE X Q.Flexural wave propagation in single-walled carbon nanotubes[J].Journal of Computational and Theoretical Nanoscience,2008,5(4):581−586.
[4]PANG M,ZHANG Y Q,CHEN W Q.Transverse wave propagation in viscoelastic single-walled carbon nanotubes with small scale and surface effects[J].Journal of Applied Physics,2015,117(2):024305.
[5]余阳,杨洋.基于非局部应变梯度欧拉梁模型的充流单壁碳纳米管波动分析[J].振动与冲击,2017,36(8):1−8.YU Yang,YANG Yang.Wave propagation of fluid-filled single-walled carbon nanotubes based on the nonlocal-strain gradient theory[J].Journal of Vibration and Shock,2017,36(8):1−8.
[6]王碧蓉,邓子辰,徐晓建.基于梯度理论的碳纳米管弯曲波传播规律的研究[J].西北工业大学学报,2013,31(5):774−778.
WANG Birong,DENG Zichen,XU Xiaojian.Modified Timoshenko beam models for flexural wave dispersion in carbon nanotubes with shear deformation considered[J].Journal of Northwestern Polytechnical University,2013,31(5):774−778.
[7]
SHEN Jun,WU Jingxiang,SONG Jie,et al.Flexural waves of carbon nanotubes based on generalized gradient elasticity [J].Physica Status Solidi(b),2012,249(1):50−57.
[8]GAFOUR Y ,ZIDOUR M,TOUNSI A,et al.Sound wave propagation in zigzag double-walled carbon nanotubes embedded in an elastic medium using nonlocal elasticity theory[J].Physica E:Low-Dimensional Systems and Nanostructures,2013,48:118−123.
[9]CAMMARATA R C.Surface and interface stress effects on interfacial and nanostructured materials[J].Materials Science and Engineering:A,1997,237(2):180−184.
[10]LI Li,HU Yujin,LING Ling.Wave propagation in
viscoelastic single-walled carbon nanotubes with surface effect under magnetic field based on nonlocal strain gradient theory[J].Physica
E:
Low-dimensional Systems
and
Nanostructures,2016,75:118−124.
[11]ZHEN
Yaxin.
Wave
propagation
in
fluid-conveying
viscoelastic single-walled carbon nanotubes with surface
and
图11考虑表面效应与不考虑表面效应情况下双壁碳纳
米管中弯曲波的相速度
Fig.11Phase velocities of flexural waves in DWCNTs
with or without surface effect
2297
第51卷中南大学学报(自然科学版)
nonlocal effects[J].Physica E:Low-Dimensional Systems and Nanostructures,2017,86:275−279.
[12]NARENDAR S,RA VINDER S,GOPALAKRISHNAN S.
Study of non-local wave properties of nanotubes with surface effects[J].Computational Materials Science,2012, 56:179−184.
[13]AIFANTIS E C.Update on a class of gradient theories[J].
Mechanics of Materials,2003,35(3/4/5/6):259−280. [14]WANG Lifeng,HU Haiyan.Flexural wave propagation in
single-walled carbon nanotubes[J].Physical Review B, 2005,71(19):195412.
[15]JR WEA VER W,TIMOSHENKO S P,YOUNG D H.
Vibration problems in engineering[M].4th ed.Chichester, United Kingdom:John Wiley&Sons Limited,1990:537−538.
[16]MILLER R E,SHENOY V B.Size-dependent elastic
properties of nanosized structural elements[J].
Nanotechnology,2000,11(3):139−147.
[17]BENOIT J M,BUISSON J P,CHAUVET O,et al.Low-
frequency Raman studies of multiwalled carbon nanotubes: Experiments and theory[J].Physical Review B,2002,66(7): 073417.
[18]YOON J,RU C Q,MIODUCHOWSKI A.Timoshenko-
beam effects on transverse wave propagation in carbon nanotubes[J].Composites Part B:Engineering,2004,35(2): 87−93.
[19]NATSUKI T,HAYASHI T,ENDO M.Wave propagation of
carbon nanotubes embedded in an elastic medium[J].Journal of Applied Physics,2005,97(4):044307.
[20]LI Xianfang,WANG Baolin,MAI Y W.Effects of a
surrounding elastic medium on flexural waves propagating
in carbon nanotubes via nonlocal elasticity[J].Journal of Applied Physics,2008,103(7):074309.
[21]DONG Ke,WANG Xi.Wave propagation in carbon
nanotubes embedded in an elastic matrix[J].Archive of Applied Mechanics,2007,77(8):575−586.
[22]WANG Qing,ZHOU Guanyou,LIN Kuochi.Scale effect on
wave propagation of double-walled carbon nanotubes[J].
International Journal of Solids and Structures,2006,43(20): 6071−6084.
[23]TANG Yugang,LIU Ying,ZHAO Dong.Viscoelastic wave
propagation in the viscoelastic single walled carbon nanotubes based on nonlocal strain gradient theory[J].
Physica E:Low-Dimensional Systems and Nanostructures, 2016,84:202−208.
[24]KAZEMI F,HAJIAHMADI A,RAJABI M.Flexural wave
propagation analysis of single-walled carbon nanotubes using molecular structural mechanics approach[J].Materials Research Express,2019,6(9):095048.
[25]王鹏.考虑表面和非局部效应的载流碳纳米管波动分析
[D].浙江:浙江大学建筑工程学院,2017:47−50.
WANG Peng.Wave analysis of current-carrying carbon nanotubes considering surface and non-local effects[D].
Zhejiang:Zhejiang University.College of Architecture and Engineering,2017:47−50.
[26]ANSARI R,SAHMANI S.Small scale effect on vibrational
response of single-walled carbon nanotubes with different boundary conditions based on nonlocal beam models[J].
Communications in Nonlinear Science and Numerical Simulation,2012,17(4):1965−1979.
(编辑陈灿华)
2298。