1. 引言
弹性波在固体中传播时,当固体中存在夹杂等缺陷时周围会发生弹性波散射的现象,从而引起动应力集中 [1] 。弹性波散射问题中最典型的就是SH波的散射。宏观情况下半空间内可移动夹杂对SH波散射问题已有了一些成果 [2] [3] ,更多的研究成果可以查阅专著 [4] 。随着纳米材料和纳米元器件广泛的研究和应用,人们对材料性能的要求越来越高,纳米量级的材料和纳米元器件更是得到了广泛的科学研究和技术应用,比如在材料科学、医疗器械、电子装置、涂料器材等行业。由于在纳米尺度下,材料表面面积与体积之比明显增大,从而导致表面效应显著,使得纳米材料与纳米元器件的力学性能不同于宏观尺度下的力学性能。因此,当原子直径远小于材料直径时,表面效应可以不考虑,但是当材料尺寸缩小到纳米尺度时,其对材料的性能有非常大的影响,此时就必须考虑表面效应。在1906年,Gibbs [5] 首次提出了固体表面张力。经过多年的发展,1975年Gurtin和Murdoch [6] [7] 建立了研究材料表面力学行为的数学模型,利用力和力矩平衡建立了表面应力的张量性质,研究了表面与基体应力的关系,形成了表面弹性理论。近年来人们利用表面弹性理论研究了在剪切载荷作用下纳米接触问题 [8] ,两个孔洞(夹杂)相互作用对P波的散射 [9] [10] ,以及SH波诱发的阵列纳米孔洞周围的弹性波散射和动应力集中 [11] 。本文研究了弹性半空间内浅埋纳米可移动夹杂对SH波的散射。利用设虚源法 [4] 得到预先满足半空间水平边界处应力自由的散射波函数,结合复变函数法、波函数展开法和Graf加法公式,得到了固体中的应力场,讨论了半空间内浅埋纳米移动夹杂对稳态入射平面SH波的散射影响,计算了夹杂周围的应力强度因子。
2. 理论模型与分析
2.1. 问题描述
如图1所示,弹性半空间介质内含有一个密度
和剪切模量
的可动刚性夹杂(
),夹杂的半径为R,掩埋深度为H。假设基体的剪切模量和质量体密度分别为
,
。一列稳态入射平面SH波
(入射波函数)在介质内传播,由于可动夹杂的存在,它将与半空间边界作用产生一列反射波
,与夹杂相互作用产生散射波
。引入复变量
,
,在复平面
下,它们都满足如下Helmholtz方程
(1)
Figure 1. Scattering of SH waves by shallow moving inclusions in semi-space at nanoscale
图1. 纳米尺度下半空间内浅埋可动夹杂对SH波的散射问题
其中
为波数,
为圆周频率,
为波速。文中忽略了时间因子
。
在复变坐标系下,经典弹性力学中应力和位移的关系为
,
. (2)
其中
,
,
,i为虚数单位,
和
分别是反平面问题中的径向和环向的切应力。
入射SH波的波函数可表示为 [3]
(3)
反射波的波函数可表示为
(4)
其中
为入射波振幅。
由于入射波
,反射波
与可动夹杂及其边界共同作用产生了散射波
,利用设虚源法,得到在复平面下满足半空间水平边界处应力自由的散射波函数为
(5)
其中
是待定的系数,
表示n阶第一类Hankel函数。
利用Graf加法公式,
(6)
将(6)代入(5),散射波函数可以表示为
(7)
将公式(3),(4)和(7)分别代入(2),入射波,反射波,散射波对应的应力场分别为
(8)
(9)
(10)
(11)
(12)
(13)
其中
,
,
分别为入射波,反射波,散射波对应的径向应力。
,
,
分别为入射波,反射波,散射波对应的环向应力。
利用叠加原理,介质内总应力场
、
可以表示为
(14)
(15)
介质内的总位移场为
(16)
2.2. 边界条件
不计表面残余张力,表面应力与表面应变的关系为
(17)
其中
和
为材料表面的Lamé常数,
为Kronecker符号。
假设弹性表面与基体并无滑移,则表面处的平衡方程为
(18)
(19)
其中
为界面的应力张量,
为表面的外法向向量,
为表面的曲率张量,希腊字母
、
取1、2、3,拉丁字母i、j均取1、2,下标中的逗号表示求导。
在
处,根据(18)式和(19)式,得到夹杂表面处边界条件为
(20)
根据表界面本构关系式(17)得到界面的切应力
为
(21)
将(21)式代入(20)式,得到边界条件为
(22)
式中
,s为纳米尺度下衡量表面效应的参数,可以看出s和夹杂半径R成反比关系。对于宏观问题,半径R非常大,那么s远远小于1,因此表面效应可以忽略不计,但缩小R到纳米尺度时,s的作用就非常显著,此时必须要考虑表面效应 [11] 。
2.3. 问题求解
由于刚性动夹杂内部产生的应力和应变场可以忽略不计,仅仅考虑其反平面运动。进一步假设动夹杂本身的反平面运动
,利用牛顿第二定律,列出动夹杂满足的运动方程为
(23)
其中单位厚度上的质量
,加速度
,则(23)变为
(24)
在上式中考虑到
,
,及
,
,则又有
(25)
根据边界条件(22)和牛顿第二定律,我们有
(26)
由于夹杂和基体之间是密切接触的,因此基体内的总位移场在夹杂处应该与夹杂的反平面位移一致。即得到下面这个等式
(27)
将(3)、(4)和(5)式代入(27)式得到
(28)
式中
。
将(15)和(28)式代入(26)式,得到关于待定系数
的无穷代数方程组
(29)
其中
(30)
(31)
其中
,
。
在具体求解时,在方程(29)两边同时乘以
,并在
内积分。
(32)
在上式中
(33)
根据上面的理论推导,给出了半空间内浅埋纳米移动夹杂对SH波散射的无穷代数方程组,利用截取有限项的方法求得待定系数
。
1) 当
时
(34)
2) 当
时
(35)
3) 当
时
(36)
最后将数值结果代入(14)、(15)和(16)式可以得到介质内的总位移场
和总应力场
与
。
3. 算例与结果分析
研究可动夹杂边界处的动应力集中程度,通常采用动应力集中因子DSCF来表示,定义如下
(37)
其中
为入射应力的最大幅值。
下面给出了动应力集中因子与无量纲入射波波数kR,剪切模量
,掩埋比
,入射波波幅
,入射角
,质量体密度比
和表面参数s关系的数值结果。
图2给出了入射角
入射时表面效应对可动夹杂边界处动应力集中因子的影响,取入射波数
,掩埋比
,入射波波幅
,密度比
。当s趋近于0时得到的动应力集中因子与宏观条件下的结果吻合,这正好验证了本文的正确性。如图所示,考虑表面效应之后,由于入射波和反射波的相互干涉,曲线出现了凸起,随着表面参数s的增加,在整个区间内动应力集中因子逐渐减小。
Figure 2. The influence of surface effect on DSCF at the boundary of movable inclusion at the α = π/2
图2. α = π/2入射时表面效应对可动夹杂边界处DSCF的影响
图3给出了入射角
入射时波数kR对动应力集中因子的影响,取密度比
,表面参数
。然后分别取波数
、
和
观察无量纲波数kR的变化对动应力集中因子的影响。从图中可以得出,随着kR的增大,动应力集中因子在整个区间内都是增大的。
图4考虑了纳米尺度下质量体密度比对动应力集中因子的影响,我们固定入射角
,入射波波数kR,入射波波幅
,掩埋比
。从图中可以得出,随着质量体密度比
的不断增大,动应力集中因子在整个区间内都是增大的。
图5考虑了动应力集中因子随着入射角
的变化情况,我们固定入射波波数kR,入射波波幅
,
Figure 3. The influence of the change of wave number kR on DSCF at α = π/2 incidence at nanoscale
图3. 纳米尺度下α = π/2入射时波数kR的变化对DSCF的影响
Figure 4. Effect of bulk density ratio on DSCF at nanometer scale
图4. 纳米尺度下质量体密度比对DSCF的影响
Figure 5. The influence of the change of incident angle α on DSCF
图5. 入射角α的变化对DSCF的影响
掩埋比
,质量体密度比
。从图中可以得出,入射角
度时,动应力集中因子最大值为5,入射角
入射时,半空间内可移动夹杂周围动应力集中因子在整个区间内是关于纵轴对称的。
4. 结论
本文研究了半空间内浅埋纳米移动夹杂对SH波散射的影响。计算结果主要考虑了入射角以
入射时可动夹杂边界处动应力集中因子的变化情况,同时也讨论了其他入射角对动应力集中因子的影响。结果表明:当表面参数s趋近于0时得到的动应力集中因子与宏观条件下的结果具有一致性,考虑表面效应之后,随着表面参数s的增加,在整个区间内动应力集中因子逐渐减小,表面效应显著。不管是高频还是低频,随着表面参数s的增加,在整个范围内DSCF都会减小。当表面效应固定不变时,随着无量纲波数kR的增大,应力集中因子DSCF也会随之增大。在入射角以
入射时,DSCF在整个区间内具有对称性。
基金项目
国家自然科学基金(1106004和11362009)。