A new numerical true triaxial engine based on discrete element method accounting for rolling resistance contact is developed. By this engine, we have simulated mechanical behavior of granular materials under complex stress loading path in this study. Stress-strain responses of a kind of typical granular sand under several stress loading path in meridian and deviatoric stress space are provided. The results show that the three dimensional effects like the intermediate principal stress play an important role in the modeling processes. Theoretical analysis in strength characteristic implies the strength criteria with three parameters such as unified strength criterion and van Eekelen strength criterion are capable of describing cohesionless granular material behaviors in three dimensional stress states. Moreover, the case study for Chende sand further demonstrates the numerical true triaxial engine, is a potential tool. As compared to conventional triaxial compression test, this new developed apparatus could be widely used to "measure" elastic-plastic behavior in three dimensional stress space for finite element analysis in geotechnical problems. (C) 2019 The Society of Powder Technology Japan. Published by Elsevier B.V. and The Society of Powder Technology Japan. All rights reserved.