Frost heave modelling using porosity rate function

Frost-susceptible soils are characterized by their sensitivity to freezing that is manifested in heaving of the ground surface. While significant contributions to explaining the nature of frost heave in soils were published in late 1920s, modelling efforts did not start until decades later. Several models describing the heaving process have been developed in the past, but none of them has been generally accepted as a tool in engineering applications. The approach explored in this paper is based on the concept of the porosity rate function dependent on two primary material parameters: the maximum rate, and the temperature at which the maximum rate occurs. The porosity rate is indicative of ice growth, and this growth is also dependent on the temperature gradient and the stress state in the freezing soil. The advantage of this approach over earlier models stems from a formulation consistent with continuum mechanics that makes it possible to generalize the model to arbitrary three-dimensional processes, and use the standard numerical techniques in solving boundary value problems. The physical premise for the model is discussed first, and the development of the constitutive model is outlined. The model is implemented in a 2-D finite element code, and the porosity rate function is calibrated and validated. Effectiveness of the model is then illustrated in an example of freezing of a vertical cut in frost-susceptible soil. Copyright © 2006 John Wiley & Sons, Ltd.