To date, the SFEM has only been investigated for bilinear and Wachspress approximations and is limited to linear reproducing conditions. The goal of this paper is to extend the strain smoothing to higher order elements and to investigate numerically in which condition strain smoothing is beneficial to accuracy and convergence of enriched finite element approximations. We focus on three widely used enrichment schemes, namely: (a) weak discontinuities; (b) strong discontinuities; (c) near-tip linear elastic fracture mechanics functions. The main conclusion is that strain smoothing in enriched approximation is only beneficial when the enrichment functions are polynomial (cases (a) and (b)), but that non-polynomial enrichment of type (c) lead to inferior methods compared to the standard enriched FEM (e.g. XFEM).

No abstract

No abstract

Clayey rocks have a complex microstructure with multiple characteristic lengths. Deformation under mechanical loading generally induces damage by microcracking, which essentially concerns the scale of mineral inclusions embedded in the clay matrix. The modelling of these materials is considered within the framework of a double scale approach, by numerical homogenisation, of the squared finite element method type. This allows a heterogeneous microstructure of the material to be taken into account and a distribution of morphological properties to be introduced. Emphasis is placed on the generation of microstructures satisfying experimental observations, and keeping a certain simplicity to fit into the framework of double scale modelling. The material characteristics and behaviour are defined at the grain scale: the mineralogical properties include the mineral phase proportions and the grain morphology, while the material constituents are represented by elastic grains separated by damageable cohesive crack models. Then, the overall microscale behaviour of the material under solicitation is derived from equilibrated elementary area (EA) configuration and computational homogenisation. The variability of the material response is studied with regard to small-scale aspects as microstructure variability, microstructure size, grain angularity, and properties of grain contacts. Deformation anal-

Summary Double‐scale numerical methods constitute an effective tool for simultaneously representing the complex nature of geomaterials and treating real‐scale engineering problems such as a tunnel excavation or a pressuremetre at a reasonable numerical cost. This paper presents an approach coupling discrete elements (DEM) at the microscale with finite elements (FEM) at the macroscale. In this approach, a DEM‐based numerical constitutive law is embedded into a standard FEM formulation. In this regard, an exhaustive discussion is presented on how a 2D/3D granular assembly can be used to generate, step by step along the overall computation process, a consistent Numerically Homogenised Law. The paper also focuses on some recent developments including a comprehensive discussion of the efficiency of Newton‐like operators, the introduction of a regularisation technique at the macroscale by means of a second gradient framework, and the development of parallelisation techniques to alleviate the computational cost of the proposed approach. Some real‐scale problems taking into account the material spatial variability are illustrated, proving the numerical efficiency of the proposed approach and the benefit of a particle‐based strategy.

The article presents a numerical finite element study of fluid leakage in concrete. Concrete cracking is numerically modelled in the framework of a macroscopic probabilistic approach. Material heterogeneity and the re- lated mechanical effects are taken into account by defining the elementary mechanical properties according to spatially uncorrelated random fields. Each finite element is consid- ered as representative of a volume of heterogeneous ma- terial, whose mechanical behaviour depends on its own volume. The parameters of the statistical distributions defining the elementary mechanical properties thus vary over the computational mesh element-by-element. A weak hydro-mechanical coupling assumption is introduced to represent the influence of cracking on the variation of transfer properties: it is assumed that the mechanical cracking of a finite element induces a loss of isotropy of its own permeability tensor. At the elementary level, an ex- perimentally enhanced parallel plates model is used to re- late the local crack permeability to the elementary crack aperture. A Monte Carlo-like approach allows to statisti- cally validate the numerical method. The self-consistency of the proposed modelling strategy is finally explored through the numerical simulation of the hydro-mechanical splitting test, recently proposed by authors to evaluate the real-time evolution of the transfer properties of a concrete sample under loading

SUMMARYIn this work, the consequences of using several different discrete element granular assemblies for the representation of the microscale structure, in the framework of multiscale modeling, have been investigated. The adopted modeling approach couples, through computational homogenization, a macroscale continuum with microscale discrete simulations. Several granular assemblies were used depending on the location in the macroscale finite element mesh. The different assemblies were prepared independently as being representative of the same material, but their geometrical differences imply slight differences in their response to mechanical loading. The role played by the micro-assemblies, with weaker macroscopic mechanical properties, on the initiation of strain localization in biaxial compression tests is demonstrated and illustrated by numerical modeling of different macroscale configurations.

Summary This paper presents a multiscale model based on a FEM×DEM approach, a method that couples discrete elements at the microscale and finite elements at the macroscale. FEM×DEM has proven to be an effective way to treat real‐scale engineering problems by embedding constitutive laws numerically obtained using discrete elements into a standard finite element framework. This proposed paper focuses on some numerical open issues of the method. Given the nonlinearity of the problem, Newton's method is required. The standard full Newton method is modified by adopting operators different from the consistent tangent matrix and by developing ad‐hoc solution strategies. The efficiency of several existing operators is compared, and a new and original strategy is proposed, which is shown to be numerically more efficient than the existing propositions. Furthermore, a shared memory parallelization framework using OpenMP directives is introduced. The combination of these enhancements allows to overcome the FEM×DEM computational limitations, thus making the approach competitive with classical FEM in terms of stability and computational cost.