Next Article in Journal
Potentials of Polyacrylonitrile Substitution by Lignin for Continuous Manufactured Lignin/Polyacrylonitrile-Blend-Based Carbon Fibers
Previous Article in Journal
Remote-Controlled Activation of the Release through Drug-Loaded Magnetic Electrospun Fibers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Simulation of Engineering Cementitious Composite Beams Strengthened with Fiber-Reinforced Polymer and Steel Bars

by
Nadim I. Shbeeb
1,2,†,
Wasim S. Barham
1,* and
Wala’a Alyahya
3
1
Department of Civil Engineering, Faculty of Engineering, Jordan University of Science and Technology, P.O. Box 3030, Irbid 22110, Jordan
2
Department of Civil Engineering, Faculty of Engineering, Fahad Bin Sultan University, P.O. Box 15700, Tabuk 71454, Saudi Arabia
3
Civil Engineering Department, Jordan University of Science and Technology, P.O. Box 3030, Irbid 22110, Jordan
*
Author to whom correspondence should be addressed.
On sabbatical leave from Jordan University of Science and Technology.
Submission received: 4 April 2024 / Revised: 27 May 2024 / Accepted: 5 June 2024 / Published: 17 June 2024

Abstract

:
In this paper, the flexural performance of the Engineering Cementitious Composite (ECC)-concrete composite beam hybrid reinforced by steel and Fiber Reinforced Polymer (FRP) bars is assessed using nonlinear finite element analysis. The concrete damage plasticity model is used to model the nonlinear behavior of ECC and concrete materials. A perfect bond is assumed at the interface surface between the ECC and concrete. The validity of the numerical model is established through comparison with a previously published experimental study (overall error of about 5.4%). Consequently, the developed model is utilized to consider the effect of hybrid (FRP/steel) tensile reinforcement ratio, thickness of the ECC layer, type of FRP bars, and compressive strength of concrete on the flexure performance. It was evident from the results that the ratio of hybrid (FRP/steel) tensile reinforcement should be carefully chosen to achieve an adequate balance between ductility and carrying load capacity. Additionally, the thickness of the ECC layer plays a crucial role in controlling the hybrid reinforcement’s tensile ratio to prevent rapid failure following the yielding of steel rebars within the ECC layer. Furthermore, the type of FRP bars used in the hybrid reinforcement has influenced the flexural behavior of the composite beam. Conversely, increasing the compressive strength of the concrete has minimal impact on enhancing the mechanical characteristics of the beams, even when considering a change in the type of FRP bars.

1. Introduction

The use of Fiber-Reinforced Polymer (FRP) is on the rise in reinforced concrete (RC) structures. To reduce rust and corrosion problems, the use of FRP bars as an alternative to steel bars is a prime area of research. Compared to steel, FRP bars possess higher tensile strength, a lower weight, higher moisture resistance, higher corrosion resistance, a lower elastic modulus, and lower ductility [1,2,3,4]. The disadvantage of using FRP bars as a prime reinforcement in concrete structures is due to their low ductility, which leads to brittle failure in concrete structures. To overcome the brittle failure of concrete structures reinforced with FRP, some researchers proposed the use of both FRP and steel bars as reinforcement. As a result, enhanced ductility, durability, and corrosion resistance were observed [5,6,7,8,9]. Another alternative to enhance ductility and lower crack width is the use of engineering cementitious composite (ECC) in the tension zone around the reinforcement [10]. ECC properties such as high resistance to erosion, freezing-thawing, and sulfate environments enhance its desirability for use in the construction industry [11,12,13]. Other applications for ECC include repairing bridge decks and aging infrastructure [14,15,16] and retrofitting walls that exhibit seismic loads [17]. Ref. [18] conducted an experimental program to study the flexural performance of concrete beams reinforced by a hybrid system consisting of steel and basalt FRP bars (BFRP). The BFRP showed the same bond strength with concrete as steel bars. In terms of deflection, crack width, and flexural capacity, the hybrid steel and BFRP beams are higher than the RC beams with steel only. The authors also concluded that for structures with high requirements for flexural capacity and low requirements for deflection, hybrid RC beams can be used. Ref. [19] studied the effects of layered ECC-concrete beams on flexural performance. The authors found that the increase in both flexural strength and ductility is proportional to the increase in the ECC thickness in the tension zone. Ref. [20] experimentally and theoretically investigated the effect of using BFRP on the flexural properties of ECC beams and ECC/concrete composite beams. The experimental results showed that using an ECC matrix reinforced by BFRP instead of steel bars improved deformation capacity, strength, and energy dissipation and showed effective resistance to shear failure compared to traditional concrete beams. Refs. [21,22] studied the behavior of concrete/ECC hybrid composite beams reinforced with steel and FRP bars in terms of ductility, energy dissipation, ultimate strength, and stiffness. It was concluded that there is an improvement in the ultimate moment, yield moment, and stiffness of hybrid composite beams when compared to RC beams with the same reinforcement ratio. Also, the energy dissipation of the ECC beam is higher than that of the traditional RC beam and concrete/ECC beam. As the previous studies show, strengthening and repairing the flexural members using a composite of ECC and concrete material reinforced by a combination of steel and FRP bars is a relatively functional mechanism to enhance the flexural behavior of concrete structures. However, since the utilization of this technique is relatively new, there are limited numerical studies on this technique. Therefore, due to the high cost of the ECC and fiber materials for conducting studies, more research is needed to better understand the member’s performance and confirm its effectiveness. The use of non-linear finite element analysis is needed to simulate the behavior of composite hybrid reinforced structural members. The objectives of this research are to develop a non-linear finite element model of ECC-concrete composite beams reinforced with a hybrid FRP and steel bars system that was experimentally studied by [22] for verification purposes and extend the model to study the effect of concrete compressive strength, beam height, steel reinforcement ratio, FRP reinforcement ratio, and type of FRP bars on the composite hybrid beams ductility and load capacity.

