Processing math: 100%

内嵌声学黑洞薄板振动特性数值模拟方法研究

王博涵, 杨德庆, 夏利福

王博涵, 杨德庆, 夏利福. 内嵌声学黑洞薄板振动特性数值模拟方法研究[J]. 中国舰船研究, 2019, 14(4): 30-39. DOI: 10.19693/j.issn.1673-3185.01366
引用本文: 王博涵, 杨德庆, 夏利福. 内嵌声学黑洞薄板振动特性数值模拟方法研究[J]. 中国舰船研究, 2019, 14(4): 30-39. DOI: 10.19693/j.issn.1673-3185.01366
Wang Bohan, Yang Deqing, Xia Lifu. Study on numerical simulation method for vibration characteristics of shell embedded with acoustic black hole[J]. Chinese Journal of Ship Research, 2019, 14(4): 30-39. DOI: 10.19693/j.issn.1673-3185.01366
Citation: Wang Bohan, Yang Deqing, Xia Lifu. Study on numerical simulation method for vibration characteristics of shell embedded with acoustic black hole[J]. Chinese Journal of Ship Research, 2019, 14(4): 30-39. DOI: 10.19693/j.issn.1673-3185.01366
王博涵, 杨德庆, 夏利福. 内嵌声学黑洞薄板振动特性数值模拟方法研究[J]. 中国舰船研究, 2019, 14(4): 30-39. CSTR: 32390.14.j.issn.1673-3185.01366
引用本文: 王博涵, 杨德庆, 夏利福. 内嵌声学黑洞薄板振动特性数值模拟方法研究[J]. 中国舰船研究, 2019, 14(4): 30-39. CSTR: 32390.14.j.issn.1673-3185.01366
Wang Bohan, Yang Deqing, Xia Lifu. Study on numerical simulation method for vibration characteristics of shell embedded with acoustic black hole[J]. Chinese Journal of Ship Research, 2019, 14(4): 30-39. CSTR: 32390.14.j.issn.1673-3185.01366
Citation: Wang Bohan, Yang Deqing, Xia Lifu. Study on numerical simulation method for vibration characteristics of shell embedded with acoustic black hole[J]. Chinese Journal of Ship Research, 2019, 14(4): 30-39. CSTR: 32390.14.j.issn.1673-3185.01366

内嵌声学黑洞薄板振动特性数值模拟方法研究

基金项目: 

国家自然科学基金资助项目 51479115

详细信息
    作者简介:

    王博涵, 男, 1996年生, 硕士生

    杨德庆, 男, 1968年生, 博士, 教授, 博士生导师

    夏利福, 男, 1994年生, 硕士生

    通讯作者:

    杨德庆

  • 中图分类号: U661.44

Study on numerical simulation method for vibration characteristics of shell embedded with acoustic black hole

