Discussion on the dynamic mechanism of Great North China area based on the observed stress data

Introduction The Great North China, located at longitude 106E to 124E and latitude 31N to 42N, in-cludes three secondary active tectonic blocks, Ordos, Yanshan and North China plain (Figure 1). The geological tectonics of these three secondary blocks is much different from each other. As a stable block with high rigidity, the Ordos block is mostly surrounded by down-faulted basins with an inactive interior since Cenozoic, although the fault zones along its boundary are strongly active wi…     

Preliminary study on variation characteristics of ocean tide dynamic stress in crust and its relationship with earthquakes

~~Preliminary study on variation characteristics of ocean tide dynamic stress in crust and its relationship with earthquakes@陆明勇$China Center for Earthquake Disaster Emergency and SAR,Beijing 100049,China @郑文衡$Huazhong University of Science and Technology,Wuhan 430074,China①Strong earthquake catalogues in China. Earthquake catalogues reported monthly by Center for Analysis and Prediction, CEA. ①YI Zhi-gang. Institute of Crustal Dynamics, China Earthquake Administrati…     

Analysis of pressuremeter geometry effects in clay using critical state models

The conventional interpretation methods of pressuremeter testing effectively approximate pressuremeter membranes as infinitely long. As a result, the effects of the two‐dimensional geometry of pressuremeters are ignored, leading to an overestimation of soil shear strength by pressuremeter testing, as demonstrated in several previous studies. This paper presents results of a numerical study of two‐dimensional geometry effects on self‐boring pressuremeter tests in undrained clay. The results are obtained using critical state soil models with an effective stress formulation. This is in contrast to most (if not all) existing studies on pressuremeter geometry effects, which were based on perfectly plastic soil models (e.g. Yu (Cavity expansion theory and its application to the analysis of pressuremeters. DPhil Thesis, The University of Oxford, 1990), Yeung and Carter (Proc. 3rd Int. Symp. on Pressuremeters, 1990), and Houlsby and Carter (Géotechnique, 1993; 43 (4):567–576)). The present study suggests that the overestimation of soil strength due to the neglect of finite pressuremeter length is significantly affected by the soil model used in the calculations. It is found that for clays with a high overconsolidation ratio (OCR) the strength overestimation predicted using critical state soil models could be considerably smaller than that predicted using perfectly plastic soil models. The main conclusion of this numerical study is that care must be exercised before directly applying any numerically determined pressuremeter geometry correction factors in practice. Copyright © 2005 John Wiley & Sons, Ltd.     

A multiphase approach for evaluating the horizontal and rocking impedances of pile group foundations

This paper advocates the use of a multiphase model, already developed for static or quasi‐static geotechnical engineering problems, for simulating the behaviour of piled raft foundations subject to horizontal as well as rocking dynamic solicitations. It is shown that such a model, implemented in a FEM code, yields appropriate predictions for the foundation impedance characteristics, provided that shear and bending effects in the piles are taken into account, thus corroborating the findings of the asymptotic homogenization theory. Besides, it is notably pointed out that such a multiphase‐based computational tool makes it possible to assess the dynamic behaviour of pile groups in a much quicker way than when using direct numerical simulations, which may face oversized problems when large pile groups are concerned. Copyright © 2016 John Wiley & Sons, Ltd.     

Transient seepage analysis in zoned anisotropic soils based on the scaled boundary finite‐element method

The scaled boundary finite‐element method, a semi‐analytical computational scheme primarily developed for dynamic stiffness of unbounded domains, is applied to the analysis of unsteady seepage flow problems. This method is based on the finite‐element technology and gains the advantages of the boundary element method as well. Only boundary of the domain is discretized, no fundamental solution is required and singularity problems can be modeled rigorously. Anisotropic and non‐homogeneous materials satisfying similarity are modeled with no additional efforts. In this study, firstly, formulation of the method for the transient seepage flow problems is derived followed by its solution procedures. The accuracy, simplicity and applicability of the method are demonstrated via four numerical examples of transient seepage flow – three of them are available in the literature. Homogenous, non‐homogenous, isotropic and anisotropic material properties are considered to show the versatility of the technique. Excellent agreement with the finite‐element method is observed. The method out‐performs the finite‐element method in modeling singularity points. Copyright © 2014 John Wiley & Sons, Ltd.     

Computational model coupling mode II discrete fracture propagation with continuum damage zone evolution

We propose a numerical method that couples a cohesive zone model (CZM) and a finite element‐based continuum damage mechanics (CDM) model. The CZM represents a mode II macro‐fracture, and CDM finite elements (FE) represent the damage zone of the CZM. The coupled CZM/CDM model can capture the flow of energy that takes place between the bulk material that forms the matrix and the macroscopic fracture surfaces. The CDM model, which does not account for micro‐crack interaction, is calibrated against triaxial compression tests performed on Bakken shale, so as to reproduce the stress/strain curve before the failure peak. Based on a comparison with Kachanov's micro‐mechanical model, we confirm that the critical micro‐crack density value equal to 0.3 reflects the point at which crack interaction cannot be neglected. The CZM is assigned a pure mode II cohesive law that accounts for the dependence of the shear strength and energy release rate on confining pressure. The cohesive shear strength of the CZM is calibrated by calculating the shear stress necessary to reach a CDM damage of 0.3 during a direct shear test. We find that the shear cohesive strength of the CZM depends linearly on the confining pressure. Triaxial compression tests are simulated, in which the shale sample is modeled as an FE CDM continuum that contains a predefined thin cohesive zone representing the idealized shear fracture plane. The shear energy release rate of the CZM is fitted in order to match to the post‐peak stress/strain curves obtained during experimental tests performed on Bakken shale. We find that the energy release rate depends linearly on the shear cohesive strength. We then use the calibrated shale rheology to simulate the propagation of a meter‐scale mode II fracture. Under low confining pressure, the macroscopic crack (CZM) and its damaged zone (CDM) propagate simultaneously (i.e., during the same loading increments). Under high confining pressure, the fracture propagates in slip‐friction, that is, the debonding of the cohesive zone alternates with the propagation of continuum damage. The computational method is applicable to a range of geological injection problems including hydraulic fracturing and fluid storage and should be further enhanced by the addition of mode I and mixed mode (I+II+III) propagation. Copyright © 2016 John Wiley & Sons, Ltd.     

A model for moisture and heat flow in fractured unsaturated porous media

The present study investigates propagation of a cohesive crack in non‐isothermal unsaturated porous medium under mode I conditions. Basic points of skeleton deformation, moisture, and heat transfer for unsaturated porous medium are presented. Boundary conditions on the crack surface that consist of mechanical interaction of the crack and the porous medium, water, and heat flows through the crack are taken into consideration. For spatial discretization, the extended finite element method is used. This method uses enriched shape functions in addition to ordinary shape functions for approximation of displacement, pressure, and temperature fields. The Heaviside step function and the distance function are exploited as enrichment functions for representing the crack surfaces displacement and the discontinuous vertical gradients of the pressure and temperature fields along the crack, respectively. For temporal discretization, backward finite difference scheme is applied. Problems solved from the literature show the validity of the model as well as the dependency of structural response on the material properties and loading. Copyright © 2016 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.