A Modified Compression Field Theory Based Analytical Model of RC Slab-Column Joint without Punching Shear Reinforcement

RC slab–column structures are widely used because of the advantages of small space occupation for horizontal members, easy construction and good economy. However, slab–column joints are prone to punching shear failures, which deteriorates structural safety. This paper provides an analytical model to predict the punching shear capacity of the RC slab–column joint. A database of 251 test results is established for the shear punching capacity of slab–column joints without punching shear reinforcement. The performance of existing design codes in predicting the shear resistance of slab–column joints is investigated and compared based on the database. Then, based on the modified compression field theory (MCFT) model, an equation for calculating the punching shear resistance of slab–column joints without punching shear reinforcement is established. The prediction results of the analytical model are enhanced by using the regression analysis method. The model proposed in this paper is based on both reliable theoretical and the summary of a large number of test results, which has higher prediction accuracy than the design codes.


Introduction
The RC slab-column structure has no beams between the columns, and thus, is a form transfer load directly from slabs to columns [1,2]. The slab-column structure also has the advantages of flexible layout, full use of height space, fast construction speed and low cost. It is widely used in warehouses, underground garages, bridges and ports [3,4]. However, the RC slab column structure may undergo punching shear at the slab column joints under bending and shear loads. With small slab-column deformation, sudden punching shear failure can occur which leads to serious problems [5][6][7][8]. Brittle punching shear failure of RC slabs is usually caused by high shear stresses around the slab-column joint [9]. Due to the mechanisms of brittle failure at slab-columns joints, it is easy to cause progressive collapse, which affects the integrity of the structures [10]. On 24 June 2021, a condominium building collapsed [11] in which 98 people were killed in the accident (Figure 1). Scholars have carried out a lot of research in the hope of improving the understanding of the punching shear [9,12,13]. However, the mechanism of punching shear of slab-column structures is complex; therefore, the punching shear problem is still open for investigation. Many researchers started with theoretical approaches to find mechanical explanations for punching shear and predicted the punching shear resistance of slabcolumns by using various theoretical models. Bazant and Cao [14] introduced the fracture mechanics model to study the punching shear process of slab-column and Hallgren [15] also used nonlinear fracture mechanics to study the effect of size effect on high-strength concrete. Yankelevsky et al. [16], Nielsen [17], Johansen [18], Silva [19], Cai and Lin [20], studied the punching shear of slab-column based on the theory of plasticity. Based on seismic mitigation strategies and cost-benefit criteria, the deformation and forces at slabcolumn joints need to be analyzed and quantified [21]. However, the above models give poor consideration to the deformation capacity of slab-columns [21]. Based on elastic thinplate theory, Long [22] presented an analytical model for punching shear capacity and summarized the criteria for different failure types of slab-column joints. Kinnunen and Nylander [23] proposed a model to calculate the ultimate rotation of the slab to study the punching shear capacity of RC slabs. The object of study for this model is sector elements in the radial flexure cracks. This model considers that when the radial inclined compression stress and the tangential compression strain near the column reach critical values the punching occurs. Then, the size effect was taken into account and a solution for the ultimate rotation was derived by Broms [12,24] to make the model more accurate. According to nonlinear fracture mechanics, Hallgren [15] proposed a model that the height of the concrete compression zone was used to derive the effects of size effect for the punching capacity of high-strength RC slab-column joints. Based on the conclusion of Walraven [25], assuming that the width of a critical shear crack is proportional to the rotation of the slab, Ruiz and Muttoni [2] proposed critical shear crack theory (CSCT). The effects of crack width and aggregate interlocking on the punching resistance of the slabcolumn joints were considered in the model. The accuracy of the relationship between the rotation angle ( ϕ ) and the internal forces in CSCT was verified by experiments. Then, the aggregate was considered and CSCT was improved by Guandalini et al. [26]. The Many researchers started with theoretical approaches to find mechanical explanations for punching shear and predicted the punching shear resistance of slab-columns by using various theoretical models. Bazant and Cao [14] introduced the fracture mechanics model to study the punching shear process of slab-column and Hallgren [15] also used nonlinear fracture mechanics to study the effect of size effect on high-strength concrete. Yankelevsky et al. [16], Nielsen [17], Johansen [18], Silva [19], Cai and Lin [20], studied the punching shear of slab-column based on the theory of plasticity. Based on seismic mitigation strategies and cost-benefit criteria, the deformation and forces at slab-column joints need to be analyzed and quantified [21]. However, the above models give poor consideration to the deformation capacity of slab-columns [21]. Based on elastic thin-plate theory, Long [22] presented an analytical model for punching shear capacity and summarized the criteria for different failure types of slab-column joints. Kinnunen and Nylander [23] proposed a model to calculate the ultimate rotation of the slab to study the punching shear capacity of RC slabs. The object of study for this model is sector elements in the radial flexure cracks. This model considers that when the radial inclined compression stress and the tangential compression strain near the column reach critical values the punching occurs. Then, the size effect was taken into account and a solution for the ultimate rotation was derived by Broms [12,24] to make the model more accurate. According to nonlinear fracture mechanics, Hallgren [15] proposed a model that the height of the concrete compression zone was used to derive the effects of size effect for the punching capacity of high-strength RC slab-column joints. Based on the conclusion of Walraven [25], assuming that the width of a critical shear crack is proportional to the rotation of the slab, Ruiz and Muttoni [2] proposed critical shear crack theory (CSCT). The effects of crack width and aggregate interlocking on the punching resistance of the slab-column joints were considered in the model. The accuracy of the relationship between the rotation angle (ϕ) and the internal forces in CSCT was verified by experiments. Then, the aggregate was considered and CSCT was improved by Guandalini et al. [26]. The theoretical models mentioned above study the punching shear behavior of slab-columns from multiple perspectives. The prediction of punching shear becomes more accurate, but the models and equations are relatively complex.
In RC slabs without punching shear reinforcement, the brittle punching failure has a higher probability [6], and thus, the slab-column joints without punching shear reinforcement are the research subjects in this paper. As far as we are concerned the punching shear (within vertical load) of RC slab-column joints without punching shear reinforcement are similar to shear diagonal tension failure of RC beams with a large shear-to-span ratio, and this failure is caused by the direct stress and shear stresses in the concrete. Therefore, based on the Modified Compression Field Theory (MCFT), a punching shear equation is proposed by transforming the three-dimensional mechanics into two-dimensional plane mechanics. With detailed specimen parameters, the equation has a relatively simple form and can achieve relatively accurate prediction.
The purpose of this paper is to establish an analytical model for punching shear capacity of RC slab-column joint that is based on justified theoretical deduction and the experimental database. In Section 1, a database that includes 251 slab-column joints tests are established. The database collects and collates the specific parameters of the slabcolumn joints without punching shear reinforcement. In Section 2, a total of five design codes are selected in this paper. The performance of existing design codes in predicting the shear resistance of slab-column joints are investigated and compared with the experimental results from the database. In Section 3, the Modified Compression Field Theory (MCFT) is applied to give the mechanical basis of punching shear behavior in interior slab-column joints without punching shear reinforcement. An equation for the punching shear resistance of slab-column joints without punching shear reinforcement is established. In order to improve the proposed model, parameters in the model are determined by regression analysis. The proposed model is compared with the existing methods provided by other design codes, verifying the correctness and effectiveness of the model.

