Crack Propagation Analysis

Defining Initially Bonded Crack Surfaces in Abaqus/Standard

Potential crack surfaces are modeled as secondary and main contact surfaces (see About Contact Pairs in Abaqus/Standard). Any contact formulation except the finite-sliding, surface-to-surface formulation can be used. The predetermined crack surfaces are assumed to be initially partially bonded so that the crack tips can be identified explicitly by Abaqus/Standard. Initially bonded crack surfaces cannot be used with self-contact.

Define an initial condition to identify which part of the crack is initially bonded. You specify the secondary surface, the main surface, and a node set that identifies the initially bonded part of the secondary surface. The unbonded portion of the secondary surface will behave as a regular contact surface. Either the secondary surface or the main surface must be specified; if only the main surface is given, all of the secondary surfaces associated with this main surface that have nodes in the node set will be bonded at these nodes.

If a node set is not specified, the initial contact conditions will apply to the entire contact pair; in this case, no crack tips can be identified, and the bonded surfaces cannot separate.

If a node set is specified, the initial conditions apply only to the secondary nodes in the node set. Abaqus/Standard checks to ensure that the node set defined includes only secondary nodes belonging to the contact pair specified.

By default, the nodes in the node set are considered to be initially bonded in all directions.

Input File Usage

INITIAL CONDITIONS, TYPE=CONTACT

Bonding Only in the Normal Direction

For fracture criteria based on the critical stress, critical crack opening displacement, or crack length versus time, it is possible to bond the nodes in the node set (or the contact pair if a node set is not defined) only in the normal direction. In this case the nodes are allowed to move freely tangential to the contact surfaces. Friction (Frictional Behavior) cannot be specified if the nodes are bonded only in the normal direction.

Bonding only in the normal direction is typically used to model bonded contact conditions in Mode I crack problems where the shear stress ahead of the crack along the crack plane is zero.

Input File Usage

INITIAL CONDITIONS, TYPE=CONTACT, NORMAL

Activating the Crack Propagation Capability in Abaqus/Standard

The crack propagation capability must be activated within the step definition to specify that crack propagation may occur between the two surfaces that are initially partially bonded. You specify the surfaces along which the crack propagates.

If the crack propagation capability is not activated for partially bonded surfaces, the surfaces will not separate; in this case the specified initial contact conditions would have the same effect as that provided by the tied contact capability, which generates a permanent bond between two surfaces during the entire analysis (see Defining Tied Contact in Abaqus/Standard).

Input File Usage

DEBOND, SECONDARY=secondary_surface_name, MAIN=main_surface_name

Propagation of Multiple Cracks

Cracks can propagate from either a single crack tip or multiple crack tips. The crack propagation capability in Abaqus/Standard requires that the surfaces be initially partially bonded so that the crack tips can be identified. A contact pair can have crack propagation from multiple crack tips. However, only one crack propagation criterion is allowed for a given contact pair. Crack propagation along several contact pairs can be modeled by specifying multiple crack propagation definitions.

Defining and Activating Crack Propagation in Abaqus/Explicit

In Abaqus/Explicit potential crack surfaces are modeled as bonded general contact surfaces (see About General Contact in Abaqus/Explicit) in the context of surface-based cohesive behavior (see Contact Cohesive Behavior). The capability is implemented using a pure main-secondary formulation. As is the case in Abaqus/Standard, the predetermined crack surfaces are assumed to be initially partially bonded so that the crack tips can be identified explicitly.

To identify which pair of surfaces determine the crack and which part of the crack is initially bonded, you must define and assign a contact clearance (see Contact Initialization for General Contact in Abaqus/Explicit). You first define a contact clearance to specify the node set that is initially bonded, and then you assign this contact clearance to a pair of two single-sided surfaces that define the crack. The unbonded portion behaves as a regular contact surface. The nodes in the node set are considered to be initially bonded in all directions.

The crack tip is identified only from the specified two surfaces and the node set. No attempt is made to determine a crack tip from all surfaces included in the general contact domain. Consequently, to be able to identify the crack tip, the surface including the specified node set must extend past the node set. Otherwise, the surfaces will not debond, and the crack cannot propagate.

You complete the definition of the crack propagation capability by defining a fracture-based cohesive behavior surface interaction. You activate the crack propagation by assigning it to the pair of surfaces that are initially partially bonded. If the fracture criterion is met, crack propagation occurs between these two surfaces. Cohesive behavior is also used to specify the elastic behavior of the bonds (see Contact Cohesive Behavior).

If a fracture-based surface interaction is not assigned to a pair of surfaces, the crack definition is incomplete. Unlike Abaqus/Standard where the identified nodes will stay bonded if the crack is not activated, in Abaqus/Explicit the nodes identified by the contact clearance definition will separate without generating any interface stress.

Similar to Abaqus/Standard, cracks can propagate from single or multiple crack tips for the same pair of surfaces.

Input File Usage

Use the following options:

CONTACT CLEARANCE, NAME=clearance_name, SEARCH NSET=bonded_nset_name
**
SURFACE INTERACTION, NAME=interaction_name
COHESIVE BEHAVIOR
FRACTURE CRITERION
..**
CONTACT
CONTACT CLEARANCE ASSIGNMENT
secondary_surface, main_surface, clearance_name
CONTACT PROPERTY ASSIGNMENT
secondary_surface, main_surface, interaction_name

Specifying a Fracture Criterion

You can specify the crack propagation criteria, as discussed below. Table 1 shows which criteria are supported by Abaqus/Standard and Abaqus/Explicit. Only one crack propagation criterion is allowed per contact pair even if multiple cracks are present.

Table 1. Fracture criterion.
Crack propagation criterion	Abaqus/Standard	Abaqus/Explicit
Critical stress	Yes	No
Critical crack opening displacement	Yes	No
Crack length versus time	Yes	No
VCCT	Yes	Yes
Enhanced VCCT	Yes	No
Fatigue crack growth	Yes	No

Crack propagation analysis is carried out on a nodal basis. The crack-tip node debonds when the fracture criterion, f, reaches the value 1.0 within a given tolerance:

f_{L L} \leq f \leq f_{U L},

where $f_{U L} = 1 + f_{t o l}$ and $f_{L L} = 1$ for VCCT, enhanced VCCT, and fatigue crack growth criteria or $f_{L L} = 1 - f_{t o l}$ for other fracture criteria. You can specify the tolerance $f_{t o l}$ . In Abaqus/Standard, if $f > 1 + f_{t o l}$ , the time increment is cut back such that the crack propagation criterion is satisfied except in the case of an unstable crack growth problem where multiple nodes at and ahead of a crack tip are allowed to debond without the cut back of increment size in one increment. The default value of $f_{t o l}$ is 0.1 for the critical stress, critical crack opening displacement, and crack length versus time criteria and is 0.2 for the VCCT and enhanced VCCT criteria.

Input File Usage

