Bridging scale method for discrete particle assembly-Cosserat continuum modeling
Received:January 14, 2011  Revised:September 21, 2011
View Full Text  View/Add Comment  Download reader
DOI:10.7511/jslx20123013
KeyWord:granular material  bridging scale method  discrete element method  cosserat continuum  multi-scale  non-reflecting boundary condition
     
AuthorInstitution
万柯 大连理工大学 工业装备结构分析国家重点实验室, 大连
李锡夔 大连理工大学 工业装备结构分析国家重点实验室, 大连
Hits: 1980
Download times: 1195
Abstract:
      On the basis of the bridging scale method (BSM) initially proposed for molecular dynamics-Cauchy continuum modeling ,a new version of the BSM that couples the discrete particle assembly model using the discrete element method (DEM) and the Cosserat continuum model using the finite element method (FEM) in fine and coarse scales respectively is presented.The present BSM applies the DEM only to limited local regions of the whole computational domain for the purpose of accurate simulation of material failure with discontinuous deformation characteristics in microscopic scale,and meantime applies the FEM that costs much less both computational time and storage space to the whole domain.In addition,different time step sizes are allowed to the time integration schemes used to the coarse and fine scales respectively,as a consequence,both computational accuracy and efficiency of the present BSM is greatly enhanced.With the coarse and fine scale decomposition of translational and rotational displacements,in light of the principle of virtual work applied to the FEM nodes of Cosserat continuum and the particle centers of discrete particle assembly respectively,two decoupling sets of equations of motion of the combined coarse-fine scale system are derived.The interfacial condition between coarse and fine domains in the case of quasi-static loading,and the non-reflecting boundary condition (NRBC),which is capable of effectively eliminating spurious reflective waves at the interface between coarse and fine domains under dynamic loads are presented and discussed.The numerical results for 2D example problems illustrate the applicability and advantages of the present BSM,as well as the availability of the proposed interfacial condition for dynamical response simulations of granular materials.