幂指数棱台声子晶体对薄板振动弯曲波的调控特性研究

2024-01-10张振华

振动工程学报 2023年6期

金星，张振华

（海军工程大学舰船与海洋学院，湖北武汉 430033）

引言

板壳类结构在船舶、航空、车辆等领域有着广泛的应用，其减振降噪问题一直是工程界的重要关注点。声学超材料概念的提出为薄板结构的减振降噪提供了新的技术途径。其中声学黑洞（ABH）［1］的结构厚度设计成幂指数函数（h（x）=ε|x|m（m≥2））的形式，可使沿该方向传播的振动波能流集中在结构的尖端位置，再通过附加阻尼材料吸收振动能量以达到减振效果，其在梁或板的减振［2］和能量调控［3］等方面具有广阔的应用前景。此外，周期性结构产生的带隙也打开了结构减振降噪的另一条思路。

近年来国内外学者对带有声学黑洞的周期性结构进行了大量研究。Zhu 等［4］采用平面波展开法和有限元法对嵌入声学黑洞薄板结构进行了研究，发现其可对弯曲波产生双折射效果。Zhao 等［5］提出了一种改进的声学黑洞板结构，并通过数值模拟与仿真验证其可对弯曲波产生准直与聚焦效应。Tang 等［6］利用小波分解能量法对嵌有多个声学黑洞的欧拉-伯努利梁进行了研究，发现其因声学黑洞效应产生了低频段的弯曲波带隙。Gao 等［7］提出了一种中心嵌有组合声学黑洞的梁结构，发现其在1200 Hz 频率下存在两个带隙，并指出第一个带隙是由纵向和横向弯曲振动的耦合作用引起的。Deng 等［8］利用高斯展开法计算了环形声学黑洞的频散曲线，通过分析等频线，发现其可在特定频率对弯曲波产生自准直效应。他还提出一种嵌有声学黑洞的圆柱壳结构，并在此基础上通过添加加强筋来增强壳体刚度，研究表明这种结构可有效阻隔弯曲波在壳体的轴向传播［9］。同时还利用高斯展开法对嵌有十字型声学黑洞凹槽的板结构进行了研究，利用结构的带隙特性抑制了板的弯曲振动［10］。Tang 等［11］提出了一种嵌有双叶型声学黑洞的梁结构，并结合局域共振效应与布拉格散射效应产生了宽频带隙。Zhou 等［12］通过数值模拟研究了一种复合型声学黑洞梁的动态和静态特性，结果表明复合型声学黑洞梁在结构强度与抑制振动方面都优于传统的声学黑洞梁结构。O’Boy等［13］利用瑞利-里茨变分能量法对声学黑洞板的固有频率进行了研究，并发现在凹槽中将适量阻尼比附加于整个板具有更好的减振效果。Ji 等［14］提出了一种可贴附于结构上的圆形ABH 吸振器结构，并通过有限元法与试验验证其具有较好的减振效果。

声子晶体的缺陷态也是目前的一个研究热点。声子晶体周期性的破坏称为声子晶体的缺陷。Sigalas［15］发现在以铅圆柱为散射体，环氧树脂为基体的二维声子晶体结构中，改变散射体半径而形成的缺陷可对弹性波产生局域作用。赵寰宇等［16］研究了以水为基体、钢柱为散射体的声子晶体的带隙特性，与实验测得的缺陷局域模态频率对比具有较好的一致性。江卫锋等［17］发现将声子晶体的点缺陷引入声子晶体的周期性结构中可以扩大声子晶体的带隙范围。徐驰等［18］验证了声子晶体的线缺陷具有聚声效应。

虽然目前对声学黑洞已开始了相关研究，但利用周期性声学黑洞的带隙特性实现结构减振的研究比较匮乏，且提出的构型大多只存在方向带隙［19］，尚少有关于利用幂指数结构周期性排列形成完全带隙从而抑制板弯曲振动的研究。传统的声学黑洞构型的结构厚度通常以幂指数函数的形式减小，会削弱结构强度。本文借鉴声学黑洞构型，提出了一种幂指数棱台声子晶体，可通过胶接的方式将轻质材料制成的幂指数棱台结构贴附于薄板上，利用声子晶体的聚焦效应与带隙特性抑制板的弯曲振动。研究发现该声子晶体具有三个宽频弯曲波完全带隙。结合其模态振型和频散曲线，对带隙产生的机理进行了分析，并研究了带隙的影响因素；通过数值仿真和试验对结构的带隙特性进行了验证；最后利用声子晶体的线缺陷与带隙特性实现了弯曲波波导。

