留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

胶结型含可燃冰砂土剪切特性的离散元模拟

王璇 徐明

王璇, 徐明. 胶结型含可燃冰砂土剪切特性的离散元模拟[J]. 工程力学, 2021, 38(2): 44-51. doi: 10.6052/j.issn.1000-4750.2020.03.0174
引用本文: 王璇, 徐明. 胶结型含可燃冰砂土剪切特性的离散元模拟[J]. 工程力学, 2021, 38(2): 44-51. doi: 10.6052/j.issn.1000-4750.2020.03.0174
Xuan WANG, Ming XU. DISCRETE ELEMENT SIMULATION OF THE SHEAR BEHAVIOR OF CEMENTED METHANE HYDRATE-BEARING SANDS[J]. Engineering Mechanics, 2021, 38(2): 44-51. doi: 10.6052/j.issn.1000-4750.2020.03.0174
Citation: Xuan WANG, Ming XU. DISCRETE ELEMENT SIMULATION OF THE SHEAR BEHAVIOR OF CEMENTED METHANE HYDRATE-BEARING SANDS[J]. Engineering Mechanics, 2021, 38(2): 44-51. doi: 10.6052/j.issn.1000-4750.2020.03.0174

胶结型含可燃冰砂土剪切特性的离散元模拟

doi: 10.6052/j.issn.1000-4750.2020.03.0174
基金项目: 国家自然科学基金项目(51978382);国家973计划项目(2014CB047003)
详细信息
    作者简介:

    王璇(1997−),女,安徽人,硕士生,主要从事含可燃冰砂土的离散元数值模拟研究(E-mail: wangxuan19@mails.tsinghua.edu.cn)

    通讯作者: 徐明(1974−),男,湖北人,副教授,博士,博导,主要从事岩土力学试验研究及理论分析(E-mail: mingxu@mail.tsinghua.edu.cn)

