RC Elements Under Cyclic Loading State of The Art Report
RC Elements Under Cyclic Loading State of The Art Report
RC ELEMENTS
UNDER CYCLIC
LOADING
STATE OF THE ART REPORT
I Thomas Telford
Published by Thomas Telford Publications, Thomas Telford Services Ltd, 1 Heron
Quay, London E14 4JD
First published 1991 as CEB Bulletin d'Information No. 210 Behaviour and
analysis of reinforced concrete structures under alternate actions inducing
inelastic response - volume 1: general models
Thomas Telford edition published 1996
A catalogue record for this book is available from the British Library
Classification
Availability: Unrestricted
Content: Subject area review
Status: Committee guided
Users: Designers, civil engineers
ISBN: 978-0-7277-3548-5
©CEB, 1996
All rights, including translation reserved. Except for fair copying, no part of this
publication may be reproduced, stored in a retrieval system of transmitted in any
form or by any means, electronic, mechanical, photocopying or otherwise, without
the prior written permission of the Publisher, Books, Publications Division, Thomas
Telford Services Ltd, Thomas Telford House, 1 Heron Quay, London El4 4JD.
This book is published on the understanding that the author is solely responsible for
the statements made and opinions expressed in it and that its publication does not
necessarily imply that such statements and/or opinions are or reflect the views or
opinions of the publishers.
Preface
During the second half of the 1980s, the CEB, with its new future Model
Code already under development, felt that the time had arrived to adjourn
and make known its ideas on the broad and rapidly evolving field of non-
linear analysis of structures. In fact, while the number of applications of
non-linear analyses for important and unconventional structures was
increasing fast, and recourse to such analyses was explicitly allowed by
the Model Code, no guidance was available to the designer for the selection
of the models more suited to each specific case. On the other hand, the
existence and the accessibility of a number of sophisticated finite element
codes for the non-linear analysis of reinforced concrete structures, was a
cause more for concern than for relief, because of their lack of validation
over a sufficient spectrum of different applications.
The terms of reference given by the CEB to the new group to be formed
stated that it should concern itself with those models and methods of
analysis adequate for dealing with general load and stress histories. This
decision was obviously of great consequence, in that it involved the
exclusion of a number of classical and well proven models and methods, but
it was certainly a necessary one: the assumption of a proportional increase
of all loads acting on a structure is not only unrealistic but is also not
conservative in assessing the likelihood of attaining a given limit-state.
Furthermore, since a fairly large number of real loads fluctuate with time,
and some of them may also reverse their action, it was finally decided to
characterize the scope of the group's work with the widest of its attributions,
by titling it 'Behaviour and analysis of reinforced concrete structures under
alternate actions inducing inelastic response'.
From the early discussions on the plan of activity, the members concurred
that the first step had to be the preparation of a document of fundamental
nature. The document would contain a review of the pertinent existing
knowledge in the area of constitutive modelling of concrete, steel, bond and
of their interaction. It should also discuss the problems encountered in
assembling the various elements with the purpose of constructing the model
of an element made of reinforced concrete. In selecting the material for
inclusion, the main criterion would be its ability to reproduce, adequately if
crudely, non-proportional stress-strain paths, with reversed loading.
Whether physically or empirically based, very simple or sophisticated,
long-established or brand new, the models would be presented in framework
as rational as possible, certainly accompanied by comments on their
advantages and limitations, but without attempts to rank them on the basis
of their relative merits.
This programme of work, with the addition of a final part dealing with the
solution strategies for the numerical integration of non-linear equations, has
materialized in the present document. Admittedly distant from the
immediate needs of professional practice, it aims at being a basic, com-
prehensive and relatively easy-to-read document, hopefully useful as a
reference for non-specialists wishing a more conscious approach to the use
of the available sophisticated computational tools.
The material in this book has been diffused initially as CEB Bulletin 210.
It is now being published integrally in an editorially revised form.
Paolo E. Pinto
Rome, March 1996
Acknowledgements
1. Concrete in compression 1
1.1. Introduction, 2
1.2. Experimental behaviour of concrete under general
multiaxial stress histories, 5
1.3. Models derived from the theory of elasticity, 9
1.4. Models based on the theory of plasticity, 16
1.5. Models based on the mechanism of microcracking
or elastic damage, 26
1.6. Models based on the association of plasticity and
microcracking mechanics, 31
2. Concrete in tension 42
2.1. Introduction, 42
2.2. The tensile behaviour of concrete, 42
2.3. The fictitious crack model, 43
2.4. Models for the bulk behaviour, 43
2.5. Models for the stress-crack width relation, 43
2.6. Rheological material model for the cr-w relation, 47
2.7. Continuous function model, 54
3. Reinforcing steel 58
3.1. Introduction, 58
3.2. Giuffre, Menegotto and Pinto model, 60
3.3. Simple non-linear stress-strain model, 63
5. Interface behaviour 98
5.1. Introduction, 98
5.2. Aggregate interlock, 102
5.3. Dowel action, 122
5.4. Concluding remarks, 132
6. Finite element modelling of reinforced concrete 134
6.1. Introduction, 134
6.2. Finite element modelling of the reinforcement and
of its interaction with concrete, 135
6.3. Modelling of the concrete component, 136
6.4. Composite steel/concrete modelling, for two-dimensional
homogeneously reinforced elements, 161
References 179
1. Concrete in compression
In the stress space, the locus of all the stress combinations for which a
proportionally loaded concrete specimen reaches its maximum load-
carrying capacity is usually called the 'failure surface' or, a better term,
the 'ultimate strength surface'.
It is a convenient simplification to assume that an ultimate strength
surface constructed for proportional loading may still be used for defining
failure conditions for more general load cases, provided the individual
components of the stress tensor do not deviate too much from monotonic
path.
In dealing with the modelling of the behaviour of concrete over its
complete range of response, the knowledge of the ultimate strength surface
is important, since it allows identification of regions of stress states beyond
which complete failure is reached according to different mechanisms.
The different failure modes are shown schematically in the following
figure (adapted from Pramono and Willam1), with the corresponding
representative stress points in a ultimate strength surface plotted in the
Rendulic plane. For example, experiments indicate the existence of a
transition point (TP) which separates brittle softening behaviour from a
ductile regime, this latter occurring when the lateral confinement ratio
increases. The hatched area in Fig. 1 corresponds to softening post-peak
behaviour which, however, exhibits markedly different characteristics
depending on whether the specimen is under one principal compressive or
tensile stress.
For the case of predominant tensile stresses, in which failure occurs along
a well-defined direction and in the form of a localized crack, fracture
mechanics approaches have long since been introduced to derive stress-
displacement relationships of the cracked zone. This approach is still the
most commonly used although it has to be complemented with some criteria
for 'smearing' the discrete phenomenon over an equivalent continuum. All
other cases, from the compression failure modes to the mixed ones, are
treated in the classical context of the constitutive laws for continuous media.
This abrupt change from a discrete to a continuous type of modelling is
obviously not satisfactory from a physical standpoint.
There are in fact some recent schools of thought2 supporting the view that
Macroscropic models
Models derived from Models based on Models based on
theory of elasticity theory of plasticity progressive damage
Equivalent Invariant- Classical Bounding surface Fracturing Continuum
uniaxial based damage variables damage
Elwi- Stankowski- Vermeer- Fardis- Dougill17 Krajcinovic18
Murray8 Gerstle10 De Borst12 Chen14
Lemaitre19
Buyukozturk- Shafer- Han- Chen- Resende-
Shareef9 Ottosen" Chen13 Buy ukozturk15 Martin20
Yang et al. Krajcinovic-
Fonseka21
Mazars22
Coupled models
Micromechanics models
Microplanes Ba2ant-Prat
of the two first stress and strain invariants (octahedral values). The
underlying assumption of isotropy is not substantially faulted, at a
macroscopic level, up to stress states close to failure.
• The direct relationships between volumetric and deviatoric stresses
and strains are governed by the evolution of two independent moduli
K and G.
• A coupling effect between octahedral shear stress and volumetric
strain is systematically observed, and can be described by
introducing a coupling modulus H.
According to the findings above, three stress and strain dependent moduli
appear to be sufficient for describing the behaviour of concrete under
monotonically increasing loads. Within the same restricted range of
response, a discussion on the features of concrete behaviour which are
universally agreed, together with some guidance on the possibility of
analytical modelling, are given in the Synthesis Report prepared by CEB in
1983.17
Concurrently with the described efforts aiming at screening and
systematizing the available experimental knowledge, computer codes for
the analysis of structures subject to general load types were being developed
at a fast pace, thus creating the need for constitutive laws able to follow
non-monotonic, non-proportional loadings and unloadings. This necessity
was readily recognized by the research field, which was prompt to respond
with a spectrum of advanced analytical models overcoming the previous
limitations, sometimes even beyond the boundaries of previously
established knowledge. These models have been the subject of a review,
last updated in the early 1980s.6
The broadening of the experimental knowledge to cover at least the most
basic non-proportional loading paths is a very recent process, and much
work remains to be done. Only a few of the laboratories active in concrete
research have made such types of tests, the largest amount of data coming
from the University of Colorado,25 Imperial College in London,26 the
Bundesanstalt fur Material Priifung in Berlin,28 and the Air Force Weapons
Laboratory in the New Mexico State University.29 The data obtained by
these institutions constitute the basis for validating all the models that will
be presented below.
= O3 o2 =
(a) (b)
obtained, as well as the coupling between the deviatoric stress TQ and the
volumetric strain eo-
From the first type of test one can observe that for continuously
increasing a0 the curve <r0 — eo starts with an initial maximum value of the
tangent bulk modulus K = Ko; the slope of the curve subsequently decreases
in an almost linear fashion with increasing cr0, up to a certain threshold
beyond which K remains approximately constant with a residual value of
about 15-20% of Ko. Unloadings and reloadings up to the same stress point
show an almost linear elastic behaviour with a modulus close to Ko. No
significant dependency of K from the amount of deviatoric stress is
observed.
The coupling between a0 and 70 is shown in Fig. 3. For a given increase
of the hydrostatic stress <T0, the amount of deviatoric strain 70 taking place is
an increasing function of TQ. This deformation appears to be irreversible
upon the unloading of CTQ, and it is also insensitive to a subsequent reloading
of cr0, at least up to its previously attained maximum value.
The direct relationship between TQ and 70 as obtained from the test type of
Fig. 2(b) is shown in Fig. 4. The three curves refer to different <TO values:
their common feature is a softening behaviour. The increase of cr0 also
increases the maximum shear deformation attainable prior to failure. As for
the volumetric quantities, in this case also unloadings and reloadings occur
almost elastically and with a modulus close to the initial one, at least for
unloading points sufficiently far from the ultimate stress.
Finally, in Fig. 5 the curves To — eo for different values of the hydrostatic
pressure are represented. There are two main notable features: the coupling
-225
r = 0-69 r = ;>07
O
-1-50
-0-75
I \
CU
i
10
/ •
= 2-28
— -
/
.
-
Jo = 3 0 8
o0 = 5-36
1
0-5
ICL
12 16
1-5
5-36
10
7/
o = 2 28
0-5
-8
CONCRETE IN COMPRESSION
14 - 1
12 -
a0 = 8 ksi 10
6- 8
/ jf 3/
V
6 /2
4 / 1
1/ fr
J
2
cu
0 ' . 1 . 1 . 1 , 1 . 1 . 1 . 1 1 1 . 1 . 1 . 1 . 1 . 1 . 1 . 1
10 12 14
Fig. 6. Non-proportional loading tests: (a) stress paths; (b) curves <7j — e* for each path
The state of stress is first increased from zero to a hydrostatic value of o-0
= 560 kg/cm2, which is henceforth held constant. It is noted that the initial
state of stress is such as to induce inelastic response. Path 1 consists of
increasing az, while ax and ay are decreased so as to keep 00 constant.
It can be seen in Fig. 6(b) that the curve az — ez proceeds inelastically
(softening behaviour), while the two unloadings are approximately elastic.
When az is subsequently decreased to its original value, all the three cr, — e,
branches are linear elastic. In path 2 the stresses az and ax exchange their
role, but the increase of ax in this case is such as to remain below the value
previously attained by az.
It can be seen that the curve ax — ex proceeds inelastically, while the two
unloadings of o~z and ay are approximately elastic. This observation is
already in contrast with the |To|max criterion, according to which path 2
should have been entirely within the elastic range.
The further paths shown in Fig. 6(a) and 6(b) support the statement that
an adequate criterion for distinguishing loading can be based on the prin-
cipal stresses: loading occurs (and inelastic behaviour with it) whenever one
(or more) principal stress exceeds its previously attained maximum value.
In concluding this brief overview of the available experimental
information on the behaviour of concrete under multiaxial stress states, it
has to be stressed that the tests carried out up to now have not explored
some of the phenomena whose existence and importance are well
recognized in the uniaxial case, and which must be presumed to play a
no less vital role in the three-dimensional case. These phenomena include,
among other aspects, the post-peak behaviour and the stiffness degradation
due to repeated loadings and unloadings. The more advanced mathematical
models that will be presented in the following paragraphs are in fact able to
describe these phenomena: their calibration, however, cannot be based at
present other than on mono- or bi-dimensional results, so that they must be
cautiously looked upon as reasonable extrapolations from known data.
1.3. Models derived The distinctive aspect of this voluminous and still widely used category of
from the theory of models is that the incremental relationship between strains and stresses
elasticity maintains the same form as for the elastic solids, i.e.
day — Djju deu (1)
RC ELEMENTS UNDER CYCLIC LOADING
where D,/W is the usual stiffness tensor, but its elastic parameters are now
made functions of the current stress state (tangent values).
Since the physical behaviour of concrete subjected to non-proportional
stress paths in the range approaching failure is substantially different from
that of an elastic body, the formal structure of equation (1) is a priori
inadequate for the general case. Despite this intrinsic limitation, the
assiduous work completed over the past decade to improve the realism and
enlarge the scope of these models, coupled with the firm experimental basis
from which they have emerged and against which they have been carefully
calibrated, have earned these models a stable place within the class of the
so-called 'empirical' models.
The discriminating feature for these models to be considered in the
present context is the inclusion of a criterion for distinguishing loading from
unloading and reloading, with different behaviours for each of these cases:
otherwise, the observed strong non-holonomy of the stress-strain relation-
ship could not be reproduced. Other aspects, such as, for example, an
accurate description of the coupling between deviatoric and volumetric
components of stresses and strains, although certainly important in some
cases, are not deemed equally essential.
The number of parameters entering the matrix D,//W is potentially very
large so that some initial hypotheses have to be introduced to get an
acceptable balance between accuracy, ease of correlation with experimental
results, and simplicity. Two main ideas have been pursued so far.
(a) To consider inelastic concrete as an orthotropic material, with
behaviour which, along each of the principal directions, is governed
by an 'equivalent' uniaxial stress-strain law: Elwi and Murray,8
Buyukozturk and Shareef.9 This model requires the definition of a
minimum of four state parameters.
(b) To describe the behaviour of concrete through a relationship between
the stress and strain invariants (Stankowski and Gerstle,10 Shafer and
Ottosen"). Depending on the number of invariants used, the two
mentioned approaches require the definition of four and eight
parameters, respectively.
A brief account of the two categories of model is given below, with closer
reference to the papers of Elwi and Murray and Stankowski and Gerstle,
respectively.
0 0 0
(A* 12^3 I + 0 0 0
symmetrical 0 0 0
(2)
<f>G\2 0 0
<j>G2i 0
G31
10
CONCRETE IN COMPRESSION
0-, = (3)
where
=
RE EQI
e (c cr,f e,c
?ff- 1 1
R =RE
(Re -
and the physical meaning of the parameters aic, e/c, cr,f, e,? is illustrated in
Fig. 7.
The current (tangent) value of the modulus Et is now obtained upon
differentiation of equation (3): E, = d<r,/de,u, which yields an expression of
the type
Ei — Ei(eiu/aic, e,c, an, e,f) (4)
and the current value of e m is the integral along the loading path
£. ^3 (5)
Fig. 7. Stress-equivalent
strain curve
RC ELEMENTS UNDER CYCLIC LOADING
For the Poisson's ratio v two choices have been made. One is to consider
v as a constant,9 which is a reasonable approximation only to 70-80% of
peak stress, beyond which v must increase in order to reproduce the well-
known dilatancy effect in the range approaching failure. For v = constant,
the matrix D of equation (2) becomes
h) v{\ + 0 0 0
(1 - v2)E2 i/(l + 0 0 0
symmetrical (1 — 0 0 0
<j)G\2 0 0 (6)
SGiT. 0
where / c ' and ecu are the uniaxial compressive strength and the uniaxial
ultimate compressive strain, respectively.
The remaining parameters in equation (3) are the coordinates (o/f, e,f of
a point sufficiently far away on the softening branch of the curve. It is
known that experimental determination of the softening branch, especially
at large strains, is test-dependent and often even questioned on theoretical
grounds. The empirical proposal made by Elwi and Murray consists of
taking e,f = 4e,c, aic = 4ai{.
12
CONCRETE IN COMPRESSION
\
[l/Y
where Aao, Ae0 are the hydrostatic stress-strain increments, and Aro, A70
the deviatoric stress-strain increments.
Notice that in assigning a generic stress increment: Acr,(/ = 1,2,3), the
corresponding Acr0 and ATO are uniquely defined, and so are the increments
Aeo, A70 from equation (7). These two latter quantities, however, are not in
general sufficient for the determination of the principal strain tensor
Ae,(/— 1,2,3), unless the stress-strain state is axisymmetric. For the
general case, it is assumed that the increments of the deviatoric strain
increments are parallel to the corresponding deviatoric stress increments,
although it is recognized that this is only a rough approximation of the real
behaviour.
1.3.2.1. Loading criterion and loading surface. Loading is defined as
the case in which one or more of the principal stresses increases beyond its
previously attained highest value. For stress states other than loading (i.e.
unloading and reloading) deformations are assumed to vary elastically
according to the tangent modulus at the origin for the 'direct' components of
strains (i.e. Aeo =/(A(7o) and A70 = / ( A r o ) ) , while the 'coupled'
components of strains (i.e. Aeo —/(Aro) and A70 =/(Aoo)) are assumed
to be completely unrecoverable.
Although the model being described is of a purely empirical origin and
could therefore be presented and used as such, it is possible to find certain
links with classical flow theories which are useful in providing a better
insight into its physical bases.
To this purpose, the concept of a loading surface well known in strain
hardening plasticity is recalled, and the model under consideration is
attributed a loading surface which changes its shape according to the
evolution of the stress state. More precisely, for stress states below 60-70%
of the peak stress, the loading surface is assumed as normal to the direction
of the increasing 07 (principal stress loading criterion), while for higher
stresses it tends to coincide with the failure surface. The mechanism of this
evolution will be clarified below.
1.3.2.2. Expressions of the moduli. The expressions of the bulk
modulus K and of the shear modulus G are given directly based on a large
number of experimental results.
K = Ko{\ - Ck <T0)
G = Go( 1 - ro/Tou) [GO = Gin( 1 - Cga0)} (8)
Equation (8) shows that both K and G display a softening behaviour with
increasing cr0 and To, respectively, and also that Go is a decreasing function
of the hydrostatic stress a. (rou is the current ultimate deviatoric strength for
the angle of similarity defined by the current stress state.)
As already mentioned, unloading and reloading proceed elastically with
moduli KQ and G,n, respectively. Expressions for the two coupling moduli
H= YY= —
Ae0 A70
have been derived for general load histories by Stankowski and Gerstle,
based on the results of ad hoc tests. If recourse is made to the loading
surface, however, the use of these empirical expressions is no longer
necessary.
Consider the stress state represented by point A in Fig. 8. Since there is
still distance between point A and the failure state, the loading surface is
13
RC ELEMENTS UNDER CYCLIC LOADING
Aeg = F(TO/TQU)
P
A7 0 V2
obeying the following conditions
14
CONCRETE IN COMPRESSION
• * * l fiTl no,
{ B2 B 3\ \
L
AwJ
with the following definitions for the quantities Ae and Aw
Ae = SjjAejj/3T0 Aw = skmsIC&Aeik/'iTl
which show Ae and Aw to bear similarities with the second and third
incremental strain invariants, respectively.
After splitting the two 'diagonal' moduli A\ and B2 into an elastic and a
'corrective' part
=
B2 B\ + B2 = 2Go + B
and with the aid of some analytical manipulations it is possible to transform
equation (10) into a complete tensorial relationship of the form
ACT,,- = DijU Aekt (11)
containing eight parameters to be defined.
1.3.2.3. Loading criterion and loading surface. Characteristic of the
model is the presence of two separate loading functions, one for deviatoric
and one for volumetric loading. Bearing in mind the form of the constitutive
relationship equation (10), the loading functions are cast in the form
volumetric loading AL, = (A', Ae0 + A'2Ae + A'3Aw)
deviatoric loading ALd = (B\ Ae0 + B'2Ae + B'3Aw) (12)
where the primed quantities are suitable combinations of the unprimed ones.
Unloading occurs when AL, and/or ALd < 0, and the separate effects are
directly accounted for in the stress-strain relationship, since this latter can
be expressed as a function of ALy and ALd a s given by equation (12).
It is important to note that since the loading criteria are formulated in
terms of strain rate quantities, they are applicable also for the softening
range.
1.3.2.4 Determination of material parameters. A detailed explanation
of the procedure for calibrating the eight parameters to match the behaviour
15
RC ELEMENTS UNDER CYCLIC LOADING
= 0 = yield surface).
Geometrically, the condition above states that during plastic loading the
stress increments can only move on a plane tangent to the yield surface. A
negative value of the loading index implies that the stress point is moving
toward the interior of the yield surface, and no plastic strains occur in this
case.
A fundamental postulate of perfect plasticity is that plastic flow takes
place along the direction of the normal to the yield surface (normality rule).
The magnitude of the plastic flow is undetermined.
Perfect plasticity is unsuited in many respects for describing the
behaviour of concrete. One of the reasons for this is the gradual
16
CONCRETE IN COMPRESSION
/ f a ) = *(ep) [ep
The consistency condition now writes
one obtains
whence for the single components deL taking into account the normality
assumption, one has
df
. P= J day Qf_ ,
damn damn)
17
RC ELEMENTS UNDER CYCLIC LOADING
18
CONCRETE IN COMPRESSION
Uniaxial compressive
Fig. 9. Non-uniform " loading path
hardening plasticity model Failure surface
(from Han and Chen13)
Loading surface
19
RC ELEMENTS UNDER CYCLIC LOADING
Tension-tension
,= 0
Still, other no less important features remain beyond reach of this model,
K(om)
20
CONCRETE IN COMPRESSION
The concept of a bounding surface has been used for modelling the
behaviour of such diverse materials as steel,36 soils,37 and concrete.1 l6 ' 38
Depending on the specific application, the concept may require adaptations
and additions with respect to those basic features which have been
mentioned. For example, metal plasticity requires an enclosed surface to
define yielding of the material, while in some other models a second
internal surface has been introduced with a meaning and a purpose
analogous to that of the loading surface in hardening plasticity. The general
formulation of bounding surface plasticity is briefly summarized in the
following. The bounding surface is given by the expression
where the bar over the a-^ indicates that they belong to F = 0. The state
variables qn accounts for the effects of the load history and are functions of
the accumulated plastic strains. Following Dafalias and Herrmann,37 the
plastic strains rate is given by
where L is the loading index (positive if plastic strains are to occur) and Ry
is the gradient of an appropriate potential function. The loading index L is
taken as proportional to the internal product of the gradient of a loading
surface by the vector of the stress increments
L N d
The proportionality factor is the inverse of the plastic modulus: hP. The
definition of L given above applies to perfect as well as to hardening
plasticity and has the advantage of correctly indicating 'loading' when the
material is in the softening range. (In this case Ny d<7,j becomes negative,
but since ffl is also negative their product is positive.)
Note that L is not only a loading index, but gives also the modulus of the
strain increment vector.
21
RC ELEMENTS UNDER CYCLIC LOADING
22
CONCRETE IN COMPRESSION
The components of the plastic strain rates, for each one of the two loading
surfaces, are expressed, according to the general formulation previously
given, as
A$ = LNl (13)
where L, loading index, (>0 for plastic loading), is given by
L = ±N»d*u (14)
with Njj = Q£-, with / = 0 = loading surface and H = generalized plastic
modulus.
L also gives the intensity of the plastic strain increment, the direction
being given by Njj which may be different from the normal to the loading
surface Nj- (non-associated flow rule).
Damage variables are included in the formulation to represent the effect
of progressive damage, mainly due to microcracking. These internal
variables may be directly identified or are connected to a scalar ep, which is
a measure of an effective accumulated deviatoric plastic strain
4^ '/2
The elastic part of the total strain is expressed as
(15)
«-<*;
where ^(a/y, &), the complementary free energy density, depends on the
damage variable k, which, in turn, is related to ep.
The total strain rate may be written as
de;y = d 4 + de? = de)/ + d6^ + de? (16)
where the elastic part has been divided into a reversible contribution
d9
= Cijlm dalm, Cijlm{k) = (17)
ooij oai
23
RC ELEMENTS UNDER CYCLIC LOADING
(19)
L
When plastic loading is active but the two surfaces are not in contact the
stress point belongs only to the loading surface. An image stress point on the
bounding surface, atj, is defined through a proper mapping rule, together
with a distance S between 07, and its image <7,y.
The plastic modulus in equation (14) is then defined as
H = H +f{8) (20)
H is still given by equation (19), computed at al}, and/(5) is a positive
decreasing function of 8 such that/(O) = 0 (which implies H = H when the
two surfaces are in contact).
Of the two plastic mechanisms used to formulate the concrete model,
only the first one uses the concept of bounding surface.
The loading surface/! = 0 and the bounding surface F = 0 are in this case
two similar cones in the space of the principal stresses, with the same origin
and same axis, that is the hydrostatic axis. Equivalently, the three stress
invariants / = an, J = (1/2 sy sp) and the Lode angle 6 in the n plane, may
be used to define them (Figs 12 and 13).
The loading surface always passes through the current stress point. When
plastic loading associated to /1 = 0 occurs, the computation of the plastic
strain rate requires the definition of the mapping rule between the current
stress point and its image on the bounding surface; this simply consists of a
F= 0
It 0
Fig. 12. Loading surface and bounding surface in Fig. 13. Loading surface and bounding surface
the I-J plane'6 in the it plane is
24
CONCRETE IN COMPRESSION
f, = 0
25
RC ELEMENTS UNDER CYCLIC LOADING
1.5. Models based These models aim at reproducing, at a phenomenological macrolevel, the
on the mechanism of internal mechanism of microcracking and growth of defects and micro-
microcracking or cavities which arise in brittle materials, such as concrete, as deformation
elastic damage progresses. The weakening of material in these models is revealed by a
decrease of stiffness without permanent deformation when the material is
unloaded. This mechanism may also be termed elastic damage.
This idea has been implemented in two formally distinct ways, which are
conceptually very much related. The first, proposed by Dougill,23 and
known as the theory of progressively fracturing solids, has been developed
in strong analogy with the theory of plasticity. The second, based on the
explicit introduction of damage variables in the constitutive relation,
belongs to the vast field of investigation known as continuous damage
mechanics.
26
CONCRETE IN COMPRESSION
In order to obtain the flow rule for the material, consider the incremental
version of equation (21)
dehk k thk (23)
This expression suggests that the stress increment day is the sum of an
elastic component do? given by the first term in equation (23) and a
component — da1- = dCyhk ehk, caused by the change in stiffness and termed
the fracture stress decrement. With reference to a unixial case the two terms
are sketched in Fig. 15.
Application of Il'iushin's postulate of plasticity and of the restriction that
the total work done in a small cycle of deformation should be non-negative
leads to the requirement that
AW{ = ]j de,:/ > 0 (24)
The flow rule then follows by analogy with the theory of hardening
plasticity
da ij = d A — (25)
oe-y
and the consistency condition of the loading surface may be used to
determine the scalar function dA. Assuming, for example, the simple form
of F, F = F(e,y) — H(Wl), which depends on a single scalar H, expressed
as a function of the fracturing work W1, we get
(26)
dey dWl
where dWf = 1/2 da\: ey is the shaded area indicated in Fig. 15.
Using all the previous relations finally we get
—2 dehk (27)
dH (dF/demn)e
which is the general constitutive equation for a progressively fracturing
solid.
An analytical expression of the function dVfldH has been deduced by
Spooner and Dougill on the basis of data obtained by uniaxial tests. From
the previous equations, and the additional postulate that the change in
stiffness caused by an increment of deformation is independent of the
deformation path, Dougill has obtained an expression for the stiffness
27
RC ELEMENTS UNDER CYCLIC LOADING
e
U =fij(a'hk) (30)
the same constitutive law may hold for the damaged material by the simple
substitution of an effective stress tensor aeg in place of a.
Assuming that linear elastic behaviour governs the effective stress-strain
relationship
= Ee (31)
then the nominal stress-strain relationship reduces to the simple secant
relationship
a = E(l -d)e = E'e (32)
where 0<d< 1.
This shows a very simple and direct relationship between damage and
stiffness reduction; such a relationship may be useful for identifying the
damage level through experimental measurements of elastic moduli at
unloading. These simple considerations introduce the simplest form of an
elastic damage constitutive relationship. This needs to be completed with
the definition of an evolutionary law for the damage parameter as a function
of the stress and strain history.
The outlined version of an elastic damage model, with a single damage
parameter, is the simplest constitutive relationship of this category. In spite
of its simplicity it has been implemented by various authors, Resende and
Martin,20 Krajcinovic and Fonseka,21 Mazars,22 with some success in
describing the behaviour of concrete.
Also evident is the close relationship of this elastic damage formulation
with the progressively fracturing material introduced by Dougill;23 the only
difference being that here no fracture surface or flow rule, reminiscent of
28
CONCRETE IN COMPRESSION
e=- ^ (33)
GQ(\-d)e-G0Aee e>0
[ }
G0{l-d) e<0
As for the volumetric strain rate, ev, it is assumed to depend both on &m
and on the deviatoric strain, to make it possible to reproduce the known
behaviour of concrete dilation under shear deformation; it is expressed as
the sum of an elastic component and of an inelastic damage component
e v - eev + ^ (37)
The elastic component e® is a non-linear function of hydrostatic stress in
order to obtain the effect of progressive incompressibility of the material as
am increases; in incremental terms it is written as
am =f(eev, ev)eev (38)
The damage component is in turn expressed as
29
RC ELEMENTS UNDER CYCLIC LOADING
Dilation Increasing
hydrostatic
pressure
is larger the larger is the value of the damage variable d prior to unloading.
At large values of d the unloading stiffness appears negative, as also
experimentally observed (see Fig. 18).
The final volumetric relation is obtained upon substitution in equation
(38) of the difference ev - e^, in place of e®
om = f (
and hence
30
CONCRETE IN COMPRESSION
1.6. Models based With reference to the actual behaviour exhibited by real materials such as
on the association of concrete, a positive feature of models based on progressive fracturing or
plasticity and elastic damage lies in their natural ability to reproduce both the hardening
microcracking and the softening range, together with stiffness degradation. In spite of these
mechanics qualities, the complete reversibility proper to these models makes them
unsuited for describing the behaviour of concrete, which is known to present
unrecoverable deformation. Since this phenomenon is well treated by
plasticity theory, the idea has arisen of combining the plastic and the
fracturing mechanisms into a single composite or mixed model. This has the
additional appeal of reproducing at macrolevel the underlying mechanical
processes of plastic slippage and microcracking both present in the real
behaviour. The coupled mechanism is illustrated qualitatively in Fig. 19.
Figure 19 represents in (a) a typical elastoplastic-hardening relationship
with unloading stiffness equal to initial stiffness; in (b) an elastic damage
relationship with stiffness degradation as deformation progresses, but no
residual strains; and in (c) a more realistic relationship where both stiffness
degradation and residual strains are present.
The first authors to exploit this possibility were Ba2ant and Kim.42 Many
more proposals have followed. A very brief account of these is given here
with reference to Fig. 20 which attempts a unitary frame in which such
proposals can be accommodated. The boxes relative to elastoplastic and
continuous damage models have already been commented on in specific
parts of this Chapter. Here a short comment is given on the most recent
proposals of the mixed or composite type.
In the pioneering work by BaZant and Kim a classical elastoplastic
mechanism with a Drucker and Prager type surface and mixed type
hardening is coupled with a progressively fracturing mechanism described
according to the Dougill23 theory. Two loading surfaces are introduced, one
of which controls the plastic deformation and the other the fracturing stress
relaxation.
A similar idea has been used by Han and Chen.43 They associate to a
Dougill type fracturing mechanism an elastoplastic model which has the
peculiarity of being defined in strain space as an alternative to the more
usual formulation in stress space. This allows the adoption of a single
loading function for the components of the stress increments of the two
types.
In recent times classical plasticity has found some effective association
with the continuous damage theory. Among the papers devoted to the
formalization of the different categories which can be generated by this
approach a special mention has to be given to papers by Krajcinovic ' and
Lemaitre.44
A specific proposal which refers to concrete has been presented by
31
RC ELEMENTS UNDER CYCLIC LOADING
(Bazant-Bhat [7],
Han-Chen) [43]
(Resende) (Lubliner ef a/.)
Progressive [40] [2]
fracturing
Continuum (Krajcinovic,
damage theory - Lemaitre,...)
[41.44]
Continuum damage
mechanics
32
CONCRETE IN COMPRESSION
33
RC ELEMENTS UNDER CYCLIC LOADING
CM = CK = CK(1 - (45)
where Cy, CD, CM, CK represent the secant moduli Cy = Fy(ev)/ev,
CD = F D (e D )/e D , CM = FM(eM)/eM, CK = FK(eK)/eK; C v , C&, CM, C£
are the initial values of Cy, CD, CM, CK; and cjy, CJD, WM, UK are the
volumetric, deviatoric and shear damage on the level of the microplane. The
shear microplane moduli are equal for both directions (CM = CK). The
damage microplane functions for the virgin loading are approximated as
34
CONCRETE IN COMPRESSION
behaviour there is no damage (u;y = 0) and the response for virgin loading
is described by
(47)
A.3. Unloading, In the previous version of microplane model45 it was possible to represent
reloading and cycling only the first unloading, and even that with considerable errors compared to
loading experiments. To model unloading, reloading and cycling loading and do so
even for general triaxial stress states, more complex rules on the microplane
level are needed. After a number of trials, the following unloading-
reloading rules, which are different for each microplane strain component,
have been chosen and verified.
In contrast to virgin loading, the stress-strain relation must be in the
incremental form and is as follows
day = Cydey; dot) = (48)
where Cy, Co, CM, CK represent unloading-reloading tangent moduli
which are defined as follows (see Fig. 22)
\ -a)a/(e-ei)
e > eP; = e P - aP/C0 + {3(e - e P ) (49)
e < ep; = 0
where, <7p and ep denote the positive or negative peak stress and the corres-
ponding strain for each microplane component, using values <rp~, ep" and dp
ep for the positive and negative peaks; a and /? are empirically chosen
constants with values between 1 and 0. The general shape of the unloading-
reloading rules for the microplane components is shown in Fig. 22.
"D
op 7
(b)
Op
CM.K
Fig. 22. Stress-strain
relations on microplane:
(a) volumetric tension-
compression; (b) deviatoric;
(c) volumetric compression-
tension; (d) shear (C)
35
RC ELEMENTS UNDER CYCLIC LOADING
A.4. Incremental Since on each microplane the kinematic constraint is used, the microplane
macroscopic stress- stresses cannot be exactly equal to the resolved components of the
strain relations macrostresses (static constraint). The micro-macro equivalence of the
stresses can be enforced only approximately, which is done by means of the
principle of virtual work. This yields an incremental macroscopic stress-
strain relation of the form
= CijrsAers - (50)
in which
- -z- -
1 / kn/i r« - CD)
7s L ->
- ( msnr)C u
A.5. Rate effect Experimental evidence indicates that concrete stiffness, strength and
ductility are rather sensitive to the deformation rate. This is known as the
rate effect. This effect is no doubt caused by creep in the bulk of test
specimen as well as time-dependent rupture of bonds in the fracture process
zone, which both cause stress relaxation. A simple rheological model,
UMIUII
- I1IIUIII
IIUUUI
consisting of one elastic spring coupled in series with a parallel coupling of
a damage unit and viscous dashpot (Fig. 23), is introduced for each
microplane component — volumetric, deviatoric and shear.49 It can be
shown that the viscosity of the dashpot must be highly stress-dependent. It
41 should be also noted that a similar rheological model for the rate effect in
damage (but without the stress-dependence of viscosity) has recently been
Fig. 23. Rheological model
introduced by de Borst.50 Dropping the subscripts V, D, M and K, the stress-
for each microplane
component strain relation for each microplane component may be written as
(53)
where e is the total microplane strain including the flow (creep) strain, the
upper dots denote the time derivatives, Ct and Co are material parameters.
A.6. Numerical In every iteration of the load or time step r, based on the kinematic
algorithm in each load constraint (equation (41)), the known macrostrains e//r and their known
step increments Ae/,r can be used to calculate the strains and strain increments
on each microplane. Furthermore, the known values of CN = ey + CD, CM,
CK, ACN = Aey + ACD, ACM and ACK, are used to calculate the stresses on
36
CONCRETE IN COMPRESSION
A.7. Numerical The behaviour of the general microplane model is demonstrated on the
examples material level, using only one uniformly strained finite element, loaded in
three different ways — tensions, compression and shear. Cyclic behaviour
of the three-point bend and compression specimens in plane stress is also
simulated. It should be noted that the plane stress state is a relatively
complex state to simulate with the microplane model. The microplane
model is a full three-dimensional model and in plane stress finite elements
the lateral strains (out-of-plane strains) need to be calculated from the
condition that the lateral stresses are zero. The basic material parameters
used are: initial Young's modulus E = 20000MPa, Poisson's ratio // = 0-18.
The microplane material parameters are chosen as follows a = 0005, b =
0-043, p = 0-75, q = 200, «, = 000007, e2 = 0-0020, e3 = 00020, e4 = 0,m =
0-85, n = 2-25, k = 2-25. Most of these values are the same as in Bazant and
Prat,45 except m, n and k (which have been adjusted so that descending
2-5- 40-
Displacement control: Displacement control:
cte/df = infinite dt/df = infinite
dt/df = 005 s" 1 de/df = 0-1E-4/0-1E-4 (1/s)
2-0-
de/df = 0025 s"' 30-
20-
10-
0-5-
37
RC ELEMENTS UNDER CYCLIC LOADING
1-5-
i
0-5-
-0-5-
Monotonic loading
-1-5
00002 00004 00002 00004
Strain Strain
Displacement control:
2-5- df/df = 0025 s '
Monotonic loading
00002 00004
Fig. 25. Calculated concrete behaviour in uniaxial
Strain cyclic tension with three different rates
40
Displacement control: Displacement control:
no rate effect df/df = 1 0 s" 1
Monotonic loading
38
CONCRETE IN COMPRESSION
16
Displacement control: Displacement control
no rate effect df/df = 1 0 s '
8-
w 0
55
-8"
-16
-0010 -0005 0 0005 0010 -0010 - 0 005 o 0005 0010
Strain Strain
Fig. 27. Calculated concrete behaviour in cyclic shear with two different rates
The present results indicate that the cyclic microplane model can
realistically predict the cyclic behaviour of concrete in different stress-
strain conditions using the same material parameters. The results shown are
only one among many possible material responses that can be modelled.
The microplane model parameters, the parameters that control loading and
unloading on each microplane in each direction as well as the rheological
model parameters are here kept constant. Varying these parameters as
material functions of the strain state and history it would be possible to
obtain a number of different kinds of responses.
Furthermore, the results of the analysis indicate that the rate effect has
significant influence on the stress-strain response. If the rate of loading
decreases, the peak stress also decreases while the post-peak descending
stress-strain curve becomes less steep. The present model can predict the
drop of stresses after repeated unloading-reloading cycles in post-peak
strain softening. As will be demonstrated later, this effect is significant
when simulating the structural behaviour.
Finally, to demonstrate the response of the present model in finite
element applications, the behaviour of the three-point bend (Fig. 29(a)) and
compression specimens (Fig. 29(b)) is analysed. Four-point isoparametric
quadrilateral plane stress finite elements with four integration points are
used. In both cases, symmetric response is assumed (Fig. 29). The material
x = Elem
(b)
Fig. 28. Comparison of cyclic uniaxial compression between: (a) test data (Sinha et al: ) and (b) calculated results
39
RC ELEMENTS UNDER CYCLIC LOADING
F/2
0 Q
3| -
1 E
g
II
Thickness = 100 mm
1L
1/2 ,
1 ?i (i r. r
ft/2 = 150 mm j
Thickness = 1 0 mm
(a) (b)
Fig. 29. Geometry and finite element mesh for (a) three-point bend specimen, (b) compression specimen
40 32
Displacement control:
do/df = 62-5 mm/s
Fu = 25-46 kN
0-94 Fu
30- 24
16-
40
CONCRETE IN COMPRESSION
A.8. Conclusions The conclusions reached in this Appendix are listed below.
(a) The present general microplane model is capable of realistically
predicting the behaviour of plain concrete under a broad range of
strain states and histories using the same material parameters. The
model is fully three-dimensional. Together with the non-local strain
concept, the model can be effectively used in finite element
simulation of the failure process in the concrete structures under
rather general loading types and histories.
(b) The rate effect can be implemented in the microplane model by
combining damage with Maxwell's rheological model. The model is
able to give an approximately correct hysteretic loop area and a steep
initial unloading slope. For shear, it exhibits the pinched form of
hysteretic loops. Generally, increasing the rate increases the concrete
strength and the post-peak descending branch of the stress-strain
curve becomes steeper.
(c) The examples indicate that the cyclic microplane model, together
with rate effect, predicts material damage due to repeated unloading-
reloading cycles quite realistically.
41
2. Concrete in tension
2.1. Introduction There is a substantial difference in the tensile stress-strain behaviour of
concrete and, for example, steel. To understand the behaviour of concrete,
this difference has to be examined. Steel, having passed the yield point!
shows a strain-hardening behaviour. Referred to the real cross-section of the
material, the strain-hardening behaviour lasts up to the ultimate load.
Concrete shows no — or almost no — strain hardening. Strain softening
starts immediately after the elastic limit is reached. The stress-strain
behaviour is governed by micro- and macrocracking within a certain crack
band.
2.2. The tensile The tensile behaviour of concrete is illustrated in Fig. 32. A deformation-
behaviour of concrete controlled centric tension test may be performed on the specimen shown on
the left of Fig. 32. The specimen is loaded with the force P. The total
deformation is measured over the length /. At the left side of Fig. 32 the
specimen is plotted for the load steps A, B and C. The load steps are also
marked in the load-deformation curve at the right side. Load step A is
before peak load, load step B at peak load and load step C after peak load in
the descending branch of the load-displacement curve.
Already before the peak load is reached, some microcracking occurs (Fig.
32(a)). As the microcracking is uniformly distributed at the macrolevel, a
uniform strain over the length of the specimen may be assumed. The strain e
is plotted over the length of the specimen in Fig. 32 right next to the
specimen.
Immediately before the peak load, an accumulation of microcracks occurs
at the weakest part of the specimen. At the macrolevel this leads to an
additional strain over the length h of this weak part. A crack band of width h
develops (see Fig. 32(b)).
Having passed the peak load, the crack band localizes more and more.
The crack band width diminishes, and the deformation within the crack
band increases. The final failure occurs due to one single crack.
The total deformation of the specimen may be split up in the bulk
deformation — which is almost linearly elastic up to the peak — and the
deformation of the crack band.
W = Al-el
p
Fig. 32. Tensile behaviour of *
concrete (a) (b) (c)
42
CONCRETE IN TENSION
2.3. The fictitious Hillerborg58 introduced the fictitious crack model, where he collected the
crack model deformation of the crack band into the crack width w of one single
'fictitious' crack. The relation between the crack width of the 'fictitious'
crack and the stress is the a-w relation. Two mechanisms contribute to the
stress transfer over the crack. As mentioned, it is actually a fictitious crack:
the crack width is the collected deformation of a band of microcracks.
Within this crack band material bridges transfer the load. After the
formation of a real single crack, the stress transfer is possible due to
aggregate interlock. In most cases a crack will run along the interface
between the aggregate grains and the cement matrix. The grains are pulled
out of the matrix and due to this, friction forces between grains and matrix
occur. The grains act like friction blocks and transfer friction forces over the
crack.
2.4. Models for the The bulk behaviour is described by the a — e relation of Fig. 32(c). This
bulk behaviour relation is linear almost up to the ultimate load. The modulus of elasticity is
equal to the initial tangential modulus of elasticity in compression.
Poisson's ratio is in the order of v - 015-0-25. Just before the ultimate
load is reached, the a — e relation bends off from the linear behaviour. For
the state of loading where the non-linear behaviour starts, the available
experimental results show great differences. In some experiments the non-
linearity starts at about half of the ultimate load, in other experiments the
curve is linear up to the peak load. These differences in the experimental
results are caused by different boundary conditions. Any source of non-
uniformity, like internal bending due to non-uniform cracking, eigenstresses
due to differential shrinkage and temperature, notch effects, etc., causes
non-linearities. Theoretical considerations59 show that even when the
material behaves linearly elastic up to the ultimate load, a specimen may
behave non-linearly before the ultimate load is reached. The non-linear
behaviour of test specimens is more a structural behaviour than a material
behaviour. Because of this it seems acceptable to assume linear-elastic
material behaviour up to the ultimate load, for both the loading and
unloading part.
2.5. Models for the A large number of experimental results are available for the a — w relation
stress - crack width in the case of monotonic loading. Fig. 33 shows some suggestions for the
relation formulation of a — w relations.5 ~63
Except for the linear one, which is often chosen because it is very simple, all
a — w relationships show the same behaviour, qualitatively. The first branch
of the relation is very steep, whereas the second branch is more flat. Since most
practical applications are not very sensitive to the exact shape of the a — w
relation, the choice of any of them is mainly a matter of practical convenience.
Contrary to the monotonic loading, only few models are available for the
case of cyclic loading. Fig. 34 shows the models of Rots et al.,64 Gylltoft,65
Reinhardt et al.66 and Yankelevsky and Reinhardt.67 These models can be
compared with the experimental results from Reinhardt et al. shown in Fig.
35. In spite of some differences in details, the models presented in Fig. 34
show some common features. For the description of the hysteresis loop they
use polygonal shapes, which are independent of the shape of the envelope.
Rots et al., Gylltoft, and Yankelevsky and Reinhardt use a a — e relation
(crack band model), Reinhardt et al. use a stress-crack width relation
(fictitious crack model).
The simple model used by Rots et al. (and others) connects the point on
the envelope where unloading starts straight with the origin of the
coordinate system. Unloading and reloading run on the same line, i.e.
there is no hysteresis loop. After running back and forth the envelope
43
RC ELEMENTS UNDER CYCLIC LOADING
function is valid again. This and similar models approximate the real
behaviour very roughly. They should be used only if cyclic effects have
negligible relevance in the overall behaviour.
Gylltoft uses a set of lines parallel to the initial slope of the curve to
describe the cyclic behaviour. The path O A B C D in Fig. 34 describes
monotonic loading. The microcracking is considered to start when the
tensile strength is reached at point A. If unloading is to be performed after
this stage, e.g. from point E, the gradient of the unloading path is chosen to
be equal to the linear elastic gradient also used for the initial loading path.
la-
44
CONCRETE IN TENSION
When zero stress is reached the remaining strain is denoted ef. When
changing over to compressive stresses, partial closing of the existing
microcracks occurs, resulting in the strain-gap FG. The size of the gap is
chosen to be linearly proportional to ef(FG = gcet)- The parameter gc is
considered to be governed by the properties of the concrete. The microcrack
will not close completely. This constraint can comprise effects due to
imperfect matching of irregular crack surfaces, particles within the
microcracks which have come loose, and irreversible deformations. The
path GH in compression is parallel to path OA.
When reloading into tensile stress takes place, a partial opening of the
existing microcracks occurs, resulting in a strain-gap GF' = gogcei; go is
considered to be a material parameter. Being back in a state of tension
again, the path F'E' is parallel to path OA. When approaching the virgin
curve once more, the remaining part of the stress-strain curve is lowered in
order not to exceed the stress level at which the virgin curve was left earlier.
Point E' and E are situated on the same stress level. After this it is proposed
that path E'C should be followed until ultimate fracture occurs. Path E'C' is
parallel to path EC.
The model of Gylltoft is a better approximation to the real material
behaviour than the first simple model, but it must be pointed out here that
the shifting of the envelope line EC to the line E'C contradicts experimental
results. After running through the hysteresis loop, the load-displacement
curve of cyclic tension tests reaches the envelope curve of monotonic
loading again.
The model of Reinhardt et al.66 uses the a-6 envelope curve of equation
(54) suggested by Cornelissen et a/.63 which was already shown in Fig. 33
C, = 3
C 2 = 6-93
wc = 160/zm
w =6
Loading and unloading steps are modelled as follows.
(a) After the peak value, the stress follows the envelope curve up to a\,
where unloading starts.
(b) Since the subsequent path is partly determined by the reloading
displacement £3, first 03 and £3 have to be calculated. The stress
45
RC ELEMENTS UNDER CYCLIC LOADING
drop, ACT, can be obtained from Fig. 36. The diagram presents three
lines obtained from different experiments. Type b represents tensile/
tensile loading with the lower stress almost zero. Type c and d belong
to loadings where the lower stress is compressive —015/ct or —/ct,
respectively. The stress drop is dependent on <J\,8\ and the lower
stress oi- The stress 03 has to be calculated as 03 = o\ — ACT. The
crack opening 63 must be calculated fromCT3and equation (54) by an
approximate procedure.
(c) Unloading takes place at constant <5| until line A is reached. Further
unloading follows the lines A, B and C until equilibrium stress 01 and
crack opening 82 are reached. Line A is given by equation (55) and
line B by equation (56).
-=1f--0 (55)
/ct 6 \S3 )
CT
- 1 (56)
It 0-3&
46
CONCRETE IN TENSION
47
RC ELEMENTS UNDER CYCLIC LOADING
The distribution function 3>a(w) must describe the rapid failure of the
primary bearing mechanism. The following formulation was chosen
'w» (58)
48
CONCRETE IN TENSION
100
Fig. 38. cr-w relation fitted to '
experimental results Experimental
Material model
0-75-
0-50- <
0-25-
friction blocks and two springs as shown in Fig. 39. The symbols used are
defined as follows.
F\ = first spring
F2 = second spring
Ri = first friction block
R2 = second friction block
w = total displacement (= crack width)
M>\ = displacement of friction block Rx
w2 = displacement of friction block R2
P(w) = total load (= load in F{)
P\ = transferable friction force in Rx
P2 = transferable friction force in R2
K\ = spring stiffness of F\
K2 = spring stiffness of F2
If the mechanism shown in Fig. 39 is forced to a certain displacement, a
relation between load and deformation of the type shown in Fig. 40 is
obtained.
1-0-1
I hwwvH h/vm P
— - w2 — - w, w
49
RC ELEMENTS UNDER CYCLIC LOADING
Wl = Wx--± (59)
K2
Since a crack does not open until the tensile strength of concrete is reached,
the mechanism must be able to transfer the maximum load for w = 0. This is
achieved by adjusting the displacements wx and w2 at point 1 in such a way
that both friction blocks R\ and R2 transfer exactly the maximum friction
force P\ or P2, respectively. At point 1 the following equations are valid
P =Pi+P2
w =0
L.
w2 = w , - ^ (60)
K2
After reversal of the deformation direction at point 2, spring Fi is
compressed until the friction force at /?i is equal to —Pi. When this occurs,
point 3 is reached. Between point 2 and point 3 both friction blocks are
motionless. Between point 2 and point 3 the following equations are valid
/.,, ... \ is-
— \W — W\)l\.\
W\ = W,
w2 = w\ (61)
where w\ - displacement W\ at point 2
w\ = displacement w2 at point 2
From point 3 on, R { slides in a negative direction. The stiffness of spring
F2 is now adjusted in a way that on achievement of a crack width w = 0, the
point 4* on the load axis will be reached. Point 4* has the ordinate — 2P\.
This leads to
p(w = 0) = - 2 P ,
w =0
from which
50
CONCRETE IN TENSION
P =(w-wl)Kl
w2K2
Kl+K2
W =
2
At point 4 the deformation direction is changed again. Spring Fx is now
elongated until friction block R\ starts sliding in a positive direction at point
5. Both friction blocks Rx and R2 are motionless up to point 5. Between
point 4 and point 5 the following equations are valid
P ={w-w\)K\
W\ =w\
(64)
From point 5 on, R\ slides in a positive direction until point 6 is reached,
which is the starting point 2 of the hysteresis loop. The displacement wx
again results from equilibrium conditions of block R{. Between point 5 and
point 6 the following equations are valid
p = (w- wx)Kx
W2K2 + wKi - />,
K{ + K 2
1
w2 (65)
From point 6 on, both friction blocks slide in a positive direction. Now the
same relations as between point 1 and point 2 are valid again. If the
deformation direction will be changed once more, another hysteresis loop
will result. From the definition of the spring stiffness K2 in equation (62),
point 4* with the coordinates (w = 0, P = —2P\) will be reached in every
hysteresis loop, if the deformation w = 0 is achieved.
51
RC ELEMENTS UNDER CYCLIC LOADING
"b2,0 \ o b1.0 \
width. With increasing opening of the crack, the elements fail and will not
be reactivated when the crack is closed again. The distribution functions
must, therefore, refer to the maximum crack width (= wmax) ever reached in
the deformation history. The following equations are proposed for the two
functions <&a(w) and
_ , x /'"'max
$ a ( w ) = e v wa
= e (66)
and therefore
max
—ac 0 ( — 1) for negative deformation direction
<Tc(w) = ' \ w )
max
—<*c 0 ( — 1) C2 for positive deformation direction
' V w /
(68)
Group A
0-4
0-2
Group B
30 10 20 30
w.ftm
Fig. 42. Hysteresis loop of -0-2J
element groups A and B
52
CONCRETE IN TENSION
± -0-4-
-0-8-
30
53
RC ELEMENTS UNDER CYCLIC LOADING
(b)
2.7. Continuous The model gives a description for the stress-crack-opening relation, while
function model crack opening is defined according to the fictitious crack model by
Hillerborg et al.5S It consists basically of continuous functions and is,
therefore, called the 'continuous function model'.
54
CONCRETE IN TENSION
Aa
unloading curve (I), the gap in the envelope curve (II), and the reloading
curve (III) (see Fig. 47). It has been decided to use only characteristic points
in the a — w relation (/",, wc, weu, creu, 0L) as variables in the expressions.
The parameter wc is the crack opening where stress can no longer be
transferred and is defined as 5,14GF//t. The expressions chosen are based on
a close inspection of the experimental results. The softening relation is
described by the following equation
/ UA31 / W\ W 3
(69)
(ci —] expl-c 2 —I ( 1 + q ) exp(-c2)
0-014 fa - - -0-57, 1 —
(70)
When reloading starts from a lower stress level CTL, the gap in the envelope
curve (II) is known by the expression for winc
Win
• = 0-l^i fa 1+3 (71)
The coordinates of the returning point at the envelope curve (wer, <rei) can
now be found with
= Weu + (72)
and equation (69). The reloading curve (III), starting from the point at the
lower stress level (H>L, <TL) up to point (wer, o~eT) at the envelope curve, is
determined by
1 / w - wh w-wL
1 - 1 -
w fir -
(73)
C 3 +l
55
RC ELEMENTS UNDER CYCLIC LOADING
£3 = 3 3 1 - I—
ft
-1
c4 = 2 3, / t - ' + 0-5
/t
(75)
= 0
Envelope
— . ^ curve
56
CONCRETE IN TENSION
to
Q.
S
55
-1- I, = 3 MPa
G F = 100 J/m2
n-c = 171 10 6 m
-2-
10 15 20 25
6
Crack opening w: 10 m
(77)
57
3. Reinforcing steel
This Chapter presents a brief review of reinforcing steel models which can
be used in the finite element analysis of reinforced concrete structures
which are subjected to large alternate actions. Since a short summary of
models proposed before 1983 was presented in Bulletin d'Information
161,72 the present report mostly addresses models published since then and
gives only brief evaluations about earlier models.
3.1. Introduction Numerous models of the hysteretic stress-strain behaviour have been
published.73"77 Most of these are phenomenological in that they attempt to
describe the behaviour by parameters which cannot be directly derived from
the physical behaviour at the microscopic level. Even though such physical
models have been published, the large computational effort associated with
their use in the finite element analysis of large reinforced concrete structures
excludes them from consideration in the present report.
It is of interest to note that most analytical models proposed to date have
attempted to address the hysteretic behaviour of steel in general, rather than
reinforcing steel in particular. The range of cyclic strain history, to which
reinforcing bars are likely to be subjected, differs significantly from that of
structural steel members in that compressive strains are not as large as
tensile strains. This is caused by an interplay of bond deterioration and crack
closure, which prevents reinforcing bars from excessive yielding in com-
pression before spalling of concrete cover has occurred. As long as a section
is uncracked or the crack previously formed is closed and the concrete cover
has not spalled, compressive forces are largely carried by concrete, while, if
the crack is open, bond deterioration in the vicinity of the crack prevents
any large force build-up in reinforcing bars. This results in small negative
steel strain increments and consequently in small compressive steel stress
increments. In either case bars under compression hardly ever enter into the
inelastic compressive strain range before the concrete cover has spalled.
Thus cyclic strain histories of reinforcing bars tend to centre about an ever-
increasing tensile plastic strain. This fact can be used in the simplification of
the proposed models of the hysteretic behaviour of reinforcing bars.
Another important consideration in the search for an appropriate model is
its numerical efficiency. Many evaluations of the steel stress—strain relation
will be required in the non-linear finite element analysis of large reinforced
concrete structures. This suggests that the model should be as simple as
possible. Counteracting the need for simplicity is the requirement of
accuracy, which is based on the observation that cracks running through the
depth of a RC member can remain open during load reversals, causing the
hysteretic response of the member to be entirely controlled by the behaviour
of reinforcing steel. In particular, the Bauschinger effect has to be accounted
for. Even though multilinear models have been proposed in the past, it
seems that simple non-linear models have a clear advantage. In any case, the
opposing requirements of numerical efficiency, i.e. simplicity on the one
hand and accuracy on the other, call for a compromise between the two
extremes.
Some of the characteristic features of the hysteretic behaviour of
reinforcing steel are shown in Fig. 50. To retain accuracy the proposed steel
models should account for the non-linear monotonic envelope, in particular,
the onset of strain hardening, the Bauschinger effect and the isotropic strain
hardening under plastic strain reversals. Several non-linear analytical
models have been proposed to date. To a large extent they are devised to
58
REINFORCING STEEL
Increase in strength
beyond initial yield stress
after reversal:
isotropic strain hardening
Degradation in steel
modulus after reversal:
cyclic strain softening
— = Ep = f(a) (81)
where ep is the plastic strain, e denotes the steel strain, a denotes the steel
stress and B? denotes the plastic modulus. Models in categories (a) and (b)
are the most popular to date, because of their simplicity. The most
successful in category (a) is the model proposed by Giuffre, Menegotto and
Pinto, while in category (b) most models are based on modifications of the
model originally proposed by Ramberg and Osgood. There appears to be
some advantage in using formulation (a) in finite element analysis. This
advantage stems from the fact that, usually, strains are used as primary
variables in section analyses or are derived first from strain-displacement
relations. In this context it is advantageous not to have to solve a non-linear
equation of type (b), but to be able to directly determine the steel stresses
from the strains. Because of this fact, attention is focused in this report on
two particular models which follow the formulation (a). The first model is a
modification of the original model of Giuffre, Menegotto and Pinto. 8081
59
RC ELEMENTS UNDER CYCLIC LOADING
The second is a new simplified model which follows both formulations (a)
and (b) simultaneously and thus can be used with great efficiency in
explicitly calculating the steel stress from a given strain (formulation (a)) or
in explicitly calculating the steel strain from a given stress (formulation (b)).
The latter problem occurs often in connection with the solution of the stress
transfer problem based on bond stress-slip relations.
3.2. Giuffre - The original model of Giuffre, Menegotto and Pinto takes on the form
Menegotto - Pinto
model
where
e - er
€ =• (83)
and
a—
a = (84)
aa — ar
Equation 82 represents a curved transition from a straight line asymptote
with slope EQ to another asympote with slope Ex (lines (a) and (b),
respectively, in Fig. 51). &0 and e0 are stress and strain at the point where the
two asymptotes of the branch under consideration meet (point A in Fig. 51);
similarly, aT and er are stress and strain at the point where the last strain
reversal with stress of equal sign took place (point B in Fig. 51); b is the
strain-hardening ratio, that is the ratio between slope E\ and EQ and R is a
parameter which influences the shape of the transition curve and allows a
good representation of the Bauschinger effect. As indicated in Fig. 51
(eo, cr0) and (er, at) are updated after each strain reversal.
R is considered dependent on the strain difference between the current
asymptote intersection point (point A in Fig. 52) and the previous load
reversal point with maximum or minimum strain depending on whether the
corresponding steel stress is positive or negative (point B in Fig. 52). The
expression for R takes the form
600
a.
60
REINFORCING STEEL
-10 10
Normalized steel strain t'
(85)
61
RC ELEMENTS UNDER CYCLIC LOADING
(c) the descending lower branch curve originating at the reversal point
with largest e value
(d) the current curve originating at the most recent reversal point.
Due to the above restrictions, reloading after partial unloading takes place
along curve (b) rather than curve (a) as shown in Fig. 53. However, the
discrepancy between the adopted analytical model and the actual behaviour
seems acceptable, since the case depicted in Fig. 53 is the most un-
favourable one.
The above implementation of the model corresponds to its simplest form:
elastic and yield asymptotes are assumed to be straight lines, the position of
the limiting asymptotes corresponding to the yield surface is assumed to be
fixed at all times and the slope Eo remains constant (Fig. 51). More complex
asymptotes have been proposed for the original model by Stanton and
McNiven.82
In spite of the simplicity in formulation, the model is capable of
reproducing well experimental results with strain histories typical of
structural steel, i.e. strain excursions of equal magnitude in tension and
compression. Its major drawback stems from its failure to allow for
isotropic strain hardening. This fact can be of importance when modelling
cyclic behaviour of reinforcing bars in RC members.
The model proposed by Filippou et al.79 retains the bilinear asymptotes of
the original model in the interest of simplicity. It introduces, however, a
shift of these asymptotes to account for isotropic strain hardening. This is
accomplished by shifting the position of the yield asymptote before
computing the new asymptote intersection point following a strain reversal.
The shift is effected by moving the initial yield asymptote through a stress
shift <rst parallel to its direction. This idea was firtst introduced by Stanton
and McNiven. The imposed stress shift of the yield asymptotes ideally
depends on several parameters of strain history. Stanton and McNiven have
suggested that the main parameter is the sum of the absolute values of
plastic strains up to the most recent strain reversal. Filippou et al. decided to
choose the maximum plastic strain as the main parameter on which the yield
asymptote shift depends. The proposed relation takes the form
where emax is the absolute maximum strain at the instant of strain reversal,
ey, (Ty are, respectively, strain and stress at yield, and 0:3 and 0:4 are
experimentally determined parameters. It is interesting to note that in a
recent study it has been found that the model of Giuffre, Menegotto and
Pinto as modified by Filippou et al. is even capable of accurately
representing the hysteretic behaviour of steel members.
Comparison of analytical with experimental results obtained from cyclic
tests on reinforcing bars84 is shown in Fig. 54. Fig. 54(b) depicts a stress-
strain history typical of top reinforcing bars where yielding in compression
is very limited, while Fig. 54(a) shows a stress-strain history typical of
bottom reinforcing bars. The accuracy achieved with the model which
includes isotropic strain hardening is seen to be very satisfactory. Since no
isotropic strain hardening of the tension envelope was introduced in the
model of Fig. 54(a) the discrepancy between analytical and experimental
results is quite significant. This problem can be solved by also allowing for
istropic strain hardening of the tension envelope. It should be noted that the
model which includes isotropic strain hardening is almost as simple and
computationally efficient as the one which does not allow for isotropic
strain hardening.
62
REINFORCING STEEL
100
50
-50
Analytical
Experimental
-100
-001 001 002 003 004 005 -001 001 002 003 004 005
Steel strain: in/in Steel strain: in/in
(a) (b)
Fig. 54. Comparison of model with experimental results: (a) bottom reinforcing bar; (b) top reinforcing bar
3.3. Simple non- The hysteretic behaviour of the model is described by monotonic envelope
linear stress-strain curves in tension and compression (Fig. 55), a load reversal curve for
model unloading and reloading in the opposite direction and a set of rules for
general strain or stress histories. The ingredients of the proposed model are
presented in detail below.
63
RC ELEMENTS UNDER CYCLIC LOADING
Linear functions are used to describe the behaviour of the model in the
linear elastic and the yield plateau range, while two options are available in
the strain-hardening region. In the first option the strain-hardening behav-
iour is described by a non-linear exponential equation. This equation is
selected so that the steel stress can be expressed as an explicit function of
strain and vice versa. Since the same equation is also adopted for the load
reversal curve of the model, further discussion is deferred until the
description of the load reversal curve. In the second option the strain-
hardening behaviour is approximated by a linear function.
The trilinear monotonic envelope is defined by the following parameters:
the yield stress ay\, the initial elastic slope Eo, the slope E{ of the plastic
yield plateau, which is always assigned a very small value, the secant slope
E2 in the strain-hardening range and the strain at the onset of strain
hardening esh. These parameters are determined from experimental results
of monotonic tension tests on reinforcing bars.
The yield strain ey and the yield stress at the onset of strain hardening ay2
are determined from these parameters
The steel stress can be expressed as a function of strain using the following
relations
a = EQ€ for e < ey (87)
a = <7yl + (e — ey)E\ for ey < e < ey2 (88)
G = ayl + (esh - ey)Ei + (e - esh)E2 for e > esh (89)
Equations (87)-(89) can be inverted to obtain the strain as a function of stress
f o r
= 7T cr<^yi ( 90 )
E
e = (^-^-) + ey for <ryi < a < ay2 (91)
l.O (95)
64
REINFORCING STEEL
= ~{hc - (96)
Ae.* = -(Ac - (97)
where
(98)
and emax is the maximum previous tensile strain at the instant of strain
reversal. The constants au a2, a3 and a 4 in equations (95) and (98) are
determined from cyclic loading tests of reinforcing bars.
a = and e* =
S(0,0)
65
RC ELEMENTS UNDER CYCLIC LOADING
Point 5 is the last point of load reversal and point E is the point with
maximum previous strain in the same direction of loading as the current
load reversal curve (Fig. 56). Under istropic strain hardening point E is
located on the shifted envelope curve. If no previous yielding has taken
place, point E is the yield point of the shifted envelope curve.
The elastic modulus of the normalized stress-strain relation takes the
form
The simplicity of equations (99) and (100) permits the use of only three
parameters a, b and A. These parameters need to satisfy two obvious
conditions, namely, that the load reversal curve passes through points 5 and
E in Fig. 56. These conditions imply that a* = 0 for e* = 0 and a* = 1 for
e* = 1. From these conditions a and b can be determined
a=b =
-e~x
Substituting the values of a and b in equations (99) and (100) gives
l — e~*
^ =^ 1 (103)
Setting the left hand side of equation (103) equal to E* results in a non-
linear equation for parameter A
- ^ - , = E* (104)
66
REINFORCING STEEL
1 -e- A c -
for a > (105)
-e~x
where
and e* =
E* =E2
(rj u — cr sh )
esh is the strain and ay2 the stress at the onset of strain hardening,
respectively. eu and au are the strain and stress corresponding to the tensile
strength of a reinforcing bar under monotonic loading. These values need to
be specified as input along with the other five parameters described earlier
for the trilinear case. In this case, however, the instantaneous tangent
modulus Et at the onset of strain hardening replaces the secant slope E2 as a
parameter (Fig. 55).
Comparison of analytical results with experimental data from cyclic tests
on reinforcing bars conducted by Ma et a/.84 are shown in Figs 57 and 58.
The analytical predictions in Fig. 57 are based on a trilinear monotonic
100
50-
-50-
• Proposed model
-——- Filippou's model
-100
-001 001 002 003 004 005 -001 001 002 003 004 005
Steel strain: in/in Steel strain: in/in
(a) (b)
Fig. 57. Comparison of simple model with trilinear envelope curve with experimental results: (a) bottom reinforcing
bar; (b) top reinforcing bar
67
RC ELEMENTS UNDER CYCLIC LOADING
100
50-
8
I 0
-50-
Analytical
Experimental
-100'
-001 001 002 003 004 005 -001 001 002 0-03 004 005
Steel strain: in/in Steel strain: in/in
Fig. 58. Comparison of model with non-linear envelope curve with experimental results: (a) bottom reinforcing bar;
(b) top reinforcing bar
100
50-
S o
-50-
Analytical
Experimental
-100-
-001' 001 002 003 004 005 -001 001 002 003 004 005
Steel strain: in/in Steel strain: in/in
(a) (b)
Fig. 59. Comparison of proposed model with modified model of G-M-P: (a) bottom reinforcing bar; (b) top
reinforcing bar
150
68
REINFORCING STEEL
69
4. Bond between concrete and steel
The bond behaviour under alternating loading is analysed following a short
description of the behaviour under monotonically increasing slip and
unidirectional cyclic loading, respectively. Various proposals for modelling
are discussed and evaluated.
4.1. Introduction The behaviour of reinforced concrete structural elements, subjected to large
alternate actions, is highly dependent on the interaction between steel and
concrete. In earthquake-resistant design of structures, economical require-
ments usually lead to the need for large seismic energy input absorption and
dissipation through large but controllable inelastic deformations. The
structural and non-structural damage incurred from an earthquake should
not lead to collapse of the structure or to the endangerment of human life.
Therefore the sources of potential structural brittle failure must be
eliminated and degradation of stiffness and strength under repeated loadings
must be minimized or delayed long enough to allow sufficient energy to be
dissipated through stable hysteretic behaviour.
In reinforced concrete, one of the sources of brittle failure is the sudden
loss of bond between reinforcing bars and concrete, which has been the
cause of severe local damage to, and even collapse of, many structures
during strong earthquakes. Even if no anchorage failures occur, the
hysteretic behaviour of reinforced concrete structures, subjected to alternate
actions, is highly dependent on the interaction between steel and concrete.
Cyclic loadings are generally divided into two categories. The first
category is the so-called iow-cycle' loading, or a load history containing
few cycles (less than one hundred) but having very large bond stress ranges.
Low-cycle loadings commonly arise in seismic and high wind loadings. The
second category is the so-called 'high-cycle' or fatigue loading, which is a
load history containing many cycles (typically thousands or millions), but at
a low bond stress range. Bridge members, offshore structures, and members
supporting vibrating machinery are often subjected to high-cycle or fatigue
loading. High-cycle loadings are considered a problem at service load
levels, while low-cycle loadings produce problems at the ultimate limit
state.
The bond behaviour under cyclic loading can further be subdivided
according to the type of stress applied. The first is repeated or unidirectional
loading, which implies that the bar stress does not change sense (tension to
compression) during a load cycle, the usual situation for fatigue loading.
The second is reversed cycle loading, where the bar is subjected
alternatively to tension and compression. Stress reversals are the typical
cases for seismic loading. In this Chapter, the behaviour of bond under low-
cycle loading is described.
4.2. Bond behaviour Bond behaviour is a combination of adhesion, bearing of lugs and friction.
under monotonically Adhesion is related to the shear strength of the steel-concrete interface and
increasing slip is the result primarily of chemical bonding. Bearing forces perpendicular to
the lug faces arise as the bar is loaded and tries to slide. In this phase
microcracking and microcrushing 89 of concrete in front of lugs are
produced. Friction is produced by bearing force on the interactional surface
and by shearing off the concrete between the lugs on the cylindrical
concrete surface at the tip of lugs.
70
BOND BETWEEN CONCRETE AND STEEL
100
71
RC ELEMENTS UNDER CYCLIC LOADING
Si 0-6 mm 0-6 mm 10 mm
s2 0-6 mm 0-6 mm 30 mm
a 0-4 0-4
Slips T
f 0-15rmax 0-40rmal
Fig. 62. Analytical bond stress-slip relationship for 'Failure by splitting of the concrete
monotonic loading (CEB93) f Failure by shearing of the concrete between the ribs.
01
1 0
0-40
0- 1
001
4.3. Bond behaviour Cyclic loading produces a progressive deterioration of bond. The
under unidirectional deterioration can be observed as a slip increment.94"96 This process can
cyclic loading lead to failure at a cyclic stress level lower than the ultimate stress under
monotonic loading.
4.3.1. Behaviour
yo
According to the pull-out test results of Rehm and Eligehausen,96 with
constant amplitude cyclic loading, if no fatigue failure of bond occurs the
slip against number of load cycles relationship in double logarithmic scale is
approximately linear (Fig. 63). Rehm and Eligehausen also concluded that if
a fatigue failure of bond does not occur, then previous load cycles do not
negatively affect the bond stress-slip behaviour near ultimate load
compared with the monotonic loading behaviour (Fig. 64). Otherwise the
comparison of test results with cyclic and sustained loading show that cyclic
loading can be considered as a time accelerator compared with a sustained
load.
Based on pull-out test results with cyclic load Balazs97'98 distinguished
three different phases of the bond fatigue process (Fig. 65). During the
72
BOND BETWEEN CONCRETE AND STEEL
beginning of the load cycles (first phase) the slip rate is decreasing then
(second phase) keeps constant and is finally followed by a rapid increase
(third phase) producing pull-out failure (Fig. 65(a)). The slip versus number
of load cycles relationship in normal scale contains concave, linear and
convex portions (Fig. 65(b)). The intermediate linear portion ends at S(T U ),
i.e. the slip at monotonic bond strength which marks the beginning of the
failure branch. Therefore, S(TU) is proposed97 as a safe fatigue failure
criterion.
Load history has a significant effect on slip. For instance93 a periodical
increase of the maximum value of cyclic force produces a higher slip rate of
the intermediate linear portion (Fig. 66).
The accumulation of bond damage is supposed to be caused by the
progressive growth of microcracks and concrete crushing in front of the
ribs. Their effect is observed by increasing slips. Thus, increase of slip can
demonstrate the damage accumulation of bond.
Test results by Balazs show three different phases of slip against number
of load cycles for constant amplitude cyclic bond tests and not a linear
relationship as the Palmgren-Miner law assumes (Fig. 67). During the first
and third phase the fatigue damage proceeds more quickly than predicted by
the hypothesis. In the intermediate linear stage the damage is slower.
73
RC ELEMENTS UNDER CYCLIC LOADING
20 40 60
Number of load cycles
Pull-out
failure
20 40 60
Number of load cycles n Bond stress: MPa
(b) (c)
74
BOND BETWEEN CONCRETE AND STEEL
• Palmgren-Miner hypothesis
Loaded end (test result)
Unloaded end (test result)
I Pull-out
I failure
100-
Load history
75-
Slip:
93%
(62%)
0-6-
-002
75%
(50%)
s
0-4- <d
Q
60%
(40%)
r' 25-
0-2-
Fig. 66. Slip increase under periodically increasing Fig. 67. Bond damage accumulation with respect
cyclic load, <j>8 def. bar, lh = 12<j>, C16 (Baldzs97) to the Palmgren and Miner hypothesis, test result:
4>8 def. bar, lh = 12<f>, C16 (Baldzs97)
75
RC ELEMENTS UNDER CYCLIC LOADING
D D D D
ac
= o2 = o3 = ac
NSC
P
(a)
— —— + --rR(x)
(c) (d)
i
(e)
Fig. 68. Model by Fehling": (a) stress-strain curves of concrete of different strength under 3-axial hydrostatic
compression; (b) modelling of local bond behaviour; (c) assumed stress distribution near lug ofrebar; (d) strut-and-
tie model; (e) unloading and reloading in plastic zone (model and simplification)
76
BOND BETWEEN CONCRETE AND STEEL
1
DK152: experiment
I
|
/ /!
i / A MmL/ /// i4
1' m%///!mm
// //Iif!ml t ¥ I'M
/ iJmAwo
M mm
If
Sir
r ' _
i
1
005 0-10 0-15 0-20 0-25
Slip: mm
(a)
DK152: model
1 / /t iff
1 r^ i mil
. — •
1 Ilkr' h
it
////
m
D-
{
UJiL
1 w
JiMWlJiw
w
Fig. 69. Bond stress-slip -1
relationships; Fehling9 :
(a) test result; (b) result of 005 0-10 0-15 0-20 0-25
calculation using model in Slip: mm
Fig. 68 (b)
77
RC ELEMENTS UNDER CYCLIC LOADING
10-
12-
ir ~-^v ^ Monotonic loading
1 8-
After 10
cycles
4-
J
i
-10 -8 -6 -4 -2 4 6 8 10 12
s: mm
--4
—8
1
--12
Series 2-4
Monotonic loading - ^ cycling between
— 16 s = ± 0-44 mm
(a)
. Monotonic loading
After 1 cycle
-10 -8 -6
Series 2-8
cycling between
Fig. 70. Bond stress-slip s = ± 4-57 mm
curves for cycling at different
maximum slips (Eligehausen
et al92)
78
BOND BETWEEN CONCRETE AND STEEL
4.4. Bond behaviour In the case of reversed cyclic loading the type of loading (load or slip
under alternating controlled), and the rate of loading (frequency) are important parameters for
loading bond strength.102103 Seismic loading, the most common case of reversed
cyclic loading, represents an intermediate step between load and slip
controlled cycles and is characterized by a wide frequency content.
Previous state-of-the-art reports closely linked to the present field are
available for reference.91'103106 Extensive studies considering several
parameters are also available. 78 ' 7992107 " 115 The main characteristics and
some modelling possibilities are summarized below.
4.4.1. Behaviour
Cycles with reversed loading produce degradation of bond strength and
bond stiffness that is more severe than for the same number of cycles with
unidirectional loading. 106112 Degradation primarily depends on the peak
slip in either direction reached previously (Fig. 70).9 Other significant
parameters are the number of cycles and the difference between the peak
values of slip between which the bar is cyclically loaded. Under otherwise
constant conditions the largest deterioration will occur for full reversals of
slip. Whenever the load cycles are limited to produce slip in one direction
only, there is not very significant degradation of the bond strength (Fig.
71 (a)).110
If the peak bond stress during cycling does not exceed 70-80% of the
monotonic bond strength r, the ensuing bond stress-slip relationship at first
loading in the reverse direction and at slip values larger than the one at
which the specimen was cycled, is not significantly affected by up to ten
repeated cycles (Fig. 70(a)).92 The bond resistance at peak slip deteriorated
moderately with increasing number of cycles.
Loading to slip values inducing a bond stress larger than 80% of the
monotonically obtained bond strength in either direction leads to
degradation in the bond stress-slip behaviour in the reverse direction
(Fig. 70(b) and (c)). As the peak slip increases, deterioration of bond
resistance is increased. Deterioration also increases with the number of
cycles, and is more extensive for full reversals of slip than for half-cycles
(Fig. 71).
-2 15
6000
4000-
2000-
o
m
-2000-
-40
-100 100 200 300 400 500 -600 -400 -200 0 200 400 600
Slip: x10~ 3 in Slip: x10 3
in
(a) (b)
Fig. 71. Comparison of monotonic and reversed cyclic bond stress-slip curves: (a) without change of slip sign;
(b) with change of slip sign; Hawkins et al."°
79
RC ELEMENTS UNDER CYCLIC LOADING
n
1+
1-2+2-3+3-4+4-5+5-
> (b)
^^ 1 i/i
Slip: mm
/
/
life 15
J 1/4 ~ -|
1
--5
4
10
r 7
/ 20-
I—10
If n
4+
3~4 4"
/
/ —15
V
(a)
80
BOND BETWEEN CONCRETE AND STEEL
0-8
— f
0-6/1/
II 11 iT/7
4
w
wV
1 M
-UL TO 111
L
- - 0-6"
/ k
j1
/ I °' U-' ' g
- 4 - If 1Iff n)flF- —p
>w>
-ii b
Ml
Ml 111
ml
it
tut,
n mi
[it n 0-
— 10
\\
TAT
— 15
/ s
3- rti
0-
T
reach the bottom of the adjacent lug (point D) the bond stress begins to
drop, and by the time the slip has reached the lug spacing only the frictional
component remains (point E).
Under load reversals (Fig. 76), the initial loading follows the monotonic
curve, but the cyclic load behaviour is sensitive to the level of slip at which
the reversal occurs. Three possible qualitative models have been proposed,
depending on whether or not inclined cracks have formed. In the first case
(Fig. 76(a)) imposing a slip reversal at an arbitrary slip value below the
level where inclined cracking occurs, results in a stiff unloading branch
(path AF) because only a small part of the slip is caused by inelastic
concrete deformations. As soon as slip in the opposite direction is imposed,
the friction branch is reached (path FH). The slope of this portion of the
curve is small because the surface of the concrete surrounding the bar is
si
I-oad h story F- k
!0-
r: N/mm^
IP Mon(,lonic
• env slotje
-15 + 2.-
y
V'1/1/1/1/ " • 1J
7
y
\/,
0-
/\ / y
/ , ^-' T
//
/
i
/T
if
•5
/
/
- — •
-M
^ ^
H}
1 UJ
pi
0-
--5
mm
25
-X
4-
>
+ T^.Sliprm m
1 ]
Monotonic
/-envelope
jjT~ ^ Reversed
i
Fig. 74. Bond deterioration /
-" i
f 10 X; ; •*<•-cyclic
under increasing reversed
cyclic slip values, test result:
r
7 /
/
s*
envelope
s
lb ru)
<p8 def. bar, lb = 6<j>, C25;
Baldzs10*
81
RC ELEMENTS UNDER CYCLIC LOADING
_± 1
Gaps A Concrete •-
Bond
crack
(a) (b)
smooth. As soon as the cracks close, the stiffness differs little from that of
the monotonic envelop (point I).
Unloading from point I, where the slip in the two directions is about
equal, the curve (path IKL) is very similar to that from the initial unloading
curve (path AFH). The major difference is that due to previous cracking and
crushing of the concrete in front of the ribs, the point where the bond
stresses begin to pick up again (point L) will be shifted to the right of the
origin. The lug will not be bearing fully until point M is reached. Further
loading follows the bond-slip curve up to the monotonic envelope.
If unloading occurs after the inclined cracks have formed (Fig. 76(b)),
and therefore near the slip at which ultimate bond stress has been attained,
the unloading path is similar to that of the first case up to point F. Since
there is more damage to the concrete, a higher frictional resistance is
mobilized (point G). When the loading is reversed the lug presses against a
key whose resistance has been lowered by inclined cracks over a part of its
length that were induced by the first half-cycle. The splitting cracks created
in the first half-cycle close at a higher load than those of the first case (point
l i l t i I 1 1 1 t t
C
r Closed old r Monotonic
Gaps Gap sF, Old crack loading
crack
M
Kf- TL
10
'Crushed concrete New crushing/ C r u s n i n 9
/ / f r o m first half cycle' /I Monotonic and new / / / s h e a r cracks < /
Sorne_crushed concrete /// loading . Monotonic
shear cracks ' f r o m loadingOC
loading
(a) (b)
i 1 i i I I 1 t
Monotonic
loading Fig. 76. Bond mechanism under reversed cyclic
loading; Eligehausen et al.
82
BOND BETWEEN CONCRETE AND STEEL
s L = (sB + Sj)/2
IN = <"H su = (Sj + s P )/2
= ar B TQ = arp
83
RC ELEMENTS UNDER CYCLIC LOADING
i = 1 - Vn
Monotonic envelope
84
BOND BETWEEN CONCRETE AND STEEL
Monotonic
loading
Monotonic '
loading
T
i = TE, rn = 0-9 i> r D ' = a0 r 0 To; = «b
= rE s M = 0-4 SG = aE r E = aE
Fig. 79. Reversed cyclic bond »L TE- TE;
85
RC ELEMENTS UNDER CYCLIC LOADING
a
TD (106)
where
SSJ, Ssji: sum of the peak slip values having a magnitude larger than
in the previous cycles for loading in tension or compression,
respectively;
sD: slip at point D;
C*TD) A-D> £TDI PD'. constants, evaluted from test results.
(d) The frictional bond resistance is assumed to be equal to TE of the
monotonic envelope and independent of the number of cycles.
(e) The bond stress-slip relationships for the reloading branch (path
MRN) and for additional cycles between fixed slip limits are very
similar to those proposed by Morita and Kaku. l08
The above model is a major improvement, because it takes into account
several features observed in experiments and it is approximately valid for
cycling between arbitrary slip values. However, in spite of being rather
complicated, it is not general. Some 20 parameters are needed to describe
the bond stress-slip relationship for cyclic loading, which have no clear
physical meaning and must be evaluated from test results. Furthermore, the
assumptions on which the calculation of the reduced envelope is based need
improvement. For example, an arbitrary number of cycles (> 1) in well-
confined concrete between s max = 2sD and smin = —2sD reduces TQ,
independent of the number of cycles, by 13%. TO, however, is reduced to
zero after eight cycles between almost the same peak slip values if only the
value of 5max is slightly increased in each cycle.
4.4.2.4. Model by Hawkins et a/.110 The monotonic envelope is
idealized by a trilinear bond stress-slip diagram. The characteristic values
of slips, bond stresses and tangents are obtained by statistical evaluation of
the test data yielding to three linear equations.
During load reversals, deterioration in the loading envelope does not
begin until the slip range exceeds s'o, where s'o is the slip for r max for
monotonic loading. The deteriorated capacity in the positive direction
Bond
r
stress
86
BOND BETWEEN CONCRETE AND STEEL
follows the line M'X^N' where that line is derived from the line M'XN" that
parallels line MN (Fig. 80). The post-peak response envelope for the
negative direction mirrors that for the positive direction. Rules for deriving
the post-peak response envelope for the positive direction are
% = V = SN" = °-5(slug + W)
sx = sx, = 0-45slug = 0 - 5 ^ - W/2)
TX, = 045r x but TX> > 0-l/c' or 2-8 N/mm2
sY = 0-50siug; r Y > 0-1 f'c or 2-8 N/mm2
T = Ty — 0-5(s — Sy), Sy < S < Sff
87
RC ELEMENTS UNDER CYCLIC LOADING
Load
histories:
Slip
I A1"
IM'lll
Monotonic envelope
Reversed cyclic
envelope
x, = 0-9 to 1 0
= 0-7to 1 0
= 0-5to10
89
RC ELEMENTS UNDER CYCLIC LOADING
relationship is like that of a rigid plastic model, the only difference being
that frictional bond resistance decreases with an increasing number of
cycles. A similar behaviour as described is followed if the slip is reversed
again at point S (path STU). To complete the illustration of the model,
details about the different branches referred to in the above overall
description are given below.
The description of 'monotonic envelope' is given in section 4.2 and by
Figs 61 and 62 respectively. 'Reduced envelopes' are obtained from the
monotonic envelope by decreasing the characteristic bond stresses T\ and r?
(see Figs 62 and 82) through reduction factors, which are formulated as a
function of one parameter, called herein the 'damage parameter a". For no
damage, d = 0, the reloading branch reaches the monotonic envelope. For
full damage, d= \, bond is completely destroyed (r = 0). The rationale for
this assumption is given by Fig. 83, which shows that reloading curves for
similar specimens subjected to different loading histories appear to form a
parametric family of curves.
The deterioration of the monotonic envelope seems to depend on the
16
- Monotonic loading
• After N = 1 cycle
12- Cycling between
±Smax
1 8-
90
BOND BETWEEN CONCRETE AND STEEL
91
RC ELEMENTS UNDER CYCLIC LOADING
10 1-5
(T3 in Fig. 82). The relationship found in the tests is shown in Fig. 85.
However, if cycling is done between fixed values of slip (e.g. between fixed
5 max and smin in Fig. 82), n is reduced more rapidly than the ultimate r 3 of
the corresponding reduced envelope. Therefore, the analytical function in
Fig. 84 is used only for the calculation of the frictional resistance for the
first slip reversal. For subsequent cycles, T{ (e.g. Tf in Fig. 82) is reduced
from this initial value by multiplying it with an additional reduction factor
which depends on the energy dissipated by friction alone. If unloading is
done from a larger slip value than the peak slip in the previous cycle (path
STU), the new frictional bond resistance (rfu) is linearly interpolated
between two values. The first value is related to r^ of the corresponding new
reduced envelope using the analytical function given in Fig. 85, and the
second value is the n reached in the last cycle (Tf in Fig. 82). This
interpolation is done in order to have a smooth transition in the values Tf.
Note that the concept of relating damage to one scalar quantity, like the
normalized dissipated energy, has provided a basis for a relatively easy
generalization of the bond behaviour for random excitations. The bond
model selected can easily be extended to cover bond of bars under
conditions different from those reported herein, such as different bar
diameter, pattern of deformation (lugs), concrete strength, degree of
confinement, effect of transverse pressure, etc. This requires that the
pertinent experimental data necessary for computing the different
parameters, in particular the monotonic envelope, be obtained. If these
are not available, the suggestions given by Eligehausen et al.9 could be
used for choosing the required parameters.
The bond conditions in a joint vary along the embedment length. For an
interior joint, three different regions have been identified (see Fig. 86).
They show differences both in the shape of the monotonic envelopes,
different for positive and negative slip, and in the rate at which degradation
occurs. Of course, there is a gradual variation in the behaviour proceeding
from an unconfined region to a confined one.
The possibility of extending the analytical model shown here for confined
concrete to the unconfined regions, using information presented by
Viwathanatepa et a/.113 has been detailed in Reference 92. The analytical
92
BOND BETWEEN CONCRETE AND STEEL
bond model can be generalized by two points. First, instead of only one, two
different monotonic envelopes are specified, one for positive and one for
negative slip values (compare Fig. 86). Secondly, the normalizing energy,
Eo, used in the computation of damage is chosen as the larger one between
EQ and EQ. These quantities define the corresponding areas under the
monotonic envelopes for positive and negative slip values up to slip value
s3. To take into account different rates of damage in the two directions of
loading, the pertinent dissipated energy, E, used for computing the reduced
envelopes, is multiplied by an amplification factor, b, which is different for
the upper and lower curve. The factors b+ and b~ are specified as input
values. Similar rules for the computation of damage apply to the friction
part of the curves.
More details regarding the quantification of the various parameters
involved, and of their distribution along the anchorage length in an interior
beam-column joint, are available.78'92 A comparison between experimental
and analytical results concerning a bar with a bond length of 25 bar
diameters is presented in Fig. 87.
4.4.2.7. Model by Pochanart and Harmon.115 The following rules
define the hysteresis rules for bond under generalized excitations (Fig.
88(c)).
(a) Under monotonic loading, the bond-slip relationship follows the
monotonic envelope.
300- i 18E
12F 15E
9D S - 21E
^-
JL—
A*
* 0- • " ^
1
e
(L-
-6K
21L
18L.
-300- ' 15L
0 -2
i: mm
(a)
Fig. 87. Comparison between (a) experimental and (b) analytical results; Ciampi et al. 78
93
RC ELEMENTS UNDER CYCLIC LOADING
10
Bond stress
Monotonic envelope
Reduced envelope Reduced
envelope
94
BOND BETWEEN CONCRETE AND STEEL
and by r{, the current value of the developed frictional stress. Also needed is
an offset value that tells where the strength component of the reduced
envelope begins. In Fig. 88(b), this offset value is set to one-quarter of ex.
T[ is composed of two components: the strength component and the
friction component. The strength component is obtained from Fig. 88(a)
given the total excursion. The friction component is determined from the
friction model shown in Fig. 88(a) by keeping track of the developed
friction lost during load cycling. The slip s\ corresponding to T[ can be
obtained easily since experimental results show that s{ minus the positive
excursion is a constant percentage of (, the undamaged length of the
concrete key which equals ls-ex. r 3 and rf are obtained by keeping track of
the developed friction and the loss of friction due to load cycling.
The deterioration of the frictional bond stress can be determined by
assuming that the frictional bond stress reduces by about 18% from the
frictional bond stress of the previous cycle, assuming no new frictional bond
resistance is developed during the cycle; that is, that the test is performed
under slip control. The deterioration of the bond strength is obtained given
the maximum total excursion and the strength model of Fig. 88(a). These
two rules define the damage law for bond.
The model appears to be able to predict a wide range of bond-slip
behaviour, including the results of both slip-controlled and load-controlled
tests from 100 to 40000 cycles (the limits were given by the authors). The
behaviour under cyclic loading can be predicted reasonably well by
repeated applications of a simple damage law applied to the reduced
envelope. The shifting of the ordinate axis of the bond-slip co-ordinate
system is a visual interpretation of the slip increase without an increase of
bond stress during the frictional bond behaviour; however, the value of
shifting and the shape of the ascending branch may need more specific
studies. The basic difference of the models proposed by Ciampi,
Eligehausen, Botero and Popov and by Pochanart and Harmon is the
definition of the bond damage during cycling.
95
RC ELEMENTS UNDER CYCLIC LOADING
AA. B E M F C
Fig. 89. Assumed bond stress distribution (a) and hysteretic bond-slip relationship (b); Filippou79'"4
end of the bar, point D, and point C which is located at the transition
between the confined and the unconfined concrete region.
(c) The bond stress distribution in the confined concrete region between
points B and C is assumed piecewise linear and is established
iteratively satisfying the equilibrium and compatibility conditions
using the bond model shown in Fig. 89(b).
80
Finite element model
Proposed model
60
-1
100 200 300 400 500 100 200 300 400 500
x: mm
(c)
Fig. 90. Distribution of (a) steel strain; (b) steel stress; (c) slip and (d) bond stress (in frame joint with increasing
magnitude of end slip); Filippou"4
96
BOND BETWEEN CONCRETE AND STEEL
4.5 Conclusions Cycles with reversed loading produce degradation of bond strength and
bond stiffness that is more severe than for the same number of load cycles
with unidirectional loading. Degradation primarily depends on the peak slip
in either direction reached previously. Other significant parameters are rib
pattern, concrete strength, confining effects, number of load cycles and peak
value of slip between which the bar is cyclically loaded. Under otherwise
constant conditions, the largest deterioration will occur for full reversals of
slip. The accumulation of bond damage is supposed to be caused by the
progressive growth of microcracks and concrete microcrushing in front of
the protruding lugs.
Although the models simplify the real behaviour, they intend to take into
account parameters which appear to control the behaviour observed in the
experiments, these being: deterioration of related slip, increase of slip under
load-controlled reversals, decrease of bond stress under slip-controlled
reversals, deterioration of frictional bond stress and deterioration of residual
bond strength. Using the described rather sophisticated bond models, the
behaviour of bars embedded in concrete under cyclic excitations can be
predicted with sufficient accuracy for practical purposes.
97
5. Interface behaviour
This Chapter presents a review of the models that can be used in finite
element analysis of RC structures, for the description of concrete-to-
concrete and steel-to-concrete interfaces. The attention is focused on
aggregate interlock and dowel action. Whenever possible, both cyclic and
monotonic behaviours are considered. Whenever necessary, test data are
recalled and commented on.
|B| Bending
. Crack
Vdu = 2MJL =
[¥] Shear
\K\ Kinking
98
INTERFACE BEHAVIOUR
Debonding
• • • • - . • ' • ' • •. }
A
Xbonding
••'.'•• i t l Friction
w
1
(a) (b)
Fig. 93. Through crack with: (a) aggregate debonding; (b) partially formed crack with contour lines of longitudinal
displacement in process zone; Cedolin et al."6
99
RC ELEMENTS UNDER CYCLIC LOADING
f,
• - /
1
1 [] (a)
^Bars
(b) (c) (d)
Fig. 94. Dowel action: (a) working load situation; (b) ultimate load situation; (c) dowel action and aggregate
interlock at collapse (vertical and hoop rebars); (d) dowel action in side-by-side elements
values of a few tenths of a millimetre, well beyond the values that are
generally examined with regard to the process zone at the tip of a
developing crack (a few thousandths or hundredths of a millimetre, Fig.
93(b)).
Before examining the general philosophy of the different approaches that
can be adopted in the formulation of the constitutive relationships, it is
necessary to look at the different effects that variable loads have on the
various interface mechanisms. Bond is very often required to withstand
either alternate forces (such as in beam-column joints subjected to seismic
loads) or pulsating forces (such as in a beam subjected to dead weight and to
live loads): consequently, bond degradation due to variable loads is of
paramount importance, with regard to strength and stiffness.
In monolithic structures (which may consist of different parts cast at
different times) dowel action is generally required to transmit pulsating
forces: this situation occurs both in the working load stage (Fig. 94(a)) and
at collapse (Fig. 94(b)). Alternate loads may occur in a few cases, such as in
RC secondary containment shells of nuclear reactors, subjected to seismic
loads and to internal pressurization (Fig. 94(c)): in these cases dowel action
and aggregate interlock are the only means for shear transfer across
horizontal and vertical cracks, if the shell is only orthogonally reinforced
(the former cracks run mostly along the horizontal construction joints, the
latter cracks are due to the hoop stresses). In side-by-side structures
connnected by dowels, such as in RC slabs resting on a deformable soil
(Fig. 94(d)), alternate forces are the rule, rather than the exception.
Since aggregate interlock and dowel action are very often activated by the
same forces and are strictly related to the same crack pattern, the
considerations made so far for dowel action also hold for aggregate
interlock. One must, however, bear in mind that for small crack widths
aggregate interlock dominates over dowel action, while at increasing
interface displacements (as is generally the case with load cycling) the role
of dowel action becomes prevalent, also because interface deterioration
builds up more quickly than concrete deterioration under the dowels.
Finally, the interaction between the two shear transfer mechanisms is even
more complex because of the variable restraint stiffness provided by the
reinforcement crossing the crack (bond-induced tension stiffening
deteriorates under load cycles).
As to aggregate debonding, which is related to highly deteriorated
concrete in tension, only monotonic loads are relevant, because variables
100
INTERFACE BEHAVIOUR
Fig. 95. Tests for the study of aggregate interlock: (a) at constant confinement stiffness; (b) at variable confinement
stiffness; (c) at constant crack opening (infinite confinement stiffness); (d) at constant confinement action (zero
confinement stiffness); (e) at constant crack dilatancy
Onset of cracking
• • Test results
Gradual
removal of
Restraining bars restraining
Fig. 96. Example of direct bars
tension test for the study of Fa - ^ ^ F°.
concrete behaviour in
tension, with aggregate
debonding; Giuriani and
Rosati134
101
RC ELEMENTS UNDER CYCLIC LOADING
(a)
P/2 P/2
5.2. Aggregate The shear transfer mechanism based on aggregate interlock has been known
interlock for a long time in its behavioural aspects, owing to the many test results
obtained in the 1960s, 1970s and early 1980s. The different test results shed
light on the many parameters involved, either kinematic (crack opening and
slip, crack dilatancy, initial crack width), mechanical (confinement across
the crack plane, aggregate and matrix strengths) or morphological (overall
and local roughnesses, maximum aggregate size and aggregate type, i.e.
rounded, crushed, light-weight). Unfortunately, scant attention has been
devoted so far to path dependency, either in the displacement field or in the
stress field, to cyclic behaviour and to the three-dimensional aspects of a
crack, while the statistical aspects of aggregate distribution along a crack
have been explored by very few researchers.
Remarkable experimental effort made it possible to understand such basic
102
INTERFACE BEHAVIOUR
concepts as crack dilatancy (i.e. the coupling between shear stress and crack
opening, which also involves the coupling between normal stress and crack
slip) and shear-confinement interaction, but no comparable effort was
devoted to the formulation of rational constitutive laws until recently,
becuase only recently have the properties of aggregate interlock been
recognized as material properties. In fact, the formulation of suitable
constitutive laws for cracked concrete became necessary after the concept of
smearing the cracks over an entire element was shown to be highly suitable
for finite element analysis.
So far, the greatest attention has been devoted to planar cracks subjected
to monotonic loading, in order to formulate the incremental stiffness matrix
of a crack, which may be generally assumed in the form (see Fig. 98 for the
symbols)
Bat\ fd8n\
where Bnn, Bnt, Btn, Btt are crack stiffness coefficients, which depend on
8n, 8t, o^n, a^t and possibly also other state parameters. The stress-
displacement relations in equation (113) are analogous to the stress-strain
relations of incremental plasticity. Note that crack stiffness matrix is neither
symmetric, nor positive definite, so that crack response tends to be unstable.
However, the response is usually stabilized by the restraint provided by the
reinforcement.
In order to formulate the stiffness coefficients of matrix B, test results at
constant crack opening (Fig. 95(c)) and at constant crack slip should be
available; since tests at constant crack slip have scant inherent meaning,
tests at constant confinement (Fig. 95(d)) have been performed, but these
tests are far less numerous than the tests of constant crack opening and at
constant (or variable) confinement stiffness, the latter being of very limited
use in the formulation of matrix B, since all parameters (cr^, o£ t , 8n, 8t)
evolve during each test.
Due to the relative scarcity of test data, simpler and less general
formulations have been adopted so far, according to the total deformation
theory. Generally the functional relationships have been expressed in a
direct stress-displacement form or in a mixed form
*) (114)
n, fit) (115a)
fin) fit = A(<4, «„) (115b)
Of course the weakness of a total deformation approach is the path-
independency of the response, while path-dependency must normally be
expected of inelastic behaviour.
In finite element analysis, once an element becomes cracked, it is
convenient to adopt a 'smeared crack' approach in order to introduce the
properties of cracks, aggregate interlock included. By replacing the actual
103
RC ELEMENTS UNDER CYCLIC LOADING
Bars
py sin 0
(a) (b)
Fig. 99. Cracked reinforced concrete Fig. 100. Equivalent steel ratio pn
(117)
where c is the concrete cover, <f> is the bar diameter and p is the geometrical
steel ratio; k\, k2 and ^3 are constants, related to bond efficiency, type of
deformed bars and type of loading.
Equation (117) is valid only for one array of parallel bars at right angles
to the crack pattern, which is seldom the case, since generally two arrays of
bars are present, neither of them at right angles to the cracks (Fig. 100(a)).
In these cases an 'equivalent' steel ratio, pn, (Fig. 100(b)) must be worked
out and introduced into equation (117). With regard to this, a criterion for
establishing the equivalence between pn and (px, py) should be introduced,
but so far nothing is available, except in the limit case with completely
yielded bars; in this case a simple equilibrium condition leads to the
expression for pn
2
9 + pyfsySm20 i.e. pn = pxcos20 + pysin20 (118)
where 6 is the angular deviation between the x-bars and the axis at right
angles to the cracks.
In order to work out the stiffness matrix, Dc, of the concrete as a whole
(cracked concrete + solid concrete between the cracks, called in the
104
INTERFACE BEHAVIOUR
(120)
Since the stresses and the strains in the solid concrete are usually small, a
of
1 — Vc 0
F sc
= 1 1 0 (121)
-Vc
3(1 - 2 v c
0 0 3(1 - 2vc)Kc I2G
where Kc = Ec/[3(\ - 2vc)].
Should the simplification of linear elasticity be dropped, one could
determine the flexibility matrix F s c by differentiating some total stress-
strain relations for concrete, or just assume incremental isotropy and
evaluate the tangent shear modulus, GCT, and the bulk modulus (or modulus
of volume expansion), KCT, as functions of octahedral normal and shear
strains, using some of the models presented in Chapter 1.
Once the tangent stiffness of plain concrete has been evaluated, the
influence of the reinforcement can be introduced, by adding the tangent
stiffness matrix of the reinforcement, Ds, to Dc. Because of steel-to-concrete
compatibility at the macrolevel, the average strains of the reinforcment are
assumed to be equal to those of plain concrete; consequently, the stress-
strain relations for the reinforcement in the crack reference system (n,t)
(Fig. 98) can be written as follows
dcf = D s de where de = de8 = dec (122)
Note that the effects of tension stiffening are taken care of in the stress-
strain relation of the steel of each array, and ultimately are implied in the
formulation of the stiffness coefficients Df-.
Since crack orientation does not generally coincide with the orientation of
either of the bar arrays, the matrix Ds must be obtained by transformation
P2 Q2 PQ
D*R,, R, = Q2 P2 -PQ
i=l,2
-2PQ 2PQ P2-Q2_
""• B / Y f,
0 0
0 0 0 (123)
0 0 0
with Q = sinl?,, P = cos0(, / = 1 for x-bars (0| = 9, Fig. 100(a)) and / = 2
for y-bars (02 = \ir + 0).
OABYR: with tension stiffening
OYR: without tension stiffening
In D^, px and py are the geometrical steel ratios of the two arrays, and
Fig. 101. Stress-strain law of Es(ex), £s(cy) are the tangential stiffnesses of the steel, tension stiffening
steel with and without tension included (Fig. 101). Note that Es(ej) is never smaller than the elastic
stiffening Young's modulus of the steel, and can even be a few times larger.
105
RC ELEMENTS UNDER CYCLIC LOADING
w Upper
" block
106
INTERFACE BEHAVIOUR
(a) (b)
of either material internal friction and/or local interface roughness, the latter
being subjected to degradation due to possible cyclic loading. The shape of
the general roughness is introduced by means of suitable, prefixed functions
of the coordinate x, such as arcs of second order parabolae.
In its original form, the general roughness model is only a 'functional'
model and not a 'predictive' model, since the many geometrical parameters
involved are a priori unknown and a posteriori unmeasurable. These
parameters may be determined by means of a multiple regression analysis,
starting from a sufficiently large data base. As regards qualitative predic-
tions, the model seems to work well (Fig. 103)135'136 and to be suitable to
cyclic loading.
With reference to the shear stiffness of a reinforced concrete block
cracked in the direction x (Fig. 103(a)), the following expression can be
derived for the relationship between the shear stress and the strains
1 C\
7yx (125)
1 1
+ _
UC Cty
where: Gc = elastic shear modulus of the concrete; c\ - constant; ay, (3y =
functions of the extensional and dowel stiffnesses of the cracks, and of crack
spacing. Expressions for G011 (shear secant modulus of cracked concrete)
are given by Leombruni et a/.136 Equations (125) was worked out by means
of a multiple regression analysis performed on test data by Gergely,
Jimenez-Perez, Laible and White at Cornell University.137'138
l39 142
5.2.2. Rough crack model ~
The rough crack model is based on some general properties that are to be
expected for rough cracked surfaces, and that can be deduced by
considering a few simple micromechanical models (Fig. 104(a) and (b)):
the interface stresses are assumed to be mostly dependent on the
displacement ratio r = 6t/6n (wedge effect); for small crack widths, the
confinement stress vanishes (aggregate debonding); for large crack widths,
the stresses tend to vanish due to concrete microcracking and crushing, and
eventually to the loss of contact between the crack faces. The general
roughness is disregarded and only the local roughness is introduced;
consequently, since the number of contact points along the crack interface
may be regarded as infinite, the resulting stress-displacement relations may
be considered continuous and smooth.
The interface stresses are formulated in the following way
<t = Fl(6n/dtk)F2(r) (126a)
107
RC ELEMENTS UNDER CYCLIC LOADING
+a2\r\
<t
x0-25
(126b)
where r u = r o [l - \/(28n/d&)].
The rough crack model can also be extended to cyclic loads, provided
that suitable stress displacement relationships are formulated both for
reloading and unloading.
As regards pulsating loads, the tests at constant crack opening by Paulay
and Loeber121 show that the shear stiffness of a crack is actually greater in
the following cycles than in the first cycle, because of some contact between
the general roughness of the crack faces (Fig. 104(a)), after the deterioration
of the local roughness in the first cycle. Though different from the first
cycle, the following cycles are practically coincident (at least up to / = 20, i
being the number of load repetitions).
Since test data are really limited and even non-existent for confinement, a
suitable and reliable constitutive law can hardly be formulated. A tentative
formulation for pulsating loads is presented by Gambarova;140 for loading
and reloading the constitutive laws are very similar to equation (126), while
completely different formulations are presented for unloading. Typical
hysteretic cycles under pulsating loads are shown in Fig. 106.
Aggregate particles
Crushed mortar
Fig. 104. Crack morphology
(F = local interface forces):
(a) wedge effect; (b)
aggregate debonding; (c)
microc racking and concrete
crushing
108
INTERFACE BEHAVIOUR
50\
Equation (127); f'c = 31 MPa,
da (maximum aggregate size)
i 2
• - —
= 19 mm n
0-5 10 1-5 20
r = d,ldn
(b)
ru = 7-83
17
A <°Sr>" = 5-80
An = 0-125 mm
I -'A • I J\ dn = constant
I / /
/
(d,)Ud,y,
Shear displacement d,
(a)
121
Fig. 106. (a) Typical hysteretic cycles at constant crack opening under pulsating loads (b) fit of test results
with the loading and unloading equation presented by Gambarova'40; f'c = 31 MPa, da (maximum aggregate size)
= 19 mm
109
RC ELEMENTS UNDER CYCLIC LOADING
(c)
<„ = - (128)
in which A, and An are very complex functions of the relative displacements
6n, 6t), maximum aggregate size (da) and relative aggregate volume (pk =
ratio between the total volume of the aggregates and the total volume of the
concrete: in most cases, the value 0 • 75 was adopted by the authors)" 8 " 119
Ax and An have an integral formulation, and the integrations have to be
performed numerically, with respect to the aggregate diameter which is
considered as a variable related to the statistics of the contact points at the
interface and to aggregate grading. A linearization of equation (128) is
proposed by Walraven et al. The agreement between the values predicted by
the two-phase model and the test results is fairly good (Fig. 109).
The extension of the two-phase model to cyclic loads (either pulsating or
alternate) and to sustained loads has been recently attempted144 taking
advantage of the inherent capacity of the model to follow and to quantify
the damage process (Fig. 108): for any couple of values 6n, St), both the
contact areas and the contact forces can be directly evaluated (Fig. 108(a)).
The fit of some of Laible's test results is presented in Fig. 109:138 a
(a) reduction of the friction coefficient and of the contact areas had to be
introduced because load recycling makes crack interfaces smoother,
particularly in concretes with poor aggregate strength.
r= IV = 1 N = 15
\
E
10:
//
/'
| h 10 10-
05: / 05- E 0 5-
-10
/,'J--
-0-5
/- "b-5 ' ' Vo 1 5 10 15
f-
• / / ' -0-5
d
' -0-5- - -0-5-
N
f/ \ -
(b)
-10 - Test
-10- -10-
/ -- - Calculation
Fig. 108. (a) Definition of
new points along the crack
interface; (b) various stages Fig. 109. Fit ofLaible et al. 's test results (RC specimens with constant confinement
of the crack damaging stiffness)138; f'c = 26MPa, da = 32 mm, /x = 0-20, apu = 30MPa, £° = 0-75mm
process (initial crack opening)
110
INTERFACE BEHAVIOUR
where Ko = initial shear stiffness, function of 6° (initial crack width) and / c '
(concrete compressive strength); b = empirical parameter, function of Ko, f[
and CT^JJ (confinement stress).
As for the second constitutive law, 6n = W(8t, o^n), the formulation is
still empirical and based on the authors' test data
ill
RC ELEMENTS UNDER CYCLIC LOADING
dA dF
dSt = (135)
dF
kn (crack normal stiffness) = ——
-i
= -rr^- (dilatancy ratios)
d6Q
j[dot
daC
nt\ A
\daj or
1 -
WJ ~ 1 d6n
(136)
112
INTERFACE BEHAVIOUR
o Numerical [139]
• Tests [121, 138, 146]
Fig. 111. Incremental slip-confinement model:146 (a) shear stress; (b) shear stiffness plotted against crack slip
(shear displacement); (c) frictional coefficient plotted against crack opening (normal displacement)
= 39 MPa 10 y.
= 0-3 mm
<>„
Al
a Jw It
I 5
\d -0-4 0
' '/ ' 0-2 04
A,', mm
- 5 - — Experiment
— Analysis
-10-
(a) (b)
1 1 1
fc = 23 MPa bn = 1-Ov^
Analysis 4 ' - 4
Fig. 112. Contact density
model — cyclic shear
- - unloading
— loading j t A0
0
1
1
- 2
transfer at constant crack 1
1 d
113
RC ELEMENTS UNDER CYCLIC LOADING
tT 0 for
nt=°'nn = © < 0, (no local contact) (142)
where At is the total area corresponding to the nominal unit area of crack
surface (A, is considered as an interface — or material — constant) and k is
a reduction factor which introduces the effects of crack displacements on
the contact area; a and k are functions of crack displacements 6n and 6t and
of contact orientation 8 ; (At p(Q)) is the effective contact area along the
planes whose orientation falls in the interval 6 ± dQ/2.
Assuming that friction is negligible ( 6 is independent of crack dis-
placements), that crack displacements are very small with respect to crack
roughness (k = 1), adopting a reasonable formulation for />(= 0-5 cos 8 )
and for/l t (= 20 500f^'^/E c ), giving a a simple formulation based on crack
displacements via geometrical considerations, equation (141) can be worked
out, as well as the crack stiffness matrix B (equation (113)), by derivation of
the same equation (141). By integration, according to any prefixed
displacement path, the crack response curves (a^v o£n) can be obtained.
In this way, several test results (either at constant or at variable crack width)
were successfully fitted, under monotonic as well as cyclic loads (pulsating
or alternate loads, Fig. 112(b) and (c)).
The proposed model is very attractive for its ingenuity and simplicity;
some basic ideas go back to the microplane model148 (contact planes,
contact stresses) and the two-phase model (plastic behaviour in the contact
zones at the interface).
114
INTERFACE BEHAVIOUR
Fig. 113. Generalized rough crack model — idealizations of crack morphology:149 (a) castellated crack; (b) rigid
saw-tooth crack; (c) deformable saw-tooth crack
cracks, the generalized rough crack model has been limited so far to plane
problems, though extensions to three-dimensional behaviour is possible.
Since this is based on crack smearing, the transformation from crack
displacements (which appear in the constitutive laws of the models) to
'equivalent' strains requires the introduction of a 'reference' length, le,
which is closely related to crack spacing.
Constitutive crack models (based on constitutive relationships), as well as
physical models (based on direct modellization of a rough interface (Fig.
113)) can be introduced into the generalized rough crack model: so far, the
so-called deformable tooth model (Fig. 113(c)) has been adopted, while the
solid concrete has been modelled by means of linear elasticity or by the
endochronic theory. The stiffness approach was used by Riggs and Powell
to calculate the matrix D. 149
In the deformable tooth model the initial shape of the crack is a saw-tooth
and under loading a quasi-static sliding of the two surfaces occurs, with
friction and surface deformations (one face of the crack is assumed to be
deformable, while the other is assumed to be stiff). Since the model is
characterized by nine constants, which depend on the mechanical and
geometrical parameters, and the determination of the crack state is complex
in itself, the implementation of the model in a pre-existing finite element
program proved to be rather difficult. However, Paulay and Loeber's test
results are well fitted, while the fitting of cracked concrete cyclic behaviour,
although promising, requires further improvements (Fig. 114).
Slip: mm Slip: mm
0-5 10 1-5 0 05 10 1-5
1 i i
- 20
S. — Experiment
— Analysis 15 |
_
1 A
i i
\A i
10 2
^5
1
/
Fig. 114. Cracked concrete
II
/ hy i
i - 0-5
cyclic response: (a)
analytical results for 15
it r
002 004 006 002 004 006
cycles; (b) fit of Laible et Slip: in Slip: in
al. 's test results'38 (a) (b)
115
RC ELEMENTS UNDER CYCLIC LOADING
Microplanes
Paulay and Loeber's tests
—- 75% of test results
0-4 0-8
<V mm
(c) (d)
subjected to triaxial or biaxial stresses, 150151 the microplane model has been
applied also to describe shear transfer across blunt cracks, by modelling the
cracks as 'crack bands', with a width equal to maximum aggregate size.148
Within a crack band, the microplanes stand for the thin mortar layers
(Fig. 115(a)) which transmit either compressive stresses (contact planes)
between two contiguous aggregate particles, or tensile stresses (planes with
strain softening); Bazant and Gambarova give the microplane constitutive
law a hypoelastic formulation, with different branches for loading,
unloading and reloading in compression and tension148
a = Fie )e (143)
The tangent stiffness matrix of the microplane system is obtained by
means of an integration performed over a unit-sphere (spatial cracks) or
over a unit-circle (planar cracks), representing all possible microplane
orientations in the space or at right angles to the reference plane (Fig.
115(b)). In the latter case
2
(144)
with aijkm = nitijtikUm, where n are the director cosines of each microplane
(1,2 = references axes, x and y in Fig. 115(b)).
Simple asymptotic or exponential laws have been adopted for the
microplane constitutive relationship (Xn(en), both for loading and unloading
(Fig. 115(c)). The model can describe the path-dependency of the concrete
(even if the constitutive relationship is of the total stress/total strain type)
and the reorientation of the principal stresses with respect to principal
strains.
116
INTERFACE BEHAVIOUR
T P : Two-phase model
(Walravenefa/. 118120 )
Most available test results have been successfully fitted (see for instance
Fig. 115(d)), and the extension to cyclic loads — though not yet attempted
— may be easily achieved, because of the inherent simplicity of the model,
which is based on a single, many-branched relationship, valid for each
microplane (equation (143)).
Recently the shear stiffness of the microplanes has been introduced for a
better description of concrete non-linear behaviour, and the model has been
extended to cyclic loading (cyclic compression).152
Three of the aforementioned models have very recently153 been intro-
duced into the DIANA computer code and a very comprehensive comparison
has been carried out in order to assess the performance of the models with
regard to the description of a few tests performed at constant crack opening
and at constant restraining stiffness, and to the numerical stability once the
models are implemented in a finite element program.
The study regards the rough crack model in two versions,139'141 the two-
phase model and its linearization,118119 and the contact density model.147
The three models perform fairly well, on the whole, but the contact density
model seems to have the edge, particularly from the numerical point of
view. An example of the thorough analysis performed by the Dutch
researchers is shown in Fig. 116.
117
RC ELEMENTS UNDER CYCLIC LOADING
such relevant aspects as crack roughness and aggregate shape and size
(which are principally responsible for crack dilatancy), concrete strength,
aggregate-to-mortar friction are ignored.
Furthermore, there is general agreement among researchers that in many
structural problems the value of the shear-retention factor (either constant or
variable) is not critical, provided that a sufficiently large value is adopted
(a > 0 • 1). Otherwise the system of equations in finite element analysis
becomes ill-conditioned and the iterative solution has to be stopped because
of numerical instability (divergence). Whether divergence is just a numer-
ical problem or has a physical meaning is still an important question to be
answered, as observed by Rots et al.64
In the earliest period (1970-1975) the factor a was generally given a
constant value (between 0, smooth cracks, and 1, fully locked-up cracks
(Fig. 117(a)); 0 - 4 - 0 - 5 in most cases). Later on, a two-step or a
continuously decreasing formulation was adopted (see for instance Cedolin
and Dei Poli154), after crack opening was recognized to be a major factor in
reducing crack shear-stiffness (Fig. 117(d)).
In most formulations the crack opening is replaced by the equivalent
strain e ^ ( = Sn/s, where s is the crack spacing) at right angles to the cracks.
Very often the solid concrete contribution to normal strain is neglected:
e m = e^ 1 . For the literature on the shear-retention factor see Walraven and
Reinhardt,119 Rots et al.,64 Walraven and Keuser.155
Among the most recent formulations,155'156 can be quoted
-1 ir
7fr,n = 0 • 002(n -
118
INTERFACE BEHAVIOUR
where
= cycle number
<$t = factional response and crack slip at the end of the loading
branch of the nth cycle
= frictional response in the first cycle, for a crack slip equal to
= confinement stress acting on crack interface
6tn = crack slip corresponding to the maximum mobilized shear
stress under monotonic actions.
By definition the secant shear moduli are
s/ \s
and the following expression can be worked out for an/a\
^ = 1 - (145e)
QCR
1
For n = 5-15, aclf'c = 0-25-0-5, St/6tu =0-5-1-5 the values of an/a\ fall
within the range 0-4-0-8, according to equation (145e)
-cl Esp\
(1-A.) +•
2(l+v 2 1 ) 2(1+v 8 )
+ (1- + • (146)
Cracked
sections
2(l+v 1 2 ) 2(1+v 8 )
where Ai, A2 (bond-slip moduli) are the ratios of the debonded lengths of
Fig. 118. Orientation of the rebars to crack spacing (Ai and A2 depend on crack width); Eci and Ec2
orthotropic axes 1, 2, 3, are the tangent moduli of concrete; p\ and p2 are the 'resolved' steel ratios;
defined by initial
cracking;1'8 x, y, z global the reinforcement is suppoed to be elastic. No comparison with test data is
coordinate system presented. No extension to cyclic loading is performed. Aggregate interlock
119
RC ELEMENTS UNDER CYCLIC LOADING
id
— —
—
1 2 3
Applied shear stress: MPa
Rebars
(b)
and dowel action are neglected. In most practical cases the value of a varies
between 0-75 and 1.
5.2.10.2. Duchon's model.159 Reinforced concrete elements subjected
to membrane forces are modelled as a system of vertical, horizontal and
two-way inclined bars, plus a system of concrete struts, after concrete
cracking (Fig. 119). By writing equilibrium and compatibility equations
along the cracks (where the internal forces are absorbed entirely by the
reinforcement), the stresses and the strains in the steel and concrete, as well
as the moduli of the resultant material, can be worked out. No closed-form
expression for G01* is presented, but an expression similar to the one given
by Collins160 can be derived. Aggregate interlock and dowel action are
neglected. No extension to cyclic loading is performed. The values
predicted for G 0 * are very low (CfR/Gc = 0-1-0-2).
5.2.10.3. Collins' model.160 Reinforced concrete and plain concrete
planar elements subjected to prevailing shear are modelled by means of
diagonal compression fields (Fig. 120), after concrete cracking. Considering
the equilibrium requirements (both in the longitudinal and transverse
direction) and the compatibility requirements (the directions of the principal
tensile strain and stress are assumed to be coincident) the behaviour of a
member in shear is analysed, after the average stress-average strain
relationships for both the steel (bars and stirrups) and the concrete struts are
given suitable formulations. The angle a^ of the diagonal compression field
turns out to be a function of the steel ratios in the horizontal and vertical
(a)
Strut
(d)
(b)
Steel
o, = 0, a2 * 0
Fig. 120. ad, ad: diagonal Fig. 121. Saw-tooth idealization of single crack: (a) crack opening; (b)
compression field; 1, 2: principal crack opening and slip (aggregate interlock and dowel action); (c)
directions in cracked concrete, at behaviour of diagonal concrete struts; (d) compressive stress-strain
onset of cracking relationship for concrete in strut
120
INTERFACE BEHAVIOUR
directions. As a result, the stresses and the strains, as well as the effective
shear modulus are evaluated prior to yield of the stirrups, while the ultimate
shear stress is evaluated after yield of the stirrups. Dowel action is
neglected.
The effective shear modulus, CFR, is as follows
s h e a r
? ) (147)
where n = EJEC is the modular ratio and pt, p\ are the steel ratios in the
transverse and longitudinal directions. No application to cyclic loading has
been attempted so far.
5.2.10.4. Perdikaris and White's model}61 A saw-tooth type
idealization is adopted for the cracks (Fig. 121 (a)), which are assumed to
run parallel to the two-way reinforcement, if the tensile stresses prevail
(Tky 5: 0-7 MPa with x,y = reinforcement directions). Furthermore, diagonal
cracking is assumed to be the prevalent cracking mode for sufficiently large
shear stresses (T^ > 0-7 MPa (Fig. 121(b)). Based on the local equilibrium
and compatibility requirements, and on largely empirical formulations for
the shear stiffnesses, KlSr (interlock shear transfer) and KDA (dowel action),
three expressions are worked out for CJ~R
1 I
1 + F, + F2
with F, = (SnEc + 1)-', F 2 = -SnEcvc{SnEc + 1)"' (149)
where Ec, Gc, vc are the elastic moduli, and 5 n is a function of the 'tension
stiffening factors' Ax and Ay:
K
S - '
On = — "
121
RC ELEMENTS UNDER CYCLIC LOADING
5.3. Dowel action Even more than aggregate interlock, dowel action has been known for a
long time as an effective shear-transfer mechanism, which is present in
different types of connections and whose importance goes beyond RC and
PC structures. Stud and pin connectors, anchor bolts and reinforcing bars
develop dowelling forces in the surrounding concrete, which produce highly
localized compressive stresses accompanied by less localized tensile
stresses. The latter, combined with the wedging action of the surface
deformations of the steel, produce splitting cracks, particularly if concrete
cover, bar free interspacing and transverse reinforcement are not properly
designed.
Contrary to aggregate interlock which is essentially a material property
(crack roughness depends on aggregate type and shape, on cement paste
strength and on mutual adhesion), dowel action is essentially a structural
property, because its effectiveness is strongly related to the detailing of the
reinforcement, to the shape of the section and even to the loads or to the
constraints, since transverse pressure can markedly reduce or increase
dowel stiffness and strength.
Because of the many parameters involved, the majority of the tests
performed so far regard the ultimate capacity of dowels, for different values
of the geometrical parameters of the section and reinforcement, while
several basic aspects such as cyclic loading, fatigue, concrete deterioration
and also mathematical modelling have received a great deal less attention.
On the other hand, in RC members dowel action becomes active when
cracking is of major importance, but this stage is reached often in the ulti-
mate load situation, when the ultimate capacity and not the other afore-
mentioned aspects are relevant. Moreover, dowel action is a very localized
phenomenon in itself, unless a regular one-way or two-way reinforcement is
accompanied by a regular system of cracks, which is seldom the case,
except in planar elements, with regularly distributed loads and no diffusive
processes.
With reference to RC and PC members, dowel action has several specific
aspects. The local bending and shear of the bars is always accompanied by a
complex triaxial state of stress in the concrete. Depending on the geometry
of the concrete mass interacting with the bar, tensile stresses (Fig. 123(a)) or
compressive-shear stresses (Fig. 123(c)) prevail. In both cases, the dowel
122
INTERFACE BEHAVIOUR
W,, A
•T
123
RC ELEMENTS UNDER CYCLIC LOADING
Dowel
capacity. These equations trace back very often to relatively simple, but
effective limit anlaysis models. 127164l65(b) Finally, the long-established
rational model by Johnston and Zia166 is still the most comprehensive
attempt to describe the different resistant mechanisms, which are activated
as soon as the concrete cover cracks and the stirrups (if any) come into
action. An extensive review of the literature on dowel action pushing
against concrete cover is presented in a paper by Dei Poli et al. together
with new test data.130
124
INTERFACE BEHAVIOUR
Fig. 126. Regularly cracked and reinforced field with one system of cracks at an angle to the y-bars: (a) forces
developed in the bars; (b), (c) forces and displacements in the crack reference system; (d), (e) effects of crack
skewness on dowel action (x-bars): (d) negative and (e) positive (Nx, Ny, V^,, Vyx = normal and shear forces per unit
length of a crack)
stiffening and dowel action are present: the former is introduced through
E*(exR) ^ Es and the latter through the foundation modulus k(T^R).
In equation (150), n is the elastic modular ratio (=EJEC) and a is the
parameter governing the elastic foundation theory
kdh \
a =
4EsJhJ
where Jb and db are the moment of inertia and the diameter of the section of
the bars.
Whenever the cracks have another direction with respect to the bars (Fig.
126(a)), not only the tensile and shear forces developed by the bars would
have two resolved components in the crack reference system (Fig. 126(b)),
but for each bar array the foundation modulus k would be different,
depending on the skewness of the crack plane (tension stiffening in the bars
is little affected, or totally unaffected by the orientation of the crack plane).
For instance, for 6t > 0 (Fig. 126(c)) the x-bars push against a limited cover
due to the skewness of the crack (Fig. 124(d)), while for 6t > 0 (Fig. 126(c))
the same bars push against an unlimited cover (Fig. 126(e)) and the
foundation modulus is larger than in the previous case (see also section 3.2).
In equation (150), the modulus El(e^R) should be given one of the
'tension-stiffening' formulations available in the literature.16? As for dowel
action, the non-linear behaviours of the concrete and of the reinforcement
are lumped into £ ( 7 ^ ) , although the two non-linearities should in principle
be dealt with separately. The equivalent strain 7 ^ is the smeared
component of the crack displacements, at right angles to the x-bars (Fig.
126)
(5tcos0)
"xy
VcosflJ
125
RC ELEMENTS UNDER CYCLIC LOADING
A\\ mm
(b)
crushed zone underneath the bar and close to the shear plane. Once the
concrete flake gets detached, the non-linearity of the embedment becomes
even more pronounced, since a non-negligible lever arm accompanies the
dowel force (Fig. 127(a) and (e)* = 0-4-0-6 db). Typical force-displacement
curves are shown in Fig. 127(b):129 as a rule, the response curves are
bilinear, with an elasto-plastic behaviour, but a limited softening for small
diameter bars (db= 14 mm) and a limited hardening for large diameter bars
(db = 24 mm) is often found in the tests. With regard to this point, high
strength concretes with basaltic aggregates exhibit a definitely elasto-
softening behaviour </c' > 80 MPa).131
In order to describe dowel action at increasing loads (monotonic loading)
the model of the beam resting on a cohesionless foundation is still
acceptable, but different formulations should be adopted depending on the
purpose of the modellization. Should the purpose be the description of the
load-displacement curve limited to the loaded section of the dowel (i.e. in
the shear plane), a relatively simply unidimensional formulation for k would
be sufficient; even a constant value, if the load level is sufficiently low with
respect to the ultimate capacity (service loads).
Five different formulations are cited below and the corresponding force-
displacement curves are compared with some of the results shown by Dei
Poli et al. (Fig. 128).129 The equations are presented for only two
formulations, since in the other cases the details are easily found in the
literature referenced.
5.3.2.1. Walraven and Reinhardt.119 The force-displacement
response traces back to the theory of the beam resting on an elastic
foundation (see also Section 5.3.2.4): for the foundation modulus an
expression already worked out by Paulay and Loeber is adopted (k = 188
JVC0'85), with some adjustments in order to introduce the effects of concrete
stiffness in a more consistent way and to improve the fitting of the authors'
own test results on dowel action.
126
INTERFACE BEHAVIOUR
V
• ^ . 0-2 0-4 0-6 0-8 10
Gambarova; B - Brenna
et al. db = 24 mm
^ ^ -
\\B>
40 40 -
Ibffx30 v ^>•"(3 < S_
fa s = 18 mm
20
Bv 10 w db ~ 14 mm
20 -
results
W 0 i i 1
i i
1 2 1 2
Ay. mm / ] , : mm
= — Wi with /3 = , a = (151)
4EsJh
where keq is the so-called 'equivalent foundation modulus', which is a non-
linear elastic modulus instrumental in the description of the force-
displacement response of the loaded end-section of the dowel, but unable
to describe the behaviour of the dowel in the other sections. The formulation
proposed by Di Prisco and Gambarova131 is based on the ratio V/Vu (which
is considered here as a 'damage index'), where Vu can be evaluated by
means of the available equations (see for instance Dulacksa126); as for the
roles of f and db, the proposal by Soroushian et al. is adopted165(a)
with ks = 127VU0AC6
V
for — < 04 7 = 2-12
V
for — > 04 7 = /0-544 + 0-026 cos h 8 — -
Ml
127
RC ELEMENTS UNDER CYCLIC LOADING
-Dulacska
-Soroushian
f'c = 30 MPa ':'••': : ' • • ' • • ' ' • • )
. -•:•:•:•.••••• Walraven — Reinhardt (pa = 001)
200 - • Walraven - Reinhardt (eM : 002)
1 Tests A and B
••>.•' :•••'.'-.•
I Walraven — Reinhardt (oa 003)
Fig. 129. Ultimate dowel N
capacity with no tensile • Slanted shear plane 6 = 45° J^
160 - N
force (a) and with tensile
force (b): (a) fits of test
120 "
data by Dei Poli et al.;'29
(b) failure envelopes:
V*u = ultimate capacity 80 -
according to Dulacska ; Soroushian
Vu(Soroushianetal.163) **
40 "
= 0-82 V*u; V (Walraven Dulacska f' = 30 MPa fsv = 440 MPa
and Reinhardt'w)= 0-72
ii i i
V* • N * = A u f •
12 16 20 24 28 32 10
psi = steel ratio of the db'. mm
longitudinal reinforcement (a)
128
INTERFACE BEHAVIOUR
Specimen A5 - d b = 18 mm
fc= 31-2MPa
50 I Cycles
Ay > 0
40
30
20
10
Fig. 130 (right). Dowel response under reversed transverse
displacement (Vintzeleou and Tassios164): empirical model
Fig. 131 (far right). Typical hysteresis loops for pulsating loads
(Dei Poli et al.v). dh = 18 mm, fc = 31 -2 MPa
i /29i
Section B
(t f . . . • • • • . . : . •
' < • • •
* > •
51)
\
(a)
Section A
*.
(b)
129
RC ELEMENTS UNDER CYCLIC LOADING
plane; fc - 24MPa,fsy =
440 MPa, s - distance 40
: . : : • ; • .
11 / O'b = 18 mm
-" D7
24
16 D1
— —
• —
D3
8 \// db= 14 mm
n
f Cover
f H
2 3
A,: mm
(a)
130
INTERFACE BEHAVIOUR
40 - / s =db
° 40
>
J 30 30 -
= 3d b
S
20 20
10 10 - S 50-85%
I 1
12 16 20 24 28 32 12 16 20 24 28 32
da: mm d^. mm
(a) (b)
0-83 10- 4
y — (153)
W
• Taylor (see reference 130)
Vn = 9078 + 0-1 blU. (154)
Vu = 37 *„{/(£) (156)
V u = 1-64 (157)
• Jimenez-Perez er a/.170
3-79c
Vu = M d 3-25 + (158)
131
RC ELEMENTS UNDER CYCLIC LOADING
/>•
Tests D IO Concrete cover c = 2d0
(refs 130, 131) I • Concrete cover c = da
• Krefeld and Thurston
(•*)-
60 • Taylor
• Soroushian er al.
A Houde and Mirza
• Baumann and Rusch
50 Jimenez-Perez ef al.
(*)
40
i
30
Fig. 135. Ultimate dowel capacity: (a) assumptions and symbols adopted to make the predictions of equations (153)
to (158) comparable to test results by Dei Poll et al.;'3 (b) effects of the position of the first stirrup on dowel
capacity
The predictions obtained with equations (153) to (158) are shown in Fig.
135(b), where the test results by Dei Poli et al. are also presented.130 As
regards cyclic loading, the same general considerations made in section 3.2
(dowel action against concrete core) hold for dowel action pushing against
concrete cover.
Under fully reversed transverse displacements (Fig. 136)133'164 the dowel
response becomes asymmetrical. Although Vintzeleou and Tassios do not
commit themselves, equation (152) seems to hold also for dowel action
pushing against concrete cover (with no stirrups).
In Fig. 137 the curves of stiffness versus number of cycles are shown:
under repeated loading the degradation turns out to be less severe than under
full reversals.
132
INTERFACE BEHAVIOUR
n th cycle
3rd
Fig. 136. Typical hysteresis loops for fully reversed Fig. 137. Degradation of normalized stiffness due to
transverse displacements: I + A l m a x | = | - Ai m a x | cycling:133 RL = repeated loading (pulsating loads); AL
(Vintzeleou and Tassios)13'' IM (dimensions in mm) = alternate loading (full reversals); AL* = alternate
loading after a series of cycles at low force levels
Summing up, only the joint effort of the people involved in basic research
on each single mechanism and of the analysts who introduce the different
models into their comprehensive computer programs, trying to fit the
various structural behaviours, will allow the code-makers to accomplish a
synthesis, and to help the professional engineers to design better and safer
structures, which is the ultimate scope of the whole chain of events
discussed.
133
6. Finite element modelling of reinforced
concrete
6.1. Introduction At the macroscale, composite materials are typically modelled as homo-
geneous, isotropic or not, with effective or equivalent properties. This is the
case for plain concrete (a composite of aggregates and cement paste), and
for fibre reinforced concrete. In reinforced concrete the mechanical
behaviour of its two constituents is very different, and the distribution of
reinforcement is typically non-uniform. So, a homogeneous description of
reinforced concrete is rarely used, and if so, only in regions where the
reinforcement is uniformly distributed. Instead, the two materials are
usually modelled individually. Such a separate representation of the two
materials, which fits very well with the finite element method, rests on the
(possibly questionable) assumption that the sum of the two constituent
materials describes reliably and without bias the behaviour of the actual
composite.
Models for the individual materials, concrete and steel (like the ones
previously presented in Chapters 1-3), for their interaction through bond
(Chapter 4) and for the behaviour of their interface (Chapter 5), can be used
as building blocks to construct fairly sophisticated models of reinforced
concrete at the microlevel. In models of this type separate finite elements
are used for the concrete (three- or two-dimensional ones, depending on the
dimensionality of the analysis) and separate ones for the steel (mono-
dimensional 'truss' elements representing a certain length of a single bar, if
the problem is analysed in three dimensions, or of a layer of bars with the
same two-dimensional position, if it is analysed in two dimensions).
Moreover, bond interaction between steel and concrete is modelled through
special elements, which connect nodes on steel elements to nodes on the
adjacent concrete elements, with the same initial position in space. Finally,
discrete cracks may be inserted between neighbouring concrete elements, by
splitting each node on their interface into two, and connecting the latter
through special contact or crack link elements, which model the behaviour
of the interface.
Integration of all the individual models into a comprehensive microlevel
model of reinforced concrete, capable of reliably predicting its behaviour up
to and beyond ultimate strength, under all possible loadings, including
cyclic, is not an easy task. Moreover, even if such a powerful tool is
available, its applicability to realistic reinforced concrete structures may be
limited by its prohibitively large requirements in computer time and
memory. For this reason, in recent years considerable effort has been
devoted in developing macrolevel composite steel/concrete models for
reinforced concrete. Models of this type treat reinforced concrete as a
continuum, with properties and behaviour which are the integrated outcome
of those of its constituent materials and of their interaction. Development of
macrolevel models has been promoted also by the belief of many leading
researchers that the interaction between steel and concrete is so strong and
complex that our real material is reinforced concrete, a composite with
properties which are not simply a superposition of those of its constituent
materials.
In the following, we will first present an overview of non-linear finite
element modelling of reinforced concrete with separate elements for the
concrete and the steel, and we will then examine 'macrolevel' composite
steel/concrete elements for uniform distribution of the reinforcement.
134
FINITE ELEMENT MODELLING OF REINFORCED CONCRETE
6.2. Finite element For discrete modelling of single bars (or of groups of bars with the same
modelling of the location, as for example the top and the bottom reinforcement of a beam) in
reinforcement and of two-dimensional or three-dimensional analyses, we use one-dimensional
its interaction with truss elements. These elements are typically two-noded, but if compatibility
concrete of displacement expansions with higher order concrete elements (e.g. eight-
noded isoparametric ones in two-dimensional analyses) is desired, higher
order one-dimensional elements (e.g. three-noded ones) must be used for the
steel bars. In axisymmetric analyses single reinforcing bars or groups of
concentrated bars in the hoop (circumferential) direction are modelled by
special single-node hoop elements, whereas constant density layers of hoop
and/or meridional steel within a surface of revolution, are modelled using
two- (or three-) noded axisymmetric membrane elements. In two-dimen-
sional or three-dimensional analyses, uniform density distributions of the
reinforcement (e.g. closely spaced stirrups in a beam, or reinforcement
meshes in a shear wall) may be modelled with two-dimensional or three-
dimensional elements, respectively, which are superimposed to (i.e. have
the same nodes as) the concrete elements. Such smeared reinforcement
elements can, of course, co-exist with one-dimensional elements for
concentrated reinforcement.
In all cases above, the constitutive relations of the steel elements are
derived in a rather trivial manner, considering that the normal stress and
strain in the direction of the reinforcing steel are related through the uniaxial
models presented in Chapter 3 for reinforcing steel.
When the overall macroscopic behaviour is mainly of interest, the nodes
of the steel and the concrete elements which have the same coordinates are
also assigned the same degrees of freedom. Then the interaction between the
two materials, such as bond-slip and dowel action effects, are modelled
implicitly, by appropriate, though phenomenological, modification of the
constitutive relations of steel and/or concrete. Incorporating tension stiff-
ening in the steel or concrete model is an example of such an implicit
representation of interaction effects. If the detailed and explicit description
of the interaction response at the microlevel is required, then nodes which
are connected with the steel and the concrete elements and have the same
(or practically the same) coordinates are assigned different degrees of
freedom, and special contact elements are inserted between them. Such
contact elements can be used only in association with the discrete, not the
smeared, representation of reinforcement.
The early version of these contact elements are the dimensionless 'bond
link' elements, connecting a single concrete node to the corresponding steel
node. The constitutive relation of these elements is equivalent to two non-
linear springs (three for three-dimensional analyses), relating lumped nodal
forces to the corresponding relative displacements between these two nodes.
One spring is in the direction of the steel bar(s) and models the bond-slip
effects. Any of the cyclic bond-slip models presented in section 4.4.2. of
Chapter 4, can be used for this purpose. The other spring(s) is normal to the
bar, and may model dowel action effects. If modelling of these latter effects
is not considered necessary, then infinite stiffness is assigned to these latter
spring(s), or the degrees of freedom of the concrete and steel nodes normal
to the bar are constrained to be the same. Coupling between the concrete-
steel interaction in the longitudinal and transverse to the bar directions, can
also be included in the stiffness matrix of the bond link element. Although
dimensionless, bond link elements require information regarding the
orientation of the steel bars to which they are connected. In this respect,
it is worth mentioning that bond link elements can very conveniently be
composed of standard two-noded truss elements, one of them along the
direction of the bar, and another (or two, if a three-dimensional analysis is
135
RC ELEMENTS UNDER CYCLIC LOADING
performed) normal to it. Each of these truss elements is connected with one
of the two nodes to be linked (the steel or the concrete node), and with an
extra node rigidly connected to the other one of the two linked nodes. So,
without requiring more degrees of freedom than a special bond link element,
the truss elements provide automatically the required orientation informa-
tion. One only needs to implement in the uniaxial constitutive relation of the
truss element the appropriate cyclic bond-slip or dowel action model.
To avoid the artificial lumping of interaction forces at the linked concrete
and steel nodes, and to improve the flexibility and the efficiency of
modelling the steel-concrete interaction, continuous contact elements have
been developed.171"174 These elements are arranged along the entire steel-
concrete interface, and they are typically isoparametric, with order com-
patible with that of the concrete and steel elements which they connect.
One-dimensional and two-dimensional (surface) contact elements in three-
dimensional space are used to connect one-dimensional discrete steel
elements or two-dimensional surface ones, respectively, to the surrounding
concrete. The elements have double nodes with the same coordinates, one
connected to a concrete element and the other to a steel element. The
constitutive model relating the three stress components in the contact
element (for the general case of three-dimensional analysis) to the three
deformation measures, offers the possibility of modelling the effect of
normal stress on the steel-concrete interface upon the bond-slip behaviour.
It is sometimes convenient to combine the continuous contact element with
the corresponding steel element into a single element, which incorporates
both the reinforcing steel and its interaction with the concrete.174
Nearly all widely used general purpose non-linear finite element
programs, which have a reinforced concrete modelling capability along
with their capabilities for metal plasticity, geomaterial and rubber
modelling,175"1 7 do not have special bond link or continuous contact
elements. So, in effect, they can account for bond between concrete and
steel only through modelling the tension stiffening effect. It is noteworthy
that Mehlhorn and his co-workers have implemented their contact element
into ADINA. 1 7 2 1 7 3
6.3. Modelling of the 6.3.1. Distinction between 'concrete in compression' and 'concrete in
concrete component tension'
It has become clear in Chapters 1 and 2, that the behaviour of plain concrete
under 'predominantly compressive' stresses exhibits fundamental differ-
ences from that under 'predominantly tensile' ones. For relatively low stress
levels, under which the behaviour of uncracked plain concrete is essentially
linear-elastic, such differences do not exist. However, when stresses ap-
proach and reach the failure criterion, the failure mode, and more impor-
tantly, the post-failure behaviour differ significantly: under predominantly
tensile loading, attainment of the failure criterion leads to cracking and
localization of extensional deformations along a well defined crack plane.
From then on the behaviour becomes strongly anisotropic: normal to the
crack plane the material strain-softens in the manner described and
modelled in Chapter 2, whereas parallel to the crack plane the behaviour is
essentially that of uncracked concrete with the dimensionality of the
problem reduced by one, so under triaxial stress conditions it becomes
biaxial, and under biaxial, it becomes uniaxial. (As discussed in detail in
section 6.4.3., in reinforced concrete the behaviour parallel to a crack plane
is not unaffected by what happens normal to the crack plane, due to the
tensile stress component which develops in this latter direction because of
the transfer of forces from the reinforcement crossing the crack through
bond, as well as other effects.) On the contrary, under predominantly
136
FINITE ELEMENT MODELLING OF REINFORCED CONCRETE
(fct is the uniaxial tensile strength of concrete), and in the two tensions-one
compression region
0"ci > o"c2 > 0 > crc3
by the plane
137
RC ELEMENTS UNDER CYCLIC LOADING
The failure surface in this latter stress subspace is 180 a Drucker-Prager cone,
described by the simple expression
K(tors) =
= ffgT D (tors)
138
FINITE ELEMENT MODELLING OF REINFORCED CONCRETE
(b) To compute the (total) stress vector ac, which corresponds to the
strain vector e, calculated at the integration point of interest from the
most recently computed nodal displacements. The total stress vector
is used then to form the vector of internal nodal forces of the element
T
<7cV (164)
For cyclic loading, or, in general, for load paths which are not monotonic/
proportional, differential (rate) type of constitutive relations only, such as
the ones reviewed in Chapter 1, are appropriate. Constitutive models of the
total (secant) type, such as those based on non-linear elasticity, should be
limited to monotonic, nearly proportional loading, as they can describe only
path-independent behaviour. It should be mentioned, though, that secant
type formulations allow full error control via equilibrium iterations, as the
total stress vector ac, and hence the vectors of internal nodal forces and of
unbalanced loads, is calculated exactly through the usually algebraic
expressions giving ac in terms of e. (Another advantage of the secant
formulations is that they can treat softening behaviour without numerical
problems and without resorting to special techniques, such as the arc-length
method). On the contrary, in formulations of the rate type, computation of
the total stress vector requires step-by-step numerical integration of the
(path-dependent) constitutive equations. Therefore, equilibrium iterations
through the unbalanced load vector can correct the solution only up to the
accuracy of the operator for the integration of the incremental constitutive
relation. This means that, unless a very fine step is used for this integration
operation, significant linearization errors may accumulate, leading to
considerable drift in the solution.
Among the incremental (rate) type of stress-strain models for concrete in
compression, the hypoelastic ones (either the orthotropic or the invariant-
based models) reviewed in section 1.3. seem most convenient, as far as
numerical implementation is concerned. A word of caution is due at this
point: as the loading/unloading criterion of such models is rather arbitrary
and lacking a sound theoretical basis, these models may violate the
continuity requirement for nearly neutral loading. However, models based
on the theory of plasticity and employing convex yield and loading surfaces
and an associated flow rule, satisfy Drucker's stability postulate which
guarantees uniqueness of solution. Nevertheless, as noted in section 1.4.1.,
experimental evidence and theoretical arguments regarding the non-linear
behaviour of frictional materials, suggest that validity of the normality rule
for concrete is neither a realistic nor a necessary assumption. The
alternative, i.e. non-associated plasticity (as applied in DIANA for
example), 182 gives rise to non-symmetric rigidity and stiffness matrices,
and requires a non-symmetric equation solver.
The criteria that seem to govern selection of the constitutive model for
concrete in compression to be implemented in a non-linear finite element
code, especially a general purpose three-dimensional one, are the following.
The model should be simple and numerically stable, yet it should be
theoretically sound and capable of reflecting the important features of the
139
RC ELEMENTS UNDER CYCLIC LOADING
140
FINITE ELEMENT MODELLING OF REINFORCED CONCRETE
141
RC ELEMENTS UNDER CYCLIC LOADING
6l
4 1
J10I
ir
NODF-
I 7 8\
Steel element Crack element
/
ELF
A Dowel-action element
All these variants of the discrete cracking approach require pairs of nodes
at opposite sides of a crack. So, either new nodes have to be added to the
model upon formation of a crack, or double nodes have to be provided from
the beginning along all interelement boundaries where cracks may poten-
tially develop, and allowed to separate when the crack forms and opens up.
In the first case every time a new crack is formed or an existing one
propagates, the finite element topology changes, the band width has to be
reminimized, etc. In the second case, a large number of extra nodes, and
corresponding linkage or contact elements along interelement boundaries,
need to be added to the model from the very beginning.
The discrete cracking approach of the two-dimensional special purpose
program MICRO developed at Delft,176178 deserves special mention. The
program uses special hybrid triangular elements for uncracked concrete,
with linear interpolation of stresses over the element and linear interpolation
of displacements along the element sides. (The stress and displacement
distributions are related through certain conditions that guarantee
equilibrium within the element and minimize lack of interelement equili-
brium along the boundaries.) When the cracking criterion is reached at the
centre of gravity of an element, a discrete crack is introduced through this
latter point. A discontinuous displacement interpolation is added to the
linear ones along the two sides intersected by the crack, amounting to three
new degrees of freedom: one slip translation along the crack, and a crack
width that varies linearly along the crack (Fig. 140(a)). The same is done to
the stress field in the element, which also becomes discontinuous (Fig.
140(b)). Appropriate relations are introduced between the new displacement
and the new stress degrees of freedom across a crack according to
Hillerborg's fictitious crack model,184 see section 2.3, and to Walraven's
model for aggregate interlock.'18 This is equivalent to introducing a contact
element along the discrete crack, which connects the two faces of the crack
through the normal stress-crack width (see Figs 33, 34, 45 and 46 of Chapter
2) and the shear stress-shear slip-crack width relations. Only one point of
intersection with cracks is allowed along each side of a concrete element.
So, if such a point of intersection exists already, a new crack, in the same or
in the neighbouring element, is forced to pass through the same point. In this
way a multilinearly continuous crack may form, crossing several elements
(Fig. 140(c)). The occurrence of cracking imposes no mesh rearrangement
or disconnection between elements. Moreover the degrees of freedom are
arranged in two sets: the first one, containing the original degrees of
freedom prior to cracking, is fixed, while the other, with the new ones, is
augmented every time a new crack forms. The two matrix equations for the
two sets of degrees of freedom are coupled through the unbalanced load
142
FINITE ELEMENT MODELLING OF REINFORCED CONCRETE
(b)
(c)
Fig. 140. Discrete crack modelling in program MICRO: (a) discontinuous displacement interpolation along crack;
(b) discontinuity in stress fields, for one and two cracks; (c) possible configuration of cracking
143
RC ELEMENTS UNDER CYCLIC LOADING
mentioned above 174176 employ such crack contact elements, with crack
cohesive and shear stresses dependent on crack width and shear slip.
The discrete crack approach lends itself better to the modelling of the
interface shear phenomena, and allows direct calculation of crack width. So,
whenever the behaviour is governed by one or by a few dominant cracks,
and/or the local behaviour is of interest, the discrete crack approach is the
natural way to go. However, in its general form, in which cracking is
allowed to form anywhere in the finite element mesh, this approach is
numerically inconvenient, as it requires redefinition of the mesh and
introduction of new degrees of freedom. From the computational point of
view, this is a serious drawback. In addition, up to the present time the
discrete crack approach has been applied only to two-dimensional problems
under monotonic loading, and its extension to three-dimensional conditions
and to cyclic, or generalized loading is not immediately obvious. Regarding
cyclic loading, it should be noted that although cyclic versions of the
individual constituent models of the discrete crack approach are available
(see the cyclic fictitious crack model presented in Chapter 2, and the
interface shear models with cyclic capabilities in Chapter 5), finite element
implementation of a general cyclic discrete crack model still has to wait for
the solution of important issues regarding closing of cracks and opening of
new ones in other directions, reopening of cracks, etc.
Because of its computational inconvenience and its essentially fixed-
crack character, the discrete crack approach enjoys today limited acceptance
and application. Nevertheless, when a major discontinuity in the geometry
(e.g. at the junction of two different structural elements) predetermines the
location of dominant cracks the best way of modelling the problem seems to
be to combine discrete cracks at the predefined location with a smeared
crack modelling for the rest of the structure. It is worth mentioning that a
special interface contact element for cyclic loading has been developed by
combining cyclic models of interface shear along the crack and of bond-slip
along the inclined bars crossing it.186 This special element has been
successfully used,186 at the junction of a cyclically loaded shear wall with
its foundation, the rest of the wall and the foundation being modelled using
composite steel/concrete elements with smeared cracking. The reader
should be cautioned at this point that continuous interface contact elements,
although apparently more general than the lumped ones inserted between
nodes in the early applications of the discrete crack concept, have been
found to introduce artificial fluctuations in the fraction profile along the
interface, when the crack is closed (a situation simulated by means of a very
high dummy stiffness of the interface element, in the direction normal to the
crack).187 So, in this respect the more primitive lumped interface contact
elements seem to be superior.
The discrete crack approach fits the physical concept of fracture, as
cracks are introduced via discontinuities in the displacement fields. For
problems with distributed fracture, however, as for example in shear walls
with densely distributed reinforcement and diffuse cracking, the smeared
crack approach seems closer to physical reality. According to this latter
approach, cracks are assumed to be distributed (smeared) over a concrete
element, or over the tributary volume of an integration point, and the local
displacement discontinuity at crack locations is smoothed out over this latter
volume, through the displacement interpolation functions of the finite
elements. Only the concrete tangent rigidity matrix at the (cracked)
integration point or element has to be modified upon and after cracking,
without changing the finite element topology. Due to this latter feature, and
as it imposes no restrictions on the crack direction, the smeared crack
approach has gained much wider acceptance than the discrete crack
144
FINITE ELEMENT MODELLING OF REINFORCED CONCRETE
145
RC ELEMENTS UNDER CYCLIC LOADING
taken to end at 20-80 times/ ct /£, l79-190>191 with higher values corresponding
to heavier reinforcement. Therefore, for the values of h commonly used in
modelling, tension stiffening completely overshadows softening of
concrete.
The (uniaxial) a-e curve of concrete in tension, with its mesh size-
dependent falling branch for mesh objectivity, is used in the multiaxial a-e
model of concrete in tension. It is illustrative to refer once more, for
example purposes, to the two major general purpose three-dimensional
programs, ADINA and ABAQUS.175'177 In the latter's associated plasticity
model with isotropic hardening, the user is expected to provide as input a
multilinear a-e curve of the concrete in uniaxial tension, including the
falling branch. This user-supplied curve is used then by the program to
establish internally the dependence of plastic modulus on strain, for the part
of the model that applies for predominantly tensile stresses. ADINA, on the
other hand, constrains the falling branch of the cr-e curve of concrete in
uniaxial tension to be linear, with a user-supplied terminal strain. After
cracking, the tangent modulus of concrete normal to the crack is reduced to
a very small positive value, in the orthotropic tangent rigidity matrix
equation (6) in Chapter 1, and equilibrium iterations are performed using the
stress value computed from the falling branch of the a-e curve. This
approach is preferred over using the negative stiffness of the falling branch
in the tangent rigidity matrix, to avoid numerical problems. In both ADINA
and ABAQUS, tensile unloading and reloading from a point on the falling
branch is along a straight line from the origin to the point of maximum
previous tensile strain.
A word of caution is due at this point: as neither of these two smeared
crack-based programs include any bond link or contact elements to simulate
the bond-slip interaction of steel with concrete, and as both of them use a
bilinear a-e curve for steel in uniaxial tension, which does not lend itself to
modification to account for tension stiffening, this latter effect can be taken
into account only through the falling branch of concrete in tension.
Therefore, in these two programs this latter falling branch should be used to
model the tension stiffening effect of reinforced concrete, rather than the
mesh size-dependent softening of the concrete component. This means that
the terminal strain of the user-supplied falling branches should be assigned a
fairly high value.
The smeared crack approach has been criticized as incapable of pre-
dicting localization of fracture. However, it produces proper strain localiza-
tions for mixed mode cracking in both plain and reinforced concrete
elements,64 regardless of whether one adopts a fixed or a rotating crack
model.
146
FINITE ELEMENT MODELLING OF REINFORCED CONCRETE
147
RC ELEMENTS UNDER CYCLIC LOADING
148
FINITE ELEMENT MODELLING OF REINFORCED CONCRETE
(167
»
in which Go is the shear modulus of uncracked concrete, the sum in the last
term extends over the cracked directions (to direction / only, or to both / and
j , if there is cracking only normal to direction i, or in both directions,
respectively) and the shear slip modulus C s is the ratio of 77, times the crack
opening wh over the resulting shear slip Sy, along the crack (assuming that
such a proportionality between 77/w, and Sy exists).
It is worth mentioning that even though in many cases response predictions
may seem relatively insensitive to the value of the shear retention factor, a
non-zero shear retention factor must be included in a fixed crack model, to
avoid numerical instabilities and spurious kinematic modes associated with
the displacement degrees of freedom in the direction normal to a crack, if all
concrete elements connected to a node crack in the same direction.
149
RC ELEMENTS UNDER CYCLIC LOADING
1 -1 -cot 20
(o-cl — crc2) sin220
-1 1 cot 20 (168)
8(ex - €yy) 2
-cot 20 cot 20 cot 20
in which ex and ey are normal strains in global directions x and y, and 0 is the
angle from the x axis to the instantaneous direction 1, measured counter-
clockwise.
The value of the effective shear modulus is less important than in the
fixed crack approach, as the shear strains are zero in the instantaneous
directions of principal strain. Its continuously changing value is implicit in
the tangent rigidity relation between da and de in the global system. As
shown by equation (166), for two-dimensional stress conditions the
coincidence of the principal stress and strain directions yields an effective
shear modulus equal to (crcl — crc2)/2(ei — e-i). It has been shown187 that if
this latter value of the effective shear modulus is used, the threshold angle is
taken equal to 0°, and elastic unloading of the previous cracks is assumed,
then the multidirectional non-orthogonal crack model approach produces as
a limit the rotating crack model. In comparison to the general form of the
former, the rotating crack approach is much simpler and less demanding in
computer memory, at the expense, of course, of the potential for more
realistic modelling of the interface shear behaviour, offered by the
multidirectional non-orthogonal crack approach. Comparative applications
of the rotating crack model and of the non-orthogonal fixed crack one, have
been made on plain and reinforced concrete examples. In the fixed crack
model, various values of the threshold angle from 0° to 90° and of the shear
retention factor, /?, have been adopted. This factor has been taken to
decrease with the ratio a = ef/e™, of the strain normal to the crack to the
falling branch terminal strain in the fictitious crack normal stress-strain
model, as (J = (1 - a)2/[\ - (1 - a) 2 ]. The results have shown that the
rotating crack approach produces a slightly less stiff response than the non-
150
FINITE ELEMENT MODELLING OF REINFORCED CONCRETE
orthogonal fixed crack model with shear straining, even becoming negative.
In addition, the stress normal to the crack decreases with shearing. So, the
crack rotation introduces implicit shear softening and shear-normal stress
coupling in the rotating crack model, without the need to resort to the
complicated interface shear modelling required by fixed crack models.187
These same examples have shown the inadequacy of the fixed orthogonal
and non-orthogonal crack model with large values of the threshold angle.
So, overall the rotating crack approach seems to achieve the best
combination of simplicity, accuracy and convergence characteristics.
The smeared crack approach seems to suffer from stress-locking in the
vicinity of dominant discrete cracks. This problem is less serious if one uses
the rotating crack approach, or the fixed crack approach with a zero shear
retention factor ((3 = 0). In all other cases significant stresses are predicted
in directions non-parallel to the crack. As shown by the more realistic
predictions of discrete crack models for the same problem, these latter
stresses are spurious, and are responsible for the excessive stiffness of the
response. The problem of stress-locking seems to be inherent in the smeared
cracking approach, and is attributed to the imposed displacement continuity
between softening cracked elements and neighbouring uncracked ones: as
strains in the former elements increase normal to the crack, neighbouring
elements are forced to follow and their stresses increase spuriously, leading
even to the artificial extension of cracking into them, which works against
strain localization and is also undesirable. The stress-locking problem
can be reduced by employing the rotating crack approach (the fixed crack
approach with /3 — 0 leads to spurious kinematic modes and instability
problems), but it can be avoided only through use of a discrete crack model.
The fixed crack and the rotating crack approaches have been proposed
and used, the former far less successfully than the latter, for monotonic
loading, proportional or not. For cyclic loading in two dimensions, two
variants of these approaches have been developed, to deal with the
formation of cracks in two nearly orthogonal directions, and for their
opening and closing due to the reversals of loading and response. The
former approach, proposed by Okamura and his co-workers,203'2 is a fixed
crack approach with provisions for formation of a second (fixed) crack and
for shifting the directions of orthotropy to the latter (with direction 1 normal
to the new crack), if the smeared strain normal to this second crack becomes
greater than that normal to the first.
The shear stiffness of the Okamura model is based on a non-linear
interface transfer model with the following monotonic curve, which serves
as an envelope for cyclic response
/3 (169a)
TT^
in which a = 712/ei is the ratio of shear strain to that normal to the crack,
and fc is in MPa. The compressive stress component developed normal to
the crack due to the interface shear is equal to:
(169b)
151
RC ELEMENTS UNDER CYCLIC LOADING
at the stress reversal, from where reloading starts on a straight line passing
through a point on the vertical axis at an ordinate equal to 010 of the value
of T\2 at reversal. Reloading follows this line up to its intersection with the
envelope curve in the opposite direction of loading, for first loading in this
latter direction, or up to the last unloading branch from it, for subsequent
cycles. From then on, the envelope curve or this latter branch is followed.
Finally, the value a — 0-5 is specified as a criterion for shear failure.
The second approach for cyclic loading is that proposed by Stevens et
al. as a modification of the rotating crack approach followed in the
numerical implementation of the modified compression field theory,202 to
cover cyclic loading in two dimensions as well. This approach has been
developed on the basis of the results of cyclic shear tests on reinforced
concrete panels, which have shown that during certain parts of the cyclic
response the principal directions of strain lag behind those of stress by
approximately 90°. This observation instigated the alternative assumption
that the principal directions of the concrete stress increment tensor da
coincide with those of the strain increment tensor de of the reinforced
concrete continuum. In the finite element implementation of the model, at
each step of the incremental analysis the strain increment tensor de is
computed at each integration point from the displacement increments of the
step, and the corresponding principal directions of the strain increment
tensor are determined. The stress and strain tensors computed at the end of
the previous step of the analysis are transformed to the local coordinate
system defined by the aforementioned principal directions. Then the normal
stress increments in these two directions are calculated from the corres-
ponding normal strain increments and from the values of normal stress and
strain in these two directions at the end of the previous step. (The shear
stress increment is zero, as these two directions are considered as principal
directions of the stress increment tensor.) Calculation of these two normal
stress increments, and determination of those terms of the tangent rigidity
matrix D c which relate the normal stress and strain increments, is performed
independently in each direction, on the basis of the cyclic law adopted for
the cracked concrete, described later on. The shear term on the diagonal of
Dc is continuously changing and equal to (dcrci — d<7c2)/ 2(de, — dt2) (see
the corresponding term for the rotating crack approach based on total
stresses and strains. For a linear-elastic isotropic material, this shear term of
both approaches equals one-half the modulus of elasticity). However, as this
shear term refers to principal directions of the stress and strain increments,
its exact value is immaterial. Finally, the tangent rigidity matrix D c is
transformed from the coordinate system of the instantaneous principal stress
and strain increments to the global coordinate system xy.
The ability of both approaches described above for cycling loading to
reproduce cyclic test results on two-dimensional reinforced concrete panels
has been verified.203"205 However this capability seems to be largely due to
features of the models other than their provisions for crack rotation.
152
FINITE ELEMENT MODELLING OF REINFORCED CONCRETE
case results reflect only the trend of the variation of average bond, steel and
concrete stresses along the bar, to which such bond-link or contact elements
are connected. Local variation between cracks is, obviously, unaccounted
for, since the model cannot provide any information regarding the location
of individual cracks. For smeared cracking, it is computationally more
efficient to assume perfect bond between steel and concrete and account
implicitly for the effect of bond-slip on the average stresses and strains in
the steel and the concrete, through appropriate, yet empirical, modification
of the properties of these two materials. So, at the expense of some accuracy
in the description of the variation of average steel and concrete stresses,
significant savings can be effected, as the extra degrees of freedom of the
bond-link or contact elements at the steel-concrete interface can be
eliminated, and these special elements do not need to be developed and
added to the finite element library. The modification of steel or concrete
properties for this purpose, is effected through appropriate incorporation of
the tension stiffening effect in their constitutive models.
At this point it seems opportune to remember the essentials of the tension
stiffening effect in reinforced concrete, on the basis of a prismatic axial
member, cracked due to concentric axial tension. At the cracked cross-
sections, concrete stress is zero and steel stress is maximum. Due to transfer
of stresses from the bars to the concrete through bond, tensile stresses
develop in the concrete between the cracks, attaining their maximum value
at crack mid-distance. So, the average concrete stress along the element is
non-zero. Because of the transfer of tensile stresses to the concrete, steel
stress decreases with distance from the crack, becoming minimum near mid-
crack. Equilibrium requires that at any cross-section steel stress is equal to
that at the crack minus the concrete stress divided by the steel ratio, p. The
smeared (i.e. average) axial stress in any cross-section of the reinforced
concrete element is equal to the sum of the corresponding concrete stress
plus the corresponding steel stress times p, and remains constant along the
element. The smeared axial strain of the element is equal to the average
value of the steel strain between two cracks, and is considerably less than
the steel strain at the cracked cross-sections. Actually, yielding of the steel
takes place only at these cross-sections and in their immediate vicinity.
Immediately after yielding at a cracked cross-section, the reinforcement
reaches strain-hardening, although at a short distance from the crack bond
may have reduced steel stresses and strains to below the yield values.
As the smeared axial strain is considerably less than the maximum steel
strain, which corresponds, through the stress-strain curve of steel, to its
stress at the cracked cross-section, the plot of this latter steel stress against
the smeared axial strain lies above the stress-strain curve of bare steel. On
the contrary, the plot of the average steel stress along the bar against smeared
strain lies below this latter curve. Finally, the plot of the average concrete
stress along the element against smeared strain exhibits a long post-cracking
tail, far above and beyond the falling branch of plain concrete in tension. All
these effects are demonstrated in the experimental results of Fig 142, taken
from Shima et al.,206 and constitute the well-known tension stiffening effect
in cracked reinforced concrete. The term is due to the stiffening of the bare
steel, due to the participation of concrete in tension between the cracks. Its
underlying reason is the bond between steel and concrete. The lower the
reinforcement ratio, the more pronounced the tension stiffening effect is, and
the more it affects the post-cracking deformations of reinforced concrete.
Extensive experimental studies of the tension stiffening effect have been
performed, under monotonic, repeated and reversed loading. The works of
Hartl,207 Guenther208 and Fehling209 from which the results in Fig. 143 are
drawn, are good examples.
153
RC ELEMENTS UNDER CYCLIC LOADING
Strain distribution
A A A y
55
154
FINITE ELEMENT MODELLING OF REINFORCED CONCRETE
80000
60000
40000
20000
same end result, as both reproduce the desired total stress-smeared strain
curve. However, for concrete reinforced in two (usually orthogonal)
directions, loaded biaxially and cracked at an angle to the reinforcement, the
two-dimensional extensions of the two approaches are based on different
(influential) assumptions, and therefore yield different results. For the two-
dimensional problem Gupta and Maestrini183'210 have developed an 'exact'
analytical solution. Their work is based on the assumptions of
155
RC ELEMENTS UNDER CYCLIC LOADING
J & ("2)
Equation (172) has been fitted to a large number of experimental points,
which exhibit considerable scatter. An alternative to equation (172) has
been proposed more recently by Collins and Mitchell,212 yielding smaller
and more realistic values of <rci for large values of ei
f r e >
° ^ <173>
in which a, = 1-0 or 0-7 for high bond or smooth bars, respectively, and
a2 — 1-0 or 0-7 for short-term monotonic loading, or for sustained or
repeated loading, respectively.
In DIANA and according to Crisfield and Wills191 the falling branch is of
the form of equation (171) with k2 = 10 and with the value 0005 in the
denominator replaced by a strain parameter. In the application of DIANA by
Wang et a/.179 the value of this strain ranges from 0001 to 0003 and is
limited by the yield strain of steel. Crisfield and Wills state the terminal
strain as equal to 0-004. It is also to be remembered, that the user-specified
156
FINITE ELEMENT MODELLING OF REINFORCED CONCRETE
fore!>/ct/£c (174)
(1 +1000 c, (0/9O)15)
in which 0 is the angle in degrees between direction 1 (the normal to the
crack) and that of the reinforcement with the higher steel ratio. It is
noteworthy that for 6 approaching zero equation (174) yields stress values
close to or equal to fct, regardless of the strain value. These unrealistically
high values of CTCI are limited, however, by the requirement that the total
stress in the direction of any reinforcement layer should not artificially
increase beyond that corresponding to yielding of this layer. This
requirement is satisfied if
i/ (175)
in which asi is the average steel stress in reinforcement layer /, and ^, is the
angle between direction 1 and the bars of layer /. A similar limitation is also
imposed on the relatively high values of ac\ obtained from equations (172)
and (173) for large values of ei. 205 This is accomplished through a more
complicated procedure of checking whether shear stresses can be transferred
by the cracked interface, without exceeding the yield stress of steel.
All the aforementioned proposals have limited provisions for cyclic
loading. Unloading from the falling branch of the concrete ac\ — ei curve is
linearly oriented towards the origin, and reloading is again from the origin
towards the point of maximum e\ value in the previous response. 175177191
(It is noteworthy that in the rotating crack approach by Crisfield and
Wills,191 this latter maximum ever value of e\ may be associated with a
different direction of the principal tensile strain, e\.) Moreover, the above
proposals include no reduction of the average steel stress-smeared strain
curve below that of bare steel, due to the incorporation of the entire tension
stiffening effect in the stress-strain curve of concrete. On the contrary, the
two models described below are more complete in both these respects.
Okamura and his co-workers203'204 have proposed a relation between ac\
and ei which consists of a monotonic curve that serves as an envelope for
the cyclic response, and of rules for unloading and reloading. The
monotonic curve comprises the familiar ascending branch, equation (170),
a horizontal branch thereafter up to a strain of 2fct/Ec, and a falling branch of
the following form
For the exponent c the values 0 • 2 and 0 • 4 are suggested for welded wire
mesh or deformed reinforcement, respectively. The curves of unloading
from the envelope and for reloading to it consist of a superposition (through
addition of the stress-ordinates) of two stress-strain curves. The first is
intended for modelling the interaction with the steel through bond, and has
parabolic unloading from the point a c i e — ei>e on the envelope to a point on
the vertical axis at an ordinate of —0 • 0016£ei)e, at which the parabola has
horizontal tangent. Reloading from that latter envelop point, or from any
other intermediate point on the unloading parabola, is linear, directed
towards the envelope point <7cie — ei,e- The second component of the
157
RC ELEMENTS UNDER CYCLIC LOADING
Jy / \PxJ
/30N
\/c /
(|7g)
in which 6 is the angle of the reinforcement layer with ratio px with respect
to direction 1 (the normal to the crack) and / c , / y are in MPa. For welded
wire mesh coefficient 100 is replaced by 72.
The other model which associates tension stiffening with concrete but
covers cyclic loading as well, is the aforementioned model by Stevens et
al.205 In this model the envelope for concrete average stress-smeared strain
in tension, which is identical to the curve for monotonic tension, consists of
the linearly-elastic branch up to tensile strength, equation (170), and a
falling branch thereafter, that includes the entire effect of tension stiffening
and is given by
158
FINITE ELEMENT MODELLING OF REINFORCED CONCRETE
Region 1
Region 2
V7\ Region 3
Compressive
envelope
Ideal reloading curve
L \ N Cracks open
Y/ X Cracks closed Hypothetical unloading path
/ A(c c a n , fc
I A'fcm. U)
Cracking point
'cm Compressive
envelope E
) /
Hyperbolic
unloading path,
Equation 2 1 A— Ideal reloading curve
fen
Equation 1 v
Tensile
•
y
TV/
/A fen
Cracking point
envelope
envelope
Fig. 144. Rules for unloading—reloading from compression to tension in cracked reinforced concrete according to
Stevens et al.20
159
RC ELEMENTS UNDER CYCLIC LOADING
f2h
£ = ^ (181a)
160
FINITE ELEMENT MODELLING OF REINFORCED CONCRETE
161
RC ELEMENTS UNDER CYCLIC LOADING
D = T 1 D c T c + E P / T* D s/ TS(- (182)
162
FINITE ELEMENT MODELLING OF REINFORCED CONCRETE
163
RC ELEMENTS UNDER CYCLIC LOADING
f; = PA (184a>
and the corresponding strain at the peak, e*o, to
164
FINITE ELEMENT MODELLING OF REINFORCED CONCRETE
<xo =
(185a)
in which the values fc and eco > 0 are those corresponding to uniaxial
compression of concrete, and the coefficient fi equals
1
<1 (186a)
0-85-0-27-
and
165
RC ELEMENTS UNDER CYCLIC LOADING
The main feature of Vecchio and Collins proposal, which is the reduction
of compressive strength in the direction of ei with increasing value of the
transverse tensile strain e\, according to equations (184a) and (186a) or
(184b) and (186b), is not in major disagreement with the European
proposals for the reduction in the value of concrete compressive strength, in
the usual range of ei/e 2 values found in practice. For values of this latter
ratio below - 4 0 and down to - 2 0 , Vecchio and Collins' proposal gives
very large reductions of compressive strength, of the order of 80%, which
may have to be viewed with some caution. Nevertheless, such low values of
e\/e2 are very rarely associated in practice with ultimate strength, especially
when high bond bars are used.
Cervenka has also proposed equations (184a) (185a) and (187)—(189) but
with the following coefficient (3, for the reduction of strength and strain at
the peak
166
FINITE ELEMENT MODELLING OF REINFORCED CONCRETE
tne
Unloading from a strain €2,e> o° envelope is parabolic, towards a point
on the horizontal axis at a strain
€pi = -2-85[1 - exp(0-35e2,e)] (193)
and reloading is linear, towards the previous envelope point.
167
7. Solution strategies for non-linear structural
equations
7.1. Non-linear The discrete non-linear equations which express, at a given time t,
equilibrium equations equilibrium of a structure, modelled through finite element space discret-
ization, may be written as
= 0 (194)
F e x t contains the externally applied nodal forces, (equivalent to body,
surface and concentrated forces), and F i n t the internal nodal forces, that is
the structural reactions associated with the present configuration; more
explicitly F i n t = S[(u(t)] is a function of the evolution of the unknown nodal
displacements u. Time t has the meaning of physical time, in dynamic
problems, or of a parameter associated to load levels, in static problems. In
dynamic analysis the vector F e x t , besides given external loads P, would also
contain inertia and damping forces, that is
= P(t)-Mu(t)-Cu(t) (195)
At)-u(t)]=Sn
= «n (198)
168
SOLUTION STRATEGIES FOR NON-LINEAR STRUCTURAL EQUATIONS
( ^2
4
(«n+l -Un)-—Un-
4
«
and subsequently, using the already given linearization of S[u(t + At)],
K*Au=R* (197b)
where K* and R* are, in this case, respectively defined by
169
RC ELEMENTS UNDER CYCLIC LOADING
V At2
and
up to the point where the residue vector is made as small as desired. The
described methodology implies that at every iteration, within a step, the
tangent matrix K(l) and the residue vector /? (l) be evaluated. Such a strategy
is referred to as 'full Newton-Raphson method'.
Since the updating and factorizing, for solution, of the tangent stiffness at
every iteration is a computationally expensive operation, different alter-
natives have been proposed; these may be considered as modifications of
the original Newton-Raphson method.
One such alternative, called 'initial stiffness method', is to use the initial
stiffness matrix KQ throughout the steps and the iterations, by operating
always with the equations
K o Au (/) = R{i)
In this case the stiffness matrix is formed and factorized only once during the
analysis; the method corresponds to a linearization of the response about the
initial configuration. For problems with significant non-linearities this can
lead to very slow convergence or even divergence, if the system stiffens.
In order to accelerate convergence or avoid divergence it may be
effective to use the so-called 'modified Newton-Raphson strategy', where a
new tangent stiffness matrix is calculated, from time to time, in order to
follow the evolution of non-linearities; within a step such matrix is always
kept constant during iterations. Without any a priori knowledge of non-
linearities it may be effective to reformulate stiffness at the beginning of
every step.
Figures 145-147 give a graphical representation, in the case of a one
degrees of freedom system, of the different alternative schemes considered.
170
SOLUTION STRATEGIES FOR NON-LINEAR STRUCTURAL EQUATIONS
Pn + 1 xT-^
Pn+
/?W is chosen, on the basis of a line sesarch algorithm, so that the out-of-
balance loads, R{l+l) = /?(M (0 + /?W A«W) have a nearly zero component in
the direction AM. This condition is satisfied to a given tolerance.
Formulas similar to equation (200a) are used, for updating displacements,
in other acceleration schemes, for example the so-called Aitken acceleration
scheme,217"219 with the only difference that (3 is no longer a scalar, but a
diagonal matrix; this is continuously updated on the basis of information
acquired during iterations.
171
RC ELEMENTS UNDER CYCLIC LOADING
(201)
where K is the last formed and factored stiffness matrix. Quasi-Newton
trix. Quasi-N
updates are often used with line searches, so that from
V (202)
The BFGS vectors v and w are defined
v = Att('- 1 7(A«('- 1 ) T A*'1'"1)) w = -A*''"1) + aR«-l\
where a = [(-0 AR^T Au^~• >)/(*('-'> T Au('-'>)] l/2
Recalling that K /? ( / - 1 ) = A« ( '~'), it is not difficult to show that K, as
implicitly defined by equation (201), satisfies the quasi-Newton equation
(202). Proper recursive application of equations (199) and (201) enables the
search direction to be computed at different iterations using only vector
operations, that is, back substitutions and forward reductions, without
explicit computation of updated matrices or costly factorizations; the
factored K is never altered during such calculations (see Matthies and
Strang220).
7.2. Static analysis So far external loading has been assumed to be a given function of time.
(special problems) This means that, in static analysis, the analyst prescribes the load levels for
which the equilibrium configurations are to be calculated; moreover, during
equilibrium iterations, such load levels have been considered to remain
constant. If load-displacement processes have to be traced in which there is
no a priori knowledge of the load-time function, special path-following
methods are needed. This is important, for example, when the structural
response about and after limit (collapse) points, (Fig. 149), is of interest.
Among these methods a short account is given here of the displacement
control methods and the arc-length methods. Both methods provide for
automatic step selection and non-constant load iterations.
AA W )P - S(« w )
172
SOLUTION STRATEGIES FOR NON-LINEAR STRUCTURAL EQUATIONS
Limit point
Pn + 1
Fig. 148. Quasi Newton methods; updated K Fig. 149. Analysis beyond limit (peak) point
which can also be written, using the residue vector computed at previous
iteration, JfW = ® W
K(0 (203)
K w AMJ° = P
and expressed by
A H W = AAW A H I 0 + . .(')
,221
In the displacement control method suggested by Batoz and Dhatt "" a
single displacement component is selected as a controlling parameter during
a step. At the beginning of the step that component is given the prescribed
value AMJ0^ = 6 and this constraint determines the first value of the load
increment A\(°\ In subsequent iterations the same displacement component
is always given zero increments, A«| () = 0, and this determines subsequent
value of the load increment corrections AA^. Iteration is continued, at
varying load levels, until a new equilibrium position is found (R = 0).
In the arc length method similar conceptual lines are followed, with the
main difference being that the constraint which determines AA^ during
iterations is no longer a single displacement parameter, but a more complex
function of an arc length of the solution curve in force displacement space.
Both methods allow the use of different strategies for the stiffness
reformulation and solution of equations (equation (203)), (full Newton-
Raphson, modified Newton-Raphson, BFGS updates), and may include line
search or other acceleration schemes in the unknowns updating
stage.219'222'223
7.3. Dynamic In non-linear dynamic analysis of structures equilibrium equation (194) may
analysis be more explicitly rewritten, using equation (195) as
M u(t) + C ii(t) + S(u(t)) = P(t) (204)
where S(u(t)), the vector of structural reactions, contains the non-linear
173
RC ELEMENTS UNDER CYCLIC LOADING
2 +
«n = ( )
and
( 2 o 6 )
174
SOLUTION STRATEGIES FOR NON-LINEAR STRUCTURAL EQUATIONS
from which it is possible to solve for the unknowns « n+1 , assuming that
displacements, velocities and accelerations at previous time steps are
available.
It is a common characteristic of explicit methods that they contain in the
expression of velocities and accelerations, at time n, the values of
displacements at time n+l. This permits the writing of equilibrium
equations at time n, in which the unknowns un+i only appear in the
expressions of inertia and damping forces and not in the terms which refer
to the structural reaction vector. Since the mass matrix is often diagonal and
damping forces are either completely neglected or simplified with the
assumption of a diagonal damping matrix, the solution may be advanced
without solving any system of equations. If this is the case, since no
procedure of triangularization of a matrix is required, there is no need to
assemble a matrix like K, but only force vectors representing structural
reactions, and the solution can be essentially carried out at the element
level, with important savings of high-speed storage.
A disadvantage of explicit methods is, on the other hand, that they require
the time step At to be smaller than a critical value which can be calculated
from the mass and stiffness properties of the complete structure. More
specifically this critical value is a fraction of the smallest period of the finite
element assemblage with m degrees of freedom. Integration schemes which
require the use of a time step smaller than a critical value are said to be
conditionally stable. Since the effective use of conditionally stable methods
is limited to certain problems, it is important to consider schemes which are
unconditionally stable. These can be found only among the implicit schemes
which belong to the category illustrated next.
7.3.1.2. Implicit methods. An important and simple example of this
category, already used as an example in 7.1.1., is the so called trapezoidal
rule or constant average acceleration rule. All the implicit methods are such
that displacements at time n+l depend not only on displacements, velocities
and accelerations at previous steps, but also on the quantities at step n+l.
In order to advance solution, equilibrium equations need to be written at
time n+l, where all the terms (inertia, damping and structural reaction
forces) depend on the unknowns un+i. This implies that at each step a
stiffness matrix has to be assembled, and a system of equations has to be
solved. The trapezoidal rule is defined by the following formulas
M+ c
A7 ) " " + M ""
Numerical integration can now be advanced solving the above linear
175
RC ELEMENTS UNDER CYCLIC LOADING
176
SOLUTION STRATEGIES FOR NON-LINEAR STRUCTURAL EQUATIONS
be measured in terms of period elongation and amplitude decay (in the exact
solution there is no decay, since there is no physical damping), as a function
of the time ratio At/T.
The analysis for various integration schemes shows that all methods are
accurate, that is, they show negligible period elongation and amplitude
decay, for At/T smaller than about 0 0 1 . The differences come out for
larger At/T.
The characteristics of the integration errors exhibited by the test problem
for different At/T may be used in the discussion of the simultaneous
integration of m equation (208) which is equivalent to the integration of
equation (205). We observe that the equations for which the time step-to-
period ratio is small are integrated accurately, but that the response in the
equations for which At/T is large is necessarily totally inaccurate. Using
explicit conditionally stable methods we are forced to use At smaller than a
critical value which is small enough to obtain accuracy in the integration of
practically all m equations. Using unconditionally stable schemes the time
step can be much larger and should only be small enough that the response
in all modes which contribute significantly to the total structural response is
calculated accurately. The other model components can not be accurate but
the errors are unimportant because the response provided by those
components is negligible. However, it is important to stress again that this
applies only if the integration scheme is unconditionally stable, i.e. if it may
be assumed that the amplitudes in the component which are negligible do
not spuriously grow.
Note that numerical damping can be considered a good feature of implicit
unconditionally stable integration schemes because it may be used to damp
out and practically suppress the response of those modes for which, since
the time step-to-period ratio is large, the response cannot be calculated
accurately anyway. A method with good damping characteristics is, for
example, presented by Hilber et al.
177
RC ELEMENTS UNDER CYCLIC LOADING
section 7.1., that the problem becomes very similar to the corresponding
static problem (compare again equations (197), (197a), (199) and (199a)).
Instead of a tangent stiffness matrix, an 'effective tangent stiffness matrix',
K*, has to be considered, as, for example, in equation (207) or (197a), with
contribution also from the inertia and damping matrices; additional
contribution from these matrices appears also in the expression of the
residue vector R*. Expressions of both K* and R* depend on the particular
implicit method which is used; but it is important to notice that the
contribution of the inertia term in K* is always in the form (c/A/ 2 )M.
When using the same solution strategies, the inertia of the system renders
its dynamic response, in general, smoother than its static counterpart, and
convergence of the iterations can be expected to be quicker. This improved
convergence characteristic in dynamic analysis depends on the contribution
of the mass matrix to the effective tangent stiffness. Indeed this
contribution, being inversely proportional to At2, becomes dominant when
the time step is small.
The first solutions of non-linear dynamic analysis using implicit
integration schemes were performed without equilibrium iterations. It has
been subsequently recognized that the iterations are very important, since
any error allowed in the incremental solution at a particular time directly
affects, in a path-dependent manner, the solution at any subsequent time.
Indeed, because of the high path dependency of any dynamic response, the
analysis of a non-linear dynamic problem requires iterations, at each time
step, more stringently than static analysis does.
178
References
1. Pramono E. and Willam K. Fracture energy-based plasticity formulation of
plain concrete J. EM, Am. Soc. Civ. Engrs, 115, June 1989.
2. Lubliner J. et al. A plastic-damage model for concrete. Int. J. Solids Struct.,
25, No. 3, 1989.
3 Ba2ant Z. P. and Oh B. H. Microplane model for progressive fracture of
concrete and rock. /. EM, Am. Soc. Civ. Engrs, 111, 1985.
4. Bazant Z. P. and Prat P. C. Microplane model for brittle-plastic material: I
theory, II verification. / EM, Am. Soc. Civ. Engrs, 114, Oct. 1988.
5 Chen W. F. Plasticity in reinforced concrete. McGraw-Hill, New York,
1982.
6. Finite element analysis of reinforced concrete. State of the art report.
American Society of Civil Engineers, New York, 1982.
7. Bazant Z. P. and Bhat P. Endochronic theory of inelasticity and failure of
concrete. J. EM, Am. Soc. Civ. Engrs, 102, Aug. 1976.
8. Elwi A. A. and Murray D. W. A 3D hypoelastic concrete constitutive
relationship. J. EM, Am. Soc. Civ. Engrs, 105, Aug. 1979.
9. Buyukozturk O. and Shareef S. S. Constitutive modeling of concrete in
finite element analysis. Computers and Structures, 21, No. 3, 1985.
10. Stankowski T. and Gerstle K. H. Simple formulation of concrete behavior
under multiaxial load histories. ACI J., 82, No. 2, Mar.-Apr. 1985.
11. Shafer G. S. and Ottosen N. S. An invariant-based constitutive model.
Structural Research Series 8506, Dept of Civil Environmental and
Architectural Engineering, University of Colorado, Boulder, 1985.
12. Vermeer P. A. and de Borst R. Non-associated plasticity for soils, concrete
and rocks. Heron Journal, 29, 1984.
13. Han D. J. and Chen W. F. A nonuniform hardening plasticity model for
concrete materials. J. Mech. Mat, 4, 1985.
14. Fardis M. N. and Chen E. S. A cyclic multiaxial model for concrete.
Computational Mechanics, 1, 1986, 301-315.
15. Chen E. S. and Buyukozturk O. Constitutive model for concrete in cyclic
compression. J. EM, Am. Soc. Civ. Engrs, 111, 1985.
16. Yang B. et al. A bounding surface plasticity model for concrete. J. EM, Am.
Soc. Civ. Engrs, 111, Mar. 1985.
17. Comite Euro-International du Beton. Concrete under multiaxial states of
stress. Constitutive equations for practical design. CEB, Lausanne, 1983,
Bulletin d'Information 156.
18. Krajcinovic D. Continuous damage mechanics. Appl. Mech. Rev., 37, 1984.
19. Lemaitre J. How to use damage mechanics. J. Nucl. Eng. & Design, 80,
1984.
20. Resende L. and Martin J. B. A progressive damage continuum model for
granular materials. Comp. Meth. Appl. Mech. Engng, 42, 1984.
21. Krajcinovic D. and Fonseka G. U. The continuous damage theory of brittle
materials. J. Appl. Mech., Am. Soc. Mech. Engrs, 48, 1981.
22. Mazars J. Description of the multiaxial behavior of concrete with an elastic
damaging model. Proc. RILEM Symp. on concrete under multiaxial
conditions. INSA-OPS Toulouse, May 1984.
23. Dougill J. W. On stable progressively fracturing solids. Z. Angew. Math.
Phys., 27, 1976.
24. Gerstle K. H. et al. Behavior of concrete under multiaxial stress states. J.
EM, Am. Soc. Civ. Engrs, 106, Dec. 1980.
25. Scavuzzo R. et al. Stress-strain curves under multiaxial load histories.
179
RC ELEMENTS UNDER CYCLIC LOADING
180
REFERENCES
181
RC ELEMENTS UNDER CYCLIC LOADING
182
REFERENCES
183
RC ELEMENTS UNDER CYCLIC LOADING
184
REFERENCES
185
RC ELEMENTS UNDER CYCLIC LOADING
186
REFERENCES
187
RC ELEMENTS UNDER CYCLIC LOADING
188
REFERENCES
189
RC ELEMENTS UNDER CYCLIC LOADING
190