Constitutive and numerical framework for modeling joints and faults in rock

A three‐dimensional constitutive model for joints is described that incorporates nonlinear elasticity based on volumetric elastic strain, and plasticity for both compaction and shear with emphasis on compaction. The formulation is general in the sense that alternative specific functional forms and evolution equations can be easily incorporated. A corresponding numerical structure based on finite elements is provided so that a joint width can vary from a fraction of an element size to a width that occupies several elements. The latter case is particularly appropriate for modeling a fault, which is considered simply to be a joint with large width. For small joint widths, the requisite equilibrium and kinematic requirements within an element are satisfied numerically. The result is that if the constitutive equation for either the joint or the rock is changed, the numerical framework remains unchanged. A unique aspect of the general formulation is the capability to handle either pre‐existing gaps or the formation of gaps. Representative stress–strain plots are given to illustrate both the features of the model and the effects of changes in values of material parameters. Copyright © 2015 John Wiley & Sons, Ltd.     

Generalized plasticity state parameter‐based model for saturated and unsaturated soils. Part 1: Saturated state

The behavior of granular materials is known to depend on its loose or dense nature, which in turns depends both on density and confining pressure. Many models developed in the past require the use of different sets of constitutive parameters for the same material under different confining pressures. The purpose of this paper is to extend a basic generalized plasticity model for sands proposed by Pastor, Zienkiewicz and Chan by modifying the main ingredients of the model flow—rule, loading–unloading discriminating direction and plastic modulus—to include a dependency on the state parameter. The proposed model is tested against the available experimental data on three different sands, using for each of them a single set of material parameters, finding a reasonably good agreement between experiments and predictions. Copyright © 2010 John Wiley & Sons, Ltd.     

An adaptive finite element model for wellbore breakout simulation

A finite element formulation is proposed and implemented for analysing the stability of excavated wells using the DiMaggio-Sandler constitutive elastoplastic model with a typical carbonate reservoir configuration. The quality of the finite element approximation is ensured by applying smooth curved elements adapted to the wellbore geometry, and h − p adaptive finite element meshes in the plastic zone. General purpose procedures are defined to transfer the elastoplastic deformation history to newly created integration points. A breakout damage criterion is proposed based on the second invariant of the deviatoric plastic deformation tensor. This damage criterion is used to apply a mesh movement algorithm to represent material collapse. The automatic successive application of the breakout damage criterion results in elliptical realistically looking geometries obtained in experiments reported in the literature.     

Elastoplastic model with loading memory surfaces (LMS) for monotonic and cyclic behaviour of geomaterials

A new elastoplastic model called loading memory surface based on the critical state concept and the multi‐surface framework is proposed for geomaterials. The model uses a hypoelastic formulation and two plastic mechanisms. The formulations of the model are made in three‐dimensional stress–strain space and work under both monotonic and cyclic loadings. A newly introduced formalism makes it possible to obtain the cyclic response directly from the monotonic loading one. This formalism gives a three‐dimensional generalization of the well‐known Masing rule. The model has been validated against test results of Hostun sand under several conditions: monotonic and cyclic, drained and undrained, tests in compression and in extension, and at different confining pressures and different densities. Copyright © 2014 John Wiley & Sons, Ltd.     

A numerical study of mechanical behavior of a cement paste under mechanical loading and chemical leaching

Based on relevant experimental data of a petroleum cement paste under mechanical loading and chemical leaching, an elastic‐plastic model is first proposed by taking into account plastic shearing and pore collapse. The degradation of mechanical properties induced by the chemical leaching is characterized by a chemical damage variable which is defined as the increase of porosity. Both elastic and plastic properties of the cement paste are affected by the chemical damage. The proposed model is calibrated from and applied to describe mechanical responses in triaxial compression tests respectively on sound and fully leached samples. In the second part, a phenomenological chemical model is defined to establish the relationship between porosity change and calcium dissolution process. The dissolution kinetics is governed by a diffusion law taking into account the variation of diffusion coefficient with calcium concentration. The chemical model is coupled with the mechanical model, and both are applied to describe mechanical response of cement paste samples subjected to progressive chemical leaching and compressive stresses. Comparisons between experimental data and numerical results are presented.     

