DEM–FEM analysis of soil failure process via the separate edge coupling method

The concurrent multiscale method, which couples the discrete element method (DEM) for predicting the local micro‐scale evolution of the soil particle skeleton with the finite element method (FEM) for estimating the remaining macro‐scale continuum deformation, is a versatile tool for modeling the failure process of soil masses. This paper presents the separate edge coupling method, which is degenerated from the generalized bridging domain method and is good at eliminating spurious reflections that are induced by coupling models of different scales, to capture the granular behavior in the domain of interest and to coarsen the mesh to save computational cost in the remaining domain. Cundall non‐viscous damping was used as numerical damping to dissipate the kinetic energy for simulating static failure problems. The proposed coupled DEM–FEM scheme was adopted to model the wave propagation in a 1D steel bar, a soil slope because of the effect of a shallow foundation and a plane‐strain cone penetration test (CPT). The numerical results show that the separate edge coupling method is effective when it is adopted for a problem with Cundall non‐viscous damping; it qualitatively reproduces the failure process of the soil masses and is consistent with the full micro‐scale discrete element model. Stress discontinuity is found in the coupling domain. Copyright © 2017 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.     

Three‐dimensional responses of buried corrugated pipes and ANN‐based method for predicting pipe deflections

This study investigated localized responses, such as circumferential stresses, on corrugation and pipe deflections. Also, this study examined the effect of corrugation geometry on the overall and localized response of corrugated pipes with refined three‐dimensional modeling of the entire soil–pipe interaction system, including corrugation. To investigate the availability of the traditional two‐dimensional method, the results from the three‐dimensional finite element method (FEM) were compared with those from the two‐dimensional FEM. The soil–pipe modeling techniques of this study were verified by comparing the FEM results by Utah State University and analytical results. An artificial neural network (ANN)‐based model to predict vertical deflections of buried corrugated pipes was developed to overcome the shortcomings of existing methods and obtain results that are close to the level of accuracy of FEM results. In order to train an ANN, analyses on a large amount of data were executed with various standardized pipe geometries and burial depths regulated by the Korea Highway Corporation using the two‐dimensional FEM verified in this study. The widely used back propagation algorithm was adopted. The ANN‐based model developed in this study was shown to be an effective tool by comparing the results with test data and sensitivity analyses were executed based on the data from the developed ANN. Copyright © 2010 John Wiley & Sons, Ltd.     

基于二阶锥规划有限元增量加载法的条形浅基础极限承载力分析

地基极限承载力分析是土力学研究中的一个经典课题。基于Hellinger-Reissner混合变分原理和有限元方法,将岩土体弹塑性问题构造成基于有限元框架的二阶锥规划（second-order cone programming,SOCP）问题,进而提出一种基于二阶锥规划理论的增量有限元法,即FEM-SOCP法。将岩土体弹塑性问题构造成二阶锥规划的数学优化问题,可以避免采用传统弹塑性计算中复杂的应力点积分等算法和屈服面棱角的平滑处理。此外,对于二阶锥规划问题,可以采用具有原始?对偶内点求解法的标准数学规划求解器MOSEK进行求解。将增量加载FEM-SOCP法应用于经典的基底粗糙的条形浅基础地基极限承载力分析中,分别考虑了关联和非关联塑性条件下的Mohr-Coulomb屈服准则。数值结果表明：所提出的增量加载FEM-SOCP法获得的地基承载力系数及地基承载力与传统FEM计算结果基本一致,而与常规有限元计算结果相比,基于增量加载的FEM-SOCP法所获得的屈服区更加平滑。     

Elasto-plastic analysis of an underground opening by FEM and coupled FEBEM

The elasto-plastic analysis of a circular underground opening using the finite element method (FEM), the coupled finite element and the boundary element method (FEBEM) is presented. The coupling procedure and elasto-plastic formulation are discussed. The effect of in situ stress ratio on yielded zone, displaced shape and principal stresses is presented. The results of FEBEM analysis are compared with those obtained from FEM. The computation time and the number of iterations as required by FEBEM and FEM are compared.     

A 3D composite finite element used to model a vertical drain and its surrounding remoulded soil

