Glass FRP-Reinforced Geopolymer Based Columns Comprising Hybrid Fibres: Testing and FEA Modelling

This study seeks to evaluate the effectiveness of glass-FRP-reinforced geopolymer concrete columns integrating hybrid fibres (GFGC columns) and steel bar-reinforced geopolymer concrete columns incorporating hybrid fibres (SFGC columns) under eccentric and concentric loadings. Steel fibre (SF) and polypropylene fibres (PF) are two types of fibres that are mixed into hybrid fibre-reinforced geopolymer concrete (HFRGC). Eighteen circular concrete columns with a cross-section of 300 mm × 1200 mm were cast and examined under axial loading up to failure. Nine columns were cast with glass-FRP rebars, whereas the other nine were cast with steel rebars. Using ABAQUS, a nonlinear finite element model was established for the GFGC and SFGC columns. The HFRGC material was modelled using a simplified concrete damage plasticity model, whereas the glass-FRP material was simulated as a linear elastic material. It was observed that GFGC columns had up to 20% lower axial strength (AST) and up to 24% higher ductility indices than SFGC columns. The failure modes of both GFGC and SFGC columns were analogous. Both GFGC and SFGC columns revealed the same effect of eccentricity in the form of a decline in AST. A novel statistical model was suggested for predicting the AST of GFGC columns. The outcomes of the experiments, finite element simulations, and theoretical results show that the models can accurately determine the AST of GFGC columns.


Introduction
The processing of Portland cement releases carbon dioxide as a byproduct. The latest studies have used green concrete, called "geopolymer concrete", to minimize the environmental impact of concrete structures [1,2]. Even though many studies have looked into the durability and mechanical properties of geopolymer, the integration of fibre-reinforced polymer (FRP) rebars in geopolymer columns have yet to be thoroughly researched. The critical parameters of reinforced concrete (RC) samples such as ductility and strength may be compromised due to the corrosion of steel and decreased serviceability of concrete structures. FRP rebars are potential substitutes for steel reinforcement. Glass-FRP rebars are gaining regulatory approval and credibility due to their high corrosion resistance, thermal conductivity, better tensile strength, and low density. FRP reinforcement is more effective in corrosive environments since it increases service life, reduces repair and overall construction costs [3,4]. Before failure, FRPs are anisotropic materials with linear elastic properties in compressive as well as flexural members [5][6][7][8][9]. FRPs along with steel rebars augmented the AST and ductility of the concrete columns as compared to plain concrete columns [10,11]. To overcome the brittleness of the concrete material, steel and polypropylene fibres were integrated into geopolymer concrete and glass-FRP rebars in this work.
Researchers in structural engineering have recently become interested in the excellent performance of FRP rebars in reinforced concrete construction. In comparison to steel-reinforced columns, glass-FRP-reinforced concrete columns show higher lateral deflection [12]. Glass-FRP stirrups usually fail due to bond-slip in the columns [13]. Due to their linear elastic nature, the interaction graphs of FRP-RC columns do not depict the balancing points until failure [14]. Furthermore, the procedures for determining the minimum reinforcement ratio of FRP-RC columns have been provided in the literature to avoid brittletension failure of FRP rebars on the tension side [15]. The structural response of glass-FRP and steel-reinforced concrete columns is the same except for AST. The AST of glass-FRP reinforced columns is 7% less compared to steel reinforced columns [16]. The applications of glass-FRP rebars in the compression members perform better with concrete than that in the tensile members, so using glass-FRP rebars as a compression reinforcement in concrete is preferable [17]. In circular columns, the exchange of steel reinforcements for an equivalent number of glass-FRP rebars reduces AST and bending moment while marginally increasing specimen ductility under several loading scenarios [18]. Columns reinforced with glass-FRP rebars must have enough lateral confinement to bear maximum loads equal to or greater than steel-reinforced columns [19]. The failure efficiency of full-scale glass-FRP-reinforced columns is different depending on lateral volumetric reinforcement. The test results show that large-scale glass-FRP-reinforced columns fail at 0.7% due to longitudinal bar buckling, and at 2.7%, volumetric lateral reinforcement fails due to core concrete crushing and lateral confinement rupture [20]. Different empirical methods for estimating the compressive strength of glass-FRP-reinforced columns have been developed by taking into account the contribution of glass-FRP reinforcement [16,18,20]. The literature has shown that the moment and AST of glass-FRP-reinforced circular samples had low AST and moment as compared to steel-reinforced, otherwise identical samples [18]. The number and size of longitudinal glass-FRP rebars were augmented to improve the AST, axial strain, and resistant effectiveness of glass-FRP-reinforced hollow columns [21,22].
Hadi et al. [23] investigated the structural performance of geopolymer compressive samples confined by BFRP tubes and/or basalt-FRP rebars. According to the researchers, steel-reinforced geopolymer compressive samples had an AST that was 27% to 34% higher than that of steel-reinforced ordinary Portland cement (OPC) compressive samples, while their ductility was reduced by 16-27%. However, the geopolymer compressive samples reinforced and confined with basalt-FRP rebars and tubes augmented ductility by 4-7% while the AST was reduced by 5-19% compared to OPC columns confined with the same material. Danda et al. [24] investigated the effect of NaOH molarity using GGBS and found that increasing NaOH molarity enhanced axial strength (AST) and reduced axial deflection of compressive samples. The activation method used for the manufacturing of geopolymers results in the emission of greenhouse gases and more energy consumption, both of which offer a safety risk. Additionally, the curing temperature and duration have a considerable impact on geopolymer synthesis [25,26]. Saranya et al. [27] examined the dynamics of steel fibre-reinforced dolomite-GGBS geopolymer short compressive samples under compressive loading and discovered that they had greater AST than steel fibre-reinforced OPC compressive samples. Furthermore, dolomite-GGBS geopolymer compression samples without ductile detailing to provide appropriate ductility and toughness exhibited similar performance to OPC compression samples, resulting in a lower cost-to-strength ratio. Maranan et al. [28] examined the response of glass-FRP reinforced geopolymer compressive samples, indicating that these compressive samples exceeded their OPC counterparts in terms of compressive strength.
Due to reducing the time and costs of research, the most important methods for a technical aspect to characterize the interaction process, the diversity of the failure mechanism of FRPs, and the unique response between FRPs and concrete are finite element models (FEMs) [29]. Excessive deformations and local buckling regulate the structural integrity of glass-FRP-reinforced samples, according to tests and finite element simulations [30,31]. Numerical modelling can accurately describe the structural efficiency of glass-FRP columns under buckling load [32,33]. Finite element simulations were used to simulate the efficiency of glass-FRP reinforced short concrete columns with varied slenderness ratios [34,35]. Initial failure mode and post-buckling efficiency test results are very comparable to finite element analysis. The numerical models have a high accuracy in evaluating the load-deflection response and failure modes of glass-FRP-reinforced columns [4,36,37]. Even though some authors have attempted to compare the response of glass-FRP-reinforced hybrid fibrereinforced geopolymer concrete (HFRGC) columns to plain concrete, none of them have looked at the efficiency of such elements integrating hybrid fibres. The primary goal of this study is to use investigations and nonlinear finite element simulations to investigate the structural response of GFGC columns under eccentric and concentric stresses.
The SFGC columns were compared to GFGC columns in terms of structural efficiency. Different reinforcing lateral tie spacing, hybrid fibres, loading eccentricity, and loading circumstances were all examined. FEM was given using a modified concrete damaged plasticity (CDP) model. A novel theoretical model for the assessment of the AST of GFGC columns was also proposed. When highly corrosive and light glass-FRP composite materials are utilized in concrete columns, the findings of this study may be valuable to the concrete industry.  Table 1 illustrates many characteristics of coarse aggregates. Coarse aggregates have a maximum size of 10 mm. Lawrencepur sand, which is abundantly available in the area, is used as fine aggregate. The fineness modulus and apparent particle density are 2.25 and 2586 kg/m 3 respectively. The granular analysis of coarse aggregates and sand used in this research is illustrated in Figure 1. The diameters of SF and PF were 0.55 mm and 24 µm, and their lengths were 25 mm and 12 mm respectively. To enhance the geopolymer mix's workability, a superplasticizer named Sika ViscoCrete ® -3425 was used. Using a trial testing approach, the mix design of geopolymer (given in Table 2) with a nominal density of 2400 kg/m 3 was obtained. In GFGC, two readily available waste materials were used as binders, i.e., 45% class F fly ash and 55% GGBS. As an activator, a mass ratio of 1:2.5 was used for a mixture of NaOH (14 M molarity) and Na 2 SiO 3 . The workability test performed according to ASTM C143 [38] reported a 125 mm slump of fresh GFGC. The setting time of the geopolymer was ninety minutes according to ASTM C807-13 [39]. Six concrete cylinders of cross-section 150 × 300 mm were cast. Compressive strength was 32 MPa of the same mix after 28 days.   The longitudinal reinforcement and transverse stirrups for SFGC specimens made from #3 (9.5 mm) and #4 (12.7 mm) steel rebars, respectively. These diamete preferably selected in low-seismicity zones [40,41]. The same number of rebars w plied as a longitudinal and lateral reinforcement for the manufacturing of GFGC mens as in SFGC. E-glass fibres coated with additives, fillers, and thermosetting vin ter resin were used to make the painted glass-FRP rebars. Ligatures were employed investigation, with SupAnchor ® fibres accounting for 80% of the volumetric quantit engineering properties of steel and glass-FRP rebars are summarized in Table 3.

