用各向材料实现宽带声波功能器件设计

1、.用各向同性材料实现宽带声波功能器件设计利用坐标变换的方法并结合保角映射技术,本文介绍了一种利用常规的各向同性材料来设计声波器件的方法。基于此理论,设计了二维声波隐身斗篷,并进展有限元模拟,证明了该器件的有效性,另外由于设计中没有利用材料共振的性质,所以器件是宽带有效的。该方法将有助于拓宽声波功能器件的设计,并为实验验证声波器件提供了可能。:保角坐标变换方法,声波器件设计,隐身斗篷0 引言声波传播和散射问题的深化研究是一个有趣的科学命题,蕴含着重要的应用背景。但是鉴于数学描绘的复杂性,人们对于声波传播的认知还主要集中于均质和各向同性介质中。对于非均质、各向异性介质中波传播的特性,人们认知还比较

2、少。不过真实世界中的介质往往是非均匀且各向异性的。在这样的条件下波的传播问题更加复杂,但也更加精彩隐身斗篷,因此值得深化的讨论。假如说给定材料和几何边界的情况下去求解波场分布情况是一个正问题,那么根据给定外界波场反推材料属性和相关的几何边界情况就是一个波的反问题。正问题的研究已经非常困难,反问题的研究那么更具有挑战性,但这样的研究却是很有意义的,因为它解决的其实是一个波的控制问题。本文中声波器件实际上就是一些经过精细设计的特殊波导构造;在构造外形和材料分布的控制下,以到达改变波的传播特性的目的。与传统的波导构造相比,它具有许多新的特点,对于波场的控制有了极大的灵敏性,应用前景广泛。当前许多研究

3、领域,都对电磁波和声波的控制问题给予了很高的重视,也获得了很大的进展隐身斗篷,追根溯源,除去新材料的开展,可以不夸大的说,一个很重要的原因在于所谓的坐标变换方法的运用。该方法由Pendry和Leonhardt 1,2最先提出,本文的设计也是基于此方法,再进展改进后得到的。针对电磁波的控制问题,Pendry等人提出的坐标变换方法的核心思想是利用Maxwell方程组在坐标变换下的形式不变性进展材料设计,近年来得到了许多有趣的电磁波控制器件,如电磁波平面透镜、弯波器、以及隐身斗篷等,其中一些器件逐渐被随后的实验所验证3,4。实际上该方法并不仅仅限于电磁波,把它用于声波控制的研究也得到了越来越多的重视

4、5隐身斗篷,这一研究分支被称为变换声学。但是根据目前的大多数变换声学方法所设计出声波器件都具有很强的各向异性的性质6,例如要求动态质量密度为二阶张量形式等。为了到达这种要求往往需要在部分共振意义上进展动力学等效处理,这样不但使设计很难在工程中实现,而且也会受到工作带宽的限制。针对这样的问题,本文将介绍一种利用保角变换的方法来加以改进的思路,使最终的设计仅用通常的各向同性工程材料就可以实现,而且整个设计过程中没有部分共振的要求,从而可以实现宽带效果。1 保角坐标变换方法下面首先给出此问题的公式推导7。不同于已有理论,在公式推导中,没有假设材料参数为二阶张量形式,即认为所有材料参数仍保持为标量。标

5、量形式的声波的控制方程为:1式中是声波压力, 为声介质密度, 为声介质体积模量。对于任意给定的坐标变换隐身斗篷,将原空间的坐标映射到像空间中,同时要求,可以证明方程1的左边将变成:,2式中。在保角变换时,由于雅可比矩阵可以表示为,其中是一个正的标量,而 是一个正交矩阵。再将代入方程2,得到。因此方程1可以表示为:,3式中 , 与分别是和的标量缩放因子。进一步考虑阻抗匹配的条件:4并要求公式1与3具有一样的表达形式,就可以得到最终器件设计所需的材料参数:,5这样得到的器件完全可以通过使用通常的各向同性材料实现隐身斗篷,从而使得到的声学器件更加便于工程实现。2 声波隐身斗篷实现利用上节的理论公式,

6、本节希望进展声学器件的设计。其中最具代表性的应该是隐身斗篷。它可以看作是目前电磁波隐身斗篷在声波领域的一个拓展。本文所给出的例子在文献8中曾被讨论过,其根本思 如图1所示:黄色梯形区域代表隐身区域,当入射波从左侧入射的时候,由于外部隐身器件显示为灰色的区域的作用,导致声波发生弯曲。根据精细设计的隐身斗篷会改变声波的传播途径,使其光滑地绕过隐身区域,之后在其前方又恢复成原来的途径,同时声波不会在隐身区域的前方发生反射。这样无论在隐身区域前方或者是前方进展声波探测,其所得的波场与隐身区域不存在时探测的情况将是一样的,从而最终实现了对区域隐身的目的。遗憾的是隐身斗篷,由于Pendry所用方法的限制,

