A comparative study of the hypoplasticity and the fabric‐dependent dilatant double shearing models for granular materials

In this paper, we consider the mechanical response of granular materials and compare the predictions of a hypoplastic model with that of a recently developed dilatant double shearing model which includes the effects of fabric. We implement the constitutive relations of the dilatant double shearing model and the hypoplastic model in the finite element program ABACUS/Explicit and compare their predictions in the triaxial compression and cyclic shear loading tests. Although the origins and the constitutive relations of the double shearing model and the hypoplastic model are quite different, we find that both models are capable of capturing typical behaviours of granular materials. This is significant because while hypoplasticity is phenomenological in nature, the double shearing model is based on a kinematic hypothesis and microstructural considerations, and can easily be calibrated through standard tests. Copyright © 2006 John Wiley & Sons, Ltd.     

A study on mechanical properties and fracturing behavior of Carrara marble with the flat-jointed model

Overcoming two inherent limitations presented by the standard bonded-particle model (BPM), the flat-jointed model still lacks extensive application. This study hereby conducts a comprehensive investigation on the application of the flat-jointed model in rock mechanics study. First, after a careful examination on the mechanical behavior of the flat-joint contact, a systematic calibration has been conducted to build numerical models matching the mechanical properties of Carrara marble. Second, the validated model is used to simulate the fracturing behavior of marble specimens containing a single or en-echelon flaws. Finally, the appearance of particle flow code (PFC) models has been compared with the results obtained from optical observation in order to establish a possible correlation between PFC models and real marble regarding microcracking behavior. Discussions suggest that to remove the particle size effect, the model resolution (ratio of sample dimension to the average particle diameter) should be larger than 30 in the flat-jointed model and 150 in the BPM, respectively. Results from the above investigations suggest that the flat-jointed model is capable of matching both the mechanical properties and fracturing behavior of Carrara marble. The numerical model is more reliable to predict the fracturing path at the instant of initiation than at failure. Development of initial fractures is not only controlled by the magnitude of tension force but also influenced by the extent of distribution of tension force. The appearance and components of a fracture in flat-jointed models probably reveal not only its property (tensile or shear) but also the corresponding microcracking behavior of real rocks.     

基于CLoE与Gudehus-Bauer亚塑性模型的颗粒材料应变局部化对比研究

基于CLoE与Gudehus-Bauer亚塑性模型数值模拟了平面应变条件下Hostun砂的应变局部化现象。从侧向压力和初始缺陷两个方面对比研究了两种模型所预测应变局部化的产生及演化模式。结果表明：（1）两种模型均能反映Hostun砂刚度随着侧向压力提高而增大的现象。（2）相比Gudehus-Bauer亚塑性模型,CLoE亚塑性模型所得出的应变局部化形态与试验结果更加一致。（3）CLoE亚塑性模型能够反映随着荷载增加,砂的体积先膨胀后缩小的特点。（4）相比Gudehus-Bauer亚塑性模型,CLoE亚塑性模型所得到的应变-应力曲线能够更明显地反映应变局部化带中单元的软化现象。（5）CLoE亚塑性模型能够更好地模拟由初始缺陷导致的不均匀应变。总的来说,所得的数值结果表明,CLoE亚塑性模型能够较好地模拟侧向压力和初始缺陷对应变局部化的影响,在模拟应变局部化现象方面较Gudehus-Bauer更有优势。然而,现有CLoE亚塑性模型无法考虑孔隙比,也未包含颗粒材料内尺度变量,有待进一步完善。     

Numerical modeling to predict the spread of landslide of schist’s area under climatic event of Ain El Hammam (Algeria)

The several reactivations of the landslide of Ain El Hammam (AEH) after each important weather event compel us to look closely at its triggering factors and predict its mechanisms and its evolution at the longer term. In this sense, the prediction of the slope behavior becomes necessary. This paper presents a numerical model of the AEH landslide using Plaxis® software. This model considers hydraulic effects such as precipitation and pore pressure even in the unsaturated parts. Soil and rock behaviors are described with proper elasto-plastic models named Hardening Soils and Jointed Rocks. The first model takes into account hardening on isotropic and deviatoric mechanism as well as a non-associated flow rule. The second model considers a non-isotropic elasticity with perfect plasticity along with given sliding directions. The hydraulic and mechanical models are coupled with an effective stress concept. To detect unstable areas in the landslide, we developed a Matlab® program to take into account the Hill’s bifurcation criterion, which is based on sign of the second-order work. It has been proved that this criterion allows detecting all failure modes that can appear in rate-independent materials and especially the ones that develop before the plasticity limit criterion. From such computations, we can predict the shape and position of slip surface responsible of the actual ground movement of the slope. To validate the numerical results, analysis of field measurement is included. We use high resolution of electrical tomography to delineate the geometry and position of failure surface and approve our results.     

