摘要
在滚石冲击破坏建(构)筑物灾害中,滚石形状往往具有不确定性,建(构)筑物表面形状也不尽相同。当滚石冲击结构物时,由于滚石形状的不同导致其所表现出来的冲击动力响应特征也不同。采用ABAQUS有限元软件中的Explicit动力有限元方法,分别开展不同形状滚石冲击不同表面形状结构物的数值模拟研究,并对不同缓冲层的缓冲耗能特性进行对比分析。结果表明:在质量、初速度相同的初始条件下,形状不同的滚石冲击平面墙状的结构物时,冲击力峰值大小依次为正方体>圆柱体(底面)>圆柱体(侧面)>球体>三棱锥;当滚石冲击面形状由“尖”变“钝”时,结构物受滚石冲击处被破坏程度将变小,但其整体变形区域增大;对不同表面形状的结构物承受滚石冲击时,冲击力峰值大小依次为平面状>凸面状>凹面状>波浪面状>斜面状,相比而言,斜面状结构物具有较好的抗滚石冲击性能;对在结构物前方设置不同缓冲层时对滚石冲击的缓冲耗能特性进行对比分析,得出EPS-砂土复合垫层具有较好的缓冲耗散性能,其在实际的滚石灾害防护工程中具有广阔的应用前景。
我国疆域辽阔,地形多样,其中山地面积约占总面积的2/3,众多的山区地形孕育了大量的地质灾害,如滚石、崩塌、滑坡、泥石流等。尤其在我国的西南地区,自2008年5月12日汶川地震后加剧了山体破碎,孤石遍布,诱发了大规模的滚石地质灾害,导致部分房屋、车辆被损毁,道路、桥梁等基础设施被破坏,交通干道被迫中断,给人民生命财产、基础设施运营等造成了严重的损失和安全隐

图1 典型滚石冲切破坏建(构)筑物地质灾
Fig.1 Typical geohazards of rockfalls impact where
buildings or structures were destroyed
当前,关于滚石地质灾害的主要研究方法包括理论研究、物理模型试验、现场原位试验及数值模拟研究。其中,理论研究虽能描述滚石运动前后速度及能量的变化,但是难以预测滚石运动轨迹;物理模型试验因存在缩尺效应,不能完全重现滚石灾害过程,因而与实际相比存在较大误差;现场原位试验由于成本高、周期长,因此很少有人开展现场原位试验研究;数值模拟成本低、时间短,可模拟复杂工况,作为试验的补充手段,备受科研工作者的青睐。
何思明
滚石防护主要可分为主动防护与被动防护,但由于山区地形条件限制,主动防护往往难以实施,故工程上常采用棚洞、拦石墙和帘式网等被动防护措施来减轻或避免滚石灾害对灾害易发区沿线公路、铁路等构筑物造成损害。许多学者针对不同的防护结构开展了试验研究,Kishi N.
现有研究大多是对滚石冲击结构物机理的探究,但对于滚石形状和受冲击结构物形状对滚石冲击过程的影响方面的研究相对较少。本文将滚石和结构物表面分别设置为4种不同的形状,采用数值模拟研究方法探究滚石冲击结构物过程,并对比分析不同缓冲层对滚石冲击的耗能效果。
滚石是指各种石块因某种原因从地质体表面失稳脱落后经下落、回弹、跳跃、滚动或滑动等运动方式中的一种或几种的组合沿着坡面向下快速运动,最后在较平缓的地带或障碍物附近停止下来的一个动力演化过

图2 滚石运动示意
Fig.2 Schematic diagram of rockfall movement
由
本文采用ABAQUS有限元软件中的Explicit动力有限元方法进行数值模拟计算。由于滚石形状具有不确定性,不同形状的滚石对结构物的冲击响应特征差异较大。因此,本文研究的对象是不同形状的滚石对平面墙状构筑物的冲击动力响应特征,所以建立了如

图3 4种不同形状滚石模型的结构示意
Fig.3 Modeled structures of four different
shaped rockfalls
本文的另一研究对象是同一形状滚石对具有不同表面形状结构物的冲击响应特征,设置了如

图4 4种不同表面形状结构物示意
Fig.4 Schematic diagram of four structures
with different surface shapes
当不同形状的滚石冲击不同表面的结构物时,它们之间相互接触的工况条件也不同。如

图5 不同滚石冲击构筑物原理
Fig.5 Schematic diagram of impact of structures by different rockfalls
此外,为探究不同缓冲层对滚石冲击的缓冲耗能效果,建立了如

图6 滚石冲击结构物的数值计算模型
Fig.6 Numerical calculation model of the structures
impacted by rockfalls

图7 不同形状滚石的冲击力时程曲线
Fig.7 Time history curve of impact force of the rockfalls with various shapes
以结构物中心受滚石冲击后的单元位移变形情况来表征其承受滚石冲击点的破坏程度。

