1976 CCR Hillerborg p773 PDF
1976 CCR Hillerborg p773 PDF
1976 CCR Hillerborg p773 PDF
(Communicated by Z. P. Bazant)
(Received August 24, 1976)
ABSTRACT
Une m6thode est present,e, par laquelle la m6chanique des ruptures est
i n t r o d u i t e dans l'analyse des ~l~ments f i n i s ~ l ' a i d e d'un module, oO
les contraintes sont suppos6es d'op6rer sur les c6t6s d'une fissure
tant que cette f i s s u r e est ~ t r o i t e .
Cette hypoth~se peut etre consid~r~e comme un moyen d'exprimer l ' a b -
sorption Gc d'6nergie en usant l'approche de l ' ~ q u i l i b r e d'6nergie.
Cette hypoth6se est aussi j u s t i f i 6 e par les r 6 s u l t a t s des essais de
tension.
Pour en prouver la v a l i d i t Y , cette m~thode a ~t~ appliqu~e au fl6chisse-
ment d'une poutre non arm~e et f o u r n i t une explication de la diff6rence
entre la r~sistance au moment de flexion et la r~sistance ~ l ' e f f o r t de
tension, ainsi que de la variation de la r~sistance au moment de flexion
en fonction de la profondeur de la poutre.
773
774 Vol. 6, No. 6
A. H i l l e r b o r g , M. Mod4er, P-E. Petersson
Proposed approach
There are many methods to choose from f r a c t u r e mechanics, e.g.
the stress i n t e n s i t y f a c t o r approach
the energy balance approach
the " s t r i p - y i e l d " model according to Dugdale
the cohesive force model according to Barenblatt.
The d i f f e r e n t methcds are known to give coherent r e s u l t s .
In the stress i n t e n s i t y f a c t o r approach the stresses near the
crack t i p are studied. These stresses t h e o r e t i c a l l y a ag_j~roach i n f i n i t y
at the crack t i p according to the expression o = K/V2~r, where r is
the distance from the crack t i p and K is a c o e f f i c i e n t , the stress i n -
t e n s i t y f a c t o r , depending on the load, the crack dimensions, etc. When
K reaches a c r i t i c a l value Kc, the crack propagates.
The stress i n t e n s i t y f a c t o r approach has been used a great deal
in FEM analysis. The d i r e c t method requires a FEM mesh with very small
elements close to the crack t i p , which l i m i t s i t s a p p l i c a b i l i t y to
Vol. 6, No. 6 775
FRACTURE MECHANICS, ANALYSIS, CONCRETE
0
and c o r r e s p o n d s to t he area between the curve and the c o o r d i n a t e
a x i s ' in Fig. 3.
776 Vol, 6, No. 6
A. H i l l e r b o r g , M. Modeer, P-E. Petersson
Cr
Fig. 3 Assumed v a r i a t i o n of
stress o with crack
ft width w, general case
0
I
w! W
71dw = G (I)
C
0
which means that the energy absorbed per newformed u n i t crack area is
the same as in the energy balance approach. The model of Fig. 2 may
thus be looked upon as a way of expressing the energy balance approach.
At the same time the assumption of Fig. 2 may be looked upon as
a r e a l i t y . Stresses may be present in a microcracked zone as long as
the corresponding displacement is small. This has been c l e a r l y demon-
strated in tension t e s t s , using a very r i g i d t e s t i n g equipment, e . g . ,
by Evans and Marathe / 4 / ; cf. Fig. 5.
By the a p p l i c a t i o n of the proposed model the curve ~(w) may be
chosen in d i f f e r e n t ways, e.g. according to Figs. 4a, b or c, which
a l l show simple mathematical r e l a t i o n s . For t y p i c a l y i e l d i n g materi-
a l s , l i k e mild s t e e l , Fig. 4a seems to be the best choice. I t corre-
sponds e x a c t l y to the Dugdale model with f t = ~y.and ~I = COD at i n i -
t i a t i o n of crack growth. The d i s c o n t i n u i t y may glve r l s e to some prob-
lems by the a p p l i c a t i o n in FEM, but they are not serious.
~dw = ftwl/2
0
or from ( I ) ,
wI = 2Gc/f t (2)
For ordinary concrete Gc/f t seems to be of the order 0.005 -
Vol. 6, No. 6 777
FP~ACTURE MECHANICS, ANALYSIS, CONCRETE
O.Ol mm, cf. / I / , and thus wI of the order 0.01 - 0.02 mm. In the
a p p l i c a t i o n we f u r t h e r assume that the concrete is l i n e a r - e l a s t i c
u n t i l f t is reached.
Fig. 5 shows a comparison between our assumptions with Gc/f t
= 0.01 mm, ~ / f t = I0 000 and a tension t e s t from / 4 / with a gage
length of 1 (25 mm). This corresponds to a t h e o r e t i c a l average
elongation over thergage length when ~ reaches O, i . e . w = Wl, of
2.0.01/25 = 800-I0 - . The assumptions seem to agree reasonably with
the t e s t r e s u l t . A lower value of E/f t would have improved the agree-
ment, but from the point of view of the energy balance approach the
E-value corresponding to unloading is most important and t h i s j u s t i -
f i e s the choice E/f t = I0 000.
ft
Fig. 5 A test result
from a t e n s i l e
t e s t according
to Evans & Mara-
the / 4 / , compared
to a corresponding
assumed r e l a t i o n
by the analysis
or , , , ~ .
0 200 ZOO 600 BOO Strain x 106
(1" gage length)
C M
~
Fig. 6 Bent rectangu-
4J
42 l a r beam and
41
40 a corresponding
J9 FEM representa-
3'8
J7 tion
.16
35
33
32
3t
1,4
NIt
Mo
Fig. 7 Calculated bending moment
M versus crack depth
0 crack
i depth
i node
point
As the r e a l t i o n
EGc = K2
C
holds for plain stress and approximately for plain s t r a i n , we may also
write 2
I c = (Kc/ft) (4)
Bendin~ strength
Tensile strength
Fig. 8 Theoretical v a r i a -
t i o n of r a t i o be- 2,0 -
tween bending and
t e n s i l e strength !
with beam depth H
and I c = ( K c / f t ) ~ =
: EGc/f~ 1.5 J
.~oO t
1.0
~,E ~
T .a~ ~....~
0.5
c~ j
,, .....<
r
Fig. 9 Assumed d i s t r i -
bution of s h r i n k - o L
age s t r a i n s o 2 3 4 H/I c
Relative bending
strenq th
1,0 - e- : = ~ ,
Fig. I0 Test r e s u l t s of
bending s t r e n g t h
0.5 - - versus beam depth,
summarized by Mey-
er / 3 / , compared
to t h e o r e t i c a l cur-
ve f o r 1 = I00 nTn
C
0
0 100 200 300 500 H mm
Meyer 1963 ~8 V a L u e s
o wOt
o dried O.O0"J'~ * H 2 P,~urs
~' drded O,OtS. H 2 h ~ r s
1,5
I Fig. I I Test r e s u l t s o f
I
i i
,.o bending s t r e n g t h
versus beam depth
f o r one q u a l i t y
o f concrete, tested
wet and d r i e d in
45% RH, compared to
o, b L
t h e o r e t i c a l curves
f o r 1 = I00 mm
C
, 1 i"
o
!
i ! II-
loo 2OO 300 H mm
Vol. 6, No. 6 781
FRACTURE MECHANICS, ANALYSIS, CONCRETE
Conclusion
The proposed method of combining fracture mechanics and f i n i t e
element analysis seems to y i e l d r e a l i s t i c results regarding crack f o r -
mation and propagation as well as regarding f a i l u r e even i f a coarse
element mesh is used. This opens up the p o s s i b i l i t y of studying comp-
licated problems with a limited amount of computer work.
References