有限元仿真模拟在《理论声学》教学中的应用：利用散射栅实现声子

1、有限元仿真模拟在?理论声学?教学中的应用：利用散射栅实现声子聲波在周期性复合介质即声子晶体中的传播行为是?理论声学?课程中非常重要的一块教学内容。声波在声子晶体中传播时，受到单体麦氏散射和周期性布拉格散射的共同作用，形成独特的能带构造。如何在能带构造对应的相空间实现声波调控是教学中的一大难点。比方，学生通常很难理解声子晶体形式在外表栅构造散射下波矢量平行分量的增加效应。文章提出在声子晶体教学中引入计算机仿真模拟1，2，基于SLultiphysis软件开展有限元数值计算，以带有外表散射栅的声子晶体的透射性质为例，向学生明晰展示了通过添加外表栅构造实现的声波相空间的有效调控，并以此为根底在声子晶体

2、带隙频率范围内实现了声波的定向激发，为新型声学器件的设计提供了有意的理论指导3-9。在这种仿真模拟教学环境下，学生能很方便的自主设置各种参数并查看结果，有利于加深学生对声波在复合介质中折射性质的理解和掌握，具有较好的教学效果。二、仿真教学模型的设计及结果分析我们考虑一个由水圆柱四方排列于水银基体中而形成的二维声子晶体，水圆柱的半径被设定为0.4，晶体排列的周期被设定为a=1。首先我们利用SLultiphysis有限元软件数值计算了该声子晶体的带构造，图1给出了计算所得的第一条能带。在计算中所用的材料参数如下：水银的密度为p=1.35104kg/3、纵波波速为=1.45103/s；水的密度为p=

3、1.0103kg/3、纵波波速为=1.49103/s。从图中可以看到，该声子晶体在归一化频率0.2到0.28之间存在着一个沿X方向的方向带隙。根据作者多年的教学经历，此时，在参照?理论声学?中有关声波带隙的知识点晶体中的带隙是声波被制止传播的频率范围后，多数学生很容易武断地得出结论，认为位于此频率范围的声波在任何情况下都不能沿着晶体的X方向传播。事实上，通过引入外表栅构造，入射声波波矢的平行分量将获得增量，这种相空间的调制作用将使晶体中本来制止声波传播的方向变为允许传播的状态，从而实现由禁带向导带的转变。在传统的教学中，带有外表栅构造的声子晶体其色散计算比拟复杂，很难给出解析结果，通常只能通过

4、光学中的光栅方程作简单的类比，于是进一步向学生分析以上由禁带向导带的转变将变得比拟困难。在此，我们提出利用基于SLultiphysis软件的有限元方法数值计算带有外表栅构造的声子晶体色散性质，利用该软件的本征模态分析模块可以计算出声子晶体本征形式相空间的等频色散线图。该图保存了声子晶体相空间中的所有形式信息，并可利用SL软件非常方面地调出每一个单独形式进展分析，从而进一步确定每种形式对应的物理图像。详细过程如下，我们在一块沿着X方向排列的声子晶体板系统中引入如图2a所示的散射栅构造。该散射栅构造由位于声子晶体板左右两侧对称分布的两排水圆柱组成。散射栅的半径被设定为0.24，栅构造沿着平行于声子

5、晶体板方向的周期为2等于两倍晶体周期，栅构造与声子晶体板之间的间隔 被设定为0.15。考虑归一化频率为0.2694的声波从图2a中左面入射。在没有散射栅构造时，入射声波将遇到晶体-基体沿X方向的界面。图2b中深黑色的虚线半圆给出了归一化频率为0.2694的声波在水银基体中的等频色散线，图2b中实线圆那么代表该频率下声子晶体中对应的等频色散线。从图2b中我们可以看出，由于此时频率较低，声波在水银基体中的等频色散圆沿X方向与晶体中等频色散线没有任何交点，这意味着当没有散射栅时，在该频率下从基体中以任意角度入射的声波都将在声子晶体板界面被全反射，这一分析结果验证了我们从图1中获得的带隙结论此时的分析

6、结果与学生在课本中学习到的带隙的概念是一致的；当存在散射栅时，由于栅构造在平行晶体板方向具有周期性，栅构造带来的周期性散射使得入射波矢的平行分量获得一个增量？啄k=n2？仔/T，此处T=2a为散射栅沿平行晶体方向的周期，n为任意整数。对于一阶情形n=1，当频率为0.2694时，从基体中以任意角度入射的声波首先将在左侧散射栅的作用下获得一个平行于晶体界面的波矢增量？啄k=0.52？仔/a，这在图2b中相当于左侧的深黑色虚线半圆集体向上平移了0.52？仔/a的间隔 ，于是本来该全反射的入射声波将按照广义的折射定律参照图2b中的晶体等频色散线实线圆发生折射。根据作者的教学经历，通常分析至此时，学生已

7、经可以对由散射栅导致的禁带向导带的转变有了一个较为清楚的认识。我们进一步注意到，在该频率下实线圆半径相对于深黑色虚线半圆来说比拟小，而且此时晶体中等频色散线随频率增加逐渐向点会聚形成以点为中心的负折射色散态，也就是说，此时入射声波将受到晶体的微量负折射作用。于是，假如在图2a中复合体系的左侧放置一个频率为0.2694的点源，在声子晶体和散射栅的共同作用下，在复合体系的右侧我们将得到一束具有准直效果的出射声波10。为本文由论文联盟搜集整理了验证以上分析，我们利用有限元方法数值模拟了图2a中复合体系对频率为0.2694的输入点源的作用效果，结果由图2给出。从图中可以看到，在散射栅声子晶体复合系统的作用下，我们得到了较好的声定向激发效果。以上的数值分析模拟过程可以使学生对本节教学内容产生形象的认识，从而加深他们对外表栅的相空间调控这一知识难点的理解和掌握。三、結束语计算机仿真教学能有效帮助学生形成对知识点的形象认