A mesoscopic framework is developed to quantify the thermal-mechanical-chemical responses of polymer-bonded explosive (PBX) samples under impact loading. A mesoscopic reactive model is developed for the cyclotetramethylenetetranitramine (HMX) crystal, which incorporates nonlinear elasticity, crystal plasticity, and temperature-dependent chemical reaction. The proposed model was implemented in the finite element code ABAQUS by the user subroutine VUMAT. A series of three-dimensional mesoscale models were constructed and calculated under low-strength impact loading scenarios from 100 m/s to 600 m/s where only the first wave transit is studied. Crystal anisotropy and microstructural heterogeneity are responsible for the nonuniform stress field and fluctuations of the stress wave front. At a critical impact velocity (>= 300 m/s), a chemical reaction is triggered because the temperature contributed by the volumetric and plastic works is sufficiently high. Physical quantities, including stress, temperature, and extent of reaction, are homogenized from those across the microstructure at the mesoscale to compare with macroscale measurements, which will advance the continuum-level models. The framework presented in this study has important implications in understanding hot spot ignition processes and improving predictive capabilities in energetic materials. (C) 2016 Elsevier B.V. All rights reserved.