基于PFC3D的滑坡堰塞坝堆积过程与形态模拟

发布时间：2021-11-19 21:07

　　基于PFC3D离散元软件,开展雅砻江唐古栋滑坡1967年失稳堆积堰塞坝的反演模拟,并基于反演参数对强变形A区失稳堆积进行预测模拟。结果表明:①通过开展唐古栋滑坡1967年失稳堆积堰塞坝的反演分析,校核得到一组合理的细观力学参数;②强变形A区失稳并整体下滑约1 600 m后在河谷堆积形成堰塞坝,堰塞坝形态呈中部高两边低的梯形形态分布,坝顶宽437.91 m,坝底宽994.39 m,坝高92.65 m;③堰塞坝在3种不同来流量下发生溃坝,其洪峰流量均远大于楞古水电站拟选坝址校核洪水流量。通过模拟发现,PFC3D软件对于模拟滑坡堰塞坝堆积过程及堆积形态有较好的适用性,可以获取堰塞坝堆积的三维形态和较准确的坝高。 

【文章来源】：人民长江. 2020,51(04)北大核心

【文章页数】：8 页

【部分图文】：

唐古栋滑坡全貌及分区

基本计算原理是运用显示中心差分法进行动态松弛,计算过程是在每个时步内交替利用力-位移定律与牛顿第二定律,时步迭代并遍历整个颗粒集合体(见图2)。材料的本构行为通过接触点上的本构模型来模拟,颗粒间的接触本构模型有3种:接触刚度模型、滑动模型和接触连接模型。

表3 唐古栋滑坡PFC3D数值模拟细观参数Tab.3 Mesoscopic parameters of Tanggudong landslide for PFC3D simulation 颗粒密度/(kg·m-3) 颗粒最小半径Rmin/m 颗粒半径比/Rmax/Rmin 法向刚度/(N·m-1) 切向刚度/(N·m-1) 摩擦系数 连接半径/m 黏结模量Ec/pa 黏结刚度比(kn/ks) 法向强度/(N·m-2) 切向强度/(N·m-2) 摩擦系数 2020 3 1.67 1×106 1×106 0.324 1.0 2.5×107 1.2 2×107 1.5×107 0.152.3.2 堰塞坝形态反演分析

【参考文献】：
期刊论文
[1]金沙江“10·10”白格堰塞湖溃坝洪水反演分析[J]. 陈祖煜,张强,侯精明,王琳,马利平.  人民长江. 2019(05)
[2]金沙江白格堰塞体结构形态与溃决特征研究[J]. 蔡耀军,栾约生,杨启贵,徐复兴,张胜军,石裕,易杜靓子.  人民长江. 2019(03)
[3]金沙江雅鲁藏布江堰塞湖应急处置回顾与思考[J]. 金兴平.  人民长江. 2019(03)
[4]基于PFC2D的山区中小型水库滑坡堵江模拟研究[J]. 李子隆,王浪.  人民珠江. 2019(03)
[5]金沙江白格堰塞湖溃决洪水预报误差与改进[J]. 王敏,卢金友,姚仕明,朱勇辉,周建银.  人民长江. 2019(03)
[6]基于PFC2D的岸坡覆盖层失稳堵江模拟研究[J]. 李子隆,许模,赵勇,王子忠,郭健,李明月.  水力发电. 2018(01)
[7]唐古栋滑坡成因机制研究[J]. 易志坚,黄润秋,吴海燕,何顺勋,杨建,周龙寿.  工程地质学报. 2016(06)
[8]滑坡型堰塞坝的形成条件与过程分析[J]. 庞林祥,莫大源,李爱华.  人民长江. 2016(11)
[9]堰塞坝及其溃决模拟研究评述[J]. 王光谦,王永强,刘磊,王大宇.  人民黄河. 2015(09)
[10]唐古栋滑坡变形破坏机制及岸坡稳定性研究[J]. 肖华波,王刚,郑汉淮,夏万洪.  长江科学院院报. 2014(11)

博士论文
[1]唐家山高速短程滑坡堵江及溃坝机制研究[D]. 罗刚.西南交通大学 2012

硕士论文
[1]川藏交通廊道冰碛物滑坡堵江风险及线路对策研究[D]. 梁承洋.西南交通大学 2015
[2]楞古水电站唐古栋巨型滑坡成因机制及稳定性研究[D]. 易志坚.成都理工大学 2010
[3]白龙江中游滑坡堵江问题研究[D]. 原俊红.兰州大学 2007



本文编号：3505870