图8 结构物承受不同滚石冲击的单元位移变形时程曲线
Fig.8 Time history curve of the unit displacement of the structures impacted by various rockfalls
为了更加清晰地观察、对比分析不同形状滚石冲击结构物中部位置时的动力响应特征,将各模型进行中心剖分,相应的模型位移云图如

图9 不同形状滚石冲击结构物时的中心剖分模型位移云图
Fig.9 Displacement contours of the central subdivision models when rockfalls with different shapes impact structures
建立的滚石冲击不同表面形状结构物的数值计算模型如

图10 不同表面形状结构物工况下的滚石冲击力时程曲线
Fig.10 Time history curve of rockfall impact force on the structures with various shapes

图11 不同表面结构物受滚石冲击时单元位移变形时程曲线
Fig.11 Time history curve of the unit displacement of the structures with various shapes impacted by rockfalls

图12 正方体滚石冲击不同表面形状结构物时的中心剖分模型位移云图
Fig.12 Displacement contours of the central division models when the cube rockfall impacts
the structures with different surface shapes
综上所述,对于表面形状为凹面、凸面、波浪面和斜面状的结构物,斜面状的结构物具有更好的抵抗滚石冲击性能,虽然当其遭受滚石冲击后的位移变形量(约为0.47 mm)略大于平面状的结构物,但在上述4种表面形状的结构物中,斜面状具有冲击力小、结构简单、便于施工、节省材料及遭受滚石冲击时形成的整体变形区域很小等众多优点。
为减轻滚石冲击对结构物造成严重损害,现有研究多是在被保护结构(或滚石被动防护结构)前设置各种材质的缓冲层以降低滚石冲击力对结构物的损伤。将目前常用的EPS泡沫垫层、砂土垫层、泡沫铝垫层及EPS-砂土复合垫层4种缓冲材料的缓冲性能进行比较分析,对于EPS-砂土复合垫层, EPS置于结构物和砂土之间。采用的数值计算模型如

图13 不同缓冲层下的滚石冲击力时程曲线
Fig.13 Time history curve of rockfall impact force for various buffer layers
因此,泡沫铝和砂土缓冲层对滚石冲击力的缓冲效果比另外2种缓冲层差,EPS和EPS-砂土复合缓冲层的缓冲效果较优。

