3.7. Modeling Interface Delamination with Contact Elements (Debonding)

Interface delamination with contact elements is referred to as debonding. Debonding is modeled with contact elements which are bonded and have a cohesive zone material model defined. There are several advantages to using debonding to model interface delamination. Existing models with contact definitions can be easily modified to include debonding, and standard contact and debonding can be simulated with the same contact definitions.

Also see Modeling Interface Delamination with Interface Elements.

3.7.1. Analyzing Debonding

Modeling debonding with contact elements involves the same steps as any other contact analysis. (For details about setting up a contact analysis, see Surface-to-Surface Contact (Pair-Based).) If you are familiar with setting up a contact analysis, you can easily include debonding in your model. Simply add a bonded contact option and a cohesive zone material model for the contact elements.

3.7.2. Contact Elements

Debonding can be defined in any model that includes the following types of contact:

The contact elements must use bonded contact (KEYOPT(12) = 2, 3, 4, 5 or 6) with the augmented Lagrangian method or pure penalty method (KEYOPT(2) = 0 or 1). Debonding is activated by associating a cohesive zone material model (TB,CZM) with the contact elements.

3.7.3. Material Definition

3.7.3.1. Material Characteristics

The bilinear material model developed exclusively for contact elements (TB,CZM with TBOPT = CBDD or CBDE) is the recommended bilinear model for debonding. The material behavior, defined in terms of contact stresses (normal and tangential) and contact separation distances (normal gap and tangential sliding), is characterized by linear elastic loading followed by linear softening. The slope of the curve depends on contact stiffness and a debonding parameter which is defined in terms of material constants.

Two other cohesive zone material models that were developed for use with the interface elements can also be used to model contact debonding: the bilinear material model (TB,CZM with TBOPT = BILI) and the exponential material model (TB,CZM with TBOPT = EXPO).

Debonding allows three modes of separation:

  • Mode I debonding for normal separation

  • Mode II debonding for tangential separation

  • Mixed mode debonding for normal and tangential separation

Debonding is also characterized by convergence difficulties during material softening. Artificial damping is provided to overcome these problems. An option for tangential slip under compressive normal contact stress for mixed-mode debonding is also provided.

After debonding is completed, the surface interaction is governed by standard contact constraints for normal and tangential directions. Frictional contact is used if friction is specified for the contact elements.

3.7.3.2. Bilinear Material Behavior - Contact (TBOPT = CBDD and CBDE)

The cohesive zone material model with bilinear behavior (TB,CZM with TBOPT = CBDD or CBDE ) is defined as:

where:

P = normal contact stress (tension)
τy = tangential contact stress in y direction
τz = tangential contact stress in z direction
Kn = normal contact stiffness
Kt = tangential contact stiffness
un = contact gap
uy = contact slip distance in y direction
uz = contact slip distance in z direction
d = debonding parameter

To model bilinear material behavior with tractions and separation distances, use TB,CZM with TBOPT = CBDD. You also input the following material constants with the TBDATA command:

Constant Symbol Meaning
C1σmax Maximum normal contact stress [1]
C2 Contact gap at the completion of debonding
C3τmax Maximum equivalent tangential contact stress [1]
C4 Tangential slip at the completion of debonding
C5ηArtificial damping coefficient
C6βFlag for tangential slip under compressive normal contact stress; must be 0 (off)  or 1 (on)
  1. For contact elements using the force-based model (see the description of KEYOPT(3) for CONTA175 and CONTA177 ), input a contact force value for this quantity.

To model bilinear material behavior with tractions and critical fracture energies, use TB,CZM with TBOPT = CBDE. You also input the following material constants with the TBDATA command:

Constant Symbol Meaning
C1σmax Maximum normal contact stress [1]
C2Gcn Critical fracture energy density (energy/area) for normal separation [2]
C3τmax Maximum equivalent tangential contact stress [1]
C4Gct Critical fracture energy density (energy/area) for tangential slip [2]
C5ηArtificial damping coefficient
C6βFlag for tangential slip under compressive normal contact stress; must be 0 (off)  or 1 (on)
  1. For contact elements using the force-based model (see the description of KEYOPT(3) for CONTA175 and CONTA177 ), input a contact force value for this quantity.

  2. For contact elements using the force-based model (see the description of KEYOPT(3) for CONTA175 and CONTA177 ), this quantity is critical fracture energy.

Example 3.14: Defining a Cohesive Zone Material

TB,CZM,,,,CBDD      ! bilinear behavior with tractions and separation distances
TBDATA,1,C1,C2,C3,C4,C5,C6

For more information about the cohesive zone material model, see Cohesive Zone Material for Contact Elements in the Material Reference.

3.7.3.3. Bilinear Material Behavior - Interface (TBOPT = BILI)

This model follows a bilinear law for traction-separation and differs slightly from the bilinear material behavior for contact. See Cohesive Zone Material (CZM) Model in the Mechanical APDL Theory Reference for details on the differences.

To define this material, use the TB,CZM command with TBOPT = BILI. Specify the material constants as data items C1 through C6 on the TBDATA command as described in Material Constants -- Bilinear Law.

3.7.3.4. Exponential Material Behavior (TBOPT = EXPO)

This model follows an exponential law for traction separation. To define this material, use the TB,CZM command with TBOPT = EXPO. Specify the material constants as data items C1, C2, and C3 on the TBDATA command as described in Material Constants -- Exponential Law.

3.7.3.5. Post-Debonding Behavior at the Contact Interface

For the cohesive zone materials with bilinear material behavior (TBOPT = CBDD, CBDE or BILI on the TB command), you can specify that the cohesive zone interface be “healed” if the surfaces come into contact again after debonding.

To activate this option, use the TBFIELD,CYCLE command to define the CZM material as a function of healing cycle number. For more information, see Post-Debonding Behavior of Cohesive Zone Material in the Contact Technology Guide.

3.7.4. Result Output

All applicable output quantities for contact elements are also available for debonding:

Output Quantities Symbol Meaning
PRESPNormal contact stress [1]
SFRICτt Tangential constant stress [1]
TAUR and TAUSτy and τz Components (tangential constant stress) [1]
GAPun Contact gap
SLIDEut Tangential slip
TASR and TASSuy and uz Components (tangential slip)
  1. For contact elements using the force-based model (see the description of KEYOPT(3) for CONTA175 and CONTA177 ), this output item is reported as a force quantity.

Debonding specific output quantities are also available and are output as NMISC data. The output quantities vary based on the CZM model used, as outlined in the tables below.

TB,CZM with TBOPT = CBDD or CBDE (Bilinear Law for Contact)
Output QuantitiesSymbolMeaning
DTSTART(no symbol)Debonding time history
DPARAMdn, dt, or dm Debonding parameter
DENERI and DENERIIGn and Gt Fracture energies
TB,CZM with TBOPT = BILI (Bilinear Law for Interface)
Output QuantitiesSymbolMeaning
DTSTART(no symbol)Debonding time history
DPARAMdn, dt, or dm Debonding parameter
DENER Gtotal Total debonding energy
TB,CZM with TBOPT = EXPO (Exponential Law)
Output QuantitiesSymbolMeaning
DENERGtotal Total debonding energy