Abstract
The evolution of the microstructure of an assembly of cohesionless granular materials with associated pores, which carry the overall applied stresses through frictional contacts is a complex phenomenon. The macroscopic flow of such materials take place by the virtue of the relative rolling and sliding of the grains on the microβscale. A new discrete element method for biaxial compression simulations of random assemblies of oval particles with mixed sizes is introduced. During the course of deformation, the new positions of the grains are determined by employing the static equilibrium equations. A key aspect of the method is that, it is formulated for ellipse crossβsectional particles, hence desirable inherent anisotropies are possible. A robust algorithm for the determination of the contact points between neighbouring grains is given. Employing the present methodology, many aspects of the behaviour of twoβdimensional assemblies of oval crossβsectional rods have been successfully addressed. The effects of initial void ratio, interparticle friction angle, aspect ratio, and bedding angle on the rolling and sliding contacts are examined. The distribution of normals to the rolling and sliding contacts have different patterns and are concentrated along directions, which are approximately perpendicular to one another. On the other hand, the distribution of all contact normals (combined rolling and sliding) are close to that of rolling contacts, which confirm that rolling is the dominant mechanism. This phenomenon becomes more pronounced for higher intergranular friction angle. Characteristics of the rolling and sliding contacts are also discussed in the context of the force angle, which is the inclination of contact force with respect to the contact normal. Copyright Β© 2003 John Wiley & Sons, Ltd.