1. Introduction
As urbanization continues to advance, the demand for sand in construction is steadily increasing. While China accounts for 35% of the world’s total sand and gravel production [
1], its natural river sand reserves are diminishing year by year. Faced with the tightening supply of river sand, finding reliable alternatives has become an urgent priority.
China boasts extensive coastal regions with exceptionally rich reserves of sea-sand resources. Numerous investigations have demonstrated that seawater and sea sand, carrying corrosive ions such as chloride ions, exert an impact on concrete [
2,
3,
4,
5]. The presence of chloride ions in seawater and sea sand can expedite the setting and hardening processes of cement, accelerate hydration, and promote the formation of more hydration products [
6,
7,
8,
9]. Refs. [
10,
11,
12] have shown that sea-sand concrete (SSC) exhibits higher compressive strength at 3 d, 7 d, and 28 d compared to ordinary concrete. Furthermore, SSC demonstrates a greater drying shrinkage rate and improved impermeability. Moreover, due to a relatively low gradient in the chloride ion concentration between the pore solution and external solution, SSC displays lower chloride ion permeability than ordinary concrete does. Therefore, the use of sea sand to produce SSC can bring huge economic benefits to the construction industry.
With the modern human activity expanding to the ocean, there are more and more infrastructures constructed and serviced in marine harsh environments. To achieve an economical and sustainable development, the feasibility of using economical construction materials and the mechanical performance of structures after long-term application in harsh environments, especially the materials made from marine resources and structures serviced in marine environments, need to be extensively investigated. This growing trend is accompanied by the widespread adoption of fiber-reinforced polymer (FRP) bars, owing to their numerous advantages, such as high strength, excellent corrosion resistance, and eco-friendliness [
13]. FRP is widely used for its advantages of a low carbon, lightweight, high strength, acid and alkali corrosion resistance, and fatigue resistance [
14]. Using FRP to reinforce concrete can not only enhance the overall stiffness of the structure and improve the durability of the structure but also reduce the building environmental pollution and carbon emissions. Substituting FRP bars for conventional steel bars in SSC components offers a fundamental solution to the issue of steel corrosion.
Hua et al. [
15] conducted four-point bending tests on SSC beams reinforced with FRP bars. The findings suggest that as the reinforcement ratio increases, the failure mode of the specimens shifts from tension failure to equilibrium failure before ultimately transitioning to compression failure. Chen et al. [
16] conducted experimental and numerical research on the flexural performance of SSC beams reinforced with FRP bars exposed to a marine environment. The results demonstrate that the ultimate load-bearing capacity of the specimens increased by 28.0% in a salt-spray environment but decreased by 13.0% in a tidal environment. As the exposure time lengthened, the ultimate load-bearing capacity of the specimens exhibited an initial increase, followed by a decrease, with a simultaneous reduction in both the initial stiffness and ductility. Zhou et al. [
17] conducted experimental research on the shear performance of ultra-high-performance seawater and sea-sand concrete beams reinforced with FRP bars. The results indicate that the load–deflection curve of the beams is bilinear, displaying pronounced brittle characteristics upon failure. A higher elastic modulus of the FRP bars corresponds to a greater ultimate load-carrying capacity of the beams and reduced deflection. In summary, the failure mode of FRP-reinforced SSC beams is primarily characterized by the FRP bar fracture, representing a typical brittle failure. While a hybrid reinforcement of steel and FRP bars can transition the beam’s failure mode from brittle to ductile in nature, it is important to note that a steel bar exposed to chloride environments is susceptible to corrosion.
Steel-continuous-fiber composite bars (SFCBs) are a novel structural reinforcement material with a steel core wrapped in FRP, offering advantages such as high strength, large elastic modulus, and exceptional corrosion resistance [
18]. Additionally, an SRCB exhibits commendable secondary stiffness. Utilizing SFCBs in conjunction with FRP bars not only enhances the flexural performance and failure mode of all-FRP-reinforced SSC beams but also addresses the corrosion issue in SSC beams featuring a hybrid reinforcement of steel and FRP bars. To date, research on SFCBs combined with FRP bars in SSC beams is scarce. Existing studies have been constrained by limitations in the number of specimens tested, resulting in an incomplete exploration of the parameters influencing beam flexural behavior. Therefore, it is imperative to conduct a parametric study on SSC beams with hybrid reinforcement of SFCBs and FRP bars.
In this work, three-dimensional refined FE models were developed for SSC beams with hybrid reinforcement of SFCBs and FRP bars. The model accounts for FRP fracture behavior and the bond-slip behavior between the tensile bars and concrete. Based on the validated model, parametric studies were conducted to analyze the effects of the SSC strength, sectional steel ratio and yield strength, as well as the out-wrapped FRP elastic modulus and ultimate tensile strength of SFCBs, on the flexural performance of the beams. In addition, three failure modes were defined based on reasonable assumptions, and formulae for the bearing capacity were proposed.
4. Model Verification
4.1. Experimental Phenomena
Beams G2SF1 and G1SF2 were used as examples to conduct a comparative analysis between the FE model and experimental observations, as illustrated in
Figure 9. Strains were employed to depict the crushing and cracking phenomena in SSC beams. In the concrete, the red color represents the element strains reaching the ultimate tensile strain (
εtu = 0.00010), while black represents strains reaching the ultimate compressive strain (
εcu = −0.0035). In the reinforcement framework, black indicates element strains reaching the FRP bar fracture strain (
εfu = 0.017).
The figure reveals that when the load reaches the ultimate capacity, the tensile reinforcement elements in beam G2SF1 experience strains that reach the fracture strain, the concrete undergoes crushing, and the out-wrapped FRP reinforcement also fractures. In contrast, the increased SFCB reinforcement ratio in beam G1SF2 prevents the tensile reinforcement from fracturing, resulting in a failure mode characterized by concrete crushing alone. The strain profiles and experimental observations of these two specimens closely match. The FE model, based on the concrete plastic damage behavior and the SFCB fracture model, accurately captures the failure phenomena of both the concrete and tensile reinforcement.
Observing the strain distribution profiles of all the specimens, it is evident that the failure mode in all cases is characterized by flexural failure. Beam S3 exhibits a failure mode marked by concrete crushing, while beams G1S2, G2S1, and G3 experience failure due to FRP bar fracture. In contrast, beam G2SF1 demonstrates a failure mode featuring both concrete crushing and FRP bar fracture, while beam G1SF2 experiences failure solely due to concrete crushing. Importantly, the FE results align consistently with the experimental results. In summary, the FE model developed in this study accurately predicts the failure modes of SSC beams with hybrid reinforcement.
4.2. Load–Deflection Curve
Figure 10 shows the load–deflection curves (
P-
δ) for all the specimens. The label “Test” refers to the experimental results, while “FE” represents the FE results.
The figures clearly show that the load–deflection curve of the steel-reinforced SSC beam (S3) demonstrates a pronounced yielding phase. The load–deflection curves of the SSC beams reinforced with steel and FRP bars (G1S2 and G2S1) display a trilinear response due to the steel bar yielding. In contrast, the load–deflection curve of the FRP-reinforced SSC beam (G3) follows a bilinear pattern, reflecting the linear-elastic properties of the FRP bar. Similarly, the load–deflection curves of the SSC beams reinforced with SFCBs and FRE bars (G2SF1 and G1SF2) also exhibit a trilinear behavior, mainly attributed to the yielding of the SFCB core steel bar.
A comparison of beams S3, G1S2, G2S1, and G3 reveals that these four beams exhibit nearly identical flexural stiffness during the elastic stage (prior to cracking). However, as the beams transition into the elastic–plastic stage (after cracking), an increase in the FRP reinforcement ratio leads to a decrease in flexural stiffness, with beam G3 exhibiting the lowest value. Consequently, as the actual reinforcement ratio of the beams increases, the ultimate load capacity also increases, albeit with a reduction in the ultimate deflection.
A comparison between beams G2SF1 and G1SF2 reveals that with an increase in the SFCB reinforcement ratio, the ultimate load and deflection of beam G1SF2 have increased by 7.9% and 6.7%, respectively, in comparison to beam G2SF1. This is attributed to the increased SFCB reinforcement ratio, which results in an enhancement of the beam’s flexural stiffness. Despite the higher FRP reinforcement ratio in beam G2SF1, the lower elastic modulus of the FRP bar compared to SFCB leads to a reduced flexural stiffness compared to beam G1SF2.
The experimental results of all the specimens closely align with the FE results. Therefore, the refined FE model, which accounts for bond-slip behavior and reinforcement material fracture, can accurately predict the flexural capacity and deformations of SSC beams with hybrid reinforcement of SFCBs and FRP bars.
5. Parametric Analyses
To address the limited number of specimens in the experimental research, this section employs a validated FE model to perform a parametric analysis of SSC beams with hybrid reinforcement of SFCBs and FRP bars. In each parametric analysis, beam G1SF2 serves as the reference, with a single parameter changed at a time. The primary focus of this analysis centers on the impact of the SSC strength
fc, SFCB sectional steel ratio
as (ratio of SFCB sectional steel to the total cross-sectional area), SFCB core steel bar yield strength
fy, SFCB out-wrapped FRP elastic modulus
Esf,f, and SFCB out-wrapped FRP ultimate tensile strength
fsf,u on the flexural performance of the beams. Detailed parameter configurations are provided in
Table 3.
A comparative analysis was conducted on the analysis results, including the characteristic loads (cracking load Pcr, yield load Py, and ultimate load Pu), corresponding deflections (cracking deflection δcr, yield deflection δy and ultimate deflection δu), ultimate energy dissipation Eu, ductility coefficient μ, and failure modes. The analysis of the failure patterns reveals that flexural failure in SSC beams reinforced with hybrid reinforcement of SFCBs and FRP bars predominantly occurs in three categories: concrete crushing, FRP bar fracture, and SFCB fracture.
Due to the absence of a yield point in an all-FRP reinforcement concrete beam, the conventional definition of the ductility factor is no longer applicable to reinforced concrete beams with steel bars. In accordance with ACI 440.1R-15 [
27], the ductility factor
μ is defined as the ratio of energy dissipation at the ultimate state
Eu to the energy dissipation at normal service condition
Escr, as shown in Equation (12). Upon reaching the deflection limit, the area enclosed by the load–deflection curve and the coordinate axes represents energy dissipation under normal service conditions, denoted as
Escr. Similarly, the enclosed area represents energy dissipation at the ultimate state, denoted as
Eu. According to code GB 50608-2010 [
28], normal service conditions are defined when the mid-span deflection limit
δscr is 1/200 of the total span length
l.
Figure 11 shows the compressive and tensile stress–strain curves of SSC with an axial compressive strength
fc of 19.5, 34.2, 41.8 and 49.4 MPa, respectively.
Figure 12 depicts the uniaxial tensile stress–strain curves of an SFCB under various parametric designs.
5.1. SSC Strength
Figure 13 illustrates the influence of the SSC strength
fc on the flexural performance of the beams, while
Table 4 presents the FE model results.
Combining
Figure 13 with
Table 4, it can be observed that as the SSC strength increases, the cracking load, yield load, and ultimate load of the beams all significantly increase, while the corresponding deflection decreases. This is attributed to the heightened SSC strength, resulting in an increased total bending moment. This is due to the combination of compressive forces in the concrete’s compression zone and tensile forces in the reinforcement of the tension zone, ultimately leading to a higher characteristic load. Because FRP bars have a relatively low elastic modulus, they reach the ultimate tensile strain before the concrete undergoes crushing, causing a reduction in the deflection. This leads to a transition in the failure mode from concrete crushing to FRP bar fracture.
The energy dissipation under the normal service condition is essentially the same for all the cases, but due to a reduced ultimate deflection, the beams with higher SSC strength exhibit lower ultimate energy dissipation, resulting in a corresponding decrease in ductility. Therefore, in the case of hybrid reinforcement using SFCBs and FRP bars, it is advisable to increase the reinforcement ratio of the FRP bars appropriately to prevent premature fracture of the FRP bars, which can result in the wastage of concrete performance.
5.2. Sectional Steel Ratio of SFCBs
Figure 14 illustrates the influence of the SFCB’s sectional steel ratio on the flexural performance of the beams within the SFCB.
Table 5 presents the FE model results.
Combining
Figure 14 and
Table 5 reveals that as the SFCB’s sectional steel ratio increases, the load–deflection curve of the beam undergoes a transition from a bilinear behavior to a trilinear one. Beams with
as of 0 (all-FRP beam) and
as of 0.39 (a hybrid reinforcement of SFCBs and FRP bars) exhibit the same ultimate load capacity, with the former experiencing a 6.4% increase in the ultimate deflection compared to the latter. This is due to the lower flexural stiffness of the all-FRP beam compared to the hybrid reinforcement beam; under the same load, the deformation of the all-FRP beam is greater than that of the SFCBs and FRP bars hybrid-reinforced beam.
As the SFCB’s sectional steel ratio increases from 0.39 to 0.77 and 1.0, the elastic modulus and yield strength of the SFCB increase accordingly. However, the secondary stiffness and ultimate tensile strength decrease, while the beam’s flexural stiffness notably improves. The failure mode shifts from concrete crushing to FRP bar fracture. In the compressed concrete zone, the ultimate compressive strain has not been reached, resulting in a significant reduction in both the ultimate load and deflection. As the SFCB’s sectional steel ratio increases, the energy dissipation at the normal service condition increases with the rise in initial stiffness. However, the ultimate energy dissipation significantly decreases as the ultimate deflection reduces, resulting in a corresponding decrease in ductility.
5.3. Core Steel Bar Yield Strength of SFCBs
Figure 15 illustrates the impact of the SFCB’s core steel bar yield strength
fy on the flexural performance of beams, while
Table 6 presents the FE model results.
Combining
Figure 15 and
Table 6, it is evident that an increase in the SFCB’s core steel bar yield strength
fy leads to constant values for the beam’s cracking load and cracking deflection. Concurrently, the yield and ultimate loads increase, resulting in improved flexural stiffness. However, the corresponding deflection experiences a slight reduction. This is because an increase in the SFCB’s core steel bar yield strength does not affect the elastic modulus of the SFCB but only increases its tensile strength. The failure mode of the beam is determined by concrete crushing, and as the tensile strength of the SFCB increases, the ultimate load of the beam also increases. The energy dissipation under the normal service condition increases accordingly, while the ultimate energy dissipation remains nearly unchanged, leading to a decrease in ductility.
When the SFCB’s core steel bar yield strength increases from 300 MPa to 414 MPa, 500 MPa, and 600 MPa, the reference price per ton of reinforcement materials increases by 8.1%, 16.2%, and 24.3%, respectively. However, the ultimate load only increases by 1.5%, 2.6%, and 3.6%, respectively. Therefore, considering cost-effectiveness, increasing the yield strength of SFCB’s core steel bar is not economically viable.
5.4. Out-Wrapped FRP Elastic Modulus of SFCBs
Figure 16 depicts the influence of the SFCB’s out-wrapped FRP elastic modulus
Esf,f on the flexural performance of the beams, while
Table 7 presents the FE model results.
Combining
Figure 16 and
Table 7, it is evident that as the SFCB’s out-wrapped FRP elastic modulus
Esf,f increases, the cracking, yield, and ultimate loads of the beam all significantly increase, while the corresponding deflection decreases. This is because increasing the SFCB’s out-wrapped FRP elastic modulus results in the higher initial stiffness of the SFCB, consequently enhancing the flexural stiffness of the beam. At the same strain level, the SFCB can withstand larger concrete compressive strains, leading to increased tensile forces in the beam’s cross-section tension zone and the corresponding balancing compressive forces in the compression zone. The concrete in the compression zone reaches the ultimate compressive strain more rapidly, resulting in a higher load-bearing capacity in the beam at the same deformation and reduced ultimate deflection. As the SFCB’s out-wrapped FRP elastic modulus increases, the beam’s flexural stiffness increases during both the elastic and elastic–plastic phases. This leads to a rise in the energy dissipation under the normal service condition. However, due to a reduction in the ultimate deflection, the ultimate energy dissipation significantly decreases, resulting in reduced ductility.
Figure 16 reveals that when
Esf,f is 50, 60, and 70 GPa, the beam failure mode is concrete crushing. However, when
Esf,f is 80 GPa, the beam failure mode is the fracture of the SFCB. This is because the fracture strain of the SFCB decreases with an increase in the
Esf,f. The fracture strain (0.0125) of the SFCB with an
Esf,f of 80 GPa is the smallest and is less than the fracture strain (0.0197) of the FRP bar. Consequently, the SFCB fractures prior to the FRP bar.
In summary, increasing the SFCB’s out-wrapped FRP elastic modulus significantly enhances the beam’s load-bearing capacity and flexural stiffness, reducing mid-span deformations. However, this improvement comes at the cost of reduced ductility. Therefore, when designing the SFCB’s out-wrapped FRP elastic modulus, caution should be exercised to avoid designing it excessively high, which could lead to premature SFCB fracture.
5.5. The Ultimate Tensile Strength of SFCB’s Out-Wrapped FRP
Figure 17 illustrates the impact of the SFCB’s out-wrapped FRP ultimate tensile strength
fsf,u on the flexural performance of beams, and
Table 8 presents the FE model results.
Combining
Figure 17 with
Table 8, it is evident that an increase in the SFCB’s out-wrapped FRP ultimate tensile strength does not result in any change in the beam’s cracking and yield loads, corresponding deflections, or energy dissipation under the normal service condition. However, for beams with an
fsf,u of 750 MPa, premature fracture of the SFCB leads to lower ultimate loads, deflections, and ultimate energy dissipation compared to the other specimens. Additionally, these beams exhibit relatively lower ductility. This is because the SFCB’s out-wrapped FRP ultimate tensile strength has no influence on the initial stiffness of the SFCB but only affects the fracture strain. The smaller the SFCB’s out-wrapped FRP ultimate tensile strength, the smaller the fracture strain of the SFCB. Under the same load, a smaller fracture strain of the SFCB results in earlier fracture of the SFCB. Beams with an
fsf,u of 997 MPa, 1250 MPa, and 1500 MPa exhibit identical ultimate loads and deflections. This is because the failure of these beams is determined by concrete crushing, and excessively increasing the SFCB’s out-wrapped FRP ultimate tensile strength has no impact on the flexural performance of the beams.
6. Theoretical Analyses
6.1. Basic Assumptions
The following assumptions were considered: (1) the perfect bond between the reinforcement and SSC has been observed in previous experimental investigations conducted by Xiao et al. [
19]; (2) each cross-sectional plane perpendicular to the beam axis maintains its planarity under loading, and (3) the tensile strength of concrete is disregarded.
6.2. Failure Modes and Boundary Failure
The boundary state of the failure modes is shown in
Figure 18, where the red color of the reinforcement is expressed as the SFCBs, and the yellow color is expressed as the FRP bars. According to the constitutive models of the materials, the failure modes of SSC beams with hybrid reinforcement of SFCBs and FRP bars can be divided into three situations, where
εc and
εcu are the compressive strain and ultimate compressive strain of the SSC, respectively;
εsf,
εsf,y and
εsf,u are the tensile strain, yield strain and ultimate tensile strain of the SFCBs, respectively;
εf and
εfu are the tensile strain and ultimate tensile strain of the FRP bars, respectively;
εhu is the minimum ultimate tensile strain of the FRP bars and SFCBs,
εhu = min {
εfu,
εsf,u};
xc1 and
xc2 are the height of the compressive zone corresponding to boundary failures 1 and 2, respectively; and
ht1 and
ht2 are the height of the cross-section part in tension corresponding to boundary failures 1 and 2, respectively.
- (1)
Failure mode 1: εc = εcu, εf = εsf < εsf,y.
In failure mode 1, the SFCBs do not undergo yielding, yet the maximum compressive strain attains the ultimate strain of the SSC. This mirrors the characteristics of over-reinforced concrete beams, a condition impractical due to its propensity for brittle failure.
- (2)
Failure mode 2: εc = εcu, εsf,y < εf = εsf < εhu.
In failure mode 2, the SFCBs undergo yielding but the strain value falls short of reaching the ultimate tensile strain of both the FRP bars and SFCBs. Failure ensues when the compressive strain of the SSC reaches its ultimate limit. This outcome is anticipated in practical structures due to their inherently ductile nature.
- (3)
Failure mode 3: εc < εcu, εsf,y < εf = εsf = εhu.
In failure mode 3, the SFCBs undergo yielding, and the strain value attains the ultimate tensile strain of both the FRP bars and SFCBs. However, the compressive strain of the SSC falls short of the ultimate compressive strain. This scenario is impractical due to the lesser reinforcement.
Assuming the relative actual compressive height coefficient ξc = xc/h0, the boundary relative actual compressive height coefficient ξcb can be defined according to the plane section assumption.
Depending on the relative actual compressive height coefficient ξc, the failure modes can be categorized as follows. If ξc > ξcb1, failure mode 1 occurs; if ξcb1 ≤ ξc ≤ ξcb2, failure mode 2 occurs; and if ξc < ξcb2, failure mode 3 occurs.
6.3. Simplified Method for Calculating the Ultimate Moment
This section develops a simplified calculation for the section ultimate moment. It primarily relies on a simplified rectangular stress block model for concrete compression. Three distinct cases are proposed based on different failure modes.
6.3.1. Failure Mode 1
Figure 19 illustrates the simplified stress–strain distribution in the cross-section when the over-reinforced SSC beam fails.
According to the Chinese code GB 50010-2010 [
21], the upper steel bars are not considered, because the upper steel bars are arranged as erecting bars. Equation (15) can be obtained according to the force equilibrium of the cross-section, where
xc and
x are the actual and height of compressive concrete, respectively,
x = xcβc; and
αc and
βc are the equivalent rectangular stress factors. The standard values are
αc = 1.0,
βc = 0.8.
Combining Equations (15) and (16) allows for the calculation of the compressive concrete height
x. Substituting the compressive concrete height
x into Equation (17), the flexural capacity of the SSC beams with hybrid reinforcement of SFCBs and FRP bars can be determined.
6.3.2. Failure Mode 2
Figure 20 illustrates the simplified stress–strain distribution in the cross-section when the appropriate hybrid reinforced SSC beam fails.
Equation (18) can be obtained according to the force equilibrium of the cross-section. Combining Equations (16) and (17), the flexural capacity of the SSC beams with hybrid reinforcement of SFCBs and FRP bars can be calculated.
6.3.3. Failure Mode 3
Figure 21 illustrates the simplified stress–strain distribution in the cross-section when the under-reinforced hybrid reinforced SSC beam fails.
Equation (19) can be obtained according to the force equilibrium of the cross-section. The concrete strain
εc can be obtained by Equation (20). Combining Equations (19) and (20), the flexural capacity of the SSC beams with hybrid reinforcement of SFCBs and FRP bars can be calculated by Equation (21).
The ultimate moment of the FRP bar and SFCB hybrid-reinforced SSC beam can be determined through checking the different cases, and the flow chart for predicting the ultimate moment is shown in
Figure 22.
6.4. Verification of the Ultimate Load
Figure 23 illustrates the errors between the FE and predicted ultimate loads, where
Pu,FE and
Pu,P represent the FE and predicted results of all the specimens analyzed in
Section 5. As shown in
Figure 23, the mean and the coefficient of variation of
Pu,P/
Pu,F are 0.91 and 0.04, respectively. Thus, the predicted results agree well with the FE results, and the simplified method for calculating the ultimate moment of SSC beams with hybrid reinforcement of SFCBs and FRP bars can be provided for reference in structural design.