Abstract
A method of multiscale analysis based on mathematical homogenization theory has been developed for quasi-static equilibrium problems of granular media. The micro-scale problem is analyzed by a discrete numerical model assuming elastic and frictional contact between rigid particles. This two-scale analysis enables us to obtain the macroscopic/phenomenological inelastic deformation response of a representative volume element (RVE). To examine the macroscopic deformation properties of the assumed RVE, a series of numerical experiments involving pure rotation of the principal stress axes are carried out. The necessity of incorporating the non-coaxiality induced by the tangent effect and the anisotropy in the yield condition is revealed in the phenomenological constitutive description of the deformation under principal stress axes rotation.
