Abstract: Accurate prediction of thermal distribution and thermal damage inside skin tissue subjected to laser irradiation is a difficult problem in the biomedical field. To explore this problem, a strategy to derive an analytical solution using Laplace and Hankel transforms was proposed. A three-dimensional non-Fourier heat conduction model was established, and temperature-rise function was introduced to simplify the problem. The analytical solution was derived through integral transforms and their inversions. The prediction results of skin thermal distribution and thermal damage under different irradiation durations and radii were obtained. And the effects of laser irradiation duration and beam radius on skin tissue were researched. Furthermore, a comprehensive comparison was made between the skin tissue damage prediction outcomes achieved through non-Fourier heat conduction models and through classical Fourier heat conduction models. The analytical results indicate distinct differences in predicted thermal damage depth between the non-Fourier and classical Fourier heat conduction models. Specifically, under short-duration irradiation, the non-Fourier model predicts shallower damage depth than the Fourier model, while the opposite trend is observed under prolonged irradiation.