2. Materials and Methods

Ref. [22] tested 32-hybrid reinforced composite ECC beams under four-point loading. The beam specimens were divided into eight sequences according to the reinforcement ratio of BFRP and steel bars. Each sequence has four specimens with different ECC layer thicknesses. The strain in the composite beam was measured by placing strain gauges along its height, while its deflection at various points was determined using a displacement transducer at the load point, mid-span point, and support point, as shown in Figure 1 (a more detailed schematic can be found in [22]). The results of [22] are used to verify the nonlinear finite element model developed in this study.

2.1. Material Properties

2.1.1. Concrete

In this study, the compressive behavior of concrete is modeled according to Carreira and Chu [23], and the tension behavior is modeled according to Aslani and Jowkarmeimandi [24], as those models were the most suitable for the data used in the verification section (Section 2.3). The cubic compressive ( f c ) strength is equal to 30 MPa, Poisson’s ratio (ν) is equal to 0.2, and the ultimate strain (εu) is equal to 0.003. The cylinder’s compressive strength is established by multiplying the cube’s compressive strength by 0.8 [25]. The damage in concrete is modeled using the concrete damage plasticity model (CDPM) introduced by [26]. The concrete plastic damage parameters are shown in Table 1.

2.1.2. ECC Material

Ref. [22] tested three prismatic specimens of size 40×40×160 mm to obtain the compressive stress–strain curves of ECC. The compressive stress of ECC recorded was 31.4 MPa, and its corresponding strain was 0.0036. To use the CDPM in [27], the compressive strength must be obtained through testing cylinder specimens; thus, [28] is used to obtain the cylinder compressive strength (27 MPa). As for the tension properties, the tensile strength at first cracking was recorded at 2 MPa with a corresponding strain of 0.00023, and the ultimate tensile strength was recorded at 2.4 MPa with a corresponding strain of 0.025. The tensile stress–strain behavior is modeled as an elastic-strain hardening material.

2.1.3. Steel Reinforcement and Steel Plates

The steel properties used in [22], such as the yield stress (fy) and ultimate stress (fu), are summarized in Table 2. The Poisson ratio is equal to 0.3. The reinforcement steel is modeled as an elastic-strain hardening material (bi-linear). The steel plates (used to distribute the loads in the experiment) are defined as an elastic material with a modulus of elasticity of 200 GPa and a Poisson’s ratio of 0.3.

2.1.4. FRP Reinforcement

In this study, four FRP types are used: Aramid FRP (AFRP), Glass FRP (GFRP), Carbon FRP (CFRP), and Basalt FRP (BFRP). BFRP is used for verification of the proposed model. Per [22], BFRP is a linear elastic material with an ultimate tensile strength of 1250 MPa and a modulus of elasticity of 50 GPa. The stress–strain curve of all types of FRP bars is shown in Figure 2. Table 3 shows the properties of the FRP reinforcement.

2.2. Nonlinear Finite Element Model (NLFEM)

In this paper, the finite element software ABAQUS/Standard ver. 6.9 [27] is used. The convergence and mesh sensitivity analysis carried out in the verification step (Section 2.3) yielded a medium-mesh size of 25 mm. The experimental study [22] showed no delamination or separation between the ECC layer and concrete; thus, the interface surface between the ECC layer and concrete is assumed to be perfectly bonded. The surface between the support/load steel plates and the beam is tied together to eliminate any relative motion between them, this is achieved by using the tie constraint described in [27]. A perfect bond is assumed between the reinforcement (steel bars, stirrups, FRP bars) and either concrete or ECC. This was characterized using the embedded region constraint in [27]. Both the ECC and concrete are assumed to be the host region and the reinforcement, regardless of its type, as the embedded region. The displacement control method is utilized to apply the load to the simply supported beam. The simply supported boundary conditions are applied along the support’s plates. Details of the finite element model are shown in Figure 3.

2.3. Model Verification

