Behavior of Non-Shear-Strengthened UHPC Beams under Flexural Loading: Influence of Reinforcement Depth

This study was carried out in order to study the flexural behavior of fiber-reinforced ultra-high-performance concrete (UHPC) containing hybrid microsteel straight fibers and natural fine aggregates under four-point flexural loading. The experimental results revealed that the fiber pullout mechanism had a progressive pullout (collapse) mode. A highly flexural crack developed when the fiber pulling mechanism was explicitly triggered, leading to the failure of most beams. The maximum load in beams reinforced by 1.2, 1.6, and 2.0% exceeded that in beams without longitudinal reinforcement by 56, 73, and 94%, respectively. Further, bar reinforcements at 125, 115, 95, 85, and 75 mm depths led to increases of 56, 55, 73, 96, and 94% in beam load capacity, respectively. In addition, bar reinforcement at 115, 95, 85, and 75 mm depths reduced the beams’ ductility by 40, 23, 35, and 39% compared to those with 125 mm depth. All studied UHPC beams had an uncracked phase that extended to a curvature of about 7.5 × 10−6 rad, which occurred at about 10 kNm. The use of the design of experiments was exploited in this investigation to develop a prediction model for the ultimate moment capacity of UHPC beams. This prediction model took into account the sectional and material properties of UHPC beams. To carry out this analysis, a database of 25 beams, developed by other investigators, as well as the present authors, was utilized. With a mean prediction-to-test ratio of 0.92, this prediction model had a reasonable performance capacity. In turn, this model was used to generate isoresponsive surface contours that could be used for UHPC beam design.


Background
As of today, concrete is one of the principal materials used in construction. Normal concrete, however, has some demerits (i.e., its poor ductility and tensile strength) that make it unsuited for use in advanced infrastructure [1]. The drawbacks of conventional concrete have been largely overcome through the development of steel-fiber-reinforced high-performance concrete (also known as ultra-high-performance concrete (UHPC). In the 1990s, the first class of UHPC (known as reactive powder concrete) emerged [2], and the properties of this concrete have been extensively investigated [3][4][5]. The structural performance of the ductility, bearing, and shear capacities of normal and high-strength concrete beams have been significantly enhanced by the inclusion of steel fibers [6][7][8][9][10]. The introduction of steel fibers in high-strength concrete (FRC) beams also enhanced their deformability, moment capacity, and flexural rigidity [11,12].
At the material level, abundant information is available on the special rank properties of UHPC compared to normal concrete. For instance, the development of cement-hydration indexes and their relationships with the strength characteristics of UHPC were previously studied [13]. There is, however, a gap in the literature regarding the preparation of UHPC mixes that contain ternary hybrids of microsteel fibers of various aspect ratios in the presence of a binary mixture of desert sands at ambient temperatures. These sands are found in the Arabian Peninsula, which is considered the global tank for desert sands. In the literature, a single type of desert sand available in Saudi Arabia has been proposed for use in manufacturing UHPC mixes with single types of steel fibers [14]. Further, a research team in Algeria developed a UHPC mix that incorporates dune sand to strengthen a multipurpose structure [15]. In another study [16], steel fibers, macro-polypropylene fibers, and dura fibers were combined to create high-performance concrete (HPC). Additionally, the flexural, self-healing, and sound-release properties of UHPC were reported [17]. The modified Andreasen and Andersen packing model, D-optimal design, and Artificial Neural Networks have recently been employed to optimize UHPC mixtures [18]. In the last decade, studies on the mechanical behavior of UHPC beams have gradually broadened at the structural level, and full-scale bar-reinforced UHPC beams with various fiber types and reinforcement percentages were experimentally investigated [19]. Additionally, the validity of current approaches to predict the loading capacity of UHPC beams was investigated by comparison with the experimental results of small-scale specimens [20].
In practical applications, the superiority of UHPC (e.g., its high strength-to-weight, ductility, and durability properties) has contributed to the emergence of new prospects for concrete. The use of UHPC, for example, can help reduce the weight of a structure by up to 70% [21]. The source for the superior ductility of UHPC is its discontinuous fibrous system, which controls cracks and can minimize the need for continuous bar reinforcement [1]. Moreover, the use of a low water-to-binder percentage and high content of fine-constituent materials has resulted in the higher durability to UHPC, allowing this concrete to tolerate harsh environmental conditions [22]. Due to these unique properties, UHPC has been effectively employed in infrastructure (e.g., pedestrian [23,24] and road [25,26] bridges), buildings (e.g., [27][28][29]), and the retrofitting of pre-existing structures [30,31]. Notably, encouraging findings were reported on the use of UHPC in coastal facilities [32,33], earthquake-prone areas [34], and composite piling systems [35]. More recently, the structural responses of UHPC beams with different bar areas and fiber dosages to four-point loading tests were reported [36]. In the present literature, there is no information on the behavior of hybrid-fiber-reinforced UHPC beams with various effective depths. Additionally, few practical methods are available to calculate the moment capacity of bar-reinforced UHPC beams incorporating hybrid fibers, although a few models are recommended by some engineering institutes to estimate the loading strength of UHPC beams containing steel fibrous systems [37][38][39]. To broaden the applicability of UHPC, additional research aiming at more accurate methods for calculating the moment capacity of reinforced UHPC beams is required.

Research Goal and Objectives
The current investigation aimed to investigate the structural behavior of hybrid-fiber reinforced UHPC beams containing steel bars. The basic variable used in the experimental investigations was the steel-bar depth. The corresponding bar-reinforcement percentages ranged between 1.2 and 2.0%. The developed UHPC contained about 2.5% hybrid microsteel straight fibers and natural desert aggregates. The findings based on the four-point loading test exhibited the mechanical features of the hybrid-fiber reinforced UHPC beams. These results include the modes of failure, load, deflection, ductility, strain, cracking, and moment-curvature properties. Moreover, a nonlinear analytical model to predict the ultimate moment capacity of the UHPC beam's response surface considering the beam's sectional and material properties was developed. Using the proposed model, the iso-response level plots of the normalized ultimate moments of reinforced UHPC beams with different concrete compressive strengths, reinforcements, and effective-to-overall depth ratios were plotted. These contour plots could be used to design fiber-reinforced UHPC beams. In the current research, ordinary Portland cement (PC) that satisfied the ASTM C 150 requirements was used to develop the UHPC, while silica fume (SF) and fly ash (FA, class F) were used as binders. The physicochemical properties of these types of cement are presented in Table 1. Moreover, two kinds of locally available natural desert sands identified as red sand (RS) and white sand (WS) were employed as the UHPC fine aggregates. Table 2 lists the physical characteristics of the RS and WS, and Figure 1 provides the particle size distribution of PC, SF, FA, RS, and WS. overall depth ratios were plotted. These contour plots could be used to design fiberreinforced UHPC beams.

Materials
In the current research, ordinary Portland cement (PC) that satisfied the ASTM C 150 requirements was used to develop the UHPC, while silica fume (SF) and fly ash (FA, class F) were used as binders. The physicochemical properties of these types of cement are presented in Table 1. Moreover, two kinds of locally available natural desert sands identified as red sand (RS) and white sand (WS) were employed as the UHPC fine aggregates. Table 2 lists the physical characteristics of the RS and WS, and Figure 1 provides the particle size distribution of PC, SF, FA, RS, and WS.  Additionally, microstreel smooth fibers with three different lengths (F1, F2, and F3) were employed to reinforce the UHPC. Figure 2 and Table 3 present the geometrical and physicomechanical properties of these fibers, respectively. It is worth noting that the cementitious matrix is bonded strongly to the nonsmooth fibers. Their presence will, however, be limited as an increased amount normally hinders workability due to the snatching effect with the fresh matrix, which adversely affects rheology. As a consequence, the sharpness of the protrusion on the surface has a significant effect on the rheological properties of UHPC. A feasible alternative here would be to use smooth microsteel fibers.
were employed to reinforce the UHPC. Figure 2 and Table 3 pre physicomechanical properties of these fibers, respectively. It i cementitious matrix is bonded strongly to the nonsmooth fibe however, be limited as an increased amount normally hinders snatching effect with the fresh matrix, which adversely consequence, the sharpness of the protrusion on the surface has rheological properties of UHPC. A feasible alternative here w microsteel fibers.  In this investigation, UHPC was designed by applying an op concrete's packing constituents. This process was performed b Andreassen model [Equation (1)] [40][41][42][43] along with experime fresh and hardened properties of UHPC-namely, flowability, c strength tests. In Equation (2), CPFT is the cumulative (volume) is the particle size, is the smaller particle size, is the maxim is the modulus of distribution (an exponent that generally ra depending upon the workability conditions [44]). Notably, the int cause workability issues. Therefore, self-compacting UHPC w purpose, a polycarboxylate ether-based superplasticizer (comme Glenium 51, SP) was utilized. The unit weight and water content o and 65.19%, respectively. Table 4 presents the optimized consti UHPC. Moreover, the SP's water content and aggregate absorp into consideration when calculating the UHPC's mixing water. CPFT = 100 − −   In this investigation, UHPC was designed by applying an optimization process to the concrete's packing constituents. This process was performed by applying the modified Andreassen model [Equation (1)] [40][41][42][43] along with experimental observations of the fresh and hardened properties of UHPC-namely, flowability, compressive, and flexural strength tests. In Equation (2), CPFT is the cumulative (volume) percentage finer than, d is the particle size, d o is the smaller particle size, D is the maximum particle size, and q is the modulus of distribution (an exponent that generally ranges from 0.21 to 0.37, depending upon the workability conditions [44]). Notably, the introduction of fibers could cause workability issues. Therefore, self-compacting UHPC was developed. For this purpose, a polycarboxylate ether-based superplasticizer (commercially known as Master Glenium 51, SP) was utilized. The unit weight and water content of the SP were 1080 kg/m 3 and 65.19%, respectively. Table 4 presents the optimized constituent components of the UHPC. Moreover, the SP's water content and aggregate absorption of water were taken into consideration when calculating the UHPC's mixing water.  Figure 3 illustrates the procedure for preparing the UHPC in a special concrete mixer. This concrete was developed in four phases over 10-18 min. These steps involved first mixing the dry constituents until achieving reasonable homogeneity and then adding mixing water that included SP. After obtaining the concrete mixture with practical flowability, hybrid fibers were added in the third phase. This fibrous mixture was then continuously mixed until a suitable fiber dispersion was attained. Finally, the developed UHPC was cast in 50 mm layers into the molds of the beams that were prepared with the designed continuous tensile reinforcement. Moreover, 50 mm cubic specimens and 40 × 40 × 160 mm prisms were collected to evaluate the compressive and flexural behaviors of the various UHPC batches. Notably, all concrete specimens were kept in water tanks (21 ± 2 • C temperature and 100% relative humidity) until the testing date (28 days after casting).
Appl. Sci. 2021, 11, x FOR PEER REVIEW 5 of 21 mixing the dry constituents until achieving reasonable homogeneity and then adding mixing water that included SP. After obtaining the concrete mixture with practical flowability, hybrid fibers were added in the third phase. This fibrous mixture was then continuously mixed until a suitable fiber dispersion was attained. Finally, the developed UHPC was cast in 50 mm layers into the molds of the beams that were prepared with the designed continuous tensile reinforcement. Moreover, 50 mm cubic specimens and 40 × 40 × 160 mm prisms were collected to evaluate the compressive and flexural behaviors of the various UHPC batches. Notably, all concrete specimens were kept in water tanks (21 ± 2 °C temperature and 100% relative humidity) until the testing date (28 days after casting).

Mechanical Properties of the UHPC
Following ASTM C109, a 3000 kN universal testing machine was used to test the UHPC specimens under compression with a loading rate of 0.2 MPa/s. Additionally, the flexural strength of the fiber-reinforced UHPC was determined using a 30 kN universal testing machine (INSTRON, Model 3367 at 0.2 mm/min rate) according to ASTM C1609. Figure 4 shows the flexural test setup. As shown in this figure, two Linear Variable Differential Transformers (LVDTs, Tokyo Sokki, model FDP 50A, with 300 × 10 −6 strain/mm sensitivity) were utilized to acquire the deflection at the specimen's midspan. Moreover, the behavior of the beam's reinforcing bar specimens (600 mm length) under uniaxial tension was investigated using a 3000 kN universal testing machine (at a rate of 1.12 mm/min), and the ASTM A370 guidelines were followed. The average mean yielding strength and elasticity modulus of the high-tensile-strength rebars were 520 MPa and 210 GPa, respectively. The previously mentioned material properties were evaluated as the average of three typical replicas. The concrete specimens were then taken in six different batches.

Mechanical Properties of the UHPC
Following ASTM C109, a 3000 kN universal testing machine was used to test the UHPC specimens under compression with a loading rate of 0.2 MPa/s. Additionally, the flexural strength of the fiber-reinforced UHPC was determined using a 30 kN universal testing machine (INSTRON, Model 3367 at 0.2 mm/min rate) according to ASTM C1609. Figure 4 shows the flexural test setup. As shown in this figure, two Linear Variable Differential Transformers (LVDTs, Tokyo Sokki, model FDP 50A, with 300 × 10 −6 strain/mm sensitivity) were utilized to acquire the deflection at the specimen's midspan. Moreover, the behavior of the beam's reinforcing bar specimens (600 mm length) under uniaxial tension was investigated using a 3000 kN universal testing machine (at a rate of 1.12 mm/min), and the ASTM A370 guidelines were followed. The average mean yielding strength and elasticity modulus of the high-tensile-strength rebars were 520 MPa and 210 GPa, respectively. The previously mentioned material properties were evaluated as the average of three typical replicas. The concrete specimens were then taken in six different batches.
Appl. Sci. 2021, 11, x FOR PEER REVIEW  Table 5 summarizes the mechanical properties of the UHPC batches. The pro include the compressive, flexural, and tensile strengths ( , , and , respective the elasticity modulus ( ). In this investigation, and were directly me however, was calculated using Equation (2) (following the French code in [37]) equation, ℎ (in mm) is the height of the flexural test prism (40 mm). However, to e the of fine-grained UHPC, the model proposed by Ma et al. in [45] was adop given by Equation (3)

Details of Specimens
In the current experimental program, 12 beams (two samples of six various set investigated, and two typical beams were constructed for each set. The first set w control set (i.e., without bar reinforcement); however, the remaining sets were rein by two 12 mm bars. The developed beams had the dimensions of 150 mm (width mm (height) × 600 mm (length). The main variable in this study was the effective ( ) of the tensile reinforcement bars while keeping the shear span-to-effective dep at a constant / ratio (1.8) for all beams. Figure 5 outlines the dimensio  Table 5 summarizes the mechanical properties of the UHPC batches. The properties include the compressive, flexural, and tensile strengths ( f c , f b , and f t , respectively) and the elasticity modulus (E c ). In this investigation, f c and f b were directly measured; however, f t was calculated using Equation (2) (following the French code in [37]). In this equation, h (in mm) is the height of the flexural test prism (40 mm). However, to evaluate the E c of fine-grained UHPC, the model proposed by Ma et al. in [45] was adopted, as given by Equation (3):

Details of Specimens
In the current experimental program, 12 beams (two samples of six various sets) were investigated, and two typical beams were constructed for each set. The first set was the control set (i.e., without bar reinforcement); however, the remaining sets were reinforced by two 12 mm bars. The developed beams had the dimensions of 150 mm (width) × 150 mm (height) × 600 mm (length). The main variable in this study was the effective depth (d) of the tensile reinforcement bars while keeping the shear span-to-effective depth ratio at a constant a/d ratio (1.8) for all beams. Figure 5 outlines the dimensions and reinforcements of the investigated specimens, while Table 6 lists the properties of the UHPC beams.

Details of the Testing System
The structural responses of UHPC beams (Table 6) to four-point loading were studied in the current research. The pertinent test details are presented in Figure 6, and the test was conducted using a hydraulic universal testing machine (AMSLER, 2000 kN capacity with a 0.4 mm/min loading rate). The test instruments included five different sets. First, three 60 mm concrete-strain gauges (C-1, C-2, and C-3) were attached at various locations ( Figure 6) to evaluate the concrete strains. Secondly, two 5 mm electrical-sensitive strain gauges were fixed on the B to F beams' steel rebars (before casting) to obtain their strain responses. Thirdly, two vertically mounted LVDTs were attached beneath the beam's midspan to achieve vertical displacement. Additionally, two horizontally located LVDTs were attached near the beam's centerline to acquire the crack width response. Lastly, the applied loads were measured by load-cells bonded to the machine's actuator. A data acquisition device (Tokyo Sokki, model TDS-630 with a speed of 1000 channels in 0.1 s) was employed to synchronize the data of the above-mentioned instruments. After testing each beam, high-quality photos were captured to assess the cracking patterns and modes of failure.  The structural responses of UHPC beams (Table 6) to four-point loading were studied in the current research. The pertinent test details are presented in Figure 6, and the test was conducted using a hydraulic universal testing machine (AMSLER, 2000 kN capacity with a 0.4 mm/min loading rate). The test instruments included five different sets. First, three 60 mm concrete-strain gauges (C-1, C-2, and C-3) were attached at various locations ( Figure 6) to evaluate the concrete strains. Secondly, two 5 mm electrical-sensitive strain gauges were fixed on the B to F beams' steel rebars (before casting) to obtain their strain responses. Thirdly, two vertically mounted LVDTs were attached beneath the beam's midspan to achieve vertical displacement. Additionally, two horizontally located LVDTs were attached near the beam's centerline to acquire the crack width response. Lastly, the applied loads were measured by load-cells bonded to the machine's actuator. A data acquisition device (Tokyo Sokki, model TDS-630 with a speed of 1000 channels in 0.1 s) was employed to synchronize the data of the above-mentioned instruments. After testing each beam, high-quality photos were captured to assess the cracking patterns and modes of failure. Appl. Sci. 2021, 11, x FOR PEER REVIEW 8 of 21  Figure 7 shows the cracking responses of the UHPC beams to flexural testing (0). In the elastic stage of this test (Figure 8), no cracks were detected. Microcracks were typically initiated at the beam's soffit surface and propagated vertically toward its top surface. The pullout of fibers from the concrete caused these microcracks [20]. When the load on a single fiber is greater than the fiber-matrix bonding resistance, fiber pullout ensues, which increases the load that the other fibers must overcome [1]. This mechanism made the fibers undergo a progressive pullout (collapse) mode. Figure 7 presents the enhanced deformability of the bar-reinforced UHPC (B-F) beams compared to the control beam (A); however, all beams (except for E1) failed due to one major flexural crack formation. This cracked section developed when the fiber-pullout phenomenon explicitly occurred. This mechanism typically took place when the beam approached its ultimate load, followed by steel bar yielding.

Modes of Failure and Crack Pattern
In the E1 case, two major flexural cracks extending from the beam's bottom surface to the loading points caused the failure. The best deformability was observed for E2, where tightly spread flexural cracks diffused along the beam's pure flexural zone. However, no clear-cut relationship between the reinforcement depth and cracking intensity and distribution was observed. The observed discrepancy in the failure patterns of the beam replicas could be attributed to their variable fiber distribution (dispersion). Therefore, the combined mechanism of fiber pullout and bar yielding likely precipitated the failure of the UHPC beams.  Figure 7 shows the cracking responses of the UHPC beams to flexural testing (0). In the elastic stage of this test (Figure 8), no cracks were detected. Microcracks were typically initiated at the beam's soffit surface and propagated vertically toward its top surface. The pullout of fibers from the concrete caused these microcracks [20]. When the load on a single fiber is greater than the fiber-matrix bonding resistance, fiber pullout ensues, which increases the load that the other fibers must overcome [1]. This mechanism made the fibers undergo a progressive pullout (collapse) mode. Figure 7 presents the enhanced deformability of the bar-reinforced UHPC (B-F) beams compared to the control beam (A); however, all beams (except for E1) failed due to one major flexural crack formation. This cracked section developed when the fiber-pullout phenomenon explicitly occurred. This mechanism typically took place when the beam approached its ultimate load, followed by steel bar yielding.

Modes of Failure and Crack Pattern
In the E1 case, two major flexural cracks extending from the beam's bottom surface to the loading points caused the failure. The best deformability was observed for E2, where tightly spread flexural cracks diffused along the beam's pure flexural zone. However, no clear-cut relationship between the reinforcement depth and cracking intensity and distribution was observed. The observed discrepancy in the failure patterns of the beam replicas could be attributed to their variable fiber distribution (dispersion). Therefore, the combined mechanism of fiber pullout and bar yielding likely precipitated the failure of the UHPC beams.

Load-Deflection Response
The measured load-deflection relations for the UHPC beams are presented in Figure 9, while the data related to these curves are summarized in Table 7. The basic data in this table are defined in Figure 8. These data include the elasticity limit (P cr ), yield (P y ), and ultimate (P u ) loads, as well as their corresponding midspan displacements (i.e., ∆ cr , ∆ y , and ∆ u ). Notably, P y was estimated based on the loads corresponding to the steel's yielding strains (0.0025) measured at the steel bar surfaces. Additionally, ∆ m (Table 7) was taken as the beam's maximum deflection at total failure. to the loading points caused the failure. The best deformability was observed for E2, where tightly spread flexural cracks diffused along the beam's pure flexural zone. However, no clear-cut relationship between the reinforcement depth and cracking intensity and distribution was observed. The observed discrepancy in the failure patterns of the beam replicas could be attributed to their variable fiber distribution (dispersion). Therefore, the combined mechanism of fiber pullout and bar yielding likely precipitated the failure of the UHPC beams.   As shown in Figure 8, the typical load-deflection behavior of the bar-reinforced and unreinforced UHPC beams started with a linear relation and ended with the development of flexural cracks that decreased the flexural stiffness of the beams (decreasing the slope of the load-deflection curve). The end of this second linear stage of the load-deflection response was characterized by bar-reinforcement yielding. Afterward, the beam's deflection increased nonlinearly up to the peak load.

Load-Deflection Response
As listed in Table 7, the plain UHPC beam A2 had the minimum P cr (85 kN), while the reinforced beam C2 (with 115 mm effective depth) had the maximum P cr (186 kN). The average P cr (a measure of the first cracking load for all beams induced by the UHPC matrix behavior) was 136.6 kN, with a standard deviation and coefficient of variation of 34.4 kN and 18.5%, respectively. Additionally, P u was dependent on the existence of bar reinforcements and their effective depths. As expected, the lowest average load-carrying capacity of beams was observed for the plain UHPC beams (172.8 kN). The maximum load of the UHPC beams was notably increased by the introduction of the bar reinforcement. Compared to the A-type beam, the use of bar reinforcements at 125, 115, 95, 85, and 75 mm enhanced P u by 56, 55, 73, 96, and 94%, respectively. However, the maximum P u (339.3 kN) was recorded for the bar-reinforced UHPC beam with an effective depth of 85 mm. The ductility analysis of UHPC beams is discussed in Section 3.3.   The measured load-deflection relations for the UHPC beams are presented in Figure  9, while the data related to these curves are summarized in Table 7. The basic data in this table are defined in Figure 8. These data include the elasticity limit ( ), yield ( ), and ultimate ( ) loads, as well as their corresponding midspan displacements (i.e., ∆ , ∆ , and ∆ ). Notably, was estimated based on the loads corresponding to the steel's yielding strains (0.0025) measured at the steel bar surfaces. Additionally, ∆ (Table 7) was taken as the beam's maximum deflection at total failure.    Based on the data given in Table 7, some robust correlations with a correlation coefficient (R 2 ) greater than 70% were observed ( Figure 10). As expected, linear relations with negative slopes were obtained for P y , P u , and ∆ m with d. This finding indicates that the yielding and ultimate loads linearly decreased as the effective depth of bar reinforcement increased. This finding is consistent with that reported in [46] for high-performance concrete beams with overall depths of 610-900 mm. Moreover, the trendline of µ u was found to be nonlinear for d (third-order polynomial). This result indicates that the UHPC beam's ultimate ductility could be maximized effectively at various effective depths. Therefore, the typical design code recommendation that limits the zero-sectional stain depths of concrete beams to ensure satisfactory ductility is likely invalid for UHPC beams. This result contradicts the findings of Bernardo and Lopes [47] for high-strength concrete beams.

Ductility Analysis
The ductility of a structural system is a measure of its plasticity-to-elasticity deformations and can be used to describe a system's toughness. This structural parameter is generally measured by the deflection ductility index (µ) [1,8,36,48]. The amount of the apparent mid-span deflection is, therefore, required for the evaluation of µ. In the current investigation, two deflection ductility indices (i.e., µ u and µ m ) were calculated to quantify the ductility properties of UHPC beams. Here, µ u and µ m were obtained by dividing the peak load and ultimate deflections by those corresponding to the bar-reinforcement yielding (∆ y ), as given by Equations (4) and (5), respectively: The results of the ductility indices are presented in Table 7. In this table, the B-type beam was taken as the control beam for comparing the ductility properties of the barreinforced UHPC beams, as the unreinforced beam clearly had inferior ductility. In Table 7, the percentage differences between the ductility parameters and the control ones are given in parentheses. From both calculated indices, it can be concluded that the ductility of the UHPC beam decreased as its effective depth decreased. Based on µ m , the ductility of the beam with a 75 mm effective depth (F-type) was surprisingly comparable to that of µ m with a 125 mm effective depth (B-type). Thus, µ u was found to be more practical for evaluating the ductility of UHPC beams. According to this property, the best tensile reinforcement effective depth after 125 mm was 105 mm. Based on this particular investigation, a 125 mm effective depth appears to provide the highest capability for fibers to control cracks. of was found to be nonlinear for (third-order polynomial). This result indicates that the UHPC beam's ultimate ductility could be maximized effectively at various effective depths. Therefore, the typical design code recommendation that limits the zero-sectional stain depths of concrete beams to ensure satisfactory ductility is likely invalid for UHPC beams. This result contradicts the findings of Bernardo and Lopes [47] for high-strength concrete beams.

Ductility Analysis
The ductility of a structural system is a measure of its plasticity-to-elasticity deformations and can be used to describe a system's toughness. This structural parameter is generally measured by the deflection ductility index ( ) [1,8,36,48]. The amount of the apparent mid-span deflection is, therefore, required for the evaluation of . In the current investigation, two deflection ductility indices (i.e., and ) were calculated to quantify the ductility properties of UHPC beams. Here, and were obtained by dividing the peak load and ultimate deflections by those corresponding to the barreinforcement yielding (Δ ), as given by Equations (4) and (5), respectively: The results of the ductility indices are presented in Table 7. In this table, the B-type beam was taken as the control beam for comparing the ductility properties of the barreinforced UHPC beams, as the unreinforced beam clearly had inferior ductility. In Table  7, the percentage differences between the ductility parameters and the control ones are given in parentheses. From both calculated indices, it can be concluded that the ductility of the UHPC beam decreased as its effective depth decreased. Based on , the ductility of the beam with a 75 mm effective depth (F-type) was surprisingly comparable to that of with a 125 mm effective depth (B-type). Thus, was found to be more practical for evaluating the ductility of UHPC beams. According to this property, the best tensile

Load-Strain Response
In the current research, the load-strain responses were directly measured by strain gauges attached to the concrete and steel bars, namely, C-1, C-2, C-3, and S strain gauges, as shown in Figure 6. As an example of these responses, the load-strain curves for beam B2 are presented in Figure 11. In this figure, all strain measurements were compressive (negative) up to 20 mm below the beam's top surface, while the mean strains were tensile (positive) on the steel rebars (at 125 mm below the top surface). The concrete strain at C-1 progressively increased even after the peak load; however, strains on C-2 and C-3 were likely constant before reaching the beam's end-of-elasticity moment and tended to decrease after the peak load. This reduction in concrete strains indicated the propagation of tensile cracks to the positions of the associated strain gauges. Notably, the ultimate concrete strains at failure were about 0.0038. Additionally, after the ultimate load, the strains of C-3 were close to zero, which indicated that the sectional neutral axis was touching the location of C-3.
likely constant before reaching the beam's end-of-elasticity moment and tended to decrease after the peak load. This reduction in concrete strains indicated the propagation of tensile cracks to the positions of the associated strain gauges. Notably, the ultimate concrete strains at failure were about 0.0038. Additionally, after the ultimate load, the strains of C-3 were close to zero, which indicated that the sectional neutral axis was touching the location of C-3. The load-strain responses of steel bars were similar to the beams' momentcurvatures (discussed in Section 3.6). This response was typically linear at various stages, with reduced slopes as the beam's flexural stiffness decreased after achieving the beam's elasticity limit and bar yielding. However, the tensile strains of steel were almost constant (about 0.0085) after reaching the beam's ultimate load. Figure 12 shows the load-crack width responses for the UHPC beams. This figure illustrated that the introduction of tensile reinforcement bars in the beam notably increased the beam's cracking capacity before failure, as the maximum crack widths of the bar-reinforced concrete beams (10.8-18.4 mm) were significantly less than those of the control beam (7.4-8.0 mm), as outlined in Table 7. As an established fact, the maximum crack width increased as the reinforcement effective depth decreased. After the beam reached its maximum load, as the effective depth increased, the rate of crack propagation decreased (shown in Figure 12). This preferable behavior indicated that the efficient design of bar reinforcements in UHPC beams could help control the development of cracking and thus play an important role in beam serviceability and durability. The load-strain responses of steel bars were similar to the beams' moment-curvatures (discussed in Section 3.6). This response was typically linear at various stages, with reduced slopes as the beam's flexural stiffness decreased after achieving the beam's elasticity limit and bar yielding. However, the tensile strains of steel were almost constant (about 0.0085) after reaching the beam's ultimate load. Figure 12 shows the load-crack width responses for the UHPC beams. This figure illustrated that the introduction of tensile reinforcement bars in the beam notably increased the beam's cracking capacity before failure, as the maximum crack widths of the barreinforced concrete beams (10.8-18.4 mm) were significantly less than those of the control beam (7.4-8.0 mm), as outlined in Table 7. As an established fact, the maximum crack width increased as the reinforcement effective depth decreased. After the beam reached its maximum load, as the effective depth increased, the rate of crack propagation decreased (shown in Figure 12). This preferable behavior indicated that the efficient design of bar reinforcements in UHPC beams could help control the development of cracking and thus play an important role in beam serviceability and durability.

Moment-Curvature Relations
In the current study, the moment-curvature responses of the UHPC beams were evaluated for the three characterized phases (uncracked, cracked, and nonlinear phases (Figure 13a)) using the loading and strain measurements. The readings of the strain gauges were employed to determine the sectional strain distribution over the effective depth under the assumption of plane sections remaining planar after bending. Thus, the linear distribution of these strains was applied to evaluate the beam's curvature. Gauge measurements were employed as long as their values had reasonable consistency, while the values associated with crack propagation in concrete surfaces were ignored. Therefore, a general evaluation of representative moment-curvature curves for the UHPC beams is shown in Figure 13b. These curves demonstrate that all beams had an uncracked phase that extended to a curvature of about 7.5 × 10 −6 rad, which occurred at about 10 kNm (at the first cracking moment). Beyond this curvature, the momentcurvature behavior can be assumed as linear. In this linear cracked phase, the curvature quickly increased for the unreinforced beam (A-type); however, the bar reinforcement helped control the curvature deformations of the UHPC (B-F type) beams. The maximum observed curvatures (about 100 × 10 −6 rad) were observed for the beams with bar-effective depths of 125 and 115 mm. These findings suggested that increasing the bar-effective depth of the UHPC beam would likely enhance its curvature capacity.

Moment-Curvature Relations
In the current study, the moment-curvature responses of the UHPC beams were evaluated for the three characterized phases (uncracked, cracked, and nonlinear phases (Figure 13a)) using the loading and strain measurements. The readings of the strain gauges were employed to determine the sectional strain distribution over the effective depth under the assumption of plane sections remaining planar after bending. Thus, the linear distribution of these strains was applied to evaluate the beam's curvature. Gauge measurements were employed as long as their values had reasonable consistency, while the values associated with crack propagation in concrete surfaces were ignored. Therefore, a general evaluation of representative moment-curvature curves for the UHPC beams is shown in Figure 13b. These curves demonstrate that all beams had an uncracked phase that extended to a curvature of about 7.5 × 10 −6 rad, which occurred at about 10 kNm (at the first cracking moment). Beyond this curvature, the momentcurvature behavior can be assumed as linear. In this linear cracked phase, the curvature quickly increased for the unreinforced beam (A-type); however, the bar reinforcement helped control the curvature deformations of the UHPC (B-F type) beams. The maximum observed curvatures (about 100 × 10 −6 rad) were observed for the beams with bar-effective depths of 125 and 115 mm. These findings suggested that increasing the bar-effective depth of the UHPC beam would likely enhance its curvature capacity.

Prediction of the Ultimate Moment Capacity of UHPC Beams
In this research, the ultimate moment capacity of UHPC beams was predicted using the Minitab programming package [49,50]. Here, the prediction approach relied upon the design of an experimental model that correlates the biased impacts of individual independent variables to the nonlinear observed response surfaces. This numerical modeling approach can yield continuous formulas for measured responses, taking into consideration its input limitations. Additionally, the iso-responsive contours (surfaces) can be assessed in terms of the designed investigational properties by employing this method.
In this investigation, four input parameters were used to predict the normalized ultimate moment capacities of beams (i.e., /( )), where and are the width and effective depth of the beam, respectively, and is the compressive strength of the concrete. These parameters were , , , and , which represent (MPa), (the concrete's tensile strength in MPa), (the percentage of reinforcement, %), and /ℎ (the effective-to-overall depth ratio), respectively. The normalized ultimate-moment-model

Prediction of the Ultimate Moment Capacity of UHPC Beams
In this research, the ultimate moment capacity of UHPC beams was predicted using the Minitab programming package [49,50]. Here, the prediction approach relied upon the design of an experimental model that correlates the biased impacts of individual independent variables to the nonlinear observed response surfaces. This numerical modeling approach can yield continuous formulas for measured responses, taking into consideration its input limitations. Additionally, the iso-responsive contours (surfaces) can be assessed in terms of the designed investigational properties by employing this method.
In this investigation, four input parameters were used to predict the normalized ultimate moment capacities of beams (i.e., M u / bd 2 f c ), where b and d are the width and effective depth of the beam, respectively, and f c is the compressive strength of the concrete. These parameters were x 1 , x 2 , x 3 , and x 4 , which represent f c (MPa), f t (the concrete's tensile strength in MPa), ρ (the percentage of reinforcement, %), and d/h (the effective-to-overall depth ratio), respectively. The normalized ultimate-moment-model prediction equation and coefficients are given by Equation (6) and presented in Tables 8 and 9, respectively.  Table 9. Normalized-ultimate-moment model interaction coefficients.
β 12 12 12 β 13 13 13 β 14 14 14 β 23  A wide-ranging database (Table 10) based on the present experiments and independent data was prepared for the development of this model. The collected data points encompass 25 bar-reinforced UHPC concrete beams. Notably, the previous data were carefully selected to match the current data. The primary matching criteria included concrete's comparable compressive strength and tensile-bar arrangement, a lack of transverse reinforcement, and a similar test setup under flexural loading. Accordingly, some of the test results from the four datasets [1,20,51,52] are included in Table 10. The prediction capability of the ultimate moment capacity of the bar-reinforced UHPC beam model (Equation (6)) is summarized in Table 10 and Figure 14. Table 10 shows that the suggested prediction model offered reasonable performance for the majority of the investigated data points, as the ratio of predictions to test results generally approached unity. Exceptions were observed for only four data points: R13-1, R14-1 [1], U-B3, and U-B5 [52]. This discrepancy was likely attributable to the size effect, casting methods, and testing configurations of these beams. Excluding the outlier data points (R13-1, R14-1, U-B3, and U-B5), the mean ratio was 0.92, with a standard deviation of 0.09; however, the overall average was 0.83, and the associated standard deviation was 0.23 (Table 10).  Figure 14 indicates that the majority of the predicted /( ) values using Equation (6) were close to the line of equality (LOE (prediction = test results)). However, few data points were underneath the LOE, while the majority of the data points were above this line. This finding indicates that the proposed model (Equation (6)) is likely to insignificantly overestimate the moment capacity of the UHPC beam. In conclusion, the proposed prediction model with a rational professional safety factor can be used to design UHPC beams in an adequate manner. Most of the predicted results occurred within an error band of ±10%, although a few data points were close to the ±20% model error band. This analysis demonstrated that the prediction formula (Equation (6)) can be efficiently employed to calculate the ultimate moment capacity of non-shear-strengthened and barreinforced UHPC beams under flexural loading conditions. Figure 15 presents the predicted contours of / for reinforced UHPC beams with different concrete compressive strengths ( ), reinforcement ratios ( ), and effectiveto-overall depth ( /ℎ) ratios using the proposed model (Equation (6)). In Figure 15a, the concrete's tensile strength ( ) and were assumed to be 17 MPa and 1.5%, respectively. However, and /ℎ were assumed to be 30 MPa and 0.55, respectively, as shown in Figure 15b. This figure could be used as a design aid for the UHPC beams. Notably, Figure  15a shows that at an of up to about 170 MPa, the /ℎ had little influence on the beam's moment capacity, as increasing the /ℎ at a constant did not significantly change the /( ). Figure 15b instead illustrates that the moment capacity of the UHPC beam was almost constant when was in a range of 150-190 MPa, and was in a range of up to 2%. Additionally, this figure shows that increasing the by more than 2% at a constant will likely decrease the beam's moment capacity. This result is expected, as over-reinforcement in beams will lead to brittle failure, obviating the benefits of the increased reinforcement ratios to increase the beam's loading capacity [53].  Figure 14 indicates that the majority of the predicted M u / bd 2 f c values using Equation (6) were close to the line of equality (LOE (prediction = test results)). However, few data points were underneath the LOE, while the majority of the data points were above this line. This finding indicates that the proposed model (Equation (6)) is likely to insignificantly overestimate the moment capacity of the UHPC beam. In conclusion, the proposed prediction model with a rational professional safety factor can be used to design UHPC beams in an adequate manner. Most of the predicted results occurred within an error band of ±10%, although a few data points were close to the ±20% model error band. This analysis demonstrated that the prediction formula (Equation (6)) can be efficiently employed to calculate the ultimate moment capacity of non-shear-strengthened and bar-reinforced UHPC beams under flexural loading conditions. Figure 15 presents the predicted contours of M u /bd 2 f c for reinforced UHPC beams with different concrete compressive strengths ( f c ), reinforcement ratios (ρ), and effectiveto-overall depth (d/h) ratios using the proposed model (Equation (6)). In Figure 15a, the concrete's tensile strength ( f t ) and ρ were assumed to be 17 MPa and 1.5%, respectively. However, f t and d/h were assumed to be 30 Mpa and 0.55, respectively, as shown in Figure 15b. This figure could be used as a design aid for the UHPC beams. Notably, Figure 15a shows that at an f c of up to about 170 Mpa, the d/h had little influence on the beam's moment capacity, as increasing the d/h at a constant f c did not significantly change the M u / bd 2 f c . Figure 15b instead illustrates that the moment capacity of the UHPC beam was almost constant when f c was in a range of 150-190 Mpa, and ρ was in a range of up to 2%. Additionally, this figure shows that increasing the ρ by more than 2% at a constant f c will likely decrease the beam's moment capacity. This result is expected, as over-reinforcement in beams will lead to brittle failure, obviating the benefits of the increased reinforcement ratios to increase the beam's loading capacity [53].

Conclusions
In the current research, the structural behavior of hybrid-fiber-reinforced UHPC beams containing steel bars under four-point flexural loading was studied. The basic variable in the experimental investigations was the steel bar depths. The corresponding bar reinforcement percentages ranged between 1.2 and 2.0%. The developed UHPC contained about 2.5% hybrid microsteel straight fibers and fine desert aggregates. The following conclusions were drawn based on this investigation.
The microcracks in the UHPC beams were typically initiated at the beam's soffit surface and propagated vertically toward its top surface. The fiber's pullout from the concrete caused these microcracks. When the load on a single fiber exceeded the fibermatrix bonding resistance, the pullout of the fiber ensued, thereby increasing the load, which the other fibers had to overcome. This mechanism made the fiber undergo a progressive pullout (collapse) mode.
The enhanced deformability of the bar-reinforced UHPC beams was noted compared to the control beam; however, the majority of beams failed due to one major flexural crack. This cracked section developed when fiber pullout explicitly occurred. This phenomenon was typically caused when the beam approached its ultimate load, followed by bar yielding.
The typical load-deflection behaviors of the bar-reinforced and -unreinforced UHPC beams started with linear relations and ended with the development of flexural cracks that decreased the flexural stiffness of the beams (decreasing the slope of the loaddeflection curve). The end of this second linear stage of the load-deflection response was characterized by bar-reinforcement yielding. Afterward, the beam's deflection increased nonlinearly up to the peak load.
The plain UHPC beam experienced the minimum cracking load and the lowest average load-carrying capacity of the beams. The maximum load of UHPC beams was notably increased by the introduction of the bar reinforcement.
The ductility of the UHPC beam decreased as the beam's effective depth decreased. Moreover, the deflection ductility parameter based on the ultimate defection was more practical in evaluating the ductility of UHPC beams, as it provided positive expected results for all the investigated beams.
The introduction of tensile reinforcement bars in the UHPC beam notably increased the beam's cracking capacity before failure. Additionally, the maximum crack width increased as the reinforcement effective depth decreased. After reaching the maximum load, as the effective depth increased, the rate of crack propagation decreased. This favorable behavior indicated that the efficient design of the bar reinforcement in an UHPC

Conclusions
In the current research, the structural behavior of hybrid-fiber-reinforced UHPC beams containing steel bars under four-point flexural loading was studied. The basic variable in the experimental investigations was the steel bar depths. The corresponding bar reinforcement percentages ranged between 1.2 and 2.0%. The developed UHPC contained about 2.5% hybrid microsteel straight fibers and fine desert aggregates. The following conclusions were drawn based on this investigation.
The microcracks in the UHPC beams were typically initiated at the beam's soffit surface and propagated vertically toward its top surface. The fiber's pullout from the concrete caused these microcracks. When the load on a single fiber exceeded the fiber-matrix bonding resistance, the pullout of the fiber ensued, thereby increasing the load, which the other fibers had to overcome. This mechanism made the fiber undergo a progressive pullout (collapse) mode.
The enhanced deformability of the bar-reinforced UHPC beams was noted compared to the control beam; however, the majority of beams failed due to one major flexural crack. This cracked section developed when fiber pullout explicitly occurred. This phenomenon was typically caused when the beam approached its ultimate load, followed by bar yielding.
The typical load-deflection behaviors of the bar-reinforced and -unreinforced UHPC beams started with linear relations and ended with the development of flexural cracks that decreased the flexural stiffness of the beams (decreasing the slope of the load-deflection curve). The end of this second linear stage of the load-deflection response was characterized by bar-reinforcement yielding. Afterward, the beam's deflection increased nonlinearly up to the peak load.
The plain UHPC beam experienced the minimum cracking load and the lowest average load-carrying capacity of the beams. The maximum load of UHPC beams was notably increased by the introduction of the bar reinforcement.
The ductility of the UHPC beam decreased as the beam's effective depth decreased. Moreover, the deflection ductility parameter based on the ultimate defection was more practical in evaluating the ductility of UHPC beams, as it provided positive expected results for all the investigated beams.
The introduction of tensile reinforcement bars in the UHPC beam notably increased the beam's cracking capacity before failure. Additionally, the maximum crack width increased as the reinforcement effective depth decreased. After reaching the maximum load, as the effective depth increased, the rate of crack propagation decreased. This favorable behavior indicated that the efficient design of the bar reinforcement in an UHPC beam can help control the development of cracks, which plays an important role in the serviceability and durability of beams.
All the studied UHPC beams had an uncracked phase that extended to a curvature of about 7.5 × 10 −6 rad, which occurred at about 10 kNm during the first cracking moment. Beyond this curvature, the moment-curvature behavior was shown to be linear. In this linear cracked phase, the curvature quickly increased for the unreinforced beam; however, the bar reinforcement helped control the curvature deformations of the UHPC beams. The bar-effective depth of the UHPC beam will likely enhance the beam's curvature capacity.
The design of the experimental numerical modeling approach was utilized to develop a prediction model for M u / bd 2 f c in terms of f c , f t , ρ, and d/h and prepare a wideranging database based on the present experiments and independent data. The suggested prediction model offered reasonable performance for the majority of the investigated data points, with a mean prediction-to-test ratio of 0.92. The conclusions drawn here are perhaps limited to lab-sized beams; however, normalized (size-free) design formulas have been proposed. Prospective research would investigate a size-based model of UHPC beams.