Abstract: Based on Latin Hypercube Sampling (LHS) and the Spearman rank correlation, a parameter sensitivity analysis method for constitutive model is established. Using this method, we implement the holistic parameter sensitivity analysis for viscoelastic damage constitutive model of PBX. Taking the compression of a cylindrical PBX specimen as example, we analyze the parameter sensitivity of viscoelastic damage constitutive model of PBX. The sensitivities of all constitutive parameters to end-face loads and their orders are determined. The impacts of the sample sizes and the loading rate on the sensitivity of constitutive parameters are discussed. The comparison with results reported in existing literatures shows that the method presented in this paper is feasible for the sensitivity analysis of multi-parameters. It can be used to predict the sensitivity of parameters to experimental results in advance, and provides a theoretical basis for reducing parameters to be identified and determining appropriate simulation experiments.