Phase field models have become an effective tool for predicting complex crack configurations including initiation, propagation, branching, intersecting and merging. However, several computational issues have hindered their utilisation in engineering practice, such as the convergence challenge in implicit algorithms, numerical stability issues in explicit methods and significant computational costs. Aiming to providing a more efficient numerical algorithm, this work integrates the explicit integral operator with the recently developed neighbored element method, for the first time, to solve the coupled governing equations in phase field models. In addition, the damage irreversibility can be ensured automatically, avoiding the need to introduce extra history variable for the maximum driving force in traditional algorithms. Six representative fracture benchmarks with different failure modes are simulated to verify the effectiveness of the proposed method, including the multiple cracks in heterogeneous concrete at mesoscale. It is found that this semi-explicit numerical algorithm yields consistent crack profiles and load capacities for all examples to the available experimental data and literature. In particular, the computational cost is significantly reduced when compared to the traditional explicit modelling. Therefore, the presented numerical algorithm is highly attractive and promising for phase-field simulations of complicated 3D solid fractures in structural-level engineering practices.
Read full abstract