• ISSN 1001-1455  CN 51-1148/O3
  • EI Compendex、CA收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

超材料混凝土的带隙特征及对冲击波的衰减效应

张恩 路国运 杨会伟 曹瑞东 陈鹏程

引用本文:
Citation:

超材料混凝土的带隙特征及对冲击波的衰减效应

    作者简介: 张 恩(1991- ),男,博士研究生,18734862993@163.com;
    通讯作者: 路国运, luguoyun@tyut.edu.cn
  • 中图分类号: O382

Band gap features of metaconcrete and shock wave attenuation in it

    Corresponding author: LU Guoyun, luguoyun@tyut.edu.cn ;
  • CLC number: O382

  • 摘要: 借鉴超材料的研究思路,在混凝土中引入谐振骨料,设计出具有消波特性的超材料混凝土。首先,通过结构动力学方法计算超材料混凝土的有效质量,从而建立了超材料混凝土带隙起始频率及截止频率的简化模型,并给出了带隙起始频率及截止频率的理论表达式。然后,分析了涂层弹性模量、芯柱密度、基体密度、骨料体积占比和芯柱边长与软涂层厚度比对超材料混凝土带隙特征的影响。最后,采用数值模拟的方法,对比了超材料混凝土和普通混凝土对冲击波的衰减效应。研究结果表明:(1)低弹性模量涂层能够形成低频带隙,但带隙宽度较窄,而高弹性模量涂层能够形成较宽的带隙,但带隙起始频率较高;(2)通过选择高密度芯柱材料和低密度基体材料,可以得到低频、宽带隙特征;(3)通过增大骨料体积占比和芯柱边长与软涂层厚度比可以实现扩宽带隙的目的;(4)与普通混凝土相比,超材料混凝土对冲击波具有更好的衰减作用。
  • 图 1  本文中提出的超材料混凝土

    Figure 1.  Metaconcrete designed in this paper

    图 2  超材料混凝土单胞及一维晶格系统

    Figure 2.  A unit cell and one-dimensional lattice systems for metaconcrete

    图 3  有效质量与激振频率的关系

    Figure 3.  Relationship between effective mass and excitation frequency

    图 4  带隙起始频率及截止频率的简化模型[22]

    Figure 4.  Simplified models for band gap start and cutoff frequencies[22]

    图 5  涂层弹性模量对带隙特性的影响

    Figure 5.  Effects of coating elastic modulus on band gap

    图 6  芯柱密度对带隙特征的影响

    Figure 6.  Effects of core column density on band gap

    图 7  基体密度对带隙特征的影响

    Figure 7.  Effects of matrix density on band gap

    图 8  γ=10时,骨料体积占比对带隙特征的影响

    Figure 8.  Effect of volume fraction of aggregate on band gap when γ=10

    图 9  β=0.36时,γ对带隙特征的影响

    Figure 9.  Effect of γ on band gap when β=0.36

    图 10  有限元模型俯视图

    Figure 10.  The top views of finite element models

    图 11  冲击波压力-时间曲线

    Figure 11.  Pressure-time curve of shock wave

    图 12  各截面平均应力时程曲线

    Figure 12.  Average stress-time curve at each section

  • [1] WEGENER M. Metamaterials beyond optics [J]. Science, 2013, 342(6161): 939–940. DOI: 10.1126/science.1246545.
    [2] LANDY N, SMITH D R. A full-parameter unidirectional metamaterial cloak for microwaves [J]. Nature Materials, 2013, 12(1): 25–28. DOI: 10.1038/nmat3476.
    [3] LI B, TAN K T, CHRISTENSEN J. Heat conduction tuning by hyperbranched nanophononic metamaterials [J]. Journal of Applied Physics, 2018, 123(20): 205105. DOI: 10.1063/1.5023487.
    [4] LI B, TAN K T, CHRISTENSEN J. Tailoring the thermal conductivity in nanophononic metamaterials [J]. Physical Review B, 2017, 95(14): 144305. DOI: 10.1103/PhysRevB.95.144305.
    [5] WANG Y Z, LI F M, WANG Y S. Influences of active control on elastic wave propagation in a weakly nonlinear phononic crystal with a monoatomic lattice chain [J]. International Journal of Mechanical Sciences, 2016, 106: 357–362. DOI: 10.1016/j.ijmecsci.2015.12.004.
    [6] BANERJEE B, NAGY P B. An introduction to metamaterials and waves in composites [J]. Journal of the Acoustical Society of America, 2012, 131(2): 1665. DOI: 10.1121/1.3672699.
    [7] AL BA'BA'A H B, NOUH M. Mechanics of longitudinal and flexural locally resonant elastic metamaterials using a structural power flow approach [J]. International Journal of Mechanical Sciences, 2017, 122: 341–354. DOI: 10.1016/j.ijmecsci.2017.01.034.
    [8] LIU X N, HU G K, HUANG G L, et al. An elastic metamaterial with simultaneously negative mass density and bulk modulus [J]. Applied Physics Letters, 2011, 98(25): 251907. DOI: 10.1063/1.3597651.
    [9] LI J, CHAN C T. Double-negative acoustic metamaterial [J]. Physical Review E, 2004, 70(5): 055602(R). DOI: 10.1103/PhysRevE.70.055602.
    [10] LIU Z Y, ZHANG X X, MAO Y W, et al. Locally resonant sonic materials [J]. Science, 2000, 289(5485): 1734–1736. DOI: 10.1126/science.289.5485.1734.
    [11] LIU Z Y, CHAN C T, SHENG P, et al. Elastic wave scattering by periodic structures of spherical objects: theory and experiment [J]. Physical Review B, 2000, 62(4): 2446–2457. DOI: 10.1103/PhysRevB.62.2446.
    [12] LIU Z Y, CHAN C T, SHENG P. Three-component elastic wave band-gap material [J]. Physical Review B, 2002, 65(16): 165116. DOI: 10.1103/PhysRevB.65.165116.
    [13] LIU Z Y, CHAN C T, SHENG P. Analytic model of phononic crystals with local resonances [J]. Physical Review B, 2005, 71(1): 014103. DOI: 10.1103/PhysRevB.71.014103.
    [14] YANG Z, DAI H M, CHAN N H, et al. Acoustic metamaterial panels for sound attenuation in the 50−1 000 Hz regime [J]. Applied Physics Letters, 2010, 96(4): 041906. DOI: 10.1063/1.3299007.
    [15] YANG Z, MEI J, YANG M, et al. Membrane-type acoustic metamaterial with negative dynamic mass [J]. Physical Review Letters, 2008, 101(20): 204301. DOI: 10.1103/PhysRevLett.101.204301.
    [16] 吴健, 白晓春, 肖勇, 等. 一种多频局域共振型声子晶体板的低频带隙与减振特性 [J]. 物理学报, 2016, 65(6): 064602. DOI: 10.7498/aps.65.064602.
    WU J, BAI X C, XIAO Y, et al. Low frequency band gaps and vibration reduction properties of a multi-frequency locally resonant phononic plate [J]. Acta Physica Sinica, 2016, 65(6): 064602. DOI: 10.7498/aps.65.064602.
    [17] 张印, 尹剑飞, 温激鸿, 等. 基于质量放大局域共振型声子晶体的低频减振设计 [J]. 振动与冲击, 2016, 35(17): 26–32. DOI: 10.13465/j.cnki.jvs.2016.17.005.
    ZHANG Y, YIN J F, WEN J H, et al. Low frequency vibration reduction design for inertial local resonance phononic crystals based on inertial amplification [J]. Journal of Vibration and Shock, 2016, 35(17): 26–32. DOI: 10.13465/j.cnki.jvs.2016.17.005.
    [18] BRÛLÉ S, JAVELAUD E H, ENOCH S, et al. Experiments on seismic metamaterials: molding surface waves [J]. Physical Review Letters, 2014, 112(13): 133901. DOI: 10.1103/PhysRevLett.112.133901.
    [19] 郜英杰, 范华林, 张蓓, 等. 超材料消波混凝土板在二维平面波作用下的削波效应研究 [J]. 振动与冲击, 2018, 37(20): 39–44. DOI: 10.13465/j.cnki.jvs.2016.17.005.
    GAO Y J, FAN H L, ZHANG B, et al. Wave attenuation of super-material wave absorbing concrete panel subjected to two-dimensional plane wave [J]. Journal of Vibration and Shock, 2018, 37(20): 39–44. DOI: 10.13465/j.cnki.jvs.2016.17.005.
    [20] LI Q Q, HE Z C, LI E, et al. Design of a multi-resonator metamaterial for mitigating impact force [J]. Journal of Applied Physics, 2019, 125(3): 035104. DOI: 10.1063/1.5029946.
    [21] HUANG H H, SUN C T, HUANG G L. On the negative effective mass density in acoustic metamaterials [J]. International Journal of Engineering Science, 2009, 47(4): 610–617. DOI: 10.1016/j.ijengsci.2008.12.007.
    [22] SHENG X, ZHAO C Y, YI Q, et al. Engineered metabarrier as shield from longitudinal waves: band gap properties and optimization mechanisms [J]. Journal of Zhejiang University: Science A: Applied Physics and Engineering, 2018, 19(9): 663–675. DOI: 10.1631/jzus.A1700192.
    [23] MITCHELL S J, PANDOLFI A, ORTIZ M. Metaconcrete: designed aggregates to enhance dynamic performance [J]. Journal of the Mechanics and Physics of Solids, 2014, 65: 69–81. DOI: 10.1016/j.jmps.2014.01.003.
  • [1] 孙伟王建国姜绍飞 . 减小高速列车荷载对地面结构损害的数值分析. 爆炸与冲击, 2009, 29(2): 155-161. doi: 10.11883/1001-1455(2009)02-0155-07
    [2] 黄家蓉刘瑞朝何翔孙桂娟徐鹏 . 侵彻过载测试信号的数据处理方法. 爆炸与冲击, 2009, 29(5): 555-560. doi: 10.11883/1001-1455(2009)05-0555-06
    [3] 赵秦晋李伯松 . 冲击响应谱频率特征值法研究. 爆炸与冲击, 1989, 9(1): 43-48.
    [4] 刘爱文俞言祥付长华吕红山陈鲲赵纪生周正华王伟 . 50t科学爆破振动的衰减特征及地形效应. 爆炸与冲击, 2010, 30(1): 21-26. doi: 10.11883/1001-1455(2010)01-0021-06
    [5] 刘波杨黎明李东杰欧阳科 . 侵彻弹体结构纵向振动频率特性分析. 爆炸与冲击, 2018, 38(3): 677-682. doi: 10.11883/bzycj-2016-0282
    [6] 刘义佳卢文波陈明严鹏 . 结构爆破振动响应的频率与持续时间依赖性分析. 爆炸与冲击, 2019, 39(8): 085203-1-085203-11. doi: 10.11883/bzycj-2019-0142
    [7] 卢强王占江丁洋刘晓新郭志昀吴玉蛟 . 线黏弹性球面发散应力波的频率响应特性. 爆炸与冲击, 2017, 37(6): 1023-1030. doi: 10.11883/1001-1455(2017)06-1023-08
    [8] 武旭张云鹏郭奇峰 . 台阶地形爆破振动放大与衰减效应研究. 爆炸与冲击, 2017, 37(6): 1017-1022. doi: 10.11883/1001-1455(2017)06-1017-06
    [9] 周佩杰王坚陶钢周杰 . 泡沫材料对冲击波的衰减特性. 爆炸与冲击, 2015, 35(5): 675-681. doi: 10.11883/1001-1455(2015)05-0675-07
    [10] 李铮朱瑞赓胡再龙陈宜明杨静岫蔡金祥 . 爆炸地震波振速的特征系数与衰减指数的研究. 爆炸与冲击, 1986, 6(3): 221-229.
    [11] 高启栋卢文波杨招伟严鹏陈明 . 水平光面爆破激发地震波的成分及衰减特征. 爆炸与冲击, 2019, 39(8): 085201-1-085201-13. doi: 10.11883/bzycj-2018-0280
    [12] 李景云 . 关于起始爆轰波形的测试问题. 爆炸与冲击, 1984, 4(4): 46-53.
    [13] 胡时胜王悟潘艺周璟 . 泡沫材料的应变率效应. 爆炸与冲击, 2003, 23(1): 13-18.
    [14] 严成欧卓成段卓平黄风雷 . 脆性材料动态强度应变率效应. 爆炸与冲击, 2011, 31(4): 423-427. doi: 10.11883/1001-1455(2011)04-0423-05
    [15] 田杰胡时胜蒋家桥黄西成 . 泡沫金属材料的孔径效应. 爆炸与冲击, 2004, 24(4): 347-351.
    [16] 巫绪涛廖礼 . 脆性材料中应力波衰减规律与层裂实验设计的数值模拟. 爆炸与冲击, 2017, 37(4): 705-711. doi: 10.11883/1001-1455(2017)04-0705-07
    [17] 张磊胡时胜陈德兴张守保余泽清刘飞 . 钢纤维混凝土的层裂特征. 爆炸与冲击, 2009, 29(2): 119-124. doi: 10.11883/1001-1455(2009)02-0119-06
    [18] 孟庆昌张志宏李启杰 . 高速射弹超空泡流动的重力和压缩性效应. 爆炸与冲击, 2016, 36(6): 781-788. doi: 10.11883/1001-1455(2016)06-0781-08
    [19] 王永刚胡时胜王礼立 . 爆炸荷载下泡沫铝材料中冲击波衰减特性的实验和数值模拟研究. 爆炸与冲击, 2003, 23(6): 516-522.
    [20] 梁龙河王政曹菊珍 . 长杆弹对混凝土的侵爆效应. 爆炸与冲击, 2008, 28(5): 415-420. doi: 10.11883/1001-1455(2008)05-0415-06
  • 加载中
