Cosmic ray research field has entered a brave new world of precise direct measurements in the multi-TeV energy domain thanks to the success of the newest space-borne instruments. DAMPE is a satellite astroparticle detector aimed at probing cosmic ray electrons and gamma rays from few GeV to 10 TeV and cosmic ray nuclei between 10 GeV and few hundred TeV with un- precedented energy resolution. It was successfully launched in December 2015, providing, in particular, the first direct observation of a break in a cosmic ray electron plus positron spectrum around 0.9 TeV. DAMPE presents a unique opportunity of directly measuring the cosmic ray spectrum and composition close to the knee region. However, one of the most challenging limita- tions on the precision of such measurement arises from the simulation of hadronic interactions in the detector. In spite of the utmost thick and finely segmented calorimeter onboard DAMPE, the cosmic ray proton and nuclei energy measurements still rely significantly on the hadronic Monte Carlo simulation. In this contribution, we discuss the challenges of the hadronic Monte Carlo detector simulation for TeV--PeV particles and review the existing solutions. We also present an in-house software tool which allows connecting the high-energy hadronic event generators, used in CORSIKA and implemented in the CRMC package, to the GEANT4 tool kit. Finally, we dis- cuss preliminary simulation results with DAMPE setup based on this tool. Comparison with the alternative FLUKA simulation is presented as well.