Research has been undertaken to investigate the deformation properties of cigarette packages subjected to different loads. Numerical simulation enables a hypothetical package produced from a hypothetical paperboard to be tested. A finite element model was implemented in the Abaqus system in order to calculate stresses and strains in the deformed packages. The model behaves as expected up to deformations with limited damage. The calculated initial slopes of the reaction force responses were generally in good agreement with the corresponding measured slopes, although the model had a nonlinear force response that was not found in the experimental data. Incorporating cohesive behaviour into the model would allow the elements to disconnect from each other. Statistical analysis of the couplings between the initial process responses and the material properties of the paperboards revealed that only the initial yield stress parameter significantly affected the slope.