上海交通大学学报

上海交通大学学报（自然版）, 2021, 55(11): 1453-1458 doi: 10.16183/j.cnki.jsjtu.2020.242

韧带型声子晶体狄拉克锥特性研究

陈炉云1,2, 王健,1,2, 崔益烽3, 孔慧1,2

1.上海交通大学 海洋工程国家重点实验室, 上海 200240

2.上海交通大学 海洋智能装备与系统教育部重点实验室, 上海 200240

3.陆军装备部驻南京地区军事代表局驻上海地区第二军事代表室, 上海 200082

Dirac Cone Characteristics of Hexachiral Phononic Crystal

CHEN Luyun1,2, WANG Jian,1,2, CUI Yifeng3, KONG Hui1,2

1. State Key Laboratory of Ocean Engineering, Shanghai Jiao Tong University, Shanghai 200240, China

2. Key Laboratory of Marine Intelligent Equipment and System of the Ministry of Education, Shanghai Jiao Tong University, Shanghai 200240, China

3. Shanghai Second Militaryx Representative Office, Nanjing Military Representative Office of the Ministry of Army Equipment, Shanghai 200082, China

通讯作者: 王 健,男,助理研究员,电话(Tel.):021-34207165;E-mail:nsms_sjtu@sjtu.edu.cn.

责任编辑: 孙伟

收稿日期: 2020-07-29  

Received: 2020-07-29  

作者简介 About authors

陈炉云(1975-),男,湖南省郴州市人,博士生,主要从事结构减振降噪研究.

摘要

声子晶体的能带结构特性是评估声学超构材料减振降噪性能的重要指标.以二维六韧带手性声子晶体为例,对其能带结构和狄拉克锥特性进行数值分析,在声子晶体布里渊区中心处获得四重偶然简并狄拉克点.通过调整韧带的设计参数,实现对双狄拉克锥的打开,形成新的定向带隙.研究韧带参数与定向带隙宽度的关系,并进一步讨论能带结构反转问题.研究成果可为手性声子晶体在弹性波操控和声学拓扑绝缘体的应用提供技术支持.

关键词： 声子晶体; 能带结构; 狄拉克锥; 线性频散; 方向带隙

Abstract

The band structure properties of phononic crystal is important to evaluate the vibration and noise reduction of acoustic metamaterials. Taking the 2D hexachiral phononic crystal as an example, the band structure and Dirac cone properties were investigated by numerical analysis, and the four-fold accidental degenerate Dirac point was obtained in the center of Brillouin zone. By adjusting the design parameters of ligament structure, a double Dirac cone was broken and a novel directional band gap was formed. The influence of geometric parameters on the directional band gaps width was investigated, and the band structure inversion problem was further discussed. This research can provide support for the application of hexachiral phononic crystal in elastic wave manipulation and acoustic topological insulator.

Keywords： phononic crystal; band structure; Dirac cone; linear dispersion; directional band gap