Experimental Database
According to a literature review, a database of 251 tests is established. The following criteria are used in collecting tests data: (1) The selected tests are RC slab-column joints (within vertical load), and the shape of the specimen is rectangular, which can exclude the influence caused by the shape; (2) The flexure reinforcement is arranged at the bottom of the slab, and there is no other punching resistance element; (3) The position of the column is in the middle of the RC slab.
The parameters recorded in the database include the side length of the slab (B); the side length of the load area or column (c); the effective height of slab (h 0 : the distance from the top of the ultimate compression zone to the center of flexure reinforcement); the ratio of punching-span (λ: ratio of half of the span to the effective height); the cylindrical compressive strength of concrete ( f c ); the axial tensile strength of concrete ( f t ); the reinforcement ratio (ρ); the yield strength of reinforcement ( f y ); the experimental punching shear capacity (V t ). The concrete compressive strength of some tests is 150 × 150 mm cube compressive strength f cu,k . Thus, it needs to be converted to the cylindrical compressive strength f c through Equation (1). Figure 2 demonstrates the distribution of the above mentioned parameters. As depicted in Figure 2, the range of the side length is mainly distributed in 125-1925 mm. The range of side length of the loading area or column is mainly distributed in 150-280 mm. It clearly shows that most specimens are not full scale tests. The effective height of the specimen is around 115 mm. Moreover, 75% of the specimens exhibit the slenderness ratio h/B = 0.07-0.10. This value conforms to international standards and ensures that the slab's flexural behavior is acceptable [27]. In order to study the effects of size effect, another 25% of the specimens which the slenderness ratio h/B are more than 0.10 or less than 0.07 are also collected in the database. Zaghlool and Paiva [28] reported that if the rotation constraint is imposed on the boundary of the slab-column joints, the punching resistance capacity can be increased at least 10%. In order to reduce the effects of boundary constraints, the tests selected are mostly simply-supported. The range of f c is 11-78 MPa, mostly less than 60 MPa. The range of f t is 1.21-6.01 MPa. The specimens selected in the database are RC slab-column joints, and the proportion of high-strength concrete is small. Some scholars [3,8,29] have found that the reinforcement ratio has a great influence on the punching capacity of the slab-column joints. Therefore, the range of ρ selected is wide in the database, which is in the range 0.33% to 2.73% and mostly concentrated in 0.6−1.8%. The range of V t is 105-1041 kN, where 91% of the specimens are less than 700 kN.
The another 25% of the specimens which the slenderness ratio h/B are more than 0.10 or less than 0.07 are also collected in the database. Zaghlool and Paiva [28] reported that if the rotation constraint is imposed on the boundary of the slab-column joints, the punching resistance capacity can be increased at least 10%. In order to reduce the effects of boundary constraints, the tests selected are mostly simply-supported. The range of ' c f is  MPa, mostly less than 60 MPa. The range of t f is 1.21-6.01 MPa. The specimens selected in the database are RC slab-column joints, and the proportion of high-strength concrete is small. Some scholars [3,8,29] have found that the reinforcement ratio has a great influence on the punching capacity of the slab-column joints. Therefore, the range of ρ selected is wide in the database, which is in the range 0.33% to 2.73% and mostly concentrated in 0.6%−1.8%. The range of t V is 105-1041 kN, where 91% of the specimens are less than 700 kN. The

