Buildings 12 02092
Buildings 12 02092
Buildings 12 02092
Article
Flat Slabs in Eccentric Punching Shear: Experimental Database
and Code Analysis
Daniel Vargas 1 , Eva O. L. Lantsoght 1,2, * and Aikaterini S. Genikomsou 3
Abstract: Eccentric punching shear can occur in concrete slab–column connections when the connec-
tion is subjected to shear and unbalanced moments. Unbalanced moments occur in all floor slabs
at the edge and corner columns. As such, this problem is of practical relevance. However, most
punching experiments in the literature deal with concentric punching shear at internal columns.
This paper presents a developed database of 128 experiments of flat slabs under eccentric punching
shear, including a summary of the testing procedure of each reference and a description of the slab
specimens. Additionally, a linear finite element analysis of all the specimens is included to determine
the relevant sectional shear forces and moments. Finally, the ultimate shear stresses from the database
experiments are compared to the shear capacities determined with ACI 318-19, Eurocode 2 and
the Model Code 2010. The comparison shows that the Model Code 2010 is the most precise in the
predictions with an average tested to predicted ratio of 0.82 and a coefficient of variation of 29.63%.
It can be concluded that improvements to the current design methods for eccentric punching shear
Citation: Vargas, D.; Lantsoght, are necessary.
E.O.L.; Genikomsou, A.S. Flat Slabs
in Eccentric Punching Shear: Keywords: code provisions; database; eccentric punching shear; experiments; flat slab; linear finite
Experimental Database and Code
element models; punching; reinforced concrete; shear; unbalanced moments
Analysis. Buildings 2022, 12, 2092.
https://doi.org/10.3390/
buildings12122092
The first comprehensive studies on punching shear were performed in the 1960s by
Kinnunen and Nylander [5], but their mechanical models resulted in complicated expres-
sions, which code makers found impractical for use [6]. Instead, empirical expressions
based on the available test results were created for the development of the code provisions.
Given that there is a lack of experimental information on eccentric punching shear on large-
scale flat slabs, it became difficult to provide a satisfactory design expression [2]. To account
for the eccentric loading, ACI 318-19, Eurocode 2 EN 1992-1-1:2005 and the fib Model Code
2010 models considered the shear stress distribution on the critical perimeter [7–9], assum-
ing either a linear or plastic stress distribution. The punching perimeter is at a distance d/2
from the column in ACI 318-19, and it is at the same distance, but with rounded corners,
in the fib Model Code 2010. Eurocode 2 uses a perimeter at a distance 2d from the column,
and the perimeter has rounded corners as well.
Despite the efforts undertaken by investigators through the years, the current design
methods cannot accurately predict the punching shear strengths when unbalanced mo-
ments act on slabs. Nowadays, the advances in materials and new analysis methodologies,
such as nonlinear finite element analysis and better instrumentation techniques for ex-
perimental campaigns, have helped researchers start proposing a reshaping of the design
codes to best meet the real performance of the slab–column connections under eccentric
loading [10].
This work aims to present a wider view of the problem by compiling and analyzing
experiments on eccentric punching shear from the literature. The analysis of the compiled
experiments can be used to examine the performance of the currently available build-
ing codes and identify which types of experiments would be a valuable contribution to
the body of knowledge. Additional experiments could be used to refine and improve
the existing models. In addition, the developed database can serve those who are work-
ing on mechanical models of punching shear to check the performance for the case of
eccentric punching.
This article compiles 128 experiments on flat slabs in eccentric punching shear. Vertical,
horizontal and combined loading setups are reported in the literature. Both slabs with and
without shear reinforcement are included in the developed database. The internal forces
of the slabs for the maximum applied load, i.e., at the onset of punching shear failure, are
typically not available in the references. To complete the missing information, a linear finite
element model of each experiment is constructed. The experimental shear capacities from
the database are then compared to the strengths predicted by the design expressions found
in ACI 318-19 [7], Eurocode 2 NEN-EN 1992-1-1:2005 [8] and the fib Model Code 2010 [9].
2. Methods
2.1. Overview of Code Provisions
2.1.1. ACI-318-19
The punching shear provisions from ACI 318-19 are based on empirical equations
derived from test results by Moe [11] and analyses by the ACI-ASCE Committee 426 [12].
The ACI 318-19 method is based on the maximum shear stress vu on the critical perimeter
bo of the slab, which is located at 0.5 d from the face of the column, where d is the average
slab effective depth. The maximum shear stress vu should not exceed the nominal shear
strength of the slab vn . Figure 1 is a sketch of the shear stresses produced by axial load and
moment transfer on an internal slab–column connection [1].
Buildings 2022, 12, 2092 3 of 49
Figure 1. Shear stress produced by applied load and moment transfer, modified from Ref. [1]:
(a) transfer of unbalanced moments to column; (b) shear stress caused by direct shear; (c) shear stress
caused by unbalanced moments; (d) total shear stress: sum of (b,c).
Vu γv Mu c
vu = ± (1)
bo d Jc
where Vu is the factored shear being transferred from the slab to the column acting on the
centroid of the critical section; c is the distance from the centroid of the critical section to
Buildings 2022, 12, 2092 4 of 49
the point where the shear stress is calculated; Jc is the polar moment of inertia of the critical
section; and γv Mu is the fraction of moment transferred by the shear, with γv as follows:
γv = 1 − γ f (2)
1
γf = q b1 (3)
2
1+ 3 b2
where b1 is the total width of the critical section measured perpendicular to the axis around
which the moment acts, and b2 is the total width parallel to the axis [1]. Figure 2 shows a
sketch of the critical perimeter of an interior, edge and corner slab–column connection.
Figure 2. Critical perimeter of an interior, edge and corner slab–column connections, modified
from Ref. [1]: (a) interior slab–column connection; (b) edge slab–column connection; (c) corner
slab–column connection.
Buildings 2022, 12, 2092 5 of 49
Av f yt
vs = (8)
bo s
where Av is the sum of the area of all legs of reinforcement on the peripheral line, which
is geometrically like the perimeter of the column section; fyt is the yield strength of the
transverse reinforcement; and s is the spacing of transversal reinforcement [7].
Section 22.6.6.1 [7] indicates that the value of vc for shear-reinforced slabs shall not
exceed the following: p
vc = 0.17λs λ f c0 (9)
p
vc = 0.25λs λ f c0 (10)
0.33 p
vc = 0.17 + λs λ f c0 (11)
β
0.083αs d p
vc = 0.17 + λs λ f c0 (12)
bo
Equation (9) is used for stirrup reinforcement, and Equations (10)–(12) are used for
headed shear stud reinforcement, where the least of them shall be taken. When shear
reinforcement is used, the critical perimeter bo shall be taken outside the reinforced section,
as illustrated in Figure 3 [7].
Figure 3. Critical perimeter for a shear-reinforced interior, edge and corner slab–column connection,
modified from Ref. [7]: (a) interior slab–column connection; (b) edge slab–column connection;
(c) corner slab–column connection.
Buildings 2022, 12, 2092 7 of 49
Figure 4. Basic control perimeter for an interior, edge and corner slab–column connection, modified
from Ref. [8]: (a) interior slab–column connection; (b) edge slab–column connection; (c) corner
slab–column connection.
Punching shear is evaluated based on the following stresses: vRd,c —the design value
of the punching shear resistance of a slab without punching shear reinforcement; vRd,s —the
value of the punching shear resistance of a slab with punching shear reinforcement; and
vEd —the maximum acting shear stress along the control section. If vEd ≤ vRd,c , then punch-
ing shear reinforcement is not necessary. If the support reaction is eccentric with respect to
the control perimeter, the maximum shear stress is
VEd
v Ed = β EC (13)
U1 d
Buildings 2022, 12, 2092 8 of 49
MEd U1
β EC = 1 + k c (14)
VEd W1
where W 1 represents the shear distribution on the control perimeter; VEd is the design value
of the sectional shear force; MEd is the design value of the sectional bending moment; and kc
is a coefficient on the ratio between the column dimensions given by Table 6.1 of NEN-EN
1992-1-1:2005 [8]. A few values of kc are 0.6 for a c1 /c2 ratio of 1.0 and 0.70 for a c1 /c2 ratio
of 2.0, where c1 and c2 are the dimensions of the column (see Figure 5). W 1 is calculated as
Z U
i
W1 = |e|dl (15)
0
where Ui is the length of the control perimeter under consideration; dl is a length increment
of the perimeter; and e is the distance of dl from the axis around which the moment MEd
acts [8]. Figure 5 shows shear distribution due to an unbalanced moment at a slab–column
connection, indicating that the Eurocode approach assumes a fully plastic distribution of
the shear stresses.
For an internal rectangular column where the loading is eccentric to both orthogonal
axes, βEC shall be calculated as follows:
s 2 2
ey
ex
β EC = 1 + 1.8 + (16)
bx by
where ey and ex are the eccentricities MEd /VEd along the axes y and x, respectively, and bx
and by are the dimensions of the control perimeter. For edge slab–column connections,
where the eccentricity is perpendicular to the slab edge toward the interior, and there
is no eccentricity parallel to the edge, the control perimeter may be reduced to U1 *, as
illustrated in Figure 6a. For corner slab–column connections, where the eccentricity is
toward the interior of the slab, the control perimeter may be reduced to U1 *, as illustrated
in Figure 6b [8]. This approach was adopted for safety reasons, considering that when
there is a moment around an axis parallel to the slab edge or a moment at a corner column,
the experimental results showed that punching failure is typically preceded by torsional
cracking at the edge of the slab [14].
Buildings 2022, 12, 2092 9 of 49
Figure 6. Reduced basic control perimeter, modified from Ref. [8]: (a) edge slab–column connection;
(b) corner slab–column connection.
For edge slab–column connections, if there are eccentricities in both orthogonal direc-
tions, βEC shall be calculated as
U1 U
β EC = + k c 1 e par (17)
U1∗ W1
where epar is the eccentricity parallel to the slab edge. For edge and corner column connec-
tions, where the eccentricity is toward the interior of the slab, βEC shall be calculated as
U1
β EC = (18)
U1∗
If the eccentricity is toward the exterior, βEC shall be calculated using Equation (16).
The punching shear resistance of slabs without shear reinforcement vRd,c is calcu-
lated as 1
v Rd,c = CRd,c k (100ρl f ck ) 3 ≥ vmin (19)
with vRd,c taken as 0.18/γc , with γc the material factor for concrete (γc = 1.5) and k the size
effect factor, calculated with the following expression, with d in [mm]
r
200
k = 1+ ≤2 (20)
d
Buildings 2022, 12, 2092 10 of 49
The reinforcement ratio is the geometric average of the reinforcement ratio in the y (ρly )
and x (ρlx ) direction:
p
ρl = ρlx .ρly (21)
The lower bound of the shear capacity is a nationally determined parameter, with a
recommended expression for vmin as
1/2
vmin = 0.035k3/2 f ck (22)
where Asw is the area of one perimeter of shear reinforcement around the column; sr is the
radial spacing of perimeters of shear reinforcement; fywd,ef is the effective design strength of
the punching shear reinforcement; and α is the angle between shear reinforcement and the
horizontal plane of the slab.
Figure 7. Basic control perimeter, modified from Ref. [9]: (a) interior column; (b) edge slab–
column connection.
Then, for calculating the punching shear resistance of the slab, a control perimeter b0
is used. This perimeter accounts for the non-uniform distribution of shear forces along
b1,MC , which can be caused by concentrations of the shear forces due to moment transfer
between the slab and the supported area because of eccentricities in the load application [9].
Figure 8 illustrates the eccentricity of the resultants [9].
Buildings 2022, 12, 2092 11 of 49
b0 = k e b1,MC (24)
1
ke = (25)
1 + beuu
where eu is the eccentricity of the resultant shear forces with respect to the centroid of b1,MC ,
and bu is the diameter of a circle with the same area as the region inside b1,MC .
The punching shear resistance VRd is calculated as
The design shear resistance attributed to the concrete is calculated using the following
expression, with the compressive strength of the concrete, fck , in [MPa]:
p
f ck
VRd,c = k ψ b0 dv (27)
γc
kψ is a parameter that depends on the rotations of the slab and shall be calculated as
1
kψ = ≤ 0.6 (28)
1.5 + 0.9k dg ψd
where d is the mean value of the effective depth of the slab for x and y directions, and kdg
shall be calculated as follows, with dg in [mm]:
32
k dg = ≥ 1.15 (29)
16 + d g
Buildings 2022, 12, 2092 12 of 49
where ∑Asw is the sum of the area of all the shear reinforcement acting on the zone between
0.35dv and dv , which has a length of 0.65dv (see Figure 9) [9].
Es ψ
σswd = ≤ f ywd (31)
6
where fywd is the yield strength.
The load-rotation behavior of the slab is calculated as follows:
rs f yd msd 1.5
ψ = 1.5 (32)
d Es m Rd
where rs is the distance from the column axis to the line of contra-flexure of the radial
bending moments; fyd is the yield strength of the flexural reinforcement; Es is the modulus
of elasticity of the flexural steel; msd is the average moment per unit length for calculating
flexural reinforcement in the support strip; and mRd is the average flexural strength per
unit length in the support strip [9]. The values of the mechanical parameters in the formula
can be calculated with different levels of approximation (LoA), where increasing levels of
approximation indicate increasing precision but also increasing computational time and
effort [6].
LoA I assumes that msd = mRd , which implies that the strength of the slab will be
governed by its bending moment capacity. For regular slabs with a long-to-short span
length ratio 0.5 ≤ Lx /Ly ≤ 2.0, rs can be estimated as follows:
Figure 10. Slab dimensions and support strip dimensions, modified from Ref. [9].
Figure 11. Slab geometries and test slab–column connection: (a) square interior slab–column connec-
tion [19]; (b) rectangular edge slab–column connection [17]; (c) square corner slab–column connec-
tion [18].
For Refs. [2,17,18,22,24,25,28,29], the age of the specimens at the time of testing was not
provided, and thus, it is assumed to be 28 days. For Refs. [2,11,18,21–23,25–30], the tensile
strength of the concrete fct was not reported by the authors. To complete this information,
the expression developed by Sarveghadi [32] was used:
p
f ct = 0.76 f c (34)
For internal slab–column connections, Refs. [2,11,19] presented the ultimate load ap-
plied to the slab–column connection and its eccentricity; on the other hand, Refs. [21–24]
reported the ultimate moment applied to the slab–column connection. For the edge
slab–column connections, Refs. [17,19,22] presented the ultimate load applied to the slab–
column connection and its eccentricity, and Refs. [20,23,25–27] reported the ultimate mo-
ment applied to the slab–column connection. Finally, for corner slab–column connections,
Refs. [18,20,26,30] reported the ultimate moment applied to the slab–column connection.
In all the works, the test setup caused this moment to act diagonally on the slab. Figure 11c
illustrates this type of loading. On the other hand, Refs. [28,29] reported the ultimate load
applied to the slab–column connection and its eccentricity. For the database, the diagonally
applied moment was divided into its components in the x and y directions. All values
in the database are presented in SI units. The information from Refs. [11,20–23,30] was
converted from US customary to SI units.
Figure 13. Example of finite element analysis on an Albuquerque [17] specimen: (a) Applied loads
in the model and support conditions (side view); (b) Rigid line links on the column x and y axis,
punching perimeter defined by ACI 318-19 code [7] (plan view); (c) Mesh generation (diagonal view);
(d) Internal shear stress on the critical perimeter, calculated from the applied load and the self-weight
of the specimen.
The internal, edge and corner specimens followed the same modeling process. Figure 14
shows the typical shear stress distributions around the control perimeter considered; note
that the loads applied are displayed to show the loading eccentricity.
Figure 14. Example of typical shear stress distributions around the control perimeter defined by ACI
318-19 code [7]: (a) Internal slab–column connection; (b) Edge slab–column connection; (c) Corner
slab–column connection.
Figure 15 shows the distribution of the most important parameters in the database.
Figure 15a shows that the majority of the slabs are made of normal strength concrete.
The developed database cannot be used to gain insight into the eccentric punching shear
capacity of high-strength concrete slab–column connections. Figure 15b shows that a tensile
reinforcement ratio in the range of 1.25–1.50% was commonly used in the tested slabs.
Typical slab designs use reinforcement ratios of 0.6–0.8%. None of the experiments in
the database used these practical values, with most slabs being over-reinforced in flexure
to achieve a punching shear failure. The distribution of the average effective depth of
the slabs is presented in Figure 15c. This plot shows that at least half of the specimens
had an effective depth d in the range 100 mm–125 mm, and another large portion of the
specimens had an effective depth close to 75 mm. The reported specimens are small-scale
specimens that do not give us insights regarding the size effect for eccentric punching shear.
Figure 15d shows the ratio between the shear span and the average effective depth a/d.
The range of a/d in the experiments covers only situations in which no direct load transfer
can occur; as such, for this database, there is a consistency in the range of a/d. Figure 15e
shows the maximum aggregate size values reported in the literature. The values reported
are consistent with the values shown in Figure 15c; relatively small maximum aggregate
sizes are used for the fabrication of specimens with small depths.
Figure 15. Distribution of the most important parameters in the database: (a) concrete compressive
strength fc ; (b) tensile reinforcement ratio ρl ; (c) effective depth d; (d) shear span to average effective
depth ratio a/d; (e) maximum aggregate size dg .
3. Results
3.1. Parameters Studied
The raw data from the database are used to analyze the effect of different experi-
mental parameters on the sectional shear stress at failure. The ACI 318-19 [7] expression
(Equation (1)) is used for determining the shear stress on the perimeter vu . Normalized
shear stresses are used to discard the influence of the concrete compressive strength fc .
An analysis of the shear stress normalized to the square root and to the cubic root of the
concrete compressive strength is carried out first. Figure 16a,b show the relation between
the normalized shear strength and fc , and, as can be seen, for the experimental results
Buildings 2022, 12, 2092 21 of 49
studied, normalizing the shear strength to the square root of the concrete compressive
strength is preferable. A similar observation was made for the shear capacity of steel fiber
reinforced concrete beams [34].
Figure 16. Shear stresses normalized to the concrete compressive strength: (a) normalized to the
square root; (b) normalized to the cubic root.
Figure 17. Parameter studies based on the normalized shear stress at failure of all entries in the
database: (a) effective depth d; (b) longitudinal reinforcement ratio ρl ; (c) shear span to depth
ratio a/d.
in Tables 2–4. Figure 19 shows the comparison between the SCIA Engineer [33] FEM
results vFEM and the predicted shear capacities vpred for ACI 318-19 [7], with the statistical
properties of this comparison in Table 5. The FEM results of the shear stress were compared
only to ACI 318-19 [7], as NEN-EN 1992-1-1:2005 [8] and fib Model Code 2010 [9] assume a
plastic stress distribution on the punching perimeter. ACI 318-19 [7], on the other hand,
assumes a linear stress distribution. The results of the linear finite element analysis were
compared to check the alignment of the assumptions of linear behavior. The results for
all the entries of the database are presented in Tables A8–A10. Some experiments only
use the moment on the slab-column connection and do not use a load on the slab. For
these referenced works, the NEN-EN 1992-1-1:2005 [8] and fib Model Code 2010 [9] models
were not evaluated. Equation (14) from NEN-EN 1992-1-1:2005 [8] uses the value of shear
force applied to the slab–column connection for calculating the enhancement factor for
eccentric shear, βEC , so that the shear stress caused by unbalanced moment only cannot
be determined. The same problem arises when applying Equations (24) and (25) from the
fib Model Code 2010 [9]; the eccentricity eu is calculated from the resultant shear forces
applied to the slab–column connection.
Figure 18. Comparison between experimental vtest and predicted shear capacities vpred for three
design methods from existing codes.
The validation of the spreadsheet used for calculating the code predictions is available
in the public domain [35].
In Table 2, the statistical results are first presented in general and are then presented
separately by the type of slab–column connection being analyzed. The number of experi-
ments evaluated varies from code to code because specimens that were not tested under
direct shear could not be evaluated by NEN-EN 1992-1-1:2005 [8] and Model Code 2010 [9].
Buildings 2022, 12, 2092 24 of 49
Table 2. Statistical results from the comparison between the tested and predicted capacities, Part I.
The number of specimens used for the evaluation is shown within brackets. First, all results are shown
together. Then, the results are subdivided into interior slab–column connections with unbalanced
moment, edge slab–column connections and corner slab–column connections.
All Results
AVG STD COV (%)
ACI (128) 1.65 0.58 35.39
EC2 (122) 1.52 0.69 45.38
MC2010 (122) 0.82 0.24 29.63
Internal slab–column connections
AVG STD COV (%)
ACI (37) 1.41 0.31 22.27
EC2 (37) 1.15 0.22 19.00
MC2010 (37) 0.81 0.17 20.47
Edge slab–column connections
AVG STD COV (%)
ACI (55) 1.79 0.61 33.78
EC2 (51) 1.70 0.88 51.93
MC2010 (51) 0.79 0.21 26.20
Corner slab–column connections
AVG STD COV (%)
ACI (36) 1.67 0.69 41.18
EC2 (34) 1.66 0.54 32.75
MC2010 (34) 0.88 0.34 38.70
Table 3. Statistical results from the comparison between the tested and predicted capacities, Part
II. The number of specimens used for the evaluation is shown within brackets. First, all results
are shown together. Then, the results are subdivided into interior slab–column connections with
unbalanced moment, edge slab–column connections and corner slab–column connections.
Table 4. Statistical results from the comparison between the tested and predicted capacities, Part
III. The number of specimens used for the evaluation is shown within brackets. First, all results
are shown together. Then, the results are subdivided into interior slab–column connections with
unbalanced moment, edge slab–column connections and corner slab–column connections.
Shear-Reinforced Slabs
AVG STD COV (%)
ACI (18) 1.64 0.61 37.32
EC2 (18) 1.59 0.75 47.54
MC2010 (18) 0.93 0.20 21.35
Internal slab–column connections
AVG STD COV (%)
ACI (9) 1.30 0.40 30.64
EC2 (9) 1.18 0.21 17.65
MC2010 (9) 0.93 0.16 17.14
Edge slab–column connections
1.89 0.57 30.02
ACI (7) 2.10 0.99 46.99
EC2 (7) 0.83 0.14 17.27
MC2010 (7) 1.89 0.57 30.02
Corner slab–column connections
AVG STD COV (%)
ACI (2) 2.29 0.90 39.09
EC2 (2) 1.61 0.29 18.30
MC2010 (2) 1.25 0.26 20.88
Figure 19. Comparison between the SCIA Engineer [33] FEM results and the predicted shear capaci-
ties for ACI 318-19 [7].
Buildings 2022, 12, 2092 26 of 49
Table 5. Statistical results from the comparison between SCIA Engineer [33] FEM results for the
acting shear stress and the predicted shear capacities for ACI 318-19 [7].
In Table 3, the statistical results are first presented in general and are then presented
separately by the type of slab–column connection being analyzed. The number of experi-
ments evaluated varies from code to code because specimens that were not tested under
direct shear could not be evaluated by NEN-EN 1992-1-1:2005 [8] and Model Code 2010 [9].
In Table 4, the statistical results are first presented in general and are then presented
separately by the type of slab–column connection being analyzed. The number of experi-
ments evaluated varies from code to code because specimens that were not tested under
direct shear could not be evaluated by NEN-EN 1992-1-1:2005 [8] and Model Code 2010 [9].
As can be observed in Table 2, ACI 318-19 [7] and NEN-EN 1992-1-1:2005 [8] tend to
be on the conservative side in terms of the average tested to predicted shear stresses. The
Model Code 2010 [9] predicted shear capacities that are on average below the tested shear
stress at failure (average = 0.82). Considering the overall results, the Model Code 2010 [9]
also has the smallest coefficient of variation (COV = 29.63%) on the tested to predicted
shear stresses. These results could be considered generally unsatisfactory. Nevertheless,
all three codes performed better when evaluating only internal slab–column connections
without shear reinforcement, with COV values under 20% for all. The tested to predicted
values using NEN-EN 1992-1-1:2005 [8] show the lowest maximum value for one entry
(0.24); however, the tested to predicted value using the Model Code 2010 [9] for this entry is
small as well (0.58). The entry analyzed is named C/C/4, from Ref. [30], and it is a corner
slab–column connection, unreinforced in shear (see Table A8 in Appendix B).
All three models evaluated performed differently when predicting strengths on internal,
edge and corner slab–column connections. ACI 318-19 [7] and NEN-EN 1992-1-1:2005 [8]
tend to render more conservative results for edge and corner slab–column connections, as
shown in Tables 3 and 4. NEN-EN 1992-1-1:2005 [8] showed the largest scatter when evalu-
ating edge slab–column connections (COV = 51.93%), which is considered unacceptable.
As can be observed in Figure 19 and Table 5, the ACI 318-19 [7] assumption of a
linear elastic model distribution leads to an overestimation of the real performance of the
specimens. Replacing the shear stresses calculated assuming a linear stress distribution with
the results presented in Tables A8–A10, which are the maximum shear internal forces of the
slab at failure, measured on the punching perimeter described in the ACI 318-19 code [7]
divided by the effective depth of the slab d in the comparison to the shear capacity from ACI
318-19 (see Table 5), results in larger (i.e., more overly conservative) values for the tested
to predicted shear. At the same time, the COV on the tested to predicted values increases,
indicating that when using ACI 318-19, both the shear stress from the code provisions
and the shear capacity should be used together. Using the results from the FEM models,
increasing the level of precision of the assumption of a linear elastic model distribution,
which was demonstrated to lead to an overestimation of the shear capacity of the slab,
tends to enlarge this overestimation.
As can be observed in Figure 20 and Table 6, the ACI 318-19 [7] assumption of a
linear elastic model distribution does not hold up for finer calculations using the SCIA
Engineer [33] FEM results. Although the average value for the comparison between the
shear capacities for ACI 318-19 [7] and SCIA Engineer [33] FEM results for the acting shear
stress is relatively close to 1.0, the COV on these results is not acceptable (COV = 54.37) (see
Tables A11–A13).
Buildings 2022, 12, 2092 27 of 49
Figure 20. Comparison between the SCIA Engineer [33] FEM results and the shear stress for ACI
318-19 [7], calculated using Equation (1).
Table 6. Statistical results from the average value for the comparison between the acting shear stress
for ACI 318-19 [7] and SCIA Engineer [33] FEM results for the acting shear stress.
Figure 21. Parameters studied based on the comparison between experimental vtest and predicted
shear capacities vpred for three design methods from existing codes: (a) concrete compressive strength
fc ; (b) effective depth d; (c) longitudinal reinforcement ratio ρl ; (d) shear span to depth ratio a/d;
(e) eccentricity.
As can be observed in Figure 21a, ACI 318-19 [7] tends to underestimate the influence
of the compressive strength of the concrete fc . On the other hand, the fib Model Code [9]
and NEN-EN 1992-1-1:2005 [8] correctly take the influence of the concrete compressive
strength into account.
When considering the influence of the effective depth d on the tested to predicted ratios
(Figure 21b), we can see that NEN-EN 1992-1-1:2005 [8] and the fib Model Code [9] take the
size effect correctly into account for the specimens tested. In contrast, for ACI 318-19 [7],
no influence of size is considered. Additionally, we observe an important increase in the
Buildings 2022, 12, 2092 30 of 49
conservativism of NEN-EN 1992-1-1:2005 [8] as the effective depth increases. For the tested
sizes, the value of d only influences the size and location of the punching perimeter. None
of the specimens were of a thickness at which the size effect factor from the code equations
starts to play a role.
As can be seen in Figure 21c, NEN-EN 1992-1-1:2005 [8] and the fib Model Code
2010 [9] consider that the influence of the longitudinal reinforcement is larger than that
observed in the experiments. The influence is taken into account more realistically in
ACI 318-19 [7]. As the specimens in the database contained large reinforcement ratios,
experiments on slabs with lower reinforcement ratios are necessary to further study the
effect of the reinforcement ratio on the tested to predicted ratios determined with the
studied codes.
Figure 21d shows that, for the Model Code 2010 [9] and NEN-EN 1992-1-1:2005 [8], as
the shear span to depth ratio a/d increases, the values of vtest /vpred tend to become more
conservative. On the other hand, for ACI 318-19 [7], the values of vtest /vpred tend to remain
relatively constant as the shear span to depth ratio increases. None of the studied specimens
are in the range of small a/d values for which direct transfer of the load to the support
plays a role.
Figure 21e shows that, for the Model Code 2010 [9] ACI 318-19 [7], the values of
vtest /vpred remain constant despite the increase in the absolute value of eccentricity. Nonethe-
less, for NEN-EN 1992-1-1:2005 [8], these values drastically decrease and become unsafe
when eccentricities exceed 1000 mm. This observation can be explained by Equation (18) for
most edge and corner slab–column connection entries due to the direction of the eccentricity.
This equation does not consider the magnitude of the eccentricity being evaluated, just its
direction; thus, for large eccentricities, this approach becomes unsafe.
4. Discussion
High-strength concrete slabs are not included in the study due to the lack of experi-
ments on high-strength concrete slab–column connections. Future studies may investigate
the behavior of this type of slabs in comparison with the ones considered in this study.
Ngo [36] presented a study on concentric experiments on high-strength concrete slabs
without shear reinforcement and concluded that the use of high-strength concrete improves
the punching shear resistance. However, as the aggregate interlock capacity decreases for
higher strength concrete, further research on this topic is warranted.
Real-scale slab–column experiments on punching shear are not commonly found in
the literature. None of the entries in this database is considered a realistic size slab, as none
of these has an effective depth over 200 mm. Since a size effect occurs in punching shear
for concentric slab–column connections [37], experimental research on larger slabs under
eccentric punching shear is necessary.
In most cases, the results found in the literature indicate that there is an important
reduction in the punching capacity when unbalanced moments occur in the slab–column
connection. Nevertheless, most experiments and research works focus on concentric
punching shear. The code provisions are based on empirical equations, which include the
effect of eccentricities by different methods, such as critical perimeter reduction or increase
in the applied shear stress, but a mechanics-based model that is practical enough to be
implemented in the building codes is lacking. Mechanics-based models, such as the critical
shear crack theory used in the Model Code 2010 [9], are developed for the case of concentric
punching shear and use simplified methods for the extension to eccentric punching shear.
For this database, the empirical methods showed large scatter in the results of the
tested to predicted capacities, represented by the high coefficients of variation. This
observation may be explained by the fact that all methods under consideration were
originally developed for concentric punching shear and validated with concentric punching
shear tests and were then extended to the use for eccentric punching shear. For the future
development of design codes, more attention should be paid to eccentric punching shear
and the mechanical basis of the problem. Lower bound plasticity-based models have been
Buildings 2022, 12, 2092 31 of 49
shown to lead to good results in an exploratory study [38] and are a promising path toward
better code models for all cases of punching shear.
5. Conclusions
The lack of understanding regarding eccentric punching shear presents a practical
problem because local forces typically control the slab’s design. The transfer of unbalanced
moments from the slab to the column causes an increase in the resulting shear stresses.
When this effect is not well understood, it may lead to a punching failure of the slab–column
connection and a possible collapse of the building. This study evaluates the available code
provisions against 128 experimental results reported in the literature.
Analyzing the available experimental results from the database resulted in the follow-
ing conclusions:
• There is a lack of experiments on eccentric punching shear.
• All experiments are carried out on slabs of under 200 mm in depth. As such, the
experiments cannot be used to evaluate the size effect in shear.
• Most specimens have large reinforcement ratios to avoid a flexural failure before
reaching the punching shear capacity of the slab and are not representative of typical
floor slabs.
• All specimens are cast using normal-strength concrete.
From the comparison between the experimental shear capacities and the capacities
predicted by the available codes, the following conclusions result:
• The presence of unbalanced moments has a large impact on the capacity of the codes
to accurately predict the performance of slab–column connections, especially for edge
and corner connections.
• The closest to 1.0 average value for tested to predicted shear capacity is obtained
with the Model Code 2010 [9] provisions, although the results are on the unsafe side
(average of tested/predicted = 0.82).
• Evaluating all experiments, the coefficient of variation of the tested to predicted shear
capacities is lower for the expressions of the Model Code 2010 [9], based on the critical
shear crack theory, than for the empirical expressions of NEN-EN 1992-1-1:2005 [8]
and ACI 318-19 [7].
• In general, the coefficient of variation of the tested to predicted shear capacities is
lower for the experiments without shear reinforcement than for the experiments
with shear reinforcement. This observation indicates that further studies should
address the distribution of shear stresses for shear-reinforced slab–column connections
with unbalanced moments, as well as the contribution of the various shear-carrying
mechanisms for the case with shear reinforcement.
• The NEN-EN 1992-1-1:2005 [8] results for tested to predicted shear capacity showed a
smaller scatter for internal and corner slab–column connections than the results of the
Model Code 2010 [9] and ACI 318-19 [7].
• The Model Code 2010 [9] has the lowest coefficient of variation for the tested to
predicted shear stress for edge slab–column connections, whereas NEN-EN 1992-1-
1:2005 [8] has the smallest COV for internal and corner slab–column connections.
• The fib Model Code generally results in the best performance in the parameter studies
of the tested to predicted ratios, indicating that this mechanics-based model has a
better representation of the parameters than the empirical models in NEN-EN 1992-1-
1:2005 [8] and ACI 318-19 [7].
• For eccentricities over 1000 mm, the use of NEN-EN 1992-1-1:2005 [8] becomes unsafe
for inward eccentricities at edge slab–column connections.
From the comparison between the SCIA Engineer [33] FEM results and the predicted
shear capacities for ACI 318-19 [7], the following conclusions result:
• Using a more refined calculation method to determine the acting shear stresses does
not improve the results in terms of tested to predicted shear stresses when using
Buildings 2022, 12, 2092 32 of 49
the ACI 318-19 [7] shear capacities. This observation indicates that the shear stress
distribution that follows from a linear finite element analysis does not correspond to
the slab at failure when cracking and redistribution of stresses occur.
• The ACI 318-19 [7] shear capacity should be combined with the method for determin-
ing the shear stress on the punching perimeter described in ACI 318-19 [7].
• Further research is necessary on the capacity methods that can be combined with a
linear finite element analysis.
A better understanding of eccentric punching shear and further experiments on
deeper slabs, slabs with high-strength concrete and carefully instrumented slabs with shear
reinforcement are necessary to obtain safe designs, optimize the design of building floors
and develop better tools for the assessment of existing building slabs.
Author Contributions: Conceptualization, E.O.L.L.; methodology, D.V. and E.O.L.L.; software, D.V.;
validation, A.S.G. and E.O.L.L.; formal analysis, D.V.; investigation, D.V. and E.O.L.L.; resources,
E.O.L.L.; data curation, D.V. and E.O.L.L.; writing—original draft preparation, D.V. and E.O.L.L.;
writing—review and editing, A.S.G.; visualization, D.V.; supervision, E.O.L.L.; project administration,
E.O.L.L.; funding acquisition, E.O.L.L. All authors have read and agreed to the published version of
the manuscript.
Funding: This research is part of the program of Collaboration Grants 2019 of Universidad San Fran-
cisco de Quito. The APC was funded by the open access initiative of Delft University of Technology.
Data Availability Statement: The data and calculations are available in the public domain through
https://doi.org/10.5281/zenodo.7317297 (accessed on 3 April 2022).
Acknowledgments: The authors would like to thank the program of Collaboration Grants 2019 of
Universidad San Francisco de Quito for the financial support.
Conflicts of Interest: The authors declare no conflict of interest.
Nomenclature
Appendix A
Reference Name Lx (mm) Lx (mm) c1 (mm) c2 (mm) h (mm) d (mm) a (mm) av (mm)
L1 2280 2280 305 305 178 151 1026 874
L3 2280 2280 305 305 178 151 1026 874
L4 2280 2280 305 305 178 151 1026 874
Narayani [19]
L5 2280 2280 305 305 178 151 1026 874
L6 2280 2280 305 305 178 151 1026 874
L10 2280 2280 305 305 178 151 1026 874
P16A 3000 3000 300 300 150 121 1375 1225
Krüger [2] P30A 3000 3000 300 300 150 121 1375 1225
PP16B 3000 3000 300 300 150 121 1375 1225
M2A 1829 1829 305 305 152 114 787 635
M4A 1829 1829 305 305 152 114 787 635
M2 1829 1829 305 305 152 114 787 635
M3 1829 1829 305 305 152 114 787 635
Moe [11] M6 1829 1829 254 254 152 114 787 660
M7 1829 1829 254 254 152 114 787 660
M8 1829 1829 254 254 152 114 787 660
M9 1829 1829 254 254 152 114 787 660
M10 1829 1829 254 254 152 114 787 660
B.3 1524 1524 203 203 102 76 737 635
B.4 1524 1524 203 203 102 76 737 635
Anis [21]
B.5 1524 1524 203 203 102 76 737 635
B.6 1524 1524 203 203 102 76 737 635
C/I/1 914 914 127 127 76 56 419 355
C/I/2 914 914 127 127 76 56 419 355
C/I/3 914 914 127 127 76 56 419 355
C/I/4 914 914 127 127 76 56 419 355
Stamenkovic [23]
C/Ir/1 914 914 152 76 76 56 419 343
C/Ir/2 914 914 152 76 76 56 419 343
C/Ir/3 914 914 152 76 76 56 419 343
C/Ir/4 914 914 152 76 76 56 419 343
A12 2134 1219 152 152 76 62 1067 991
Hanson [22] B16 2134 1219 152 305 76 62 1067 991
C17 2134 1219 305 152 76 62 1067 914
S3 2500 2500 300 300 180 145 1250 1100
S4 2500 2500 300 300 180 143 1250 1100
Pina Ferreira [24]
S6 2500 2500 300 300 180 144 1250 1100
S8 2500 2500 300 300 180 144 1250 1100
Buildings 2022, 12, 2092 35 of 49
Reference Name ρ (%) fy (MPa) Es (GPa) fc (MPa) Age (Days) fct (MPa) dg (mm)
L1 1.78% 398 188 32.80 28 2.70 19.00
L3 1.78% 398 188 33.10 28 2.66 19.00
L4 1.78% 398 188 45.80 28 3.45 19.00
Narayani [19]
L5 1.78% 398 188 35.00 28 2.98 19.00
L6 1.78% 398 188 42.10 28 3.07 19.00
L10 1.78% 398 188 41.80 28 2.47 19.00
P16A 1.30% 480 200 35.00 28 4.50 16.00
Krüger [2] P30A 1.30% 480 200 35.00 28 4.50 16.00
PP16B 1.30% 480 200 35.00 28 4.50 16.00
M2A 1.50% 481 196 15.51 25 2.99 38.10
M4A 1.50% 481 196 17.65 23 3.19 38.10
M2 1.50% 481 196 25.72 22 3.85 38.10
M3 1.50% 481 196 22.72 20 3.62 38.10
Moe [11] M6 1.34% 328 196 26.48 26 3.91 38.10
M7 1.34% 328 196 24.96 24 3.80 38.10
M8 1.34% 328 196 24.61 24 3.77 38.10
M9 1.34% 328 196 23.24 22 3.66 38.10
M10 1.34% 328 196 21.10 25 3.49 38.10
B.3 2.19% 331 205 38.06 28 4.69 9.53
B.4 2.19% 331 205 37.23 28 4.64 9.53
Anis [21]
B.5 2.19% 331 205 36.20 28 4.57 9.53
B.6 2.19% 331 205 39.16 28 4.76 9.53
C/I/1 1.17% 413 192 45.02 7 5.10 9.53
C/I/2 1.17% 413 192 37.09 7 4.63 9.53
C/I/3 1.17% 413 192 31.92 7 4.29 9.53
C/I/4 1.17% 413 192 31.37 7 4.26 9.53
Stamenkovic [23]
C/Ir/1 1.17% 413 192 28.27 7 4.04 9.53
C/Ir/2 1.17% 413 192 36.54 7 4.59 9.53
C/Ir/3 1.17% 413 192 35.71 7 4.54 9.53
C/Ir/4 1.17% 413 192 33.23 7 4.38 9.53
A12 1.36% 372 192 33.23 28 4.38 9.53
Hanson [22] B16 1.36% 341 192 30.41 28 4.19 9.53
C17 1.36% 341 192 35.99 28 4.56 9.53
S3 1.46% 540 213 50.30 28 4.30 9.50
S4 1.48% 540 213 49.20 28 4.40 9.50
Pina Ferreira [24]
S6 1.47% 540 213 50.10 28 4.90 9.50
S8 1.47% 540 213 48.40 28 4.00 9.50
Buildings 2022, 12, 2092 36 of 49
Reference Name Lx (mm) Lx (mm) c1 (mm) c2 (mm) h (mm) d (mm) a (mm) av (mm)
L1 2350 1700 300 300 180 147 2000 1850
L2 2350 1700 300 300 180 146 2000 1850
L3 2350 1700 300 300 180 146 2000 1850
L4 2350 1700 300 300 180 146 2000 1850
L5 2350 1700 300 300 180 146 2000 1850
L6 2350 1700 300 300 180 146 2000 1850
Albuquerque [17] L7 2350 1700 300 300 180 146 2000 1850
L8 2350 1700 300 300 180 146 2000 1850
L9 2350 1700 300 300 180 146 2000 1850
L10 2350 1700 300 300 180 146 2000 1850
L11 2350 1700 300 300 180 146 2000 1850
L12 2350 1700 300 300 180 146 2000 1850
L13 2350 1700 300 300 180 146 2000 1850
ES1 1295 2280 305 305 178 151 1029 876
ES2 1295 2280 305 305 178 151 1029 876
ES3 1295 2280 305 305 178 151 1029 876
Narayani [19] ES4 1295 2280 305 305 178 151 1029 876
ES5 1295 2280 305 305 178 151 1029 876
ES6 1295 2280 305 305 178 151 1029 876
ES7 1295 2280 305 305 178 151 1029 876
Z—IV (1) 965 1829 178 178 152 121 870 781
Z—V (1) 965 1829 267 267 152 121 826 692
Z—V (2) 965 1829 267 267 152 121 826 692
Z—V (3) 965 1829 267 267 152 118 826 692
Zaghlool [20]
Z—V (4) 965 1829 267 267 152 121 826 692
Z—V (5) 965 1829 267 267 152 121 826 692
Z—V (6) 965 1829 267 267 152 121 826 692
Z—VI (1) 965 1829 356 356 152 121 781 603
M(T)/E/1 914 914 127 127 76 56 812 749
M(T)/E/2 914 914 127 127 76 56 812 749
C(T)/E/1 914 914 127 127 76 56 812 749
C(T)/E/2 914 914 127 127 76 56 812 749
C(T)/E/3 914 914 127 127 76 56 812 749
C(T)/E/4 914 914 127 127 76 56 812 749
Stamenkovic [23]
V/E/1 914 914 127 127 76 56 812 749
M(II)/E/1 914 914 127 127 76 56 812 749
C(II)/E/1 914 914 127 127 76 56 812 749
C(II)/E/2 914 914 127 127 76 56 812 749
C(II)/E/3 914 914 127 127 76 56 812 749
C(II)/E/4 914 914 127 127 76 56 812 749
Buildings 2022, 12, 2092 37 of 49
Reference Name Lx (mm) Lx (mm) c1 (mm) c2 (mm) h (mm) d (mm) a (mm) av (mm)
Hanson [22] D15 1143 1219 152 152 76 62 1067 991
Ritchie [25] 1 1350 1900 250 250 150 122 1175 1050
E1 711 1219 203 203 140 105 610 508
E2 711 1219 203 203 140 105 610 508
E4 711 1219 203 203 140 105 610 508
E1-1 1571 1719 203 203 140 105 1220 1118
E1-2 1571 1719 203 203 140 105 1220 1118
Sudarsana [26]
E1-3 1571 1719 203 203 140 105 1220 1118
E1-4 1571 1719 203 203 140 105 1220 1118
E2-1 961 2940 203 203 140 105 610 508
E2-2 961 2940 203 203 140 105 610 508
E2-3 961 2940 203 203 140 105 610 508
E2-4 961 2940 203 203 140 105 610 508
ZJESSS 1060 1770 250 250 150 119 850 725
Zaghloul [27]
ZJES 1060 1770 250 250 150 119 850 725
Reference Name ρ (%) fy (MPa) Es (GPa) fc (MPa) Age (Days) fct (MPa) dg (mm)
L1 1.00% 558 192 46.80 28 3.40 9.50
L2 1.30% 558 192 44.70 28 3.00 9.50
L3 1.30% 558 192 45.10 28 3.10 9.50
L4 1.30% 558 192 46.00 28 3.30 9.50
L5 1.30% 558 192 51.40 28 4.10 9.50
L6 1.30% 558 192 52.10 28 4.30 9.50
Albuquerque [17] L7 1.50% 558 192 50.00 28 3.70 9.50
L8 1.40% 558 192 50.50 28 3.90 9.50
L9 1.50% 558 192 57.60 28 3.20 9.50
L10 1.50% 558 192 59.30 28 3.60 9.50
L11 1.50% 558 192 43.10 28 3.10 9.50
L12 1.50% 558 192 43.60 28 3.30 9.50
L13 1.50% 558 192 44.10 28 3.40 9.50
ES1 0.88% 398 188 33.80 28 2.70 19.00
ES2 0.88% 398 188 32.80 28 2.70 19.00
ES3 0.88% 398 188 51.30 28 3.10 19.00
Narayani [19] ES4 0.88% 398 188 50.00 28 3.31 19.00
ES5 0.88% 398 188 37.60 28 2.63 19.00
ES6 0.88% 398 188 40.40 28 2.23 19.00
ES7 0.88% 398 188 45.80 28 3.38 19.00
Buildings 2022, 12, 2092 38 of 49
Reference Name ρ (%) fy (MPa) Es (GPa) fc (MPa) Age (Days) fct (MPa) dg (mm)
Z—IV (1) 1.23% 476 207 27.34 28 2.99 19.05
Z—V (1) 1.23% 474 207 34.34 28 3.52 19.05
Z—V (2) 1.65% 474 207 40.47 28 3.61 19.05
Z—V (3) 2.23% 475 207 38.75 28 3.79 19.05
Zaghlool [20]
Z—V (4) 1.23% 475 207 35.03 28 4.10 19.05
Z—V (5) 1.23% 476 207 35.16 28 3.58 19.05
Z—V (6) 1.23% 476 207 31.30 28 3.63 19.05
Z—VI (1) 1.23% 476 207 25.99 28 2.83 19.05
M(T)/E/1 1.17% 413 192 30.34 7 4.19 9.53
M(T)/E/2 1.17% 413 192 33.09 7 4.37 9.53
C(T)/E/1 1.17% 413 192 38.47 7 4.71 9.53
C(T)/E/2 1.17% 413 192 32.41 7 4.33 9.53
C(T)/E/3 1.17% 413 192 33.99 7 4.43 9.53
C(T)/E/4 1.17% 413 192 34.34 7 4.45 9.53
Stamenkovic [23]
V/E/1 1.17% 413 192 35.85 7 4.55 9.53
M(II)/E/1 1.17% 413 192 36.20 7 4.57 9.53
C(II)/E/1 1.17% 413 192 34.82 7 4.48 9.53
C(II)/E/2 1.17% 413 192 35.51 7 4.53 9.53
C(II)/E/3 1.17% 413 192 34.89 7 4.49 9.53
C(II)/E/4 1.17% 413 192 36.54 7 4.59 9.53
Hanson [22] D15 1.15% 365 192 31.10 28 4.24 9.53
Ritchie [25] 1 1.12% 432 192 26.20 28 3.89 9.50
E1 0.90% 420 183 43.62 28 5.02 10.00
E2 0.90% 420 183 42.41 28 4.95 10.00
E4 0.90% 420 183 43.62 28 5.02 10.00
E1-1 1.60% 420 183 52.80 28 5.52 10.00
E1-2 1.60% 420 183 52.80 28 5.52 10.00
Sudarsana [26] E1-3 1.60% 420 183 55.00 28 5.64 10.00
E1-4 1.60% 420 183 52.80 28 5.52 10.00
E2-1 0.80% 420 183 52.80 28 5.52 10.00
E2-2 0.80% 420 183 52.80 28 5.52 10.00
E2-3 0.80% 420 183 55.00 28 5.64 10.00
E2-4 0.80% 420 183 55.00 28 5.64 10.00
ZJESSS 1.40% 400 192 42.00 28 4.93 10.00
Zaghloul [27]
ZJES 1.40% 400 192 42.00 28 4.93 10.00
Buildings 2022, 12, 2092 39 of 49
Reference Name Lx (mm) Lx (mm) c1 (mm) c2 (mm) h (mm) d (mm) a (mm) av (mm)
Z—I (1) 1067 1067 178 178 152 121 965 1067
Z—II (1) 1067 1067 267 267 152 121 921 1067
Z—II (2) 1067 1067 267 267 152 121 921 1067
Zaghlool [20] Z—II (3) 1067 1067 267 267 152 118 921 1067
Z—II (4) 1067 1067 267 267 152 121 921 1067
Z—II (6) 1067 1067 267 267 152 121 921 1067
Z—III (1) 1067 1067 356 356 152 121 876 1067
C5 711 711 305 305 140 105 559 711
C6 711 711 305 305 140 105 559 711
Sudarsana [26]
C7 711 711 305 305 140 105 559 711
C8 711 711 305 305 140 105 559 711
S101 530 530 100 100 100 80 480 530
S201 530 530 100 100 100 80 480 530
S301 530 530 100 100 100 80 480 530
Desayi [28]
S102 530 530 100 100 100 80 480 530
S202 530 530 100 100 100 80 480 530
S302 530 530 100 100 100 80 480 530
SC1 1525 1525 300 300 125 100 1375 1525
SC2 1525 1525 300 300 125 100 1375 1525
SC3 1525 1525 300 300 125 100 1375 1525
SC4 1525 1525 220 220 125 100 1415 1525
SC5 1525 1525 220 220 125 100 1415 1525
Walker [29]
SC7 1525 1525 220 220 125 100 1415 1525
SC8 1000 1000 160 160 80 64 920 1000
SC9 1000 1000 160 160 80 64 920 1000
SC11 1000 700 160 160 80 60 920 1000
SC12 1000 700 300 300 80 60 850 1000
C/C/1 914 914 127 127 76 56 813 914
C/C/2 914 914 127 127 76 56 813 914
Stamenkovic [30]
C/C/3 914 914 127 127 76 56 813 914
C/C/4 914 914 127 127 76 56 813 914
NH1 1075 1075 250 250 150 114 910 1075
NH2 1075 1075 250 250 150 114 910 1075
Ghali [18] NH3 1075 1075 250 250 150 114 910 1075
NH4 1075 1075 250 250 150 114 910 1075
NH5 1075 1075 250 250 150 114 910 1075
Buildings 2022, 12, 2092 40 of 49
Reference Name ρ (%) fy (MPa) Es (GPa) fc (MPa) Age (Days) fct (MPa) dg (mm)
Z—I (1) 1.23% 379 207 32.68 28 4.34 19.05
Z—II (1) 1.23% 389 207 33.03 28 4.37 19.05
Z—II (2) 1.65% 405 207 33.44 28 4.39 19.05
Zaghlool [20] Z—II (3) 2.23% 451 207 27.72 28 4.00 19.05
Z—II (4) 1.23% 389 207 30.75 28 4.21 19.05
Z—II (6) 1.23% 381 207 33.58 28 4.40 19.05
Z—III (1) 1.23% 379 207 33.65 28 4.41 19.05
C5 1.11% 420 183 44.40 28 5.06 10.00
C6 1.11% 420 183 44.40 28 5.06 10.00
Sudarsana [26]
C7 1.11% 420 183 44.40 28 5.06 10.00
C8 1.11% 420 183 44.40 28 5.06 10.00
S101 0.53% 720 192 45.00 28 5.10 9.50
S201 0.80% 720 192 45.00 28 5.10 9.50
S301 1.07% 720 192 25.00 28 3.80 9.50
Desayi [28]
S102 0.53% 720 192 31.00 28 4.23 9.50
S202 0.80% 720 192 34.00 28 4.43 9.50
S302 1.07% 720 192 28.00 28 4.02 9.50
SC1 1.14% 595 192 43.30 28 5.00 20.00
SC2 1.11% 595 192 47.90 28 5.26 20.00
SC3 1.13% 595 192 37.40 28 4.65 20.00
SC4 1.14% 595 192 40.80 28 4.85 20.00
SC5 1.71% 595 192 46.50 28 5.18 20.00
Walker [29]
SC7 1.71% 595 192 43.80 28 5.03 20.00
SC8 1.37% 595 192 37.40 28 4.65 20.00
SC9 1.24% 595 192 34.30 28 4.45 20.00
SC11 1.27% 595 192 27.20 28 3.96 20.00
SC12 1.18% 595 192 40.70 28 4.85 20.00
C/C/1 1.17% 413 192 38.06 28 4.69 9.50
C/C/2 1.17% 413 192 35.37 28 4.52 9.50
Stamenkovic [30]
C/C/3 1.17% 413 192 32.27 28 4.32 9.50
C/C/4 1.17% 413 192 38.27 28 4.70 9.50
NH1 1.45% 440 200 41.50 28 4.90 9.50
NH2 1.45% 440 200 42.20 28 4.94 9.50
Ghali [18] NH3 1.45% 440 200 36.40 28 4.59 9.50
NH4 1.45% 440 200 36.90 28 4.62 9.50
NH5 1.45% 440 200 33.20 28 4.38 9.50
Buildings 2022, 12, 2092 41 of 49
Appendix B
Hanson [22] B16 2.21 1.82 1.22 1.12 1.24 0.90 1.98 3.31 0.60
C17 1.94 1.98 0.98 1.21 1.32 0.92 1.81 3.60 0.50
S3 4.63 3.16 1.47 2.45 1.72 1.42 1.87 1.74 1.08
S4 4.83 2.70 1.79 2.57 2.34 1.10 1.33 1.18 1.13
Pina Ferreira [24]
S6 3.52 2.34 1.51 1.86 1.51 1.23 3.50 3.95 0.89
S8 5.62 3.05 1.84 3.03 2.32 1.31 1.71 1.42 1.21
Table A11. Internal slab–column connections—SCIA (FEM) [33] results compared with shear stress
according to ACI 318-19, Equation (1), and capacity results for ACI 319-19 [7].
Table A12. Edge slab–column connections—SCIA (FEM) [33] results compared with shear stress
according to ACI 318-19, Equation (1), and capacity results for ACI 319-19 [7].
Table A13. Corner slab–column connections—SCIA (FEM) [33] results compared with shear stress
according to ACI 318-19, Equation (1), and capacity results for ACI 319-19 [7].
References
1. Wight, J.; MacGregor, J. Analysis, and Design. In Reinforced Concrete Mechanics & Design, 6th ed.; Pearson: London, UK, 2010;
Volume 1.
2. Krüger, G. Résistance au Poinçonnement Excentré des Planchers-Dalles; EPFL: Écublens, Switzerland, 1999.
3. Muttoni, A.; Ruiz, M.F.; Simões, J.T. The theoretical principles of the critical shear crack theory for punching shear failures and
derivation of consistent closed-form design expressions. Struct. Concr. 2018, 19, 174–190. [CrossRef]
4. Mirzaei, Y. Post punching behaviour of reinforced concrete slab-column connections. In Proceedings of the 7th International FIB
PhD Symposium, Stuttgart, Germany, 11–13 September 2008.
5. Kinnunen, S.; Nylander, H. Punching of concrete slabs without shear reinforcement. Elander 1960, 158, 112.
6. Muttoni, A.; Ruiz, M.F. The levels-of-approximation approach in MC 2010: Application to punching shear provisions. Struct.
Concr. 2012, 13, 32–41. [CrossRef]
7. ACI Committee 318. ACI CODE-318-19: Building Code Requirements for Structural Concrete and Commentary; American Concrete
Institute: Farmington Hills, MI, USA, 2019.
8. European Committee for Standardization and British Standards Institution. Eurocode 2: Design of Concrete Structures; British
Standards Institution: London, UK, 2005.
9. The International Federation for Structural Concrete. fib Bulletin 55. Model Code 2010 First Complete Draft; fib: Lausanne,
Switzerland, 2010; Volume 1. [CrossRef]
10. Park, H.; Choi, K. Improved Strength Model for Interior Flat Plate-Column Connections Subject to Unbalanced Moment. ASCE J.
Struct. Eng. 2006, 694, 694–704. [CrossRef]
11. Moe, J. Shearing Strength of Reinforced Concrete Slabs and Footings under Concentrated Loads; Development Department Bulletin No.
D47; Portland Cement Association: Skokie, IL, USA, 1961.
12. ACI-ASCE Committee 426. The Shear Strength of Reinforced Concrete Members. Chapter 5, Shear Strength of Slabs. Proc. ASCE
J. Struct. Div. 1974, 100, 1543–1591. [CrossRef]
13. Mast, P.E. Stresses in Flat Plates Near Columns. J. Proc. 1970, 67, 761–768.
14. The International Federation for Structural Concrete. fib Bulleting 2. Textbook on Behavior, Design and Performance Updated Knowledge
of the CEB/FIP Model Code 1990; fib: Lausanne, Switzerland, 1999; Volume 2. [CrossRef]
Buildings 2022, 12, 2092 49 of 49
15. Muttoni, A. Punching Shear Strength of Reinforced Concrete Slabs without Transverse Reinforcement. ACI Struct. J. 2008,
105, 440–450.
16. Fernández Ruiz, M.; Muttoni, A. Application of Critical Shear Crack Theory to Punching of Reinforced Concrete Slabs with
Transverse Reinforcement. ACI Struct. J. 2009, 106, 485–494.
17. Albuquerque, N.G.B.; Melo, G.S.; Vollum, R.L. Punching Shear Strength of Flat Slab-Edge Column Connections with Outward
Eccentricity of Loading. ACI Struct. J. 2016, 113, 1117–1129. [CrossRef]
18. Hammil, N.; Ghali, A. Punching Shear Resistance of Corner Slab-Column Connections. ACI Struct. J. 1994, 91, 697–707. [CrossRef]
19. Narayani, N. Shear Reinforcement in Reinforced Concrete Column Heads. Ph.D. Thesis, Faculty of Engineering, Imperial College
of Science and Technology, London, UK, 1971.
20. Zaghlool, E.E.-D.R.F. Strength and Behaviour of Corner and Edge Column-Slab Connections in Reinforced Concrete Flat Plates.
Ph.D. Thesis, University of Calgary, Calgary, AB, Canada, 1971.
21. Anis, N.N. Shear Strength of Reinforced Concrete Flat Slabs without Shear Reinforcement. Ph.D. Thesis, Imperial College London,
London, UK, 1970.
22. Hanson, N.W.; Hanson, J.M. Shear and Moment Transfer between Concrete Slabs and Columns. J. PCA Res. Dev. Lab. 1968,
11, 2–16. Available online: https://www.concrete.org/publications/internationalconcreteabstractsportal/m/details/id/19463
(accessed on 3 April 2022).
23. Stamenkovic, A. Local Strength of Flat Slabs At Column Heads. Ph.D. Thesis, Imperial College London, London, UK, 1970.
24. De Pina Ferreira, M.; Oliveira, M.H.; Melo, G.S.S.A. Tests on the punching resistance of flat slabs with unbalanced moments. Eng.
Struct. 2019, 196, 109311. [CrossRef]
25. Ritchie, M.; Ghali, A.; Dilger, W.; Gayed, R.B. Unbalanced Moment Resistance by shear in Slab-Column Connections: Experimental
Assessment. ACI Struct. J. 2006, 103, 74–82.
26. Sudarsana, I.K. Punching Shear in Edge and Corner Column Slab Connections of Flat Plate Structures. Ph.D. Thesis, Department
of Civil Engineering, University of Ottawa, Ottawa, ON, Cananda, December 2001.
27. Zaghloul, A. Punching Shear Strength of Interior and Edge Column-Slab Connections in CFRP Reinforced Flat Plate Structures
Transferring Shear and Moment. Ph.D. Thesis, Department of Civil and Environmental Engineering, University of Ottawa,
Ottawa, ON, Canada, February 2007.
28. Desayi, P.; Seshadri, H.K. Punching shear strength of flat slab corner column connections. In Part 1. Reinforced Concrete Connections;
Department of Civil Engineering, Indian Institute of Science: Bangalore, India, 1997.
29. Walker, P.R.; Regan, P.E. Corner column-slab connections in concrete flat plates. ASCE J. Struct. Eng. 1987, 113, 704–720. [CrossRef]
30. Stamenkovic, A.; Chapman, J.C. Local strength at column heads in flat slabs subjected to a combines vertical and horizontal
loading. Proc. ICE 1974, 57, 205–232.
31. Vargas, D.; Lantsoght, E.; Genikomsou, K. Spreadsheet for Flat Slabs in Eccentric Punching Shear: Experimental Database and Analysis;
Zenodo: Quito, Ecuador, 2022. [CrossRef]
32. Sarveghadi, M.; Gandomi, A.H.; Bolandi, H.; Alavi, A.H. Development of prediction models for shear strength of SFRCB using a
machine learning approach. Neural Comput. Appl. 2019, 31, 2085–2094. [CrossRef]
33. Nemetschek Group. SCIA Downloads. 17 May 2022. Available online: https://www.scia.net/en/scia-engineer/downloads
(accessed on 3 April 2022).
34. Lantsoght, E. Database of Shear Experiments on Steel Fiber Reinforced Concrete Beams without Stirrups. Materials 2019, 12, 917.
[CrossRef] [PubMed]
35. Vargas, D.; Lantsoght, E.; Genikomsou, K. Spreadsheet Validation for Flat slabs in Eccentric Punching Shear Experimental Database and
Analysis; Zenodo: Quito, Ecuador, 2022. [CrossRef]
36. Ngo, T. Punching shear resistance of high-strength concrete slabs. Electron. J. Struct. Eng. 2001, 1, 52–59. Available online:
https://ejsei.com/EJSE/article/view/14 (accessed on 3 April 2022). [CrossRef]
37. Guandallini, S.; Burdet, O.L.; Muttoni, A. Punching Tests of Slabs with Low Reinforcement Ratios. ACI Struct. J. 2009, 106, 87–95.
38. Carrera, B.; Lantsoght, E.O.L.; Alexander, S.D.B. Application of Strip Model to Edge Column-Lab Connections Loaded with Outward
Eccentricity; ACI: Farmington Hills, MI, USA, 2022; Volume 353, pp. 124–141. Available online: https://www.concrete.org/
publications/internationalconcreteabstractsportal.aspx?m=details&id=51737115 (accessed on 3 April 2022).