知识共享许可协议
内嵌声学黑洞薄板振动特性数值模拟方法研究王博涵,采用知识共享署名4.0国际许可协议进行许可。
  • 摘要:
      目的  针对内嵌二维声学黑洞薄板的数值模拟方法进行对比研究,
      方法  提出使用板壳单元对声学黑洞薄板进行离散建模的方法。首先,设计内嵌二维声学黑洞薄板算例,应用体单元建立该薄板的有限元模型,以验证其声学黑洞效应。然后,采用相同尺度的板壳单元来模拟该二维声学黑洞薄板,并对2种建模方法的模态和振动特性进行数值对比验证。最后,讨论声学黑洞区域厚度阶梯变化尺度以及幂律对整体结构振动特性的影响。
      结果  研究表明,使用板壳单元离散内嵌声学黑洞薄板计算其振动特性的精度是可靠的,板壳单元离散方法相比体单元离散法具有计算规模小、计算效率高的优点。同时,合理调节声学黑洞区域的厚度阶梯变化尺度以及幂律会改变局部模态对应的固有频率。
      结论  研究结果可为声学黑洞结构的设计与应用提供指导。
    Abstract:
      Objectives  To study a numerical simulation method for a plate embedded with a two-dimensional acoustic black hole,
      Methods  a method for discrete modeling of the structure with acoustic black hole by a shell element is proposed. Firstly, an example of a plate embedded with a two-dimensional acoustic black hole is designed in this paper, a finite element model of the plate is established with a solid element, and the acoustic black hole effect of the plate is verified; secondly, the plate embedded with two-dimensional acoustic black hole is simulated with a shell element with the same size, and the modes and vibration responses obtained by the two modeling methods are numerically compared and verified; and finally, the influence of step change scale and power law of acoustic black hole region on the vibration characteristics of the whole structure are discussed.
      Results  The results show that it is reliable to calculate the accuracy of the vibration characteristics of the shell with acoustic black hole by means of shell element dispersion, and the shell element discrete method can reduce the computation scale and improve the computation efficiency compared with the solid element discrete method. At the same time, the reasonable adjustment of the step change scale and power law of the acoustic black hole region will influence the natural frequency of local modes.
      Conclusions  The research results can provide guidance for the design and application of the structure with acoustic black hole.
  • 声学黑洞(Acoustic Black Hole,ABH)效应是指在通过合理控制面板厚度且遵循幂函数变化逐渐减小为零的过程中,面板内弯曲波传播速度逐渐减小,在理想情况下,当面板厚度为0时,波速衰减为0而波幅趋于无穷的现象。声学黑洞的概念及基本理论由Pekeris[1]在1946年提出,他发现在靠近海岸线处海水深度逐渐减少,海水中声波的相速度也会随之逐渐衰减为0,因此不发生反射。随后,Mirnov[2]在楔形结构中发现固体结构内部的弯曲波存在类似的现象。Krylov等[3-6]将声学黑洞概念应用于梁、板壳结构的减振中,并在理论、数值与实验3个层面对梁、平板内的声学黑洞效应进行了系统性研究。现有研究表明,声学黑洞区域的厚度变化若遵循幂函数衰减规律,且幂函数的指数不小于2,则弯曲波传播至声学黑洞区域时,由于结构厚度减少,波速会逐渐减小,波长变短。当声学黑洞区域接近厚度为0的位置时,累计相位趋于无穷大,弯曲波被永远局限在声学黑洞区域内。声学黑洞结构通过控制厚度变化来达到聚集结构弯曲波的效果,在船舶、航天器和车辆等结构的振动能量控制、结构减振降噪等方面有潜在应用价值。

    目前,在声学黑洞研究中,数值仿真方法的应用最为广泛。黄薇等[7]使用有限元软件建立二维声学黑洞模型,研究弯曲波在声学黑洞区域的时域传播过程,结合有限元数值结果与振动功率流的结果,分析弯曲波能量聚集过程。Conlon等[8]将多个二维声学黑洞结构以周期性规律布置在薄板结构中,结合均匀板的响应特征,计算并综合分析该薄板结构的振动和声辐射特性。Jia等[9]将二维声学黑洞板与动力吸振器相结合,在板的声学黑洞区域连接一个弹簧质量系统,分析整个板的振动响应特征。上述数值仿真模型均采用体单元进行离散。

    对于二维声学黑洞板,其厚度远小于长、宽方向的尺度。在使用体单元进行离散时,为保证单元在3个方向的尺度相差不是太大,需要在结构的长、宽方向上划分大量单元。因此,使用体单元进行数值模拟存在单元数量多、计算效率低的缺点。而二维声学黑洞区域板厚度与另外2个方向的尺度相比属于薄板的范畴,从力学特性角度分析,合理选择划分网格的尺度,使用板壳单元离散得到的计算结果也应是可靠的。

    因此,本文拟对内嵌二维声学黑洞薄板采用板壳单元模拟声学黑洞效应的方法进行深入研究,并与体单元离散建模的计算结果进行对比,验证使用板壳单元模拟声学黑洞效应的可靠性。最后讨论声学黑洞区域厚度阶梯变化尺度以及幂律对其振动性能的影响。

    二维声学黑洞结构示意图如图 1所示,其中声学黑洞区域的半径由rABH表示。二维声学黑洞区域xOz截面的厚度遵循幂函数式:

    图  1  二维声学黑洞结构示意图
    Figure  1.  The structure of two-dimensional acoustic black hole
    h(x)=εxm (1)

    式中:h(x)为距中心点x处声学黑洞结构(板)的厚度;ε为常数;幂律m为正有理数,被称为黑洞效应指数因子。Mironov[2]已证明,当结构的厚度变化满足幂函数曲线且指数不小于2时,就能满足结构声学黑洞的基本要求,将传导至声学黑洞区域的弯曲波全吸收。

    弯曲波在二维声学黑洞结构中传播的控制方程为[10]

    2x2[D(x,y)(2wx2+μ2wy2)]+2(1μ)2xy[D(x,y)2wxy]+2y2[D(x,y)(2wy2+μ2wx2)]ω2ρh(x,y)w=0 (2)

    式中:D(x,y)=Eh3(x,y)12(1μ2),为板的弯曲刚度,E为弹性模量,μ为泊松比;ρ为结构材料的密度;w为挠度;ω为圆频率。式(2)的解可以写为

    w(x,y)=A(x,y)eikpφ(x,y) (3)

    式中:A(x, y)为振幅;φ(x, y)为相位;kp为波数。将式(3)代入式(2),若使等式成立,则等式右侧的实部和虚部必须为0。由实部等于0并舍弃关于A(x, y)和φ(x, y)]的高阶导数项,即可得到

    k(x,y)=121/4k1/2p[h(x,y)]1/2 (4)

    式中,k(x, y)为与位置有关的波数。

    二维声学黑洞是由一维声学黑洞的截面h(x)=εxmm≥2)以x=0点为中心旋转1周形成的,因此,一维声学黑洞的特性可以一定程度上反映二维的情况。对于xOz截面,式(4)可写为

    k(x)=121/4k1/2p(εxm)1/2 (5)

    根据式(5),当x→0时,弯曲波波数将趋于正无穷,这意味着弯曲波波速将趋于0,波幅则趋于无限大,振动将被局限在声学黑洞中心区域。因此,理想的声学黑洞在理论上可以100%吸收板内弯曲波,取得显著的减振降噪效果。

    在实际加工制造过程中,板厚无法严格遵循幂函数连续变化至0,总会在接近黑洞区域中心存在一个截断厚度。当板厚小于截断厚度时不再继续衰减,而是维持在常量h0。这导致声学黑洞区域无法“吞噬”全部弯曲波,而会对入射波产生一定反射。虽然截断厚度不可避免,但如果将阻尼材料附加于声学黑洞区域,就能有效降低黑洞区域的弯曲波反射率[11]

    基于文献[7]中的二维声学黑洞板算例,设计了3种四边简支矩形板模型。

    模型1:四边简支均匀厚度矩形薄板,板长a=600 mm,宽b=400 mm,板厚h=5 mm,如图 2所示。板的材料为钢,密度为7 850 kg/m3,弹性模量为210 GPa,泊松比为0.3,结构阻尼系数为0.02。

    图  2  四边简支均匀厚度矩形薄板
    Figure  2.  The simply supported rectangular uniform thin plate

    模型2:在模型1基础上,板心处开出半径r=100 mm的圆孔,如图 3所示。

    图  3  中心处开孔的均匀厚度矩形薄板
    Figure  3.  The rectangular uniform thin plate with a hole in the center

    模型3:在模型1基础上,板心处制造满足幂律m=2的半径r2=100 mm的声学黑洞区域,如图 4所示。黑洞区域板厚度变化遵循幂函数hx)=(h/r22)×x2 =[5/(100)2x2。考虑加工原因,无法使中心区域的厚度无限趋近于0,这里取r1=20 mm范围内板厚为一固定小值,最小厚度h0=0.2 mm。

    图  4  中心处变厚度的二维声学黑洞板
    Figure  4.  The two-dimensional ABH plate with variable thickness in the center

    在点(-260 mm,0,0)处施加1 N的正弦激励力,激励频率为1~8 000 Hz,频率步长为10 Hz。选取图 2~图 4中所示的4个振级评价点,输出各点的振动速度级并进行对比。其中振动速度级的参考值为10-9 m/s。评价点1~4的坐标分别为(160 mm,0,0)、(230 mm,0,0)、(154.445 mm,102.014 mm,0)和(230 mm,160 mm,0)。

    采用体单元对图 4所示的声学黑洞板(模型3)进行离散,得到如图 5所示的有限元模型。该模型共有86 000个单元,108 105个节点。网格划分时,黑洞区域体单元长度为2 mm,沿板厚度方向划分4层体单元。四边简支边界条件施加在板厚度方向的中间面块体节点处。在激励点一侧选取评价点5,坐标为(-230 mm,0,0),比较载荷作用端同侧的评价点5与黑洞另一侧评价点1~4的振动速度级,以此判断声学黑洞是否能吸收弯曲波,进而减小振动响应。

    图  5  体单元离散的二维声学黑洞板有限元模型
    Figure  5.  The finite element model of two-dimensional ABH plate with solid elements

    使用MSC/Patran有限元软件对声学黑洞板进行频响分析,激励力频率范围为1~8 000 Hz,频率步长为10 Hz,计算5个评价点的振动速度级,结果如图 6所示。

    图  6  5个评价点振动速度级对比结果
    Figure  6.  Comparison of vibration speed level between five evaluation points

    图 6可见,5个评价点的振动速度级在低频段相差不大,这说明低频时声学黑洞效应并不明显。在大部分高频段,评价点1~4的振速小于评价点5的振速,这说明弯曲波在经过黑洞后,由于结构的声学黑洞效应,导致板右侧的振动减小。算例中的声学黑洞不是理想黑洞,在黑洞中心处存在截断(中心等厚度区域),不可能吸收所有的弯曲波,但是根据Krylov等[11]的研究,在黑洞中心截断区域表面敷设少量阻尼会有很好的吸收效果。

    根据弹性力学和有限元法关于板壳力学的特性理论,在计算效率和计算精度2个方面,采用板壳单元离散声学黑洞板计算应该优于等量单元数情况下采用体单元离散板结构后的相关计算。本文采用板壳单元离散声学黑洞板,并讨论其与采用体单元离散声学黑洞板的计算结果的等效性。首先介绍有限元模型的网格划分方式,然后对结构模态以及在1~8 000 Hz单位激励下板的频率响应进行对比。

    对4种有限元模型进行讨论:

    1)有限元模型(a),采用板壳单元模拟模型1所示的均匀厚度薄板(共有21 500个单元、21 621个节点);

    2)有限元模型(b),采用板壳单元模拟模型2所示中心处开孔的均匀厚度薄板(共有6 000个单元、6 220个节点);

    3)有限元模型(c),采用板壳单元(单元长度为2 mm)模拟模型3所示声学黑洞板(共有21 500个单元、21 621个节点);

    4)有限元模型(d),采用体单元(单元长度为2 mm)模拟模型3所示声学黑洞板(共有86 000个单元、108 105个节点)。

    有限元模型(a)和(c)如图 7所示,两者的区别仅在于板厚度参数的设置不同。对于有限元模型(c),黑洞区域有限元模型被离散成阶梯厚度板,单元宽度为2 mm,厚度依照幂律变化。有限元模型(b)如图 8所示。有限元模型(d)如图 5所示。

    图  7  有限元模型(a)和(c)
    Figure  7.  Finite element model (a)and (c)
    图  8  有限元模型(b)
    Figure  8.  Finite element model (b)

    根据四边简支矩形薄板固有频率的解析解公式,均匀厚度薄板前5阶固有频率为

    ωij=π2(i2a2+j2b2)Dρh (6)

    式中:ij分别为薄板在长度和宽度方向上的半波数。

    使用MSC/Patran有限元软件对上述4种模型结构进行模态计算,其前5阶固有频率如表 1所示。板壳单元离散的声学黑洞板的前5阶振型图如图 9所示,体单元离散的声学黑洞板的前5阶振型图如图 10所示。

    表  1  4种模型的前5阶固有频率
    Table  1.  First five order natural frequencies of four models
    振型 薄板解析解/Hz 薄板数值解/Hz 开孔板数值解/Hz 声学黑洞板(板壳单元)/Hz 声学黑洞板(体单元)/Hz
    (1,1) 110.98 110.56 113.61 101.21 110.81
    (2,1) 213.42 212.43 191.18 194.46 199.00
    (1,2) 341.47 340.59 297.17 290.58 302.25
    (3,1) 384.15 382.52 423.44 373.60 391.94
    (2,2) 443.91 441.51 411.09 416.46 423.93
    下载: 导出CSV 
    | 显示表格
    图  9  板壳单元离散的声学黑洞板前5阶振型图
    Figure  9.  The first five modes of ABH plate with shell elements
    图  10  体单元离散的声学黑洞板前5阶振型图
    Figure  10.  The first five modes of ABH plate with solid elements

    表 1可见,薄板前5阶固有频率的解析解和数值解相吻合,误差小于1%;此外,声学黑洞板的固有频率介于薄板与开孔板两者之间或与两者接近。

    对比使用板壳单元和体单元模拟声学黑洞板的固有频率计算结果,可以看到,两者非常接近,板壳单元离散声学黑洞板的固有频率计算结果普遍小于同阶次体单元离散的计算结果,相对误差为5%~9%。由图 9图 10可以看出,2种模拟方法的前5阶振型图基本相同。至于(3,1)阶振型的情况,本文关注的是结构的整体模态,体单元建模振型图的中心区域与板壳单元的不同,是因为此时钢板的整体模态图形显示受到了局部模态图形显示比例尺度的影响。总之,计算结果验证了板壳单元建模在模态计算方面的合理性,且计算误差较小。

    施加前文所述激励力进行频响分析,激励力频率范围为1~8 000 Hz,频率步长为10 Hz。对有限元模型(a),(c)和(d)进行频响计算,计算结果如图 11所示。

    图  11  不同模拟方法的4个评价点振动速度级
    Figure  11.  Vibration velocity level of four evaluation points in the case of different simulation methods

    图 11可见,采用板壳单元和体单元模拟声学黑洞板,在低频段,板壳单元和体单元离散的2种模型评价点的计算结果较吻合;在高频段,因固有频率存在一定差异,两者响应计算结果的峰值位置有一定差异,但峰值十分接近。这说明采用板壳单元模拟声学黑洞板的计算结果也是精确可靠的。

    图 11还反映了声学黑洞板的减振性能。采用两种方式建模的声学黑洞板与薄板(模型(a))相比,三者在低频段的振动速度级基本吻合,这说明声学黑洞板在低频段减振效果不明显。对于高频段,可以看到声学黑洞板的减振效果很明显,与文献[10]的研究结果相同,说明采用板壳单元模拟声学黑洞板时,其频响计算结果也体现了黑洞效应。

    对黑洞区域离散时,板厚度阶梯变化尺度分别取5,10和16 mm进行建模,分析黑洞区域厚度阶梯变化尺度对声学黑洞板振动性能的影响。

    声学黑洞区域板厚分布不均匀,满足幂律变化规律,厚度在板中心处逐步变薄。因此,模态分析时在黑洞区域会出现局部模态,在某些特定频率,黑洞中心区域会发生共振,吸收能量。图 12给出了厚度阶梯变化尺度为2 mm时黑洞区域的一阶局部模态。

    图  12  厚度阶梯变化尺度为2 mm时黑洞区域一阶局部模态振型
    Figure  12.  The first order local mode with 2 mm step change scale

    研究发现,厚度阶梯变化尺度的改变会导致黑洞区域局部共振频率发生变化。厚度阶梯变化尺度分别为2,5,10和16 mm的一阶局部模态对应的固有频率如表 2所示。

    表  2  不同厚度阶梯变化尺度一阶局部模态对应的固有频率
    Table  2.  Natural frequencies of first order local mode with different step change scales
    厚度阶梯变化尺度/mm 一阶局部固有频率/Hz
    2 625.95
    5 568.59
    10 452.94
    16 313.09
    下载: 导出CSV 
    | 显示表格

    表 2可见,厚度阶梯变化尺度对局部模态对应的固有频率产生较大影响,且随着厚度阶梯变化尺度的增大,局部模态对应的固有频率会大幅降低。在实际应用中,局部模态可以起到吸收振动能量的作用,因此合理调整厚度阶梯变化尺度,可以将声学黑洞板的局部模态对应的固有频率调整到预期减振频率。

    施加前文所述激励力,对黑洞区域厚度阶梯变化尺度分别为2,5,10和16 mm的声学黑洞板进行频响分析,计算结果如图 13所示。

    图  13  不同厚度阶梯变化尺度情况下4个评价点的振动速度级
    Figure  13.  Vibration velocity level of four evaluation points in the case of different step change scales

    分析图 13可知,低频时不同厚度阶梯变化尺度声学黑洞板的速度响应吻合较好,与无声学黑洞薄板振速相比,没有明显的减振效果。高频时不同厚度阶梯变化尺度的声学黑洞板均具有声学黑洞效应,有明显的减振效果,且其减振频率范围不同。由此说明,合理调整厚度阶梯变化尺度可以对目标频段进行减振设计。

    在幂律m=3,4的情况下,黑洞区域厚度阶梯变化尺度取2 mm,对板壳单元和体单元离散声学黑洞板进行模态分析,并与m=2声学黑洞板的计算结果进行比较。幂律变化对声学黑洞板整体模态特性的影响如表 3所示。

    表  3  幂律变化对声学黑洞板前5阶整体固有频率的影响
    Table  3.  Influence of power law variation on the first five order natural frequencies of ABH plate
    振型 m=2 m=3 m=4
    板壳单元离散/Hz 体单元离散/Hz 板壳单元离散/Hz 体单元离散/Hz 板壳单元离散/Hz 体单元离散/Hz
    (1,1) 101.21 110.81 102.93 111.62 104.46 111.97
    (2,1) 194.46 199.00 190.87 196.31 189.52 195.05
    (1,2) 290.58 302.25 283.95 297.35 283.18 295.98
    (3,1) 373.60 391.94 382.88 398.43 393.64 402.66
    (2,2) 416.46 423.93 412.31 421.37 410.43 419.92
    下载: 导出CSV 
    | 显示表格

    表 3可以看出,幂律变化对声学黑洞板整体固有频率的影响很小,无论采用板壳单元还是体单元离散声学黑洞板,固有频率的计算结果基本一致。在幂律相同的情况下,使用体单元离散的声学黑洞板的整体固有频率要高于板壳单元离散的情况。

    幂律变化对黑洞区域局部模态特性的影响如表 4所示。其中,m=2时,板壳单元离散模型的前3阶局部模态振型如图 14所示,对于m=3,4的情况,其前3阶局部模态振型与m=2的情况类似,在此不再赘述,仅列出具体的固有频率值。

    表  4  幂律变化对声学黑洞板前3阶局部固有频率的影响
    Table  4.  Influence of power law variation on the first three order local natural frequencies of ABH plate
    阶数 m=2 m=3 m=4
    板壳单元离散/ Hz 体单元离散/ Hz 板壳单元离散/ Hz 体单元离散/ Hz 板壳单元离散/ Hz 体单元离散/ Hz
    1 625.95 762.14 168.69 199.79 37.04 50.415
    2 938.44 1 011.7 328.04 372.57 74.802 100.55
    3 981.78 1 055.3 328.49 372.68 74.802 100.56
    下载: 导出CSV 
    | 显示表格
    图  14  板壳单元离散模型的前3阶局部模态振型图(m=2)
    Figure  14.  The first three order local mode shapes with shell elements (m=2)

    以上计算结果表明,当幂律越大时,由于声学黑洞板黑洞区域厚度变化较为剧烈,黑洞中心区域的板厚相比邻近区域更薄,导致黑洞区域的局部共振模态在频率很低的情况下就会出现。由表 4可见,m=4时,一阶局部模态在37.04 Hz出现;m=3时,一阶局部模态在168.69 Hz出现。对于高幂律的情况,与m=2相比,低频的局部模态十分密集。模态计算结果还表明,采用板壳单元和体单元模拟局部模态的结果存在一定差异,板壳单元模拟的局部模态结果要低于体单元模拟的情况。

    对幂律m=2,3,4的声学黑洞板分别进行频响分析,网格划分时黑洞区域厚度阶梯变化尺度为2 mm,计算结果如图 15所示。

    图  15  不同幂律下4个评价点的振动速度级
    Figure  15.  Vibration velocity level of four evaluation points in the case of different power laws

    图 15所示的频响分析结果表明,低频段时三者相差不大,黑洞效应不明显。在高频段,不同幂律声学黑洞板的减振频段不同,例如m=4的声学黑洞板在3 000~4 000 Hz频段的减振效果优于m=2的声学黑洞板。

    针对内嵌二维声学黑洞薄板,本文对采用板壳单元模拟声学黑洞效应的计算方法进行数值实验计算,并与体单元的模拟结果进行了对比,验证板壳单元模拟方法的有效性。探讨了声学黑洞板厚度阶梯变化尺度以及幂律对其振动性能的影响,为声学黑洞板的设计与实际应用提供指导。研究表明:

    1)采用板壳单元建模方法,声学黑洞板整体模态以及频率响应计算结果与体单元建模方法的基本相同,模态计算结果误差小于10%。板壳单元建模大大减少了模型单元数量,提高了计算效率。这证明了采用板壳单元模拟声学黑洞板的合理性和有效性。

    2)声学黑洞区域的厚度阶梯变化尺度以及幂律的变化可以调节黑洞区域局部模态对应的固有频率,影响具体频段的减振效果。在实际应用中,可以利用该特点对目标频段进行减振设计。

  • 图  1   二维声学黑洞结构示意图

    Figure  1.   The structure of two-dimensional acoustic black hole

    图  2   四边简支均匀厚度矩形薄板

    Figure  2.   The simply supported rectangular uniform thin plate

    图  3   中心处开孔的均匀厚度矩形薄板

    Figure  3.   The rectangular uniform thin plate with a hole in the center

    图  4   中心处变厚度的二维声学黑洞板

    Figure  4.   The two-dimensional ABH plate with variable thickness in the center

    图  5   体单元离散的二维声学黑洞板有限元模型

    Figure  5.   The finite element model of two-dimensional ABH plate with solid elements

    图  6   5个评价点振动速度级对比结果

    Figure  6.   Comparison of vibration speed level between five evaluation points

    图  7   有限元模型(a)和(c)

    Figure  7.   Finite element model (a)and (c)

    图  8   有限元模型(b)

    Figure  8.   Finite element model (b)

    图  9   板壳单元离散的声学黑洞板前5阶振型图

    Figure  9.   The first five modes of ABH plate with shell elements

    图  10   体单元离散的声学黑洞板前5阶振型图

    Figure  10.   The first five modes of ABH plate with solid elements

    图  11   不同模拟方法的4个评价点振动速度级

    Figure  11.   Vibration velocity level of four evaluation points in the case of different simulation methods

    图  12   厚度阶梯变化尺度为2 mm时黑洞区域一阶局部模态振型

    Figure  12.   The first order local mode with 2 mm step change scale

    图  13   不同厚度阶梯变化尺度情况下4个评价点的振动速度级

    Figure  13.   Vibration velocity level of four evaluation points in the case of different step change scales

    图  14   板壳单元离散模型的前3阶局部模态振型图(m=2)

    Figure  14.   The first three order local mode shapes with shell elements (m=2)

    图  15   不同幂律下4个评价点的振动速度级

    Figure  15.   Vibration velocity level of four evaluation points in the case of different power laws

    表  1   4种模型的前5阶固有频率

    Table  1   First five order natural frequencies of four models

    振型 薄板解析解/Hz 薄板数值解/Hz 开孔板数值解/Hz 声学黑洞板(板壳单元)/Hz 声学黑洞板(体单元)/Hz
    (1,1) 110.98 110.56 113.61 101.21 110.81
    (2,1) 213.42 212.43 191.18 194.46 199.00
    (1,2) 341.47 340.59 297.17 290.58 302.25
    (3,1) 384.15 382.52 423.44 373.60 391.94
    (2,2) 443.91 441.51 411.09 416.46 423.93
    下载: 导出CSV

    表  2   不同厚度阶梯变化尺度一阶局部模态对应的固有频率

    Table  2   Natural frequencies of first order local mode with different step change scales

    厚度阶梯变化尺度/mm 一阶局部固有频率/Hz
    2 625.95
    5 568.59
    10 452.94
    16 313.09
    下载: 导出CSV

    表  3   幂律变化对声学黑洞板前5阶整体固有频率的影响

    Table  3   Influence of power law variation on the first five order natural frequencies of ABH plate

    振型 m=2 m=3 m=4
    板壳单元离散/Hz 体单元离散/Hz 板壳单元离散/Hz 体单元离散/Hz 板壳单元离散/Hz 体单元离散/Hz
    (1,1) 101.21 110.81 102.93 111.62 104.46 111.97
    (2,1) 194.46 199.00 190.87 196.31 189.52 195.05
    (1,2) 290.58 302.25 283.95 297.35 283.18 295.98
    (3,1) 373.60 391.94 382.88 398.43 393.64 402.66
    (2,2) 416.46 423.93 412.31 421.37 410.43 419.92
    下载: 导出CSV

    表  4   幂律变化对声学黑洞板前3阶局部固有频率的影响

    Table  4   Influence of power law variation on the first three order local natural frequencies of ABH plate

    阶数 m=2 m=3 m=4
    板壳单元离散/ Hz 体单元离散/ Hz 板壳单元离散/ Hz 体单元离散/ Hz 板壳单元离散/ Hz 体单元离散/ Hz
    1 625.95 762.14 168.69 199.79 37.04 50.415
    2 938.44 1 011.7 328.04 372.57 74.802 100.55
    3 981.78 1 055.3 328.49 372.68 74.802 100.56
    下载: 导出CSV
  • [1]

    Pekeris C L. Theory of propagation of sound in a half-space of variable sound velocity under conditions of formation of a shadow zone[J]. The Journal of the Acoustical Society of America, 1946, 18(2):295-315. doi: 10.1121/1.1916366

    [2]

    Mironov M A. Propagation of a flexural wave in a plate whose thickness decreases smoothly to zero in a finite interval[J]. Soviet Physics Acoustics-USSR, 1988, 34(3):318-319. https://www.researchgate.net/publication/284775645_Propagation_of_a_flexural_wave_in_a_plate_whose_thickness_decreases_smoothly_to_zero_in_a_finite_interval

    [3]

    Krylov V V. Laminated plates of variable thickness as effective absorbers for flexural vibrations[C]//Proceedings of the 17th International Congress on Acoustics. Rome, Italy: [s.n], 2001.

    [4]

    Krylov V V, Tilman F J B S. Acoustic 'black holes' for flexural waves as effective vibration dampers[J]. Journal of Sound and Vibration, 2004, 274(3/4/5):605-619. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=65a6d33de221a56fcca0fea617a2e651

    [5]

    Krylov V V. New type of vibration dampers utilizing the effect of acoustic 'black holes'[J]. Acta Acustica United with Acustica, 2004, 90(5):830-837. http://cn.bing.com/academic/profile?id=c37231ce6e501515780873b21e9da23d&encoded=0&v=paper_preview&mkt=zh-cn

    [6]

    Bowyer E P, Krylov V V. Experimental investigation of damping flexural vibrations in glass fiber composite plates containing one-and two-dimensional acoustic black holes[J]. Composite Structures, 2014, 107:406-415. doi: 10.1016/j.compstruct.2013.08.011

    [7] 黄薇, 季宏丽, 裘进浩, 等.二维声学黑洞对弯曲波的能量聚集效应[J].振动与冲击, 2017, 36(9):51-57. http://d.old.wanfangdata.com.cn/Periodical/zdycj201709008

    Huang W, Ji H L, Qiu J H, et al. Energy focusing effect of two-dimensional acoustic black hole on flexural waves[J]. Journal of Vibration and Shock, 2017, 36(9):51-57(in Chinese). http://d.old.wanfangdata.com.cn/Periodical/zdycj201709008

    [8]

    Conlon S C, Fahnline J B, Shepherd M R, et al. Vibration control using grids of acoustic black holes: how many is enough?[C]//Inter-noise and Noise-con Congress and Conference Proceedings. San Francisco, USA: Institute of Noise Control Engineering, 2015.

    [9]

    Jia X X, Du Y, Zhao K M. Vibration control of variable thickness plates with embedded acoustic black holes and dynamic vibration absorbers[C]//Proceedings of ASME 2015 Noise Control and Acoustics Division Conference at InterNoise 2015. San Francisco, California, USA: American Society of Mechanical Engineers, 2015.

    [10] 季宏丽, 黄薇, 裘进浩, 等.声学黑洞结构应用中的力学问题[J].力学进展, 2017, 47:201710. doi: 10.6052/1000-0992-16-033

    Ji H L, Huang W, Qiu J H, et al. Mechanics problems in application of acoustic black hole structures[J]. Advances in Mechanics, 2017, 47:201710(in Chinese). doi: 10.6052/1000-0992-16-033

    [11]

    Krylov V V, Winward R E T B. Experimental investigation of the acoustic black hole effect for flexural waves in tapered plates[J]. Journal of Sound and Vibration, 2007, 300(1/2):43-49. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=042862cb75b2b04cc0fb2b59a2d5c9cb

  • 期刊类型引用(9)

    1. 魏海婴,温华兵,康钦伟,黄惠文,史自强,吕世金. 约束阻尼对二维声学黑洞薄板减振特性的影响. 噪声与振动控制. 2023(02): 253-258+291 . 百度学术
    2. 余勇花,吴天歌,杨淇,叶锦啸,夏兆旺. 基于声学黑洞的复合隔声结构声学特性研究及试验. 声学技术. 2023(04): 509-514 . 百度学术
    3. 魏建红,向阳,廖琳,张春吉. 声学黑洞在圆柱壳上的应用及结构优化研究. 噪声与振动控制. 2023(05): 59-65 . 百度学术
    4. 赵晓宇,吴卫国,林永水. 波导吸振器的低频减振设计和应用. 力学学报. 2023(11): 2636-2646 . 百度学术
    5. 赵楠,王禹,陈林,庞福振. 分布式声学黑洞浮筏系统隔振性能研究. 振动与冲击. 2022(13): 75-80 . 百度学术
    6. 李敬,万志威,李天匀,朱翔. 周期声学黑洞结构弯曲波带隙与振动特性. 噪声与振动控制. 2021(02): 21-27 . 百度学术
    7. 颜伏伍,张家铭,李良栋,朱亚伟,刘志恩. 基于二维声学黑洞结构的面板振动噪声控制理论. 塑性工程学报. 2021(07): 184-192 . 百度学术
    8. 曾周末,李广智,黄新敬,李健,杜立普. 基于倒置声学黑洞的弹性波汇聚检测研究. 仪器仪表学报. 2021(09): 151-159 . 百度学术
    9. 赵业楠,杨德庆,王博涵. 声学黑洞俘能器在气垫船舱室噪声控制中的应用研究. 中国造船. 2020(03): 58-67 . 百度学术

    其他类型引用(13)

图(15)  /  表(4)
计量
  • 文章访问数:  538
  • HTML全文浏览量:  201
  • PDF下载量:  102
  • 被引次数: 22
出版历程
  • 收稿日期:  2018-07-29
  • 网络出版日期:  2021-05-07
  • 刊出日期:  2019-08-20

目录

/

返回文章
返回