. Steel and Glass-FRP Rebars
The longitudinal reinforcement and transverse stirrups for SFGC specimens were made from #3 (9.5 mm) and #4 (12.7 mm) steel rebars, respectively. These diameters are preferably selected in low-seismicity zones [40,41]. The same number of rebars was applied as a longitudinal and lateral reinforcement for the manufacturing of GFGC specimens as in SFGC. E-glass fibres coated with additives, fillers, and thermosetting vinyl ester resin were used to make the painted glass-FRP rebars. Ligatures were employed in this investigation, with SupAnchor ® fibres accounting for 80% of the volumetric quantity. The engineering properties of steel and glass-FRP rebars are summarized in Table 3.

Specimen Details
Eighteen (18) circular HFRGC columns were cast in this study. Nine of those were used for glass-FRP rebars and the other nine were for used steel rebars. Stirrups were used for lateral reinforcement for all 18 samples. The goal of the specimen testing was to see how the axial efficiency of columns was affected by hybrid fibres, reinforcing type, tie spacing, and loading eccentricity. The SFGC columns were cast as control specimens. All the columns had a cross-section of 300 mm × 1200 mm. These columns can fit into a testing machine but are sizable enough to be referred to as full-size samples. Figure 2 presents the details of the manufactured specimens. Six 12.7 mm-long rebars were used in each GFGC column, resulting in a reinforcement ratio of 1.10, which is ideal for non-seismic low-rise buildings. As a lateral reinforcement, #3 glass-FRP stirrups were used, positioned 75, 150, and 250 mm apart, resulting in volumetric ratios of 1.10%, 0.55%, and 0.32%, respectively. Correspondingly, six #4 long rebars were used to reinforce the SFGC columns, yielding a 1.10% reinforcement ratio. For transverse reinforcement, #3 steel stirrups were employed at a position of 75 mm, 150 mm, and 250 mm, giving volumetric ratios of 1.10%, 0.55%, and 0.32%, respectively. The concrete cover provided to all columns was 20 mm. Table 4 lists the geometric and testing details of the samples. The column framework was fabricated with PVC pipes with an inner diameter of 300 mm and a wall thickness of 5 mm. After the installation of the reinforcement cages, the HFRGC was poured into the mold and vibrated continuously with a poker. The samples were cured at room temperature.
for lateral reinforcement for all 18 samples. The goal of the specimen testing was to see how the axial efficiency of columns was affected by hybrid fibres, reinforcing type, tie spacing, and loading eccentricity. The SFGC columns were cast as control specimens. All the columns had a cross-section of 300 mm x 1200 mm. These columns can fit into a testing machine but are sizable enough to be referred to as full-size samples. Figure 2 presents the details of the manufactured specimens. Six 12.7 mm-long rebars were used in each GFGC column, resulting in a reinforcement ratio of 1.10, which is ideal for non-seismic low-rise buildings. As a lateral reinforcement, #3 glass-FRP stirrups were used, positioned 75, 150, and 250 mm apart, resulting in volumetric ratios of 1.10%, 0.55%, and 0.32%, respectively. Correspondingly, six #4 long rebars were used to reinforce the SFGC columns, yielding a 1.10% reinforcement ratio. For transverse reinforcement, #3 steel stirrups were employed at a position of 75 mm, 150 mm, and 250 mm, giving volumetric ratios of 1.10%, 0.55%, and 0.32%, respectively. The concrete cover provided to all columns was 20 mm. Table 4 lists the geometric and testing details of the samples. The column framework was fabricated with PVC pipes with an inner diameter of 300 mm and a wall thickness of 5 mm. After the installation of the reinforcement cages, the HFRGC was poured into the mold and vibrated continuously with a poker. The samples were cured at room temperature.

Testing Arrangement
The columns were analysed employing compression testing equipment with a load capacity of 5000 kN. All of the specimens were loaded at a rate of 0.003 mm/s. A plaster capping with a wall thickness of 10 mm and a width of 100 mm was provided at the top and bottom to prevent end crushing and non-uniform load distribution on the cross-section during sample testing. A consistent compressive force was delivered to the top surface of the concentrically loaded sample. In an eccentrically loaded sample, as shown in Figure 3, a load was applied at the appropriate eccentricity via a steel rod [41][42][43]. Three LVDTs were installed at 120 • apart. The mechanical properties, such as longitudinal strain, longitudinal deflection, and compressive loads, were documented automatically by using a data logger of the testing machine.

Testing Arrangement
The columns were analysed employing compression testing equipment wit capacity of 5000 kN. All of the specimens were loaded at a rate of 0.003 mm/s. A capping with a wall thickness of 10 mm and a width of 100 mm was provided at and bottom to prevent end crushing and non-uniform load distribution on the cr tion during sample testing. A consistent compressive force was delivered to the face of the concentrically loaded sample. In an eccentrically loaded sample, as sh Figure 3, a load was applied at the appropriate eccentricity via a steel rod [41][42][43] LVDTs were installed at 120° apart. The mechanical properties, such as long strain, longitudinal deflection, and compressive loads, were documented autom by using a data logger of the testing machine.

