Adaptive phase field analysis with dual hierarchical meshes for brittle fracture. (September 2019)

• Record Type: Journal Article
• Authors: Goswami, Somdatta; Anitescu, Cosmin; Rabczuk, Timon
• Is Part Of: Engineering fracture mechanics. Volume 218(2019)
• Publication Date: 2019-09
• Subjects: Brittle fracture -- Phase field -- Dual mesh -- PHT-splines -- Adaptive; Fracture mechanics -- Periodicals; Rupture, Mécanique de la -- Périodiques; Fracture mechanics; Periodicals; 620.112605
• DOI: 10.1016/j.engfracmech.2019.106608
• Languages: English
• ISSNs: 0013-7944

Abstract

Highlights: An adaptive phase field model based on PHT-splines is proposed. Within this adaptive framework, a new and efficient data transfer algorithm is proposed. The mechanical mesh is refinement driven by a recovery based error estimator. The efficiency is demonstrated through several numerical examples. Abstract: We present a novel dual-mesh based adaptive phase field method for solving fracture problems. As the name suggests, the proposed scheme is solved using two meshes with different characteristic element sizes ( h ); a coarser mesh for the elastic field and a finer mesh for the phase field. To facilitate the exchange of information between the meshes, an efficient data transfer algorithm is proposed. The efficiency of the formulation is enhanced by implementing an adaptive h -refinement scheme using polynomial splines over hierarchical T-meshes (PHT-splines). The inherent hierarchical nature of PHT-Splines andtheir localization property is exploited in the proposed data transfer algorithm as well as in implementing adaptivity. Independent refinement strategies have been proposed for both meshes. To illustrate the performance of the proposed approach, six numerical examples have been presented. For most examples, the critical load and computational efficiency (in terms of CPU time) have been reported for different combinations of element sizes of both meshes. The present work concludes that the proposed approach with h = 0.5 l 0 for the elastic mesh and h = 0.25 l 0 for the phase field mesh yields results as accurate as results obtained from the conventional phase field model with h = 0.25 l 0, where l 0 denotes length scale parameter. Moreover, the elastic analysis being performed on a coarser mesh makes the proposed approach faster and cheaper than the conventional approach, indicating potential for future applications in other phase field problems.