7、如前面所说,其所需材料要具有较强的各向异性的性质,这给工程实现带来了很多困难。图1 隐身斗篷示意图Fig.1 Scheme of the tunnel cloak为了说明本文方法的有效性,本文重新设计了该种隐身斗篷。首先需要构造恰当的保角变换,这里采用的是Schwarz-Christoffel S-C变换。S-C变换是复变分析中常见的用于产生保角变换的方法,根本思路可以参见文献9。这里我们将图2a映射到图2b,即将半平面映射为在线段和射线, 上的半平面。图2 利用S-C变换得到变换关系Fig.2Construction of the transformation这样可以直接得到映射方程:,式中

8、6从图2b不难发现,只要在此根底上将其沿轴作镜像复制就可以得到最终的隐身斗篷。在斗篷的外部形式之后,通过公式5和6,不难进一步得到其中材料分布的情况隐身斗篷,如图3所示。其中图3a为密度缩放因子,图3b为体积模量因子。图3 隐身斗篷内的材料分布情况Fig.3 Distributionprofile of the material parameters a mass density and b bulk modulus为了进一步说明所提出的方法的有效性,本文还利用大型商用有限元软件Abaqus对设计得到的隐身斗篷进展了数值模拟。图4给出了有隐身斗篷与没有隐身斗篷的情况下的声压场分布情况,入射波为

9、平面压力波,入射波频率。图中虚线给出了斗篷的外边界。从比照的情况不难看出,在4b有隐身斗篷的时候,声波传播的途径按照预期发生了弯曲,之后又恢复了其原有的途径;当4a没有隐身斗篷时,入射声波与内部边界发生了明显的散射,外场呈现无规律的分布。在这种情况下,可以说位于梯形隐身区域内部的任何物体对于声压场是没有影响的。从另一个角度讲隐身斗篷,位于隐身区域的物体对于外界观察者来说似乎并不存在一样,到达了预期的隐身效果。这说明本文中所给出的方法是实在有效的。图4 波场分布情况比照Fig.4The pressure fields awithout cloak and b with cloak at由于本设计

10、方法没有利用材料的部分共振性质,所以自然回避了共振隐身斗篷带宽过窄的缺点。当入射波频率为时,从图5不难发现该斗篷的隐身效果仍然比较好。图5波场分布情况比照Fig.5 The pressure fields a without cloak andb with cloak at3 结论本文利用保角变换的方法,给出了设计各向同性声波器件的思路,并首次设计了楔形的声学隐身斗篷。通过有限元模拟,说明该方法设计的隐身斗篷是有效的。另外由于本方法回避了部分共振材料的使用,使得设计的器件具有宽带的性质。该方法虽然目前限于研究弹性波中较简单的声波问题,但其思路对于讨论固体波器件是有借鉴价值的10。固体波斗篷和其

11、它器件设计中保角变换的应用是我们下一步工作研究的重点,作者相信利用保角变换降低器件设计中对材料各向异性的要求是很有前景的。1Pendry J B, Schurig D, Smith D R. Controlling electromagnetic fields J. Science, 2019, 312: 17801782.2Leonhardt U. Optical conformal mappingJ. Science, 2019, 312: 17771780.3Schurig D, Mock J J, Justice B J, Cummer S A, Pendry J B, Starr A

12、 F, Smith D R.Metamaterial electromagnetic cloak at microwave frequenciesJ. Science, 2019,314: 977980.4Liu R, Ji C,Mock J J, Chin J Y, Cui T J, Smith D R. Broadband ground-plane cloakJ.Science, 2020, 323: 366369.5Chen H Y, ChanC T. Acoustic cloaking and transformation acousticsJ. Journal of Physics

13、D:Applied Physics, 2019, 43, 113001.6Cummer S A, Schurig D. One path to acoustic cloakingJ. New Journal of Physics,2019, 9, 45.  7Ren C Y, Xiang Z H, Cen Z Z. Design of acoustic devices with isotropic materialvia conformal transformation J. Applied Physics Letters, 2019, 97, 0441014.8Li J, Pendry JB. An acoustic metafluid: realizing a broadband acoustic cloakJ. New Journal of Physics, 201