A Course On Damage Mechanics-Springer-Verlag Berlin Heidelberg (1996) Professor Jean Lemaitre (Auth.)
A Course On Damage Mechanics-Springer-Verlag Berlin Heidelberg (1996) Professor Jean Lemaitre (Auth.)
A Course On Damage Mechanics-Springer-Verlag Berlin Heidelberg (1996) Professor Jean Lemaitre (Auth.)
A Course on
Damage Mechanics
With a Foreword by
Professor H. Lippmann
Springer
Professor Jean Lemaitre
Laboratoire de Mecanique et Technologie
E.N.S. de Cachan, C.N.R.S.
U niversite Paris 6
61, avenue du President Wilson
F-94235 Cachan Cedex, France
ISBN 978-3-540-60980-3
This work is subject to copyright. All rights are reserved, whether the whole or part of the
material is concerned, specifically the rights of translation, reprinting, reuse of illustrations,
recitation, broadcasting, reproduction on microfilm or in any other way, and storage in data
banks. Duplication of this publication or parts thereof is permitted only under the provisions
of the German Copyright Law of September 9, 1965, in this current version, and permission
for use must always be obtained from Springer-Verlag. Violations are liable to prosecution
under the German Copyright Law.
Springer-Verlag Berlin Heidelberg 1996
Originally published by Springer-Verlag Berlin Heidelberg in 1996
The use of general descriptive names, registered names, trademarks, etc. in this publication
does not imply. even in the absence of a specific <tatement, that such names are exempt from
the relevant protective laws and regulations and therefore free for general use.
Typesetting: Fotosatz-Service Kohler OHG, Wiirzburg
Cover design: de'blik, Berlin
SPIN: 10710372 55/3012 - 5 4 3 2 1 - Printed on acid-free paper
Foreword
A new branch of science usually develops thus. Somebody publishes the basic
ideas. Hesitatingly at first, then little by little, other original contributions appear,
until a certain threshold is reached. Then, overview articles are printed, conferences
are held, and a first mention is made in textbooks, until specialized monographs
are written. Continuum damage mechanics has reached that status now.
To analyze or, if possible, to predict the failure of machine parts or other structures
is one of the main goals of engineering science. Consequently fracture mechanics
became one of its leading branches. It was based on the analysis of existing cracks.
However, especially under conditions of cyclic loading, this might be too late to
prevent a disaster. Therefore, the question regarding the precursory state, that is,
the evolution of internal damage before macrocracks become visible, was then
posed. One of the successful approaches to the problem was Weibull's theory which
examined, in a statistical manner, the " weakest link " in the material volume under
consideration. Unfortunately it proved too difficult mathematically to be applied to
complicated parts or structures. Therefore it was highly appreciated by the scientific
community when L.M. Kachanov published in 1958 a simple model of material
damage which subsequently could be extended to brittle elastic, plastic or viscous
materials under all conditions of uniaxial or multiaxial, simple or cyclic loadings,
so that it may be considered nearly universal. As usual in a living science many
questions are still open, and the limitations of the model need further research.
Jean Lemaitre, the author of the present work, has become one of the most dis-
tinguished representatives of damage mechanics. His monograph highlights all the
aspects of the subject, from the basic ideas, to theoretical and thermodynamical
perspectives, up to numerical or experimental applications. It also devotes suf-
ficient space to the problems associated with parameter identification. I like the
style of presentation, which is always vivid or humorous and never dry, so that
the understanding even of the more difficult passages is facilitated . An excellent
book has finally appeared, most valuable to any engineer or scientist dealing with
modem problems of strength of materials.
I wish it every success.
Notation xv
Bibliography 222
The choice of notation is always a difficult problem in a new field because all letters
have already been assigned a meaning. I have tried to be as classical as possible,
but I could not completely avoid using the same letter for different meanings. Be
careful and please forgive me.
Operators
X a scalar
xi components of a vector x
x ij components of a second-order tensor x
[xl matrix
X time derivative of x(x = dx/dt)
Xi ,j gradient of x
divergence of x
Xij,j
d,o,D differential operators
H Heaviside function: H(x) = if x < 0, H(x) I if x >
Ln Napierian logarithm
Symbols
a crack length
aijkl elastic moduli
A Gurson-Tvergaard material parameter
A crack area
d length of microcells
D damage variable
Ds,D" damage variables in shear and tension
Dij , D components and second-order anisotropic damage tensor
D ijkl components of anisotropic fourth-order damage tensor
Dc critical damage at crack initiation
DIe critical damage in pure tension
f yield function
l body force
F force
F function potential of dissipation
FD damage potential function
strain-stress function
vector of strain rate discontinuity
strain energy release rate
energy toughness of materials
Z intensity of electricity
Ikl unit second-order tensor components
Ikl
'J .
unit fourth-order tensor components
factor
specific kinetic energy rate
stress intensity factor
cyclic plasticity material parameter
creep material parameter
viscosity material parameter
stress concentration coefficient
Notation XVII
m mass
M material point
M strain hardening material exponent
n number of defects
n viscosity material exponent
ii unit normal vector
N number of cycles
N Norton's creep material exponent
NR number of cycles to rupture
T resistivity
T strain of isotropic hardening
R isotropic hardening stress variable
RXJ isotropic hardening material parameter
Ro,Ra stiffness
Rv triaxiality function
s specific entropy
s crack area
S sectional area
SD damaged sectional area
S damage strength material parameter
t time
T temperature
i1 displacement vector
v wave speed
V electrical potential difference
V volume
coordinate
unidimensional and components of tensorial kinematic hardening
stress variable
kinematic hardening material parameter
amplitude
dissipation
p mass density
stresses at microscale
(5H hydrostatic stress
D D deviatoric stress tensor
(5 ij , ( j
The damage of materials is the progressive physical process by which they break.
The mechanics of damage is the study, through mechanical variables, of the mecha-
nisms involved in this deterioration when the materials are subjected to loading. At
the microscale level this is the accumulation of microstresses in the neighborhood
of defects or interfaces and the breaking of bonds, which both damage the material.
At the mesoscale level of the representative volume element this is the growth and
the coalescence of microcracks or microvoids which together initiate one crack.
At the macroscale level this is the growth of that crack. The two first stages may
be studied by means of damage variables of the mechanics of continuous media
defined at the mesoscale level. The third stage is usually studied using fracture
mechanics with variables defined at the macroscale level.
When studying engineering materials such as metals and alloys, pofymers and com-
posites, ceramics, rocks, concrete, and wood, it is very surprising to see how such
materials, which have different physical structures, are similar in their qualitative
mechanical behavior. All show elastic behavior, yielding, some form of plastic or
irreversible strain, anisotropy induced by strain, cyclic hysteresis loops, damage
by monotonic loading or by fatigue, and crack growth under static or dynamic
loads. This means that the common mesoscopic properties can be explained by a
few energy mechanisms that are similar for all these materials. This is the main
reason it is possible to explain material behavior successfully with the mechan-
ics of continuous media and the thermodynamics of irreversible processes, which
model the materials without detailed reference to the complexity of their physical
microstructures.
lines of dislocations where atoms are missing. If a shear stress is applied, the
dislocations may move by the displacement of bonds, thus creating a plastic strain
by slip without any debonding as shown in Figure 1.1.
If the dislocation is stopped by a microdefect or a microstress concentration, it
creates a constrained zone in which another dislocation may be stopped. This
second process cannot occur without a debonding damage as shown in Figure 1.2.
Several arrests of dislocations nucleate a microcrack. Other damage mechanisms
in metals are intergranular debonding and decohesion between inclusions and the
matrix.
All these mechanisms create plastic microstrains.
- In polymers, damage occurs by the breakage of bonds that exist between the
long chains of molecules.
- In composites damage is the debonding between the fibers and the polymeric
matrix.
- In ceramics it is mainly microdecohesions between the inclusions and the matrix.
- In concrete, the early damage mechanism is also a decohesion between aggre-
gates and the cement with the complex influence of water.
- In wood, the weak point where damage occurs is the bonding of the celulosic
cells.
In all cases elasticity is directly influenced by the damage, since the number
of atomic bonds responsible for elasticity decreases with damage. This coupling,
I I
.L
I I
which occurs at the level of the state of the material defined here by the elastic
strain and the damage, is called a "state coupling" .
Another important property to consider is that the damage is always much more
localized than the strain. Remember that the damage, or debonding of atoms, is
4 Phenomenological Aspects of Damage
Brittle damage
The damage is called brittle when a crack is initiated at the mesoscale without a
large amount of plastic strain. Just to give an order of magnitude, let us say that
the ratio of plastic strain to elastic strain is below unity:
cp
-<
ce
This means that the cleavage forces are below the forces that could produce slips
but are higher than the debonding forces. The degree of localization is high. An
example is given in Figure 1.3.
Ductile damage
On the other hand, the damage is called ductile when it occurs simultaneously with
plastic deformations larger than a certain threshold PD' It results from the nucleation
of cavities due to decohesions between inclusions and the matrix foHowed by their
growth and their coalescence through the phenomenon of plastic instability. As a
consequence, the degree of localization of ductile damage is comparable to that of
plastic strain. An example is given in Figure 1.4.
Creep damage
When a metal is loaded at elevated temperature, for instance a temperature above
1/3 of the melting temperature, the plastic strain involves viscosity; that is, the
material may be deformed at constant stress. When the strain is large enough, there
1.1 Physical Nature of the Solid State and Damage 5
MPa cr
600
400
200
Fig. 1.3. a Fractograph showing microcracks damage around an inclusion in a silicon nitride ceramic;
b Tensile stress-strain curve up to rupture (after F. Hild)
are intergranular decohesions which produce damage and an increase of the strain
rate through the period of tertiary creep. As for ductile damage, the gradients of
creep damage are similar to the visco-plastic strain gradients. An example is given
in Figure 1.5.
a
o'(MPa) Damage
500
400
300
200
100
e
b 0 10 20 30 40 50 60 70 80 >10- 2
Fig. 1.4. a Fractograph showing microcavities damage at the grain boundaries of a steel (after L. Engel
and H. Klingele); b Tension stress-strain curve of A!S! 1010 annealed steel (after J. Dufailly)
damage localization is higher than for ductile or creep damage. Because of the high
values for the stress, the low cycle fatigue is characterized by low values of the
number of cycles to rupture, N R:
If the material is strain loaded, the damage induces a drop of the stress amplitude
as shown in Figure 1.6 for two stress-strain loops corresponding to the stabilized
cycle and a cycle close to the rupture.
1. 1 Physical Nature of the Solid State and Damage 7
a
-2 ep
.10
20
10
t (h)
b o 10000 15000
Fig. 1.5. a Fractograph showing intergranular cavitation damage formed at 800 C in a nickel-based
alloy (after L. Engel and H. Klingele); b Creep curve under constant tension stress of A 542 stainless
steel at 550C (after A. Benallal)
For metals, the damage can be either intergranular or trans granular microcracking
following slip-band arrests.
t- - - - -- - - - - - - - - - - --
- 0.2
----- ----~-------------------------
b
Fig. 1.7. a Fractograph of transcrystalline fatigue fracture in a nickel alloy (after L. Engel and
H. Klingele); b Cyclic tension compression curves for high cycle fatigue of A 316 stainless steel
Note that for brittle damage and high cycle fatigue damage, a stress-strain curve ob-
tained from a classical tension-compression test at the meso-scale usually does not
represent the "true" behavior for strain and damage because the space localization
induces microplastic and damage zones much smaller than those of the specimens.
Nevertheless, it is used because mechanical experiments at the microscale are diffi-
cult to perform; but please bear in mind that the results are averages of nonuniform
quantities over a mesovolume. The microhardness test may help to characterize a mi-
crovolume as it involves a size of the order of microns but its state of stress is complex.
1.1 Physical Nature of the Solid State and Damage 9
Fig. 1.6. a Fractograph showing ductile cycle striations of fatigue in an 18:8 austenitic steel (after
L. Engel and H. Klingele); b Cyclic tension compression curves for low cycle fatigue of A 316 L
stainless steel (after J. Dufailly)
the micro level where transgranular microcracking occurs only on some planes and
most often at the surface of the specimen by the mechanism of intrusion-extrusion.
The number of cycles to failure may be very large:
NIl> 100000 cycles.
As a consequence, the localization of damage is high and the drop of stress at the
mesolevel for an experiment similar to that of Figure 1.6 occurs later on Figure 1.7.
10 Phenomenological Aspects of Damage
Look carefully at the pictures of materials in Figure 1.8 taken with a microscope
and find the mechanism of strain damage or rupture for each case. The answers
are, in random order:
Fig. 1.8. Micrograph pictures taken from An Atlas of Metal Damage by L. Engel and H. Klingele and
also from D. Fran,>ois and Y.Berthaud for the two last pictures
1.2 Mechanical Representation of Damage II
Hint: if the answers are correctly matched to the pictures, the sequence of capital
reference letters should spell the name of a most famous wine.
It follows from the preceding section that damage may be interpreted at the mi-
croscale as the creation of microsurfaces of discontinuities: breaking of atomic
bonds and plastic enlargement of microcavities. At the mesoscale, the number of
broken bonds or the pattern of microcavities may be approximated in any plane by
the area of the intersections of all the flaws with that plane. In order to manipu-
late a dimensionless quantity, this area is scaled by the size of the representative
volume element. This size is of primary importance in the definition of a variable
continuous .in the sense of continuum mechanics. At one point, it must be the rep-
resentative effect on failure of microdefects over the mesoscale volume element.
It is similar to plasticity where the plastic strain Ep represents, at one point, the
average of many slips.
Let us consider a damaged body and a Representative Volume Element (RVE) at
a point M oriented by a plane defined by its normal ii and its abscissa x along the
direction ii Figure 1.9.
- Let 85 be the area of the intersection of the plane with the RVE;
- Let 85 D", be the effective area of the intersections of all microcracks or micro-
cavities which lie in 85;
- The value of the damage D(M, ii, x) attached to the point M in the direction
ii and at the abscissa x is:
-)
D(M ,n,x 8S D x
=~.
In order to define a continuous variable over the RYE for its deterioration to failure
in two parts, one must look at all the planes varying with x and consider that which
is most damaged:
D(M ,ii) = Max[D(M
(x )
ii x )] ,
, ,
Ir-D-(-M-'ii-)-=-~-:S-D--'I
It follows from this definition that the value of the scalar variable D (which depends
upon the point and the direction considered) is bounded by 0 and 1:
0::::; D::::; 1;
D = 0 undamaged RYE material;
-+
D =1 fully broken RYE material in two parts.
-+
In fact, the failure occurs for D < I through a process of instability (see Sec-
tions 1.2.4 and 3.1.3.).
Consideration of the simple one-dimensional case of a homogeneous damage pic-
tured in Figure 1.10 leads to the simple definition of damage as the effective surface
density of microdefects:
D_ SD
- S
t:F
.. +! .
50
...
-, I ~ 5
If all the defects are open in such a way that no microforces are acting on the
surfaces of microcracks or microcavities represented by 5 D , it is convenient to
introduce an effective stress 0- related to the surface that effectively resists the
load, namely (5 - 5 D ):
F
F
or L. -a_
I_1=D - ___a_--,I
A way to avoid a micromechanical analysis for each type of defect and each type
of mechanism of damage is to postulate a principle at the mesoscale.
In thermodynamics, the method of local state assumes that the thermomechanical
state at a point is completely defined by the time values of a set of continuous
state variables depending upon the point considered. This postulate applied at the
microscale imposes that the constitutive equations for the strain of a microvolume
element are not modified by a neighboring microvolume element containing a
microcrack. Extrapolating to the mesoscale, this means that the strain constitutive
equations written for the surface 85 - 85 D are not modified by the damage or that
the true stress loading on the material is the effective stress 0- and no longer a.
The following principle results:
"Any strain constitutive equation for a damaged material may be derived in the
same way as for a virgin material except that the usual stress is replaced by the
effecti ve stress"
14 Phenomenological Aspects of Damage
Elasticity law
This is direct state coupling through the concept of effective stress:
Undamaged material I Damaged material
D=O O < D<l
a a
ce = E ce = E( I - D)
Contraction for isotropic damage:
c22 = c33 = -vce;
E is Young's modulus of the undamaged material;
v is Poisson's ratio.
The elasticity modulus of the damaged material defined by the ratio E=
E= E(I -D).
Figure 1.11 shows an example of the variation of the elasticity modulus as ductile
damage progresses. It can also be evaluated from the microcrack pattern by means
of micromechanics (8. Budiansky, R. O'Connell, 1976).
Plasticity
This is a kinetic coupling on the evolution of plastic strain which has to be written
in the plasticity criterion used to derive the kinetic constitutive equations.
In order to model plasticity two kinds of strain hardening are usually considered:
- the isotropic hardening related to the density of dislocations or flow arrests;
- the kinematic hardening related to the state of internal microstress concentra-
tions. The corresponding back stress defines the center of the elastic domain in
tension compression (or in three dimensions).
1.2 Mechanical Representation of Damage 15
I
E \
I 0.5
100 99000 \
I
e
oOL.--2~0--'-~4~0--6-r0-~8~0-.Ll0~0---l'1~0- 2 0 0'---~20-~4~0--60
~-~8~
0 --I~OO----C'1~0 -2
Fig. 1.11. Variation of the elasticity modulus with damage for copper (99.9%) (after J. Dufailly. 1975)
If a y is the yield stress, R the stress due to isotropic hardening and X the back
stress, both functions of the plastic strain, the one-dimensional plasticity criterion
defining the current threshold of yield limit, and represented in Figure 1.12 is
a=ay+ R + X
or f = la-XI-R-a y = O.
f is the yield function from which the kinetic constitutive equation for plastic strain
is derived (see Section 2.1.3):
f = 0 f < 0
:lip -::j:. 0 if { ~nd { or
f = 0 j < 0
To achieve this, let us write the total strain as
-------------
-- - ~
-
~...,..::;;;:,;=
Oy r - - - - i - -__ OR
p
O~--_+---~------_4-----+
pD I
D I
t ~I[k
o
~
L..-_ _ _"-~_ _ _ _ ___I._ _ -+.
p Fig. 1.12. Schematic evolution of the plas-
ticity yield criterion with damage
16 Phenomenological Aspects of Damage
When damage occurs, according to the principle of equivalence, the yield function
f must be written as:
f = II~D-XI-R-O-y = 0
0- = (o-y+R+X)(I-D)
show that the damage equally decreases the yield stress, the isotropic strain hard-
ening stress and the back stress as shown in Figure 1.12.
Rupture criterion
The rupture at the mesoscale is a crack initiation which occupies the whole surface
of the RYE; that is, D = I. In many cases this is caused by a process of instability
which suddenly induces the decohesion of atoms in the remaining resisting area.
It corresponds to a critical value of damage Dc' which depends upon the material
and the conditions of loading.
The final decohesion of atoms is characterized by a critical value of the effective
stress acting on the resisting area. Let us call it 0- 00 because it is the maximum
stress that could be applied to the material (see Figure 1.12):
_ 0-
0-=---=0-
I - Dc 00
Damage threshold
Before the microcracks are initiated, creating the damage modelled by D, they must
nucleate by the accumulation of microstresses accompanying incompatibilities of
microstrains or by the accumulation of dislocations in metals. This corresponds in
1.2 Mechanical Representation of Damage 17
the pure tension case to a certain value EpD of the plastic strain below which no
damage by microcracking occurs:
Ep < EpD ----+ D = O.
Finally, the four main relations which comprise the basis of damage mechanics
are:
(J ..
Ep = for elasticity;
E ( l-D )
Il-D
-(J- - X 1- R - (Jy
.
= 0 as the plastic yield criterion;
where k is the correction factor due to the stress concentration in the vicinity of
the crack tip line to define the effective damaged area.
-u
Fig. 1.13. Single cracked meso volume element
18 Phenomenological Aspects of Damage
Ra = Ro(l - D) ,
with R = EI2 = EI
~"O I '
E being the Young's modulus of the material.
laWI 1 F2 dR
G = -2" aA F=const = -2" R~ 2Jrada;
and the stress intensity factor for a circular crack in an infinite body loaded in
mode I:
2 2F
K = -a00,j1i(i = - -12 ,j1i(i.
Jr Jr
For a case of generalized plane strain (1/ being the Poisson 's ratio),
K2
G = e(l-l}).
1- 1/ 2 2
- e a da ,
which can be integrated from a = 0 --+ Ra = Ro, to the current state (a, R,J.
The result is
R = [~ 16(1-1/2)a3 ]-'
a Ro + 3EI4
R = R ( _ 16( I - 1/ 2 )a3 R )
a H{) I 3EI4 o'
k -_ 16(137r
- V2) ~ '" ~
I - 1.55 I .
This means that for this particular case and for a/I I, the value of the damage
is always smaller than the crack density.
If the damage consists of microcracks which are not perfectly homogeneous in the
RYE, there is almost no chance of obtaining even one crack in any plane of inter-
section! Therefore the damage can only be evaluated by means of the lines which
indicate the intersection of the microcracks with the plane of observation. An equiv-
alent isotropic damage may be obtained as follows : Consider a micromechanical
element of n cells of size d x d x d, each containing a crack of different appar-
ent size a i (Figure 1.15). Assuming square cracks for simplicity, and a correction
factor equal to I, the damage in the plane of the crack of each cell i is
aT
)
D = OSD
1 as d2 '
20 Phenomenological Aspects of Damage
Fig. 1.14. Direct evaluation of damage (after Engel, KlingeJe and C. Lienard)
The equivalent isotropic damage in the plane of observation of the n cells is taken as
the mean value of the damage in each cell. This is the simplest way to homogenize
discrete quantities in order to define a continuous variable
n a2
~-t
D = i=1 d
n
If a micrograph of size 12 is considered with n cells corresponding to n cracks, the
surface of the micrographic may be expressed as
12 = n. d2 .
~a2
t
D
y'
An example of damage is shown in Figure 1.15., where grain boundary precipitates
of titanium carbide caused intercrystalline failure in a martensitic NiCoMo steel
subjected to triaxial stresses.
1.3 Measurement of Damage 21
~ /
, ~ )
~ ......
/'
~
Fig. 1.15. Micromechanic model and direct evaluation of damage (after L. Engel, H. Klingele and
C. Lienard)
D -'r ~a2
0.88
(J
cF =
E(I - D)'
This destructive method requires the machining of specimens in order to run me-
chanical tests. It assumes uniform homogeneous damage in the specimen gauge
section.
If E = E( I - D) is considered as the effective elasticity modulus of the damaged
material, the values of the damage may be derived from measurements of E,
provided that Young's modulus E is known:
This very useful method requires accurate strain measurements. Strain gauges are
commonly used and E is most accurately measured during unloadings. An example
22 Phenomenological Aspects of Damage
() (MPal
soo
400
300
a
Ol
Ol a
200 co
Ol
co
C1)
a
"
0"
Ol 'I
co 0I C
Ol
100 Ol I I
~
"
UJ UJ iii iii
00 10 20 30 SO 60 70 80 90 100 110
o
0.8 0.85 - - - - - - -- - - - - - - - - - -
I
0.6 I
I
I
0.4 I
I
I
I
0.2 I
I 2
I ep .10
OL-_.---r--~_,--_.--._--.__.--_r--r_~~
o 20 40 60 1 80 100
Fig. 1.16. Measurement of ductile damage on 99.9% copper at room temperature (after J. Dufailly)
is given in Figure 1.16 for ductile damage of copper induced by large deformations,
where a v is the true stress a v = a( 1 + E)
This technique may be used for any kind of damage as long as the damage is
uniformly distributed in the volume on which the strain is measured which is the
main limitation of the method. If the damage is too greatly localized, as for high
cycle fatigue of metals, for example, another method must be used.
Some other precautions related to non-linearities have to be taken. At the beginning
and at the end of the unloading paths in the plane (a, E) there are small nonlin-
earities, owing to viscous or hardening effects and also due to the experimental
devices as well. It is best to ignore them and to identify E in the range
It is most important to always use the same procedure to evaluate E and the
evolution of E.
For ductile or low cycle fatigue damage in metals, the procedure could be somewhat
disturbed by an early decrease of E at low strain levels or during the first cycles.
This is due to microplasticity related to reversible movements of dislocations, and
1.3 Measurement of Damage 23
/
0.5
/
II
0.4
II
/
r<)
co
-I'
0.3 X3
/
o
0.2
!!
III~
II
!II
I
0.1
6~
o 2 4
- a'n (MPa)
Fig. 1.17. Evolution of brittle damage in a shear direction of a three-dimensional carbon-carbon com-
posite (after M. Poss)
24 Phenomenological Aspects of Damage
Concrete
. Xz
25c
Equivalent strain 'i = .r;r:;:;r
1 2 = f:1 X f2
t
28~
.
,
'"Q 26
.,
)( 24 \.
0
0.. 22
::iE
20
w'"
18
16
O ~--~----~----~--~--~---~
0.5 1.0 1.5 2.0 2 .5
xl0 3
0.5
Fig. 1.18. Evolution of brittle damage in compression for a concrete (au =40 MPa) (after J. Mazars)
For frequencies higher than 200 kHz the longitudinal wave speed v L and the
transversal wave speed vT in a linear isotropic elastic cylindrical medium are:
2 E 1- v 2 E I
vL = P(I + v)( I - 2v) VT = P2(1 + v) '
where E is the Young's modulus, p the density and v Poisson's ratio. A measure
of the longitudinal wave speed of a damaged material gives
E I-v
p (I + v) (I - 2v)'
1.3 Measurement of Damage 25
where E and p are the actual damaged elastic modulus and density. Poisson's ratio
does not vary with damage if elasticity is isotropic , and damage a scalar.
The damage is calculated by:
The method considered here is destructive: indeed, to measure the speed v L or the
time a wave takes to propagate through a certain thickness, one needs to limit that
thickness by two surfaces. If the space distribution of the damage is not uniform,
the thickness must be of an order of magnitude coherent with that of the RYE,
therefore that requires the body to being analyzed to be cut into parts.
The present limitations of this method are related to the size, to be analyzed, which,
for metals, is too small in comparison with ultrasonic transducers size and accuracy
of time measurements. Nevertheless the method gives good results for concrete.
The ranges of frequency to be used are:
I_ a
I- D
-xl-R-a II
=0
01
tot / toto
0.6
0
0.4
,
I
0.2
o 0.5
(1' / (1'r
Fig. 1.19. Evolution of brittle damage from the propagation lime ~t of ultrasonic waves in concrele
(after Y. Berthaud).
H = k'a,.
This threshold corresponds to the actual yield stress,
as = (all + R + X) (I - D).
then
H = k'(all+R+X)(I-D).
In fact, the microhardness test itself increases the strain hardening by an amount
which corresponds. to a plastic strain e~ of the order of 5 to 8%. H is then always
related to (el' + e~), eli being the current plastic strain.
1.3 Measurement of Damage 27
H(*o P +oH)
Jl
= k'(a'I+Rx
+X x ) = const.
and let H be the microhardness at the point where the strain hardening has to be
evaluated:
H = k'(a y + R+ X)(I - D)
H
or H = --...1t.(a y +R+X)(I - D) ;
ay
28 Phenomenological Aspects of Damage
H
/
3000
1000
900
800
700
600
SOO
400
300 o'y=320 D=I-.t!.
0.5 H*
200
lOa Fig. 1.20. Evolution of ductile damage by
0.1 microhardness on AIS! 1010 at room tem-
00 10 20 30 40 perature (after R. Billardonj
then
or
Several other methods, which are based upon the influence of damage on some
physical or mechanically measurable properties, may also be used.
1.3 Measurement of Damage 29
0.3
0.2
0.1
Fig. 1.21. Damage field on a plastic deformed built notched specimen as obtained from microhardness
measurements on A)S) 1010 (after J. Dufailly, R. Billardon)
Variation of density
In the case of pure ductile damage, the defects are cavities which can be assumed
to be roughly spherical; this means that the volume increases with damage. The
corresponding decrease of density is measurable with apparatuses based on the
Archimedean principle.
If (p - p) / p is the relative variation of density between the damaged state p and
the initial nondamaged state p, it easy, by means of micromechanics, considering a
spherical cavity of radius l' in a spherical R.V.E.of initial radius R and mass m, to
derive the following relationships between the surface damage D and the variation
of density or porosity, assuming no residual microstress:
m
p = 4 '
-7rRJ
3
m
p-p RJ _rJ
- I
p RJ + rJ R' + 1'3 '
bS D 7r1'2 ( rJ ) 2/3
D
bS 7r(RJ + rJ)2fJ RJ + rJ
30 Phenomenological Aspects of Damage
D ( 1--p)2/3
p
The effective intensity of the electrical current can be defined in the same way as
the effective stress was defined, using the surface definition of damage:
--: 1
1 = I-D'
I being the intensity which actually exists in the cohesive parts of a damaged
volume element.
Given the potential difference V, Ohm's law for a non damaged element of length
I, area s and resistivity r is written as
V = r-l;
s
0.15
~ 0.10
a.
-.
a.
"'? 0.05
Q
14
12
.
M
52 10 c
0 8
6
4
2
eeq
0 0.5 1.0
Fig. 1.22. Evauation of ductile damage on steel for different strain hardenings (after F. Moussy)
1.3 Measurement of Damage 31
- l -
V = f-i,
8
where f is the resistivity affected by the damage by means of the change of volume
only (p - p)1 p (p being the density). Bridgman's law gives f as:
-~
V r 8
- -
f
V
81-D
D = I- ~~ with ~ = 1+ K D3 /2.
Vr r
For small values of D, the correction term (fir) due to the volume change is close
to I (for instance, D c::: 0.1 ~ f i r = 1.064); then
V
D"'I--
- V
This method is known as the "potential drop method".
Some comparisons between measurements of damage obtained with this method
and from elasticity modulus variations D = I - EI E or methods described
hereafter are shown in Figure 1.23.
As the damage measurements are made during tension tests, a correction of these
examples is required to account for the length change l( I + f), the area change
8(1 -2VCe -lOp) and the volume change due to elastic strain: 8l(1 -2v)c e . The
agreement between the results obtained by the different methods is sufficient for
the practical accuracy involved.
The influence of damage on plasticity may be used to measure the low cycle fatigue
damage.
The one-dimensional law of cyclic plasticity at stabilization may be written as a
power relationship between the amplitude of stress 6..(7 and the amplitude of plastic
strain 6..c p (see Section 3.4.1.):
32 Phenomenological Aspects of Damage
!I
0.2 0.2
!I
0 308 MPa I 0 0' = 240 MPa !I ~
<.>
E
I <.>
!I
I ~it
0.1 ,
~ 0.1 !I I E.2 G>
>
~ / G>
i_ -~
~
0.5 0.5
a N/NR c tltR
0.1 0.1
b d
Fig. 1.23. Comparison between damage measurements by variation of the electrical potential and by
a the elasticity modulus change on stress-controlled fatigue of A 316 stainless steel, b the elasticity
modulus change on strain-controlled fatigue of A 316 stainless steel, c the tertiary creep method on
creep damage of IN 100 superalloy at IOOOoC. d the cyclic plasticity response on stress controlled
fatigue of IN 100 superalloy at IOOOC (after G. Cailletaud)
[K,(~~D) 1 M
1.3 Measurement of Damage 33
12 -
10.23- 10 r-- ilo'*
0
0...
:E 8 cp =It 0.3 'I.
N
S2 6 -
0
<J
4
2
I I I I
0 5 10 15 20
1.0 N'10 2
0.8 -
0.6 I-
0
0.4
'I
0.2 :1_ one mesocrack
I
I I develops
0 5 10 15 20
N .10 2
Fig. 1.24. Evolution of low cycle fatigue damage on AISI 1010 deep drawn steel
34 Phenomenological Aspects of Damage
Assuming that the damage process begins at the end of secondary creep, during
tertiary creep one may write:
(!:r,
from which one derives
D = 1-
~ 2
a.
w
e -*
o
t (h)
. t
0.8
I:
/:
o
0.6
r
(
0.4
I :
I :
0.2 +- one meso crack
~ develops
:1
o 20 30
t (h)
Fig. 1.25. Evaluation of creep damage on IN 100 superalloy at lOOOee (after H. Policella)
1.3 Measurement of Damage 35
This method yields good results which are in accordance with those obtained
by measuring the variation of the elasticity modulus. An example is given in
Figure 1.25.
Acoustic emission
It is a good method for detecting the location of the damaged zone, but the results
remain qualitative as far as the values of the variable D are concerned.
By way of conclusion, Figure 1.26 offers, in the form of a wine selection chart (!),
some advice for choosing the proper method of damage measurement, depending
on the kind of damage involved.
Low High
cycle cycle
Damage Brittle Ductile Creep fatigue fatigue
8S D
Micrography D= * ** ** * *
8S
Density D= ( pr/
1- -
1
** * *
P
Ultrasonic {i2
D= I - --.b.
V2
*** ** ** * *
waves L
Micro-hardness D= 1- -
H ** *** ** *** *
H*
Electrical V
D= 1 - -:::- * ** ** * *
resistance V
three stars *** means "very good" like "Bourgogne Chambertin 1985"
two stars ** means "good" like " Bourgogne Pommard 1983 "
one stars * means "try to see" like " Beaujolais 1989"
no stars means "do not try" like ???
A votre sante I
When the strain hardening is saturated and the damage still zero,
then aM
I-D
Fig. 1.27. Cyclic strain controlled test on MAR M 509 stainless steel.
0
1.0
0.8 -
0.6
0.4
0.2
Fig. 1.28. Evolution of the damage
0 measured by the stress amplitude drop
0 0.2 1.0 and by the elasticity modulus change
1.3 Measurement of Damage 37
From Figure 1.27b, draw the evolution 0/ D as a/unction o/the numher o/cycles N
and compare in Figure 1.28 with the evolution of D calculated Fom the elasticity
modulus change measured on Figure 1.27a:
Chapter 2
Thermodynamics and Micromechanics of Damage
- Internal variables:
* " the elastic strain tensor of component cij associated with the stress tensor
0';
* p the plastic strain tensor of component cfj. It turns out that its associated
variable is -0' due to the partitioning of strains:
._(2. )I/2
P -
p .p
"jCijCi j ,
Fig. 2.1. Plastic yield locus in the principal stress space (The trisectrice is normal to the page!)
2. 1 Three-Dimensional Analysis of Isotropic Damage 41
to be isotropic, it has the same value in all directions, and the scalar
D = 8S D ~
8S \in
characterizes completely the three-dimensional state of damage in the RYE at
the point considered.
Let us call Y the associated variable of D which will be derived from the
state potential. As D is dimensionless and since the product - Y b is the power
involved in the process of damage, we already see that - Y is a volume energy
density.
Figure 2.2 is the chart of variables involved in the continuous mechanics of thermo-
elasticity, plasticity and damage.
Let us postulate that the state laws are derived from a state potential: a continuous
scalar function, concave with the temperature, convex with the other state variables
and containing the origin. Taking the Helmholtz free energy,
For elastoplasticity or elastoviscoplasticity, the strains act only through their dif-
ference e - eP = ee:
eP ..... . . .. . . . . - {j
r ..... .. . . .. R
9 0 ..xD
D ..... ... . .. .. . .... Y
1- D 'J 1- D '
since the damage is isotropic, it has the same influence on all the components
of the Cauchy stress tensor. For anisotropic damage, the variable is no longer a
scalar and the effective stress needs the introduction of a more complex operator
(see Section 2.2). The difference of the behavior in tension and compression is
not taken into account here either, it will be, however, in Section 2.3.4.
- no state coupling neither between plasticity and elasticity nor between damage
and plasticity,
- classical expressions for isotropic hardening and kinematic hardening giving
saturation of hardening for large plastic strain (see Section 2.1.4),
- isothermal processes:
a ij
a'lj; =
= P~ e (
a 'i j klEkl 1-
D) ,
UE ij
where a is the fourth order elasticity stiffness tensor or, by inversion for the
isotropic case,
e _ I +v ~ _ ~ akk
Eij - E 1 - D E I _ D bij ,
E being Young's modulus, v Poisson's ratio and bij the Kronecker delta. Values
for E and v are given for some materials in Section 3.5.
The isotropic strain hardening scalar stress is expressed as:
R oo and b are two parameters which characterize the isotropic strain hardening for
each material (see Section 5).
The kinematic strain hardening tensorial stress is expressed as
D a'lj; 2
X ij = p~ = -3 X -x/yo:,j.
uO:' J
Xoo and I are two parameters of nonlinear kinematic hardening written with the
factor 2/3 to ensure a simple expression in one dimension. They must be identified
for each material (see also Section 3.5).
The associated variable for D is defined as
Y a'lj; I ~~
= p aD = -iaijklCij Ckl
In order to work with a positive quantity, let us write
Y = -Yo
Let us now find the relationship between Y and the elastic strain energy density we.
By definition,
dW e = (Jijd<j.
Integrating with the law of elasticity, and assuming no variation of damage, that
is, D = const, yields.
IY ~I
Y is also equal to one half the variation of the strain energy density corresponding
to a variation of damage at constant stress: d(J ij = 0 (J.L. Chaboche 1976).
Starting with the law of elasticity,
d(Jij = aij hl [( 1 - D)dck:l - ck:ldDj = 0,
d ee dD
or ci.:l = ch:[ 1_ D '
together with
dweliT=const
or dW e I
dD iT=const
and from the definition of Y = - y
y = ~ dW e I
2 dD a=const
This last equation allows Y to be called the strain energy density release rate.
This is the energy released by loss of stiffness of the RYE in which the damage
occurred.
44 Thennodynamics and Micromechanics of Damage
1-211 ~
E I-D
2
Y =
We
I-D = 2
lJ e q
2E(I-D)2 ["3(1+11)+3(1-211) ( lJIJHeq ) 2] .
H
() e q
is the triaxiality ratio, which plays a very important role in the rupture of
Rv = ~(I
3
+11)+3(1-211) H)2 (IJ
lJ eq
2.1 Three-Dimensional Analysis of Isotropic Damage 45
~l '
-a* 0 0
0 3
[ a' I * I *
[a] = ~ 0 aH 3 a , [aD] 0 --a 0
0 3
0 0 I *
--a
3
::' :
ap'l
~' }
3
W
e 2E(1 - D) '
Rv = I
Writing the equality wf'(a*) = wp (O') = Y(I - D) ,
a;",R,/
2E(1 - D) 2E(1 - D)
defines the damage equivalent stress as
Ira-*-=-a-p-q-R-~/-2-'1
This differs from the von Mises equivalent stress by the triaxiality function, which
is in accordance with the physics described in Section 1.1: plasticity is mainly due
to slips which do not depend upon the hydrostatic stress; damage is debonding
influenced by the hydrostatic stress or the triaxiality ratio. This phenomenon is
influenced by Poisson's ratio, which governs the elastic volume change; the factor
RI/ increases in most cases when the Poisson 's ratio decreases.
As an example, let us consider the damage equivalent stress for a case of plane
stress:
[~
46 Thermodynamics and Micromechanics of Damage
Figure 2.3 represents this expression in the plane (~, ( 2 ) for different values of
a* a *
Poisson's ratio.
Having defined all the state and associated variables, a second potential will give
the kinetic constitutive relations to describe the evolution of the phenomena.
First of all, we must write the second principle of thermodynamics to ensure the
validity of the model chosen. Let us start from the Clausius-Duhem inequality,
. . Ti
a i/: ij - p( 'Ij! + sT) - qi T : : 0,
where if is the heat flux vector associated with the temperature gradient
grad T (= 8x =8T
i
T .)
"
8'1j! ) .e
( a ij - P8E;'fj Ci j -
( 8'1j! )
P s + 8T
.p -
t + a ijCij 8'1j! . 8'1j! .
P 8r r - P8a ij a ij
8'1j! . Ti
- p - D - q """:'" > 0
8D ' T -
11=0
0.25
0.28
0.32
Fig. 2.3. Representation of the damage criterion
0.5 in the case of plane stress (after D.Baptiste)
2.1 Three-Dimensional Analysis of Isotropic Damage 47
aii~j - Ri' - X ij
D
aij -
-'
YD - q.i;T 2: 0.
-yiJ 2: 0.
The laws of plasticity or viscoplasticity coupled with damage are derived from this
potential by means of a scalar multiplier which is always positive. This ensures the
normality condition of yielding for plasticity, which does not depend upon time
explicitly, and the viscosity property for the viscoplasticity. Within the so-called
"normality rule of generalized standard materials,"
of
-A
ou
of
-oRA (Visco) Plasticity constitutive equations;
of .
- oX DA
of /\\
D = - -= ---7
.. I
Kmetlc f
aw 0 damage '
evolutIOn.
oY
In order to define the yield condition, we need a loading function which generalizes
in three-dimensions the one-dimensional plasticity criterion in tension introduced
in Section 1.2.4, namely,
I-1-I J
D
-XI-R-IJ
y
=0
Let us work with the von Mises criterion which states that the (visco) plastic strains
are governed by the elastic shear (or deviatoric) energy density, the so-called "J2
theory". The corresponding equivalent stress is
lJ eq = (~1J31J3) 1/2 ,
1
(TD being the stress deviator:1J3 = lJ ij - 3lJkkDij' Together with kinematic hard-
ening, the von Mises criterion is applied to define the size of the yield locus
regardless of the translation X D defined in Section 2.1.1 and shown in Figure 2.1.
Therefore, it acts upon the difference between (TD and XD.
Furthennore, in the presence of damage, the coupling between the damage and the
plastic strain is written in accordance with the principle of strain equivalence (see
Section 1.2.3). The yield criterion is written in the same way as for a non dam aged
material except that the stress is replaced by the effective stress, which, for isotropic
damage is:
_ (T
(T = - - .
I-D
Then, the loading function f is written as
If = (a-D-XD)eq-R-lJ y I
with (a- D _ XD)
eq
= [~( 1J3 _XI)) ( 11J3
2 1- D - D
_XI)) ]1/2
Z] Z]
Remember that XD is a deviator and IJ y is the yield stress of each material which
depends upon the temperature (see Section 3.5).
__ . _ ___ . _ _. ____..--f:
.. 0 ...........
. /Ioading
/1<0 f: 0
unloading "-- f < 0 /
/
/
o Fig. 2.5. Loading-unloading condition in one di-
mension for plasticity and damage
Viscoplasticity
The plastic strain becomes a time-dependent phenomenon as the temperature in-
creases. For metals, this plastic viscosity is important if the loading occurs at a
temperature above about a third of the melting temper~ture T m:
I 0
T 2 3Tm =K.
For the uniaxial state of stress without any damage, Norton's law expresses the
relation between the steady state plastic strain rate and the viscous stress:
(1v
-- K vC pI IN ,
where K" and N are two material parameters. This relation is usually applied for
low plastic strain rates and does not account for the saturation of viscosity where
the strain rate is very high. Then for damage mechanics dealing with conditions
close to failure , it is more appropriate to work with an expression like:
(1" = K oo [I - exp ( - ~ )] or i p = Ln (I _; : )- n ,
where K ex:, and n are two material parameters (see Section 3.5)
(i P --+ 00) => ((1" --+ K ex:, = const).
The stress a is now the sum of the yield stress, the strain hardening stresses Rand
X as for plasticity, and a viscous stress (1,, :
a = (1y + R + X + (1v'
which means that the loading function f may be positive:
f = la-XI-R-(1y = a" > O.
Referring back to the multiplier ~ of the three-dimensional case and using the von
Mises function ,
ifj of ~ , f
oa = [(u D -X D )eq-R-(1y] > o for visco plasticity;
ij
50 Thermodynamics and Micromechanics of Damage
by extension of the uniaxial case, ). is taken such that (see also Section 2.1.4)
.
p = 1 - D = Ln
~ (f )-n
1 - K 00
I) As already mentioned in Section 1.2.4, below a certain value of the plastic strain,
a threshold cpD' no damage occurs. This allows us to identify the (visco) plastic
constitutive equations in their (large) domain of validity without considering any
coupling with damage.
2) The kinematic back stress X , defined as the locus of the center of the elastic
domain,
- increases with the plastic strain;
cr
-;h;--
tR ,
cry /
cry i E<E
E ,_-
I
X
-~--
/
I
I
\
\
I
0
I 0
1 1 I
1
I , I
-~
I
..J
~ ,-"
Those equations, valid for one-dimensional tension, serve as a guide to writing the
three-dimensional potential of dissipation where the kinetic coupling acts only by
the effective stress deviator aD /(1 - D). In the framework of the "State Kinetic
Coupling theory", it corresponds to:
3
F = ( a- D -X D) eq-R-ay + 4X
D D
X ij X ij +FD (Y:(r.D))
00
The first three terms correspond to the loading function f. Working with a potential
equal to the loading function F = f gives rise to the so-called "associate plasticity
theories." The fourth term _3_
4Xoo
XB XB is responsible for the nonlinear kinematic
hardening, including the factor 3/4 to ensure a simple expression in one dimension.
The last term F D is the part of the potential from which the kinetic law of damage
evolution is derived. Experiments show that it does not depend explicitly upon a,
X or R:
iJ = _ BF ). = BFD ).
BY BY
The analytical expression for F D will be discussed in detail in Chapter 3.
Now let us derive the general constitutive equations valid for plasticity and vis-
coplasticity:
.p
Eij
52 Thermodynamics and Micromechanics of Damage
I
3 (a&
2.. 2 I=D - X i j I=D -
D) I oX& .
a;;:;
2. (~_XD)
I-D
A,
eq
.p
3 ((jlJ-XlJ)
'J IJ
~
p = (~ififj )'/2 ,
which is, for the one-dimensional case, simply the absolute value of the strain rate
ip:
p= lipl
Let us calculate p from the expression of ifj:
(
2 ) 1/2 fi 3 ~
P= :3ififj = V:3 2. ( (jD )
I_D-XD (I-D)
eq
a& - XiD)
. [( I=D j
(I=D
a& - XiD)]1
j
/2,
which, together with the definition of
( ~_XD)
I-D eq
= [ ~2 ( I afJ- D -XlJ) ( I afJ- D -XlJ)] 1/2
'J IJ
yields
~
p= I-D
The constitutive equation for the isotropic strain hardening is
of .
r= - ORA = A.
This shows that
Ir = ~ = p( I - D) I
When no damage occurs, the state variable for isotropic hardening is the accumu-
lated plastic strain,
D = 0 -+ r = p.
For incremental calculations it is more convenient to work with a constitutive
equation for the associated variable R. From the state potential,
R = Roo [l-exp(-br)]
or R = R oobexp( -br)r
or R= b(Roo - R)i
2.1 Three-Dimensional Analysis of Isotropic Damage 53
(1E _ XD)
( I-D 8( -Xm
3 ZJ 8XIJ 3
XDij A,
__ -XD
2-'--(:--cr""""D==-----'-:--)--'-- + 2X
(X . ZJ
ZJ
00
I-D eq
_
(xij - Cij I -
.p ( D) 3
- 2Xoo ij A ,
xD ;
or, with the definition of X D from the state potential in order to have a constitutive
equation for this associated variable,
D 2 D 2 .
X ZJ = 3X oo l'(Xij or X ij = 3Xoo l'(Xij'
Plasticity (f = 0)
8f
8(~-XD)
I-D eq
3 (~I-D
-XIJ)
ZJ
- X D)
2 ( -cr-D
8(1ZJ I - D'
\-D eq
54 Thermodynamics and Micromechanics of Damage
aD)
aj
a(~_XD)
I-D 3
(t=D- ij D
X ij
eq
aX8 -2(~_XD) ,
I-D eq
aj
aD
a(5-
aD
x t = ~ (~[ a8
aD 2 1- D
-XIJ ] [a8
'J 1- D
-XIJ]) 'J
1/2
3 a ijD ) a ij D
-2
12
( ---x
1- D
D
ij (I - D)2
- 2 ( ___
uD XD )
I-D eq
aj = -I
aR .
Then:
. 3
( 1~jD ) ( .
D - X8 a a ijDD )
XD)
ij . D .
j = 2 (~ _ 1- D - X ij + (I _ D)2 - R = o.
I-D eq
D
X ij =
2.1 Three-Dimensional Analysis of Isotropic Damage 55
( ~_XD)
I-D
(I-D) .
~q
D)
( I=DIJij D
- Xij
. XOO-Y -2(~_XD)
3
8Y
8FD
( D D
IJ ij )
-yX ij +(I_D)2 +b(Roo -R)
I-D eq
D
The damage influences the plasticity by the effective deviatoric stress afl ~
J-D
and by the partial derivative of the damage potential 8:-:.
This case is much simpler, as the loading function gives the viscosity effect by the
viscosity law:
p= I ~D = Ln ( I - /xJ -n
A = Ln [ 1- (iT -X )e q. - R - I J y
D D ]-n
J-D Kx
in Section 3.1.4)
-D XD .
.p 3 (Tij - ij A
eij - 2 (iT D - XD)eq 1 - D'
R = b(Roo - R)~,
. D _ "I
Xij -
[23" X .p (
ooeij 1 - D) - XDJ
ij A ,
. BFD
D = BY A,
(T ..
- tJ
with (Tij = 1 - D'
[23 ( 1 - D
D -
) (D
(Tij
1_ D -
. D
Xij
(Tij D
X ij
)]' /2 ,
D 1
(Tij = (Tij - 3"(Tkk 8ij'
Plasticity (f :S 0)
if f = 0 and j = 0
3 -D D (Tij BFD D
-2((Tij - X ij ) (
-D
1 _ D BY + "IXij ) ]
if f < 0 or j = 0
~ =0
Viscoplasticity (f > 0)
.
A [D
(iT -X D)e -R-(T ]-n
- - = L n \- q y
I-D Koo
2.1 Three-Dimensional Analysis of Isotropic Damage 57
-~J '
0
~ [~ ~l ~ " [~ ~ 1+ [~
0 0 1
--
[a] a 0 lEI -/.I
'p 2
0 0 -/.I 0 0
2 2
- 0 0 - 0 0
3 3
I I
[aD] = a 0 -- 0 [xD] = X 0 - - 0
3 3
I 1
0 0 -- 0 0 --
3 3
From the constitutive equations obtained previously with D = 0 and iJ = 0,
f = la-XI-R-a y ,
C = c e + cP '
a
c e = E'
a-X .
Cp = la _ XI A = Asign(a - X),
R= b(Roo - R)~ ,
2 2 2
-X
3 = -X
3 ool'VE p - -'VXA
3 ' ,
~ = r = p = IEpl.
32
X)a 23(a -
For plasticity, A = --------'=-=-----=-3-=2----'
la - XI [Xx; ,+ b(Roo - R)]- 23(a - X),X
or
For viscoplasticity,
C = Ce +c P '
a
Ce = E'
R= b(Roo - R) ICpl '
58 Thermodynamics and Micromechanics of Damage
i p = Ln [1 -
10 - XI- R -
Koo
0
Y
]-n for viscoplasticity
For monotonic loading (Iipl = ip) the two strain-hardening equations may be
integrated with the initial condition lOp = 0 -+ R = 0 and X = 0:
- Xoo and /, for the kinematic hardening. The partition of stresses is made as
shown in Figure 2.6, the plot of X as a function of lOp allows us to identify X x
and /, either on a semi-Log diagram or by numerical identification through the
least square method, for example;
- Rand b for the isotropic hardening by the same method.
A schematic example is given in Figure 2.7.
For viscoplasticity two more parameters need to be identified. The most convenient
method is to derive them from a relaxation test at constant strain for a large value
of the plastic strain when the strain hardening is almost saturated (X ~ X x' R ~
R x ) and when the damage has not yet occurred:
o
10 = const -+ i = 0 -+ ip +i p = 0 -+ ip -_.
E'
if = Ln (0 -
1_ Xx - Rx - 0)-n Y
E Kx
If the stress is recorded as a function of the time and if (X x + R x + 0 Y)is known
from the above method for a very low strain rate tension compression test, one may
calculate i p and plot it as a function of 0 ; then n and K x are the values which
2.2 Analysis of Anisotropic Damage 59
fit best with the experimental points. A systematic method consists in a numerical
identification of both coefficients nand K :xl using the least square method.
o'(MPa)
+0'
e r (%)
~
0.05 0.1 e(%)
O~----~----,-.
-3 3 e(%)
0.05
Fig. 2.8. Cyclic test result on AISI 316 L stainless steel at room temperature (after J.L. Chaboche)
Nevertheless all materials under special loading conditions and most brittle mate-
rials develop anisotropic damage for which the damage variable can no longer be
a scalar.
[a] = [
0'1
0
0
0'2
0
0
1 if 0'1 < 0'2 < 0'3'
o 0 ~
I-D
area of the strain-equivalent undamaged configuration. Let us develop this idea for
anisotropic damage.
Consider a RYE of a damaged material unloaded in its current configuration and a
sectional area defined by its orientation ii, its surface 8S and its shape (Figure 2.9).
Consider also the strain-equivalent undamaged configuration in the same coordinate
system. The R.Y.E. has changed and particularly the sectional area has been trans-
formed into a section of different orientation ii due to the anisotropy of damage
and Of the smaller area 85 (Figure 2.9). Let us assume that the shape remains the
same; later on we shall see that this hypothesis corresponds to orthotropic damage.
ii8S is transformed into fi8S
By analogy with the one-dimensional case, the second-order tensor which trans-
forms the vector ii8S into the vector fi8S characterizes the damage and is denoted
by (I - D) where 1 is the unit second order tensor
(I - D) . ii8S = fi8S or (Iij - Dij )n j 8S = ni8S
D is the second-order tensor damage variable. It is symmetric and has the three
principal values D l , D 2 , D3 in the principal coordinates (Xl ' X2 , X 3) which cor-
responds to the reduction of the areas of a tetrahedron drawn in the Xl' X 2 , X3
0' 1- D i1
s -
s
one dimension
I
_['1<
e
L:2 x1 L::: x1
3 dimensions
Fig. 2.9. Isotropic damage in one dimension and anisotropic damage in three dimensions.
62 Thermodynamics and Micromechanics of Damage
TiDS = Ti DS.
Introducing the effective stress tensor ij as
we can write
R
.....
n
'"'-'
OA, = (1-0, lOA,
x2
Q
'"'-'
OA) = ( 1- 0) I iiA)
'"'-'
0A 2 = (1-02 10A 2
but this leads to a nonsymmetric tensor. As only the symmetric part accounts for
the constitutive equations of elasticity and plasticity, let us consider
iT = 2I [u(I-D) - 1+(I-D)-I u ].
When the principal directions of stress and damage coincide (see Exercise 2.2.4.),
0"1
I-DI
~
[aJ 0
l-D2
0
0"3
0 0
l-D3
The definition of D as a symmetric second order tensor and of its associated
effective stress are the basis for developing a complete theory of the orthotropic
damage mechanics in the same maner as for the isotropic case. It is only "a bit"
more complicated! In particular, the associate variable also becomes a second-order
tensor Y.
The Helmholtz free energy state potential for linear elasticity coupled with damage
is written as a scalar in a form similar to that in Section 2.1.1:
or
By definition of the effective stress, (ai jkl . iikl~s) is the operator of damage. It is
a fourth-order tensor, by analogy with the one-dimensional case let us call it:
(/ D ) -,
ijkl -
--I
i jkl = a i jrsarsk l'
where D i jkl are the components of the fourth order tensor damage variable D, and
I ijkl are the components of the fourth-order unit tensor I for symmetrical fourth
order tensors:
D ij kl = I ijkl -
- aijr s ar skl '
-I
or D=I-a : a- 1 ,
where the symbol (:) means the contracted product upon two indices.
The associated variable Y is also a fourth order tensor:
I -
.1.
'f/ = 2p a i jklCij Ckl
e e
+ .'f/p
1.
'
Y
i j kl = "2aklr
l ee
sCtJCr s'
This theory corresponds to the general case of anisotropy. We shall see in Sec-
tion 2.3.3 that micromechanics of a general pattern of microcracks leads to damage
representation as a combination of a scalar, a second-order tensor and a fourth-order
tensor.
independent scalar variables which characterize the damage that influences the
elastic shear energy and the elastic hydrostatic energy:
- D s acts on the deviatoric stress components;
- D n acts on the hydrostatic stress.
The damage may be isotropic or anisotropic if the material is not subjected to
damage in certain material directions.
For the general case of anisotropic elasticity, but without any damage, the comple-
mentary energy density is
1S
we = 2" ij klai jakl ,
where S is the compliance fourth-order tensor S = a-I. It may be split into two
parts as in Section 2.1.2 by means of the stress deviator (j D and the hydrostatic
stress a H :
aH = 3akk '
a[)
tJ
= a t J - aHDtJ'
1 D D 1 2
2" S ijk/aij a kl + 2"SijklaHDijDk l'
~~
shear energy hydrostatic energy
H _OW'
__e OW
a {Y s or Ys = oDe
oDs s
- the two variables associated to D s and Dn' _ OW' ,
_ OW'
__e or Yn = __
Yn e
oDn oDn
the damage dissipation being
YsD s + YnDn '
With this theory it is also possible to consider anisotropic damage for which
D n i= 0 in certain material directions, and D n = 0 in the non damageable
directions. To illustrate this point, let us develop the whole theory for an elementary
ply of a fiber composite made of (0. Allix 1989)
- parallel fibers considered as linearly elastic, anisotropic and nondamageable up
to brittle failure;
- a matrix whose behavior is governed by anisotropic linear elasticity coupled
with damage by progressive cracking in the direction of the fibers;
66 Thermodynamics and Micromechanics of Damage
- an interface which is linearly elastic and damageable in shear and elastic and
brittle in transverse tension perpendicular to the fibers.
Assume plane stress conditions in the ( X l' X 2 ) plane. Take EI as Young's modulus
in the X I fiber direction, E2 as Young's modulus in the transverse x 2 direction, v l2
and v21 as the corresponding Poisson's ratios, and G I2 as the shear modulus. As
the damage does not affect the behavior in the X I direction, must be equal to Dn
in the x 2 direction and equal to Ds in (1,2) shear, then the complementary energy
density must be written as
The plane strain components of the law of elasticity coupled with this particular
damage are derived from this elastic stress potential (0: and (3 being indices taking
the values I and 2):
ow'~
oD,
ow'
1'" --
~
oD"
Demonstrate that the principal directions oj the stresses remain constant with time
(but may possibly be different at each point):
[a] =
SI
[ 0
0
S2
0
0
1 T(t).
o 0 S3
Demonstrate that the triaxiality ratio remains constant with time:
I
SH = :3 SH :'
2I [Sik ( 1- D) kj+
- I (
1-
D) -I S ]
ik kj T(t ) =
E ,.
l+v Cij+( I+v)(I_2v)cu,u;j,
vE e i:
and show that the elastic strain tensor and the strain energy release rate tensor
are also "proportional":
c
C'i j(M.t) = eij(M) T(t ),
y l ee y 2
ij = 2a i lllUCj'ITIEU = ;j(M)' T(t ) .
- Assuming that the damage rate depends linearly on the strain energy density
release rate (see Chapter 3), show that the damaf?e tensor is also' 'proportional"
Jor the initial condition Di j (t = 0) = 0:
The main mechanism of brittle damage is the nucleation and growth of microcracks
and their coalescence to initiate a crack at the mesoscale. Let us derive from the
microcrack pattern the relation between the isotropic damage variable D defined
at the mesoscale, and the surface density of microcracks.
Consider for simplicity a R.Y.E. at the mesoscale as a cube of dimension l l l
and define the microscale of cells as the cube d d . d in which a microcrack of
any area Si and any orientation may lie (Figure 2.11). The number of cells is [3 / d 3
although the number of cracks, is n(n ::; l3 / d3 ).
The geometry having been defined, modeling consists here in writing the balance of
the dissipated energy calculated by classical fracture mechanics and the dissipated
energy calculated by continuous damage mechanics.
/ /
/ /
L
MESO SCALE RVE
~
V
/ ~ / Micro scale cell
- .
-~ Id
\ ~
/ - ~ Microcrack area Sj
/
- /
~ GiSi
" .="~ Y;Di
'3 d .
I I
Assuming that brittle growth of microcracks occurs at G = G e = cons!., and that
brittle damage occurs at Y = 1";: = cons!.,
". 3'"
n n
Furthermore, when s; =
~
I I
or D =
e I
The term Gjy)3 may be calculated from the rupture criterion for mesocrack
initiation in one plane or here in a set of cells occupying the flat volume [ .[ . d.
Due to the localization of damage, other microcracks may be neglected.
If one assumes that the mesocrack is initiated when a part k of this flat volume
has microcracks of size d2 , corresponding to D = Dc at the mesoscale, then
n
LSi = kP
I
This allows us to write D as
n
D =
2:>
I 1 Dc
-[2-T '
which shows that D is the surface density of microcracks multiplied by a correction
factor (here, De/ k) as was mentioned for the definition of the damage variable in
Section 1.2.1.
If the simplest fracture criterion is considered,
"
70 Thennodynamics and Micromechanics of Damage
D=
This result was also found by the purely geometrical considerations described in
Section 1.3.1.
By the way, this calculation gives the order of magnitude of a characteristic length
which permits correlation between fracture mechanics and damage mechanics,
namely l, the size of the R.Y.E.
Taking the simple fracture criterion k = I, Dc = I,
Gc
YP
c
[2 ;
then It ~ G,
Y c
I
for most metallic materials from the table in Section 3.5:
light alloys steels and alloys
0.005 <
- Gr: <- 0.05 MPa.m
2 -< Yc <- 10 MPa
2.5 10- 3 :::;l:::; 5 10- 3 m
and for concrete in tension: G c :::::' 3 10- 5 MPa.m, ~: :::::' 1.5 10- 4 MPa: l
2 10- 1 m.
This shows that the size of the physical RYE must be of the order of millimeters for
metals and of the order of decimeters for concrete as pointed out in Section 1.1.3,
where slightly smaller values were involved.
As a specific example, let us derive the kinetic damage evolution law at the
mesoscale which corresponds to fatigue microcrack growth at the microscale of
Figure 2.11.
With Dc = k, the damage was found to be
I: si
yI .
The surface growth rate of one crack Si may be expressed as a function of the
strain energy release rate of one cell G i by means of the Paris law of fatigue crack
growth.
2.3 Micromechanics of Damage 71
This example shows that the damage rate is an increasing function of the strain
energy density release rate and is quasi-proportiona~ to it for most materials for
which 'T/ '::::: 4. It is also proportional to its rate Y. This fact will be used as
a guideline for the derivation of a general kinetic law of damage evolution in
Chapter 3.
The main mechanisms of ductile damage are the nucleation, growth and coalescence
of microcavities by large local plastic deformations. As in the previous section, let
us derive a relation between the isotropic damage variable D at the mesoscale and
the density of microvoids.
Let us again consider a RYE at the mesoscale as a cube ll l with n cavity cells
of dimension d3 (Figure 2.12).
On this very simple geometry let us again write the balance of energy calculated
from the growth of cavities and from the damage concept.
According to the Gurson model, the porosity at the mesoscale P is equal to the
1
hydrostatic part of the plastic strain E~ = 3E~k due to the growth of voids. For
the geometrical model under consideration,
nd3
P=P'
where n is the number of cavities.
Writing the equality for the rates at the mesoscale,
p = i~
The total power density dissipated at the mesoscale for the homogenized stress aij
and plastic strain rate ifj is
/ /
/ /
/ /
RVE
~ /
Micro cavity
~/
I ~
'/~ ~ I~ I~ Id
~ /
~, /
/
which can be split into two parts by means of the deviatoric and the hydrostatic
quantities:
In the Gurson model, the porosity rate is also the sum of two terms accounting
for nucleation and growth. For nucleation the Tvergaard kinetic law is used:
P = Aaeq + BaH'
where A and B are material parameters.
Assuming for simplicity a sudden nucleation of cavities of a fixed size d,
. d3
P = nr'
. .1 1
D = Pd = d (Aaeq+Ba H)
or b ~ ~ff" (A+B!:)
It is convenient to express the damage as a function of the accumulated plastic
p= aeq .
Er'
(JH (JH
and also
a eq (J eq ,
~ET (A+B::) P
b) Growth by enlargement of cavities of a fixed number n:
. dd
D = 2n r
d2
or with D = n(2'
The problem of crack growth has received much attention in the past 20 years. An
essential result for this problem is offered by the McClintock and Rice and Tracey
analysis, which derives the rate of growth of a cylindrical or a spherical cavity
of volume V in a perfectly plastic infinite body as a function of the accumulated
plastic strain rate p and the triaxiality ratio (J H/ (J eq:
V = 0.85Vpexp (~(JH)
2 (Jeq
or, with V = d3 ,
b = O.57DjJexp (~(JH)
2 (Jeq
.
These two examples of cavity nucleation and growth show that the damage rate is
- proportional to the accumulated plastic strain rate ;
- an increasing function of the triaxiality ratio (J H/ (J eq;
- a function of the state through ET or D.
These properties will also serve as a guideline for the derivation a general kinetic
law of damage evolution in Chapter 3.
Consider a cell at the microscale containing a circular crack defined by Euler's an-
e,
gles if; of the orientation of its plane, by its radius a and its area A (Figure 2.13).
The matrix is isotropic linear elastic and the response is perfectly brittle.
X2
e)
"'-\
[rrl'j
Micro Ell x,
radius a
[
- Crack area A
~j
X3
MESO[ Volume V
Let aM and EM be the stress and elastic strain tensor fields in the cell, decomposed
in
aM = aO +a c ,
EM = EO + EC,
where aO and EO are the stress and elastic strain "far fields", which are constant
in the cell when no crack is involved, and a C and E C are the disturbances in stress
and strain due to the crack. These fields verify the classical equations of continuum
mechanics, assuming a static case and no body forces:
aO .. = 0,
tJ,J
aijnj -- Fi'
Eij
= Soijklakl,
where So is the elastic compliance tensor of the matrix.
The total stress and strain also verify that
E~j + E~j = Sfj kd a k1 + akl),
and the stress disturbance verify that
C . . = 0,
a tJ,}
C
a ij n j -- 0 ,
The important role in the loss of stiffness of the cell due to the crack is played
by the crack opening displacement vector b defined by the discontinuity of the
displacement across the opened crack,
bi = [[u~ + udl.
The relation existing between b and the "far field" stress (70 is somewhat difficult
to obtain. Using Eshelby's result for an ellipsoidal inclusion of vanishing thickness
and Mura's computations, the result with respect to the crack coordinates, where
e~ is normal to the crack in its center, is
1 1
The total stress averaged over the volume of the unit cell is
, IJ
'J
=~
V
(lJ o
'J
+ IJ'Jc. ) dV = 1J'J0 + ~
V
IJc
'J
dV.
V v
To perform the integration by parts let us write the equation together with
xj, k = Djk as
IJ
'J = IJ'oJ + V1- 1 ,. x J,, k
IJk
C dV
v
or
IJ 'J = 1J0
'J +V ~ [I IJcknk
'J
x dA - 1 IJ"" X
'.' ,.k J. dV]
A v
The equilibrium equations for the stress disturbance require the two integrals to
vanish and, as expected; we obtain
o
lJij = lJ ij
The average of the strain may be deduced from the stress and from the contribution
of the crack-opening displacement:
Eij = V1 I( 0
Eij + Eij C )
dV = 0
Sijkl VI I( 0
IJkl + IJkl
C ) *,
dV + Eij
V V
Eij o llJkl + Eij
= sijk *
The determination of cij first requires the calculation of the average of the crack-
opening displacement over the crack area:
1
A
b; dA ;
then, the transformation from local (e;, e2'eJ ) to global (X I ,X2, x3 ) coordinates by
means of the rotation matrix gij:
[Bill
B22 = 8 (1 - v 2 ) [ 22 . 1
B33 1TE(2-v) (2-v)
E and v are Young's modulus and Poisson's ratio, and f is Coulomb's coefficient
of sliding friCtion between the two surfaces of the crack.
78 Thennodynamics and Micromechanics of Damage
The final expression for the homogenized stress strain relation in a single crack
cell loaded in tension is
iknjnl(Jkl
+Oilnjnk(Jkl
+Ojknlni(Jkl
+Ojlnink(Jkl
-2vn i n j n k n l(Jkl
Consider now a set of cells each containing a crack of different orientation. Let us
calculate the strain energy of the whole.
According to Taylor's hypothesis we consider a reasonably dilute concentration of
microcracks which allows us to neglect the interactions between the cracks. The
above analysis may be used considering a random distribution of orientations and
sizes of the crack (Figure 2.14). The strain energy is simply the superposition of
the contribution of each cell.
For one active crack cell, the strain energy density is calculated from
32 (,\ + 2J-l)J-l a3
-- -
3 3'\ +4J-l V
/' - - , \
a'11
- - / /
-
[
-" ~ Micro e It
I / \ /
'- \ I / -
, / - - /'
A. _ vE E
f.1.=
- (I + v) (I - 2v) , 2(1 + v)"
3
The measure of the damage is here ~, its relation with the D variable used all
along the course is obvious.
Considering a cubic cell of length [ : V = [3 ,
and
V
If we now consider a set of cells with a random distribution of cracks, the calcu-
lation of the mean value of the energy requires us to consider the following terms
as random variables:
a3
vn i nj '
a3
V n i njnkn/.
Calling p ( ~ , ii) the probability density associated with the relative void area of
a3
the cracks V and their orientations ii, let us introduce the Leckie and Onat tensors:
v
-Wijkl"2cijCkl
By analogy with the isotropic case, this expression shows that the description of
the anisotropic damage is achieved by three variables:
- a scalar W as in the isotropic case ;
- a second-order tensor w ij as may be obtained by a purely geometric analysis
(see Section 2.2.1);
- a fourth-order tensor wijkl as it is obtained by a pure thermodynamical analysis
(see Section 2.2.2).
/ /
/ /
Crack area s
Damaged zone
Let us write the overall equilibrium equations in the plane of the crack.
s
or with the simple definition of damage D [2 and the principle of strain
equivalence,
I-D
- In compression, P = P - , the damage as a state variable is still D but its effect
is different. The resisting area [2 may be decomposed in two parts: ([2 - s),
for which the elasticity modulus is E; and s, for which the elasticity modulus
is one of the damaged zone with the crack closed E = E(l - DJ. Then the
equilibrium equation is written as
P- = J
l' -s
arM) dS+ J arM) dS
a- = Ec-[l- D + (I - DJD]'
a- = Ec-(I - DDc)
The principle of strain equivalence requires that ij = Ec- and
_ a
a - -:-----..,..
-(I-D.Dc)"
This simple example shows that the crack closure parameter h is of the order of
magnitude of Dc' a quantity which depends upon the material and the loading:
Ih ~ Dc I
In practice, h is considered to be a constant, in order to avoid a tremendous
amount of algebra! A value of h = 0.2 often gives results close to the experi-
ments. It is also possible to identify its value from the measurements of elasticity
modulus in tension and in compression on a damaged material as explained in
Section 3.1.6.
(x ) = x if x 2 0
or with the Macauley brackets { (x) = 0 if x < 0
or u = (u ) - (-u ).
2.3 Micromechanics of Damage 83
00 ..
a'J
(-21((ji) ((jij ) ) = ",0 . .
u(j.'J
(~ ((ji~,)(akj )8ij) = (aij ),
o~ (~ ((jkk)2)
'.J
= o~ . (~ ((jk:k)((jH,))
1.J
= ((jkk )8ij ,
I+ v v I +v
E(l - D) (a ij ) - E(l - D) (akk)8 ij - E(I _ Dh) (-(jij)
v
+ E(l - Dh) (-(jkk )8ij
84 Thennodynamics and Micromechanics of Damage
_ (-akl)+~((-arS)Drs-(-arr))Dkl] )
1 - Dh Dkl Dij
Care should be taken with respect to the order of the operations; for example,
([(arr)DkllDkl)Dij = (a rr )3Dij ,
and also tr(a) :j: (tr(a))
or (akl)b kl :j: (a kk ),
v v
(a ij ) + ~ ((akl)Dkr-<akk) )Dij (-a ij ) + ~(-akl)Dkl - (-akk)D ij
aij 1- D 1- Dh
Let us introduce more convenient notation:
1 .
a H = :3 a kk the hydrostatic tress;
+ _ 1 1
a H - :3 (akl)D kl , aH = :3(-akl)Dk1 ;
a + -_ (
a) + -3v- (a + - (a H)) D .
'J 'J 1 - 2v H 'J'
3v
aij = (-aij ) + 1_2)aH-(-aH))Dij ;
+ -
~-~
1- D 1- Dh
I
-a
3
+
at = 0
I
a H = -30' - (positive)
ali - (-a H ) =0
a
I-Dh
The damage mechanics in unilateral conditions may be developed as it has been
done all through this chapter by simply using the new expression for the effective
stress in the potentials.
+ 2pE(I I - Dh) (I
[
+ v) (-aij )( -ajj ) -
2
v( -akk ) ] ,
-Y = -P aD
eN; = - 2E(1 1_ D)2 [
(I + v) (aij )(aij ) - v(akk )
2]
2E(1 - D) '
86 Thennodynamics and Micromechanics of Damage
Writing the equality with the expression of w; for the three-dimensional case,
eJ* = {(l + lI)(eJij )(eJ ij ) - lI ( eJkk ) 2 +
I D } 1/2
+ I ~Dh [(I + 1I)(-eJij ) (-a;) -1I(-akk)2]
_ ( 1- D ) 1/2 _
for compression (CJ -), eJ* - I _ Dh leJ I
where H is the tangent stiffness tensor, a function of the state of stress or elastic
and plastic strains and of damage as a softening process.
Since this is a problem of continuum mechanics, the displacement and the stress
vector are continuous even across the surface (n). With the notation of the discon-
tinuity of a quantity X across a surface,
[llill = 0,
[aijnJ = O.
Using Maxwell's compatibility equation for the general shape of the discontinuity
of a second order tensor gradient, = ~ [grad ii + (grad ii)T], which must be
written as
. I
[cijll = 2 (ginj +nigj ) ,
where is a vector to be determined:
- = 0, no discontinuity;
- i= 0, sufficient condition for the discontinuity to exist.
Let us see for which property this corresponds to H, starting with
Hkl
'J = H '--lk
J -' H 'J kl = H JI-
- kl ,
then, H i jklgknl = H i jklnk gl
1det( n . H . n) = 0 I
This is a sufficient condition for the existence of a surface of discontinuity of strain,
which can be used for damage localization if H depends upon D.
If an orientation n can be found in an R.Y.E. such that det(n H n) = 0, then a
mesocrack will be initiated in the plane perpendicular to n.
As an exercise, let us work out the simple tension of the plane sheet problem
depicted in Figure 2.17:
a
[a] = [ 0 0 0 ,
0 0] o
-v *c; o ].
000 o -v* c;
88 Thennodynamics and Micromechanics of Damage
The constitutive equations of the material do not need to be introduced. Let us only
consider a constant contraction coefficient v* and the general law
a= H""i-v*i(H1l22+H'133)'
0= H2211i-v*i(H2222+H2233)'
0= H3311i-v*i(H3322+H3333).
The components of the matrix [H] may be calculated analytically from the consti-
tutive equations of each material.
Nevertheless, this particular case may be solved by considering the discontinuities
of the strain rates. The shear strain '2 being zero in the tension problem con-
sidered, the discontinuity of its rate is assumed to be zero also. This corresponds
to a particular mode, the principal one, among several others also possible. From
Maxwell's compatibility equation,
. I
[[E'2 ]] = 2" (g,n2 +n,92) = o.
For the problem in planes (x l ,X2) or (Xllx3)' the components of the unit normal
vector ii to a possible discontinuity line defined by its angle () with the x I axis are
n [:~~l.
g is a vector in the (x" x 2 ) plane; let us call it
9 [~:~~l
Then [e 12] = ~ (-y sin X cos () + ')' sin () cos X) = 0
and we know from the strain matrix in tension (Figure 2.17) that E22
Assuming this relation is still valid locally,
This leads to
cos 2 B = v* sm
. 2 B or tan B = ~- .
v*
This important result shows that the orientation of the mesocrack that will be
initiated does not depend upon the entire constitutive equation, but only upon the
contraction coefficient.
Two cases are of practical interest:
B = Arc tan v-rs = 54~ is the classic angle of the Luder's bands which may
be easily observed by experiments;
* v * = v, Poisson's ratio, corresponding to elastic damageable materials:
v = 0.2 is the value for concrete: B = +65.9~ ;
v = 0.3 is the value for cast iron B = +61.3~;
v = 0.4 is the value for plexiglass B = + 57 . 7~.
The stress strain condition for damage localization is deduced from the equilib-
rium condition, namely, no discontinuity upon the stress vector, together with the
constitutive equations:
This localization condition corresponds for this tension problem to the peak of the
stress strain curve:
dO"II = 0.
d C II
This result is a necessary and sufficient condition for damage localization; it does
not give the condition for crack initiation, which may occur later by the rupture
criterion D = Dc'
To summarize, the mesocrack:
- is initiated for D = Dc;
- growths in the direction normal to n defined by det( n. H . n) = 0;
- has a size of the order of magnitude of the RYE, also determined by the corre-
lation between damage mechanics and fracture mechanics, l = ~e .
e
2.3.6 Exercise on the Fiber Bundle System (1 . HuIt 1982, Chrzanowski 1986)
,!//////// / ////P///,ij-
Cross section
A
~ ~
i~ I ~~.
a
~
h b
from EE =a -+D=EE(I-a)-O'R(I-b)
(x= D) R(x=D) 2(0' Rb - Eca) .
Determine the relation between the load F and the damage D in correspondence
with the three particular cases given in Figure 2.19:
F = 7J A(I-D+aD-aD 2)(I-b+2bD).
R l -a +2aD
The first case is interesting because it shows the possibility of gradual rupture from
F = Fo until of/aD = 0 in a way somewhat similar to ductile damage failure
(but here without any plastic strain).
From of/aD = 0, deduce the following equation from which the smallest root
gives Dc (This is a good exercise for practicing algebra!):
8a 2bD3 + (2a 2 + IOab - 12a2b)D2 + (2a + 4b - 2a 2 - IOab + 6a 2b)D +
+ I - 3b + a2 + 2ab - a2b = O.
Note that Dc is much smaller than I.
Creep rupture
The fibers are now considered to be linearly viscoplastic, with the constitutive
equation:
a( x, t)
i (t) = K N(x ) where K N = K N [I + c(2x -I)] .
92 Thennodynamics and Micromechanics of Damage
2
I+a <b<1
3- 2a+ a 2
F F F
o o
o Dc o o
Progressive Instantaneous brittle rupture
ductile rupture
Fig. 2.19. Motion of the rupture damage front D in the fiber bundle model
The elastic strain is neglected but the ultimate tensile stress is still
aR = aR[1 +b(2x-I)]'
with o< b < I, - 1 ::; c ::; 1 and b > c.
As the creep strain progresses, the volume remains constant.
F F
Determine the true stress A as a function of the engineering stress - where Ao
Ao
is the initial cross-section area,
dL F F
dE = _. LA = L A ---+ - = -expE.
L ' 0 0 A Ao
The load F is constant ; it induces creep and failure of the fibers. Determine
the creep rate E: as a function of time before failure, with the initial condition
t = 0 ---+ E = O.
From the equilibrium equation
J
I
F
a(x ,t)dx ---+E:
A
o
and by integration,
it follows that
2.3 Micromechanics of Damage 93
* if
F I-b . f' .
* if -A ::; aR-- -> to IS mIte,
o l- c
there is delayed rupture of an incubation time to > O. These two cases are
represented in Figure 2.20.
The failure propagation is governed by two equations:
J
I
* equilibrium a (x, t ) dx = :
D
Jc\
I
F
or t) K N [I + c(2x - 1) J dx = Aexp C:(t)
()
Deduce the following relation between the damage D( t) and the creep strain i :
D=O { i= ~exp c: }
AoK N
.
-aR I - b
c: = = - - - -
b = (l- c)(I-b) aR .
2(b- c) - (l- c)2((I-b) K N
KN I- c
Note that the damage rate does not depend upon the load.
94 Thermodynamics and Micromechanics of Damage
D D
.E..~o: l-b ...~o: l-b
AD Rl-e AD R l-c
1+c 2 1+e 2
b~ 3-2c+c 2 b> 3-2c+c 2
Immediate Immediate
instantaneous gradual
failure failure
D D
Delayed Delayed
instantaneous gradual Fig. 2.20. Creep rupture of the
failure failure fiber bundle model
if b::;
3 - 2c + c2 '
* gradual development of the damage:
iJ finite,
1+ c2
if
b > 3 _ 2c+ c2'
These two cases appear in Figure 2.20.
Courageous readers can also calculate the time to rupture and the critical value
of the damage x = Dc which produces the final instantaneous failure (some
approximations must be introduced)!
Chapter 3
Kinetic Laws of Damage Evolution
The physical nature of damage has been defined in the first chapter. The defini-
tion of the damage variable as the effective surface density of microcracks in a
Representative Volume Element associated with the effective stress concept and
the principle of equivalence has given rise to methods of damage measurement
through changes in elasticity or plasticity. These concepts have been generalized
to the three-dimensional case in the second chapter by means of two potentials.
Accordingly, state coupling occurs between elastic strain and damage, and kinetic
coupling takes place between plastic strain and damage, which allowing us to cal-
culate strains and damage up to failure if the constitutive equations for the damage
are known.
Many models have been proposed in the past; their basic properties are contained
in a unified formulation of the potential of dissipation, which will be studied in
detail in this chapter. After the kinetic law of isotropic damage is derived, it is
applied to several classical cases of loading giving rise to different kinds of dam-
age such as brittle, quasi-brittle, ductile, and low cycle fatigue or high cycle fa-
tigue, with special emphasis on the particular properties of metals, ceramics, poly-
mers, composites and concrete. The common main feature is the proportionality
of the damage rate to the strain energy density release rate and to the accumu-
lated plastic strain rate beyond a plastic strain threshold and up to a critical value
of the damage variable. The accumulated pla~ic strain which governs the dam-
age is defined on the meso R.V.E. or at the microscale when the damage is very
localized. The case of perfect plasticity is of great interest for further develop-
ments.
D ijD+FD (-
F = ((7-D -X D) f'q-R-a y + 4X3 XiJX Y; (r,D )) ,
ex:
with no analytical expression for the damage term F D from which the damage rate
b is derived
b = _BF~ = BFD~
BY BY'
with ~ = r= zi(\ - D),
96 Kinetic Laws of Damage Evolution
and
The choice of this function F D is of course the key to representing the damage
evolution. It is a matter of logic through the "State Kinetic Coupling theory", of
physical observations, of experimental data and of micromechanical mechanism
modeling to quantify the effects. Like Hooke 's law for elasticity or the Prandtl-
Reuss law for plasticity, a unique kinetic law of damage evolution may be de-
veloped to model the general trends of all kinds of damage in many kinds of
materials.
D= 0 if p < PD'
This allows us to introduce a step function in the potential F D which has to be
a function of the state variable r taken as a parameter:
I if r 2:: PD
Oifr < PD
In monotonic uniaxial loading, PD may be identified with the uniaxial damage
threshold EpD but in fatigue or creep processes, PD is a function of the applied
stress as will be established in Section 3.1.2.
- On the basis of a thermodynamical analysis, the main causal variable for the
damage is the strain energy density release rate Y, as it is the dual variable of
3.1 Unified Fonnulation of Damage Laws, O. Lemaitre 1987) 97
or
I'D-'-=-~-p- I,
. H-(-P--P-D-")
with the rupture condition for crack initiation,
D = Dc'
- Three material parameters are introduced to characterize the damage evolution:
* S --t the energy strength of damage;
* p D --t the damage threshold function of the material and the loading (see
Section 3.1.2.)
* Dc --t the critical damage, a function of the material and the loading (see
Section 3.1.3).
The effects of the temperature T are taken into account by the variation of these
coefficients with T and by the accumulated plastic strain rate p, also a function of
the temperature.
Several properties not directly introduced in the formulation will be observed as
consequences:
- the nonlinear accumulation of damage;
- the effect of mean stress in fatigue;
- the nonlinear interaction of different kinds of damage;
- the brittleness of materials as the applied stress or the triaxiality increase.
strain hardenings and no damage (D =
-+. --
'I' (JijCij P - X ijD a' ij >
.p - R' _ ,
which may be calculated from :
F = (uD-XD)eq-R-(Jy+_3_xBxi~;
4Xoo
- its associated normality flow rule,
.p of . of
Cij = o(J . A, P = -oRA A,
'J
- and the yield criterion (Section 2.1.3)
f = (u D - XD) eq - R - (Jy = 0,
oF of D of ) .
= ( (Jij 0(J . + R oR + Xij oXD P 2 0,
'J 'J
J J
t t
t DDt
j ~2 CY CY pdt - j CY ypdt
CYijij
eq
o 0
Ws = j(CYeq -CYy)dp.
o
This fonnula may be used; however, it is still easier to assume a perfectly plastic
material of plastic threshold CY eq 2: CY y' whose stored energy is a function of
the difference between CYeq and the fatigue limit CY j' that is the stress, still in the
conventional elastic range below which no damage may occur. It is identified as
the stress amplitude which corresponds to a very large number of cycles to failure
in a tension-compression fatigue test: NRc::,- 106 to 107 cycles. Some values are
given in the chart in Section 3.5.
Then, as CYeq = const.
Ws = (CY eq - CYj)p.
The characteristic value of that energy corresponding to microcrack initiation
(D = 0, b > 0) is taken as that of the pure tension reference case having a
damage threshold C PD and a plastic threshold CY u' the ultimate stress.
(CY u -CYj)CPD'
Writing the equality of these two energies gives the value of the damage threshold
PD for any kind of loading:
(CY eq -CYj)PD = (CY u -CYj)CPD'
ues for several materials may be found in the table in Section 3.5). We then
3.1 Unified Formulation of Damage Laws, (J. Lemaitre 1987) 101
obtain
.
We aijklci/:kl
e 'e (I -
D) -
1
"2aijklCijCk/
e e D'
Physically, this means that an instability may occur if the energy released by loss
of stiffness becomes equal to or greater than the energy needed for the creation of
new damage surfaces of decohesions.
In fact, during this stage of instability prior to failure, the strain energy density
release rate is almost constant. For example, if a perfectly plastic material of plastic
threshold as is submitted to a proportional loading for which Rv = const. (see
Section 2.2.4),
a eq
- - - as = O,
1- D
a;Rv
2E = const.
The instability may occur at any time during the damage growth when some per-
turbation slightly increases the quantity Y. Then this instability condition cannot
simultaneously be the rupture criterion. From both the practical and theoretical
points of view, a physical condition of rupture must be added.
To this end, let us assume that the final stage of damage rupture by atomic decohe-
sion, (that is mesocrack initiation) is governed by the amount of energy dissipated
in damage growth:
jY
Dc
dD = const. at failure .
o
This quantity may be identified from the uniaxial case in tension, still taken as a
reference:
a eq = a R }
Rv = 1 Y
D = Die
and
In three dimensions, the case of perfect plasticity in proportional loading is, for
simplicity, considered as an approximation:
Dc Dc 2 Dc
a;Rv D
j YdD = j 2E(I-D)2
a eqRv dD = ja;R vdD
2E 2E e
0 0 0
a*
as a-* a* a eq R vI / 2 an d a eq = a .. (1 - D).
I-D'
Then the rupture criterion gives the critical value of the damage as a function of its
value for the one-dimensional reference case Di e' the effective damage equivalent
stress ij* and the ultimate stress a" with the limiting value Dc = I:
ij*2 a2
2EfD" = 2f;;D lc '
Dc = Die -:2
a
a
2
: :; 1
3. 1 Unified Fonnulation of Damage Laws, (J. Lemaitre 1987) 103
This formula shows that the critical value of the damage for mesocrack initiation
decreases as the effective damage equivalent stress a* increases either by the
stress a eq' the damage D or by the triaxiality ratio a H / a eq contained in the Rv
expression.
DI e is the critical value of the damage of the case of reference under tension. Any
other one-dimensional case (like fatigue, for example) may have a different value
of De' depending upon the applied effective stress a.
It is interesting to study the different possible cases in the graph ( Dc, a*) of
2
Figure 3.1, where the rupture relation Dc Die ~
'a
:2 corresponds to a line which
defines:
- a safe domain below which a stress-damage path does not produce crack initi-
ation,
- a rupture domain upon which the material is certainly broken.
Rupture line
Safe
domaine
o 0
2
o o
3 4
Fig. 3.2. Damage rupture paths. I) Fatigue !::!.U const; 2) Fatigue !::!.C: const; 3) Ductile case;
4) creep u = cons!.
This mesocrack initiation is completed by the localization criterion giving the ori-
entation of the crack developed in Section 2.3.5.
The entire set of equations governing damage evolution can be summarized as
follows:
b
o if P < PD
In most cases Y is the strain energy density release rate defined in Section 2.1.1. In
order to take into account the micro-crack closure effect responsible for different
damage behaviors in tension and compression, the only change to be made is in
the definition of Y derived in Section 2.3.4 as a function of the crack closure
parameter h.
3.1 Unified Fonnulation of Damage Laws, (1. Lemaitre 1987) 105
P > PD
Taking 0'8 = O'u ' the maximum value of the stress on a stress strain curve in
tension, where dO' = 0, is often a good choice (Figure 3.3 for example).
dE
Let us write the corresponding kinetic laws for plastic strain and damage taken
from sections 2.1.4 and I. I of the present chapter.
of .
-A if {f = I ~~ - 0'8 =0
OO'ij j = 0 with ~ = p(1 - D),
O';qRv .
ifp > PD with Y = 2E(I_D)21
aeq +0'8D = 0;
. 0'2
then: D = 2IlSRvp,
106 Kinetic Laws of Damage Evolution
.p _ 3ES a iJ
Cij - 71-DR'
s v
This shows that for perfectly plastic behavior when P > PD' the plastic strain is
directly related to the damage.
Another interesting result coming from the consistency condition is
iJ = _ a eq .
as
Notice that the equivalent stress rate is negative, owing to the softening damage
process.
The last interesting result comes from the plasticity criterion written in tenns of
the elastic strain: the law of elasticity coupled to the damage written in tenns of
the deviatoric and the hydrostatic stress and strain is
e eD + cHe fj ij' c eH
I e
Cij = Cij 3' c kk '
I
a i j = aB+aHfjij' aH 3' a kk '
c eD I +v D
I
1-2v
c eH
tJ E(I_D)a ij , E(I_D)a H
It follows that
and
or, with f = 0,
eD e D 1/ 2 _
(Cij Ci j ) - V~3' (I +E)
v \
a 8
_
- const .,
which makes the second invariant of the elastic strain deviator constant for any
state of stress verifying the plasticity criterion.
In tension, where .
a eq
I-D
EC e - a s = 0,
as
Ce = E = const.
This property is shown schematically in Figure 3.3 for a classical tension test and
for a cyclic tension-compression test, strain controlled at constant strain amplitude.
It proves that the measure of the damage by the variation of the elasticity modulus
and by the stress amplitude drop are equivalent.
3.1 Unified Formulation of Damage Laws. O. Lemaitre 1987) 107
I+ v_ v_
E aij - Eakkti;j '
-D
3 a ij . . _
.p
Eij --P If a eq -a~ = 0,
2 as
2
D= 2isRj) if P 2: PD'
In order to take into account the micro crack closure effect explained in Sec-
tion 2.3.4, the assumption of perfect unilateral conditions may be made; it corre-
sponds to a zero value of the crack closure parameter:
h = O.
The strain energy density release rate from Section 2.3.4 reduces to
I +v < aij >< aij > v < au: >2
y = 2E (I-D)2 -2E (I-D)2
and, with this expression, the kinetic law of damage remains unchanged:
. y
D = -p if P 2: PD
S
a b
Fig. 3.3. Constant elastic strain. a tension test; b cyclic tension compression test
108 Kinetic Laws of Damage Evolution
Viscoplasticity
The case of nonhardening behavior reduces the set of elastic perfect viscoplasticity
constitutive equations from those of Section 2. 104 to
Cij = e
Cij
+ Cij'
p
.p
c1,)
.
2
2 3 ( ECH )
with -(I+v)+-- -_-
3 1- 2v a eq
Damage threshold
au -af
PD = cPD
a eq - af
For perfectly plastic materials this becomes
and
Rupture criterion
Then
and
D .
== --1. If a == au.
Rv S
3.1 Unified Formulation of Damage Laws, (1. Lemaitre 1987) 109
The detennination of the coefficients 5, [PD' and Di e that characterize the damage,
together with a f and au related to plasticity, must be worked out for each material
and temperature, from experiments perfonned on those cases in which the damage
is easiest to measure. This involves:
. a2
D = 2E5(1 - D)2 Ii p IH(" - PD) .
- [p D is the plastic strain below which there is no appreciable damage (see
Figure 3.4).
- au is the classical ultimate stress as shown in Figure 3.4. It is also interesting
to know the yield stress a y as it may help to make a better choice of the plastic
threshold as when needed.
- a f is the fatigue limit taken, as already mentioned in Section 3.1 .2, as the
stress amplitude corresponding to a number of cycles to failure in a tension
O'c - -- - - - - - - -
dO
Ep
o~----~----------~
- The last parameter S is determined from the slope of the curve: damage D
versus the plastic strain cp :
0"2
D- .
- 2ES(1 - D)2c p
or
dD 0"2
dc p 2ES(1 - D)2'
At each point of the curve, D is known, 0" is known from the stress strain curve,
dD / dc p is estimated and E is known from a previous identification:
0"2
S-
- 2E( I - D)2 dD .
dc p
Several points may be considered in order to obtain S as the best average.
The main difficulty involved in this identification lies in obtaining a good stress
strain curve in the softening range where necking occurs. To avoid any instability, a
"soft" machine is required with an excellent feedback system, which can accurately
impose the strain in the damaged region of the specimen. This strain must be
measured locally by a small strain gauge (see Section 1.3).
with -
0" = -0"- un d er tension
.
I-D
and I _0"Dh under compression (see Section 2.3.4),
with the condition for crack initiation being,
2
D Ie ~
0-*2
=D Ie'
O"eq
as I and 1- D = O"u
N cycles
0-10
200
iif!~f 2000
2300
2400
2500
2550
Fig. 3.5_ Low cycle fatigue stress strain pattern, for AISI 316 L stainless steel (after J. Dufailly)
h = 1 - ~- / E = 0.2
\- E+/E '
.tE
1.0 t---~=======::::::~~~=~
_______ E- -
0.8 '\
E'
0.6
0.4
0.2
N
taken as the mean value for several values of N. Also deduce the damage as a
function of N by
t+ E-]
D = 1-- or D = h
I [ 1- .
E'
P = 2 LAcPi'
i=1
PD = Max(p(D=O)) :::: 6,
Dc = Max(D(p)) :::: 0.15.
3) Determine the parameter S as the average value of several computations such
as:
S -- 0"7w
dD ' S:::: 7M Pa
2E(I-D)2-
dp
o
1.0 1 - - - - - - - - - - -_ _- ,
0.8
0.6
0.4
I
0.2 D1c =0.15 ,I
r -- ---~
p
o 2 4 6 8 10
Fig. 3.7. Damage evolution for a low cycle fatigue, on AISI 316 stainless steel
3.2 Brittle Damage of Metals, Ceramics, Composites and Concrete 113
Elastic (plastic)
E(o'y' o'u.o'f)
Furthennore, the weakness also comes from the fatigue limit aj assumed to be
reduced in the same proportion as the plastic threshold:
a~
J1
at = at-
ay
The complete inclusion problem will be solved numerically in Section 4.3. by "lo-
cally coupled analysis". Here, however, some approximations allow us to derive
the rupture conditions without resolution of the complete set of constitutive equa-
tions. Starting with the kinetic damage law for the inclusion,
. yJ1
D - - _ s p' J1 ,
we wish to express yJ1 and pJ1 as functions of macroscopic quantities such as the
elastic strain and the stress (T.
According to the Lin-Taylor hypothesis, we may assume that the inclusion is sub-
jected to the state of strain (or strain rate) of the matrix, which is taken to be
unifonn:
Neglecting the elastic strain e J1 in comparison to the plastic strain PJ1 in the
inclusion allows us to write:
P
'J1 _
- (
32 Cij
'PJ1 'PJ1)
Cij
1/2 _
-
(23 .J1D .J1D
Cij Cij
)
1/ 2
The inclusion being perfectly plastic, then, from the yield criterion,
all
eq _ II
l-D -as'
The triaxiality factor may also be expressed as a function of the stress at the
mesoscale:
R~ 2
3(I+v)+3(1-2v) ( ;~J1)2
a~q (I - D)a~;
3.2 Brittle Damage of Metals, Ceramics, Composites and Concrete 115
cH -ea
1-2v
H from pure elasticity at the mesoscale;
aJl a
then a~ = (I - D)a H and ff
a eq
~.
as
au-a! au-a!
The damage threshold PD = cP D becomes PD = cP D Jl Jl because
a eq - a! as - a!
for the inclusion, a eq = a~ and a! = aj, the fatigue limit remaining a! for the
case of reference.
It is not possible to write the yield criterion exactly in terms of macroscopic
quantities. The following approximation is made:
D= 0 if (J eq < (J J .
with
c - co R +- D - D - D __ a u_
2
Ceq - Ceq - " - Ie Jl2 R .
as v
This last condition is, strictly speaking, a crack initiation at the microscale but can
be shown to also correspond to the brittle fracture of the matrix ! The strain energy
release rate at the mesoscale corresponding to a microcrack of surface d 2 is
G = _~ 8WI
2 8A <1=consl.
I
But -28W1<1=consl is also equal to the energy dissipated in the inclusion by the
damaging process:
8A
Assuming a constant strain energy density release rate Y Yr and a critical value
of the damage Dr = I yields
G(DA) = Y,d.
Using the result of Section 2.3.1 in which a simple relation was found between the
critical value of the strain energy density release rate Y" and the dimension of the
116 Kinetic Laws of Damage Evolution
In practical problems of engineering the triaxiality ratio varies from 0 (pure shear)
to 5 or 6 at the tips of very sharp notches, a H = ~ being the value for pure
a eq 3
tension.
For pure brittle damage , plot the von Mises equivalent stress at failure divided
by the one-dimensional damage equivalent stress as a function of the triaxiality
ratio. Take Poisson's ratio as a parameter: v = 0.2, v = 0.3, v = 0.4
(Figure 3.9).
0.8
0.6
0.4
0.2
o 1/3 2 3 4 5 6
Recalling Section 1.1.4, the damage process is called ductile or creep when it
occurs simultaneously with large permanent strains. At the microscale it mainly
involves the nucleation and growth of cavities causing ductile damage in metals
and polymers; in metals submitted to elevated temperature, it causes creep damage
on the form of grain boundary sliding and decohesion.
This is a particular case of loading which often occurs in structures for which
integration is simple.
118 Kinetic Laws of Damage Evolution
R
Ceg- PD =RV(I)
CR-CPD Rv(~
V=O.2
o 2 3 4 5 6
Fig. 3.10. Strain criterion for quasi-brittle failure
The following results were obtained in Section 2.2.4. Starting with the definition
of a proportional loading in a structure:
it follows that
O"H
const .,
O"eq
Rv = const.
Assuming as in Section 3.1.4, that the ductile damage occurs when the strain-
hardening is saturated or almost saturated, then from the plasticity criterion:
3.3 Ductile and Creep Damage of Metals and Polymers 119
This approximation makes the strain energy density release rate a constant:
(J~qRv (J2 R
Y = 2E (1 _ D)2 ~Ev = const.
and the integration of the damage constitutive equation becomes obvious:
J~~~
t
~'E~ Jdp;
2 P
D = p(t)H(p - PD)
o Po
This makes it possible to determine the value of the damage at each point of a
structure as a function of the accumulated plastic strain, provided that the loading
is proportional.
In one dimension and monotonic loading:
Rv = I,
p = Ep,
(J2
Dc
.
= Dl c ~
(Js
;
(J ;
2ES(cPR - cpo )'
This allows us to write the one-dimensional evolution of damage as
E -E
D = D P PD,
cE - E
PR Po
which is the equation of a straight line lying between the two points (D 0,
Ep = Ep o ) and (D = D r.,E p = E p R ) as shown in Figure 3.11.
Dc - - - - - - - - - - - - -
Figure 3.12, which is similar to Figure 3.10. The one-dimensional case is of course
represented by the point
These limit curves depend upon the type of material by means of PD and slightly
by means of Poisson's ratio if 0.25 :S v :S 0.33. They strongly depend upon
the triaxiality ratio. The ductility measured by the maximum accumulated plastic
strain at crack initiation decreases substantially with the triaxiality ratio; it is well
known that a "sharp notch makes the material more brittle". (In fact it is not the
material itself but the effect of the hydrostatic loading on it).
3.3 Ductile and Creep Damage of Metals and Polymers 121
rr
- neglecting the kinematic hardening in the plastic constitutive equation.
1.2
0.8
0.6
0.4
0.2
o 1/3 2 3 4 5 6
The plane stress condition imposes for the principal stress components that:
o
a2
0]
0 ,
o 0
from which the plastic incompressible condition is written for tr.([c]) = 0:
101 ~ [~ ; -(o,~+J
then
2 (2 2 )\/2
PR = J3 c\ +c2 +c\c2 .
The triaxiality ratio a H / a eq may also be calculated as a function of the strain from
the plastic constitutive equations.
From Section 2.1.4,
-D D) A.
. p _ 3 (aij - X ij
Cij - 2 (u D - XD) eq 1- D'
Assuming no kinematic hardening (X = 0) and proportional loading, it is straight-
forward to demonstrate the proportionality between strains and components of the
deviatoric stress:
A plot of e21 eRas a function of elle R is the classical failure limit in deep drawing
of sheets (Figure 3.14). As for the master curve of ductile fracture, the curves
depend slightly upon Poisson's ratio but the result for each particular case depends
strongly upon the type of material by means of its characteristic e R'
The application of the kinetic damage law to creep damage is also straightforward.
It differs from the case of ductile damage only in that the accumulated plastic strain
rate comes from a visco plastic constitutive equation:
in which Y may possibly take into account the microcrack closure effect (see
Section 3.1.4) Remember that the material parameter 5 depends upon temperature.
From Section 2.1.4,
. - L 1-
(iT D - XD )e q - R - (J
y
]-n
P- n[ K '
x
with p= I-D
~0.8eR
"" \ '
\ 0.6 v = 0.33
" 0.54 v = 0.27
0.45 \ 0.4 /
0.64
'" \ /
"" 0.2
,,\ '\ /
/
/1
./
//
///
/./
Fig. 3.14. Limit curves of deep
- 0.1 o 0.1 drawing
124 Kinetic Laws of Damage Evolution
For simplicity, let us consider a simple perfect visco plastic Norton's law with
neither strain hardening nor yield stress (see Section 2.1.3). Then the visco plastic
multiplier reduces to:
~ . [ (Je q ] N
I- D = P= K" (I - D)
where K" and N are material parameters which can be found for some materials
from the table in Section 3.5. Then
. (J;qRV [ (Je q ] N
D = 2ES(1 _ D)2 K,,(I _ D) H(p - PD)
(J~+2Rv
or D = 2ESKt'(l-D)N+2H(P-PD).
The damage becomes a time dependent phenomenon like the visco plastic strain.
In the one-dimensional case this constitutive equation reduces to the Kachanov's
earlier model:
b =[ (J ] N +2 H( )
A(l- D) P - PD ,
I
with A = (2ESK::) N+2.
Let us calculate the evolution of the damage as a function of the time in a simple
creep process where (J = const.
Let t* be the time needed to reach cp = PD by creep without any damage,
which leads to
The time t* must be calculated from a visco plastic model. For simplicity let us
again take the Norton's perfect visco plastic law:
or
7 (;y j
o
dE,
0
dt , t' = PD (;J- N
(J
E =
e E'
Ep [Kv(I(J_D)] N
D [A(I~D)] N+2 H(Ep-PD) ;
For Ep < PD and with the initial condition t = 0 ~ Ep = 0:
D = 0 ~ E
p
= (~)
Kv N t.'
o
or
The shape of the graph of the total strain C = c e + cp as a function of time for
different values of (J' is given in Figure 3.16.
Fig. 3.16. Evolution of the strain up to rupture in creep processes, 0"3 < 0"2 < 0"1
3.4 Fatigue Damage 127
For Ep < PD the relation between the stress, the plastic strain and the time
calculated in Section 3.3.3 is
Ep = (; )N tor t = Ep (.!!..-)-N
v Kv
Then, it is easy to plot t as a function of (J taking Ep as a parameter (Figure 3.17).
For Ep 2: PD' the relation between the stress, the damage and the time is:
D 1-[I-(N+3)mN+2(t-t'r~3
with t* = PD (;,J- N
PD I - (I _D)N+3
(;J +
or t = --=---=N N +2 .
(N+3)(~)
Taking D as a parameter one may plot t as a function of (J (Figure 3.17).
At rupture, the critical value of the damage Dc is
Dc = I,
which gives the time to rupture to be plotted on the same graph (Figure 3.17).
STRESS (MPa)
700
~~ - D=O.Ol
""'" ~"
.~ - D=0.02
- D=0.05
~~ ~ ::::; ... ~ "!!o.. - D=O.l
......... ~ "I
~ r....~ . - D=0.2
:::::r-. ~ ........
~
~ .;:<~ - D=l
r'" .........;;
::::
.
~ ~~ . t:--. ~ ~~
~ f:::r-.r- ~~, ~
.....
- Ep=O.ool ~ ::::~ tS:; r:::
r-...
)D
- Ep=0.002
~
i""
r...r--
~ .....-...,
- Ep=O.005
~~
100 -
-
-
Ep=O.Ol
Ep=O.02 ~ ............. ) Ep
80 I I I II Ill!
10 100 1000 10000 100000
TIME (hours)
Let us start with the general kinetic damage law written for the one-dimensional
case defined by the stress a(t) and the strain c(t) as functions of time t:
. y
D = Sp H(p-pv),
a~qRv a2
y = 2E(I-D)2 - 2E(I-D)2'
P= Ic'pl,
2
~c = const . ~ ~a = const. ,
~cp = const.
8p
8N
J
I cycle
lipl dt = 2~cp'
Then
au - at
No = -PD- wIth
'
PD = cpD----"-
2~cp' a M -at
The calculation of N D requires the integration of the kinetic damage law over one
cycle and then over N cycles:
0'2
D 2ES(1 _ D)2 lipl
- The material is first assumed to be perfectly plastic during one cycle with a
threshold as = ~a, and the variation of the damage is neglected for the
integration over one cycle (Figure 3.19).
- The strain and damage behaviors are assumed to be identical both under tension
and compression, this corresponds to the crack closure parameter h being equal
to I (Section 2.3.4). This allows us to calculate the damage per cycle as:
A cyclic relationship between flc p and fla is assumed. It can be derived numer-
ically from the equations of Section 2.1.4 and has strong experimental support:
Kp and M are material parameters which can be found for some materials in
the table in Section 3.5.
It follows that:
oD K~ flc(M+2)/M
oN 4ES p ,
D.Ep being constant throughout the damage process, the damage is linear with N:
D = J
N
8N 4ES p
2
8D 8N = Kp D.E(M+2) /M(N - N ).
0
N"
D =
c
K~
4ES
D.E(M+2)/M (N
p R
_~)
2D.E '
P
4ESDc A - (M+2) /M
and N R - No = K2 UEp .
P
Fig. 3.20. Manson-Coffin curve of low cycle fatigue of a strainhardened steel: AISI 10 10 at temperature
T = 20C
132 Kinetic Laws of Damage Evolution
o 0
~-------------- ~
I
I
1
I
I N -NO N
1---
=---------______ ~-~
- L_ _
~
o 0
Then in practice a single curve exists in the diagram (D, N - No ), but the
NR-No
damage evolution curves differ by No in the (D, ~) diagram (Figure 3.21).
An example of real low cycle damage evolution measured by the variation of the
elasticity modulus is given in Figure 3.22.
where:
- N R is the number of cycles to failure of the creep fatigue process corresponding
to some value of Dc;
3.4 Fatigue Damage 133
1.0 D
0.8
0.6
A -2
) tiE. = 1.8,10
-2
0.4 ,. M:=!,IO
" -2
I, 6.E. = 0.3,10
0.2
N
NR
Fig. 3.22. Damage evolution in low cycle fatigue of strain-hardened AISI 1010 steel
102 < N R < l<f cycles
(J' (J'
B C B C
~eFp .le~
A 0 ep
F E
~ep
- N'R is the number of cycles to failure of a pure cyclic creep process at constant
stress amplitude a = aM for which the critical value of damage is D~ ;
- Nt; is the number of cycles to failure of a pure fatigue process without any
hold time,
NR - ~cp
_ -I [PD
2 +
ES [1 - (1 - Dc?]
2
1
'
3a M
2) The cycle in Figure 3.23 is partitioned in
- (AB + DE) considered as a pure fatigue process of plastic strain amplitude
~cF ,
p'
-(Be + EF) considered as a pure cyclic creep process of plastic strain ampli-
tude ~c~ ,
The damage per cycle of the pure fatigue process is
8D a~ F
8N = ES(1 _ D)2 ~cp ,
8D a~ ~ c
8N = ES(1 _ D)2 c p
and
It follows that
3.4 Fatigue Damage 135
N N
The expression ~ + ---.!i may be smaller or greater than I depending upon
NR N'R
the relative values of Dc' D; and D~ . This property is known as the nonlinear
creep fatigue interaction. As 6.c p 6.c: + 6. c~, the sum of the two ratios is
equal to I if Dc = D; = D~ :
NR NR
N F + NC = I,
R R
which is the Taira rule of linear creep fatigue interaction.
At very high temperature, the superalloys show a strong nonlinearity in their creep
fatigue interaction in the sense that a small amount of creep may considerably re-
duce the number of cycles to failure as shown in the bottom curve of the interaction
diagram in Figure 3.24.
If the amplitude of the loading is low, the amplitude of the plastic strain may be
very small, even negligible at the mesoscale in comparison with the amplitude of
the elastic strain. This corresponds to high values of the number of cycles to failure .
For instance,
N R > 100000 cycles.
Another feature which makes the damage analysis of high cycle fatigue difficult is
its high degree of localization. Very often only a very small microelement damaged
at the free surface of the body gives rise to one microcrack by slips (stage I) which
later propagates perpendicularly to the loading (stage 2). Then, for high cycle fa-
tigue, materials may be considered as quasi-brittle and modeled by a damageable
microinclusion embedded in an elastic mesoelement. Its complete numerical anal-
ysis, including high cycle fatigue, is performed in Section 4.3. Nevertheless, the
linear interaction
main features of high cycle fatigue may be derived from the simplified analysis of
the quasi-brittle damage presented in Section 3.2.2.
The following kinetic law of damage evolution
2
. af .
D = 2ESRv (~reqH(ceq - PD) WIth PD =
r1 r1
o~ _ _ _ .'
and
Then
A first integration over one cycle gives the damage per cycle.
da,
JD =2a7
-
IN
-
E2 S
[ (l+v) (!la
--af) +--
2
1 - 2 V [(!la)3
9a7
- -af
2
3]]
A second integration as a function of the number of cycles gives the evolution of
the damage for the initial condition,
where No is the number of cycles for Ceq to reach the threshold PD:
N = PD
o 2.:lc
or
The graph of the amplitude of stress .:la as a function of the number of cycles to
failure is known as the Woehler curve. An example is given in Figure 3.26.
Using the expression for (N R - No) ' it is possible to express the evolution of the
damage in a simpler way:.
138 Kinetic Laws of Damage Evolution
M(MPa)
1200
+
800
+
400
( cycles)
Fig. 3.26. Woehler curve of high cycle fatigue of AISI 316 stainless steel
This shows that the damage is a linear function of the number of cycles. With
the approximation De = Di e' a unit straight line represents the evolution of the
damage in the graph (D, ~--~o)' but as No decreases with the amplitude of
the stress, different lines represent the evolution of the damage as a function of
the number of cycles. These types of damage evolution are shown schematically
in Figure 3.27.
An example of real evolution measured by means of elasticity changes is shown
in Figure 3.28.
o 0
: N -No N
: NR-N o NR
o ~------------------~I--~ 0
Fig. 3.27. Evolution of high cycle fatigue damage as modeled by the simplified theory
3.4 Fatigue Damage 139
1M =600MPa
0.2
/ ,1'.0'=500
/ ,
0.1 / 1/1'.0'=480
/ /I N
Fig. 3.28. Fatigue damage evolu-
~ tion for AISI 316 stainless steel
o 0.2 0.4 0.6 0.8 1.0 104 < N R < 106 cycles
-2 -1 /0 / 2 0
./
/
/
/ / -1
/ /
I
I
/
-2
/
-3
-4
I
/
which reduces to
Y in tension a > 0;
2E(1 - D)2
Y = 0 in compression a < O.
D = YII..~.q 1f' {a eq ~ af
5 Ceq ~ CD
.
With the results obtained previously and the new conditions, it becomes
. af2
D = 2E 2 5 iF in tension
b =0 in compression.
The integration of the damage rate concerns only the positive stress taking the
micro plasticity threshold equal to zero for simplicity
bD
bN =
J
aM
a2 a2
2E;5 da = 2E;5 (aM-af)
l1 f
Other parameters may act on the phenomenon of fatigue. The influence of temper-
ature is taken into account by the material parameters which are functions of the
temperature. In high cycle fatigue the influence of frequency is negligible; in the
model it disappears when using pure plasticity at the microscale and pure elastic-
ity at the mesoscale. The influence of environment can be evaluated by chemical
means used to define appropriate internal variables; it is, however beyond, the
scope of the damage mechanics.
NR = ~ni'
~~ = 1.
NRi
This "rule" is in fact no longer a rule if a kinetic damage law is used; it is only
a consequence of the mathematical properties of the model used. Let us see under
which conditions the model developed in Section 3.4.3 has the property of linear
or nonlinear accumulation.
Consider two levels of a one-dimensional loading defined in Figure 3.30.
Using the high cycle damage evolution model written as
N - No(f':.a)
D = Die ,
N R - No(f':.a)
(j
J~M~ I I t
mn
I
0 l'1
\ I
;
n1 n2
NR
Fig. 3.30. Two-level fatigue case
142 Kinetic Laws of Damage Evolution
Find for which condition the Palmgreen-Miner's rule of linear accumulation would
be valid
iJ = II (D)I2(x) with x = a or E or p
;2
but it may give wrong results if the amplitudes differ by more than 10 or 20%.
To quantify these differences please draw the graphs as a function of ...!2:L
R2 NRI
corresponding to the two-level fatigue sequence for the conditions in Figure 3.31 :
3.5 Damage of Interfaces 143
\
,..,
.
" .... ...... ,
\"
~
\ '- ...........
0.8 ''''
\
\
...
\
\
0.6
0.4
Then, three components of the stress in the interface are known from the stresses
in layers A or B:
144 Kinetic Laws of Damage Evolution
, (1- .
I,j
ou A ou A
~dx}+~dx2
UX I UX2
Then, three components of the strain in the interface, complementary to the known
stresses, are known from the strains on layers A or B (plus a rotation condition
which does not play any role in classical continuum mechanics):
3.5 Damage of Interfaces 145
The missing terms in [0'/] and [1: / ] are of course those which cannot be defined due
to the surface character of the interface: no area to support loads in the x 1 and x 2
directions, no length to define strains in the X3 direction.
and the same for the state variable ei3 to which 0'i3 is associated
olf;
0'1 = oe = k,,: ee(1- D), ef3 have the dimensions of length,
olf;
r = -on = k n n 1/1' , r has the dimensions of stress,
ji = ~~ = - (~' : 1: 1 : 1: 1 +~ : e:e).
146 Kinetic Laws of Damage Evolution
Y =:= - ji is the strain surface energy release rate which can be written as a function
of the variables known by a classical structure calculation el and (JI
2 I~'I
Q = ~ eeqi>eq
2 St
Where 1] and St are material constants and eD' also a material constant, is the
damage threshold.
- A plastic-type potential F governed by a criterion/written as a function of an
equivalent stress (J!q taking into account the unilateral character of interfaces
which do not show any damage in compression
The ratio of the ultimate normal stress (JR to the ultimate shear stress r R models
the large difference of strength which often exists between tension and shear in
interfaces. <>are the Macauley brackets, II is the absolute value.
(J!q = 1133 in pure tension,
l1!q = 0 in pure compression,
The coupling with damage is neglected in the plastic criterion and potential:
f = (J;q - r - (Jy = 0,
ji2
F=f+ 28 ,
"
(Jy is a yield stress and S" another parameter characteristic of the interface.
The kinetic laws of evolution derive from these potentials by the normality rule
an _
oQ
SaIl - ;;-:- ,
ueafJ
of. .
n=--A=A
or '
D = _ oQ _ of i = y (e;: i;q + ~)
oji oji Se S,,'
i ~ 0 is the plastic multiplier determined by the consistency condition J = 0,
.I .
(Jeq-r= 0,
D = y [e;: i ;q(if eeq ~ eD) + ((J;q~ (Jy Y-1:S" (if {j: ~)1
Finally D = 0 if not,
D = Dc -+ debonding of a Representative Surface Element.
The rheological parameters of the interface introduced are k", ke' Se, 1] , (JR, !R' eD'
(Jy, k", 11, Dc: too many! In fact some of them may be identified as group parame-
ters and a simplified model restricts them to four only.
Last assumption to make the model as simple as possible! The yield stress, always
difficult to determine, is considered to be zero. Then
if {f=
1=0
0
or
Qualitative identification
It is possible to relate the parameters ka' k~, k" and Sa to engineering parameters
such as (1 R or 'R and others easier to measure.
- Consider the case of mode I debonding in pure tension
A mesodebonding occurs for D = Dc when the tensile stress reaches the ultimate
tensile stress (1R . These conditions are obtaned by integration of the damage
differential equation with the initial condition of an undamaged state:
0"33=0 -+ D=O
from which
J1 1 - (1- Dd 3 J1 +2
2kak~Sa (1~+2 3
from' which
Quantitative identification
For each nature of interface, four coefficients need to be determined:
aR: the ultimate stress in pure tension,
, R: the ultimate stress in pure shear,
Dc: the critical damage at debonding,
fJ.: the exponent characteristic of the nonlinearity
Even in the case of glue the identification tests must be performed on a multimate-
rial that is the interface and its two adjacent layers because an isolated interface,
glue for example, has different properties than a built-in interface.
- A pure tensile test allows the determination of a R .
- A pure shear test with unloading allows the determination of 'R and Dc by
extrapolation of the Damage evolution measured by
E G
D=l - :E=l- G,
Gbeing the damaged shear modulus and G its value in the undamaged state. This
formula using G = /2(1 + v) assumes that the Poisson's ratio is not affected by
the damage which is exact in the case of isotropic damage.
- The best way to obtain the exponent fJ. is to perform fatigue tests in tension
(aM = 0, aM) or in shear and to deduce fJ. from the slope of the Woehler curve
giving the logarithm of the stress amplitude as a function of the logarithm of
the number of cycles to debonding.
The equation of the Woehler curve is obtained by integration of the above damage
law over one cycle first neglecting the variation of D:
oD
oN
150 Kinetic Laws of Damage Evolution
(a
aD = ~ 1- (1- Dd 3 M)Il+ 2
oN 3 (1-D)2 aR '
and then as a function of the number of cycles with the conditions
N=O ~ D=O,
N=NR +- D=D c ,
Having obtained the Woehler curve from experiments we calculate the exponent
Jl from
or
- Show that this criterion is the exact solution given by the damage law when
integrated in the case of proportional loading in shear and tension:
- In the case of shear and compression, the integration is more complex due to
the sign of a 33 < 0 but the numerical values are not far from those given by
I _
a eq - aRo
3.5 Damage of Interfaces 151
For example show that for a constant value of 0" 33/0" R = -1,0"23 = 0 and Il = 4 the
integration of the damage law gives 0" d, R = 1.371 and the simple criterion
O"13/'R = 1.414.
Damage law:
_~1-(1-Dd3(
dD - 3 ( 2 1+
O"i3)(_
2
O"i3)3/2(_
1+ 2
O"i3)- 1/20"13
1+ 2 2 dO"13'
i-D) 'R 'R 'R 'R
Integration:
Criterion:
0"3310"331 + 0"13 = 1
O"~ ,~ ,
0"33> 0 --+ half a circle,
0" 33 < 0 --+ two branches of an hyperbola.
4 'to
cr
ll
r----r----r----+----+----+----+----4--0-+---- 1----~---
cr,
-8 -7 -6 -5 -4 -3 -2 -1 2
-2
-3
-4
Fig. 3.33. Debonding criterion (from J.P. Sermage)
152 Kinetic Laws of Damage Evolution
Users should be aware that all numbers of the following table are to be considered as
orders of magnitude only, because each set corresponds to a specific material which
may show slightly different properties from another one identically denominated.
When performing a structural calculation, use available characteristics chosen with
the name of the material and be prepared to perform calculation of corrections with
the numbers identified, if possible in situ, on the real material of the real structure.
Furthermore, keep in mind that each parameter is related to a specific model.
For example, the value of the yield stress of a given material may be different
for a perfectly plastic model, a linear kinematic hardening model or an isotropic
hardening model.
Dear "identify-or" , each time you perform an identification please write the results
in the following chart and send it to me! Thank you.
Material Temps Elasticity Plasticity Strain hardening Viscoplasticity DAMAGE Crack
= C closure
E v af ay au X oo -y R oo b K oo n S pD Die h
MPa MPa MPa MPa MPa MPa MPa S-1 MPa
Steel
AISI 20 190 000 0.28 200 320 700 - - 2.4 0.44 0.2 0.2
1010
Stainless 20 200 000 0.30 180 260 700 200 2 300 I - - 7 0.10 0.15 0.2
steel 600 140000 0.32 3 6 - 150 12 0.2 0 0.5 0.5
AISI316
Copper 20 100 000 0.33 100 190 300 - - 0.4 0.35 0.85 0.6
Alumi-
nium 20 72 000 0.32 250 300 500 - - 1.7 0.03 0.23 0.2
alloy
2024
Ceramic 20 400 000 0.2 300 306 500 - - - - - - 6 0.10 0.99 0.2
Alumina
at micro
scale
Concrete
in 20 30000 0.2 1.5 2.5 3 - - - - - - 2.510- 7 0 0.2 0.2
tension ;.>
V>
at micro
~e ____ L - _____ ,
o
-
3
~
o
...,
5"
&
;;;>
~
V>
\.>
-
Chapter 4
Analysis of Crack Initiation in Structures
The fundamental reason for studing damage mechanics is to understand why and
how materials break. Together with physics, metallurgy and chemistry this knowl-
edge allows us to improve the mechanical properties of materials and to design
new multimaterials. The practical reason for studing damage mechanics is to pre-
dict when materials, as they are currently made, will break upon submission to
mechanical and thermal loading. This involves the analysis of real components of
structures in real or presumed situations.
This analysis requires two main steps: modeling behavior of the material and solv-
ing the boundary value problem of stress, strain and damage fields in the structure.
The first step has been studied in detail in the previous three chapters (the con-
stitutive equations for plasticity and viscoplasticity are given in Section 2.1.4; the
constitutive equations for damage are given in Chapter 3). In the present chapter
we shall deal with the second step, by considering how these constitutive equations
may be integrated in the classical framework of continuum mechanics, the basic
variables of which are stress, strain, temperature and time. In the first section, we
review the classical equations and methods. Then, several levels of computation
are developed, whereby the coupling between damage and strain is: neglected for
a rough approximation; taken into account only for the most damaged points in
small-scale damage zones; or, taken into account everywhere for a fully coupled
analysis of structure in large-scale damaged zones.
Point M
Structure S - j - - -
Boundary 6S
Due to the complexity of the partial differential equations and the nonlinearities
introduced by the constitutive equations, there is no analytical method for solving
this type of problem, in the general case. It can only be solved by means of
approximations in some cases and by numerical timelike incremental analysis,
usually with the aid of computers.
Figure 4.2 gives an example of a small circular hole in a large elastic plate. On
the boundary of the hole where x 2 = 0 the stress is one-dimensional.
Many particular. geometries have been analyzed either by approximate analytical
methods, photoelasticity or the numerical finite element method. Some classical
results taken from handbooks listed in the literature references are given below.
All concern isotropic elasticity.
h I
20
--t1.
I.
Fig. 4.3.
{
fla > lO
"Infinite" plate ---+
hlb > lO
Stress concentration factor at A and in the direction of O"~
o a
K T -- 0"22
0 1+2-
0"2 b
I
o
- -~~- 0'1
Fig. 4.4.
aA or B or C
22
KT aO
2
~
KT
0'0
2
4
L1
3
C
2
B
C A B
- 0
3 4 5 I/r
-1
-2
-3
-4
Fig. 4.5.
4.1 Stress-Strain Analysis 159
{ Kf
K~ = 3.5
1
Example: - = 2 ....... =0
r
Kfj = 3.15
\
J
\
i\ .........
r-- K2
I
I
0.2 o
d/b
Fig. 4.6.
Fig. 4.7.
Bending -- KT
Example: hlr = 4 -- KT = 5
Fig. 4.8.
4.1 Stress-Strain Analysis 161
~\
Tension 0.5 1.70 1.60 1.53 1.47 1.39 1.21
A
K - a 22 T
- 0.05 0.10 0.20 0.27 0.50 1.0
T - 0
a 22A
~\
0.5 1.61 1.49 1.39 1.34 1.22 1.07
Fig. 4.9.
162 Analysis of Crack Initiation in Structures
-Square shoulder with fillet in circular shaft subjected to torsion (Figure 4.10)
C--;c
Fig. 4.10.
K -
T -
A
a R8
a 0R8A
r 0.005 0.01 0.02 0.03 0.04 1.06 0.08 0.10 0.12
~\
Bending 2.00 3.0 2.25 2.00 1.82 1.65 1.51 1.44 1.39
1.20 3.00 2.5 2.00 1.75 1.62 1.50 1.40 1.34 1.30
1.09 2.20 1.88 1.53 1.40 1.30 1.20 1.16 1.15 1.15
A
K - a R8 h
a\r
T - 0 0.5 3 5 9
a R8
0:
Fig_ 4_11.
When plasticity occurs, even locally, it is much more complicated. Neuber's method
allows us to calculate the elastoplastic stress concentration from a purely elastic
calculation. It applies for all notch problems at the root of which there is a stress
concentration.
The basic heuristic hypothesis, checked for thousands of applications, postulates
that for a one-dimentional elastic problem, the product of the stress by the strain
a . at the most loaded point is equal to the same product for the same problem
analyzed in pure elasticity a E . E
I", ~ aea E ~~ I
Assuming the elastic problem to be solved, a E and E are known. The second
equation with which to calculate the two variable a and is the elastoplastic
constitutive equation. For the simple case of monotonic loading, let us take some
stress function representing the stress strain curve:
Then
a2
agc(a) = : = const.
is the equation with which to find the stress at the root of an elasto plastic notch.
A simple geometrical construction follows from the above two equations.
In the stress strain coordinate system;
a2
.-K
a
E
is an hyperbola containing the point (a = a E' aE )
E .
164 Analysis of Crack Initiation in Structures
is the classical elastic strain hardening stress strain curve. The solution is obviously
given by the intersection of these two curves (Figure 4.12).
This method may be generalized straighforwardly to cyclic loading by using the
cyclic elasto plastic stress strain curve 9 c to relate the amplitude of stress and strain:
(~O"E?
~O" . ~c = ~O" E . ~c E = E
~c = 9c (~0").
The generalization to a three-dimensional state of stress is also easy if we restrict
the application to proportional loadings.
The basic heuristic hypothesis is
O"ijc'fj = 2w e , the elastic strain energy which has been calculated in Section 2.1.2
as:
2(I+v) 2 3(1-2v) 2
2we = 3E O"eq+ E O"H
with D = 0,
\
\
\
\
\
dE ------- Elastic solution
.
since p .IS a d
Cij ' '8i jCij
eVlator. P-O-
For monotonic proportional loading the plastic strain may be written as (see Sec-
tion 2.1.4)
3 aIJ
1.)
-2 g ( a eq ) - ,
a eq
then
a;; being known from the elastic calculation, a eq may be determined if a~ (which
cannot be calculated without the complete analysis) is considered close to (a~)2
to make the last term negligible.
This method allows us to numerically solve the set of equations from problems
in mechanics as described in Section 4.1. It consists in replacing the problem of
partial differential equations by a problem of linear algebraic equations in which
the unknowns are the displacements of the nodes defining the finite elements into
which the structure is devided. In each element the continuum mechanical equa-
tions are solved for a class of simple displacement fields: linear, quadratic or of
higher order with space coordinates. By consequence of the piecewise approxi-
mation of the displacement, and of the strain and stress fields, we are obliged
to considerer a mesh size small enough to discretize the gradients of stress or
strain.
The algebric system to be solved is linear for elasticty, the only difficulty being
the size of the banded matrix for 1000, 10000, 100000 or even more degrees of
freedom. When plasticity and damage occur, the problem becomes nonlinear and
the only way to solve it is by time like or load incremental linearization.
166 Analysis of Crack Initiation in Structures
It is not the purpose of this course to describe the finite element method in detail;
many good books already do. But as we are concerned with damage, two dificulties
arise: the choice of the meshes and the choice of the time or load increments.
Convergence regarding these two parameters is difficult to ensure and only heuristic
procedures may be advised:
- The algorithms used to integrate the highly nonlinear constitutive equations must
be "strong" relative to stability.
- Implicit integration schemes associated with Newton's iterative method are
preferred
- The local integration of the constitutive equations is performed by using an
"elastic predictor" to calculate the first increment as elastic and a "plastic
corrector" to ensure, by iterations, that the plasticity criterion and kinetic laws
are satisfied up to the required accuracy.
- In order to compute the tangent stiffness matrix of the structure at each it-
eration, a tangent modulus J consistent with the discretization, calculated as
baij = Jijklbckl is used instead of aij = LijklEkl. This ensures a much better
convergence.
- As to the space discretization is concern, the convergence regarding the mesh
size is of the same order of difficulty as for classical elastoplastic calculations
up to crack initiation, that is just before the localization of damage. Afterwards
the problem is no longer elliptic; its solution involves strain rate discontinuities
which cannot be modeled by classical finite elements.
- When a periodic loading is considered, as for fatigue, it is not possible to
perform the tremendous number of increments needed for all of the thousands,
or millions of cycles. A method for saving computer time involves the use of
a "jump-in-cycles" procedure which avoids calCulating the stress strain cycles
when the process is a quasi-steady state process.
Assume a piecewise periodic loading and consider a number of cycles N i for
which the damage is known to be Di. The structural calculation of this cycle (i)
is performed several times until the stabilization is reached at each point. The
coupled constitutive equations are such that for constant damage, a periodic
input implies an asymptotic periodic output. Only two or three calculations are
sufficient.
This gives (~~) i at each Gauss point of the elements of the structure, this
rate will be assumed to be constant during a certain number of cycles !:1N i,
which may be very large.
The criterion for choosing !:1N i, is a certain amount of damage !:1D at the most
damaged point M*. !:1D must not be so large as to violate the coupling, but large
enough to obtain a reasonable computer time. For applications !:1D = D lc /50
4.1 Stress-Strain Analysis 167
tlN = tlp
8p ,
8N
where tlp is the best compromise between accuracy and rapidity. CR being the
order of magnitude of the strain to rupture in pure tension, tlp = ~~ gives
good results. Nevertheless, be careful when using those heuristic procedures!
Consider a plate with a central hole of small radius in comparison to the dimensions
of the plate, which is loaded in tension by a "far" stress 0'22 = O'~ (Figure 4.13).
Assuming the material to be linear elastic, find the stress concentration factor at
the boundary of the hole (XI = r,x2 = 0) from Section 4.1.1;
K T -- 0'22
0
-
-
3.
0'2
Considering an increase of the external load O'~, plasticity may occur at the point
of stress concentration. Assuming the power strain function,
where O'y is the yield stress Ko and M are other parameters, determine the stress
concentration factor KTfor 0'22 2:: O'y by Neuber's method;
0'22 c 22 = O'E c E ,
168 Analysis of Crack Initiation in Structures
0"22)M _ (30"g)2
0"22 ( K - E '
E
= 3 (SL)
30"g
~:;::
0"0
Draw the graph of KT as a function of -1. for M = 3 and for the limiting case
O"y
of perfect plasticity M = 00 (Figure 4.14).
3 + - - - -....
,..* -,..
v -
RI /2 , Rv =
V eq v ~(I+V)+3(1_2V)(aH)2
3 a eq
The critical time t R at which a crack is initiated is reached when the damage
itself reaches its critical value given by the rupture criterion (see Section 3.1.3)
(12
Dc = Dlc (1-:2 : :;
1.
Let us take here for simplicity
Dc = 1.
4.2 Uncoupled Analysis of Crack Initiation 171
I
t R is given by the expression
J
tR
J
tR
p(t) dt
to
The strain history is imposed on the mesovolume element considered; this allows
us to take p( t) as a given function, which is particularly simple to calculate if the
elastic strain is neglected:
2 ) 1/ 2
p(t) = ( "3 Ei j (tAj (t )
- the same strain history is imposed and the elastic strain is again neglected: the
function p( t} is the same.
The critical time for crack initiation in the coupled case t'n is deduced from the
same kinetic damage law:
. a;qR v .
D = 2ES(1 _ D)2P(t}H(p-PDl
or
iJ .( }H
= a;Rv
2ES P t (p-PDl'
where (I - D)2 disappears due to the coupling in the plasticity criterion. The
integration is obvious:
t
2ES J(}d
D = a;Rv P t t.
to
Taking again D = 1 as the critical condition:
D = 1~ t = t'n,
J
t fR
2ES
p(t) dt ~R
as v
to
Comparison with the uncoupled case shows that
J J
t~ tR
2 ES ES
tR = to+-3~R
as vp as vP
t'n = to+2~R.
The uncoupled calculation is always conservative. It is a lower bound, but far from
the coupled solution. In other terms, component design may benefit from coupled
calculations, whkh can prove the enhanced safety of components or indicate how
light-weight economical components could be built.
The scheme of the uncoupled method is sketched in Figure 4.15.
Elasto-(plastic) Kinetic
constitutive equations damage law
0:0
2 0'0
2
Ad,ot,
vL
-r
I
. x1 --
rfy
o'M
Po
P D
.".
,-
,
./ _:~---
I
....,I
I I
0 Ny NO NR
Lla~ = 8N.
2a y
N =-
y 38
Step 2: Ny ::; N < No.
The material is elastoplastic at the free boundary of the hole (Xl = r, x 2 = 0).
Then the stress concentration factor is now a function of the stress at infinity and by
consequence, a function of the number of cycles. For a cyclic loading (Section 4.1.2
and Exercise 4.1.4),
KT = 3(3~;1) :::;: ,
M-I
2a ) M+I 2
Lla = KTLla~ = 38 ( 3: N M+I .
Using the cyclic strain-hardening power function already mentioned in Section 3.4.1
(here without any damage),
J
No
2 LlCp(N) dN = PD = cPD
Ny
from which
If you work with numbers for a specific material this formula will not look so bad!
4.2 Uncoupled Analysis of Crack Initiation 175
Step 3: No ~ N ~ N R"
A low cycle fatigue damage occurs with a varying amplitude of stress. Take the
damage per cycle derived in Section 3.4.1:
8D K2 M2
= P LlE M
8N 4ES P
The uncoupled hypothesis assumes the plastic strain amplitude to be related to the
stress amplitude without any consideration of damage:
38(~)~ M
NMI.
2M
Kp
Quite often, the damage is so localized that the volume of the damaged material
is small in comparison to the macroscale of the structural component and even
to the RYE. This allows us to perform an uncoupled analysis at the macroscale
and to consider the coupling between strain and damage only on the RYE of the
critical point as shown schematically in Figure 4.17. This is the case of small-scale
damage.
The uncoupled analysis differs from the coupled calculation at the level of the
constitutive equations of the critical point. The latter gives a more accurate result,
but remember that the uncoupled analysis is a lower bound on the lifetime.
The method of locally coupled analysis may often be used with good accuracy for
brittle and fatigue types of damage.
Damage localization results from stress concentration, of course, but also occurs,
because some weakness always exists at the microscale. Let us generalize, for any
kind of damage, what has been said for quasi-brittle materials and for high cycle
fatigue in Section 3.2.2 and 3.4.3. The mechanical model was, a two-scale volume
element, elastic or elastoplastic at the mesoscale and elastoplastic and damageable
at the microscale (Figure 3.8).
The only material characteristics which differ in the matrix and in the inclusion is
the yield stress of the inclusion a~ , which is assumed to be equal to the fatique
limit a f of the material, and its fatigue limit, assumed to be reduced in the same
proportion:
Elasto-(plastic) Coupled
constitutive equations elasto-plastic and damage
constitutive equations
This takes into account the microinternal stresses and the weak defects always
existing everywhere in all materials. As the fatigue limit is smaller than the yield
stress, it allows for plasticity which induces damage at the microscale, whereas the
matrix remains elastic or elastoplastic and undamaged for the same loading.
The second assumption which simplifies calculations is the Lin-Taylor strain com-
patibility hypothesis, which states that the state of strain at the microscale is equal
to the state of strain at the mesoscale as derived from the classical structural cal-
culation,
g'" = g
Then, there is no boundary value problem to be considered. Only the set of coupled
constitutive equations must be solved for the given history of the strain.
As in Section 3.2.2, the crack initiation condition at the microscale D = Dc is
also the brittle crack initiation at the mesoscale G = G c.
As the input of the locally coupled method is the strain at some point, it can be
set up as a postprocessor to any finite element code of structural calculations. The
most sophisticated constitutive equations to be used are:
- the set of equations developed in Section 2.1.4 for linear elasticity and for
isotropic and kinematic plastic hardening,
- associated with the kinetic damage law developed in Section 3.1;
- in the case of unilateral damage microcrack closure conditions as explained in
Section 2.3.4.
(J eq
1 _ D - (J~ = 0 and f = 0; if not, ifj = 0,
IL2
D=~R . if not, D = 0;
2ES vP
178 Analysis of Crack Initiation in Structures
2
Crack initiation if De = Dl e f -
(J~ Rv
These equations may be used for piecewise perfect plasticity by considering several
values of the plastic threshold (J~ as the loading or the timelike parameter vary.
This allows us to take some cases of high values of strainhardening and the cyclic
stress strain curve for multilevel fatigue processes into consideration.
Then, the material parameters must be considered as follows:
f = aeq - (J~ = 0,
4.3 Locally Coupled Analysis 179
with >. = Ev E
(I - 2v)(1 + v) J.l = 2(1 + v)'
Efjn is the plastic strain at the beginning of the increment, t:::.p is accumulated
plastic strain increment, and the other quantities are the values at the end of the
increment.
The iteration operates on
J+ ::'CfJ
'J
= 0,
&h &h
hij 'J 'J C = 0 ,
C kl+-&
lI
+-&- p
akl P
where CfJ and C p may be explicitly calculated. The resulting expressions are
C =
J -Nh
'J 'J 'th N .
p 3J.l ' WI 'J
DAMAGE 90 is a friendly computer code with fewer than 600 FORTRAN instruc-
tions; it calculates the evolution of damage up to crack initiation with the method
described in Section 4.3.2. It may be used in an interactive way either as a postpro-
cessor of the strain output of a finite element code or as an autonomous program.
It is written in FORTRAN 77 as available on a Digital V AX computer. The CPU
time for one complete execution is a few seconds. It distinguishes between two
loading cases:
- general loading history where the history of loading is defined by the values
of the strain components at given timelike parameter values. DAMAGE 90
interpolates linearly between these values;
- piecewise periodic loading for which the loading is a certain number of blocks
of cycles defined by the two consecutive maximum and minimum set of strain
components and the number of cycles in each block. The strain interpolations are
also linear. For a large number of cycles, the jump-in-cycles procedure (which
alows "jumps" of cycles to save computer time as explaned in Section 4.1.3)
may be used. It is capable of accounting for initial values of damage Do and
plastic strain Po'
t 80 Analysis of Crack Initiation in Structures
The questions asked by DAMAGE 90 to introduce the input data and the charts of
the output data as they appear on the screen of the terminal are given below. (The
subscript f..l is omitted a!' == as).
The complete listing is also given hereafter for those of you who wish to play with
DAMAGE 90.
UI DAMAGE91) u,
wI II give you the damage gl'owth up to crack I nit I at Ion
***111111111111111111111111111.11111111111111111111111
u ELASTICITY . Give
YOUNG's modulus
200 . E3
POISSOW s I'at. 10
0.3
I I PERFECT PLASTICITY: plastic threshold SIGs given with loading. Give
Fatigue limit. SIGf
21)1).
Yield st.ress SIGy
31)1).
Ultlmat.e str'ess SIGu
51)0.
DAMRGE EVOLUTION dD (Y/ S ) dp . Give : ~;
13.El6
u DAI1AGE THRESHOLD dD =13 If p(pD. Do you know I'D? 'Y' or 'Ii'
Y
Give t.he value of pD :
10.E-2
CRACK INITIATION : D'Dc . Do you know Dc ? 'Y' or ' N'
Y
Give the va I ue 0 f Dc (remember : 13 ( Dc ( 1) :
0.99
U : INITIAL CONDITIONS Give :
The value of po
0.
The value of Do
D.
H LOADIIiG
Is t.he stress state uniaxial? ' Y' or ' N'
Y
Is t.he st.raln history cyclic? 'Y' or 'N'
Ii
4.3 Locally Coupled Analysis 181
Itl I I II.II.IIIII.II.I.III 1
Give male,' i cd conslanls and lhe sl,'a ins h i slo,'~!
II ELASTICITY. Give
YOUNG ' e modu 1us
72.E ... :::
POISSOW s "'al io
0.32
I I PERFECT PLASTICITY: plastic lhr eshold SIGs given wilh loading. Gille
Fatigue limil SIGf
303 .
Yield stress SIGy
3136.
Ultimate slress SIGu
500.
:1:1 DAI'1AGE EVOLUTIOli dD (Y/ S) dp . Give : S
6.
u Dl'li'1AGE THRESHOLD dDO if p<pD. Do you know pO? ''I'' or' "1'1 '
Ii
Give EpD in lens ion: pDEpDI(SIGu -S IGf) / [SIGs- IS IGfl.2/ SIGyIJ
10. [-2
II CRACK INITIATION : DDc . Do you know Dc ? 'Y' or ' N'
Ii
Give D1e in lension Dc D1e I [(SIGu/ SIGs).12J / Rv
13 . 99
INITIAL CONDITIONS Give:
The value of po
0.
The ....J'due of Do
0.
U : LOADIIiG
Is lhe sl,'ee.s st.ale uniCl xia l ? ""1" or' ' 1'1 '
N
Is lhe sl,'a in h i 5lo,'~! eye 1 ic ? ",// or- ' Ii'
v
182 Analysis of Crack Initiation in Structures
Listing of DAMAGE 90
c**********************"'*"*""',.
C'" DAMAGE 90
C Elastic-perfectly plastic law coupled to a ductile damage model.
C'l' Fully implicit integration scheme.
C'2' Jump in cycles procedure.
C'" Written by Issam Doghri
C U. C. SANTA BARBARA, DEPT . Of MECHANICAL ENGINEERING
C'" Version : May 1990
C"*'It*'ltlt""'***"*"'*""*"'*"*'**"""'*""** **'*""'*'*'"
IMPLICIT REAL'S (A-H,o-Z)
CHARACTER'l UNIAXI,ANsPD,ANsDC,CYCLIC,JUMP,CONVER,sTAB,COUPL
CHARACTER'20 flLE1,flLE2
CHARACTER' 6 COMMENT
PARAMETER(NTENs=6,NSTATV=ll)
REAL'8 sTRS(NTENS) ,sTATEV(NsTATV) ,STATEVI (NsTATV),
sTRAN(NTENs) ,DsTRAN(NTENS) ,sIGHI (6) ,SIGHf(6),
INCUSE(50),TIM(200),HIsT(6,200),EPsB(6,50),EPss(6,50),
PERDIV(6),sYIEL(50)
INTEGER ISLOPE(6),IPER(6),INDEX(6),NBCYCL(50)
COMMON/ETIQ1jEO,XNU,sIGS,SMPD,SO
COMMON/ETIQ2/0NE,TWO,THREE
DATA INDEX 111,22,33,12,13,231
C'" READ DATA
MS=6
MB=5
NAXI=6
NCYCLE=O
sMPD=l.D+IO
DCRIT=O . 99DO
WRITE (MS , ')'
WRITE (MS , ')'
WRITE(MS,"')' ft. It". **10 It **, ***,. * ,** **11 *It.,. *"" **'" ,," ft** "It, #tit.
WRITE(MS,')
WRITE(MS, *)' It* ELASTICITY . Give
WRITE(MS,') , YOUNG"s modulus
READ(MB,' )EO
WRITE(MS,') , POISSON"s ratio
READ(MB,')XNU
WRITE(MS,')'" PERfECT PLASTICITY plastic threshold sIGs
. given with loading . Give :'
WRITE(MS,'), fatigue limit SIGf
READ(MB,' )SIGf
WRITE(MS,'), Yield stress sIGy
READ(MB,' )SIGY
WRITE(MS,'), Ultimate stress SIGu .,
READ(MB,')sIGU
WRITE(MS,')'" DAMAGE EVOLUTION dD = (Y/S) dp . Give S'
READ(MB,')SO
FORMAT (A)
WRITE(MS,')'" DAMAGE THRESHOLD dD=0 if p<pD . Do you know pD?
, ,y" or ' 'N'"
READ(MD,l)ANsPD
If( (ANSPD.EQ. 'y') .OR. (ANSPD.EQ. 'Y') )THEN
WRITE(MS,')' Give the value of pD :'
READ(MD,')sHPD
If(SMPD.EQ.0.)SMPD=1.D-6
ELSE If( (ANsPD.EQ. 'n') .OR. (ANSPD.EQ. 'N') ITHEN
186 Analysis of Crack Initiation in Structures
CALL OtrrPtrr (NUl. NU2 TIME. STnAN. STRS STATEV NTENS NSTA'1'V
STRSB.STAR)
~IILE(TIME+DTIME.LE . TFIN)
DO IST=l. NAXI
DSTnAN( IST)=DTIME' (IIIST( 1ST. IGIV+l )-IIIST( 1ST. IGIV) )
/(TIM(IGIV+l)-TIM(IGIV))
END DO
CONVER='y'
DO IST=l.NSTATV
STATEV(IST)=STATEVI(IST)
END DO
SIGS=SYIEL(IGIV)
IF(SIGS . LT.SIGF.OR.SIGS . GT.SIGU)STOp .. PLEASE REVIE:W TilE
VALUES OF SIGf SIGu and SIGs .'
CALL INTEGR(STnAN.DSTnAN.NTENS.NSTATV.UNIAXI.COUPL .
& CONVER.STRS.STATEV.STRSB.STAR.RNU)
C'" IF NO CONVERGENCE. DIVIDE TIME INC. BY 2
IF(CONVER.EQ.y . AND . STATEV(ll).GT.O.)THEN
ENERG2=ENERG3+(SIGS-SIGf'(SIGF/SIGY))'STATEV(11)
IF(ANSDC . EQ.n)DCRIT=DIC((SIGU/SIGS)2)!RNU
IF(ILOOKO. EQ. 0 .AND. (ANSPD. EQ. 'y' .OR.ANSPD. EQ . 'Y' )
.AND.STATEV(2).GT.l . OSSMPD.OR.
ILOOKO . EQ.O .AND. (ANSPD. EQ . 'n' . OR .ANSPO. EQ . N' ) .111'10.
ENERG2 .GT . 1. OSENERG1 . OH .
STATEV(3).GT . l . OSDCHIT)TIIEN
CONVER='n'
END IF
END IF
IF(CONVER.EQ.n)THEN
IPASS=IPASS+l
DTIME=DTIME/2.
C'" IF CONVERGENCE
ELSE
IF( (ANSDC.EQ . 'n' . OR . ANSDC . EQ . 'N')
.AND . DCRIT.GT.O.99DO)DCRIT=O.99DO
IPASS=O
DO IST=l.NSTATV
STATEVI(IST)=STATEV(IST)
END DO
DO IST=1.6
STRJIN ( 1ST) =STHIIN ( 1ST) +DSTRJIN ( 1ST)
END DO
TIME=TIME+DTIME
CALL OtrrPtrr (NUl. NU2 TIME. STRAN. STRS. STATEV NTENS NS1ATV
STRSB.STAR)
IF(STATEV(3).GE . DCRIT)THEN
WH1TE(MS.)
WHITE(MS.) CRACK INITIATION.'
GO TO 308
END IF
IF(STATEV(3).GE.l.l)STOp DAMAGE CANNOT ~XCEED 1 . '
C'" FIND pD
IF(ILOOKO.EQ.O .AND . STATEV(ll) . GT.O. )THEN
ENERG3=ENERG2
IF( (ANSPD.EQ. 'n' .OR.ANSPD.EQ.W)
.AND. ENERG3 .GE.ENERGl )TIlEN
SMPD=STATEV(2)
ILOOKO=l
ELSE IF( (ANSPD.EQ. 'y') .OR. (IINSPD.EQ. 'Y') .IIND.
STATEV(2).GE.SMPD)THEN
ILOOKO=l
END IF
END IF
IF(DAllS(TIME-TIM(IGIV+l)) .LE . DTMIN)TIIEN
IGIV=IGIV+l
END If
END IF
DTIME=DTIME1.1
IF(TIME+DTIME.GT.TIM(IGIV+l)+DTMIN.AND . IGIV+l.LE.NGIV)TIlEN
DTIME=TIM(IGIV+l)-TIME
END IF
4.3 Locally Coupled Analysis 189
IF(IPASS.EQ.13.0R . DTIME.LT.DTMIN)TIIEN
WRITE(MS,')"" NO CONVERGENCE'
WRITE(MS,')'DTlME=' ,DTlME,'IPASS=' ,IPASS
GO TO 308
END IF
END DO
END IF
C'" IF TIlE LOADING IS CYCLIC
IF( (CYCLIC.EQ. 'y') .OR. (CYCLIC . EQ. 'Y') )TIIEN
TPER=1.
DTMIN=TPER/40000 .
TlME=O.
ICYCLE=1
NRCYCL=ICYCLE
DIITlME=O.
DO 1=1,6
SIGIII(I)=O.
END DO
SMPI=SMPO
0111=00
IF( (JUMP.EQ. 'y') .OR . (JUMP.EQ. 'Y') )THEN
COUPL='n'
OPEN (UNIT=80,name='fatigue .out' ,status='unknown')
WRITE (80, ' ) 'REAL TIME TIME D p'
WRITE (80, #I ) , #I " " #I" #I #I" 11"" #I "'. #I"" #I *" .. *""" .. I< #I ,
#I
C IF CONVERGENCE
ELSE
11'( (ANSDC.EQ. 'n' .OR.ANSDC.EQ. 'N')
.AND.DCRIT . GT.O. 99DO )OCRIT=O. 99DO
IPJlSS=O
DO IST=1. NSTATV
STATEVI(IST)=STATEV(IST)
END DO
DO IST=1.6
STRAN(IST)=STRAN(IST)+DSTRAN(IST)
END DO
TIME=TIME+DTHIE
RTIME=TIME+DHTIME
CALL OUTPUT(NU1.NU2.RTIME.STRAN.STRS.STIITEV.NTENS.NSTATV.
STRSIl.STIIR)
IF (STIITEV( 3) .GE .DCRIT)TIlEN
WRITE(MS. *).
WRITE(MS.) CRACK INITIATION .'
GO TO 300
END IF
IF (STATEV( 3) .GE . 1.1 )STOI' .. DAMAGE CANNOT EXCEED 1.'
DO IST=1.6
11'( DABS (TIME-IPER( 1ST) *PERDIV( 1ST) ) .LE. DTMIN)TIIEN
IPER(IST)=IPER(IST)+1
DTIME=1./INCUSE(IBL)
END IF
11'( ISLOPE( 1ST) .GT. 0 .AND. DABS (STRAN( IST)-EPSB( IS1'. IIlL) ) . LE.
EPSMIN.OR.ISLOPE(IST).LT.O.AND.DABS(STRAN(IST)-EPSS(IST.IBL))
. .LE.EPSMIN)THEN
ISLOPE(IST)=-ISLOPE(IST)
END IF
END DO
C'" END OF A CYCLE
IF(DABS(TIME-TPERICYCLE) . LE . DTMIN)THEN
IDELTN=O
DO 1=1,6
SIGHF ( I )=S'rRS ( I )
END DO
SMPF=STATEV(2)
DHF=STATEV(3)
DELTP=SMPF-SMPI
DELTD=DHF-DIII
C'" FIND pD
IF(ILOOKO.EQ.O)THEN
ENERG3=ENERGO+(SIGS-SIGF*(SIGF/SIGY))'
DELTP'(NRCYCL-NCYSUM)
11'( (IINSPD.EQ. 'n' .OR.ANSPD.EQ. 'N')
.AND.ENERG3 . GE . ENERG1)THEN
SMPD=SMPF
NCYCLO=NRCYCL-NCYSUM
ILOOKO=l
ELSE 11'( (ANSPD.EQ. 'y' .OR.ANSPD.EQ. "';') .AND.
ST~TEV(2).GE . SMPD)THEN
NCYCLO=NRCYCL-NCYSUM
ILOOKO=l
END IF
END IF
C'" JUMP IN CYCLES PROCEDURE
11'( (JUMP . EQ . 'y') .OR. (JUMP . EQ. 'Y') )TIIEN
WRITE (00, .) END OF TIlE CYCLE'
WRITE ( 00.400 )RTIME. TIME. Dill' .S~IPF
DDMIIXO=DCRIT/50.
11'( (IINSDC.EQ. 'n') .OR. (IINSDC . EQ.N'))
DDMllXO=DlC/50 .
11'( (UNIAXI .EQ. 'y' ) .OR . (UNIAXLEQ . 'Y') )TIIEN
YYY=SIGSSIGS/ 2./EO
ELSE
TRDEPS=O .
DO IST=l.3
TRDEPS=TRDEPS+EPSIl(IST.IBL)-EI'SS(IST.IIlL)
END DO
4.3 Locally Coupled Analysis 191
YYl=(l.+XNU)*SIGS'SIGS/3./E0
YY2=EO*TRDEPS*TRDEPS/(l.-2.*XNU)/6.
YYY=YYl+YY2
END If
DPMAXO=SO'DDMAXOjYYY
DPMAX=DPMAXO
PSLOPE=NBCYCL(IDL)*DELTP
If(DPMAX.GT.PSLOPE)DPMAX=PSLOPE
WRITE(60,')'ILOOKO=' ,ILOOKO,'NCYCLO=' ,NCYCLO,'pD=' ,SMPD
If(COUPL . EQ . 'n' .OR . STATEV( 2) . LT. SMPD) Til EN
STAB='y'
DO KK=l,6
If(DABS(SIGllf(KK)-SIGUI (KK) ) .GT. Eill ) STNI= , n'
END DO
WRITE(80,*)'ENERGl=' ,ENEHGl,'ENERG3=' , ENERG3
If(STAB.EQ. 'y') TIIEN
WRITE(60,') 'STABILIZED CYCLE'
If (STATEV( 2) .LT .SMPD)T1IEN
C** * JUMP Of CYCLES BEfORE DAMAGE GRowrll
IDELTN=IDINT(DPMAX/DELTP)
If(NRCYCL+IDELTN.GT.NCYSUM+NBCYCL(IBL))
IDELTN=NCYSUM+NBCYCL(IBL)-NRCYCL
ENERG2=ENERGO+(SIGS-SIGY) *DELTP' (NRCYCL+IDELTN-NCYSUM)
If(ENERG2 . GT.ENERGl'l . 05)
IDELTN=(ENERGl-ENEllG3)/(SIGS-SIGY)/DELTP
WRITE(60,*)'*** JUMP Of ',IDELTN,' CYCLES'
WRITE(60,')'DPMAX=' ,IDELTN'DELTP, 'YYY=' ,YYY
ELSE
WRITE(60,') 'COUPLED CO~WUTATION fOil NEXT CYCLE
COUPL='y'
END If
END If
I::LSI::
C'" JUMP Of CYCLES ,\FTI::H DNlAGI:: GHOWfIi
DOMJ\X=DDNA.,(O
DSLOPE=(NBCYCL(IBL)-NCYCLO)'DELTD
If(DDMAX . GT.DSLOPE)DDMAX=DSLOPE
NCYCLO=O
IDELTN=IDINT(DDMAX/DELTD)
If(NRCYCL+IDELTN.GT . NCYSUM+NBCYCL(IDL))
IDELTN=NCYSUM+NBCYCL(IBL)-NRCYCL
If(Dllf+DELTD'IDELTN.GT.DCRIT'l.05)
IDELTN=IDINT((DCRIT-DHf)/DELTD)
If(IDELTN.GT.IDINT(DPMAX/DELTP))IDELTN=IDINT(DPMAX/DELTP)
WRITE (80, * ) , ** * JUMP Of ',IDELTN,' CYCLES'
WRITE(80,')'DDMAX=',IDELTN*DELTD, 'DPMAX=' ,IDELTN'DELTP
COUPL='n'
END If
If(DELTI'.GT . DPMAX)WRITE(80,')' .. PHOIJL"M
DPMAXO IS TOO SMALL . '
If(DELTD.GT . DDMAX)WHITE(80 , ')' ... PHODLEM
DDMAXO IS TOO SMALL . '
If(DELTD . GT . DDNN( .OH . DELTP .GT . DPMAX) IDEL1'N=0
If(IDELTN . LT.O.)STOP"'* PROBLEM: NEGATIVE JUMP Of CYCLES . '
END If
DO J=l,6
SIGIII(J)=SIGllf(J)
END DO
SMPI=SMPf+DELTP'IDELTN
DUI=Dllf+DELTD' I DELTN
STATEVI(2)=SMPI
STATEVI(3)=DIII
D1ITlME=DIITIME+TPJ::H'IDELTN
NRCYCL=NRCYCL+IDELTN
R'fIME=TlME+DIITIME
If( (JUM.I'.EQ . 'y') . OIL (JUll1' . EQ . ', ' ) )THEN
W1\ITE (80, ' ) ,-------- -----------------,
WRITE(60,') 'BEGINNING Of TilE CYCLE'
WRITE(80,')'CYCLE (REAL) =' ,NHCYCL+l
WRITE(80,')'CYCLE (MACHINE) =' , ICYCLE+l
WRITE(80,400)RTINE,TIME,DHI , SMl'I
END If
192 Analysis of Crack Initiation in Structures
SUBROUTINE INTEGR(STRAN.DSTRAN.NTENS.NSTATV.UNIAXI.COUPL.
CONVER.STRS.STATEV.STRSB.STAR.RNU)
CHITERE= . FALSE .
ELSE
C'" PLIISTIC CORRECTIONS
CALL VTRANVl(N,KESTHS,XNKES,NISS)
C. . . Correction over p
DENOM= 3 . 'XMU
CSMP=(F-XNKES)jDENOM
C.. .. Correction over effective stress
DENOM=DENOM'DSMP jESTHSIJ+ (1. )
DO I=l,NISS
CESTRS(I)=(2.j3.)'(XNKES-DENOM'F)'N(I)-KESTRS(I)
CESTRS(I)=CESTHS(I)/DENOM
END DO
C'" END OF THE ITERATIVE TEST ON THE YIELD CONDITION
END IF
C'" END OF THE LOOP ON THE PLIISTIC CORRECTIONS
END DO
C'" END OF TEST ON THE YIELD CONDITION
END IF
500 CONTINUE
IF(CONVER.EQ . 'y' )TIIEN
IF( (UNIIIXI .EQ. 'y') .OR. (UNIIIXI .EQ. 'Y' ) )THEN
DSTRAN(2)=STRN(2)-STRAN(2)
DSTRAN(3)=DSTRAN(2)
END IF
C.... Plastic strain increment
DO I=l,NISS
DPSTRN(I)=N(I)'DSMP
END DO
C. ... Damage increment
DO=O .
CALL DAMAGE(THEPS,ESTRSIJ,Y,RNU)
IF(DSMP.GT . O. . AND.SMP+DSMP.GE.SMPD.AND .
COUPL . EQ . 'y' )TIIEN
DO=Y'DSMPjSO
END IF
STAR=(ONE-D-DD)'DSQRT(2.'EO'Y)
C ... plastic multiplier increment
DSMR=(ONE-D-DD)'DSMP
C.... Compute stresse s
DO I=l,NISS
STRS(I)=ESTRS(I)'(ONE-D-DD)
END DO
STHSIJ=ESTHSIJ'(ONI::-D-DD)
C.... STORE STATEV AT THE END OF THE INCHEMENT
STATEV(l)=SMR+DSMR
STATEV(2)=SMP+DSMP
STATEV( 3 )=D+DD
STATEV(4)=DD
DO I=l,NISS
STATEV( 1+4 )=PSTIlN ( I ) +DPSTRN ( I)
END DO
4.3 Locally Coupled Analysis 195
STATEV(ll)=DSMP
END IF
RETURN
END
SUBROUTINE INVAR(V,~I,VB,NISS)
C*~*********************************
C . . . 1st and 2nd stress invariant s
IMPLICIT REAL'8(A-H,O-Z)
REAL'O V(6)
COMMON/ETIQ2/0NE ,TWO ,THREE
C
VH=(V( 1 )+V( 2 )+V( 3) )/TIIREE
CONST=O.
00 I=l,NISS
IF(I.LE . 3)CONST=CONST+(V(I)-VH)'(V(I)-VH)
IF(I.GE.4)CONST=CONST+V(I)'V(I)
END 00
VB=DSQI\T(TIIREE'CONSTjTWO)
RETURN
END
SUBROUTINE VTRANVI(V,V1,VTV1,NISS)
c*******************************
C . . . Inner product of 2 symmetric 2nd order tensor s
IMPLICIT REAL'O (A-H,O-Z)
REAL'O V(6),V1(6)
C
VTV1=0.
00 I=l,NISS
VTV1=VTV1+V( I) 'VI (I)
END 00
RETURN
END
(j (MPa)
SOO
400
300
200
0.6
100 0.4
0.2
0 0
0 5 10 15 20 Fig. 4.18. Numerical simulation of a pure
e ('!.) tension test at the microscale
4.3 Locally Coupled Analysis 197
and find for each numerical test the set of strains (E 1\ , E22 ) which corresponds to
crack initiation. You will obtain the so-called "damage forming limit curve" in
plane strain, as shown in Figure 4.19.
Perform calculations for the case of an elastic tension case at the mesoscale: Ell
given with E22 = E33 = -VE I \' with a plastic threshold a s = 200 MPa, and plot
the result at the mesoscale (an elastic-brittle material) and at the microscale. as in
Figures 4.20 and 4.21.
300 *
200
100
200
*
100 o
0.10
0.05
O--'r-~---r--~~O~~-r--~~ Fig. 4.21. Stress strain and damage at the
o 2 3 4 5 microscale when pure tension is applied
at the mesoscale
Consider a pure tension case at the microscale in which E II is the main input, the
other components E22 and E33 being calculated by DAMAGE 90 as
The cyclic strain hardening is taken into account by using piecewise perfect plas-
ticity according to the cyclic stress strain curve:
Perfonn some numerical simulations to obtain the number of cycles to crack ini-
tiation as a function of the stress amplitude, and plot the results together with the
number of cycles No to crack nucleation corresponding to P = PD (Figure 4.22).
Consider the case Ell = 3.5% (~E = 7%). The number of cycles to failure is
40. Compare to an hypothetical case of pure tension at the mesoscale where the
4.3 Locally Coupled Analysis 199
input is
E" = 3.51O- 2 , E22 = E33 = -VEil'
Due to the effect of the triaxiality the number of cycles to failure would be only 8.
Plot the graphs of stresses, strains and damage for the pure tension case at the
microscale as in Figure 4.23.
t. E11 ("!o)
8.0 \
\
\
\
4.0 \
\
\ No NR
\
2.0 ,
"-
.......
' ,,-
1.0 . . . . . .......
-. -
- ....,---
~""'
0.5 t - - - - . , - - - - , - - - , - - - . - - - - , - - - -
10 10~
N Fig. 4.22. Fatigue rupture curve
e" = 0.47% during 3088 cycles and E" = +0.425% after ~ N R 52400
cycles;
Ell = 0.425 % during 43828 cycles and Ell = +0.47% after ~ N R 50000
cycles;
Check that for
Multiaxial fatigue
You may also play with DAMAGE 90 to obtain contours of the number of cycles to
failure in biaxial fatigue as it is shown for the same aluminum alloy as previously:
- for plane strain in Figure 4.25:
cil = x, c22 = y, c33 = OJ
- for tension and shear in Figure 4.26:
c II = x, c 12 = y, all other components = OJ
"
0.04 (J",,(MPa)
600
400
0.02
200
0
N 0
"
-0.02 -200
-400
-0.04
-600
-0.06 -800
0 10 20 30 40 -0.06 - -0.02 0 0.04
(J"eq (MPa)
SOO 5
p
400 4
300 3
200 2
100
D
0 0
0 10 20 30 40 50 0 10 20 30 40 so
N N
Fig. 4.23. Very low cycle fatigue simulation in tension at the microscale
4.4 Fully Coupled Analysis 201
~
NR
2
1.0
8"
0.8
N
0.6
0.2
O+----,r----.-----,----,---~~---
o 0.2 0.4 0.6 0.8 1.0
Fig. 4.24. Accumulation diagram for two-level simulation in tension at the mesoscale
dE22 %
2
.7
.6
.5
.4
NR=7800 cycles
.3
.2
.1
L\12 %
-2-
.4
.3
.2
.1
.1 .2 .3 .4 .5 .6 L\1I %
-2- Fig. 4.26. Biaxial fatigue envelope for the case of ten-
sile and shear strains imposed
The classical way to perform this calculation is to implement the coupled constitu-
tive equations from Section 2.1.4 and the kinetic damage law from Section 3.1 in
a finite element code as explained in Section 4.1.3. The difficulty for convergence
of iterations and for the computer time of calculations are of the same order as
for classical elastoplastic analysis. Such routine is available in the French code
SYSTUS and possible with the American code ABAQUS.
What makes this method interesting is that it allows us to take into account pre-
existing strain hardening and damage as initial fields in some part of the structural
component to be analyzed. It is always of great interest to evaluate the remaining
strenght or lifetime of a component after an accident which has induced plastic
deformation and damage.
Three quantities may be introduced as material initial conditions at the Gauss points
of the finite element meshes:
- Ro, the isotropic strain hardening stress as an initial value for its differential
constitutive equatio.n,
Coupled
elasto-plastic and damage
constitutive equations
or, in the absence of damage, Po' an initial value of the accumulated plastic
strain since in this case, there is a one-to-one relationship between p and R:
R = Roo[l- exp( -bp)] ;
- X{?, the kinematic strain hardening stress tensor as an initial value for its dif-
ferential constitutive equation,
XB = I' [~Xooifj(I-D)-XB~];
- Do an initial value of damage for its kinetic law,
b Y.
= SP.
One more quantity which often exist as a characteristic of the structure, but not
of the material, is a state of residual stresses, after welding, for example. A
self-equilibrated residual stress field can also be introduced at the Gauss points
as an initial condition:
o
aij'] = oIn S }
a?j such that { Figure 4.1
a?j,n j = 0 on 8S
The first three material quantities are not difficult to introduce in a finite element
code, but the main problem is to identify them. As far as isotropic strain hardening
and damage are concerned, if the damaged structure is available for testing, one
can evaluate Ro and Do, at least on the surface, by microhardeness measurements
as explained in Section 1.3.3. As far as kinematic strain hardening is concerned its
initial value may be evaluated by means of phenomenological considerations(!!)
or neglected; the latter is not an approximation for the stabilized cycle which does
not depend upon the back stress X{?:
Note that those initial conditions may be considered in the uncoupled analysis as
a constant reduction of the elasticity modulus and the strain hardening variables.
It is out of the scope of this course to describe in detail the finite element procedure
for elastoplasticity or viscoplasticity coupled with damage calculations. The set
of constitutive equations of Section 2.1.4 together with the kinetic damage law
identified in Section 3.1.4 may be implemented in any classical finite element code
having plasticity routines. Best results are obtained with the algorithm described
briefly in Section 1.3 of the present chapter.
The practical case considered here schematically is an underground gas pipeline
which has been damaged by the shock of a machine digging in the soil. The
pressure inside the pipe varies periodicallly twice a day due to gas consumption. It
induces fatigue of the material especially where the shock occured. The question
is: What is the remaining strength or lifetime of the damaged pipe before it breaks?
Obviously, it is of primary importance to decide when and how to make a repair.
204 Analysis of Crack Initiation in Structures
The pipe and the loading being considered are defined in Figure 4.28.
A simple calculation of the undamaged cylinder unconstrained at the ends by
the expansion joints gives the state of stress everywhere. Assuming the material
to be elastic, one can deduce the displacements everywhere and particularly at
the boundaries A and B of the substructure containing the damage that will be
studied in more detail. Figure 4.29 shows the geometrical model submitted to a
one-dimensional displacement within the hypothesis of plane strain. Note that this
substructuring neglects the redistribution of stress due to plasticity that will occur
at the root of the notch as a result of the loading 71(t).
The material is a low-carbon steel. Its elastoplastic and damage parameters have
been determined from several tension tests according to procedures described in
Sects.2.1.5 and 3.1.5:
- E, v for elasticity;
- a y' Roo b, X oo ' 'Y for plasticity;
- S, cPD ' Die for damage together with af and au
More interesting is the strain hardening and the damage induced by the shock in
the idented zone; these must be considered as initial conditions, as explained in
Section 4.4.1. These damage and strain hardening fields have been identified from
microhardness experiments on a polished section of a plate being subjected to a
shock similar to that which damaged the pipe. The method used is described in
Section 1.3.3. The results are shown in Figure 4.30.
A I B p
o~----------------------------~
N
u
42mm
1----- - - - - - - - ------------
The values of Ro and Do are introduced at each Gauss point of the finite ele-
ment meshes represented in Figure 4.30. The rest of the meshing is indicated in
Figure 4.29. The results of the structural calculation corresponding to a mono-
tonic loading performed in order to determine the remaining strength, are shown
in Figures 4.31 and 4.32.
From the evolution of the damage in Figure 4.31, the displacement which will
cause crack initiation is u c:::' 0.12 mm. It turns out to be about 50% of the elastic
displacement corresponding to the load of the pipe when the von Mises equivalent
stress is equal to the yield stress ay. This particular damage has reduced the strenght
of the pipe by a factor of about 50%.
If the initial strain hardening and damage are not taken into account the damage
evolution shows a critical displacement of about 0.20 mm, that is, an error of 66%.
As far as crack initiation by fatigue due to the cyclic pressure in the pipe is con-
erned, the same calculation is performed with a periodic displacement as input.
To avoid a considerable amount of calculation, the procedure of integration by
damage increments (instead of time or cycle increments) described in Section 4.13
is used. For the case under consideration, the number of cycles to crack initiation
of the indented substructure (or the pipe) was about 100 times less than the fatigue
life of the initial pipe. As the nominal number of cycles to failure of the undam-
aged pipe line was of the order of 105 , this corresponds to 1000 cycles, that is,
more than a year to organize the repair without cutting off gas delivery. Thank
you.
206 Analysis of Crack Initiation in Structures
o'eq(MPa) 0
600 1.0
0.8
400
0.6
0.4
200
0.2
u u
0
0 0.02 0.04 0.06 0.08 0.10
mm
0
0 0.02 0.04 0.06 0.08 0.10 mm
Fig. 4.31. Evolution of the von Mises equivalent stress and the damage at the critical point M *
o'eq (MPa)
600
u =0.12mm
400
u=O.OSmm
200
O.S
0.4
u=0.12mm
0.2
z
OL-------~~----r_------._--_.
o 2 3 (mm)
Fig. 4.32. Von Mises equivalent stress and damage along the z axis
Then, to cover the whole process of fracture of a component, the following steps
may be worked out (Figure 4.33):
taken into account by their mean effects in the identification of the material
parameters. But for materials such as ceramics or cast iron. The large size of some
defects imposes the introduction of a statistical effect. Furthermore the only way
to explain the size effect on fracture or the high scatter in fatigue failure for
example is to take into account the probability of having defects of a given size in
the tested sample or structure.
d2
D~[2 .
Then, the statistical characteristics of the surface of initial defects may be directly
translated into the statistics of initial damage Do provided the definition of the size
a representative volume element be possible. Practically, this size I may be the
mean distance between the centers of adjacent microcracks, that is the inverse of
their linear density. These defects may be randomly distributed in size, orientation,
shape ... but we restrict ourselves to the distribution in size only given by the
probability density f(Do} of the equivalent initial damage Do
Prob (Do :s;; Do :s;; Do + dD o) = f(Do} dD o .
The function f(Do} may be deduced from measurements of defects through obser-
vations with a microscope. We shall see that some statistical characteristics may be
obtained more easily from mechanical characteristics at measoscale.
The probability density characterizes the distribution of the defects or the initial
damage in a structure but it is also the probability of having an initial damage of
a given value in each representative volume element.
J
(a1;)2 = (a1;(Do) - a1;)2 f(Do)dDo,
o
1
(a1;)2 = a; J(150 - Do)2 f(Do) dD o ,
o
Do = D (Do -15 0 )2 f(D o)dD o J/ 2 being the standard deviation of the initial
damage distribution.
As an example, consider an uniform distribution of initial damage
212 Analysis of Crack Initiation in Structures
the mean value and the standard deviation of the failure stress are
where (Jf is the fatigue limit of the material, S the damage strength, Rv the
triaxiality function.
Assuming a proportional loading for which Rv = const.
D(t) = Do + 2is(J2
Rvceq(t)'
Do becoming a random variable with a probability density f(Do), the mean value
of the accumulated stain to crack initiation is
or with
4.5 Statistical Analysis with Microdefects 213
(e:q)2
Dc
= J(e:q(D o ) -
o
e:q) f(Do) dD o ,
2
-
where Do =
[DC! (Do _DO)2 f(Do)dDoJ1 /2 is the standard deviation of the initial
damage distribution.
the mean value of the number of cycles to failure NR and its standard deviation NR
immediately follow
- NR -
NR = Do
Dc-Do
explains a well-known experimental result which shows that the standard devia-
tion of the number of cycles to failure NR , or the scatter in fatigue, increases with
the number of cycles to failure NR . An example is given in Figure 4.35.
Another interesting consequence is related to the identification of initial damage
or defects of a given material for which the Woehler curve is known by the mean
value of the number of cycles to failure NR and its standard deviation NR .
Solving the above equations for Do and Do gives the mean value and the standard
deviation of initial damage
214 Analysis of Crack Initiation in Structures
Experiments
2 .--~...,..,-~,---~~-,.-j -- PF =50 %
-- - --PF = 90 %
\
'- ,
......... P = 10 %
.... F
1.5 ___ ~----------r---~
----------- -
-E
CI3
1
. .... ..
.. ,....
.
............... ..... ......... ......... _....
zo 0.5 L.4~~~.........L-
5 ~~~'--'--'----6~~~....... 7
10 10 10 10
Number of Cycles to Failure
Fig. 4.35. Woehler curve of a sa cast iron with its scatter (from A. S. Beranger)
2 A
- aff, -
Do = D ---NR
C ES '
at ~ au(l- Dc) = aR
The probability of survival of the structure considered as the set of the R VE links
is given by the product
but
!
where PFO = J f(Do} dD o ,
Dc
1
then Ln(1- PF) = - JLn(1- PFo)dV,
Vo v
This formula explains the size effects: Vo being fixed by the nature of the material,
if the volume of the structure V increases, the probability of failure PF also
increases.
The well-known Wei bull formula may be derived for a three-dimensional state of
stress with two additional hypothesis.
with
with
then
4.5 Statistical Analysis with Microdefects 217
v
- --+ 00 =>
Vo
The main difficulty in a volume effect analysis is the identification of the elemen-
tary representative volume Vo. Practically, Vo maybe taken as the volume of the
uniform stressed part of the specimens used to identify the constitutive equations
of the material under consideration, but it is difficult to avoid a certain degree of
subjectivity!
A better identification would be to determine the probability density of defects
f(Do} which allows for the calculation of PFo if Dc (1*) is known. Then from tests
on a real structure PF may be derived and Vo calculated from
PF = 1- (1- PFo}v/vo .
>-
::l
'"
~o
DM =0,1 ...,
~ "-"-'-"-"-' .. _.. _.. _.. _..-.. _..-.. -.. _.. _. --------_.._.._._.._.._.. _..-.. _.._.._.._.._.._.._.._.'-"-"-"-"-"-"-"-'"
-~.--.-.-- - - . - - - - - - - - . - - - _ _ _ . _ _ _ _. __ 0' ._ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ __ ~
~
0.6 D .. =0,5 ~
z
a
0.4
0.2
V
o ~1--~4-~-4--~4--+~--~4-~-4--~4-~-4--~4-~-4--+-4--+~--+-~-+~~+-~
Vo
o 5 10 15 20 25 30
Fig. 4.36. Volume effect on the mean damaged equivalent stress at fracture
4.5 Statistical Analysis with Microdefects 219
bility of failure for the same distribution of defects is larger than for a uniform
stress field involving the same strain energy.
To model this effect, it is convenient to introduce the concept of effective volume
V eff defined by an equivalence in elastic strain energy as it governs the phenomenon
of damage.
If We is the elastic strain energy density
W~AX its maximum value at the stress concentration point of the structure,
Jwe(M)dV = VeffW~AX,
V
or
Then, the probability of failure formulas are written with the effective volume
calculated by an elastic calculation instead of the real volume of the structure,
PF = 1- exp ( - Va [O'~AXJm)
Veff
~ .
L
..
Fig. 4.37. Beam in bending fatigue
220 Analysis of Crack Initiation in Structures
12M
L1e = Ebh 2 .
- The number of cycles to failure from Section 4.5.3 is
show that:
- The probability of crack initiation is
DM
PFo = J f(Do)dDo ,
(Do)e
NR aJL1e
(Do), = D, - ES .
Taking into account the effects fo volume and stress heterogeneity linear with X2
show that:
- The effective volume is
Lbh
Veff = -3-
Notice that the effective volume is reduced to Lbh/6 if pure unilateral conditions
are considered that is no damage in compression.
- The volume ratio is
Veff 1
Vo 3'
4.5 Statistical Analysis with Microdefects 221
JE
Bibliography
Strain Localization and Size Effects Due to Cracking and Damage, CNRS-
NSF Workshop, Cachan, France, 1988.
Fracture and Damage of Concrete and Rock, FDCR 2, Vienna (Austria), 1992
Inelastic Deformation, Damage and Life Analysis, Int. Symposium. Hawal (USA),
1995
Micromecbanics of Plasticity and Damage of Multipbase Materials, IUTAM Sym-
posium, Paris (France), 1995
Mecbanical Bebaviour of Damaged Solids, Mecamat Workshop, Fontainebleau
(France), 1995
Author and Subject Index
Ladeveze, 64, 83
Gauss point, 166 Lame's constants, 79
generalized standard materials, 47 Leckie, viii, ix, 60, 79
Germain, viii Lee Hao, vii
Geymonat, 86 Legendre-Fenchel, 85
Gibbs specific enthalpy, 83 Lemaitre, v, 11, 13, 39, 41,95,176
Goodman's diagram, 139 Lienard, 20, 21, 168
Gunawardena, ix Lin-Taylor hypothesis, 81 , 114, 136, 177
Gurson model, 72, 74 linear accumulation, 142
linear creep fatigue interaction, 135
Lippmann, v, vii, viii, ix
listing of DAMAGE, 90, 185
Hayhurst, 60, 206 loading function , 48 , 49
Heaviside, xv, xvi loading-unloading condition, 49
Helmholtz free energy, 41 localization, 86, 176
Herath, ix locally coupled analysis, 176
high cycle fatigue, 7, 27, 135, 197 low cycle fatigue, 5, 27, 127, 197
high-alloy steel, 19 lower bound, 171
Hild, 5, 128,209 Luder's bands, 89
homogenization, 69
H ult, viii, 90
Hutchinson, viii, 208 Macauley, xv
hydrostativ stress, 65 macrocracks, 205
macroscale, 4
Manderscheid, ix
identification, 57, 109, 110 Manson-Coffin law, 131
implicit integration, 166 MAR M 509, 36
implicit scheme, 178 Marquis, viii, ix, 39
IN 100 superalloy, 32 master curves, 120
inclusions, 4 material parameters, 57, 98, 104, 109, 152
incremental linearization, 165 Maxwell's compatibility equation, 87
initial conditions, 155 Mazars,24
initial damage, 179, 202 McClintock, 74
initial plastic strain, 179 McMeeking, viii
initial strain hardening, 202 mean stress, 138
input data, 180 measurements, 19
instability, 86 mechanism of damage, 10
interface, 66, 208 mesoscale, 4, 116
intergranular decohesions, 5 metal forming, 119
intergranular microcracking. 7 metals, 1, 3, 117
internal variables, 40 micro crack closure effect, 107
isochronous creep damage curves, 126 micro-meso element, 68, 72
isotropic damage, 39 microcavities, 72
isotropic hardening, 14, 40, 42, 51 microcrack, 2, 19, 68
Author and Subject Index 227
Weibull, v, 216
wine selection, 35