当前位置:首页 >> 能源/化工 >>

FRACTURE LECTURE OF ABAQUS


Basic Concepts of Fracture Mechanics
Lecture 1

L1.2

Overview
? Introduction

? Fracture Mechanisms
? Linear Elastic Fracture Mechanics

/>? Small Scale Yielding
? Energy Considerations

? The J-integral
? Nonlinear Fracture Mechanics ? Mixed-Mode Fracture ? Interfacial Fracture ? Creep Fracture

? Fatigue

Modeling Fracture and Failure with Abaqus

L1.3

Overview
? This lecture is optional.

? It aims to introduce the necessary fracture mechanics concepts and quantities that are relevant to the Abaqus functionality that is presented in the subsequent lectures.
? If you are already familiar with these concepts, this lecture may be omitted.

Modeling Fracture and Failure with Abaqus

Introduction

L1.5

Introduction
? Fracture mechanics is the field of solid mechanics that deals with the behavior of cracked bodies subjected to stresses and strains. ? These can arise from primary applied loads or secondary selfequilibrating stress fields (e.g., residual stresses).

Modeling Fracture and Failure with Abaqus

L1.6

Introduction
? Objective of fracture mechanics

? The objective of fracture mechanics is to characterize the local deformation around a crack tip in terms of the asymptotic field around the crack tip scaled by parameters that are a function of the loading and global geometry.

Modeling Fracture and Failure with Abaqus

Fracture Mechanisms

L1.8

Fracture Mechanisms
? For engineering materials, such as metals, there are two primary modes of fracture: brittle and ductile. ? Brittle fracture ? Cracks spread very rapidly with little or no plastic deformation.

? Cracks that initiate in a brittle material tend to continue to grow and increase in size provided the loading will cause crack growth.
? Ductile fracture

? Three stages: void nucleation, growth, and coalescence. ? The crack moves slowly and is accompanied by a large amount of plastic deformation. ? The crack typically will not grow unless the applied load is increased.
Modeling Fracture and Failure with Abaqus

L1.9

Fracture Mechanisms
? Brittle fracture in polycrystalline materials displays either cleavage (transgranular) or intergranular fracture. ? This depends upon whether the grain boundaries are stronger or weaker than the grains .

Cleavage fracture

Modeling Fracture and Failure with Abaqus

L1.10

Fracture Mechanisms
? Ductile fracture has a dimpled, cup-and-cone fracture appearance .

? Ductile fracture surfaces have larger necking regions and an overall rougher appearance than a brittle fracture surface.

Modeling Fracture and Failure with Abaqus

L1.11

Fracture Mechanisms
? Fracture process zone

? The fracture process zone is the region around the crack tip where dislocation motions, material damage, etc. occur.
? It is a region of nonlinear deformation. ? The fracture process zone size is characterized by

? a number of grain sizes for brittle fracture or
? either inclusion or second phase particle spacings for ductile fracture. ? Different theories have been advanced to describe the fracture process in order to develop predictive capabilities ? LEFM

? Cohesive zone models
? EPFM

? Etc.
Modeling Fracture and Failure with Abaqus

Linear Elastic Fracture Mechanics

L1.13

Linear Elastic Fracture Mechanics
? Fracture modes

? Linear Elastic Fracture Mechanics (LEFM) considers three distinct fracture modes: Modes I, II, and III
? These encompass all possible ways a crack tip can deform. ? Mode I: ? The forces are perpendicular to the crack, pulling the crack open. ? This is referred to as the opening mode.

Modeling Fracture and Failure with Abaqus

L1.14

Linear Elastic Fracture Mechanics
? Mode II:

? The forces are parallel to the crack.
? One force pushes the top half of the crack back and the other pulls the bottom half of the crack forward, both along the same line. ? This creates a shear crack: the crack slides along itself.

? This is referred to as the in-plane shear mode. ? The forces do not cause out-ofplane deformation.

Modeling Fracture and Failure with Abaqus

L1.15

Linear Elastic Fracture Mechanics
? Mode III:

? The forces are transverse to the crack.
? This causes the material to separate and slide along itself, moving out of its original plane ? This is referred to as the out-of-plane shear mode. ? The objective of LEFM is to predict the critical loads that will cause a crack to grow in a brittle material.

Modeling Fracture and Failure with Abaqus

L1.16

Linear Elastic Fracture Mechanics
? Stress intensity factor

? For isotropic, linear elastic materials, LEFM characterizes the local crack-tip stress field in the linear elastic (i.e., brittle) material using a single parameter called the stress intensity factor K.

? K depends upon the applied stress, the size and placement of the
crack, as well as the geometry of the specimen.

? K is defined from the elastic stresses near the tip of a sharp crack
under remote loading (or residual stresses).

? K is used to predict the stress state ("stress intensity") near the tip
of a crack.
? When this stress state (i.e., K) becomes critical, a small crack grows ("extends") and the material fails. ? This critical value is denoted KC and is known as the fracture toughness (it is a material property; discussed further later).

Modeling Fracture and Failure with Abaqus

L1.17

Linear Elastic Fracture Mechanics
? Asymptotic crack tip solutions

? The stress and strain fields in the vicinity of the crack tip are expressed in terms of asymptotic series of solutions around the crack tip.
? They are valid only is a small region near the crack tip. ? This size of this region is quantified by small scale yielding assumptions (discussed later). ? The stress intensity factor is the parameter that relates the local crack-tip fields with the global aspects of the problem.

Modeling Fracture and Failure with Abaqus

L1.18

Linear Elastic Fracture Mechanics
? The leading-order terms of the asymptotic solution are:

? ij (r ,? ) ?
where

KI K K fijI (? ) ? II fijII (? ) ? III fijIII (? ), 2? r 2? r 2? r

x2

r

?

x1

r is the distance from the crack tip,
? = atan(x2/x1),

KI is the Mode I (opening) stress intensity factor,
KII is the Mode II (in-plane shear) stress intensity factor, KIII is the Mode III (transverse shear) stress intensity factor, and the

fija define the angular variation of the stress for mode a.

Modeling Fracture and Failure with Abaqus

L1.19

Linear Elastic Fracture Mechanics
? Crack-tip singularity

? The predicted stress state at the crack tip in a linear elastic (brittle) material possesses a square-root singularity:

??

1 . r

? In reality, the crack tip is surrounded by the fracture process zone where plastic deformation and material damage occur. ? Inside this zone, the LEFM solution is not valid.

? Outside of this zone (i.e., sufficiently "far" from the fracture process zone), the LEFM is accurate provided the plastic/damage zone is “small enough.”
? This is called small-scale yielding (discussed further later).

Modeling Fracture and Failure with Abaqus

L1.20

Linear Elastic Fracture Mechanics
? Some comments on fracture toughness

? Fracture toughness is strongly dependent on temperature.
Fracture toughness

Temperature

? The brittle-ductile transition temperature range depends on the material.

? For many common metals it may lie within the reasonable operating temperature range for the design, so the temperature dependence of the fracture toughness must be considered.

Modeling Fracture and Failure with Abaqus

L1.21

Linear Elastic Fracture Mechanics
? Experimentally, the fracture toughness KC is a function of specimen thickness. ? Since plane strain gives the practical minimum value of KC …

? The plane strain value is usually the value that is determined experimentally. ? However, if the application is fracture of thin sheets of material, KC values somewhere between the plane stress and plane strain values may be appropriate.
Fracture toughness

KC
Thickness



Modeling Fracture and Failure with Abaqus

L1.22

Linear Elastic Fracture Mechanics
? Aside from temperature and thickness, the fracture toughness is also a function of the crack extension. ? The fracture toughness as a function of crack extension is called the resistance curve (shown below).
ductile

Variation in fracture toughness with crack growth is Kr(Da):

Kr(0)= KC
brittle

? The resistance curve is used to predict crack growth stability.

Modeling Fracture and Failure with Abaqus

L1.23

Linear Elastic Fracture Mechanics
? Crack growth and stability

? The condition for continued crack growth for a crack length a + Da is

Kapplied ? K R (Da).
? The condition for stable continued crack growth is

?K applied ?a
load

?

dK R . d Da

Modeling Fracture and Failure with Abaqus

Small-Scale Yielding

L1.25

Small-Scale Yielding
? Small-scale yielding (SSY) means the region of inelastic deformation at the crack tip is contained well within the zone dominated by the LEFM asymptotic solution. ? For LEFM to be valid, there must be an annular region around the crack tip in which the asymptotic solution to the linear elasticity problem gives a good approximation to the complete stress field.

Plastic zone

K-dominated zone

Transition zone

Modeling Fracture and Failure with Abaqus

L1.26

Small-Scale Yielding
? The size of the process zone and the plastic region must be sufficiently small so that this is true. Typical shapes of plastic zones follow:

plane strain

plane stress (diffuse)

plane stress (Dugdale)

Modeling Fracture and Failure with Abaqus

L1.27

Small-Scale Yielding
? We can estimate the plastic zone size, rp, by setting ?22 = ?0 in the LEFM asymptotic solution, where ?0 is the yield stress. This gives (for Mode I)

rp ?

1 2?

? KI ? 1 ? KI ? ? ? ? ? ? . 6 ? ?0 ? ? ?0 ?

2

2

? Since the tractions across the boundary of the plastic zone have no net force or moments (St. Venant’s principle), the effect on the elastic field surrounding the plastic zone decays rapidly with distance from the boundary, becoming negligible at ~3rp.
? LEFM predicts infinite stress at the crack tip —obviously this is unrealistic. ? But we can use LEFM results if the region of inelastic deformation near the crack tip is small enough that there is a finite zone outside this region where the LEFM asymptotic solution is accurate.

Modeling Fracture and Failure with Abaqus

L1.28

Small-Scale Yielding
? If a is a characteristic dimension in the problem, such as remaining ligament size or thickness or crack length, then, to have a finite zone rK in which the K-field dominates, we need

1 ? K IC ? a / 5 ? rK ? 3rp ? ? ? 2 ? ?0 ?
or

2

? K IC ? a ? 2.5 ? ? . ? ?0 ?

2

ASTM Standard for validity of LEFM

? This is the limit on specimen size in ASTM Standard E-399 for a valid KIC test.

? KIC is KC (the fracture toughness) in Mode I.
? The fracture toughness represents the critical value of K required to initiate crack growth.
Modeling Fracture and Failure with Abaqus

L1.29

Small-Scale Yielding
? For some typical metal materials rp is calculated by matching the yield stress to the Mises stress of the K field and the minimum characteristic length is calculated using the ASTM standard limit. ? For materials with high fracture toughness the size of the specimen for a valid fracture test is very large. T
(? C)

Material
A061-T651 (Al) A075-T651 (Al)

?0
(MPa)

KIC
(MN/m3/2)

rp
(mm)

Characteristic dimension (mm)

20 20

269 620

33 36

5 0.35

38 8.4

AISI 4340 (Steel)
A533-B (Steel)

0
93

1500
620

33
200

0.05
11

1.2
260

Modeling Fracture and Failure with Abaqus

Energy Considerations

L1.31

Energy Considerations
? Energy principles play an important role in studying crack problems.

? This is motivated by the fact that crack propagation always involves dissipation of energy. Sources of energy dissipation include:
? Surface energy, plastic dissipation, etc. ? By considering fracture from an energetic point of view, crack growth criteria can be postulated in terms of energy release rates. ? This approach offers an alternative to the K-based fracture criteria discussed earlier and reinforces the connection between global and local fields in fracture problems.
? The energy release rate is a global parameter while the stress intensity factor is a local crack-tip parameter.

Modeling Fracture and Failure with Abaqus

L1.32

Energy Considerations
? The energy available to grow a crack is defined as

G?-

? ( PE ) , ?a Loads

where PE is the potential energy and G is the Energy Release Rate. ? We consider the difference in the energy for two essentially identical specimens, one with crack length a, the other with crack length a + Da. ? The area under the loaddisplacement curve gives -PE for elastic materials.

Modeling Fracture and Failure with Abaqus

L1.33

Energy Considerations
? For isotropic linear elastic materials, one can show that

1 - v2 2 G? K for plane strain E
and

K2 G? for plane stress. E
? In a three-dimensional body under general loading that contains a crack with a smoothly changing crack-tip line, the energy release rate (assuming linear elasticity) per unit crack front length is

? Thus, we see the stress intensity factors are directly related to the energy release rate associated with infinitesimal crack growth in an isotropic linear elastic material.
Modeling Fracture and Failure with Abaqus

1 - v2 2 1 2 2 G? ( K I ? K II )? K III . E 2G

L1.34

Energy Considerations
? Initiation of crack growth in SSY

? The necessary condition for crack growth expressed in terms of the energy release rate is G ? GC.

? GC is a material property and represents the energy per unit crack
advance going into:

? the formation of new surfaces,
? the fracture process, and

? plastic deformation.
? As noted earlier, for linear elastic materials, G and K are related.

? This leads to an alternative condition for K ? KC.
? Recall KC is the fracture toughness of the material.

Modeling Fracture and Failure with Abaqus

The J-integral

L1.36

The J-integral
? The J-integral is used in rate-independent quasi-static fracture analysis to characterize the energy release associated with crack growth. ? It can be related to the stress intensity factor if the material response is linear. ? As will become apparent in the next section, it also has the advantage that it provides a method for analyzing fracture in nonlinear materials.

Modeling Fracture and Failure with Abaqus

L1.37

The J-integral ? J is defined as follows:

? ? ?u J ? ? Wn1 - i ? ij n j ? ds ?x1 ? ??

?

x2

x1
? It is path independent when contours are taken around a crack tip.
? The definition of J assumes: ? The material is homogeneous in the crack direction.

? The material is elastic.
? For linear elastic materials, the value of J is equal to the energy release rate associated with crack advance:

J ?G
Modeling Fracture and Failure with Abaqus

L1.38

The J-integral ? J in small-scale yielding
? Choose ?, the contour for J, to fall entirely within the annular region in which the K fields dominate.

3rp
? The integrand for J can be evaluated directly in terms of the (known) K fields. Direct calculation for Mode I in a linear elastic material gives

1 - v2 2 J ?G ? K I for plane strain and E 1 J ? G ? K I2 for plane stress. E
Modeling Fracture and Failure with Abaqus

Nonlinear Fracture Mechanics

L1.40

Nonlinear Fracture Mechanics
? LEFM applies when the nonlinear deformation of the material is confined to a small region near the crack tip. ? For brittle materials, it accurately establishes the criteria for failure. ? However, severe limitations arise when the region of the material subject to plastic deformation before a crack propagates is not negligible. ? Nonlinear fracture mechanics attempts to extend LEFM to consider inelastic effects.

? The theory is sometimes called Elastic-Plastic Fracture Mechanics (EPFM). ? However, the theory is not based on an elastic-plastic material model, but rather a nonlinear elastic material.

Modeling Fracture and Failure with Abaqus

L1.41

Nonlinear Fracture Mechanics

n

? Consider a material that has a power-law hardening form,

?? ? e ?a? ? , e0 ? ?0 ?
where ?0 is the effective yield stress, e0 = ?0 / E is the associated yield strain, E is Young's modulus, and a and n are chosen to fit the stressstrain data for the material.
Modeling Fracture and Failure with Abaqus

n

L1.42

Nonlinear Fracture Mechanics
? For such a material, Hutchinson, Rice, and Rosengren (extended to mixed mode loading by Shih) showed that the near-tip fields have the form
Loading parameter is J

? ij ? ? 0 ? e ij ? e 0 ?
?

?

? a? e I r ? 0 0 n ? J

J

1 ? n?1

? ij (? ), ?

? a? e I r ? 0 0 n ?

n ? n?1

?ij (? ), e
n ? n?1

? J ?i ? ae 0 r ? ui - u ? a? e I r ? 0 0 n ?

?i (? ). u

?i is the displacement relative to the displacement of the crack ? Here ui - u ?i. These fields are commonly referred to as the HRR crack-tip fields. tip, u
Modeling Fracture and Failure with Abaqus

L1.43

Nonlinear Fracture Mechanics
? Why not elastic-plastic?

? The HRR field assumes a nonlinear elastic power law material:
?? ? e ?a ? ? e0 ??0 ? ? Under monotonic loading, this nonlinear elastic material can be matched to the behavior of an elastic-plastic material whose hardening behavior is accurately modeled by a power law.
n

? Thus, evaluating J allows us to characterize the strength of the singularity in the crack-tip region in an elastic-plastic material subjected to monotonic loading.
Modeling Fracture and Failure with Abaqus

L1.44

Nonlinear Fracture Mechanics
? In unloading situations, the HRR fields do not describe the state around the crack tip, and hence J does not characterize the strength of the stress state ahead of a crack tip for plastic materials. Use caution when: ? The loading is not monotonic and an incremental plasticity material is used ? Crack growth occurs under monotonic loading (individual material particles may unload even when the overall structure is being loaded).

? The HRR solution:
? Gives the leading term in an asymptotic expansion of the deformation around the crack tip for a power law material; and ? Does not take into account finite-strain effects.

Modeling Fracture and Failure with Abaqus

L1.45

Nonlinear Fracture Mechanics
? Some comments on the HRR fields

? The HRR fields, thus, describe the near-tip crack fields in terms of J.

? J gives the strength of the near-tip singularity in any power-law material
(nonlinear elastic or plastic) solid
? Recall that in LEFM K plays this role in linear elastic materials.

? J-based fracture mechanics is applied in much the same way as LEFM.
? Crack growth initiates when J reaches a critical value: J ? JC . ? To apply the theory, must ensure conditions for J-dominance are satisfied (discussed next).

Modeling Fracture and Failure with Abaqus

L1.46

Nonlinear Fracture Mechanics ? J-dominance

? J-dominance refers to situations when J can be used as a method of
predicting fracture.

? In general, J is an adequate characterization when there exists a state of high triaxial tension (high triaxiality) ahead of the crack tip. ? High triaxiality ahead of the crack tip leads to low fracture toughness. ? Examples: states of small-scale and well-contained yielding (where the plastic zone is surrounded by an elastic zone):
? Deeply notched bend specimen c ?d
d c

Modeling Fracture and Failure with Abaqus

L1.47

Nonlinear Fracture Mechanics
? In some situations the crack-tip stress field does not exhibit high triaxiality. ? Example: large-scale yielding (the plastic zone extends to the free boundaries of the body): ? Fully plastic flow of single-edge cracked specimens under tension loading

? Shallow cracks under bending
? Center-cracked panel

? A two-parameter approach can be used to extend the fracture characterization to such cases (discussed next).

Modeling Fracture and Failure with Abaqus

L1.48

Nonlinear Fracture Mechanics
? Two-parameter fracture mechanics ? The Williams’ expansion of the Mode I stress field about a sharp crack in a linear elastic body with respect to r, the distance from the crack tip, is

? ij (r ,? ) ?

KI fij (? ) ? T ?1i?1 j ? O(r1/2 ). 2? r

? The T-stress thus represents a stress parallel to the crack faces. ? The magnitude of the T-stress affects the size and shape of the plastic zone and the region of tensile triaxiality ahead of the crack tip. ? For positive T-stress, J-dominance exists and a single parameter J can be used for a fracture criterion. ? For negative T-stress, a two-parameter approach (J, T) is required to characterize the stress fields.

Modeling Fracture and Failure with Abaqus

Mixed-Mode Fracture

L1.50

Mixed-Mode Fracture
? Under general loading almost all theories for the direction of crack growth assume or predict that the continued crack growth will be with KII = 0. ? Can assume that macroscopic cracks growing with continuously turning tangents will advance straight ahead, presumably under Mode I conditions.

? The crack curvature will evolve in such a way as to maintain this in response to the loading. ? If the loading changes such that the local crack-tip stress field experiences a large change in local stress intensities, mixed-mode fracture will occur.

Modeling Fracture and Failure with Abaqus

L1.51

Mixed-Mode Fracture
? Different criteria for homogeneous, isotropic linear elastic materials have been proposed, including:

? The maximum tangential stress criterion. ? The maximum energy release rate criterion. ? The KII = 0 criterion.
? Although all three imply that KII = 0 as the crack extends, they predict slightly different angles for crack initiation.
Comparison of predictions of crack propagation direction for different ratios of KII / KI

Modeling Fracture and Failure with Abaqus

Interfacial Fracture

L1.53

Interfacial Fracture
? Many engineering applications involve bonded materials.

? Examples:
? adhesive joints;

? protective coatings;
? composite materials;

? etc.
? Engineers must be able to predict the strength of the bond. ? Interfacial fracture mechanics provides a method by which to do this. ? It extends LEFM to predict the behavior of cracks between two linear elastic materials.

Modeling Fracture and Failure with Abaqus

L1.54

Interfacial Fracture
? Once a crack has started to grow in an isotropic, homogeneous material, it generally does so in an opening mode; that is, in Mode I. ? A crack lying on an interface can kink off the interface and grow under Mode I conditions, or it can grow along the interface under mixed mode conditions. ? Whether the crack kinks off the interface or propagates along it is frequently determined through energy considerations.

Modeling Fracture and Failure with Abaqus

L1.55

Interfacial Fracture
? If the crack kinks off the interface, the fact that there is an interface is important only in how it influences the stress and strain fields. ? If the crack grows along the interface, it grows under mixed mode conditions due to material asymmetry and possibly (though not necessarily) under mixed remote loading conditions. ? In such situations the conditions for crack growth depend on the interface properties. It is not sufficient to define crack initiation and growth criterion based on the conventional fracture toughness, KC.

? Specifically KC = KC (?).
? Toughness depends strongly on the mode mixity ?.

Modeling Fracture and Failure with Abaqus

L1.56

Interfacial Fracture
? Asymptotic fields

? The asymptotic stress field for an interfacial crack between linear elastic materials is given by

? K * ie ? ? ij (? , e ) ? ? ij ? Re ? r ? ? 2? r ? where K* = K1 ? iK2 is the complex stress intensity factor (i.e., it has real ? ij ?? , e ? is a complex function of the angle and imaginary parts) and ?
and material mismatch parameter e :

e?

? (? - 1) - ?2 (?1 - 1) 1 1- ? log , where ? ? 1 2 , and 2? 1? ? ?1 (? 2 ? 1) ? ?2 (?1 ? 1)
for plane stress for plane strain, axi, 3D

? 3 -? ? ? ? ? 1 ?? ? ?3 - 4?

Modeling Fracture and Failure with Abaqus

L1.57

Interfacial Fracture
? The complex exponent rie indicates that the stresses will oscillate near the crack tip:

? Both the stresses and crack opening displacements will oscillate wildly as the crack tip is approached. ? At some distance ahead of the crack tip, the fields settle down.
? The fracture criterion should be measured at this point. Provided the location of this point is the same in different specimens, a fracture criterion is valid.

Modeling Fracture and Failure with Abaqus

Creep Fracture

L1.59

Creep Fracture
? High-temperature fracture

? For temperatures above 0.3?M (where ?M is the melting temperature on an absolute scale), metals will typically creep. ? In plastics creep can occur even at room temperature.
? There are typically two mechanisms that are active in creep fracture:

? Blunting of the crack tip due to a relaxing stress field.
? This tends to retard crack growth.

? Accumulation of creep damage (microcracks, void growth, and coalescence). ? This enhances crack growth.
? Steady-state creep crack growth occurs when the two effects balance one another.

Modeling Fracture and Failure with Abaqus

L1.60

Creep Fracture
? The stress state around a crack tip in a material that can creep is more complicated than for the corresponding plasticity problem. ? Because of the time-dependent effects there is no one parameter that can characterize the stress state around the crack tip for all possibilities.

