›› 2015, Vol. 36 ›› Issue (11): 3307-3314.doi: 10.16285/j.rsm.2015.11.036

• 数值分析 • 上一篇    下一篇

基于土体细观结构重构技术的渗流场数值模拟

申林方1,王志良1,李邵军2   

  1. 1. 昆明理工大学 建筑工程学院,云南 昆明650500;2. 中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室,湖北 武汉 430071
  • 收稿日期:2014-04-10 出版日期:2015-11-11 发布日期:2018-06-14
  • 通讯作者: 王志良,男,1982年生,博士,主要从事岩土工程渗流方面的研究工作。E-mail: wangzhiliangtj@126.com E-mail:shenlinfang@kmust.edu.cn
  • 作者简介:申林方,女,1982年生,博士,主要从事土体渗流数值模拟方面的研究工作。
  • 基金资助:

    国家自然科学基金资助(No. 51179187);昆明理工大学人才培养基金资助(No. KKSY201306023)。

Numerical simulation for seepage field of soil based on mesoscopic structure reconfiguration technology

SHEN Lin-fang1, WANG Zhi-liang1, LI Shao-jun2   

  1. 1. Faculty of Civil Engineering and Architecture, Kunming University of Science and Technology, Kunming, Yunnan 650500, China; 2. State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan, Hubei 430071, China
  • Received:2014-04-10 Online:2015-11-11 Published:2018-06-14
  • Supported by:

    Project supported by the National Natural Science Foundation of China (Grant No. 51179187) and the Personnel Training Fund of Kunming University of Science and Technology (Grant No. KKSY201306023).

摘要: 真实土体的细观结构由许多个大小不一的土颗粒团组成,传统四参数随机生长法(QSGS)构建的土体结构土颗粒团比较均匀,与实际情况存在较大的差异。为弥补这一缺陷,考虑土体孔隙率及自相关函数的影响,对传统的四参数随机生长法进行改进,实现了更接近于真实土体的细观结构重构。在此基础上基于格子Boltzmann方法,采用D2Q9模型,通过设置模型入口、出口边界为非平衡态外推格式,左右边界及土颗粒边界为标准反弹格式的边界条件,建立了模拟重构土体细观渗流场的二维模型。同时,针对一算例编制了相应的计算程序,研究了恒定流速入渗情况下重构土体的细观渗流场。研究表明:土体的渗流方向优先选择连通性较好孔隙所形成的通道,流速受控于通道整体连通性的优劣。整体贯通型的通道流动速度较快,部分连通的孔隙中其流动速度相对较慢。即使局部孔隙空间较大,其渗流速度仍取决于是否位于贯通型通道上。

关键词: 四参数随机生长法, 土体细观结构, 格子Boltzmann方法, 渗流场, 数值模拟

Abstract: The mesoscopic structure of real soil consists of numerous particles with different sizes, and yet the traditional reconstruction method of quartet structure generation set (QSGS) produces relatively uniform soil particles,which has significant discrepancy with the actual mesoscopic structure of soil. In order to remedy this deficiency, the QSGS is improved by considering the influences of porosity and autocorrelative function of soil, so that the mesoscopic structure reconfiguration closer to real soil is obtained. Based on the reconfiguration model combined with Lattice Boltzmann method , the D2Q9 model is applied to construct a two dimensional model for simulating mesoscopic seepage field of reconfiguration soil by setting the non-equilibrium extrapolation scheme at the inlet and outlet boundaries, and the bounce-back scheme at the soil particles’ boundary as well as the left and right boundaries. Meanwhile, according to a case study, a corresponding program is developed to simulate mesoscopic seepage field of reconstructed soil with constant inlet seepage velocity. The results show that the pore fluid preferentially flows through the channel with good connectivity. Meanwhile, the flow velocity is controlled by the whole connectivity of the channel. The velocity is much faster in whole-through-channel than that in partial-through-channel. Even if there exists a large pore space, the velocity still depends on whether or not it is on the channel with integral connectivity.

Key words: quartet structure generation set (QSGS), mesoscopic structure of soil, lattice Boltzmann method(LBM), seepage field, numerical simulation

中图分类号: 

  • O 357.3
