Pultrusion is a continuous and cost-effective process for a production of composite structural components with a constant cross-sectional area. In most cases the numerical models of pultrusion processes are firmly connected with the final product and applied technology, and therefore no more than one algorithm is utilized for a simulation purpose. For this reason in the present investigation, two different approaches for a numerical simulation of pultrusion processes are proposed, compared and discussed. Each of them is constructed by using the general-purpose FE software that results in considerable savings in development time and costs, and also makes available various modeling features of the FE package. The first procedure is developed in ANSYS Mechanical environment and based on the mixed time integration scheme and nodal control volumes method to decouple the coupled energy and species equations. The second procedure is performed by using ANSYS CFX software and presents quite new approach not used widely for a pultrusion modeling. The cure reaction in this case is introduced as an additional variable. The developed procedures have been validated by the results of other authors that additionally gave the possibility to compare algorithms in their applicability for an accurate thermo-chemical simulation of pultrusion processes. Good agreement between present finite element results and published numerical-experimental results has been observed for the temperature and degree of cure fields in all test problems studying the pultrusion of cylindrical rod, flat plate and I-beam profiles.