• ISSN 0258-2724
  • CN 51-1277/U
  • EI Compendex
  • Scopus 收录
  • 全国中文核心期刊
  • 中国科技论文统计源期刊
  • 中国科学引文数据库来源期刊

注浆渗透扩散的多物理场耦合数值分析

陈锋 杨杰 张冲 余祯 刘先峰

陈锋, 杨杰, 张冲, 余祯, 刘先峰. 注浆渗透扩散的多物理场耦合数值分析[J]. 西南交通大学学报. doi: 10.3969/j.issn.0258-2724.20220763
引用本文: 陈锋, 杨杰, 张冲, 余祯, 刘先峰. 注浆渗透扩散的多物理场耦合数值分析[J]. 西南交通大学学报. doi: 10.3969/j.issn.0258-2724.20220763
CHEN Feng, YANG Jie, ZHANG Chong, YU Zhen, LIU Xianfeng. Numerical Analysis of Multiphysics Coupling of Grout Penetration[J]. Journal of Southwest Jiaotong University. doi: 10.3969/j.issn.0258-2724.20220763
Citation: CHEN Feng, YANG Jie, ZHANG Chong, YU Zhen, LIU Xianfeng. Numerical Analysis of Multiphysics Coupling of Grout Penetration[J]. Journal of Southwest Jiaotong University. doi: 10.3969/j.issn.0258-2724.20220763

注浆渗透扩散的多物理场耦合数值分析

doi: 10.3969/j.issn.0258-2724.20220763
基金项目: 中国铁道科学研究院基金(2020YJ105)
详细信息
    作者简介:

    陈锋(1980—),男,研究员,研究方向为路基工程及地基处理,E-mail:chenfeng7@163.com

    通讯作者:

    刘先峰(1980—),男,教授,博士,研究方向为特殊岩土力学及工程应用,E-mail:Xianfeng.liu@swjtu.edu.cn

  • 中图分类号: O242.2

