基于梯度塑性模型的多孔介质应变局部化分析的基本理论

讨论了引入梯度塑性模型情况下的含液多孔介质双重内尺度律的特征以及相互作用问题。在引入梯度塑性本构模型的情况下,一方面材料应变局部化分析数值结果的正则性得到保证;而另一方面,含液多孔介质会出现双重内尺度律参数的共同作用问题。本文给出了此时内尺度律预测的一个基本方法,并对稳定性问题进行了分析,讨论了不同情况下实波速存在的条件,给出了对于给定的渗透系数情况下实波速存在的波数区间,并对相关现象进行了解释。     

基于热-流-固耦合模型的石油钻井施工过程数值分析

岩土介质多场耦合问题需考虑诸多因素,温度、渗流及应力之间的耦合关系复杂,试验条件不易控制,且难以实现,因此,辅以数值模拟手段具有重要的意义。基于混合物理论,推导出岩土介质温度、渗流和应力耦合的数学模型及其控制方程,提出该数学模型的求解方法,以MATLAB语言为平台,将Abaqus程序作为一个模块嵌入迭代算法程序中,编制了多场耦合分析程序,并给出了2个典型算例验证该方法的有效性和实用性。然后,将建立的多场耦合模型和计算程序应用于石油钻井施工过程的模拟,重点分析井壁围岩内温度场、渗流场和应力场的变化规律,以及钻井液温度的变化对井壁稳定性的影响。研究成果对我国地下石油、核废料储存等工程设计和施工具有一定的指导意义。     

Contaminant transport with non‐equilibrium processes in unsaturated soils and implicit characteristic Galerkin scheme

A non‐equilibrium sorption—advection—diffusion model to simulate miscible pollutant transport in saturated–unsaturated soils is presented. The governing phenomena modelled in the present simulation are: convection, molecular diffusion, mechanical dispersion, sorption, immobile water effect and degradation, including both physical and chemical non‐equilibrium processes. A finite element procedure, based on the characteristic Galerkin method with an implicit algorithm is developed to numerically solve the model equations. The implicit algorithm is formulated by means of a combination of both the precise and the traditional numerical integration procedures. The stability analysis of the algorithm shows that the unconditional stability of the present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results illustrate good performance of the present algorithm in stability and accuracy, and in simulating the effects of all the mentioned phenomena governing the contaminant transport and the concentration distribution. Copyright © 2000 John Wiley & Sons, Ltd.     

Propagation of Rayleigh waves in fluid‐saturated non‐homogeneous soils with the graded solid skeleton distribution

The propagation characteristic of Rayleigh waves in a fluid‐saturated non‐homogeneous poroelastic half‐plane is addressed. Based on Biot's theory for fluid‐saturated media, which takes the inertia, fluid viscosity, mechanical coupling, compressibility of solid grains, and fluid into account, the dispersion equations of Rayleigh waves in fluid‐saturated non‐homogeneous soils/rocks are established. By considering the shear modulus of solid skeleton variation with depth exponentially, a small parameter, which reflects the relative change of shear modulus, is introduced. The asymptotic solution of the dispersion equation expressing the relationship between the phase velocity and wave number is obtained by using the perturbation method. In order to analyze the effects of non‐homogeneity on the propagation characteristic of Rayleigh waves, the variation of the phase velocity with the wave number is presented graphically and discussed through numerical examples. Copyright © 2016 John Wiley & Sons, Ltd.     

Development of a technique for modelling clay liner desiccation

This paper presents a model for the analysis of clay liner desiccation in a landfill barrier system due to temperature effects. The model incorporates consideration of fully coupled heat‐moisture‐air flow, a non‐linear constitutive relationship, the dependence of void ratio and volumetric water content on stress, capillary pressure and temperature, and the effect of mechanical deformation on all governing equations. Mass conservative numerical schemes are proposed to improve the accuracy of the finite element solution to the governing equations. The application of the model is then demonstrated by examining three test problems, including isothermal infiltration, heat conduction and non‐isothermal water and heat transport. Comparisons are made with results from literature, and good agreement is observed. Copyright © 2003 John Wiley & Sons, Ltd.     

Non‐local damage modelling of concrete: a procedure for the determination of model parameters

This paper presents a procedure for the determination of parameters of non‐local damage models. This is to assure a consistent response of a non‐local damage model, as choice of the internal length and other parameters of the model are varied. Correlations between the internal length and other parameters governing the local constitutive behaviour of the model are addressed and exploited. Focus is put on the relationship between the internal length of the non‐local model and the width of the fracture process zone. Numerical examples are used to demonstrate the rigour of the proposed method. Copyright © 2006 John Wiley & Sons, Ltd.     

