层状复合材料具有诸多有趣和有用的波动特性，因此在相关领域受到普遍重视。在冲击波物理领域，对层状复合材料中的应力波传播特性研究至少有半个多世纪的历史，人们发现通过改变各层材料的属性、厚度、界面结合强度和加载条件等，可以在层状结构中产生十分丰富的应力波现象。例如：Lundergan等^[1]研究了应力波在10层环氧树脂和10层不锈钢交替组成的结构中的色散现象；Oved等^[2]研究了有机玻璃和钢交替组成的多层结构中的应力波共振现象；Benson等^[3]研究了铜/铝、钛/铝多层结构中的应力波异常衰减现象；Zhuang等^[4]研究了聚碳酸酯分别与铝、不锈钢及玻璃交替组成的多层结构中的应力波，发现应力波速度低于其组成材料；Franco Navarro等^[5]研究了铝/钨层状材料中的应力波传播，发现通过改变加载脉冲的宽度与结构内在特征时间的比值，可以产生3种不同类型的应力波，即振荡型冲击波、一串孤立波、单一孤立波。由于具有对应力波的衰减作用及其他功能（如提升断裂韧性），某些层状复合结构常用作冲击波衰减器和复合装甲^[6-8]。

在振动和声学领域，声子晶体^[9-10]和声学超材料^[10-11]研究于20世纪末21世纪初陆续兴起，层状复合材料属于其中的一维情形。狭义而言，声子晶体和声学超材料有所不同，前者的关键特性之一是由布拉格散射（Bragg scattering，BS）形成频率禁带从而禁止某频率范围的波通过，后者的关键特性是由局域共振（Local resonance，LR）机制产生负密度、负模量、负折射等负参数效应^[12-13]。但在文献中，两种名称常常混用，即声子晶体也可以称为声学超材料，而声学超材料如果具有周期性结构也可以称为声子晶体。声子晶体和声学超材料的波动特性使其在减振降噪中具有重要的应用价值。为方便起见，本研究统一使用声子晶体的名称，但区分LR和BS两种类型。

本研究将报道我们在一维黏弹性声子晶体色散与耗散关系研究中的初步进展。黏弹性是层状复合结构中常用的基础材料（即高分子材料）的基本力学属性，但是从现有文献上看，无论在冲击波物理相关领域还是在振动和声学领域，黏弹性对层状复合结构波动特性的影响都没有得到充分的研究。在声子晶体研究中，大多数都不考虑黏弹性，或者仅用十分简化的黏弹性本构进行描述^[9-10]。直到近年来才将可真实反映高分子材料黏弹性的广义麦克斯韦本构^[14]逐渐应用于声子晶体的相关研究，如利用平面波展开法研究二维黏弹性声子晶体的禁带和耗散特性^[15-17]。

从高强度应力波防护到减振降噪，一维黏弹性声子晶体的色散与耗散关系都值得深入研究。正如本研究将展示的：一方面，黏弹性改变了纯弹性声子晶体的波动衰减特性，甚至改变了其禁带的存在性；另一方面，声子晶体的周期性结构也改变了无结构均质黏弹性材料的波动衰减特性（总的来看是促进了衰减）。

具体而言，本研究将首先建立离散形式的一维黏弹性LR和BS型声子晶体模型，其中离散形式作为连续介质的简化近似^[18-21]，黏弹性由广义麦克斯韦本构模型表示，进而从运动方程出发，推导求解模型的色散与耗散关系。本研究结果有助于进一步推动基于高分子材料的声子晶体研究，也可为冲击波物理相关领域层状复合结构中的应力波研究提供参考。

1.   离散形式的一维黏弹性声子晶体模型

黏弹性材料在任意时刻的行为都与其经历的全部历史有关，其应力-应变关系由遗传积分表示^[14]

式中：$\sigma$为应力，$\varepsilon$为应变，E(t)为模量，s为时间，t为当前时刻。E(t)的具体形式取决于本构模型，对于广义麦克斯韦本构模型^[14]，E(t)可表示为

式中：E_l、$\eta$_l、$\tau$_l分别表示第l个麦克斯韦单元的弹性模量、黏性系数和松弛时间，$\tau$_l = $\eta$_l/E_l；N为麦克斯韦单元的数目。