FRACTURE CRITERION, TOLERANCE= $f_{t o l}$ , TYPE=type

Critical Stress Criterion

This criterion is available only in Abaqus/Standard.

If you specify a critical stress criterion at a critical distance ahead of the crack tip, the crack-tip node debonds when the local stress across the interface at a specified distance ahead of the crack tip reaches a critical value.

This criterion is typically used for crack propagation in brittle materials. The critical stress criterion is defined as

f = \sqrt{{(\frac{{\hat{σ}}_{n}}{σ^{f}})}^{2} + {(\frac{τ_{1}}{τ_{1}^{f}})}^{2} + {(\frac{τ_{2}}{τ_{2}^{f}})}^{2}}, {\hat{σ}}_{n} = max (σ_{n}, 0),

where $σ_{n}$ is the normal component of stress carried across the interface at the distance specified; $τ_{1}$ and $τ_{2}$ are the shear stress components in the interface; and $σ^{f}$ and $τ_{1}^{f}$ are the normal and shear failure stresses, which you must specify. The second component of the shear failure stress, $τ_{2}$ , is not relevant in a two-dimensional analysis; therefore, the value of $τ_{2}^{f}$ need not be specified. The crack-tip node debonds when the fracture criterion, f, reaches the value 1.0.

If the value of $τ_{1}^{f}$ is not given or is specified as zero, it will be taken to be a very large number so that the shear stress has no effect on the fracture criterion.

The distance ahead of the crack tip is measured along the secondary surface, as shown in Figure 1. The stresses at the specified distance ahead of the crack tip are obtained by interpolating the values at the adjacent nodes. The interpolation depends on whether first-order or second-order elements are used to define the secondary surface.

Input File Usage

FRACTURE CRITERION, TYPE=CRITICAL STRESS, DISTANCE=n

Critical Crack Opening Displacement Criterion

This criterion is available only in Abaqus/Standard.

If you base the crack propagation analysis on the crack opening displacement criterion, the crack-tip node debonds when the crack opening displacement at a specified distance behind the crack tip reaches a critical value. This criterion is typically used for crack propagation in ductile materials.

The crack opening displacement criterion is defined as

f = \frac{δ}{δ_{c}},

where $δ$ is the measured value of crack opening displacement and $δ_{c}$ is the critical value of the crack opening displacement (user-specified). The crack-tip node debonds when the fracture criterion reaches the value 1.0.

You must supply the crack opening displacement versus cumulative crack length data. In Abaqus/Standard the cumulative crack length is defined as the distance between the initial crack tip and the current crack tip measured along the secondary surface in the current configuration. The crack opening displacement is defined as the normal distance separating the two faces of the crack at the given distance.

You specify the position, n, behind the crack tip where the critical crack opening displacement is calculated. The value of this position must be specified as the length of the straight line joining the current crack tip and points on the secondary and main surfaces (Figure 2).

Abaqus/Standard computes the crack opening displacement at that point by interpolating the values at the adjacent nodes. The interpolation depends on whether first-order or second-order elements are used to define the secondary surface. An error message will be issued if the value of n is not within the end points of the contact pair.

Input File Usage

FRACTURE CRITERION, TYPE=COD, DISTANCE=n

Modeling Symmetry

In problems where the debonding surfaces lie on a symmetry plane, you can specify that Abaqus/Standard should consider only half of the user-specified crack opening displacement values. In this case the initial bonding must be in the normal direction only (see Bonding Only in the Normal Direction above).

Input File Usage

FRACTURE CRITERION, TYPE=COD, DISTANCE=n, SYMMETRY

Crack Length Versus Time Criterion

This criterion is available only in Abaqus/Standard.

To specify the crack propagation explicitly as a function of total time, you must provide a crack length versus time relationship and a reference point from which the crack length is measured. This reference point is defined by specifying a node set. Abaqus/Standard finds the average of the current positions of the nodes in the set to define the reference point. During crack propagation the crack length is measured from this user-specified reference point along the secondary surface in the deformed configuration. The time specified must be total time, not step time.

The fracture criterion, f, is stated in terms of the user-specified crack length and the length of the current crack tip. The length of the current crack tip from the reference point is measured as the sum of the straight line distance of the initial crack tip from the reference point and the distance between the initial crack tip and the current crack tip measured along the secondary surface.

Referring to Figure 3, let node 1 be the initial location of the crack tip and node 3 be the current location of the crack tip. The distance of the current crack tip located at node 3 is given by

l_{3} = l_{1} + Δ l_{12} + Δ l_{23},

where $l_{1}$ is the length of the straight line joining node 1 and the reference point, $Δ l_{12}$ is the distance between nodes 1 and 2, and $Δ l_{23}$ is the distance between nodes 2 and 3 measured along the secondary surface.

The fracture criterion, f, is given by

f = \frac{l - (l_{3} - Δ l_{23})}{Δ l_{23}},

where l is the length at the current time obtained from the user-specified crack length versus time curve. Crack-tip node 3 will debond when the failure function f reaches the value of 1.0 (within the user-defined tolerance).

If geometric nonlinearity is considered in the step (Defining an Analysis), the reference point may move as the body deforms; you must ensure that this movement does not invalidate the crack length versus time criterion.

Abaqus/Standard does not extrapolate beyond the end points of your crack data. Therefore, if the first crack length specified is greater than the distance from the crack reference point to the first bonded node, the first bonded node will never debond and the crack will not propagate. In this case Abaqus/Standard will print warning messages in the message (.msg) file.

Input File Usage

FRACTURE CRITERION, TYPE=CRACK LENGTH, NSET=name

VCCT Criterion

This criterion is available in both Abaqus/Standard and Abaqus/Explicit.

The Virtual Crack Closure Technique (VCCT) criterion uses the principles of linear elastic fracture mechanics (LEFM) and, therefore, is appropriate for problems in which brittle crack propagation occurs along predefined surfaces.

VCCT is based on the assumption that the strain energy released when a crack is extended by a certain amount is the same as the energy required to close the crack by the same amount. For example, Figure 4 illustrates the similarity between crack extension from i to j and crack closure at j.

In Figure 5 nodes 2 and 5 will start to release when

f = \frac{G_{I}}{G_{I C}} = \frac{1}{2} (\frac{v_{1, 6} F_{v, 2, 5}}{b d}) \frac{1}{G_{I C}} \geq 1.0,

where $G_{I}$ is the Mode I energy release rate, $G_{I C}$ is the critical Mode I energy release rate, b is the width, d is the length of the elements at the crack front, $F_{v, 2, 5}$ is the vertical force between nodes 2 and 5, and $v_{1, 6}$ is the vertical displacement between nodes 1 and 6. Assuming that the crack closure is governed by linear elastic behavior, the energy to close the crack (and, thus, the energy to open the crack) is calculated from the previous equation. Similar arguments and equations can be written in two dimensions for Mode II and for three-dimensional crack surfaces including Mode III.