The model is verified by comparing the following load-deflection curves: ultimate deflection, ultimate load, crack pattern, and failure mode, between the experiment and NLFEM. Three groups of beams are modeled; each group has four tested beams from the published work of [22]. The first group is composed of steel-reinforced beams, the second is composed of hybrid-reinforced beams, and the third is composed of FRP-reinforced beams. All groups have the same geometry (150 mm × 200 mm × 1500 mm) with the same top steel reinforcement (two steel bars- size 10 mm diameter) and the same shear reinforcement (8 mm stirrup at 100 mm spacing). Table 4 shows the reinforcement detail and thickness of the ECC layer. Note that zero ECC thickness indicates a fully concrete beam, while 200 mm ECC thickness indicates a fully ECC beam.
Figure 4, Figure 5 and Figure 6 show the load-deflection curves between the experiment (EXP) and the finite element analysis (FEA). Based on the loading deflection curves, suitable agreement is observed between the experimental data and the NLFEM response. Both Group I (steel reinforcement only) and II (hybrid reinforcement) exhibit ductile failure, while Group III (BFRP reinforcement only) exhibits brittle failure, which is in agreement with the literature.
The ultimate load and corresponding deflection from the experimental and NLFEM are listed in Table 5. The error percentage between the experimental and NLFEM responses is within acceptable levels, overall error is 5.25% and 3.96% for the ultimate deflection and the ultimate load, respectively, which is less than the usual acceptable error (10%).
Figure 7, Figure 8, Figure 9 and Figure 10 show the comparison between the experimental failure modes in [22] and the NLFEM for group II (composite hybrid reinforced beams). It is shown that the dominant failure is ductile for all beams. For all beams, the steel reinforcement yielded, while the axial stress in BFRP reinforcement ranged from 46% to 66% of its ultimate tensile strength (the axial stress for reinforcement is only shown in Figure 7b). It is observed that increasing the ECC layer decreased the formation of diagonal shear cracks, conforming to the findings in [31,32,33,34]. In all the beams, the cracks extended to the compression zone. Based on the validation study, the NLFEM developed was used to conduct the parametric study in the next section.

3. Results and Discussion

3.1. Parametric Study Parameters

The parameters used in this study to investigate the mechanical characteristics of ECC beams hybrid reinforced with FRP and steel bars are longitudinal tensile steel reinforcement ratio (ρs), longitudinal tensile FRP reinforcement ratio (ρfrp), ECC layer thickness (tECC), FRP bar type, beam height, and compressive strength of concrete, as listed in Table 6. The mechanical characteristics evaluated in this study are the ultimate load (Pu), the ultimate deflection (ΔU), the ductility index (μ, defined as the ratio between ultimate deflection and yield deflection), the energy absorption (area under the load-deflection curve), and the stiffness (slope of the second linear portion). The deflection load response for all developed models is composed of three stages: the first is from the initial to the point where the first crack appeared, known as the pre-cracking or un-cracked stage. The second one is from the cracking point to the point at which the steel rebars are yielding, known as the yield point. Finally, the third stage is from the yield point to the last point, representing the failure point.

3.2. Effect of FRP Type and Concrete Compressive Strength

In studying the effect of the FRP type and compressive strength of concrete the following parameters are kept constant: beam height (200 mm), tECC (100 mm), ρs (0.0043), and ρFRP (0.0038). As shown in Figure 11, regardless of the concrete compressive strength and FRP type, the initial slope (pre-cracking stage) is approximately the same. As expected, the highest load-carrying capacity is for the CFRP reinforcement and 40 MPa compressive strength of concrete, followed by AFRP, BFRP, and GFRP. Due to the higher tensile strength of CFRP reinforcement when compared to other FRP reinforcements (Figure 2).
In Figure 12, the ultimate load, ductility index, energy absorption, and stiffness were normalized to the GFRP values shown for the different concrete compressive strengths (24 MPa, 30 MPa, and 40 MPa). It is evident that the same trend holds irrespective of concrete compressive strength, and the influence of the FRP type is more pronounced at lower concrete compressive strengths than expected. Regarding the ductility index, there is no discernable difference between GFRP and BFRP. Despite BFRP being stiffer and possessing a higher tensile strength than GFRP, GFRP fails at a higher strain value than BFRP, as shown in Figure 2.

3.3. Effect of Hybrid Tensile Reinforcement Ratio and tECC

