Non-human primate (NHP) spinal cord injury experiments exhibit high intersubject variability in biomechanical parameters even when a consistent impact protocol is applied to each subject. Optimizing impact parameters to reduce this variability through experiments is logistically challenging in NHP studies. Finite element models provide a complimentary tool to inform experimental design without the cost and complexity of live animal studies. A morphologically variable virtual population (N= 10) of NHPs quantified the interaction of morphological variability and different impact conditions in a unilateral cervical contusion, including impactor size (4 and 5 mm) and mediolateral alignment over the cord midline (0.5 and 1 mm). We explored the effect of these conditions on the magnitude and intersubject variability of impact force and cord lateral slippage. The study demonstrated that a 1-mm mediolateral alignment maximized peak forces and minimized lateral slippage. A 5-mm impactor was beneficial in increasing peak forces, whereas a 4-mm impactor reduced lateral slippage. Comparatively, intersubject variability in peak forces and lateral slippage were minimized with a 0.5-mm mediolateral alignment. The study highlights that impact parameters selected based on peak forces may not be beneficial in reducing variability. The study also showed that morphology was an important contributor to variability. Integrating morphology variability through a virtual population in an injury simulation to investigate mechanistic research questions will more effectively capture the heterogeneity of experiments and provide better insights for effective experimental design.