Crystallization in certain polymers, like natural rubber, is characterized by the specific geometrical arrangement of atoms in macromolecules caused by high strains. Starting from crystallization nuclei, polymer chains leave their natural entangled structure, stretch out, fold back and stack. Eventually, they build regions with a regular structure, also called lamellae. The process must be taken into consideration when planning manufacturing processes since it significantly influences mechanical and thermal properties of the final product. The present contribution deals with the thermomechanical model for crystallization of unfilled polymers, which involves displacements and temperature as global degrees of freedom, and the degree of network regularity as an internal variable. The mechanical part of the model uses the dissipation potential with two special features: Firstly, the thermodynamically consistent framework is developed to simulate the reduction of the network regularity during the unloading phase. Secondly, the microstructure evolution under the cyclic tensile load is visualized. The thermal part of the model is based on the solution of the heat equation. The resulting, coupled thermomechanical problem is solved in a monolithic way. Finally, selected numerical examples are compared with experimental data of natural rubber without fillers.