On the linear elastic responses of the 2D bonded discrete element model     

Gao-Feng Zhao  Qiuyue Yin  Adrian R. Russell  Yingchun Li  Wei Wu  Qin Li 《国际地质力学数值与分析法杂志》2019,43(1):166-182

The bonded discrete element model (DEM) is a numerical tool that is becoming widely used when studying fracturing, fragmentation, and failure of solids in various disciplines. However, its abilities to solve elastic problems are usually overlooked. In this work, the main features of the 2D bonded DEM which influence Poisson's ratio and Young's modulus, and accuracy when solving elastic boundary value problems, are investigated. Outputs of numerical simulations using the 2D bonded DEM, the finite element method, a hyper elasticity analysis, and the distinct lattice spring model (DLSM) are compared in the investigation. It is shown that a shear interaction (local) factor and a geometric (global) factor are two essential elements for the 2D bonded DEM to reproduce a full range of Poisson's ratios. It is also found that the 2D bonded DEM might be unable to reproduce the correct displacements for elastic boundary value problems when the represented Poisson's ratio is close to 0.5 or the long-range interaction is considered. In addition, an analytical relationship between the shear stiffness ratio and the Poisson's ratio, derived from a hyper elasticity analysis and applicable to discontinuum-based models, provides good agreement with outputs from the 2D bonded DEM and DLSM. Finally, it is shown that the selection of elastic parameters used the 2D bonded DEM has a significant effect on fracturing and fragment patterns of solids.  相似文献   

Improved algorithms for generalized finite element simulations of three-dimensional hydraulic fracture propagation     

Nathan Shauer  Carlos Armando Duarte 《国际地质力学数值与分析法杂志》2019,43(18):2707-2742

This paper reports improvements to algorithms for the simulation of 3-D hydraulic fracturing with the Generalized Finite Element Method (GFEM). Three optimizations are presented and analyzed. First, an improved initial guess based on solving a 3-D elastic problem with the pressure from the previous step is shown to decrease the number of Newton iterations and increase robustness. Second, an improved methodology to find the time step that leads to fracture propagation is proposed and shown to decrease significantly the number of iterations. Third, reduced computational cost is observed by properly recycling the linear part of the coupled stiffness matrix. Two representative examples are used to analyze these improvements. Additionally, a methodology to include the leak-off term is presented and verified against asymptotic analytical solutions. Conservation of mass is shown to be well satisfied in all examples.  相似文献   

格尔木河水质现状及生态风险评估          下载免费PDF全文

李雷明  武君  徐娟  温小虎  闵秀云  郭效瑛  柴小丽  李廷伟 《盐湖研究》2019,27(1):89-98

为了解经济迅速发展背景下格尔木河流域水体中痕量元素分布、水环境质量现状和生态风险,对格尔木河9个采样点水样21种痕量元素的分布及相关性进行了分析,以Cu、Zn、Hg、Cd、Cr和Pb为目标重金属,评估了格尔木河水质现状及生态风险。结果表明:重金属Pb浓度在所有水样中均低于检出限,其它痕量元素浓度在0.001(Cd)～6297.013(Sr)μg·L~(-1)范围变化。根据相关性分析,推测格尔木河水体中Li、Sc、Ti、V、Cu、Ge、Rb、Sr、Mo、Cd来源可能相同。采样点水样6种重金属浓度均达到地表水环境质量I类的标准。采用5种方法评价了格尔木河水质现状,研究区水质处于无污染的理想状态。格尔木河重金属引发的潜在生态风险处于低风险水平(生态风险指数变化范围为0.35～0.68)。  相似文献   

橡胶砂弹性动力学参数的弯曲-伸缩元试验研究     

姚玉文  刘方成  补国斌  景立平 《岩土力学》2020,41(7):2369-2379

橡胶砂作为一种廉价环保的土工材料，应用前景十分广阔，关于其力学特性的研究具有重要意义。采用弯曲–伸缩元法对橡胶砂进行了剪切波和压缩波联合测试，分别采用离散频率扫描法和初达波法确定了剪切波和压缩波的传播时间，进而得到其弹性动力学参数：初始剪切模量、侧限模量和泊松比，分析了橡胶含量和围压对橡胶砂弹性动力学参数的影响。结果表明：在同样围压下，随着橡胶含量的增加，橡胶砂初始剪切模量和侧限模量逐渐减小，初始泊松比逐渐增大；在同样橡胶含量下，随着围压的增大，橡胶砂初始剪切模量和侧限模量逐渐增大，初始泊松比逐渐减小。最后，在此基础上进行了两种因素耦合效应的分析以及相关力学机制的探讨。  相似文献   

黄陵断穹核部鳞片石墨矿地球化学特征与成因研究     

程林 《地质与勘探》2020,56(4):745-758

黄陵断穹核部石墨矿是我国鳞片(晶质)石墨矿的主要产区,对该区石墨矿的研究具有重要理论和实际意义。该区已有石墨矿研究成果多侧重于矿床地质方面,对矿床成因分析缺乏地球化学支撑。为此,本文从地球化学层面查明石墨矿原岩类型,碳质富集岩相、变质作用和晶质石墨的形成,为矿床成因和成矿作用提供地球化学依据。通过系统采集区内主要石墨矿样品,进行主量元素、微量元素、稀土元素和同位素分析,根据矿床地球化学原理和分析方法,阐明石墨矿的物质来源,成矿元素集散行为和矿床形成机制。根据主量元素含量和组合,推断石墨的原岩为碳质粉砂质泥岩。微量元素的双重性指示石墨矿的物质有两个来源:中太古代野马洞岩组和中太古代东冲河片麻杂岩。独特的稀土元素配分型式说明物源转化为原岩时存在着强烈的地球化学反应。石墨矿的碳同位素组成支持碳质来自有机物的观点。  相似文献   