Hypoplastic model for sands with loading surface

Although the hypoplastic models for sands have exhibited good predictive capability in monotonic loading, they are not able to reproduce memory effects and predict excessive plastic accumulation under cyclic loading. To overcome these issues, a loading surface has been incorporated into a hypoplastic model. This surface is capped and has two hardening variables. Notions from the bounding surface plasticity were borrowed in order to formulate the hardening functions. With this novel model, some salient features can be described: the model can account for the accumulation of plastic deformation, a memory effect is provided by the new surface, and stress-induced anisotropy effects observed in sands are successfully simulated. A short calibration guide of the parameters is given, and some simulations for Hostun RF loose sand and Toyoura sand are presented.     

A hypoplastic model for gas hydrate‐bearing sandy sediments

Gas hydrate‐bearing sediments (GHBSs) have been considered as a potential energy resource. In this paper, the mechanical properties of GHBS are firstly investigated by the integrated test apparatus for synthesis of GHBS using silty sand as skeleton. Triaxial tests indicate an obvious transition of stress‐strain relationship from strain hardening under low hydrate saturation and strain softening under high hydrate saturation. The hypoplastic models coupled with Drucker‐Prager criterion and the Mohr‐Coulomb criterion are proposed to analyze the stress‐strain relationship of GHBS with considering the effective porosity because of the hydrate filling in the pores of GHBS. The strain hardening and softening behaviors are well predicted with less material parameters compared with the classical models. Compared with the test results, the proposed hypoplastic models are verified to be capable of capturing the salient features of the mechanical behaviors of GHBS under the conditions of little temperature change and no hydrate dissociation.     

Mechanical and hydraulic properties of rocks related to induced seismicity

Witherspoon, P.A. and Gale, J.E., 1977. Mechanical and hydraulic properties of rocks related to induced seismicity. Eng. Geol., 11(1): 23–55.The mechanical and hydraulic properties of fractured rocks are considered with regard to the role they play in induced seismicity. In many cases, the mechanical properties of fractures determine the stability of a rock mass. The problems of sampling and testing these rock discontinuities and interpreting their non-linear behavior are reviewed. Stick slip has been proposed as the failure mechanism in earthquake events. Because of the complex interactions that are inherent in the mechanical behavior of fractured rocks, there seems to be no simple way to combine the deformation characteristics of several sets of fractures when there are significant perturbations of existing conditions. Thus, the more important fractures must be treated as individual components in the rock mass.In considering the hydraulic properties, it has been customary to treat a fracture as a parallel-plate conduit and a number of mathematical models of fracture systems have adopted this approach. Non-steady flow in fractured systems has usually been based on a two-porosity model, which assumes the primary (intergranular) porosity contributes only to storage and the secondary (fracture) porosity contributes only to the overall conductivity. Using such a model, it has been found that the time required to achieve quasi-steady state flow in a fractured reservoir is one or two orders of magnitude greater than it is in a homogeneous system. In essentially all of this work, the assumption has generally been made that the fractures are rigid.However, it is clear from a review of the mechanical and hydraulic properties that not only are fractures easily deformed but they constitute the main flow paths in many rock masses. This means that one must consider the interaction of mechanical and hydraulic effects. A considerable amount of laboratory and field data is now available that clearly demonstrates this stress-flow behavior. Two approaches have been used in attempting to numerically model such behavior: (1) continuum models, and (2) discrete models. The continuum approach only needs information as to average values of fracture spacing and material properties. But because of the inherent complexity of fractured rock masses and the corresponding decrease in symmetry, it is difficult to develop an equivalent continuum that will simulate the behavior of the entire system. The discrete approach, on the other hand, requires details of the fracture geometry and material properties of both fractures and rock matrix. The difficulty in obtaining such information has been considered a serious limitation of discrete models, but improved borehole techniques can enable one to obtain the necessary data, at least in shallow systems. The possibility of extending these methods to deeper fracture systems needs more investigation. Such data must be considered when deciding whether to use a continuum or discrete model to represent the interaction of rock and fluid forces in a fractured rock system, especially with regard to the problem of induced seismicity. When one is attempting to alter the pressure distribution in a fault zone by injection or withdrawal of fluids, the extent to which this can be achieved will be controlled in large measure by the behavior of the fractures that communicate with the borehole. Since this is essentially a point phenomenon, i.e., the changes will propagate from a relatively small region around the borehole, the use of a discrete model would appear to be preferable.     

