搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

各向异性三维非对称双锥五模超材料的能带结构及品质因数

蔡成欣 陈韶赓 王学梅 梁俊燕 王兆宏

各向异性三维非对称双锥五模超材料的能带结构及品质因数

蔡成欣, 陈韶赓, 王学梅, 梁俊燕, 王兆宏
PDF
HTML
导出引用
  • 五模超材料是一类可以解除剪切模量的人工固体微结构, 具有类似流体的性质, 在声波调控中有着潜在的应用. 声学变换作为声波调控的一种重要手段, 在超材料声学器件的设计中被广泛使用. 声学变换的引入会压缩均匀各项同性五模超材料. 因此, 需要研究各向异性对三维非对称双锥五模超材料带隙及品质因数的影响. 本文利用有限元方法, 对各项异性三维非对称双锥五模超材料的能带结构及品质因数进行了研究. 研究结果表明, 三维非对称双锥五模超材料单模区域的相对带宽随着各向异性的增大而减小, 第一带隙的相对带宽增大到123%, 品质因数增大到6.9倍. 本研究可为声学变换在三维非对称双锥五模超材料中的应用提供指导.
      通信作者: 蔡成欣, cxcai2018@haut.edu.cn
    [1]

    Milton G W, Cherkaev A V 1995 J. Eng. Mater. Technol. 117 483

    [2]

    Kadic M, Bückmann T, Stenger N, Thiel M, Wegener M 2012 Appl. Phys. Lett. 100 191901

    [3]

    陈毅, 刘晓宁, 向平, 胡更开 2016 力学进展 46 201609

    Chen Y, Liu X N, Xiang P, Hu G K 2016 Adv. Mech. 46 201609

    [4]

    王兆宏, 蔡成欣, 楚杨阳, 刘广顺 2017 光电工程 44 34

    Wang Z H, Cai C X, Chu Y Y, Liu G S 2017 Opto-Electron. Eng. 44 34

    [5]

    Milton G W, Briane M, Wills J R 2006 New J. Phys. 8 248

    [6]

    Norris A N 2008 Proc. R. Soc. A 464 2411

    [7]

    Scandrett L C, Boisvert J E, Howarth T R 2010 J. Acoust. Soc. Am. 127 2856

    [8]

    Scandrett L C, Boisvert J E, Howarth T R 2011 Wave Motion 48 505

    [9]

    Boisvert J E, Scandrett L C, Howarth T R 2016 J. Acoust. Soc. Am. 139 3404

    [10]

    Schittny R, Bückmann T, Kadic M, Wegener M 2013 Appl. Phys. Lett. 103 231905

    [11]

    Gokhale N H, Cipolla J L, Norris A N 2012 J. Acoust. Soc. Am. 132 4

    [12]

    Kadic M, Buckmann T, Schittny R, Gumbsch P, Wegener M 2014 Phys. Rev. Appl. 2 054007

    [13]

    Cai C X, Wang Z H, Li Q W, Xu Z, Tian X G 2015 J. Phys. D: Appl. Phys. 48 175103

    [14]

    Huang Y, Lu X G, Liang G Y, Xu Z 2016 Phys. Lett. A 380 1334

    [15]

    Wang G, Jin L, Zhang L, Xu Z 2017 AIP Adv. 7 025309

    [16]

    Tian Y, Wei Q, Cheng Y, Xu Z, Liu X J 2015 Appl. Phys. Lett. 107 221906

    [17]

    Sun Z Y, Jia H, Chen Y, Wang Z, Yang J 2018 J. Acoust. Soc. Am. 143 1029

    [18]

    Chen Y, Liu X N, Hu G K 2015 Sci. Rep. 5 15745

    [19]

    Chen J G, Liu J H, Liu X Z 2018 AIP Adv. 8 085024

    [20]

    张向东, 陈虹, 王磊, 赵志高, 赵爱国 2015 物理学报 64 134303

    Zhang X D, Chen H, Wang L, Zhao Z G, Zhao A G 2015 Acta Phys. Sin. 64 134303

    [21]

    陆智淼, 蔡力, 温激鸿, 温熙森 2016 物理学报 65 174301

    Lu Z M, Cai L, Wen J H, Wen X S 2016 Acta Phys. Sin. 65 174301

    [22]

    Chen H Y, Chan C T 2007 Appl. Phys. Lett. 91 183518

    [23]

    Chen H Y, Chan C T 2010 J. Phys. D: Appl. Phys. 43 113001

    [24]

    Cai C X, Han C, Wu J F, Wang Z H, Zhang Q H 2019 J. Phys. D: Appl. Phys. 52 045601

    [25]

    Wang Z H, Cai C X, Li Q W, Li J, Xu Z 2016 J. Appl. Phys. 120 024903

    [26]

    Bückmann T, Schittny R, Thiel M, Kadic M, Milton G W, Wegener M 2014 New J. Phys. 16 033032

  • 图 1  (a)各向同性五模超材料晶胞与(b)−(e)各向异性五模超材料晶胞结构示意图

    Fig. 1.  The unit cell structure of isotropy (a) and (b)−(e) anisotropic pentamode materials.

    图 2  P点沿空间对角线方向偏移0.25倍对角线长度时的能带结构

    Fig. 2.  The band structure of pentamode material with ${O_{\rm{s}}}P/\sqrt 3 a = 0.25$.

    图 3  模型1的第一带隙与单模区域的(a)上下界频率; (b)相对带宽

    Fig. 3.  (a) The upper and lower edges and (b) relative bandwidth of the first phononic band gaps and single mode area of model 1.

    图 7  各向异性对非对称双锥五模超材料品质因数的影响

    Fig. 7.  The influence of anisotropy on the figure of merit of asymmetric double-cone pentamode materials.

    图 4  模型2的第一带隙与单模区域的(a)上下界频率; (b)相对带宽

    Fig. 4.  (a) The upper and lower edges and (b) relative bandwidth of the first phononic band gaps and single mode area of model 2.

    图 5  模型3的第一带隙与单模区域的(a)上下界频率; (b)相对带宽

    Fig. 5.  (a) The upper and lower edges and (b) relative bandwidth of the first phononic band gaps and single mode area of model 3.

    图 6  模型4的第一带隙与单模区域的(a)上下界频率; (b)相对带宽

    Fig. 6.  (a) The upper and lower edges and (b) relative bandwidth of the first phononic band gaps and single mode area of model 4.

  • [1]

    Milton G W, Cherkaev A V 1995 J. Eng. Mater. Technol. 117 483

    [2]

    Kadic M, Bückmann T, Stenger N, Thiel M, Wegener M 2012 Appl. Phys. Lett. 100 191901

    [3]

    陈毅, 刘晓宁, 向平, 胡更开 2016 力学进展 46 201609

    Chen Y, Liu X N, Xiang P, Hu G K 2016 Adv. Mech. 46 201609

    [4]

    王兆宏, 蔡成欣, 楚杨阳, 刘广顺 2017 光电工程 44 34

    Wang Z H, Cai C X, Chu Y Y, Liu G S 2017 Opto-Electron. Eng. 44 34

    [5]

    Milton G W, Briane M, Wills J R 2006 New J. Phys. 8 248

    [6]

    Norris A N 2008 Proc. R. Soc. A 464 2411

    [7]

    Scandrett L C, Boisvert J E, Howarth T R 2010 J. Acoust. Soc. Am. 127 2856

    [8]

    Scandrett L C, Boisvert J E, Howarth T R 2011 Wave Motion 48 505

    [9]

    Boisvert J E, Scandrett L C, Howarth T R 2016 J. Acoust. Soc. Am. 139 3404

    [10]

    Schittny R, Bückmann T, Kadic M, Wegener M 2013 Appl. Phys. Lett. 103 231905

    [11]

    Gokhale N H, Cipolla J L, Norris A N 2012 J. Acoust. Soc. Am. 132 4

    [12]

    Kadic M, Buckmann T, Schittny R, Gumbsch P, Wegener M 2014 Phys. Rev. Appl. 2 054007

    [13]

    Cai C X, Wang Z H, Li Q W, Xu Z, Tian X G 2015 J. Phys. D: Appl. Phys. 48 175103

    [14]

    Huang Y, Lu X G, Liang G Y, Xu Z 2016 Phys. Lett. A 380 1334

    [15]

    Wang G, Jin L, Zhang L, Xu Z 2017 AIP Adv. 7 025309

    [16]

    Tian Y, Wei Q, Cheng Y, Xu Z, Liu X J 2015 Appl. Phys. Lett. 107 221906

    [17]

    Sun Z Y, Jia H, Chen Y, Wang Z, Yang J 2018 J. Acoust. Soc. Am. 143 1029

    [18]

    Chen Y, Liu X N, Hu G K 2015 Sci. Rep. 5 15745

    [19]

    Chen J G, Liu J H, Liu X Z 2018 AIP Adv. 8 085024

    [20]

    张向东, 陈虹, 王磊, 赵志高, 赵爱国 2015 物理学报 64 134303

    Zhang X D, Chen H, Wang L, Zhao Z G, Zhao A G 2015 Acta Phys. Sin. 64 134303

    [21]

    陆智淼, 蔡力, 温激鸿, 温熙森 2016 物理学报 65 174301

    Lu Z M, Cai L, Wen J H, Wen X S 2016 Acta Phys. Sin. 65 174301

    [22]

    Chen H Y, Chan C T 2007 Appl. Phys. Lett. 91 183518

    [23]

    Chen H Y, Chan C T 2010 J. Phys. D: Appl. Phys. 43 113001

    [24]

    Cai C X, Han C, Wu J F, Wang Z H, Zhang Q H 2019 J. Phys. D: Appl. Phys. 52 045601

    [25]

    Wang Z H, Cai C X, Li Q W, Li J, Xu Z 2016 J. Appl. Phys. 120 024903

    [26]

    Bückmann T, Schittny R, Thiel M, Kadic M, Milton G W, Wegener M 2014 New J. Phys. 16 033032

  • [1] 庄飞, 肖三水, 何江平, 何赛灵. 二维正方各向异性碲圆柱光子晶体完全禁带中缺陷模的FDTD计算分析和设计. 物理学报, 2002, 51(9): 2167-2172. doi: 10.7498/aps.51.2167
    [2] 宋顾周, 马继明, 王奎禄, 周鸣. 厚针孔射线成像品质因数的研究. 物理学报, 2012, 61(10): 102902. doi: 10.7498/aps.61.102902
    [3] 周建华, 刘虹遥, 罗海陆, 文双春. 各向异性超常材料中倒退波的传播研究. 物理学报, 2008, 57(12): 7729-7736. doi: 10.7498/aps.57.7729
    [4] 庄飞, 何赛灵, 何江平, 冯尚申. 大带隙的二维各向异性椭圆介质柱光子晶体. 物理学报, 2002, 51(2): 355-361. doi: 10.7498/aps.51.355
    [5] 张利伟, 赵玉环, 王勤, 方恺, 李卫彬, 乔文涛. 各向异性特异材料波导中表面等离子体的共振性质. 物理学报, 2012, 61(6): 068401. doi: 10.7498/aps.61.068401
    [6] 庄 飞, 沈建其. 双轴各向异性负折射率材料光纤中光子波函数几何相位研究. 物理学报, 2005, 54(2): 955-960. doi: 10.7498/aps.54.955
    [7] 李安华, 董生智, 李卫. 烧结Sm2Co17型永磁材料的力学性能及断裂行为的各向异性. 物理学报, 2002, 51(10): 2320-2324. doi: 10.7498/aps.51.2320
    [8] 凌瑞良, 冯进, 冯金福. 三维各向异性耦合谐振子体系的量子化能谱与精确波函数. 物理学报, 2010, 59(12): 8348-8358. doi: 10.7498/aps.59.8348
    [9] 蔡 力, 韩小云, 温熙森. 长波条件下二维声子晶体中的弹性波传播及各向异性. 物理学报, 2008, 57(3): 1746-1752. doi: 10.7498/aps.57.1746
    [10] 王 刚, 温激鸿, 韩小云, 赵宏刚. 二维声子晶体带隙计算中的时域有限差分方法. 物理学报, 2003, 52(8): 1943-1947. doi: 10.7498/aps.52.1943
    [11] 钟会林, 吴福根, 姚立宁. 遗传算法在二维声子晶体带隙优化中的应用. 物理学报, 2006, 55(1): 275-280. doi: 10.7498/aps.55.275
    [12] 齐共金, 杨盛良, 白书欣, 赵 恂. 基于平面波算法的二维声子晶体带结构的研究. 物理学报, 2003, 52(3): 668-671. doi: 10.7498/aps.52.668
    [13] 王 刚, 温激鸿, 刘耀宗, 郁殿龙, 温熙森. 大弹性常数差二维声子晶体带隙计算中的集中质量法. 物理学报, 2005, 54(3): 1247-1252. doi: 10.7498/aps.54.1247
    [14] 卢敏, 黄惠莲, 余冬海, 刘维清, 魏望和. 不同晶面银纳米晶高温熔化的各向异性. 物理学报, 2015, 64(10): 106101. doi: 10.7498/aps.64.106101
    [15] 穆全全, 刘永军, 胡立发, 李大禹, 曹召良, 宣 丽. 光谱型椭偏仪对各向异性液晶层的测量. 物理学报, 2006, 55(3): 1055-1060. doi: 10.7498/aps.55.1055
    [16] 杨宏伟, 袁 洪, 陈如山, 杨 阳. 各向异性磁化等离子体的SO-FDTD算法. 物理学报, 2007, 56(3): 1443-1446. doi: 10.7498/aps.56.1443
    [17] 翁紫梅, 陈 浩. 单离子各向异性影响下的一维铁磁链中的孤子. 物理学报, 2007, 56(4): 1911-1918. doi: 10.7498/aps.56.1911
    [18] 孟繁义, 吴 群, 傅佳辉, 顾学迈, 李乐伟. 三维各向异性超常媒质交错结构的亚波长谐振特性研究. 物理学报, 2008, 57(10): 6213-6220. doi: 10.7498/aps.57.6213
    [19] 王志军, 王锦程, 杨根仓. 各向异性作用下合金定向凝固界面稳定性的渐近分析. 物理学报, 2008, 57(2): 1246-1253. doi: 10.7498/aps.57.1246
    [20] 孟繁义, 吴 群, 傅佳辉, 杨国辉. 各向异性超常媒质矩形波导的导波特性研究. 物理学报, 2008, 57(9): 5476-5484. doi: 10.7498/aps.57.5476
  • 引用本文:
    Citation:
