Abstract
A computational approach for analyzing boundary value problem with granular materials is developed based on the two-scale homogenization method. While the microscopic problem leads to the variational inequality which reflects microscopic friction-contact responses, the homogenized structure reveals the nonlinear response stemming directly from the averaged microscopic behaviors. Although the macroscopic problem is analyzed by the continuum-based FEM, the microscopic one is remodelled by rigid grains with spring and friction devices. In addition to the investigation of the effects of grains' configurations in a unit cell, the applicability of this two-scale modeling is discussed via representative numerical examples. In particular, the bi-axial compression tests on a plane specimen is simulated to illustrate the feasibility of the proposed method.
