Comparison between a μFE model and DEM for an assembly of spheres under triaxial compression

This paper presents a simple case of a Face Centred Cubic (FCC) array of 2,000 spheres under triaxial compression to compare the results obtained using the Discrete Element Method (DEM) and a micro finite element model (μFE). This μFE approach was developed so that the internal structure of the soil can be obtained using x-ray computed tomography and converted into a numerical fabric. The individual grains are represented as continuum deformable bodies and the inter-granular interaction based on the defined contact laws. In order to demonstrate the simple contact constitutive behaviour used in this μFE model, the response for two contacting elastic spheres is compared with theoretical equations. The strength at failure of the packing of 2,000 spheres is seen to yield similar values for DEM, μFE and the analytical solution. When comparing the evolving void ratio, a good agreement between the two numerical models was observed for very small strains but as the strain increases, the values start to diverge, which is believed to be related with the rigidity of the grains used in DEM.