Explicit FEA Cracks

From OSUPDOCS
Jump to navigation Jump to search

FEA calculations can model explicit cracks.

Explicit Cracks

Feacrack.jpg

Cracks in FEA calculations are set up by configuring the mesh to include a crack. A crack can be internal to the mesh or, for symmetric loading, can be on an edge of the mesh and defined with the help of boundary conditions. These two crack geometries are illustrated on the right. In the first form, the two element edges to the right of the crack tip are initially aligned on the same coordinates. When the crack is loaded in tension, the edges separate and the point where they are joined is the crack tip. The static FEA run in NairnFEA does not include contact and therefore only opened cracks can be fully modeled.

The second form of crack is for symmetric loading. The crack plane is on the bottom of the object and the displacement boundary conditions determine the location of the crack tip

Crack Closure Calculations

In crack closure calculations, nodal forces and displacements are used to find mode I and mode II energy release rate for crack growth. The calculations are done after the analysis and NairnFEAMPM is configured to do the calculations. A recent paper on crack closure calculations for a variety of elements and allowing for crack surface tractions is in Nairn (2011).[1] This reference was the basis for the NairnFEAMPM calculations. It pointed out that for best results, the crack plane edges in front of and behind the crack tip must initially be collinear and of be of equal length (Δa).

Crack closure calculations can be done for any element type allowed in NairnFEA. In addition, one special crack tip element available and it is configured by using the CrackTip command.

Crack Tip Command

Qtrpt.jpg

The CrackTip command coverts all elements at the crack tip to be quarter point elements. Quarter point elements are 6 or 8-node quadratic elements where the midside nodes near the crack tip are 1/4 of the way from the tip node to the next edge node rather than half way as in conventional quadratic elements (see figure on the right). Quarter point elements at the crack tip are created by first creating a mesh with 8-node quadrilaterals or 6-node triangles and then adding the scripted command:

   CrackTip (keyid)

or

   CrackTip (x),(y)

In XML files, the corresponding commands are

   <CrackTip keypt='(keyid)'/>

   <CrackTip x='(x)' y='(y)'/>

where

  • (keyid) specifies a previously defined keypoint located at the crack tip.
  • Alternatively, the crack tip can be specified by the ((x),(y)) coordinates (in mm). The node closest to that point will be the crack tip node.

In theory, quarter point elements represent the singular stress field of a crack tip better than conventional elements. They were first proposed by Barsoum (1974).[2] Quarter point element are available in NairnFEA for research use, but they have drawbacks, and therefore are not recommended. Recent experience with crack-closure calculations shows that they converge more slowly to the correct energy release rate than conventional quadratic elements with midside nodes.[1] Furthermore, prior literature has errors in quarter-point element crack-closure equations, although NairnFEAMPM has the new corrected equations from Nairn (2011).[1] Finally, visualization of strain from FEA results is done by differentiation of the displacement field. Since quarter-point elements have singular strains at the crack tip, the strain in those elements will be undefined. These elements can still plot stresses.

References

  1. 1.0 1.1 1.2 J.A. Nairn, "Generalized Crack Closure Analysis for Elements with Arbitrarily-Placed Side Nodes and Consistent Nodal Forces," Int. J. Fracture, 171, 11-22 (2011). (PDF version)
  2. R. S. Barsoum, "Application of Quadratic Isoparametric Finite Elements in Linear Fracture Mechanics," Int. J. Fracture, 10, 603-605 (1974).