有限元仿真模拟在《理论声学》教学中的应用：利用散射栅实现声子晶体相空间模式调控.doc

有限元仿真模拟在《理论声学》教学中的应用：利用散射栅实现声子晶体相空间模式调控聲波在周期性复合介质中的传播行为是《理论声学》课程中非常重要的一块教学内容声波在声子晶体中传播时，受到单体麦氏散射和周期性布拉格散射的共同作用，形成独特的能带结构如何在能带结构对应的相空间实现声波调控是教学中的一大难点比如，学生通常很难理解声子晶体模式在表面栅结构散射下波矢量平行分量的增加效应文章提出在声子晶体教学中引入计算机仿真模拟[1，2]，基于 COMSOL Multiphysics 软件开展有限元数值计算，以带有表面散射栅的声子晶体的透射性质为例，向学生清晰展示了通过添加表面栅结构实现的声波相空间的有效调控，并以此为基础在声子晶体带隙频率范围内实现了声波的定向激发，为新型声学器件的设计提供了有意的理论指导[3-9]在这种仿真模拟教学环境下，学生能很方便的自主设置各种参数并查看结果，有利于加深学生对声波在复合介质中折射性质的理解和掌握，具有较好的教学效果 二、仿真教学模型的设计及结果分析 我们考虑一个由水圆柱四方排列于水银基体中而形成的二维声子晶体，水圆柱的半径被设定为，晶体排列的周期被设定为 a=1mm。

首先我们利用 COMSOL Multiphysics 有限元软件数值计算了该声子晶体的带结构，图 1 给出了计算所得的第一条能带在计算中所用的材料参数如下：水银的密度为 p=×104kg/m3、纵波波速为c=×103m/s；水的密度为 p=×103kg/m3、纵波波速为cw=×103m/s从图中可以看到，该声子晶体在归一化频率到之间存在着一个沿 X 方向的方向带隙根据作者多年的教学经验，此时，在参照《理论声学》中有关声波带隙的知识点“晶体中的带隙是声波被禁止传播的频率范围” 后，多数学生很容易武断地得出结论，认为位于此频率范围的声波在任何情况下都不能沿着晶体的 X 方向传播事实上，通过引入表面栅结构，入射声波波矢的平行分量将获得增量，这种相空间的调制作用将使晶体中原本禁止声波传播的方向变为允许传播的状态，从而实现由禁带向导带的转变在传统的教学中，带有表面栅结构的声子晶体其色散计算比较复杂，很难给出解析结果，通常只能通过光学中的光栅方程作简单的类比，于是进一步向学生分析以上“由禁带向导带”的转变将变得比较困难 在此，我们提出利用基于 COMSOL Multiphysics 软件的有限元方法数值计算带有表面栅结构的声子晶体色散性质，利用该软件的“本征模态分析” 模块可以计算出声子晶体本征模式相空间的等频色散线图。

该图保留了声子晶体相空间中的所有模式信息，并可利用 COMSOL 软件非常方面地调出每一个单独模式进行分析，从而进一步确定每种模式对应的物理图像具体过程如下，我们在一块沿着 X方向排列的声子晶体板系统中引入如图 2 所示的散射栅结构该散射栅结构由位于声子晶体板左右两侧对称分布的两排水圆柱组成散射栅的半径被设定为，栅结构沿着平行于声子晶体板方向的周期为 2mm，栅结构与声子晶体板之间的距离被设定为考虑归一化频率为的声波从图 2 中左面入射在没有散射栅结构时，入射声波将遇到晶体-基体沿 X 方向的界面图 2 中深黑色的虚线半圆给出了归一化频率为的声波在水银基体中的等频色散线，图 2 中实线圆则代表该频率下声子晶体中对应的等频色散线从图 2中我们可以看出，由于此时频率较低，声波在水银基体中的等频色散圆沿 X 方向与晶体中等频色散线没有任何交点，这意味着当没有散射栅时，在该频率下从基体中以任意角度入射的声波都将在声子晶体板界面被全反射，这一分析结果验证了我们从图 1 中获得的带隙结论；当存在散射栅时，由于栅结构在平行晶体板方向具有周期性，栅结构带来的周期性散射使得入射波矢的平行分量获得一个增量？啄 k=n，此处 T=2a 为散射栅沿平行晶体方向的周期， n 为任意整数。

对于一阶情形，当频率为时，从基体中以任意角度入射的声波首先将在左侧散射栅的作用下获得一个平行于晶体界面的波矢增量？啄 k=，这在图 2 中相当于左侧的深黑色虚线半圆集体向上平移了的距离，于是本来该全反射的入射声波将按照广义的折射定律参照图 2 中的晶体等频色散线发生折射根据作者的教学经验，通常分析至此时，学生已经能够对由散射栅导致的“禁带向导带” 的转变有了一个较为清楚的认识我们进一步注意到，在该频率下实线圆半径相对于深黑色虚线半圆来说比较小，而且此时晶体中等频色散线随频率增加逐渐向 M 点汇聚，也就是说，此时入射声波将受到晶体的微量负折射作用于是，如果在图 2 中复合体系的左侧放置一个频率为的点源，在声子晶体和散射栅的共同作用下，在复合体系的右侧我们将得到一束具有准直效果的出射声波[10]为联盟了验证以上分析，我们利用有限元方法数值模拟了图 2 中复合体系对频率为的输入点源的作用效果，结果由图 2 给出从图中可以看到，在散射栅——声子晶体复合系统的作用下，我们得到了较好的声定向激发效果以上的数值分析模拟过程能够使学生对本节教学内容产生形象的认识，从而加深他们对表面栅的相空间调控这一知识难点的理解和掌握。

三、結束语  　　计算机仿真教学能有效帮助学生形成对知识点的形象认识，是一项攻克知识难点的强有力的教学手段文章以声子晶体表面栅结构对声波相空间的调控作用这一知识难点为例，探讨了仿真模拟教学在《理论声学》课程中的应用，通过将难以理解的知识点转化为模拟结果中生动形象的仿真结果，得到了较好的教学效果。