Determination of fiber/matrix interface debond growth parameters from cyclic loading of single fiber composites Andrejs Pupurs 1,2 , Janis Varna 1 , Povl Brøndsted 3 , Stergios Goutianos 3 1 Luleå University of Technology, S-97187 Luleå, Sweden 2 Swerea SICOMP, SE-94126 Piteå, Sweden 3 Technical University of Denmark, Risø Campus, DK-4000 Roskilde, Denmark 6 th International Conference on Composites Testing and Model Identification (CompTest 2013) 22-24 April 2013, Aalborg, Denmark.
18
Embed
Determination of fiber/matrix interface debond growth ... · E r G f m mech th m z f f z II ' fH fDD 2 1 2 d Q d E k f z m f zr m » ¼ º « ¬ ª f QQ f mz m f m z f f r f zr th
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Determination of fiber/matrix interface debond growth parameters from cyclic loading of
single fiber composites Andrejs Pupurs1,2, Janis Varna1, Povl Brøndsted3, Stergios Goutianos3
1Luleå University of Technology, S-97187 Luleå, Sweden 2Swerea SICOMP, SE-94126 Piteå, Sweden
3Technical University of Denmark, Risø Campus, DK-4000 Roskilde, Denmark
6th International Conference on Composites Testing and Model Identification (CompTest 2013) 22-24 April 2013, Aalborg, Denmark.
2
Fatigue damage in UD composites
Failure of composite similar to static failure. Random fiber
breaks and failure in few load cycles.
N=1
Applied strain level is too low (below fatigue
limit). N=∞
Random fiber breaks during first cycle, N=1.
Further damage development in form of
growing fiber/matrix interface debond cracks.
After certain amount of load cycles, aligned matrix
cracks appear and propagate as debond
cracks. Bridged cracks followed by new fiber
cracks
2
Failure of composite similar to static failure. Random fiber
breaks and failure in few load cycles.
N=1
Applied strain level is too low (below fatigue
limit). N=∞
Random fiber breaks during first cycle, N=1.
Further damage development in form of
growing fiber/matrix interface debond cracks.
After certain amount of load cycles, aligned matrix
cracks appear and propagate as debond
cracks. Bridged cracks followed by new fiber
cracks
c
FL
f c
FL f
3
Objective
• To focus on fiber/matrix interface debond growth in tension-tension fatigue
• To use fracture mechanics (energy release rate) for the debond growth analysis
• To calculate the energy release rate and to determine its dependence on debond length and to describe it by simple fitting function
• Using experimental data and the calculated energy release rate to evaluate whether a power law can describe the debond length growth rate in cyclic loading
• To determine material constants in this power law, if it is applicable
4 Debond growth in tension-
tension fatigue
ld
mIIGB
dNdA
ldn ld / rf
mII
dn GBdNdl *
B* B 2 rf2where
dA – increase of the crack surface area
GII G( max ) G( min )
max 0min 0
dN – increase of number of applied load cycles
mKBdNdAmax
min
Paris law:
)(Nfld
B, m – material constants
K – stress intensity factor range
5
Determination of power law parameters
• Proper values of Paris law parameters m and B can only be found if experimental data of ld=f(N) are available.
• Experimental measurements of debond growth in fatigue is complicated for UD composites.
• Single fiber (SF) composite with transparent matrix, is a convenient tool to measure debond crack growth in cyclic loading.
• Analytical and FEM calculations of GII for SF composites are necessary to account for the stress state.
• Proper vam and B cexperimeavailable.
• Experimegrowth in composite
• Single fibtranspareto measucyclic load
• Analyticalfor SF comaccount fo
mII
dn GBdNdl *
6
Concentric cylinder model
Lf
ld
Lf
z z
fiber break
debonded fiber region
z
T
F M
7
Long and short debonds
• Long debonds • Debond growth in self-
similar manner • Analytical model for GII
calculation
d) Plateau stress-state in bonded fiber region.
c) Debond tip region with stress singularity.
b) Plateau stress-state in debonded fiber region.
a) Very complex stress state at the fiber break.
Stress state regions
• Short debonds • Strong interaction
between stress states in regions a) and c)
• Numerical methods for GII calculation
ld
F M
8
Analytical solution. Energy release rate for long debonds
• J.A. Nairn, Y.C. Liu, On the use of energy methods for interpretation of results of single-fiber fragmentation experiments, Composite Interfaces, vol.4, pp. 241-267, 1996.
2
4Tkk
rEG f
zmthmechmf
fz
II
21
2dQ
dE
k fz
mf
zrm f
zmm
fzm
frf
z
fzr
th dQ
dE
k2
12
m
mf
r
fr
EEd
111 1
2
22 d
Ed f
z
fzr
1
221dE
Q fz
fzr
Em = 3 [GPa] Em = 3.5 [GPa] Em = 4 [GPa]
m = 0.3 0.997427 0.997011 0.996598
1.014151 1.016443 1.018716 mk
thk
9
Numerical calculation of energy release rate for short debonds
• Energy release rate calculation using Virtual Crack Closure Technique (VCCT)
• Calculations performed using FEM software ANSYS
uz z
r
fiber break
sym
met
ry a
xis
T
oute
r bou
ndar
y
GII (ld ) limdld 0
12dld
uldld
ld dld
z dld rzld z dz
Lf
Symmetry b.c.
10
Magnification of energy release rate for short debonds
BEM results taken from: • E. Graciani, V. Mantič, F. París, J. Varna, Numerical analysis of debond propagation
in the Single Fibre Fragmentation Test, Composites Science and Technology, vol.69, pp. 2514-2520, 2009.
1.0
1.4
1.8
0 5 10k m
*
ldn
FEM
BEM
1.0
1.5
2.0
2.5
0 10 20 30
k j*
ldn
CF, k m
CF, k th
GF, k m
GF, k th
11
General expression for energy release rate
• Long debonds:
2
4Tkk
rEG f
zmthmechmf
fz
II
km km* (ldn)km kth kth
* (ldn)kth
)()()( 0**
dndnthdnm lklklk
30 ;1
301 ;1 321
20
dn
dndndn
l
llaala
k
• Short debonds:
• Fitting the whole range:
GII (ldn)rf Ez
f k02(ldn)4
( mechmax )2 ( mech
min )2 2( m zf ) T( mech
maxmechmin )
• Energy release rate change in one load cycle:
12
Experimental measurements
• Fiber was pre-stressed with 7g of weight suspended at each end of fiber • First a static load was applied until first break was initiated in fiber ( 1st) • Tension-tension cyclic loading with frequency f=2Hz, R=0.1 was applied
under load control. Measurements of debond length ld were performed after selected number of load cycles.
• Two different values of maximal fatigue strain level were studied: Samples A, B: max = 1.76%; Sample C: max = 1.32%. The strain levels approximately correspond to 80% and 60% of the stress 1st respectively.
[GPa] [-] [1/°C] [GPa] [-] [1/°C]
70 0.20 4.70·10-6 3.17 0.33 91.00·10-6
Ezf
zrf
zf Em m m
1st fiber break
Glass fiber/Epoxy
13
Experimental results. Debond growth in cyclic loading
• The strain energy release rate in Mode II for debond growth was analyzed combining analytical solution for long debonds and FEM solution for short debonds.
• Using the quantified debond length versus number of cycles data it was shown that the power law with respect to the strain energy release rate change is applicable for debond growth characterization in tension-tension fatigue.
• Simulations showed that the obtained parameters give acceptable predictions for cases, when the debond grows as well as when it does not grow.
• Being material properties the power law parameters determined in this study using single fiber composites can also be applied for the case of UD composites made of the same material system
17
Related papers
• A. Pupurs, S. Goutianos, P. Brøndsted, J. Varna, Interface debond crack growth in tension-tension cyclic loading of single fiber polymer composites, Composites Part A, vol.44, pp. 86-94, 2013.
• A. Pupurs, J. Varna, FEM modeling of fiber/matrix debond growth in tension-tension cyclic loading of unidirectional composites, International Journal of Damage Mechanics, In Press, 2013. (available online)