Numerical Analysis of Multiphysics Coupling of Grout Penetration

  • 摘要:

    非饱和地层中的注浆渗透扩散是一个复杂的多物理场过程,为更加精确分析浆液在饱和与非饱和地层中的扩散特性,并估算注浆渗透扩散范围和注浆挤密区域,以混合理论为基础,建立一个非饱和多孔介质中的多物理场耦合模型. 通过ABAQUS二次开发,构建一种新型八节点五自由度四边形轴对称Serendipity单元,实现对注浆过程中土体变形、土体孔隙率、孔隙压力和浆液浓度分布的数值求解,以及对土体饱和度、渗透系数等状态变量的实时更新;结合一个三维轴对称注浆算例,分析浆液水灰比、注浆压力、土体初始干密度以及土体初始含水率对粉砂地层注浆效果的影响,并得到浆液水平和竖向扩散距离随不同因素变化的拟合曲线. 研究结果表明:浆液扩散范围受水灰比影响最显著,受注浆压力影响次之,受含水率和干密度影响最小;浆液扩散范围随水灰比增加而增长,水灰比大于1.0时增长显著;注浆管壁周围会形成挤密区域,浆液扩散区域内土体同时受到注浆压力的挤压和孔隙压力的支撑作用;随着远离注浆管壁,土体孔隙率在挤密区域内逐渐减小,在挤密区域外逐渐恢复,且挤密区域随注浆压力的增加而增大;研究成果可为土体注浆加固范围计算提供理论指导.

     

  • 图 1  ABAQUS调用UEL子程序计算流程

    Figure 1.  Calculation flow chart of ABAQUS calling UEL subroutine

    图 2  新型五自由度轴对称单元示意

    Figure 2.  Illustration of novel five-degree-of-freedom axisymmetric element

    图 3  数值模型的网格和几何尺寸(单位:m)

    Figure 3.  Mesh and geometry of numerical model (unit: m)

    图 4  工况6的浆液相对浓度分布

    Figure 4.  Relative concentration of grout in case 6

    图 5  工况6的孔隙液体饱和度分布

    Figure 5.  Fluid saturation of pore in case 6

    图 6  工况6的孔隙压力分布

    Figure 6.  Pore pressure in case 6

    图 7  不同水灰比浆液的浓度分布

    Figure 7.  Concentration of the grout in the surrounding soil for different water cement ratio

    图 8  不同水灰比浆液的最大扩散范围

    Figure 8.  The maximum diffusion range of grout with different water cement ratio

    图 9  水平路径上的浆液浓度分布随时间变化

    Figure 9.  Grout concentration distribution along horizontal path at different time

    图 10  不同注浆时间下的浆液扩散范围

    Figure 10.  Grout penetration range at different grouting time

    图 11  浆液扩散距离随水灰比、注浆压力、干密度和含水率变化的拟合曲线

    Figure 11.  Fitting curves of grout penetration distance with water-cement ratio, grouting pressure, dry density, and water content

    图 12  注浆土体的孔隙率变化及分布特征

    Figure 12.  Variation and distribution characteristics of porosity of soil during grouting

    表  1  数值模拟主要试验变量

    Table  1.   Main test variables for numerical simulation

    工况 注浆压力/MPa W/C 黏度时变[16]/
    (MPa·s)
    初始干密
    度/(g·cm−3
    初始含水率/%
    1 1.0 0.5 $ {\mu _{0.5}} = 407.41{{\mathrm{e}}^{0.0005t}} $ 1.6 10
    2 1.0 0.7 $ {\mu _{0.7}} = 43.864{{\mathrm{e}}^{0.0088t}} $ 1.6 10
    3 1.0 0.8 $ {\mu _{0.8}} = 28.183{{\mathrm{e}}^{0.0114t}} $ 1.6 10
    4 1.0 1.0 $ {\mu _{1.0}} = 15.632{{\mathrm{e}}^{0.0021t}} $ 1.6 10
    5 1.0 1.5 $ {\mu _{1.5}} = 10.867{{\mathrm{e}}^{0.0017t}} $ 1.6 10
    6 1.0 3.0 $ {\mu _{3.0}} = 6.8812{{\mathrm{e}}^{0.0009t}} $ 1.6 10
    7 0.3 1.0 $ {\mu _{1.0}} = 15.632{{\mathrm{e}}^{0.0021t}} $ 1.6 10
    8 0.6 1.0 $ {\mu _{1.0}} = 15.632{{\mathrm{e}}^{0.0021t}} $ 1.6 10
    9 0.9 1.0 $ {\mu _{1.0}} = 15.632{{\mathrm{e}}^{0.0021t}} $ 1.6 10
    10 1.2 1.0 $ {\mu _{1.0}} = 15.632{{\mathrm{e}}^{0.0021t}} $ 1.6 10
    11 1.0 1.0 $ {\mu _{1.0}} = 15.632{{\mathrm{e}}^{0.0021t}} $ 1.5 10
    12 1.0 1.0 $ {\mu _{1.0}} = 15.632{{\mathrm{e}}^{0.0021t}} $ 1.4 10
    13 1.0 1.0 $ {\mu _{1.0}} = 15.632{{\mathrm{e}}^{0.0021t}} $ 1.3 10
    14 1.0 1.0 $ {\mu _{1.0}} = 15.632{{\mathrm{e}}^{0.0021t}} $ 1.6 5
    15 1.0 1.0 $ {\mu _{1.0}} = 15.632{{\mathrm{e}}^{0.0021t}} $ 1.6 15
    16 1.0 1.0 $ {\mu _{1.0}} = 15.632{{\mathrm{e}}^{0.0021t}} $ 1.6 20
    下载: 导出CSV

    表  2  不同水灰比下的水泥浆液密度

    Table  2.   Density of grout with different water-cement ratios

    W/C 0.5 0.7 0.8 1.0 1.5 3.0
    水泥浆密度/
    (g·cm−3
    1.823 1.662 1.603 1.512 1.372 1.204
    下载: 导出CSV

    表  3  粉质砂土的水土特征[23]

    Table  3.   Soil and water characteristic curve of silty sand

    饱和度/% 100 94.9647 78.3324 28.0630 23.5915 20.0685 18.6214 17.3607
    基质吸力/kPa 0 15.2868 19.6129 34.2297 53.6164 91.5849 156.062 219.9930
    下载: 导出CSV

    表  4  各因素对浆液扩散距离影响的拟合曲线方程

    Table  4.   Fitting equations of grout penetration distance with various factors

    因素 最大扩散半径 最大扩散深度
    水灰比 $ y = {{ - }}0.16 + 0.55x - 0.09{x^2} $ $ \;y = {{ - }}0.09 + 0.44x - 0.08{x^2} $
    注浆
    压力
    $ y = 0.12 + 0.55x - 0.22{x^2} $ $ y = 0.16 + 0.38x - 0.14{x^2} $
    干密度 $ y = 1.21{{ - }}0.54x $ $ y = 0.93{{ - }}0.37x $
    含水率 $ y = 0.35{{ - 0}}{\text{.01}}x $ $ y = 0.34{{ - 0}}{\text{.15}}x $
    下载: 导出CSV
  • [1] 叶飞,朱合华,何川. 盾构隧道壁后注浆扩散模式及对管片的压力分析[J]. 岩土力学,2009,30(5):1307-1312. doi: 10.3969/j.issn.1000-7598.2009.05.020

    YE Fei, ZHU Hehua, HE Chuan. Back-filled grouts diffusion model and its pressure to segments of shield tunnel[J]. Rock and Soil Mechanics, 2009, 30(5): 1307-1312. doi: 10.3969/j.issn.1000-7598.2009.05.020
    [2] 周晓敏,徐衍,刘书杰,等. 金矿超深立井含水围岩注浆加固的应力场和渗流场研究[J]. 岩石力学与工程学报,2020,39(8):1611-1621.

    ZHOU Xiaomin, XU Yan, LIU Shujie, et al. Research on stress and seepage fields of surrounding rock grouting of ultra-deep shafts of gold mines[J]. Chinese Journal of Rock Mechanics and Engineering, 2020, 39(8): 1611-1621.
    [3] 张忠苗,邹健,何景愈,等. 考虑压滤效应下饱和黏土压密注浆柱扩张理论[J]. 浙江大学学报(工学版),2011,45(11):1980-1984.

    ZHANG Zhongmiao, ZOU Jian, HE Jingyu, et al. Cavity expansion theory of compaction grouting in saturated clay considering pressure filtration[J]. Journal of Zhejiang University (Engineering Science), 2011, 45(11): 1980-1984.
    [4] 刘奇,陈卫忠,袁敬强,等. 基于渗流-侵蚀理论的岩溶充填介质注浆加固效果评价[J]. 岩石力学与工程学报,2020,39(3):572-580.

    LIU Qi, CHEN Weizhong, YUAN Jingqiang, et al. Evaluation of grouting reinforcement effect for karst filling medium based on seepage-erosion theory[J]. Chinese Journal of Rock Mechanics and Engineering, 2020, 39(3): 572-580.
    [5] 马明杰,杨新安,周建,等. 非达西渗流作用下饱和黏土压密注浆扩孔理论分析[J]. 哈尔滨工业大学学报,2022,54(3):32-40.

    MA Mingjie, YANG Xin’an, ZHOU Jian, et al. Theoretical analysis of cavity expansion for compaction grouting in saturated clay under non-Darcy seepage[J]. Journal of Harbin Institute of Technology, 2022, 54(3): 32-40.
    [6] YIN Z Y, YANG J, LAOUAFA F, et al. A framework for coupled hydro-mechanical continuous modelling of gap-graded granular soils subjected to suffusion[J]. European Journal of Environmental and Civil Engineering, 2020, 27(1): 1-22.
    [7] YANG J, YIN Z Y, LAOUAFA F, et al. Analysis of suffusion in cohesionless soils with randomly distributed porosity and fines content[J]. Computers and Geotechnics, 2019, 111: 157-171. doi: 10.1016/j.compgeo.2019.03.011
    [8] YANG J, YIN Z Y, LAOUAFA F, et al. Hydromechanical modeling of granular soils considering internal erosion[J]. Canadian Geotechnical Journal, 2020, 57(2): 157-172. doi: 10.1139/cgj-2018-0653
    [9] YANG J, YIN Z Y, LAOUAFA F, et al. Three-dimensional hydromechanical modeling of internal erosion in dike-on-foundation[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2020, 44(8): 1200-1218. doi: 10.1002/nag.3057
    [10] RAATS P A C. Applications of the theory of mixtures in soil science[J]. Mathematical Modelling, 1987, 9(12): 849-856. doi: 10.1016/0270-0255(87)90003-0
    [11] DE BOER R. Contemporary progress in porous media theory[J]. Applied Mechanics Reviews, 2000, 53(12): 323-370 doi: 10.1115/1.3097333
    [12] UZUOKA R, BORJA R I. Dynamics of unsaturated poroelastic solids at finite strain[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2012, 36(13): 1535-1573. doi: 10.1002/nag.1061
    [13] ZHOU Z L, ZANG H Z, WANG S Y, et al. Filtration behaviour of cement-based grout in porous media[J]. Transport in Porous Media, 2018, 125(3):435-463. doi: 10.1007/s11242-018-1127-x
    [14] LIU X X, SHEN S L, XU Y S, et al. A diffusion model for backfill grout behind shield tunnel lining[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2021, 45(4): 457-477. doi: 10.1002/nag.3164
    [15] REVIL A, CATHLES L. Permeability of shaly sands[J]. Water Resources Research, 1999, 35(3): 651-662. doi: 10.1029/98WR02700
    [16] 雷进生. 碎石土地基注浆加固力学行为研究[D]. 武汉:中国地质大学,2013.
    [17] LIANG Y, ZHANG J, LAI Z S, et al. Temporal and spatial distribution of the grout pressure and its effects on lining segments during synchronous grouting in shield tunnelling[J]. European Journal of Environmental and Civil Engineering, 2020, 24(1): 79-96. doi: 10.1080/19648189.2017.1364299
    [18] Hibbitt, Karlsson, Sorensen. ABAQUS/Standard: User's manual[M]. Providence: [s.n.], 1997.
    [19] LADYZHENSKAYA O A. The mathematical theory of viscous incompressible flow[M]. New York: Gordon & Breach, 1969.
    [20] BREZZI F. On the existence, uniqueness and approximation of saddle-point problems arising from Lagrangian multipliers[J]. Revue Française D’automatique, Informatique, Recherche Opérationnelle Analyse Numérique, 1974, 8: 129-151.
    [21] 王勖成,邵敏. 有限单元法基本原理和数值方法[M]. 2版. 北京:清华大学出版社,1997.
    [22] YANG J. Numerical analyses of the multi-physics problem of sinkholes in the vicinity of a dike or a linear geo-structure[D]. Nantes: École centrale de Nantes, 2019.
    [23] 李文. 路基沉陷劈裂注浆处治试验研究及力学参数计算[D]. 长沙:长沙理工大学,2013.
    [24] 林鸿州,李广信,于玉贞,等. 基质吸力对非饱和土抗剪强度的影响[J]. 岩土力学,2007,28(9):1931-1936. doi: 10.3969/j.issn.1000-7598.2007.09.031

    LIN Hongzhou, LI Guangxin, YU Yuzhen, et al. Influence of matric suction on shear strength behavior of unsaturated soils[J]. Rock and Soil Mechanics, 2007, 28(9): 1931-1936. doi: 10.3969/j.issn.1000-7598.2007.09.031
  • 加载中
图(12) / 表(4)
计量
  • 文章访问数:  52
  • HTML全文浏览量:  22
  • PDF下载量:  6
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-08-03
  • 修回日期:  2023-01-20
  • 网络出版日期:  2024-03-12

目录

    /

    返回文章
    返回