A generalized anisotropic deformation formulation for geomaterials

In this paper, the combined finite-discrete element method (FDEM) has been applied to analyze the deformation of anisotropic geomaterials. In the most general case geomaterials are both non-homogeneous and non-isotropic. With the aim of addressing anisotropic material problems, improved 2D FDEM formulations have been developed. These formulations feature the unified hypo-hyper elastic approach combined with a multiplicative decomposition-based selective integration for volumetric and shear deformation modes. This approach is significantly different from the co-rotational formulations typically encountered in finite element codes. Unlike the co-rotational formulation, the multiplicative decomposition-based formulation naturally decomposes deformation into translation, rotation, plastic stretches, elastic stretches, volumetric stretches, shear stretches, etc. This approach can be implemented for a whole family of finite elements from solids to shells and membranes. This novel 2D FDEM based material formulation was designed in such a way that the anisotropic properties of the solid can be specified in a cell by cell basis, therefore enabling the user to seed these anisotropic properties following any type of spatial variation, for example, following a curvilinear path. In addition, due to the selective integration, there are no problems with volumetric or shear locking with any type of finite element employed.

[1]  O. K. Mahabadi,et al.  Verification of the Implementation of Rock-Reinforcement Elements in Numerical Analyses Based on the Hybrid Combined Finitediscrete Element Method (FDEM) , 2015 .

[2]  A. Munjiza,et al.  Fracture and fragmentation of thin shells using the combined finite–discrete element method , 2013 .

[3]  Hari S. Viswanathan,et al.  Fracture-permeability behavior of shale , 2015 .

[4]  Joseph P. Morris,et al.  Dynamic simulations of geological materials using combined FEM/DEM/SPH analysis , 2008 .

[5]  O. K. Mahabadi,et al.  Numerical modelling of a Brazilian Disc test of layered rocks using the combined finite-discrete element method , 2009 .

[6]  Petr Krysl,et al.  Locking‐free continuum displacement finite elements with nodal integration , 2008 .

[7]  M. Ortiz,et al.  Computational modelling of impact damage in brittle materials , 1996 .

[8]  Joseph P. Morris,et al.  Simulations of Fracture and Fragmentation of Geologic Materials using Combined FEM/DEM Analysis , 2005 .

[9]  M. Gurtin,et al.  An introduction to continuum mechanics , 1981 .

[10]  Antonio Munjiza,et al.  Numerical comparison of some explicit time integration schemes used in DEM, FEM/DEM and molecular dynamics , 2004 .

[11]  Antonio Munjiza,et al.  MRCK_3D contact detonation algorithm , 2010 .

[12]  R. P. Swift,et al.  The Combined Finite-Discrete Element Method applied to the Study of Rock Fracturing Behavior in 3D , 2011 .

[13]  A. Munjiza The Combined Finite-Discrete Element Method , 2004 .

[14]  Davide Elmo,et al.  An Integrated Numerical Modelling–Discrete Fracture Network Approach Applied to the Characterisation of Rock Mass Strength of Naturally Fractured Pillars , 2010 .

[15]  Esteban Rougier,et al.  Computational Mechanics of Discontinua , 2011 .

[16]  Giovanni Grasselli,et al.  ROCK IMPACT MODELLING USING FEM/DEM , 2010 .

[17]  Chunning Ji,et al.  Large scale simulation of red blood cell aggregation in shear flows. , 2013, Journal of biomechanics.

[18]  Ted Belytschko,et al.  Triangular composite finite elements , 2000 .

[19]  A. Munjiza,et al.  Large Strain Finite Element Method: A Practical Course , 2015 .

[20]  A. Munjiza The Combined Finite-Discrete Element Method: Munjiza/Discrete Element Method , 2004 .

[21]  Antonio Munjiza,et al.  Validation of a three-dimensional Finite-Discrete Element Method using experimental results of the Split Hopkinson Pressure Bar test , 2014 .

[22]  Antonio Munjiza,et al.  Combined single and smeared crack model in combined finite-discrete element analysis , 1999 .