Punching Shear Design Provisions
As for the design codes of RC structures, most of them determine punching shear failure by the nominal stress in the critical section. Although various aspects are considered in the design codes, such as the shape of punching shear critical section, size of the critical section, nominal stress calculation, reinforcement ratio and other controlling factors. There are five design codes shown in Table 1. In summary, the design codes selected in this paper can be split into two categories as follows: (1) design codes that do not take into account the effect of reinforcement ratio, including "Code for

Punching Shear Design Provisions
As for the design codes of RC structures, most of them determine punching shear failure by the nominal stress in the critical section. Although various aspects are considered in the design codes, such as the shape of punching shear critical section, size of the critical section, nominal stress calculation, reinforcement ratio and other controlling factors. There are five design codes shown in Table 1. In summary, the design codes selected in this paper can be split into two categories as follows: (1) design codes that do not take into account the effect of reinforcement ratio, including "Code for Design of Concrete Structures

Design Codes
Punching Shear Strength The section height influence coefficient β h = 1.0; η 1 is the influence coefficient of the shape of loaded area; η 2 is the influence coefficient (u m /h 0 ); u m is the critical perimeter; (The critical section does not consider the angle). β s is the ratio of long side to short side (2 ≤ β s ≤ 4 for GB50010 [49]); α s is the column position influence coefficient, which is taken as 40, 30, 20 for the inner, edge and corner column; α se is taken as 1.15, 1.4, 1.5 for the inner, edge and corner column. α se is taken as 4, 3, 2 for the inner, edge and corner column; ξ is the size effect factor, which is taken 1 (normal weight concrete) or 0.75 (lightweight concrete); k is the size effect coefficient. Table 2 clearly shows the difference in the parameters included in the punching shear design. The main influencing parameters include concrete strength f c and f t , flexure reinforcement ratio ρ, the critical section u m , effective height h 0 , the position of the column and size effect. The punching shear capacity is related to the concrete tensile strength directly. Thus, the tensile strength of concrete f t is used in the equation of punching shear strength in GB50010 [49]. ACI318 [50], EC2 [52], CSA [51] and JSCE [53], take 0.33 f c , ( f c ) 1/3 , 0.38 f c and 0.2 f c as the tensile strength of concrete in the equation of punching shear strength. In addition, only EC2 [52] and JSCE [53] include the effect of flexure reinforcement ratio, and the punching shear capacity is proportional to ρ 1/3 . All design codes hope to increase the punching shear capacity by increasing the effective height h 0 of the slab. Compressive strength of concrete Figure 3 shows the ratio of the experimental punching shear capacity (V t ) to the value calculated by the design codes (V e ) for RC slab-column joints. The tests are divided by the ratio of punching-span (λ), and the ordinate is V t /V e . The data analyzed include average value, standard deviation, coefficient of variation, the maximum and minimum of V t /V e and the percentage of exceeding the average (Table 3).

