广义多极子技术在含任意夹杂的二维声子晶体标量波能带结构计算中的应用.pdfVIP

广义多极子技术在含任意夹杂的二维声子晶体标量波能 带结构计算中的应用 史志杰1 汪越胜1 张传(增)2 1北京交通大学工程力学研究所，100044 2锡根 大学 土木 工程 系 ，锡 根D-57068,德 国 中 文 摘 要 = 摘 要 基 于 广 义 多 极 子 技 术 提 出 了 用 于 计 算 具 有 任 意 形 状 夹 杂 的 二 维 声 子 晶 体 标 量 波 带 结 构 的 多 重 单 极 法 。 该 方 法 中 除 了 引 入 用 以 展 开 波 场 的 点 源 外 ， 在 适 当 位 ( 置 ) 又 额 外 引 入 一 个 单 极 点 源 作 为 激 励 。 改 变 激 励 的 频 率 ， 特 征 频 率 则 以 一 个 适 当 选 择 的 响 应 函 数 极 值 的 形 式 确 定 。 通 过 逦 历 所 要 计 算 的 频 率 范 围 及 第 一 不 可 约 布 里 渊 区 ( F B Z ) 的 边 界 即 可 得 到 带 结 构 。 通 过 与 有 限 元 法 ( F E M ) 的 结 果 比 较 表 明 该 方 法 可 以 准 确 地 计 算 非 圆 夹 杂 的 固 ／ 固 体 系 声 子 晶 体 的 带 结 构 。 该 方 法 具 有 计 算 速 度 快 ， 建 模 简 单 ， 模 型 参 数 与 材 料 参 数 无 关 以 及 可 以 计 算 任 意 区 间 频 率 的 特 点 。 中 文 关 键 词 ＝ 关 键 词 ： 广 义 多 极 子 技 术 ; 多 重 单 极 法 ; 带 结 构 ; 特 征 值 问 题 根据广义多极子技术，基体和夹杂中的位移司分别用F式逼近： (1) 甜r(x)=“。(x)+∑竺。q，。联”(尼。lx—x。，。1)，x=(x，y)∈Q。， (2) “男“(x)=∑竺．口：，。Hjl’(尼：Ix—x：，。1)，x=(x，y)∈Q：， 式(1)中“。(x)表示由外部激励引起的位移场，x抽表示第，价多极源。“野“及“男” 在界质面和单胞外边界上分别满足连续性条件及B10ch条件。在介质面及单胞外边界上 根据有最小二乘法可 取配置点可得超定矩阵方程：A，(缈，t，七，)．a=b(国，t，七，)。 得最优权值a掣=(ar，a罗)2。 由此可得石(∞，也，砖)=√1／Ⅳf∑羔。卜：’(x“)12，式 中x“为基体Q，内随机分布的检测点。 另外可定义误差范数： 五(国，t，砖)=0A，(国，t，砖)·a叫一b(缈，戤，足，)¨根据z和厶可定义特征值寻找函 数：石(∞，t，砖)=(1／石(∞，t，砖))”p／(厶(国，t，b))p。给定不可约第一布里渊区 边界上任意一点(t，砖)，遍历所要计算的频率范围即可得到特征值寻找函数 六(∞，颤，七，)的极小值点。这些极小值点就是我们所要寻找的和(颤，七，)相对应的特征 值。遍历不可约第一布里渊区的边界即可得到完整的带结构。为使特征值寻找函数的 极信点更加突出，可引入两种正则化方法，见．文献【21。 三、数值结果与结论 图l给出了多重单极模型，图2给出了剪切模态能带结构。通过与有限元结果的 比较可看出两种方法的结果吻合的非常好。 这充分证明了该方法对于有棱角的夹杂体 系也能给出非常精确的结果。这种夹杂体系对于多重散射法及狄利克雷．纽曼映射法等 常规波函数展开法是困难的。 此外该方法具有计算速度快，建模简单，模型参数与材 料参数无关以及可以计算任意区间频率的特点。 。．…，’”： ‘芝==!互==五￡二’ ■ I!．．． ．。一一一7ii■一j了二i： 1性。．茹≯—1