1. 引言
在我国西部山区交通中,不可避免遇到落石等灾害问题。针对小型落石的防护中,棚洞占有极其重要的地位 [1] [2]。棚洞属于明洞的范畴,即棚式明洞。指为抵挡危岩落石、塌方等自然灾害对已有道路的影响,在隧道进出口或线路特殊地段独立修建的以吸收落石为主要功能的棚式建筑物 [3]。与其他措施相比,棚洞明显减少对环境的破坏,照明通风成本减少,长远来看,棚洞更具有优越性 [4]。
自阳友奎 [3] 在2009年发明了由钢拱架与金属网组合形成的隔离防护飞石或落石的柔性棚洞的实用性专利以来,我国的各位专家学者对柔性棚洞不断有新的突破。
在针对柔性棚洞的使用效果的研究中,汪敏 [4] 、刘成清等 [2] 分别利用模型实验、数值模拟等方法研究了普通钢筋混凝土棚洞与柔性棚洞在实践中的优缺点和使用限制。在实际工程运用中,罗章波等 [5] 展示了一个较为完善的研究思路。刘成清等 [6] 将柔性钢棚洞的结构转化为简单的力学模型,研究了柔性钢棚洞的耗能与极限承载力的问题,提出针对柔性棚洞的设计理论。
本文以九路寨某隧道出口的柔性钢棚洞为例,在对落石运动分析的基础上,使用有限元数值模拟软件模拟落石撞击钢棚洞的过程,分析钢棚洞结构各构件能量、轴力的变化,利用柔性钢棚洞设计理论对已有结构进行验证。
2. 工程概况
九路寨生态旅游区位于保康县西南部歇马镇境内。通往本旅游区某公路隧道为双向单车道,宽6米,高5米,隧道出口上方为险峻的陡崖,陡崖高约30米。
2.1. 地形地貌和地质情况简介
本地为构造剥蚀中低山地貌区,地势险峻,山体陡直,多为断岩绝壁,深沟峡谷。山体坡表多为第四系崩积物所覆盖,厚度较薄。基岩大面积裸露,顶部被密集植被覆盖(图1)。
该地区地层由新到老分别为耕表土(Qpd),松散至稍密状,厚约0.5 m;粉质黏土(Q4el+dl),硬塑状,厚约0.8 m;下伏基岩为古生代奥陶系上统五峰组(O3w)泥质灰岩、白云质灰岩等。
该地区岩层产状平缓,倾向19˚,倾角8˚~10˚。发育两条主要节理,节理1倾向146˚,倾角80˚;节理2倾向285˚,倾角76˚ (图2)。根据隧道区周边地质调绘,沟底和路边落石多呈类球状,大小不一,直径0.05~0.36 m。
![](//html.hanspub.org/file/1-2750877x9_hanspub.png)
Figure 1. The situation of road tunnel exit
图1. 公路隧道出口情况
2.2. 柔性钢棚洞简介
该工程拟采用钢拱架加柔性网的形式进行设计。钢拱架一般使用H型钢。柔性网可分为环形网,支撑绳,加密钢丝网三个部分,环形网与支撑绳主要来承受落石的撞击,吸收落石下落的能量。加密钢丝绳用来阻止小型落石的进入。
根据道路的宽度6 m以及隧道的高度5 m,同时考虑棚顶与隧道最高处应有一定的安全空间,一般为1.5 m,设置棚洞的基本断面尺寸。
柔性钢棚洞的钢结构基本采用的是H 型钢桁架 [7],选用的尺寸为H500 mm × 300 mm。在效能和经济两方面的考虑,对防护网中的环形网采用RXI-025 (ROCCO)环形网 [8] [9]。以带支座的支撑绳纵向五跨,横向一跨,最终形成8个2.5 m × 3 m的网格。防护网分为上下两层,上层为环形网与支撑绳,每个环与周边四个环相互连接,圆环直径为0.3 m;下层为加密钢丝网,加密钢丝网直径为0.05 m × 0.05 m。钢棚洞棚顶防护网平面示意,见图3,整体示意图如图4所示。
![](//html.hanspub.org/file/1-2750877x11_hanspub.png)
Figure 3. The plane diagram of roof protection net
图3. 棚顶防护网平面示意图
![](//html.hanspub.org/file/1-2750877x12_hanspub.png)
Figure 4. The overall schematic diagram of flexible steel shed-tunnel model
图4. 柔性钢棚洞模型整体示意图
3. 边坡落石运动分析
3.1. 落石产生原因
落石的形成与其工程地质条件有极其紧密的关系。主要有地形地貌条件,气候、地层岩性,地质构造,风化作用及人类工程活动等 [10] [11]。
3.2. 落石下落的运动状态分析
当崩塌落石在重力作用下,沿山体坡面滑动、滚动或跳动,当与坡面接触后发生反弹,根据坡面接触点的法向弹性系数和切向摩擦系数的不同,落石的弹跳高度亦不同 [12] [13]。此时,接触坡面对落石产生消能作用,致使落石的动能衰减,直至落石停留动能为零 [14]。整个过程遵循能量的转化和守恒定律。
3.3. 落石运动模型计算
3.3.1. 模型参数
导入工程地质剖面图的CAD文件,建立模型。对不同的地层设置不同的法向弹性系数、切向弹性系数(表1)。根据散落在道路中的落石及边坡节理调查,将道路边坡单块落石方量取0.2 m³,按平均密度2500 kg/m³,设置落石质量为500 kg。只在自重作用下沿着斜坡面滑动的状况进行50次模拟。落石模拟轨迹中,边坡上的落石基本全部掉落在公路内。
![](Images/Table_Tmp.jpg)
Table 1. The setting parameters of slope formation
表1. 边坡地层设定参数
3.3.2. 计算结果及落石冲击力分析
在落石的能量方面,在19 m到25 m处能量较高(图5,图6),最高达到了135 kJ,公路在17 m到23 m的范围,正是在公路靠近陡崖的一侧,设计的棚洞须能够承受住这些冲击力。
![](//html.hanspub.org/file/1-2750877x13_hanspub.png)
Figure 5. The simulation of rock fall trajectory
图5. 对落石轨迹的模拟
![](//html.hanspub.org/file/1-2750877x14_hanspub.png)
Figure 6. The kinetic energy variation of falling rock at horizontal position
图6. 水平位置落石动能变化图
4. 落石对柔性钢棚洞冲击数值模拟
本文拟用有限元数值模拟软件模拟落石对柔性网的冲击过程,对已有的结构进行显式非线性动力分析,从而获取模拟后的各构件能量及轴力等结果。
4.1. 落石对柔性钢棚洞冲击受力分析
落石冲击力是棚洞结构设计所依据的主要荷载之一 [15],冲击力代表着落石下落对棚洞造成的作用力。
冲击力计算公式应用较为广泛的主要有公路路基规范和日本道路公团 [16] 等方法。将落石假设为直径0.36 m的球体,密度为2500 kg/m³,质量为500 kg。路基土层为黏土层,填土重度为18.7 k N/m³,内摩擦角为25°,拉梅常数取1000 kN/m²,泊松比为0.3。使用公路路基设计规范求得冲击力为184.01 kN,使用日本道路公团的计算冲击力为P = 161.98 kN。取两者的较大值即184.01 kN。
4.2. 模型建立
Link160作柔性网中环形网、支撑绳的模拟单元 [17]。采用Beam161梁单元模拟柔性棚洞的拱形钢、挂网支架及横撑。Sold164单元可模拟落石。根据已有的构件的各项基本参数(表2),设置相应的单元参数(表3)。在柔性钢棚洞的模型中,采用非弹性非线性的Plastic Kinematic (塑形随动)材料模型 [18]。落石一般模拟成球形,同体积球型落石与片状、立方体状等落石相比,在位移、冲击力、能量方面破碎性最大 [19]。为简化模型将环形柔性网概化为同面积的正方形进行模拟(图7)。
![](//html.hanspub.org/file/1-2750877x15_hanspub.png)
Figure 7. The schematic diagram of simulation mode
图7. 模拟模型示意图
![](Images/Table_Tmp.jpg)
Table 3. The model element set table
表3. 模型单元设置表
4.3. 数值模拟结果分析
从各构件能量、变形、轴力等方面进行分析。
4.3.1. 能量分析
在落石下落中,由于重力作用,与柔性网的摩擦撞击,会引起总能量变大变小的现象。
能量的途径是落石的动能转化为柔性网的内能。由图8、图9可知,落石最高总能量为120 kJ。
![](//html.hanspub.org/file/1-2750877x16_hanspub.png)
Figure 8. The time-history curve of component total energy
图8. 构件总能量时程曲线
![](//html.hanspub.org/file/1-2750877x17_hanspub.png)
Figure 9. The time-history curve of each component energy
图9. 各构件能量时程曲线
4.3.2. 构件变形分析
选取落石竖向位移最大时刻,绘制各构件的位移云图(图10~12)。
位移数值向下为负,对于环形网的最大位移为1.22 m,安全距离取1.5 m是合适的。从中心到四周,位移量逐渐减少。在支撑绳中的最大位移发生的位置与环形网基本一致,都为落石落入接触的中心位置附近。最大位移为1.18 m。
钢拱架发生最大位移是由于落石下落冲击引起,因落石对柔性网的挤压使得连接柔性网的两个横撑发生相对移动。最大位移的位置为两个钢拱架之间的横杆中,受到挤压力发生位移。最大位移为0.51 m。
4.3.3. 构件轴力分析
选取环形网、支撑绳及钢拱架中轴力较大的几个单元,以单元为研究对象,分析轴力与时间的变化关系。
![](//html.hanspub.org/file/1-2750877x18_hanspub.png)
Figure 10. The diagram of maximum displacement of ring network
图10. 环形网最大位移图
![](//html.hanspub.org/file/1-2750877x19_hanspub.png)
Figure 11. The diagram of maximum displacement of support rope
图11. 支撑绳最大位移图
![](//html.hanspub.org/file/1-2750877x20_hanspub.png)
Figure 12. The diagram of maximum displacement of steel arch
图12. 钢拱架最大位移图
环形网与支撑绳只受到压力的作用,环形网所受轴力较小但起伏较大,最大超过了15 kN将近20 kN (图13),落石下落直接接触点为支撑绳,因此支撑绳最大超过300 kN,达到了350 kN,起到了主要承受拉力的作用(图14)。
钢拱架既受压力也受拉力的作用(图15),轴力为正时为拉力,负时为压力。钢拱架的拉力最大不超过120 kN,但压力最大不超过220 kN,因此对钢拱架来说主要受到压力的作用。压力的位置主要在两个钢拱架之间的横杆中。拉力主要在于挂网支架之间的横撑之间。
![](//html.hanspub.org/file/1-2750877x21_hanspub.png)
Figure 13. The time-history curve of axial force of ring network element
图13. 环形网单元轴力时程图
![](//html.hanspub.org/file/1-2750877x22_hanspub.png)
Figure 14. The time-history curve of axial force of support rope
图14. 支撑绳轴力时程图
![](//html.hanspub.org/file/1-2750877x23_hanspub.png)
Figure 15. The time-history curve of axial force of steel arch shelf
图15. 钢拱架轴力时程图
4.4. 数值模拟对结构尺寸修正
在能量分析中,根据软件模拟的落石最大能量为120 kJ,这数据与数值模拟的135 kJ较为接近,取其最大值135 kJ。这也从侧面证明了数值模拟的一个相对准确程度。
根据最大轴力可以修正最小截面。最小截面面积如表4所示。
柔性网根据截面尺寸,选用ROCCO R5/3/300型即RXI-025防护网,钢拱架由于经济因素,可以选用尺寸为500 mm × 200 mm型号的H型钢。
![](Images/Table_Tmp.jpg)
Table 4. The minimum section area of each member
表4. 各构件最小截面面积
5. 基于能量法对柔性钢棚洞的验证
5.1. 理论假设
利用西南交大刘成清等人提出的能量法理论 [6] 对柔性钢棚洞进行验证,该设计理论从钢拱架以及柔性钢绳索开始,针对两者的屈服机制,分别求得在两种情况下的钢拱架、刚绳索以及防护网最终的吸收能量值,加在一起即为整个钢棚洞的耗能值。理论假设:1) 落石下落位置取最不利的情况,一般为两个钢拱架中间的刚绳索中。2) 落石撞击的力,由两条支撑绳承受。
5.2. 理论验算
在设计中,柔性网是由环形网及支撑绳共同承担,但在能量法中,只有两个刚绳索承担作用力,因此需要求得环形网与支撑绳等效为两根刚绳索的半径。
通过对比钢拱架顶部与根部屈服分析可知钢拱架的根部先于顶部出现塑性铰。钢拱架结构极限荷载计算中,在钢结构根部先出现塑性铰的情况下,会经历弹性阶段、2个塑性铰阶段、3个塑性铰阶段、最终达到极限状态。通过验算,最终极限状态时极限承载力为448.92 kN。而刚绳索的极限承载力为143.80 kN。刚绳索比钢拱架先屈服。最终的总耗能为钢拱架、刚绳索以及柔性网耗能的总和,为206.54 kJ。
6. 结语
本文以九路寨公路隧道出口落石防治设计为例,通过对落石的研究,研究落石轨迹、冲击能量以及冲击力。使用有限元数值模拟软件模拟落石对棚洞的碰撞作用,得到构建内部能量轴力等数据,最后基于能量法对已有钢棚洞进行检验。
在本次设计中得到的结论如下:1) 在柔性钢棚洞对落石的冲击中,钢拱架受到的最大压力为拉力的2倍左右,因此要注意钢拱架的抗压性能;2) 在对落石能量吸收中,柔性网中支撑绳吸收了大部分的能量,因此要重视对支撑绳的选取;3) 在利用能量法理论计算时,由于能量法理论模型为两条纵线与实际结构的环形网有较大出入,能量法偏于保守。