计量
  • 文章访问数:  679
  • PDF下载量:  42
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-03-12
  • 修回日期:  2020-04-28
  • 上网日期:  2020-05-14
  • 刊出日期:  2020-07-01

各向异性三维非对称双锥五模超材料的能带结构及品质因数

  • 1. 河南工业大学粮食信息处理与控制教育部重点实验室, 郑州 450001
  • 2. 河南工业大学河南省粮食光电探测与控制重点实验室, 郑州 450001
  • 3. 河南工业大学信息科学与工程学院,郑州 450001
  • 4. 西安交通大学电子科学与工程学院, 物理电子与器件教育部重点实验室, 西安 710049
  • 通信作者: 蔡成欣, cxcai2018@haut.edu.cn

摘要: 五模超材料是一类可以解除剪切模量的人工固体微结构, 具有类似流体的性质, 在声波调控中有着潜在的应用. 声学变换作为声波调控的一种重要手段, 在超材料声学器件的设计中被广泛使用. 声学变换的引入会压缩均匀各项同性五模超材料. 因此, 需要研究各向异性对三维非对称双锥五模超材料带隙及品质因数的影响. 本文利用有限元方法, 对各项异性三维非对称双锥五模超材料的能带结构及品质因数进行了研究. 研究结果表明, 三维非对称双锥五模超材料单模区域的相对带宽随着各向异性的增大而减小, 第一带隙的相对带宽增大到123%, 品质因数增大到6.9倍. 本研究可为声学变换在三维非对称双锥五模超材料中的应用提供指导.

