基于COMSOLMultiphysics的有限元仿真模拟在《理论声学》教学中的应用：缺陷态的分析-最新教育文档.docx

基于 COMSOL Multiphysics 的有限元仿真模拟在《理 论声学》教学中的应用：缺陷态的分析 概述 声子晶体是物理类专业研究生课程 《理论声学》 中非常重要 的一块教学内容。 在声子晶体中引入缺陷将导致缺陷态、 波导态 以及界面态等本征模态的出现。 如何有效计算缺陷声子晶体的色 散关系并对这些模态进行本征分析是教学中的一大难点。比如， 学生通常很难从一般的色散计算结果中有效地区分波导态以及 界面态。文章提出在教学中引入计算机仿真模拟 [1-3] ，利用基 于 COMSOMLultiphysics 软件的有限元方法数值计算缺陷声子晶 体的色散关系。 在这种全矢量的色散计算方式中， 学生能很方便 的调取所有模态的本征信息， 从而能全面而精准地掌握缺陷声子 晶体丰富多彩的色散性质，具有较好的教学效果。 有限元仿真教学模式的设计 以铁柱四方排列在空气中形成的二维声子晶体为例。 声子晶 体的填充率被设定为0.577。通过在晶体中沿着 X方向移除一排 铁柱，从而引入一条沿 X方向的线缺陷。利用基于 COMSOL Multiphysics 软件的有限元方法，我们数值计算了该缺陷声子 晶体结构的色散关系。 图 1 给出了计算所得归一化频率介于 0 到 1 之间的带结构， 其中阴影部分代表没有缺陷的完美声子晶体的 X方向投影带结 构。此处我们采用的是超胞算法， 图 1 中右下插图给出了计算中 所用的超胞结构示意图；在求解过程中选取的边界条件为 FIoquet periodicity 边界条件。 计算中所用材料参数为： 铁的密度=7.67 X 103kg/m3 纵波波速 CI=6.01 X 103m/s 横波波速 Ct=3.32X103m/s 空气的密度 =1.29X10-3kg/m3 纵波波速 C=0.34X103m/s 从图 1 中我们可以很清晰的看到， 在声子晶体的完全带隙中 频率大约介于 0.4 和 0.7 之间）出现了三条导带， 分别被标记 为带 1、带 2 以及带 3。很明显，这三条位导带位于完整声子晶 体的带隙频率范围内， 其产生的机理一定与晶体中线缺陷结构的 引入有关。 根据作者多年教学的经验， 大多数学生此时很容易先 入为主地做出判断，认为带 1 、带 2 以及带 3 均为由缺陷引入的 波导态。事实上，在这三条带中，带 2和带 3 是波导态，而带 1 则是与缺陷结构引入的表面所对应的表面态。在传统的教学中， 声子晶体的带结构通常由多重散射或者时域有限差分法计算得 出，而这两种方法一般很难在最终计算结果中保留所有模式的本 征信息，于是进一步向学生分析这三条导带的本征性质变得比较 困难。 文章中图 1 的色散关系由基于 COMSOML uItiphysics 软件的 有限元方法， 这款软件在计算中保留了所有模式的本征信息并且 能够非常方面的调出分析， 从而进一步确定每种模式对应的物理 图像。 在此我们提出如下教学思路： 首先，通过计算具有不同缺陷尺寸声子晶体带结构的方式向学生展示图1 中带 1 在此我们提出如下教学思路： 首先， 通过计算具有不同缺陷 尺寸声子晶体带结构的方式向学生展示图 1 中带 1 、带 2 以及带 3 随缺陷尺寸的变化趋势，所计算结果由图 2 （a）表示。 图 2（ a ）中不同形状的数据点代表具有不同缺陷尺寸的声 子晶体结构，从该图可以看出，随着缺陷尺寸的改变，带 2 和带 的位置变化非常明显，而带 1 的位置几乎不变。 这说明带 2 和带 3 与缺陷性质密切相关（由此我们推断带 2 和带 3 是线缺陷中的波导态），但带 1和缺陷的关系并不大；其 次，通过在图 1 中引入水线（此处即为空气中的声速线，由图中 直线表示）不难发现，带 1 始终位于体系水线以外，这说明带 1 中的模式在空气中始终是衰逝的，考虑到这些模式同时也在晶体 中的模式在空气中始终是衰逝的， 考虑到这些模式同时也在晶体 内衰逝（因为带 1 位于晶体的带隙范围内），我们不难推断出带 1 其实是与缺陷两侧晶体表面所对应的表面态。 根据作者的教学经验， 通常分析到这一步时， 学生已经对这 三条导带的性质有了一个较为清楚的认识；接着，我们调出 COMSOL Multiphysics 中保存的本征模式形态以便于学生对以上 分析结果形成感性认识。 图2 （b）和（c）分别给出了带3和带2上任意一点在缺陷 中的本征声压分布， 至此我们不难看出， 带 3 和带 2 分别对应着 波导模式中的对称以及反对称模式 [4] ；最后，为了进一步形象 化以上讨论结果，我们从带 1、带 2 以及带 3 中各自选取一个频 率，将工作在这三个频率下的点源分别放置在由图 3(a)、( b)、 (c)所示的体系中，我们利用 COMSOL Multi physics