1 研究模型

如图1 所示，声子晶体原胞由幂指数棱台与薄板组成。所述幂指数棱台是指：上端面和下端面为正方形，棱面呈幂指数曲线向内侧弯曲，上端面中心设置有柱形空洞的棱台，如图2 所示。幂指数棱台底座的正方体边长为a，边缘厚度为hA，总高度为HA，凹陷宽度为rA。幂指数棱台在x与y方向的厚度变化为棱台中间的圆柱形空洞半径为rA1，空洞的深度为HA1。薄板的长宽与棱台下端面一致，厚度为hB。

图1 声子晶体周期性结构Fig.1 The periodic structure of phononic crystals

图2 声子晶体晶胞Fig.2 The crystal cell of phononic crystal

2 幂指数棱台声子晶体带隙分析

2.1 声子晶体中弹性波基本方程

根据弹性力学理论，薄板内弯曲波的控制方程可写为：

式中t为时间；r为位置矢量；w(r，t)为板的横向位移；ρ(r)为密度；D(r)=E(r)h3(r)/{12[1-ν2(r)]}为板的弯曲刚度；E(r)为板的杨氏模量；h(r)为板的厚度；ν(r)为泊松比。

2.2 频散曲线计算

根据Bloch 定理，声子晶体板的位移场可表示为：

式中k=(kx，ky)为第一布里渊区的Bloch 波矢；R为平移矢量；ω为圆频率；wk(r)为与材料参数具有相同周期的周期函数。

第一布里渊区如图3 阴影所示，根据原胞的平移周期性，通过改变第一布里渊区的波矢k，利用有限元软件求解特征值，可得到频散曲线。声子晶体的几何参数与材料参数如表1 所示，其中薄板材料为钢，幂指数棱台材料为聚碳酸酯（PC）。由表1 参数计算出幂指数棱台声子晶体的频散曲线如图4所示。

表1 声子晶体的几何参数与材料参数Tab.1 Geometrical parameters and material parameters of phononic crystals

图3 声子晶体第一布里渊区Fig.3 The first Brillouin zone of phononic crystal

图4 声子晶体频散曲线Fig.4 Frequency dispersion curves of phononic crystals

在有限厚度的板结构中，一般具有弯曲波、纵波、水平剪切波三种基本波模式。三者可分别表示为A 模式、S 模式和SH 模式，且与频率阶数（n≥0）相关，可以分别写为An，Sn和SHn（n=0，1，2，…）［20］。为分析频散曲线对应的波模式，图5 列出了频散曲线中点A～E 本征模态在x，y和z方向上的位移分量。其中，点A，B，C，D 处本征模态在z方向上的振动幅值远大于其他两个方向，而点E，F 处本征模态则分别主要在x，y方向产生较大的位移，这表明点A～D 所在频散曲线（图4 黑点）对应弯曲波模式，点E，F 所在曲线（图4 红点）则对应其他波模式［21］。因薄板结构主要在垂向产生较大的振动，故本文重点关注频散曲线中的弯曲波模式，若隐去其他类型的波模式，可以发现幂指数棱台声子晶体具有三个弯曲波完全带隙，其频段如表2 所示。

表2 声子晶体完全带隙范围Tab.2 Range of complete band gap of phononic crystals

图5 位移场分布Fig.5 Distribution of the displacement field

为更好地揭示带隙机理，计算了自由边界条件下声子晶体原胞的固有频率，其第一阶固有频率为3475 Hz，对应第一带隙频段。由图6 可知，原胞模态振型显示波的能量主要聚集于棱台的斜边处，这是由于棱台厚度的变化导致弯曲波的波长减小，振幅增大，降低了弯曲波的群速度与相速度，导致弯曲波能量被限制于幂指数斜面中，即能量聚焦效应。这表明第一带隙是因聚焦效应引起的局域共振产生的［21］。声子晶体原胞的第五阶固有频率（6806 Hz）对应第二带隙；第八阶固有频率（8749 Hz）对应第三带隙。图6 中原胞固有频率的模态振型因聚焦效应出现了能量在幂指数斜面聚集的现象，这表明幂指数棱台产生的三个完全带隙都是因聚焦效应引起的局域共振产生的。