Finite Element Modelling
This section explains how to forecast the structural performance of GFGC an columns under axial loads using detailed finite element modelling of the com damaging efficiency of HFRGC and glass-FRP rebars. ABAQUS 6.14 was used to finite element simulations of samples. The columns' ultimate failure mechanism, u strength, initial stiffness, crack propagation efficiency, and post-failure stiffness ef are all taken into account in the proposed FEM. The FEM was standardized u

Finite Element Modelling
This section explains how to forecast the structural performance of GFGC and SFGC columns under axial loads using detailed finite element modelling of the complicated damaging efficiency of HFRGC and glass-FRP rebars. ABAQUS 6.14 was used to perform finite element simulations of samples. The columns' ultimate failure mechanism, ultimate strength, initial stiffness, crack propagation efficiency, and post-failure stiffness efficiency are all taken into account in the proposed FEM. The FEM was standardized using the findings of a control specimen's experiment (GFGC75-0). Three-dimensional solid stress sections were used to model the HFRGC. Three-dimensional deformable truss members were utilized to simulate the glass-FRP and steel reinforcing rebars. The complex efficiency of HFRGC was examined by using an improved concrete damaged plasticity model. The compression samples were restrained at the bottom end. The 'embedded region' constraint, which connects the compatible degrees of freedom (DOF) of the wire elements to the needed DOF of the 3D stress elements of concrete, was used to model the contact between the HFRGC and reinforcing material [4]. A uniformly distributed compressive force was supplied to the upper region of the concentrically loaded columns. The displacement control method was utilized to impose a line load of the same magnitude at the required eccentricity, validating the identical loading circumstances used in the experiments. The 'tie' constraint in ABAQUS was applied to fasten 50 mm-thick steel plates at the top and bottom surfaces of the samples for the application of boundary conditions. The boundary conditions that were applied to the simulated specimens were stated in Figure 4. contact between the HFRGC and reinforcing material [4]. A uniformly distributed compressive force was supplied to the upper region of the concentrically loaded columns. The displacement control method was utilized to impose a line load of the same magnitude at the required eccentricity, validating the identical loading circumstances used in the experiments. The 'tie' constraint in ABAQUS was applied to fasten 50 mm-thick steel plates at the top and bottom surfaces of the samples for the application of boundary conditions. The boundary conditions that were applied to the simulated specimens were stated in Figure 4.