This work addresses an enrichment technique for the three-dimensional (3D) finite element (FE) analysis of a vertical drain foundation because (1) 1D and 2D simulations are insufficient to integrally describe the consolidation behaviour and (2) drains are small both in spacing and size, resulting in enormous computational costs for a traditional 3D FE analysis. Based on the idea of the semi-analytical finite element method (FEM), which combines general FEM with the high accuracy of a closed-form solution, a new spatial element that contains a drain well and its neighbouring smear zone is presented. This new combined element is depicted by eight global independent nodes and two local dependent nodes, and a classical analytical theory is introduced to set up the relationship between the two kinds of nodes. Because permeability diversity between the drain and the smear zone is considered, both the effects of smearing and well resistance are taken into account with the composite element method (CEM). A detailed derivation of the CEM is performed using the weighted residual method. The accuracy of the proposed method is validated with a totally penetrating, single-drain ground analysis for seven calculation conditions. Additionally, the proposed CEM saves 1/4–1/2 mesh elements and helps to avoid slender elements for the FEM analysis of the drained foundation.     

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.     

A three‐dimensional integral equation model for calculating poro‐ and thermoelastic stresses induced by cold water injection into a geothermal reservoir

Poro‐mechanical and thermo‐mechanical processes change the fracture aperture and thus affect the water flow pattern in the fracture during the cold water injection into enhanced geothermal systems (EGS). In addition, the stresses generated by these processes contribute to the phenomenon of reservoir seismicity. In this paper, we present a three‐dimensional (3D) partially coupled poro‐thermoelastic model to investigate the poroelastic and thermoelastic effects of cold water injection in EGS. In the model, the lubrication fluid flow and the convective heat transfer in the fracture are modeled by the finite element method, while the pore fluid diffusion and heat conductive transfer in the reservoir matrix are assumed to be 3D and modeled by the boundary integral equation method without the need to discretize the reservoir. The stresses at the fracture surface and in the reservoir matrix are obtained from the numerical model and can be used to assess the variation of in situ stress and induced seismicty with injection/extraction. Application of the model shows that rock cooling induces large tensile stresses and increases fracture conductivity, whereas the rock dilation caused by fluid leakoff decreases fracture aperture and increases compressive total stresses around the injection zone. However, increases in pore pressure reduce the effective stresses and can contribute to rock failure, fracture slip, and microseismic activity. Copyright © 2009 John Wiley & Sons, Ltd.     

Coupled HM analysis using zero‐thickness interface elements with double nodes. Part I: Theoretical model

In recent years, the authors have proposed a new double‐node zero‐thickness interface element for diffusion analysis via the finite element method (FEM) (Int. J. Numer. Anal. Meth. Geomech. 2004; 28 (9): 947–962). In the present paper, that formulation is combined with an existing mechanical formulation in order to obtain a fully coupled hydro‐mechanical (or HM) model applicable to fractured/fracturing geomaterials. Each element (continuum or interface) is formulated in terms of the displacements (u) and the fluid pressure (p) at the nodes. After assembly, a particular expression of the traditional ‘u–p’ system of coupled equations is obtained, which is highly non‐linear due to the strong dependence between the permeability and the aperture of discontinuities. The formulation is valid for both pre‐existing and developing discontinuities by using the appropriate constitutive model that relates effective stresses to relative displacements in the interface. The system of coupled equations is solved following two different numerical approaches: staggered and fully coupled. In the latter, the Newton–Raphson method is used, and it is shown that the Jacobian matrix becomes non‐symmetric due to the dependence of the discontinuity permeability on the aperture. In the part II companion paper (Int. J. Numer. Anal. Meth. Geomech. 2008; DOI: 10.1002/nag.730 ), the formulation proposed is verified and illustrated with some application examples. Copyright © 2008 John Wiley & Sons, Ltd.     

The 2D hydro‐mechanically coupled response of a rock mass with fractures via a mixed BEM–FEM technique

This paper describes the essential features of a numerical technique for the simulation of the coupled fluid flow and deformation in a 2D assembly of poroelastic blocks and transmissive fractures. The boundary element method (BEM) is applied to each block to reduce Navier and diffusion equations to a set of integral equations involving block boundary terms, whereas a Galerkin weighted‐residuals finite element method (FEM) is applied to the fracture diffusion equations. In addition, fracture local equilibrium is rendered through spring‐like equations relating the stresses to the relative displacements of the fracture walls. A time‐marching process is implemented leading to an algebraic system where the right‐hand side vector is built based on the collected solutions of the previous time steps. The technique requires the meshing of the fracture network only. The accuracy of the results is adequate even with relatively coarse meshes without the resort to small time steps at the beginning of the simulation. It furnishes outputs that focus only on the salient features of the response. The efficiency of the technique is demonstrated through the illustration of the results of three examples. Copyright © 2007 John Wiley & Sons, Ltd.     

