ABSTRACT
A numerical method based on a damage-viscoplasticity model for rock implemented with the polygonal finite element method is presented in this paper. As many rocks display polygonal mineral texture, polygonal elements are naturally suited for modelling these rocks. A heterogeneity description based on random clusters of elements representing the constituent minerals of the rock is included. The numerical simulations of dynamic indentation of rock with a spherical indenter are carried out to demonstrate the method. The effect of surface roughness is accounted for in the simulation showing that it adds to the reliability of the results.
