A thermodynamics-based mathematical model for the frost heave of unsaturated freezing soils