A finite element algorithm for parameter identification of material models for fluid saturated porous media

In this contribution an algorithm for parameter identification of geometrically linear Terzaghi–Biot‐type fluid‐saturated porous media is proposed, in which non‐uniform distributions of the state variables such as stresses, strains and fluid pore pressure are taken into account. To this end a least‐squares functional consisting of experimental data and simulated data is minimized, whereby the latter are obtained with the finite element method. This strategy allows parameter identification based on in situ experiments. In order to improve the efficiency of the minimization process, a gradient‐based optimization algorithm is applied, and therefore the corresponding sensitivity analysis for the coupled two‐phase problem is described in a systematic manner. For illustrative purpose, the performance of the algorithm is demonstrated for a slope stability problem, in which a quadratic Drucker–Prager plasticity model for the solid and a linear Darcy law for the fluid are combined. Copyright © 2001 John Wiley & Sons, Ltd.     

An FEM/VOF hybrid formulation for fracture grouting modelling

A numerical model using a hybrid formulation of a finite element method (FEM) coupled with the volume of fluid (VOF) technique to simulate the fracture grouting processes in soils is described. The numerical model considered the couplings of the stress distribution, with two-phase fluid flows, and the mesh element damage. The hardening of grout in soil is described by a time-dependent Young’s modulus and viscosity. Crack initiation, branching, propagation, and grout vein growth in homogeneous and heterogeneous soils can be numerically reproduced. Although the method is developed particularly for simulating fracture grouting, the processes of compaction grouting and permeation grouting can also be numerically simulated. Some grouting cases have been simulated with results similar to the experimental results. This further confirms the adequacy and the power of the numerical approach.     

Implementation of a Locally Conservative Numerical Subgrid Upscaling Scheme for Two-Phase Darcy Flow

We present a locally mass conservative scheme for the approximation of two-phase flow in a porous medium that allows us to obtain detailed fine scale solutions on relatively coarse meshes. The permeability is assumed to be resolvable on a fine numerical grid, but limits on computational power require that computations be performed on a coarse grid. We define a two-scale mixed finite element space and resulting method, and describe in detail the solution algorithm. It involves a coarse scale operator coupled to a subgrid scale operator localized in space to each coarse grid element. An influence function (numerical Greens function) technique allows us to solve these subgrid scale problems independently of the coarse grid approximation. The coarse grid problem is modified to take into account the subgrid scale solution and solved as a large linear system of equations posed over a coarse grid. Finally, the coarse scale solution is corrected on the subgrid scale, providing a fine grid representation of the solution. Numerical examples are presented, which show that near-well behavior and even extremely heterogeneous permeability barriers and streaks are upscaled well by the technique.     

多尺度有限单元法求解非均质多孔介质中的三维地下水流问题

应用多尺度有限单元法模拟非均质多孔介质中的三维地下水流问题。与传统有限单元法相比，多尺度有限单元法的基函数具有能反映单元内参数变化的优点，所以这种方法能在大尺度上抓住解的小尺度特征获得较精确的解。在介绍多尺度有限单元法求解非均质多孔介质中三维地下水流问题的基本原理之后，对参数水平方向渐变垂直方向突变的非均质多孔介质中的三维地下水流和Borden实验场的三维地下水流分别用多尺度有限单元法和传统等参有限单元法进行了计算，结果表明在模拟高度非均质多孔介质中的三维地下水流问题时，多尺度有限单元法比传统有限单元法有效，既节省计算量又有较高的精度；在模拟非均质性弱的多孔介质中的三维地下水流问题时，多尺度有限单元法虽然也能在大尺度上获得较为精确的解，但效果不明显。     

一种基于有限元的岩石圈长期变形数值计算方法