In the general case involving Mode I, II, and III the fracture criterion is defined as

f = \frac{G_{e q u i v}}{G_{e q u i v C}} \geq 1.0,

where $G_{e q u i v}$ is the equivalent strain energy release rate calculated at a node, and $G_{e q u i v C}$ is the critical equivalent strain energy release rate calculated based on the user-specified mode-mix criterion and the bond strength of the interface. The crack-tip node will debond when the fracture criterion reaches the value of 1.0.

Abaqus provides three common mode-mix formulae for computing $G_{e q u i v C}$ : the BK law, the power law, and the Reeder law models. The choice of model is not always clear in any given analysis; an appropriate model is best selected empirically.

BK Law

The BK law model is described in Benzeggagh (1996) by the following formula:

G_{e q u i v C} = G_{I C} + (G_{I I C} - G_{I C}) {(\frac{G_{I I} + G_{I I I}}{G_{I} + G_{I I} + G_{I I I}})}^{η},

G_{e q u i v} = G_{I} + G_{I I} + G_{I I I} .

To define this model, you must provide $G_{I C}, G_{I I C},$ and $η$ . This model provides a power law relationship combining energy release rates in Mode I, Mode II, and Mode III into a single scalar fracture criterion.

Input File Usage

FRACTURE CRITERION, TYPE=VCCT, MIXED MODE BEHAVIOR=BK

Power Law

The power law model is described in Wu (1965) by the following formula:

\frac{G_{e q u i v}}{G_{e q u i v C}} = {(\frac{G_{I}}{G_{I C}})}^{a_{m}} + {(\frac{G_{I I}}{G_{I I C}})}^{a_{n}} + {(\frac{G_{I I I}}{G_{I I I C}})}^{a_{o}} .

To define this model, you must provide $G_{I C}, G_{I I C}, G_{I I I C}, a_{m}, a_{n},$ and $a_{o}$ .

Input File Usage

FRACTURE CRITERION, TYPE=VCCT, MIXED MODE BEHAVIOR=POWER

Reeder Law

The Reeder law model is described in Reeder (2002) by the following formula:

G_{e q u i v C} = G_{I C} + (G_{I I C} - G_{I C}) {(\frac{G_{I I} + G_{I I I}}{G_{I} + G_{I I} + G_{I I I}})}^{η} +

(G_{I I I C} - G_{I I C}) (\frac{G_{I I I}}{G_{I I} + G_{I I I}}) {(\frac{G_{I I} + G_{I I I}}{G_{I} + G_{I I} + G_{I I I}})}^{η},

G_{e q u i v} = G_{I} + G_{I I} + G_{I I I} .

To define this model, you must provide $G_{I C}, G_{I I C}, G_{I I I C},$ and $η$ . The Reeder law is best applied when $G_{I I C} \neq G_{I I I C}$ . When $G_{I I C} = G_{I I I C}$ , the Reeder law reduces to the BK law. The Reeder law applies only to three-dimensional problems.

Input File Usage

FRACTURE CRITERION, TYPE=VCCT, MIXED MODE BEHAVIOR=REEDER

Defining Variable Critical Energy Release Rates

You can define a VCCT criterion with varying energy release rates by specifying the critical energy release rates at the nodes.

If you indicate that the nodal critical energy rates will be specified, any constant critical energy release rates you specify are ignored, and the critical energy release rates are interpolated from the nodes. The critical energy release rates must be defined at all nodes on the secondary surface.

Input File Usage

Use both of the following options:

FRACTURE CRITERION, TYPE=VCCT, NODAL ENERGY RATE
NODAL ENERGY RATE

Enhanced VCCT Criterion

This criterion is available only in Abaqus/Standard.

The enhanced VCCT criterion is very similar to the original VCCT criterion described above. As in the original VCCT criterion, the fracture criterion in the general case involving Mode I, II, and III is defined as

f = \frac{G_{e q u i v}}{G_{e q u i v C}} \geq 1.0.

The crack-tip node debonds when the fracture criterion reaches the value of 1.0. However, unlike the original VCCT criterion, you can specify two different critical fracture energy release rates: $G_{C}$ for the onset of a crack and $G_{C}^{P}$ for the growth of a crack. When the enhanced VCCT criterion is used in the general case involving Mode I, II, and III fracture, the amount of energy dissipated associated with the release of the debonding force is controlled by the critical equivalent strain energy release rate required to propagate the crack, $G_{e q u i v C}^{P}$ , rather than by the critical equivalent strain energy release rate required to initiate the crack, $G_{e q u i v C} .$ The formulae for calculating $G_{e q u i v C}^{P}$ are identical to those used for $G_{e q u i v C}$ for different mixed-mode fracture criteria.

Input File Usage

FRACTURE CRITERION, TYPE=ENHANCED VCCT

Fatigue Crack Growth Criterion

This criterion is available only in Abaqus/Standard.

If you specify the fatigue crack growth criterion, progressive delamination growth at the interfaces in laminated composites subjected to sub-critical cyclic loadings can be simulated. This criterion can be used only in the general fatigue crack growth approach (Linear Elastic Fatigue Crack Growth Analysis) or a low-cycle fatigue analysis using the direct cyclic approach (Low-Cycle Fatigue Analysis Using the Direct Cyclic Approach). The onset and delamination growth are characterized by using the Paris law, which relates the relative fracture energy release rate to crack growth rates as illustrated in Figure 6. An alternative form to the original Paris law, which better accounts for the mixed mode fatigue crack growth, is also available as shown in Figure 7. The fracture energy release rates at the crack tips in the interface elements are calculated based on the above mentioned VCCT technique.

The Paris regime is bounded by the energy release rate threshold, $G_{t h r e s h}$ , below which there is no consideration of fatigue crack initiation or growth, and the energy release rate upper limit, $G_{p l}$ , above which the fatigue crack will grow at an accelerated rate. $G_{C}$ is the critical equivalent strain energy release rate calculated based on the user-specified mode-mix criterion and the bond strength of the interface. The formulae for calculating $G_{C}$ have been provided above for different mixed mode fracture criteria. You can specify the ratio of $G_{t h r e s h}$ over $G_{C}$ and the ratio of $G_{p l}$ over $G_{C}$ . The default values are $\frac{G_{t h r e s h}}{G_{C}} = 0.01$ and $\frac{G_{p l}}{G_{C}} = 0.85$ .

Input File Usage

FRACTURE CRITERION, TYPE=FATIGUE

Onset of Delamination Growth

The onset of delamination growth refers to the beginning of fatigue crack growth at the crack tip along the interface. In a fatigue crack growth analysis the onset of the fatigue crack growth criterion is characterized by $Δ G$ , which is the relative fracture energy release rate when the structure is loaded between its maximum and minimum values. The fatigue crack growth initiation criterion is defined as

f = \frac{N}{c_{1} Δ G^{c_{2}}} \geq 1.0,