? This makes measuring the relevant parameters more difficult.
? Hence, creep fracture is not as well established as elastic-plastic fracture.
Initially, the crack-tip field is the elastic field. ? cr ) ? O(e ? el ) around the Stationary crack: O(e ? el ) ? O(e ? cr ) crack tip (RR field); around this field O(e (K field). ? el ) ? O(e ? cr ) Growing crack: region develops where O(e (HR field), which is in turn surrounded by the RR field. Eventually the HR field envelops the RR field (which ultimately disappears).

Modeling Fracture and Failure with Abaqus

L1.61

Creep Fracture
? Contour integrals

? The contour integral for creep fracture is called the C(t)-integral.
? It plays an analogous role to the J-integral in the context of timedependent creep fracture. ? Its development assumes a power law creep material:

? ?e ? el ? e ? cr e

? The C(t)-integral is proportional to the rate of growth of the crack-tip creep zone for a stationary crack under small-scale creep conditions:

?? ? ? ? ? e0 ? ? E ? ?0 ? ? ?

n

? Under steady-state creep conditions, when creep dominates throughout the specimen, C(t) becomes path independent and is known as C*.

?j ? ?u ? n ?ij n1 - ni ? ij C (t ) ? ? ije ds. ? ? ? r ?0 ? n ? 1 ?x1 ?

?

Modeling Fracture and Failure with Abaqus

L1.62

Creep Fracture
? Asymptotic fields for stationary crack

? The near tip stress and strain fields were obtained by Riedel and Rice in terms of C(t). They are known as the RR fields and are analogous to the HRR fields in power law hardening plasticity.
C(t) acts like a time-dependent
loading parameter

? ij ? ? 0 ?
cr ?ij e

? C (t ) ? ? e ? I r ? 0 0 n ?

1 ? n?1

? ij (? , n) ? ?ij (? , n) e

? C (t ) ?0 ? ?e ? ? e ? I r ? 0 0 n ?

n ? n?1

Crack tip fields are similar to those for an elastic-plastic material

?ij (? , n) is approximately Here In is a function of n and the magnitude of ? 1.

Modeling Fracture and Failure with Abaqus

L1.63

Creep Fracture
? Small-scale vs. extensive creep

? For the case of no crack growth the loading parameters that characterize the crack-tip fields are reasonably well understood.
? Under small-scale creep conditions with secondary creep, K is the loading parameter characterizing the crack-tip field. ? For extensive secondary creep C* is a loading parameter characterizing the crack-tip field upon which a fracture criterion may be based. ? Suitable criteria for crack extension that will predict an initiation time for crack growth for general cases are not yet available.
Modeling Fracture and Failure with Abaqus

??

K ? (? ) ? r

Small-scale creep

creep zone

Extensive creep

Fatigue

L1.65

Fatigue
? Fatigue is a special kind of failure in which cracks gradually grow under a prolonged period of subcritical loading. ? It is the single most common cause of failure in metallic structures.

Damage at the ball grid array (BGA) in a solder joint after 2700 thermal loading cycles

? The Paris Law can be used to predict crack growth as a function of cycles (or time):

da ? C (DK ) n , where dN DK ? K max - K min

Modeling Fracture and Failure with Abaqus

L1.66

Fatigue
? Abaqus offers a direct cyclic low-cycle fatigue capability based on the Paris Law. ? Models progressive damage and failure both in bulk materials and at material interfaces for a structure subjected to a sub-critical cyclic loading. ? For more advanced fatigue analysis capabilities, consult www.safetechnology.com. ? fe-safe is a suite of fatigue analysis software that has a direct interface to Abaqus.

Modeling Fracture and Failure with Abaqus

Modeling Cracks
Lecture 2

L2.2

Overview
? Crack Modeling Overview ? Modeling Sharp Cracks in Two Dimensions ? Modeling Sharp Cracks in Three Dimensions ? Finite-Strain Analysis of Crack Tips ? Limitations Of 3D Swept Meshing For Fracture ? Modeling Cracks with Keyword Options

Modeling Fracture and Failure with Abaqus

Crack Modeling Overview

L2.4

Crack Modeling Overview
? A crack can be modeled as either

? Sharp
? Small-strain analysis

? Singular behavior at the crack tip
? Requires special attention

? In Abaqus, a sharp crack is modeled using seam geometry ? Blunted
? Finite-strain analysis

? Non-singular behavior at crack tip
? In Abaqus, a blunted crack is modeled using open geometry ? For example, a notch

Modeling Fracture and Failure with Abaqus

L2.5

Crack Modeling Overview
? Mesh refinement

? Crack tips cause stress concentrations.
? Stress and strain gradients are large as a crack tip is approached.

? The finite element mesh must be refined in the vicinity of the crack tip to get accurate stresses and strains. ? The J-integral is an energy measure; for LEFM, accurate J values can generally be obtained with surprisingly coarse meshes, even though the local stress and strain fields are not very accurate.
? For plasticity or rubber elasticity, the crack-tip region has to be modeled carefully to give accurate results.

Modeling Fracture and Failure with Abaqus

L2.6

Crack Modeling Overview
? The crack-tip singularity in small-strain analysis

? For mesh convergence in a small-strain analysis, the singularity at the crack tip must be considered.

? J values are more accurate if some singularity is included in the
mesh at the crack tip than if no singularity is included. ? The stress and strain fields local to the crack tip will be modeled more accurately if singularities are considered. ? In small-strain analysis, the strain singularity is:
? Linear elasticity ? ? r -? ? Perfect plasticity ? ? r -1 ? Power-law hardening ? ? r -n/(n+1)

Modeling Fracture and Failure with Abaqus

Modeling Sharp Cracks in Two Dimensions

L2.8

Modeling Sharp Cracks in Two Dimensions
? In two dimensions…

? The crack is modeled as an internal edge partition embedded (partially or wholly) inside a face.
? This is called a seam crack ? The edge along the seam will have duplicate nodes such that the elements on the opposite sides of the edge will not share nodes. ? Typically, the entire 2D part is filled with a quad or quad-dominated mesh. ? At the crack tip, a ring of triangles are inserted along with concentric layers of structured quads. ? All triangles in the contour domains must be represented as degenerated quads.
Modeling Fracture and Failure with Abaqus

L2.9

Modeling Sharp Cracks in Two Dimensions
? Example: Slanted crack in a plate

? In Abaqus/CAE a seam is defined by through the Crack option underneath the Special menu of the Interaction module.
? The seam will generate duplicate nodes along the edge.

Seam

Create face partition to represent the seam; assign a seam to the partition.

Modeling Fracture and Failure with Abaqus

L2.10

Modeling Sharp Cracks in Two Dimensions
? To define the crack, you must specify ? Crack front and the crack-tip ? Normal to the crack plane or the direction of crack advance ? The crack advance direction is called the q vector.

Select the vertex at either end as the crack front. (Repeat for the other end.)

Crack tip same as crack front in this case

The crack extension direction (q vector) defines the direction in which the crack would extend if it were growing. It is used for contour integral calculations.

Modeling Fracture and Failure with Abaqus

L2.11

Modeling Sharp Cracks in Two Dimensions
? Other options for defining the crack front and crack tip
Crack front for a geometric instance

Crack tip for an orphan mesh

Crack front may be:
Vertex/Node Edges/Element edges Faces/Elements Geometric Instances Orphan Mesh Crack tip may be:

Vertex/Node

Geometric Instances

Orphan Mesh

Modeling Fracture and Failure with Abaqus

L2.12

Modeling Sharp Cracks in Two Dimensions
? Example: crack on a symmetry plane

? If the crack is on a symmetry plane, you do not need to define a seam.
? This feature can be used only for Mode I fracture.

Crack normal

Crack tip

Modeling Fracture and Failure with Abaqus

L2.13

Modeling Sharp Cracks in Two Dimensions
? Modeling the crack-tip singularity with second-order quad elements

? To capture the singularity in an 8-node isoparametric element:
? Collapse one side (e.g., the side made up by nodes a, b, and c) so that all three nodes have the same geometric location at the crack tip. ? Move the midside nodes on the sides connected to the crack tip to the ? point nearest the crack tip.

Modeling Fracture and Failure with Abaqus

L2.14

Modeling Sharp Cracks in Two Dimensions
? If nodes a, b, and c are free to move independently, then

A B ?? ? as r ? 0 r r
everywhere in the collapsed element.
? If nodes a, b, and c are constrained to move together, A = 0: ? The strains and stresses are square-root singular (suitable for linear elasticity). ? If nodes a, b, and c are free to move independently and the midside nodes remain at the midsides, B = 0 : ? The singularity in strain is correct for the perfectly plastic case.

? For materials in between linear elastic and perfectly plastic (most metals), it is better to have a stronger singularity than necessary.
? The numerics will force the coefficient of this singularity to be small.

Modeling Fracture and Failure with Abaqus

L2.15

Modeling Sharp Cracks in Two Dimensions
? Usage:

The crack tip nodes are independent: r -1 singularity

Quarter-point midside nodes on the sides connected to the crack tip The crack tip nodes are constrained: r -? singularity

3

4
1,2,3,4

1, 2

2 3

1

1,1,2,3
Modeling Fracture and Failure with Abaqus

L2.16

Modeling Sharp Cracks in Two Dimensions
? Aside: Controlling the position of midside nodes for orphan meshes ? Singularity controls cannot be applied to orphan meshes. ? Use the Mesh Edit tools to adjust their position.

Modeling Fracture and Failure with Abaqus

L2.17

Modeling Sharp Cracks in Two Dimensions
? If the side of the element is not collapsed but the midside nodes on the sides of the element connected to the crack tip are moved to the ? point: ? The strain is square root singular along the element edges but not in the interior of the element. ? This is better than no singularity but not as good as the collapsed element.

nodes moved to ? points
Modeling Fracture and Failure with Abaqus

L2.18

Modeling Sharp Cracks in Two Dimensions
? Angular resolution

? We need enough elements to resolve the angular dependence of the strain field around the crack tip.
? Reasonable results are obtained for LEFM if typical elements around the crack tip subtend angles in the range of 10? (accurate) to 22.5? (moderately accurate).

? Nonlinear material response usually requires finer meshes.

Modeling Fracture and Failure with Abaqus

L2.19

Modeling Sharp Cracks in Two Dimensions
? Modeling the crack-tip singularity with first-order quad elements

? Collapsing the side of a first-order quadrilateral element with independent nodes on the collapsed side gives

A ? ? as r ? 0. r

Modeling Fracture and Failure with Abaqus

L2.20

Modeling Sharp Cracks in Two Dimensions
? Example: Slanted crack in a plate

? To enable the creation of degenerate quads, you must create swept meshable regions around the crack tips (using partitions) and specify a quad-dominated mesh.
24 elements around crack tip: 15? angles

Quarterpoint nodes

CPE8R elements; typical nodal connectivity shows repeated node at crack tip: Quad-dominated mesh + swept technique for the circular regions surrounding the crack tips

Quadratic element type assigned to part

8, 8, 583, 588, 8, 1969, 1799, 1970

All crack-tip elements repeat node 8 in this example (nodes are constrained).

Modeling Fracture and Failure with Abaqus

L2.21

Modeling Sharp Cracks in Two Dimensions
? Example (cont’d): Alternate meshes ? No degeneracy:

? Degenerate with duplicate nodes:

With swept meshable region: CPE6M elements at crack tip — cannot be used for fracture studies in Abaqus.

With arbitrary mesh, singularity only along edges connected to crack tip.

CPE8R elements at crack tip but no repeated nodes:
1993, 1992, 583, 588, 2016, ...

Coincident nodes located at crack tip

Modeling Fracture and Failure with Abaqus

L2.22

Modeling Sharp Cracks in Two Dimensions
? Example (cont’d): Deformed shape

Focused mesh; deformation scale factor = 100

Arbitrary mesh; deformation scale factor = 100

Modeling Fracture and Failure with Abaqus

Modeling Sharp Cracks in Three Dimensions

L2.24

Modeling Sharp Cracks in Three Dimensions
? In three dimensions…

? The seam crack is modeled as a face partition that is either partially or totally embedded into a solid body.
? This can be done by partitioning or using a cut (Boolean) operation.

? The face along the seam will have duplicate nodes such that the elements on the opposite sides of the face will not share nodes.
? Wedge elements must be created along the crack front. ? Generally, this will require partitioning.

Penny-shaped seam crack: Full model

Quarter model

Wedge elements

Meshed model

Modeling Fracture and Failure with Abaqus

L2.25

Modeling Sharp Cracks in Three Dimensions
? Options for defining the crack front and crack line
Crack front for a geometric instance
Crack line for an orphan mesh

Crack front may be:
Edges/Element edges Faces/Element faces Cells/Elements Geometric Instances Orphan Mesh Geometric Instances Orphan Mesh

Crack line may be:
Edges/Element edges

Modeling Fracture and Failure with Abaqus

L2.26

Modeling Sharp Cracks in Three Dimensions
? Specifying the crack growth direction in three dimensions

? In 3D you can specify either the
? normal to the crack plane (only when the crack is planar)

or the
? virtual crack extension direction (the q vector).

? Only a single q vector can be defined for geometric instances.
? The implications of this will be discussed shortly.

Modeling Fracture and Failure with Abaqus

L2.27

Modeling Sharp Cracks in Three Dimensions
? Modeling the crack-tip singularity in three dimensions

? 20-node and 27-node bricks can be used with a collapsed face to create singular fields.
C3D20(RH)
midplane edge plane

2 nodes collapsed to the same location

crack line midside nodes moved to ? points

3 nodes collapsed to the same location

Modeling Fracture and Failure with Abaqus

L2.28

Modeling Sharp Cracks in Three Dimensions
? On an edge plane (orthogonal to the crack line):
Double-edge notch specimen (symmetry model)

??

A as r ? 0 r

??

A B ? as r ? 0 r r

??

B as r ? 0 r

Crack line

Edge plane nodes displace independently

Edge plane nodes displace together

Modeling Fracture and Failure with Abaqus

L2.29

Modeling Sharp Cracks in Three Dimensions
? On a midplane for 20-node bricks:

? If the two nodes on the collapsed face at the midplane can displace independently, ? ? r -1 at the midplane (i.e., element interior).

? If on each plane there is only one node along the crack line, no singularity is represented within the element. ? In either case the interpolation is not the same on the midplane as on an edge plane. ? This generally causes local oscillations in the J-integral values along the crack line.
Modeling Fracture and Failure with Abaqus

L2.30

Modeling Sharp Cracks in Three Dimensions
? On a midplane for 27-node bricks with all the extra nodes on the element faces:
midplane

C3D27(RH)

edge plane

3 nodes collapsed to same location

centroid

crack line
3 nodes collapsed to same location

Modeling Fracture and Failure with Abaqus

L2.31

Modeling Sharp Cracks in Three Dimensions
? If all midface nodes and the centroid node are included and moved with the midside nodes to the ? points, the singularity can be made the same on the edge planes and midplane. ? Abaqus does not allow the centroid node to be moved from the geometric centroid of the element. ? Therefore, the behavior at the midplane will never be the same as at the edge planes. ? This usually causes some small oscillation of the crack fields along the crack line. ? The midface node marked “A” is frequently omitted.

? This creates differences in interpolation between the midplane and the edge planes and, hence, causes further oscillation in the cracktip fields.
? These oscillations are minor in most cases.

Modeling Fracture and Failure with Abaqus

L2.32

Modeling Sharp Cracks in Three Dimensions
? Example: Conical crack in a halfspace ? A conical crack in an infinite halfspace is considered. ? Only the aspects related to the geometric modeling are considered here. ? The results of this analysis (J-integral values, etc) will be considered in the next lecture. ? The modeling procedure is outlined next.

Modeling Fracture and Failure with Abaqus

L2.33

Modeling Sharp Cracks in Three Dimensions
1 Example (cont’d): Create the basic geometry

? Because of symmetry, only a quarter model is created

a = 15
r = 10

q

= 45?

Large solid block (300 × 300 × 300) used to represent the half-space.

Conical shell of revolution (revolved 90? ); this will be used to cut the block.

Modeling Fracture and Failure with Abaqus

L2.34

Modeling Sharp Cracks in Three Dimensions
2 Example (cont’d): Merge the block and cone

? This will create the edges and surface necessary to define the seam and the crack.
Instance and merge the two parts to create a new part. The instance must be independent.

Modeling Fracture and Failure with Abaqus

L2.35

Modeling Sharp Cracks in Three Dimensions
3 Example (cont’d): Define the seam and the crack front/line

Only one q vector can be defined for geometry. The q vectors will be adjusted at the end of the modeling process by editing an orphan mesh.
Modeling Fracture and Failure with Abaqus

L2.36

Modeling Sharp Cracks in Three Dimensions
4 Example (cont’d): Partition the block for meshing
The regions surrounding the crack front are partitioned to permit structured meshing.

A small curved tube is centered at the crack tip; this region is meshed with a single layer of wedge elements. This mesh is swept along the length of the tube.

Modeling Fracture and Failure with Abaqus

L2.37

Modeling Sharp Cracks in Three Dimensions
? Aside: Why is the small curved tube needed?
The swept meshing technique sweeps a mesh through a cross section. For the curved tube, this implies the sweep direction is along its length. In order for Abaqus to automatically create a focused mesh at the crack tip, however, it would need to sweep around the circumference. To overcome this, two concentric tubes are used; the smaller one is meshed with a single layer of wedge elements (which is then swept along the length of the tube). If only a single curved tube was created (shown at right), the mesh around the crack tip would be arbitrary—not focused (wedge elements not created).

Modeling Fracture and Failure with Abaqus

L2.38

Modeling Sharp Cracks in Three Dimensions
? Aside: What about the seam?

? After all the partitions are created for meshing purposes, the definition of the seam remains intact.

Mesh seam

Modeling Fracture and Failure with Abaqus

L2.39

Modeling Sharp Cracks in Three Dimensions
5 Example (cont’d): Mesh the part

? Specify appropriate edge seeds to create a focused mesh around the crack front with minimal mesh distortion.

Modeling Fracture and Failure with Abaqus

L2.40

Modeling Sharp Cracks in Three Dimensions
6 Example (cont’d): Adjust the

q vectors

? As noted earlier, only a single q vector can be defined for geometry. As seen in the figure, the vector that was defined is only accurate at the left end of the crack line. ? Individual q vectors can be defined on an orphan mesh, however. Thus, either…
? Create a mesh part (Mesh module)
To take advantage of the input file approach, define a set that contains the conical region before writing the input file. Then you will be able to easily create a display group based on this set when manipulating the orphan mesh.

or
? Write an input file and import the model ? This approach has the advantage that it preserves attributes (sets, loads, etc).

Modeling Fracture and Failure with Abaqus

L2.41

Modeling Sharp Cracks in Three Dimensions
? For the orphan mesh, adjust each vector individually

To redefine this particular vector, select these nodes as the start and end points of the vector.

Modeling Fracture and Failure with Abaqus

L2.42

Modeling Sharp Cracks in Three Dimensions
? For all elements, the singularities are modeled best if the element edges are straight. ? In three dimensions the planes of the element perpendicular to the crack line should be flat. ? If they are not, when the midside nodes are moved to the ? points, the Jacobian of the element at some integration points may be negative. ? One way to correct this is to move the midside nodes slightly away from the ? points toward the midpoint.

Modeling Fracture and Failure with Abaqus

L2.43

Modeling Sharp Cracks in Three Dimensions
? Example: Conical crack model

Modeling Fracture and Failure with Abaqus

Finite-Strain Analysis of Crack Tips

L2.45

Finite-Strain Analysis of Crack Tips
? Finite-strain analyses:

? Singular elements should not be used (normally).
? The mesh must be sufficiently refined to model the very high strain gradients around the crack tip if details in this region are required. ? Even if only the J-integral is required, the deformation around the crack tip may dominate the solution and the crack-tip region will have to be modeled with sufficient detail to avoid numerical problems.

Modeling Fracture and Failure with Abaqus

L2.46

Finite-Strain Analysis of Crack Tips
? Physically, the crack tip is not perfectly sharp, and such modeling makes it difficult to obtain results. ? Instead, we model the tip as a blunted notch, with a suggested radius ? 10-3rp. ? Here, rp is the size of the plastic zone (discussed in Lecture 1).

? The notch must be small enough that under the applied loads, the deformed shape of the notch no longer depends on the original geometry.
? Typically, the notch must blunt out to more than four times its original radius for this to be true.

Modeling Fracture and Failure with Abaqus

L2.47

Finite-Strain Analysis of Crack Tips
? Geometric modeling of blunt cracks

? In 2D, the geometry of a blunted (or open) crack is modeled as a cut having a significant thickness.
? Meshing is done in the usual way. ? A very fine mesh is required at the crack tip. ? This can be achieved by simply assigning small element sizes to the notch.

Modeling Fracture and Failure with Abaqus

L2.48

Finite-Strain Analysis of Crack Tips
? 3D open cracks can be created in Abaqus/CAE in one of two ways: ? Adding a Cut feature in the Part module. ? Subtracting a flaw from the original part with a Boolean operation in the Assembly module. ? Hex meshing more difficult due to irregular geometry. ? Creating a fine mesh at the crack front generally requires many partitions.
Quarter model Penny shaped open crack: Full model

Meshed model

Partitions to control mesh
Modeling Fracture and Failure with Abaqus

Refined mesh

L2.49

Finite-Strain Analysis of Crack Tips
? The size of the elements around the notch must be about 1/10 th the notch-tip radius. Biased edge seeds can
reduce the size of the mesh by focusing small elements towards the crack tip. SEN specimen

crack-tip mesh

rnotch
10% of rnotch
Modeling Fracture and Failure with Abaqus

L2.50

Finite-Strain Analysis of Crack Tips
? For J-integral evaluation, the region on the surface of the blunted notch should be used to define the crack front.

Crack tip region

q vector
Crack surface is detected automatically

The blunted notch surface is the crack front region

Symmetry plane

? For the J- and Ct-integrals to be path independent, the crack surfaces must be parallel to one another (or parallel to the symmetry plane). ? If this is not the case, Abaqus automatically generates normals on the crack surface. ? If the notch radius shrinks to zero, all nodes that would be at the crack tip should be included in the crack-tip node set.
Modeling Fracture and Failure with Abaqus

L2.51

Finite-Strain Analysis of Crack Tips
? If the mesh is so coarse that the integration points nearest the crack tip are far from the tip, most of the details (accurate stresses and strains) of the finite-strain region around the crack tip will be lost. ? However, accurate J values may still be obtained if cracks are modeled as sharp.

Modeling Fracture and Failure with Abaqus

L2.52

Finite-Strain Analysis of Crack Tips
? Example: SEN specimen

Deformed shape

Moderate blunting
Undeformed shape

Severe blunting

Deformed vs Undeformed Shapes
Modeling Fracture and Failure with Abaqus

Contours of PEEQ

L2.53

Finite-Strain Analysis of Crack Tips
? In situations involving finite rotations but small strains, such as the bending of slender structures, a small keyhole around the crack tip should be modeled.

crack-front region

? The region defining the crack front for the contour integral consists of the region on the keyhole.

? The elements should not be singular.
Modeling Fracture and Failure with Abaqus

Limitations Of 3D Swept Meshing For Fracture

L2.55

Limitations Of 3D Swept Meshing For Fracture
? For curved regions cannot generate wedges at the center using a hexdominated approach and then sweep along the length of the region. ? This was discussed earlier in the context of the conical crack problem. ? To create a focused mesh in this case, embed a small tube within a larger concentric tube. Mesh the smaller tube with a single layer of wedge elements; the surrounding regions are meshed with hex elements.
Sweep direction

Modeling Fracture and Failure with Abaqus

L2.56

Limitations Of 3D Swept Meshing For Fracture
? Partition for a penny-shaped crack

? Illustrates the limitation that the path for the partition must be perpendicular to its bounding surfaces; thus, cannot properly partition along the arc of a circle as shown in this example:

Tangent direction of arc
arc (not a semi-circle as in previous example)

Cross-sectional view of block

Partition by sweeping circular edge along arc

Modeling Fracture and Failure with Abaqus

L2.57

Limitations Of 3D Swept Meshing For Fracture
? The workaround is to partition the face with circular arcs, and then partition the cell using the n-sided patch technique.
n-sided patch

Face partition

Note that the cross-sectional area of the swept region is not constant along its length because the tangents at the ends are not perpendicular to the block (generalized sweep meshing)

Resulting mesh around the crack front using wedge elements

Modeling Fracture and Failure with Abaqus

Modeling Cracks with Keyword Options

L2.59

Modeling Cracks with Keyword Options
? Defining a crack with keyword options:

? The *CONTOUR INTEGRAL option is used to define both, the crack itself and the fracture output, in an Abaqus input ( .inp) file.
? In this section, we focus solely on the crack-specific parameters of this option. ? These include:
*CONTOUR INTEGRAL, SYMM, NORMAL

? In the next lecture, we discuss the output-specific parameters of this option. ? As noted earlier, the main requirements in defining a crack are: ? Defining the crack front

? Defining the crack extension direction

Modeling Fracture and Failure with Abaqus

L2.60

Modeling Cracks with Keyword Options
? Crack symmetry
*CONTOUR INTEGRAL, SYMM

? The crack lies on a plane of symmetry and only half the structure is being modeled ? This feature should only be used for Mode I problems.

Modeling Fracture and Failure with Abaqus

L2.61

Modeling Cracks with Keyword Options
? Crack extension
*CONTOUR INTEGRAL, NORMAL

? The NORMAL parameter is used to define the normal to the crack plane when the crack is planar. ? Usage:
*contour integral, normal nx, ny, nz nodeSet1, nodeSet2, ...

? In this case, give a list of the node set names defining the crack front from one end to the other end, in sequential order, without missing any points on the crack line. ? In two-dimensional cases, only one node set is needed.

These sets define the crack front; the first node in each set defines the crack tip node for that set.

(An optional CRACK TIP NODES parameter is available to specify the crack tip nodes directly).

Modeling Fracture and Failure with Abaqus

L2.62

Modeling Cracks with Keyword Options
? Example: Penny-shaped crack in an infinite space
*Contour integral, symm, normal, ... 0.0, 1.0, 0.0 Crack-Front-1, Crack-Front-2, Crack-Front-3, ...

Crack-Front-1

Modeling Fracture and Failure with Abaqus

L2.63

Modeling Cracks with Keyword Options
? If the NORMAL parameter is omitted, we must give the crack-tip node set name, and the crack propagation direction q, at each node set defining the crack front. ? Usage:
*contour integral, ... nodeSet1, (qx)1, (qy)1, (qz)1 nodeSet2, (qx)2, (qy)2, (qz)2 :

? Data must start with the node set at one end and be given for each node set defining the crack line sequentially until the other end of the crack is reached.
? The first node in each set is the crack tip node for that set unless the CRACK TIP NODES parameter is used.

? This format allows nonplanar cracks to be analyzed.

Modeling Fracture and Failure with Abaqus

L2.64

Modeling Cracks with Keyword Options
? Example: conical crack in an infinite half-space
*Contour integral, ... Crack-Front-1, 0.707107, -0.707107, 0. Crack-Front-2, 0.705994, -0.707107, 0.0396478 Crack-Front-3, 0.702661, -0.707107, 0.0791708

Crack-Front-1

Modeling Fracture and Failure with Abaqus

L2.65

Modeling Cracks with Keyword Options
? Generating a focused mesh with keyword options

? Example: DEN specimen
? The focused mesh shown in the figure will be generated with the use of keyword options. ? The options include

*NODE
*NGEN

*NFILL
*ELEMENT

*ELGEN

Modeling Fracture and Failure with Abaqus

L2.66

Modeling Cracks with Keyword Options
? Node definitions
*node 1, 16001, 101, 4101, 12101, 16101, *ngen, 1, *ngen, 101, 4101, 12101, 0.0125, 0.0000 0.0125, 0.0000 0.0250, 0.0000 0.0250, 0.0125 0.0000, 0.0125 0.0000, 0.0000 nset=tip 16001, 1000 nset=outer 4101, 1000 12101, 1000 16101, 1000
12101
8101 4101

14101

2101

16101

101

tip *NGEN generates nodes incrementally between any two previously defined nodes. In this example, 17 crack-tip nodes are created (contained in the set tip); the 17 nodes on the outer boundary are contained in set outer.

Start node

End node

Increment in node number

Modeling Fracture and Failure with Abaqus

L2.67

Modeling Cracks with Keyword Options
? Quarter-point nodes
*nfill, singular=1 tip, outer, 10, 10
This parameter generates quarterpoint nodes; the 1 indicates the first bound represents the crack tip

Start set: first bound

End set: second bound

Node number increment

8021 4021 2021 1021 21

Number of intervals between bounding nodes

4011 11

31

*NFILL generate nodes for a region of a mesh by filling in nodes between two bounds. In this example, 10 rows of nodes are generated between each tip node and its corresponding outer node.
Modeling Fracture and Failure with Abaqus

L2.68

Modeling Cracks with Keyword Options
? Element definitions
*element, type=cps8r 1, 1, 21, 2021, 2001, 11, 1021, 2011, 1001 *elgen, elset=plate 1, 5, 20, 10, 8, 2000, 1000
First row of elements Total number of rows Nodes 1, 1001, and 2001 are coincident

1

2021 1021

*ELGEN generates elements incrementally.

21

11
1

In this example, 5 elements form the first row (extending radially outward from the tip); a total of 8 rows of elements (based on the first row) are created around the crack tip.

Modeling Fracture and Failure with Abaqus

L2.69

Modeling Cracks with Keyword Options
? Crack-tip nodes

? If the crack-tip nodes are permitted to behave independently, the strength of the strain-field singularity is ? ? r -1. ? The crack-tip nodes can be constrained using equations, multi-point constraints, using repeated nodes in the element definition, etc. For example, to constrain the crack-tip nodes with a multi-point constraint:
*nset, nset=constrain, generate 1, 15001, 1000 *mpc tie, constrain, 16001

? Only node 16001 is independent in this case.
? The strain-field singularity is ? ? r -?.

Modeling Fracture and Failure with Abaqus

Fracture Analysis
Lecture 3

L3.2

Overview
? Calculation of Contour Integrals ? Examples ? Nodal Normals in Contour Integral Calculations

? J-Integrals at Multiple Crack Tips
? Through Cracks in Shells
? Mixed-Mode Fracture

? Material Discontinuities
? Numerical Calculations with Elastic-Plastic Materials

? Workshop 1
? Workshop 2

Modeling Fracture and Failure with Abaqus

Calculation of Contour Integrals

L3.4

Calculation of Contour Integrals
? Abaqus offers the evaluation of J-integral values, as well as several other parameters for fracture mechanics studies. These include: ? The KI, KII, and KIII stress intensity factors, which are used mainly in linear elastic fracture mechanics to measure the strength of local crack tip fields; ? The T-stress in linear elastic calculations; ? The crack propagation direction: an angle at which a preexisting crack will propagate; and ? The Ct-integral, which is used with time-dependent creep behavior.

? Output can be written to the output database ( .odb), data (.dat), and results (.fil) files.

Modeling Fracture and Failure with Abaqus

L3.5

Calculation of Contour Integrals
? Domain representation of J

? For reasons of accuracy, J is evaluated using a domain integral. ? The domain integral is evaluated over an area/volume contained within a contour surrounding the crack tip/line.
? In two dimensions, Abaqus defines the domain in terms of rings of elements surrounding the crack tip. ? In three dimensions, Abaqus defines a tubular surface around the crack line.

Modeling Fracture and Failure with Abaqus

L3.6

Calculation of Contour Integrals
? Different contours (domains) are created automatically by Abaqus. ? The first contour consists of the crack front and one layer of elements surrounding it.

? Ring of elements from one crack surface to the other (or the symmetry plane).
? The next contour consists of the ring of elements in contact with the first contour as well as the elements in the first contour.

Contour 1

Contour 2

? Each subsequent contour is defined by adding the next ring of elements in contact with the previous contour.

Contour 3

Contour 4

Modeling Fracture and Failure with Abaqus

L3.7

Calculation of Contour Integrals
? The J-integral and the Ct-integral at steady-state creep should be path (domain) independent. ? The value for the first contour is generally ignored. ? Examples of contour domains:

2nd contour

1st contour crack-front nodes
Crack-tip node

2nd contour

1st contour

Crack-tip node

Modeling Fracture and Failure with Abaqus

L3.8

Calculation of Contour Integrals
? Usage:
*CONTOUR INTEGRAL, CONTOURS= n, TYPE={J, C, T STRESS, K FACTORS}, DIRECTION = {MTS, MERR, KII0}

Specifies the number of contours (domains) on which the contour integral will be calculated

This is the output frequency in increments

Note: In this lecture, we focus on the output-specific parameters of the *CONTOUR INTEGRAL option. The crack-specific parameters SYMM and NORMAL were discussed in the previous lecture.

Modeling Fracture and Failure with Abaqus

L3.9

Calculation of Contour Integrals
? Usage (cont’d):
*CONTOUR INTEGRAL, CONTOURS= n, TYPE={J, C, T STRESS, K FACTORS}, DIRECTION = {MTS, MERR, KII0}

? J for J-integral output,
? C for Ct-integral output. ? T STRESS to output T-stress calculations ? K FACTORS for stress intensity factor output

Modeling Fracture and Failure with Abaqus

L3.10

Calculation of Contour Integrals
? Usage (cont’d):
*CONTOUR INTEGRAL, CONTOURS= n, TYPE={J, C, T STRESS, K FACTORS}, DIRECTION = {MTS, MERR, KII0}
Three criteria to calculate the crack propagation direction at initiation

? Use with TYPE=K FACTORS to specify the criterion to be used for estimating the crack propagation direction in homogenous, isotropic, linear elastic materials:
? Maximum tangential stress criterion (MTS) ? Maximum energy release rate criterion (MERR) ? KII = 0 criterion (KII0)
Modeling Fracture and Failure with Abaqus

L3.11

Calculation of Contour Integrals
? Output files
*CONTOUR INTEGRAL, OUTPUT

? Set OUTPUT=FILE to store the contour integral values in the results (.fil) file.
? Set OUTPUT=BOTH to print the values in the data and results files. ? If the parameter is omitted, the contour integral values will be printed in the data (.dat) file but not stored in the results (.fil) file.

Modeling Fracture and Failure with Abaqus

L3.12

Calculation of Contour Integrals
? Loads

? Loads included in contour integral calculations:
? Thermal loads.

? Crack-face pressure and traction loads on continuum elements as well as those applied using user subroutines DLOAD and UTRACLOAD. ? Surface traction and crack-face edge loads on shell elements as well as those applied using user subroutine UTRACLOAD. ? Uniform and nonuniform body forces.
? Centrifugal loads on continuum and shell elements. ? Not all types of distributed loads (e.g., hydrostatic pressure and gravity loads) are included in the contour integral calculations. ? The presence of these loads will result in a warning message.

Modeling Fracture and Failure with Abaqus

L3.13

Calculation of Contour Integrals
? Other loads not included in contour integral calculations:

? Contributions due to concentrated loads are not included.
? If needed, modify the mesh to include a small element and apply a distributed load to the element. ? Contributions due to contact forces are not included.

? Initial stresses are not considered in the definition of contour integrals.

Modeling Fracture and Failure with Abaqus

Examples

L3.15

Examples
? Penny-shaped crack in an infinite space

? Model characteristics
? The mesh is extended far enough from the crack tip so that the finite boundaries will not influence the crack-tip solution. ? The radius of the penny-shaped crack is 1.

? Two types of loading are considered: ? Uniform far-field loading
? Nonuniform loading on the crack face: p = Ar n.

Modeling Fracture and Failure with Abaqus

L3.16

Examples
? Different mesh characteristics:
20

? Axisymmetric or three-dimensional
? Fine or coarse focused meshes

? With or without ? point elements
? Various element types used:
20

? First- and second-order
? With and without reduced integration

Axisymmetric model

Crack tip
Focused mesh around crack tip
Modeling Fracture and Failure with Abaqus

L3.17

Examples
? Fine mesh vs. coarse mesh (axisymmetric and 3D models)

0.0004

0.08

The fine mesh is shown to the left; the coarse mesh above. The length perpendicular to crack line of the crack-tip elements are indicated.

~0.08
Modeling Fracture and Failure with Abaqus

L3.18

Examples
? Axisymmetric model: geometry

Symmetry planes

Close up of crack tip region for coarse mesh model (identical for fine mesh model—only the inner semicircular region is smaller)

Model geometry

Modeling Fracture and Failure with Abaqus

L3.19

Examples
? Axisymmetric model: crack definition
Crack tip with extension direction

Set to 0.5 to use midpoint rather than ? point elements

Modeling Fracture and Failure with Abaqus

L3.20

Examples
? 3D model: geometry and mesh
Fine 3D mesh

? A 90? sector is modeled because of symmetry.
Additional partition required for swept mesh

Symmetry planes

On planes perpendicular to the crack front, the mesh is very similar to the axisymmetric mesh Partitions used for coarse mesh model (identical for fine mesh model—only the inner semicircular region is smaller) In the circumferential direction around the crack line, 12 elements are used.

Modeling Fracture and Failure with Abaqus

L3.21

Examples
? Why is the additional partition required?

? Without the additional partition, the region shown below would require irregular elements at the vertex located on the axis of symmetry.
? This is not supported by Abaqus.
Irregular elements required here because revolving about a point

A 7-node element is an example of an irregular element.

Modeling Fracture and Failure with Abaqus

L3.22

Examples
? 3D model: crack definition

? Orphan mesh created to edit q vectors.

Modeling Fracture and Failure with Abaqus

L3.23

Examples
? Contour integral output requests (axisymmetric and 3D)

Separate output requests are required for J, K-factors, and the T-stress.

Modeling Fracture and Failure with Abaqus

L3.24

Examples
? Loads (axisymmetric and 3D)
The far-field load is suppressed.

Modeling Fracture and Failure with Abaqus

L3.25

Examples
? Results

? MISES stress shown below for the axisymmetric fine mesh.

J analytical ? J numerical J analytical

?100%

Deformation scale factor = 250

Analytical

Contour 1

Contour 2

Contour 3

Contour 4

Contour 5

5.796E-02

5.8169E-02
Contour 6 5.8064E-02

5.8095E-02
Contour 7 5.8044E-02

5.8121E-02
Contour 8 5.8024E-02

5.8104E-02
Contour 9 5.8005E-02

5.8084E-02
Contour 10 5.7985E-02

Modeling Fracture and Failure with Abaqus

L3.26

Examples J values from meshes with ? point elements (reduced integration)
3-D Loading Uniform far field Uniform crack face Nonuniform crack face (n = 1) Nonuniform crack face (n = 2) Analytical result C3D20R Coarse Fine Axisymmetric CAX8R Coarse Fine

.0580
.0580 .0358 .0258 .0201

.0578
.0578 .0356 .0256 .0199

.0580
.0580 .0357 .0260 .0206

.0579
.0579 .0356 .0256 .0200

.0581
.0581 .0358 .0258 .0202

Nonuniform crack face (n = 3)

? Abaqus values are based on the average of contours 3 ?5 in each mesh.
Modeling Fracture and Failure with Abaqus

L3.27

Examples J values from meshes with ? point elements (full integration)
3-D Loading Uniform far field Uniform crack face Nonuniform crack face (n = 1) Nonuniform crack face (n = 2) Analytical result C3D20 Coarse Fine Axisymmetric CAX8 Coarse Fine

.0580
.0580 .0358 .0258 .0201

.0577
.0577 .0355 .0255 .0198

.0572
.0572 .0352 .0253 .0197

.0578
.0578 .0356 .0255 .0199

.0580
.0580 .0358 .0258 .0201

Nonuniform crack face (n = 3)

? Abaqus values are based on the average of contours 3 ?5 in each mesh.
Modeling Fracture and Failure with Abaqus

L3.28

Examples J values from meshes without ? point elements (reduced integration)
3-D Axisymmetric

Loading Uniform far field
Uniform crack face Nonuniform crack face (n = 1) Nonuniform crack face (n = 2) Nonuniform crack face (n = 3)

Analytical result

C3D20R
Coarse Fine

C3D8R
Coarse

CAX8R
Coarse Fine

CAX4R
Coarse

.0580

.0574

.0580

.0563

.0574

.0581

.0562

.0580 .0358
.0258 .0201

.0574 .0350
.0250 .0193

.0580 .0357
.0260 .0206

.0563 .0336
.0234 .0177

.0574 .0350
.0250 .0193

.0581 .0358
.0258 .0202

.0562 .0337
.0236 .0179

? Abaqus values are based on the average of contours 3 ?5 in each mesh.
Modeling Fracture and Failure with Abaqus

L3.29

Examples J values from meshes without ? point elements (full integration)
3-D Axisymmetric

Loading Uniform far field
Uniform crack face Nonuniform crack face (n = 1) Nonuniform crack face (n = 2) Nonuniform crack face (n = 3)

Analytical result

C3D20
Coarse Fine

C3D8
Coarse

CAX8
Coarse Fine

CAX4
Coarse

.0580

.0573

.0572

.0552

.0574

.0580

.0557

.0580 .0358
.0258 .0201

.0573 .0350
.0249 .0193

.0572 .0352
.0253 .0197

.0552 .0329
.0229 .0172

.0574 .0350
.0250 .0193

.0580 .0358
.0258 .0201

.0557 .0333
.0232 .0175

? Abaqus values are based on the average of contours 3 ?5 in each mesh.
Modeling Fracture and Failure with Abaqus

L3.30

Examples
? Conclusions

? 3D fine meshes with second-order elements are more sensitive to the choice of integration rule when determining J. ? The results are still very accurate (within 2% of analytical value).
? The inclusion of the singularity helps most in the coarser meshes.

? For mesh convergence in small strain, the singularity must be included.

Modeling Fracture and Failure with Abaqus

L3.31

Examples
? Conical crack in a half-space

? At each node set along the crack front, the crack propagation direction is different.

Modeling Fracture and Failure with Abaqus

L3.32

Examples
? Three-dimensional model

? Displaced shape and Mises stress distribution of full threedimensional model.

Deformation scale factor = 1.e6

Modeling Fracture and Failure with Abaqus

L3.33

Examples ? J values of three-dimensional mesh
? There is some oscillation between J values evaluated at corner nodes compared to J values evaluated at midside nodes.

Variation of J with angular position
1.338E-07
J-integral

1.336E-07
1.334E-07

3D contour 5 3D contour 4 3D contour 3 3D contour 2 0 45 Angle (degrees) 90

1.332E-07
1.330E-07 1.328E-07

Modeling Fracture and Failure with Abaqus

L3.34

Examples
? Axisymmetric model and results

Contours 3-5 have converged

Axisymmetric results are used as reference results.

Modeling Fracture and Failure with Abaqus

L3.35

Examples
? Comparison of axisymmetric and 3D results
Variation of J with angular position Contour 1
1.380E-07 J -integral 1.360E-07 1.340E-07 J -integral 3D AXI

Variation of J with angular position Contour 2
1.334E-07 1.333E-07 1.332E-07 1.331E-07 1.330E-07 1.329E-07 0 45 Angle (degrees) 90

3D AXI

1.320E-07
1.300E-07 0 45 Angle (degrees) 90

Variation of J with angular position Contour 3
1.336E-07 J -integral 1.334E-07 1.332E-07 1.330E-07 3D AXI J -integral

Variation of J with angular position Contour 5
1.338E-07 1.336E-07 1.334E-07 1.332E-07 1.330E-07 1.328E-07

3D AXI

1.328E-07 0 45
Angle (degrees)

90

0

45
Angle (degrees)

90

Modeling Fracture and Failure with Abaqus

L3.36

Examples
? Since the three-dimensional mesh is quite coarse around the axis of symmetry, these results are considered to be good —the error is less than 0.5% for all but the first contour.

% difference in J between AXI and 3D results
3.5 3.0 2.5 2.0 1.5 1.0 0.5 0.0 0 45 Angle (degrees) 90 Contour 1

% difference

Contour 2
Contour 3

Contour 4
Contour 5

Modeling Fracture and Failure with Abaqus

L3.37

Examples
? Submodeling

? We can use submodeling to create two meshes that are significantly smaller than the full threedimensional model.
? The top-right figure is the coarse mesh global model in the vicinity of the crack.

? The bottom-right figure shows the refined submodel mesh overlaid on the global model mesh.

Modeling Fracture and Failure with Abaqus

L3.38

Examples ? J values of submodel:
% difference

% difference in J between AXI and 3D results
4.5 4.0 3.5 3.0 2.5 2.0 1.5 1.0 0.5 0.0 0 45 Angle (degrees) 90

? Inaccuracies are introduced by the coarser mesh used in the global model. ? Errors in J are less than 1%. ? CPU time was reduced by a factor of 3.
Variation of J with angular position
1.326E-07 J -integral 1.324E-07 1.322E-07 1.320E-07 1.318E-07 0 45 90 3D contour 5

Contour 1 Contour 2 Contour 3 Contour 4 Contour 5

Variation of J with angular position Contour 5
1.335E-07 J -integral 1.330E-07 1.325E-07 1.320E-07 1.315E-07 0 45 Angle (degrees) 90 3D

3D contour 4 3D contour 3 3D contour 2

AXI

Angle (degrees)

Modeling Fracture and Failure with Abaqus

L3.39

Examples
? Compact Tension Specimen

? This is one of five standardized specimens defined by the ASTM for the characterization of fracture initiation and crack growth.
? The ASTM standardized testing apparatus uses a clevis and a pin to hold the specimen and apply a controlled displacement.

Modeling Fracture and Failure with Abaqus

L3.40

Examples
? Model details ? Plane strain conditions assumed. ? The initial crack length is 5 mm. ? Elastic-plastic material

Prescribed load line displacement

Crack seam

? Low alloy ferritic steel

q-vector

1/√r singularity modeled in the crack-tip elements

Modeling Fracture and Failure with Abaqus

L3.41

Examples
? Results

Small strain analysis

Finite strain analysis

Modeling Fracture and Failure with Abaqus

L3.42

Examples

At small to moderate strain levels, the small and finite strain models yield similar results.

Finite strain effects must be considered to represent this level of deformation and strain accurately.

Modeling Fracture and Failure with Abaqus

Nodal Normals in Contour Integral Calculations

L3.44

Nodal Normals in Contour Integral Calculations
? Sharp curved cracks

? For sharp cracks, if the crack faces are curved, Abaqus automatically determines the normal directions of the nodes on the portions of the crack faces that lie within the contour integral domains.
? This improves the accuracy of the contour integral estimation. ? The normal is not used at the crack-tip node, however.

Normals to top crack surface nodes

n (normal to
crack plane) Normals to bottom crack surface nodes

q

Modeling Fracture and Failure with Abaqus

L3.45

Nodal Normals in Contour Integral Calculations
? Example: sharp curved crack

Contour #
J without normals J with normals

1
3.363 3.600

2
2.980 3.602

3
2.475 3.605

4
1.888 3.605

5
1.283 3.605

Modeling Fracture and Failure with Abaqus

L3.46

Nodal Normals in Contour Integral Calculations
? Blunt cracks and notches

? All nodes on the notch should be included in the crack-tip node set.
? The J-integral results are more accurate since the q vector is parallel to the crack surface in this case, as illustrated below.
Paths for contour integrals

Crack surface

Crack surface

n q
Single node in crack-tip node set; normals calculated on nodes of blunted surface; q not parallel to crack surface. All nodes on blunted surface in crack-tip node set; q parallel to crack surface.

q

Modeling Fracture and Failure with Abaqus

J-Integrals at Multiple Crack Tips

L3.48

J-Integrals at Multiple Crack Tips
? Abaqus can calculate J (or Ct ) at multiple crack tips ? Abaqus/CAE: multiple crack tips and history output requests ? Input file: repeated use of the *CONTOUR INTEGRAL option.

? If the domain for one crack tip envelopes the other crack tip, the J value will go to zero (as it should).

Modeling Fracture and Failure with Abaqus

Through Cracks in Shells

L3.50

Through Cracks in Shells
? Second-order quadrilateral shell elements must be used if contour integral output is requested. ? Sides of S8R elements should not be collapsed. If a focused mesh is used, the crack tip must be modeled as a keyhole whose radius is small compared to the other dimensions measured in the plane of the shell.

Shell mesh

Crack-tip mesh for S8R elements

Modeling Fracture and Failure with Abaqus

L3.51

Through Cracks in Shells
? S8R5 elements can be collapsed and midside nodes moved to the 1/4 points.

Shell mesh

Crack-tip mesh for S8R5 elements

? The q vector must lie in the shell surface.

? It should be tangent to the surface.

Modeling Fracture and Failure with Abaqus

L3.52

Through Cracks in Shells
? Example: Circumferential through crack under axial load

? Mean radius R = 10.5 in ? Wall thickness t = 0.525 in ? Crack half-angle q = p / 4

? Longitudinal membrane stress = 100 psi

Modeling Fracture and Failure with Abaqus

L3.53

Through Cracks in Shells
? Model details ? Axial load is applied using a shell edge load ? Symmetry used to reduce mode size

Edge loads

symmetry

Modeling Fracture and Failure with Abaqus

L3.54

Through Cracks in Shells
? Modeling a crack with a keyhole

Crack front

q vector

Crack tip

Modeling Fracture and Failure with Abaqus

L3.55

Through Cracks in Shells
? Results

Deformed shape—axial loading

J values—axial loading

Modeling Fracture and Failure with Abaqus

L3.56

Through Cracks in Shells
? In shell element meshes, mechanical loads which act normal to the shell surface and are applied within the contour integral domain are not taken into account in the calculation of the contour integral. ? For example, pressure loads are not considered because they act normal to the shell surface ? Conversely, axial edge loads are considered because they act in the shell surface. ? Two workarounds exist:

? Run successive shell models with differing crack lengths and numerically differentiate the potential energy ? Use solid elements (if the response is membrane dominated)

Modeling Fracture and Failure with Abaqus

L3.57

Through Cracks in Shells
? Using numerical differentiation to obtain J:

? ( PE ) J =? ?a Constant Load =? PE a ?Da ? PE a Da

Potential energy:

PE = ALLSE ? ALLWK
.

Constant Load

? The PE values should be obtained from two separate analyses, with crack lengths differing by Da. ? The values of PE in the Abaqus data (.dat) file are generally not printed to a sufficient number of figures to be useful for this calculation and must be read from the results ( .fil) file. ? A similar technique can be used to get Ct at long times.

Modeling Fracture and Failure with Abaqus

L3.58

Through Cracks in Shells
? Using solid elements:

? If membrane deformation is dominant, the shell can be modeled with a single layer of 20-node bricks since these solid elements include loading contributions to contour integrals.

Modeling Fracture and Failure with Abaqus

L3.59

Through Cracks in Shells
? To obtain accurate values of J through the shell thickness with solid elements, more than one element should be used in the thickness direction.

J values will show significant path dependence unless
averaged. ? If only one element is used through the thickness, the values can be averaged by thinking of J as a force per unit length: ? The average is calculated as if the J values were equivalent nodal forces:

J

shell

=

J A ? 4J B ? JC . 6

A B C

Modeling Fracture and Failure with Abaqus

L3.60

Through Cracks in Shells
? Aside: Generating a solid element mesh from a shell mesh.

? A shell mesh can easily be converted to a solid one using the ―Offset Mesh‖ tool. ? Creates solid layers from a shell mesh.

Modeling Fracture and Failure with Abaqus

L3.61

Through Cracks in Shells
? Example: Circumferential through crack in an internally pressurized, closed-end pipe ? The same pipe discussed earlier, now subjected to 10 psi internal pressure + axial load (which simulates the closed end).

? Comparison of J values using one layer of C3D20R elements through the thickness :
CONTOUR
1 At Node A At Node B At Node C Averaged 2

J values ? 100
3 4 5

A
B C

2.0965 3.7396 5.0226 3.6796

2.1317 3.6992 5.0501 3.6631

2.1505 3.7004 5.0813 3.6722

2.1557 3.6968 5.1471 3.6817

2.1697 3.6904 5.2373 3.6948

Modeling Fracture and Failure with Abaqus

L3.62

Through Cracks in Shells
? Example: Circumferential through crack under axial load revisited

? Now we revisit the problem in which the pipe is subjected to an axial load.
? Comparison of J values using one layer of C3D20R elements through the thickness:
CONTOUR

J values ? 100
1 2 3 4 5

At Node A At Node B
At Node C Averaged Analytical

2.2122 3.7629
4.9560 3.7033 3.7181

2.2524 3.7202
4.9893 3.6871

2.2700 3.7212
5.0175 3.6954

2.2740 3.7184
5.0737 3.7036

2.2850 3.7136
5.1492 3.7148

Modeling Fracture and Failure with Abaqus

L3.63

Through Cracks in Shells
? Comparing these results with the shell element results presented earlier: ? Errors with respect to the analytical solution for the 3D model are less than 1%.

? Much closer agreement because transverse shear effects are considered in the 3D model.
? Only in-plane stress and strain terms are included in the Abaqus J calculations for shells. ? Transverse shear terms are neglected.

Modeling Fracture and Failure with Abaqus

Mixed-Mode Fracture

L3.65

Mixed-Mode Fracture
? Abaqus uses interaction integrals to compute the stress intensity factors. ? ? This approach accounts for mixed-mode loading effects. Note that the J- or Ct-integrals do not distinguish between modes of loading.

?

Usage:
*CONTOUR INTEGRAL, TYPE=K FACTORS

?

Stress intensity factors can only be calculated for linear elastic materials.

Modeling Fracture and Failure with Abaqus

L3.66

Mixed-Mode Fracture
? Example: Center slant cracked plate under tension
?
22.5?
22.5? Element type

?? ??
0.185 (?2.9%)*
0.185 (?2.9%)

? ?? ??
0.403 (?0.2%)
0.403 (?0.2%)

CPE8
CPE8R

67.5?
67.5?

CPE8
CPE8R

1.052 (?3.6%)
1.053 (?3.8%)

0.373 (?1.0%)
0.374 (?1.3%)

K0 = ? p a
*Values enclosed in parentheses are percentage differences with respect to the reference solution. See Abaqus Benchmark Problem 4.7.4 for more information.

? = 22.5?

? = 67.5?
Modeling Fracture and Failure with Abaqus

Material Discontinuities

L3.68

Material Discontinuities
? The J-integral will be path independent if the material is homogeneous in the direction of crack propagation in the domain used for the contour integral calculation. ? If there is material discontinuity ahead of the crack in this region, the *NORMAL option can be used to correct the calculation of J so that it will still be path independent.

? The normal to the material discontinuity line must be specified for all nodes on the material discontinuity that will lie in a contour integral domain.

n

Modeling Fracture and Failure with Abaqus

L3.69

Material Discontinuities
? Example: J-integral analysis of a two material plate ? As an example, the figure shows a single-edge notch specimen made from two materials in which the material interface runs at an angle to the sides of the specimen.

? The material containing the crack (left) has a Young’s modulus of 2 ? 105 MPa and a Poisson’s ratio of 0.3.
? The uncracked material (right) has Young’s modulus of 2 ? 104 MPa and a Poisson’s ratio of 0.1. ? The specimen is stretched by uniform displacement at its ends.

Modeling Fracture and Failure with Abaqus

L3.70

Material Discontinuities ? J-integral analysis of a two material plate (cont’d)
? Along the material discontinuity, the normal to the discontinuity is given using the *NORMAL option.
? The normal needs to be defined on both sides of the discontinuity.
*NORMAL LEFT, NORM, 1.0, 0.125, 0.0 RIGHT, NORM, -1.0, -0.125, 0.0

Modeling Fracture and Failure with Abaqus

L3.71

Material Discontinuities
? The calculated J-integral values for 10 contours are as follows:
Contour 1 J (N/mm)

Without normals
55681

With normals
55681

2
3

57085
57052

57085
57052

4
5

57058
35188

57058
57116

6
7

31380
27536

57114
57114

8
9

23512
19172

57113
57116

10

14181

57094

? The need for the normals on the interface (contours 5 –10) is clear.
Modeling Fracture and Failure with Abaqus

Numerical Calculations with Elastic-Plastic Materials

L3.73

Numerical Calculations with Elastic-Plastic Materials
? For Mises plasticity the plastic deformation is incompressible.

? The rate of total deformation becomes incompressible (constant volume) as the plastic deformation starts to dominate the response. ? All Abaqus quadrilateral and brick elements suitable for use in J-integral calculations can handle this rate incompressibility condition except for the ―fully‖ integrated quadrilaterals and brick elements without the ―hybrid‖ formulation.
? Do not use CPE8, CAX8, C3D20 elements with these materials. They will ―lock‖ (become overconstrained) as the material becomes more incompressible.

Modeling Fracture and Failure with Abaqus

L3.74

Numerical Calculations with Elastic-Plastic Materials
? Second-order elements with reduced integration (CPE8R, C3D20R, etc.) work best for stress concentration problems in general and for crack tips in particular. ? If the displaced shape plot shows a regular pattern of deformation, this state is an indication of mesh locking. ? Locking can be seen in quilt contour plots of hydrostatic pressure for first-order elements—the pressure shows a checkerboard pattern.

? Change to reduced integration elements if you are using fully integrated elements. ? Increase the mesh density if you already using reduced integration elements. ? If these steps do not help, use hybrid elements.
? Hybrid elements must be used for fully incompressible materials (such as hyperelasticity, linear elasticity with n = 0.5).

Modeling Fracture and Failure with Abaqus

L3.75

Numerical Calculations with Elastic-Plastic Materials
? Results with elastic-plastic materials (and nonlinear materials in general) are more sensitive to meshing than for small-strain linear elasticity. ? Meshes adequate for linear elasticity may have to be refined. ? The more complex the solution, the more J values tend to be path dependent. ? A lack of path dependence can be an indication of a lack of mesh convergence; however, path independence of J does not prove mesh convergence.

Modeling Fracture and Failure with Abaqus

Workshop 1

L3.77

Workshop 1
? Crack in a three-point bend specimen

? Two-dimensional geometry
? Mesh sensitivity study

? Focus vs. unfocused mesh
? Quarter-point vs. mid-side nodes

Modeling Fracture and Failure with Abaqus

Workshop 2

L3.79

Workshop 2
? Crack in a helicopter airframe component

? Three-dimensional geometry
? Create mesh and evaluate response for cracks at different locations

Modeling Fracture and Failure with Abaqus

Material Failure and Wear
Lecture 4

L4.2

Overview
? Progressive Damage and Failure

? Damage Initiation for Ductile Metals
? Damage Evolution

? Element Removal
? Damage in Fiber-Reinforced Composite Materials

? Failure in Fasteners
? Material Wear and Ablation

Modeling Fracture and Failure with Abaqus

Progressive Damage and Failure

L4.4

Progressive Damage and Failure
? Abaqus offers a general capability for modeling progressive damage and failure in engineering structures ? Material failure refers to the complete loss of load carrying capacity that results from progressive degradation of the material stiffness. ? Stiffness degradation is modeled using damage mechanics. ? Progressive damage and failure can be modeled for: ? Ductile materials

? Continuum constitutive behavior
? Fiber-reinforced composites

? Interface materials
? Cohesive elements with a traction-separation law

? Damage and failure of cohesive elements are discussed in the next lecture.

Modeling Fracture and Failure with Abaqus

L4.5

Progressive Damage and Failure
? Two distinct types of ductile material failure can be modeled with Abaqus ? Ductile fracture of metals ? Void nucleation, coalescence, and growth ? Shear band localization ? Necking instability in sheet-metal forming ? Forming Limit Diagrams ? Marciniak-Kuczynski (M-K) criterion ? Damage in sheet metals is not discussed further in this seminar.

Modeling Fracture and Failure with Abaqus

L4.6

Progressive Damage and Failure
? Components of material definition

? Undamaged constitutive behavior (e.g., elastic-plastic with hardening)
? Damage initiation (point A) ? Damage evolution (path A–B) ? Choice of element removal (point B)
Keywords

?

Undamaged response

A
Damaged response

*MATERIAL
*ELASTIC *PLASTIC *DAMAGE INITIATION,CRITERION=criterion *DAMAGE EVOLUTION *SECTION CONTROLS, ELEMENT DELETION=YES
Modeling Fracture and Failure with Abaqus

B
Multiple damage definitions are allowed

?

Typical material response showing progressive damage

Damage Initiation Criteria for Ductile Metals

L4.8

Damage Initiation Criteria for Ductile Metals
? Damage initiation defines the point of initiation of degradation of stiffness ? It is based on user-specified criteria

? Ductile or shear
? It does not actually lead to damage unless damage evolution is also specified ? Output variables associated with each criterion ? Useful for evaluating the severity of current deformation state

? Output
DMICRT DMICRT > 1 indicates damage has initiated
Ductile Shear

Different damage initiation criteria on an aluminum double-chamber profile

Modeling Fracture and Failure with Abaqus

L4.9

Damage Initiation Criteria for Ductile Metals
? Ductile criterion:

? Appropriate for triggering damage due to nucleation, growth, and coalescence of voids
? The model assumes that the equivalent plastic strain at the onset of damage is a function of stress triaxiality and strain rate.
Pressure stress

? Stress triaxiality h = - p / q
Mises stress

? The ductile criterion can be used with the Mises, Johnson-Cook, Hill, and Drucker-Prager plasticity models, including equation of state.

Ductile criterion for Aluminum Alloy AA7108.50-T6 (Courtesy of BMW)

Modeling Fracture and Failure with Abaqus

L4.10

Damage Initiation Criteria for Ductile Metals
? Usage:

? Specify the equivalent plastic strain at the onset of damage as a tabular function of
? Stress triaxiality ? Strain rate
*DAMAGE INITIATION, CRITERION=DUCTILE

? pl , T , f ? pl , h , ? i
Equivalent fracture strain at damage initiation Temperature and field variable dependence optional The criterion for damage initiation is met when wD = 1.

? Output:
DUCTCRT (wD)

Modeling Fracture and Failure with Abaqus

L4.11

Damage Initiation Criteria for Ductile Metals
? Shear criterion:

? Appropriate for triggering damage due to shear band localization
? The model assumes that the equivalent plastic strain at the onset of damage is a function of the shear stress ratio and strain rate. ? Shear stress ratio defined as:

qs = (q + ks p) /tmax
? The shear criterion can be used with the Mises, Johnson-Cook, Hill, and Drucker-Prager plasticity models, including equation of state.
ks = 0.3

Shear criterion for Aluminum Alloy AA7108.50-T6 (Courtesy of BMW)

Modeling Fracture and Failure with Abaqus

L4.12

Damage Initiation Criteria for Ductile Metals
? Usage:

? Specify the equivalent plastic strain at the onset of damage as a tabular function of
? Shear stress ratio ? Strain rate
*DAMAGE INITIATION, CRITERION=SHEAR, KS=ks

? pl , T , f ? pl , q s , ? i
Equivalent fracture strain at damage initiation Temperature and field variable dependence optional

ks is a material parameter

? Output:
SHRCRT (wS)

The criterion for damage initiation is met when wS = 1.

Modeling Fracture and Failure with Abaqus

L4.13

Damage Initiation Criteria for Ductile Metals
? Example: Axial crushing of an aluminum double-chamber profile

Cross section

Quasi-static buckling mode

Modeling Fracture and Failure with Abaqus

L4.14

Damage Initiation Criteria for Ductile Metals
? Model details

? Steel base:
? C3D8R elements

? Enhanced hourglass control
? Elastic-plastic material

Rigid plate with initial downward velocity Aluminum chamber

? Aluminum chamber:
? S4R elements ? Stiffness hourglass control ? Rate-dependent plasticity ? Damage initiation

? General contact
? Variable mass scaling
Steel base: bottom is encastred.

Modeling Fracture and Failure with Abaqus

L4.15

Damage Initiation Criteria for Ductile Metals
? Material definition : Keywords interface
*MATERIAL, NAME=ALUMINUM T6 (Courtesy of BMW) 7 *DENSITY strain rate=0.001/s 6 2.70E-09 strain rate=250/s 5 *ELASTIC 4 7.00E+04, 0.33 *PLASTIC,HARDENING=ISOTROPIC,RATE=0 3 : 2 *DAMAGE INITIATION, CRITERION=DUCTILE 1 5.7268, 0.000, 0.001 0 0 0.2 0.4 0.6 4.0303, 0.067, 0.001 stress triaxiality 2.8377, 0.133, 0.001 : ? pl Strain rate, ? 4.4098, 0.000, 250 2.5717, 0.067, 250 Stress triaxiality, h 1.5018, 0.133, 250 : Equivalent fracture strain at damage initiation, pl
strain at damage initiation

Ductile criteria for Aluminum Alloy AA7108.50-

?

Modeling Fracture and Failure with Abaqus

L4.16

Damage Initiation Criteria for Ductile Metals
? Material definition : Keywords interface (cont'd)
*MATERIAL, NAME=ALUMINUM : *DAMAGE INITIATION, CRITERION=DUCTILE 5.7268, 0.000, 0.001 4.0303, 0.067, 0.001 : *DAMAGE INITIATION, CRITERION=SHEAR, KS=0.3 0.2761, 1.424, 0.001 0.2613, 1.463, 0.001 0.2530, 1.501, 0.001 : 0.2731, 1.424, 250 Strain rate, ? pl
Shear criteria for Aluminum Alloy AA7108.50-T6 (Courtesy of BMW)
strain at damage initiation

0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 1.6 1.7

strain rate=0.001/s strain rate=250/s 1.8 1.9 2

shear stress ratio

0.3025, 1.463, 250 0.3323, 1.501, 250 :

?

Shear stress ratio,

qs

Equivalent fracture strain at damage initiation, ? pl

Modeling Fracture and Failure with Abaqus

L4.17

Damage Initiation Criteria for Ductile Metals
? Material definition : Abaqus/CAE interface

: *DAMAGE INITIATION, CRITERION=DUCTILE 5.7268, 0.000, 0.001 4.0303, 0.067, 0.001 2.8377, 0.133, 0.001 : 4.4098, 0.000, 250 2.5717, 0.067, 250 1.5018, 0.133, 250 :

Modeling Fracture and Failure with Abaqus

L4.18

Damage Initiation Criteria for Ductile Metals
? Material definition : Abaqus/CAE interface (cont'd)

: *DAMAGE INITIATION, CRITERION=SHEAR, KS=0.3 0.2761, 1.424, 0.001 0.2613, 1.463, 0.001 0.2530, 1.501, 0.001 : 0.2731, 1.424, 250 0.3025, 1.463, 250 0.3323, 1.501, 250 :

Modeling Fracture and Failure with Abaqus

L4.19

Damage Initiation Criteria for Ductile Metals
? Results (without damage evolution)

Ductile
Quasi-static response

Shear

Modeling Fracture and Failure with Abaqus

Damage Evolution

L4.21

Damage Evolution
? Damage evolution defines the post damage-initiation material behavior. ? That is, it describes the rate of degradation of the material stiffness once the initiation criterion is satisfied. ? The formulation is based on scalar damage approach:
Stress due to undamaged response

? = (1 - d )?

? The overall damage variable d captures the combined effect of all active damage mechanisms. ? When damage variable d = 1, material point has completely failed. ? In other words, fracture occurs when d = 1.

Modeling Fracture and Failure with Abaqus

L4.22

Damage Evolution
? Elastic-plastic materials
Undamaged response

? For a elastic-plastic material, damage manifests in two forms
? Softening of the yield stress ? Degradation of the elasticity

?
? y0

?
(d = 0)
- d?
Degradation of elasticity softening

? The strain softening part of the curve cannot represent a material property.
? The above argument is based on ? Fracture mechanics considerations ? Mesh sensitivity

?0
E

E

(1 - d ) E

? 0pl

? fpl

?

Schematic representation of elastic-plastic material with progressive damage.

Modeling Fracture and Failure with Abaqus

L4.23

Damage Evolution
? To address the strain softening issue, Hillerborg’s (1976) proposal is adopted. ? The fracture energy to open a unit area of crack, Gf , is assumed to be a material property. ? The softening response after damage initiation is characterized by a stress-displacement response (rather than a stress-strain response) ? This requires the introduction of a characteristic length L associated with a material point.

Modeling Fracture and Failure with Abaqus

L4.24

Damage Evolution
? The fracture energy is written as

Gf =

??

? fpl
pl 0

? pl = L? y?

?

u fpl 0

? pl ? yu

where u pl is the equivalent plastic displacement.
? The characteristic length L is computed automatically by Abaqus based on element geometry. ? Elements with large aspect ratios should be avoided to minimize mesh sensitivity. ? The damage evolution law can be specified either in terms of fracture energy (per unit area) or in terms of the equivalent plastic displacement. ? Both approaches take into account the characteristic length of the element. ? The formulation ensures that mesh-sensitivity is minimized.

Modeling Fracture and Failure with Abaqus

L4.25

Damage Evolution
? Displacement-based damage evolution

d
1
0
(a) Tabular

d
1
0

d
1
0

u

pl

u

pl f

u

pl

(b) Linear

(c) Exponential

u fpl u pl

*DAMAGE EVOLUTION,TYPE=DISPLACEMENT, SOFTENING={TABULAR,LINEAR,EXPONENTIAL}

Modeling Fracture and Failure with Abaqus

L4.26

Damage Evolution
? Procedure for generating d vs u pl table from tensile test data 1. Plot true stress, ? vs. total displacement u measured over the gauge length L For stress values in the softening branch (i.e. beyond damage initiation), compute damage parameter d from the expression ? = (1 - d )? Compute the corresponding plastic displacement u pl as shown in the schematic.
Undamaged response

?
? y0

d = 0; u = 0
pl
u
pl f

?
- d?
softening

2.

?0
E L E L
(1-d ) E L

3.

u
pl pl f

u pl

u

pl f

d = 1; u = u

4.

In the absence of intermediate data, choose linear softening and provide value of

Schematic representation of tensile test data in stress – displacement space for elastic-plastic materials

Modeling Fracture and Failure with Abaqus

L4.27

Damage Evolution
? Energy-based damaged evolution
?y ? y0
Gf
(a) Linear

u fpl =

2G f

? y0

?y ? y0
Gf

NOTE: The response is linear or exponential only if the undamaged response is perfectly plastic

u fpl u pl

(b) Exponential

u pl

*DAMAGE EVOLUTION,TYPE=ENERGY, SOFTENING={LINEAR,EXPONENTIAL}

Modeling Fracture and Failure with Abaqus

L4.28

Damage Evolution
? Example: Tearing of an X-shaped cross section
Tie constraints

Fix this end

Pull and twist this this end

Failure modeled with different mesh densities

*damage initiation, criterion=fld 0.20, *damage evolution, type=displacement, softening=tabular 0.0, 0.0 1.0, 0.003 damage-plastic displacement data pairs
Modeling Fracture and Failure with Abaqus

L4.29

Damage Evolution
? Comparison of reaction forces and moments confirms mesh insensitivity of the results.

Modeling Fracture and Failure with Abaqus

L4.30

Damage Evolution
? Example: Axial crushing of an aluminum double-chamber profile

? Dynamic response with damage evolution

*Material, name=Aluminum : *Damage initiation, criterion=Ductile : *Damage evolution, type=displacement 0.1, *Damage initiation, criterion=Shear, ks=0.3 : *Damage evolution, type=displacement 0.1,

Modeling Fracture and Failure with Abaqus

L4.31

Damage Evolution
? With damage evolution, the simulation response is a good approximation of the physical response.

Simulation without damage evolution

Aluminum double-chamber after dynamic impact

Simulation with damage evolution

Modeling Fracture and Failure with Abaqus

Element Removal

L4.33

Element Removal
? Abaqus offers the choice to remove the element from the mesh once the material stiffness is fully degraded (i.e., once the element has failed). ? An element is said to have failed when all section points at any one integration point have lost their load carrying capacity. ? By default, failed elements are deleted from the mesh.

Modeling Fracture and Failure with Abaqus

L4.34

Element Removal
? Removing failed elements before complete degradation

? The material point is assumed to fail when the overall damage variable D reaches the critical value Dmax. ? You can specify the value for the maximum degradation Dmax.
? The default value of Dmax is 1 if the element is to be removed from the mesh upon failure.

Modeling Fracture and Failure with Abaqus

L4.35

Element Removal
? Usage:
*SECTION CONTROLS, NAME=Ec-1, ELEMENT DELETION=YES, MAX DEGRADATION=0.9

:
** Refer to the section controls by name on the element section definition. *SOLID SECTION, ELSET=Elset_1, CONTROLS=Ec-1, MATERIAL=Material_1 :

Modeling Fracture and Failure with Abaqus

L4.36

Element Removal
? Retaining failed elements

? You may choose not to remove failed elements from the mesh.
*SECTION CONTROLS, ELEMENT DELETION = NO

? In this case the default value of Dmax is 0.99, which ensures that elements will remain active in the simulation with a residual stiffness of at least 1% of the original stiffness. ? Here Dmax represents
? the maximum degradation of the shear stiffness (three-dimensional), ? the total stiffness (plane stress), or ? the uniaxial stiffness (one-dimensional).

? Failed elements that have not been removed from the mesh can sustain hydrostatic compressive stresses.
Modeling Fracture and Failure with Abaqus

L4.37

Element Removal
? Output

? The output variable SDEG contains the value of D. ? The output variable STATUS indicates whether or not an element has failed. ? STATUS = 0 for failed elements ? STATUS = 1 for active elements ? Abaqus/Viewer will automatically remove failed elements when the output database (.odb) file includes STATUS.

Failed elements removed by default when STAUS output is available

failed elements

Deactivate status variable to view failed elements
Modeling Fracture and Failure with Abaqus

Damage in Fiber-Reinforced Composite Materials

L4.39

Damage in Fiber-Reinforced Composite Materials
? Abaqus offers a general capability for modeling progressive damage and failure in fiber-reinforced composites. ? Material failure refers to the complete loss of load carrying capacity that results from progressive degradation of the material stiffness. ? Stiffness degradation is modeled using damage mechanics. ? The model must be used with elements with a plane stress formulation (plane stress, shell, continuum shell, and membrane elements) ? Four different modes of failure are considered: ? fiber rupture in tension; ? fiber buckling and kinking in compression; ? matrix cracking under transverse tension and shearing; and ? matrix crushing under transverse compression and shearing

Common damage types in composite laminates

Modeling Fracture and Failure with Abaqus

L4.40

Damage in Fiber-Reinforced Composite Materials
? User interface

? Damage Initiation

*DAMAGE INITIATION, CRITERION=HASHIN, ALPHA=<alpha> XT,XC,YT,YC,SL,ST

Modeling Fracture and Failure with Abaqus

L4.41

Damage in Fiber-Reinforced Composite Materials
? Damage Evolution
*DAMAGE EVOLUTION, TYPE=ENERGY, SOFTENING=LINEAR

Gft,Gfc ,Gmt,Gmc

? Viscous Regularization
*DAMAGE STABILIZATION
ηft, ηfc, ηmt, ηmc

Modeling Fracture and Failure with Abaqus

L4.42

Damage in Fiber-Reinforced Composite Materials
? Output

? Initiation Criteria Variables
? HSNFTCRT – tensile fiber Hashin’s criterion

? HSNFCCRT – compressive fiber Hashin’s criterion
? HSNMTCRT – tensile matrix Hashin’s criterion

? HSNMCCRT – compressive matrix Hashin’s criterion
? Damage Variables ? DAMAGEFT – tensile fiber damage ? DAMAGEFC – compressive fiber damage ? DAMAGEMT – tensile matrix damage

? DAMAGEMC – compressive matrix damage

Modeling Fracture and Failure with Abaqus

L4.43

Damage in Fiber-Reinforced Composite Materials
? Output (cont'd)

? Status
? STATUS – element status (1 – present, 0 – removed)

?

Energies
? Damage energy (ALLDMD,DMENER,ELDMD,EDMDDEN)

? Viscous regularization (ALLCD, CENER, ELCD, ECDDEN)

Modeling Fracture and Failure with Abaqus

L4.44

Damage in Fiber-Reinforced Composite Materials
? Example: Analysis of blunt notched fiber metal laminate

? Fiber metal laminates (FMLs) are composed of:
? laminated thin aluminum layers

? Intermediate glass fiber-reinforced epoxy layers

Modeling Fracture and Failure with Abaqus

L4.45

Damage in Fiber-Reinforced Composite Materials
? Geometry of blunt notched fiber metal laminate (Glare 3 3/2 –0.3)
1/8 part model

Aluminum core and exterior

a through-thickness hole

? Through-thickness view of the laminate:

glass fiber-reinforced epoxy layers

Example Problem 1.4.6, "Failure of blunt notched fiber metal laminates”
Modeling Fracture and Failure with Abaqus

L4.46

Damage in Fiber-Reinforced Composite Materials
? Results

damage in matrix

and

damage in fibers Net blunt notch strength (MPa)
Test (De Vries, 2001) Abaqus 446 453

for one of glass fiber-reinforced epoxy layers

Modeling Fracture and Failure with Abaqus

L4.47

Damage in Fiber-Reinforced Composite Materials
? Abaqus allows the import of the damage model for fiber-reinforced composites from Abaqus/Explicit to Abaqus/Standard. ? Details of the import capability will not be covered in this lecture (please refer to ―Importing and transferring results,‖ Section 9.2 of the Abaqus Analysis User’s Manual). ? One typical application is the analysis of Barely Visible Impact Damage (BVID) in composite structures used in aerospace applications. ? Non-visible damage to composite structures is a significant concern in the aerospace industry.

from McGowan, D.M., and Ambur, D.R., NASA TM-110303

Damage-Tolerance Characteristics of Composite Fuselage Sandwich Structures With Thick Facesheets

Modeling Fracture and Failure with Abaqus

Damage in Fasteners

L4.49

Damage in Fasteners
? Connection methodologies—point fasteners

? Fastener (spot weld) compliance and failure are available in Abaqus.

multiple layers

attachment points

radius of influence

Modeling Fracture and Failure with Abaqus

L4.50

Damage in Fasteners
? Fastener failure

? Model combines plasticity and progressive damage
– Response depends on loading angle (normal/shear) – Stages: ? Rigid plasticity with variable hardening ? Damage initiation
F
0
?

S ? 0?

N ? 90?

Spot weld
Plasticity + Damage

? Progressive damage evolution using fracture energy

45?
90?
Plasticity
damage initiation boundary

Schematic representation of the predicted numerical response
Modeling Fracture and Failure with Abaqus

u pl

L4.51

Damage in Fasteners
? Example

? Spot-welded hat section of three layers of sheet metals subjected to severe compressive loading
Deformable fastener still holding

Failed fasteners

Rigid spot welds

Compliant spot welds with damage

Modeling Fracture and Failure with Abaqus

Material Wear and Ablation

L4.53

Material Wear and Ablation
? Material wear/erosion in Abaqus/Standard

? Many applications require the modeling of wear/erosion of material at one or more surfaces
? Capability enables modeling of material wear/erosion on the surface of the body ? Idea is to erode material while receding mesh away from surface (with same number and topology of elements) ? Involves remeshing, state mapping—handled through an Arbitrary Lagrangian-Eulerian (ALE) technique ? User interface takes advantage of existing adaptive meshing framework to define mesh motion

Adaptive mesh domain for modeling material wear. Wear extent/velocity applied as mesh constraints

Modeling Fracture and Failure with Abaqus

L4.54

Material Wear and Ablation
? Applications

? Geotechnical
? Well bore sand production

? Plastic strain, fluid velocity
? Aerospace

? Rocket motor ablation
? Pyrolysis, char formation ? Solid propellants ? Automotive ? Tire wear

? Disk brake wear
? Manufacturing
Fluid velocity dependent wear of a well bore

? Machining
Modeling Fracture and Failure with Abaqus

L4.55

Material Wear and Ablation
? User interface
*Adaptive mesh, elset=... *Adaptive mesh constraint, type=[velocity|displacement], User *Adaptive mesh controls

? Adaptive mesh constraints define mesh motion (wear extent or velocity)

? Wear criterion
? General descriptions possible through user subroutine UMESHMOTION

? User access to solution variables
? Nodal

? Material
? Contact

? A local surface coordinate system is provided

Modeling Fracture and Failure with Abaqus

L4.56

Material Wear and Ablation
? Example of wear criterion

? Tire wear
? Use of CSLIP, CSHEAR, CPRESS

? =?E ? h
Rate of recession of tread

Rate of frictional energy dissipation Proportionality constant

h

Modeling Fracture and Failure with Abaqus

L4.57

Material Wear and Ablation
? Example: erosion of material from oil bore hole perforation tunnel ? Setup consists of bore hole with perforations, loaded by weight of material above

? Pore pressure gradient leads to flow into perforation ? Material wear rate controlled by fluid flux, transport concentration, porosity, sand production coefficient, and the local plastic deformation ? Optimum design to minimize wear rate ? Example Problem 1.1.22

Perforation tunnel

Bore hole
Geometry of oil well

Courtesy of Exxon

Modeling Fracture and Failure with Abaqus

L4.58

Material Wear and Ablation
? Analysis steps

? Geostatic
? Model change removal of well bore and casing (drilling operation)

? Apply pore pressure; establish steady state conditions
? Transient soils consolidation (during which the erosion occurs)

? Ablation relation:
V = 10 × (PEEQ - 0.028)

Erosion velocity

Modeling Fracture and Failure with Abaqus

L4.59

Material Wear and Ablation
? Adaptive mesh constraints
*Adaptive mesh, elset=Adaptive-Zone, Freq=1, Mesh=40 *Adaptive mesh constraint, constraint type=Lagrangian

Lag
*Adaptive mesh constraint, type=velocity, user

Rock-Perf, 1, 1, 1.0

Lag: Nodes on back face of adaptive domain Rock-Perf

Adaptive-Zone

Cut section of the adaptive mesh domain showing the perforation tunnel

Modeling Fracture and Failure with Abaqus

L4.60

Material Wear and Ablation
? User subroutine
c
include 'aba_param.inc' c parameter parameter dimension dimension dimension dimension dimension (zero=0.d0, ten=10.d0, peeqCrit=0.028d0) (nelemmax=100) array(1000) ulocal(*) jgvblock(*),jmatyp(*) alocal(ndim,*) jelemlist(nelemmax),jelemtype(nelemmax) subroutine umeshmotion(uref,ulocal,node,nndof,lnodetype,alocal, $ ndim,time,dtime,pnewdt,kstep,kinc,kmeshsweep,jmatyp,jgvblock,lsmooth)

locnum = 0 jtyp = 1 peeq = zero nelems = nelemmax call getNodeToElemConn(node,nelems,jelemlist, $ jelemtype,jrcd,jgvblock) call getVrmAvgAtNode(node, jtyp, 'PE', array, jrcd, $ jelemlist, nelems, jmatyp, jgvblock) peeq = array(7)

When NDIM=3 the 3-direction is normal to the surface

if (peeq .gt. peeqCrit) then ulocal(ndim) = ulocal(ndim)- ten*(peeq - peeqCrit) end if return ulocal passed in as the value end the mesh smoothing algorithm

determined by

Modeling Fracture and Failure with Abaqus

L4.61

Material Wear and Ablation
? Results

Material wear at bore hole/perforation junction

Total volume lost due to erosion is available with history output variable VOLC

Modeling Fracture and Failure with Abaqus

L4.62

Material Wear and Ablation
? Mesh smoothing

? Two options
? Original configuration projection method ? Smoothing performed according to the original configuration Original-configuration ? Volume-based smoothing
smoothing

? Either method can include a geometric-based enhancement

Volumetric smoothing

Modeling Fracture and Failure with Abaqus

L4.63

Material Wear and Ablation
? Smoothing permitted in conjunction with UMESHMOTION constraints

? Enables UMESHMOTION to describe normal mesh motions, while the smoothing algorithm handles the tangential mesh motions.

Modeling Fracture and Failure with Abaqus

L4.64

Material Wear and Ablation
? Limitations

? Available for a subset of continuum elements
? Available only for following procedures using geometric nonlinearity

? Static
? Soils

? Coupled Temperature-Displacement
? Tracer particles not supported

Modeling Fracture and Failure with Abaqus

Element-based Cohesive Behavior
Lecture 5

L5.2

Overview
? Introduction

? Element Technology
? Constitutive Response

? Viscous Regularization
? Modeling Techniques

? Examples
? Workshop 3 (Part 1) ? Workshop 4

Modeling Fracture and Failure with Abaqus

L5.3

Overview
? Historical perspective

? The concept of a cohesive zone has been around for some time:
? Dugdale (1960) and Barenblatt (1962) were the first to apply the concept of a cohesive stress zone to fracture modeling. ? Many extensions since then.

? For example, Needleman (1987) recognized that cohesive elements are particularly attractive when interface strengths are relatively weak compared to the adjoining materials.
? Examples: composite laminates and parts bonded with adhesives

Modeling Fracture and Failure with Abaqus

Introduction

L5.5

Introduction
? Cohesive behavior is useful in modeling adhesives, bonded interfaces, and gaskets. ? Models separation between two initially bonded surfaces

? Progressive failure of adhesives
? Delamination in composites

? Idealize complex fracture mechanisms with a macroscopic “cohesive law,” which relates the traction across the interface to the separation.
? The cohesive behavior can be: ? Element-based ? Modeled with cohesive elements ? Surface-based ? Modeled with contact pairs in Abaqus/Standard and general contact in Abaqus/Explicit

T-peel analysis: Cohesive elements are used for modeling adhesive patches

Failed adhesive is red
(CSDMG = 1)

Rail crush: Cohesive surfaces
Modeling Fracture and Failure with Abaqus

L5.6

Introduction
? Element-based cohesive behavior—cohesive elements

? Cohesive elements allow very detailed modeling of adhesive connections, including
? specification of detailed adhesive material properties, direct control of the connection mesh, modeling of adhesives of finite thickness, etc. ? Cohesive elements in Abaqus primarily address two classes of problems:

? Adhesive joints
? Adhesive layer with finite thickness ? Typically the bulk material properties are known ? Delamination ? Adhesive layer of “zero” thickness

? Typically the bulk material properties are not known
Modeling Fracture and Failure with Abaqus

L5.7

Introduction
? The constitutive modeling depends on the class of problem:

? Based on macroscopic properties (stiffness, strength) for adhesive joints ? Continuum description: any Abaqus material model can be used
? Modeling technique is relatively straightforward: cohesive layer has finite thickness; standard material models (including damage). ? The continuum description is not discussed further in this lecture.

? Based on a traction-separation description for delamination
? Linear elasticity with damage

? Modeling technique is less straightforward: typical applications use zero-thickness cohesive elements; non-standard constitutive law ? This application is the primary focus of this lecture

Modeling Fracture and Failure with Abaqus

L5.8

Introduction
? In addition, the uniaxial response of a laterally unconstrained adhesive patch can also be modeled ? This represents the behavior of a gasket. ? Limited capability for modeling gaskets with cohesive elements: ? The complexity of the response in the thickness direction is not as rich as with gasket elements available in Abaqus/Standard. ? Compared to gasket elements, however, cohesive elements:

? are fully nonlinear (can be used with finite strains and rotations);
? can have mass in a dynamic analysis; and ? are available in both Abaqus/Standard and Abaqus/Explicit.

? The use of cohesive elements for modeling gaskets is not discussed further in this lecture.

Modeling Fracture and Failure with Abaqus

L5.9

Introduction
? Surface-based cohesive behavior—cohesive surfaces ? This is a simplified and easy way to model cohesive connections, using the traction-separation interface behavior. ? It offers capabilities that are very similar to cohesive elements modeled with the traction-separation constitutive response. ? However, it does not require element definitions.

? In addition, cohesive surfaces can bond anytime contact is established (“sticky” contact) ? It is primarily intended for situations in which interface thickness is negligibly small. ? It must be defined as a surface interaction property. ? Damage for cohesive surfaces is an interaction property, not a material property.
? The kinematics of cohesive surfaces is different from that of cohesive elements. ? By default, the initial stiffness of the interface is computed automatically.
Modeling Fracture and Failure with Abaqus

L5.10

Introduction
? Cohesive elements are the focus of this lecture.

? Cohesive surfaces are discussed in the next lecture.
? A workshop exercise will allow you to compare and contrast the two cohesive modeling techniques in the context of a simple problem.

Modeling Fracture and Failure with Abaqus

Element Technology

L5.12

Element Technology
? Element types*
Top face

? 3D elements
? COH3D8

? COH3D6
? 2D element

Bottom face

? COH2D4
? Axisymmetric element ? COHAX4 ? These elements can be embedded in a model via ? shared nodes or ? tie constraints.
*Cohesive pore pressure elements are also available.
Modeling Fracture and Failure with Abaqus

L5.13

Element Technology
? Element and section definition
*ELEMENT, TYPE = COH3D8 *COHESIVE SECTION, ELSET =..., RESPONSE = {TRACTION SEPARATION, CONTINUUM, GASKET },

THICKNESS = { SPECIFIED, GEOMETRY},
MATERIAL = ...

Specify thickness in dataline (default is 1.0)

Modeling Fracture and Failure with Abaqus

L5.14

Element Technology
? Default thickness of cohesive elements

? Traction-separation response:
? Unit thickness

? Continuum and gasket response
? Geometric thickness based on nodal coordinates

Modeling Fracture and Failure with Abaqus

L5.15

Element Technology
? Output variables

? Scalar damage (i.e., degradation) variable
? SDEG

? Variables indicating whether damage initiation criteria met or exceeded
? Discussed shortly

? Element status flag
? STATUS

Modeling Fracture and Failure with Abaqus

L5.16

Element Technology
? Import of cohesive elements

? The combination of Abaqus/Standard and Abaqus/Explicit expands the range of applications for cohesive elements.
? For example, you can simulate the damage in a structure due to an impact event then study the effect of the damage on the structure's load carrying capacity.

Modeling Fracture and Failure with Abaqus

Constitutive Response

L5.18

Constitutive Response
? Delamination applications
T
N

? Traction separation law
? Typically characterized by peak strength (N) and fracture energy (GTC) ? Mode dependent ? Linear elasticity with damage ? Available in both Abaqus/Standard and Abaqus/Explicit

GT C
?
Typical traction-separation response
7 6 5

? Modeling of damage under the general framework introduced earlier ? Damage initiation
? Traction or separation-based criterion ? Damage evolution
GTC

Shear mode
Normal mode

4 3 2 1 0

? Removal of elements

0

0.2

0.4

0.6

0.8

1

Mode Mix

Dependence of fracture toughness on mode mix
Modeling Fracture and Failure with Abaqus

L5.19

Constitutive Response
? Linear elasticity with damage

? Linear elasticity
? Defines behavior before the initiation of damage ? Relates nominal stress to nominal strain ? Nominal traction to separation with default choice of unit thickness ? Uncoupled traction behavior: nominal stress depends only on corresponding nominal strain

? Coupled traction behavior is more general

*ELASTIC, TYPE = { TRACTION, COUPLED TRACTION }

Modeling Fracture and Failure with Abaqus

L5.20

Constitutive Response
? The elastic modulus for the traction separation law should be interpreted as a penalty stiffness.

N
N max

? For example, for the opening mode:

Kn = Nmax / ?ninit

? In Abaqus, nominal stress and strain quantities are used for the traction separation law.
? If unit thickness is specified for the element, then the nominal strain corresponds to the separation value. ? Elastic response governed by Kn. ? If you specify a non-unit thickness for the cohesive element, you must scale your data to obtain the correct stiffness Kn. Example on next slide.

Kn
1
init ?n

? nfail

?n

Displacement at damage initiation in normal (opening) mode

Modeling Fracture and Failure with Abaqus

L5.21

Constitutive Response
? Example: Peel test model
N = En? n = K n? n
Abaqus evaluates this… …which is equivalent to this

? n = ? n / heff ? K n = En / heff
init Assume separation at initiation ? n = 1e-3 and Nmax = 6.9e9.

A
En=Kn?heff

For model A: use geometric thickness
init heff = hgeom =1e-3?? n = ? n /heff = 1; init

Nmax = En = 6.9e9 ? Kn = 6.9e12 For model B: specify unit thickness heff = 1 ?? n = ? n / heff =1e-3;
init init

Nmax = 6.9e9 ? En = Kn = 6.9e12

B
Geometric thickness (based on nodal coordinates) of the adhesive hgeom = 1e-3
Modeling Fracture and Failure with Abaqus

L5.22

Constitutive Response
? Damage initiation ? Mixed mode conditions ? Maximum stress (or strain) criterion:

? ?n ?t ? s ? MAX ? , , ? =1 ? N max Tmax Smax ?

?n

?? n for ? n ? 0 =? ?0 for ? n ? 0

? Output:

? MAXSCRT
? MAXECRT
* DAMAGE INITIATION, CRITERION = { MAXS, MAXE }

Modeling Fracture and Failure with Abaqus

L5.23

Constitutive Response
? For example, for Mode I (opening mode) the MAXS condition implies damage initiates when ?n = Nmax.

N N max

Damage initiation point

*Damage initiation,criterion=MAXS 290.0E6, 200.0E6, 200.0E6

?n
Nmax
Tmax Smax

Modeling Fracture and Failure with Abaqus

L5.24

Constitutive Response
? Quadratic stress (or strain) interaction criterion:

? ?n ? ? ?t ? ? ? s ? ? ? ?? ? ?? ? =1 ? N max ? ? Tmax ? ? Smax ?
? No damage initiation under pure compression ? Output:
? QUADSCRT

2

2

2

? QUADECRT
* DAMAGE INITIATION, CRITERION = { QUADS, QUADE }

Modeling Fracture and Failure with Abaqus

L5.25

Constitutive Response
? Summary of damage initiation criteria
Maximum nominal stress criterion Maximum nominal strain criterion

? ? ? ? ? MAX ? n , s , t ? = 1 ? N max Smax Tmax ?
*DAMAGE INITIATION, CRITERION=MAXS

? ?s ?t ? ? ?n ? MAX ? max , max , max ? = 1 ?s ?t ? ? ?? n ?
*DAMAGE INITIATION, CRITERION=MAXE
max max ?n , ?s , ? tmax

N max , S max , Tmax

Quadratic nominal stress criterion

Quadratic nominal stress criterion

? ?n ? ? ? s ? ? ?t ? ? ? ?? ? ?? ? =1 ? N max ? ? Smax ? ? Tmax ?
*DAMAGE INITIATION, CRITERION=QUADS

2

2

2

? ? n ? ? ? s ? ? ?t ? ? max ? ? ? max ? ? ? max ? = 1 ? ? n ? ? ? s ? ? ?t ?
*DAMAGE INITIATION, CRITERION=QUADE
max max ?n , ?s , ? tmax

2

2

2

N max , S max , Tmax
?n: nominal stress in the pure normal mode ?s: nominal stress in the first shear direction ?t: nominal stress in the second shear direction
Note : ? n =

?n: nominal strain in the pure normal mode ?s: nominal strain in the first shear direction ?t: nominal strain in the second shear direction

?n
To

, ?s =

?s
To

, ?t =

?t
To

where ?n,

?s, and ?t are components of relative displacement between the top and bottom of the cohesive element; and To

is the original thickness of the cohesive element. Modeling Fracture and Failure with Abaqus

L5.26

Constitutive Response
? Damage evolution ? Post damage-initiation response defined by:
?
?
-d?
(1 - d )?
K0
(1 - d ) Κ 0

? = ?1 - d ??
? d is the scalar damage variable

d = 0: undamaged
d = 1: fully damaged

?
K0
Typical damaged response

d monotonically increases

Modeling Fracture and Failure with Abaqus

L5.27

Constitutive Response
? Damage evolution is based on energy or displacement ? Specify either the total fracture energy or the post damage-initiation effective displacement at failure ? May depend on mode mix

N
N max
Area under the curve is the fracture energy

GT C

? Mode mix may be defined in terms of energy or traction

?n
Displacement at failure ? n in normal (opening) mode
fail

Modeling Fracture and Failure with Abaqus

L5.28

Constitutive Response
? Displacement-based damage evolution ? Damage is a function of an effective displacement:

Traction
Linear postinitiation response

?=

?n

2

? ? s2 ? ? t2

? The post damage-initiation softening response can be either ? Linear ? Exponential ? Tabular

? init

? fail ?

Modeling Fracture and Failure with Abaqus

L5.29

Constitutive Response
? Keywords interface for displacement-based damage evolution
*DAMAGE EVOLUTION, TYPE = DISPLACEMENT, SOFTENING = { LINEAR | EXPONENTIAL | TABULAR }, MIXED MODE BEHAVIOR = TABULAR

? For LINEAR and EXPONENTIAL softening:

? Specify the effective displacement at complete failure ?fail relative to the effective displacement at initiation ?init. ? For TABULAR softening:
? Specify the scalar damage variable d directly as a function of ? –?init. ? Optionally specify the effective displacement as function of mode mix in tabular form.

? Abaqus assumes that the damage evolution is mode independent otherwise.

Modeling Fracture and Failure with Abaqus

L5.30

Constitutive Response
? Abaqus/CAE interface for displacement-based damage evolution

Modeling Fracture and Failure with Abaqus

L5.31

Constitutive Response
? Energy-based damage evolution

? The fracture energy can be defined as a function of mode mix using either a tabular form or one of two analytical forms:
? Power law
? ? ?

? GI ? ? GII ? ? GIII ? ? ? ?? ? ?? ? =1 ? GIC ? ? GIIC ? ? GIIIC ?
? BK (Benzeggagh-Kenane)
?
For isotropic failure (GIC = GIIC), the response is insensitive to the value of ?.

? Gshear ? GIC ? ? GIIC - GIC ? ? ? = GTC ? GT ? where Gshear = GII ? GIII GT = GI ? Gshear
Modeling Fracture and Failure with Abaqus

L5.32

Constitutive Response
? Keywords interface for energy-based damage evolution
*DAMAGE EVOLUTION, TYPE = ENERGY, SOFTENING = { LINEAR | EXPONENTIAL}, MIXED MODE BEHAVIOR = { TABULAR | POWER LAW | BK }, POWER = value

? Specify fracture energy as function of mode mix in tabular form, or ? Specify the fracture energy in pure normal and shear deformation modes and choose either the POWER LAW or the BK mixed mode behavior

Modeling Fracture and Failure with Abaqus

L5.33

Constitutive Response
? Abaqus/CAE interface for energy-based damage evolution

Modeling Fracture and Failure with Abaqus

L5.34

Constitutive Response
? Example
Normal (opening) mode:
N max
Traction (nominal stress) Cohesive material law: Traction, Damage Evolution

? The preceding discussion was very general in the sense that the full range of options for modeling the constitutive response of cohesive elements was presented. ? In the simplest case, Abaqus requires that you input the adhesive thickness heff and 10 material parameters:
*Elastic, type=traction

Kn =
Kn
1
init ?n

GIC
(area under entire curve)

En heff

Separation

? nfail

En, Et, Es
*Damage initiation, criterion = maxs

What do you do when you only have 1 property and the adhesive thickness is essentially zero?
Diehl, T., "Modeling Surface-Bonded Structures with

Nmax, Tmax, Smax
*Damage evolution, type=energy, ABAQUS Cohesive Elements: Beam-Type Solutions," mixed mode behavior=bk, power=? ABAQUS Users' Conference, Stockholm, 2005. GIC, GIIC , GIIIC
Modeling Fracture and Failure with Abaqus

L5.35

Constitutive Response
? Example (cont’d)

? Common case: you know GTC for the surface bond.
? Assume isotropic behavior

GIC = GIIC = GIIIC = GTC
? For MIXED MODE BEHAVIOR = BK, this makes the response independent of ? term, so set ? = any valid input value (e.g., 1.0)

? Bond thickness is essentially zero
? Specify the cohesive section property thickness heff = 1.0
? Nominal strains = separation; elastic moduli = stiffness

? Isotropy also implies the following:

En = Et = Es = Eeff

(=Keff since we chose heff = 1.0)

Nmax = Tmax = Smax = Tult
Modeling Fracture and Failure with Abaqus

L5.36

Constitutive Response
? Example (cont’d)

? Introduce concept of damage initiation ratio:

?ratio= ?init /?fail, where 0 ? ?ratio ? 1.
? Use GC and equation of a triangle to relate back to Keff and Tult :

Keff =

2 GTC

? ratio ? 2 fail

Tult =

2 GTC

? fail

? The problem now reduces to two penalty terms: ?fail and ?ratio.

? Assume ?ratio = ?.
? Choose ?fail as a fraction of the typical cohesive element mesh size.

? For example, use ?fail = 0.050 ? typical cohesive element size as a starting point.

Modeling Fracture and Failure with Abaqus

L5.37

Constitutive Response
? Example (cont’d)

? Thus, after choosing the two penalty terms, a single (effective) traction-separation law applies to all modes (normal + shear):

Effective properties:
Tult
Traction (nominal stress) Cohesive material law: Traction, Damage Evolution

K eff =
1

Eeff heff

K eff GTC
(area under entire curve)

? init
Separation

? fail

*Cohesive section, thickness=SPECIFIED, ... 1.0, : : *Elastic, type=TRACTION Keff, Keff, Keff *Damage initiation, criterion = MAXS Tult, Tult, Tult *Damage evolution, type=ENERGY, mixed mode behavior=BK, power=1 GTC, GTC , GTC

Modeling Fracture and Failure with Abaqus

L5.38

Constitutive Response
? Example (cont’d)

? What if the response is dynamic? What about the density?
? The density of the cohesive layer should also be considered a penalty quantity. ? For Abaqus/Explicit, the effective density should not adversely affect the stable time increment. Diehl suggests the following rule:

?eff = Eeff

? D tstable ?? ? ft 2 D heff ?

? ? ? ?

2

? Dtstable = stable time increment
without cohesive elements in the model ? ft2D = 0.32213 (for cohesive elements whose original nodal coordinates relate to zero element thickness)

? The Abaqus Analysis User’s Manual provides additional guidelines for determining a cohesive element density that minimizes the effect on the stable time increment in Abaqus/Explicit.
Modeling Fracture and Failure with Abaqus

Viscous Regularization

L5.40

Viscous Regularization
? Cohesive elements have the potential to cause numerical difficulties in the following cases ? Stiff cohesive behavior may lead to reduced maximum stable time increment in Abaqus/Explicit ? Potentially addressed through selective mass scaling ? Unstable crack propagation may lead to convergence difficulties in Abaqus/Standard ? Potentially addressed through built-in viscous regularization option specific to cohesive elements

Modeling Fracture and Failure with Abaqus

L5.41

Viscous Regularization
? Viscous regularization

? Material models with damage often lead to severe convergence difficulties in Abaqus/Standard
? Viscous regularization helps in such cases ? Helps make the consistent tangent stiffness of softening material positive for sufficiently small time increments ? Similar approach used in the concrete damaged plasticity model in Abaqus/Standard

? = ?1 - dv ??
? = 1 ?d - d ? d v v

?

Modeling Fracture and Failure with Abaqus

L5.42

Viscous Regularization
? Consistent material tangent stiffness

D = ?1 - d ? K 0 - f

?d ? ?? ??

K0 is the undamaged elastic stiffness

f is a factor that depends on the details of the damage model
? Viscous regularization ensures that when

Dt

?

? 0 , D = (1 - d ) K 0

? “Offending” second term is eliminated when the analysis cuts back drastically

Modeling Fracture and Failure with Abaqus

L5.43

Viscous Regularization
? User interface for viscous regularization
*COHESIVE SECTION, CONTROLS = control1 *SECTION CONTROLS, NAME = control1, VISCOSITY = factor

? Add-on transverse shear stiffness may provide additional stability *COHESIVE SECTION
*TRANSVERSE SHEAR STIFFNESS ? Output ? Energy associated with viscous regularization: ALLCD

Modeling Fracture and Failure with Abaqus

L5.44

Viscous Regularization
? Example: Multiple delamination problem (Alfano & Crisfield, 2001) – Industry standard AlfanoCrisfield nonsymmetric delamination examples ? Plies are initially bonded with predefined cracks, then peeled apart in a complex sequence
12 layers Initial cracks 2 layers Interface elements

10 layers a
1

? Example done in Abaqus/Standard and Abaqus/Explicit ? Effect of viscous regularization is investigated

a
2

a
2

L

Modeling Fracture and Failure with Abaqus

L5.45

Viscous Regularization

? = 5.e - 4

? = 1.e - 3

?=0
? = 1.e - 4

? = 2.5e - 4

Modeling Fracture and Failure with Abaqus

L5.46

Viscous Regularization
? Effect of viscous regularization on convergence of multiple delamination problem: ? Significant improvements with small regularization factor

Viscous regularization factor

Total number of increments

0.

375

1.0e-4
2.5e-4 1.0e-3

171
153 164

Modeling Fracture and Failure with Abaqus

Modeling Techniques

L5.48

Modeling Techniques
? Model problem: double-cantilever beam

? Alfano and Crisfield (2001)
? Pure Mode I

? Displacement control
? Analyzed using

u

? 1D (B21),
? 2D (CPE4I), and ? 3D (C3D8I) elements

-u
Initial crack

? Delamination assumed to occur along a straight line ? Beams: Orthotropic material

? Cohesive layer: Traction-separation with damage

Modeling Fracture and Failure with Abaqus

L5.49

Modeling Techniques
? One-dimensional model

? Use tie constraints between the cohesive layer and the beams
? Require distinct parts for the beam and cohesive zone geometry

? Geometry

Modeling Fracture and Failure with Abaqus

L5.50

Modeling Techniques
? One-dimensional model (cont’d)

? Assembly

Create 2 instances of the beam; one of the cohesive zone

Position the parts to leave gaps between them; this will later facilitate picking surfaces

Modeling Fracture and Failure with Abaqus

L5.51

Modeling Techniques
? One-dimensional model (cont’d)

? Tie constraints
coh-top beam-top

beam-bot

coh-bot Define tie constraints between mating surfaces.

The cohesive side should be the slave surface (because it is a softer material) This approach is required when quadratic displacement elements are used.

Modeling Fracture and Failure with Abaqus

L5.52

Modeling Techniques
? One-dimensional model (cont’d)

? Properties: beam

Modeling Fracture and Failure with Abaqus

L5.53

Modeling Techniques
? One-dimensional model (cont’d)

? Properties: adhesive

Modeling Fracture and Failure with Abaqus

L5.54

Modeling Techniques
? One-dimensional model (cont’d)

? Meshing

1 Cohesive elements can only
be assigned to sweep meshable regions
Sweep path must be aligned with thickness direction

3 Assign seeds and mesh

Only one element through the thickness

2 Assign cohesive element
type to the swept region

Modeling Fracture and Failure with Abaqus

L5.55

Modeling Techniques
? One-dimensional model (cont’d)

? Meshing (cont’d)

4

Edit the nodal coordinates of each part instance so that they all have the same 2-coordinate

Toggle this off; otherwise, nodes will project back to their original positions

Final mesh

Modeling Fracture and Failure with Abaqus

L5.56

Modeling Techniques
? Two-dimensional model

? All geometry is 2D and planar
? Properties, attributes, etc. treated in a similar manner to the 1D case presented earlier ? Modeling options include: ? Shared nodes

? Tie constraints
? Similar to the 1D model

Modeling Fracture and Failure with Abaqus

L5.57

Modeling Techniques
? Two-dimensional model (cont’d)

? Shared nodes
1

Define a finite thickness slit in the beam as shown below

? Use the actual overall thickness of the DCB ? The center region represents the cohesive layer
2

Mesh the part:

Modeling Fracture and Failure with Abaqus

L5.58

Modeling Techniques
? Two-dimensional model (cont’d)

? Shared nodes (cont’d)
3

Edit the coordinates of the nodes along the interface

Modeling Fracture and Failure with Abaqus

L5.59

Modeling Techniques
? Two-dimensional model (cont’d)

? Tie constraints
1

Create two instances of the beams and position them as shown below.

? Suppress the visibility of the instances to facilitate picking surfaces, etc.
2

Create a finite thickness cohesive layer, position it appropriately in the horizontal direction, define surfaces, etc. ? After meshing, adjust the coordinates of all the nodes in the cohesive layer so that they lie along the interface between the two beams.

Modeling Fracture and Failure with Abaqus

L5.60

Modeling Techniques
? Three-dimensional model

? All geometry is 3D
? Solid geometry for beams

? Solid or shell geometry for cohesive layer
? Modeling options include

? Shared nodes
? Tie constraints

Modeling Fracture and Failure with Abaqus

L5.61

Modeling Techniques
? Three-dimensional model (cont’d)

? Shared nodes

1 Partition the geometry and
define a mesh seam between these two faces

Modeling Fracture and Failure with Abaqus

L5.62

Modeling Techniques
? Three-dimensional model (cont’d)

? Shared nodes (cont’d)

2

Mesh the part with solid (continuum elements).

3

Create a orphan mesh

Mesh→Create Mesh Part

Modeling Fracture and Failure with Abaqus

L5.63

Modeling Techniques
Tip 1: Remove elements from top region with display groups (select by angle)

4 Create a single zero-thickness
solid layer by offsetting from the midplane (selected by angle) of the orphan mesh created in the previous step

Tip 2: Use the selection options tools to facilitate picking. In particular, select from interior entities.

Create a set for the new layer so you can easily assign element type and section properties.
Modeling Fracture and Failure with Abaqus

L5.64

Modeling Techniques
? Three-dimensional model (cont’d)

? Shared nodes (cont’d)

5 Assign section properties and
the element type to the set created in the previous step

Modeling Fracture and Failure with Abaqus

L5.65

Modeling Techniques
? Three-dimensional model (cont’d)

? Tie constraints
? The cohesive region can be defined as ? Solid (with finite thickness)
? Edit nodal coordinates of cohesive elements as in previous examples

? Shell geometry
? Mesh geometry then create orphan mesh ? Offset a zero-thickness layer of solid elements from the orphan mesh
Define surfaces automatically to facilitate tie constraints

Modeling Fracture and Failure with Abaqus

L5.66

Modeling Techniques
? Three-dimensional model (cont’d)

? Tie constraints (cont’d)

When defining the tie constraints, query the mesh stack direction to determine when the “top” and “bottom” surfaces should be used

Brown = top

Purple = bottom

Modeling Fracture and Failure with Abaqus

L5.67

Modeling Techniques
? What if I don't use Abaqus/CAE?

?

In this case do the following in the preprocessor of your choice:
1. Generate the mesh for the structure and cohesive layer (temporarily assigning an arbitrary element type to the cohesive layer) 2. Position the layer of cohesive elements over the interface 3. Define surfaces on the structure and cohesive layer

4. Write the input file

Surface top-beam

Surface top-coh

Surface bot-beam

Surface bot-coh

Modeling Fracture and Failure with Abaqus

L5.68

Modeling Techniques
? Edit the input file: 5. Change the element type assigned to the cohesive layer
*element, elset=coh, type=coh2d4

6. Assign cohesive section properties
*cohesive section, elset=coh, material=cohesive, response=traction separation, stack direction=2, controls=visco 1.0, 0.02 : *material, name=cohesive *elastic, type=traction 5.7e+14, 5.7e+14, 5.7e+14 *damage initiation, criterion=quads 5.7e7, 5.7e7, 5.7e7 *damage evolution, type=energy, mixed mode behavior=bk, power=2.284 280.0, 280.0, 280.0

Modeling Fracture and Failure with Abaqus

L5.69

Modeling Techniques
? The stack direction defines the thickness direction based on the element isoparametric directions. ? ? Set STACK DIRECTION = { 1 | 2 | 3 } to define the element thickness direction along an isoparametric direction. 2D example (extends to 3D):
?2
201 202 201

?1
202

?1
101
102 Thickness direction

?2
101 102

Element connectivity: 101, 102, 202, 201
Stack direction = 2

Element connectivity: 102, 202, 201, 101 Stack direction = 1

Modeling Fracture and Failure with Abaqus

L5.70

Modeling Techniques
? Edit the input file (cont'd): 7. Define tie constraints between the surfaces
Cohesive surface is the slave *tie, name=top, adjust=yes, position tolerance=0.002 top-coh, top-beam *tie, name=bot, adjust=yes, position tolerance=0.002 bot-coh, bot-beam The position tolerance should be large enough to contain the slave nodes when measured from the master surface. In this case the overclosure is equal to 0.001 on either side of the interface so a position tolerance of 0.002 is sufficient to capture all slave nodes. 0.001

Setting adjust=yes will force Abaqus to move the slave (cohesive element) nodes onto the master surface. By adjusting both the top and bottom cohesive surfaces in this way, a zero-thickness cohesive layer is produced.

Modeling Fracture and Failure with Abaqus

L5.71

Modeling Techniques
? Results

Modeling Fracture and Failure with Abaqus

L5.72

Modeling Techniques
? Effect of viscous regularization

Viscous regularization factor

Total number of increments

1.e-5 2.5e-5
5.0e-5 1.0e-4

636 163
129 90

Modeling Fracture and Failure with Abaqus

L5.73

Modeling Techniques
? Effect of mesh refinement

? Typically, you will need to use a much finer mesh (for both the stress/displacement and cohesive elements) than may be necessary for a problem without cohesive elements.

Modeling Fracture and Failure with Abaqus

L5.74

Modeling Techniques
? Non-planar geometry

? The technique for embedding a layer of solid elements into an orphan mesh is not restricted to planar geometry.
? As an example, consider the following fiber-matrix pullout model

matrix
Orphan mesh

fiber

Modeling Fracture and Failure with Abaqus

L5.75

Modeling Techniques
? Failure driven by mismatch in CTEs

View cut of the matrix-fiber interface at 100% of the applied load (magnified 5×)

Failure levels at 38% of the applied load

Modeling Fracture and Failure with Abaqus

L5.76

Modeling Techniques
? Cohesive elements on a symmetry plane

? The traction-separation law is based on N max the separation between the top and bottom faces of the cohesive element.
? On a symmetry plane, however, the separation that is computed is ? the actual value. ? To account for this, specify:

N

GC area = 2
2Kn
1
init ?n

? 2? the cohesive stiffness that would be used in a full model. ? ? the fracture toughness that would be used in a full model. ? Linear equations between the nodes on the top and bottom faces in the lateral directions.

? nfail
2
2 En En = heff heff / 2

?n

2
2Kn =

Modeling Fracture and Failure with Abaqus

L5.77

Modeling Techniques
? Symmetry example
Symmetric model (top) overlaid on full model

Constraint on lateral displacements

Symmetric model

Full model

Constitutive thickness is same as for the full model so double the elastic modulus to double the cohesive stiffness

Modeling Fracture and Failure with Abaqus

Examples

L5.79

Examples
? Composite components in aerospace structures (Courtesy: NASA) ? Stress concentrations around stiffener terminations and flanges

? Residual thermal strains at the interface at room temperature
? Analysis of the effects of residual strains on skin/stiffener debonding ? Delamination initiation and propagation
Beginning of separation After separation

Abaqus/Standard simulation of skin/stiffener debonding

Example Problem 1.4.5

Modeling Fracture and Failure with Abaqus

L5.80

Examples

Abaqus/Standard simulation of skin/stiffener debonding
Modeling Fracture and Failure with Abaqus

L5.81

Examples
? Electronic packaging (Courtesy: INTEL) ? Solder to motherboard fracture due to static overload ? Experiments to assess integrity of solder joints under various loading conditions (e.g., board bending) ? Strain in motherboard at which solder joint fails

Ball grid array
Modeling Fracture and Failure with Abaqus

L5.82

Examples
Debonded solder balls

Damage severity in cohesive layer between motherboard and solder balls
Modeling Fracture and Failure with Abaqus

L5.83

Examples
? Delamination of a composite

? This model is a representative of composite delamination.
? It comprises 3 layers of composite with adhesive layers applied between composite layers. ? The composite delaminates under the impact of a heavy mass displayed in light greenish shade in the animation.

Cohesive layers

Modeling Fracture and Failure with Abaqus

L5.84

Examples
? Impact of moving mass with a stationary wall

? Brick wall modeled with adhesives applied to each face of each brick.
? Simulating damage of the (stationary) wall from high velocity impact with a heavy mass ? Analysis performed in Abaqus/Explicit. ? This model is a representative of several problems that can be modeled using cohesive elements

? Hydroplaning
? Machining ? Oil Drilling ? Excavation

Section of the model illustrating the application of cohesive layers around the bricks.

? Effect of explosion on a building.
Modeling Fracture and Failure with Abaqus

L5.85

Examples
? Deformation sequence

Modeling Fracture and Failure with Abaqus

Workshop 3 (Part 1)

L5.87

Workshop 3 (Part 1)
? Crack growth in a three-point bend specimen using element-based cohesive behavior ? Generate cohesive element mesh ? Define/assign traction-separation behavior and damage properties

Layer of cohesive elements

Modeling Fracture and Failure with Abaqus

Workshop 4 (Optional)

L5.89

Workshop 4 (Optional)
? Crack growth in a helicopter airframe

? Use the mesh offset tool to create a layer of cohesive elements
? Impose symmetry conditions on the cohesive elements using linear equations

Cohesive element thickness shrunk to zero

Modeling Fracture and Failure with Abaqus

Surface-based Cohesive Behavior
Lecture 6

L6.2

Overview
? Surface-based Cohesive Behavior

? Element- vs. Surface-based Cohesive Behavior
? Workshop 3 (Part 2)

Modeling Fracture and Failure with Abaqus

Surface-based Cohesive Behavior

L6.4

Surface-based Cohesive Behavior
? Surface-based cohesive behavior provides a simplified way to model cohesive connections with negligibly small interface thicknesses using the traction-separation constitutive model.

? It can also model “sticky” contact (surfaces can bond after coming into contact). ? The cohesive surface behavior can be defined for general contact in Abaqus/Explicit and contact pairs in Abaqus/Standard (with the exception of the finite-sliding, surface-to-surface formulation).
? Cohesive surface behavior is defined as a surface interaction property. ? To prevent overconstraints in Abaqus/Explicit, a pure master-slave formulation is enforced for surfaces with cohesive behavior.

Modeling Fracture and Failure with Abaqus

L6.5

Surface-based Cohesive Behavior
? User interface

Abaqus/CAE

Abaqus/Standard
*SURFACE INTERACTION, NAME=cohesive *COHESIVE BEHAVIOR ... *CONTACT PAIR, INTERACTION=cohesive surface1, surface2

Abaqus/Explicit
*SURFACE INTERACTION, NAME=cohesive *COHESIVE BEHAVIOR ... *CONTACT *CONTACT PROPERTY ASSIGNMENT surface1, surface2, cohesive

Modeling Fracture and Failure with Abaqus

L6.6

Surface-based Cohesive Behavior
? The formulae and laws that govern surface-based cohesive behavior are very similar to those used for cohesive elements with traction-separation behavior: traction ? linear elastic traction-separation, ? damage initiation criteria, and ? damage evolution laws.
GC
separation

? However, it is important to recognize that damage in surface-based cohesive behavior is an interaction property, not a material property. ? Traction and separation are interpreted differently for cohesive elements and cohesive surfaces:
Cohesive elements
Relative displacement (?) between the top and bottom of the cohesive layer separation Nominal strain (?) = Initial thickness (To) traction Nominal stress (?) Cohesive surfaces

Contact separation (?)
Contact stress (t) =

Contact force (F) Current area (A) at each contact point

Modeling Fracture and Failure with Abaqus

L6.7

Surface-based Cohesive Behavior
? Linear elastic traction-separation behavior

? Relates normal and shear stresses to the normal and shear separations across the interface before the initiation of damage.
? By default, elastic properties are based on underlying element stiffness. ? Can optionally specify the properties.

? Recall this specification is required for cohesive elements.
? The traction-separation behavior can be uncoupled (default) or coupled.

*COHESIVE BEHAVIOR, TYPE= { UNCOUPLED, COUPLED} Optional data line to specify Knn, Kss, Ktt

Modeling Fracture and Failure with Abaqus

L6.8

Surface-based Cohesive Behavior
? Controlling the cohered nodes

? The slave nodes to which cohesive behavior is applied can be controlled to define a wider range of cohesive interactions: Can include:
? All slave nodes ? Only slave nodes initially in contact

? Initially bonded node set
1 ? Applying cohesive behavior to all slave nodes (default)

? Cohesive constraint forces potentially act on all nodes of the slave surface. ? Slave nodes that are not initially contacting the master surface can also experience cohesive forces if they contact the master surface during the analysis.
*COHESIVE BEHAVIOR, ELIGIBILITY = CURRENT CONTACTS

Modeling Fracture and Failure with Abaqus

L6.9

Surface-based Cohesive Behavior
2 Applying cohesive behavior only to slave nodes initially in contact

? Restrict cohesive behavior to only those slave nodes that are in contact with the master surface at the start of a step.
? Any new contact that occurs during the step will not experience cohesive constraint forces. ? Only compressive contact is modeled for new contact.

*COHESIVE BEHAVIOR, ELIGIBILITY = ORIGINAL CONTACTS

Modeling Fracture and Failure with Abaqus

L6.10

Surface-based Cohesive Behavior
3 Applying cohesive behavior only to an initially bonded node set

(Abaqus/Standard only)
? Restrict cohesive behavior to a subset of slave nodes defined using *INITIAL CONDITIONS, TYPE=CONTACT. ? All slave nodes outside of this set will experience only compressive contact forces during the analysis. ? This method is particularly useful for modeling crack propagation along an existing fault line.

*COHESIVE BEHAVIOR, ELIGIBILITY = SPECIFIED CONTACTS

Modeling Fracture and Failure with Abaqus

L6.11

Surface-based Cohesive Behavior
? Example: Double cantilever beam (DCB)

? Analyze debonding of the DCB model using the surface-based cohesive behavior in Abaqus/Standard.
? To model debonding using surface-based cohesive behavior, ? you must define:
1 ? contact pairs and initially bonded crack surfaces;

2 ? the traction-separation behavior;
3 ? the damage initiation criterion; and 4 ? the damage evolution.

? You may also
5 ? specify viscous regularization to facilitate solution convergence u in Abaqus/Standard.

? Note: Steps 3, 4, and 5, will be covered later in this lecture.

-u
Initial crack Cohesive interface

Note: Only the Keywords interface is illustrated in the example; the Abaqus/CAE interface is illustrated in the workshop exercise.
Modeling Fracture and Failure with Abaqus

L6.12

Surface-based Cohesive Behavior
1 ? Define contact pairs and initially bonded crack surfaces

? The initially bonded portion of the slave surface (i.e., node set bond) is identified with the *INITIAL CONDITIONS, TYPE=CONTACT option.
*NSET, NSET=bond, GENERATE 1, 121, 1 *SURFACE, NAME=TopSurf _TopBeam_S1, S1 *SURFACE, NAME=BotSurf _BotBeam_S1, S1 *CONTACT PAIR, INTER=cohesive TopSurf, BotSurf *INITIAL CONDITIONS, TYPE=CONTACT TopSurf, BotSurf, bond slave surface master surface
a list of slave nodes that are initially bonded

bond

TopSurf

BotSurf

Note: Frictionless contact is assumed.

Modeling Fracture and Failure with Abaqus

L6.13

Surface-based Cohesive Behavior
2 ? Define traction-separation behavior t
Kn (Ks , Kt)

? In this model, the cohesive behavior is only enforced for the node set bond.
? Use the ELIGIBILITY=SPECIFIED CONTACTS parameter to enforce this behavior.

1

? Recall the default elastic properties are based on underlying element stiffness. Here we specify the properties.
bond

?
Kn, Ks, and Kt: normal and
tangential stiffness components

TopSurf

BotSurf

... *CONTACT PAIR, INTER=cohesive TopSurf, BotSurf *INITIAL CONDITIONS, TYPE=CONTACT TopSurf, BotSurf, bond *SURFACE INTERACTION, NAME=cohesive *COHESIVE BEHAVIOR, ELIGIBILITY=SPECIFIED CONTACTS 5.7e14, 5.7e14, 5.7e14 Optional

Modeling Fracture and Failure with Abaqus

Kn

Ks

Kt

L6.14

Surface-based Cohesive Behavior
? Damage modeling for cohesive surfaces ? Damage of the traction-separation response for cohesive surfaces is defined within the same general framework used for cohesive elements. ? The difference between the two approaches is that for cohesive surfaces damage is specified as part of the contact interaction properties.
t
max max max tn ts , tt

?

?

max ?n ? smax , ? tmax

?

?

? nf ? sf , ? t f

?

?

?

max max tn , ts , and ttmax :

peak values of the contact stress
max ?n , ? smax , and ?tmax :

peak values of the contact separation

? nf , ? sf , and ? t f :
separations at failure

Modeling Fracture and Failure with Abaqus

L6.15

Surface-based Cohesive Behavior
? User interface Abaqus/Standard
*SURFACE INTERACTION, NAME=cohesive *COHESIVE BEHAVIOR *DAMAGE INITIATION *DAMAGE EVOLUTION *CONTACT PAIR, INTERACTION=cohesive surface1, surface2

Abaqus/CAE

Abaqus/Explicit
*SURFACE INTERACTION, NAME=cohesive *COHESIVE BEHAVIOR *DAMAGE INITIATION *DAMAGE EVOLUTION *CONTACT *CONTACT PROPERTY ASSIGNMENT surface1, surface2, cohesive

Modeling Fracture and Failure with Abaqus

L6.16

Surface-based Cohesive Behavior
? Damage initiation criteria
Maximum stress criterion Maximum separation criterion

? ts tt ? ? tn ? MAX ? max , max , max ? ?1 ts tt ? ? tn ? ?
*DAMAGE INITIATION, CRITERION=MAXS
max max max tn , ts , tt

? ?s ?t ? ? ?n ? MAX ? max , max , max ? ?1 ?s ?t ? ?? n ? ?
max ?n , ? smax , ?tmax

*DAMAGE INITIATION, CRITERION=MAXU

Quadratic stress criterion

Quadratic separation criterion

? tn ? ? ts ? ? tt ? ? max ? ? ? max ? ? ? max ? ? 1 ? tn ? ? ts ? ? tt ?
*DAMAGE INITIATION, CRITERION=QUADS
max max max tn , ts , tt

2

2

2

? ? n ? ? ? s ? ? ?t ? ? max ? ? ? max ? ? ? max ? ? 1 ? ? n ? ? ? s ? ? ?t ?
max ?n , ? smax , ?tmax

2

2

2

*DAMAGE INITIATION, CRITERION=QUADU

tn: normal contact stress in the pure normal mode ts: shear contact stress along the first shear direction tt: shear contact stress along the second shear direction

?n: separation in the pure normal mode ?s: separation in the first shear direction ?t: separation in the second shear direction

Note: Recall the damage initiation criteria for the cohesive elements: if the initial constitutive thickness To = 1, then ? = ?/To = ?. In this case, the separation measures for both approaches are exactly the same. Modeling Fracture and Failure with Abaqus

L6.17

Surface-based Cohesive Behavior
? Example: Double cantilever beam
3 ? Define the damage initiation criterion

? The quadratic stress criterion is specified for this problem.
... *CONTACT PAIR, INTER=cohesive TopSurf, BotSurf *INITIAL CONDITIONS, TYPE=CONTACT TopSurf, BotSurf, bond *SURFACE INTERACTION, NAME=cohesive *COHESIVE BEHAVIOR, ELIGIBILITY=SPECIFIED CONTACTS 5.7e14, 5.7e14, 5.7e14 *DAMAGE INITIATION, CRITERION=QUADS 5.7e7, 5.7e7, 5.7e7
max tn max ts

bond

TopSurf

BotSurf

ttmax

Modeling Fracture and Failure with Abaqus

L6.18

Surface-based Cohesive Behavior
? Damage evolution

? For surface-based cohesive behavior, damage evolution describes the degradation of the cohesive stiffness.
? In contrast, for cohesive elements damage evolution describes the degradation of the material stiffness. ? Damage evolution can be based on energy or separation (same as for cohesive elements). ? Specify either the total fracture energy (a property of the cohesive interaction) or the post damage-initiation effective separation at t failure.

? May depend on mode mix
? Mode mix may be defined in terms of energy or traction

max max max tn ts , tt

?

?
GTC

max ?n ? smax , ? tmax

?

?

? nf ? sf , ? t f

?

?

?

Modeling Fracture and Failure with Abaqus

L6.19

Surface-based Cohesive Behavior
? Separation-based damage evolution

? Damage is a function of an effective separation:

t
max max max tn ts , tt

?

?

Linear postinitiation response

??

?n

2

? ? s2 ? ? t2
max ?n ? smax , ? tmax

? As with cohesive elements, the post damage-initiation softening response can be either:
? Linear ? Exponential ? Tabular

?

?

? nf ? sf , ? t f

?

?

?

Modeling Fracture and Failure with Abaqus

L6.20

Surface-based Cohesive Behavior
? Separation-based damage evolution (cont’d)

? Usage:

*DAMAGE EVOLUTION, TYPE = DISPLACEMENT, SOFTENING = { LINEAR | EXPONENTIAL | TABULAR }, MIXED MODE BEHAVIOR = TABULAR

Modeling Fracture and Failure with Abaqus

L6.21

Surface-based Cohesive Behavior
? Energy-based damage evolution

? As with cohesive elements, the energy-based damage evolution criterion can be defined as a function of mode mix using either a tabular form or one of two analytical forms:
Power law Benzeggagh-Kenane (BK)
?

? GI ? ? GII ? ? GIII ? ? ? ?? ? ?? ? ?1 ? GIC ? ? GIIC ? ? GIIIC ?

?

?

? Gshear ? GIC ? ? GIIC - GIC ? ? ? ? GTC G ? T ? where Gshear ? GII ? GIII GT ? GI ? Gshear

?

Modeling Fracture and Failure with Abaqus

L6.22

Surface-based Cohesive Behavior
? Energy-based damage evolution (cont’d)

? Usage:

*DAMAGE EVOLUTION, TYPE = ENERGY, SOFTENING = { LINEAR | EXPONENTIAL}, MIXED MODE BEHAVIOR = { TABULAR | POWER LAW | BK }, POWER = value

Modeling Fracture and Failure with Abaqus

L6.23

Surface-based Cohesive Behavior
? Example: Double cantilever beam
4 ? Define damage evolution

? The energy-based damage evolution based on the BK mixed mode behavior is specified.
?G ? - GIC ? ? shear ? ? GTC ? GT ?
?

GIC ? ? GIIC

bond

TopSurf

BotSurf

... *CONTACT PAIR, INTER=cohesive TopSurf, BotSurf *INITIAL CONDITIONS, TYPE=CONTACT TopSurf, BotSurf, bond *SURFACE INTERACTION, NAME=cohesive *COHESIVE BEHAVIOR, ELIGIBILITY=SPECIFIED CONTACTS 5.7e14, 5.7e14, 5.7e14 *DAMAGE INITIATION, CRITERION=QUADS 5.7e7, 5.7e7, 5.7e7 *DAMAGE EVOLUTION, TYPE=ENERGY, MIXED MODE BEHAVIOR=BK, POWER=2.284 280.0, 280.0, 280.0

?

GIC

GIIC

GIIIC

Modeling Fracture and Failure with Abaqus

L6.24

Surface-based Cohesive Behavior
? Viscous regularization

? Can be specified to facilitate solution convergence in Abaqus/Standard for surface-based cohesive behavior when stiffness degradation occurs.
? Output: ? Energy associated with viscous regularization: ALLCD

*DAMAGE STABILIZATION

Modeling Fracture and Failure with Abaqus

L6.25

Surface-based Cohesive Behavior
? Example: Double cantilever beam
5 ? Specify a viscosity coefficient for

the cohesive surface behavior

bond

TopSurf

BotSurf

... *CONTACT PAIR, INTER=cohesive TopSurf, BotSurf *INITIAL CONDITIONS, TYPE=CONTACT TopSurf, BotSurf, bond *SURFACE INTERACTION, NAME=cohesive *COHESIVE BEHAVIOR, ELIGIBILITY=SPECIFIED CONTACTS 5.7e14, 5.7e14, 5.7e14 *DAMAGE INITIATION, CRITERION=QUADS 5.7e7, 5.7e7, 5.7e7 *DAMAGE EVOLUTION, TYPE=ENERGY, MIXED MODE BEHAVIOR=BK, POWER=2.284 280.0, 280.0, 280.0 *DAMAGE STABILIZATION 1.e-5

viscosity coefficient,

?

Modeling Fracture and Failure with Abaqus

L6.26

Surface-based Cohesive Behavior
? Example: Double cantilever beam

? Summary of the input for the traction-separation response
Cohesive elements
*COHESIVE SECTION, MATERIAL=cohesive, RESPONSE=TRACTION SEPARATION, ELSET=coh_elems, CONTROLS=visco , 0.02 *MATERIAL, NAME=cohesive *ELASTIC, TYPE=TRACTION 5.7e14, 5.7e14, 5.7e14 *DAMAGE INITIATION, CRITERION=QUADS 5.7e7, 5.7e7, 5.7e7 *DAMAGE EVOLUTION, TYPE=ENERGY, MIXED MODE BEHAVIOR=BK, POWER=2.284 280.0, 280.0, 280.0 *SECTION CONTROLS, NAME=visco, VISCOSITY=1.e-5

Cohesive surfaces
*SURFACE INTERACTION, NAME=cohesive *COHESIVE BEHAVIOR, ELIGIBILITY=SPECIFIED CONTACTS 5.7e14, 5.7e14, 5.7e14 *DAMAGE INITIATION, CRITERION=QUADS 5.7e7, 5.7e7, 5.7e7 *DAMAGE EVOLUTION, TYPE=ENERGY, MIXED MODE BEHAVIOR=BK, POWER=2.284 280.0, 280.0, 280.0 *DAMAGE STABILIZATION 1.e-5

Modeling Fracture and Failure with Abaqus

L6.27

Surface-based Cohesive Behavior
? Results

u2 = 0.006

Cohesive elements

Failed cohesive elements

u2
u2 = 0.006

Cohesive surfaces

u2

Modeling Fracture and Failure with Abaqus

Element- vs. Surface-based Cohesive Behavior

L6.29

Element- vs. Surface-based Cohesive Behavior
Preprocessing

? Cohesive elements
? Gives you direct control over the cohesive element mesh density and stiffness properties. ? Constraints are enforced at the element integration points.
? Refining the cohesive elements relative to the connected structures will likely lead to improved constraint satisfaction and more accurate results.

? Cohesive surfaces

Integration points on an 8-node cohesive element

? Are easily defined using contact interactions and cohesive interaction properties. ? A pure master-slave in formulation is used.

? Constraints are enforced at the slave nodes.
? Refining the slave surface relative to the master surface will likely lead to improved constraint satisfaction and more accurate results .

Modeling Fracture and Failure with Abaqus

L6.30

Element- vs. Surface-based Cohesive Behavior
Initial configuration:

? Cohesive elements
? Must be bonded at the start of the analysis.

? Once the interface has failed, the surfaces do not re-bond.
? Cohesive surfaces

? Can bond anytime contact is established (i.e., “sticky” contact behavior). ? Cohesive interface need not be bonded at the start of the analysis. ? You can control whether debonded surfaces will stick or not stick if contact occurs again. ? By default, they do not stick.

Modeling Fracture and Failure with Abaqus

L6.31

Element- vs. Surface-based Cohesive Behavior
Constitutive behavior:

? Cohesive elements
? Allow for several constitutive behavior types:

? Traction-separation constitutive model
? Including multiple failure mechanisms

? Continuum-based constitutive model
? For adhesive layers with finite thickness ? Uses conventional material models

? Uniaxial stress-based constitutive model
? Useful in modeling gaskets and/or single adhesive patches

? Cohesive surfaces ? Must use the traction-separation interface behavior. ? Intended for bonded interfaces where the interface thickness is negligibly small.

? Only one failure mechanism is allowed.
Modeling Fracture and Failure with Abaqus

L6.32

Element- vs. Surface-based Cohesive Behavior
Influence on stable time increment (Abaqus/Explicit only):
? Le ? ?t ? ? ? c d ? ?

? Cohesive elements
? Often require a small stable time increment.

? Cohesive elements are generally thin and sometimes quite stiff.
? Consequently, they often have a stable time increment that is significantly less than that of the other elements in the model. ? Cohesive surfaces ? Cohesive surface behavior with the default cohesive stiffness properties is formulated to minimally affect the stable time increment. ? Abaqus uses default contact penalties to model the cohesive stiffness behavior in this case. ? You can specify a non-default cohesive stiffness values. ? However, high stiffnesses may reduce the stable time increment.

Modeling Fracture and Failure with Abaqus

L6.33

Element- vs. Surface-based Cohesive Behavior
Mass:

? Cohesive elements
? The element material definitions include mass.

? Cohesive surfaces
? Do not add mass to the model.

? Indented for thin adhesive interfaces; thus, neglecting adhesive mass is appropriate for most applications. ? However, nonstructural mass can be added to the contacting elements if necessary.

Modeling Fracture and Failure with Abaqus

L6.34

Element- vs. Surface-based Cohesive Behavior
Summary:

? Cohesive elements
? Are recommended for more detailed adhesive connection modeling.

? Additional preprocessing effort (and often increased computational cost) is compensated for by gaining:
? Direct control over the connection mesh ? Additional constitutive response options
? E.g., model adhesives of finite thickness

? Cohesive surfaces ? Provides a quick and easy way to model adhesive connections. ? Negligible interface thicknesses only ? Surfaces can bond anytime contact is established (“sticky” contact)
? Model contact adhesives, Velcro, tape, and other bonding agents that can stick after separation.

Modeling Fracture and Failure with Abaqus

Workshop 3 (Part 2)

L6.36

Workshop 3 (Part 2)
? Crack growth in a three-point bend specimen using surface-based cohesive behavior ? Repeat the element-based exercise using surface-based behavior ? Use default traction-separation elastic properties ? Compare with element-based results

Modeling Fracture and Failure with Abaqus

Virtual Crack Closure Technique (VCCT)
Lecture 7

L7.2

Overview
? Introduction

? VCCT Criterion
? Output

? VCCT Plug-in
? Comparison with Cohesive Behavior

? Examples
? Workshop 5

Modeling Fracture and Failure with Abaqus

Introduction

L7.4

Introduction
? Motivation is aircraft composite structural analysis ? To reduce the cost of laminated composite structures, large integrated bonded structures are being considered.

? In primary structures, bondlines and interfaces between plies are required to carry interlaminar loads.
? Damage tolerance requirements dictate that bondlines and interfaces carry required loads with damage.
Modeling debonding along skin-stringer interface

Modeling Fracture and Failure with Abaqus

L7.5

Introduction
? Analysis requirements for composite damage

? Apply Linear Elastic Fracture Mechanics (LEFM) to bondlines and interfaces
? 2D and 3D delaminations ? Propagation

? Mode separation
? Multiple cracks

? Non-linear behavior (e.g., postbuckling)
? Composite structure

? Practical (CPU time, minimum set of models)

Modeling Fracture and Failure with Abaqus

L7.6

Introduction
? VCCT uses LEFM concepts Pure Mode I
Modified VCCT

? Based on computing the energy release rates for normal and shear crack-tip deformation modes.
? Compare energy release rates to interlaminar fracture toughness.
? See Rybicki, E. F., and Kanninen, M. F., "A Finite Element Calculation of Stress Intensity Factors by a Modified Crack Closure Integral," Engineering Fracture Mechanics, Vol. 9, pp. 931-938, 1977.

Node numbers are shown

Nodes 2 and 5 will start to release when: 1 v1,6 Fv,2,5 ? GI ? GIC Mode II treated 2 bd similarly where GI ? mode I energy release rate GIC ? critical mode I energy release rate b ? width Fv ,2,5 ? vertical force between nodes 2 and 5 v1,6 ? vertical displacement between nodes 1 and 6

Modeling Fracture and Failure with Abaqus

VCCT Criterion

L7.8

VCCT Criterion
? The debond capability is used to perform the crack propagation analysis for initially bonded crack surfaces. ? The crack propagation analysis allows for five types of fracture criteria :
1 ? Critical stress criterion

2 ? Crack opening displacement criterion
3 ? Crack length vs. time criterion 4 ? VCCT criterion 5 ? Low-cycle fatigue criterion

? Defining case 4, “VCCT criterion,” is the subject of this lecture.

? The details of cases 1, 2, and 3 are not discussed here. Please consult the Abaqus Analysis User’s Manual for more details. ? The details of case 5 will be discussed later in Lecture 8 “Low -cycle Fatigue.”

Modeling Fracture and Failure with Abaqus

L7.9

VCCT Criterion
? When using VCCT to model crack propagation,

? you must:
1 ? define contact pairs for potential crack surfaces;
2 ? define initially bonded crack surfaces; 3 ? activate the crack propagation capability; and 4 ? specify the VCCT criterion.

? you also may: ? define spatially varying critical energy release rates; ? use viscous regularization, contact stabilization, and/or automatic stabilization to overcome convergence difficulties for unstable propagating cracks;

? use a linear scaling technique to accelerate convergence for VCCT.

Modeling Fracture and Failure with Abaqus

L7.10

VCCT Criterion
? Defining the VCCT criterion is not currently supported in Abaqus/CAE.

? However, the VCCT plug-in is available and allows you to interactively define the debond interface(s).
? The details of the VCCT plug-in will be discussed later in this lecture. ? Downloaded from “VCCT plug-in utility,” SIMULIA Answer 3235.

Modeling Fracture and Failure with Abaqus

L7.11

VCCT Criterion
? Example: Double cantilever beam (DCB)

? Analyze debonding of a DCB model using the VCCT criterion.
? Steps required for setting up the model include:
? Define slave (TopSurf) and master (BotSurf) surfaces along the debond interface.
? Define a set (bond) containing the initially bonded region (part of TopSurf in this example).

? The Keywords interface is illustrated in this example.
bond

TopSurf BotSurf

Modeling Fracture and Failure with Abaqus

L7.12

VCCT Criterion
1 ? Define contact pairs for potential crack surfaces

? Potential crack surfaces are modeled as slave and master contact surfaces.
? Any contact formulation except the finite-sliding, surface-to-surface formulation can be used. ? Cannot be used with self-contact.
bond
*NSET, NSET=bond, GENERATE 1, 121, 1 *SURFACE, NAME=TopSurf _TopBeam_S1, S1 *SURFACE, NAME=BotSurf _BotBeam_S1, S1 *CONTACT PAIR, INTER=... TopSurf, BotSurf slave surface
Modeling Fracture and Failure with Abaqus

TopSurf

BotSurf

Note: The frictionless interaction property is assumed.
master surface

L7.13

VCCT Criterion
2 ? Define initially bonded crack surfaces

? The initially bonded contact pair is identified with the *INITIAL CONDITIONS, TYPE=CONTACT option.

bond

TopSurf

BotSurf

*NSET, NSET=bond, GENERATE 1, 121, 1 *SURFACE, NAME=TopSurf _TopBeam_S1, S1 *SURFACE, NAME=BotSurf _BotBeam_S1, S1 *CONTACT PAIR, INTER=... TopSurf, BotSurf *INITIAL CONDITIONS, TYPE=CONTACT TopSurf, BotSurf, bond slave surface master surface

a list of slave nodes that are initially bonded

Modeling Fracture and Failure with Abaqus

L7.14

VCCT Criterion
? The unbonded portion of the slave surface will behave as a regular contact surface. ? If the node set that includes the initially bonded slave nodes is not specified, the initial contact condition will apply to the entire contact pair. ? In this case, no crack tips can be identified, and the bonded surfaces cannot separate. ? For the VCCT criterion, the initially bonded nodes are bonded in all directions.

Modeling Fracture and Failure with Abaqus

L7.15

VCCT Criterion
3 ? Activate the crack propagation capability

? The ?DEBOND option is used to activate crack propagation in a given step.
? The SLAVE and MASTER parameters identify the surfaces to be debonded. *NSET, NSET=bond, GENERATE
1, 121, 1 *SURFACE, NAME=TopSurf _TopBeam_S1, S1 *SURFACE, NAME=BotSurf _BotBeam_S1, S1 *CONTACT PAIR, INTER=... TopSurf, BotSurf *INITIAL CONDITIONS, TYPE=CONTACT TopSurf, BotSurf, bond *STEP, NLGEOM *STATIC ... *DEBOND, SLAVE=TopSurf, MASTER=BotSurf

bond

TopSurf

BotSurf

Modeling Fracture and Failure with Abaqus

L7.16

VCCT Criterion
4 ? Specify the VCCT criterion

? The BK law model is used in this example.
BK law:

GequivC ? GIC ? ? GIIC

? GII ? GIII ? ? GIC ? ? ? ? GI ? GII ? GIII ?

?

bond

TopSurf

BotSurf

*NSET, NSET=bond, GENERATE 1, 121, 1 *SURFACE, NAME=TopSurf _TopBeam_S1, S1 *SURFACE, NAME=BotSurf _BotBeam_S1, S1 *CONTACT PAIR, INTER=... TopSurf, BotSurf *INITIAL CONDITIONS, TYPE=CONTACT TopSurf, BotSurf, bond *STEP, NLGEOM *STATIC ... *DEBOND, SLAVE=TopSurf, MASTER=BotSurf *FRACTURE CRITERION, TYPE=VCCT, MIXED MODE BEHAVOIR=BK 280.0, 280.0, 0.0, 2.284

GIC

GIIC

GIIIC

?

Modeling Fracture and Failure with Abaqus

L7.17

VCCT Criterion
? The crack-tip node debonds when the fracture criterion, f,

f ?

Gequiv GequivC

,

reaches the value 1.0 within a given tolerance, ftol:

1 ? f ? 1 ? ftol .
where

Gequiv is the equivalent strain energy release rate, and GequivC is the critical equivalent strain energy release rate calculated
based on the user-specified mode-mix criterion and the bond strength of the interface.
? For the VCCT criterion, the default value of ftol is 0.2. ? Use following option to control ftol:
*FRACTURE CRITERION, TYPE=VCCT, TOLERANCE=ftol
Modeling Fracture and Failure with Abaqus

L7.18

VCCT Criterion
? In the DCB model, the tolerance is set to 0.1.
*NSET, NSET=bond, GENERATE 1, 121, 1 *SURFACE, NAME=TopSurf _TopBeam_S1, S1 *SURFACE, NAME=BotSurf _BotBeam_S1, S1 *CONTACT PAIR, INTER=... TopSurf, BotSurf *INITIAL CONDITIONS, TYPE=CONTACT TopSurf, BotSurf, bond *STEP, NLGEOM *STATIC ... *DEBOND, SLAVE=TopSurf, MASTER=BotSurf *FRACTURE CRITERION, TYPE=VCCT, MIXED MODE BEHAVOIR=BK, TOLERANCE=0.1 280.0, 280.0, 0.0, 2.284

bond

TopSurf

BotSurf

Modeling Fracture and Failure with Abaqus

L7.19

VCCT Criterion
? In addition to the BK law model, Abaqus/Standard also provides two other commonly used mode-mix criteria for computing GequivC: the Power law and the Reeder law models. ? An appropriate model is best selected empirically. ? Power law
Gequiv GequivC ? G ? ?? I ? ? GIC ?
am

? G ? ? ? II ? ? GIIC ?

an

? G ? ? ? III ? ? GIIIC ?

ao

*FRACTURE CRITERION, TYPE=VCCT, MIXED MODE BEHAVIOR=POWER GIC, GIIC, GIIIC, am, an, ao

? Reeder law ? Applies only to three-dimensional problems
? ? GIII GequivC ? GIC ? ? GIIC ? GIC ? ? GIIIC ? GIIC ? ? ? GII ? GIII ? ? ? ? GII ? GIII ? ?? Gi ? ?? ? ? ? ? ?
?

?

*FRACTURE CRITERION, TYPE=VCCT, MIXED MODE BEHAVIOR=REEDER GIC, GIIC, GIIIC, ?
Modeling Fracture and Failure with Abaqus

L7.20

VCCT Criterion
? Spatially varying critical energy release rates

? The VCCT criterion can be defined with varying energy release rates by specifying the critical energy release rates at all nodes on the slave surface.
? In this case, the critical energy release rates should be interpolated from the critical energy release rates specified at the nodes with the *NODAL ENERGY RATE option. ? However, the exponents (e.g., ?) are still read from the data lines under the *FRACTURE CRITERION option.
*NODAL ENERGY RATE node ID1, GIC, GIIC, GIIIC model data node ID2, GIC, GIIC, GIIIC ... *STEP *STATIC ... *FRACTURE CRITERION, TYPE=VCCT, MIXED MODE BEHAVIOR=BK, NODAL ENERGY RATE GIC, GIIC, GIIIC, ?
Modeling Fracture and Failure with Abaqus

L7.21

VCCT Criterion
? Viscous regularization for VCCT ? Can be used to overcome some convergence difficulties for unstable propagating cracks. ? Example: DCB
*NSET, NSET=bond, GENERATE 1, 121, 1 *SURFACE, NAME=TopSurf _TopBeam_S1, S1 *SURFACE, NAME=BotSurf _BotBeam_S1, S1 *CONTACT PAIR, INTER=... TopSurf, BotSurf *INITIAL CONDITIONS, TYPE=CONTACT TopSurf, BotSurf, bond *STEP, NLGEOM *STATIC ... *DEBOND, SLAVE=TopSurf, MASTER=BotSurf, VISCOSITY=0.1 *FRACTURE CRITERION, TYPE=VCCT, MIXED MODE BEHAVOIR=BK, TOLERANCE=0.1 280.0, 280.0, 0.0, 2.284

? Set the value of the viscosity coefficient to 0.1.
bond

TopSurf

BotSurf

Modeling Fracture and Failure with Abaqus

L7.22

VCCT Criterion
? In addition, contact and automatic stabilization that are not specific to VCCT can be also used to aid convergence. ? They are built into Abaqus/Standard and are compatible with VCCT. ? Note that the crack propagation behavior may be modified by the damping forces. ? Therefore, monitor the damping energy (ALLVD or ALLSD) and compare it with the total strain energy in the model (ALLSE) to ensure that the results are reasonable in the presence of damping.

? ALLVD stores the damping energy generated from viscous regularization. ? ALLSD stores the damping energy generated from contact stabilization and automatic stabilization.

Modeling Fracture and Failure with Abaqus

L7.23

VCCT Criterion
? Linear scaling to accelerate convergence for VCCT

? Abaqus provides a linear scaling technique to quickly converge to the critical load state. This reduces the solution time required to reach the onset of crack growth.
? This technique works best for models in which the deformation is nearly linear before the onset of crack growth. ? Once the first crack-tip node releases, the linear scaling calculations will no longer be valid and the time increment will be set to the default value.

? Usage:
*CONTROLS, LINEAR SCALING

?

where ? is the coefficient of linear scaling.

? For details of linear scaling to accelerate convergence for VCCT, see “Crack propagation analysis,” Section 11.4.3 of the Abaqus Analysis User’s Manual.
Modeling Fracture and Failure with Abaqus

L7.24

VCCT Criterion
? Tips for using the VCCT criterion

? Crack propagation problems using the VCCT criterion are numerically challenging.
? To help you create a successful model, several tips for using the VCCT criterion are provided: ? The master debonding surfaces must be continuous. ? The tie MPCs should NOT be used for the slave debonding surface to avoid overconstraints. ? A small clearance between the debonding surfaces can be specified to eliminate unnecessary severe discontinuity iterations during incrementation as the crack begins to progress. …… ? Note: More tips are provided in “Crack propagation analysis,” Section 11.4.3 of the Abaqus Analysis User’s Manual.

Modeling Fracture and Failure with Abaqus

Output

L7.26

Output
? The following output options are provided to support the VCCT criterion:

? Abaqus/CAE supports the surface output requests for VCCT.

*OUTPUT, FIELD, FREQUENCY=freq *CONTACT OUTPUT, MASTER=master, SLAVE=slave *OUTPUT, HISTORY, FREQUENCY=freq *CONTACT OUTPUT, [(MASTER=master, SLAVE=slave)|(NSET=nset)]

Modeling Fracture and Failure with Abaqus

L7.27

Output
? The following bond failure quantities can be requested as surface output:

DBT
DBSF

The time when bond failure occurred
Fraction of stress at bond failure that still remains

DBS
OPENBC

Stress in the failed bond that remains
Relative displacement behind crack.

CRSTS
ENRRT EFENRRTR BDSTAT

Critical stress at failure.
Strain energy release rate. Effective energy release rate ratio. Bond state (=1.0 if bonded, 0.0 if unbonded)

? All of the above variables can be visualized in Abaqus/Viewer.

? The initial contact status of all of the slave nodes is printed in the data (.dat) file.

Modeling Fracture and Failure with Abaqus

L7.28

Output
? Example: DCB

? Request surface output:

bond

... *INITIAL CONDITIONS, TYPE=CONTACT TopSurf TopSurf, BotSurf, bond BotSurf *STEP, NLGEOM *STATIC ... *DEBOND, SLAVE=TopSurf, MASTER=BotSurf, VISCOSITY=0.1 *FRACTURE CRITERION, TYPE=VCCT, MIXED MODE BEHAVOIR=BK, TOLERANCE=0.1 280, 280, 280, 2.284 ... *OUTPUT, FIELD, VAR=PRESELECT field output *CONTACT OUTPUT, SLAVE=TopSurf, MASTER=BotSurf DBT, DBS, OPENBC, CRSTS, ENRRT, BDSTAT *OUTPUT, HISTORY history output *CONTACT OUTPUT, SLAVE=TopSurf, MASTER=BotSurf, NSET=bond DBT, DBS, OPENBC, CRSTS, ENRRT, BDSTAT *NODE OUTPUT, NSET=tip U2, RF2 *END STEP
Modeling Fracture and Failure with Abaqus

L7.29

Output
? Results

VCCT

Modeling Fracture and Failure with Abaqus

VCCT Plug-in

L7.31

VCCT Plug-in
? VCCT plug-in

? provides an interactive interface to define the debond interface(s).
? supports the following keyword options required for VCCT analysis:
*INITIAL CONDITIONS, TYPE=CONTACT *DEBOND, SLAVE=slave, MASTER=master, OUTPUT=[fil|dat|both], VISCOSITY= ? *FRACTURE CRITERION, TYPE=VCCT, MIXED MODE BEHAVIOR=[BK|POWER|REEDER], TOLERANCE= ftol, NODAL ENERGY RATE *NODAL ENERGY RATE *CONTROLS, LINEAR SCALING

? For details please refer to “VCCT plug-in utility,” SIMULIA Answer 3235.

Modeling Fracture and Failure with Abaqus

L7.32

VCCT Plug-in
? Example: Double Cantilever Beam (DCB)

? The VCCT plug-in is discussed in the context of the Keywords interface presented earlier.

bond

TopSurf

BotSurf

initially bonded region master surface slave surface

Modeling Fracture and Failure with Abaqus

L7.33

VCCT Plug-in
1 ? Define contact pairs for potential crack surfaces

? Frictionless contact is assumed.

*SURFACE INTERACTION, NAME=IntProp-1 1. *FRICTION 0.0 *CONTACT PAIR, INTERACTION=IntProp-1 TopSurf, BotSurf

bond

TopSurf

BotSurf

Modeling Fracture and Failure with Abaqus

L7.34

VCCT Plug-in
2 ? Define the VCCT criterion 2a ? Select the fracture criterion, viscosity

coefficient, and cutback tolerance.
... *STEP, NLGEOM *STATIC ... *DEBOND, SLAVE=TopSurf, MASTER=BotSurf, VICOSITY=0.1 *FRACTURE CRITERION, TYPE=VCCT, TOLERANCE=0.2, MIXED MODE BEHAVOIR=BK 280, 280, 280, 2.284
bond

TopSurf

BotSurf

Modeling Fracture and Failure with Abaqus

L7.35

VCCT Plug-in
2b ? Specify critical strain energy release rates

... *STEP, NLGEOM *STATIC ... *DEBOND, SLAVE=TopSurf, MASTER=BotSurf, VICOSITY=0.1 *FRACTURE CRITERION, TYPE=VCCT, TOLERANCE=0.2, MIXED MODE BEHAVOIR=BK 280, 280, 280, 2.284
bond

TopSurf

BotSurf

Modeling Fracture and Failure with Abaqus

L7.36

VCCT Plug-in
? The VCCT plug-in also supports defining spatially varying critical energy release rates. ? Click mouse button 3 to manage the table.

*NODAL ENERGY RATE node ID1, GIC, GIIC, GIIIC node ID2, GIC, GIIC, GIIIC ... *STEP *STATIC ... *FRACTURE CRITERION, TYPE=VCCT, MIXED MODE BEHAVIOR=BK, NODAL ENERGY RATE GIC, GIIC, GIIIC, ?

Modeling Fracture and Failure with Abaqus

L7.37

VCCT Plug-in
3 ? Define the VCCT bonded interface

? Select the initially bonded region, the crack propagation output file and frequency, and the debond initiation step.

? Note: The VCCT plug-in allows specification of linear scaling.
*INITIAL CONDITIONS, TYPE=CONTACT TopSurf, BotSurf, bond *STEP, NAME=Step-1 *STATIC, NLGEOM ... *DEBOND, SLAVE=TopSurf, MASTER=BotSurf, VISCOSITY=0.1 *FRACTURE CRITERION, TYPE=VCCT, MIXED MODE BEHAVOIR=BK 280, 280, 280, 2.284

Modeling Fracture and Failure with Abaqus

L7.38

VCCT Plug-in
? The relevant keywords will be generated when Abaqus/CAE writes the input file.

surface interaction initial contact conditions

debond

fracture criterion

field output

history output
Modeling Fracture and Failure with Abaqus

Comparison with Cohesive Behavior

L7.40

Comparison with Cohesive Behavior
? VCCT and cohesive behavior are very similar in their application and formulation. ? Both theories ? are used to model interfacial shearing and delamination crack propagation and failure, ? use an elastic damage constitutive theory to model the material's response once damage has initiated, and ? dissipate the same amount of fracture energy between damage initiation and complete failure.

Modeling Fracture and Failure with Abaqus

L7.41

Comparison with Cohesive Behavior
? The fundamental difference between VCCT and cohesive behavior is in the way crack propagation is predicted. ? In VCCT an existing flaw is assumed. ? VCCT is appropriate for brittle crack propagation problems. ? However, cohesive behavior can model damage initiation.

? Damage initiation in cohesive behavior is based strictly on the predefined ultimate (normal and/or shear) stress/strain limit. ? Cohesive behavior can be used for both brittle and ductile crack propagation problems.

Modeling Fracture and Failure with Abaqus

L7.42

Comparison with Cohesive Behavior
? VCCT may be viewed as more fundamentally based on fracture mechanics. ? The damage initiation and damage evolution are both based on fracture energy, whereas cohesive behavior use the fracture energy only during damage evolution. ? Applicability of VCCT is limited to “self-similar” crack propagation analyses. ? This implies a steady-state running crack.

? Difficult to reproduce in practice.

Modeling Fracture and Failure with Abaqus

L7.43

Comparison with Cohesive Behavior
? Summary: Complementary techniques for modeling of debonding
VCCT
Use the debond framework (surface based) Assumes an existing flaw Brittle fracture using LEFM occurring along a well defined crack front Requires GI, GII, and GIII Crack propagates when strain energy release rate exceeds fracture toughness Crack surfaces are rigidly bonded when uncracked.

Cohesive behavior
Interface elements (element based) or contact (surface based) Can model crack initiation Ductile fracture occurring over a smeared crack front modeled with spanning cohesive elements or cohesive contact Requires E, σmax, GI, GII, and GIII Crack initiates when cohesive traction exceeds critical value and releases critical strain energy when fully open Crack surfaces are joined elastically when uncracked.

Available only in Abaqus/Standard

Available in Abaqus/Standard and Abaqus/Explicit

? Both are needed to satisfy general fracture requirements
Modeling Fracture and Failure with Abaqus

Examples

L7.45

Examples
? Verification problems

? DCB
? SLB

? ENF
? Alfano-Crisfield

? Alfano, G., and M. A. Crisfield, “Finite Element Interface Models for the Delamination Analysis of Laminated Composites: Mechanical and Computational Issues,” International Journal for Numerical Methods in Engineering, vol. 50, pp. 1701–1736, 2001.
? Also available as Abaqus Benchmark Problem 2.7.1 with cohesive elements ? NASA Panel ? Reeder, J.R., Song, K., Chunchu, P.B., and Ambur, D.R., “Postbuckling and Growth of Delaminations in Composite Plates Subjected to Axial Compression,” AIAA 2002 -1746.
Modeling Fracture and Failure with Abaqus

L7.46

Examples
? Compression Buckling/Delamination Single Disbond (Unreinforced)
Multiple crack tips Buckling driven delaminations

30000 Euler buckling 25000 20000

Load (lb)

15000 10000 5000 0 0 0.01 0.02 0.03 0.04 0.05 Displacement (in)

FEA closed form

Modeling Fracture and Failure with Abaqus

L7.47

Examples

Modeling Fracture and Failure with Abaqus

L7.48

Examples
? Compression Buckling/Delamination Multiple Disbonds (Unreinforced)

Multiple cracks can also be addressed

Modeling Fracture and Failure with Abaqus

L7.49

Examples

Modeling Fracture and Failure with Abaqus

L7.50

Examples
? T-Joint Pull–off Model

Modeling Fracture and Failure with Abaqus

L7.51

Examples
? Postbuckling Behavior of Skin-Stringer Panels

? VCCT can be applied to determine the global strength and failure mode for typical aerospace composite structures like this skin/stringer panel

Courtesy Boeing

Modeling Fracture and Failure with Abaqus

L7.52

Examples

Displacement imposed at corner nodes

Contact surfaces defined for region of fracture

Modeling Fracture and Failure with Abaqus

L7.53

Examples

Crack tip

Initially debonded nodes Initially bonded nodes

Modeling Fracture and Failure with Abaqus

L7.54

Examples

The Abaqus Tech Brief on skin/stringer bonded joint analysis can be downloaded from www.simulia.com

Modeling Fracture and Failure with Abaqus

L7.55

Examples

Modeling Fracture and Failure with Abaqus

Workshop 5

L7.57

Workshop 5
? Crack growth in a three-point bend specimen using VCCT

? Repeat the cohesive-based exercises using VCCT and compare results

Modeling Fracture and Failure with Abaqus

Low-cycle Fatigue
Lecture 8

L8.2

Overview
? Introduction

? Low-cycle Fatigue in Bulk Materials
? Low-cycle Fatigue at Material Interfaces

Modeling Fracture and Failure with Abaqus

Introduction

L8.4

Introduction
? Low-cycle fatigue analysis is a quasi-static analysis of a structure subjected to sub-critical cyclic loading. ? It can be associated with thermal as well as mechanical loading. ? In Abaqus can simulate low-cycle fatigue in: ? bulk ductile materials

? material interfaces

Modeling Fracture and F

相关文章:
Abquas中vcct与debond对比
Application of the She... 10页 免费 VCCT分层损伤分析过程 12页 7下载券 ...Abaqus裂纹扩展分析 5页 免费 XFEM分析 17页 免费 Fracture Mechanics wit... ...
abaqus结构分析单元类型
of numbers,enjoy~ ;updated in 20130116, based on the wordfile "abaqus_67...*FRACTURE *FRAME *FREQUENCY *FRICTION *GAP *GASKET *GEL *GEOSTATIC *GLOBAL...
abaqus XFEM 方法介绍
abaqus XFEM 方法介绍_机械/仪表_工程科技_专业资料...否 则就提示某些单元的 fracture damage 模型未能...the compressive behavior of the crack surfaces has...
abaqus材料子程序
abaqus材料子程序_材料科学_工程科技_专业资料。各向...(2)) C Critical values of stresses XT=PROPS(3...Fracture energy in matrix ETA=0.001 C Current ...
abaqus帮助文档_shear damage
Section 23.2.2 of the Abaqus Analysis User's Manual. To define shear ...5. Enter damage parameters in the Data table: Fracture Strain Equivalent ...
cohesive模拟例子
abaqus 的关于断裂那部分培训资料里面有一些例子,附着很详细的建立过程,挺不错的...选中 strains 中的 SDEG 参数、 failure\fracture 中的 DMICRT 参数、 state ...
edu_ecologychuanke1477656347
09 lecture8-动态分析2 12分钟 10 lecture9-动态接触 20分钟 11 lecture10 33分钟 12 lecture11 17分钟 第5章 abaqus初级教学视频01...
ABAQUS 裂纹
IS REACHED CHECK ENERGY FRACTURE CRITERION , C---...OF MAIN CODE SUBROUTINES 感觉掌握好以上这些方法这些对 ABAQUS 的断裂分析应该算是比较完善了,个人见解 ...
求救
Simulation Module Indicates End of Simulation *** INFO: End of Automatic Remeshing...无法和 ABAQUS 相比,尤其是在考虑模具断裂的时候,无法进行 fracture mechanic...
关于断裂能量
fracture energy,ABAQUS 有自带的 linear,exponential sofening law, 可以定义其他...然后选第三个选项 do not allow the number of elements to change 不准 seed...
更多相关标签:
fracture of teeth | fracture | fracture mechanics | fracture toughness | engineering fracture | fracturefx 如何安装 | fracturefx | brittle fracture |