图14 设置不同缓冲层后滚石冲击结构物时的中心剖分模型位移云图
Fig.14 Displacement contours of the central subdivision models when rockfalls impact structures
with different buffer layers
通过运用ABAQUS有限元软件中的Explicit动力有限元方法对滚石冲击结构物动力响应特征进行数值模拟研究,就本文研究假设的4种形状和几何尺寸的滚石建立的单一几何体模型而言,可得出以下主要结论:
(1)不同形状滚石对平面墙状结构物的冲击力峰值大小依次是正方体>圆柱体(底面)>圆柱体(侧面)>球体>三棱锥;滚石冲击面形状越尖锐,结构物受滚石冲击处的破坏程度越严重,但结构物内部整体变形区域越小;当滚石冲击面形状由“尖”变“钝”时,结构物受滚石冲击处的破坏程度变小,但整体变形区域变大。
(2)正方体滚石冲击不同表面形状结构物时,各形状结构物所受滚石冲击力峰值大小依次是平面>凸面状>凹面状>波浪面状>斜面状;其受冲击区域破坏程度大小依次是波浪面>斜面>凹面>凸面>平面;凹面状和凸面状结构物在受滚石冲击时的整体变形区域较大,而波浪面状和斜面状结构物形成的整体变形区域较小;相比而言,斜面状结构物具有较好的抗滚石冲击性能。
(3)对不同缓冲层的滚石冲击缓冲耗散性能进行对比分析,得出EPS-砂土复合缓冲层和EPS缓冲层的缓冲性能较优,但在工程实际应用中EPS缓冲层易被滚石冲切破坏,还会造成环境污染,因而EPS-砂土缓冲层为一种工程实际中更有广阔应用前景的缓冲垫层型式。
(4)本文提出的滚石模型均为规则样式,得到的相关结果和结论可为实际的滚石灾害防护提供一定的指导依据,但在实际的滚石地质灾害中,其形状较为复杂多样且不规则,使得其运动和冲击特性更为复杂,在后续研究中应进一步加强对非规则滚石冲击结构物动力响应特征的探究。
参考文献(References)
张路青,杨志法,许兵.滚石与滚石灾害[J].工程地质学报,2004(3):225-231. [百度学术]
ZHANG Luqing, YANG Zhifa, XU Bing. Rock falls and rock fall hazards[J]. Journal of Engineering Geology, 2004(3):225-231. [百度学术]
刘洋.滚石冲击棚洞防护结构动力响应及作用机理研究[D].成都:成都理工大学,2017. [百度学术]
LIU Yang. Dynamic response and interaction mechanism of the rock shed subject to rockfall impact[D]. Chengdu: Chengdu University of Technology, 2017. [百度学术]
何思明,吴永,杨雪莲.滚石坡面冲击回弹规律研究[J].岩石力学与工程学报,2008,27(S1):2793-2798. [百度学术]
HE Siming, WU Yong, YANG Xuelian. Study of rock‑fall motion on slope[J]. Chinese Journal of Rock Mechanics and Engineering, 2008,27(S1):2793-2798. [百度学术]
叶四桥,巩尚卿.落石碰撞法向恢复系数的模型试验研究[J]. 中国铁道科学,2015,36(4):13-19. [百度学术]
YE Siqiao, GONG Shangqing. Research on normal restitution coefficient of rockfall collision by model tests[J]. China Railway Science, 2015,36(4):13-19. [百度学术]
黄润秋,刘卫华.基于正交设计的滚石运动特征现场试验研究[J].岩石力学与工程学报,2009,28(5):882-891. [百度学术]
HUANG Runqiu, LIU Weihua. In‑situ test study of characteristics of rolling rock blocks based on orthogonal design[J]. Chinese Journal of Rock Mechanics and Engineering, 2009,28(5):882-891. [百度学术]
裴向军,刘洋,王东坡.滚石冲击棚洞砂土垫层耗能缓冲机理研究[J].四川大学学报(工程科学版),2016,48(1):15-22. [百度学术]
PEI Xiangjun, LIU Yang, WANG Dongpo. Study on the energy [百度学术]
dissipation of sandy soil cushions on the rock‑shed under rockfall impact load[J]. Journal of Sichuan University (Engineering Science Edition), 2016,48(1):15-22. [百度学术]
王东坡,何思明,吴永,等.滚石防护棚洞EPS垫层结构缓冲作用研究[J].振动与冲击,2014,33(4):199-203. [百度学术]
WANG Dongpo, HE Siming, WU Yong, et al. Cushioning effect of rock sheds with EPS cushion on rock‑falls action[J]. Journal of Vibration and Shock, 2014,33(4):199-203. [百度学术]
王东坡,何思明,李新坡,等.冲击荷载下EPS垫层棚洞耗能减震作用研究[J].四川大学学报(工程科学版),2012,44(6):102-107. [百度学术]
WANG Dongpo, HE Siming, LI Xinpo, et al. Study on the dissipating effects of shed with EPS cushion under impact load[J]. Journal of Sichuan University (Engineering Science Edition), 2012,44(6):102-107. [百度学术]
王静峰,赵鹏,袁松,等.复合垫层钢棚洞抵抗落石冲击性能研究[J].土木工程学报,2018,51(S2):7-13. [百度学术]
WANG Jingfeng, ZHAO Peng, YUAN Song, et al. Numerical study on impact resistance of steel shed gallery with composite cushion[J]. Journal of Civil Engineering, 2018,51(S2):7-13. [百度学术]
Kishi N., Konno H., Ikeda K. Prototype impact tests on ultimate impact resistance of PC rock‑sheds[J]. International Journal of Impact Engineering, 2002,27(9):969-985. [百度学术]
Lambert S., Gotteland P., Nicot F. Experimental study of the impact response of geocells as components of rockfall protection embankments[J]. Natural Hazards and Earth System Sciences, 2009,9(2):459-467. [百度学术]
Buzzi O., Spadari M., Giacomini A., et al. Experimental testing of rockfall barriers designed for the low range of impact energy[J]. Rock Mechanics and Rock Engineering, 2013,46(4):701-712. [百度学术]
Peila D., Oggeri C., Castiglia C. Ground reinforced embankments for rockfall protection: Design and evaluation of full scale tests[J]. Landslides, 2007,4(3):255-265. [百度学术]
OUYANG Chaojun, LIU Yang, WANG Dongpo, et al. Dynamic analysis of rockfall impacts on geogrid reinforced soil and EPS absorption cushions[J]. KSCE Journal of Civil Engineering, 2019,23(1):37-45. [百度学术]
张晓强.桥隧相连落石防护新型结构动力响应分析[D].北京:北京交通大学,2018. [百度学术]
ZHANG Xiaoqiang. Dynamic response analysis for a new type of rockfall protection structure with connection of the bridge and the tunnel[D]. Beijing: Beijing Jiaotong University, 2018. [百度学术]
胡卸文,梅雪峰,杨瀛,等.落石冲击荷载作用下的桩板拦石墙结构动力响应[J].工程地质学报,2019,27(1):126-136. [百度学术]
HU Xiewen, MEI Xuefeng, YANG Ying, et al. Dynamic response of pile‑plate rock retaining wall under impact of rockfall[J]. Journal of Engineering Geology, 2019,27(1):126-136. [百度学术]