A numerical model that accounts for fully coupled long-term large strain consolidation and heat transfer provides a more realistic analysis for various applications, including geothermal energy storage and extraction, buried power cables, waste disposal, groundwater tracers, and landfills. Despite its importance, existing models rarely simulate fully coupled large-strain long-term consolidation and heat transfer effectively. To address this research gap, this study presents a numerical model, called Consolidation and Heat Transfer (i.e., CHT), designed for one-dimensional (1D) coupled large-strain consolidation and heat transfer in layered soils, with the added capability to account for thermal creep. The model employs a piecewise-linear approach for the coupled long-term finite strain consolidation and heat transfer processes. The consolidation algorithm extends the functionality of the CS-EVP code by incorporating thermally induced strains. The heat transfer algorithm accounts for conduction, thermomechanical dispersion, and advection, assuming local thermal equilibrium between fluid and solid phases. Heat transfer is consistent with the spatial and temporal variation of void ratio and seepage velocity in the long-term consolidating layer. This paper details the development of the CHT model, presents verification checks against existing numerical solutions, and demonstrates its performance through several simulations. These simulations illustrate the effects of seepage velocity, thermal boundary conditions, and layered soil configurations on the coupled heat transfer and consolidation behavior of saturated compressible soils.