Abstract

Regarding the complicated behavior of continuous deep beams, a research program including three parts was conducted. First part: three continuous concrete deep beams with different shear span-to-depth ratios (a/h) were tested. The effects of varying a/h ratio on ultimate strength and failure modes were investigated. Second part: the nonlinear finite element (FE) analyses were performed to simulate the experimental specimens and 21 large-scale continuous deep beams. The main parameters investigated were a/h ratio from 0.33 to 2 and considered 40 MPa, 60 MPa, and 80 MPa. Third part: the strut-and-tie modeling of different design codes and indeterminate strut-tie method were studied for continuous deep beams. Regardless of the a/h ratio, all beam specimens failed in shear mode with main diagonal cracks. Although EC2 load prediction was conservative for all beam models, the ACI and CSA predictions for concrete deep beams with high compressive strength were unsafe. The indeterminate truss model showed closer results to FE analysis in comparison with ACI, EC2, and CSA strut-and-tie method.

1. Introduction

Deep beams are used in shear walls, transfer girders, silos, dams, bridges, and tall buildings. The shear behavior of a simply supported reinforced concrete (RC) beam is complex because of the large number of parameters involved [1]. It becomes even more sophisticated for continuous RC deep beams due to their small shear span-to-depth ratio and statically indeterminate nature.

In deep beams, the structural response is determined by the nonlinear strain distribution even in the elastic area, in contrast to ordinary beams. Besides, since shear deformation is not negligible, ordinary beam theory cannot predict the deep beams’ load distribution. Furthermore, in the continuous concrete deep beams, maximum negative moment and maximum shear occurred at the same point. The contra flexure’s points in continuous beams might have occurred next to the critical sections of shear. Therefore, the experimental load capacity equations prepared for simply supported RC deep beams are not suitable for continuous deep beams [2].

The deep beam is specified as having a clear span (l) equal to or less than four times the overall depth (h) or a beam in which its concentrated loads are within a space of 2 h from the face of the support, by “The American Building Code” [3]. Alternatively, a deep beam is defined as having l/h ratio less than 3 by “The Euro Code 2” (EC2) [4], when “The Canadian Code” [5] defines a deep beam with l/h ratio less than 2.

To analyze deep beams, it is commonly recommended to use “strut-and-tie modeling.” This model is created on the basic characteristic of direct transfer of shear forces from the load points to the supports through compression stresses (arch action) [3, 57]. An indeterminate and determinate strut-tie modeling has been mentioned by FIB (2010) for simply supported deep beams. This model recommends an arch load transfer mechanism for beams with a/z 0.5 (a is the shear span and z is the distance between the top and bottom longitudinal reinforcement). The model suggests the truss load transfer mechanism for beams with a/z 2. It also proposed a mixture of an arch and truss load transfer mechanisms as an indeterminate strut-tie model for beams with 0.5 a/z 2 [8].

Several research types have been done on single-span deep beams; however, the deep beams used in engineering structures are mostly continuous ones [915].

Chen et al. [16] revealed that the stress distribution and crack patterns of continuous deep beams and simple deep beams are different. The continuous deep beam’s diagonal strut had a higher degree of deterioration than a simple deep beam with the same parameters and characteristics.

Oh and Shin’s study [14] revealed that deep beams with low a/d ratio (a is the shear span and d is effective depth) and high strength concrete suddenly failed without any prognosis. Regardless of concrete compressive strength, deep beams’ final failure modes depend on the a/d ratio. Singh et al. [2] showed that continuous deep beams follow the general trend of increased shear strength with a decrease in a/h ratio, similar to simply supported deep beams. Yang et al. [17] indicated that beams with a low a/h ratio had failed with main diagonal cracks, and a/h ratio has a major effect on the growth of diagonal cracks. They also revealed that the effect of section depth (h) on ultimate strength is more noticeable in continuous deep beams and high strength concrete beams than simply supported beams.

Ashour and Morley’s study [18] showed that web reinforcement of continuous deep beams is related to the a/h ratio, and the ultimate loads of continuous deep beams are influenced by the web reinforcement. Thus, it is concluded that the final load capacity of continuous deep beams is dependent on the a/h ratio. Yang and Ashour [19] reported that, with increasing concrete compressive strength, the ultimate load capacity of concrete continuous beams would be increased. It is also mentioned that, with an increasing a/h ratio, the load capacity will be decreased.

Arabzadeh’s study [20] showed that the bottom reinforcements in continuous deep beams have an important influence on failure modes and ultimate load-bearing capacity. He also revealed that the midspan deflection is the maximum deflection in deep beams regardless of boundary conditions.

Aguilar et al. [9] reported that ACI building code, either in shear design or in strut-and-tie design, has a conservative estimate of single-span beams’ shear strength. Khatab et al. [21] revealed that the ultimate capacity of continuous deep beams with low and medium strength concrete has a good agreement with strut-and-tie predictions of different design codes.

Chen et al. [22] proposed a cracking strut-and-tie model for continuous deep beams. This model considers the effect of web-shear crack on the STM of these beams. They showed that the cracking strut-and-tie model correctly predicts the effect of a/h ratio, reinforcement ratio, effective depth, and concrete compressive strength on these beams’ shear strength.

In this research, an experimental study of continuous deep beams was investigated along with the tested specimens’ numerical study. The FE analysis results and earlier experimental data were compared. Also, the numerical study of 21 full-scale continuous deep beams was conducted. The results were evaluated with strut-and-tie modeling of three different design codes and the indeterminate strut-tie method.

2. Experimental Study

2.1. Test Specimens

In this research, three continuous reinforced concrete deep beams were prepared and tested. The overall geometric dimensions and details of reinforcement for specimens are revealed in Tables 1 and 2 and Figure 1. All specimens had the same overall length, clear span, and width equal to L = 1500 mm, l = 500 mm, and b = 100 mm, respectively. The main parameter to be studied was the shear span-to-overall depth ratio (a/h). To develop three different a/h ratios, the overall depth (h) was considered h = 250 mm for BT1, h = 500 mm for BT2, and h = 750 mm for BT3. According to Librescu and Song [23], the beam can be categorized as a slender, thin beam when b/h ≤ 0.1 and h/l ≤ 0.1; otherwise, it is a thick beam. The experimental specimens were thick beams and did not need to be maintained from the out-of-plane stability.

Details of reinforcement for specimens are shown in Figure 1. Additional reinforcement cages were provided at each loading and reaction points to prevent premature local failure. The material properties of concrete and steel are presented in Tables 1 and 2.

2.2. Test Setup

A 3000 kN capacity testing frame was used to test the specimens under five-point loading, shown in Figure 2. Some steel plates were positioned between the loading points and the deep beams to avoid local failure in loading points. The lengths of steel plates were defined based on the ACI strut-and-tie model. The bearing stresses applied by the steel plates on the underlying nodes should be less than the allowable bearing stresses for these nodes. Side support plates and loading plates were considered with an equal length of 150 mm. To prevent the concrete crushing at the middle support zone, the plate had a length of 200 mm. In other words, due to the bigger dimensions of the steel plate at the middle support, the crushing of concrete does not happen at an early stage. A steel beam on top was used to divide the load into two loading points, as displayed in Figure 2. For measuring vertical deflection at midspans, two linear variable differential transformers (LVDTs) were mounted. A data logger was used to obtain the test results from LVDTs automatically.

2.3. Experimental Results
2.3.1. Cracking Behavior and Failure Modes

In beam specimen BT1, the first crack was shaped diagonal in the left span from loading plate to side support plate under 120 kN of loading. By increasing the applied load, diagonal cracks were formed at both spans between the side support plate and loading plate. Due to proper bending reinforcement at the top and bottom of the beam, no significant bending cracks were observed. Afterward, with the gradual increase in applied load, shear cracks appeared between the loading plate and the middle support. Finally, the beam failed at the load of 364 kN with a midspan deflection of 9 mm in the shear mode, as shown in Figure 3. In beam specimen BT2, the first diagonal crack was gradually formed from the side support toward the beam’s inner part under 351 kN of loading. Then, some fine cracks were shaped parallel to the initial diagonal crack. At the load of 530 kN, cracks were created around the loading and side support plates. As increasing the applied load, shear cracks appeared between the loading and the middle support plate. Finally, the beam failed at the load of 760 kN with a midspan deflection of 7.2 mm with brittle shear failure, as illustrated in Figure 3. The main failure modes in both BT1 and BT2 specimens were shear failure regardless of (a/h) ratio, same as failure modes presented in past studies for continuous concrete deep beam [17, 24]. In beam specimen BT3, the first crack was created near the side support plate at the load of 470 kN. By increasing the load to 670 kN, diagonal shear cracks were formed between the side support plates and loading plates. With increasing the applied load, local cracks were created under the loading plates. Then, the bending cracks, caused by negative moment, were formed at the top of the beam. Afterward, the BT3 specimen failed at a load of 1070 kN and midspan deflection of 6.5 mm with major diagonal cracks along with concrete crushing and cracks nearby loading and support plates. The local support failure shows that the additional reinforcement provided near the loading and bearing plates was not enough to prevent premature failure. The cracked form of specimen BT3 is illustrated in Figure 3. Therefore, the failure mode was shear cracking and support failure. As it is shown in Table 3, the initial diagonal cracks were found at 32%, 46%, and 43% of ultimate strength in specimens BT1, BT2, and BT3, respectively.

2.3.2. Load-Deflection Curves

Figure 4 displays the deflections at left midspan points of specimens with different a/h ratios. As the left and right midspan deflections were so close to each other, all of the deflections illustrated in curves are related to the left midspan. As the a/h ratio decreases, the initial stiffness and final load capacity increase, and the midspan deflection decreases. As the a/h ratio increases and the beam height decreases, the beam behavior approaches the normal beam flexural behavior with high ductility. Therefore, by decreasing the a/h ratio and increasing height, the structural function approaches the shear behavior, and the beams’ ductility decreases. As shown in Figure 4, the BT3 specimen with a/h = 0.33 is the stiffest specimen, and the BT1 specimen with a/h = 1 has the most ductility. Like Rogowsky et al. and Ashour’s research, the ductility of beams depends on the a/h ratio [24, 25].

3. Finite Element Analyses

3.1. Modeling Assumptions

In this study, several models of the continuous concrete deep beams were analyzed using ABAQUS software [26].

Three-dimensional 8-noded brick first-order fully integration continuum elements (C3D8—Bricks) were adopted to model the continuous concrete deep beams and bearing and loading plates. The longitudinal reinforcing bars and the steel stirrups were discretized with three-dimensional 2-node first-order linear truss elements (T3D2). To define the interaction of steel and concrete, all longitudinal and shear reinforcements were embedded in the concrete region.

The tie constraint is used to constrain each of the nodes on the bearing and loading plates to have the same value of displacement with the point on the concrete that it contacts. There is no relative displacement between two different surfaces, which are tied together by a tie constraint.

The mesh size is optimized to satisfy both computational time and accuracy. At the point of applying the load, to utilize an acceptable amplitude curve, the load was applied gradually as displacement offered on a steel plate. The FE model of continuous deep beams and the typical mesh is presented in Figure 5. For all beams, the same methods of modeling and mesh generation are applied.

The outline of either damage or plasticity or even both can be used to model the concrete nonlinearity [27]. Plasticity is characterized by the formation of unrecoverable deformation when loads are removed, and the lessening of elastic stiffness characterizes the damage. Since the concrete damage plasticity (CDP) model has been shown to work properly in similar studies, this model was applied for modeling concrete response [28, 29]. Important parameters to define the CDP model are shown in Table 4.

is the concrete Poisson’s ratio, Ψc (deg) is the angle of dilation, is equibiaxial to uniaxial compressive stress, and K is the ratio of the second stress invariant on the tensile meridian to that on the compressive meridian at initial yield.

The stress-strain behavior under tension and compression is also needed to describe the CDP model. The concrete stress-strain under compression was modeled based on Hognestad et al. [30]. The model is expressed in equations (1) and (2) as follows and shown in Figure 6.

For ,

For ,where is the concrete compressive stress in MPa, is concrete compressive strength in MPa, is the strain at peak compressive strength , is the compressive strain of concrete, is the ultimate concrete compressive strain, and is the elasticity modulus of concrete.

The concrete stress-strain diagram in tension is assumed to be linear up to peak tensile stress, and after cracking, it has a linear part up to zero stress at a total strain of 0.001, as shown in Figure 6. The total strain is usually assumed to be ten times the cracking strain; however, calibration of this value in finite element modeling of concrete has been recommended [26].

In Figure 6, is the tensile stress, is the concrete tensile stress, is the elasticity modulus of concrete, is the strain at peak tensile strength, and is the ultimate concrete tensile strain.

Concrete damage plasticity assumes that a variable, , characterizes the degradation of elastic stiffness. This variable is the function of two damage variables, and , which are uniaxial compressive and tensile damage variables as presented in equation (3). These two variables ( and ) are assumed to be functions of their corresponding plastic strains, and , as presented in equations (4) and (5) [31]. The damage variables can be from zero to one. Zero value is for undamaged material and one is for the total loss of strength.

and are the compressive and tensile stiffness recovery, and are the tensile and compressive stresses, and is the initial elastic modulus [26]. In these equations, and have a constant value of 0.7 and 0.1, respectively, suggested by Birtel and Mark [31].

To model the concrete with three different compressive strengths, the tensile strength and the modulus of elasticity were considered as equations (6) and (7), according to ACI318-14 [3].

The steel reinforcement behavior was considered elastoplastic. The yield and ultimate strength and elastic modulus of reinforcement were described in Table 2. Furthermore, Poisson’s ratio was assumed at 0.3.

3.2. FE Analyses of Tested Specimens and Comparison with Experimental Study

The numerical analyses of test specimens BT1, BT2, and BT3 were conducted. The resulting load-deflection diagrams are compared to experimental results, and the outcomes are presented in Figure 7. Based on numerical and experimental studies, the yield and ultimate loads are satisfactorily close to each other. As the load-displacement diagrams show, the simulated model graphs are stiffer than the experimental tests’ outcomes. This issue may come from the fact that microcracks exist in real concrete; however, FE models do not simulate them, though they reduce the stiffness of the RC models [32]. Furthermore, some ambient factors may have affected the stiffness of the experimental specimen.

The numerical failure modes of specimens BT1, BT2, and BT3 are shown in Figure 8. It should be noted that in this figure, the tensile damage of concrete is illustrated to show the concrete cracking. Comparing diagrams in Figure 8 with experimental crack patterns, shown in Figure 3, indicates that both numerical and experimental studies show the same failure modes.

3.3. Numerical Analysis

In this research, the parametric study of 21 continuous RC deep beam models was conducted. Table 5 displays the overall geometric dimensions and the reinforcement information for all the beam models used.

The parametric study was conducted by varying the shear span-to-overall depth ratio (a/h) and the studied beams’ concrete compressive strength. First, the shear span-to-overall depth ratio (a/h) was varied by changing the beams’ overall depth (h) to come up with different a/h ratios. Second, the concrete compressive strength () was considered 40 MPa, 60 MPa, and 80 MPa, to compare high strength concrete with medium and normal ones. The beam models beginning with the letter C have a 4000 mm clear span, 2000 mm shear span, and 500 mm width. These beams’ overall depth is assumed to be between 1000 and 2000 mm for various a/h ratios. The beam models beginning with the letter B have a 3000 mm clear span, 1500 mm shear span, and 300 mm width. The height of these beams is assumed to be between 1500 and 4500 mm for different a/h ratios. Each beam has been evaluated with 3 different compressive strengths of concrete. The beam models, designated only by C or B, have a compressive concrete strength of 40 MPa. Beam models labeled CM or BM have 60 MPa compressive strength, and beams named CH or BH have 80 MPa compressive strength. Details of these beam models are shown in Table 5.

3.4. Numerical Analysis Load-Deflection Response

The load-deflection diagram of continuous deep beams with different a/h ratios and is presented in Figure 9. Based on the results, a decrease in a/h ratio causes a rise in the stiffness and ultimate load capacity of the beam models similar to Singh [2] and Yang and Ashour’s [19] studies. Beams with the same a/h ratio and same overall depth show different initial stiffness as their varies from 40 MPa to 80 MPa. The deflection of midspan points shows that less deformation and ductility occur when beams have a lower a/h ratio. This can be explained by the fact that beams having large a/h ratios act similarly to the flexural critical structure. Beams with a low a/h ratio are acting similarly to the shear critical structure with an arch load transfer mechanism. As expected, the bearing capacity of beams and initial stiffness increases with adding similar to Yang and Ashour’s study [19] and confirmed in Figure 9.

3.5. Numerical Results Compared with Design Codes

In ACI 318-14 for analyzing concrete deep beams, the strut-and-tie model (STM) is suggested. Figure 10 shows a graphic STM for continuous reinforced concrete deep beams with two-point symmetrical top loads [17, 19].

The ultimate load-bearing capacity of continuous concrete deep beams (which shows the failure of struts of concrete, can be obtained from equations (8)–(10):

shows the angle between the strut and the longitudinal axis, is the shear span, is the overall section depth, is the cover of bottom reinforcement, is the cover of top reinforcement, is the concrete cylinder compressive strength, is the beam width, and . As illustrated in Figure 10, , , and are the widths of loading plates, exterior, and interior support plates.

Each design code utilizes a different effectiveness factor, for concrete. Table 6 displays the effectiveness factors indicated in EC2, CSA-A23.3-04, and ACI 318-14.

is the effectiveness factor, is the strut coefficient to account for the effect of cracking and crack-control reinforcement on the effective compressive strength of the concrete, is the angle between the reinforcement and a strut, is the total area, is bars spacing, is the beamwidth, is the concrete cylinder compressive strength, is the shear span, and is the effective beam depth.

Table 7 indicates a comparison between the numerical analysis result and the design codes STM calculations. The FE analyses show the effect of section depth on shear strength of deep beams, while the strut-and-tie model of ACI, CSA, and BC2 did not take into account this effect. The closest result for the current analysis is the ACI STM. The EC2 ultimate load capacity prediction resulted in conservative values for all continuous deep beams. The final load-bearing prediction of ACI and CSA was conservative for beams with normal and medium strength concrete (). Moreover, ACI’s ultimate load predictions underestimate the result of CH2 and BH1, and CSA underestimates the result of BH4, which are all made of high strength concrete (). With increasing from 40 MPa to 60 MPa and 80 MPa, the ACI and CSA final load prediction becomes closer to the same beams’ finite element analysis. Also, Table 7 shows that as the a/h ratio decreases, the prediction of CSA becomes closer to the FE analysis. The CSA prediction in beams C1, CM1, CH1, C2, CM2, and CH2 is too conservative, which may come from the fact that “The Canadian Code” [5] defines deep beam with l/h ratio less than 2, and these beams are not considered deep, according to this standard.

In the same shear span-to-depth ratio, with doubling from 40 MPa to 80 MPa, the failure load of C1, C2, C3, B1, B2, B3, and B4 increases between 1.67 and 1.83 times. An increase of from 40 MPa to 60 MPa, in the same beams with similar a/h ratios increase the final load of C1, C2, C3, B1, B2, B3, and B4 between 1.32 and 1.46 times.

4. Indeterminate Strut-Tie Model for Two-Span Continuous Deep Beams

In deep beams, a fraction of the load is transferred by the arch function. Considering this load transfer mechanism, ACI318-14 recommends the use of STM [3].

However, the recommended strut-and-tie method only considers the arch function. FIB (2010) has proposed an indeterminate and determinate strut-tie model for deep beams. This model recommended the arch load transfer mechanism for deep beams with a/z ≤ 0.5 and truss load transfer mechanism for beams with a/z ≥ 2. It also suggested a combination of truss and arch load transfer mechanism as an indeterminate truss for deep beams with 0.5 < a/z < 2 according to Figure 11 [8]. Like FIB’s strut-tie model, an indeterminate truss for single-span deep beams with 1 ≤ a/z ≤  has been presented by Foster and Gilbert [35]. Their model suggested the arch function for deep beams with a/z ≤ 1 and truss load mechanism transfer for deep beams with a/z . It recommended a combination of truss and arch load transfer mechanism for deep beams with 1 ≤ a/z ≤ .

Due to the nonlinear behavior in continuous deep beams, the indeterminate external and internal model, which reflects the final load and the complex nonlinear structural behavior of these beams, is presented in Figure 12 for beams with a/d ≥ 1 [34]. This indeterminate-type model represents the combination of arch function and truss mechanism [13]. Past studies showed that the role of horizontal shear reinforcement in shear strengthening of beams that do not have too many horizontal reinforcements is not significant. Therefore, in this model, the effect of horizontal shear reinforcement is not considered [24, 25].

In this model, the load distribution ratio, determined as the part of applied load transferred by a vertical tie, was proposed as equations (11) and (12) to help the designing of continuous concrete deep beams [34].

In these equations, is the load distribution ratio, is the shear span, is the effective depth, is the concrete compressive strength, is the flexural reinforcement ratio, and is the balanced flexural reinforcement ratio.

In the first step, to design a beam with the indeterminate strut-tie method, it is necessary to change the indeterminate truss into a determinate one. The conversion of an indeterminate truss to an externally determinate is done by the load distribution ratio, specified in equation (13) [34].where is the reaction distribution ratio, is the shear span, and is the effective depth.

In the second step, the load distribution ratio is used to transfer load with arch and truss functions. As the load increases, one of the arch or truss mechanisms fails. After initial failure, the indeterminate strut-tie model becomes a determinate one that can continue to transfer load by another mechanism. The second load transfer mechanism will fail by increasing the applied load, and the beam load-bearing ends. The important point is that the indeterminate strut-tie method is simultaneously reflecting the truss and arch load transfer mechanisms. The final load-bearing capacity of deep beams with shear failure is based on both mechanisms’ consecutive failure.

To show the calculation procedure and effectiveness of the indeterminate STM, beam B1 with a/h = 1 was selected and analyzed as an indeterminate strut-tie model, shown in Figure 12. The values of and for beam B1 were calculated as 0.269 and 0.375, respectively, according to equations (11)–(13).

The support reactions (, P = external load applied vertically) and the cross-sectional force of the vertical steel tie () were calculated. Maximum widths of struts and ties were computed by reaction distribution ratio, load distribution ratio, and nodal zone shapes. Required cross-sectional widths (areas) of struts and ties at the first failure were computed. The first failure of beam B1 with the indeterminate strut-tie model was happened by concrete strut as an arch mechanism failure at a load of 4905 kN (74% of numerical failure load). After initial failure, the indeterminate strut-tie model became determined and capable of transferring the load with struts and truss action ties. The remaining capacity of strut and ties after the first failure was computed. Adding load, the second failure of the model occurred as a truss action failure. When the second failure happened, the truss became unstable, which could not transfer any more loads. The final load of 6674 kN was estimated with this method. There is a 1% difference in the estimated final load of FEA and the indeterminate STM method, which is quite acceptable.

To evaluate the performance of numerical studies conducted in this research, beams C1, C2, C3, CM1, CM2, CM3, CH1, CH2, and CH3 were selected and studied with the indeterminate strut-tie model of Figure 12. The results of the indeterminate strut-tie model, according to Table 8, have a good agreement with finite element analyses. The difference in the results of FE analyses and indeterminate truss is up to 8%, which is quite acceptable. However, the difference between FEA and ACI318-14 is up to 25%. In comparing FEA with indeterminate truss, less difference percentage reflects the simultaneous operation of arch and truss function in these beams.

5. Conclusions

Based on this study, the FE analyses and experimental tests were conducted. The parameters investigated were a/h ratio and . The results were compared to the strut-and-tie and indeterminate truss model. The adjacent agreement between the finite element prediction and test results according to load-deflection curves, crack patterns, and ultimate loads indicates the model’s reliability. Below are some conclusions obtained based on the research discussed in this paper:(i)The ultimate strength of concrete beams increased when the a/h ratio decreased, and the midspan deflection decreased with decreasing a/h ratio.(ii)The place between the support and load plates encountered main diagonal cracks, which caused all tested specimens’ failure regardless of the a/h ratio.(iii)The ultimate strength of continuous deep beams was decisively controlled by concrete compressive strength. Where the concrete compressive strength of two beams with the same a/h ratio increased, the ultimate strength increased remarkably.(iv)The EC2 ultimate load predictions were conservative for all two-span concrete deep beams. The CSA and ACI final load predictions were unsafe for some of the beam models with high strength concrete. By comparing design codes, the ACI building code put forth the nearest prediction to the analysis compared to the EC2 and the CSA.(v)The prediction of CSA becomes closer to the FE analysis as the a/h ratio decreased, and with increasing a/h ratio, the FEA/CSA ratio increases(vi)Comparing the indeterminate truss model with the ACI strut-and-tie model, the indeterminate truss showed closer results with FEA in beams with a/d ≥ 1, which indicate the concurrent function of arch and truss in these beams.

Notations

L:Overall length
a:Shear span
h:Overall section depth
l:Clear span
d:Effective depth of the beam
:Vertical displacement
:Concrete cylinder compressive strength
:Concrete tensile strength
:Elasticity modulus
:Concrete initial elastic modulus
:Damage variable
:Uniaxial compressive damage variable
:Uniaxial tensile damage variable
:Compressive plastic strain
:Tensile plastic strain
:Compressive stiffness recovery
:Tensile stiffness recovery
νc:Concrete Poisson’s ratio
Ψc:Angle of dilation
:Ratio of equibiaxial to uniaxial compressive stress
K:Ratio of the second stress invariant on the tensile meridian to that on the compressive meridian at an initial yield
:Concrete compressive stress
:Strain at peak compressive strength
:Compressive strain of concrete
:Ultimate concrete compressive strain
:Tensile stress
:Compressive stress
:Concrete tensile stress
:Strain at peak tensile strength
:Ultimate concrete tensile strain
:Angle between the longitudinal axis and the strut
:Cover of bottom reinforcement
:Cover of top reinforcement
 = b:Beamwidth
:Width of loading plates
:Width of exterior support plates
:Width of interior support plates
:Effectiveness factor for concrete
:Angle between the reinforcement and a strut
:Total area of reinforcement
:Spacing of reinforcement
:Load distribution ratio
:Flexural reinforcement ratio
:Balanced flexural reinforcement ratio
:Reaction distribution ratio
:Support reaction
P:External load applied vertically
:Cross-sectional force of the vertical steel tie.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.