Pore‐scale modeling of fluid‐particles interaction and emerging poromechanical effects

A micro‐hydromechanical model for granular materials is presented. It combines the discrete element method for the modeling of the solid phase and a pore‐scale finite volume formulation for the flow of an incompressible pore fluid. The coupling equations are derived and contrasted against the equations of conventional poroelasticity. An analogy is found between the discrete element method pore‐scale finite volume coupling and Biot's theory in the limit case of incompressible phases. The simulation of an oedometer test validates the coupling scheme and demonstrates the ability of the model to capture strong poromechanical effects. A detailed analysis of microscale strain and stress confirms the analogy with poroelasticity. An immersed deposition problem is finally simulated and shows the potential of the method to handle phase transitions. Copyright © 2013 John Wiley & Sons, Ltd.     

Isogeometric analysis of THM coupled processes in ground freezing

An isogeometric analysis (IGA) based numerical model is presented for simulation of thermo-hydro-mechanically (THM) coupled processes in ground freezing. The momentum, mass and energy conservation equations are derived based on porous media theory. The governing equations are supplemented by a saturation curve, a hydraulic conductivity model and constitutive equations. Variational and Galerkin formulation results in a highly nonlinear system of equations, which are solved using Newton-Raphson iteration. Numerical examples on isothermal consolidation in plane strain, one-dimensional freezing and heave due to a chilled pipeline are presented. Reasonably good agreements were observed between the IGA based heave simulations and experimental results.     

A discretized analytical solution for fully coupled non‐linear simulation of heat and mass transfer in poroelastic unsaturated media

Mathematical simulation of non‐isothermal multiphase flow in deformable unsaturated porous media is a complicated issue because of the need to employ multiple partial differential equations, the need to take into account mass and energy transfer between phases and because of the non‐linear nature of the governing partial differential equations. In this paper, an analytical solution for analyzing a fully coupled problem is presented for the one‐dimensional case where the coefficients of the system of equations are assumed to be constant for the entire domain. A major issue is the non‐linearity of the governing equations, which is not considered in the analytical solution. In order to introduce the non‐linearity of the equations, an iterative discretized procedure is used. The domain of the problem is divided into identical time–space elements that cover the time–space domain. A separate system of equations is defined for each element in the local coordinate system, the initial and boundary conditions for each element are obtained from the adjacent elements and the coefficients of the system of equations are considered to be constant in each step. There are seven governing differential equations that should be solved simultaneously: the equilibrium of the solid skeleton, mass conservation of fluids (water, water vapor and gas) and energy conservation of phases (solid, liquid and gas). The water vapor is not in equilibrium with water and different phases do not have the same temperature. The governing equations that have been solved seem to be the most comprehensive in this field. Three examples are presented for analyzing heat and mass transfer in a semi‐infinite column of unsaturated soil. Copyright © 2009 John Wiley & Sons, Ltd.     

Coupling of finite and hierarchical infinite elements: application to a non‐homogeneous cross‐anisotropic half‐space subjected to a non‐uniform circular loading

A method is presented for coupling cubic‐order quadrilateral finite elements with the finite side of a new coordinate ascent hierarchical infinite element. At a common side shared by a hierarchical infinite element and an arbitrary number of finite elements, the displacements are minimized in the least square sense with respect to the degrees‐of‐freedom of the finite elements. This leads to a set of equations that relate the degrees‐of‐freedom of the finite and hierarchical infinite elements on the shared side. The method is applied to a non‐homogeneous cross‐anisotropic half‐space subjected to a non‐uniform circular loading with Young's and shear moduli varying with depth according to the power law. A constant mesh constructed from coupled finite and hierarchical infinite elements is used and convergence is sought simply by increasing the degree of the interpolating polynomial. The displacements and stresses produced by conical and parabolic circular loads applied on the surface are obtained. The efficiency of the proposed method is demonstrated through convergence and comparison studies. New results produced by a frusto‐conical circular load applied on the surface of a half‐space made up of heavily consolidated London clay are provided. The non‐homogeneity parameter and degree of anisotropy are shown to influence the soil response. Copyright © 2012 John Wiley & Sons, Ltd.     

Thermo-hydro-mechanical modeling of unsaturated soils using isogeometric analysis: Model development and application to strain localization simulation

