Finite Elements in Analysis and Design: E. Syroka-Korol, J. Tejchman, Z. Mróz
Finite Elements in Analysis and Design: E. Syroka-Korol, J. Tejchman, Z. Mróz
Finite Elements in Analysis and Design: E. Syroka-Korol, J. Tejchman, Z. Mróz
art ic l e i nf o
a b s t r a c t
Article history:
Received 29 July 2013
Received in revised form
30 April 2014
Accepted 12 May 2014
The paper presents results of FE analysis of mechanical size effects in longitudinally reinforced concrete
slender beams without shear reinforcement failing in shear mode. The simulations were performed
under plane stress conditions for three beams of different sizes and a xed shape (height/length ratio).
The attention was focused on deterministic and statistical size effects related to the nominal beam shear
strength. Concrete was assumed as an isotropic elasto-plastic material exhibiting non-local softening.
The bond strength between concrete and reinforcement was assumed to depend on interface slip with
both stable and softening responses. Statistical simulations were performed for spatially correlated
Gaussian random elds of tensile strength using a stratied sampling reduction method. The FE
numerical results were compared with the respective own experimental test results.
& 2014 Elsevier B.V. All rights reserved.
Keywords:
Elasto-plasticity
Non-local softening
Random elds
Reinforced concrete beams
Size effect
Strain localization
1. Introduction
The size effect phenomenon in quasi-brittle structures is
related to a transition from a ductile behaviour of small specimens
to a totally brittle response of large ones. Thus, the nominal
strength sN decreases with increasing characteristic specimen
dimension D. The reasons for this behaviour are: (a) intense strain
localization regions with a certain volume (i.e. micro-cracked
damage regions called also fracture process zones, FPZ) which
precede discrete macro-cracks; their size related to D contributes
to a deterministic size effect and (b) a spatial variability/randomness of local material properties contributing to a statistical size
effect that becomes dominant with increasing D.
A strong size effect also occurs in reinforced concrete beams
without shear reinforcement wherein diagonal sheartensile fracture takes place in concrete. It was experimentally observed
among others by Leonhardt and Walther [1], Kani [2,3], Bhal [4],
Taylor [5], Walraven [6], Chana [7], Iguro et al. [8], Bazant and
Kazemi [9], Shioya et al. [10], Kim and Park [11], Grimm [12],
Ghannoum [13], Kawano and Watanabe [14], PodgorniakStanik [15], Yoshida [16], Angelakos et al. [17], Lubell et al. [18]
and Syroka-Korol [19]. The diagonal cracks at failure had in
Corresponding author.
E-mail addresses: esyroka@pg.gda.pl (E. Syroka-Korol),
tejchmk@pg.gda.pl (J. Tejchman), zmroz@ippt.gov.pl (Z. Mrz).
http://dx.doi.org/10.1016/j.nel.2014.05.005
0168-874X/& 2014 Elsevier B.V. All rights reserved.
26
Fig. 1. Size effect curve of type II with large cracks or notches by Bazant [20,22]
(N nominal strength, D characteristic specimen size, LEFM linear elastic
fracture mechanics); (a) linear scale and (b) loglog scale.
2. Constitutive models
2.1. Concrete
The concrete deformational response was simulated by assuming an isotropic elasto-plastic constitutive model with a nonlocal softening which was used in our previous calculations
[19,26,3436]. This relatively simple isotropic model for concrete
27
Fig. 4. Bond relationship: radial normal stress sr,rs and bond slip versus radial
normal strain r,rs between concrete and reinforcement during splitting failure [51].
28
5
with
k
Er r;rs; max
sr;rs; max
and
r;rs
:
r;rs; max
;
7
Er E
cef f 0:5b 2 0:252b
where b is the bar diameter, ceff is the effective concrete cover and
E and are the modulus of elasticity of concrete and Poisson's
ratio, respectively.
The limit radial normal stresses and strains were dened as
follows:
cef f 0:88
cef f 1:08
f
sr;rs; max 2f t
; r;rs; max 4:2 t
8
b
E
b
r;rs;res
ft
E
2cef f c0
;
b
b
29
related to the radial normal stress sr,rs by the ctitious friction rule
computed as follows:
tan b ; ;
r:rs
0:5b
10
where b is the cone angle between the bar axis and cone-shaped
cracking surface starting from ribs. The ctitious value was
assumed as b[o] 0.1fc [MPa]. In turn, the bond stress was
sr;rs cot b
11
Cx1 ; x2 i f i x1 f i x2 ;
13
i1
Table 1
Geometry of reinforced concrete slender beams of Figs. 7 and 8.
Dimension Small-size beam
SL20
Medium-size beam
SL40
Large-size beam
SL80
L [mm]
Leff [mm]
H [mm]
D [mm]
a [mm]
b [mm]
3000
2700
400
360
1080
540
6000
5620
800
750
2250
1120
1500
1200
200
160
480
240
The spatially correlated random elds H(x,) (here with the zero
mean and unit variance) were dened according to the Karhunen
Love expansion [54,55] by an innite linear combination of
orthogonal functions with random coefcients
1 p
Hx;
15
i i f i x;
i1
Fig. 7. Cross-section of slender reinforced concrete beams with horizontal bars used in calculations and experiments: (a) small-size SL20, (b) medium-size SL40 and
(c) large-size SL80.
30
Fig. 8. Geometry and loading scheme of reinforced concrete slender beams without stirrups.
Fig. 9. Experimental crack pattern on front (solid lines) and back (dashed lines) side of medium-size reinforced concrete beam.
Fig. 10. FE mesh for small-size beam SL20 (a) and large-size beam SL80 (b) (note that beams are not proportionally scaled).
j; k A N;
17
i x 2j x k;
18
resulting error was minimized by sorting eigenvalues in a decreasing order and controlling the sum of i for i 1,, M.
An analytical solution of Eq. (17) is available with e.g. an
exponential autocorrelation function (Ghanen and Spanos [56]).
In other cases, Eq. (17) has to be solved numerically. In our study,
the wavelet-Galerkin method proposed by Phoon et al [57,58] was
used wherein conventional bases (e.g. trigonometric or polynomials) were replaced by Haar wavelets. A family of the orthogonal
and
f k x di i x T xDk ;
k
19
i0
31
Fig. 12. Contours of non-local tensile softening parameter 1 in small-size reinforced concrete beam SL20 at failure with various bond-slip models: (a) perfect bond, (b) by
Do rr [51], (c) by CEB-FIP [39] and (d) by Akkermann [52] from FE analyses compared to experimental crack pattern (single lines).
4. FE input data
The FE-analyses of longitudinally reinforced concrete slender
beams without shear reinforcement were performed with 3 different beam sizes of a similar geometry from experiments by Syroka
Korol [19]. The beams had the same dimensions H L as in
the experiments: 200 1500 mm2 (small-size beam SL20),
400 3000 mm2 (medium-size beam SL40) and 800 6000 mm2
(large-size beam SL80) (Table 1). All specimens had the constant
thickness t200 mm and constant reinforcement percentage
1% (Fig. 7). The beams were subjected to 4-point bending
(Fig. 8) under the constant shear span to the effective depth ratio
a/D3. The numerical calculations were carried out under plane
stress conditions since experimental crack patterns were very
similar on both sides of all beams (Fig. 9 shows these patterns in
a medium-size beam).
The specimens were cast from concrete C35/45 of a maximum
aggregate size 32 mm (the characteristic compressive strength
fck 35 MPa, the characteristic tensile strength fctk 2.2 MPa,
Young modulus E 34 GPa and Poisson's ratio 0.2). The internal
friction angle was 141 with rs
bc 1.2 (Eq. (A2) and (A4)) and the
dilatancy angle 81 [34]. The calculations were carried out with
the characteristic length lc 5 mm based on the experimentally
measured mean width of localized tensile zones wloc by means of
the DIC technique equal to 20 mm. The non-locality parameter
was m 2 [42]. The ultimate non-local softening parameter in
tension was mainly u1 0.0207 (Gf gf4lc E 180 N/m) and in compression u2 0.0057 (Gc gc4lc E2700 N/m) based on initial calculations. The steel behaviour was specied by the yield stress
fy 500 MPa, the elastic modulus Es 200 GPa and the Poisson's
ratio s 0.3.
The bar diameters were b 10 mm, 16 mm and 20 mm in
the beams SL20, SL40 and SL80, respectively, and b 10 mm for
the second reinforcement layer in the large-size beam SL80. In the
model by Do rr [51] Eq. (7), the slip displacement was u 0.06 mm
(as originally proposed by Do
rr [51]). When using the bond model
by CEB-FIP [39], for the unconned concrete and good bond
conditions case, the material parameters were: 1 0.6 mm,
2 0.6 mm, 3 1 mm, 0.4, max 2.0fck 11.83 MPa and
f 0.15max 1.77 MPa. In the bond model by Akkermann [52]
(Eqs. (5)(11)), we assumed ceff 22 mm, b 0.2 [39] and
c0 0.18 mm (based on own initial FE simulations).
32
Fig. 13. Evolution of non-local tensile softening parameter 1 in small-size reinforced concrete beam SL20 with bond-slip denition by Akkermann [52] from FE analysis at
different vertical resultant force: (a) P 30 kN, (b) P 50 kN, (c) P 85 kN, (d) Pmax 91.4 kN and experimental crack pattern after failure [19] (e). (For interpretation of the
references to color in this gure, the reader is referred to the web version of this article.)
FE meshes of quadrilateral elements composed of four diagonally crossed triangles were used (Fig. 10). The nite elements
had mainly the area of 7.5 10 mm2 (width height) which was
equal to (1.52) lc to obtain mesh-independent results [24,26].
In the case of the large-size beam, the ne mesh covered the beam
mid-part of 4410 mm only. Totally, 130 6001690 344 nite elements
were used depending upon the beam size. The steel bars were
modelled as 2D elements with the width of 1.5 mm and height of
2.02.5 mm depending on the bar diameter b. Since our FE
analyses were two-dimensional, the bar height was taken as
0.5 b to obtain a similar contact surface as in experiments
(2 0.5b). The inuence of the mesh size on FE results was
checked by comparing 2 different meshes in the small-size
reinforced concrete beam SL20 assuming the bond model by
Akkermann [52].
33
Fig. 14. Comparative deterministic FE results with coarse (element size 15 20 mm2) and ne mesh (element size 7.5 10 mm2) for small-size reinforced concrete beam
SL20: (A) load-deection diagram and (B) contours of non-local tensile softening parameter 1 at failure (a) coarse mesh and (b) ne mesh.
Fig. 15. Evolution of diagonal sheartensile localized zone with marked points: (a) at P 76 kN, (b) at P 86 kN, (c) at peak and (d) at failure (small-size reinforced concrete beam).
mesh are similar (Fig. 14A and B). The failure force for the coarse
mesh (element size 15 20 mm2) was 93 kN while for the ne
mesh (element size 7.5 10 mm2) was 91 kN (Fig. 14A). Based on
these results, we assumed that there was no need to decrease
more the size of nite elements (less than 7.5 10 mm2).
34
35
Fig. 17. Contours of non-local tensile softening parameter 1 (a) compared to experimental crack pattern (b) and its distribution along beam bottom (c) and along mid-height
(d) from deterministic simulations with small-size reinforced concrete beam.
and
0:5f
nn
"r
#
"r
#
a x2 a x
a x2 a x
1
1
0:5s c f ck
D
D
D
D
20
21
36
Fig. 18. Contours of non-local tensile softening parameter 1 (a) compared to experimental crack pattern (b) and its distribution along beam bottom (c) and along mid-height
(d) from deterministic simulations with reinforced concrete medium-size beam.
Fig. 19. Contours of non-local tensile softening parameter 1 (a) compared to experimental crack pattern (b) and its distribution of along beam bottom (c) and along midheight (d) from deterministic simulations with large-size beam.
37
Table 2
Calculated normal un and tangential ut displacement increments between points of Fig. 15 across diagonal localized zone during its evolution (small-size reinforced
concrete beam).
Vertical force
u [mm]
Distance ab
Distance cd
Distance ef
Distance gh
Distance ik
Distance ln
Distance mo
Distance pr
76 kN
un
ut
0.01
0.12
0.01
0.12
0.02
0.13
0.04
0.13
0.09
0.11
0.01
0.09
0.02
0.12
0.15
0.12
86 kN
un
ut
0.03
0.15
0.05
0.15
0.10
0.17
0.15
0.18
0.22
0.14
0.11
0.07
0.04
0.13
0.27
0.12
At peak
un
ut
0.09
0.16
0.14
0.17
0.23
0.21
0.31
0.22
0.40
0.18
0.29
0.01
0.16
0.11
0.40
0.09
After failure
un
ut
0.16
0.19
0.26
0.19
0.39
0.25
0.49
0.28
0.61
0.20
0.46
0.06
0.29
0.11
0.48
0.07
Fig. 20. Contours of non-local tensile softening parameter 1 from deterministic simulations compared to experimental crack pattern (lines) with: (a) small-size,
(b) medium-size and (c) large-size beam.
Fig. 21. Exemplary random elds of tensile strength ft in: (a) small-size, (b) medium-size and (c) large-size beam used in FE analyses (note that beams are not proportionally
scaled).
with
3:5
1
c 1:6p0:27 1 p 0:15 0:58; ;
D
f ck
22
38
Considering statistical simulations, all mean values underestimate the experimental shear strength by 10% on average (this
result is probably caused by the assumption of a deterministic
bond model in statistical computations or a sampling region
limited to beam shear regions). The mean shear strength reduction
in statistical computations with respect to deterministic results
was close to the assumed coefcient of variation covft of the tensile
strength. Substituting into the experimentally calibrated SEL formula (Eq. (1)) the reduced value of o, expressed by o(1 covft)
2.87(1 0.12) 2.53 MPa, a perfect match between the new SEL
curve and the statistical results was obtained (Fig. 27). Therein, the
brittleness number dened by the ratio D/D0 was not modied.
This indicates that the size effect on the nominal shear strength is
predominantly of a deterministic type.
With respect to FE numerical results on a size effect by Sato
et al. [31], our deterministic results indicate a stronger reduction of
the nominal shear strength with similar increasing beam size (38%
between the small- and large-size beam against 20%).
7. Concluding remarks
The following conclusions can be drawn from our plane stress
non-linear FE analyses with a constant and spatially correlated
random distribution of the local tensile strength in geometrically
similar longitudinally reinforced concrete beams without shear
reinforcement failing in shear:
Fig. 22. Vertical force P deection u diagrams (12 curves) for reinforced concrete
beams from statistical FE simulations (solid lines) compared to one deterministic
curve (dashed line): (a) small-size beam, (b) medium-size beam, and (c) largesize beam.
local softening enhanced by a characteristic length of microstructure was able to realistically reproduce both a very
complex experimental failure mode in beams (combined
diagonal sheartensile failure and horizontal bond splitting
failure) and a pronounced reduction of the nominal shear
strength with increasing beam size by taking a bond model
into account.
The size effect on the nominal shear strength of slender RC
beams without shear reinforcement was of a deterministic type
only. The nominal shear strength from deterministic simulations strongly decreased with increasing specimen size. The
deterministic size effect was caused by localized zones with a
constant width and a linearly varying height. It was in agreement with the energetic size effect by Bazant. Thus, it was not
caused by non-symmetric cracking combined with the unintended out-of-plane beam action. The mean shear strengths in
statistical computations were solely lower by 10% than corresponding deterministic outcomes. The mean statistical shear
strength was proportional to the mean tensile strength reduced
by the standard deviation.
The effect of the bond denition between concrete and reinforcement was pronounced and led to the different failure
mode, load capacity, beam stiffness and number of localized
zones. The application of the bond model by Akkermann [52]
contributed to a sheardiagonal tensile failure mechanism in
numerical simulations as in experiments.
A position of a critical diagonal localized zone in deterministic
calculations was close to experimental outcomes. In statistical
calculations, a critical diagonal localized zone was located at a
higher distance from the beam support, hence the mean
bearing capacity was lower than the corresponding deterministic one. The height of a cracked section in experiment was
successfully reproduced by both deterministic and statistical
analyses. The height of a critical diagonal localized zone was
slightly larger in statistical calculations than in deterministic
ones. The width of localized zones was similar in deterministic
and statistical calculations.
39
Fig. 23. Contours of non-local tensile softening parameter 1 from statistical simulations of small-size reinforced concrete beam corresponding to: (a) highest failure force,
(b) mean failure force, and (c) lowest failure force.
Fig. 24. Contours of non-local tensile softening parameter 1 from statistical simulations of medium-size reinforced concrete beam corresponding to: (a) highest failure
force, (b) mean failure force, and (c) lowest failure force.
Fig. 25. Contours of non-local tensile softening parameter 1 from statistical simulations of large-size reinforced concrete beam corresponding to: (a) highest failure force,
(b) mean failure force, and (c) lowest failure force.
40
1
f 2 sij ; 2 q p tan c 2 q p tan 1 tan sc 2 r 0
3
A2
r
3
s s ;
q
2 ij ij
tan
Fig. 26. FE results of nominal shear strength 1.5 V/(tD) against effective beam
depth D from deterministic analyses for reinforced concrete beams without shear
reinforcement as compared with experiments, size effect law (Eq. (1)) and upper
bound plastic theory results (Eqs. (20) and (21)).
1
p skk ;
3
A3
31 r sbc
;
1 2r sbc
A4
g1 f 1
A5
g 2 q p tan
A6
st 1 f t 1 A1 1 3 exp A2 1 A3 1
A7
with
A1
Fig. 27. Mean statistical nominal shear strength 1.5 V/(tD) against effective
beam depth D from FE analyses for reinforced concrete beams without shear
reinforcement compared to modied size effect law SEL of type II (Eq. (1)) by
expression o(1-covft).
Acknowledgements
The research work has been carried out within the project:
Innovative ways and effective methods of safety improvement
and durability of buildings and transport infrastructure in the
sustainable development nanced by the European Union
(POIG.01.01.02-10-106/09-01) and the project Experimental and
numerical analysis of coupled deterministic-statistical size effect
in brittle materials nanced by National Research Centre NCN
(UMO-2013/09/B/ST8/03598).
Appendix
In order to describe the concrete behaviour, two failure criteria
were assumed. In tension regime, the Rankine criterion was used
with the yield function f1 using isotropic softening and associated
ow rule and in compression, the DruckerPrager yield surface f2
with isotropic hardening/softening and non-associated ow rule
was used [36]
f 1 si ; 1 maxfs1 ; s2 ; s3 ;g st 1 r 0;
A1
c1
;
u
A2
c2
u
and
A3
1
1 c31 exp c2 ;
u
A8
A10
[3] G.N.J. Kani, How safe are our large reinforced concrete beams? ACI J. 64 (1967)
128142.
[4] N.S. Bhal, ber den Einuss der Balkenho
he auf die Schubtragfhigkeit von
einfeldrigen Stahlbetonbalken mit und ohne Schubbewehrung (Dissertation),
Stuttgart Universitt, 1968.
[5] H.P.J. Taylor, Shear strength of large beams, J. Struct. Div., ASCE 98 (1972)
24732489.
[6] J.C. Walraven, The Inuence of Depth on the Shear Strength of Lightweight
Concrete Beams without Shear Reinforcement, Stevin Laboratory Report no.
5-78-4, Delft University of Technology, 1978.
[7] P.S. Chana, Some aspects of modelling the behavior of reinforced concrete
under shear loading, Tech. Rep., Cement and Concrete Assoc., Wexham
Springs, 1981.
[8] M. Iguro, T. Shioya, Y. Nojiri, H. Akiyama, Experimental studies on shear
strength of large reinforced concrete beams under uniformly distributed load,
Concr. Library Int. JSCE 5 (1985) 137146.
[9] Z. Bazant, M.T. Kazemi, Size effect on diagonal shear failure of beams without
stirrups, ACI Struct. J. 88 (1991) 268276.
[10] T. Shioya, H. Akiyama, Application to design of size effect in reinforced
concrete structures, in: H. Mihashi, H. Okamura, Z.P. Bazant (Eds.), Size Effect
in Concrete Structures, Spon, London, 1994, pp. 409416.
[11] J.K. Kim, Y.D. Park, Shear strength of reinforced high-strength concrete beams
without web reinforcement, Mag. Concr. Res. 46 (1994) 716.
[12] R. Grimm, Einu bruchmechanischer Kenngro
en auf das Biege- und
Schubtragverhalten hochfester Betone (Ph.D. thesis), Fachb. Konstr. Ingenieurbau der TH Darmstadt und DafStb H.477, Beuth Verlag GmbH, Berlin, 1997.
[13] W.M. Ghannoum, Size effect on shear strength of reinforced concrete beams
(M.A.Sc thesis), MC Gill University, Department of Civil Engineering and
Applied Mechanics, Montreal, Canada, 1998.
[14] H. Kawano, H. Watanabe, Shear strength of reinforced concrete columns
effect of specimen size and load reversal, in: Proceedings of the 2nd Italy
Japan Workshop on Seismic Design and Retrot of Bridges, Rome, Italy, 1997,
pp. 141154.
[15] B.A. Podgorniak-Stanik, The inuence of concrete strength, distribution of
longitudinal reinforcement, amount of transverse reinforcement and member
size on shear strength of reinforced concrete members (M.A.Sc. thesis), Dept.
of Civil Engineering, University of Toronto, 1998.
[16] Y. Yoshida, Shear reinforcement for large lightly reinforced concrete members
(M.A.Sc. thesis), Department of Civil Engineering, University of Toronto, 2000.
[17] D. Angelakos, E.C. Bentz, M.P. Collins, Effect of concrete strength and minimum
stirrups on shear strength of large members, ACI Struct. J. 98 (2001) 290300.
[18] A. Lubell, T. Sherwood, E. Bentz, M. Collins, Safe shear design of large, wide
beams, Concr. Int. 26 (2004) 6678.
[19] E. Syroka-Korol, Theoretical and experimental study on size effect in concrete
beams reinforced with steel or basalt bars (Ph.D. thesis), Gdansk University of
Technology, 2012.
[20] Z. Bazant, Size effect in blunt fracture: concrete, rock, metal, J. Eng. Mech.
ASCE, 110, 1984518535.
[21] Z. Bazant, Q. Yu, Designing against size effect on shear strength of reinforced
concrete beams without stirrups: I. Formulation, J. Struct. Eng. 131 (2005)
18771885.
[22] Z. Bazant, J. Planas, Fracture and Size Effect in Concrete and Other Quasi-brittle
Materials, CRC Press LLC, Boca Raton, US, 1998.
[23] Q. Yu, Size effect and design safety in concrete structures under shear (Ph.D.
thesis), Northwestern University, 2007.
[24] L. Skarynski, E. Syroka, J. Tejchman, Measurements and calculations of the
width of the fracture process zones on the surface of notched concrete beams,
Strain 47 (s1) (2011) e319e332.
[25] J. Bobiski, J. Tejchman, J. Grski, Notched concrete beams under bending
-calculations of size effects within stochastic elasto-plasticity within non-local
softening, Arch. Mech. 61 (2009) 125.
[26] E. Syroka-Korol, J. Tejchman, Z. Mrz, FE calculations of a deterministic and
statistical size effect in concrete under bending within stochastic elastoplasticity and non-local softening, Eng. Struct. 48 (2013) 205219.
[27] E. Syroka-Korol, J. Tejchman, Experimental investigations of size effect in
reinforced concrete beams failing by shear, Eng. Struct. 58 (2013) (2014)
6378.
[28] M.D. Kotsovos, M.N. Pavlovic, A possible explanation for size effects in
structural concrete, Arch. Civ. Eng., Pol. Acad. Sci. XL (1994) 243261.
[29] M.D. Kotsovos, M.N. Pavlovic, Size effects in structural concrete: a numerical
experiment, Comput. Struct. 64 (1997) 285295.
[30] F.J. Vecchio, W. Shim, Experimental and analytical re-examination of classic
concrete beam testsJ. Struct. Eng., ASCE (2004)460469.
[31] Y. Sato, T. Tadokoro, T. Ueda, Diagonal tensile mechanism of reinforced
concrete beams, J. Adv. Concr. Technol. 2 (2004) 327341.
[32] A. Slobbe, A.M. Hendriks, J. Rots, Sequentially linear analysis of shear critical
reinforced concrete beams without shear reinforcement, Finite Elem. Anal.
Des. 50 (2012) 108124.
[33] Q. Yu, Z. Bazant, Can stirrups suppress the size effect on shear strength of RC
beams? J. Struct. Eng., ASCE 137 (2011) 607617.
41