Pure Mechanics Crack Model For Shear Stress Transfer in Cracked Reinforced Concrete
Pure Mechanics Crack Model For Shear Stress Transfer in Cracked Reinforced Concrete
An analytical model is presented that is capable of predicting or the main shear crack can either be computed as 11.57 MPa
assessing the response of cracks in reinforced concrete elements (1678 psi) according to the Contact Density Model,10 or as
subjected to monotonic, cyclic, or reversed cyclic in-plane shear 2.76 MPa (387 psi) employing the model of Walraven and
and normal stresses. The model is formulated in terms of global Reinhardt.9
and local equilibrium, compatibility, and stress-strain relation-
Interface shear-stress transfer across a crack is a relatively
ships. The small number of required empirical parameters are
complex process depending on parameters such as loading
determined from the test results of 14 pre-cracked reinforced
concrete panels under uniform shear and uniaxial tensile stresses. history, crack width, crack slip, compressive stress acting
The proposed model is shown to capture the important aspects of across the crack, concrete strength, aggregate size, crack
the response of cracks subjected to complex loading conditions, roughness, and tensile stress in the reinforcement crossing
resulting in improved simulations of crack behavior. In addition, it the crack. Given the many mechanical and geometrical
is shown that the model can be employed to conduct crack-based parameters involved, crack models normally rely on a signif-
structural assessments that estimate the stress state of a struc- icant number of empirical equations with constants obtained
ture based on crack displacement information collected as part of by fitting to available results.
site-monitoring activities. This paper presents a simple rational crack behavior
model for cracks subjected to cyclic and/or reversed cyclic
Keywords: aggregate interlock modeling; crack behavior; cyclic loading;
nonlinear response; reinforced concrete; shear strength; structural shear and axial loads. The model, referred to as the “Pure
assessment. Mechanics Crack Model” (PMCM), is formulated in terms
of equilibrium, compatibility and constitutive consider-
INTRODUCTION ations. The only empirical equation employed is a constitu-
After cracking, reinforced concrete (RC) members carry tive equation, successfully employed over the course of the
a significant portion of the applied shear loads across cracks past three decades.2,14 A significant feature of the proposed
through aggregate interlock action. For example, it has been model is that it can be used to perform crack-based struc-
shown that RC beams not containing shear reinforcement tural assessments (that is, estimating the stress state of
carry up to 76% of the total shear load across existing shear the structure from crack opening, sliding, and orientation)
cracks, with the remaining portion of the shear being resisted without needing empirical parameters and requiring only
by the uncracked portion of the cross section.1 Because of a small number of relatively easily obtainable properties
this, to better model the behavior of RC structures subjected of the structure.
to shear loads, it is necessary to accurately understand and
simulate the response of the cracks.2,3 RESEARCH SIGNIFICANCE
The shear-transfer behavior of cracked concrete inter- The model presented in this paper can be used to predict
faces has been a subject of interest to many researchers over the behavior of cracks under complex loading conditions or
the past few decades.4-15 Based on available experimental to assess the stress state and the reserve capacity of the crack
results, a number of these investigators proposed crack based on measurable crack information.
behavior models to predict the shear response of cracked
interfaces subjected to monotonic and cyclic loading condi- OVERVIEW AND CONCEPTUAL MODEL
tions.9-11 However, all the proposed models show limitations A recent experimental study conducted at the University of
depending on the context they have been used in and none Toronto involved 14 precracked RC panel elements subjected
have been shown to perform consistently better than the to cyclic or reversed cyclic shear and axial tension.4,15 The
others.4 In this regard, it is of concern that the shear stress elements were tested using the Panel Element Tester and
transfer predicted by the different models for particular crack were subjected to uniform tensile stress fx and uniform shear
widths and crack slips can vary by factors of up to five.4 stress v, as shown in Fig. 2. Each element was reinforced in
To illustrate the variability of existing crack models, the two orthogonal directions, x and y, and was crossed by a
consider the crack outlined in Fig. 1, which formed in concrete
with a strength of 38 MPa (5510 psi) and a maximum aggre- ACI Structural Journal, V. 114, No. 2, March-April 2017.
gate size of 14 mm (0.55 in.). If it is assumed that the crack MS No. S-2016-157, doi: 10.14359/51689460, received April 18, 2016, and
reviewed under Institute publication policies. Copyright © 2017, American Concrete
width is 1.25 mm (0.05 in.) and the crack slip is 2.5 Institute. All rights reserved, including the making of copies unless permission is
obtained from the copyright proprietors. Pertinent discussion including author’s
mm (0.10 in.), the amount of shear stress transferred across closure, if any, will be published ten months from this journal’s date if the discussion
is received within four months of the paper’s print publication.
w
ε sx = − ε sh (1)
Leff
(L )
2
eff + w − ε sh Leff + s 2 − Leff
ε sx = (2)
Leff
Note that the only difference between Eq. (1) and (2) is
that the latter accounts for second-order effects. For this
reason, Eq. (2) may be more appropriate for analysis of “ulti-
mate strength” cases. The term εsh represents the strain due
to concrete shrinkage and can be taken as 0.3 × 10–3 in the
Fig. 1—Shear-critical reinforced concrete beam. (Adapted absence of more accurate measurements.
from Sherwood et al.1)
w
tan (α ) = (3)
s
In a similar approach to that of Walraven and Reinhardt,9
it is assumed that at the time a crack opens, the crack surface
consists of a series of spherical aggregate particles. To
further simplify the problem, it is assumed that the crack
plane intersects the most critical of these particles at their
middepth. The crack slip can then be computed from geom-
etry alone as
a g − a g 2 − 4 ⋅ w2
s= (4)
2
where ag represents the maximum aggregate size.
Note that Eq. (4) is only valid at the initial stages of the
life of a crack. As local crushing phenomena begin to occur,
the crack slip can no longer be computed just from geom-
etry, as will be discussed later.
Fig. 4—Free-body diagram of aggregate at contact point
Equilibrium conditions (note that steel inclination is intentionally exaggerated).
As already mentioned, vertical equilibrium of the free body
shown in Fig. 2(b) requires that the shear stress v applied to the right face of the element and the crack plane, these two
the right face of the panel must be exactly balanced by shear global equilibrium equations for the element can be formu-
stresses on the crack interface, vci. Similarly, horizontal lated as
equilibrium requires that the tensile force in the longitudinal
reinforcement at the crack must be in equilibrium with the fci = ρx ∙ fsx,cr – n ∙ v (5)
applied tensile force fx and the compressive stress fci that
develops on the crack plane. Assuming a unit area for both
vci = v (6)
Unloading phase
The crack system subjected to a given combination of
shear and tensile stresses displaces until an equilibrium posi-
tion has been reached. The problem at hand is now to deter-
mine the crack response under decreasing applied loads. The
local stress state at a contact point evolves from the scenario
depicted in Fig. 4(a) to that outlined in Fig. 4(b) and 4(c).
Over the course of the unloading phase, two possible
scenarios can occur:
1. The displacement configuration remains essentially
unchanged; or
2. Crack width and crack slip tend to decrease somewhat
proportionally to the reduction in applied load magnitude.
Crack width and crack slip are insensitive to stress reduc-
tions if
vci
≥ tan α (9)
f ci
If vci/fci < tanα, crack displacements can be either affected
or unaffected by a reduction in applied loads. More specifi-
cally, the crack motion tends to reverse if
f ci ⋅ ( µ ⋅ cos α − sin α )
vact ,1 = (12)
(cos α + µ ⋅ sin α )
The shear stress that produces loss of contact between the
opposing faces of the crack can be computed as
ρsx ⋅ f sx ⋅ cos α
vact , 2 = (13)
(sin α + n cos α )
The activation stress is then taken as the value whose
absolute value is the smaller of that given by Eq. (12)
and (13).
Constitutive relations
Constitutive relationships are required to link stresses
and strains for both the steel reinforcement and the aggre-
gate-paste matrix at the crack location. Note that the
compatibility and equilibrium equations introduced so far
are not sufficient in number to make the problem at hand
statically determinant.
In regards to the steel reinforcement crossing the crack
plane, the stress fsx,cr is assumed to depend only on the axial
strain in the reinforcement εsx,cr. In addition, it is assumed
that the dowel action of the reinforcement can be neglected.
Consequently, axial stress and axial strain can be related
through the usual bilinear uniaxial stress-strain relationships
such as