Analytical and numerical analysis for frost heaving stress distribution within rock joints under freezing and thawing cycles

Water-bearing joints within rock engineering in cold areas are often subjected to frost heaving force in cold season due to water–ice phase transition. To evaluate the damage and stability of rock mass in cold regions, a 3D model that considers moisture migration loss during freezing and thawing was established to study the characteristics of frost heaving force within joints. Then, the numerical simulation of cyclic freeze-thawing of water-bearing joints was carried out through equivalent expansion coefficient and particle flow calculation methods. The distribution of frost heaving force in and around the joints was obtained. According to the results of the numerical tests and theoretical calculations, the frost heaving force in joints is basically stable, the tensile stress concentration area appears at the joint tip, and the frost heaving force decreases gradually away from the jointed rock mass area. The frost heaving force decreases considerably with increasing cycle number and moisture migration loss but it increases with increasing mechanical strength and joint geometric size of rock and ice. The comparison between the numerical solution of the equivalent expansion coefficient method and the theoretical solution shows that the force size and distribution law of frost heaving for the two methods are consistent.