English Abstract

    • 五模超材料是Milton与Cherkaey[1]提出的一类具有固体形态和“流体”特性的人工微结构极值材料. 通过解除压缩波和剪切波的耦合, 使其周期性结构很难被压缩, 但却很容易发生剪切形变[2]. 由这样的结构单元周期排列构成的声学超材料, 可以实现类似“流体”的性质. 五模超材料可调的模量各向异性、较小的填充率、宽频等特点赋予五模超材料优越的水下声波调控能力[3,4].

      2006年, Milton等[5]通过对传统弹性动力学方程在曲线变换下的变化规律, 提出了五模超材料用于弹性波隐身的可行性. 2008年, Norris[6]系统分析了惯性声隐身与五模超材料声隐身衣, 提出了变换声学理论在五模超材料中应用的可行性. 2010年, Scandrett等[7,8]和Boisvert等[9]利用分层化的五模超材料提出了声隐身衣的设计方法. 理想的五模超材料是由窄直径很小的双锥单元连接构成的面心立方结构, 结构稳定性较差. 为了能够得到稳定性相对较好的三维五模超材料, Kadic等[2]和Schittny等[10]加入实心小球来替代点连接, 并利用侵入式激光直写光刻技术与3D打印技术制作出了微米及毫米量级的三维实物结构. 此后, 研究人员开始对五模超材料的声学与力学特性进行研究[11-15], 设计出声学超表面[16]、声学波导[17]与声学隐身衣[18,19]等多种二维五模超材料声波调控器件.

      基于变换声学理论, 五模超材料声学隐身衣能够避免惯性声学隐身衣壳层内部质量奇异点的问题, 这为五模声学隐身衣的设计和制作提供了基础[6]. 由于在二维情况下, 声波方程可以与麦克斯韦方程组一一对应, 因此对五模超材料的研究多集中在二维结构器件的设计与实验. 2015年, 张向东等[20]采用分层近似的方法, 对圆柱形分层五模超材料声学隐身衣的散射声压场进行了理论分析, 并通过数值计算验证了理论推导的正确性. 2016年, 陆智淼等[21]对基于五模超材料的声学变换方程进行了研究, 分析了影响五模超材料声学隐身衣性能的因素及规律, 结果表明合适的坐标变换方程能够有效地改善隐身效果. 同年, 陈毅等基于准对称映射梯度算法获得了任意形状的二维五模声学隐身衣. 此外, Chen等[19]还提出基于五模超材料的隐身地斗篷.

      由于三维五模超材料结构的复杂性, 目前对三维五模超材料的研究还比较缺乏. Chen和Chan[22,23]利用直流电导方程在坐标变换下的不变性, 将声波方程和直流电导方程建立了一一对应的关系, 首次导出一般的三维变换声学方程, 给出了声学坐标变换前后三维空间等效参数需要满足的对应关系, 为三维五模超材料声学隐身衣的设计提供了可行的研究方法和理论基础. 在设计三维五模超材料隐身衣时, 除了需要对外部声波进行调控之外, 还需要考虑对隐身衣内部声波的抑制[13]. 声学坐标变换下, 五模超材料的周期性结构会被不同程度地压缩, 引起双锥单元的结构参数发生变化, 并对其能带结构、品质因数和填充率产生作用, 进而影响五模超材料对其内部声波的控制[24,25]. 因此, 本文基于有限元方法研究各向异性对三维非对称双锥五模超材料能带特性及品质因数的影响, 为变换声学在三维五模超材料声学器件的设计及应用提供参考.

    • 三维各向同性非对称双锥五模超材料晶胞结构如图1(a)所示, 由十六个非对称双锥相交构成面心立方晶格. 非对称双锥的宽直径为D, 窄直径分别为d1d2, 高度为H, 晶格常数为a. 构造三维各向异性五模超材料有很多方式, 其中最直接的方法就是: 选取一个原胞的窄直径交点P为参考点, 分别沿四个不同方向移动交点P来构造三维各向异性五模超材料, 对应的各向异性晶胞结构如图1(b)图1(e)所示, 分别命名为模型1(X轴)、模型2(Y轴)、模型3(Z轴)、模型4(体对角线). 当点P (0.25a, 0.25a, 0.25a)如图1(a)所示时, 三维五模超材料为各向同性的, 其各向异性程度最低; 当P点移动时, 距离点(0.25a, 0.25a, 0.25a)的位置越远, 其各向异性程度就越高.

      图  1  (a)各向同性五模超材料晶胞与(b)−(e)各向异性五模超材料晶胞结构示意图

      Figure 1.  The unit cell structure of isotropy (a) and (b)−(e) anisotropic pentamode materials.

    • 各向异性的引入会对组成三维非对称双锥五模超材料原胞的基元结构参数H产生影响, 使得原胞内原本相同的四个基元变得不同, 从而影响原胞内部的本征振动形态. 因此, 各向异性会影响三维非对称双锥五模超材料的带隙特性、单模区域、压缩模量及剪切模量. 为了系统地研究各向异性对三维非对称双锥五模超材料的能带结构及五模特性的影响, 利用有限元仿真软件Comsol Multiphysics, 在布洛赫边界条件下对各向异性原胞进行数值计算. 原胞结构参数为 a = 37.3 mm, H = 16.15 mm, D = 3 mm, d1 = 0.6 mm, d2 = 0.3 mm. 组成材料的质量密度为1190 kg/m3、泊松比为0.4、杨氏模量 3 GPa.

      图2P点沿空间对角线方向偏移0.25倍对角线长度时所对应的能带结构图, 其中横纵坐标分别对应整个简约布里渊区边界与频率. 第一带隙(黑色区域)的下边界(fl)与上边界(fu)频率分别为 9.017与10.466 kHz, 对应的相对带宽$\left(\dfrac{{\Delta \omega }}{{{\omega _{\rm{g}}}}}\right. = \left.\dfrac{{{f_{\rm{u}}} - {f_{\rm{l}}}}}{{{{({f_{\rm{u}}} + {f_{\rm{l}}})} / 2}}}\right)$为0.149, 频率落在带隙范围内的压缩波和剪切波的传播均被抑制. 单模区域(灰色区域)的下边界与上边界频率分别为 0.457与2.333 kHz, 相对带宽为1.345. 在单模区域频率范围内, 压缩波与剪切波将被解耦合, 即只有压缩波可以传播, 剪切波将被抑制.

      图  2  P点沿空间对角线方向偏移0.25倍对角线长度时的能带结构

      Figure 2.  The band structure of pentamode material with ${O_{\rm{s}}}P/\sqrt 3 a = 0.25$.

      理想的三维五模超材料是难压缩、易形变的固体人工微结构, 其体积弹性模量B与剪切模量G具有较大比值. 为了描述单模区域对应压缩波与剪切波的解耦合能力, 将体积弹性模量B与剪切模量G的比值定义为品质因数(figure of merit). 品质因数越大, 五模超材料解除压缩波与剪切波耦合的能力就越好. 对于各向异性三维非对称双锥五模超材料, 弹性模量B与剪切模量G的值可由以下式子得到[26]:

      ${C_{11}} = 2\rho {\left(\upsilon _{110}^{T,xy}\right)^2} + {C_{12}}, $

      ${C_{12}} = \rho {\left(\upsilon _{110}^{L,xy}\right)^2} - {C_{44}} - \rho {\left(\upsilon _{110}^{T,xy}\right)^2}, $

      ${C_{44}} = \rho {\left(\upsilon _{110}^{T,z}\right)^2}, $

      $G = {C_{44}}, $

      $B = \left({C_{11}} + 2{C_{12}}\right)/3, $

      式中${C_{11}}$, ${C_{12}}$${C_{44}}$为三维五模超材料弹性矩阵的3个独立弹性系数, $\upsilon _{110}^{T, xy}$, $\upsilon _{110}^{L, xy}$$\upsilon _{110}^{T, z}$分别为压缩波与剪切波沿[110] ($\varGamma K$)方向的相速度, $\rho = f{\rho _0}$是三维五模超材料的等效密度, $f$${\rho _0}$分别为三维五模超材料的填充率与组成材料的质量密度. 因此, 图2中三维五模超材料的品质因数约为251.

    • 基于有限元方法, 对于四种类型的三维各向异性非对称双锥五模超材料的原胞进行能带结构数值计算, 并对其第一带隙、单模区域、品质因数的变化规律进行了系统研究, 结果如图3图7所示.

      图  3  模型1的第一带隙与单模区域的(a)上下界频率; (b)相对带宽

      Figure 3.  (a) The upper and lower edges and (b) relative bandwidth of the first phononic band gaps and single mode area of model 1.

      图  7  各向异性对非对称双锥五模超材料品质因数的影响

      Figure 7.  The influence of anisotropy on the figure of merit of asymmetric double-cone pentamode materials.

      对于模型1, 当P点沿X轴方向移动时, 其第一带隙与单模区域随移动距离的变化如图3所示. 当${O_x}P/a = 0.25$时, 五模超材料的原胞为各向同性的三维结构, 这时P点的各向异性程度最小; 当${O_x}P/a = 0.05$或0.45时, 各向异性程度最大. 随着各向异性程度的增加, 第一带隙的下界频率从8.845 kHz上升到11.753 kHz, 对应的上界频率从9.538 kHz上升到14.846 kHz, 其相对带宽从0.148增加到0.233, 大约增加了57.4%. 单模区域的下界频率在0.457与0.601 kHz之间变化, 对应的上界频率在1.938与2.333 kHz之间变化, 其相对带宽从1.344减小到1.053, 大约减小了21.6%.

      对于模型2, 当P点沿Y轴方向移动时, 各向异性五模超材料的第一带隙与单模区域随移动距离的变化如图4所示. 随着各向异性程度的增加, 第一带隙的下界频率在8.017—11.27 kHz之间变化, 对应的上界频率在8.421—14.852 kHz之间变化, 其相对带宽从0.147增加到0.328, 大约增加了123%. 单模区域的下界频率从0.457 kHz上升到0.632 kHz, 对应的上界频率从2.333 kHz减小到1.448 kHz, 其相对带宽从1.344减小到0.881, 大约减小了34.4%. 单模区域相对带宽在各向异性程度最大处明显下降的主要原因是由于第一带隙相对带宽增大, 单模区域的频率范围受到压缩.

      图  4  模型2的第一带隙与单模区域的(a)上下界频率; (b)相对带宽

      Figure 4.  (a) The upper and lower edges and (b) relative bandwidth of the first phononic band gaps and single mode area of model 2.

      对于模型3, 当P点沿Z轴方向移动时, 各向异性五模超材料的第一带隙与单模区域随移动距离的变化如图5所示. 随着各向异性程度的增加, 第一带隙的下界频率从8.839 kHz上升到11.764 kHz, 对应的上界频率从9.527 kHz上升到14.847 kHz, 其相对带宽从0.148增加到0.236, 大约增加了59.4%. 单模区域的下界频率在0.438—0.562 kHz之间变化, 对应的上界频率在1.941—2.333 kHz之间变化, 其相对带宽从1.344减小到1.101, 大约减小了18.1%.

      图  5  模型3的第一带隙与单模区域的(a)上下界频率; (b)相对带宽

      Figure 5.  (a) The upper and lower edges and (b) relative bandwidth of the first phononic band gaps and single mode area of model 3.

      对于模型4, 当P点沿对角线方向移动时, 各向异性五模超材料的第一带隙与单模区域随移动距离的变化如图6所示. 随着各向异性程度的增加, 第一带隙的下界频率在4.701—10.034 kHz之间变化, 对应的上界频率在4.992—17.794 kHz之间变化, 其相对带宽从0.147增加到0.256, 约增加了74.1%. 单模区域的下界频率在0.322—0.923 kHz之间变化, 对应的上界频率在1.217—2.333 kHz之间变化, 其相对带宽从1.344减小到0.449, 大约减小了66.6%. 在${O_{\rm{s}}}P/\sqrt 3 a = 0.05$处, 第一带隙的拓宽对单模区域的上界产生了较强的影响, 使单模区域的上界频率下降, 同时其下界频率上升, 从而单模区域相对带宽在此处下降比较明显.

      图  6  模型4的第一带隙与单模区域的(a)上下界频率; (b)相对带宽

      Figure 6.  (a) The upper and lower edges and (b) relative bandwidth of the first phononic band gaps and single mode area of model 4.

      “流体”特性是非对称双锥五模超材料的主要特性之一, 而品质因数是衡量这种性质的重要参数. 各向异性对五模超材料品质因数的影响如图7所示, 由于三维非对称双锥五模超材料的各向异性变化是沿着四个不同的方向, 因此对原胞内四个双锥基元高度的影响是不同的, 即不同方向上的各向异性对三维非对称双锥五模超材料能带结构的影响也不相同. 当P点沿X轴与Y轴方向变化时, 其品质因数最大值约为251; 当P点沿Z轴变化时, 其品质因数最大值约为1730; 当P点沿对角线方向变化时, 其品质因数最大值约为327. 相对于各向同性三维非对称双锥五模超材料, 各向异性三维非对称双锥五模超材料的品质因数可以增加近6.9倍.

    • 通过调节基元窄直径连接点P的位置, 构建了四种各向异性三维非对称双锥五模超材料模型. 由于不同方向上的各向异性变化对非对称双锥结构基元参数的影响不尽相同, 其能带结构与品质因数也会受到不同程度的影响, 本文给出了各向异性程度与三维非对称双锥五模超材料带隙特性、单模区域与品质因数的关系. 与各向同性三维非对称双锥五模超材料相比, 各向异性三维非对称双锥五模超材料第一带隙的相对带宽可被有效地拓宽, 品质因数可增大6.9倍. 研究结果可为基于三维非对称双锥五模超材料的声学器件设计提供参考. 由于各向异性的引入, 使得三维非对称双锥五模超材料的结构更为复杂, 样件加工的需求进一步提高, 整体稳定性也有所降低, 因此对具有较高稳定性且易加工的新型三维五模超材料原胞结构还需进一步研究探索.

参考文献 (26)

目录

    /

    返回文章
    返回