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

青海湖流域植被盖度时空变化研究     

高黎明  张乐乐 《地球信息科学学报》2019,21(9):1318-1329

高寒区植被变化一直是气候和生态学领域关注的热点问题。本研究基于MODIS NDVI数据计算的植被覆盖度数据和高分辨率气象数据,分析了青海湖流域2001-2017年植被覆盖度分布格局及动态变化,探讨了其对气候变化、人类活动和冻土退化的响应。结果表明：① 近十几年青海湖流域植被覆盖度整体表现为增加趋势,不同植被类型增幅存在差异性,草地增幅最大,达到6.1%/10a,其它植被类型增幅在2%~3%/10a之间;② 流域局部地区仍存在植被退化现象,研究期植被退化面积表现为先增加后减小的变化趋势。2006-2011年重度退化区集中在青海湖东岸,2011-2017年重度退化区集中在流域的西北部,这些区域是青海湖流域荒漠分布区,植被覆盖度较低,是今后生态恢复需重点关注的区域;③ 气候变化是流域植被覆盖度变化的主导因素,气候变化对青海湖流域主要植被类型覆盖度变化的贡献率为84.21%,对草原、草甸和灌丛植被覆盖度变化的贡献率分别为81.84%、87.47%和75.96%;④ 人类活动对流域主要植被类型覆盖度变化的贡献率为15.79%,对草原、草甸和灌丛植被覆盖度变化的贡献率分别为18.16%、12.53%和24.04%,环青海湖地区人类活动对植被恢复有促进效应,在青海湖流域北部部分地区人类活动的破坏力度仍大于建设力度;⑤ 冻土退化对青海湖流域草甸和灌丛植被覆盖度变化影响很小,主要影响草原植被覆盖度变化,冻土退化造成草原植被覆盖度增长速率减小了1.2%/10a。  相似文献   

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock          下载免费PDF全文

N. Zolfaghari  E. Dontsov  A. P. Bunger 《国际地质力学数值与分析法杂志》2018,42(4):587-617

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.  相似文献   

Micromechanical approach to effective viscoelastic properties of micro‐fractured geomaterials          下载免费PDF全文

Cássio B. Aguiar  Samir Maghous 《国际地质力学数值与分析法杂志》2018,42(16):2018-2046

The aim of this paper is to formulate a micromechanics‐based approach to non‐aging viscoelastic behavior of materials with randomly distributed micro‐fractures. Unlike cracks, fractures are discontinuities that are able to transfer stresses and can therefore be regarded from a mechanical viewpoint as interfaces endowed with a specific behavior under normal and shear loading. Making use of the elastic‐viscoelastic correspondence principle together with a Mori‐Tanka homogenization scheme, the effective viscoelastic behavior is assessed from properties of the material constituents and damage parameters related to density and size of fractures. It is notably shown that the homogenized behavior thus formulated can be described in most cases by means of a generalized Maxwell rheological model. For practical implementation in structural analyses, an approximate model for the isotropic homogenized fractured medium is formulated within the class of Burger models. Although the approximation is basically developed for short‐term and long‐term behaviors, numerical applications indicate that the approximate Burger model accurately reproduce the homogenized viscoelastic behavior also in the transient conditions.  相似文献   

黑龙江哈尔滨白渔泡国家湿地公园沼泽、林地和农田土壤物理、化学和生物性质差异     

丁俊男  于少鹏  史传奇 《湿地科学》2020,18(1):77-84

为了揭示黑龙江哈尔滨白渔泡国家湿地公园沼泽、林地和农田土壤物理、化学和生物性质的差异,于2018年7月25日~8月2日,在湿地公园内,在天然芦苇(Phragmites australis)沼泽、林地、旱田和水田中设置采样地,采集不同深度(0~10 cm、10~20 cm和20~30 cm)的土壤样品,测定土壤样品的物理、化学和生物指标。研究结果表明,白渔泡国家湿地公园不同采样地土壤指标存在差异;与天然芦苇沼泽土壤相比,其它采样地土壤的含水量明显偏低,土壤全氮、全磷、碱解氮和有机质含量都明显偏小,水田土壤速效磷含量偏大;天然芦苇沼泽土壤脲酶、硝酸还原酶、纤维素酶、蛋白酶和β-葡萄糖苷酶活性都高于林地和农田土壤,水田0~10 cm和10~20 cm深度土壤的硝酸还原酶活性显著高于旱田和林地;与天然芦苇沼泽土壤相比,旱田土壤小于0.25 mm的小团聚体含量偏大,而其它采样地土壤的各粒级团聚体的比例变化较小,水田土壤团聚体平均重量直径比天然芦苇沼泽和旱田土壤低。  相似文献   

