To understand the physiology of tissue burns for successful clinical treatment, it is important to investigate the thermal behavior of human skin tissue subjected to heat injury. In this paper, a fundamental solution-based hybrid¯nite element formulation is proposed for numerically simulating steady-state temperature distribution inside a multilayer human skin tissue during burning. In the present approach, since only element boundary integrals are involved, the computational dimension is reduced by one as the fundamental solutions used analytically satis¯es the bioheat governing equation. Further, in multi-layer skin modeling, the burn is applied via a heating disk at constant temperature on a part of the epidermal surface of the skin tissue. Numerical results from the proposed approach are¯rstly veri¯ed by comparing them with exact solutions of a simple single-layered model or the results from conventional¯nite element method. Thereafter, a sensitivity analysis is carried out to reveal the e®ect of biological and environmental parameters on temperature distribution inside the skin tissue subjected to heat injury.