二维三角陀螺声子晶体的波调控研究*

王伟能 杨晓东 张伟

（北京工业大学机电学院，北京 100124）

引言

超材料是指一些具有超常物理性质的人造材料，一般是自然界中所没有的.声子晶体作为一种新型超材料，在近些年受到广泛关注.声子晶体具有能带特性，波在能带结构通带范围内可以传播，在禁带内则被抑制［1］.振动和声音都可以通过弹性波的形式传播，因此声子晶体为振动控制和噪声抑制提供了新的方法手段，问题也可转化为对声子晶体中弹性波的调控研究.

声子晶体的概念是在周期结构的研究基础上类比光子晶体提出的.Sigalas［2］将球形散射体置于一种基体材料中并组成三维周期结构，发现弹性波在其中的传播存在带隙现象.Kushwaha［3］等在研究镍铝二维周期复合材料时首次明确提出了声子晶体的概念，并且给出了结构的弹性波带隙.对声子晶体能带结构的传统调控方法，主要是改变结构的几何构型，或使用能与其他物理场发生耦合的材料.Wang等［4］用有限元仿真对三维周期弹性体进行了研究，通过施加外力使材料变形来改变结构的带隙范围，材料的弹性性质使得带隙的改变具有往复性.Bergamini［5］设计了一种压电陶瓷声子晶体梁，改变通过压电陶瓷的电流频率，即可获得不同的色散关系.Zhang［6］等人对力和磁场作用下的二维磁致伸缩声子晶体板进行了研究，发现增大磁场强度会拓宽原有带隙宽度并且产生新的带隙，而增加预紧力会抵消磁场对能带结构的影响.这些调控方法依赖于材料属性或其它物理场，有的调节带隙范围窄，有的调节具有延时性.而通过其它方式，如参激调节、陀螺调节，则具有实时调控的能力，且能获得较大的禁带.Yang等［7］人研究了一维弹簧质量声子晶体，将高频振动作为参数激励提供伪刚度，实现了采用参数激励调节结构带隙.邸士杰［8］提出了一种陀螺和质量块交替连接的一维声子晶体模型，分析了陀螺转速、连接刚度等对结构带隙的影响.Brun等［9，10］将陀螺引入二维周期晶格，发现陀螺的参数可以控制结构的色散性质，并提出用这种材料来设计屏蔽装置和指定方向的波导，但其所推导的动力学方程并不能全面地体现陀螺的影响.Zhou［11］将离散陀螺周期结构的研究拓展到了陀螺连续体，并通过连续介质力学中的Hermit质量密度张量建立了声学陀螺连续介质的表面方程，解释了单向波产生的原因.李夏临［12］等研究了一种无约束的二维陀螺声子晶体，发现结构的参数可以改变带隙.

本文对二维陀螺声子晶体的波进行调控研究，用拉格朗日方程推导了结构的动力学方程，根据Bloch定理求解色散关系，分析陀螺转速对色散关系的影响，指出波的传播方向随频率的变化趋势，并对理论带隙进行了数值验证.

1 二维三角陀螺声子晶体

本文所研究的二维三角陀螺声子晶体模型如图1所示，单元体为一个上端连接质点、下端铰接的陀螺，各个单元体通过质点间的弹性连接组成平面正三角排列周期结构，在质点运动之前陀螺转轴垂直于质点组成的平面.每个质点的质量为M0，质点间的弹性连接长度为l、刚度为c、质量忽略不计.假设陀螺章动角较小，质点只在平面内平移，不随陀螺自转而转动.

图1 三角陀螺声子晶体Fig.1 Triangular gyroscope phononic crystals

定义图中处于x方向第m列（列的平面方向与x轴成60度角）、y方向第n行的质点（m，n），在x方向的位移为 xm，n、y方向的位移为 ym，n.取陀螺-质点单元进行运动分析，对章动角采取线性化假设，由拉格朗日方程可推得结构的动力学方程为：

2 色散性质分析

假设质点的运动是时谐的，即：

其中，M=dig[M，M]为质量矩阵，G=iGR为陀螺矩阵，C为刚度矩阵：

下面分析陀螺项系数G的取值对能带结构的影响，取M=1、c=1、l=1，图2为G分别取0.5、1.5和2时的色散曲面以及对应的沿不可约布里渊区边界的色散曲线，不可约布里渊区示意图位于G取0.5的色散曲线图右上角.色散图中有两支色散曲面，当G取值较小时，两支曲面具有频率重叠部分，当G增大时，两支曲面逐渐分离且变得扁平，曲面间出现带隙.由于陀螺项系数G正比于陀螺转速，在陀螺的转动惯量、转轴长度等固有属性确定的情况下，可以通过调节陀螺转速，控制结构禁带是否出现，以及出现时的带宽，从而通过或抑制指定频段的波.

图2 G不同取值的色散性质Fig.2 Dispersion properties of different value of G

图3是G取0.5时色散曲面的等频率截线，其中（a）对应下支，（b）对应上支，曲线的法向为对应频率的波的传播方向.当波的频率较低时，曲线为圆形，法线向四周各向同性散开，波可向各个方向同等传播.当波的频率较高时，截线法线集中于六个方向，表示波也集中在这六个方向传播.但上下曲面的高频范围的六个传播方向角度不同.

图3 G=0.5色散曲面截线Fig.3 G=0.5 Contours of dispersion surfaces

3 带隙的数值验证

二维声子晶体中波的传播方向有无数个，因此我们选定一个方向进行带隙验证，图4为G取4时的x轴方向色散曲线，两支色散曲线中间有较宽的带隙.对方程（1）使用MATLAB进行求解，在结构中间质点施加一个圆形轨迹位移激励，激励频率范围从0到6，观测x轴方向上与激励点有一定距离的一个质点在一段时间内的最大位移.由于结构是有限的，为避免波的反射，观测时间范围不超过波到达边界的时间.

图4 G=4 x方向色散曲线Fig.4 G=4 Dispersion curves along x direction

图5是不同频率的激励时观测到的最大位移与激励的幅值比r的变化曲线，当波的频率处在禁带（包括带隙）范围内时，波的传播受到了明显抑制，而通带内波的传播效果明显优于禁带.从数值验证的结果来看，通过调节三角陀螺声子晶体的陀螺转速，可获得能带带隙以抑制波的传播，且波受到抑制的频率范围与理论获得的带隙吻合良好.

图5 不同频率的波响应曲线Fig.5 Response of waves of different frequencies

4 结论

本文对一种三角陀螺声子晶体进行了力学建模，得到了结构的动力学方程及色散关系，发现陀螺的存在会对二维弹簧质量声子晶体结构起到一个附加质量的作用，并对质点施加了陀螺力.存在两支色散曲面，通过对陀螺转速的调节，可以控制其间的带隙，从抑制指定频率范围的波.色散曲面等频截线表明，对应于两支曲面的波，当频率较低时，波向各个角度同等传播，当频率较高时，波有六个集中的传播方向，但上下曲面的传播方向角度不同.数值验证结果显示，波在一定频段内受到了明显抑制，频率范围与理论获得的带隙吻合良好，证明了理论获得的带隙确实存在，且精度较高.