where $c_{1}$ and $c_{2}$ are material constants and $N$ is the cycle number. The interface elements at the crack tips will not be released unless the above equation is satisfied and the maximum fracture energy release rate, $G_{m a x}$ , which corresponds to the cyclic energy release rate when the structure is loaded up to its maximum value, is greater than $G_{t h r e s h}$ . If you do not specify the onset criterion, Abaqus/Standard assumes that the onset of fatigue crack growth is satisfied automatically.

Fatigue Delamination Growth Using the Paris Law

The rate of the delamination growth per cycle is given by the Paris law if $G_{t h r e s h} < G_{m a x} < G_{p l},$

\frac{d a}{d N} = c_{3} Δ G^{c_{4}},

where $c_{3}$ and $c_{4}$ are material constants.

Ratcliffe and Johnston (2014) and Deobald et al. (2017) proposed the following alternate form of the fatigue law, which better accounts for mixed-mode fatigue crack growth:

\frac{d a}{d N} = c_{3} G_{T M a x}^{c_{4}} .

In the above expression, $G_{T M a x}$ is the total maximum strain energy release rate (as opposed to the strain energy release rate change over a cycle used in the original form of the Paris law), while $c_{3}$ and $c_{4},$ are material parameters that depend on mode-mix and stress ratios. Abaqus does not support the above form of the crack growth rate equation directly, but instead allows specification of $d a / d N$ as a tabular function of $G_{T M a x}$ , the mode-mix ratio, and the stress ratio.

In addition, user subroutine UMIXMODEFATIGUE provides a general capability for implementing a user-defined fatigue crack growth law.

At the end of cycle $N$ , Abaqus/Standard extends the crack length, $a_{N}$ , from the current cycle forward over an incremental number of cycles, $Δ N$ to $a_{N + Δ N}$ by releasing at least one element at the interface. Given the Paris law parameters, $c_{3}$ and $c_{4},$ combined with the known node spacing $Δ a_{N}_{j} = a_{N + Δ N} - a_{N}$ at the interface elements at the crack tips, the number of cycles necessary to fail each interface element at the crack tip can be calculated as $Δ N_{j}$ , where j represents the node at the jthe crack tip. The analysis is set up to release at least one interface element after the loading cycle is completed. The element with the fewest cycles is identified to be released, and its $Δ N_{m i n} = m i n (Δ N_{j})$ is represented as the number of cycles to grow the crack equal to its element length, $Δ a_{N}_{m i n} = m i n (Δ a_{N}_{j})$ . The most critical element is completely released with a zero constraint and a zero stiffness at the end of the completed cycle. As the interface element is released, the load is redistributed and a new relative fracture energy release rate must be calculated for the interface elements at the crack tips for the next cycle. This capability allows at least one interface element at the crack tips to be released after each completed cycle and precisely accounts for the number of cycles needed to cause fatigue crack growth over that length.

If $G_{m a x} > G_{p l}$ , the interface elements at the crack tips will be released by increasing the cycle number count, $d N$ , by one only.

For information on how to accelerate the fatigue crack growth analysis and to provide a smooth solution for the crack front, see Controlling Element Fracture.

Input File Usage

Use one of the following options to define the fatigue delamination onset and growth criterion based on the original form of the Paris law:

FRACTURE CRITERION, TYPE=FATIGUE,
MIXED MODE BEHAVIOR=BK (default)

FRACTURE CRITERION, TYPE=FATIGUE,
MIXED MODE BEHAVIOR=POWER

FRACTURE CRITERION, TYPE=FATIGUE,
MIXED MODE BEHAVIOR=REEDER

Use the following option to define the fatigue delamination onset and growth criterion based on the alternative form of the fatigue law:

FRACTURE CRITERION, TYPE=FATIGUE,
MIXED MODE BEHAVIOR=TABULAR

Use the following option to define a fatigue delamination onset and growth criterion using user subroutine UMIXMODEFATIGUE:

FRACTURE CRITERION, TYPE=FATIGUE,
MIXED MODE BEHAVIOR=USER

Specifying How a Debonding Force Is Released after a Fracture Criterion Is Met in Abaqus/Standard

After debonding, the traction between two surfaces is initially carried as equal and opposite forces at the secondary node and the corresponding point on the main surface. The debonding force is released as the crack opens and advances. Once complete debonding has occurred at a point, the bond surfaces act like standard contact surfaces with associated interface characteristics. There are several ways to release the debonding force, depending on the fracture criterion that you specify.

Specifying a Debonding Amplitude Curve

When you use the critical stress, critical crack opening displacement, or crack length versus time fracture criteria, you can define how this force is to be reduced to zero with time after debonding starts at a particular node on the bonded surface. You specify a relative amplitude, a, as a function of time after debonding starts at a node. Thus, suppose the force transmitted between the surfaces at secondary node N is ${T^{N} |}_{0}$ when that node starts to debond, which occurs at time ${t^{N} |}_{0}$ . Then, for any time $t > {t^{N} |}_{0}$ the force transmitted between the surfaces at node N is ${a (t - {t^{N} |}_{0}) T^{N} |}_{0}$ . The relative amplitude must be 1.0 at the relative time 0.0 and must reduce to 0.0 at the last relative time point given.

The best choice of the amplitude curve depends on the material properties, specified loading, and the crack propagation criterion. If the stresses are removed too rapidly, the resulting large changes in the strains near the crack tip can cause convergence difficulties. For large-strain problems severe mesh distortion can also occur. For problems with rate-independent materials a linear amplitude curve is normally adequate. For problems with rate-dependent materials the stresses should be ramped off more slowly at the beginning of debonding to avoid convergence and mesh distortion difficulties. To reduce the likelihood of convergence and mesh distortion difficulties, you can reduce the value of the debond stress by 25% in 50% of the time to debond. The solution should not be strongly influenced by the details of the unloading procedure; if it is, this usually indicates that the mesh should be refined in the debond region.

Input File Usage

DEBOND, SECONDARY=secondary, MAIN=main
Data lines to define debonding amplitude curve

Ramping Down Debonding Force for the VCCT and the Enhanced VCCT Criteria

For the VCCT and the enhanced VCCT criteria, when the energy release rate exceeds the critical value at a crack tip, you can either release the traction between the two surfaces at the crack tip immediately during the following increment or release the traction gradually during succeeding increments with the reduction of the magnitude of the debonding force being governed by the critical fracture energy release rate. The latter approach is sometimes recommended to avoid sudden loss of stability when the crack tip is advanced. The enhanced VCCT criterion is meaningful only when used in conjunction with the latter approach. When the former approach is used, the results obtained by using the enhanced VCCT criterion are identical to those obtained by using the original VCCT criterion.

Input File Usage

Use the following option to release the traction immediately:

DEBOND, SECONDARY=secondary, MAIN=main, 
DEBONDING FORCE=STEP

Use the following option to release the traction gradually:

DEBOND, SECONDARY=secondary, MAIN=main, 
DEBONDING FORCE=RAMP

Ramping Down Debonding Force for the Fatigue Crack Growth Criteria

At the end of the loading cycle, the interface element that takes the fewest cycles, $Δ N_{m i n}$ , to fracture is released completely with zero constraint and a zero debonding force immediately after fracture. All of the other interface elements satisfied with the accuracy of the damage extrapolation condition (as defined in Controlling Element Fracture) are fractured. The debonding force is removed immediately after fracture or ramped down gradually according to a scalar damage variable defined as:

D_{N + Δ N} = D_{N} + \frac{\frac{d a}{d N} Δ N_{\min}}{a_{N}},

where $\frac{d a}{d N}$ is the crack growth rate, and $a_{N}$ is the element length. The debonding force becomes zero when the scalar damage variable reaches the value of 1.0 within a given tolerance, f:

D_{N + Δ N} \geq 1 - f .

The default value of f is 0.05.

Input File Usage

Use the following options to release the traction immediately:

DEBOND, DEBONDING FORCE=STEP (default) 
FRACTURE CRITERION, TYPE=FATIGUE

Use the following options to release the traction gradually:

DEBOND, DEBONDING FORCE=RAMP
FRACTURE CRITERION, TYPE=FATIGUE, TOLERANCE= $f$

Releasing Multiple Nodes in One Increment in Abaqus/Standard

For an unstable crack growth problem, sometimes it is more efficient to allow multiple nodes at and ahead of a crack tip to debond in one increment without cutting back the increment size when the VCCT or the enhanced VCCT fracture criterion is satisfied. This capability is activated automatically if you specify an unstable growth tolerance, $f_{t o l}^{u}$ . In this case if the fracture criterion, f, is within the given unstable growth tolerance:

1 + f_{t o l} \leq f \leq 1 + f_{t o l}^{u},

where $f_{t o l}$ is the tolerance described earlier, the time increment size by default is immediately reduced automatically to a very small value, $β Δ t_{\min}$ . This allows more nodes to debond until $f < 1$ for all the nodes ahead of the crack tip.

However, you can optionally specify the maximum number of cutbacks allowed, $n$ , to be controlled by the regular tolerance, $f_{t o l}$ , prior to the activation of the unstable growth tolerance in an increment. After this limit is reached, the time increment size is automatically recovered to a larger value, $α Δ t_{p r e}$ , where $Δ t_{\min}$ is the minimum time increment allowed; $Δ t_{p r e}$ is the time increment size prior to the unstable crack growth; and $α$ , $β$ , and $n$ are scaling parameters. The default values of $α$ , $β$ , and $n$ are 0.5, 2.0, and 0, respectively.

If you do not specify a value for the unstable growth tolerance, the default value is infinity. In this case the fracture criterion, f, for unstable crack growth is not limited by any upper-bound value in the above equation. The forces at those debonded nodes can be completely released immediately during the following increment or ramped down gradually during succeeding increments with the reduction of the magnitude of the debonding force being governed by the critical fracture energy release rate. If the forces at the debonded nodes are ramped down gradually, an unsymmetric matrix storage and solution scheme (Matrix Storage and Solution Scheme in Abaqus/Standard) is used by default to improve convergence.

Input File Usage

Use one of the following sets of options:

DEBOND, DEBONDING FORCE=RAMP
FRACTURE CRITERION, TYPE=VCCT, UNSTABLE GROWTH TOLERANCE= $f_{t o l}^{u}$ 
CONTROLS, TYPE=NO CUTBACK SCALING

or

DEBOND, DEBONDING FORCE=STEP
FRACTURE CRITERION, TYPE=VCCT, UNSTABLE GROWTH TOLERANCE= $f_{t o l}^{u}$

Procedures

Crack propagation analysis can be performed for static or dynamic overloadings using the following procedures:

It can also be performed for sub-critical cyclic fatigue loadings using the following procedures:

Controlling Time Incrementation during Debonding in Abaqus/Standard

When automatic incrementation is used for any criteria other than VCCT, enhanced VCCT, or fatigue crack growth, you can specify the size of the time increment used just after debonding starts. By default, the time increment is equal to the last relative time specified. However, if a fracture criterion is met at the beginning of an increment, the size of the time increment used just after debonding starts will be set equal to the minimum time increment allowed in this step.

For fixed time incrementation the specified time increment value will be used as the time increment size after debonding starts if Abaqus/Standard finds it needs a smaller time increment than the fixed time increment size. The time increment size will be modified as required until debonding is complete.

Input File Usage

DEBOND, SECONDARY=secondary, MAIN=main, TIME INCREMENT=t

Viscous Regularization for VCCT in Abaqus/Standard

The simulation of structures with unstable propagating cracks is challenging and difficult. Nonconvergent behavior may occur from time to time. While the usual stabilization techniques (such as contact pair stabilization and static stabilization) can be used to overcome some convergence difficulties, localized damping is included for VCCT or enhanced VCCT by using the viscous regularization technique. Viscous regularization damping causes the tangent stiffness matrix of the softening material to be positive for sufficiently small time increments.

Input File Usage

Use one of following options:

FRACTURE CRITERION, TYPE=VCCT, VISCOSITY= $μ$

FRACTURE CRITERION, TYPE=ENHANCED VCCT, VISCOSITY= $μ$

Linear Scaling to Accelerate Convergence for VCCT in Abaqus/Standard

For most crack propagation simulations using VCCT or the enhanced VCCT criterion, the deformation can be nearly linear up to the point of the onset of crack growth; past this point the analysis becomes very nonlinear. In this case a linear scaling method can be used to effectively reduce the solution time to reach the onset of crack growth.

Suppose that an applied “trial” load at increment $t = t_{i}$ is just a fraction of the critical load at the onset time of crack growth, $t = t_{c r i t}$ . The following algorithm is used in Abaqus/Standard to quickly converge to the critical load state:

Δ t_{i + 1} = (β_{i} \sqrt{\frac{G_{e q u i v C}}{G_{e q u i v}}} - 1) t_{i},

where initially $β_{i}$ would be set between 0.7 and 0.9 depending on the degree of nonlinearity (the default value is 0.9). When $Δ t_{i + 1}$ becomes smaller than 0.5% (indicating that the load is within 0.5% of its critical value), the next $β_{i}$ is automatically set to 1.0 to cause the most critical crack-tip node to precisely reach the critical value at the next increment. After the first crack-tip node releases, the linear scaling calculations are no longer valid and the time increment is set to the default value. Cutback is then allowed.

Input File Usage

CONTROLS, TYPE=VCCT LINEAR SCALING

Nonlocal Smoothing of the Crack Tangential Directions Using Least-Squares Approximation

