Comprehensive analysis of mode-I cracking in ice: Exploring full-range rate dependency

Ice, as a crystalline material, demonstrates a unique response in its mode-I fracture toughness to variations in loading rates, characterized by an initial decrease and subsequent increase in toughness. This phenomenon has been explored in limited studies. In this work, an extensive numerical analysis of mode-I ice cracking is conducted by using the distinct lattice spring model (DLSM). Norton-Bailey Drucker-Prager DLSM (NB-DP-DLSM) is initially employed to investigate the classical explanations of creep and stress relaxation for the anomaly observed in ice’s fracture toughness at lower loading rates, but this approach does not successfully replicate the experimentally observed strain rate dependency. Then, two rate-dependent constitutive models are introduced to further examine the mode-I fracture mechanics of ice. Our numerical simulations of three-point bending tests show that the relationship between fracture toughness and strain rate at lower levels is more accurately captured by rate-dependent models. For higher strain rates, our numerical modeling of notched semi-circular bending tests indicates that a rate-independent constitutive model can replicate the loading rate dependency of ice’s mode-I fracture toughness. In conclusion, these observations suggest that a reverse-stage-like dynamic constitutive model for DLSM can potentially capture the full range of loading rate dependencies observed in ice.