1. Introduction
The various advantages of the flat slabs, such as user-friendliness, ideal surface finish, and absence of pendant joists, along with reducing the construction time, made this system an appropriate structural system regarding the roof applications [
1]. Despite these benefits, the greatest challenge to this design is how to improve the resistance against punching shear. Brittleness and being followed by an unexpected reduction in load-bearing without any previous warning are the main concerns for this type of failure. Therefore, such an unacceptable failure should be prevented [
2,
3]. Various methods have been suggested for rehabilitating reinforced concrete (RC) flat slabs to enhance their punching shear capacity. These approaches include stirrup shear reinforcement, involving the post-installation of steel bolts in the slab [
4,
5,
6,
7]; shear studs, employing vertical bars welded at their tops to square anchor heads and at the bottom to a steel strip [
8,
9,
10]; shear heads, welded to the column and inserted between the layers of the reinforcement [
11,
12,
13]; shear bolts, similar to shear stirrups [
14,
15,
16]; fiber reinforced polymer (FRP) fan, involving FRP rods drilled into the slab’s thickness [
17,
18,
19,
20,
21]; and FRP sheets, attached to the bottom of the slab [
22,
23,
24,
25].
Failure owing to the punching shear in slab–column connections has a complex nature and becomes more complicated when the unbalanced moment, due to span discontinuity or lateral load, needs to be transferred from the slab to the column, which cannot be circumvented at edge slab–column connections. An imbalanced moment results in an asymmetrical distribution of shear forces, affecting the punching shear strength of connections between edge slabs and columns [
26]. Most building codes incorporate a reduction in punching shear capacity, which is typically achieved by decreasing the basic control perimeters, using constant load factors, or assuming a linear or plastic shear distribution along a control perimeter [
27]. Although the punching behavior of strengthened reinforced concrete (RC) slabs is studied widely in other research, only limited information is accessible for strengthened exterior or edge slab–column connections. Furthermore, the code provisions are mainly based on experimental results of moment transfer on interior columns. Exterior and edge slab–column connections are assessed in the ACI code using the eccentric shear method. However, the load-bearing capacity of the exterior slab–column connections is notably undervalued using this method [
28,
29,
30,
31,
32]. Hence, the ACI code adopts an improved strength model, which is derived from experiments conducted on exterior connections. In this model, the unbalanced moment capacity of the exterior connection is determined solely based on the flexural moment capacity of the designated slab width, without taking into account the impact of eccentric shear [
33]. Furthermore, the effect of strengthening these types of connections has not been considered.
This paper aims to investigate the behavior of RC flat slab strengthening against punching shear using FRP strips. Furthermore, the current lack of accessible information on the strengthened exterior or edge slab–column connection is addressed. In addition, the FRP configuration is studied experimentally in order to gain optimum performance. The results of the experimental test are compared with the ACI code to determine the accuracy of the proposed equations for the estimation of the RC slab punching shear capacity. In addition, the FE analysis is performed for a detailed investigation of the response of the RC slab under various loading locations. The nonlinear finite element model of RC slabs is generated using ABAQUS software (
https://en.wikipedia.org/wiki/Abaqus). The numerical analysis is compared with the results obtained by experimental outcomes. Results show a satisfactory correlation with experiments, and finally, a parametric study is performed to determine the effect of loading location on the response of the RC slab.
3. Finite Element Model and Validation
To ensure the accuracy of the numerical results obtained from the finite element software, an experimental slab specimen subjected to static loading is replicated and assessed using the commercial finite element analysis (FEA) software ABAQUS Explicit. The geometry of the RC slab is identical to the experimental slab and is meshed in 3D using
solid elements. This element is a 3D eight-node linear brick element with reduced integration. The “R” in C3D8R stands for “Reduced Integration”, which means that this element uses reduced integration to improve its performance in certain situations [
42,
43]. The steel rebar is modeled using wire elements, and
elements are used for meshing them. This element is a three-node linear triangular element that is used for modeling two-dimensional (2D) structural and plane stress problems. It is often employed in plane stress or plane strain analyses where the out-of-plane deformation is negligible [
44]. In the finite element analysis, the steel material exhibits a constitutive behavior characterized by a bi-linear relationship, adhering to an elastic, perfectly plastic model in both tension and compression. The interaction between the steel bar and the surrounding concrete is treated as a constraint within the “embedded region” [
45].
A concrete damaged plasticity (CDP) model is adopted to represent the uniaxial compressive and tensile response of the concrete. This model accounts for stiffness degradation due to tensile and compressive failure, as well as concrete damping. One of the distinctive features of the CDP model is its ability to model damage evolution in concrete. The model incorporates a damage variable that evolves with loading, representing the initiation and progression of cracks in the material [
46]. The stress–strain relationships are governed by a scalar damaged elasticity equation, as described in Equation (
1), where
represents the modulus of elasticity of concrete, and
,
,
, and
d represent (compressive or tensile) stress, strain, plastic strain, and the damage parameter, respectively [
47].
Certain parameters, including the compressive strength and tensile strength of concrete, as well as the tensile strength of steel rebar, were determined through experimental tests. Concrete exhibited a compressive strength of 45 MPa, a tensile strength of 2.9 MPa, and a modulus of elasticity of 24,400 MPa. The steel rebar had a yielding stress of 400 MPa and an ultimate strength of 600 MPa. However, other parameters essential for modeling the RC slab in ABAQUS were derived from theoretical values based on previous research. The shape factor, denoted as
, represents the ratio of the second stress invariant for tension and compression at the same hydrostatic stress and governs the shape of the yield surface. For normal concrete, the default value of
is 0.667. Another crucial parameter is the ratio of biaxial to uniaxial compression stress (
), which describes the material state during biaxial compression. In the CDP model, the dilation angle (
) and the flow potential eccentricity (
) play a significant role in defining the potential plastic flow of concrete under a three-dimensional stress state. The dilation angle represents the change in volumetric strain during plastic deformation. The flow potential eccentricity signifies the rate at which the flow function approaches the asymptote. The viscosity parameter is an additional factor employed to define concrete in ABAQUS, aimed at enhancing the convergence rate during simulations [
26,
48]. For this parameter, it is recommended to use a small value to improve the convergence rate [
49] or use 0 [
50].
As reported by Daneshvar et al. [
51], they specified the dilation angle (
) as 56, plastic potential eccentricity (
) as 0.1, the ratio of biaxial stress to uniaxial stress (
) as 1.16, the ratio of the second stress invariant on the tensile meridian, which is known as shape parameter (
), as 0.66, and viscosity parameter (
) as 0. The essential input parameters for modeling concrete in ABAQUS are summarized in
Table 3. Given that the specimen supports were constructed using high-strength steel and exhibited no noticeable deformation during testing, they were modeled as rigid body elements in the numerical models [
52]. The boundary conditions for the model were set to simulate simple support at all four sides. Furthermore, the loading plate was represented using rigid elements and was not included in the slab model itself. The simulation of the loading plate and concrete contact problem involves employing surface-to-surface interaction with the hard contact property, utilizing a penalty contact formulation. Using the hard contact, the stress and the resulting displacement are fully transmitted to the slab [
53].
Given that the predominant failure modes observed in the experimental specimens were concrete crushing and punching shear, the CFRP sheets were represented in the model using shell elements, specifically the
element. This element is a 4-node element with six degrees of freedom at each node—three displacements and three rotations—utilizing a transverse shear strain field in its formulation [
54]. To connect the CFRP sheets to the RC slabs in ABAQUS, a “Tie” constraint was employed. In ABAQUS, the “Tie” constraint is a type of constraint that is used to simulate the bonding or connection between two surfaces or sets of nodes. It is often employed in contact analysis to model interactions between different parts of a structure.
Validation
To ensure the accuracy of the numerical simulation, a comparison between numerical and experimental results was conducted.
Figure 4 illustrates the load–deflection curves, revealing a satisfactory agreement between the experimental specimens and numerical models. Discrepancies between the numerical and experimental findings can be attributed to factors such as the finite element model’s characteristics, the degree of concrete heterogeneity, assumptions made in the analysis, and the chosen boundary conditions. The reasonable agreement enhances confidence in the reliability of the numerical models, potentially allowing for an expanded scope of models and parameters in a parametric study.
Table 4 provides a summary of errors in predicting the ultimate load-bearing capacity of RC slabs using numerical models. The average error across all cases is 1%, with the highest error observed in S2-A numerical models, reaching nearly 10% in predicting the ultimate load-bearing capacity.
5. Comparison with Current ACI Code Requirements
Because the ACI 318 code offers design methods for predicting the punching shear capacity of RC flat slabs, the results of the punching shear tests were compared with the values estimated using the punching shear principles outlined in ACI 318 [
64]. In the ACI code, punching shear strength can be regarded as the product of the shear strength and the area at the critical section, which is at a specified distance from the column face. The shear strength of the critical section can be represented as Equation (
2).
In this equation,
stands for the perimeter length of the critical section at a distance equal to half of the slab section from the column face,
d is slab thickness,
indicates the vertical component of all the effective prestress forces crossing the critical section, and
is the ratio of the long side over the short side of the column and needs to be the smaller of 3.5 or
. For an interior, edge, and corner column,
is 4, 3, and 2, respectively. Equation (
2) can be implemented for concrete with a compressive strength of less than 35 MPa; otherwise, Equation (
3) must be used.
The comparison of experimental test results and design values calculated using the ACI code is listed in
Table 7. In this table, the ACI code presents a conservative estimation for RC slabs and is more dominant in the B series, which strengthened the RC slabs with CFRP sheets. In addition, the ACI code methods are able to estimate the punching shear of the unstrengthened RC slab with reasonable accuracy. However, the accuracy of these methods declines with distancing the column from the centerline. In other words, ACI equations cannot precisely estimate the punching shear capacity of an RC slab with an edge column. The amount of error in predicting the punching shear of unstrengthened edge and corner slab–column connections are 42% and 14%, respectively. The accuracy of the ACI code prediction increased with utilizing CFRP sheets.
6. Parametric Study on the Location of Loading
By conducting experimental studies and establishing benchmarks for validating finite element models, a pathway has been cleared for subsequent investigations into the punching behavior of strengthened RC slabs. In order to evaluate the influence of loading location on the response of RC slabs, a series of numerical specimens were generated, leading to the development of 35 distinct finite element models. These models were analyzed, and the combined variations in loading location on peak load and ductility were determined for each slab. It’s important to highlight that these models exclusively pertain to the present investigation and should not be extrapolated to other comparable components.
Figure 10 shows the numerical specimens in a parametric study, where the variation in the location of the loading region is evaluated. The location of the loading region was treated as a variable in both the
x and
y directions. The variable’s range, measured from the center point of the loading plate, spans from 50 to 950 mm, with increments of 50 mm.
Figure 11a indicates that moving away from supports causes a reduction in peak load. Loading closer to the supports creates a leverage effect, allowing the slab to distribute the applied load more efficiently. Moreover, when the load is applied near the midspan of the slab, it generates a higher bending moment at that location. This increased bending moment requires the slab to undergo greater flexural deformation, potentially leading to an earlier onset of failure. In addition, when a load is applied in the middle, the shear forces are not effectively distributed, leading to higher localized shear stresses and, consequently, punching shear. This type of failure can lead to a reduced load-bearing capacity, especially if the slab is not adequately reinforced around the loading region. Therefore, the expected load-bearing capacity of the RC slab is reduced.
As can be seen in
Figure 11b the ductility increases toward the center of the slab, which has the least out-of-plane stiffness. The maximum ductility is about 10 and decreases to less than 3.5 when approaching the loading location toward the RC slab supports. This means that although the peak load and load-bearing capacity of the slab increase with approaching the loading location toward the RC slab supports, the ductility of the slab reduces, and a brittle failure of RC slabs can be expected. Brittle failure in an RC slab refers to a sudden and catastrophic failure with little or no warning before collapse. Unlike ductile failure, where structures show significant deformation and warning signs before failure, brittle failure occurs with minimal deformation, and the structure may fail abruptly. This type of failure results in a rapid loss of load-carrying capacity. This type of failure is typically characterized by the development of multiple, closely spaced cracks that lead to the structural elements’ inability to support loads, ultimately resulting in collapse.