A model for pore‐fluid‐sensitive rock behavior using a weathering state parameter

Chalk and other porous rocks are known to behave differently when saturated with different pore fluids. The behavior of these rocks varies with different pore fluids and additional deformation occurs when the pore fluid composition changes. In this article, we review the evidence that behavior in porous rocks is pore‐fluid‐dependent, present a constitutive model for pore‐fluid‐dependent porous rocks, and present a compilation of previously published data to develop quantitative relationships between various pore fluids and mechanical behavior. The model proposed here is based on a state parameter approach for weathering and has similarities to models previously proposed for weathering‐sensitive rocks in that the values for parameters that characterize material behavior vary as a function of weathering. Comparisons with published experimental data indicate that the model is capable of reproducing observed behavior of chalk under a variety of loading conditions and changes in pore fluid composition. Copyright © 2008 John Wiley & Sons, Ltd.     

http://www.sciencedirect.com/science/article/pii/S1674987114001649

Hydrocarbon production from shale has attracted much attention in the recent years. When applied to this prolific and hydrocarbon rich resource plays, our understanding of the complexities of the flow mechanism(sorption process and flow behavior in complex fracture systems- induced or natural) leaves much to be desired. In this paper, we present and discuss a novel approach to modeling, history matching of hydrocarbon production from a Marcellus shale asset in southwestern Pennsylvania using advanced data mining, pattern recognition and machine learning technologies. In this new approach instead of imposing our understanding of the flow mechanism, the impact of multi-stage hydraulic fractures, and the production process on the reservoir model, we allow the production history, well log, completion and hydraulic fracturing data to guide our model and determine its behavior. The uniqueness of this technology is that it incorporates the so-called "hard data" directly into the reservoir model, so that the model can be used to optimize the hydraulic fracture process. The "hard data" refers to field measurements during the hydraulic fracturing process such as fluid and proppant type and amount, injection pressure and rate as well as proppant concentration. This novel approach contrasts with the current industry focus on the use of "soft data"(non-measured, interpretive data such as frac length, width,height and conductivity) in the reservoir models. The study focuses on a Marcellus shale asset that includes 135 wells with multiple pads, different landing targets, well length and reservoir properties. The full field history matching process was successfully completed using this data driven approach thus capturing the production behavior with acceptable accuracy for individual wells and for the entire asset.     

非饱和土水-力本构模型及其隐式积分算法

在已有工作基础上建立了水力-力学耦合的非饱和土本构模型,在硬化方程中考虑饱和度的影响,同时在土水特征曲线中考虑了塑性体变的影响,从而使模型可以反映非饱和土中的毛细现象与土中弹塑性变形现象的耦合行为。采用隐式积分方法,建立了非饱和土耦合模型的数值模型,并推导了得到了水力-力学耦合的非饱和土的一致切线模量。利用该算法编制了本构模型计算的子程序,使其能向外输出切线刚度矩阵,用于有限元计算。为了验证该算法和程序的正确性,用所编制程序对不同路径下的土体行为进行了预测。通过预测结果与试验结果相对比,表明程序预测结果与试验数据相吻合,模型可以较好地模拟土体的水力-力学耦合行为特性。     

Evaluation of a pumping test of the Snake River Plain aquifer using axial-flow numerical modeling

The Snake River Plain aquifer in southeast Idaho is hosted in a thick sequence of layered basalts and interbedded sediments. The degree to which the layering impedes vertical flow has not been well understood, yet is a feature that may exert a substantial control on the movement of contaminants. An axial-flow numerical model, RADFLOW, was calibrated to pumping test data collected by a straddle-packer system deployed at 23 depth intervals in four observation wells to evaluate conceptual models and estimate properties of the Snake River Plain aquifer at the Idaho National Engineering and Environmental Laboratory. A delayed water-table response observed in intervals beneath a sediment interbed was best reproduced with a three-layer simulation. The results demonstrate the hydraulic significance of this interbed as a semi-confining layer. Vertical hydraulic conductivity of the sediment interbed was estimated to be about three orders of magnitude less than vertical hydraulic conductivity of the lower basalt and upper basalt units. The numerical model was capable of representing aquifer conceptual models that could not be represented with any single analytical technique. The model proved to be a useful tool for evaluating alternative conceptual models and estimating aquifer properties in this application. Electronic Publication     

