Fuel burnup analysis requires a high computational cost for full core calculations, due to the amount of the information processed for the total reaction rates in many burnup regions. Indeed, they reach the order of millions or more by a subdivision into radial and axial regions in a pin-by-pin description. In addition, if multi-physics approaches are adopted to consider the effects of temperature and density fields on fuel consumption, the computational load grows further. In this way, the need to find a compromise between computational cost and solution accuracy is a crucial issue in burnup analysis. To overcome this problem, the present work aims to develop a methodological approach to implement a Reduced Order Model (ROM), based on Proper Orthogonal Decomposition (POD), in fuel burnup analysis. We verify the approach on 4 years of burnup of the TMI-1 unit cell benchmark, by reconstructing fuel materials and burnup matrices over time with different levels of approximation. The results show that the modeling approach is able to reproduce reactivity and nuclide densities over time, where the accuracy increases with the number of basis functions employed.