Numerical Evaluation of the Opening Effects on the Reinforced Concrete Slab Structural Performance

A finite element method was used in this study to investigate the effects of openings on the resistive behavior of concrete slabs.&e presented modeling procedure is used to conduct numerical analyses on the response of reinforced concrete slab subjected to inplane monotonic loads in X (perpendicular to the beam) and Z (parallel to the beam) directions. Initially, the developed numerical model was calibrated and compared with laboratory results. In building this three-dimensional model, it is attempted to accurately model all nonlinear properties of steel and concrete materials as well as the interactions between them. &en, the behavior of bilaterally concrete slabs under different loads was investigated and used as a reference. Finally, the effect of openings under different loads on the strength of concrete slabs was studied. &e results confirm that openings have a great influence on the change of hardness, ductility, initiation and cracking path, and stress distribution under shear and gravitational loading. Moreover, by adding an opening inside the diaphragm, not only did beam and block flooring showmore fragile behavior, but also its strength and resistance decreased against lateral load. Given the results of the parametric study of the effect of layout, generally, its place became critical at the state that opening disturbed transmission of shear stresses to the collector beams. By adding the area of the opening and loading in X direction, the concentration of the tensile stresses (equivalent to main maximum stresses) was at the tensile edge as well as at the middle of the flooring around the opening. It is worth noting that an increase in the opening’s area caused the number of tensile stresses to be increased at the middle of the flooring. Meanwhile, the concentration of maximum compressive stresses which is equivalent to the main minimum stresses was at the compressive edge, started at the area of the collectors, and stretched to the edge of the opening. Among different layouts, X-1 and Z-3 states were more critical than other states. Considering openings with different layouts, X-1 and Z-3 have the most stiffness deteriorating and strength in such a way that stiffness deteriorating and strength were 39.93% and 37.89%, respectively, for Z-3model and 38.68% and 43.33%, respectively, for Z-3 model.