土水特征曲线(SWCC)的滞回特性模拟研究

为了模拟非饱和土的水力、力学特性在多次降雨、蒸发过程中的变化特性,本文以传统域模型的基本原理为基础,推到得到了一个能够模拟SWCC滞后性的计算模型,该模型计算方法简单,易于通过程序实现,且通过与试验数据和其他模型的计算结果对比,验证了该模型的合理性,同时,这个模型可为研究在复杂吸力变化状态下非饱和土的渗流特性、强度特性以及本构模型的研究打下一定的基础。     

Directional response of a reconstituted fine‐grained soil—Part II: performance of different constitutive models

In this paper, the performance of different advanced constitutive models for soils is evaluated with respect to the experimentally observed behaviour of a soft reconstituted clay subject to a wide range of loading directions, see (presented in the companion paper). The models considered include a three‐surface kinematic hardening elastoplastic model; the CLoE hypoplastic model; a recently proposed K‐hypoplastic model for clays, and an enhanced version of the same model incorporating the concept of intergranular strain. A clear qualitative picture of the relative performance of the different models as a function of the loading direction is obtained by means of the incremental strain response envelopes. The definition of suitable error measures allows to obtain further quantitative information in this respect. For the particular initial conditions and loading programme considered in this study, the kinematic hardening and the enhanced K‐hypoplastic models appear to provide the best performance overall. Copyright © 2006 John Wiley & Sons, Ltd.     

Coupled modelling of hydro‐mechanical behaviour of unsaturated compacted expansive soils

The Barcelona basic model cannot predict the mechanical behaviour of unsaturated expansive soils, whereas the Barcelona expansive model (BExM) can only predict the stress–strain behaviour of unsaturated expansive soils without the water‐retention behaviour being incorporated. Moreover, the micro‐parameters and the coupling function between micro‐structural and macro‐structural strains in the BExM are difficult to determine. Experimental data show that the compression curves for non‐expansive soils under constant suctions are shifted towards higher void ratios with increasing suction, whereas the opposite is true for expansive soils. According to the observed water‐retention behaviour of unsaturated expansive soils, the air‐entry value increases with density, and the relationship between the degree of saturation and void ratio is linear at constant suction. According to the above observation, an elastoplastic constitutive model is developed for predicting the hydraulic and mechanical behaviour of unsaturated expansive soils, based on the existing hydro‐mechanical model for non‐expansive unsaturated soil. The model takes into consideration the effect of degree of saturation on the mechanical behaviour and that of void ratio on the water‐retention behaviour. The concept of equivalent void ratio curve is introduced to distinguish the plastic potential curve from the yield curve. The model predictions are compared with the test results of an unsaturated expansive soil, including swelling tests under constant net stress, isotropic compression tests and triaxial shear tests under constant suction. The comparison indicates that the model offers great potential for quantitatively predicting the hydraulic and mechanical behaviour of unsaturated expansive soils. Copyright © 2011 John Wiley & Sons, Ltd.     

Linear Elastic and Cohesive Fracture Analysis to Model Hydraulic Fracture in Brittle and Ductile Rocks

Hydraulic fracturing technology is being widely used within the oil and gas industry for both waste injection and unconventional gas production wells. It is essential to predict the behavior of hydraulic fractures accurately based on understanding the fundamental mechanism(s). The prevailing approach for hydraulic fracture modeling continues to rely on computational methods based on Linear Elastic Fracture Mechanics (LEFM). Generally, these methods give reasonable predictions for hard rock hydraulic fracture processes, but still have inherent limitations, especially when fluid injection is performed in soft rock/sand or other non-conventional formations. These methods typically give very conservative predictions on fracture geometry and inaccurate estimation of required fracture pressure. One of the reasons the LEFM-based methods fail to give accurate predictions for these materials is that the fracture process zone ahead of the crack tip and softening effect should not be neglected in ductile rock fracture analysis. A 3D pore pressure cohesive zone model has been developed and applied to predict hydraulic fracturing under fluid injection. The cohesive zone method is a numerical tool developed to model crack initiation and growth in quasi-brittle materials considering the material softening effect. The pore pressure cohesive zone model has been applied to investigate the hydraulic fracture with different rock properties. The hydraulic fracture predictions of a three-layer water injection case have been compared using the pore pressure cohesive zone model with revised parameters, LEFM-based pseudo 3D model, a Perkins-Kern–Nordgren (PKN) model, and an analytical solution. Based on the size of the fracture process zone and its effect on crack extension in ductile rock, the fundamental mechanical difference of LEFM and cohesive fracture mechanics-based methods is discussed. An effective fracture toughness method has been proposed to consider the fracture process zone effect on the ductile rock fracture.     