Modelling of HFRGC
To anticipate the structural efficiency of concrete elements, a precise definition of HFRGC in FEM is necessary. The elastic efficiency of HFRGC was described in this study using a Poisson's ratio of 0.2 [44] and an elastic modulus of 4734 f c [45]; where f c is the compressive strength of HFRGC, which is 32 MPa. By combining isotropic damaged elasticity, isotropic compressive plasticity, and isotropic tensile plasticity, the CDP model can accurately reflect the nonlinear and irreversible degradation of plain concrete. The CDP model is applied to model the complicated damaging efficiency of HFRGC in the current investigation, with minor adjustments as suggested by Chi et al. [46]. The use of the Chi et al. [46] model for numerical simulations of HFRGC is certified with equitable accuracy after careful review and evaluation between the calculations and experimental results.
The viscosity parameter (µ), yielding of the surface (K (ψ) can all be defined in the CDP model to model the plasticity efficiency of HFRGC. As stated in the validation section of the presented work, the parameter was calibrated to ensure optimum results in comparison to the test data, and the other parameters in the CDP model for HFRGC were changed. For both plain and HFRGC, the default value for the parameter 'e' is always 0.1. The parameter σ h f bo /σ h f co for HFRGC in the CDP model can be represented by Equation (1) [46].
where k c and k t are the constants for compression and tension meridians, respectively. Because of the scarcity of experimental studies in this area, the constant k c has yet to be established for HFRGC. Therefore, the correlation for these constants can be represented by Equations (2) and (3), respectively [44].
In these expressions, λ s f and λ p f are the fibre reinforcement values for SF and PF, respectively which can be defined by the mentioned expressions.
where V s f shows the volumetric fraction of SF, V p f shows the volumetric fraction of PF, l s f shows the length of SF, l p f shows the length of PF, d s f shows the diameter of SF, and d p f shows the diameter of PF. The ratio of length to the diameter of fibres is known as the aspect ratio. After incorporating all the values of several parameters in Equation (1), the obtained value of σ h f bo /σ h f co that was considered in the present work was 1.48. It is should be noted that the value of σ bo /σ co = 1.16 is obtained when k c = 1 and k t = 1 are used to demonstrate the plain concrete.
As a result, the constants' relationships can be represented. The shape factor (K c ) for plain concrete is typically in the range of 0.64 to 0.80 [47]. When plain concrete is subjected to low hydrostatic stresses, the strength experimental measurements are shown to be very close. Furthermore, the value of 0.70 provides a more acceptable correlation when the concrete is subjected to high hydrostatic stresses [46]. As a result, for HFRGC, the relationship for K h f c can be written in terms of the compressive and tensile meridians' enhancing constants [44,46].
Because the fibres are more resistant to crack opening and spreading, the tension meridian will be enhanced more than the compression meridian, i.e., k t k c ≥ 1. Consequently, the inequality K h f c ≥ K c will always be true for the HFRGC. The ψ is an important parameter that controls the flow rule. The value drops in the case of HFRGC, resulting in less plastic volumetric strain in the concrete matrix that effectively confines it. An equation for the dilation angle of HFRGC (ψ h f ) can be expressed in the form of λ s f and λ p f [46].
According to the elastoplasticity hypothesis, which depicts the nonlinearity and irreversible deformation of concrete, the overall strain of concrete (ε) may be separated into two categories: elastic strain (ε el ) and plastic strain (ε pl ).
Concrete nonlinearity can be attributable to either plasticity or damage, alone or a combination of both. Consideration of unloading stiffness is important for the simulation of the non-linearity of concrete. The deterioration of stiffness increases as strain increases, implying that the concrete damaging evolution is linked to the strain elements i.e., ε el and ε pl [46]. By involving two variables known as the uniaxial compression (d c ) and the uniaxial tension parameter (d t ), the CDP model in FEA defines the concrete damage. From Figure 5, the correlation for the compressive stress of plain concrete can be depicted by Equation (9).
where E o is the elastic modulus of concrete defined according to ACI 318-95 [45], ε c is the concrete strain under compressive loading and ε pl c is the plastic portion of the concrete strain under compressive loading. The parameter d c can be stated by Equation (10) as recommended by Wang and Chen [48].
Because the fibres are more resistant to crack opening and spreading, the tension meridian will be enhanced more than the compression meridian, i.e., ≥ 1. Consequently, the inequality ≥ will always be true for the HFRGC. The is an important parameter that controls the flow rule. The value drops in the case of HFRGC, resulting in less plastic volumetric strain in the concrete matrix that effectively confines it. An equation for the dilation angle of HFRGC ( ) can be expressed in the form of and [46].
According to the elastoplasticity hypothesis, which depicts the nonlinearity and irreversible deformation of concrete, the overall strain of concrete ( ) may be separated into two categories: elastic strain ( ) and plastic strain ( ).
Concrete nonlinearity can be attributable to either plasticity or damage, alone or a combination of both. Consideration of unloading stiffness is important for the simulation of the non-linearity of concrete. The deterioration of stiffness increases as strain increases, implying that the concrete damaging evolution is linked to the strain elements i.e., and [46]. By involving two variables known as the uniaxial compression ( ) and the uniaxial tension parameter ( ), the CDP model in FEA defines the concrete damage. From Figure 5, the correlation for the compressive stress of plain concrete can be depicted by Equation (9).
where is the elastic modulus of concrete defined according to ACI 318-95 [45], is the concrete strain under compressive loading and is the plastic portion of the concrete strain under compressive loading. The parameter can be stated by Equation (10) as recommended by Wang and Chen [48].
where m c is the parameter that supervises the speed of compressive damage evolution having a value of 0.1 for the normal concrete [49], ε in c, norm is the compressive inelastic concrete strain that could be expressed as: ε in c /ε in cu . Here, ε in cu is the compressive strain corresponding to the inelastic concrete strain with a value of 0.033 [49]. Therefore, the modified compressioncontrolling parameter for HFRGC (m h f c ) can be expressed by Equation (11) [46].
where a m1 and b m1 are the constants linked to the characteristic parameter of fibres with magnitudes of 0.452 and 0.054, respectively [46].  (12) and (13), respectively.
The ultimate compressive stress (σ h f co ) and corresponding compressive strain (ε h f co ) of HFRGC can be expressed by the following relations [46].
A strain-hardening branch of the stress-strain curve characterizes the post-cracking efficiency of concrete. Strain-hardening represents the pre-ultimate efficiency of the stressstrain curve through plastic flow, while strain-softening represents the post-ultimate efficiency [50]. Figure 6 represents the tensile efficiency of concrete. The equation for tensile stress of concrete (σ t ) can be presented as: where ε t shows the concrete strain under tension and ε pl t is the plastic region of the concrete strain under tension. The parameter d t can be represented by Equation (17) as proposed by Wang and Chen [48].
where m t is the variable that controls the tensile damaging evolution speed, having a magnitude of 0.05 for plain concrete [49], and ε ck t, norm is the normalized inelastic concrete strain under tension that might be expressed as: ε ck t /ε ck tu . Here, ε ck tu is the strain corresponding to the normalized inelastic concrete strain under tension having a value of 0.1ε in cu = 0.0033 [49]. As the speed of crack initiation and propagation is reduced in HFRGC, the modified tension controlling parameter for HFRGC (m h f t ) can be expressed by Equation (18) [46].
In this relationship, and are the variables linked to the characteristic parameter of PF and SF, with values of 0.156 and 0.628, respectively [46]. Thus, for HFRGC, the tension damage parameter and tensile stress ( ) can be expressed by Equation (19) and Equation (20), respectively. In this relationship, a m2 and b m2 are the variables linked to the characteristic parameter of PF and SF, with values of 0.156 and 0.628, respectively [46]. Thus, for HFRGC, the tension damage parameter d The ultimate tensile stress (σ h f to ) and corresponding tensile strain (ε h f to ) of HFRGC can be expressed by the following relations [46].
By considering the above relationship, the yield stress and the inelastic strain of HFRGC can be calculated for introduction into the FEA models.

Modelling of Glass-FRP and Steel Rebars
The glass-FRP and steel rebars were simulated as 3D truss elements. Table 3 shows the elastic moduli, yield strengths, and cross-sectional details of the reinforcements. For steel reinforcement, the Poisson's ratio was 0.3 [5,51], while for glass-FRP reinforcement it was 0.25 [41]. The response of glass-FRP rebars was simulated up to failure by treating them as a linear elastic material [52][53][54]. Steel rebars were treated as an elastic-plastic material with comparable compression and tension properties to simulate their efficiency [51,55] as expressed in Figure 7. Compressive stress was found to be 55% of tensile strength in glass-FRP rebars [56].

Calibration of FEA Model
The experimental results were used to check the given FEM of the referenced specimen (GFGC75-0). The proposed FEM was scaled for different end conditions, HFRGC's viscosity parameter (µ), multiple element types, and element mesh sizes to generate finite element assumptions that were close to the experiments. A similar relationship with the load-deflection efficiency of the reference sample was seen when all the DOF at the lower

Calibration of FEA Model
The experimental results were used to check the given FEM of the referenced specimen (GFGC75-0). The proposed FEM was scaled for different end conditions, HFRGC's viscosity parameter (µ), multiple element types, and element mesh sizes to generate finite element assumptions that were close to the experiments. A similar relationship with the loaddeflection efficiency of the reference sample was seen when all the DOF at the lower end of the samples was controlled while the upper end was able to translate and rotate on any axis. [57]. The effect of the parameter on the load-deflection response was observed to be remarkable, as shown in Figure 8a. The value of 0.0048 for µ demonstrated close agreement with the test results. efficiency of steel rebars.

Calibration of FEA Model
The experimental results were used to check the given FEM of the referenced specimen (GFGC75-0). The proposed FEM was scaled for different end conditions, HFRGC's viscosity parameter (µ), multiple element types, and element mesh sizes to generate finite element assumptions that were close to the experiments. A similar relationship with the load-deflection efficiency of the reference sample was seen when all the DOF at the lower end of the samples was controlled while the upper end was able to translate and rotate on any axis. [57]. The effect of the parameter on the load-deflection response was observed to be remarkable, as shown in Figure 8a. The value of 0.0048 for demonstrated close agreement with the test results.  The mesh size was modified to provide the most optimum finite element calculations in a suitably limited time of FEA to minimize disruption in the circulation of strains and stresses in the HFRGC and glass-FRP rebars. The sensitivity of the load-deflection curve as a function of mesh size was shown in Figure 8b. Mesh sizes from 20 to 70 mm were investigated in this study. A similar agreement between the FEM model and experimental results was noted for a 40 mm screen size. The load-deflection response was used to investigate the effect of several kinds of the element of HFRGC and glass-FRP rebars in the control model. The studied kinds of the element of HFRGC were C3D8R and C3D20R for hexahedral elements, C3D4H and C3D10H for tetrahedral elements, and C3D6H and C3D15H for triangular elements, with C3D8R showing the highest accuracy as reported in previous works [4,37,46]. The glass-FRP and steel rebars were examined for T3D3H and T3D2H elements. A close correlation for the test results was observed when T3D2H elements were used for the reinforcing rebars. The effect of several element types is demonstrated in Figure 9. The mesh size was modified to provide the most optimum finite element calculations in a suitably limited time of FEA to minimize disruption in the circulation of strains and stresses in the HFRGC and glass-FRP rebars. The sensitivity of the load-deflection curve as a function of mesh size was shown in Figure 8b. Mesh sizes from 20 to 70 mm were investigated in this study. A similar agreement between the FEM model and experimental results was noted for a 40 mm screen size. The load-deflection response was used to investigate the effect of several kinds of the element of HFRGC and glass-FRP rebars in the control model. The studied kinds of the element of HFRGC were C3D8R and C3D20R for hexahedral elements, C3D4H and C3D10H for tetrahedral elements, and C3D6H and C3D15H for triangular elements, with C3D8R showing the highest accuracy as reported in previous works [4,37,46]. The glass-FRP and steel rebars were examined for T3D3H and T3D2H elements. A close correlation for the test results was observed when T3D2H elements were used for the reinforcing rebars. The effect of several element types is demonstrated in Figure 9.

Failure Modes and Cracking Patterns
All the columns were noted to fail in the upper half. Figure 10 portrays the failure modes of some specimens after the test. Both groups of samples, i.e., GFGC and SFGC specimens, show similar failure modes and processes. Fracture in the longitudinal reinforcement is the major reason for GFGC column failure, while SFGC columns fail due to

Failure Modes and Cracking Patterns
All the columns were noted to fail in the upper half. Figure 10 portrays the failure modes of some specimens after the test. Both groups of samples, i.e., GFGC and SFGC specimens, show similar failure modes and processes. Fracture in the longitudinal reinforcement is the major reason for GFGC column failure, while SFGC columns fail due to the buckling of long rebars. For the columns with a smaller distance between stirrups, the failure of the specimens was understood as the rupturing of stirrups and damage to the core. Conversely, the buckling of longitudinal rebars caused damage to the samples with a larger distance between stirrups. Up to 80% of ultimate load specimens behaved elastically after being loaded due to transverse confinement. The cracking of HFRGC started with a small vertical hairline fracture that was propagated from the upper corner of the columns after a further increase in load. Before reaching the ultimate load, cracks start propagating and their width expands. Cracks started to propagate along with the height of the specimen, and deflection augmented progressively after reaching the ultimate load. At this moment, the concrete cover was expanded, and the lateral stirrups were activated. When the AST of specimens decreased up to about 60% of the ultimate load, the lateral confinement reached the limit of its strength, triggering a rupture sound in the GFGC samples. The SFGC specimens triggered no rupturing sound of the steel rebars. The geopolymer concrete and hybrid fibres demonstrated a good bonding effect with the reinforcement rebars. When the load is applied eccentrically on GFGC and SFGC samples, crushing of the concrete during failure starts from the compression side, while minor cracks appear on the tension side. Due to tensile forces acting normal to the plane of cracking and the bridging effect of fibres, concrete remains uncrushed from the tension side and no serious damage occurs in the concrete cover [58].

Ductility Index
The rotation, axial deformation, curvature, strain, absorbed energy, and dissipated energy of elements are some of the variables that can be used to express ductility. The ductility value of the GFGC and SFGC samples was calculated utilising the following equation in this research [42,59,60]: where ∆ is the value of the axial deflection of the specimen in the elastic section of the load-deflection curve corresponding to 75% of the peak load, and ∆ is the area under the axial load-deflection curve up to 75%. As illustrated in Figure 11, ∆ is the area under the axial load-deflection curve up to ∆ , where ∆ is the value of the axial deflection of the columns in the post-ultimate collapse curve of the load-deflection corresponding to 85% of the peak load.  When the load is applied eccentrically on GFGC and SFGC samples, crushing of the concrete during failure starts from the compression side, while minor cracks appear on the tension side. Due to tensile forces acting normal to the plane of cracking and the bridging effect of fibres, concrete remains uncrushed from the tension side and no serious damage occurs in the concrete cover [58].

Ductility Index
The rotation, axial deformation, curvature, strain, absorbed energy, and dissipated energy of elements are some of the variables that can be used to express ductility. The ductility value of the GFGC and SFGC samples was calculated utilising the following equation in this research [42,59,60]: where ∆ 75 is the value of the axial deflection of the specimen in the elastic section of the load-deflection curve corresponding to 75% of the peak load, and A ∆75 is the area under the axial load-deflection curve up to 75%. As illustrated in Figure 11, A ∆85 is the area under the axial load-deflection curve up to ∆ 85 , where ∆ 85 is the value of the axial deflection of the columns in the post-ultimate collapse curve of the load-deflection corresponding to 85% of the peak load. When the load is applied eccentrically on GFGC and SFGC samples, crushing of the concrete during failure starts from the compression side, while minor cracks appear on the tension side. Due to tensile forces acting normal to the plane of cracking and the bridging effect of fibres, concrete remains uncrushed from the tension side and no serious damage occurs in the concrete cover [58].

Ductility Index
The rotation, axial deformation, curvature, strain, absorbed energy, and dissipated energy of elements are some of the variables that can be used to express ductility. The ductility value of the GFGC and SFGC samples was calculated utilising the following equation in this research [42,59,60]: where ∆ is the value of the axial deflection of the specimen in the elastic section of the load-deflection curve corresponding to 75% of the peak load, and ∆ is the area under the axial load-deflection curve up to 75%. As illustrated in Figure 11, ∆ is the area under the axial load-deflection curve up to ∆ , where ∆ is the value of the axial deflection of the columns in the post-ultimate collapse curve of the load-deflection corresponding to 85% of the peak load.  Figure 11. Ductility index measurement for specimens. Figure 11. Ductility index measurement for specimens.
According to the literature, the GFGC specimens were found to be more ductile than identical SFGC specimens [42]. On average, the ductility values of GFGC columns were 24% greater than those of SFGC columns. The ductility values for all of the GFGC and SFGC columns that were evaluated are shown in Figure 12. The columns with a decreased stirrup distance had greater ductility indices, as expected. The average ductility indices were 2.66, 2.33, and 2.26 for the GFGC columns having a ligature distance of 75 mm, 150 mm, and 250 mm, respectively. Likewise, the ductility indices for the SFGC columns with a ligature distance of 75 mm, 150 mm, and 250 mm were 2.33, 1.87, and 1.66, respectively. The confinement of the concrete core is attributed to the reduced stirrup distance which restrained the long rebars and caused the improvement in ductility by absorbing more energy [60]. Additionally, the enhancement of the loading eccentricity of GFGC and SFGC columns is directly proportional to the ductility which can be observed from Figure 12. According to the literature, the GFGC specimens were found to be more ductile than identical SFGC specimens [42]. On average, the ductility values of GFGC columns were 24% greater than those of SFGC columns. The ductility values for all of the GFGC and SFGC columns that were evaluated are shown in Figure 12. The columns with a decreased stirrup distance had greater ductility indices, as expected. The average ductility indices were 2.66, 2.33, and 2.26 for the GFGC columns having a ligature distance of 75 mm, 150 mm, and 250 mm, respectively. Likewise, the ductility indices for the SFGC columns with a ligature distance of 75 mm, 150 mm, and 250 mm were 2.33, 1.87, and 1.66, respectively. The confinement of the concrete core is attributed to the reduced stirrup distance which restrained the long rebars and caused the improvement in ductility by absorbing more energy [60]. Additionally, the enhancement of the loading eccentricity of GFGC and SFGC columns is directly proportional to the ductility which can be observed from Figure 12.

Reinforcement Type Effect
The SFGC columns had 20% greater AST compared to the GFGC columns. The SFGC75-0 carried a peak load of 1949 kN at an axial deformation of 5.54 mm. The concentric SFGC column had an AST of 2240 kN with a center-to-center distance of 75 mm. In comparison to the concentric SFGC75 column, the identical GFGC75 column showed a decrease in AST, while its deflection was 14.9% compared to 6.6% in the SFGC75 column respectively. The AST of the concentrically loaded GFGC150 and GFGC250 exhibited values of 1651 kN and 1775 kN respectively, which were 74% and 87% of the AST of identical SFGC columns, respectively. The average axial deformation of the SFGC columns at the ultimate AST was 7% larger than that of the GFGC columns.
When compared to their GFGC counterparts, all of the eccentric SFGC columns had higher AST. The samples GFGC75-35 carry an AST of 1372 kN that was 23% lower than SFGC75-35. When loading eccentricity was augmented, GFGC columns had lower AST than SFGC identical columns. The AST of the GFGC specimen with a 150 mm stirrup distance was up to 54.5% lower while eccentricity was ameliorated from 35 mm to 70

Reinforcement Type Effect
The SFGC columns had 20% greater AST compared to the GFGC columns. The SFGC75-0 carried a peak load of 1949 kN at an axial deformation of 5.54 mm. The concentric SFGC column had an AST of 2240 kN with a center-to-center distance of 75 mm. In comparison to the concentric SFGC75 column, the identical GFGC75 column showed a decrease in AST, while its deflection was 14.9% compared to 6.6% in the SFGC75 column respectively. The AST of the concentrically loaded GFGC150 and GFGC250 exhibited values of 1651 kN and 1775 kN respectively, which were 74% and 87% of the AST of identical SFGC columns, respectively. The average axial deformation of the SFGC columns at the ultimate AST was 7% larger than that of the GFGC columns.
When compared to their GFGC counterparts, all of the eccentric SFGC columns had higher AST. The samples GFGC75-35 carry an AST of 1372 kN that was 23% lower than SFGC75-35. When loading eccentricity was augmented, GFGC columns had lower AST than SFGC identical columns. The AST of the GFGC specimen with a 150 mm stirrup distance was up to 54.5% lower while eccentricity was ameliorated from 35 mm to 70 mm. Identical SFGC columns reinforced with steel rebars showed a 39.3% reduction in AST when the eccentricity was augmented from 35 mm to 70 mm. As a result, the literature for plain concrete columns reveals that the SFGC columns surpass the GFGC columns in terms of AST [18,42].

Stirrup Spacing Effect
The literature shows that reductions in stirrup distance cause an increase in the AST of both SFGC and GFGC specimens [42,60]. The AST results for eccentric and concentric loading show good agreement with the literature. For concentrically loaded SFGC and GFGC columns, reducing the stirrup distance from 150 mm (by volume of 0.55%) to 75 mm (by volume of 1.10%) resulted in a percentage increase of 7.7% and 18% in AST, respectively. The confinement of the concrete core and longitudinal bar increases as the stirrup distance decreases, which causes the increment in AST of the specimens [60]. For SFGC columns, a percentage increase of 3.5% was seen when the stirrup distance was reduced from 250 mm (by volume of 0.32%) to 150 mm (by volume of 0.55%); however, the effect was the opposite for GFGC columns. AST and axial deflections were reduced when stirrup spaces were decreased for GFGC columns. In this research, AST was reduced by 6.9% when the distance was reduced from 250 mm to 150 mm for GFGC samples. Though the smaller stirrup distance is more preferable for glass-FRP-reinforced columns, this reduction in the AST may be ascribed to issues with testing the fabrication of the GFGC columns with a 75 mm stirrup distance.

Eccentricity Effect
The AST of all specimens was remarkably reduced as a result of the loading eccentricity. Figure 13 reflects the extent of eccentricity on the AST of SFGC and GFGC columns. When eccentricity was enriched from 0 mm to 35 mm and 70 mm, the AST of GFGC75 specimens decreased by 42% and 120%, respectively. When the eccentricity was adjusted by the same numbers, the AST of the specimens GFGC150 decreased by 35.7% and 109.7%, respectively. Similarly, raising the eccentricity from 0 to 35 and 70 mm reduced the AST of GFGC250 samples by 42.7% and 111.5%, respectively. when the eccentricity was augmented from 35 mm to 70 mm. As a result, the literature for plain concrete columns reveals that the SFGC columns surpass the GFGC columns in terms of AST [18,42].

Stirrup Spacing Effect
The literature shows that reductions in stirrup distance cause an increase in the AST of both SFGC and GFGC specimens [42,60]. The AST results for eccentric and concentric loading show good agreement with the literature. For concentrically loaded SFGC and GFGC columns, reducing the stirrup distance from 150 mm (by volume of 0.55%) to 75 mm (by volume of 1.10%) resulted in a percentage increase of 7.7% and 18% in AST, respectively. The confinement of the concrete core and longitudinal bar increases as the stirrup distance decreases, which causes the increment in AST of the specimens [60]. For SFGC columns, a percentage increase of 3.5% was seen when the stirrup distance was reduced from 250 mm (by volume of 0.32%) to 150 mm (by volume of 0.55%); however, the effect was the opposite for GFGC columns. AST and axial deflections were reduced when stirrup spaces were decreased for GFGC columns. In this research, AST was reduced by 6.9% when the distance was reduced from 250 mm to 150 mm for GFGC samples. Though the smaller stirrup distance is more preferable for glass-FRP-reinforced columns, this reduction in the AST may be ascribed to issues with testing the fabrication of the GFGC columns with a 75 mm stirrup distance.

Eccentricity Effect
The AST of all specimens was remarkably reduced as a result of the loading eccentricity. Figure 13 reflects the extent of eccentricity on the AST of SFGC and GFGC columns. When eccentricity was enriched from 0 mm to 35 mm and 70 mm, the AST of GFGC75 specimens decreased by 42% and 120%, respectively. When the eccentricity was adjusted by the same numbers, the AST of the specimens GFGC150 decreased by 35.7% and 109.7%, respectively. Similarly, raising the eccentricity from 0 to 35 and 70 mm reduced the AST of GFGC250 samples by 42.7% and 111.5%, respectively. Similarly to GFGC samples, as the eccentricity of SFGC75 was increased from 0 to 35 mm and 70 mm, the reduction in AST was 32.7% and 114.7%, respectively. As in the GFGC150 samples, eccentricity augmented from 0 to 35 mm and 70 mm, the AST of  Similarly to GFGC samples, as the eccentricity of SFGC75 was increased from 0 to 35 mm and 70 mm, the reduction in AST was 32.7% and 114.7%, respectively. As in the GFGC150 samples, eccentricity augmented from 0 to 35 mm and 70 mm, the AST of SFGC150 was reduced by 37.3% and 91.3%, respectively. Furthermore, raising the eccentricity from 0 mm to 35 mm and 70 mm lowered the AST of the sample SFGC250 by 34.6% and 108.7%, respectively. Due to the loading eccentricity, both SFGC and GFGC specimens showed nearly the same amount of percentage losses in their AST.

Control Model
After confirmation for mesh size, element types, boundary conditions, and other HFRGC parameters, the reference FEM (GFGC75-0) showed variations of 2.6% and 6.8% for maximum AST and corresponding axial deflection from the tests. The elastic behavior, fracture patterns, and post-ultimate collapse efficiency of the reference sample were precisely described by the proposed FEM. Differences in the experimental and FEA testing environments, such as geometric imperfections, boundary conditions, material strength, sample manufacturing flaws, the differences in the features of glass-FRP rebars and steel rebars, as well as testing equipment accuracy, can be held responsible for minor discrepancies in the proposed FEM. The inconsistencies could potentially be attributed to FEA assumptions, such as assuming an ideal connection between HFRGC and rebars. The proposed FEM showed a strong connection with the tests and can be applied to analyse other SFGC and GFGC columns. Table 5 summarizes the research observations and numerical indicators obtained using the proposed FEM for the GFGC and SFGC columns. The SFGC columns have larger ultimate loads than their GFGC counterparts in general. The largest discrepancy between test observations and FEM calculations was obtained for the ultimate AST of SFGC150-70 with a percentage error of 9.3%. The specimen GFGC250-70 had the largest variation in axial deflection at maximum AST with a percentage error of 24.5%. The specimens' initial geometric flaws cause large variations in results that were not accounted for in their FEA models. Even though the FEM was confirmed using the experimental outcomes of a GFGC column (GFGC75-0), it anticipated the axial efficiency of SFGC specimens. The average percentage discrepancies between the experiments and FEA calculations for ultimate strength for GFGC and SFGC columns were 3.6% and 3.6%, respectively, while axial deflection discrepancies at ultimate strength were 7.8% and 6%, respectively. The majority of the FEM results underestimated the specimens' ultimate AST. The AST and deflection of eccentric columns were accurately identified using the recommended FEM. It captures the concentric SFGC column's AST with a percentage inaccuracy of 4.5% and the eccentric GFGC column AST with a percentage inaccuracy of 4%. The tensile efficiency of HFRGC was exactly predicted in this research based on the accuracy of FEM predictions for eccentrically loaded SFGC and GFGC columns. Thus, the model of the tensile efficiency of HFRGC is needed for precise estimates of the eccentric columns. The difference between the experimental outcomes and FEM predictions for all HFRGC columns was 3.6% for ultimate AST and 6.9% for respective axial deflection. As a result of the safe under predictions for the AST of SFGC and GFGC columns, the suggested numerical model can be adopted and endorsed. Figure 14 shows the analytical and simulation load-deflection curves of eccentrically and concentrically loaded GFGC samples. The efficiency of both eccentric and concentric GFGC columns was properly estimated by the recommended FEM in both the elastic and inelastic regions, with only slight deviations. FEM outcomes of load-deflection behavior for some samples, such as GFGC75-70 and GFGC150-35, are stiffer in the elastic range. The slight differences in the post-ultimate collapse efficiency of the samples between the experiments and FEM results could be related to the complicated damage and degradation efficiency of HFRGC after cracks have been started, and the interaction between glass-FRP rebars and HFRGC needs to be adjusted. Nevertheless, the entire load-deflection behavior seems to be realistically monitored by the recommended FEA models. The AST and deflection of eccentric columns were accurately identified using the recommended FEM. It captures the concentric SFGC column's AST with a percentage inaccuracy of 4.5% and the eccentric GFGC column AST with a percentage inaccuracy of 4%. The tensile efficiency of HFRGC was exactly predicted in this research based on the accuracy of FEM predictions for eccentrically loaded SFGC and GFGC columns. Thus, the model of the tensile efficiency of HFRGC is needed for precise estimates of the eccentric columns. The difference between the experimental outcomes and FEM predictions for all HFRGC columns was 3.6% for ultimate AST and 6.9% for respective axial deflection. As a result of the safe under predictions for the AST of SFGC and GFGC columns, the suggested numerical model can be adopted and endorsed. Figure 14 shows the analytical and simulation load-deflection curves of eccentrically and concentrically loaded GFGC samples. The efficiency of both eccentric and concentric GFGC columns was properly estimated by the recommended FEM in both the elastic and inelastic regions, with only slight deviations. FEM outcomes of load-deflection behavior for some samples, such as GFGC75-70 and GFGC150-35, are stiffer in the elastic range. The slight differences in the post-ultimate collapse efficiency of the samples between the experiments and FEM results could be related to the complicated damage and degradation efficiency of HFRGC after cracks have been started, and the interaction between glass-FRP rebars and HFRGC needs to be adjusted. Nevertheless, the entire load-deflection behavior seems to be realistically monitored by the recommended FEA models. The load-deflection efficiency of all SFGC columns, as assessed by experimental studies and finite element simulations, is shown in Figure 15. The recommended FEM exactly predicted the load-deflection efficiency of all SFGC specimens in the elastic range. The samples SFGC75-35, SFGC150-35, and SFGC250-70 had some variations in the post-ultimate damage phase. Minor discrepancies in the proposed FEM can be attributed to differences in the experimental and FEA circumstances such as geometric imperfections, boundary conditions, material strength, specimen manufacturing faults, differences in the characteristics of glass-FRP and steel rebars, and testing equipment correctness. Because the proposed FEM could not successfully predict the post-ultimate efficiency of SFGC250-70, the specimen's strength did not decline when the stirrups were opened. The load-deflection efficiency of all SFGC columns, as assessed by experimental studies and finite element simulations, is shown in Figure 15. The recommended FEM exactly predicted the load-deflection efficiency of all SFGC specimens in the elastic range. The samples SFGC75-35, SFGC150-35, and SFGC250-70 had some variations in the postultimate damage phase. Minor discrepancies in the proposed FEM can be attributed to differences in the experimental and FEA circumstances such as geometric imperfections, boundary conditions, material strength, specimen manufacturing faults, differences in the characteristics of glass-FRP and steel rebars, and testing equipment correctness. Because the proposed FEM could not successfully predict the post-ultimate efficiency of SFGC250-70, the specimen's strength did not decline when the stirrups were opened.

Crack Patterns and Damage Behavior
The experimental failure and FEM crack patterns are shown in Figure 10. Cracks in HFRGC were displayed in the FEM using ultimate positive plastic strains. Since the orientation of the cracks is mostly perpendicular to these plastic strains, the cracking trends can be assessed effectively using these types of strains. Figure 10 shows that the recommended FEM satisfactorily captured the cracking efficiency of SFGC and GFGC specimens. The crack propagation and spalling of the concrete cover of the specimens began when the reinforcement began to yield after the ultimate load. As predicted by the measurements, the bulk of the specimens were damaged in the upper portion. It was also revealed that the recommended FEM accurately predicted the cracking patterns and failure modes of eccentric columns. The elastic modulus of glass-FRP rebars is lower than that of steel rebars, which is comparable to that of concrete. As a result, the glass-FRP rebars and HFRGC combination were more coordinated when the load was applied to GFGC specimens.

FRP-Confining Mechanism
Columns' strength and ductility are enhanced by using glass-FRP stirrups to confine concrete [61]. The confinement effect is minor at the start, but it becomes active at the ultimate strength with increasing lateral strains that tend to result in lateral pressure. The lateral glass-FRP reinforcement prevents concrete dilatation caused by lateral pressure after an axial load is applied. The AST of glass-FRP-reinforced columns was previously computed without taking into account the resistance efficiency of glass-FRP stirrups. The predictions are underestimated since the influence of glass-FRP-confinement on the AST of columns is overlooked [6]. The lateral confining pressure ( ) due to glass-FRP stirrups can be showed by Equation (24) where is the elastic modulus of glass-FRP stirrups, is the volumetric ratio of glass-FRP lateral stirrups, and is the confinement effectiveness constant. As illustrated in Figure 16, Mander et al. [62] suggested a correlation for this constant by incorporating a

Crack Patterns and Damage Behavior
The experimental failure and FEM crack patterns are shown in Figure 10. Cracks in HFRGC were displayed in the FEM using ultimate positive plastic strains. Since the orientation of the cracks is mostly perpendicular to these plastic strains, the cracking trends can be assessed effectively using these types of strains. Figure 10 shows that the recommended FEM satisfactorily captured the cracking efficiency of SFGC and GFGC specimens. The crack propagation and spalling of the concrete cover of the specimens began when the reinforcement began to yield after the ultimate load. As predicted by the measurements, the bulk of the specimens were damaged in the upper portion. It was also revealed that the recommended FEM accurately predicted the cracking patterns and failure modes of eccentric columns. The elastic modulus of glass-FRP rebars is lower than that of steel rebars, which is comparable to that of concrete. As a result, the glass-FRP rebars and HFRGC combination were more coordinated when the load was applied to GFGC specimens.

FRP-Confining Mechanism
Columns' strength and ductility are enhanced by using glass-FRP stirrups to confine concrete [61]. The confinement effect is minor at the start, but it becomes active at the ultimate strength with increasing lateral strains that tend to result in lateral pressure. The lateral glass-FRP reinforcement prevents concrete dilatation caused by lateral pressure after an axial load is applied. The AST of glass-FRP-reinforced columns was previously computed without taking into account the resistance efficiency of glass-FRP stirrups. The predictions are underestimated since the influence of glass-FRP-confinement on the AST of columns is overlooked [6]. The lateral confining pressure ( f l ) due to glass-FRP stirrups can be showed by Equation (24) [62].
where A f t shows the area of glass-FRP stirrups, d c shows the diameter of the core, s shows the pitch glass-FRP stirrups, and f f v shows the stress in glass-FRP stirrups at f co . This can be computed by employing the equation below [63]: where E f is the elastic modulus of glass-FRP stirrups, ρ t is the volumetric ratio of glass-FRP lateral stirrups, and k e is the confinement effectiveness constant. As illustrated in Figure 16, Mander et al. [62] suggested a correlation for this constant by incorporating a bending effect because of the confinement of transverse stirrups. The bending efficiency resembles the form of a second-degree parabola [61,62]. In this equation, is the clear distance between the longitudinal glass-FRP reb is the distance between glass-FRP stirrups, and is the reinforcement ratio of lo glass-FRP rebars in percentage. The functional transverse confining stress due to gla FRP stirrups can be described by Equation (27) as recommended by Mander et al. [62] = ( The parameters and can be calculated using Equation (24) and Equation (2 respectively. Ahmad et al. [64] and Raza et al. [65] recommended the correlations for culating the peak stress and respective strain of concrete constrained with glass-FRP s rups based on the empirical values given by Equation (28) and Equation (29) where is the ultimate stress of concrete constrained with glass-FRP stirrups, is unconfined concrete compressive strain respective to ultimate compressive load and is the maximum concrete compressive strain confined with glass-FRP stirrups. The rameters and are the confined concrete stiffness ratio and strain ratio which be calculated as follows, respectively:

Proposed AST Model
A series of simulations for measuring the compressive strength of FRP-RC columns has been recommended in previous research [16,20,[66][67][68][69][70]. Most of these studies did not include the influence of glass-FRP resistant pressure, except for [66,67] which offered models for glass-FRP-reinforced conventional concrete specimens. The ultimate longitudinal strain of concrete is assumed to be equivalent to the compressive strain of glass-FRP rebars in a few of the strength models [66,68]. The experimental data reflects that the AST of GFGC columns was underestimated when the transverse pressure related to the lateral constraint of the concrete core was ignored. According to the comparative analysis, the current recommended model accurately depicts the AST of HFRGC columns. This is accomplished by calculating the restricting effect of glass-FRP stirrups and assuming that the ultimate stress and strain of HFRGC are the same as those of concrete constrained by glass-FRP. The ultimate AST of GFGC columns (P n ) was studied using a new experimental setup represented by Equation (32). This theoretical model is suggested based on the regression analysis and curve fitting tools in MATLAB based on the experimental results of the present study. P n = 0.85 f co A g − A FRP 1 + 3.1 f le f co 0.83 + ε co E FRP A FRP 1.85 + 7.46(ρ k ) 0.71 (ρ ε ) 1.17 (32) In this relation, A g is the whole area of the sample in the cross-section, A FRP is the overall area of glass-FRP long rebars in the cross-section, and E FRP is Young's modulus of glass-FRP long rebars. For the SFGC and GFGC concentric specimens, the suggested model's assumptions were equated to those of proposed models, experiments, and FEM predictions. Figure 17 represents the comparison of testing results, theoretical estimates, and FEM estimates for the ultimate AST of specimens. The theoretical models proposed by some previous studies [16,20,[66][67][68][69][70] were used to assess and compare the accuracy of the predictions of the presently proposed models.

Proposed AST Model
A series of simulations for measuring the compressive strength of FRP-RC columns has been recommended in previous research [16,20,[66][67][68][69][70]. Most of these studies did not include the influence of glass-FRP resistant pressure, except for [66,67] which offered models for glass-FRP-reinforced conventional concrete specimens. The ultimate longitudinal strain of concrete is assumed to be equivalent to the compressive strain of glass-FRP rebars in a few of the strength models [66,68]. The experimental data reflects that the AST of GFGC columns was underestimated when the transverse pressure related to the lateral constraint of the concrete core was ignored. According to the comparative analysis, the current recommended model accurately depicts the AST of HFRGC columns. This is accomplished by calculating the restricting effect of glass-FRP stirrups and assuming that the ultimate stress and strain of HFRGC are the same as those of concrete constrained by glass-FRP. The ultimate AST of GFGC columns ( ) was studied using a new experimental setup represented by Equation (32). This theoretical model is suggested based on the regression analysis and curve fitting tools in MATLAB based on the experimental results of the present study.
In this relation, is the whole area of the sample in the cross-section, is the overall area of glass-FRP long rebars in the cross-section, and is Young's modulus of glass-FRP long rebars. For the SFGC and GFGC concentric specimens, the suggested model's assumptions were equated to those of proposed models, experiments, and FEM predictions. Figure 17 represents the comparison of testing results, theoretical estimates, and FEM estimates for the ultimate AST of specimens. The theoretical models proposed by some previous studies [16,20,[66][67][68][69][70] were used to assess and compare the accuracy of the predictions of the presently proposed models.  The current recommended model accurately reflected the AST of HFRGC-columns, according to the comparative analysis. The average differences of the estimates of the equations given by Mohamed et al. [20], Afifi et al. [16], Pantelides et al. [66], Khan et al. [67], Tobbi et al. [68], Samani and Attard [69], and Hadhood et al. [70] from the testing measure- Data Availability Statement: Some or all data will be available upon request from the corresponding author.