离散形式的一维黏弹性LR模型和BS模型分别是由图1(a)、图1(b)所示的单元周期性重复连接而成的一维无限长链。在有的文献^[10]中，二者分别被称作Mass in Mass Structure（LR型）和Mass and Mass Structure（BS型），这是从它们各自的质量块与连接体的不同排列形式来作的形象区分。从原理上说，BS结构是利用周期变化的材料特性与波的相互作用，而LR结构则是利用单个散射体（振子）与波的相互作用。设模型中单元的长度为L，质量块的质量为m_a和m_b。各质量块之间链接的模量由式（2）给出。设接触面积为1，则质量块的受力等于应力$\sigma$，由式（1）和式（2）联立确定。

图  1  离散形式的一维黏弹性声子晶体

Figure  1.  One dimensional viscoelastic phononic crystals in discrete form

下载: 全尺寸图片 幻灯片

2.   复波数实频率假定下的色散与耗散关系

2.1   LR模型的色散与耗散关系

考虑到模型的周期性，取第n个单元，对其中质量块m_a和m_b进行受力分析，得到如下两个方程

式中：u为相对位移，$\dot{u}$为位移u对时间s的微分。

代入牛顿第二定律，得到两个质量块的运动方程

式中：上标a、b区分质量块，下标区分单元。

为了得到色散与耗散关系，选取简谐形式的试探解

式中：U为振幅，$\omega$为频率，k为波数，L为单元长度。将式（7）、式（8）代入式（5）和式（6），得到关于频率$\omega$与波数k的方程

消去U_a和U_b，最终得到方程

可以看出，若使式（11）成立，波数k和频率$\omega$中需至少有一个为复数。若同为复数，则方程数目（分离式（11）得到的实部和虚部，2个方程）与未知数数目（两个复数的实部和虚部，4个未知数）不匹配，难以求解讨论。故以下仅将波数k视作复数，频率$\omega$视作实数，或者反过来，将频率$\omega$视作复数，波数k视作实数。这里先采取第一种假定。

令k = $\alpha$ + i$\,\beta$，其中$\alpha$称为真实波数，$\,\beta$是与耗散有关的衰减系数。代入式（11）中，利用欧拉公式，分离方程两边的实部和虚部得到

式中：$A=\displaystyle\frac{L}{2}{m}_{{\rm{a}}}{\omega }^{2}{}$ ，$B=\displaystyle\frac{L}{2}{m}_{{\rm{b}}}{\omega }^{2}$ ，${C}_{1\left(2\right)}= {}\displaystyle\sum\limits _{l=1}^{N}\displaystyle\frac{{E}_{1\left(2\right),l}{\tau }_{1\left(2\right),l}\omega }{1+{\omega }^{2}{{\tau }^{2}_{1\left(2\right),l}}}$ ，${D}_{1\left(2\right)}={}\displaystyle\sum\limits _{l=1}^{N}\displaystyle\frac{{E}_{1\left(2\right),l}{{\tau }^{2}_{1\left(2\right),l}}{\omega }^{2}}{1+{{\omega }^{2}}{{\tau }^{2}_{1\left(2\right),l}}}$，$\xi =\displaystyle\frac{{B}^{2}{C}_{2}}{{\left({D}_{2}-B\right)}^{2}+{{C}^{2}_{2}}}+2{C}_{1}$，$\psi = 2{D}_{1}-A-\displaystyle\frac{B\left({{D}^{2}_{2}}+{{C}^{2}_{2}}-B{D}_{2}\right)}{{\left({D}_{2}-B\right)}^{2}+{{C}^{2}_{2}}}$。

由方程组（12）可以得到

式中：$P={C}_{1}\psi -{D}_{1}\xi$，$Q={C}_{1}\xi +{D}_{1}\psi$，$R={{C}^{2}_{1}}+{{D}^{2}_{1}}$，$x={{\rm{e}}}^{\beta L}$。

由式（13）得到

将式（14）代入式（13）中的第1个方程，得到

求解这个关于${{\rm{sin}}}^{2}(\alpha L)$的一元二次方程，得到

这里由于${{\rm{sin}}}^{2}(\alpha L)$的非负性，舍去了负值解。将式（16）代入式（13）中的第1个方程，得到

至此，就得到了在黏弹性LR模型中真实波数 $\alpha$ 与频率 $\omega$ 的关系，即色散关系式（16），以及衰减系数 $\beta$ 与频率 $\omega$ 的关系，即耗散关系式（17）。

2.2   BS模型的色散与耗散关系

类似地，考虑到模型的周期性，选取第n个单元进行分析，得到单元中质量块m_a和m_b的运动方程

将试探解式（7）和式（8）代入上述运动方程，得到

消去U_a和U_b，最终得到方程

