Virtual Crack Closure Technique Fundamentals

What Is It?

The Virtual Crack Closure Technique (VCCT) is a method for computing the energy release rate at the front of an existing crack in a finite element model. It is based on Irwin's crack closure concept: the energy required to extend a crack by a small amount equals the work done by the nodal forces at the crack front in closing the newly created crack surfaces. VCCT computes this work from the nodal forces at the crack front and the nodal displacements at the corresponding location behind the crack front. The technique resolves the energy release rate into Mode I, Mode II and Mode III components, enabling mixed-mode fracture assessment.

Why It Matters

VCCT is one of the two principal methods for predicting delamination propagation in composite structures (the other being cohesive-zone modelling). It is well-established, computationally efficient and widely implemented in commercial FEA codes. VCCT is particularly suited to problems where the crack path is known (the interface) and a pre-existing crack exists or can be introduced. It computes the energy release rate directly from the FEA solution without requiring special elements at the interface. Understanding VCCT — its principles, its assumptions, its strengths and its limitations — is essential for delamination propagation assessment.

VCCT is a principal method for delamination propagation prediction. It is efficient, well-established and widely available. It computes the energy release rate from nodal forces and displacements at an existing crack front. Best suited when the crack path is known and a pre-existing crack exists.

Energy Release Rate Concept

The energy release rate G is the energy available for crack growth per unit crack area. It is the derivative of the potential energy with respect to the crack area: G = −∂Π/∂A. If G exceeds the fracture toughness G_c, the crack propagates. VCCT computes G by evaluating the work that would be done by the nodal forces at the crack front if the crack were extended by a small amount and the new surfaces were closed. This work, divided by the new crack area, gives the energy release rate. The concept is based on Irwin's crack closure integral, adapted for the finite element discretisation.

Energy release rate concept:

G = energy available per unit crack area
  = −∂Π / ∂A    (potential energy / crack area)

Propagation criterion:
  G ≥ G_c  →  crack propagates

VCCT computes G from:
  G = work to close new crack surfaces / new crack area
  = (nodal forces × nodal displacements) / (2 × crack area)

The factor of 2 accounts for the
closure of two surfaces. The nodal
forces are at the crack front; the
displacements are at the corresponding
location behind the front.

Existing Crack Requirement

VCCT requires an existing crack — a pre-defined delamination with a crack front. Unlike cohesive-zone modelling, which captures initiation naturally, VCCT cannot predict where a crack will initiate; it can only assess whether an existing crack will propagate. The crack must be explicitly modelled in the FEA mesh — the nodes along the crack surfaces are unconnected (duplicated), allowing the surfaces to separate. The crack front is the boundary between the cracked and uncracked regions. VCCT evaluates G at each node along the crack front. If the analysis needs to predict initiation, a pre-crack must be introduced at the suspected initiation location, or a different method (cohesive zones or strength-based) must be used.

VCCT requires an existing, pre-defined crack — it cannot predict initiation. The crack must be explicitly modelled with duplicated nodes. VCCT evaluates G at the crack front. If initiation prediction is needed, introduce a pre-crack or use cohesive zones instead.

Nodal Forces and Displacements

VCCT computes the energy release rate from the nodal forces and displacements in the FEA solution. The nodal forces are the internal forces at the crack front nodes — the forces that the elements ahead of the crack front transmit across the would-be crack surface. The nodal displacements are the relative displacements between the crack surfaces behind the crack front — the opening and sliding at the corresponding location. The product of the force and the displacement gives the work to close the crack; divided by the crack area, it gives the energy release rate. The forces and displacements are extracted from a single analysis — no second analysis with an extended crack is needed (this is the "virtual" closure).

• Nodal forces — internal forces at the crack front nodes (forces transmitted across the crack surface)
• Nodal displacements — relative displacements between crack surfaces behind the front
• Product of force × displacement = work to close the crack
• Divided by crack area = energy release rate
• Extracted from a single analysis — no second analysis needed (virtual closure)

Mode I Component

The Mode I (opening) energy release rate G_I is computed from the normal nodal force and the normal opening displacement. The normal force F_n is the force in the direction perpendicular to the crack plane (the peel direction). The normal displacement δ_n is the opening displacement between the crack surfaces behind the front. G_I = (F_n × δ_n) / (2 × ΔA), where ΔA is the crack area associated with the node. G_I represents the energy available for opening-mode crack growth. If G_I exceeds the Mode I fracture toughness G_Ic, the crack propagates in Mode I.

