This study involves the development of the auxiliary stress approach for producing elasticallyhomogeneous 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 uniand 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) splitcylinder 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.