与式（11）类似，若使式（22）成立，可将波数k视作复数，频率 $\omega$ 视作实数，或者反过来，将频率 $\omega$ 视作复数，波数k视作实数。这里同样先采取第一种假定。

令$k=\alpha +{\rm{i}}\beta$，代入式（22）并分离方程两边的实部和虚部，得到方程组

式中：$A=\displaystyle\frac{L}{2}{m}_{{\rm{a}}}{\omega }^{2}$，$B=\displaystyle\frac{L}{2}{m}_{{\rm{b}}}{\omega }^{2}$，${C}_{1\left(2\right)}={}\displaystyle\sum\limits _{l=1}^{N}\displaystyle\frac{{E}_{1\left(2\right),l}{\tau }_{1\left(2\right),l}\omega }{1+{\omega }^{2}{{\tau }^{2}_{1\left(2\right),l}}}$，${D}_{1\left(2\right)}={}\displaystyle\sum\limits _{l=1}^{N}\displaystyle\frac{{E}_{1\left(2\right),l}{{\tau }^{2}_{1\left(2\right),l}}{\omega }^{2}}{1+{\omega }^{2}{{\tau }^{2}_{1\left(2\right),l}}}$，$\xi =2({C}_{1}{D}_{2}+{C}_{2}{D}_{1})- (A+B) \text{·} ({C}_{1}+{C}_{2})$，$\psi =2\left({D}_{1}{D}_{2}-{C}_{1}{C}_{2}\right)-\left(A+B\right)\left({D}_{1}+{D}_{2}\right)+AB$。求解得到与式（13）形式相同的结果，区别在于这里$P= \left({C}_{1}{D}_{2}+{C}_{2}{D}_{1}\right)\psi -\left({D}_{1}{D}_{2}-{C}_{1}{C}_{2}\right)\xi$，$Q=\left({D}_{1}{D}_{2}-{C}_{1}{C}_{2}\right)\psi +\left({C}_{1}{D}_{2}+{C}_{2}{D}_{1}\right)\xi$，$R={\left({D}_{1}{D}_{2}-{C}_{1}{C}_{2}\right)}^{2}+ \left({C}_{1}{D}_{2}+ {C}_{2}{D}_{1}\right)^{2}$，$x={{\rm{e}}}^{\beta L}$。可以发现，继续求解后，最终得到形式上与式（16）和式（17）完全相同的色散与耗散关系式。

2.3   关于不存在禁带的判断

色散关系式（16）的右端显然非负，令其大于1，即

可推导得到${R}^{2}{Q}^{2}<0$，这在$\omega$为实数时显然不成立，因此证得0≤ sin²($\alpha$L) ≤ 1。

于是可以发现，无论频率$\omega$为何值，均满足sin²($\alpha$L)有解的条件，换言之，总有真实波数$\alpha$与之对应。于是可以确定LR和BS两种模型的色散关系中均不存在禁带。

值得注意的是，对于纯弹性LR与BS结构而言，色散关系中是存在禁带的。为了探究上述色散关系在黏弹性退化至纯弹性（松弛时间 $\tau$ 趋于无穷大）时是否存在禁带，绘制了如图2所示的不同松弛时间下的色散关系图，其中：N = 1，m_a = 1，m_b = 2，L = 2，E₁ = E₂ = 1。

图  2  不同 $\tau$ 下LR和BS模型的色散关系

Figure  2.  Dispersion relation of LR and BS model with different $\tau$

下载: 全尺寸图片 幻灯片

从图2可以看出，对两种结构而言，色散关系中确实不存在禁带，即对每个频率取值都可以找到与之对应的非零真实波数$\alpha$，色散关系难以退化至纯弹性情形下的“禁带”模型。将式（11）和式（22）与弹性情形下LR、BS结构的方程^[22-23]对比，可以发现，波数k仅出现在自然指数项中，且形式与弹性情形下的方程一致，而弹性情形下的弹性模量E与黏弹性情形下的“等效模量”E_eff相对应。

当松弛时间$\tau$趋于无穷大时，有${E}_{{\rm{e}}{\rm{f}}{\rm{f}}}=E$，意味着式（11）和式（22）在 $\tau$ 趋于无穷大时可以退化至纯弹性情形。但是，若将k视作复数，原先的方程被分离成实部和虚部两个方程进行后续的求解，就破坏了方程的形式，而造成复数波数假定下的色散关系无法退化至纯弹性情形。

2.4   周期性调制的衰减作用

由于在复数波数实数频率假定下，上述BS和LR模型的色散关系中不存在禁带，所以波的衰减将完全由黏性耗散和周期性调制引起。如果将单原子链视作均匀材料的离散化，那么将质量块交替排列的BS模型与之对比，可以看出由质量差异排列导致的周期性调制所带来的影响。