Introduction
Reinforced concrete (RC) slabs supported by RC columns are among the frequently used structural systems. e specific abilities of this system make them applicable alternatives for the construction industry. So, several studies about their simulation have been performed and design methods are proposed. Besides their advantages, RC slabcolumns systems have their own weakness in punching shear caused by brittle failure mode. To control punching shear failure, particular attention must be paid to design and implementation [1,2]. Many national and international construction codes have been investigated and several methods were introduced for calculating and implementing such systems [3]; fib Bulletin 14 [4]; Italian guide CNR DT200 [5]; Japan Society of Civil Engineers (JSCE) design recommendations [6]. By examining these codes, it can be stated that they are actually studying the same subject and the difference in expression is due to the different point of view [7]. Using bonded rebar in a concrete slab is a common method to increase the punching shear strength of the RC slab-column systems [8,9]. Furthermore, the inevitable presence of openings in these systems because of architectural and structural considerations makes the task of evaluating safe performance more complicated [10]. Hence, many analytical studies, numerical researches, and experimental efforts have been done on the performance evaluation of RC slab subjected to many various different loading types [11,12]. Most scholars believe that rebar-concrete interaction in RC slab can be simulated with the assumption of the perfect bond at the interface between them [13,14]. Some researchers investigate the effect of bondingdebonding in simulating the structural performance of RC slabs [15,16].
Finally, in 2016, Khajehdehi and Panahshah [17] conducted studies on the effect of openings on the in-plane structural behavior of RC floor slabs and concluded that the collapse mechanism of slabs with an opening significantly differed from the mechanism of slab collapse without opening. Furthermore, the presence of openings in floor diaphragms for architectural features, staircases, and elevator shafts is sometimes inevitable [18]. ese openings will result in diaphragm stiffness reduction and can decrease the load-carrying capacity of the member [19,20]. ese types of diaphragms are usually designed to ignore opening effects. erefore, their true response may be different from what is assumed. In other words, the presence of openings makes the behavior of floor diaphragms significantly more complicated and unpredictable [21]. A number of researchers have evaluated the structural behavior of RC slabs with openings in them. However, slabs were only subjected to out-of-plane gravity loads, and the effect of in-plane loads was not considered [22][23][24][25]. A number of researchers have evaluated the structural behavior of RC slabs with the presence of opening and out-of-plane loading. Radik et al. [26], Choi et al. [19], and Floruţ et al. [21] studied the effect of various strengthening methods on improving the load-bearing capacity of slabs with the presence and absence of openings using FRP, GFRP, and FRC [27]. Khajehdehi and Panahshah conducted a numerical study of the effect of the opening on the floor slab concluding that the presence of an opening plays an important role in determining the behavior of the slab plates [17].
In this study, debonding of rebar from concrete surfaces is taking into account. It is obvious that validated FE simulations may cause accurate results on the response of RC slabs. Also, it is obvious that, in numerical studies, one should use logical assumptions to get closer to the finite element model that the numerical models of reinforced concrete slabs are no exception. Until now, most studies have focused on the behavior of unilateral reinforced concrete slabs, and more limited studies have been conducted on the impact of openings on the structural performance of bilateral concrete slabs [28,29]. So far, many studies have been carried out on this issue due to the high capability of such a structural system. However, the number of comprehensive studies that have modeled and investigated the purely nonlinear behavior of materials coupled with the precise definition of the interaction between its components is limited. us, the 3D nonlinear finite element model of RC slab considering the rebar-concrete interaction is used to simulate such structural elements [30]. In a structure, the need for an opening in the roof system (slab) for architectural features, stairs, elevator shafts, pipelines, and utility ducts is inevitable. ese openings will in some status lead to a decrease in the stiffness of the roof system and collapse load due to their dimensions relative to the dimensions of the slab [19,31]. In designing these types of slabs, the effects of openings are usually ignored. erefore, their actual response may be different from that assumed. Radik et al. [26], Choi et al. [19], and Floruţ et al. [21] investigated the effectiveness of different strengthening methods on improving the load-carrying capacity of slabs with and without openings using GFRP, FRC, and FRP, while slabs were only subjected to out-of-plane loads and again the effect of in-plane loads was not considered. Zhang et al. [32] conducted a numerical study on the effect of openings on floor slabs and concluded that the presence of openings plays an important role in determining the in-plane behavior of the slabs. In the present research, a two-way RC slab system tested by Durucan and Anil is used to validate the numerical fully nonlinear FE model [33]. e main scope of the performed numerical this study is the simulation of direct and indirect impact of horizontal construction openings on expected structural functions of bilateral reinforced concrete slabs. In this regard, the influence of different positions and sizes of openings on the deformation and hardness of reinforced concrete slabs is studied in detail and the results are compared with each other to obtain a comprehensive view. [33]). In this study, the shear behavior of the slab interior connection in addition to FRP reinforced opening with bar was studied using a set of experiments. In this study, a sample without an opening and 8 other samples with openings were along the sides and diagonal. In all samples, slab dimensions have been 2000 × 2000 × 120 mm and bar dimensions were 200 × 200 mm as shown in Figure 1. In this research, Durucan and Anil's modeling in ABAQUS finite element software was used due to the lack of valid experimental researches on beam and block floorings. en after results' confirmation, a model was provided and studies of flooring were designed for beam and block according to Iran's regulation for finite element software [34].

