Efficient quasi-brittle fracture simulations of concrete at mesoscale using micro CT images and a localizing gradient damage model

• Published in: Computer methods in applied mechanics and engineering Vol. 400; p. 115559
• Main Authors: Huang, Yu-jie, Zhang, Hui, Zhou, Jin-jie, Xu, Shi-lang
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier B.V 01.10.2022; Elsevier BV
• Subjects: Algorithms; Brittle fracture; Composite materials; Computed tomography; Computed Tomography (CT); Computer simulation; Damage assessment; Damage localization; Equilibrium equations; Equivalence; Finite element method; Image-based modelling; Localizing gradient damage model; Mathematical models; Medical imaging; Mesoscale concrete; Mesoscale phenomena; Mortars (material); Quasi-brittle fracture; Shape functions; Stiffness matrix; Image-based modelling; Quasi-brittle fracture; Mesoscale concrete; Localizing gradient damage model; Computed Tomography (CT)
• ISSN: 0045-7825; 1879-2138
• DOI: 10.1016/j.cma.2022.115559

In this work, a localizing gradient damage model (LGDM) based on the generalized micromorphic theory is adopted to investigate the quasi-brittle fracture behaviour of concrete at mesoscale for the first time. Micro Computed Tomography (CT) image-based realistic concrete models are generated including aggregates, mortar, pores, and aggregate-mortar interfacial transition zones (ITZs). In the LGDM, an additional micro-force equilibrium equation is formulated with a damage-dependent interaction function to prevent spurious energetic interactions during the softening stage. The damage driving forces are evaluated as history variables through the maximum nonlocal equivalent strains at the integration points. The stiffness matrices and residual vectors for displacement and nonlocal equivalent strain fields are formulated through the conventional shape function and strain matrices of finite elements, by using the user-defined element (UEL) subroutine implemented with the Broyden–Fletcher–Goldfarb–Shanno (BFGS) algorithm to improve numerical performance. Typical nonlinear fracture benchmarks are presented to investigate the variability of mesoscale fracture evolution, crack trajectories, and load–displacement curves as well as numerical performance. It is found that the LGDM significantly improves the computational efficiency over the conventional phase-field cohesive zone model (PFCZM) with 60% savings of CPU time, and holds great potential for multiscale fracture evaluation of quasi-brittle composite materials, through the high-fidelity representation of real materials’ microstructures, and flexible simulation of the complicated nonlinear fracture without the need for re-meshing or mesh enrichments. •Simulate mesoscale fracture of concrete using a localizing gradient damage model (LGDM) for the first time.•Enable a robust and efficient monolithic framework by the BFGS quasi-Newton algorithm.•Successfully capture quasi-brittle fracture behaviour of mesoscale concrete using CT images with high fidelity.•Demonstrate considerable reduction (60%) in computational time over the phase-field cohesive zone model.•Widen the applicability and acceptability of the LGDM among researchers.