Mode I energy release rate:

G_I = (F_n × δ_n) / (2 × ΔA)

where:
F_n = normal nodal force at crack front (peel direction)
δ_n = normal opening displacement behind crack front
ΔA = crack area associated with the node

G_I = energy for opening-mode crack growth
Propagation when:  G_I ≥ G_Ic

Note: The exact formula depends on the
VCCT implementation (element type,
crack front geometry). The concept is
the same — force × displacement / area.

Mode II Component

The Mode II (forward shear) energy release rate G_II is computed from the shear nodal force and the shear sliding displacement in the direction of crack propagation. The shear force F_s is the force in the propagation direction. The shear displacement δ_s is the relative sliding between the crack surfaces in the propagation direction. G_II = (F_s × δ_s) / (2 × ΔA). G_II represents the energy available for forward-shear crack growth. If G_II exceeds the Mode II fracture toughness G_IIc, the crack propagates in Mode II. In most composite interfaces, G_IIc is higher than G_Ic, so the interface is tougher in shear.

• G_II = (F_s × δ_s) / (2 × ΔA)
• F_s = shear nodal force in propagation direction
• δ_s = shear sliding displacement in propagation direction
• G_II = energy for forward-shear crack growth
• Propagation when G_II ≥ G_IIc
• G_IIc typically higher than G_Ic — tougher in shear

Mode III Component

The Mode III (anti-plane shear) energy release rate G_III is computed from the shear nodal force and the shear sliding displacement in the direction perpendicular to the propagation direction (the anti-plane direction). G_III = (F_t × δ_t) / (2 × ΔA), where F_t is the anti-plane shear force and δ_t is the anti-plane sliding displacement. G_III represents the energy available for anti-plane shear crack growth. Mode III is less common than Modes I and II but can be important in torsion or complex 3D loading. G_IIIc is the least commonly measured fracture toughness and is sometimes assumed equal to G_IIc when data is unavailable.

Mode III energy release rate:

G_III = (F_t × δ_t) / (2 × ΔA)

where:
F_t = anti-plane shear nodal force
δ_t = anti-plane shear sliding displacement
ΔA = crack area associated with the node

G_III = energy for anti-plane shear crack growth
Propagation when:  G_III ≥ G_IIIc

G_IIIc is the least commonly measured
toughness. Sometimes assumed equal
to G_IIc when data is unavailable —
this assumption should be checked.

Total Energy Release Rate

The total energy release rate is the sum of the three mode components: G = G_I + G_II + G_III. The total G represents the total energy available for crack growth. For mixed-mode propagation, the total G is compared against the mixed-mode fracture toughness, which depends on the mode mixity. The mode mixity at each point along the crack front is the ratio of the mode components (e.g. G_II/(G_I + G_II)). The mode mixity may vary along the crack front, so the mixed-mode toughness and the propagation criterion must be evaluated at each point.

Total G = G_I + G_II + G_III. For mixed-mode propagation, compare G against the mixed-mode toughness G_c(mode mixity). The mode mixity may vary along the crack front — evaluate at each point. See the mixed-mode delamination criteria article for interaction laws.

Assumptions and Limitations

VCCT is based on several assumptions that limit its applicability. It assumes self-similar crack growth — the crack extends in the same direction as the existing crack, along the same interface. It requires a pre-existing crack — it cannot predict initiation. It assumes that the crack closure work is accurately represented by the nodal forces and displacements, which requires a sufficiently fine mesh at the crack front. It assumes linear elastic behaviour at the crack front — the material ahead of the crack is undamaged. For problems involving large process zones, nonlinear material behaviour, or crack path uncertainty, cohesive-zone modelling may be more appropriate.

Key Takeaways

• VCCT computes the energy release rate from nodal forces and displacements at an existing crack front
• Based on Irwin's crack closure concept — work to close new crack surfaces divided by crack area
• Requires a pre-existing crack — cannot predict initiation
• Resolves G into Mode I (G_I), Mode II (G_II) and Mode III (G_III) components
• Total G = G_I + G_II + G_III; mixed-mode propagation uses mixed-mode toughness
• Assumes self-similar growth, linear elastic behaviour, fine mesh at crack front
• Best suited when crack path is known and a pre-existing crack exists; for initiation, use cohesive zones