由耗散关系式（17）可以得到

将式（16）代入，绘制出耗散曲线（$\beta$-$\omega$关系），如图3所示。各项参数为：E₁ = E₂ = 1，$\tau$₁ = $\tau$₂ = 1，L = 2。可以看出，质量比越大，耗散系数（其绝对值）越大，耗散越明显，因此可以认为，虽然此时周期性结构中没有禁带，但周期性调制增强了对波的衰减作用。

图  3  N = 1和N = 3时BS模型的耗散关系

Figure  3.  Dissipation relation of BS model when N = 1 and N = 3

下载: 全尺寸图片 幻灯片

类似地，可以得到LR结构的耗散曲线，如图4所示。图4中也反映了与BS结构的对比情况。可以看出，LR结构产生的耗散作用只在某一频率区间内大于BS结构。比较不同质量比下的耗散曲线，可以看出，质量比越大，耗散作用越明显，这一点与BS模型一致。

图  4  N = 1和N = 3时LR模型的耗散关系

Figure  4.  Dissipation relation of LR model when N = 1 and N = 3

下载: 全尺寸图片 幻灯片

综上，虽然不存在禁带，但黏性产生的耗散作用仍会导致波的衰减，且对于两种模型而言，周期性质量差异的调制会增强这种耗散作用。

3.   复频率实波数假定下的色散与耗散关系

下面采取复频率实波数假定。选取如下试探解形式

式中：波数k为实数；令$\lambda$ = $\alpha$ + i$\beta$，其中 $\,\beta$ 为真实频率，$\alpha$ 是与耗散有关的衰减系数。

3.1   LR模型的色散与耗散关系

将式（27）代入运动方程式（5）和式（6）中，得到

消去U_a和U_b，得到

求得 $\lambda$ 后分离实部和虚部，虚部 $\,\beta$ 关于k的表达式便是色散关系式，实部 $\alpha$ 关于k的表达式是耗散关系式。这是一个关于 $\lambda$ 的高次方程。阿贝尔曾证明，高于四次的一般方程没有求根公式，因此只能得到N = 1时 $\lambda$ 的解析解，而对于N > 1的情形，由于 $\lambda$ 的次幂大于4，故无法得到解析解。但原则上可以通过数值计算方式得到给定k情形下 $\lambda$ 的解，从而绘制出色散图和耗散图。

3.2   BS模型的色散与耗散关系

将试探解式（27）代入运动方程式（18）和式（19）中，得到

消去U_a和U_b，得到

求得 $\lambda$ 后分离实部和虚部，得到色散与耗散关系式。同样地，这个关于 $\lambda$ 的方程仅在N = 1时存在解析解，当N > 1时由于 $\lambda$ 的次幂大于4而无法得到解析解。但原则上依然可以通过数值计算方式求得某一给定k情形下的 $\lambda$，即可画出色散图和耗散图。

3.3   N = 1和N = 3时色散和耗散关系的进一步讨论

为了进一步探究色散与耗散关系，本研究具体探讨N = 1和N = 3两种情形。前者对应方程有解析解，即可以得到色散和耗散关系的解析表达式；后者则需要通过数值求解得到散点图。以下情形均舍去了频率虚部为零的解，即不考虑真实频率为零的情况。

在N = 1的情形下，为了简化式（30）和式（33）且不失一般性，可令m_a = 1，m_b = 2，E₁ = 1，E₂ = 2，$\tau$₁ = $\tau$₂ = 1，L = 2。

对于LR模型，代入式（30），整理得到

求解这个四次方程，得到4个解析解

分离这4个解的虚部，在$k\in \left[-{\text{π}}/L,{\text{π}}/L\right]$区间画出色散曲线，如图5（a）所示。

图  5  N = 1时LR和BS模型的色散关系

Figure  5.  Dispersion relation of LR and BS model when N = 1

下载: 全尺寸图片 幻灯片

对于BS模型，式（33）可以化简为

求解这个四次方程，得到4个解析解

分离这4个解的虚部，在$k\in \left[-{\text{π}}/L,{\text{π}}/L\right]$区间画出色散曲线，如图5（b）所示。

可以看出，这种耗散由时间的系数体现，即振幅随时间衰减的情形下，两种模型是存在禁带的，在禁带范围内，对真实频率$\,\beta$而言没有与之对应的波数k。

