二维声子晶体的等效非局部弹性介质徐伟汪越胜（北京交通大学工程.docVIP

二维声子晶体的等效非局部弹性介质 徐 伟 汪越胜 （北京交通大学工程力学研究所，北京，100044） njtu.edu.cn） 【摘要】本文提出一种新的等效连续介质模型，即将声子晶体等效为非局部各向异性弹性介质。以二维声子晶体SH波传播为例，利用频散关系匹配的方法，并结合二元快速Fourier反演，得到等效非局部弹性介质的影响函数，进一步建立了动力学方程。 【关键词】 弹性波，频散关系，声子晶体，非局部，均匀化 1 引言 近10多年来，人们对复合介质的研究有着浓厚的兴趣。当声波或弹性波在周期性复合介质中传播时可能产生弹性波禁带，使得某些频率范围的波不能通过介质。对应于这样的周期复合介质，Kushwaha等[1]人于1993年提出了声子晶体的概念，并对Ni-Al四方排列的周期复合结构材料，利用平面波展开法计算获得了在剪切模态下的弹性波禁带。1995年Martinez-Sala等[8]对马德里的一座200年前制作的雕塑“流动的旋律”进行声学特性测试，第一次从实验上证实了弹性波禁带的存在。声子晶体作为一种声功能材料，可能通过对其周期结构进行设计，人为地调控声波的流动，因而在声滤波器、声波导、声变送器、声学透镜以及减震降噪控制方面有着广泛的应用前景，相关的研究已引起了国内外物理学家的极大关注[9]。人们研究发现，声子晶体具有很多美妙的特性[9]，如表面态、缺陷态以及负折射率等，这对于制作声表面器件、声波导器件有着重要的意义。 利用均匀化观点来研究非均匀介质的宏观平均特性，一直是力学所关注的问题。不管是具有周期性的还是完全随机的非均匀介质，考虑波动效应的等效连续均匀介质模型大多是在低频范围内建立的，其目的是用各种计算方法获得等效的波速或材料参数。一般情况下得到的参数是随频率变化的，在极低的频率内则近似为常数。尽管这样能得到等效的材料参数，但是很难得到真正意义上的动力学方程，而这正是力学所关心的问题，因为没有动力学方程，就很难用连续介质力学的方法对复杂结构的波动响应进行计算分析。 声子晶体是一种具有周期性的典型非均匀复合介质，其等效材料参数同样受到关注。如果能提出一个合理的连续均匀介质模型，并建立真正意义上的动力学方程，对声子晶体特性的研究将具有重要的意义。 经典的非局部理论由于考虑了材料微观结构的效应而在处理波长与晶格尺寸相近时的波动问题时获得了成功。这给新模型的建立提供了思路，即可以考虑把声子晶体等效为非局部各向异性弹性介质。声子晶体内部结构具有严格的周期性，也正是这种周期性造成了其特殊的频散关系。为此，可以考虑利用频散关系匹配的方法求得非局部弹性模量。文献[2]研究了一维周期平板结构（相当于一维声子晶体）的非局部等效介质模型，并用频散关系匹配的方法得到了非局部弹性模量。本文利用同样的思路，把该方法推广到二维情况，建立二维声子晶体的等效非局部弹性介质模型，而这种非局部介质是立方各向异性的。模型的建立分三个步骤：(1) 用平面波展开法计算二维声子晶体SH波频散关系；(2) 基于非局部弹性理论推导立方各向异性介质SH波频散关系；(3) 将两个频散关系匹配，得到等效非局部介质的弹性模量，进而获得动力学方程。 2 二维声子晶体SH波频散关系 考虑偏振方向平行于散射体的SH波。这种情况下波动方程写为[1]： (1) 其中，研究Ni圆柱-Mg基正方形排列声子晶体。散射体半径为R，晶格常数为a。用平面波展开法可计算得到频散关系： (2) 以简约频率作为纵坐标，简约波矢作为横坐标画出前10个能带，如图1所示。 3 非局部线弹性理论 非局部弹性理论考虑了物质内部特征结构与原子间长程影响，认为一点的应力与区域内所有点的应变有关[4-7]。对于各向异性材料，非局部本构关系可以写成下面的形式： (3) 其中， 为非局部弹性模量。为了问题研究的需要，引入非局部影响函数，这样可以把非局部弹性模量表示成： (4) 其中为经典的弹性常数。 非局部弹性理论中，几何方程、运动方程，与经典弹性理几何方程、应力运动方程一样，分别为： (5) (6) 对于SH波，只有位移分量： (7) 此时运动方程为：