A numerical method for the simulation of heat transfer occurring in thermoplastic composites thermostamping process is proposed. A reduced thermal model, named additive decomposition, is developed. It is based on the operator splitting method under thin shell assumption. A resolution algorithm using this decomposition is proposed, and developed in MATLAB. The approach is validated by comparing solutions obtained with a full 3D resolution and the presented method. Using this method, the computational time is proved to be about over 30 times faster. Eventually, prediction of temperature field is a prerequisite for the prediction of other phenomena, such as crystallization kinetics. Finally, the proposed method is implemented in the simulation software for thermostamping process Plasfib.