为了探究当 $\lambda$ 的方程为更高次时是否依然存在禁带，研究N = 3的情形。如前所述，此时的方程由于幂次大于4而没有解析解，只能通过数值计算方式得到某一给定波数k对应方程的数值解，再分离出虚部$\,\beta$。当给定的k足够多时，可以得到足够多的数据组(k, $\,\beta$)，从而得到色散关系图。

对于LR模型，为了简化方程，不妨假定连接两质量块m_a和m_b的两个广义麦克斯韦模型相同，令m_a = 1，m_b = 2，E_1,1 = E_2,1 = 1，E_1,2 = E_2,2 = 1，E_1,3 = E_2,3 = 1，$\tau$_1,1 = $\tau$_2,1→∞，$\tau$_1,2 = $\tau$_2,2 = 1，$\tau$_1,3 = $\tau$_2,3 = 2，L = 2。此时的模量为

将上述参数代入式（30）并化简，得到

这是关于$\lambda$的八次方程。首先在$\left[-{\text{π}}/L,{\text{π}}/L\right]$区间内取多个实数k，数值求解得到各自对应的方程的解，分离出虚部$\,\beta$，得到的色散关系如图6（a）所示。

图  6  N = 3时LR和BS模型的色散关系

Figure  6.  Dispersion relation of LR and BS model when N = 3

下载: 全尺寸图片 幻灯片

对BS模型作相同的假定，将参数代入式（33）并化简，得到

得到的色散关系如图6（b）所示。

从图6可以看到，当N = 3时禁带依然存在。尽管无法穷举或严格证明，我们似乎可以谨慎地推断N取任意值时都存在禁带（实际使用广义麦克斯韦模型时，N一般不超过20）。

值得注意的是，这里的色散关系与纯弹性情形类似^[22-23]，这是因为将i$\omega$视作一个复数$\lambda$，式（11）和式（22）成为关于$\lambda$的高次方程，即式（30）和式（33）。与复数波数假定不同，这里不需要分离实虚部求解，原方程的形式得到了保留，且这个高次方程在松弛时间趋于无穷大时会退化为关于$\lambda$的四次方程，其形式与弹性情形下关于$\omega$的四次方程一致。因此也可以预见，复数频率假定下的色散关系可以退化至弹性情形。

对式（34）、式（36）、式（39）、式（40）分离出解的实部，便得到耗散关系，如图7和图8所示。可以看出，N对耗散曲线形状的影响较大。当N相同时，LR和BS两种结构的耗散曲线相似，说明耗散作用受到结构的影响不大，应主要由黏性引起。

图  7  N = 1时LR和BS模型的耗散关系

Figure  7.  Dissipation relation of LR and BS models when N = 1

下载: 全尺寸图片 幻灯片

图  8  N = 3时LR和BS模型的耗散关系

Figure  8.  Dissipation relation of LR and BS models when N = 3

下载: 全尺寸图片 幻灯片

4.   结　论

研究了一维广义麦克斯韦黏弹性情形下LR和BS型声子晶体的色散与耗散关系。在求解运动方程时，采用了波数为复数、频率为实数，以及频率为复数、波数为实数的两种不同假定，两种假定给出了禁带存在与否的相反结果。

对于复波数假定，推导所得的LR和BS型声子晶体的色散关系在形式上一致，都不存在禁带，两者的耗散关系在形式上也是一致的。对于复频率假定，LR和BS型声子晶体均得到关于时间系数$\lambda$的一个高次方程，求解该方程并分离虚部和实部，得到色散和耗散关系，此时色散关系存在禁带，且二者的耗散关系类似。

鉴于复波数假定时黏性造成波幅随空间衰减，该假定适合于时间谐波传播问题^[24-25]；复频率假定时，黏性造成波幅随时间衰减，该假定适合于自由波传播问题^[24-25]。因此可以确定：在广义麦克斯韦黏弹性情形下，对于LR和BS两种声子晶体中的时间谐波传播问题，禁带不能发挥作用，波的衰减完全依赖于黏性耗散和周期性调制，且周期性调制会增强耗散作用；相反，对于这两种声子晶体中的自由波传播问题，禁带能够发挥作用，但在禁带之外，波动的衰减仍需借助于黏性耗散和周期性调制。

本研究不失普遍性的定性认识对于减振降噪和冲击波防护研究具有一定的参考价值。在实际应用中，还需进一步研究连续介质形式的一维声子晶体（或层状复合结构），使用真实材料参数，针对波动或振动能谱中能量集中的频率范围，优化设计结构，以期最大程度地发挥禁带、黏性耗散和周期性调制对波动或振动衰减的作用。