预测二维声子晶体能量禁带的FM-BEM.pdfVIP

预测二维声子晶体能量禁带的FM-BEM 1) *,2) * + 孟文辉 ，王连堂 ，崔俊芝 * (西北大学数学系，陕西 西安 710127) + (中国科学院数学与系统科学研究院，北京 100080) 摘要：声子晶体是一种具有声子带隙的周期性结构功能材料，对声子晶体能量禁带的研究是近些年声学研究的热 点问题之一。声子晶体能量禁带预测常用的计算方法有传递矩阵法、平面波展开法、时域有限差分法、多重散射 方法等，而边界元方法在该领域的应用还很少见。本文提出求解周期结构中声波传播问题的快速多极边界元方法 (FM-BEM) ，并将该方法应用于二维声子晶体能量禁带的预测中。针对散射体的周期分布方式，将整个结构划分 为很多大小形同的单位快，通过适当的变换，将问题的求解转化到该单位快内，从而减少了问题的求解规模。同 时，在应用快速多极算法求解中，提出一种周期树结构的建立方式，相比传统的树结构，该树结构在算法求解中 的求解效率更高。然后，应用以上的方法对钛合金(Ti-alloy)和碳化硅(SiC)组成的声子晶体的能量禁带进行了预测， 得到了较好的结果。最后，对算法的求解效率进行了分析，从中可以看出本文提出的求解周期结构中声波传播问 题的FM-BEM 方法能够高效准确的预测二维声子晶体的能量禁带。 关键词：FM-BEM ；周期结构；声子晶体；能量禁带 引 言 快速多极算法(FMM)最早由 Rokhlin 和 Greengard 在文献[1,2] 中提出，作为加速多粒子系统以及边界 积分方程数值求解的一种有效的算法，通过使用快速多极算法来求解，使得求解的存储量和计算量降低到 O(N) 。快速多极算法自提出到现在的二十多年之中，得到了迅速的发展，快速多极算法与边界元方法结合 形成了快速多极边界元方法(FM-BEM) ，该方法被应用于很多工程领域。文献[7-10] 中对快速多极算法的原 理做了全面系统的阐述。在算法应用方面，Liu 和其合作者在FM-BEM 的应用方面做了大量的工作，其中 [11] [12] [13] [14,15] 包括FMM 在位势问题 , Stokes 流问题 , 复合材料力学问题 以及声波相关问题 中的应用。此外， 文献[16,17]对快速多极边界元方法在二维和三维声波问题中的应用研究方面做了部分工作，Nishimura 等在 文献[18,19,20] 中研究了 FM-BEM 算法在周期结构中的应用。在固体力学领域的应用，姚和王做了大量的 相关研究工作[2 1,22,23] 。除此之外，FMM 算法在其它一些相关领域的应用研究工作在文献[24,25,26,27]中可 以看到。 周期结构复合材料中的声波或者弹性波的传播问题是近些年来在工程中经常碰到的问题。周期结构复 合介质是由两种或两种以上介质组成的复合体，整体上它是非均匀的，然而每种组元的内部是均匀的。复 合介质表现出不同于其组元的物理性质，复合介质的组成或结构不同，其物理性质也不同。近年来，复合 1）国家自然科学基金资助项目. 2 ）E-mail ：yonkey@. 介质的研究极其活跃，例如随机介质、电(磁)流变液、光子晶体以及声子晶体等都引起了人们的极大兴趣。 声子晶体是一种由两种或多种材料组成的具有周期性结构和弹性波带隙的声学功能材料和结构。声子晶体 就是具有声子带隙的周期性结构功能材料，它通过散射体弹性常数的周期性调制来实现声子带隙。对声子 晶体的研究具有更丰富的物理意义，也存在重大的应用价值，例如作为声音滤波器或者是隔音减振材料， 它可以更有效的隔离噪声和振动，同时又具有传统材料所不具有的特性。 对声子晶体能量带隙的研究是近些年声学研究的热点问题之一。近十年来，声子晶体的研究取得了较 [28] [29,32,33] 大的进步，用于声子晶体能带结构的理论计算方法有传递矩阵法 、平面波展开法 、时域有限差分 [30] [31] 法 、多重散射方法 等，而边界元方法在预测声子晶体带隙问题中的应用还很少见