声学超材料的复数动态质量密度

Complex Dynamic Mass Density in Acoustic Metamaterials

摘要: 在普通的声学材料中，质量密度通常是一个实数。能量耗散是由材料的形状或体积变化引起的，这与复数模量的虚部有关。这里，我们利用声学超材料也可以实现复数的动态质量密度。在这种情况下，能量损耗是由材料的动量变化引起的。我们分析了这种复数的动态质量密度的物理成因，并提出了计算这种有损耗的声学超材料的复数动态质量密度的理论方法。得到的有效复数质量密度通过有限元模拟进行了验证，包括传输研究和实现相干完美吸收。我们的工作展示了一种实现复数质量密度的方法，这找提高声波的吸收上有重要的应用。

Abstract: The mass density of normal acoustic materials is usually a real number. Energy dissipation is in-duced by the volume or shape change of the materials, which relates to the imaginary parts of complex moduli. Here, we show that by using acoustic metamaterials, complex dynamic mass density can also be realized. In this case, energy dissipation is induced by the change of momentum of the material. We analyze the physical origin of such complex mass density and provide a theoretical approach to calculate the effective complex dynamic mass density for acoustic metamaterials with dissipation. The obtained effective complex mass density is verified by finite element simulations, including both transmission studies and realization of coherent perfect absorption. Our work shows a way to realize complex mass density, which has important applications in enhancing absorption of acoustic waves.

1. 引言

在过去的一个世纪中，声波的吸收由于具有重要应用价值被广泛地研究。例如噪声控制，混响设计等。最近，声学超材料的出现给我们控制声波的吸收提供了一个新的机遇 [1] - [11] 。通过设计适当的局域共振结构，声学超材料可以呈现几乎任意的有效参数，包括负的质量密度 [1] [4] [5] 、负的体积模量 [3] 、负的剪切模量 [6] 、双负 [2] [6] - [9] 和其他不同寻常的密度和模量组成 [10] [11] 。拥有如此不同寻常参数的超材料，不仅可以在超透和高透镜 [12] - [15] ，隐身斗篷和声学幻像 [16] - [27] 上产生新的应用，而且可以极大地提高声波的吸收 [28] - [40] 。特别是，薄膜超材料已经在亚波长尺度实现了对声波的宽频吸收和全吸收 [28] - [30] 。声学相干完美吸收 [32] [33] 的概念和一个声学薄膜实现宽带完美的吸收 [31] 的理论已经被提出。科研工作者们也已经对声学多孔片状晶体 [34] [35] 和其他结构如声学黑洞 [36] - [39] 进行了研究。这些令人兴奋的发现导致了新的声学吸收技术的快速发展。

在普通的声学材料中，质量密度通常是一个实数。损耗与复数的模量有关，表明能量耗散是由材料的形状或体积变化引起的。这是由材料中相对的位移(形状或体积变化)引起的直观的摩擦(损耗的起源)产生的。然而，声学超材料为实现一种非直观的吸收提供了新的可能性，即损耗和复数的动态质量密度有关。这种复数的动态质量密度不仅代表了损耗的另外一种机制，而且还可以支持一些有趣的应用，如声学相干完美吸收 [32] [33] 。特别是，最近的超薄膜组成的纯虚数的质量密度与特定的分散体被发现在自由边界条件下支持声波的宽带完美吸收 [31] 。这样的宽带完美的吸收非常依赖质量密度，但对于弹性模量和速度的要求具有很大的自由度。这些研究结果表明了研究复数的质量密度的重要意义。在这项工作中，我们分析了这种复数的动态质量密度的物理成因，并提出了计算声学超材料的复数的动态质量密度的理论方法。通过对传输特性进行数值研究，实现了相干完美吸收。由此我们证明了该方法的有效性。

2. 计算声学超材料的复数的动态质量密度的理论方法

我们从一个简单的正方晶格弹性超材料出发，它是由铅柱子外包围着一层硅橡胶然后镶嵌在环氧树脂基体中。晶格常数是10 cm；铅柱的半径是3 cm；外围着的硅橡胶的半径是4 cm。环氧树脂的参数分别是，，；硅橡胶的参数分别是，，；铅的参数分别是，， (是密度，是杨氏模量，是泊松比)。对于这样一个声学超材料，我们基于有限元方法用COMSOL多物理场耦合软件计算出了它的能带结构，如图1所示。很显然，在频率大约300 Hz到700 Hz之间存在一个大的完全带隙，这是由晶格的局域偶极共振引起的。这个带隙是由偶极共振产生的负的动态质量密度导致的。

现在，我们在硅橡胶层里加入损耗，而实际情况也经常是这样。为了简单起见，我们把硅橡胶的杨氏模量变成。在这种情况下，我们照样可以计算给定波矢对应的本征频率，但是本征频率通常都是复数值。在图2(a)和(b)中，我们展示了在频率高于700 Hz的纵向偶极模态中，在的虚部为零的情况下本征频率的实部和虚部和的关系图，显然，本征频率是有虚部的。这样的虚部表示能量随着时间消逝。

在这项工作中，我们想要研究的是传输特性而不是瞬态。为了避免波衰逝在时间上，一个方法就是让它消逝在空间上。在图2(c)和(d)中，我们采用了复数的波矢而不是纯实数的波矢。波矢的实部固定为。有趣的是，当波矢的虚部逐渐增大的时候，本征频率的虚部会相应地变小。在图2(d)中，红点处对应的是，对应的本征频率的虚部为零，这表明在时间里的衰逝完全转移到空间里。这正是传输研究中的情况。在图2(c)中，我们可以看到纯实数本征频率仅仅比之前有吸收时稍微小一些。图2(e)图(e)展示的是在保证本征频率为纯实数情况下，的实部和虚部和频率的关系图。