Granular element method for three‐dimensional discrete element calculations     

Keng‐Wit Lim  José E. Andrade 《国际地质力学数值与分析法杂志》2014,38(2):167-188

This paper endows the recently‐proposed granular element method (GEM) with the ability to perform 3D discrete element calculations. By using non‐uniform rational B‐Splines to accurately represent complex grain geometries, we proposed an alternative approach to clustering‐based and polyhedra‐based discrete element methods whereby the need for complicated and ad hoc approaches to construct 3D grain geometries is entirely bypassed. We demonstrate the ability of GEM in capturing arbitrary‐shaped 3D grains with great ease, flexibility, and without excessive geometric information. Furthermore, the applicability of GEM is enhanced by its tight integration with existing non‐uniform rational B‐Splines modeling tools and ability to provide a seamless transition from binary images of real grain shapes (e.g., from 3D X‐ray CT) to modeling and discrete mechanics computations.Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

A stress and displacement discontinuity element method for elastic transversely anisotropic rock     

A. Bobet  V. García Marín 《国际地质力学数值与分析法杂志》2014,38(18):1898-1922

The paper presents closed‐form solutions for stress and displacement influence functions for stress discontinuity (SD) and displacement discontinuity (DD) elements, for a two‐dimensional plane‐strain elastic, transversely anisotropic medium. The solutions for SD elements are based on Kelvin's problem and for DD elements on the concept of dipoles. Stress and displacement influence functions are derived for the following elements: constant SD, linear SD, constant DD, linear DD, square root DD, parabolic DD, constant DD surface, and linear DD surface elements. The formulations are incorporated into FROCK, a hybridized boundary element method code, and are validated by providing comparisons between the results from FROCK and the finite element code ABAQUS. A limited parametric analysis shows the effects of slight anisotropy on the stress field around the tip of a crack and of the orientation of the crack with respect to the axes of elastic symmetry. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Simulation of drilling‐induced compaction bands using discrete element method     

H. Rahmati  A. Nouri  D. Chan  H. Vaziri 《国际地质力学数值与分析法杂志》2014,38(1):37-50

Rock failure is observed around boreholes often with certain types of failure zones, which are called breakouts. Laboratory‐scale drilling tests in some high‐porosity quartz‐rich sandstone have shown breakouts in the form of narrow localized compacted zones in the minimum horizontal stress direction. They are called fracture‐like breakouts. Such compaction bands may affect hydrocarbon extraction by forming barriers that inhibit fluid flow and may also be a source of sand production. This paper presents the results of numerical simulations of borehole breakouts using 3D discrete element method to investigate the mechanism of the fracture‐like breakouts and to identify the role of far‐field stresses on the breakout dimensions. The numerical tool was first verified against analytical solutions. It was then utilized to investigate the failure mechanism and breakout geometry for drilled cubic rock samples of Castlegate sandstone subjected to different pre‐existing far‐field stresses. Results show that failure occurs in the zones of the highest concentration of tangential stress around the borehole. It is concluded that fracture‐like breakout develops as a result of a nondilatant failure mechanism consisting of localized grain debonding and repacking and grain crushing that lead to the formation of a compaction band in the minimum horizontal stress direction. In addition, it is found that the length of fracture‐like breakouts depends on both the mean stress and stress anisotropy. However, the width of the breakout is not significantly changed by the far‐field stresses. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

Modeling deepwater seabed steady-state thermal fields around buried pipeline including trenching and backfill effects     

《Computers and Geotechnics》2014

Deepwater pipelines are designed to transport mixtures of oil and gas, and their associated impurities at wellhead temperatures that can be in excess of 149 °C (∼300 °F or 422 K) while the external temperature maybe in the range of 5 °C (∼41 °F or 278 K). Depending on the circumstances these pipelines may be buried for physical protection or for additional thermal insulation using robotic trenching equipment. This results in a complex cut and backfill geometry in the seafloor in addition to altering the thermal properties of the backfill. A two-dimensional boundary element model was developed specifically to address to investigate the local steady-state thermal field in the near field of the pipeline. The model allows one to account for the complex geometries in the near field associated with this burial technique, site-specific multi-layered soil conditions and the seawater adjacent to the seafloor. A parametric study was preformed to evaluate effects of the thermal power loss, burial depth, pipe diameter and soil thermal conductivity on the thermal field in the near field of a buried pipeline. The numerical examples illustrate the influence of the backfill thermal property on the temperature at the pipe wall, that the pipe diameter controls the required output thermal power needed to maintain the desired pipe wall temperature, and the importance of pipeline burial depth on seabed temperature distribution above the pipeline.  相似文献   

Critical Top Tension for Static Equilibrium Configuration of A Steel Catenary Riser     

Chainarong ATHISAKUL  Karun KLAYCHAM  Somchai CHUCHEEPSAKUL 《中国海洋工程》2014,(6):829-842

This paper aims to present the critical top tension for static equilibrium configurations of a steel catenary riser（SCR） by using the finite element method. The critical top tension is the minimum top tension that can maintain the equilibrium of the SCR. If the top tension is smaller than the critical value, the equilibrium of the SCR does not exist. If the top tension is larger than the critical value, there are two possible equilibrium configurations. These two configurations exhibit the nonlinear large displacement. The configuration with the smaller displacement is stable, while the one with larger displacement is unstable. The numerical results show that the increases in the riser＇s vertical distances, horizontal offsets, riser＇s weights, internal flow velocities, and current velocities increase the critical top tensions of the SCR. In addition, the parametric studies are also performed in order to investigate the limit states for the analysis and design of the SCR.  相似文献