In this section, the behavior of the composite hybrid reinforced beam is numerically modeled to discuss the effectiveness of the thickness of the ECC layer against various steel and BFRP reinforcement ratios on the ultimate load, ductility index, energy absorption, and stiffness. The steel reinforcement ratios used are 0.19%, 0.43%, 0.77%, and 1.22%, while the BFRP reinforcement ratios used are 0.38%, 0.86%, and 1.53% (Table 4 shows actual bar sizes). Both the concrete compressive strength and beam height are kept constant at 24 MPa and 200 mm, respectively. The ECC thicknesses used are 50 mm and 100 mm, placed in the tension zone.
As shown in Figure 13, using a higher cross-sectional area of steel reinforcement reduced the ultimate load to 95.3% and 91.1% by using 0.77% and 1.2% steel reinforcement ratios and a 0.38% BFRP ratio, respectively. However, as shown in Figure 14, a significant decrease in the ductility factor of the composite hybrid beams as the steel reinforcement ratio increases, reaching less than 3, which does not meet the service requirement, is observed.
Moreover, the early failure of the beams reinforced by steel reinforcement ratios of 0.77% and 1.2% led to a high decrease in the energy absorption of the beams, as shown in Figure 15. However, the beam reinforced with a steel reinforcement ratio of 0.19% and a BFRP ratio of 0.38% recorded a drop in the ultimate load and energy absorption and an improvement in ductility of about 89.9%, 88.7%, and 137.73%, respectively, compared with the beam reinforced with a 0.43% steel ratio.
In the case of the 50 mm ECC layer, it can be observed from Figure 16 that increasing the steel ratio is advantageous for improving the stiffness of the beams. However, an early failure occurred rapidly after the yielding of steel rebars at lower ultimate displacement, leading to insufficient ductility (less than 3) in the composite hybrid beams. This behavior, explained by the high stress concentration from the increasing steel ratio in the ECC layer, led to crushing in the ECC layer between the support and the loading point.
Generally, the 0.43% steel reinforcement ratio obtains a higher load-carrying capacity and sufficient ductility for the composite hybrid beam when the ECC layer is 50 mm in the tension zone. For the case of using a 100 mm ECC layer, providing a higher tensile steel reinforcement ratio significantly enhanced the ultimate load capacity, even with the BFRP ratio used, due to the increasing tensile resistance. For instance, when the BFRP reinforcement ratio was 0.38% and the composite hybrid beams were reinforced by 1.2% and 0.77% steel reinforcement ratio, the ultimate load increased by 122.4% and 108.6%, respectively, compared with the composite beam reinforced by 0.43% steel reinforcement ratio. In contrast to the reinforced composite hybrid beam with a 0.19% steel ratio, which caused a reduction in the ultimate load by about 83.9%.
The increase in the tensile steel reinforcement ratio results in increases in the yielding load and its corresponding deflection. In contrast, a slight decrease was shown in the ultimate deflection. This is because the compressive failure mode of all the beams means a fixed value of the concrete maximum compression strain, so a lower strain from the steel reinforcement can show enough tensile resistance to balance the concrete compression strain, and then the ultimate deflection is reduced. This behavior led to a decreasing ductility factor of the beam with an increasing reinforcement ratio. The beam reinforced by 1.2% and 0.77% steel ratio and 0.38% BFRP reduced the ductility factor by 72% and 78.7%, respectively, while the ductility improved by 126.6% by reducing the reinforcement ratio to 0.19%. Compared with the reference beam, the consumption of more steel reinforcement ratio improved the energy absorption and the flexural stiffness by about 118.63% and 131.2% at a steel reinforcement ratio of 1.2%, 112.5%, and 130%, respectively, using a steel reinforcement ratio of 0.77%. In comparison, the beam reinforced with a 0.19% steel reinforcement ratio showed drops in toughness and flexural stiffness, reaching up to 78.67% and 73.8%, respectively.
Generally, with the BFRP reinforcement ratio held constant, a correlation between the thickness of the ECC layer and the effective steel reinforcement ratio can be observed, leading to an appropriate load-carrying capacity and sufficient ductility of the composite hybrid beams. In composite beams with a 50 mm ECC layer in the tension zone, the 0.43% steel reinforcement ratio performs better in terms of load-carrying capacity and ductility. However, using 0.77% and 1.22% steel reinforcement ratios exhibits higher load-carrying capacity, while a 0.19% steel reinforcement ratio provides higher ductility performance in the case of composite beams with a 100 mm thickness of ECC material. Moreover, it can be noted from Figure 13 that in both cases of using 50 mm and 100 mm ECC layers with a constant steel reinforcement ratio, there was a correlation between the BFRP reinforcement ratio and the ultimate load of the beam. As evident, the peak load capacity, beam toughness, and stiffness for each specimen tend to be enhanced by utilizing a higher ratio of BFRP. On the other hand, the brittle, elastic nature of the BFRP bar material led to a reduction in the ductility of the beam. Therefore, a sufficient fiber reinforcement ratio must be chosen to meet the suitable ductility of the composite hybrid beam.

4. Conclusions

In this study, the flexural behavior of ECC-concrete beams reinforced by a hybrid combination of steel and FRP bars is numerically investigated based on a non-linear finite element method. Twelve models of beams reinforced with different configurations of steel and FRP bars and various matrix types were simulated, based on a previous experimental study, to ensure the accuracy of the numerical model in predicting the beams’ flexural behavior. Subsequently, thirty-seven models were developed to investigate the effectiveness of steel reinforcement ratio, BFRP reinforcement ratio, FRP bar type, concrete compressive strength, and beam depth on the composite hybrid beam’s flexural behavior. The general findings of this research are highlighted below:
  • The numerical simulation models accurately predict the load-deflection response with an error of less than 10% of the ultimate load, ultimate deflection, and stiffness. This provides researchers with a valuable tool for studying the behavior of ECC-concrete hybrid reinforced beams, particularly considering the higher costs associated with ECC material and FRP bars.
  • In hybrid composite beams (using 50 mm ECC), when the BFRP reinforcement ratio is kept constant, utilizing a higher steel reinforcement ratio results in a reduction in ultimate load capacity, toughness, and ultimate deflection, along with insufficient ductility (less than 3), leading to early failure immediately after yielding of steel rebars due to high stress concentration in the ECC layer.
  • In hybrid composite beams (using 100 mm ECC), when the BFRP reinforcement ratio is kept constant, using a higher steel reinforcement ratio leads to improvements in ultimate load, toughness, and flexural stiffness, but reduces beam ductility and shows a slight decrease in maximum deflection.
  • With the remaining steel reinforcement ratio constant in the hybrid reinforcement system, an increase in the ratio of BFRP enhances ultimate load capacity, toughness, and flexural stiffness, but causes a decline in maximum deflection and ductility.
  • The BFRP/steel tensile reinforcement ratio should be carefully designed to achieve the required strength, acceptable ductile performance, and improved load-carrying capacity of the composite beam.
  • Strengthening composite hybrid beams with CFRP bars instead of BFRP results in better mechanical behavior in terms of ultimate load, toughness, flexural stiffness, and plastic stiffness, albeit with a reduction in ductility index. Conversely, replacing BFRP with GFRP bars exhibits a higher enhancement in ductility performance but a reduction in other mechanical properties. Additionally, the use of AFRP bars provides moderate mechanical behavior, performing better than BFRP but still inferior to CFRP bars.
  • Using a higher grade of concrete strength in composite hybrid beams contributes minimally to improving mechanical behavior, with no noticeable effect on enhancement or reduction when using different types of FRP bars.
  • In composite hybrid reinforced beams, using a 100 mm ECC layer in the tension zone significantly impacts the load-deflection response and greatly improves mechanical performance with increasing effective beam depth.