图6 固有频率模态图Fig.6 Modal diagram of natural frequency

2.3 带隙影响因素分析

为揭示声子晶体的结构参数对带隙的影响规律，对声子晶体的参数进行研究。

2.3.1 幂指数棱台总高度HA变化对带隙的影响

幂指数棱台声子晶体的总高度对带隙的影响规律如图7 所示。其中声子晶体晶胞只改变总高度HA，其他参数与表1 一致。

图7 总高度变化对带隙的影响Fig.7 The influence of total height variation on band gaps

通过图7 可知，随着棱台总高度HA的升高，第一带隙的起始频率逐渐降低，而其他两个带隙的起始频率则升高。三个带隙的终止频率都是随着HA的增高而提高。对于带隙范围而言，第一带隙受HA的影响最大，三个带隙的带宽随着HA的升高而变宽，当高度为13 mm 时，第一带隙开始出现。

2.3.2 幂函数幂次m变化对带隙的影响

幂函数的幂次对于声子晶体带隙也有较大的影响，其影响规律如图8 所示。其中声子晶体原胞只改变幂次m，其他参数与表1 一致。

图8 幂次变化对带隙的影响Fig.8 The influence of power variation on band gaps

2.3.3 幂指数棱台声子晶体的边缘厚度hA变化对带隙的影响

幂指数棱台声子晶体的带隙对其边缘厚度十分敏感，如图9 所示。其中声子晶体原胞只改变边缘厚度hA，其他参数与表1 一致。

图9 边缘厚度对带隙的影响Fig.9 The influence of the thickness of the border on band gaps

随着边缘厚度hA的增加，所有带隙的频段逐渐变窄，且带隙的起始频率开始升高，当hA=0.9 mm时，第一带隙消失。这是由于边缘厚度的升高，聚焦效应开始减弱，从而导致局域共振效应减弱［21］。由此可知，可以通过降低棱台边界厚度的方式获得起始频率低且带宽更宽的带隙。

3 幂指数棱台声子晶体板对弯曲波的调控

3.1 幂指数棱台声子晶体板对弯曲波的阻隔

频散曲线可以预报带隙的频段，但还需通过分析传递损耗曲线去验证声子晶体的带隙特性，从而充分证明声子晶体板的减振特性。如图10 所示，声子晶体板长为0.3 m，宽为0.2 m，厚度为0.5 mm，其上贴附有4×4 个声子晶体原胞，声子晶体板四周为固支约束，且声子晶体的几何及材料参数与表1 一致。在点P1 处加载一个垂直于板的激励载荷，并拾取点P1 与P2 处的加速度用于计算传递损耗曲线。研究中通常将传递损耗曲线定义为［23］：

图10 振动传递模型Fig.10 The vibration transfer model

式中X1为点P2 处的加速度；X0为点P1 处的加速度。

图11 给出了振动传递模型在激励频率为1000～10000 Hz 作用下的传递损耗曲线。相较于平板结构，声子晶体板在3200～3800 Hz，6300～7100 Hz，8800～10000 Hz 三个频段对弯曲波具有较好的衰减效果，声子晶体板的振动损耗传递最高可达-45 dB，且能量衰减段与带隙频段（图11 蓝色虚线）基本一致。

图11 传递损耗曲线Fig.11 The curves of transfer loss

为进一步验证幂指数棱台声子晶体带隙阻隔弯曲波的有效性，根据表1 参数制备了声子晶体板试验模型，如图12 所示。声子晶体板由3D 打印机（Raise3D Pro2）加工出的聚碳酸酯棱台与尺寸为300 mm×300 mm×0.5 mm 的薄钢板组成，利用高强度AB 胶将2×5 个棱台粘结于薄板中部。

