The collision of two elastic or viscoelastic spherical shells is investigated as a model for the dynamic response of a human head impacted by another head or by some spherical object. Determination of the impact force that is actually being transmitted to bone will require the model for the shock interaction of the impactor and human head. This model is indended to be used in simulating crash scenarios in frontal impacts, and provide an effective tool to estimate the severity of effect on the human head and to estimate brain injury risks. The model developed here suggests that after the moment of impact quasi-longitudinal and quasi-transverse shock waves are generated, which then propagate along the spherical shells. The solution behind the wave fronts is constructed with the help of the theory of discontinuities. It is assumed that the viscoelastic features of the shells are exhibited only in the contact domain, while the remaining parts retain their elastic properties. In this case, the contact spot is assumed to be a plane disk with constant radius, and the viscoelastic features of the shells are described by the fractional derivative standard linear solid model. In the case under consideration, the governing differential equations are solved analytically by the Laplace transform technique. It is shown that the fractional parameter of the fractional derivative model plays very important role, since its variation allows one to take into account the age-related changes in the mechanical properties of bone.