Crack propagation analysis is carried out on a nodal basis. The crack-tip node debonds when the fracture criterion is satisfied. The crack front is formed by connecting the individual crack-tip nodes. For coarse or unstructured meshes, the crack front might be jagged (not smooth). You can use a moving least-squares approximation by polynomials to smooth out the crack tangential directions of the crack segments along the crack front. This approximation provides a smoother pseudocrack tangential direction, leading to a more accurate evaluation of the strain energy release rates under mixed-mode loading.

Input File Usage

FRACTURE CRITERION, POSITION=NONLOCAL

Specifying the Approximation Used in the Least-Squares Approximation

You can use a linear, quadratic, or cubic polynomial approximation for the moving least-squares approximation to smoothen the crack tangential directions. You specify the number of terms in the polynomial.

Input File Usage

Use one of the following options:

FRACTURE CRITERION, NPOLY=4
FRACTURE CRITERION, NPOLY=7 (default)
FRACTURE CRITERION, NPOLY=10

Specifying the Region of the Model Used for Nonlocal Smoothing of the Crack Tangential Directions

To control the range of elements used for nonlocal smoothing of the crack tangential directions along the crack front, you can specify a radius, $r_{c}$ , within which the elements around the crack tip along the crack front are included. The default radius is three times the typical element characteristic length along the crack front.

Input File Usage

FRACTURE CRITERION, R CRACK DIRECTION= $r_{c}$

Tips for Using the VCCT or Enhanced VCCT Criterion in Abaqus/Standard

Crack propagation problems using the VCCT or enhanced VCCT criterion are numerically challenging. The following tips will help you create a successful Abaqus/Standard model:

An analysis with the VCCT or enhanced VCCT criterion requires small time increments. Abaqus/Standard tracks the location of the active crack front node by node when the VCCT or enhanced VCCT criterion is used. Therefore, the crack front is allowed to advance only a single node forward in any single increment (although such an advance may take place across the entire crack front in three-dimensional problems). Because an analysis using the VCCT or enhanced VCCT criterion provides detailed results of the growth of the crack, you will need small time increments, especially if the mesh is highly refined.
Three different types of damping can be used to aid convergence for a model using the VCCT or enhanced VCCT criterion: contact stabilization, automatic or static stabilization, and viscous regularization. Contact and automatic stabilization are not specific to VCCT; they are built into Abaqus/Standard and are compatible with VCCT. Setting the value of the damping parameters is often an iterative procedure. If your VCCT model fails to converge due to unstable crack propagation, set the damping parameters to relatively high values and rerun the analysis. If the parameters are high enough, stable incrementation should return. However, the crack propagation behavior may have been modified by the damping forces and may not be physically correct. To monitor the energy absorbed by viscous damping, plot the damping energy and compare the results to the total strain energy in the model (ALLSE). When set properly, the value of the damping energy should be a small fraction of the total energy. Monitor the damping energy to ensure that the results of the VCCT simulation are reasonable in the presence of damping. When you use contact or automatic stabilization, Abaqus writes the damping energy to the variable ALLSD in the output database (.odb) file. When you use viscous regularization, Abaqus writes the damping energy to the variable ALLVD.
To maximize the accuracy of the debonding simulation, try to use matched meshes between the secondary and main surfaces of the debonding contact pair.
If you do use a mismatched mesh, you can maximize the accuracy of the simulation by using the small-sliding, surface-to-surface formulation for the contact pair (see Contact Formulations in Abaqus/Standard).
Printing contact constraint information to the data (.dat) file allows you to review the status of the debonding contact pair at the beginning of the analysis. By printing detailed contact conditions to the message (.msg) file, you can track the incremental behavior of the advancing crack front during the analysis. For more information about these output requests, see About Output.
You can add a small clearance to the initially unbonded portion of the debonding contact pair (Contact Initialization for Contact Pairs in Abaqus/Standard). The small clearance will help to eliminate unnecessary severe discontinuity iterations during incrementation as the crack begins to progress.
Do not use tie MPCs (General Multi-Point Constraints) for the secondary surface in a debonding contact pair. Abaqus is unable to resolve the overconstraint presented by the MPC and the debonded contact state.
You must have continuous main debonding surfaces.
You may be able to help the analysis converge by adding geometric nonlinearity (even if small-sliding is used for the debonding contact pair). For more information, see Geometric Nonlinearity.
For two-dimensional models with contact pairs involving higher-order underlying elements, the initially unbonded portion must extend over complete element faces. In other words, the crack tip in a two-dimensional, higher-order model must start at a corner node on the quadratic secondary surfaces. The crack tip must not start at a midside node.
When the surface-to-surface contact formulation is used, at least two rows of elements should be used behind the crack front.

Tips for Using the VCCT Criterion in Abaqus/Explicit

Crack propagation problems using the VCCT criterion analyzed in Abaqus/Explicit benefit from the robustness of the general contact algorithm in the context of an explicit time integrator. Nevertheless, as is the case in Abaqus/Standard, these analyses remain challenging given the discontinuous nature of the fracture phenomenon. The following tips will help you create a successful Abaqus/Explicit model:

Dynamic effects are of utmost relevance when assessing the results from a debonding analysis using the VCCT criterion. In most cases experimental and/or theoretical data are available in quasi-static settings. You must ensure that the Abaqus/Explicit analysis generates low ratios of kinetic energy to internal energy (1% or less). In practical terms this requirement often translates into avoiding the use of mass scaling in the model. Use smooth amplitudes to drive the loading to help reduce the kinetic energy in the model. Running the analysis over a longer period of time will not help in most cases because bond breakage is an inherently fast and localized process.
If appropriate, use damping-like behavior in the materials associated with the debonding plates to reduce dynamic vibrations. Unlike Abaqus/Standard, where a pure static equilibrium is achieved at the end of a converged increment, in Abaqus/Explicit the bond breakage at a given location is associated with a dynamic overshoot beyond the static equilibrium position. If the vibrations are significant (kinetic energy is clearly observable), the dynamic overshoot at nodes behind the crack tip may lead to premature debonding of the crack tip.
To maximize the accuracy of the debonding simulation, use quad meshes between the secondary and main surfaces of the debonding surfaces. Avoid using elements with aspect ratios greater than 2. In most cases mesh refinement will help with obtaining a realistic result.
Highly mismatched critical energy values between modes tend to induce crack propagation in continuously changing directions in a manner that may be unstable and unrealistic, particularly for modes II and III. Do not use such values unless experimental data suggest so.
Use frequent field output requests to evaluate the debonding evolution as the analysis progresses. In some cases this can point to nontrivial modeling deficiencies that are difficult to identify from a simple data check analysis.
Avoid the use of other constraints involving nodes on both surfaces of the debonding interface because the cohesive contact forces will compete with the constraint forces to achieve global equilibrium. Bond breakage might be hard to interpret in these cases.

Comparing VCCT and Cohesive Elements

Using VCCT to solve delamination problems is very similar to using cohesive elements in Abaqus. Table 2 describes the advantages and disadvantages of the two approaches.

For an example of the use of cohesive elements, see Delamination analysis of laminated composites. This example also shows the effect of viscous regularization on the predicted force-displacement response.

Table 2. Comparing VCCT and cohesive elements.
VCCT	Cohesive Elements
Simulation (mechanics)-driven crack propagation along a known crack surface.	Simulation (mechanics)-driven crack propagation along a known crack surface. However, cohesive elements can also be placed between element faces as a mechanism for allowing individual elements to separate.
Models brittle fracture using LEFM only.	Model brittle or ductile fracture for LEFM or EPFM. Very general interaction modeling capability is possible.
Uses a surface-based framework. Does not require additional elements.	Require definition of the connectivity and interconnectivity of cohesive elements with the rest of the structure. For accuracy, the mesh of cohesive elements may need to be smaller than the surrounding structural mesh and the associated cohesive zone. As a result, cohesive elements may be more expensive.
Requires a pre-existing flaw at the beginning of the crack surface. Cannot model crack initiation from a surface that is not already cracked.	Can model crack initiation from initially uncracked surfaces. The crack initiates when the cohesive traction stress exceeds a critical value.
Crack propagates when strain energy release rate exceeds fracture toughness.	Crack propagates according to cohesive damage model, usually calibrated so that the energy released when the crack is fully open equals the critical strain energy release rate.
Multiple crack fronts/surfaces can be included.	Multiple crack fronts/surfaces can be included.
In Abaqus/Standard crack surfaces are rigidly bonded when uncracked.	Crack surfaces are joined elastically when uncracked in Abaqus/Standard.
Requires user-specified fracture toughness of the bond.	Require user-specified critical traction value and fracture toughness of the bond, as well as elasticity of the bonded surface.

Measuring the Critical Strain Energy Release Properties for VCCT

You must obtain the critical strain energy release properties of the bonded surfaces for VCCT. The procedure to obtain the critical strain energy release properties is beyond the scope of this guide; however, you can refer to the following ASTM test specifications for guidance:

ASTM D 5528-94a, Standard Test Method for Mode I Interlaminar Fracture Toughness of Unidirectional Fiber-Reinforced Polymer Matrix Composites
ASTM D 6671-01, Standard Test Method for Mixed Mode I-Mode II Interlaminar Fracture Toughness of Unidirectional Fiber-Reinforced Polymer Matrix Composites
ASTM D 6115-97, Standard Test Method for Mode I Fatigue Delamination Growth Onset of Unidirectional Fiber-Reinforced Polymer Matrix Composites

These test specifications can be found in the Annual Book of ASTM Standards, American Society for Testing and Materials, vol. 15.03, 2000.

Initial Conditions

Initial contact conditions are used to identify which part of the secondary surface is initially bonded, as explained earlier.

Boundary Conditions

Boundary conditions should not be applied to any of the nodes on the main or secondary crack surfaces, but they can be used to load the structure and cause crack propagation. Boundary conditions can be applied to any of the displacement degrees of freedom in a crack propagation analysis (Boundary Conditions). In a fatigue crack growth analysis, prescribed boundary conditions must have an amplitude definition that is cyclic over the step: the start value must be equal to the end value (see Amplitude Curves).

Loads

The following types of loading can be prescribed in a crack propagation analysis:

Concentrated nodal forces can be applied to the displacement degrees of freedom (1–6); see Concentrated Loads.
Distributed pressure forces or body forces can be applied; see Distributed Loads. The distributed load types available with particular elements are described in Abaqus Elements Guide.

For a fatigue crack growth analysis each load must have an amplitude definition that is cyclic over the step: the start value must be equal to the end value (see Amplitude Curves).

Predefined Fields

The following predefined fields can be specified in a crack propagation analysis, as described in Predefined Fields:

Although temperature is not a degree of freedom in stress/displacement elements, nodal temperatures can be specified as predefined fields. The specified temperature affects temperature-dependent critical stress and crack opening displacement failure criteria, if specified.
The values of user-defined field variables can be specified. These values affect field-variable-dependent critical stress and crack opening displacement failure criteria, if specified.

The temperatures and user-defined field variables on secondary and main surfaces are averaged to determine the critical stresses and crack opening displacements.

In a fatigue crack growth analysis, the temperature values specified must be cyclic over the step: the start value must be equal to the end value (see Amplitude Curves). If the temperatures are read from the results file, you should specify initial temperature conditions equal to the temperature values at the end of the step (see Initial Conditions). Alternatively, you can ramp the temperatures back to their initial condition values, as described in Predefined Fields.

Material Options

Any of the mechanical constitutive models in Abaqus/Standard can be used to model the mechanical behavior of the cracking material. See Abaqus Materials Guide.

Elements

Regular, rectangular meshes give the best results in crack propagation analyses. Results with nonlinear materials are more sensitive to meshing than results with small-strain linear elasticity.

First-order elements generally work best for crack propagation analysis.

Line spring elements cannot be used in crack propagation analysis.

The VCCT, enhanced VCCT, and fatigue crack growth criteria not only support two-dimensional models (planar and axisymmetric) but also three-dimensional models with contact pairs involving first-order underlying elements (solids, shells, and continuum shells). In Abaqus/Standard use of the VCCT or enhanced VCCT criterion in two-dimensional models with contact pairs involving higher-order underlying elements is limited to crack fronts that are aligned with the corner nodes of the higher-order element faces. Use of the fatigue crack growth criterion with contact pairs involving higher-order underlying elements is not supported.

Output

Unless otherwise stated, the following discussions in this section are applied only to the critical stress, critical crack opening displacement, and crack length versus time criteria.

At the start of an analysis Abaqus/Standard will scan the partially bonded surfaces and identify all of the crack tips that are present in the model. The initial contact status of all of the secondary surface nodes is printed in the data (.dat) file. At this stage Abaqus/Standard will explicitly identify all the crack tips and mark them as crack 1, crack 2, etc. The secondary and main surfaces that are associated with these cracks are also identified.

The initial contact status of all of the secondary surface nodes is also printed in the data (.dat) file for the VCCT, enhanced VCCT, and fatigue crack growth criteria.

Printing Crack Propagation Information to the Data File

By default, crack propagation information will be printed to the data file during the analysis. For each crack that is identified Abaqus/Standard will print out the initial and current crack-tip node numbers, accumulated incremental crack length (distance from the initial crack tip to the current crack tip, measured along the secondary surface), and the current value of the user-specified fracture criterion used. Crack propagation information cannot be printed to the data file in Abaqus/Explicit.

For example, if the crack opening displacement criterion is used, the printed output during the analysis will appear as follows in the data file:

     CRACK TIP LOCATION AND ASSOCIATED QUANTITIES
CRACK  SECONDARY   MAIN  INITIAL  CURRENT  CUMULATIVE  CRITICAL
NUMBER SURFACE SURFACE CRACKTIP CRACKTIP INCREMENTAL COD
                       NODE #   NODE #   LENGTH
    $⋮$        $⋮$        $⋮$          $⋮$          $⋮$          $⋮$           $⋮$

Input File Usage

DEBOND, SECONDARY=secondary, MAIN=main

Writing Crack Propagation Information to the Results File

In Abaqus/Standard you can choose to write the crack propagation information to the results (.fil) file.

Input File Usage

DEBOND, SECONDARY=secondary, MAIN=main, OUTPUT=FILE

Writing Crack Propagation Information to Both the Data File and the Results File

In Abaqus/Standard you can write the crack propagation information to both the data and the results files.

Input File Usage

DEBOND, SECONDARY=secondary, MAIN=main, OUTPUT=BOTH

Controlling the Output Frequency

In Abaqus/Standard you can control the output frequency in increments. By default, the crack-tip location and associated quantities will be printed every increment. Specify an output frequency of 0 to suppress crack propagation output.

Input File Usage

DEBOND, SECONDARY=secondary, MAIN=main, FREQUENCY=f

Output Variables

The following bond failure quantities can be requested as surface output (see Output to the Data and Results Files, Abaqus/Standard Output Variable Identifiers, and Abaqus/Explicit Output Variable Identifiers) for all fracture criteria:

DBT: The time when bond failure occurred. For the VCCT, enhanced VCCT, and fatigue crack growth criteria, this is the time when debonding initiates.
DBSF: Fraction of stress at bond failure that still remains.
DBS: All components of remaining stress in the failed bond.
DBS1i: 1i component of stress in the failed bond that remains ( $i = 1, 2$ ).

For the VCCT, enhanced VCCT, and fatigue crack growth criteria, the following additional variables can be also requested as surface output (see Output to the Data and Results Files):

CSDMG: Overall value of the scalar damage variable.
BDSTAT: Bond state. The bond state varies between 1.0 (fully bonded) and 0.0 (fully unbonded).
OPENBC: Relative displacement behind crack when the fracture criterion is met.
CRSTS: All components of critical stress at failure
CRSTS1i: 1i component of critical stress at failure ( $i = 1, 3$ ).
ENRRT: All components of strain energy release rate.
ENRRT1i: 1i component of strain energy release rate ( $i = 1, 3$ ).
EFENRRTR: Effective energy release rate ratio, $\frac{G_{e q u i v}}{G_{e q u i v C}}$ .

Surface output requests provide the usual output of contact variables in addition to the above quantities. The bond failure quantities must be requested explicitly; otherwise, only the default output for contact will be given.

Contour Integrals

Contour integrals can be requested for two-dimensional crack propagation analyses performed using the critical stress, critical crack opening displacement, or crack length versus time fracture criteria. If the contours are chosen so that the crack tip passes through the contour, the contour value will go to zero (as it should). Therefore, in crack propagation analysis contour integrals should be requested far enough from the crack tip that the crack tip does not pass through the contour, which is easily done by including all nodes along the bond surface in the crack-tip node set specified. See Contour Integral Evaluation for details on contour integral output.

Input File Template

Abaqus/Standard Analysis

HEADING
…
BOUNDARY
Data lines to specify zero-valued boundary conditions
INITIAL CONDITIONS, TYPE=CONTACT (, NORMAL)
Data lines to specify initial conditions
SURFACE, NAME=secondary
Data lines to define secondary surface
SURFACE, NAME=main
Data lines to define main surface
**
CONTACT PAIR
secondary, main
** 
STEP (, NLGEOM)
STATIC or VISCO or COUPLED TEMPERATURE-DISPLACEMENT
DEBOND, SECONDARY=secondary, MAIN=main
Data lines to define debonding amplitude curve
FRACTURE CRITERION, TYPE=type, DISTANCE or NSET
Data lines to define fracture criterion
BOUNDARY
Data lines to define zero-valued or nonzero boundary conditions
CLOAD and/or DLOAD and/or TEMPERATURE and/or FIELD
Data lines to define loading
**
CONTOUR INTEGRAL, CONTOURS=n, TYPE=type
**Contour integrals can be requested in a two-dimensional crack propagation analysis
CONTACT PRINT
DBT, DBSF, DBS
EL PRINT
JK,
END STEP
**
STEP 
DIRECT CYCLIC, FATIGUE
DEBOND, SECONDARY=secondary, MAIN=main
FRACTURE CRITERION, TYPE=FATIGUE
Data lines to define material constants used in Paris law and  fracture criterion
BOUNDARY
Data lines to define zero-valued or nonzero cyclic  boundary conditions
CLOAD and/or DLOAD and/or TEMPERATURE and/or FIELD
Data lines to define cyclic loading
**
END STEP
**

Abaqus/Explicit Analysis

HEADING
…
BOUNDARY
Data lines to specify zero-valued boundary conditions
SURFACE, NAME=secondary
Data lines to define secondary surface
SURFACE, NAME=main
Data lines to define main surface
**
CONTACT CLEARANCE, NAME=clearance_name,
SEARCH NSET=initially_bonded_nodeset_name
SURFACE INTERACTION, NAME=interaction_name
COHESIVE BEHAVIOR
Data lines to specify elastic behavior
FRACTURE CRITERION, TYPE=VCCT, MIXED MODE BEHAVIOR=BK 
**
STEP
DYNAMIC, EXPLICIT
CONTACT
CONTACT CLEARANCE ASSIGNMENT
Data lines to assign a clearance name to a surface pair
CONTACT PROPERTY ASSIGNMENT
Data lines to assign a surface interaction to a surface pair
END STEP
**

References

Benzeggagh, M., and M. Kenane, “Measurement of Mixed-Mode Delamination Fracture Toughness of Unidirectional Glass/Epoxy Composites with Mixed-Mode Bending Apparatus,” Composite Science and Technology, vol. 56 439, 1996.
Deobald, L., G. Mabson, S. Engelstad, M. Rao, M. Gurvich, W. Seneviratne, S. Perera, T. O'Brien, G. Murri, J. Ratcliffe, C. Davila, N. Carvalho , and R. Krueger, “Guidelines for VCCT-Based Interlaminar Fatigue and Progressive Failure Finite Element Analysis,” NASA/TM-2017-219663, 2017.
Ratcliffe, J., and W. Johnston, “Influence of Mixed Mode I-Mode II Loading on Fatigue Delamination Growth Characteristics of a Graphite Epoxy Tape Laminate,” Proceedings of American Society for Composites 29th Technical Conference, 2014.
Reeder, J., S. Kyongchan, P. B. Chunchu, and D. R.. Ambur, “Postbuckling and Growth of Delaminations in Composite Plates Subjected to Axial Compression”43rd AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, Denver, Colorado, vol. 1746, p. 10, 2002.
Wu, E. M., and R. C. Reuter Jr., “Crack Extension in Fiberglass Reinforced Plastics,” T and M Report, University of Illinois, vol. 275, 1965.