有限单元法以其灵活性和精确性,成为固体地球科学中广为使用的数值方法。从短周期的地震活动到长周期的岩石圈变形、地幔对流,甚至行星演化,有限单元法几乎在固体地球科学的各个领域都占据着十分重要的位置。随着研究的深入,某些特定的地学问题给有限元计算带来挑战,尤其是岩石圈尺度大变形的数值计算,比如俯冲带的演化、剪切带中塑性流变导致的应力集中。基于显式有限元,尝试考虑粘弹塑性岩石圈大变形过程的数值计算。应用Marker-In-Cell（MIC）方法处理物质迁移。在描述基本原理和流程的基础上,对粘弹性变形、弹塑性变形、大变形过程及热传递过程等核心模块分别做了基准测试,而这四个模块是模拟岩石圈长期变形的关键。由测试结果和其他学者的（解析或数值）研究结果比对情况来看,受测试的核心模块基本达到了测试要求。可以预见,现有的基本算法可以满足研究岩石圈大变形的需要,进一步的具体研究工作将探讨这类问题。从科学问题层面讲,逐渐复杂的科学问题有利于数值模型的成熟。已达到基准测试的数值方法对下一步开展一些具体的地球动力学数值模拟研究有实际意义。      

Multiscale insights into classical geomechanics problems

We pay a revisit to some classical geomechanics problems using a novel computational multiscale modelling approach. The multiscale approach employs a hierarchical coupling of the finite element method (FEM) and the discrete element method. It solves a boundary value problem at the continuum scale by FEM and derives the material point response from the discrete element method simulation attached to each Gauss point of the FEM mesh. The multiscale modelling framework not only helps successfully bypass phenomenological constitutive assumptions as required in conventional modelling approaches but also facilitates effective cross‐scale interpretation and understanding of soil behaviour. We examine the classical retaining wall and footing problems by this method and demonstrate that the simulating results can be well validated and verified by their analytical solutions. Furthermore, the study sheds novel multiscale insights into these classical problems and offers a new tool for geotechnical engineers to design and analyse geotechnical applications based directly upon particle‐level information of soils. Copyright © 2015 John Wiley & Sons, Ltd.     

Frictionless contact formulation for dynamic analysis of nonlinear saturated porous media based on the mortar method

A finite element algorithm for frictionless contact problems in a two‐phase saturated porous medium, considering finite deformation and inertia effects, has been formulated and implemented in a finite element programme. The mechanical behaviour of the saturated porous medium is predicted using mixture theory, which models the dynamic advection of fluids through a fully saturated porous solid matrix. The resulting mixed formulation predicts all field variables including the solid displacement, pore fluid pressure and Darcy velocity of the pore fluid. The contact constraints arising from the requirement for continuity of the contact traction, as well as the fluid flow across the contact interface, are enforced using a penalty approach that is regularised with an augmented Lagrangian method. The contact formulation is based on a mortar segment‐to‐segment scheme that allows the interpolation functions of the contact elements to be of order N. The main thrust of this paper is therefore how to deal with contact interfaces in problems that involve both dynamics and consolidation and possibly large deformations of porous media. The numerical algorithm is first verified using several illustrative examples. This algorithm is then employed to solve a pipe‐seabed interaction problem, involving large deformations and dynamic effects, and the results of the analysis are also compared with those obtained using a node‐to‐segment contact algorithm. The results of this study indicate that the proposed method is able to solve the highly nonlinear problem of dynamic soil–structure interaction when coupled with pore water pressures and Darcy velocity. Copyright © 2015 John Wiley & Sons, Ltd.     

General coupling extended multiscale FEM for elasto‐plastic consolidation analysis of heterogeneous saturated porous media

This paper presents a general coupling extended multiscale FEM (GCEMs) for solving the coupling problem of elasto‐plastic consolidation of heterogeneous saturated porous media. In the GCEMs, the numerical multiscale base functions for the solid skeleton and fluid phase of the coupling system are all constructed on the basis of the equivalent stiffness matrix of the unit cell, which not only contain the interaction between the solid and fluid phases but also consider the time effect. Furthermore, in order to improve the computational accuracy for two‐dimensional problems, a multi‐node coarse element strategy for the GCEMs is proposed, and a two‐scale iteration algorithm for the elasto‐plastic consolidation analysis is developed. Some one‐dimensional and two‐dimensional homogeneous and heterogeneous numerical examples are carried out to validate the proposed method through the comparison with the coupling multiscale FEM and standard FEM. Numerical results show that the newly developed GCEMs can almost preserve the same convergent property as the standard FEM and also possesses the advantages of high computational efficiency. In addition, the GCEMs can be easily applied to other coupling multifield and multiphase transient problems. Copyright © 2014 John Wiley & Sons, Ltd.