Issue 31
J. Xavier et alii, Frattura ed Integrità Strutturale, 31 (2015) 13-22; DOI: 10.3221/IGF-ESIS.31.02
/ C P
in which
is the compliance. From the Timoshenko beam theory and Castigliano theorem an expression for the
E ) using an
compliance of the DCB specimen can be obtained. This equation can be solved for the flexural modulus ( f
) a as [4]
0 ( ) C and the corrected initial crack length 0 (
initial compliance
1
3
a
a
12
8
0
0
E C
(2)
f
0
3
BhG
Bh
5
LR
where accounts for root rotation effects and can be determined from finite element analysis [4], and G LR
is the shear
a , which is considered to account for the FPZ effect at the crack tip as
moduli of the material. The CBBM is based on eq
a
eq a a
[7, 10]. Finally, the application of CBBM to the DCB test yield the following expression for
given by:
FPZ
the strain energy release rate in mode I (resistance or R –curve) [4]
2 B h E h G 2
2 2 6 1 5 f a P G eq 2 LR
(3)
I
It is worth noting that this procedure is less sensitive to experimental errors. This is supported by the fact that the measurement of the crack length during the fracture test is not required. Besides, the inherent elastic properties variability among specimens is taken into account by computing a flexural modulus for each P curve Eq. (2). In mode I loading, strain energy release rate ( I G ) and CTOD ( I w ) can be related by the following expression [9]
w
I
( )d w w
G
(4)
I
I
I
I
o
I I f( ) w ) can then be directly obtained by differentiating the above equation
The cohesive law (
I G w w
I
( )
(5)
I
I
This data reduction scheme, however, requires the accurate evaluation of the I I f( ) G w relationship. Moreover, a suitable differentiation algorithm must be used to avoid noise amplification in the reconstruction of the constitutive cohesive law. In order to solve Eq. (5), it is proposed here to fit the I I G w data by a continuous function described by the following expression (logistic function)
1 I,0 1 ( / ) p A A w w 2 I
G
A
(6)
I
2
2 A , p and
1 A ,
A parameter must
I,0 w are constants to be determined in least-square sense. In this function, the 2
where
2 w A G G I lim
provide an estimation of the critical strain energy release rate:
.
Ic
I
Figure 2 : Finite element model of the DCB test (cohesive law, mesh and boundary conditions).
15
Made with FlippingBook Ebook Creator