In this research, we introduce the Finite Pointset Method (FPM) as an innovative approach for approximating transient linear bio-thermoelasticity in skin tissues, addressing the complex interplay between thermal and elastic effects under three distinct shock conditions. Our study demonstrates significant potential for clinical applications, providing a framework for precise measurement and optimization of therapeutic parameters. Through a comprehensive parametric analysis, we explore temperature distribution and displacement profiles within tissue, emphasizing the impact of coupling parameters and perfusion rates. Notably, our findings reveal the cooling effects induced by perfusion under both coupled and uncoupled scenarios and highlight the variability of thermal responses based on the coupling parameter. This research contributes to the interdisciplinary field of biomedical engineering, offering novel insights into haemodynamics problems using FPM and paving the way for advancements in medical technology and therapeutic strategies.