[1] 李翻翻, 陈卫忠, 雷江, 于洪丹, 马永尚, . 基于塑性损伤的黏土岩力学特性研究[J]. 岩土力学, 2020, 41(1): 132-140.
[2] 夏 坤, 董林, 蒲小武, 李璐, . 黄土塬地震动响应特征分析[J]. 岩土力学, 2020, 41(1): 295-304.
[3] 郭院成, 李明宇, 张艳伟, . 预应力锚杆复合土钉墙支护体系增量解析方法[J]. 岩土力学, 2019, 40(S1): 253-258.
[4] 闫国强, 殷跃平, 黄波林, 张枝华, 代贞伟, . 三峡库区巫山金鸡岭滑坡成因机制与变形特征[J]. 岩土力学, 2019, 40(S1): 329-340.
[5] 刘红岩. 宏细观缺陷对岩体力学特性及边坡稳定影响研究[J]. 岩土力学, 2019, 40(S1): 431-439.
[6] 金爱兵, 刘佳伟, 赵怡晴, 王本鑫, 孙浩, 魏余栋, . 卸荷条件下花岗岩力学特性分析[J]. 岩土力学, 2019, 40(S1): 459-467.
[7] 韩征, 粟滨, 李艳鸽, 王伟, 王卫东, 黄健陵, 陈光齐, . 基于HBP本构模型的泥石流动力过程SPH数值模拟[J]. 岩土力学, 2019, 40(S1): 477-485.
[8] 吴锦亮, 何吉, . 岩质边坡动态开挖模拟的复合单元模型[J]. 岩土力学, 2019, 40(S1): 535-540.
[9] 吴凤元, 樊赟赟, 陈剑平, 李军, . 基于不同侵蚀模型的高速崩滑碎屑 流动力过程模拟分析[J]. 岩土力学, 2019, 40(8): 3236-3246.
[10] 孙峰, 薛世峰, 逄铭玉, 唐梅荣, 张翔, 李川, . 基于连续损伤的水平井射孔-近井筒三维破裂模拟[J]. 岩土力学, 2019, 40(8): 3255-3261.
[11] 曹洪, 胡瑶, 骆冠勇. 滤管两端均不在含水层层面的承压不 完整井近似计算方法研究[J]. 岩土力学, 2019, 40(7): 2774-2780.
[12] 穆锐, 浦少云, 黄质宏, 李永辉, 郑培鑫, 刘 旸, 刘 泽, 郑红超, . 土岩组合岩体中抗拔桩极限承载力的确定[J]. 岩土力学, 2019, 40(7): 2825-2837.
[13] 金俊超, 佘成学, 尚朋阳. 基于应变软化指标的岩石非线性蠕变模型[J]. 岩土力学, 2019, 40(6): 2239-2246.
[14] 张 聪, 梁经纬, 阳军生, 曹 磊, 谢亦朋, 张贵金, . 堤坝脉动注浆浆液扩散机制及应用研究[J]. 岩土力学, 2019, 40(4): 1507-1514.
[15] 严 健, 何 川, 汪 波, 蒙 伟, . 高地温对隧道岩爆发生的影响性研究[J]. 岩土力学, 2019, 40(4): 1543-1550.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 刘清秉,项 伟,张伟锋,崔德山. 离子土壤固化剂改性膨胀土的试验研究[J]. , 2009, 30(8): 2286 -2290 .
[2] 况雨春,伍开松,杨迎新,马德坤. 三牙轮钻头破岩过程计算机仿真模型[J]. , 2009, 30(S1): 235 -238 .
[3] 杜文琪,王 刚. 土工结构地震滑动位移统计分析[J]. , 2011, 32(S1): 520 -0525 .
[4] 许振浩 ,李术才 ,李利平 ,侯建刚 ,隋 斌 ,石少帅. 基于层次分析法的岩溶隧道突水突泥风险评估[J]. , 2011, 32(6): 1757 -1766 .
[5] 温世清 ,刘汉龙 ,陈育民. 浆固碎石桩单桩荷载传递特性研究[J]. , 2011, 32(12): 3637 -3641 .
[6] 李顺群 ,高凌霞 ,柴寿喜. 冻土力学性质影响因素的显著性和交互作用研究[J]. , 2012, 33(4): 1173 -1177 .
[7] 钟 声 ,王川婴 ,吴立新 ,唐新建 ,王清远. 点状不良地质体钻孔雷达响应特征 ——围岩及充填效应正演分析[J]. , 2012, 33(4): 1191 -1195 .
[8] 孟 振,陈锦剑,王建华,尹振宇. 砂土中螺纹桩承载特性的模型试验研究[J]. , 2012, 33(S1): 141 -145 .
[9] 罗 刚 ,胡卸文 ,顾成壮 . 强震作用下顺层岩质斜坡动力失稳机制及启动速度研究[J]. , 2013, 34(2): 483 -490 .
[10] 付晓东,盛 谦,张勇慧. 水电站地下洞室群分步开挖的非连续变形分析[J]. , 2013, 34(2): 568 -574 .