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.  相似文献   

Bacterial productivity in the Prydz Bay and its adjacent waters,Antarctic     

邱雨生  黄奕普  陈敏  刘广山 《极地研究(英文版)》2004,15(1):14-27

1　IntroductionBacteriaandtheiractivitiesplayanimportantroleintheelementalbiogeochemicalcyclesandenergytransformingintheocean (Zhenetal.1 997) .DortchandPackard(1 989) proposedthatfoodwebsintheeutrophicwatersaredominatedbythebiomassofprimaryproducerswhilefoodwebsintheoligotrophicwatersaredominatedbythebiomassofmicrobes.Heterotrophicbacteriahadbeenshowntoplayanimportantroleinthedecompositionoflarge ,rapidlysinkingorganicparticleswithinandbelowtheeuphot iczone ,andfurthertoaffecttheelementaldyn…  相似文献   

Kuiper带小天体动力演化     

周济林  周礼勇  孙义燧 《天文研究与技术》2003,(1):76-81

Kuiper带是指太阳系内位于离太阳 30～ 50AU一个区域。 1 992年该区域陆续发现了一群半径在几十到几百公里的小天体。这些小天体在Kuiper带的分布是极其不均匀的。Kuiper带小天体的发现对人们认识太阳系的形成与演化有重要的意义。本文回顾了近年来国际上在Kuiper带小天体动力演化方面的研究 ,着重分析了目前国际上几种用以解释其非均匀分布的动力学机制 ,并提出目前该领域的一些尚未解决的问题。  相似文献   

Virtual Huanghe River System： Framework and Technology   总被引：2，自引：0，他引：2  

LU Heli LIU Guifang SUN Jiulin 《中国地理科学(英文版)》2006,16(3):255-259

1 Introduction Huanghe (Yellow) River basin is located in 32°–42°N, 96°–119°E. The area of the catchment is more than 752,000km2. The river is 5464km long with a drop in elevation of 4830m. Among the whole area, the moun- tainous and stone area accounts for 29%, loess and hills area 46%, sandy area 11% and plain area 14%, respec- tively. Different natural landscapes exist in this area. The Huanghe River flows through the Loess Plateau, where the soil is eroded seriously (Wang, 2002;…  相似文献   

四川石棉独立碲矿床蚀变岩的地质地球化学研究     

银剑钊  陈毓川  周剑雄  杨百川  郭成远 《地球科学与环境学报》1996,(4)

四川石棉大水沟独立碲矿床堪称世界首例，其矿体两侧的围岩蚀变较窄，且蚀变带与矿脉及未蚀变岩呈清楚而截然的接触关系。主要围岩蚀变为白云石化、白（绢）云母化、云英岩化及电气石化等。这些蚀变作用分别发生于１７７．７～１６５．１Ｍａ和９１．７１～８０．１９Ｍａ。蚀变作用过程中ＳｉＯ２大量析出，Ａｌ２Ｏ３、ＴｉＯ２、ＣａＯ、ＭｇＯ、Ｋ２Ｏ、Ｈ２Ｏ、ＣＯ２、Ｆｅ２Ｏ３、ＭｎＯ则有不同程度地带入。与此同时，ＲＥＥ在此过程中包集体迁出，Ｔｅ、Ｂｉ、ＡＳ等元素的含量则决定于蚀变作用及其强度，而与原岩无关。  相似文献   

有限元法与伪谱法混合求解弹性波动方程   总被引：6，自引：0，他引：6  

马德堂  朱光明 《地球科学与环境学报》2004,26(1):61-64

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.  相似文献   

Investigating the mechanisms leading to the deglaciation of past continental northern hemisphere ice sheets with the CLIMBER–GREMLINS coupled model     

Sylvie Charbit  Masa Kageyama  Didier Roche  Catherine Ritz  Gilles Ramstein 《Global and Planetary Change》2005,48(4):1887

A coupling procedure between a climate model of intermediate complexity (CLIMBER-2.3) and a 3-dimensional thermo-mechanical ice-sheet model (GREMLINS) has been elaborated. The resulting coupled model describes the evolution of atmosphere, ocean, biosphere, cryosphere and their mutual interactions. It is used to perform several simulations of the Last Deglaciation period to identify the physical mechanisms at the origin of the deglaciation process. Our baseline experiment, forced by insolation and atmospheric CO₂, produces almost complete deglaciation of past northern hemisphere continental ice sheets, although ice remains over the Cordilleran region at the end of the simulation and also in Alaska and Eastern Siberia. Results clearly demonstrate that, in this study, the melting of the North American ice sheet is critically dependent on the deglaciation of Fennoscandia through processes involving switches of the thermohaline circulation from a glacial mode to a modern one and associated warming of the northern hemisphere. A set of sensitivity experiments has been carried out to test the relative importance of both forcing factors and internal processes in the deglaciation mechanism. It appears that the deglaciation is primarily driven by insolation. However, the atmospheric CO₂ modulates the timing of the melting of the Fennoscandian ice sheet, and results relative to Laurentide illustrate the existence of threshold CO₂ values, that can be translated in terms of critical temperature, below which the deglaciation is impeded. Finally, we show that the beginning of the deglaciation process of the Laurentide ice sheet may be influenced by the time at which the shift of the thermohaline circulation from one mode to the other occurs.  相似文献   

X-ray synchrotron emission from the oblique shock in the jet of the powerful radio galaxy 3C 346     

D. M. Worrall  M. Birkinshaw 《Monthly notices of the Royal Astronomical Society》2005,360(3):926-934

We report the first detection, with Chandra , of X-ray emission from the jet of the powerful narrow-line radio galaxy 3C 346. X-rays are detected from the bright radio and optical knot at which the jet apparently bends by approximately 70°. The Chandra observation also reveals a bright galaxy-scale atmosphere within the previously known cluster and provides a good X-ray spectrum for the bright core of 3C 346. The X-ray emission from the knot is synchrotron radiation, as seen in lower-power sources. In common with these sources, there is evidence of morphological differences between the radio/optical and X-ray structures, and the spectrum is inconsistent with a one-component continuous-injection model. We suggest that the X-ray-bright knot is associated with a strong oblique shock in a moderately relativistic, light jet, at ∼ 20° to the line of sight, and that this shock is caused by the jet interacting with the wake in the cluster medium behind the companion galaxy of 3C 346. The general jet curvature can result from pressure gradients in the cluster atmosphere.  相似文献