To improve the modeling of fiber reinforced composites, we present in this work numerical methods able to compute both fiber-reinforced composites deformation in squeeze flow and thermal-kinetic evolution. The homogeneous orthotropic model for fiber composites described the anisotropy of the in-plane fiber. Both physics are related since the kinetic evolution as well as the temperature profile modify the rheology of the composites. The model used here to account for the kinetics of reaction is the Bailleul's model. An study case is presented, where the mold temperature is set to 150C to a composite sample at 40C. Thermal transfer begins as well as sample compression at constant speed. The thermo-rheological-kinetical coupling is done by means of the viscosity. Here the kinetic impacts is shown by the evolution of the temperature, velocity and stress during reaction. Reaction begins in the surface and moves towards the center. However, reaction in the core of the material, once started, is much quicker than in the surface. This translates into a mapping of viscosity values presented during the reaction, modifying the usual nature of the squeeze flow pattern.