Benchmark Experimental Study (Durucan and Anil
During the experiment, a hydraulic jack on the column imposes a fixed axial load on the column which is indicative of the gravity load obtained from the frame analysis [35].  Tables 1 and 2, respectively. Meanwhile, the way the experiment is conducted is indicated schematically as shown in Figure 2. In this experiment, a concrete reinforced slab was placed on IPE steel pieces and then a hydraulic jack imposed axial load to the column [36]. Table 3 shows the experimental results. As it can be seen, the displacement in the sample without opening is 40 mm toward the center which is equivalent to 193 kN [37].

Common Component Modeling of the Samples.
For all samples, model geometry including components such as concrete slab, reinforcements, and rigid supports and loading in the software separately, as shown in Figure 3, were modeled, and finally, the samples were made completely in the software by putting these components along each other. Since constraints and interactions among the components have a great influence on the results of the analysis, it was applied carefully and cautiously.

Types and Sizes of the Elements.
e concrete column was modeled by three-dimensional elements accessible at the library of numerical software [38]. e concrete slab was meshed using C3D8R elements as shown in Figure 4. Reinforcements and supportive components were modeled using T3D2 and C3D8R truss elements, respectively. Dimensions of the elements were considered as 20 mm for slab and steel reinforcements and 100 mm for support (see Figure 4).

Interaction and Constraints Conditions among Different Components.
Where concrete contacts steel, the surface-to-surface constraint was used without friction coefficient which is shown in Figure 5 (contact of supportive sheets with concrete slab). For the rigid sheets, a point that is called reference point is defined too that all degrees of freedom for these sheets are affected by this point. A Tie constraint is defined; when two surfaces adhere to each other, one should be determined as Master surface and another as Slave surface. As it is mentioned in the software manual, it is better that the surface which has finer materials be determined as a Slave surface and the dimensions of its elements be smaller than those of the Master surface in order to avoid numerical errors. erefore, when the constraint is used, points' displacement of the Slave surface is obtained through the points' displacement of the Master surface and indeed relative slip between these two surfaces is neglected.

Boundary Conditions.
In laboratory tests, the connection of the support is to fixed ground. erefore, in its modeling in the software, all rotational and transitional degrees of

Comparing Experimental Results and Numerical
Outputs. At first, finite component samples were studied and validated to be sure that intended samples are modeled correctly and the given parameters have reasonable values, it was seen that the numerical results and laboratory results are in good agreement, and then a series of parametric studies were conducted to investigate the efficacy of different components of the connection. Figure 6 shows the comparison of the force-displacement response of the finite element method with the experimental method.
As it can be seen from Figures 6(a) and 6(b), the results of numerical analysis and laboratory results are in good agreement, indicate completely identical behavior, and experience the maximum strength at nearly 40 mm, but at displacements higher than this amount, strength decreased due to concrete cracking. Figure 7 compares concrete slab cracking in the experimental study and numerical model. It can be seen that most cracks are observed in the area in the middle of the concrete slab.

FEM of Reinforced Concrete Slabs.
To provide a model for beam and block flooring, concrete slab and steel truss were modeled in ABAQUS software as shown in Figure 8. As      Shock and Vibration shown in this figure, the concrete component was modeled as solid and mashed with C3D8R elements. is element was a hexagonal 8-node one and reduced integration was used in it. Reinforcements were modeled as wire and were meshed with two-node three-dimensional truss elements (T3D2). Dimension for each element was considered 20 mm.

Loading and Boundary Condition for Numerical
Model. Since the cross section of the model was different in two X and Z directions, all analyses were done in two directions. Boundary conditions of the model were applied in the way that transitional degrees of freedom were closed at two edges of the slab which means for analysis in X direction and for analysis in Z direction in two edges of the slab, U x � U y � U z � 0. erefore, a volumetric load was applied incrementally with a linear trend. Applied load to the diaphragm was very higher than the regulatory amount. e cause for the higher load applied is to allow the diaphragm to pass linear area and enter to plastic range in order to study its behavior carefully. Boundary conditions and how models were loaded in two X and Z directions are indicated in Figure 9.

Studying the Diaphragm Behavior in the Numerical
Model without Opening. By applying lateral body force, distribution of main maximum plastic strains (see Figure 10) which is equivalent to cracking strains is in such a way that, by loading in direction in the model, crack expansion is perpendicular to the direction of loading and cracking is more in the place of the beam to slab connection. Meanwhile, the rupture mechanism is completely visible in the area where it plays the role of tensile, due to the existence of axial force. Meanwhile, cracks at the support's areas are sheering which is less dominant in Figure 10 because the section of the diaphragm of this mechanism is strong.
Meanwhile, in the model of loading in Z direction, crack distribution is centered more in the collector area which is because the diaphragm section is small in this direction. Meanwhile, after terrible cracks in the collector area, crack expansion was obliquely towards inside the diaphragm and these cracks penetrated more deeply at the collector area. It is worth mentioning that, as shown in Figure 10, in the models with loading in X and Z directions, the maximum amount for main maximum plastic strains which is equivalent to cracking strain is 0.02561 and 0.1041, respectively. It means that, in this area, cracks were opened widely due to the weakness of the diaphragm section to bear the loading direction. Of course, this amount of strain for far displacement of 15 mm is not far-fetched. Figure 11 shows the distribution of the main minimum plastic strain (equivalent to corrosion strains) from the strain in maximum strength; in the model with volumetric loading in X direction, distribution of corrosion is focused on the connection points of the beams to slab. Meanwhile, as it is clear on compression chord, strains are purely compressive. Meanwhile, cracks at the support's areas are sheering which are less dominant in Figure 11(a) because the diaphragm section is strong. Meanwhile, in the model with loading in Z direction, the distribution of corrosion strains of the concrete was more concentrated at the collector area because the section of the diaphragm is small in this direction as shown in Figure 11(b) too. It should be noted that the expansion of corrosion was less than cracking which can be seen clearly by comparing Figures 10 and 11.
Meanwhile, the distribution of main maximum and minimum stresses at displacement equivalent to maximum strength for models with loading in Z and X directions is shown in Figures 12 and 13 which indicate the stress transmission mechanism in the diaphragm in these two directions. As it is clear in these figures, the mechanism of collector bars and chord is clearly visible and the highest stress was created in these places. erefore, bars located in these areas can tolerate the forces resulting from the mechanism and collectors.

Reaction of Lateral Force-Displacement of Diaphragm without
Opening. Force applied to the diaphragm linearly and incrementally was extracted against transformation at the middle of the opening by the software. e amount of this force against the transformation of the middle point of the slab is indicated in Figure 14.
As shown in Figure 14, the behavior of diaphragm with beam and block section in Z direction in the direction of beam's axis, bean, and block flooring indicated more fragile behavior which is mainly because of the weakness of the diaphragm in stress transmission at the collector area. e amounts of maximum force, displacement equivalent to maximum force, and hardness of the diaphragm are indicated in Table 4.

Numerical Study of Opening Effects on RC Slab.
ere is no doubt that the existence of voids inside the diaphragm can affect the capacity to create lateral force in the diaphragm. e area and location of the opening are things that its critical place on the diaphragm performance and the amount of stress created inside it and especially around the opening should be determined. Particularly, investigating these parameters at beam and block flooring as the most common details in Iran has great importance. In this section, the location of the opening and its area at this kind of flooring has been studied parametrically.

Modeling Hypothesis of the Opening in Diaphragm.
In order to model the opening, the opening model was created by creating a square opening inside it in the way that both slabs and reinforcements were cut. Meanwhile, due to the complexities of modeling, the hypothesis of not slipping the reinforcements inside the concrete was considered by taking embedded region interaction into account. Meanwhile, no reinforcement was used around the opening. Similar to the previous section, loading on the diaphragm model in addition to opening was conducted and studied in two Z and X directions. For parametric studying, firstly the layout of the opening was studied, then the most critical      Shock and Vibration place for the opening was determined, and at the same critical place, the area of the opening was analyzed.

Parametric Study of the Opening Layout inside Diaphragm of Beam and Block Flooring.
To model the layout of opening at the diaphragm, a square opening of 2 m × 2 m was used which is equivalent to 16% of the diaphragm. e location of the opening and the name of the corresponding model have been shown in Figure 15. What is worth mentioning in Figure 15 is that, considering symmetry, X-4, X-5, and X-6 models can be ignored at the loading state in X direction as these models are, respectively, correspondent with X-2 and X-8 models. Meanwhile, at the state of loading in Z direction, Z-6, Z-7, and Z-8 models can be neglected as they are correspondent to Z-2, Z-3, and Z-4 due to being symmetric. Applied force to the diaphragm, which was applied linearly and incrementally, was extracted against the transformation of the middling of the opening through software. e amount of this force against the transformation of the slab's middle point is shown in Figure 16(a) for the models with loading in X direction and in Figure 16(b) for the models with loading in Z direction.
By comparing charts of Figures 16(a) and 16(b), it can be found that, by adding the opening inside the diaphragm of the beam and block flooring, not only does diaphragm show more fragile behavior but also its strength and hardness decrease against lateral force too. Meanwhile, it can be clearly seen by comparing these two figures that loading diaphragm in Z direction shows more fragile behavior. In such a way that, after reaching the diaphragm to its maximum capacity, it is accompanied by a sudden decrease in strength. Meanwhile, by comparing Figures 16(a) and 16(b) more closely, the results indicated in Tables 5 and 6 can be observed.
For more tangible and accurate study of the critical state for the place of the diaphragm, values calculated for stiffness deteriorating and strength are provided in Tables 5 and 6 as well as in the bar charts drawn in Figures 17(a) and 17(b).     As it can be seen from Figures 17(a) and 17(b), X-1 and Z-3 models have the highest stiffness deteriorating and strength which are more critical compared to other models. Meanwhile, the criticality of this opening layout from stress and strain distribution can be seen in the next session. According to the results of the parametric study, opening layout location is critical in two loading modes in X and Z directions, when the opening disrupts the shear stress transfer path to the collector beams. Meanwhile, the X-1 and Z-3 modes are more critical than the others. is is due to the high degree of deterioration and diaphragm resistance shown in Tables 5 and 6 as well as Figure 17.

Parametric Study of the Openings with Various Areas inside Diaphragm of Beam and Block
Flooring. In order to investigate the effect of diaphragm area on hardness,  strength, stress transfer mechanisms, cracking and crushing, different critical areas of opening were investigated at these critical locations. e names of these models and how they are applied to the models are shown in Figure 18. e force applied to the diaphragm linearly and incrementally was extracted from the software against its midspan deformation. e magnitude of this force versus slab midpoint deformation for models with loading in X direction is shown in Figure 19 and for models with loading in Z direction in Figure 20.
It can be seen that, by increasing the area of the opening inside the diaphragm, its resistance and hardness against lateral load decrease. It is also apparent by comparing the two shapes that the Z-diaphragm loading shows a slower and more brittle behavior. So once the diaphragm reaches its maximum capacity, there is a sudden drop in resistance. A closer look at Figures 19  and 20 shows the results in Tables 7 and 8.
In order to more accurately and precisely investigate the effect of the increase of the opening area on the final diaphragm hardness and resistance, the calculated values for the stiffness and resistance decline shown in Tables 7 and 8 are shown in Figure 21.
As can be seen from the previous figures, the diaphragm resistance and stiffness decrease with increasing opening area in numerical models. For models with loading in X direction, the hardness and strength      Model Name Z-1-3 X-1-1 X-1-2 Model Name  deterioration rates are approximately the same, but for models with loading in Z direction, the hardness deterioration is greater than the resistance deterioration.

Conclusion
In this study, the finite element method was used for three-dimensional modeling and nonlinear analysis on the diaphragm of the beam and block flooring. Meanwhile, the following conclusions are arrived at by studying this kind of diaphragm, with and without opening considering different layouts. By applying volumetric lateral force, distribution of main maximum plastic strains which are equivalent to cracking strains is in a way that, in the model with loading in X direction (perpendicular to the beam), expansion of the crack was perpendicular to the direction of loading. Meanwhile, at the areas that play the role of tensile edge, the rapture mechanism is completely visible due to the existence of the axial force. Meanwhile, cracks at the palace of the support are sheering which is less dominant because the section of the diaphragm is strong. Meanwhile, at the model with loading in Z direction (parallel to the beam), the distribution of the cracking is focused more at the collector area because the section of the diaphragm is small in this direction. Meanwhile, after terrible cracks at the collectors' area, expansion of the crack was obliquely towards inside the diaphragm. Meanwhile, at the area of tensile edge, the distribution of the cracking is completely visible in a tensile manner and the cracks penetrate deeply at the area near to the collectors. At the models with loading in X and Z directions, the maximum amount of main maximum plastic strain which is equivalent to the cracking strain is, respectively, 0.02561 and 0.1041. It means that, due to the weakness of the sections' diaphragm to bear loading in Z direction, this area is forced with the openness of the crack width at a high amount. Of course, this amount of strain is not far-fetched for far displacement of 15 mm for the diaphragm. By adding an opening inside the diaphragm, not only did beam and block flooring show more fragile behavior but also its strength and resistance decreased against lateral load. Considering the force-displacement charts for the models with opening with different layouts, in which opening is at proximity to the collectors, X-1 and Z-3 have the most stiffness deteriorating and strength in such a way that stiffness deteriorating and strength were39.93% and 37.89%, respectively, for Z-3 model and 38.68% and 43.33%, respectively, for Z-3 model. By studying the distribution of the maximum and minimum plastic strain (equivalent to cracking and corrosion) in the models with opening with different layouts, at loading in X direction, the most cracks occurred at the area in which beams were connected to slab as well as at the corners of the opening. Meanwhile, cracks became critical when the opening was placed in the collectors' area. e distribution of corrosion had been identical to cracking distribution too but its amount was less in the diaphragms. Meanwhile, for the models with loading in Z direction, when the opening was placed in the collectors' area (as the critical state), the most concentration of main maximum plastic strains which was equivalent to cracking was at the collector area. e most concentration of the main minimum plastic strains (equivalent to corrosion) in this state was similar to that at the area of the collector too where the slab tolerated high stresses. e amounts of main minimum plastic strain reached corrosion around the opening proximity to the collector area. Given the results of the parametric study of the effect of layout, generally, its place became critical at the state that opening disturbed transmission of shear stresses to the collector beams. Among them, X-1 and Z-3 states were more critical than other states.
By studying stiffness deteriorating and strength for the diaphragms with the opening with different areas, it was seen that, by increasing the opening area from 0% to 25% area of the diaphragm, stiffness deteriorating and strength increased but the rate of their raising decreased gradually. By studying the distribution of minimum and maximum plastic strain (equivalent to cracking and corrosion) at the models with opening with different layouts, by applying load in X direction to the diaphragm with opening, the most cracking concentrated at the top of the opening and after that at its lower corner by changing the area. Meanwhile, it cracked at the area where the slab was connected to the beams. But it is worth noting at this state of openness that all connections of the beam to slab around the opening are faced with many cracks and corrosion and worsen by increasing the area. e most corrosion was concentrated at the area of tensile edge and at the top of the opening and the highest concentration of main minimum plastic strains which is equivalent to corrosion is around the corners of the opening too. By applying load to the diaphragm with the opening in Z direction, at 15 mm displacement for the diaphragm, the most concentration of main maximum plastic strains, which is equivalent to cracking, was at the area of the collector which was caused due to the weakness of the slab's section at the area connected to the collector beams that slab tolerated high stresses. After that at top of the opening and its lower corner, the main maximum plastic strain arrived at high amounts. Meanwhile, the highest concentration of the main minimum plastic strains equivalent to cracking was at the area of the collectors where slabs bear high stresses. Meanwhile, parts of the opening's corner, which are in the direction of the collectors cracked. It is worth mentioning that, by increasing the opening's area, the amounts of cracking and corrosion increased, and the number of plastic strains became more critical.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request (amir.reza.eskenati@upc.edu).