以离散元理论为基础,采用颗粒流软件PFC建立落石冲击棚洞模型,在有回填土和无回填土两种工况下,模拟山体落石冲击框架棚洞,并对框架棚洞不同部位的冲击效应进行研究。研究表明:在落石冲击作用下,棚洞结构同一部位的动力响应不同,不同部位的动力响应也不相同,靠近落石冲击部位附近以剪应力为控制应力,破坏形式是剪切破坏,顶板端部则以正应力为控制应力,破坏形式是挤压破坏;有回填土情况下落石冲击效应小于无回填土,说明回填土对落石冲击具有非常好的缓冲作用;回填土的布置对顶板在落石冲击作用下的受力情况有影响。 Based on the discrete element theory, a particle model was made to simulate the process of rockfall impact to shed cave under with and without backfill soil and the impact effect of different parts of the shed structure was studied and analyzed. The results drew conclusions as below: The stress response of the same part was different and meanwhile, the dynamic response of different parts was neither different under the action of rockfall impact. The shear stress is taken as the control stress near the impact position of the falling rock, and the failure form is shear failure. At the end of the roof, the normal stress is used as the control stress, and the failure form is extrusion failure. The impact effect of backfill is less than that of non-backfill, which indicates that backfill has very good buffering effect on rockfall impact, and the layout of backfill has an effect on the stress distribution of roof under rockfall impact.
李俊杰1,李俊隆2,王玉锁3
1中国港湾工程有限责任公司,北京
2中交机电工程局有限公司,北京
3西南交通大学土木工程学院,四川 成都
收稿日期:2019年5月6日;录用日期:2019年5月21日;发布日期:2019年5月28日
以离散元理论为基础,采用颗粒流软件PFC建立落石冲击棚洞模型,在有回填土和无回填土两种工况下,模拟山体落石冲击框架棚洞,并对框架棚洞不同部位的冲击效应进行研究。研究表明:在落石冲击作用下,棚洞结构同一部位的动力响应不同,不同部位的动力响应也不相同,靠近落石冲击部位附近以剪应力为控制应力,破坏形式是剪切破坏,顶板端部则以正应力为控制应力,破坏形式是挤压破坏;有回填土情况下落石冲击效应小于无回填土,说明回填土对落石冲击具有非常好的缓冲作用;回填土的布置对顶板在落石冲击作用下的受力情况有影响。
关键词 :落石冲击,棚洞,离散单元法,数值模拟
Copyright © 2019 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
落石灾害在我国山区铁路建设、运行过程中经常发生,导致列车受损,轨道损毁,人员伤亡。业界人士针对落石对结构的冲击效应展开系统研究。在研究初期,落石冲击效应的研究多以试验模型为主,如隧道设计技术手册 [
近年来,很多学者对于棚洞结构的落石冲击效应方面的研究有了很大的进展。裴向军,刘洋,王东坡等 [
本文采用ITASCA公司开发的软件PFC来建立落石、框架棚洞、回填土等模型,模拟落石垂直下落在有回填土和无回填土两种情况下,对棚洞结构的冲击过程,通过在棚洞不同部位(如顶板、立柱、纵梁等)布置监测点来获取相应部位的应力–应变时程曲线,从而对框架棚洞的冲击效应展开研究并总结。
本次建模包括物理模型和本构模型。物理模型包括山体、棚洞结构、回填土层、落石等分模型,棚洞结构组成部分主要有顶板、盖梁、支柱、基础等;本构模型主要讲述不同物理模型之间的相互作用,比如落石与回填土、回填土与棚洞结构之间的相互作用。本次建模主要采用接触本构模型,包括颗粒与颗粒之间的粘结和滑移模型以及颗粒与墙体之间的接触刚度模型 [
1) 棚洞结构模型
PFC3D模型中颗粒间通过平行粘结方式粘结,平行粘结模式有五个参数,其中平行粘结法向刚度pb_kn和平行粘结切向刚度pb_ks,单位是Pa/m,代表颗粒间经过平行粘结后的法向和切向刚度,影响粘结模型的变形特性,而平行粘结法向强度pb_ns和平行粘结切向强度pb_ss单位是Pa,代表的是颗粒间法向和切向粘结强度,影响粘结模型的强度特性;平行粘结半径放大系数pb_rad,代表的是平行粘结半径放大系数,粘结半径等于最小颗粒半径乘以这个半径放大系数。
棚洞结构由顶板、盖梁、支柱、基础组成,全部用球体颗粒单元通过平行粘结来模拟。平行粘结主要有五个微观参数,分别是pb_kn、pb_ks、pb_ns、pb_ss和pb_rad,使得模拟材料的参数属性与钢筋混凝土接近。
2) 回填土、山体和落石模型
回填土以及落石都用球体颗粒来模拟,山体坡面采用墙体单元模拟。
本次离散元数值模拟采用接触刚度模型,回填土颗粒单元采用接触粘结本构模型,棚洞颗粒单元采用平行粘结结构模型。
采用离散元颗粒流数值模拟软件PFC建模,模拟工况是有回填土(t = 1 m)和无回填土(t = 0)两种情况,落石重量2.0 t,垂直下落冲击棚洞结构。落石冲击时,需要考虑落石重力、回填土的摩擦力、回填土对落石的粘滞阻尼力和落石与回填土以及棚洞结构的碰撞力,阻尼选用粘滞阻尼,粘滞阻尼系数取值0.2。
落石冲击棚洞结构的物理模型和数值建模1) 模型概述
棚洞结构设计图如图1所示,颗粒流PFC3D建模如图2所示。
图1. 棚洞结构设计图 [
图2. 落石、棚洞PFC3D模型
2) 离散元模型单元微观参数
模型微观参数 [
项目 | 单元类型 | 密度 | 摩擦系数 | 颗粒法向刚度 | 颗粒切向刚度 |
---|---|---|---|---|---|
Dens (kg/m3) | fric | Kn (N/m) | Ks (N/m) | ||
落石 | ball | 2600 | 0.5 | 1e8 | 1e8 |
回填土 | ball | 2000 | 0.5 | 1e8 | 1e8 |
棚洞 | ball | 2500 | 0.5 | 1e15 | 1e15 |
表1. 模型微观参数
3) 回填土参数标定
根据土力学 [
弹性模量E0 (MPa) | 压缩模量Es (MPa) | 泊松比υ | 粘聚力c( kPa) | 内摩擦角φ (˚) |
---|---|---|---|---|
22.0 | 38.1 | 0.366 | 40 | 21 |
表2. 三轴试验PFC3D数值模拟宏观参数结果 [
4) 工况组合
本次数值模拟工况组合表3。
变量 | H (m) | 坡度(˚) | 重量(t) | 回填土厚(m) |
---|---|---|---|---|
取值 | 5 | 90 | 2.0 | 0 |
5 | 90 | 2.0 | 1 |
表3. 工况组合 [
监测点布置情况见图3(a)~(b),单元三维空间应力符号说明见图4,图中所示应力均是应力正方向,其中三个正应力σxx、σyy、σzz正值表示拉应力,负值表示压应力 [
图3. (a) 顶板监测点布置图;(b) 纵梁、支柱和基础监测点布置图
图4. 单元三维空间应力符号 [
用三种方法研究落石冲击效应,三种方法介绍如下:
1) 单个监测点位应力时程曲线
在离散元数值模拟软件PFC中用fish语言可以实现用测量小球(measure sphere)读取监测点部位(例如选取监测点18冲击部位顶板中心)的六个应力分量,在模拟计算的过程中自动记录应力分量,绘制应力时程曲线,进而可以对该监测点的最大应力进行分析。
2) 条带峰值应力曲线
在框架棚洞结构选取顶板中心横向条带(监测点3,8,13,18,23,28,33)和纵梁纵向条带(36,37,38,39,40),在PFC中通过Fish语言记录、读取并导出条带上各监测点的应力值,从中找出各监测点各应力分量的最大值(正值/负值)。
3) 应力峰值云图
在PFC软件中用fish语言记录、读取、导出所有监测点的应力数值,提取出各应力分量的最大值,在origin软件中导入相应数据绘制框架棚洞的峰值应力云图,进而在应力云图上找到在落石冲击下最大应力发生的部位,即危险部位。
选取框架棚洞顶板横向条带几个监测点(3,8,13,18,23,18,33)进行分析和研究。框架棚洞顶板结构横向条带各监测点应力–时程曲线如图5所示。
图5. 框架棚洞顶板横向条带各监测点应力–时程动态变化曲线
从图5中可以得出以下结论:
1) 棚洞结构在落石冲击下,同一部位的各个应力分量变化情况不一致,例如图5(a)中显示,监测点3在落石冲击作用下Y向正应力突变最大,其他应力分量变化较小,说明该部位落石冲击效应以Y向正应力σyy为控制应力,同时为负值,如发生破坏就是受压破坏;
2) 棚洞结构在落石冲击下,不同部位的冲击效应不一致,如图5(a)、图5(c)、图5(k)、图5(m)中显示,顶板靠近端部位置(监测点3、8、28、33)都是Y向正应力σyy突变最大,而图5(e)、图5(g)、图5(i)中显示,顶板靠近中间(监测点13、18、23)落石冲击部位附近,则是τxy表现增大比较明显,说明顶板靠近冲击部位范围内,切应力τxy是控制应力,如果发生破坏则是剪切破坏,即顶板被落石砸穿或者砸破。
3) 回填土对落石冲击具有很大的缓冲作用。分别对比回填土厚土t = 0和t = 1 m各监测点的应力曲线,如图5(a)和图5(b),有回填土时的最大应力明显小于无回填土,可以说明回填土对落石冲击效应有很大的缓冲作用。
选取框架棚洞结构顶板结构横向条带几个监测点(3,8,13,18,23,18,33)进行分析和研究。框架棚洞顶板结构横向条带各监测点最大应力曲线见图6所示。
图6. 顶板横向条带各监测点最大应力分布曲线
从图6中可以得出以下结论:
1) 在落石冲击下,棚洞结构不同部位的动力响应不一致。如图6(a)、图6(c)、图6(e)中显示,正应力分量突变最大的是Y向正应力σyy突变最大,最大值达到σyy = 120 MPa,发生在顶板端部位置;而图6(g)、图6(i)、图6(k)中显示;切应力突变最大的则是XY方向切应力τxy,发生在顶板中间靠近落石冲击部位附近,最大值达到τxy = 25 MPa,破坏形式将是剪切破坏,即顶板被落石砸穿或者砸破。该结论与4.1章的结论吻合。
2) 回填土对落石冲击具有很大的缓冲作用。分别对比回填土厚土t = 0和t = 1 m各部位峰值应力分布情况,如图6(a)和图6(b),有回填土时最大应力明显小于无回填土,说明回填土对落石冲击具有很好的缓冲作用。该结论与4.1章的结论吻合。
选取框架棚洞顶板结构(监测点1-35)进行分析,顶板的应力峰值云图见图7(a)~(x)。
从PFC3D软件中导出监测点的应力,再将数据导入Origin软件中进行处理,绘制应力云图,目前这种方法还在探索阶段。如何能更高效、准确地绘制离散元颗粒流的应力云图,是今后应该努力探索一个方向。
图7. 顶板峰值应力云图
本文使用PFC3D软件建立落石冲击棚洞模型,模拟落石冲击过程,用Fish语言记录冲击过程中不同部位的应力分量,分别绘制了单个监测点位的应力时程曲线、条带应力峰值分布曲线、顶板结构峰值应力云图,在落石冲击下,研究棚洞结构的动态应力响应。可以得出以下结论:
1) 棚洞结构在落石冲击下,同一部位的正应力和切应力变化趋势不一致;不同部位的动态响应不同;具体分析见本文4.1、4.2。
2) 框架棚洞顶板在落石冲击作用下,落石冲击部位控制应力是XY方向剪切应力τxy,破坏形式是剪切破坏;顶板端部控制应力是Y向正应力σyy,破坏形式是受压破坏;具体分析见本文4.1、4.2。
3) 回填土对落石冲击棚洞结构具有非常明显的缓冲作用,具体分析见本文4.1、4.2;填筑回填土会影响棚洞顶板在落石冲击作用下的应力分布规律,见本文图7。
李俊杰,李俊隆,王玉锁. 落石冲击棚洞结构的PFC3D数值模拟研究 PFC3D Simulation Study of Rockfall Impact to Shed Tunnel Structure[J]. 土木工程, 2019, 08(03): 800-810. https://doi.org/10.12677/HJCE.2019.83093