Activity and kinematics of members of the TW Hydrae association     

Neill Reid 《Monthly notices of the Royal Astronomical Society》2003,342(3):837-850

We present high-resolution echelle spectroscopy of 20 stars in 16 systems catalogued as members of the TW Hydrae association, and 16 stars identified as possible new members. We have calibrated the range of coronal and chromospheric activity expected for such young stars as a function of spectral type by combining our observations with literature data for field and open cluster stars. We also compute space motions for TWA members and candidate members with proper motion measurements, using two techniques to estimate distances to stars lacking direct trigonometric parallax measurements. The mean space motion of the four TWA members with known parallaxes is  ( U , V , W : −10.0, −17.8, −4.6) km s⁻¹  . 14 of the candidates have properties inconsistent with cluster membership; the remaining two are potential new members, although further observations are required to confirm this possibility.  相似文献   

Evolutionary Processes in Low-Mass Binary Systems     

Peter P. Eggleton  Ludmila Kisseleva-Eggleton 《Astrophysics and Space Science》2005,296(1-4):327-336

We consider the evolution of certain low-mass binaries, incorporating models of (a) internal evolution, (b) tidal friction, (c) dynamo activity driven by an elementary α,Ω dynamo, (d) stellar wind driven by the activity, and (e) magnetic braking as a consequence of wind and poloidal dynamo-generated magnetic field. In some circumstances the stellar wind is found to remove mass on a nuclear timescale, as is necessary to explain some observed systems. We can hope that various uncertainties in the model may be clarified by a careful comparison of the models with such observed quantities as rotation periods. These are modified by processes (a), (b) and (e). Assuming that stellar evolution is slow, rotation rate should in some circumstances represent a balance between magnetic braking trying to slow the star down and tidal friction trying to spin it up. Preliminary attempts are promising, but indicate that some fine tuning is necessary. When there is a third body present, in an orbit which is inclined but not necessarily of short period, the eccentricity of a close binary can be strongly modified by ‘Kozai cycles’. We show that this may complicate attempts to account for spin rates of stars in close binaries.  相似文献   

The Relation between the Amplitude and the Period of Solar Cycles   总被引：3，自引：0，他引：3  

Zhan-Le Du  Hua-Ning Wang  Xiang-Tao He National Astronomical Observatories  Chinese Academy of Sciences  Beijing  Department of Astronomy  Beijing Normal University  Beijing 《中国天文和天体物理学报》2006,6(4):489-494

The maximum amplitudes of solar activity cycles are found to be well anti-correlated (r = -0.72) with the newly defined solar cycle lengths three cycles before (at lag -3) in 13-month running mean sunspot numbers during the past 190 years. This result could be used for predicting the maximum sunspot numbers. The amplitudes of Cycles 24 and 25 are estimated to be 149.5±27.6 and 144.3±27.6, respectively.  相似文献   

LU Vel (GJ 375): A M3.5Ve Flare and Double-Lined Spectroscopic Binary     

D. Montes  M. C. Gálvez  M. J. Fernández-Figueroa  I. Crespo-Chacón 《Astrophysics and Space Science》2006,304(1-4):367-369

High resolution echelle spectroscopic observations taken with the FEROS spectrograph at the 2.2 m telescope ESO confirm the binary nature of the flare M3.5V star LU Vel (GJ 375, RE J0958-462) previously reported by Christian and Mathioudakis (2002). Emission of similar intensity from both components is detected in the Balmer, Na i D₁&D₂, He i D₃, Ca ii H&K, and Ca ii IRT lines. We have determined precise radial velocities by cross correlation with radial velocity standard stars, which have allowed us to obtain for the first time the orbital solution of the system. The binary consists of two near-equal M3.5V components with an orbital period shorter than 2 days. We have analyzed the behaviour of the chromospheric activity indicators (variability and possible flares). In addition, we have determined its rotational velocity and kinematics.  相似文献