Interfacial thermal contact model for consolidation of multilayered saturated soils subjected to time-dependent heating and loading

• Four thermal contact models were developed to predict thermal contact resistance at the interface of multilayered saturated soils. • Semi-analytical solutions for thermal consolidation under time-dependent heating and loading were derived using Laplace transform techniques. • The general thermal contact model produces the most significant thermal gradient at the interface, while the perfect contact model overestimates deformation. • Increasing interfacial thermal contact resistance diminishes the influence of temperature on consolidation behaviour.