DISCRETE ELEMENT SIMULATION OF THE SHEAR BEHAVIOR OF CEMENTED METHANE HYDRATE-BEARING SANDS

  • 摘要: 可燃冰是一种新型清洁能源,广泛分布于深海土体和常年冻土中,对其开采需要深入认识含可燃冰土体的力学特性。该文提出了一种胶结型含可燃冰砂土离散元模型建立方法,并利用该模型完成了排水双轴剪切试验的模拟。通过与前人试验数据对比,验证了该离散元模型的准确性,结合该模型分析了剪切过程中含可燃冰砂土的宏细观变化规律。结果表明:砂土抗剪强度、应变软化特性及剪胀性均随可燃冰饱和度的增加而增加;可燃冰胶结断裂数量的演化与偏应力紧密相关,偏应力达到峰值时胶结断裂数量增长速率最快;剪切过程中剪切带内外土体的胶结断裂、颗粒运动、孔隙率变化均表现出明显的差异。通过对模型粘结断裂率、剪切带内外局部孔隙率等细观指标的分析,对土体的偏应力变化、体积变化等宏观现象的机理做出了更进一步的阐释。
  • 图  1  胶结型含可燃冰砂土离散元模型

    Figure  1.  DEM model of cemented hydrate-bearing sands

    图  2  不同颗粒粘结率Rb/可燃冰饱和度SMH下应力-应变关系

    Figure  2.  Stress-strain curves with different bonding ratios Rb/hydrate saturation SMH

    图  3  不同颗粒粘结率Rb/可燃冰饱和度SMH下体应变关系

    Figure  3.  Volumetric curves with different bonding ratios Rb/hydrate saturation SMH

    图  4  不同颗粒粘结率Rb/可燃冰饱和度SMH下内摩擦角φ及粘聚力c变化

    Figure  4.  Friction angle φ and cohesion c with different bonding ratios Rb/hydrate saturation SMH

    图  5  不同颗粒粘结率Rb下粘结断裂率Rbk-轴应变变化关系

    Figure  5.  Relationship between bond breaking ratio Rbk and axial strain with different bonding ratios Rb

    图  6  不同围压下粘结断裂率Rbk-轴应变变化关系

    Figure  6.  Relationship between bond breaking ratio Rbk and axial strain under different confining pressures

    图  7  不同轴应变下砂土粘结断裂分布(Rb=55.1%)

    Figure  7.  Bond breaking distribution at different axial strains (Rb=55.1%)

    图  8  不同轴应变下砂土颗粒转角(Rb=55.1%)

    Figure  8.  Rotation of particles at different axial strains (Rb=55.1%)

    图  9  不同轴应变下砂土局部孔隙率(Rb=55.1%)

    Figure  9.  Local porosity of sand at different axial strains (Rb=55.1%)

    图  10  剪切带内外颗粒位移对比(ε1=14%)

    Figure  10.  Particle displacement inside and outside shear band (ε1=14%)

    图  11  ε1=14%时不同颗粒粘结率Rb砂土剪切带内外孔隙率

    Figure  11.  Porosity of sands inside and outside shear band with different bonding ratios Rb at ε1=14%

    表  1  离散元模型参数取值

    Table  1.   Parameters of DEM model

    参数取值
    模型尺寸/mm160×80
    圆盘单元总数20 117
    圆盘单元密度/(kg/m3)2600
    圆盘单元直径/mm0.6~1.0
    成样时圆盘单元摩擦系数0.2
    成样后圆盘单元摩擦系数1.0
    初始模型孔隙率/(%)21
    圆盘单元切向刚度/(N/m)1.6×109
    圆盘单元法向刚度/(N/m)2.4×109
    墙单元切向刚度/(N/m)1.6×109
    墙单元法向刚度/(N/m)2.4×109
    切向粘结强度/Pa1.1×104
    法向粘结强度/Pa1.1×104
    颗粒粘结率Rb/(%)取值等于可燃冰饱和度SMH
    下载: 导出CSV
  • [1] Makogon Y F, Holditch S A, Makogon T Y. Natural gas-hydrates — A potential energy source for the 21st Century [J]. Journal of Petroleum Science & Engineering, 2007, 56(1/2/3): 14 − 31.
    [2] 马小娟, 呼方涛, 王锦华. 浅析可燃冰的研究现状与发展前景[J]. 科技信息, 2011(7): 351 − 382.

    Ma X J, Hu F T, Wang J H. Research status and development prospect of combustible ice [J]. Science & Technology Information, 2011(7): 351 − 382. (in Chinese)
    [3] Locat J, Lee H J. Submarine landslides: advances and challenges [J]. Canadian Geotechnical Journal, 2002, 39(1): 193 − 212. doi:  10.1139/t01-089
    [4] Soga K, Klar A. Characterisation and engineering properties of methane hydrate soils [C]// The Second International Workshop on Characterisation and Engineering Properties of Natural Soils. Singapore: Taylor and Francis Group, 2007: 2591 − 2642.
    [5] Masui A, Haneda H, Ogata Y. Effects of methane hydrate formation on shear strength of synthetic methane hydrate sediments [C]// Proceedings of the 15th International Offshore and Polar Engineering Conference. Seoul, Korea: The International Society of Offshore and Polar Engineers, 2005: 364 − 369.
    [6] Li Y, Liu W. Experimental Research on the Mechanical Properties of Methane Hydrate-Ice Mixtures [J]. Energies, 2012, 5: 181 − 192.
    [7] Miyazaki K, Masui A, Sakamoto Y, et al. Triaxial compressive properties of artificial methane-hydrate-bearing sediment [J]. Journal of Geophysical Research: Solid Earth, 2011, 116(B06102): 1 − 12. doi:  10.1029/2010JB008049
    [8] Pinkert S, Grozic J L H. Failure Mechanisms in Cemented Hydrate-Bearing Sands [J]. Journal of Chemical & Engineering Data, 2015, 60(2): 376 − 382. doi:  10.1021/je500638c
    [9] Hyodo M, Yoneda J, Yoshimoto N, et al. Mechanical and dissociation properties of methane hydrate-bearing sand in deep seabed [J]. Soils Found, 2013, 53(2): 299 − 314. doi:  10.1016/j.sandf.2013.02.010
    [10] Song Y, Zhu Y, Liu W, et al. Experimental research on the mechanical properties of methane hydrate-bearing sediments during hydrate dissociation [J]. Marine & Petroleum Geology, 2014, 51(Complete): 70 − 78.
    [11] Zhang X H, Luo D S, Lu X B, et al. Mechanical properties of gas hydrate-bearing sediments during hydrate dissociation [J]. Acta Mechanica Sinica, 2018, 34(2): 266 − 374. doi:  10.1007/s10409-017-0699-y
    [12] Cundall P A, Strack O D L. A discrete numerical model for granular assemblies [J]. Géotechnique, 1979, 29(1): 47 − 65. doi:  10.1680/geot.1979.29.1.47
    [13] 齐念, 叶继红. 弹性DEM方法在杆系结构中的应用研究[J]. 工程力学, 2017, 34(7): 11 − 20.

    Qi N, Ye J H. Application of elastic discrete element method in member structures [J]. Engineering Mechanics, 2017, 34(7): 11 − 20. (in Chinese)
    [14] 王蕴嘉, 周梦佳, 宋二祥. 考虑颗粒破碎的堆石料湿化变形特性离散元模拟研究[J]. 工程力学, 2018, 35(增刊 1): 217 − 222. doi:  10.6052/j.issn.1000-4750.2017.05.S043

    Wang Y J, Zhou M J, Song E X. DEM simulation of wetting deformation characteristics of rockfill considering particle breakage [J]. Engineering Mechanics, 2018, 35(Suppl 1): 217 − 222. (in Chinese) doi:  10.6052/j.issn.1000-4750.2017.05.S043
    [15] 孟凡净, 刘焜, 吴华伟. 颗粒流润滑剪切膨胀的力学机制研究[J]. 工程力学, 2018, 35(8): 236 − 244.

    Meng F J, Liu K, Wu H W. Research on mechanical mechanism of shear dilatancy for granular flow lubrication [J]. Engineering Mechanics, 2018, 35(8): 236 − 244. (in Chinese)
    [16] Brugada J, Cheng Y P, Soga K, et al. Discrete element modelling of geomechanical behaviour of methane hydrate soils with pore-filling hydrate distribution [J]. Granular Matter, 2010, 12(5): 517 − 525. doi:  10.1007/s10035-010-0210-y
    [17] 蒋明镜, 肖俞, 朱方园. 深海能源土微观力学胶结模型及参数研究[J]. 岩土工程学报, 2012, 34(9): 1574 − 1583.

    Jiang M J, Xiao Y, Zhu F Y. Micro-bond contact model and its parameters for the deep-sea methane hydrate bearing soils [J]. Chinese Journal of Geotechnical Engineering, 2012, 34(9): 1574 − 1583. (in Chinese)
    [18] 蒋明镜, 肖俞, 朱方园. 深海能源土宏观力学性质离散元数值模拟分析[J]. 岩土工程学报, 2013, 35(1): 157 − 163.

    Jiang M J, Xiao Y, Zhu F Y. Numerical simulation of macro-mechanical properties of deep-sea methane hydrate bearing soils by DEM [J]. Chinese Journal of Geotechnical Engineering, 2013, 35(1): 157 − 163. (in Chinese)
    [19] Jiang M, Zhu F, Liu F, et al. A bond contact model for methane hydrate-bearing sediments with interparticle cementation [J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2014, 38(17): 1823 − 1854. doi:  10.1002/nag.2283
    [20] Xu M, Song E, Jiang H, et al. DEM simulation of the undrained shear behavior of sand containing dissociated gas hydrate [J]. Granular matter, 2016, 18(79): 1 − 13. doi:  10.1007/s10035-016-0675-4
    [21] Hong J, Xu M. DEM study on the undrained mechanical behavior of gassy sand [J]. Acta Geotechnica, 2020, 15(8): 2179 − 2193. doi:  10.1007/s11440-019-00910-x
    [22] Itasca Consulting Group Inc. User's Manual for PFC2D [M].Minnesota: Itasca Consulting Group Inc., 2005.
    [23] 胡高伟, 李承峰, 业渝光, 等. 沉积物孔隙空间天然气水合物微观分布观测[J]. 地球物理学报, 2014, 57(5): 1675 − 1682. doi:  10.6038/cjg20140530

    Hu G W, Li C F, Ye Y G, et al. Observation of gas hydrate distribution in sediment pore space [J]. Chinese Journal of Geophysics, 2014, 57(5): 1675 − 1682. (in Chinese) doi:  10.6038/cjg20140530
    [24] Zhou W, Yang L, Ma G, et al. DEM modeling of shear bands in crushable and irregularly shaped granular materials [J]. Granular Matter, 2017, 19(25): 1 − 12. doi:  10.1007/s10035-017-0712-y
    [25] 杨智勇, 曹子君, 李典庆, 等. 颗粒接触摩擦系数空间变异性对颗粒流双轴数值试验的影响[J]. 工程力学, 2017, 34(5): 235 − 246.

    Yang Z Y, Cao Z J, Li D Q, et al. Effect of spatially variable friction coefficient of granular materials on its macro-mechanical behaviors using biaxial compression numerical simulation [J]. Engineering Mechanics, 2017, 34(5): 235 − 246. (in Chinese)
  • [1] 叶继红, 王佳.  基于GPU的杆系离散元并行算法在大型工程结构中的应用 . 工程力学, 2021, 38(2): 1-7. doi: 10.6052/j.issn.1000-4750.2020.07.ST03
    [2] 张帅, 程晓辉, 王天麟.  非等向固结砂土极小应变刚度的超弹性模型 . 工程力学, 2020, 37(1): 145-151. doi: 10.6052/j.issn.1000-4750.2019.02.0058
    [3] 王鹏鹏, 郭晓霞, 桑勇, 邵龙潭, 陈之祥, 赵博雅.  基于数字图像相关技术的砂土全场变形测量及其DEM数值模拟 . 工程力学, 2020, 37(1): 239-247. doi: 10.6052/j.issn.1000-4750.2019.02.0050
    [4] 孟凡净, 刘焜, 吴华伟.  颗粒流润滑剪切膨胀的力学机制研究 . 工程力学, 2018, 35(8): 236-244. doi: 10.6052/j.issn.1000-4750.2017.04.0282
    [5] 井国庆, 王子杰, 施晓毅.  多围压下三轴压缩试验与不可破裂颗粒离散元法分析 . 工程力学, 2015, 32(10): 82-88. doi: 10.6052/j.issn.1000-4750.2014.03.0205
    [6] 邵帅, 季顺迎.  块石空间分布对土石混合体边坡稳定性的影响 . 工程力学, 2014, 31(2): 177-183. doi: 10.6052/j.issn.1000-4750.2012.04.0248
    [7] 李学丰, 黄茂松, 钱建固.  基于非共轴理论的各向异性砂土应变局部化分析 . 工程力学, 2014, 31(3): 205-211. doi: 10.6052/j.issn.1000-4750.2012.10.0771
    [8] 曹光栩, 徐明, 宋二祥.  反映粗粒料应力路径相关性的一种应变硬化模型 . 工程力学, 2013, 30(4): 83-88. doi: 10.6052/j.issn.1000-4750.2011.12.0814
    [9] 王学滨, 杜亚志, 潘一山.  考虑一阶和二阶位移梯度的数字图像相关方法在剪切带测量中的比较 . 工程力学, 2013, 30(7): 282-287. doi: 10.6052/j.issn.1000-4750.2012.03.0197
    [10] 喻葭临, 于玉贞, 张丙印, 吕 禾.  基于扩展有限元方法的界面接触算法 . 工程力学, 2011, 28(4): 13-017.
    [11] 辛克贵, 何铭华.  分布粘聚元的系统理论研究 . 工程力学, 2011, 28(增刊Ⅱ): 109-128.
    [12] 侯 健.  考虑块体间碰撞作用的混凝土框架结构空间倒塌反应分析 . 工程力学, 2010, 27(6): 89-097.
    [13] 常 在, 杨 军, 程晓辉.  砂土强度和剪胀性的颗粒力学分析 . 工程力学, 2010, 27(4): 95-104.
    [14] 李增志;别社安;任增金.  抛石防波堤稳定性的离散单元法分析 . 工程力学, 2009, 26(增刊 I): 111-114.
    [15] 冀 宾, 陈万吉, 赵 杰.  岩土软化行为的网格依赖性问题研究 . 工程力学, 2008, 25(8): 0-041.
    [16] 熊保林, 邵龙潭.  考虑主应力轴旋转Gudehus-Bauer亚塑性模型的改进 . 工程力学, 2008, 25(1): 0-132.
    [17] 王学滨.  岩样单轴压缩峰后泊松比理论研究 . 工程力学, 2006, 23(4): 99-103.
    [18] 金伟良, 方韬.  钢筋混凝土框架结构破坏性能的离散单元法模拟 . 工程力学, 2005, 22(4): 67-73.
    [19] 王学滨, 潘一山, 马瑾.  剪切带内部应变(率)分析及基于能量准则的失稳判据 . 工程力学, 2003, 20(2): 111-115.
    [20] 邓子胜.  砂土振动固结抗剪强度变化规律的试验研究 . 工程力学, 2000, 17(4): 94-98.
  • 加载中
图(11) / 表 (1)
计量
  • 文章访问数:  57
  • HTML全文浏览量:  11
  • PDF下载量:  23
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-03-20
  • 修回日期:  2020-06-30
  • 网络出版日期:  2021-01-16
  • 刊出日期:  2021-01-16

胶结型含可燃冰砂土剪切特性的离散元模拟

doi: 10.6052/j.issn.1000-4750.2020.03.0174
    基金项目:  国家自然科学基金项目(51978382);国家973计划项目(2014CB047003)
    作者简介:

    王璇(1997−),女,安徽人,硕士生,主要从事含可燃冰砂土的离散元数值模拟研究(E-mail: wangxuan19@mails.tsinghua.edu.cn)

    通讯作者: 徐明(1974−),男,湖北人,副教授,博士,博导,主要从事岩土力学试验研究及理论分析(E-mail: mingxu@mail.tsinghua.edu.cn)

摘要: 可燃冰是一种新型清洁能源,广泛分布于深海土体和常年冻土中,对其开采需要深入认识含可燃冰土体的力学特性。该文提出了一种胶结型含可燃冰砂土离散元模型建立方法,并利用该模型完成了排水双轴剪切试验的模拟。通过与前人试验数据对比,验证了该离散元模型的准确性,结合该模型分析了剪切过程中含可燃冰砂土的宏细观变化规律。结果表明:砂土抗剪强度、应变软化特性及剪胀性均随可燃冰饱和度的增加而增加;可燃冰胶结断裂数量的演化与偏应力紧密相关,偏应力达到峰值时胶结断裂数量增长速率最快;剪切过程中剪切带内外土体的胶结断裂、颗粒运动、孔隙率变化均表现出明显的差异。通过对模型粘结断裂率、剪切带内外局部孔隙率等细观指标的分析,对土体的偏应力变化、体积变化等宏观现象的机理做出了更进一步的阐释。

English Abstract

王璇, 徐明. 胶结型含可燃冰砂土剪切特性的离散元模拟[J]. 工程力学, 2021, 38(2): 44-51. doi: 10.6052/j.issn.1000-4750.2020.03.0174
引用本文: 王璇, 徐明. 胶结型含可燃冰砂土剪切特性的离散元模拟[J]. 工程力学, 2021, 38(2): 44-51. doi: 10.6052/j.issn.1000-4750.2020.03.0174
Xuan WANG, Ming XU. DISCRETE ELEMENT SIMULATION OF THE SHEAR BEHAVIOR OF CEMENTED METHANE HYDRATE-BEARING SANDS[J]. Engineering Mechanics, 2021, 38(2): 44-51. doi: 10.6052/j.issn.1000-4750.2020.03.0174
Citation: Xuan WANG, Ming XU. DISCRETE ELEMENT SIMULATION OF THE SHEAR BEHAVIOR OF CEMENTED METHANE HYDRATE-BEARING SANDS[J]. Engineering Mechanics, 2021, 38(2): 44-51. doi: 10.6052/j.issn.1000-4750.2020.03.0174
  • 可燃冰的主要成分为甲烷水合物,是由甲烷分子和水分子在高压、低温条件下生成的亚稳态固体物质,广泛分布于深海和常年冻土区域[1]。作为一种燃烧热高、污染少的新型绿色能源,可燃冰在近些年来得到了人们的广泛关注。目前,许多国家已经制定了可燃冰资源开发规划,并开始大力发展可燃冰开采技术[2]

    然而,在深海开采可燃冰过程中极易引发海底坍塌、滑坡,进而诱发地震、海啸等地质灾难[3],造成巨大的人员伤亡和财产损失。为了进一步完善可燃冰开采技术,避免可燃冰开采过程中灾害的发生,对含可燃冰土体力学特性的深入研究具有非常重要的意义。

    可燃冰在砂土颗粒间的分布主要有填充、胶结、包裹、骨架等4种形态[4],不同的分布形态会使含可燃冰砂土表现出不同的力学性质。对此,学者们已经进行了有限的试验研究。其中,Masui等[5]分别向含冰砂和含水砂中通入甲烷气体,诱导生成填充型和胶结型含可燃冰砂土,并对其实施三轴试验,探究了可燃冰形态差异对于砂土抗剪强度和刚度的影响;Li等[6]通过三轴试验探究了可燃冰饱和度和围压对于填充型含可燃冰砂土内摩擦角与粘聚力的影响;Miyazaki等[7]通过三轴试验探究了可燃冰饱和度、围压、温度等对于胶结型含可燃冰砂土抗剪强度的影响;Pinkert等[8]对于剪切过程中胶结型含可燃冰砂土的胶结破坏形式进行了探究。此外,可燃冰仅能在高压、低温条件下保持稳定,在开采过程中极易发生分解生成甲烷气体,因此Hyodo等[9]、Song等[10]、Zhang等[11]通过试验研究了可燃冰分解对于含可燃冰土体力学性质的影响。上述试验初步揭示了含可燃冰砂土的宏观力学性质,但对于加载过程中土体的微观变化难以进行有效的观测,此外,真实试验受到取样制样难、试验成本高等制约。

    离散元法是一种用于反映颗粒材料力学性质和颗粒间相互作用的数值模拟方法,于1979年由Cundall等[12]提出,目前在结构、岩土、机械等多个领域得到了广泛的应用[13-15]。近年来,对含可燃冰砂土的离散元模拟已经成为一种新的研究思路。Brugada等[16]通过在砂土颗粒间生成直径更小的可燃冰颗粒,建立填充型含可燃冰砂土的离散元模型,探究了填充型可燃冰对于砂土力学性质的影响;蒋明镜等[17-19]通过分析前人试验结果并类比水泥、环氧树脂等其他胶结材料,提出了一种含可燃冰砂土微观力学胶结模型的建立方法,并在此基础上使用PFC2D软件,采用在所有颗粒接触处添加平行粘结的方式对胶结型含可燃冰砂土进行了离散元模拟;Xu等[20-21]利用离散元方法模拟了含可燃冰解离气体砂土的不排水剪切试验。然而,目前已有的成果主要集中在反映含可燃冰砂土的宏观力学性质上,而对于一些宏观现象仍缺乏相应的细微观阐释。此外,目前已有模型往往通过改变粘结强度[17-19]反映可燃冰含量的变化,然而二者之间的联系较为抽象,函数关系难以确定。

    本文在前人研究的基础上,提出一种新的胶结型含可燃冰砂土离散元模拟方法,通过在颗粒间添加接触粘结模拟胶结型可燃冰的存在。该方法通过改变添加粘结的数量反映可燃冰含量的变化,使得模型与真实土体之间的对应关系更为直观。随后,本文利用PFC2D软件对上述模型进行了排水剪切试验模拟,通过与Masui等[5]试验对比验证了模型的合理性。在此基础上,本文探讨了宏观力学特性和细观现象之间的对应机理及关联,通过分析粘结断裂、剪切带内外土颗粒运动等细观现象,对偏应力变化、体胀等宏观特性进行了进一步阐释。

    • 本文利用PFC2D软件建立可燃冰胶结模型,该模型包括砂土颗粒模型和可燃冰模型两部分。

      其中,砂土颗粒选用线性接触模型,法向刚度和切向刚度均为常数。模型不考虑砂土颗粒形状对于土体性质的影响,所有砂土颗粒均采用单个具有单位厚度的圆盘单元来模拟。模型尺寸取为160 mm×80 mm。

      对于砂土颗粒间以胶结形态存在的可燃冰,选取PFC2D软件中内置的粘结模型进行模拟。PFC2D中常用的粘结模型包括平行粘结和接触粘结[22],通过将分别添加上述两种粘结模拟所得的双轴试验曲线与前人真实试验曲线对比,确定最终选用的粘结模型为接触粘结。

      可燃冰饱和度SMH指土体中可燃冰体积占土体总孔隙的百分比,它是影响含可燃冰砂土力学性质的一项重要指标。在本胶结模型中,利用颗粒粘结率Rb来反映真实土体中可燃冰饱和度SMH的变化。颗粒粘结率Rb定义如式(1)所示:

      $$ {R_{\rm{b}}} = \frac{{{N_{{\rm{b}}}}}}{{{N_{{\rm{t}}}}}} \times 100\text{%} $$ (1)

      式中:Nt表示PFC模型的圆盘单元间接触总数量;Nb表示添加接触粘结的数量。模型中每个接触粘结的添加位置随机,以此模拟天然条件下可燃冰生成位置随机的特点[23]。由于真实土体中可燃冰含量的多少与模型中添加粘结的数量直接相关,因此令模型颗粒粘结率Rb的取值与真实土体中可燃冰饱和度Rb的大小相同,以便直观反映二者之间的对应关系。

    • 为了便于同已有试验数据进行比较,本模型在建模和参数确定过程中均选用Masui等[5]所进行的人工合成胶结型含可燃冰砂土三轴排水试验作为参照。在Masui试验中,选取初始水饱和度为30%~70%的Toyoura砂制作样品,对该含水砂土样品施加指定温度和围压,随后通入甲烷气体诱导生成甲烷水合物,诱导完成后在保持围压不变的情况下对试样实施三轴剪切[5]

      对于离散元模型,首先定义上下左右4个方向的墙单元,然后用膨胀法生成20 117个颗粒直径在0.6 mm~1.0 mm随机分布的圆盘单元,由此得到初始试样。该模型取墙单元刚度与圆盘单元刚度相等,法向刚度为切向刚度的1.5倍。确定圆盘单元摩擦系数、圆盘单元刚度等取值时,选用Masui等[5]试验所得的纯砂应力-应变曲线进行参数标定,调整模型参数使离散元纯砂模型剪切强度、刚度等与真实试验吻合。最终确定墙单元与圆盘单元的详细参数如表1所示。

      表 1  离散元模型参数取值

      Table 1.  Parameters of DEM model

      参数取值
      模型尺寸/mm160×80
      圆盘单元总数20 117
      圆盘单元密度/(kg/m3)2600
      圆盘单元直径/mm0.6~1.0
      成样时圆盘单元摩擦系数0.2
      成样后圆盘单元摩擦系数1.0
      初始模型孔隙率/(%)21
      圆盘单元切向刚度/(N/m)1.6×109
      圆盘单元法向刚度/(N/m)2.4×109
      墙单元切向刚度/(N/m)1.6×109
      墙单元法向刚度/(N/m)2.4×109
      切向粘结强度/Pa1.1×104
      法向粘结强度/Pa1.1×104
      颗粒粘结率Rb/(%)取值等于可燃冰饱和度SMH

      为模拟Masui等[5]试验中诱导含可燃冰砂土的生成步骤,在形成离散元初始试样后,利用伺服机制为试样施加指定围压,然后在颗粒间接触处添加接触粘结,用于模拟胶结型可燃冰的存在。

      接触粘结的参数包括法向和切向强度,对于本模型中的粘结强度参数,选取Masui等[5]试验中SMH=55.1%的含可燃冰砂土应力-应变曲线进行标定,调整参数使Rb=55.1%时离散元模型抗剪强度与真实试验结果相同,确定模型切向与法向粘结强度均取1.1×104 Pa。经过以上步骤所得的离散元模型示意图如图1

      图  1  胶结型含可燃冰砂土离散元模型

      Figure 1.  DEM model of cemented hydrate-bearing sands

      最后,在保持围压的情况下控制上下侧墙单元的位移,墙单元运动形成对土单元的压力,以此实施对已添加粘结试样的加载,模拟双轴剪切试验过程。离散元模型中不直接模拟水的存在,施加围压和剪切的过程直接作用在土单元骨架上,总应力等于有效应力,从而模拟排水条件。通过实时记录各墙单元的位移值,可计算出加载过程中试样的轴应变和体应变变化;通过实时记录各墙单元所受法向应力,可计算出加载过程中试样的应力变化。

      通过改变颗粒粘结率Rb、围压等参数,探究不同条件下胶结型含可燃冰砂土的宏观与细观特性。

    • 首先探究模型颗粒粘结率Rb变化对于砂土抗剪强度与体应变的影响。对前文所述的离散元模型施加围压σ3=1 MPa,模拟得不同颗粒粘结率Rb(0%、25.7%、40.7%和55.1%)条件下的胶结型含可燃冰砂土应力-应变曲线如图2所示,体应变曲线如图3所示。图中同时给出了同等试验条件下Masui等[5]三轴试验的结果对比。

      图  2  不同颗粒粘结率Rb/可燃冰饱和度SMH下应力-应变关系

      Figure 2.  Stress-strain curves with different bonding ratios Rb/hydrate saturation SMH

      图  3  不同颗粒粘结率Rb/可燃冰饱和度SMH下体应变关系

      Figure 3.  Volumetric curves with different bonding ratios Rb/hydrate saturation SMH

      可以看出,通过改变模型颗粒粘结率Rb可以较好地定性反映不同可燃冰饱和度SMH条件下含可燃冰砂土排水剪切力学特性的变化,且具体表现为:

      1) 胶结型含可燃冰砂土峰值强度与刚度随SMH/Rb的增大而增大,并在峰值强度后表现出更强烈的应变软化现象;

      2) 试样加载过程中先发生体缩后发生体胀,SMH/Rb越高,试样的剪胀性越大。

      从定量上看,该模型模拟出的土体峰值强度大小与真实试验结果较为吻合,而体应变偏小,造成差异的原因可能包括未考虑砂土颗粒形状影响、模型仅限于二维模拟等。

      进一步地,探究排水剪切过程中胶结型含可燃冰砂土内摩擦角φ与粘聚力c值的变化规律。参照Masui等[5]试验,取离散元模型颗粒粘结率Rb分别为0%、25.4%和34.25%,并对每个Rb取值下的模型分别施加1 MPa、2 MPa、3 MPa及5 MPa围压进行排水双轴试验模拟,绘制摩尔圆获得破坏包络线,计算得内摩擦角φ与粘聚力c值,试验与模拟结果对比如图4所示。

      图  4  不同颗粒粘结率Rb/可燃冰饱和度SMH下内摩擦角φ及粘聚力c变化

      Figure 4.  Friction angle φ and cohesion c with different bonding ratios Rb/hydrate saturation SMH

      对比可知,本文提出的离散元模型模拟结果与Masui试验[5]结果较为吻合,二者均反映出以下规律:胶结型含可燃冰砂土在排水剪切过程中,随着SMH/Rb的增加,内摩擦角φ大小变化不大,而粘聚力c的值呈显著的增加趋势。由此说明和纯砂相比,胶结型可燃冰主要通过提高砂土的粘聚力进而提高抗剪强度,而其对于内摩擦角的影响较小。

      考虑上述现象中可燃冰的存在带来的影响,可以作出如下解释:由于可燃冰主要作为胶结物质存在于砂土颗粒之间,因此仅能提供颗粒与颗粒之间的粘聚作用,可燃冰饱和度越大,颗粒接触处存在胶结物质的比例越高,因此粘聚力c值越高,且初始刚度越大;而可燃冰的存在对土颗粒本身的颗粒咬合、摩擦系数等影响不大,因此可燃冰饱和度变化对于内摩擦角φ的影响较小;此外,含可燃冰砂土在剪切过程中表现出比纯砂更大的剪胀性,其原因是胶结型可燃冰的存在能将多个砂土颗粒粘结成为较大的颗粒簇,在剪切错动中更容易形成较大的孔隙。

    • 研究表明,胶结型含可燃冰砂土在排水剪切过程中发生可燃冰的破碎或脱粘[8],而在本离散元模型中对应为接触粘结的断裂。接触粘结的断裂标准为颗粒粘结处受到的切向应力大于切向粘结强度,或者法向应力大于法向粘结强度。为了定量描述加载过程中模型接触粘结断裂的数量变化规律,此处定义粘结断裂率Rbk如式(2):

      $$ {R_{{\rm{bk}}}} = \frac{{{N_{{\rm{br}}}}}}{{{N_{{\rm{b}}}}}} \times 100 \text{%} $$ (2)

      式中,Nbr为发生剪切后断裂的接触粘结数量。

      图5给出了围压σ3=1 MPa,模型颗粒粘结率Rb分别为0%、25.7%、40.7%和55.1%条件下的应力-应变曲线和粘结断裂率Rbk变化曲线对比(当Rb=0%时,模型不添加粘结,故不存在对应的粘结断裂率Rbk变化曲线)。

      图  5  不同颗粒粘结率Rb下粘结断裂率Rbk-轴应变变化关系

      Figure 5.  Relationship between bond breaking ratio Rbk and axial strain with different bonding ratios Rb

      图5所示,对于不同颗粒粘结率Rb的模型,均在加载至ε1=1%时才开始发生粘结断裂,且其后加载至相同轴应变条件下,Rb越高,土体的粘结断裂率Rbk越高。此外,由图5可以看出模型中粘结的断裂情况与偏应力的变化趋势密切相关。利用该结果可以解释胶结型含可燃冰砂土剪切破坏过程中宏观与细观特性之间的对应关系,其大致可以分为3个阶段:

      1) 当轴应变小于1%左右时,接触粘结几乎没有发生破坏,粘结断裂率Rbk约等于0,即土体中可燃冰胶结物完好,此时偏应力近似随轴应变的增加而线性增大,该阶段可近似看作弹性阶段,土体尚未发生破坏。

      2) 当轴应变到达1%~8%左右时,粘结断裂率Rbk进入明显的快速增长阶段,即可燃冰的胶结作用受到大量破坏,此时偏应力先达到峰值并随后开始下降,土体表现出应变软化的特性。当偏应力达到峰值时,Rbk曲线的斜率近似达到最大值。在这一阶段,土体发生剪切破坏。

      3) 当轴应变达到约8%时,粘结断裂率Rbk开始进入缓慢上升阶段,其增长速度不断下降并最终趋向于0;与此对应,偏应力曲线在这一阶段下降缓慢,至加载结束时,偏应力几乎不再变化。

      图6给出了给定颗粒粘结率Rb=55.1%,分别施加围压σ3=1 MPa、σ3=2 MPa和σ3=3 MPa时应力-应变曲线和粘结断裂率Rbk变化曲线的对比情况。可以看出,不同围压条件下的模型破坏同样表现为前文所述的3个阶段,且施加围压越大,Rbk越高,即可燃冰胶结物的破坏数量越多。

      图  6  不同围压下粘结断裂率Rbk-轴应变变化关系

      Figure 6.  Relationship between bond breaking ratio Rbk and axial strain under different confining pressures

    • 在确定粘结断裂数量的基础上,进一步研究粘结断裂的位置分布规律。以围压σ3=3 MPa,颗粒粘结率Rb=55.1%为例进行分析。图7给出了加载至轴应变ε1=0%、ε1=3.5%、ε1=7%、ε1=14% 这4个时刻的模型粘结断裂分布情况,其中灰色短线表示尚未发生断裂的粘结,红色短线表示已经断裂的粘结。由该图可以看出,随着轴应变的增加,发生断裂的粘结数逐渐增多。通过分析不同轴应变时粘结断裂的分布情况,可以大致反映出加载过程中剪切带的形状和形成过程。对于该模型,剪切过程中形成多条呈“X”状交叉的剪切带,粘结的断裂主要发生在剪切带内部。由此可以推断,对于胶结型含可燃冰砂土,排水剪切过程中可燃冰胶结物的破坏并不是均匀分布的,其破坏主要分布于土体剪切带的内部。

      图  7  不同轴应变下砂土粘结断裂分布(Rb=55.1%)

      Figure 7.  Bond breaking distribution at different axial strains (Rb=55.1%)

      图6可知,围压σ3=3 MPa,颗粒粘结率Rb=55.1%,加载至ε1=14%时,砂土粘结断裂率Rbk仅为20%左右,而曲线增长趋势已非常缓慢。该时刻对应于图7(d),剪切带内部的绝大部分粘结已发生破坏,其中少量未破坏的粘结可使相邻的几个圆盘单元形成一个较大的颗粒簇,而剪切带外部粘结在加载过程中破坏非常少。图6图7对比可知,砂土粘结断裂率Rbk的增长主要取决于剪切带内部粘结的断裂情况,当剪切带内部粘结充分破坏时,粘结断裂率Rbk便几乎不再继续增加。

      图8给出了围压σ3=3 MPa,颗粒粘结率Rb=55.1%条件下,加载至轴应变ε1=0%、ε1=3.5%、ε1=7%、ε1=14% 这4个时刻砂土颗粒的转动示意图。可以看出,在加载至轴应变大于7%后,剪切带内外颗粒的转角出现明显差异,其中剪切带内部的颗粒转角随轴应变的增加而显著增加,而外部颗粒的转角基本保持在30°以内,未发生明显变化。已有研究表明,纯砂土在剪切过程中颗粒运动存在类似规律,即剪切带内部颗粒转角明显大于剪切带外部[24-25],结合纯砂运动规律,分析添加粘结后含可燃冰砂土剪切带内外颗粒转角存在差异的原因,可能是由于剪切带内部的粘结发生大量断裂,砂土颗粒失去约束,因此仍能发生自由转动;而剪切带外部的粘结发生断裂较少,砂土颗粒受到粘结的约束,因此颗粒转角较小。

      图  8  不同轴应变下砂土颗粒转角(Rb=55.1%)

      Figure 8.  Rotation of particles at different axial strains (Rb=55.1%)

      为了探究离散元模型加载过程中局部孔隙率的变化情况,将模型划分为30×30个大小相等的网格,计算每个网格内砂土颗粒的孔隙率,作为整体模型中该网格位置的局部孔隙率。图9给出了围压σ3=3 MPa,颗粒粘结率Rb=55.1%条件下,加载至轴应变ε1=0%、ε1=3.5%、ε1=7%、ε1=14%时模型的局部孔隙率变化情况。由图9可以看出,随着轴应变的增加,剪切带外部的孔隙率变化不明显,而剪切带内部的孔隙率有显著增加。造成这一现象的原因可能是在加载过程中,剪切带外部的砂土颗粒受粘结约束较多,呈同一方向的定向运动,因此仍能保持原先的紧密堆积;而剪切带内部的砂土颗粒由于粘结断裂而受约束减少,部分颗粒形成粒径较大的颗粒簇,且剪切过程中颗粒发生相互错动,因此可以形成较大的孔隙,如图10所示,图中箭头表示砂土颗粒位移矢量。

      图  9  不同轴应变下砂土局部孔隙率(Rb=55.1%)

      Figure 9.  Local porosity of sand at different axial strains (Rb=55.1%)

      图  10  剪切带内外颗粒位移对比(ε1=14%)

      Figure 10.  Particle displacement inside and outside shear band (ε1=14%)

      在模型划分网格的基础上,在剪切带内外分别取多个大小为3×3网格的测量单元,计算每个测量单元的孔隙率,并分别求出剪切带内部和外部测量单元孔隙率的平均值,用以探究不同颗粒粘结率Rb条件下模型剪切带内外孔隙率变化情况。取围压σ3=3 MPa,颗粒粘结率分别为Rb=0%、Rb=25.7%、Rb=40.7%和Rb=55.1%,加载至轴应变ε1=14%进行分析,结果如图11所示。由图11可知,对于纯砂土(Rb=0%)和含可燃冰砂土来说,加载至该轴应变条件下,剪切带内部孔隙率均明显高于剪切带外部;同时可以看出,加载至该轴应变时,模型颗粒粘结率Rb越高,剪切带内部孔隙率越大,而剪切带外部孔隙率变化不大。这一现象可以较好地印证前文“可燃冰饱和度SMH越高,试样的剪胀性越大”的试验结论,并进一步说明剪切带内部的孔隙率变化是造成试样体胀的主要因素。

      图  11  ε1=14%时不同颗粒粘结率Rb砂土剪切带内外孔隙率

      Figure 11.  Porosity of sands inside and outside shear band with different bonding ratios Rb at ε1=14%

    • 本文通过引入颗粒粘结率Rb这一指标来反映可燃冰饱和度SMH变化带来的影响,利用砂土颗粒间添加接触粘结的方式建立了胶结型含可燃冰砂土的离散元模型,并对该模型进行了排水双轴剪切试验的模拟。通过与Masui等[5]真实试验结果的对比,本文验证了离散元模型的可靠性,并在此基础上对于模拟结果进行了宏观和细观两个层次的分析。分析结果表明:

      (1) 含胶结型可燃冰砂土的抗剪强度、刚度、应变软化程度、以及剪胀性均随可燃冰饱和度SMH的增大而显著增大;

      (2) 可燃冰的存在对于砂土抗剪强度的影响主要体现在对粘聚力c的影响上,可燃冰饱和度SMH的增大可导致粘聚力c的增加,但内摩擦角φ的大小几乎不随SMH改变而改变;

      (3) 加载过程中,模型粘结断裂率Rbk随轴应变的增大而增大,其变化过程大致可以分为3个阶段,且与宏观上偏应力的变化相互对应,密切相关;

      (4) 在剪切过程中,胶结型含可燃冰砂土在剪切带内外的变化具有明显差异。其中剪切带内部胶结物发生大量破坏,砂土颗粒转动和局部孔隙率明显增加;剪切带外部的胶结物发生破坏较少,砂土颗粒几乎不发生转动。

      (5) 剪切带内部的孔隙率变化是造成试样体胀的主要因素,加载至相同轴应变时,可燃冰饱和度SMH越高,剪切带内部孔隙率越大。

参考文献 (25)

目录

    /

    返回文章
    返回