图(12)
计量
  • 文章访问数:  1754
  • HTML全文浏览量:  938
  • PDF下载量:  22
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-06-25
  • 录用日期:  2020-04-26
  • 网络出版日期:  2020-05-25
  • 刊出日期:  2020-06-01

超材料混凝土的带隙特征及对冲击波的衰减效应

    作者简介:张 恩(1991- ),男,博士研究生,18734862993@163.com
    通讯作者: 路国运, luguoyun@tyut.edu.cn
  • 1. 太原理工大学土木工程学院,山西 太原 030024
  • 2. 太原理工大学材料强度与结构冲击山西省重点实验室,山西 太原 030024
  • 3. 山西大学土木工程系,山西 太原 030013

摘要: 借鉴超材料的研究思路,在混凝土中引入谐振骨料,设计出具有消波特性的超材料混凝土。首先,通过结构动力学方法计算超材料混凝土的有效质量,从而建立了超材料混凝土带隙起始频率及截止频率的简化模型,并给出了带隙起始频率及截止频率的理论表达式。然后,分析了涂层弹性模量、芯柱密度、基体密度、骨料体积占比和芯柱边长与软涂层厚度比对超材料混凝土带隙特征的影响。最后,采用数值模拟的方法,对比了超材料混凝土和普通混凝土对冲击波的衰减效应。研究结果表明:(1)低弹性模量涂层能够形成低频带隙,但带隙宽度较窄,而高弹性模量涂层能够形成较宽的带隙,但带隙起始频率较高;(2)通过选择高密度芯柱材料和低密度基体材料,可以得到低频、宽带隙特征;(3)通过增大骨料体积占比和芯柱边长与软涂层厚度比可以实现扩宽带隙的目的;(4)与普通混凝土相比,超材料混凝土对冲击波具有更好的衰减作用。

