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.
The following related topics are available:
Also see Modeling Interface Delamination with Interface Elements.
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.
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.
The following material-definition topics are available for modeling interface delamination (debonding) with contact elements:
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.
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) |
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] |
C2 | Gcn | Critical fracture energy density (energy/area) for normal separation [2] |
C3 | τmax | Maximum equivalent tangential contact stress [1] |
C4 | Gct | 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) |
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.
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.
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.
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.
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.
All applicable output quantities for contact elements are also available for debonding:
Output Quantities | Symbol | Meaning |
---|---|---|
PRES | P | Normal contact stress [1] |
SFRIC | τt | Tangential constant stress [1] |
TAUR and TAUS | τy and τz | Components (tangential constant stress) [1] |
GAP | un | Contact gap |
SLIDE | ut | Tangential slip |
TASR and TASS | uy and uz | Components (tangential slip) |
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 Quantities | Symbol | Meaning |
DTSTART | (no symbol) | Debonding time history |
DPARAM | dn, dt, or dm | Debonding parameter |
DENERI and DENERII | Gn and Gt | Fracture energies |
TB,CZM with
TBOPT = BILI (Bilinear Law for
Interface) | ||
---|---|---|
Output Quantities | Symbol | Meaning |
DTSTART | (no symbol) | Debonding time history |
DPARAM | dn, dt, or dm | Debonding parameter |
DENER | Gtotal | Total debonding energy |
TB,CZM with
TBOPT = EXPO (Exponential
Law) | ||
---|---|---|
Output Quantities | Symbol | Meaning |
DENER | Gtotal | Total debonding energy |