Stacked packaging units is the main form of distribution packaging of products. Its dynamic performance is not fully understood. This paper investigated the influence of the constraint, input vibration, location and contact nonlinearity on the dynamic performance of three layers stacked packaging units. The dynamic contact force between surfaces and acceleration response of products were obtained. In sine sweep vibrations, the constraint to stacked packaging units has an obvious influence on the dynamic characteristics. The acceleration response of product is associated with the vibration mode. The force amplification factor is in general between 1.5 and 2, but it can close to 3 on top layer in the case of no fixed. In random vibrations, non‐Gaussian data of dynamic contact force appear when Gaussian data of input vibration pass through the stacked packaging units, resulting a Weibull distribution of force level‐crossing. The force level‐crossing diagram becomes more abrupt with the decreasing of input vibration level, smoother from the top contact surface to middle and bottom ones, and moves right and becomes smoother with the constraint strengthen. In the case of lower input level, Gaussian distribution of force level‐crossing appears. The force power spectral density (PSD) between bottom box and table is much larger than that between boxes, which is significantly influenced by the first resonance frequency. However, the acceleration PSD of product is significantly influenced by both the first and second resonance frequencies, and controlled by the vibration mode. It depends much on the value of input acceleration PSD around the resonance frequencies.