Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis     

Ramin Pakzad  Shanyong Wang  Scott W. Sloan 《国际地质力学数值与分析法杂志》2020,44(10):1472-1494

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.  相似文献   

Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law          下载免费PDF全文

Ahmad Hajjar  Luc Scholtès  Constantin Oltéan  Michel Antoine Buès 《水文研究》2018,32(16):2510-2525

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.  相似文献   

Implementation of a coupled plastic damage distinct lattice spring model for dynamic crack propagation in geomaterials          下载免费PDF全文

Chao Jiang  Gao‐Feng Zhao 《国际地质力学数值与分析法杂志》2018,42(4):674-693

This paper studies dynamic crack propagation by employing the distinct lattice spring model (DLSM) and 3‐dimensional (3D) printing technique. A damage‐plasticity model was developed and implemented in a 2D DLSM. Applicability of the damage‐plasticity DLSM was verified against analytical elastic solutions and experimental results for crack propagation. As a physical analogy, dynamic fracturing tests were conducted on 3D printed specimens using the split Hopkinson pressure bar. The dynamic stress intensity factors were recorded, and crack paths were captured by a high‐speed camera. A parametric study was conducted to find the influences of the parameters on cracking behaviors, including initial and peak fracture toughness, crack speed, and crack patterns. Finally, selection of parameters for the damage‐plasticity model was determined through the comparison of numerical predictions and the experimentally observed cracking features.  相似文献   

Analytical study of ground motion caused by seismic wave propagation across faulted rock masses          下载免费PDF全文

J.C. Li  N.N. Li  S.B. Chai  H.B. Li 《国际地质力学数值与分析法杂志》2018,42(1):95-109

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.  相似文献   

自然和自然贡献情景模型研究综述     

岳天祥  范泽孟  史文娇  赵娜 《地球信息科学学报》2020,22(4):743-750

本文回顾了自然和自然贡献情景模型发展的背景、历史和内容,概括总结了自然和自然贡献情景模型的发展进程以及联合国生物多样性与生态系统服务政府间科学—政策平台(IPBES)情景模型的概念框架,讨论了自然和自然贡献情景模型存在的问题和发展方向。为了在全球层面解决现有综合集成模型存在的问题,根据地球表层建模基本定理和生态环境曲面建模基本定理,提出了具有中国原创特点的自然与自然贡献情景模型概念框架。  相似文献   

Anisotropies in the redshift-space correlations of galaxy groups and clusters – I. Simulated catalogues     

N. D. Padilla  D. G. Lambas † 《Monthly notices of the Royal Astronomical Society》2003,342(2):519-531

  相似文献   

景观生态分类与制图浅议   总被引：12，自引：0，他引：12  

程维明 《地球信息科学学报》2002,4(2):61-65

本文在查阅分析大量文献和前人研究的基础上 ,对目前景观生态分类和景观制图作了详细的对比分析 ,认为景观分类需要结合实际区域现状 ,采用逐级分类的方法 ;同时利用 ETM遥感影像为数据源 ,以天山北麓为示范区.研制其土地利用土地覆盖变化的景观类型图。  相似文献   

Cosmological constraints from the local X-ray luminosity function of the most X-ray-luminous galaxy clusters     

S. W. Allen  R. W. Schmidt  A. C. Fabian  H. Ebeling 《Monthly notices of the Royal Astronomical Society》2003,342(1):287-298

  相似文献   

平衡剖面的制作流程及其地质意义   总被引：9，自引：0，他引：9  

王运所  刘亚洲  张孝义  张剑峰  尹哲 《地球科学与环境学报》2003,25(1):28-32

平衡剖面技术是地质思维和计算机技术的结晶,使对断层构造的研究提高到定量阶段,其依据是在垂直构造走向的剖面上,地层长度和面积(2D)或体积(3D)是均衡的。在此原理基础上利用数学手段对盆地的构造发育史进行正演和反演模拟,直观地再现地下构造的原始几何形态,迅速提供地震剖面的构造解释方案,并对解释结果进行检验(不平衡的剖面其解释一般有问题),为深刻认识构造发育史、分析油气运移及聚集规律提供依据,提高了工作效率。其结果也为盆地模拟、油藏模拟、定量计算构造伸缩量等地质研究打下了坚实的基础[1]。  相似文献   

Laboratory kinetics of C₂H radical reactions with ethane, propane, and n-butane at T = 96-296 K: implications for Titan     

James E. Murphy  Stephen R. Leone 《Icarus》2003,163(1):175-181

The kinetics of the reactions of C₂H radical with ethane (k₁), propane (k₂), and n-butane (k₃) are studied over the temperature range of T = 96-296 K with a pulsed Laval nozzle apparatus that utilizes a pulsed laser photolysis-chemiluminescence technique. The C₂H decay profiles in the presence of both the alkane reactant and O₂ are monitored by the CH(A²Δ) chemiluminescence tracer method. The results, together with available literature data, yield the following Arrhenius expressions: k₁(T) = (0.51 ± 0.06) × 10⁻¹⁰ exp[(−76 ± 30)K/T] cm³ molecule⁻¹ s⁻¹ (T = 96-800 K), k₂(T) = (0.98 ± 0.32) × 10⁻¹⁰exp[(−71 ± 60)K/T] cm³ molecule⁻¹ s⁻¹ (T = 96-361 K), and k₃(T) = (1.23 ± 0.26) × 10⁻¹⁰ cm³ molecule⁻¹ s⁻¹ (T = 96-297 K). At T = 296 K, k₁ is measured as a function of total pressure and has little or no pressure dependence. The results from this work support a direct hydrogen abstraction mechanism for the title reactions. Implications to the atmospheric chemistry of Titan are discussed.  相似文献