Difference between revisions of "Explicit FEA Cracks"
| (20 intermediate revisions by the same user not shown) | |||
| Line 4: | Line 4: | ||
[[File:Feacrack.jpg|right]] | [[File:Feacrack.jpg|right]] | ||
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 | 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 | 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 | ||
| Line 10: | Line 10: | ||
== Crack Closure Calculations == | == 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 | 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).<ref name='CC'>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). ([http://www.cof.orst.edu/cof/wse/faculty/Nairn/papers/CrackClosure.pdf PDF version])</ref> 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 [[FEA Elements|element type]] allowed in [[NairnFEA]]. In addition, one special crack tip element available and it is | Crack closure calculations can be done for any [[FEA Elements|element type]] allowed in [[NairnFEA]]. In addition, one special crack tip element available and it is configured by using the <tt>CrackTip</tt> command. | ||
== Crack Tip Command == | == Crack Tip Command == | ||
[[File:qtrpt.jpg|right]] | |||
The <tt>CrackTip</tt> command coverts all elements at the crack tip to be quarter point elements. Quarter point elements are [[FEA Elements|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 [[FEA Elements|8-node quadrilaterals or 6-node triangles]] and then adding the scripted command: | |||
<tt> CrackTip (keyid)</tt> | |||
or | |||
<tt> CrackTip (x),(y)</tt> | |||
In <tt>XML</tt> files, the corresponding commands are | |||
<tt> <CrackTip keypt='(keyid)'/></tt> | |||
<tt> <CrackTip x='(x)' y='(y)'/></tt> | |||
where | |||
* <tt>(keyid)</tt> specifies a previously defined [[Keypoint Command|keypoint]] located at the crack tip. | |||
* Alternatively, the crack tip can be specified by the (<tt>(x),(y)</tt>) coordinates (in [[ConsistentUnits Command#Legacy and Consistent Units|length units]]). 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).<ref name='qtr'>R. S. Barsoum, "Application of Quadratic Isoparametric Finite Elements in Linear Fracture Mechanics," <i>Int. J. Fracture</i>, <b>10</b>, 603-605 (1974).</ref> 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.<ref name='CC'/> Furthermore, prior literature has errors in quarter-point element crack-closure equations, although [[NairnFEAMPM]] has the new corrected equations from Nairn (2011).<ref name='CC'/> 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 == | == References == | ||
<references/> | <references/> | ||
Latest revision as of 13:06, 2 June 2015
FEA calculations can model explicit cracks.
Explicit Cracks
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
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 length units). 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.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)
- ↑ R. S. Barsoum, "Application of Quadratic Isoparametric Finite Elements in Linear Fracture Mechanics," Int. J. Fracture, 10, 603-605 (1974).