This study presents a thermo-hydro-mechanical (THM) model of unsaturated soils using isogeometric analysis (IGA). The framework employs Bézier extraction to connect IGA to the conventional finite element analysis (FEA), featuring the current study as one of the first attempts to develop an IGA-FEA framework for solving THM problems in unsaturated soils. IGA offers higher levels of interelement continuity making it an attractive method for solving highly nonlinear problems. The governing equations of linear momentum, mass, and energy balance are coupled based on the averaging procedure within the hybrid mixture theory. The Drucker-Prager yield surface is used to limit the modified effective stress where the model follows small strain, quasi-static loading conditions. Temperature dependency of the surface tension is implemented in the soil-water retention curve. Nonuniform rational B-splines (NURBS) basis functions are used in the standard Galerkin method and weak formulations of the balance equations. Displacement, capillary pressure, gas pressure, and temperature are four independent quantities that are approximated by NURBS in spatial discretization. The framework is used to simulate strain localization in an undrained dense sand subjected to plane strain biaxial compression under different temperatures and displacement velocities. Results show that an increase in the displacement rate leads to reduction in the equivalent plastic strain while an increase in the temperature leads to an increase in the equivalent plastic strain. The findings suggest that the proposed IGA-based framework offers a viable alternative for solving THM problems in unsaturated soils.     

Numerical solution for consolidation and desiccation of soft soils

The consolidation and desiccation behaviour of soft soils can be described by two time‐dependent non‐linear partial differential equations using the finite strain theory. Analytical solutions do not exist for these governing equations. In this paper, we develop efficient numerical methods and software for finding the numerical solutions. We introduce a semi‐implicit time integration scheme, and show numerically that our method converges. In addition, the numerical solution matches well with the experimental result. A boundary refinement method is also developed to improve the convergence and stability for the case of Neumann type boundary conditions. Interface governing equations are derived to maintain the continuity of consolidation and desiccation processes. This is useful because the soil column can undergo desiccation on top and consolidation on the bottom simultaneously. The numerical algorithms has been implemented into a computer program and the results have been verified with centrifuge test results conducted in our laboratory. Copyright © 2001 John Wiley & Sons, Ltd.     

Simulation of localization failure with strain‐gradient‐enhanced damage mechanics

Strain gradient implies an important characteristic in localized damage deformation, which can be observed in the softening state of brittle materials, and strain gradients constitute the basic behaviours of localization failure area of the materials. The most important point in strain gradient is its damaging function including an internal length scale, which can be used to express the scale effects of mechanical responses of brittle rock mass. By extending the strain gradient theory and introducing an intrinsic material length scale into the constitutive law, the authors develop an isotropic damage model as well as a micro‐crack‐based anisotropic damage model for rock‐like materials in this paper. The proposed models were used to simulate the damage localization under uniaxial tension and plain strain compression, respectively. The simulated results well illustrated the potential of these models in dealing with the well‐known mesh‐sensitivity problem in FEM. In the computation, elements with C₁ continuity have been implemented to incorporate the proposed models for failure localization. When regular rectangle elements are encountered, the coupling between finite difference method (FDM) and conventional finite element method (FEM) is used to avoid large modification to the existing FEM code, and to obtain relatively higher efficiency and reasonably good accuracy. Application of the anisotropic model to the 3D‐non‐linear FEM analysis of Ertan arch dam has been conducted and the results of its numerical simulation coincide well with those from the failure behaviours obtained by Ertan geophysical model test. In this paper, new applications of gradient theories and models for a feasible approach to simulate localized damage in brittle materials are presented. Copyright © 2002 John Wiley & Sons, Ltd.     

Spectral analysis of nonlocal regularization in two‐dimensional finite element models

Strain‐softening in geomaterials often leads to ill‐posed boundary‐valued problems (BVP), which cannot be solved with finite element methods without introducing some kind of regularization such as nonlocal plasticity. Hereafter we propose to apply spectral analysis for testing the performance of nonlocal plasticity in regularizing ill‐posed BVP and producing mesh‐independent solutions when local plasticity usually fails. The spectral analysis consists of examining the eigenvalues and eigenvectors of the global tangential stiffness matrix of the incremental equilibrium equations. Based on spectral analysis, we propose a criterion for passing or failing the test of constitutive regularization in the context of BVP. If the eigenvalues of the tangential operator are all positive then the regularization succeeds, otherwise it fails and may not prevent artificial mesh‐dependent solutions from appearing. The approach is illustrated in the particular case of a biaxial compression with strain‐softening plasticity. In this particular case, local softening plasticity is found to produce negative eigenvalues in the tangential stiffness matrix, which indicates ill‐posed BVP. In contrast, nonlocal softening plasticity always produces positive eigenvalues, which regularizes ill‐posed BVP. The dominant eigenvectors, which generate localized deformation patterns, have a bandwidth independent of mesh size, provided that the mesh is fine enough to capture localization. These mesh‐independent eigenmodes explain why nonlocal plasticity produces numerical solutions that are mesh‐independent. Copyright © 2011 John Wiley & Sons, Ltd.     

