In order to study the spherical aberration of thermal lens when the laser crystal is pumped with the pump light distributed differently and the pump light under the end-pumping condition, in this paper we establish a single-ended pump and constant temperature thermal model to analyze the working characteristics of the Nd:YVO4 crystal. The steady state heat conduction equation and Zernike polynomials are solved, and the relationship between thermal spherical aberration and distribution of pump laser is established. The model is used in simulation, and the simulation results are further analyzed theoretically, showing that under the same pump power, the spherical aberration is greatest when the pump beam is of 2-order super-Gaussian distribution. The spherical aberration decreases with the increase of pump distribution coefficient k (except the Gaussian distribution). With the increase of pump power, the influence of pump distribution coefficient k on spherical aberration is aggravated gradually, and the difference in spherical aberration caused by different values of distribution coefficient k increases gradually. The range of strongest laser power of the 2-order super-Gaussian distribution pump is analyzed and simulated. Under same pump power, the maximum range of the strongest laser power of 2-order super-Gaussian distribution pump is 0.3-0.63 times the Gauss radius. The research methods and conclusions obtained in this paper have universal applicability and can be used for quantitatively analyzing the temperature distributions, thermal deformations, optical path difference distributions, and spherical aberration distributions of other laser crystals. At the same time, this study also provides a theoretical reference for improving spherical aberration from the perspective of changing the distribution of pump light and the laser output characteristics.