Evaluation of Design Codes
to resisting punching shear will become limited when 2 % ρ > [57]. The discretene the average value of the calculation results of CSA [51] Figure 3 shows that the calculation results of the equation used in ACI318 [50] are mostly larger than the experimental results. Table 3 shows that the average value of V t /V e (1.4369) is large and the coefficient of variation is small. One of the reasons [21] for these performances is that the flexure reinforcement ratio is not considered in ACI318 [50]. The punching shear strength is typically controlled by V ACI3 [13]. However, Some researchers [54][55][56] have found that if the size of the load area is considerably greater than the height of the slab, V ACI1 and V ACI2 are more appropriate. JSCE [53] and EC2 [52] have the same characteristics, the coefficient of variation is large and the average value of V t /V e is smaller, which means the discreteness of predictions is larger than ACI318 [50], GB50010 [49] and CSA [51]. Besides, only JSCE [53] and EC2 [52] consider the effect of flexure reinforcement ratio and the contribution of the yield strength of the flexure reinforcement to resisting punching shear will become limited when ρ > 2% [57]. The discreteness and the average value of the calculation results of CSA [51] and GB50010 [49] are small. However, they both have big differences between the maximum of V t /V e and the minimum of V t /V e , which means GB50010 [49] and CSA [51] have the limitations to predict the punching resistance for different size slab-columns.

Modified Compression Field Theory
Modified compression field theory (MCFT) is a method to solve the shear problem of reinforced concrete members proposed by Vecchio and Collins [58,59]. Based on the average stress and strain of cracked RC elements, MCFT establishes balance, compatibility and constitutive equations. Punching failure of RC slab-column joints essentially is shearing failure [16]. Therefore, in this paper, the basic equations established by MCFT are used to calculate the stress of RC cracked elements under the plane shear and axial force. Then, the calculation equation for the punching shear capacity of the slab column is established.
According to the stress state of the cracked concrete element, the stress state of the reinforcement element and the average stress Mohr circle shown in Figure 4, the balance equations are established as follows: Where cx  According to the geometric deformation conditions of the cracked concrete element and the average strain Mohr circle (Figures 5 and 6), the strain compatibility equations of the cracked concrete element can be established: The stress of concrete element at the crack is as follows: where σ cx and σ cy are the average stress of the cracked concrete element in the x direction and y direction, respectively, ρ x and ρ y are the horizontal and vertical reinforcement ratios, respectively, σ sx and σ sy are transverse and longitudinal reinforcement stress, respectively, σ f is the principal tensile stress of the cracked concrete element, σ k is the principal compressive stress of the cracked concrete element, τ and τ c are the shear stress and average shear stress of the element, respectively. According to the geometric deformation conditions of the cracked concrete element and the average strain Mohr circle ( Figures 5 and 6), the strain compatibility equations of the cracked concrete element can be established: where ε 1 is the average principal tensile strain of the vertical crack, ε 2 is the principal compressive strain of parallel cracks, ε x and ε y are the average strains in the x and y directions, γ c is the average shear strain. The softening effect occurs when the concrete reaches its tensile strength, the following stress-strain relationships (Figure 7) are used [58] as: where E c is the Young's modulus concrete, ε 0 is the peak strain of the concrete, σ gmax is the peak compressive stress of softened concrete.     The softening effect occurs when the concrete reaches its tensile strength, the following stress-strain relationships (Figure 7) are used [58] as: Where c E is the Young's modulus concrete, 0 ε is the peak strain of the concrete, max g σ is the peak compressive stress of softened concrete.  The stress-strain relationship of the reinforcement (Figure 8) is elastic-perfectly plastic as follows: where E s , ε sx , f sy is the Young's modulus, strain, and yield strength of longitudinal reinforcement. The RC slab-column joint without punching shear reinforcement presents a rigidplastic failure form when punching shear occurs [60], and the failure area forms a punching cone (Figure 9). In this paper, the effect of the parts which directly contact the surface of the column (part Ⅰ) is not considered. Only consider the part of the cone connected with the column (part Ⅱ) and 4 prisms can be taken as the free-body. One of the prisms free-body is shown in Figure 10. Then integrate the stress on the surfaces (part Ⅲ) of the prisms to represent the punching shear strength.  The RC slab-column joint without punching shear reinforcement presents a rigidplastic failure form when punching shear occurs [60], and the failure area forms a punching cone (Figure 9). In this paper, the effect of the parts which directly contact the surface of the column (part I) is not considered. Only consider the part of the cone connected with the column (part II) and 4 prisms can be taken as the free-body. One of the prisms free-body is shown in Figure 10. Then integrate the stress on the surfaces (part III) of the prisms to represent the punching shear strength. The RC slab-column joint without punching shear reinforcement presents a rigidplastic failure form when punching shear occurs [60], and the failure area forms a punching cone (Figure 9). In this paper, the effect of the parts which directly contact the surface of the column (part Ⅰ) is not considered. Only consider the part of the cone connected with the column (part Ⅱ) and 4 prisms can be taken as the free-body. One of the prisms free-body is shown in Figure 10. Then integrate the stress on the surfaces (part Ⅲ) of the prisms to represent the punching shear strength.  In order to apply the model in a rational way during the analysis, three assumptions are proposed: (1) The study objects are experimental slab-column joints, which load through the column. The self-weight of the specimens is not considered, and thus, the vertical squeeze is not considered. (2) The stress distribution on the analysis surface is assumed to be uniformly distributed. (3) The pinning of the reinforcement is not considered.
Punching shear failure of the slab-column joints has a three-dimensional property. In order to transform into two-dimensional plane mechanics, take the lateral surfaces of the prism (part IV) as the analysis object ( Figure 11). Due to the assumptions mentioned before, the vertical squeeze is not considered, and thus, the vertical load at the punching crack can be neglected, and there is no vertical reinforcement in the free-body. Therefore, take cy σ and y ρ as 0.
tan f c σ τ θ = (15) The equations according to static equivalent principle can be established as follows: sin sin sin sin cos x sx f x sx k k ρ σ θ σ θ ρ σ θ σ θ τ θ + = − − (16) cos cos sin The local compressive stress k σ is so small that the possible effects can be ignored here and Equation (18) is obtained from Equations (16) and (17): Equation (19) can be obtained from Equations (4), (15) and (18): In order to apply the model in a rational way during the analysis, three assumptions are proposed: (1) The study objects are experimental slab-column joints, which load through the column. The self-weight of the specimens is not considered, and thus, the vertical squeeze is not considered. (2) The stress distribution on the analysis surface is assumed to be uniformly distributed. (3) The pinning of the reinforcement is not considered.
Punching shear failure of the slab-column joints has a three-dimensional property. In order to transform into two-dimensional plane mechanics, take the lateral surfaces of the prism (part IV) as the analysis object ( Figure 11). Due to the assumptions mentioned before, the vertical squeeze is not considered, and thus, the vertical load at the punching crack can be neglected, and there is no vertical reinforcement in the free-body. Therefore, take σ cy and ρ y as 0. The equations according to static equivalent principle can be established as follows: The local compressive stress σ k is so small that the possible effects can be ignored here and Equation (18) is obtained from Equations (16) and (17): Equation (19) can be obtained from Equations (4), (15) and (18): The equation of shear stress at the crack proposed by MCFT [58] is: The equation of shear strength of RC slab can be established as follows: The equation for punching shear strength of slab-column joint is established by integrating Equation (21) as follows: where L is the perimeter of the column. ω is the width of critical crack. a d is the size of aggregate. For ease of calculation, taking ω = 0.0005 0.9h 0 sin θ and a d = 20 mm [29]. In order to verify the accuracy of the bearing capacity equation, the database of this paper is used for testifying, as shown in Figure 12.  Figure 12 shows that there is some deviation in the results of the punching shear capacity equation established by MCFT. MCFT believes that the shear capacity of the concrete slab is provided by the shear stress transferred along the crack surface [58]. However, the contribution of the flexure reinforcement to resist punching shear is ignored, so the results are relatively small. For slab-column joints without punching shear  Figure 12 shows that there is some deviation in the results of the punching shear capacity equation established by MCFT. MCFT believes that the shear capacity of the concrete slab is provided by the shear stress transferred along the crack surface [58]. However, the contribution of the flexure reinforcement to resist punching shear is ignored, so the results are relatively small. For slab-column joints without punching shear reinforcement, the ratio of flexure reinforcement is the key index that affects the punching shear strength. Through experiments and numerical simulations, Li [61] found that when the reinforcement ratio was greater than 0.95%, the probability of punching failure occurred at the slab-column joint was the largest. Li [6] increased the reinforcement ratio by 0.4% in the experiments, but the punching shear capacity increased by 12%. Chen et al. [62] concluded that the increase in the reinforcement ratio significantly improved the punching shear capacity, ductility and lateral stiffness of the slab-column joints. EC2 [52] and JSCE [53] also considered the effects of the ratio of flexure reinforcement on the punching shear strength of the slab-column joints.
It is necessary to consider the effects of the reinforcement ratio. Therefore, the reinforcement ratio ρ is used to modify the equation of punching shear strength (Equation (22)). Moreover, referring to various design codes, the critical section u m is added to the equation. Assuming that the punching shear strength equation is: where α is the undetermined coefficient. Regression analysis can be used to study the effect of the single variable reinforcement ratio ρ for punching shear capacity, and the effects of other parameters can be eliminated at the same time. Firstly, the middle term Y = (ρ) k is defined, and the specimens in the database are grouped with similar parameters (deviation less than 3%) except ρ. For example, in the first group of specimens, ρ is different, but the values of parameters, such as f cu,k , λ, f y , and h 0 are similar (the deviation is less than 3%). The specimens in the database can be divided into m groups (i = 1, 2, . . . , m), and the number of specimens in each group is greater than or equal to 2. Each group has a total of n specimens (j = 1, 2, . . . , n).
The value of k is taken 1 at first. The ratio of the experimental punching shear capacity V t to the middle term Y = (ρ) k can be calculated. Then the average value of each group X i can be calculated as follows: The relative deviation e ij of each specimen and the average deviation e of all specimens can be obtained as follows: The coefficient k is adjusted until the average deviation e is the smallest. When k = 0.2, the average deviation e reaches the smallest value (0.05379).

Calculation Results
Based on the database in this paper, the punching shear capacity of the slab-column joints without punching shear reinforcement can be calculated using Equations (28)- (30). The calculation results are shown in Figure 13 and Table 4. It is clearly shown that the dispersion of Equation (30) is smaller than Equations (28) and (29). At the same time, referring to EC2 [52], the critical section perimeter u m is taken 4(c + 4h 0 ).  In order to verify the accuracy of Equation (31). Some specimens are selected which are out of the database in this paper. Table 5 shows the specimens parameters [63]. The punching shear strength is calculated by Equation (31) and the equations used in various design codes. The calculation results are shown in Figure 14. It can be seen that the results of the punching shear capacity calculated by the equation proposed in this paper have a lower degree of dispersion, and the calculated value is closer to the experimental value.   In the above analysis process, the solution of coefficient k can evaluate the contribution of the single parameter reinforcement ratio ρ to the punching shear capacity of slab-column joints. When k = 0.2 and the average deviation e reaches the smallest, the effect of the reinforcement ratio ρ to punching shear is best demonstrated. The revised equation of punching shear strength can be established as follows: In order to verify the accuracy of Equation (31). Some specimens are selected which are out of the database in this paper. Table 5 shows the specimens parameters [63]. The punching shear strength is calculated by Equation (31) and the equations used in various design codes. The calculation results are shown in Figure 14. It can be seen that the results of the punching shear capacity calculated by the equation proposed in this paper have a lower degree of dispersion, and the calculated value is closer to the experimental value.