A constitutive model for unsaturated soil interfaces

A constitutive model based on the disturbed state concept is presented to describe the behavior of interfaces in unsaturated soil. The model is an extension of an existing model developed for a sand–steel interface. As opposed to the original model, the modified model incorporates two independent stress variables, which are the net normal stress and matric suction. The saturated and dry state of the interface can be modeled as a special case using the constitutive model presented in this paper. The modified model is capable of capturing the main features of unsaturated interfaces observed during laboratory testing, including increasing shear strength and strain softening with increasing suction and net normal stress and increasing dilatancy with increasing suction. Laboratory tests were carried out on unsaturated interfaces in a modified direct shear test apparatus. The observed behavior of interfaces between unsaturated soil and steel plates (rough and smooth) is presented in comparison with model predictions. Copyright © 2008 John Wiley & Sons, Ltd.     

A kinematic hardening based model for unsaturated soils considering different hydraulic conditions

At present, several of the existing elastoplastic constitutive models are adapted for describing the stress–strain behavior of unsaturated soils. However, most of them present certain limitations in this field. These limitations can be related to the basic model and/or added unsaturated state variables and formulations. In this regard, inability to model the hydro‐mechanical behavior in constant water (CW) conditions is an example of these limitations. In this paper, an advanced version of CJS model is selected for adaptation to the unsaturated states. Adaptation to unsaturated states is achieved in the framework of effective stress approach. Effective stress equation and unsaturated state variables are selected based on the recent research existing in the literature. The developed model is capable of describing the complex behavior of unsaturated soil in the CW condition in addition to predicting the behavior at failure and post–failure, nonlinear elastoplastic behavior at low levels of stress and strain (by selecting a very small elastic domain), as well as wetting and collapse behaviors. In order to validate the model, results of triaxial tests in CD and CW conditions are used. The validation results indicate the good capability of the proposed model. Behavior of the unsaturated soils during wetting is an important issue. For this reason, the model is also evaluated based on the results of wetting and collapse triaxial tests. A comparison between the tests and simulation results shows that the model is able to predict the soil behavior under the wetting path. Copyright © 2016 John Wiley & Sons, Ltd.     

Capability of constitutive models to simulate soils with different OCR using a single set of parameters

Incorporation of void ratio as a state variable into constitutive models allows, in principle, to use a single set of parameters for soils with different OCRs. Two sets of experimental data on reconstituted clays are used for evaluation of three constitutive models of different complexity (Modified Cam clay model, 3SKH model, hypoplastic model for clays). Although all the models predict the influence of OCR correctly from the qualitative point of view, quantitative comparison using a suitable scalar error measure reveals merits and shortcomings of different models.     

A thermo‐mechanical model for variably saturated soils based on hypoplasticity

The paper presents a mechanical model for non‐isothermal behaviour of unsaturated soils. The model is based on an incrementally non‐linear hypoplastic model for saturated clays and can therefore tackle the non‐linear behaviour of overconsolidated soils. A hypoplastic model for non‐isothermal behaviour of saturated soils was developed and combined with the existing hypoplastic model for unsaturated soils based on the effective stress principle. Features of the soil behaviour that are included into the model, and those that are not, are clearly distinguished. The number of model parameters is kept to a minimum, and they all have a clear physical interpretation, to facilitate the model usefulness for practical applications. The step‐by‐step procedure used for the parameter calibration is described. The model is finally evaluated using a comprehensive set of experimental data for the thermo‐mechanical behaviour of an unsaturated compacted silt. Copyright © 2011 John Wiley & Sons, Ltd.     

Prediction of hydraulic conductivity for soil–bentonite mixture

Owing to its low hydraulic conductivity, soil and bentonite mixture is applied as a liner material. However, the experimental determination of hydraulic conductivity, which is controlled by various physical, chemical and mineralogical factors, requires an expensive and time-consuming setup. In the present work, multigene symbolic, genetic programming was used to model functional relationships for hydraulic conductivity. The developed model was able to generalize highly nonlinear variations in data as well as predict system behavior from experimental observations. It was found that the predictions obtained from developed model agree well with experimental observations.