采用螺栓将声子晶体板固定于刚性支架上，用于模拟固支约束条件。将激振器（DH40050）与升降台放置于模型下方，调整升降台高度，使激振器顶杆置于声子晶体板点P1 下方，以激励板产生振动。利用信号发生器（FY6900）与功率放大器（DH5872）对激振器的输出激振力进行设定。此外，在声子晶体板点P1，P2 处贴附微型加速度传感器（1A803E）用于测量两点的法向振动加速度，试验测试系统如图13 所示。

在课程进行到三分之一左右时开始发放文献阅读任务，此时学生既具有一定的知识基础，又有充分的时间完成任务。

图13 试验测试系统Fig.13 Experimental measuring system

试验中激振器的输入信号为1000～4000 Hz 的正弦波信号。利用采集仪自带软件对点P1，P2 的时域信号进行快速傅里叶变换，得到的频域信号用于计算传递损耗曲线，试验结果如图14 所示。

图14 传递损耗曲线Fig.14 The curves of transfer loss

由图14 可知，试验模型在3210～3700 Hz 频段具有明显的能量衰减段，有效抑制了板的弯曲振动响应，且能量衰减频段与带隙频段（蓝色虚线框）基本一致，验证了幂指数棱台声子晶体弯曲波带隙的减振特性。此外，试验样件在非带隙频段同样具有衰减振动的效果，这是由于声子晶体的幂指数剖面产生的能量聚集效应抑制了板的弯曲振动。试验结果表明，在带隙频段内的声子晶体板减振效果要优于在非带隙范围内的减振效果。

3.2 幂指数棱台声子晶体板对弯曲波的波导

缺陷是指对声子晶体理想周期性结构的破坏。可通过改变周期性结构中声子晶体的几何参数或是材料参数等方式生成缺陷。

对含Z 型缺陷的声子晶体板进行传递损耗计算。在长度为0.4 m，宽度为0.35 m，厚度为0.5 mm的声子晶体板结构中设置一个Z 型线缺陷，如图15所示。在点P1 处施加一个频率为1000～10000 Hz，方向垂直于板，大小为1 N 的简谐载荷，并拾取点P1与P2处的加速度，计算传递损耗曲线，结果如图16所示。相较于无缺陷的声子晶体板结构，具有线缺陷的声子晶体板的振动衰减效果有所下降。由此可见，线缺陷的引入会弱化结构的带隙特性。

图15 Z 型线缺陷Fig.15 Z-shaped line defects

图16 传递损耗曲线Fig.16 The curves of transfer loss

图17 展示了具有Z 型线缺陷的声子晶体板在载荷频率为2.6 与3.5 kHz 作用时的位移场分布。频率在带隙内的弯曲波可以沿着线缺陷进行传播，频率在带隙外的弯曲波则没有形成明显的波导效果。

图17 Z 型线缺陷声子晶体板的位移场分布Fig.17 The distribution of displacement field of phononic crystal plate of Z-shaped line defects

在声子晶体薄板结构的中心处设置一个直线型线缺陷，并在点P3 处设置一个频率为3400 Hz，垂直于板的简谐载荷，如图18 所示。当频率在带隙范围内的载荷作用于具有直线型线缺陷的声子晶体板结构时，弯曲波可平稳地通过预定的路径传播。相较于直线波导，Z 型线波导的传播距离较长，出口处的振幅明显减小。在实际应用中，需要合理安排波导传播路径的长度才能达到最佳波导效果。

图18 直线型线缺陷Fig.18 Straight line defects

利用声子晶体的波导效应可将弯曲波的大部分能量引至非重要的区域，从而在声子晶体板上生成振动屏蔽区域。如图19 所示，在长度为0.35 m，宽度为0.35 m，厚度为0.5 mm 的声子晶体板中设置T型线缺陷，在点P4（0，0）处设置频率为3700 Hz，方向垂直于板，大小为1 N 的简谐载荷。由图19 可知，带隙频段内的弯曲波可沿着T 型线缺陷进行传播，并在声子晶体板内形成振动屏蔽区域（白色虚线框），且振动屏蔽区域在带隙范围内的位移远小于T型线缺陷上的位移。上述研究表明本文提出的声子晶体结构对弯曲波具有较好的振动屏蔽效果。