PDF (3251KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

陈炉云, 王健, 崔益烽, 孔慧. 韧带型声子晶体狄拉克锥特性研究[J]. 上海交通大学学报（自然版）, 2021, 55(11): 1453-1458 doi:10.16183/j.cnki.jsjtu.2020.242

CHEN Luyun, WANG Jian, CUI Yifeng, KONG Hui. Dirac Cone Characteristics of Hexachiral Phononic Crystal[J]. Journal of shanghai Jiaotong University, 2021, 55(11): 1453-1458 doi:10.16183/j.cnki.jsjtu.2020.242

随着科技发展,声学超材料在结构的减振降噪方面得到广泛应用[1,2].与传统材料相比,通过微结构功能胞元设计,声学超材料具有负等效质量密度、负等效体积模量和负剪切模量等特性[3,4].声子晶体中周期性排列的宏观介质材料会对声波产生散射,令其形成具有周期性的能带结构[5,6].各能带之间的缝隙称为声子带隙,影响带隙参数的因素包括结构参数、材料参数、多相材料和界面参数等.改变声子晶体的设计参数可以实现对带隙结构的调整和操控,提高声子晶体减振降噪的效果[7,8,9,10].在某些声子晶体的简约布里渊区边界上,存在两条或多条能带线性地相交于一点的频散现象,此交点称为狄拉克点 [11,12].在狄拉克点处将出现赝扩散传输、边缘态、震颤和声拓扑绝缘体等物理现象 [13,14].

声子晶体的胞元形式多样,其中手性韧带型声子晶体结构具有典型的负泊松比特性,因此其能带结构具有一定的特殊性[15,16].目前,韧带型蜂窝结构的研究成果多集中在因其负泊松比而导致的抗冲击、减振等特性中[17,18,19,20].同时,因手性特性,手性韧带型声子晶体不具有C6晶体完美的对称性,这对声子晶体的声传输、声拓扑等问题有重要影响[21,22].而针对手性韧带型声子晶体声拓扑问题的研究较少,对于声子晶体频散中出现的狄拉克锥问题的研究还处于探索阶段[23,24,25,26].因此,开展手性六韧带型声子晶体的带隙特性、拓扑绝缘体特性的研究可以为声隐形、完美透镜和单向传输等物理现象提供技术支撑.

本文提出一种基于二维手性六韧带蜂窝晶格型声子晶体,通过分析其能带结构和狄拉克频散特性,获得声子晶体的四重偶然简并狄拉克点.通过调整声子晶体散射体结构参数,打破四重偶然简并并获得新的方向带隙,并对能带反转特性进行分析.研究成果可以为新型声学超构材料的设计和工程应用提供技术支持.

1 二维声子晶体基本理论

1.1 Bloch 定理

声子晶体是由相同晶胞单元在空间内无限周期性延拓而成,其周期结构的最小周期尺度为晶格常数.声子晶体的结构周期性由组成声子晶体的介质材料的弹性常数及其密度的周期性体现,即

T^Rf(r)=f(r+R)=f(r)

式中:r和R为位矢; T^R为r至r+R的点阵平移算子.

声子晶体结构的空间周期性和平移对称性,使其本征场和本征频率也具有对称性.弹性波在声子晶体中的传播符合Bloch定理,其位移场满足如下周期形式:

u(r)=eikruk(r)

式中:u(r)为与声子晶体具有相同周期性的矢量函数;eikr为相位因子;k为波矢,其仅在第一布里渊区内取值,表示当平移r时波函数只增加一个eikr.

1.2 声子晶体能带结构

利用波矢扫掠声子晶体不可约布里渊区可以得到能带结构.以波矢为横坐标,以本征值或与其相应的量(如频率)为纵坐标作图,可以得到声子晶体能带结构或频散关系.在某些频率范围内,声子晶体能带结构中不存在任何本征模式且不支持任何形式的弹性波的传播.声子晶体带隙的实质是弹性波能量随传播距离的增加而迅速呈指数衰减.带隙宽度、对弹性波传播的衰减能力等是评估声子晶体材料物理特性的重要指标.

在声子晶体能带结构中,在确定频率范围内包括q条能带和j条带隙.设在第n条能带和第n+1条能带之间存在带隙,且该带隙对应的下边界能带的最大频率值和上边界能带的最小频率值分别为max(ωn(k))和min(ωn+1(k)),并满足min(ωn+1(k))>max(ωn(k)).该带隙宽度和中心频率分别为

Δωn=min(ωn+1(k))-max(ωn(k))ωc=(min(ωn+1(k))+max(ωn(k)))/2

通过调整声子晶体的设计参数和材料组分,可以增大或减小Δωn.

2 六韧带手性蜂窝材料

六韧带手性蜂窝材料的胞元为正六边形晶格结构,单个胞元包括圆柱节点环以及与其相切连接的6条韧带,每条韧带又循环连接另一个节点环并与之相切.圆柱节点环的内部可以填充与韧带不同的组分材料,其具有明显的负泊松比特性.

声子晶体胞元示意如图1所示,其几何参数包括:韧带长度(L)、两圆柱节点环之间的距离即晶格常数(a)、节点环内径(r0)、节点环厚度(tc)、韧带厚度(t1)、相邻节点中心连线之间的夹角(β)、中心连线与韧带之间的夹角(θ)等. 图1(a)为六韧带手性蜂窝材料的二维结构图,e1和e2为声子晶体的晶格矢; 图1(b)为声子晶体胞元;图1(c)阴影区为声子晶体简约布里渊区.在图1(c)中,3个高对称点分别为Γ(0,0)、M 0,23π3a和K 2π3a,23π3a,实线为不可约布里渊区边界.

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   六韧带手性声子晶体及其布里渊区

Fig.1   Diagrams of hexachiral phononic crystal and irreducible Brillouin area


用晶格矢ei定义声子晶体结构的空间周期性.二维声子晶体的位置可以根据晶格矢的线性组合n1e1+n2e2确定,则

e1=[acosθasinθ]T=3a2a2Te2=[-acosθasinθ]T=-3a2a2T

3 数值计算

3.1 模型概述

以二维三组元六韧带手性声子晶体为例,进行固-液型声子晶体能带结构的数值计算.在COMSOL Multiphysics 5.4有限元软件平台上,依次通过确定晶格常数,实施几何建模,进行有限元划分,选择声学材料,定义周期性边界条件和选取物理场后,对特征频率进行路径扫描,可以获得声子晶体的能带结构.

声子晶体由圆柱体、韧带散射体和基体组成,且a=8 mm.声子晶体韧带和钢柱体的外包覆物材料均为环氧树脂,圆柱体材料为钢,基体材料为水.图2为二维六韧带手性声子晶体胞元的外形及其网格模型.各材料力学参数包括: 水密度为 1000 kg/m3, 声速为1480 m/s;环氧树脂密度为 1300 kg/m3, 声速为 2830 m/s;钢密度为 7850 kg/m3, 声速为 5760 m/s.

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   二维声子晶体模型和网格分布

Fig.2   Finite element model and grid distribution of 2D phononic crystal


3.2 声子晶体带隙结构计算

确定声子晶体结构参数后,在布里渊区对声子晶体特征频率进行参数化扫描可以计算声子晶体的能带结构,扫描路径为K→Γ→M→K.在数值计算中,对前8条能带进行分析.

调整声子晶体散射体结构参数r0、tc和tl,可以得到具有四重偶然简并态的双狄拉克锥,并将其作为初模型.初模型声子晶体参数:r0=1.31 mm,tc=0.35 mm, tl=0.45 mm.

图3(a)为具有四重偶然简并双狄拉克锥手性蜂窝晶格声子晶体的频散关系.声子晶体的第3条至第6条能带在频率f=218795.6 Hz 处出现闭合,形成四重偶然简并狄拉克点(点D),在布里渊区中心表现出双狄拉克锥特征.图3(b)为狄拉克点附近能带结构放大图.在四重简并点附近,频散关系呈线性,4条能带近似合并为2条线性能带.

在第1条和第2条能带之间以及第6条和第7条能带之间分别存在1条完全带隙,带隙频率参数如表1所示.

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   初模型声子晶体带隙曲线图

Fig.3   Band gap diagrams of phononic crystal for initial model


表1   声子晶体带隙频率

Tab.1  Frequency of phononic crystal of band gapsHz

模型带隙max(ωn(k))min(ωn+1(k))Δωnωc
初模型第1条71218.6152740.881522.2111979.7
第2条282668.7305206.522537.8293937.6
模型A第1条71025.2152424.081398.8111724.6
第2条282020.2301912.619892.4291966.4
模型B第1条71422.9153072.681649.7112247.8
第2条283371.1308565.725194.6295968.4

新窗口打开| 下载CSV


3.3 四重偶然简并态的打开

针对图3所示的双狄拉克锥问题,将tl定义为设计变量,分析韧带厚度参数对能带结构带隙变化的影响.以初模型的设计参数为基础,固定声子晶体晶格常数不变,增大或减小韧带厚度.数值计算表明,通过调整胞元韧带厚度可以打开双狄拉克锥.

设模型A和B的韧带厚度分别为tlA=0.4 mm和tlB=0.5 mm.分别计算模型A和B两种声子晶体的能带结构,如图4所示.其中,EAL、EAU和EBL、EBU分别为模型A和B的二重简并点.

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4   模型A和B声子晶体带隙曲线图

Fig.4   Band gap diagrams of phononic crystal for models A and B


由图4可知,增大或减小韧带厚度,能带结构中的四重简并态消失.同时,随着四重简并态的打开,第1条和第2条带隙的带隙宽度以及带隙中心频率均发生变化(见表1).在布里渊区中心,四重简并态的双狄拉克锥打开为2条二重简并能态,并形成1条方向带隙.

在模型A和B中,打开的方向带隙宽度(fv)分别为 1557.0 Hz 和 1636.8 Hz.取tl∈[0.35 mm,0.55 mm],步长 Δtl=0.01 mm.对比韧带厚度与方向带隙宽度的关系,如图5所示.

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5   韧带厚度与方向带隙宽度关系图

Fig.5   t1 versus fv


由图5可知,在tl=0.45 mm处,由于出现双狄拉克锥,所以fv=0;随着韧带厚度的增大或减小,fv近似呈线性变化,且韧带厚度的增大对fv的变化更为敏感.

4 能带反转

声子晶体设计参数的调整会破坏初模型声子晶体中的四重偶然简并特性.模型A和B声子晶体在布里渊区中心处的2个双重简并态的声压场分布如图6所示,图中数据为正、负声压相对值.

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6   声子晶体在布里渊区中心的本征模式声压场分布

Fig.6   Distributions of pressure field for eigenmodes in the center of Brillouin zone


图中,在布里渊区中心出现两种类似于电子p和d轨道的双重简并态,且均具有一对偶极模式和一对四极模式.在模型A中,p态位于低频带而d态位于高频带;在模型B中,p态位于高频带而d态位于低频带.可知,从模型A到B,声子晶体的韧带厚度增加并超过初模型的参数值, 实现了模型A与B声子晶体之间的能带反转,使声子晶体从拓扑平庸相转变为拓扑非平庸相.结合图3~5可知,初模型的韧带厚度值为一个突变点,即能带反转的临界点.这表明,狄拉克点所对应的结构设计参数对研究声子晶体能带特性具有重要意义.此外,由于六韧带声子晶体不具有完全的镜像对称性以及C6晶体对称性,所以图6中的偶极模式(px和py)和四极模式(dxy和 dx2-y2)不再关于x轴和y轴完美对称.

5 结语

以六韧带手性蜂窝晶格声子晶体为研究对象,建立求解声子晶体能带结构的数学模型.通过有限元方法计算二维声子晶体的带隙,获得蜂窝晶格声子晶体的四重偶然简并狄拉克点.固定蜂窝晶格声子晶体的晶格常数不变,调节晶体胞元的结构参数实现狄拉克点的打开并产生能带反转,验证设计参数对手性韧带声子晶体声拓扑相变的影响规律.研究成果为韧带型声子晶体在低频减振降噪中的应用提供理论依据.

参考文献

吴九汇, 马富银, 张思文, 等.

声学超材料在低频减振降噪中的应用评述

[J]. 机械工程学报, 2016, 52(13):68-78.

[本文引用: 1]

WU Jiuhui, MA Fuyin, ZHANG Siwen, et al.

Application of acoustic metamaterials in low-frequency vibration and noise reduction

[J]. Journal of Mecha-nical Engineering, 2016, 52(13):68-78.

[本文引用: 1]

阮居祺, 卢明辉, 陈延峰, 等.

基于弹性力学的超构材料

[J]. 中国科学: 技术科学, 2014, 44(12):1261-1270.

[本文引用: 1]

RUAN Juqi, LU Minghui, CHEN Yanfeng, et al.

Metamaterial based on elastic mechanics

[J]. Science China Technological Sciences, 2014, 44(12):1261-1270.

[本文引用: 1]

田源, 葛浩, 卢明辉, 等.

声学超构材料及其物理效应的研究进展

[J]. 物理学报, 2019, 68(19):7-18.

[本文引用: 1]

TIAN Yuan, GE Hao, LU Minghui, et al.

Research advances in acoustic metamaterials

[J]. Acta Physica Sinica, 2019, 68(19):7-18.

[本文引用: 1]

BARAVELLI E, RUZZENE M.

Internally resonat-ing lattices for bandgap generation and low-frequency vibration control

[J]. Journal of Sound and Vibration, 2013, 332(25):6562-6579.

DOI:10.1016/j.jsv.2013.08.014      URL     [本文引用: 1]

LIU Z Y, ZHANG X X, MAO Y W, et al.

Locally resonant sonic materials

[J]. Physica B: Condensed Matter, 2000, 338(1/2/3/4):201-205.

DOI:10.1016/S0921-4526(03)00487-3      URL     [本文引用: 1]

CUMMER S A, CHRISTENSEN J, ALÙ A.

Controlling sound with acoustic metamaterials

[J]. Nature Reviews Materials, 2016, 1(3):16001.

DOI:10.1038/natrevmats.2016.1      URL     [本文引用: 1]

JAMES R, WOODLEY S M, DYER C M, et al.

Sonic bands, bandgaps, and defect states in layered structures—Theory and experiment

[J]. The Journal of the Acoustical Society of America, 1995, 97(4):2041-2047.

DOI:10.1121/1.411995      URL     [本文引用: 1]

梁孝东, 缪林昌, 尤佺, 等.

局域共振二维声子晶体的低频带隙特性研究

[J]. 人工晶体学报, 2019, 48(7):1225-1232.

[本文引用: 1]

LIANG Xiaodong, MIAO Linchang, YOU Quan, et al.

Low-frequency band gap characteristics of locally resonant two-dimensional phononic crystal

[J]. Journal of Synthetic Crystals, 2019, 48(7):1225-1232.

[本文引用: 1]

李妍.

光子晶体和声子晶体中由偶然简并所导致的 Dirac 锥形色散关系研究

[D]. 广州: 华南理工大学, 2015.

[本文引用: 1]

LI Yan.

Research on the Dirac-cone dispersions induced by accidental degeneracy in photonic and phononic crystals

[D]. Guangzhou: South China University of Technology, 2015.

[本文引用: 1]

XIAO B, LAI K F, YU Y, et al.

Exciting reflectionless unidirectional edge modes in a reciprocal photonic topological insulator medium

[J]. Physical Review B, 2016, 94(19):195427.

DOI:10.1103/PhysRevB.94.195427      URL     [本文引用: 1]

范海燕, 夏百战.

三维声学超材料的高阶拓扑态

[J]. 科学通报, 2020, 65(15):1411-1419.

[本文引用: 1]

FAN Haiyan, XIA Baizhan.

Higher-order topological states in a three-dimensional acoustic metamaterial

[J]. Chinese Science Bulletin, 2020, 65(15):1411-1419.

[本文引用: 1]

王一鹤, 张志旺, 程营, 等.

声子晶体中的表面声波赝自旋模式和拓扑保护声传输

[J]. 物理学报, 2019, 68(22):264-271.

[本文引用: 1]

WANG Yihe, ZHANG Zhiwang, CHENG Ying, et al.

Pseudospin modes of surface acoustic wave and topologically protected sound transmission in phono-nic crystal

[J]. Acta Physica Sinica, 2019, 68(22):264-271.

[本文引用: 1]

PENG Y G, QIN C Z, ZHAO D G, et al.

Experimental demonstration of anomalous Floquet topological insulator for sound

[J]. Nature Communications, 2016, 7:13368.

DOI:10.1038/ncomms13368      URL     [本文引用: 1]

HE C, NI X, GE H, et al.

Acoustic topological insulator and robust one-way sound transport

[J]. Nature Physics, 2016, 12(12):1124-1129.

DOI:10.1038/nphys3867      URL     [本文引用: 1]

裴东亮, 杨洮, 陈猛, 等.

基于复合蜂窝结构的宽带周期与非周期声拓扑绝缘体

[J]. 物理学报, 2020, 69(2):153-161.

[本文引用: 1]

PEI Dongliang, YANG Tao, CHEN Meng, et al.

Broadband periodic and aperiodic acoustic topological insulator based on composite honeycomb structure

[J]. Acta Physica Sinica, 2020, 69(2):153-161.

[本文引用: 1]

付子义, 王晨旭, 王立国, 等.

基于COMSOL的声子晶体带结构计算新方法

[J]. 软件, 2018, 39(12):6-9.

[本文引用: 1]

FU Ziyi, WANG Chenxu, WANG Liguo, et al.

A new method for computation of phononic crystals band structure by COMSOL

[J]. Computer Engineering & Software, 2018, 39(12):6-9.

[本文引用: 1]

SPADONI A, RUZZENE M, GONELLA S, et al.

Phononic properties of hexagonal chiral lattices

[J]. Wave Motion, 2009, 46(7):435-450.

DOI:10.1016/j.wavemoti.2009.04.002      URL     [本文引用: 1]

HU L L, LUO Z R, YIN Q Y.

Negative Poisson’s ratio effect of re-entrant anti-trichiral honeycombs under large deformation

[J]. Thin-Walled Structures, 2019, 141:283-292.

DOI:10.1016/j.tws.2019.04.032      URL     [本文引用: 1]

MOUSANEZHAD D, HAGHPANAH B, GHOSH R, et al.

Elastic properties of chiral, anti-chiral, and hierarchical honeycombs: A simple energy-based approach

[J]. Theoretical and Applied Mechanics Letters, 2016, 6(2):81-96.

DOI:10.1016/j.taml.2016.02.004      URL     [本文引用: 1]

EBRAHIMI H, MOUSANEZHAD D, NAYEB-HASHEMI H, et al.

3D cellular metamaterials with planar anti-chiral topology

[J]. Materials & Design, 2018, 145:226-231.

[本文引用: 1]

XIA B Z, WANG G B, ZHENG S J.

Robust edge states of planar phononic crystals beyond high-symmetry points of Brillouin zones

[J]. Journal of the Mechanics and Physics of Solids, 2019, 124:471-488.

DOI:10.1016/j.jmps.2018.11.001      URL     [本文引用: 1]

BITTNER S, DIETZ B, MISKI-OGLU M, et al.

Observation of a Dirac point in microwave experiments with a photonic crystal modeling graphene

[J]. Physical Review B Condensed Matter, 2010, 82(1):1558-1564.

[本文引用: 1]

BITTNER S, DIETZ B, MISKI-OGLU M, et al.

Extremal transmission through a microwave photonic crystal and the observation of edge states in a rectangular Dirac billiard

[J]. Physical Review B Condensed Matter, 2012, 85(6):652-660.

[本文引用: 1]

DIEM M, KOSCHNY T, SOUKOULIS C M.

Transmission in the vicinity of the Dirac point in hexa-gonal photonic crystals

[J]. Physica B: Condensed Matter, 2010, 405(14):2990-2995.

DOI:10.1016/j.physb.2010.01.020      URL     [本文引用: 1]

TORRENT D, MAYOU D, SÁNCHEZ-DEHESA J.

Elastic analog of graphene: Dirac cones and edge states for flexural waves in thin plates

[J]. Physical Review B Condensed Matter, 2013, 87(11):269-275.

[本文引用: 1]

CHENG X J, JOUVAUD C, NI X, et al.

Robust reconfigurable electromagnetic pathways within a photonic topological insulator

[J]. Nature Materials, 2016, 15(5):542-548.

DOI:10.1038/nmat4573      URL     [本文引用: 1]

/

〈 〉