This study involves the development of the auxiliary stress approach for producing elastically-homogeneous lattice models of damage in geomaterials. The lattice models are based on random, three-dimensional assemblages of rigid-body-spring elements. Unlike conventional lattice or particle models, the elastic constants of a material (e.g., Young's modulus and Poisson's ratio) are represented properly in both global and local senses, without any need for calibration. The proposed approach is demonstrated and validated through analyses of homogeneous and heterogeneous systems under uni- and tri-axial loading conditions. Comparisons are made with analytical solutions and finite element results. Thereafter, the model is used to simulate a series of standard laboratory tests: (a) split-cylinder tests, and (b) uniaxial compressive tests of sedimentary rocks at the Horonobe Underground Research Laboratory in Hokkaido, Japan. Model inputs are based on physical quantities measured in the experiments. The simulation results agree well with the experimental results in terms of pre-peak stress-strain/displacement responses, strength measurements, and failure patterns.
