Efficient prediction of crack initiation from arbitrary 2D notches. (June 2022)

• Record Type: Journal Article
• Authors: Rettl, Matthias; Pletz, Martin; Schuecker, Clara
• Is Part Of: Theoretical and applied fracture mechanics. Volume 119(2022)
• Publication Date: 2022-06
• Subjects: Coupled criterion -- Finite element analysis -- Fracture mechanics -- Meta-modeling; Fracture mechanics -- Periodicals; 620.1126
• DOI: 10.1016/j.tafmec.2022.103376
• Languages: English
• ISSNs: 0167-8442

Abstract

Highlights: Predicting crack initiation from notches using the Coupled Criterion. Efficient implementation utilizing meta-modeling, scaling laws, and superposition. Energy release rate is computed using boundary relaxation in a submodel-like model. Computing one position takes 0.66 s (this method) instead of 2:50 min (full model). The method is easily applicable to arbitrary notches and 2D load cases. Abstract: An efficient two-scale approach for predicting mode I crack initiation from 2D notches based on the Coupled Criterion is proposed. On the scale of the local model, a voxel model containing the notch simulates the displacement field. The crack model is introduced on the smaller scale and is defined in an image space. Based on the notch curvature, the precomputed crack model can be transformed to any position on the notch surface. The displacement field of the local model is fitted at the boundaries of the transformed crack model by predefined deformation modes and results can be obtained by a superposition of precomputed crack model results. By introducing the crack in the crack model, the stiffness of this model is reduced and thus, the incremental energy release rate can be inaccurate. Therefore, a boundary relaxation approach is used to obtain more accurate energy release rates. It is shown that the method is very efficient as it requires only 3:20 min to analyze 50 positions on a notch compared to 2:21 h of a conventional approach using full FEM simulations. Thereby, the method is reliable in identifying the critical position. The predicted failure index at this position deviates by at most 10.8%. Since the crack model limits the length of initiating cracks, Irwin's length K I c 2 / σ c 2 of the material must lie below 2.53 times the radius of a circular hole under uniaxial tension. For a brittle material like A l 2 O 3, notches with a curvature radius above 31 μ m can thus be analyzed.