Two‐scale modeling of solute dispersion in unsaturated double‐porosity media: Homogenization and experimental validation

A two‐scale modeling of solute transport in double‐porosity (DP) media under unsaturated water flow conditions is presented. The macroscopic model was developed by applying the asymptotic homogenization method. It is based on theoretical and empirical considerations dealing with the orders of magnitude of characteristic quantities involved in the process. For this purpose a physical model that mimics the behavior of DP medium was built. The resulting two‐equation model relies on a coupling exchange term between micro‐ and macro‐porosity subdomains associated with local non‐equilibrium solute concentrations. The model was numerically implemented (Comsol Multiphysics^®) to simulate the macroscopic one‐dimensional physical process taking place into the porous medium of 3D periodic microstructure. A series of dispersion experiments of NaCl solution under unsaturated steady‐state flow conditions were performed. The experimental results were used first to calibrate the dispersion coefficient of the model, and second to validate it through two other independent experiments. The excellent agreement between the numerical simulations and the measurements of the time evolution of the non‐symmetrical breakthrough curves provides a proof of predictive capacity of the developed model. Copyright © 2010 John Wiley & Sons, Ltd.     

Modeling coupled THM processes of geological porous media with multiphase flow: Theory and validation against laboratory and field scale experiments

A FEM model for analysis of fully coupled multiphase flow, thermal transport and stress/deformation in geological porous media was developed based on the momentum, mass and energy conservation laws of the continuum mechanics and the averaging approach of the mixture theory over a three phase (solid–liquid–gas) system. Six processes (i.e. stress–strain, water flow, gas flow, vapor flow, heat transport and porosity evolution processes) and their coupling effects are considered, which not only makes the problem well-defined, but renders the governing PDEs closed, complete, compact and compatible. Displacements, pore water pressure, pore gas pressure, pore vapor pressure, temperature and porosity are selected as basic unknowns. The physical phenomena such as phase transition, gas solubility in liquid, thermo-osmosis, moisture transfer and moisture swelling are modeled. As a result, the relative humidity and other related variables in porous media can be evaluated on a sounder physical basis. A three dimensional computer code, THYME3D, was developed, with eight degrees of freedom at each node. The laboratory CEA Mock-up test and the field scale FEBEX benchmark test on bentonite performance assessment for underground nuclear waste repositories were used to validate the numerical model and the software. The coupled THM behaviors of the bentonite barriers were satisfactorily simulated, and the effects and impacts of the governing equations, constitutive relations and property parameters on the coupled THM processes were understood in terms of more straightforward interpretation of physical processes at microscopic scale of the porous media. The work developed enables further in-depth research on fully coupled THM or THMC processes in porous media.     

Elasto‐plastic model for cement‐treated sand

This paper presents an elasto‐plastic model for non‐linear analyses of cement‐treated sand. Various laboratory tests were systematically carried out to investigate the pre‐peak and post‐peak behaviours of a cement‐treated sand. On the basis of these experimental results, the new model was built within the framework of a relatively simple elasto‐plastic theory. Two failure criteria are employed to express tensile and shear failure characteristics observed in the experimental results of the cement‐treated sand. The proposed model can describe strain‐hardening and strain‐softening responses under both failure modes. In the strain‐softening rules, the smeared crack concept is used, and a characteristic length is considered to avoid the issue of mesh‐size dependency. Since the failure criterion and strain‐hardening/softening rules are based on the experimental evidences, the model is relatively easy to understand and the parameters used in the model have clear physical meaning. The proposed model was applied to simulate the behaviour of cement‐treated sand in various laboratory tests, allowing for a reasonable comprehensive evaluation. It was demonstrated that the proposed model is suitable for describing both the tensile and shear failure behaviours of cement‐treated sand. Copyright © 2006 John Wiley & Sons, Ltd.     

Incrementalization of a single hardening constitutive model for frictional materials

The governing equations for an elasto‐plastic constitutive model for frictional materials such as soil, rock, and concrete are presented, and the incremental form is indicated in preparation for implementation of the model in a user‐defined module for finite element calculations. This isotropic, work‐hardening and ‐softening model employs a single yield surface, it incorporates non‐associated plastic flow, and its capability of capturing the behaviour of different types of frictional materials under various three‐dimensional conditions has been demonstrated by comparison with measured behaviour, as presented in the literature. The incrementalization procedure is indicated and the resulting equations for the single hardening model are presented together with parameters for a dense sand. Following the implementation of the model, these parameters are used for evaluation of different integration schemes as presented in a companion paper by Jakobsen and Lade (Int. J. Numer. Anal. Meth. Geomech. 2002; 26 :661). Copyright © 2002 John Wiley & Sons, Ltd.