CN110598323B - 一种渗透破坏离散元模拟方法 - Google Patents
一种渗透破坏离散元模拟方法 Download PDFInfo
- Publication number
- CN110598323B CN110598323B CN201910864260.6A CN201910864260A CN110598323B CN 110598323 B CN110598323 B CN 110598323B CN 201910864260 A CN201910864260 A CN 201910864260A CN 110598323 B CN110598323 B CN 110598323B
- Authority
- CN
- China
- Prior art keywords
- simulation
- fluid
- force
- particles
- gradient
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Abstract
本发明提供了一种渗透破坏离散元模拟方法,包括:在无流体的固体颗粒模拟中,施加一种与流场水力梯度相关的力,即梯度力,用以替代流体对于颗粒的作用力,以对渗透破坏现象进行简化模拟。本发明能够在实现流固耦合效果的基础上,省略流固耦合计算过程,提高模拟计算速度。
Description
技术领域
本发明属于渗透破坏模拟技术领域,具体涉及一种渗透破坏离散元模拟方法。
背景技术
在当前的离散元软件中,多采用流固耦合方法模拟渗透破坏时的颗粒运动,该方法虽然接近实际的渗透破坏物理过程,但其计算量大,计算效率较低,模拟计算用时较长,因此需寻求合理的简化方法对该模拟进行简化,以提高渗透破坏模拟计算的效率。
发明内容
本发明的目的是提供一种渗透破坏离散元模拟方法,以解决上述技术问题。
本发明提供了一种渗透破坏离散元模拟方法,包括:
在无流体的固体颗粒模拟中,施加一种与流场水力梯度相关的力,即梯度力,用以替代流体对于颗粒的作用力,以对渗透破坏现象进行简化模拟。
进一步地,所述渗透破坏离散元模拟方法具体包括如下步骤:
步骤1,确定所要模拟的渗透破坏类型;
步骤2,确定发生渗透破坏土体的细颗粒粒径、孔隙率,以及所要模拟的流场水力梯度;
步骤3,计算阻力影响比,确定其误差是否在可接受范围之内;
步骤4,采用给定加速度场模拟流场梯度力,进行简化的渗透破坏模拟。
进一步地,步骤3中基于下述公式计算阻力影响比θ:
进一步地,步骤4中基于下述公式施加沿流线方向的加速度场,从而实现施加流场梯度力:
式中:αi为待施加的沿流线方向的加速度;γw为流体容重;γs为颗粒容重;i为待施加的水力梯度。
与现有技术相比本发明的有益效果是:
能够在实现流固耦合效果的基础上,省略流固耦合计算过程,提高模拟计算速度。
附图说明
图1是本发明一种渗透破坏离散元模拟方法的流程图;
图2为本发明一实施例中孔隙率参数为疏松时θ~r的敏感度影响曲线;
图3为本发明一实施例中孔隙率参数为普通时θ~r的敏感度影响曲线;
图4为本发明一实施例中孔隙率参数为密实时θ~r的敏感度影响曲线;
图5为本发明一实施例中细颗粒粒径为疏松时θ~i的敏感度影响曲线;
图6为本发明一实施例中细颗粒粒径为普通时θ~i的敏感度影响曲线;
图7为本发明一实施例中细颗粒粒径为密实时θ~i的敏感度影响曲线。
具体实施方式
下面结合附图所示的各实施方式对本发明进行详细说明,但应当说明的是,这些实施方式并非对本发明的限制,本领域普通技术人员根据这些实施方式所作的功能、方法、或者结构上的等效变换或替代,均属于本发明的保护范围之内。
在岩土工程常见的渗透破坏模拟中,人们主要聚焦于颗粒在流体作用下的受力,进而分析其可能发生渗透破坏的状态与特征。针对渗透破坏的流固耦合模拟,本实施例提出了一种简化的思路,并对其进行了分析验证。
当前离散元模拟中,多采用如下控制方程表征颗粒在流场中与流体耦合的受力状态:
对于该公式,本实施例将其转化为用各项力表征的公式,即为:
因此,相对于传统的纯固体离散元研究,流固耦合模拟的重点在于准确地模拟出流体对于颗粒的作用力的效果。在渗透破坏过程中,流体对任一颗粒作用力主要为沿流线方向的渗流场力受流固界面因速度差异产生的拖曳力并且在孔隙介质中,拖曳力相对较小。
根据本方案的模拟对象。依照图1所示步骤进行该模拟方法的适用性判定,当达到模拟对象的误差允许范围之内时,即可采用该简化模拟方法对渗透破坏现象进行简化模拟。
在上文的推导与分析中,利用离散元开展流固耦合模拟的目的即是正确反映流体对于颗粒的作用力。
以管涌现象为例,定义阻力影响比θ,有:
由上,θ表征了颗粒所受拖曳力在其受到流体作用力中的占比,即阻力对于流体对颗粒作用力的影响程度,即使用该简化方案所带来的相对误差大小。
针对该分析中的部分参数,其取值表如下:
为了更加简化地分析各项不确定参数与θ的响应关系,此处引入如下两个关系:
(1)Darcy公式
(2)Kozeny-Carman公式
其中,k为渗透系数,S0为颗粒单位体积的比表面积,∈为土体孔隙率。
对于离散元模拟中形状概化为标准球形的颗粒,可通过如下公式计算其单位体积的比表面积:
所以原式变为:
根据上文分析,最终联立得到阻力影响比θ的最终公式:
考虑到该方案是对传统流固耦合模拟中复杂情况的一些简化,因此需要通过对其误差进行分析来验证其有效性。
可以看出,在结合前文给出的关系,阻力影响比θ主要受如下三个参数的影响:粒径r、水力梯度i、孔隙率∈。
对于孔隙率∈,土体试样在受到流体冲蚀而发生管涌现象时,其孔隙率的不同将会导致不同的管涌破坏后果,因此针对不同孔隙率的土体试样进行研究非常有必要。根据传统渗透破坏研究的试样调查,我们在此次误差分析中分别选取孔隙率参数为0.1(密实),0.35(普通),0.5(疏松)三种状态进行分析。
1、粒径r
针对管涌现象,其为在渗流场的作用下,一定级配的无黏性土中的细小颗粒通过较大颗粒形成的孔隙发生移动而产生破坏。本实施例中着重分析无黏性土中的细小颗粒受力情况。此处主要考虑粒径在0.06-2mm的土体颗粒。
根据图2、图3、图4可以看出,在较大范围的水力梯度作用下,三种密实程度(分别为密实、正常、疏松)的试样,采用该简化方法所带来的误差,即阻力影响比θ随细颗粒粒径的变化,均保持在0.05%~0.4%之间,这符合渗透破坏现象模拟的误差要求,因此该方法具有较强的有效性。
2、水力梯度i
根据目前的实验以及数值模拟采用数据,对于具体产生管涌现象水力梯度没有较为明确的范围规定。本实施例采取一个较大的水力梯度分析范围,即i=1~20进行参数的敏感度分析。
根据图5、图6、图7可以看出,在给定细颗粒粒径的范围内,三种密实程度(分别为密实、正常、疏松)的试样,采用该简化方法所带来的误差,即阻力影响比θ随流场水力梯度的变化,亦保持在0.05%~0.4%之间,这符合渗透破坏现象模拟的误差要求,因此该方法具有较强的有效性。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。
Claims (1)
1.一种渗透破坏离散元模拟方法,其特征在于,包括:
在无流体的固体颗粒模拟中,施加一种与流场水力梯度相关的力,即梯度力,用以替代流体对于颗粒的作用力,以对渗透破坏现象进行简化模拟;
所述渗透破坏离散元模拟方法具体包括如下步骤:
步骤1,确定所要模拟的渗透破坏类型;
步骤2,确定发生渗透破坏土体的细颗粒粒径、孔隙率,以及所要模拟的流场水力梯度;
步骤3,计算阻力影响比,确定其误差是否在可接受范围之内;
步骤4,采用给定加速度场模拟流场梯度力,进行简化的渗透破坏模拟;
步骤3中基于下述公式计算阻力影响比θ:
步骤4中基于下述公式施加沿流线方向的加速度场,从而实现施加流场梯度力:
式中:αi为待施加的沿流线方向的加速度;γw为流体容重;γs为颗粒容重;i为待施加的水力梯度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910864260.6A CN110598323B (zh) | 2019-09-12 | 2019-09-12 | 一种渗透破坏离散元模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910864260.6A CN110598323B (zh) | 2019-09-12 | 2019-09-12 | 一种渗透破坏离散元模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110598323A CN110598323A (zh) | 2019-12-20 |
CN110598323B true CN110598323B (zh) | 2021-05-11 |
Family
ID=68859506
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910864260.6A Active CN110598323B (zh) | 2019-09-12 | 2019-09-12 | 一种渗透破坏离散元模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110598323B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113075106B (zh) * | 2021-03-26 | 2022-11-11 | 桂林理工大学 | 一种砾石土侵蚀阶段的水力参数及涌水量的确定方法 |
CN113237722B (zh) * | 2021-05-26 | 2022-04-22 | 中国水利水电科学研究院 | 一种土体内部侵蚀变形响应模拟方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7413570B2 (en) * | 2002-08-21 | 2008-08-19 | Kci Licensing, Inc. | Medical closure screen installation systems and methods |
CN101551830A (zh) * | 2009-03-04 | 2009-10-07 | 同济大学 | 基于流固两相流的混凝土搅拌车搅拌筒叶片设计方法 |
JP2011219224A (ja) * | 2010-04-09 | 2011-11-04 | Murata Mfg Co Ltd | 気体搬送装置 |
CN103077556A (zh) * | 2013-02-04 | 2013-05-01 | 重庆大学 | 油井出砂的三维数值模型设计 |
CN109297880A (zh) * | 2018-10-11 | 2019-02-01 | 中国电建集团华东勘测设计研究院有限公司 | 深埋水工隧洞渗透梯度模拟试验系统和试验方法 |
CN109887083A (zh) * | 2019-01-29 | 2019-06-14 | 中国石油集团测井有限公司西南分公司 | 一种裂缝-孔隙双重介质耦合渗透率模型的建立方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103970989B (zh) * | 2014-04-15 | 2017-08-25 | 昆明理工大学 | 一种基于流固界面一致条件的浸入边界流场计算方法 |
US9881108B2 (en) * | 2015-05-29 | 2018-01-30 | One Energy Enterprises Llc | Method of evaluation wind flow based on conservation of momentum and variation in terrain |
CN109902376B (zh) * | 2019-02-25 | 2021-01-15 | 北京理工大学 | 一种基于连续介质力学的流固耦合高精度数值模拟方法 |
-
2019
- 2019-09-12 CN CN201910864260.6A patent/CN110598323B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7413570B2 (en) * | 2002-08-21 | 2008-08-19 | Kci Licensing, Inc. | Medical closure screen installation systems and methods |
CN101551830A (zh) * | 2009-03-04 | 2009-10-07 | 同济大学 | 基于流固两相流的混凝土搅拌车搅拌筒叶片设计方法 |
JP2011219224A (ja) * | 2010-04-09 | 2011-11-04 | Murata Mfg Co Ltd | 気体搬送装置 |
CN103077556A (zh) * | 2013-02-04 | 2013-05-01 | 重庆大学 | 油井出砂的三维数值模型设计 |
CN109297880A (zh) * | 2018-10-11 | 2019-02-01 | 中国电建集团华东勘测设计研究院有限公司 | 深埋水工隧洞渗透梯度模拟试验系统和试验方法 |
CN109887083A (zh) * | 2019-01-29 | 2019-06-14 | 中国石油集团测井有限公司西南分公司 | 一种裂缝-孔隙双重介质耦合渗透率模型的建立方法 |
Non-Patent Citations (3)
Title |
---|
WeiChao Li等.Discrete element modeling of the Hongshiyan landslide triggered by the 2014 Ms 6.5 Ludian earthquake in Yunnan, China.《Environmental Earth Sciences (2019)》.2019, * |
散粒土管涌临界水力梯度的研究;王霜等;《水电能源科学》;20180930;第36卷(第9期);第114-117页 * |
考虑流固耦合效应的饱和砂土渗透破坏数值模拟;刘先珊等;《兰州大学学报(自然科学版)》;20131031;第49卷(第5期);第633-638、643页的第1-2节 * |
Also Published As
Publication number | Publication date |
---|---|
CN110598323A (zh) | 2019-12-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Vincens et al. | Approaches to determine the constriction size distribution for understanding filtration phenomena in granular materials | |
CN110598323B (zh) | 一种渗透破坏离散元模拟方法 | |
Alonso et al. | Numerical analysis of rapid drawdown: Applications in real cases | |
CN111709148B (zh) | 一种粘性砂土水力侵蚀破坏的离散元数值模拟方法 | |
CN106837315B (zh) | 裂缝性碳酸盐岩基质与裂缝耦合作用表征方法 | |
Poulsen et al. | Convergence of synthetic rock mass modelling and the Hoek–Brown strength criterion | |
CN112347709B (zh) | 一种基于dem-cfd耦合的渗透注浆过程模拟方法及系统 | |
Madabhushi et al. | LEAP-GWU-2015: Centrifuge and numerical modelling of slope liquefaction at the University of Cambridge | |
CN114117879A (zh) | 一种砂卵石土三轴剪切离散元模型细观参数标定方法 | |
Zhang et al. | Interpretation of cone penetration test in clay with smoothed particle finite element method | |
Moug et al. | Axisymmetric simulations of cone penetration in saturated clay | |
CN110441209A (zh) | 一种基于致密储层数字岩心计算岩石渗透率的方法 | |
Moret‐Fernández et al. | A modified multiple tension upward infiltration method to estimate the soil hydraulic properties | |
Arulanandan et al. | Numerical simulation of liquefaction-induced deformations | |
Shi et al. | Simulation of cyclic strength degradation of natural clays via bounding surface model with hybrid flow rule | |
Pells et al. | Rock mass erodibility | |
Chrisopoulos et al. | FE simulation of model tests on vibratory pile driving in saturated sand | |
Abelev et al. | Numerical investigation of dynamic free-fall penetrometers in soft cohesive marine sediments using a finite difference approach | |
CN105241911A (zh) | 基于lbm模拟低场核磁共振分析流体的方法及装置 | |
CN114611434A (zh) | 一种体现渣土改良泡沫剂渗流特性的模型构建方法及系统 | |
CN110750901A (zh) | 基于离散元模型的土体扰动范围判断方法 | |
CN110569548A (zh) | 一种各向异性岩质边坡浸润线确定方法 | |
CN112485106A (zh) | 一种控制土体状态参数的物理模型分层制备与试验方法 | |
Jiang et al. | Experimental study on the applicability of failure criteria for rockfill in three-dimensional stress conditions | |
Esmaeilzadeh et al. | Monitoring Seepage through Dam under Consolidation Using the Finite Element Method, Case Study: Karkheh Dam |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |