4.1. Ultimate Moment
Numerous mature models for calculating the flexural capacity of UHPC have been established. However, a definitive method for UHPGC has not yet been proposed. Given the strong correlation between the design theories and material properties of UHPC and UHPGC, this section will explore the applicability of existing UHPC beam flexural capacity calculation methods to UHPGC beams, aiming to inform the design and theoretical study of UHPGC components.
The French standard NF P 18-710 [
29] provides structural design methods for UHPC. NF P 18-710 necessitates iterative calculations and trial estimates to define the stress distribution within a section using UHPC cracking and ultimate tensile strengths. The NF P 18-710 method for calculating flexural capacity is illustrated in
Figure 11a. Other national standards around the world typically use simplified calculations.
In
Figure 11,
fc is the compressive strength,
fs is the strength of tensile rebar,
xt is the depth of tensile concrete,
ftd is the tensile strength, and
ftcr is the tensile cracking strength.
The method of Swiss standard MCS-EPFL for calculating the bending capacity [
30] is demonstrated in
Figure 11b. It takes into account the tensile contribution of UHPC in the tension zone, where the height of the tension zone is adjusted by a coefficient of 0.9, and simplifies the compression zone into a triangular shape.
The MCS-EPFL method can be summarized by the following formula:
The Chinese code DBJ 43/T 325-2017 [
28] simplifies the compression zone into a rectangle and does not consider the tensile strength of UHPC with stress assumptions, as depicted in
Figure 11c. In contrast, the Canadian standard CSA-S6 [
31] does not reduce the contribution of the tension zone, as shown in
Figure 11d, as same as the American design guide FHWAHIF-3-032 [
32].
Peng et al. [
20] explained that using a value of
k equal to 0.5 underestimates the flexural capacities of beams with smaller section heights (
h < 400 mm) but overestimates those of beams with
h > 600 mm. Based on flexural capacity programs and a regression of data from the literature, he proposed a reduction coefficient (
k) for the tension zone that is related to the beam height (
h) and the fiber length (
lf).
A force diagram is shown in
Figure 11e, and the calculation formula is as follows:
where
and
are the coefficients for the equivalent stress block of the compression zone, which, according to [
28], are taken as 0.92 and 0.73, respectively.
Values calculated using the above methods were compared with the experimental results, and the outcomes are presented in
Table 7.
Mt,u represents the experimental values of the five beams,
Mc,u1 represents the values calculated according to NF P18-710,
Mc,u2 represents the values calculated according to MCS-EPFL,
Mc,u3 represents the values calculated according to CSA-S6,
Mc,u4 represents the values calculated according to DBJ 43/T 325-2017, and
Mc,u5 represents the values calculated using Fei Peng’s proposed formula.
An analysis of the data presented in
Table 7 reveals that the results using DBJ43/T 325-2017 tend to be conservative, while the other calculation methods can generally predict the ultimate bending moments of the beams with reasonable accuracy, though they tend to be unsafe to varying degrees. DBJ 43/T 325-2017 does not consider the contribution of the tensile area of concrete at all, leading to significant underestimations. The NF P18-710 calculation method makes more use of the actual stress distribution and less use of simplified formulas, resulting in higher computational accuracy. However, the lack of simplification necessitates more trial calculations, which reduces its applicability. The assumption of MCS-EPFL is a triangular stress distribution in the compressed concrete area, and an estimated tensile reduction coefficient of 0.9 achieved commendable predictive results for UHPCG beams as well. Nevertheless, the Swiss standard does not consider the plastic zone of concrete during compression failure or when high reinforcement ratios are used, which leads to more pronounced underestimations for UHPCG beams with lower elastic modulus values. Fei Peng’s calculation formula achieved the smallest coefficient of variation but tended towards being unsafe.
Due to the varying degrees of overestimation caused by the above methods, this paper suggests adopting a more safety-oriented approach for the tensile reduction coefficient of UHPCG beams, selecting a value of k = 0.7. By inserting this value into Equation (4), the bending moment can be obtained as shown in
Table 8.
The method presented in this paper achieves a smaller coefficient of variation while maintaining high accuracy and conservative values, thus allowing for a relatively precise prediction of the flexural capacity of UHPGC beams.
4.2. Average Crack Spacing
The average crack spacing (
lcr) of normal concrete can be calculated using the Chinese concrete structural code GB 5001-2010 [
33].
where
lcr is the average crack spacing of the concrete,
deq is the equivalent diameter of the tensioned reinforcement,
cs is the thickness of the concrete cover,
d is the diameter of the
i-th type of tensioned reinforcement,
ni is the amount of the
i-th type of tensioned reinforcement,
h and
Ate are the effective cross-sectional height and the area around the tensioned reinforcement, respectively, and
a is the distance from the center of the tensile reinforcement to the surface of the tensile concrete.
Test results indicated that a rebar-reinforced UHPC beam could redistribute stress and generate multiple cracks prior to fiber pull-out, aligning with findings from similar research. Although the spacing between cracks decreased, the widths of the cracks were significantly reduced. Feng et al. [
16] indicated that the maximum crack widths in UHPC beams decreased by 41% in the serviceability limit states compared with those of non-fiber beams. To exactly predict the average crack spacing (
lcr) of a UHPGC beam under flexural load, Equation (13) was modified by considering the enhancement effect of steel fibers. The characteristic parameters of steel fibers were used to describe this enhancement effect, and the contribution of the steel fibers was corrected by considering the fiber influence coefficient (
ψf). Equation (19) originated from DBJ43/T 325-2017 [
28] and was validated by Feng et al. [
16].
where
ψf is the fiber influence coefficient, which can be taken as 0.05 [
16],
λf is the characteristic parameter of the steel fiber content,
lf is the fiber length, and
df is the fiber diameter.
The predicted values and experimental results of the average crack spacing for each specimen are shown in
Figure 12.
The calculation results overestimate the experimental values, with an average overestimation of 10.3%, indicating that the current crack width calculation methods are not entirely suitable for UHPGC beams. A possible reason for this is that the reduction in the elastic modulus relative to UHPC leads to a sharp increase in the number of cracks. It is worth noting that although the number of cracks has increased, the crack spacing under the serviceability limit state satisfies the requirements of DBJ43/T325-2017 very well.
4.3. Short-Term Stiffness
DBJ43/T325-2017 specifies the effective stiffness of reinforced concrete flexural members as follows:
where
I0 is the transformed section moment of inertia of the beam, which can be determined using the following formula:
where
n is the ratio of the moduli of elasticity of steel and concrete (
Es/
Ec),
h0 is the effective beam depth, and
As is the area of the tensile reinforcement in the beam.
The stiffness after cracking is calculated using the following formula:
where
βB is the coefficient of the steel fibers’ effect on the sectional short-term flexural stiffness, which can be taken as 0.2,
λf is the characteristic parameter of the steel fiber content,
ψ is the longitudinal reinforcement strain non-uniformity coefficient,
ftk is the tensile strength of concrete, and
ρte is a reinforcement ratio calculated based on the effective tensile concrete area.
σs is the longitudinal reinforcement stress. It is calculated using 0.5
Mu, 0.6
Mu, and 0.7
Mu.
η is the internal force arm coefficient, which can be taken [
34] as
.
The short-term stiffness of a flexural beam can be back-calculated using a specific formula, and the formula for calculating the deflection of a beam is as follows:
where
β is a coefficient related to the support conditions. For the experimental beams in this study, which were subjected to four-point symmetrical loading,
β is taken to be 23/216,
l0 is the calculated span of the section, and
Bs is the short-term stiffness of the section.
Transforming Equation (18) yields
Table 9 compares the predicted results for each type of flexural stiffness with the results calculated from the experiments. It was observed that the stiffness of the beams degraded significantly when the load exceeded the cracking load and continued to degrade as the load increased. The predicted values align reasonably well with the experimental values, tending to be slightly overestimated, with an average between 0.84 and 1.14. The current design specifications are suitable for calculating the short-term flexural stiffness of UHPGC beams.