A multiscale coupling approach between discrete element method and finite difference method for dynamic analysis

Summary The coupling between two distinct numerical methods presents a major challenge, especially in the case of discrete-continuum coupling for dynamic simulations. This paper presents a general multiscale framework for coupling the discrete element method (DEM) and the finite difference method (FDM). DEM has been shown to be particular suitable for capturing physical phenomena at small length scales where continuum mechanics description no longer applies. Its efficiency, however, is limited because of the computational power required. A multidomain analysis coupling DEM and FDM is thus proposed to reduce the computational cost. To couple overlapping DEM and FDM, a bridging scale term is introduced such that compatibility of dynamic behavior between the DEM-based and the FDM-based models is enforced. This multiscale method couples two commercial packages: the DEM-based code, particle flow code, and the FDM-based code, fast lagrangian analysis of continua. The new method is applied to several reference dynamic tests. Results using the proposed method compare well with benchmark simulations from FDM and DEM. Copyright © 2015 John Wiley & Sons, Ltd.

[1]  P. Cundall,et al.  A discrete numerical model for granular assemblies , 1979 .

[2]  Eugenio Oñate,et al.  Multiscale analysis using a coupled discrete/finite element model , 2008 .

[3]  Laurent Daudeville,et al.  Multidomain finite and discrete elements method for impact analysis of a concrete structure , 2009 .

[4]  Jiannong Fang,et al.  A 3D distinct lattice spring model for elasticity and dynamic failure , 2011 .

[5]  Gaël Combe,et al.  Coupling between finite and discrete element methods for the modelling of earth structures reinforced by geosynthetic , 2009 .

[6]  Xiaojing Zheng,et al.  Multiscale mechanical behaviors in discrete materials: A review , 2010 .

[7]  G. Guiochon,et al.  Coupled discrete element and finite volume solution of two classical soil mechanics problems , 2011 .

[8]  Usama El Shamy,et al.  A practical co-simulation approach for multiscale analysis of geotechnical systems , 2010 .

[9]  J. D. Doll,et al.  Generalized Langevin equation approach for atom/solid-surface scattering: General formulation for classical scattering off harmonic solids , 1976 .

[10]  Ivan Iordanoff,et al.  3D coupling approach between discrete and continuum models for dynamic simulations (DEM–CNEM) , 2013 .

[11]  Ted Belytschko,et al.  Conservation properties of the bridging domain method for coupled molecular/continuum dynamics , 2008 .

[12]  David Saiang,et al.  Stability analysis of the blast-induced damage zone by continuum and coupled continuum–discontinuum methods , 2010 .

[13]  Alain Combescure,et al.  Three dimensional SPH-FEM gluing for simulation of fast impacts on concrete slabs , 2011 .

[14]  Esteban Samaniego,et al.  AES for multiscale localization modeling in granular media , 2011 .

[15]  Guillaume Rateau,et al.  The Arlequin method as a flexible engineering design tool , 2005 .

[16]  J. Andrade,et al.  On the rheology of dilative granular media: Bridging solid- and fluid-like behavior , 2012 .

[17]  Antonio Bobet,et al.  Multiscale method for long tunnels subjected to seismic loading , 2013 .

[18]  Gregory J. Wagner,et al.  Coupling of atomistic and continuum simulations using a bridging scale decomposition , 2003 .

[19]  Pierre-Etienne Gautier,et al.  Modelling ballast behaviour under dynamic loading. Part 1: A 2D polygonal discrete element method approach , 2006 .

[20]  N. Monteiro Azevedo,et al.  Hybrid discrete element/finite element method for fracture analysis , 2006 .

[21]  Xuxin Tu,et al.  Return mapping for nonsmooth and multiscale elastoplasticity , 2009 .

[22]  Mengyan Zang,et al.  An approach to combining 3D discrete and finite element methods based on penalty function method , 2010 .

[23]  Eugenio Oñate,et al.  Combination of discrete element and finite element methods for dynamic analysis of geomechanics problems , 2004 .

[24]  Zdeněk P. Bažant,et al.  Spurious reflection of elastic waves due to gradually changing finite element size , 1983 .

[25]  Hachmi Ben Dhia,et al.  Application of the Arlequin method to some structures with defects , 2002 .

[26]  Ted Belytschko,et al.  A bridging domain and strain computation method for coupled atomistic–continuum modelling of solids , 2007 .

[27]  A. Soroush,et al.  Fluid coupled-DEM modelling of undrained behavior of granular media , 2008 .

[28]  G. Rateau Méthode Arlequin pour les problèmes mécaniques multi-échelles : applications à des problèmes de jonction et de fissuration de structures élancées , 2003 .

[29]  Ludovic Chamoin,et al.  Ghost forces and spurious effects in atomic‐to‐continuum coupling methods by the Arlequin approach , 2010 .

[30]  J. Doll,et al.  Generalized Langevin equation approach for atom/solid–surface scattering: Numerical techniques for Gaussian generalized Langevin dynamics , 1976 .

[31]  Hiroshi Morioka,et al.  FLAC/PFC coupled numerical simulation of AE in large-scale underground excavations , 2007 .

[32]  Noam Bernstein,et al.  Spanning the continuum to quantum length scales in a dynamic simulation of brittle fracture , 1998 .

[33]  Jack W. Baker,et al.  Characterization of random fields and their impact on the mechanics of geosystems at multiple scales , 2012 .

[34]  Ted Belytschko,et al.  Bridging domain methods for coupled atomistic–continuum models with L2 or H1 couplings , 2009 .

[35]  D. Owen,et al.  A combined finite‐discrete element method in transient dynamics of fracturing solids , 1995 .

[36]  Harold S. Park,et al.  Non‐reflecting boundary conditions for atomistic, continuum and coupled atomistic/continuum simulations , 2005 .

[37]  Jianqiao Ye,et al.  Dynamic simulation of crack initiation and propagation in cross-ply laminates by DEM , 2011 .

[38]  Ted Belytschko,et al.  On the L2 and the H1 couplings for an overlapping domain decomposition method using Lagrange multipliers , 2007 .

[39]  Ted Belytschko,et al.  A continuum‐to‐atomistic bridging domain method for composite lattices , 2010 .

[40]  Mengyan Zang,et al.  A contact algorithm for 3D discrete and finite element contact problems based on penalty function method , 2011 .

[41]  T. Belytschko,et al.  A bridging domain method for coupling continua with molecular dynamics , 2004 .

[42]  P. Tallec Domain decomposition methods in computational mechanics , 1994 .

[43]  J. Q. Broughton,et al.  Concurrent coupling of length scales: Methodology and application , 1999 .