English Abstract

  • 近年来关于人造具有特殊功能材料的设计和制造成为研究热点,超材料(metamaterials)是人造复合材料,具有天然材料不具备的非凡的物理性质。这些特性源于其特殊设计的微观结构,而不是化学成分。超材料的出现代表着一种崭新的材料设计理念,即人类可以在利用现有材料的基础上,按照自己的意志逆向设计微结构来获得新材料。超材料的研究主要集中在电磁超材料[1-2]、热学超材料[3-4]、光学超材料[5-6]和声学超材料[7-13]等领域。

    Liu等[10-13]提出了局域共振声子晶体,使用橡胶包裹铅球构成局域共振单元。在环氧树脂基体中周期性地排列上述局域共振单元,可以构成人工周期材料。上述周期材料突破Bragg机理的限制,使用弹性波的局域共振效应实现了在亚波长频段对低频弹性波的衰减,为声学超材料的研究开拓了新思路。Yang等[14-15]使用局域共振的思想设计了薄膜型主动声学超材料,能够对50~1 000 Hz的低频声波进行阻隔,并进一步探索出有效拓宽局域共振声学超材料带隙的方法。吴健等[16]、张印等[17]使用局域共振的思想设计了梁板类声学超材料,并对其带隙计算方法进行改进和完善,最后把该类声学超材料应用于声隐身和减振降噪等领域。Brûlé等[18]把局域共振的思想引入到防震工程中,其通过在黏土中打出中心间距为1.73 m、深度为5 m、直径为0.32 m的自稳定孔,实现了对频率为50 Hz的弹性波的阻断。

    本文中,借鉴超材料的研究思路,使用局域共振骨料代替混凝土中的普通骨料,设计一种具有冲击波衰减效应的超材料混凝土。首先,基于结构动力学计算了超材料混凝土的有效质量,并进一步建立超材料混凝土带隙起始频率、截止频率的理论模型,讨论涂层弹性模量、芯柱密度、基体密度、骨料体积占比和芯柱边长与软涂层厚度比等因素对超材料混凝土带隙起始频率、截止频率及带隙宽度的影响;其次,采用数值模拟的方法,对比超材料混凝土和普通混凝土对冲击波的衰减效应。以期研究结果可为超材料混凝土的设计及应用提供理论支持。

    • 本文中提出的超材料混凝土(metaconcrete)如图1所示,由局域共振型骨料和砂浆两部分组成。其中局域共振型骨料由高密度柱体周围包覆软涂层形成,把该骨料均匀地排列在砂浆中,形成超材料混凝土。最初,局域共振型结构多被设计成球形和圆柱形结构[10-13],但是在实际使用中,球形和圆柱形结构存在制作难度大、加工复杂等问题,为此设计了长方体形结构,其具有以下优点:芯柱制作过程简单;软涂层的工艺更易实现;便于工程应用。当冲击波从一端向另外一端传播时,由于软涂层的存在,芯柱和砂浆可以产生相对运动,起到衰减冲击波的目的[19]

      图  1  本文中提出的超材料混凝土

      Figure 1.  Metaconcrete designed in this paper

    • 图2(a)为超材料混凝土单胞,其中砂浆基体质量为m1,芯柱质量为m2,软涂层等效为刚度为k2的弹簧,砂浆和芯柱的位移分别为u1u2。超材料混凝土可视为图2(b)所示的一维晶格系统。为了计算超材料混凝土的有效质量,又可以把图2(b)假设为如2(c)所示的由刚度为k1(k1→∞)的外部弹簧连接的一维超材料混凝土晶格系统。

      图  2  超材料混凝土单胞及一维晶格系统

      Figure 2.  A unit cell and one-dimensional lattice systems for metaconcrete

      对于图2(a)所示的超材料混凝土单胞,假设涂层为均匀弹性材料,在计算软涂层的等效刚度k2时,只考虑左右两侧软涂层的拉伸和压缩作用,忽略前后两侧软涂层的剪切作用,等效刚度k2可以简化为:

      ${k_2} = \frac{{2Ebl}}{t}$

      式中:E为涂层的弹性模量,b为芯柱的边长,l为芯柱的高度,t为软涂层厚度。

      对于图2(c)所示的一维晶格系统,根据动力学平衡方程,可以给出第j个单胞的运动方程:

      $ m_1^j\ddot u_1^j + {k_1}\left( {2u_1^j - u_1^{j - 1} - u_1^{j + 1}} \right) + {k_2}\left( {u_1^j - u_2^j} \right) = 0 $

      $ m_2^j\ddot u_2^j + {k_2}\left( {u_2^j - u_1^j} \right) = 0 $

      考虑周期系统的稳态简谐振动,可设位移解[20]为:

      $ {u^{j + n}} = A{{\rm{e}}^{{\rm{i}}\left( {jqL + nqL - \omega t} \right)}} $

      式中:A为波幅,q为波数,L为单胞的边长,ω为激振频率。

      从而:

      $ \ddot u_1^j = - {\omega ^2}u_1^j,\;\ddot u_2^j = - {\omega ^2}u_2^j $

      $ u_1^{j + 1} = u_1^j{{\rm{e}}^{{\rm{i}}qL}},\;u_1^{j - 1} = u_1^j{{\rm{e}}^{ - {\rm{i}}qL}} $

      从而式(2)和式(3)可以简化为:

      $ - m_1^j{\omega ^2}u_1^j + {k_1}\left( {2u_1^j - u_1^j{{\rm{e}}^{ - {\rm{i}}qL}} - u_1^j{{\rm{e}}^{{\rm{i}}qL}}} \right) + {k_2}\left( {u_1^j - u_2^j} \right) = 0 $

      $ - m_2^j{\omega ^2}u_2^j + {k_2}\left( {u_2^j - u_1^j} \right) = 0 $

      再次使用下列恒等式:

      $ {{\rm{e}}^{{\rm{i}}qL}} + {{\rm{e}}^{ - {\rm{i}}qL}} = 2\cos \left( {qL} \right) $

      $ 1 - \cos \left( {qL} \right) = 2{\sin ^2}\left( {qL/2} \right) $

      式(7)和式(8)可以简化为:

      $ \left[ { - m_1^j{\omega ^2} + {k_2} + 4{k_1}{{\sin }^2}\left( {\frac{{qL}}{2}} \right)} \right]u_1^j - {k_2}u_2^j = 0 $

      $ - {k_2}u_1^j + \left( {{k_2} - m_2^j{\omega ^2}} \right)u_2^j = 0 $

      为了得到系统的色散关系,可以把式(11)和式(12)看成求解矩阵特征值的问题,可以简化为如下形式:

      $ \left( {{{{K}}_{\rm{r}}} - {\omega ^2}{{{M}}_{\rm{r}}}} \right){{u}} = {\bf{0}} $

      式中:KrMr分别为系统的刚度矩阵和质量矩阵,ω为频率,u为系统中自由度的位移矢量。从而,系统的运动方程可以写成下列形式:

      $ \left( {\begin{array}{*{20}{c}} { - m_1^j{\omega ^2} + {k_2} + 4{k_1}{{\sin }^2}\left( {\dfrac{{qL}}{2}} \right)}&{ - {k_2}} \\ { - {k_2}}&{ - m_2^j{\omega ^2} + {k_2}} \end{array}} \right)\left( {\begin{array}{*{20}{c}} {u_1^j} \\ {u_2^j} \end{array}} \right) = {\bf{0}} $

      当且仅当${{{K}}_{\rm{r}}} - {\omega ^2}{{ M}_{\rm{r}}} = {\bf{0}}$时,此特征值问题具有非零解。从而:

      $ m_1^jm_2^j{\omega ^4} - \left[ {{k_2}\left( {m_1^j + m_2^j} \right) + 4m_2^j{k_1}{{\sin }^2}\left( {\frac{{qL}}{2}} \right)} \right]{\omega ^2} + 4{k_1}{k_2}{\sin ^2}\left( {\frac{{qL}}{2}} \right) = 0 $

      结合式(9)~(10)和(15),可以得到:

      $\cos \left( {qL} \right) = 1 + \frac{{{\omega ^4}{m_1}{m_2} - {\omega ^2}{k_2}\left( {{m_1} + {m_2}} \right)}}{{2{k_1}{k_2} - 2{\omega ^2}{k_1}{m_2}}}$

      对于图2(d)中的一维单原子晶格系统,其等效质量meff[20-21]可以表示为:

      ${m_{e{\rm{ff}}}}{\omega ^2} = 2{k_1}\left[ {1 - \cos \left( {qL} \right)} \right]$

      结合公式(16)和(17),可以得到无量纲的meff表达式:

      $\frac{{{m_{{\rm{eff}}}}}}{{{m_{{\rm{st}}}}}} = 1 + \frac{\theta }{{1 + \theta }}\left[ {\frac{{{{\left( {{\omega / {{\omega _2}}}} \right)}^2}}}{{1 - {{\left( {{\omega / {{\omega _2}}}} \right)}^2}}}} \right]$

      式中:mst=m1+m2θ=m2/m1${\omega _2}{\rm{ = }}\sqrt {{{{k_2}} / {{m_2}}}} $

      根据式(18),当θ=1,可以给出有效质量与激振频率的关系图,如图3所示。

      图  3  有效质量与激振频率的关系

      Figure 3.  Relationship between effective mass and excitation frequency

      图3所示,当激振频率接近由弹簧-芯柱组成的内部系统的自振频率时,有效质量曲线变化特别大并且复杂,当ωω2时,达到负值。图中AB两点之间的区域,有效质量达到负值,波在该频率范围内会得到衰减,形成带隙[16-17]。因此,有效质量是激振频率的函数,正是接近内部系统自振频率的行为导致材料整体性能的变化,从而能够形成一种具有消波能力的新颖材料。

    • 结构的带隙是指不能通过该结构的波的频率范围,其特征主要包含带隙起始频率(band gap start frequency)、带隙截止频率(band gap cutoff frequency)和带隙宽度(band gap width)。对比图3,其中A点对应带隙起始频率,其主要由局域共振骨料的平动振动决定;B点对应带隙截止频率,根据式(18),求出该点对应的频率,表明该共振模式下,芯柱与基体在涂层的作用下,形成反相位,达到动态平衡。

      采用图4所示模型来进一步描述以上两种共振模式。

      图  4  带隙起始频率及截止频率的简化模型[22]

      Figure 4.  Simplified models for band gap start and cutoff frequencies[22]

      图4中,M1为单个局域共振型骨料中芯柱的等效质量,M2为单个周期内砂浆的等效质量,弹簧k为骨料涂层的等效刚度。在带隙起始频率处,芯柱M1在弹簧k的作用下发生共振;而在带隙截止频率处,芯柱M1和砂浆M2在弹簧k的连接下,以相对振动的方式发生共振,弹簧上黑点所示位置保持相对静止。

      对式(18)进行转换,可以得到带隙起始频率f1和带隙截止频率f2的函数表达式:

      ${f_1} = \frac{1}{{2{\rm{{\text{π}} }}}}\sqrt {\frac{k}{{{M_1}}}} $

      ${f_2} = \frac{1}{{2{\rm{{\text{π}} }}}}\sqrt {\frac{{k\left( {{M_1} + {M_2}} \right)}}{{{M_1}{M_2}}}} $

      由此可知:能够影响局域共振骨料等效质量或等效刚度的材料和结构参数,均对局域共振带隙的起始频率产生影响;同时,能够影响局域共振骨料等效质量、等效刚度或基体等效质量的材料和结构参数,均对局域共振带隙的截止频率产生影响。

    • 为了更加深入探究局域共振型超材料混凝土的带隙特征,本文讨论不同参数对超材料混凝土带隙特征的影响,以加深对此共振结构消波机理的认识,用于指导超材料混凝土的结构设计及材料选择。

    • 由式(19)和(20)可知,涂层弹性模量会对等效刚度产生较大的影响。为研究不同涂层材料对带隙特征的影响,在保证其他参数不变的情况下,仅改变涂层弹性模量。现选取单胞的边长L=20 cm,芯柱的边长b=10 cm,软涂层厚度t=1 cm,芯柱材料为铅(Pb,密度为11.6 Mg/m3),基体材料为砂浆(密度为2 Mg/m3)。涂层材料的相关参数参考文献[23]。

      图5给出了涂层弹性模量对带隙特征的影响。从图中可以看出,随着涂层弹性模量的增大,带隙起始频率、截止频率及宽度均增大。根据前文所述,带隙起始频率对应骨料自振,带隙截止频率对应骨料和基体相对某平衡位置产生共振,涂层起到弹簧的作用,涂层弹性模量越大,等效刚度越大,频率越高。当涂层的弹性模量为0.2 MPa时,其对应的带隙起始频率、截止频率和宽度分别为20.91、36.69和15.78 Hz;当涂层的弹性模量为1.0 MPa时,其对应的带隙起始频率、截止频率和宽度分别为46.75、82.04和35.29 Hz。可见,柔性的涂层能够形成低频带隙,但是带隙宽度较窄,而高弹性模量涂层能够形成较宽的带隙,但带隙起始频率较高。因此,可以根据所需衰减的波的频率进行涂层材料的选择。对于用作防止爆炸等高频载荷的材料,超材料混凝土必须在高频下展现共振行为,应选择高弹性模量的涂层材料。对于用作防止地震等低频荷载的材料,超材料混凝土必须在低频下展现共振行为,应选择低弹性模量的涂层材料。

      图  5  涂层弹性模量对带隙特性的影响

      Figure 5.  Effects of coating elastic modulus on band gap

    • 由式(19)和(20)可知,芯柱密度会对等效质量M1产生较大影响,因此,需要进一步研究芯柱密度对带隙特征的影响。采用3.1节的几何尺寸,涂层材料选择silicone(弹性模量为1 MPa),芯柱材料选择铝、钢、铜和铅,其对应的密度分别为2.7、7.9、8.89和11.6 Mg/m3[22]

      图6给出了带隙特征随芯柱密度变化的曲线。从图6可以看出,随着芯柱密度的增大,带隙起始频率和截止频率均降低且带隙起始频率降低的速度较快,从而带隙宽度增大。当芯柱材料为Al时,带隙起始频率、截止频率和宽度分别为96.91、118.05和21.14 Hz;当芯柱材料为Pb时,带隙起始频率、截止频率和宽度分别为46.75、82.04和35.29 Hz。可见,高密度芯柱材料能够形成低频带隙,且带隙宽度较大,而低密度芯柱材料更容易形成高频带隙,且带隙宽度较低。因此,若需要得到低频、宽带隙等特征,应选择密度较高的芯柱材料。

      图  6  芯柱密度对带隙特征的影响

      Figure 6.  Effects of core column density on band gap

    • 由式(19)和(20)可知,基体密度会对等效质量M2产生影响,因此,需要进一步研究基体密度对带隙特征的影响。采用3.1节的几何尺寸,涂层材料选择silicone(弹性模量为1 MPa),芯柱材料选择Pb(密度为11.6 Mg/m3)。

      图7给出了带隙特征随基体密度变化的曲线。从图7可以看出,随着基体密度的增大,带隙起始频率保持不变,带隙截止频率和宽度持续减小。主要原因是:带隙起始频率对应芯柱的自振,其主要由芯柱等效质量及软涂层等效刚度决定,基体密度的变化不会影响芯柱的等效质量,从而带隙起始频率不会产生变化;带隙截止频率,对应芯柱与基体动态平衡,其主要由芯柱等效质量、涂层等效刚度和基体等效质量决定,基体密度的增大会导致基体等效质量增加,从而使得带隙截止频率持续减小,又因为带隙起始频率保持不变,最终带隙宽度也会持续减小。因此,若需要得到低频、宽带隙等特征,应选择密度较低的基体材料。

      图  7  基体密度对带隙特征的影响

      Figure 7.  Effects of matrix density on band gap

    • 为了解骨料体积占比对带隙特征的影响,对于典型的正方形单元,定义骨料体积占比β=R2/a2。其中,R=b+2tb为芯柱边长,t为涂层厚度,a为周期常数。在讨论骨料体积占比对带隙特征的影响时,需要保持参数γγ为芯柱边长与软涂层厚度的比值,γ=b/t)不变。涂层材料选择silicone(弹性模量为1 MPa),芯柱材料选择Pb(密度为11.6 Mg/m3),基体密度为2 Mg/m3L=2 cm,γ=10。

      图8给出了带隙特征与骨料体积占比(β)的关系。从图8可以看出,带隙起始频率随着β的增大而持续降低,且在β从0.1到0.4时,降低速度较快,之后降低速度减缓;带隙截止频率随β的增大呈现先降低再升高的趋势,其中β从0.1到0.4时,带隙截止频率持续降低,之后,带隙截止频率随β的增大而升高,且在β从0.1到0.4时,带隙起始频率的降低速度大于带隙截止频率的降低速度,最终带隙宽度随着β的增大而持续增大。因此,可以通过增大骨料体积占比的方式来实现扩宽带隙的目的。

      图  8  γ=10时,骨料体积占比对带隙特征的影响

      Figure 8.  Effect of volume fraction of aggregate on band gap when γ=10

    • 为了研究参数γ对带隙特征的影响,需要保证骨料体积占比β不变。涂层材料选择silicone(弹性模量为1 MPa),芯柱材料选择Pb(密度为11.6 Mg/m3),基体密度为2 Mg/m3L=2 cm,骨料体积占比β=0.36。

      图9给出了骨料体积占比为36%时,参数γ对带隙特征的影响。从图9可以看出,带隙起始频率及截止频率均随着γ的增大呈现出先减小再增大的变化趋势,而带隙宽度随着γ的增大而持续增大。其中,带隙起始频率在γ=2左右取得极小值;带隙截止频率在γ=1左右取得极小值。因此,若需要得到宽带隙特征,可以通过增大参数γ来实现。

      图  9  当β=0.36时,γ对带隙特征的影响

      Figure 9.  Effect of γ on band gap when β=0.36

    • 为了说明超材料混凝土对冲击波的衰减作用,使用LS-DYNA建立如图10(a)所示的超材料混凝土有限元模型,其长度为100 mm,宽度为60 mm,厚度为20 mm。整个模型包含60 000个单元和67 771个节点,砂浆、软涂层和芯柱均使用八节点的solid单元。为了消除边界效应的影响,在上下两个侧面设置对称性边界条件;为保证混凝土板不发生位移,在前后两个侧面设置固定边界条件[19]。在模型左侧施加如图11所示冲击波。在保证计算精度的情况下,对模型进行简化。混凝土基体、涂层和芯柱均采用线弹性材料模型[23],3种材料之间采用共节点的方式进行连接。同时,为了对比超材料混凝土的优越性能,建立了如图10(b)所示的同样条件下的普通混凝土有限元模型。

      图  10  有限元模型俯视图

      Figure 10.  The top views of finite element models

      图  11  冲击波压力-时间曲线

      Figure 11.  Pressure-time curve of shock wave

    • 图10所示,对于超材料混凝土和普通混凝土,均在距离冲击端20 mm(截面A-A)、40 mm(B-B截面)、60 mm(截面C-C)和80 mm(截面D-D)处依次选取4个截面,得到各截面的平均应力时程曲线,如图12所示。从图12可以看出,对于超材料混凝土和普通混凝土,各截面的平均应力均随着时间的增加而降低,但超材料混凝土中的平均应力降低较快。以1.0 ms左右为例,超材料混凝土中各截面的平均应力在−0.2~0.2 MPa之间,而普通混凝土中各截面的应力在−0.4~0.4 MPa之间。因此,与普通混凝土相比,超材料混凝土对冲击波具有更好的衰减作用。

      图  12  各截面平均应力时程曲线

      Figure 12.  Average stress-time curve at each section

    • 借鉴超材料的研究思路,在混凝土中引入谐振骨料,设计出具有消波特性的超材料混凝土,并对其带隙特征和冲击波衰减效应进行研究,得到以下结论:

      (1)低弹性模量的涂层能够形成低频带隙,但是带隙宽度较窄,而高弹性模量涂层能够形成较宽的带隙,但带隙起始频率较高。因此,应根据所需衰减的波的频率进行涂层材料的选择。

      (2)高密度芯柱材料能够形成低频带隙,且带隙宽度较大,而低密度芯柱材料更容易形成高频带隙,且带隙宽度较小。因此,若需要得到低频、宽带隙等特征,可以选择密度较高的芯柱材料。

      (3)基体密度的变化对带隙起始频率没有影响,但基体密度的升高使带隙截止频率降低和带隙宽度减小。因此,选择低密度的基体材料可以扩宽带隙。

      (4)通过增大骨料体积占比和芯柱边长与软涂层厚度比可以实现扩宽带隙的目的。

      (5)与普通混凝土相比,超材料混凝土对冲击波具有更好的衰减作用。

参考文献 (23)

目录

    /

    返回文章
    返回