A cohesive damage-friction interface model accounting for water pressure on crack propagation G. Alfano * , S. Marfia ** , E. Sacco ** * Dipartimento di Scienza delle Costruzioni Università di Napoli, Via Claudio, 21 - 80125, Napoli, Italy e-mail: [email protected]** Dipartimento di Meccanica, Strutture, Ambiente e Territorio Università di Cassino, Via G. di Biasio 43, 03043 Cassino, Italy e-mail:[email protected], [email protected]Abstract. Aim of the work is the study of crack propagation in concrete constructions, such as dams, taking into account the water pressure effects and the damage-friction evolution with unilateral contact. With this aim, an interface model, based on the cohesive fracture, is developed. In particular, a damage-friction model based on a new multiscale approach is proposed for the interface, which is able to simulate crack propagation in mode I, mode II, and mixed mode. The criteria for crack initiation and propagation and the closure and reopening of the crack are considered for the interface. The friction and the dilatancy are characterized by a decreasing rate of friction and dilatancy angles, respectively, with the possibility of the dilatancy recovering during cyclic loading. The water pressure effect is taken into account by considering a further static pressure acting on the real crack and in the process zone. The value of this pressure is set as a function of the crack opening displacement and it varies from a prescribed initial pressure to the external water pressure. A numerical procedure is developed in order to integrate the evolutive laws governing the interface model. The interface model is implemented in finite element codes; some numerical simulations and a comparison with experimental data are performed in order to validate the proposed model and to investigate the influence of the water pressure effects on the crack propagation. Finally, the mechanical response of a concrete dam is studied. Keywords: Cohesive crack, interface, damage, friction, dilatancy, water pressure, concrete dams. 1 INTRODUCTION In the last decades, concrete dams have been the object of an increasing attention from the international scientific community. In this respect, sophisticated mechanical models, suitable numerical procedures and appropriate design methods have been developed for the evaluation of the nonlinear response of concrete dams. In particular, the presence of joints generated by a step-wise successive construction of layer concrete blocks and between the dam and the foundation can play a fundamental role in the overall dam behavior; in fact, these interfaces are the potential sites of crack growth which can lead to the weakening of the structure and to the penetration of water that exerts uplift pressure [1]. Thus, the possible crack paths are determined a priori and the study of crack propagation in dams can be based on the development of interface models suitable for concrete material. Several interface models have been proposed in literature to study crack propagation in cementitious material and at bimaterial interfaces. Among others, Carol et al. [2] developed a mixed-mode interface based on the definition of a hyperbolic activation function that provides a good fitting with experimental data. Cĕrvenka et al. [3] proposed a nonlinear elastic model based on an extension of Hilleborg’s model to study crack propagation at the rock/concrete interface. They considered the evolution of the failure function based on a softening parameter which is the norm of the inelastic displacement vector. Ruiz et al. [4] developed a three dimensional cohesive interface model that takes into account tension-shear damage and mixed-mode fracture. They developed a finite element calculations considering crack
29
Embed
A cohesive damage–friction interface model accounting for water pressure on crack propagation
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
A cohesive damage-friction interface model
accounting for water pressure on crack propagation G. Alfano
*, S. Marfia
**, E. Sacco
**
*Dipartimento di Scienza delle Costruzioni
Università di Napoli, Via Claudio, 21 - 80125, Napoli, Italy
Abstract. Aim of the work is the study of crack propagation in concrete constructions, such as dams,
taking into account the water pressure effects and the damage-friction evolution with unilateral
contact. With this aim, an interface model, based on the cohesive fracture, is developed. In particular, a damage-friction model based on a new multiscale approach is proposed for the interface, which is
able to simulate crack propagation in mode I, mode II, and mixed mode. The criteria for crack
initiation and propagation and the closure and reopening of the crack are considered for the interface.
The friction and the dilatancy are characterized by a decreasing rate of friction and dilatancy angles, respectively, with the possibility of the dilatancy recovering during cyclic loading. The water pressure
effect is taken into account by considering a further static pressure acting on the real crack and in the
process zone. The value of this pressure is set as a function of the crack opening displacement and it varies from a prescribed initial pressure to the external water pressure. A numerical procedure is
developed in order to integrate the evolutive laws governing the interface model. The interface model
is implemented in finite element codes; some numerical simulations and a comparison with experimental data are performed in order to validate the proposed model and to investigate the
influence of the water pressure effects on the crack propagation. Finally, the mechanical response of a
Initially, the damage evolution is evaluated solving equations (26); thus, the updated damage
parameter D is obtained. Then, the friction evolutionary problem is approached developing a
predictor-corrector algorithm.
The trial step is computed assuming frozen the evolution of the inelastic relative
displacement:
,1 1
,
, ,
, ,
0 0
, ,
ˆ
1ˆ 1
0
1
ˆ
di TR
n
di TR di
n
d TR di TR
s sd TR d TR
L
qTR
f i i
TR TR d TR d TR
p p e p
e
s s
σ H s s
σ σ (28)
If 0TR , the step is elastic and the trial represents the solution of the finite time step, i.e. ,di di TRs s , ,
ˆ ˆd d TR
, ,d d TR , ,d d TR , TR . Otherwise a return map algorithm is
developed in order to solve the corrector phase.
The equations to be solved are the following:
1 1
0 0
1 1 1 1 1
,
2 2 2,
1 1,
2 2,
11 1, 2 2
2
1ˆ
1ˆ
1
ˆ0
1
di
n
n
s sd dL
di did
d d TR di di
n
qdi di
n o
ddi di
n d
q
f i i
d d
di dqdi di di di
n odi d
p p e p
h s s H s s
H s s
u u e
u u
e
sv v h u u e
s
2 0
n
di
di di di
v
s u v (29)
It can be proved that when 2 2 0di diu u , i.e. when the second term of the right-hand side of
formula (29)8 does not vanish, it results 1 2 1, 2,/ /di di di di
n ns s s s .
In order to simplify the notation, in the following the Heaviside functions involved in the
model formulation are denoted as:
1 2
1 1 2 2
di di dih h s s h h u u (30)
From equations (29)3, (29)5, (29)9 and (29)10, it results:
G. Alfano, S. Marfia and E. Sacco
12
,2 1 d d TR
d
H
(31)
whose absolute value of the two members allows to obtain:
,
2
d d TR H (32)
Combining equations (31) and (32), it results:
,
,
d TRd
d d TR
(33)
Taking into account the definition (30) and formulas (29)2, (29)6 and (32), the consistency
relation (29)7 is written in the following residual form:
1 ,
1 1 1 21 1 0nq di d TR
f i iR e h H s s H
(34)
In order to evaluate at each time step, the nonlinear residual equation (34) is solved
developing a Newton-Raphson technique. In the following, the superscript k denotes the
quantities at the k th iteration and the superscript 1k those at the current iteration. The
symbol indicates the quantities variation between two consecutive iterations, i.e.
1k k .
By linearization, the equation (34) becomes:
10 k k
k
dRR R
d
(35)
with
1,
1
,1,2, 2, ,
1 1, 2,2,
1 1
1
kn
kn
qk k
f i i
di d TRq ndi k k k d TR k
n o di d TRn
R e h H
ss s h e h H
s
(36)
1,
1
,1,2, 2,
2,2,
1,
1
,1,2, 2,
1 1, ,2,
1 1
1 1
1
1
kn
kn
kn
kn
q k
f i i
k
di d TRq nk k k
o di d TRn
q k
f i
di d TRq ndi k k k
n o di d TRn
dRe h H
d
sh e h H
s
e h H
ss s h e h
s
(37)
assuming the derivatives of 1,kh and 2,kh with respect to vanishing.
G. Alfano, S. Marfia and E. Sacco
13
The plastic multiplier increment is obtained solving equation (35). Then, the iterative
procedure goes on until a convergence test is satisfied, i.e. when the values of the residual 1kR is less then a prefixed tolerance. Once is evaluated, the normal and tangential
stresses can be computed from formulas (29).
4 CONSISTENT TANGENT OPERATOR
The consistent linearization of the finite-step relationships is obtained as the derivative σ s
of the micro-macro stress equation (25):
1d
d
t
D DD D
σ σK K Ks σ
s s s s (38)
The derivative D s is equal to zero when there is no damage evolution, i.e. nD D ;
otherwise, in the case of damage growth, i.e. nD D , it is evaluated differentiating equation
(26)3:
22
11
1
D D D
s s s s s (39)
Taking into account equations (26)1 and (26)2, it results:
1
2
1
2
2
2
1
1
o
o
s
s
s
s
s (40)
1 21
2 2 14
1
2 h s sss
s
s s
(41)
For what concerns the term /d σ s in equation (38), two cases can occur in the finite time
step:
no-evolution of the inelastic relative displacement, i.e. 0 ,
variation of the inelastic relative displacement, i.e. 0 .
In the first case, the derivative of equation (27)1 is evaluated considering di di
ns s :
1 1
0 1 1
1 0
0 0
dids s di
Lp p e h s s
σH
s (42)
In the case of inelastic evolution, differentiating equation (27)1 gives:
1 2
1 2
d d
d
d d
s s
s s
σ
s (43)
with
G. Alfano, S. Marfia and E. Sacco
14
1 1
1 1
1 11 0 1 1
1 1
1 11 0 1 1
2 2
22
1 1
22
2 2
1 1
1
1
di
di
dids s di
L
dids s di
L
did
did
sh H p p e h s s
s s
sh H p p e h s s
s s
sH
s s
sH
s s
(44)
The derivatives of the inelastic relative displacements are computed differentiating equations
(29)4, (29)5 and (29)8:
1 1 1 1 1 1
1 1 1 2 2 2
1 1
1 1 2 2
1 1
1 1 2 2
, ,
2 2 2 2
, ,1 1 1 2 2 2
1 2
di di di di di di
di di
di di
di di d TR di di d TR
d TR d TR
s u v s u v
s s s s s s
u uA A
s s s s
v vB B
s s s s
s u s u
s s s s s s
C Ds s
(45)
The parameters in equation (45) are the following:
,1,2
,2,
,1
1 2,
1 2 2
1 1 1, 2
,1,2
,2,
1 1
1
1
1
1
n
di d TRn
di d TRn
d TR
d TR
q di
f i n
di d TRn
di d TRn
sA B h A
s
C H h E D H E
Eh H h A F e s s h F H
sF h
s
(46)
5 NUMERICAL RESULTS
The interface model described in Sections 2-4 has been implemented in the finite-element
codes LUSAS [27] and FEAP [28] as a new constitutive law for 2D interface elements.
Several types of numerical studies are developed in the following: first, an investigation
devoted to illustrate the main features of the proposed interface model is presented, then, a
comparison with experimental results, available in literature, is performed; finally, the
mechanical response of a concrete dam is studied.
G. Alfano, S. Marfia and E. Sacco
15
Behavior of a typical point of the interface
As a first step towards the assessment of the reliability and the efficiency of the model,
numerical simulations have been conducted to derive the mechanical response in a typical
point of the interface considering different loading-unloading histories. The simulations are
performed at level point without studying any structural problem requiring a mesh
discretization.
The following interface parameters are considered:
1 o1 1 0
2 o2 2
0.35 N/mm 0.3 MPa 0.006 mm 0
0.09 N/mm 0.7 MPa 0.010 mm 1
tan 30 0
c o
c o
o
f i
G s p
G s
(47)
From the above input data, the stiffness parameters 1K and
2K are evaluated. Furthermore, it
has been set 1 1H K and
2 2H K .
Initially, in order to investigate on the effect of the uplift water pressure and on its
dependency on the interface opening relative displacement, a mode I tensile test is performed.
The hydrostatic pressure is set equal to 0.8 MPaLp , while different values of the parameter
are considered in the computations.
In figure 4, the interface normal stress versus the crack opening is plotted for -10, 0.1,1,10,100 mm . From the numerical results different softening behaviors of the
interface can be noted. For all the analyzed values of the uplift parameter , the maximum
normal stress and the final damage relative displacement are the same of those corresponding
to 0 , i.e. 0.0Lp . When is set equal to zero a linear softening is recovered, while,
when the uplift pressure is considered, the softening curve depends on the value of the
parameter : a steeper softening is obtained increasing the value of .
Figure 4: Interface normal stress for mode I opening.
G. Alfano, S. Marfia and E. Sacco
16
Then, the shear behavior of the proposed interface model is studied. A normal stress
2 MPa is assigned and taken constant, while the tangential relative displacement is
increased from 0 to 1.2 mm.
In figure 5, the shear stress versus the tangential relative displacement 2s is plotted for
three different values of the dilatancy parameter, 0 0, 0.1, 0.5 , with -10.8 mm . It can be
noted that the interface behavior is influenced by the value of the dilatancy parameter. In fact,
increasing the value of 0 an increment of the maximum shear stress is attained, while a
steeper softening branch occurs. This special interface shear response can be interpreted
looking at the micromechanics of the REA. As described in Section 2, according to the
proposed model, the shear stress is the sum of two terms: a first one due to the elasticity of the
undamaged part of the REA and a second one due to the friction, proportional to the local
normal stress occurring on the damaged part. Because of the ditatancy, a positive normal
relative displacement occurs in the REA, which induces a decrement of the absolute value of
the normal stress acting on the undamaged part of the REA; as the overall equilibrium of the
REA has to be satisfied, the local normal stress in the cracked part of the REA increases in
absolute value, leading to an increment of the friction stress.
Figure 5: Shear stress versus tangential displacement for different dilatancy parameters 0.
In figure 6, the overall stress at the interface and the local stresses in the damaged and
undamaged part of the REA weighted for their corresponding fraction of area, i.e. 1 uD
and dD , are plotted versus the tangential relative displacement 2s , for 0 0 and 0 0.5 .
It can be noted that the normal stress in the damaged part of the REA is higher when the
dilatancy effect is considered.
G. Alfano, S. Marfia and E. Sacco
17
Figure 6: Normal stresses at in interface: overall stress, stress on the undamaged part
and stress on the damaged part.
The steeper softening behavior in presence of dilatancy is due to the development of positive
normal relative displacement in the interface which, according to formulas (7) and (8), leads
to an increase in the rate of damage parameter evolution.
Then, a numerical study of the interface subjected to cyclic shear laoding-unloading history is
performed, considering the dilatancy parameter 0 0.5 and assuming the total normal stress
2 MPa constant during the cycles. In figure 7 the shear stress is plotted versus the
relative displacement 2s . Results are obtained through a displacement control. With reference
to figure 7 and figure 8, the interface response consists in the following steps:
OA, the tangential relative displacement 2s is increased from 0 to 1.0 mm; the shear
stress reaches its maximum and, then, a softening phase occurs until the complete
damage of the interface is reached; during this first step a dilatancy effect is obtained;
AB, the tangential displacement is decreased to 0; the shear stress attains a negative
constant value; the accumulated dilatancy normal displacement is recovered;
BC, the tangential displacement is decreased to -1.0 mm; the shear stress has a negative
constant value; the dilatancy normal displacement is obtained, characterized by a
reduced value respect to the first step;
CD, the tangential displacement is increased to 0; the shear stress attains a positive
constant value; the accumulated dilatancy normal displacement is again completely
recovered.
The shear behavior of the interface subjected to several loading-unloading cycles has been
studied according to the following horizontal displacement history:
2
[s] 0 1 2 3 4 5 6 7 8 9
[mm] 0 0.125 0.25 0.375 0.5 0.625 0.75 0.875 1 1
t
s (48)
G. Alfano, S. Marfia and E. Sacco
18
Figure 7: Shear behavior for one complete loading cycle: shear stress vs horizontal
displacement.
Figure 8: Shear behavior for one complete loading cycle: dilatancy vs slip.
In figure 9, the shear stress versus the horizontal displacement is represented; a description of
the interface behavior in the time interval 0,2t is reported in detail:
OA, during the first loading phase, a partial damage occurs at the interface, with
frictional sliding at the damage part of the REA;
G. Alfano, S. Marfia and E. Sacco
19
AB, an unloading phase follows, characterized by a linear response with the initial
stiffness and for 0 a residual inelastic displacement is found; in this phase there is no
damage growth;
BC, the slope of the 2s curve changes because of the development of the inelastic
negative slip, without damage evolution;
CD, at point C, the absolute value of the horizontal displacement is equal to the one in
point A, thus the damage increases again beyond point C; the response of the interface is
characterized by inelastic slip and damage growth.
Figure 9: Shear behavior for several loading-unloading cycles: shear stress vs horizontal
displacement.
Similar considerations can be made for the subsequent cycles. The final reloading is
characterized by a flat 2s curve due only to the friction effect, as the interface is
completely damaged.
In figure 10, the dilatancy versus the slip is plotted, reporting also the time value in order to
better understand the global interface response. It can be noted that at each positive or
negative loading-unloading the dilatancy is recovered and the slope of the 1 2
di dis s curve
decreases during the cycles until the curve becomes flat, i.e. dilatancy is completely sutured.
From the numerical results concerning the mechanical response, the proposed interface model
appears able to reproduce the main features of the behavior of material discontinuities in
concrete constructions also in presence of uplift water pressure; in particular, the model can
be used to simulate the response of joint-interfaces resulting by the typical step-wise
construction process or by the interaction between rock foundation and the construction itself.
G. Alfano, S. Marfia and E. Sacco
20
Figure 10: Shear behavior for several loading-unloading cycles: dilatancy vs slip.
Comparison with experimental results
A comparison between the numerical results derived adopting the proposed interface model
with the experimental ones, obtained by Lee et al. [31], is performed. In particular, the
behavior of the interface between rock blocks of Hwangdeung granite is studied. In figure 11,
the scheme of the experiment developed in Ref. [31] is reported. Initially a normal load is
applied on the top of the upper block, then a shear loading cyclic history is performed keeping
constant the normal load. To model the experimental test a 2D analysis is carried out adopting
one element discretization for the upper block and a rigid support simulating the lower one.
The joint is modeled by the proposed interface finite element as schematically reported in
figure 11.
Figure 11: Geometry and loading of the joint.
The following interface parameters are considered:
G. Alfano, S. Marfia and E. Sacco
21
1 o1 1 0
2 o2 2
-1 -1
0
1.255 N/mm 1.97 MPa 1.25 mm 0
1.255 N/mm 1.97 MPa 1.25 mm 1
0.9 1.9 0.5 mm 0.162 0.002 mm
c o
c o
f i
G s p
G s
(49)
The joint is subjected to the loading history reported below:
1
2
[s] 0 1 2 3 4 5 6
[MPa] 0 1 1 1 1 1 1
[mm] 0 0 15 15 15 15 0
t
s
(50)
Each shear loading cycle is divided into four stages: forward advance (stage I), forward return