Conclusions
The present paper introduces the models for punching shear strength of slab-column joints without punching shear reinforcement in the design codes. A database with 251 specimens is established by searching the literature. Then, the punching shear strength is predicted by the various design codes. Finally, referring to the CMFT theory, an equation for the punching shear capacity of the slab-column joints without punching shear reinforcement is established. By summarizing the investigation of design codes and the analytical model established in this paper, the following conclusions can be obtained: (1) The database established in this paper has the characteristic of a small amount of extreme data and adequate specimens. It is suitable for the various research aspects of the slab-column joints without punching shear reinforcement in further work. (2) There are some differences in the parameters included in the design codes.

Conclusions
The present paper introduces the models for punching shear strength of slab-column joints without punching shear reinforcement in the design codes. A database with 251 specimens is established by searching the literature. Then, the punching shear strength is predicted by the various design codes. Finally, referring to the CMFT theory, an equation for the punching shear capacity of the slab-column joints without punching shear reinforcement is established. By summarizing the investigation of design codes and the analytical model established in this paper, the following conclusions can be obtained: (1) The database established in this paper has the characteristic of a small amount of extreme data and adequate specimens. It is suitable for the various research aspects of the slab-column joints without punching shear reinforcement in further work. (2) There are some differences in the parameters included in the design codes. The design codes proposed in this paper all consider the relative of punching shear capacity with the concrete strength. Where only GB50010 [49] uses tensile strength in the equation, the other design codes use compressive strength. Moreover, only EC2 [52] and JSCE [53] consider the effects of reinforcement ratio for punching shear. Other parameters, such as the position of the column, the size effect and the critical perimeter are not the same. (3) Design codes have different results in predicting punching shear for the database in this paper. ACI318 [50] has a relatively small coefficient of variation and a low degree of discreteness. But the predicting values are visibly higher than the experimental re-sults which indicate conservative results. JSCE [53] has a large coefficient of variation. EC2 [52] has a large dispersion, but the average value is the smallest of other design codes. GB50010 [49] has a small dispersion and average value, but the applicable range of the equation is small. (4) The punching shear failure criterion of RC slab-column (within vertical load) without punching shear reinforcement proposed are combined the shear diagonal tension failure of RC beams with a large shear-to-span ratio, and the failure is caused by the direct stress and shear stresses in the concrete. (5) The model of MCFT can be used to predict the punching shear is based on three assumptions, which are essential to transform the three-dimensional mechanics into two-dimensional planar mechanics. The punching shear strength is decided by integrating the stress on the punching surfaces, which is proposed by the stress state of the cracked concrete units in the critical shear crack. (6) The punching shear capacity equation established based on the MCFT model has limitations before being modified, which is conservative. To ensure the accuracy of the equation, a regression analysis is also applied according to the database. It is shown by the comparison between existing design codes and the proposed model that the model is more accurate when the reinforcement ratio and critical section are considered. The characteristic of the equation established in this paper is that both have reliable theoretical support and summaries of a large number of test results, as well as high prediction accuracy and low dispersion.