A disk-shaped domain integral method for the computation of stress intensity factors using tetrahedral meshes

A novel domain integral approach is introduced for the accurate computation of pointwise J-integral and stress intensity factors (SIFs) of 3D planar cracks using tetrahedral elements. This method is efficient and easy to implement, and does not require a structured mesh around the crack front. The method relies on the construction of virtual disk-shaped integral domains at points along the crack front, and the computation of domain integrals using a series of virtual triangular and line elements. The accuracy of the numerical results computed for through-the-thickness, penny-shaped, and elliptical crack configurations has been validated by using the available analytical formulations. The average error of computed SIFs remains below 1% for fine meshes, and 2–3% for coarse ones. The results of an extensive parametric study suggest that there exists an optimum mesh-dependent domain radius at which the computed SIFs are the most accurate. Furthermore, the results provide evidence that tetrahedral elements are efficient, reliable and robust instruments for accurate linear elastic fracture mechanics calculations

An investigation of irregular crack path effects on fracture mechanics parameters using a grain microstructure meshing technique

Electronic version of an article published as Journal of Multiscale Modeling, Vol. 4, Iss. 1, atricle 1250001, 2012, http://dx.doi.org/10.1142/S1756973712500011 © World Scientific Publishing Company, http://www.worldscientific.com/worldscinet/jmmA sub-grain size finite element modelling approach is presented in this paper to investigate variations in fracture mechanics parameters for irregular crack paths. The results can be used when modelling intergranular and transgranular crack growth where creep and fatigue are the dominant failure mechanisms and their crack paths are irregular. A novel method for sub-grain scale finite element mesh consisting of multiple elements encased in ~50–150 μm-sized grains has been developed and implemented in a compact tension, C(T), mesh structure. The replicated shapes and dimensions were derived from an isotropic metallic grain structure using representative random sized grain shapes repeated in sequence ahead of the crack tip. In this way the effects of crack tip angle ahead of the main crack path can be considered in a more realistic manner. A comprehensive sensitivity analysis has been performed for elastic and elastic-plastic materials using ABAQUS and the stress distributions, the stress intensity factor and the J-integral have been evaluated for irregular crack paths and compared to those of obtained from analytical solutions. To examine the local and macroscopic graph path effects on fracture mechanics parameters, a few extreme cases with various crack-tip angles have been modelled by keeping the macroscopic crack path parallel to the axis of symmetry. The numerical solutions from these granular mesh structures have been found in relatively good agreement with analytical solutions

The proper generalized decomposition for the simulation of delamination using cohesive zone model

The use of cohesive zone models is an efficient way to treat the damage, especially when the crack path is known a priori. This is the case in the modeling of delamination in composite laminates. However, the simulations using cohesive zone models are expensive in a computational point of view. When using implicit time integration scheme or when solving static problems, the non-linearity related to the cohesive model requires many iterations before reaching convergence. In explicit approaches, the time step stability condition also requires an important number of iterations. In this article, a new approach based on a separated representation of the solution is proposed. The Proper Generalized Decomposition is used to build the solution. This technique, coupled with a cohesive zone model, allows a significant reduction of the computational cost. The results approximated with the PGD are very close to the ones obtained using the classical finite element approach

Crack propagation in honeycomb cellular materials: a computational approach

Computational models based on the finite element method and linear or nonlinear fracture mechanics are herein proposed to study the mechanical response of functionally designed cellular components. It is demonstrated that, via a suitable tailoring of the properties of interfaces present in the meso- and micro-structures, the tensile strength can be substantially increased as compared to that of a standard polycrystalline material. Moreover, numerical examples regarding the structural response of these components when subjected to loading conditions typical of cutting operations are provided. As a general trend, the occurrence of tortuous crack paths is highly favorable: stable crack propagation can be achieved in case of critical crack growth, whereas an increased fatigue life can be obtained for a sub-critical crack propagation

Node-to-segment and node-to-surface interface finite elements for fracture mechanics

The topologies of existing interface elements used to discretize cohesive cracks are such that they can be used to compute the relative displacements (displacement discontinuities) of two opposing segments (in 2D) or of two opposing facets (in 3D) belonging to the opposite crack faces and enforce the cohesive traction-separation relation. In the present work we propose a novel type of interface element for fracture mechanics sharing some analogies with the node-to-segment (in 2D) and with the node-to-surface (in 3D) contact elements. The displacement gap of a node belonging to the finite element discretization of one crack face with respect to its projected point on the opposite face is used to determine the cohesive tractions, the residual vector and its consistent linearization for an implicit solution scheme. The following advantages with respect to classical interface finite elements are demonstrated: (i) non-matching finite element discretizations of the opposite crack faces is possible; (ii) easy modelling of cohesive cracks with non-propagating crack tips; (iii) the internal rotational equilibrium of the interface element is assured. Detailed examples are provided to show the usefulness of the proposed approach in nonlinear fracture mechanics problems.Comment: 37 pages, 17 figure