Generalized plasticity state parameter‐based model for saturated and unsaturated soils. Part II: Unsaturated soil modeling

The aim of this paper is to extend the generalized plasticity state parameter‐based model presented in part 1 to reproduce the hydro‐mechanical behavior of unsaturated soils. The proposed model is based on two pairs of stress–strain variables and a suitable hardening law taking into account the bonding—debonding effect of suction and degree of saturation. A generalized state parameter for unsaturated state is proposed to reproduce soil behavior using a single set of material parameters. Generalized plasticity gives a suitable framework to reproduce not only monotonic stress path but also cyclic behavior. The hydraulic hysteresis during a drying—wetting cycle and the void ratio effect on the hydraulic behavior is introduced. Comparison between model simulations and a series of experimental data available, both cohesive and granular, are given to illustrate the accuracy of the enhanced generalized plasticity equation. Copyright © 2010 John Wiley & Sons, Ltd.     

On the pore‐scale mechanisms leading to brittle and ductile deformation behavior of crystalline rocks

Deformation mechanisms at the pore scale are responsible for producing large strains in porous rocks. They include cataclastic flow, dislocation creep, dynamic recrystallization, diffusive mass transfer, and grain boundary sliding, among others. In this paper, we focus on two dominant pore‐scale mechanisms resulting from purely mechanical, isothermal loading: crystal plasticity and crofracturing. We examine the contributions of each mechanism to the overall behavior at a scale larger than the grains but smaller than the specimen, which is commonly referred to as the mesoscale. Crystal plasticity is assumed to occur as dislocations along the many crystallographic slip planes, whereas microfracturing entails slip and frictional sliding on microcracks. It is observed that under combined shear and tensile loading, microfracturing generates a softer response compared with crystal plasticity alone, which is attributed to slip weakening where the shear stress drops to a residual level determined by the frictional strength. For compressive loading, however, microfracturing produces a stiffer response than crystal plasticity because of the presence of frictional resistance on the slip surface. Behaviors under tensile, compressive, and shear loading invariably show that porosity plays a critical role in the initiation of the deformation mechanisms. Both crystal plasticity and microfracturing are observed to initiate at the peripheries of the pores, consistent with results of experimental studies. Copyright © 2015 John Wiley & Sons, Ltd.     

A plasticity constitutive model for unsaturated,anisotropic, nonexpansive soils

The paper describes and evaluates an incremental plasticity constitutive model for unsaturated, anisotropic, nonexpansive soils (CMUA). It is based on the modified Cam-Clay (MCC) model for saturated soils and enhances it by introducing anisotropy (via rotation of the MCC yield surface) and an unsaturated compressibility framework describing a double dependence of compressibility on suction and on the degree of saturation of macroporosity. As the anisotropic and unsaturated features can be activated independently, the model is downwards compatible with the MCC model. The CMUA model can simulate effectively: the dependence of compressibility on the level of developed anisotropy, uniqueness of critical state independent of the initial anisotropy, an evolving compressibility during constant suction compression, and a maximum of collapse. The model uses Bishop's average skeleton stress as its first constitutive variable, favouring its numerical implementation in commercial numerical analysis codes (eg, finite element codes) and a unified treatment of saturated and unsaturated material states.     

A smoothed cap model for variably saturated soils and its robust numerical implementation

This paper deals with the numerical implementation of a cap model for unsaturated soils. It provides a brief review of existing cap model approaches, based on which an improved model formulated in terms of generalised effective stress and matric suction is derived and described in detail. Although the proposed model is a multisurface plasticity model, it can efficiently be implemented using only single‐surface projections because of the smoothness of the model, which is obtained by construction. Numerical algorithms are provided for these single‐surface stress projections, using a single‐equation approach whenever possible. The robustness of the utilised single‐equation approaches is enhanced by proposing problem‐fitted start‐up procedures based on investigations of the nonlinear projection equations. A comparison of the model response with extensive material test data is used to validate the model and to demonstrate the robust application of the approach to silty sands and low to medium plasticity clays. Copyright © 2015 John Wiley & Sons, Ltd.