Molecular mechanics (MM) and molecular dynamics (MD) simulation method were applied to explore the impact of temperature (220-380 K) on the thermostability, sensitivity, and mechanical performance of RDX (1,3,5-trinitro-1,3,5-triazacyco-hexane)/HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazocane) energetic cocrystal and mixture models. The mechanical property, the maximum trigger bond length (LN-NO2), binding energy, and cohesive energy density (CED) of the pure RDX, beta-HMX crystal, the cocrystal, and mixture models were acquired and compared. The results manifest that temperature has an important impact on the binding capacity between the components of the cocrystal and mixture. The binding energies decrease as the temperature rises, and the cocrystal has larger values than those of mixture. For all the models, theLN-NO2 increases and the CEDs decrease with the rising temperature, implying that the sensitivity of the explosives increases, while the LN-NO2 values of the cocrystal are smaller than those of HMX and the CED values are between those of RDX and beta-HMX, indicating that the sensitivity has been enhanced through co-crystallization. As the temperature increases, the shear modulus (G), bulk modulus (K), and tensile modulus (E) values of all models have an evident downtrend. Simultaneously, G, K, and E values of the cocrystal model are less than those of RDX and beta-HMX, while the K/G ratio and Cauchy pressure (C-12-C-44) are larger, signifying that co-crystallization can weaken the brittleness and enhance the ductility of the pure crystals. Compared with the mixture, the cocrystal has better ductility and stability.