Author Contributions

Conceptualization, N.I.S. and W.S.B.; methodology, N.I.S. and W.S.B.; validation, W.A.; formal analysis, W.A; investigation, W.A.; resources, N.I.S. and W.S.B.; data curation, W.A.; writing—original draft preparation, N.I.S., W.S.B. and W.A.; writing—review and editing, N.I.S., W.S.B. and W.A.; visualization, W.A.; supervision, N.I.S. and W.S.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors on request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cosenza, E.; Manfredi, G.; Realfonzo, R. Behavior and modeling of bond of FRP rebars to concrete. J. Compos. Constr. 1997, 1, 40–51. [Google Scholar] [CrossRef]
  2. Gravina, R.J.; Smith, S.T. Flexural behaviour of indeterminate concrete beams reinforced with FRP bars. Eng. Struct. 2008, 30, 2370–2380. [Google Scholar] [CrossRef]
  3. Soric, Z.; Kisicek, T.; Galic, J. Deflections of concrete beams reinforced with FRP bars. Mater. Struct. 2010, 43, 73–90. [Google Scholar] [CrossRef]
  4. Qin, W.; Fernand, M.; Yang, Z.-J.; Georges, I.; Zakari, M.D. Effect of fiber volume ratio on flexural behavior of RC beam with novel HFRP rebar and steel rebar. Structures 2022, 43, 1882–1900. [Google Scholar] [CrossRef]
  5. Harris, H.G.; Somboonsong, W.; Ko, F.K. New ductile hybrid FRP reinforcing bar for concrete structures. J. Compos. Constr. 1998, 2, 28–37. [Google Scholar] [CrossRef]
  6. Aiello, M.A.; Ombres, L. Structural performances of concrete beams with hybrid (fiber-reinforced polymer-steel) reinforcements. J. Compos. Constr. 2002, 6, 133–140. [Google Scholar] [CrossRef]
  7. Leung, H.Y.; Balendran, R.V. Flexural behaviour of concrete beams internally reinforced with GFRP rods and steel rebars. Struct. Surv. 2003, 21, 146–157. [Google Scholar] [CrossRef]
  8. Lau, D.; Pam, H.J. Experimental study of hybrid FRP reinforced concrete beams. Eng. Struct. 2010, 32, 3857–3865. [Google Scholar] [CrossRef]
  9. Kara, I.F.; Ashour, A.F.; Köroğlu, M.A. Flexural behavior of hybrid FRP/steel reinforced concrete beams. Compos. Struct. 2015, 129, 111–121. [Google Scholar] [CrossRef]
  10. Maalej, M.; Li, V.C. Introduction of strain-hardening engineered cementitious composites in design of reinforced concrete flexural members for improved durability. Struct. J. 1995, 92, 167–176. [Google Scholar]
  11. Şahmaran, M.; Özbay, E.; Yücel, H.E.; Lachemi, M.; Li, V.C. Frost resistance and microstructure of Engineered Cementitious Composites: Influence of fly ash and micro poly-vinyl-alcohol fiber. Cem. Concr. Compos. 2012, 34, 156–165. [Google Scholar] [CrossRef]
  12. Özbay, E.; Karahan, O.; Lachemi, M.; Hossain, K.M.; Atis, C.D. Dual effectiveness of freezing–thawing and sulfate attack on high-volume slag-incorporated ECC. Compos. Part B Eng. 2013, 45, 1384–1390. [Google Scholar] [CrossRef]
  13. Shang, X.Y.; Yu, J.T.; Li, L.Z.; Lu, Z.D. Strengthening of RC structures by using engineered cementitious composites: A review. Sustainability 2019, 11, 3384. [Google Scholar] [CrossRef]
  14. Li, V.C. Advances in ECC research. ACI Spec. Publ. 2002, 206, 373–400. [Google Scholar]
  15. Li, V.C. On engineered cementitious composites (ECC) a review of the material and its applications. J. Adv. Concr. Technol. 2003, 1, 215–230. [Google Scholar] [CrossRef]
  16. Lim, Y.M.; Li, V.C. Durable repair of aged infrastructures using trapping mechanism of engineered cementitious composites. Cem. Concr. Compos. 1997, 19, 373–385. [Google Scholar] [CrossRef]
  17. Kesner, K.; Billington, S. Investigation of ductile cement-based composites for seismic strengthening and retrofit. In Fracture Mechanics of Concrete Structures; Swets & Zeitlinger: Lisse, The Netherlands, 2001; pp. 65–72. [Google Scholar]
  18. Ge, W.; Zhang, J.; Cao, D.; Tu, Y. Flexural behaviors of hybrid concrete beams reinforced with BFRP bars and steel bars. Constr. Build. Mater. 2015, 87, 28–37. [Google Scholar] [CrossRef]
  19. Zhang, J.; Leung, C.K.; Cheung, Y.N. Flexural performance of layered ECC-concrete composite beam. Compos. Sci. Technol. 2006, 66, 1501–1512. [Google Scholar] [CrossRef]
  20. Yuan, F.; Pan, J.; Leung, C.K.Y. Flexural behaviors of ECC and concrete/ECC composite beams reinforced with basalt fiber-reinforced polymer. J. Compos. Constr. 2013, 17, 591–602. [Google Scholar] [CrossRef]
  21. Ge, W.; Ashour, A.F.; Cao, D.; Lu, W.; Gao, P.; Yu, J.; Ji, X.; Cai, C. Experimental study on flexural behavior of ECC-concrete composite beams reinforced with FRP bars. Compos. Struct. 2019, 208, 454–465. [Google Scholar] [CrossRef]
  22. Ge, W.J.; Ashour, A.F.; Yu, J.; Gao, P.; Cao, D.F.; Cai, C.; Ji, X. Flexural behavior of ECC–concrete hybrid composite beams reinforced with FRP and steel bars. J. Compos. Constr. 2019, 23, 04018069. [Google Scholar] [CrossRef]
  23. Carreira, D.J.; Chu, K.-H. Stress-strain relationship for plain concrete in compression. J. Proc. 1985, 82, 797–804. [Google Scholar]
  24. Aslani, F.; Jowkarmeimandi, R. Stress–strain model for concrete under cyclic loading. Mag. Concr. Res. 2012, 64, 673–685. [Google Scholar] [CrossRef]
  25. Elwell, D.J.; Fu, G. Compression testing of concrete: Cylinders vs. cubes, Report FHWA/NY/SR-95/119, special report 119, Transportation research and development bureau New York state department of transportation, 1995.
  26. Lubliner, J.; Oliver, J.; Oller, S.; Oñate, E. A plastic-damage model for concrete. Int. J. Solids Struct. 1989, 25, 299–326. [Google Scholar] [CrossRef]
  27. Smith, M. ABAQUS/Standard User’s Manual; Version 6.9; Dassault Systèmes Simulia Corp.: Providence, RI, USA, 2009. [Google Scholar]
  28. Yi, S.-T.; Yang, E.-I.; Choi, J.-C. Effect of specimen sizes, specimen shapes, and placement directions on compressive strength of concrete. Nucl. Eng. Des. 2006, 236, 115–127. [Google Scholar] [CrossRef]
  29. Hosen, M.A.; Alengaram, U.J.; Jumaat, M.Z.; Sulong, N.H.R.; Darain, K.M.U. Glass Fiber Reinforced Polymer (GFRP) Bars for Enhancing the Flexural Performance of RC Beams Using Side-NSM Technique. Polymers 2017, 9, 180. [Google Scholar] [CrossRef]
  30. Zhou, F.; Zhang, J.; Song, S.; Yang, D.; Wang, C. Effect of Temperature on Material Properties of Carbon Fiber Reinforced Polymer (CFRP) Tendons: Experiments and Model Assessment. Materials 2019, 12, 1025. [Google Scholar] [CrossRef]
  31. Paegle, I.; Fischer, G. Phenomenological interpretation of the shear behavior of reinforced Engineered Cementitious Composite beams. Cem. Concr. Compos. 2016, 73, 213–225. [Google Scholar] [CrossRef]
  32. Meng, D.; Lee, C.; Zhang, Y. Flexural and shear behaviours of plain and reinforced polyvinyl alcohol-engineered cementitious composite beams. Eng. Struct. 2017, 151, 261–272. [Google Scholar] [CrossRef]
  33. Ding, Y.; Yu, K.Q.; Yu, J.T.; Xu, S.L. Structural behaviors of ultra-high performance engineered cementitious composites (UHP-ECC) beams subjected to bending-experimental study. Constr. Build. Mater. 2018, 177, 102–115. [Google Scholar] [CrossRef]
  34. OEl-Mahdy, O.; AHamdy, G.; HEl-Diasity, M.; Shalaby, Y. Performance of Reinforced Engineered Cementitious Composite Beams. Eng. Res. J.-Fac. Eng. 2022, 51, 83–101. [Google Scholar] [CrossRef]