联系到纯实数的本征频率，我们同样得到了相应的衰逝本征场 [10] 。通过分析晶胞每个边界上的应力，应变和位移，我们可以得到相应的动态有效质量密度和有效体积模量。此外，晶格常数为a的二维正方晶格的动态有效质量密度和有效体积模量可以通过以下公式得到：

(1)

和

(2)

其中有效的力和位移可以通过和得到；有效张量和有效应力可以通过，，，，，和得到。

在图3中，我们展示了得到的有效参数和。在图3中，我们可以看到有效体积模量的实部数量级是，而虚部是，虚部比实部小几个数量级。因此，尽管我们已经在硅橡胶层里加入了损耗，有效体积模量几乎仍然是纯实数。在另一方面，在图3(b)中，我们可以看到有效质量的虚部只比实部小一个数量级。特别是在频率700 Hz附近，有效质量的实部是接近于零的，同时有效质量的虚部不为零。因此，在频率700 Hz左右有效质量是接近于纯虚数的。这里的结果表明，可以通过在偶极共振中加入少量的损耗来获得复数的质量密度。通过偶极共振，在质量密度中的损耗(即虚部)可以

Figure 1. The band structure of the elastic metamaterial (inset graph)

图1. 弹性超材料的能带结构(如插图所示)

Figure 2. (a) and (b) shows the real and imaginary parts of the complex eigen-frequency as functions of the real wave vector for the longitudinal dipolar mode, respectively. Loss is added in the silicone rubber coatings by changing the modulus to be. (b) and (d) shows the real and imaginary parts of the complex eigen-frequency as functions of the imaginary part of the wave vector for a given real part of, respectively. Apparently, with increase in, the decreases gradually. At the point of (marked by a red dot in (d) and (e)), the eigen-frequency turn into a purely real number. (e) The real part and imaginary parts of that are required for eigen-frequencies f to be purely real

图2. (a)和(b)分别表示在纵向偶极模态中的虚部为零的情况下本征频率的实部和虚部和的关系图。(c)和(d)分别表示在的实部为时，本征频率的实部和虚部与的虚部的关系图。显然，本征频率的实部和虚部随着的虚部变大逐渐变小。在点 (如图中红点所示)，本征频率变成一个纯实数值。图(e)是在保证本征频率为纯实数情况下，的实部和虚部和频率的关系图

被大大提高。特别是，在特定的频率，质量密度是实部为零，从而实现纯虚数的质量密度。这样一个纯虚数的质量密度可以在提高吸收上有重要的应用 [31] 。

接下来，我们运用数值模拟(通过利用基于有限元方法的多物理场耦合软件COMSOL)来证明的有效参数的有效性。首先，我们研究了传输特性。在图4中，我们研究了透射率和反射率和样品厚度的变化关系图。我们选择了波矢对应着的频率为693.25 Hz这点，这点对应着的有效参数为有效体积模量和有效质量密度。很明显，得到的动态质量密

Figure 3. (a) The obtained for the band of dipolar longitudinal mode; (b) The obtained for the band of dipolar longitudinal mode. The longitudinal coordinates of the left are the real part of the effective parameters and the longitudinal coordinates of the right are the imaginary part of the effective parameters

图3. (a) 通过偶极纵向模态的能带得到的有效体积模量；(b) 通过偶极纵向模态的能带得到的有效质量密度。左边的纵坐标对应着有效参数的实部，右边的纵坐标对应着有效参数的虚部

Figure 4. Reflection and transmission rates as a function of the sample thickness

图4. 反射率和透射率随着样品的厚度变化关系图

Figure 5. Snapshots and amplitude distributions for displacements obtained for both the metamaterial sample and the effective medium, under two beams with opposite phases incident from the left and the right. The background medium is epoxy

图5. 左右两束反相位的波同时入射时，超材料样品和有效介质的图像和振幅分布图。背景是环氧树脂

度为复数，而且虚部比实部大得多。图4的是用有效参数得到的理论上的反射和透射曲线，而图上的点是用超材料样品仿真模拟得到。很明显，理论和仿真模拟的结果吻合的很好，说明了所得到的有效参数的准确性。

最后，我们利用复数动态质量密度的声学超材料来演示相干完美吸收。我们想到从左右两侧同时入射两束反相位的波。我们证明了只要材料具有合适的参数，入射波是会被全部吸收的 [33] 。这就意味着实现了相干完美吸收 [41] [42] 。在图5(a)和5(b)中，我们展示了超材料样品和有效介质的图像和振幅分布图。超材料样品和有效介质的宽和高相同，分别是和。在频率下，波矢，对应着的有效质量密度和体积模量分别是和。背景里的振幅几乎是一致的，这表明背景里几乎没有反射波，即了实现波的相干完美吸收。此外，在超材料样品中包络的场的形状几乎是和有效介质中的一样。

3. 结论

总之，我们研究了声学超材料复数的质量密度的物理起源，即偶极共振和耗散，提出了一种计算声学超材料复数的动态质量密度的有效方法。通过对传输特性和一个相干完美吸收的例子的数值模拟，方法的有效性已被成功验证。复数的动态质量密度意味着损耗是由材料的动量变化引起的，并且在实现声波的宽带完美吸收中有着重要的应用 [28] - [33] 。

基金项目

以上工作是由中国国家重点基础研究项目(编号2014cb360505，2012cb921501号)，国家自然科学基金(编号：11374224)支持，并通过江苏高等教育机构的优势学科建设项目资助(PAPD)。