Stochastic analysis of dynamic fracture of concrete using CT-image based mesoscale models with a rate-dependent phase field method

Concrete structures are commonly exposed to dynamic loads spanning a wide range of strain rates, and the inherent mesoscale heterogeneities complicate stochastic dynamic fracture mechanisms even more. This work develops a numerical framework using mesoscale concrete models based on micro computed tomography (CT) images to investigate such mechanisms with meaningful stochastic analyses. A rate-dependent phase field model is proposed to characterise the dynamic initiation and propagation of cracks by incorporating both micro-viscosity and macroscopic viscoelasticity, which is described by two standard Maxwell elements with different relaxation times to consider a wide range of strain rates. Moreover, the viscoelastic constitutive relation is formulated in the full strain space, which allows for a spectral decomposition of the strain tensor to determine the effective damage driving force, thus effectively addressing the issue of compressive fracture. A numerical implementation scheme is developed by combining user-defined element and material subroutines in ABAQUS/Explicit solver. Extensive Monte Carlo simulations of dynamic tension up to a strain rate of 200 s⁻¹ are performed with statistical analyses. This work reveals the intricate dynamics associated with mesoscale heterogeneities and identifies the critical transition state at 20 s⁻¹. The transition is characterised by changing modes of fracture patterns, stress wave propagation, and load-carrying capacities. A new TDIF–strain rate–standard deviation relation is also proposed and aligns well with the increasing dispersion of experimental data. The relationship between void content and tensile strength reflects the formation characteristics of crack networks, with the void content exhibiting a positive correlation with the TDIF from 20 s⁻¹ to 100 s⁻¹.