Figure 1. The schematic diagram of the composite hybrid reinforced beam.
Figure 1. The schematic diagram of the composite hybrid reinforced beam.
Fibers 12 00049 g001
Figure 2. Stress–strain curves of FRP bars [22,29,30].
Figure 2. Stress–strain curves of FRP bars [22,29,30].
Fibers 12 00049 g002
Figure 3. Finite element model: (a) Mesh; (b) Embedded constraint between reinforcement and concrete /ECC; (c) Hinge support; (d) Load application.
Figure 3. Finite element model: (a) Mesh; (b) Embedded constraint between reinforcement and concrete /ECC; (c) Hinge support; (d) Load application.
Fibers 12 00049 g003aFibers 12 00049 g003b
Figure 4. Load-deflection curves for Group I: (a) Beam B1; (b) Beam B2; (c) Beam B3; (d) Beam B4.
Figure 4. Load-deflection curves for Group I: (a) Beam B1; (b) Beam B2; (c) Beam B3; (d) Beam B4.
Fibers 12 00049 g004
Figure 5. Load-deflection curves for Group II: (a) Beam B5; (b) Beam B6; (c) Beam B7; (d) Beam B8.
Figure 5. Load-deflection curves for Group II: (a) Beam B5; (b) Beam B6; (c) Beam B7; (d) Beam B8.
Fibers 12 00049 g005aFibers 12 00049 g005b
Figure 6. Load-deflection curves for Group III: (a) Beam B9; (b) Beam B10; (c) Beam B11; (d) Beam B12.
Figure 6. Load-deflection curves for Group III: (a) Beam B9; (b) Beam B10; (c) Beam B11; (d) Beam B12.
Fibers 12 00049 g006
Figure 7. Experiment [22] vs. NLFEM cracks patterns and failure mode: (a) Beam B5; (b) Axial stress in reinforcement.
Figure 7. Experiment [22] vs. NLFEM cracks patterns and failure mode: (a) Beam B5; (b) Axial stress in reinforcement.
Fibers 12 00049 g007aFibers 12 00049 g007b
Figure 8. Experiment [22] vs. NLFEM cracks patterns and failure mode for Beam B6.
Figure 8. Experiment [22] vs. NLFEM cracks patterns and failure mode for Beam B6.
Fibers 12 00049 g008
Figure 9. Experiment [22] vs. NLFEM cracks patterns and failure mode for beam B7.
Figure 9. Experiment [22] vs. NLFEM cracks patterns and failure mode for beam B7.
Fibers 12 00049 g009
Figure 10. Experiment [22] vs. NLFEM cracks patterns and failure mode for beam B8.
Figure 10. Experiment [22] vs. NLFEM cracks patterns and failure mode for beam B8.
Fibers 12 00049 g010
Figure 11. Load mid-span deflection curve of the composite hybrid beam: (a) f’c = 24 MPa; (b) f’c = 30 MPa; (c) f’c = 40 MPa.
Figure 11. Load mid-span deflection curve of the composite hybrid beam: (a) f’c = 24 MPa; (b) f’c = 30 MPa; (c) f’c = 40 MPa.
Fibers 12 00049 g011
Figure 12. Normalized value of the mechanical properties of the composite hybrid beams with varying FRP type and grade of concrete: (a) f’c = 24 MPa; (b) f’c = 30 MPa; (c) f’c = 40 MPa.
Figure 12. Normalized value of the mechanical properties of the composite hybrid beams with varying FRP type and grade of concrete: (a) f’c = 24 MPa; (b) f’c = 30 MPa; (c) f’c = 40 MPa.
Fibers 12 00049 g012aFibers 12 00049 g012b
Figure 13. Effect of hybrid reinforcement ratio and ECC thickness on the ultimate load: (a) ρFRP = 0.38; (b) ρFRP = 0.86; (c) ρFRP = 1.53.
Figure 13. Effect of hybrid reinforcement ratio and ECC thickness on the ultimate load: (a) ρFRP = 0.38; (b) ρFRP = 0.86; (c) ρFRP = 1.53.
Fibers 12 00049 g013
Figure 14. Effect of hybrid reinforcement ratio and ECC thickness on the ductility index: (a) ρFRP = 0.38; (b) ρFRP = 0.86; (c) ρFRP = 1.53.
Figure 14. Effect of hybrid reinforcement ratio and ECC thickness on the ductility index: (a) ρFRP = 0.38; (b) ρFRP = 0.86; (c) ρFRP = 1.53.
Fibers 12 00049 g014
Figure 15. Effect of hybrid reinforcement ratio and ECC thickness on the energy absorption: (a) ρFRP = 0.38; (b) ρFRP = 0.86; (c) ρFRP = 1.53.
Figure 15. Effect of hybrid reinforcement ratio and ECC thickness on the energy absorption: (a) ρFRP = 0.38; (b) ρFRP = 0.86; (c) ρFRP = 1.53.
Fibers 12 00049 g015
Figure 16. Effect of hybrid reinforcement ratio and ECC thickness on the stiffness: (a) ρFRP = 0.38; (b) ρFRP = 0.86; (c) ρFRP = 1.53.
Figure 16. Effect of hybrid reinforcement ratio and ECC thickness on the stiffness: (a) ρFRP = 0.38; (b) ρFRP = 0.86; (c) ρFRP = 1.53.
Fibers 12 00049 g016
Table 1. Concrete plastic damage parameters, [26].
Table 1. Concrete plastic damage parameters, [26].
ParameterValue
Flow potential eccentricity (ϵ)0.1
Dilation angle (ψ)30
Compressive meridian at in initial yield (k)2/3
Viscosity (μ)0.0005
Initial biaxial compressive yield stress to initial uniaxial compressive yield stress (fb0/fc0)1.16
Table 2. Tensile properties of top and bottom steel bars [22].
Table 2. Tensile properties of top and bottom steel bars [22].
Bars TypeDiameter (mm)Elastic Modulus (GPa)fy (MPa)fu (MPa)Ultimate Strain
Top steel bar101984034950.051
Bottom steel bar121994085030.0675
Stirrups81984064850.075
Table 3. Properties of FRP reinforcement [22,29,30].
Table 3. Properties of FRP reinforcement [22,29,30].
Bars TypeUltimate Tensile Stress (MPa)Elastic Modulus (GPa)
CFRP2070156
AFRP117283
BFRP125050
GFRP108040
Table 4. Tested beam reinforcement details and ECC height [22].
Table 4. Tested beam reinforcement details and ECC height [22].
GroupBeam NotationAs (Steel)Afrp (BFRP)ECC Layer Thickness (mm)
IB12ϕ12 (226.2 mm2)_0
B22ϕ12 (226.2 mm2)_50
B32ϕ12 (226.2 mm2)_100
B42ϕ12 (226.2 mm2)_200
IIB51 ϕ12 (113.1 mm2)2ϕ8 (100.53 mm2)0
B61 ϕ12 (113.1 mm2)2ϕ8 (100.53 mm2)50
B71 ϕ12 (113.1 mm2)2ϕ8 (100.53 mm2)100
B81 ϕ12 (113.1 mm2)2ϕ8 (100.53 mm2)200
IIIB9_3ϕ8 (150.80 mm2)0
B10_3ϕ8 (150.80 mm2)50
B11_3ϕ8 (150.80 mm2)100
B12_3ϕ8 (150.80 mm2)200
Note: As: is the cross-sectional area of steel reinforcement and Afrp: is the cross-sectional area of the BFRP fiber reinforcement.
Table 5. Validation results for ultimate load and deflection.
Table 5. Validation results for ultimate load and deflection.
GroupBeam NotationUltimate Deflection (mm)Ultimate Load (kN)
Experiment [22]NLFEA%ErrorExperiment [22]NLFEA%Error
IB129.532.058.6438.636.984.2
B222.023.094.953939.832.13
B319.119.994.6639.641.284.24
B433.335.155.564545.681.51
IIB531.133.297.0442.442.610.5
B641.541.470.0753.654.792.22
B729.530.894.5147.245.593.41
B847.950.144.6854.258.738.36
IIIB931.533.416.0635.638.798.96
B1029.731.365.5940.841.391.45
B1128.129.675.5944.643.881.61
B1235.036.925.494549.08.89
Overall 5.25 3.96
Table 6. Summary of the investigated parameters.
Table 6. Summary of the investigated parameters.
Modelfʹc (MPa)Beam Height (mm)ρs(%)FRP Typeρfrp (%)tECC (mm)
M1242000.43BFRP0.38100
M2242000.43CFRP0.38100
M3242000.43AFRP0.38100
M4242000.43GFRP0.38100
M5302000.43BFRP0.38100
M6302000.43CFRP0.38100
M7302000.43AFRP0.38100
M8302000.43GFRP0.38100
M9402000.43BFRP0.38100
M10402000.43CFRP0.38100
M11402000.43AFRP0.38100
M12402000.43GFRP0.38100
M13242500.43BFRP0.38100
M14243000.43BFRP0.38100
M15242000.19BFRP0.3850
M16242000.19BFRP0.8650
M17242000.19BFRP1.5350
M18242000.43BFRP0.3850
M19242000.43BFRP0.8650
M20242000.43BFRP1.5350
M21242000.77BFRP0.3850
M22242000.77BFRP0.8650
M23242000.77BFRP1.5350
M24242001.22BFRP0.3850
M25242001.22BFRP0.8650
M26242001.22BFRP1.5350
M27242000.19BFRP0.38100
M28242000.19BFRP0.86100
M29242000.19BFRP1.53100
M30242000.43BFRP0.86100
M31242000.43BFRP1.53100
M32242000.77BFRP0.38100
M33242000.77BFPR0.86100
M34242000.77BFPR1.53100
M35242001.22BFPR0.38100
M36242001.22BFPR0.86100
M37242001.22BFPR1.53100
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shbeeb, N.I.; Barham, W.S.; Alyahya, W. Numerical Simulation of Engineering Cementitious Composite Beams Strengthened with Fiber-Reinforced Polymer and Steel Bars. Fibers 2024, 12, 49. https://0-doi-org.brum.beds.ac.uk/10.3390/fib12060049

AMA Style

Shbeeb NI, Barham WS, Alyahya W. Numerical Simulation of Engineering Cementitious Composite Beams Strengthened with Fiber-Reinforced Polymer and Steel Bars. Fibers. 2024; 12(6):49. https://0-doi-org.brum.beds.ac.uk/10.3390/fib12060049

Chicago/Turabian Style

Shbeeb, Nadim I., Wasim S. Barham, and Wala’a Alyahya. 2024. "Numerical Simulation of Engineering Cementitious Composite Beams Strengthened with Fiber-Reinforced Polymer and Steel Bars" Fibers 12, no. 6: 49. https://0-doi-org.brum.beds.ac.uk/10.3390/fib12060049

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop